У меня есть 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 составляет 1 м, а разрешение dsm составляет 0,5 м.
Я хочу рассчитать модель высоты кроны (CHM).
Самый простой способ это
CHM = dsm - dtm
Но когда я пытаюсь в R, появляется следующий код ошибки:
Error in compareRaster(e1, e2, extent = FALSE, rowcol = FALSE, crs = TRUE, :
different resolution
Есть ли простой способ игнорировать разрешение? Или я должен сделать передискретизацию данных перед дальнейшими вычислениями?
В ArcGis вы можете легко выполнить такой тип растрового расчета, потому что вам не нужно сначала повторять выборку данных.
Любые предложения будут оценены!