MODIS是目前生态学研究中很重要的一个遥感数据源，但数据下载和后期分析数据还是需要一定的遥感基础。由于官网正常的MODIS数据下载，是以景为单位，一景MODIS遥感数据的覆盖范围为10°*10°，约100km*100km公里。如果想获得某一个台站的MODIS数据产品（如GPP），就需要利用GIS分析工具根据站点经纬度来提取数值，过程稍微繁琐，特别是对于没有地学基础的研究者而言，还是有点小难度的。

      其实利用官网提供的R包，就可以实现MODIS数据的下载，很多年前官网有对应MATLAB的程序包，只是后来不知什么原因只保留R包。当然R包的使用，稍微要求一点点编程的基础，需要安装R。好在现在用R的人挺多，应该问题不大。

https://daac.ornl.gov/MODIS/MODIS-menu/modis_webservice.html 有关于 MODISTools的介绍

https://github.com/ropensci/MODISTools   是程序包的具体介绍和使用

具体步骤如下：
首先在R环境下安装官网提供的程序包-MODISTools

在R编译器中输入一下代码

1. install.packages("MODISTools")  
2. library("MODISTools")
查看MODISTools里的主要函数及使用说明
3. ？？ MODISTools  #打开MODISTools的帮助界面
4. products <- mt_products()  # 是获取MODIS的数据产品，存于products变量
   list(products)  #查看MODIS所有数据产品的具体信息，目前程序包支持的是32个产品
5. bands <- mt_bands(product = "MOD17A2H")  # 了解需获取产品包括的数据项，根据products的产品选择，注意R严格区分大小写
6. 根据需要，设置具体参数，下载站点的MODIS数据下载
subset <- mt_subset(product = "MOD11A2",  #产品编号
                    lat = 40, #纬度
                    lon = -110, #纬度
                    band = "LST_Day_1km", #产品里对应的数据项
                    start = "2004-01-01",  #开始时间
                    end = "2004-02-01",  #结束时间
                    km_lr = 1,  #上下缓冲距离单元，如果只需要站点本身，可将该值设为0
                    km_ab = 1,  #左右缓冲距离单元，与数据产品本身的空间分辨率有关
                    site_name = "testsite", #名称
                    internal = TRUE,
                    progress = FALSE)
 7. 数据的查看
    list(subset)    
 8. 数据的导出   
    write.table (subset,file='text1.csv',sep =",") # 将数据存储到工作路径下,以csv文件保存，以逗号分隔，可以通过EXCEL直接打开
    也可在file中添加具体的路径信息，如果不知道路径信息的添加，可以使用getwd()获取文件存储的位置，这样就能获得单站点的MODIS数据了。
        
 PS：更多详细的内容，可以查看官网或查看程序包。
 总体流程不算复杂，掌握了这一招，获取小范围的MODIS各种数据产品就不用愁了。