Растровая проекция до 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 * в качестве шаблона для проецирования. Таким образом, вы можете установить экстент и разрешение так, чтобы выходной растр совпадал с другими используемыми вами растровыми данными.