I would use setValues on individual layer.

Cheers,
Roman



On Wed, May 25, 2011 at 8:25 AM, Pierre Roudier <pierre.roudier@gmail.com>wrote:

> Hi Agustin,
>
> What about this:
>
> # Create a Raster object
> logo <- raster(system.file("external/rlogo.grd", package="raster"))
> # Use this rastre to create a new, void Raster object
> r <- raster(logo)
> show(r)
>
> Pierre
>
> 2011/5/25 Agustin Lobo <alobolistas@gmail.com>:
> > HI!
> >
> > I must create a layer with all missing values. I try:
> >> show(D10N2009)
> > class       : RasterBrick
> > dimensions  : 393, 606, 35  (nrow, ncol, nlayers)
> > resolution  : 0.008928571, 0.008928571  (x, y)
> > extent      : -2.004464, 3.40625, 40.49554, 44.00446  (xmin, xmax, ymin,
> ymax)
> > projection  : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
> +towgs84=0,0,0
> > values      :
> /media/Iomega_HDD/FLUXPYR/VGTFLUXPYR1/D10/FLXP_NDD10_2009.vrt
> > min values  : 0 0 0 0 0 0 0 0 0 0 ...
> > max values  : 255 255 255 255 255 255 255 255 255 255 ...
> >
> >> delme = D10N2009[[1]]*0
> >> show(delme)
> > class       : RasterLayer
> > dimensions  : 393, 606, 238158  (nrow, ncol, ncell)
> > resolution  : 0.008928571, 0.008928571  (x, y)
> > extent      : -2.004464, 3.40625, 40.49554, 44.00446  (xmin, xmax, ymin,
> ymax)
> > projection  : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
> +towgs84=0,0,0
> > values      : in memory
> > min value   : 0
> > max value   : 0
> >
> > But
> >> delme = D10N2009[[1]]*NA
> > Error in D10N2009[[1]] * NA : non-numeric argument to binary operator
> >
> > I also try
> >> NAvalue(delme) <- 0
> > but no trace of the NAs:
> >> show(delme)
> > class       : RasterLayer
> > dimensions  : 393, 606, 238158  (nrow, ncol, ncell)
> > resolution  : 0.008928571, 0.008928571  (x, y)
> > extent      : -2.004464, 3.40625, 40.49554, 44.00446  (xmin, xmax, ymin,
> ymax)
> > projection  : +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
> +towgs84=0,0,0
> > values      : in memory
> > min value   : 0
> > max value   : 0
> >
> > By now I write with NAflag=0 and read back, but there must be a
> > cleaner way of generating a layer with all NA values
> > (it's a missing layer for a particular time)
> >
> > Thanks!
> >
> > Agus
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
>
>
>
> --
> Scientist
> Landcare Research, New Zealand
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>



-- 
In God we trust, all others bring data.

	[[alternative HTML version deleted]]

