Необходимо перепроецировать растр перед мозаикой?

У меня есть список растров, которые я хотел бы сложить мозаикой. Проекция такая же, за исключением зоны 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]]) 
    }

Правильно ли я предполагаю, что сначала мне нужно перепроектировать? Спасибо за любое направление.


person Michael Davidson    schedule 26.03.2015    source источник
comment
Этот вопрос, вероятно, лучше подходит для географических информационных систем, но да, вы правы, полагая, что растры должны иметь общую CRS. Для повышения скорости вы можете рассмотреть возможность использования gdalwarp либо из командной строки, либо из пакета gdalUtils. См. Также эту ветку форума ESRI.   -  person jbaums    schedule 27.03.2015


Ответы (1)


help(mosaic) четко заявляет:

Все объекты должны иметь одинаковое начало, разрешение и систему отсчета координат.

Две разные зоны UTM не считаются имеющими одну и ту же систему отсчета (параметры различаются, если заглянуть под капот: ложное восточное положение, ложное северное положение, опорные меридианы и т. Д.).

Их следует перепроецировать на выходную CRS, которая должна быть эталонной системой, которая будет картографически действительной для большей площади, вместо того, чтобы продолжать использовать одну зону UTM, которая в конечном итоге будет показывать искажения и т. Д. При картировании больших областей.

person Benjamin    schedule 09.02.2017