我有一个具有不同分辨率的DTM和DSM。以下是每个栅格层的摘要。
> raster_dsm
class       : RasterLayer 
dimensions  : 2001, 2501, 5004501  (nrow, ncol, ncell)
resolution  : 0.5, 0.5  (x, y)
extent      : -112500.2, -111249.8, 388999.8, 390000.2  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=16.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +units=m +no_defs 
data source : D:/Test_Raster/DSM/dsm.asc 
names       : dsm 
>raster_dtm
class       : RasterLayer 
dimensions  : 1001, 1251, 1252251  (nrow, ncol, ncell)
resolution  : 1, 1  (x, y)
extent      : -112500.5, -111249.5, 388999.5, 390000.5  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=tmerc +lat_0=0 +lon_0=16.33333333333333 +k=1 +x_0=0 +y_0=-5000000 +ellps=bessel +units=m +no_defs 
data source : D:/Test_Raster/DTM/dtm.asc 
names       : dtm如你所见,dtm的分辨率是1m,dsm的分辨率是0.5m。
我想计算一个皇冠海格模型(CHM)。最简单的方法是让CHM = dtm
但是,当我在R中尝试时,会出现以下错误代码:
Error in compareRaster(e1, e2, extent = FALSE, rowcol = FALSE, crs = TRUE,  : 
  different resolution有什么简单的方法可以忽略决议吗?或者,在进一步计算之前,我必须对数据进行重新采样?在ArcGis中,您可以很容易地进行这种光栅计算,因为您不必先重采样数据。
如有任何建议,将不胜感激!
发布于 2019-05-08 04:27:26
是的,弧*会为你做这件事,但它实际上是做什么的?我认为最好避免这种模棱两可。在这种情况下,您不能使用dis/聚合,因为区段不同。所以你需要使用resample
https://stackoverflow.com/questions/56004124
复制相似问题