Hi Michelle,
The r-sig-geo list would be the ideal place for this question since
the authors of 'raster' would be likely to see it.
I haven't used the raster package before, but this seems to work:
#install.packages("raster",
repos="http://R-Forge.R-project.org")
r <- raster()
r <- setValues(r, sample(1:27, ncell(r), replace=1))
plot(r)
r[] <- r[]==1
plot(r)
If you'd rather have the non-1 cells be NA, just replace 'r[]==1'
with
'ifelse(r[]==1, 1, NA)'
hth,
Kingsford Jones
On Sat, Feb 21, 2009 at 8:46 AM, Michelle Greve
<michelle_greve at yahoo.com> wrote:> I have a raster (which I called glc), which I read
> into R as a raster with "raster.create.from.file" from the raster
package (R-forge). Values in glc range between 1 and 27 (whole numbers only).
I'd like to extract all cells with a value of 1 to create
> a new raster with only the cells that have a value of 1, or to extract
> the xy values of all raster cells with a value of 1. Could you please
> help me along? I have tried using several of the functions in the
> raster help package, e.g. "raster.get.xy.from.cell", but do not
know
> how to specify glc=1. I have tried to use e.g.:
> glc1<-raster.get.xy.from.cell(glc,glc==1)
> glc1<-raster.get.xy.from.cell(glc, glc[which(glc==1)])
>
> The reason I specifically used the raster package to read in my files is
because I have been unable to read my ascii files (which are very big) into R
using read.asciigrid because the ascii files are so big (and none of the various
suggestions I found in various R-help archives that I found for increasing
memory size were sufficient to allow me to read in my ascii files). Therefore if
alternatively somebody knows of another way to read raster or ascii files into
R, which would make it easier to conduct the extraction I mentioned in my first
paragraph, I would appreciate any suggestions.
>
> Regards,
> Michelle
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>