如何从栅格化的```SpatialPolygonsDataFrame```中获取单元尺寸？

Question

我在 R 中有 SpatialPolygonsDataFrame class，我想对其进行光栅化。我使用下面的代码栅格化多边形：

> p <- shapefile('MadaGranary')
> p
class       : SpatialPolygonsDataFrame 
features    : 1 
extent      : 100.1269, 100.5313, 5.793798, 6.446769  (xmin, xmax, ymin, ymax)
crs         : +proj=longlat +datum=WGS84 +no_defs 
variables   : 1
names       : id 
value       :  1 
> crs(p) <- NULL

> # Define RasterLayer object
> r.raster <- raster(extent(p), res = 100)
> #rasterize
> p.r <- rasterize(p, r.raster)
> print(p.r)
class      : RasterLayer 
dimensions : 1, 1, 1  (nrow, ncol, ncell)
resolution : 100, 100  (x, y)
extent     : 100.1269, 200.1269, -93.55323, 6.446769  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : layer 
values     : NA, NA  (min, max)
attributes :
 ID id
  1  1

不幸的是，这段代码已将单元格维度分配给 dimensions : 1, 1, 1 (nrow, ncol, ncell)，获取 ncols and nrows 的实际值对我来说非常重要。

而当我 print p.r 时，我需要在控制台中获取 dimensions : (nrow, ncol, ncell) 的值，以便我能够保存所有属性，以便稍后使用它进行其他分析。注意：我需要得到的resolution输出应该是表示新栅格中ncols and nrow的真实值，而不仅仅是dimensions : 1, 1, 1 (nrow, ncol, ncell)，如中所示我的代码。

谁能帮忙吗？？

Answer 1

如果你想设置 nrow 和 ncol 那么你就这样做（而不是设置分辨率，它表示每个单元格的大小）。这是一个例子：

library(raster)
#Loading required package: sp
e <- extent(0,1,0,1)
r <- raster(e, nrow=100, ncol=100)
r
#class      : RasterLayer 
#dimensions : 100, 100, 10000  (nrow, ncol, ncell)
#resolution : 0.01, 0.01  (x, y)
#extent     : 0, 1, 0, 1  (xmin, xmax, ymin, ymax)
#crs        : NA

如何从栅格化的```SpatialPolygonsDataFrame```中获取单元尺寸？

how to get cell dimensions from rasterized ```SpatialPolygonsDataFrame```?

r

raster

sp