Problem with raster
Mark,
You would not be tell from the error message, but the error occurs
because you are using an invalid CRS.
It should be:
proj4string=CRS("+proj=longlat")
instead of this:
proj4string=CRS("longlat")
Robert
On Fri, Feb 21, 2014 at 6:59 AM, Mark Payne <markpayneatwork at gmail.com> wrote:
Hi, I am having a problem with the following code - I'm trying to create a raster from a SpatialGrid, but there seems to be something going wrong along the way - the problem pops up when I try to build two rasterlayers into a brick:
halfrect.grd <- GridTopology(c(-22.75,36.75),
+ c(0.5,0.5), + c(46,55))
grd.sp <- SpatialGrid(halfrect.grd,proj4string=CRS("+longlat"))
grd.r <- raster(grd.sp)
brick(grd.r,grd.r)
class : RasterBrick dimensions : 55, 46, 2530, 1 (nrow, ncol, ncell, nlayers) resolution : 0.5, 0.5 (x, y) extent : -23, 0, 36.5, 64 (xmin, xmax, ymin, ymax) coord. ref. : +longlat This works ok so far. But then if I try to change the value, it fails
EP.r <- grd.r EP.r[] <- NA brick(EP.r,EP.r) #Doesn't
Error in x[, 2] : subscript out of bounds
values(EP.r) <- NA brick(EP.r,EP.r) #Doesn't
Error in x[, 2] : subscript out of bounds
My package details are below. Does anyone have any ideas what could be going wrong? Mark
packageDescription("raster")
Package: raster Type: Package Title: raster: Geographic data analysis and modeling Version: 2.2-12 Date: 2014-1-17 <snip>
R.version
_ platform x86_64-pc-linux-gnu arch x86_64 os linux-gnu system x86_64, linux-gnu status major 3 minor 0.1 year 2013 month 05 day 16 svn rev 62743 language R version.string R version 3.0.1 (2013-05-16) nickname Good Sport
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo