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

使用 MRT 工具在 R 中进行重投影

Reprojection in R using MRT tool

r

raster

reprojection-error