Растровая проекция до UTM в Lat-Lon в R

Я новичок в R. У меня есть rasterlayer и я хочу изменить его систему координат (с UTM на lan-lot). Я использовал следующий код:

r <-projectRaster(from = tiff,
   crs=crs("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"))

(Здесь TIFF rasterlayer перепроецироваться)

Но я получаю ошибку

"Ошибка: невозможно выделить вектор размером 439,0 Мб"

Кто-то, пожалуйста, предоставьте решение

1 ответ

Это странно. Вы прочитали все значения в tiff в память?

Проблема может уйти

  • используя аргумент имени файла в projectraster. Например projectRaster(tiff, filename="out.tif", crs= ...)

  • установка некоторых rasterOptions на более низкие значения, например, rasterOptions(chunksize=1e+06, maxmemory=1e+08)

Кроме того, это, как правило, не очень хорошая идея, чтобы использовать projectRaster с crs= аргумент. Лучше предоставить объект Raster * в качестве шаблона для проецирования. Таким образом, вы можете установить экстент и разрешение так, чтобы выходной растр совпадал с другими используемыми вами растровыми данными.

Другие вопросы по тегам