Репроекция в R с использованием инструмента MRT

Question

Я пытаюсь загрузить и перепроектировать данные MODIS в GEO (Geographic) проекцию, используя код ниже. Тем не менее, есть проблема с результатом. Размеры выходного файла составляют 1 на 1 пиксель (размеры: 1, 1, 1 (nrow, ncol, ncell) вместо 4000 на 2200.

x <- "MOD09A1"
ModisDownload(x=x,h=c(9),v=c(6),
          dates=c("2011.01.01","2011.01.01"),
          MRTpath="~\\MRT_download_Win\\bin",
          proj=TRUE,
          proj_type="GEO",
          datum="WGS84",
          resample_type="NEAREST_NEIGHBOR",
          proj_params="0 0 0 0 0 0 0 0 0 0 0 0 0 0 0",
          pixel_size=500)

Answer 1

Источником ошибки является размер пикселя. Поскольку размер пикселя MODIS составляет 500 м, я неправильно написал "pixel_size=500" в сценарии R. Это дало один пиксель 500 градусов на 500 градусов в качестве вывода. Когда тип проекции - Географический (proj_type="GEO"), размер пикселя следует указывать в градусах. Для широты плитки MODIS в этом примере (плитка h09v06) 0,004514 градусов составляют приблизительно 500 м. Следовательно, размер пикселя должен быть записан в скрипте как: pixel_size=0.004514.