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

Question

我正在尝试使用下面的代码将MODIS数据下载并重新投影到GEO（地理）投影中。但是，结果存在问题。输出文件尺寸仅为1 x 1像素（尺寸：1,1,1（nrow，ncol，ncell），而不是4000 x 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米，我错误地写了＆＃34; pixel_size = 500＆＃34;在R脚本中。这给出了一个500度，500度的单个像素作为输出。当投影类型为地理（proj_type =＆＃34; GEO＆＃34;）时，像素大小应以度为单位。对于此示例中的MODIS瓦片的纬度（瓦片h09v06），0.004514度约为500米。因此，像素大小应该在脚本中写为：pixel_size = 0.004514。