Растровая проекция до UTM в Lat-Lon в R - PullRequest
0 голосов
/ 02 июля 2018

Я новичок в 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 Ответ

0 голосов
/ 03 июля 2018

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

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

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

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

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

...