r - 使用 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 m，我在 R 脚本中错误地写了“pixel_size=500”。这给出了一个 500 度乘 500 度的单个像素作为输出。当投影类型为地理 (proj_type="GEO") 时，像素大小应以度为单位。对于本例中 MODIS 瓦片的纬度（瓦片 h09v06），0.004514 度约为 500 m。因此，像素大小在脚本中应该写成：pixel_size=0.004514。