0

有没有人有使用 MODIStp 包下载 MODIS 数据的经验？

我想以 1 公里的分辨率下载中国的 EVI 数据并按照教程https://cran.r-project.org/web/packages/MODIStsp/vignettes/MODIStsp.html创建了一个 jason 文件（最后的问题）并这样做。它运行正常，然后生成错误

  library(MODIStsp)

  MODIStsp(gui = FALSE,
           opts_file = file.path(dir_ls$base_path,'modis_spatial_file.json'))

  GDAL version in use: 3.0.4
  [Tue Mar 30 11:29:07 2021] MODIStsp --> Starting processing
  [Tue Mar 30 11:29:07 2021] Accessing http server at:  https://e4ftl01.cr.usgs.gov/MOLT/MOD13A2.006/
  [Tue Mar 30 11:29:09 2021] Retrieving list of available ` Terra ` Files for Year 2000
  [Tue Mar 30 11:29:11 2021] Downloading Terra Files for date: 2000_02_18 : 1 of:  53
  [Tue Mar 30 11:29:17 2021] Downloading Terra Files for date: 2000_02_18 : 2 of:  53
  .
  .
  .
  [Tue Mar 30 11:34:13 2021] Downloading Terra Files for date: 2000_02_18 : 52 of:  53
  [Tue Mar 30 11:34:15 2021] Downloading Terra Files for date: 2000_02_18 : 53 of:  53
  [Tue Mar 30 11:34:15 2021] [Tue Mar 30 11:34:15 2021] 53 files for date: 2000.02.18 were 
  successfully downloaded!
  [Tue Mar 30 11:34:15 2021] Processing and Reprojecting Terra EVI files for date: 2000_02_18
   Error in r[i1] - r[-length(r):-(length(r) - lag + 1L)] : 
   non-numeric argument to binary operator

这个有点神秘的错误让我感到困惑，想知道是否有人知道发生了什么？我的杰森文件如下：

{
    "selcat": "Ecosystem Variables - Vegetation Indices",
    "selprod": "Vegetation_Indexes_16Days_1Km (M*D13A2)",
     "prod_version": "6",
     "sensor": "Terra",
     "bandsel": "EVI",
     "quality_bandsel": null,
     "indexes_bandsel": null,
     "download_server": "http",
     "user": "XXXX",
     "password": "XXXX",
     "downloader": "http",
     "download_range": "Full",
     "start_date": "2000.02.18",
     "end_date": "2021.03.30",
     "spatmeth": "map",
     "start_x": 21,
     "end_x": 31,
     "start_y": 3,
     "end_y": 8,
     "bbox": [66.2891, 4.3513, 137.832, 55.6556],
     "spafile": null,
     "drawnext": "{\n  \"type\": \"FeatureCollection\",\n  \"features\": [\n    {\n      \"type\": 
     \"Feature\",\n      \"properties\": {\n        \"X_leaflet_id\": 3991,\n        \"layerId\": 
     \"3991\",\n        \"edit_id\": \"3991\"\n      },\n      \"geometry\": {\n        \"type\": 
     \"Polygon\",\n        \"coordinates\": [\n          [\n            [66.289051, 4.351296],\n            
     [66.289051, 55.65556],\n            [137.832014, 55.65556],\n            [137.832014, 
     4.351296],\n            [66.289051, 4.351296]\n          ]\n        ]\n      }\n    }\n  ]\n}",
  "out_projsel": "User Defined",
  "output_proj": "4326",
  "out_res_sel": "Native",
  "out_res": 926.6254,
  "resampling": "near",
  "reprocess": false,
  "delete_hdf": true,
  "nodata_change": false,
  "scale_val": true,
  "out_format": "GTiff",
  "compress": "LZW",
  "out_folder": "D:/modis",
  "out_folder_mod": "D:/modis_hdf",
  "MODIStspVersion": "2.0.5"
   }
4

0 回答 0