1
我無法使用柵格包河返回層不會柵格化中的R
shp<-shapefile(ZoneShape);
lcRas<-raster(lcRaster);
r<-raster(ncol=ncol(lcRas), nrow=nrow(lcRas), crs=CRS);
res(r)<-res(lcRas);
extent(r)<-extent(lcRas);
>r
class : RasterLayer
dimensions : 22610, 27959, 632152990 (nrow, ncol, ncell)
resolution : 1, 1 (x, y)
extent : 554739, 582698, 3837197, 3859807 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=15 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0
shp$GID<-1:nrow(shp);
> shp
class : SpatialPolygonsDataFrame
nfeatures : 1
extent : 554838, 582597.6, 3837297, 3859707 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=15 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
nvariables : 14
names : SP_ID, NAME, Shape_Leng, Shape_Area, GID
>zoneRas<-rasterize(shp, r, "GID")
以下錯誤光柵化shape文件使用光柵包:
trying to get slot "coords" from an object of a basic class ("NULL") with no slots.
任何人都可以看看我在想什麼/搞砸了嗎?
感謝
這是一個數據問題,因爲我能夠與另一個shapefile一起運行。但是,shapefile顯示正確,並且Arcmap中的屬性表看起來不錯。問題也在[貼在GISStackExchange上](http://gis.stackexchange.com/questions/64343/r-raster-package-rasterize-failingno-slots)。 – RichT