在R中捕捉栅格网格

Question

我正在尝试在R中对齐两个栅格网格。一旦对齐，我希望能够将它们一起添加。

我试图检查堆叠是否有效：

grid_snap <- stack(GLC2000_sdw, afriPop_sdw)

我收到以下错误：

  

compareRaster（x）出错：程度不同

栅格网格具有以下属性：

show(habi_sdw)
# class       : RasterLayer 
# dimensions  : 9187, 9717, 89270079  (nrow, ncol, ncell)
# resolution  : 0.00892857, 0.00892857  (x, y)
# extent      : -28.83706, 57.92186, -36.02464, 46.00214  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs 
# data source : C:\Users\di39\AppData\Local\Temp\R_raster_di39\raster_tmp_2015-08-12_172902_12860_17067.grd 
# names       : layer 
# values      : 0, 333707.6  (min, max)

show(Pop_sdw)
# class       : RasterLayer 
# dimensions  : 10143, 8858, 89846694  (nrow, ncol, ncell)
# resolution  : 0.008333333, 0.008333333  (x, y)
# extent      : -17.53524, 56.28143, -46.97893, 37.54607  (xmin, xmax, ymin, ymax)
# coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
# data source : C:\Users\di39\AppData\Local\Temp\R_raster_di39\raster_tmp_2015-08-12_170421_12860_12760.grd 
# names       : pop2010ppp 
# values      : 0, 128925.9  (min, max)

在alignExtent()包中使用raster似乎不是正确的方法。

我是否需要重新采样，因为分辨率略有不同？

（0.00892857 x 0.00892857）vs（0.008333333 vs 0.008333333）

Answer 1

首先使用resample(GLC2000_sdw, afriPop_sdw, method="ngb")然后使用crop (GLC2000_sdw, afriPop_sdw)以确保它们具有相同的范围。