MODIS NDVI/EVI指数下载及处理

瞬间的可爱

已于 2024-04-29 09:13:35 修改

阅读量5k

点赞数 25

分类专栏： modis 文章标签： python arcgis

于 2024-04-28 17:15:21 首次发布

本文链接：https://blog.csdn.net/m0_70398377/article/details/138278921

版权

modis 专栏收录该内容

1 篇文章

订阅专栏

纯小白，第一次接触到modis数据，直接整不会了。花了几天恶补了各种知识，整理了一下modis数据下载及处理的过程，做个备份，以防忘记。

一.了解modis遥感产品数据

遥感影像（Remote Sensing Image）是指记录各种地物电磁波大小的胶片或照片，主要分为航空影像和卫星相片。这里我们主要介绍卫星影像。卫星影像，顾名思义，在卫星上通过传感器发射电磁波并由接收器接收地面信息的影像。卫星遥感影像有Landsat系列、Modis系列、Sentinel系列、哨兵系列等遥感产品。这里，modis数据，并不是以卫星的名字命名的，modis是它所在的卫星上的传感器。
modis,全称中分辨率成像光谱仪（Moderate-resolution Imaging Spectroradiometer，缩写MODIS），是美国航空航天局研制的大型空间遥感仪器。modis仪器分别位于两颗卫星上，Terra和Aqua。

Terra & Aqua 中分辨率成像光谱仪（MODIS）

第一台MODIS仪器于1999年在Terra上发射，其赤道穿越时间为上午10：30;第二台MODIS仪器于2002年在Aqua平台上发射，其穿越赤道时间为下午1：30。Terra-和Aqua-MODIS仪器每1到2天观察一次整个地球表面，获取36个光谱波段的数据，波长从0.4 μm到14.4 μm不等。双MODIS设计旨在优化无云成像，同时最大限度地减少早晨和下午阳光产生的阴影和眩光的光学效应。

MODIS数据产品具有三种不同的分辨率（250米、500米和1公里），有助于提高我们对陆地、海洋和低层大气中发生的全球环境过程和动态的理解。MODIS数据记录有助于扩展遗产数据，例如NOAA的先进甚高分辨率辐射计，从而确保此类收藏的关键连续性，以支持我们对短期和长期全球环境变化的调查。MODIS衍生的数据产品在帮助开发和验证具有足够预测潜力的全球地球系统模型方面继续发挥着至关重要的作用，这些模型可以在政策制定者应对全球环境变化时为他们提供信息和帮助。【翻译自官网】

二.下载modis NDVI/EVI数据

1.注册earthdata的账号（第一步貌似要连外网才能成功）

2.登录官网，下载数据。

细节请参考这位博主的博客。

数据下载| LAADS DAAC网站下载植被指数数据 - 知乎 (zhihu.com)

三.理解下载的modis数据格式

下载后的数据格式是hdf。hdf，一种复杂的通用数据格式，比netcdf格式复杂。

问题来了，这些文件名是什么意思？例如第一个文件MOD13A3.A2000032.h24v04.061.2020048131447.hdf是什么意思，完全看不懂。这里，modis产品是有它自己的一套命名规则的。MOD13A3是modis的NDVI/EVI产品缩写的意思，即L3级的每月1km的SIN（正弦等积投影）的MODIS/Terra植被指数产品。

A2000032，A表示Terra卫星，2000032表示2000年第32天。

h24v04表示下载的数据在地球表面的行列号。h，表示水平方向，v,表示垂直方向。

详情请参考http://t.csdnimg.cn/dPREF

061表示产品序号。

2020048131447表示数据加工完成时间，即2020年第48天，13:14:47完成产品的加工。

然后，hdf格式还是一头雾水，应该怎么查看呢？笔者最初使用matlab的hdfinfo，hdfread函数读取这些hdf文件的，但是不太方便。这里使用HDFView软件查看hdf数据比较方便。

还是很混乱，关于modis NDVI/EVI 的hdf格式。关于具体的细节，这里我们需要查看官方文档。在modis官网植被指数那部分内容处下载即可。

详情见文档P21-P23部分内容。

关于modis植被指数，官网提供6种不同数据集。

这里，我下载的是MOD13A3，每月1km的产品。

moids产品一般是SIN（正弦等积）投影。下载是每个square tile下载的，每片每片下载的，也就是前面说的h,v序列号。

对于文件格式，是hdf格式存储，存储内容包括两部分，科学数据集和元数据集。科学数据集，就是实际的矩阵数据集，元数据集描述的是数据集产品的各种信息。hdf存储类似于windows的文件夹形式，每个子数据集分开存储。如下图，NDVI和EVI分开存储。

细节涉及太多，具体内容请查看文档。

四、处理下载的hdf文件

那怎么处理呢？网上查了一下，非常混乱，不知从何开始。还有个问题，下载的hdf数据，需要做哪些预处理？到底需不需要几何校正、辐射校正、大气校正这些步骤，一头雾水。还是只需要重投影、拼接、裁剪？后来查了一些资料，遥感影像预处理步骤是包括几何校正、辐射校正、大气校正这些步骤的，但是具体需不需要对下载的遥感影像数据处理，要看下载的产品说明，这涉及到一个产品级别的问题。MODIS VI产品是经过很多处理的，因此只需要重投影即可。

详情请参考

MODIS产品简介_中国科学院遥感与数字地球研究所-高光谱遥感应用技术研究室 (hrs-cas.com)

本来是想用matlab处理的，但是网上基本找不到matlab处理的相关代码，基本上都是python的代码，看起来非常复杂。最后决定用modis官方提供的配套的工具MRT处理，貌似非常方便。

下载、安装及使用MRT

详情请参考这位博主的博客http://t.csdnimg.cn/eIAaF

后续有时间再更。

%%%%%%%%%%%%

更新

先尝试用批处理脚本处理一个hdf数据。

1.双击ModisTool.bat，进入GUI界面

2.通过Open Input File导入数据

3.Selected Bands,这里我只需要NDVI

4.Spatial Subset保持默认

5.指定输出文件位置

6.输出类型为GEOTIFF，Resampling Type为最邻近，投影为UTM，并设置投影参数，输出像元大小

7.不要点Run，点Save Parameter File,获得参数文件

8.在D:\MRT_tool\MRT\bin（我是这么安装的）新建文本文件test01.txt,写入resample -p C:\Users\hydrology\Desktop\demo\MOD13A3.A2000.prm 注意，这里-p后面没有加=，虽然MRT使用手册上是加了=，但是我运行就报错。

然后把特test01.txt后缀改为.bat

9.打开cmd,进入D:\MRT_tool\MRT\bin，执行test01.bat,可能会报错，出现MRT Open Datum File Error的错误，去MRT安装目录data文件夹下把spheroid. txt和datum. txt 放到D:\MRT_tool\MRT\bin下或许会解决，但可能还是会报错，我最后通过重启电脑解决了这个问题。