siminair的个人博客分享 http://blog.sciencenet.cn/u/siminair

博文

MODIS单站点(小区域)数据的获取

已有 2296 次阅读 2020-11-11 17:46 |系统分类:科研笔记

     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各种数据产品就不用愁了。




https://m.sciencenet.cn/blog-3409822-1257983.html

上一篇:文章修改稿中如何将修订的地方批量自动标红
下一篇:参加大气环境遥感与协同分析青年学者论坛-线下会议

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2023-1-28 21:44

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部