Modis积雪遥感数据处理（积雪覆盖度）【时间：2001.12.01--2022.03.31】

文章以详细的方式说明在处理Modis积雪遥感数据时的 操作方法与步骤，本文仅供参考学习，不可用于商业化，谢谢。

Modis影像数据介绍

要通过Modis数据进行遥感数据的批量下载与提取，首先是需要了解一下什么是Modis数据。

EOS-MODIS 卫星

美国航空航天局（NASA）于 1999 年 12 月 18 日成功发射了地球观测系统
EOS 的第一颗先进的极地轨道环境遥感卫星 TERRA（EOS-AM1）。这颗卫星是美国国家航空航天局（NASA）地球行星任务中总数 15 颗卫星的第一颗，也是第一个提供对地球过程进行整体观测的系统。由于 TERRA 卫星每日地方时上午 10：30 时过境，因此也把它称作地球观测第一颗上午星。TERRA 卫星上共载有五个对地观测传感器，其中中分辨率成像光谱仪MODIS 是主要遥感探测仪器之一。2002 年 5 月 4 日，NASA 发射了载有另一个 MODIS 传感器的 AQUA 卫星，不过AQUA 将在地方时下午过境，这样 TERRA 与 AQUA 上的 MODIS 数据在时间更新频率上相配合，加上晚间过境数据，对于接收 MODIS 数据来说，每天可获得 2 次白天和 2 次黑夜更新的数据，这些数据的较快更新速率，对于实时地球观测（如自然灾害监测）和长期的全球系统变化的研究都有非常重要的服务价值。

MODIS 数据

目前在世界上应用广泛的中低分辨率遥感卫星数据是 NOAA-AVHRR 和
MODIS 数据。相比于 AVHRR 数据，MODIS 数据在保留了 AVHRR 的功能的同时，在数据分辨率、数据波段数及数据应用范围都有较大的改进。MODIS 数据波段范围广，比 SPOT、TM 和 AVHRR 有更高的光谱分辨率，达到 36 个波段，其中波段 1～19 和波段 26 在可见光和近红外波段，其它波段在 3～14 μm的热红外波段。MODIS 数据 1、2 波段的空间分辨率为 250m，3～7 为 500m，其余的为1000m，扫描宽度为 2330km，而 AVHRR 数据的分辨率是 1100m。但由于MODIS 数据在波段数和分辨率方面的改进，使得 MODIS 数据量增加很多。

MODIS 数据产品

NASA 根据规范和数据分级标准对 MODIS 数据按照一定空间分辨率和时间
分辨率进行连续开发了标准数据产品。MODIS 标准数据产品根据内容的不同分为0 级（未经处理的、包括全部数据信息在内的原始数据）、1 级数据产品（已赋予定标参数），在 1B 级数据产品之后，划分 2-4 级数据产品，包括：陆地标准数据产品、大气标准数据产品和海洋标准数据产品等三种主要标准数据产品类型，总计分解为 44 种标准数据产品类型（数据产品类型官方文档介绍网站）。

其他更详细的内容可以到MODIS数据介绍及下载查看

Modis数据的选择与下载方式

首先进入官网：数据下载网站链接: https://ladsweb.modaps.eosdis.nasa.gov/search/

这里可以直接搜索相对应产品的Modis数据。
可是！因为原本想处理的是MOD09A1系列的数据，奈何20年的数据量实在是太大太大了！！！！！受不了，一年限制12月份到3月份，大概120天左右，安徽省又需要三张瓦片合成，所以光一年可能就要下载360次左右，20年大概就是6000多条......
所以后来查了一下，针对积雪遥感的数据，NASA官网有现成的，是以八天为一个周期进行记录。

安徽位于东经114°54′一119°37′与北纬29°41′— 34°38′之间
从赤道开始是0°
SW（南-西）：填写初始的南纬和西经，也就是29°41′,114°54′；
NE（北-东）：填写初始的南纬和西经，也就是34°38′,119°37′。

然后再滑动滚轮到下面选择数据的开始时间和结束时间。

选择好数据之后点击+号下载数据：


警告：最好还是不要用迅雷批量下载，因为链接点进去之后需要登录才能继续进行下载任务，但是迅雷直接粘贴了链接之后链接不上下载资源，下载下来的文件完全不完整，建议还是手点下载，或者考虑IDM下载器或者DownThemAll插件，方法只能自行百度了。

上面是我已经下载好的20年的数据，范围是安徽省的，未裁剪的数据。

HDF转TIF方式（MRT工具使用）【批量处理转换等】

首先，MRT工具的安装我就不过多做详细介绍了，安装的时候注意环境变量的配置就行了。
数据准备好之后，我们先打开我们的MRT工具

1. 首先导入下载好的数据。（后缀为hdf格式的数据文件）
   这里是对三个波段瓦片进行拼接的操作。
   参考坐标系：WGS84；
   
   OK之后出来点击Run。
   
   等待加载完成即可在对应文件夹产生tif文件。
   
   这里是未裁剪的对应影像。
2. 批量处理拼接hdr数据。
   首先我们在界面里面已经设置好相对应的格式了：
   注意注意一定要保存到要处理数据的文件夹里，不然cmd命令无法直接调用。
   可以看到prm文件已经和数据在同一个目录下了。
   然后我们win+r打开运行，输入cmd调出任务窗口，依次输入命令：

（1）cd/d G:\MRT\bin
进入MRT的安装位置的bin目录下面；

（2）java -jar MRTBatch.jar -d G:\Snow\HDF\2022-2023 -p G:\Snow\HDF\2022-2023\2022-2023.prm -o G:\Snow\HDF\2022-2023\result
-d到-p之间的路径是你存储数据的路径
-p到-o之间的路径是你保存prm文件的路径
-o之后的路径是处理好之后的数据的保存路径【数据输出路径】

（3）MRTBatch.bat
然后输入命令调用工具即可。

现在就是在自动处理了，可能时间会有些长，大家可以稍作等待。

现在这是已经完成了。我们可以打开result文件夹查看一下：

27个项目，tif文件个数是9个，说明把2022-2023年的数据拼接成功了。

3. 我们用ArcMap软件打开我们的tif文件查看一下是否是正常的：
   
   可以正常打开，导入shp文件看看坐标位置是否正常即可。

积雪TIF数据处理（ArcMap处理）

针对tif文件进行栅格裁剪、重采样、符号系统颜色分级
在符号分级时遇到这样的问题，需要对栅格数据进行重采样



计算栅格面积可以用”以表格形式进行分区几何统计“，输入栅格和选择字段值，即可生成excel表格数据。