用dtm和dsm计算不同分辨率的CHM

Question

我有一个具有不同分辨率的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中，您可以很容易地进行这种光栅计算，因为您不必先重采样数据。

如有任何建议，将不胜感激！

Answer 1

是的，弧*会为你做这件事，但它实际上是做什么的？我认为最好避免这种模棱两可。在这种情况下，您不能使用dis/聚合，因为区段不同。所以你需要使用resample