概述

基于A省shp地图，获取MODIS相关数据，使用ArcMap提取出该省每个县的各土地利用类型的面积。

一、MODIS数据下载

NASA Earth Science Data Systems (ESDS)及相关介绍
https://earthdata.nasa.gov/

MODIS数据

数据介绍
https://blog.csdn.net/u010666165/article/details/81145006
1999年2月18日，美国发射了地球观测系统（EOS）的第一颗极地轨道环境遥感卫星Terra，2002年4月22日发射Aqua卫星。MODIS是搭载在TERRA和AQUA卫星上的传感器，全称为中分辨率成像光谱仪（moderate-resolution imaging spectroradiometer）。
下载地址
https://ladsweb.modaps.eosdis.nasa.gov/search/order
The Level-1 and Atmosphere Archive & Distribution System (LAADS) Distributed Active Archive Center (DAAC)
美国航空航天宇航局的戈达德航天中心用来存放数据的一个网站接口。
依照产品进行搜索，按照顺序选择传感器、时间、地点和数据。

产品选择

每一个产品都有简介，如选择MODIS Combined-Terra-Aqua 的 MCD12Q1产品。简介为：MODIS/Terra+Aqua Land Cover Type Yearly L3 Global 500m SIN Grid 土地覆盖年度产品，产品等级为L3，全球尺度500m格网。
在这里插入图片描述
详细信息：

注意: 共5种土地覆盖类型分类，选择前考虑该方法是否适合自己的专题研究。
MCD12Q1详细介绍:
https://blog.csdn.net/qq_36958801/article/details/108096958
命名规则：

时间选择跨度不能太大，否则检索时间会很长，因为是MCD12Q1是年度数据，所以时间跨度（Data Range）搜索某一年1月1号到年末最后一天即可，或（SingData）选择1月1号。年度数据一般为每年的第一天开始到第二年的第一天，即A2015001，A2016001…

下载步骤

https://blog.csdn.net/weixin_43513353/article/details/104456597
获取省级数据依照Tiles选择比较快。
在这里插入图片描述
中国相关地区的数据编号（来源如下）
https://blog.csdn.net/qq_37798011/article/details/80973188

注意: 下载前需要登录EarthData账号，注册需要gmail邮箱，手机QQ邮箱app可以注册。

二、MODIS数据处理

MRT

NASA提供MODIS Reprojection Tool（MRT），专门用于处理MODIS数据，官网提供了Linux、macos、windows的版本。

MRT下载与安装
https://blog.csdn.net/gisboygogogo/article/details/75784080
MRT数据提取
https://blog.csdn.net/gisboygogogo/article/details/79533595
MRT拼接、投影和格式转换
https://www.cnblogs.com/orchidLan/p/7210448.html

数据处理流程：下载与安装MRT→MRT数据拼接→ArcMap投影处理。
批量数据拼接分为一年之内数据和多年数据，本数据为年度数据，使用一年内数据拼接方法。

新建一个data文件夹，将".hdf"数据（即刚刚下载下来的数据）拷到该文件夹下。（注意备份，拼接后同一目录下的hdf数据将被删除）；
在DATA文件夹下新建一个Result文件夹用于存放拼接后数据;
将以下代码内容粘到一个txt文件中，更改txt文件内容，将其副本的扩展名改为bat，假如命名为MOSAIC.bat；
点击MOSAIC.bat，拼接完成。

set MRTDATADIR=D:\Modis\MRT\data
set /a DAY=2015001        rem ***开始日期***
set /a DEADLINE=2016001  rem ***结束日期***
:start
if %DAY% leq %DEADLINE% (goto ORDER) else exit
:ORDER
rem **将当天的图幅数据文件名放在一个TXT文件中**
dir *%DAY%.*.hdf/a/b/s > MOSAICINPUT.TXT
rem **拼接HDF12种数据中的第一种**
D:/Modis/MRT/bin/mrtmosaic.exe -i MOSAICINPUT.TXT -s "1 0 0 0 0 0 0 0 0 0 0 0" -o MOSAIC_TMP_%DAY%.hdf 
rem **将拼接后数据复制到已建Result 文件夹中并删除当天的数据**
copy MOSAIC_TMP_%DAY%.hdf  Result & del MOSAIC_TMP_%DAY%.hdf
del *%DAY%.*.hdf
set /a DAY= %DAY% + 365   rem ***拼接下一个时相的数据***
goto start

需要修改的地方
a. 第一行MRTDATADIR的路径为下载的MRT文件夹\bin\ModisTool.bat，右键选中编辑，“set MRT_DATA_DIR=…"语句对应路径；
b. 开始和结束时间，MCD12Q1数据的开始和结束时间相差一年，对应的倒数第二行DAY的时间为365；
c. 拼接数据（第10行），对应mrtmosaic.exe的路径，12个0代表12种数据，将0改为1则表示拼接该种数据。本数据提取第一种LC_1（采用IGBP分类方法），若采用第二种LC_2（UMD分类方法），则改为”010000000000“。

拼接完成后生成的resample.log记录了拼接过程，点开可查看出错原因。
可以利用MRT或者ENVI、ArcMap进行投影转换。

三、ArcMap 面积制表

步骤

投影后的Modis数据，通过A省shp地图进行掩膜提取，获得土地利用类型图（ArcToolbox→Spatial analysis→提取分析→按掩膜提取）；
利用面积制表，提取出每个县的土地利用类型的面积统计表（Excel）（ArcToolbox→Spatial analysis→区域分析→面积制表）。
第一行shp文件；第二行ID号或其它唯一字段；第三行土地利用数据；第四行分类类型编号Value；…