У меня есть список растров, которые я хотел бы сложить мозаикой. Проекция такая же, за исключением зоны utm. Вот координаты одного растрового стека в зоне 50 UTM:
coord. ref. : +proj=utm +zone=50 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
Вот еще один rasterStack с идентичной проекцией, за исключением зоны UTM 51:
coord. ref. : +proj=utm +zone=51 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
Когда я пытаюсь создать мозаику без предварительного перепроецирования (как это можно сделать в других ГИС), я получаю следующую ошибку:
rMosaic <- do.call(mosaic,rStacks)
Error in compareRaster(x, extent = FALSE, rowcol = FALSE, orig = TRUE,: different CRS
Это наводит меня на мысль, что мне нужно перепроецировать все растровые стеки. С projectRaster
это требует больших вычислительных затрат времени, но я считаю, что могу сделать это следующим образом:
for(i in 2:length(rStacks)){
rStacks[[i]] <- projectRaster(from=rStacks[[i]], to=rStacks[[1]])
}
Правильно ли я предполагаю, что сначала мне нужно перепроектировать? Спасибо за любое направление.
gdalwarp
либо из командной строки, либо из пакетаgdalUtils
. См. Также эту ветку форума ESRI. - person jbaums   schedule 27.03.2015