Растеризация пространственных полигонов в R дает растр со значениями NA

Question

У меня возникли проблемы с преобразованием этого пространственного полигонального кадра в растр. Когда я делаю преобразование, растр имеет NA в качестве значений. Как показано ниже:

 DL3
[1]
class       : SpatialPolygonsDataFrame 
features    : 126 
extent      : -15.04001, 46.1036, 3.759985, 31.71804  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +towgs84=0,0,0 +ellps=WGS84 
variables   : 1
names       :  LFRP 
min values  :    14 
max values  : 335.2 

Вот как я растеризую это:

##TO CONVERT TO RASTER
FunR<-function(r){
ext<-raster(extent(r))
crs(ext)<-crs(r)
D<-rasterize(r,ext,field=1,update=T)
D}

DL4<-lapply(DL3,FunR)
DL4
[1]
class       : RasterLayer 
dimensions  : 45, 40, 1800  (nrow, ncol, ncell)
resolution  : 1.52859, 0.6212901  (x, y)
extent      : -15.04001, 46.1036, 3.759985, 31.71804  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +towgs84=0,0,0 +ellps=WGS84 
data source : in memory
names       : layer 
values      : NA, NA  (min, max)

Что я могу делать неправильно? Мне нужна помощь с методом, чтобы гарантировать, что значения в кадре данных отражаются в растре, пожалуйста.

Answer 1

Принцип работает, как показано ниже.

library(raster)
SPP <- shapefile(system.file("external/lux.shp", package="raster"))
r <- raster(SPP, ncol=40, nrow=45)    
SPP2 <- rasterize(SPP, r, "ID_1")

Поскольку это не работает для вас, я предполагаю, что ваши полигоны очень малы относительно размера растровых ячеек, так что ни одна из ячеек не покрыта. Можете ли вы попробовать с гораздо меньшими ячейками сетки? Если действительно многоугольники настолько малы, возможно, имеет смысл использовать их центроиды (coordinates(SPP)) и растеризовать их (как точки).