Блог

Главная - Вопросы по программированию - Как синхронизировать растры с разными проекциями

Как синхронизировать растры с разными проекциями

#r #raster #map-projections

#r #растровые #карта-проекции

Вопрос:

Я хотел бы синхронизировать два файла geotiff с разной степенью и разной системой проекции. Я пытаюсь использовать spatial_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 Rpreqieslc_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 Rpreqieslc_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")

У меня была такая же проблема, и я решил дело. удачи

Метки: