使用MRT工具在R中重複投影

Question

我正嘗試使用下面的代碼下載並將MODIS數據重新投影到GEO（地理）投影中。但是，結果存在問題。輸出文件尺寸是1：1像素只（尺寸：1，1，1（nrow，NcoI位，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米，我在R腳本中錯誤地寫了「pixel_size = 500」，這是給一個500度乘500度的單個像素作爲輸出，當投影類型是Geographic（proj_type =「GEO」）時，像素大小應該用度數表示，對於本例中的MODIS瓦片的緯度（瓦片h09v06）0.004514度約爲500米，因此，像素尺寸應該寫在腳本中：pixel_size = 0.004514