如何用R语言批量下载MODIS13Q1数据

使用R语言批量下载MODIS13Q1数据可以通过以下步骤实现：

1. 安装和加载必要的R包：首先，确保你已经安装了`MODIS`和`raster`这两个R包。如果没有安装，可以使用以下命令进行安装：

install.packages("MODIS")
install.packages("raster")

然后，加载这两个包：

library(MODIS)
library(raster)

2. 设置MODIS下载参数：使用`MODISoptions()`函数来设置下载参数，包括数据产品、日期范围、下载路径等。例如，设置下载路径为`"C:/MODIS"`，数据产品为`"MOD13Q1"`，日期范围为从`"2019-01-01"`到`"2019-12-31"`：

MODISoptions(localArcPath = "C:/MODIS",
             outDirPath = "C:/MODIS",
             product = "MOD13Q1",
             begin = "2019-01-01",
             end = "2019-12-31")

3. 批量下载MODIS数据：使用`runGdal() `函数来批量下载MODIS数据。例如，下载NDVI数据：

runGdal(product = "MOD13Q1", 
        tileH = 18, 
        tileV = c(4, 5), 
        bands = "250m_16_days_EVI", 
        job = "MOD13Q1_download")

其中，`tileH`和`tileV`参数指定了要下载的地理瓦片，`bands`参数指定了要下载的数据集。

4. 查看下载结果：下载完成后，可以使用`listMosaicFiles()`函数来查看下载的文件列表：

listMosaicFiles(product = "MOD13Q1", 
                tileH = 18, 
                tileV = c(4, 5), 
                bands = "250m_16_days_EVI")