Как синхронизировать растры с разными проекциями
Я хотел бы синхронизировать два файла геотипов с разным экстентом и другой системой проекции. Я пытаюсь использовать Space_Sync_raster, но это не сработало. Как я должен сделать? Ниже приведены подробности файла.
class : RasterLayer
band : 1 (of 3 bands)
dimensions : 480, 480, 230400 (nrow, ncol, ncell)
resolution : 1, 1 (x, y)
extent : 0, 480, 0, 480 (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : D:\Test R\preqies\lc_2005.tif
names : lc_2005
values : 0, 255 (min, max)
class : RasterLayer
dimensions : 869, 1201, 1043669 (nrow, ncol, ncell)
resolution : 100, 100 (x, y)
extent : 771507.8, 891607.8, 9788643, 9875543 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=47 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
data source : D:\Test R\preqies\lc_2010.tif
names : lc_2010
values : 1, 22 (min, max)
Я получаю сообщение об ошибке при попытке этого кода
landuse2<-spatial_sync_raster(landuse2, landuse1, method = "ngb")
**Error in rawTransform(projfrom, projto, nrow(xy), xy[, 1], xy[, 2]) :
no arguments in initialization list**
Может кто-нибудь сказать мне, что часть, которая делает мой код не работает? Спасибо
1 ответ
Попробуйте добавить ту же проекцию:
projection(landuse2)<- "+proj=utm +zone=47 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84"
landuse2<-spatial_sync_raster(landuse2, landuse1, method = "ngb")
У меня была такая же проблема, и я решил дело. везение