[Python]Modis数据处理库PyModis ⽂章⽬录
前⾔
介绍
【NASA官⽹】
1.
2.
【PyModis简介】基于Python的开源Modis数据处理库
1. 下载功能：根据⽤户提供时间批量下载
2. 读取：提取数据元数据（包含数据质量信息）
3. 处理：镶嵌、投影、转格式
【官⽅介绍】
python处理xml文件
1. 下载⼤量modis hdf/xml ⽂件以及解析xml⽂档以获取有关HDF⽂件的信息
2. 使⽤MRT或GDAL将modis⽂件转为geotiff格式
3. 使⽤MRT或GDAL将modis⽡⽚数据拼接、重投影
4. ⽤于下载⼤量modis hdf/xml ⽂件以及解析xml⽂档以获取有关HDF⽂件的信息
5. 创建融合后影像的xml格式元数据
6. 从modis字节编码质量评估图层提取具体信息
使⽤⽅法
【两种使⽤⽅法】
1. 脚本⼯具（.py）：命令⾏调⽤
2. Pymodis库：Python代码调⽤
安装
【PyModis安装】
1. PyModis依赖于GDAL库：
2. 安装PyModis：pip install pymodis
PyModis库使⽤
【预处理】
1. createMosaicGDAL实现拼接
2. convertModisGDAL将hdf波段导出为TIFF，并可对TIFF重采样、投影变化等
downmodis
参数说明
【官⽅说明】
url        下载地址 [默认=v]
password 密码
user  ⽤户名
tiles  要下载产品的⾏列号
[格式] ⽤逗号隔开或数组形式
[default=none] 下载所有
[⽰例] "h18v03,h18v04"
[解释] 下载 18⾏3列和 18⾏4列的数据⽡⽚
path  url与产品⽬录之间的路径
product    产品名称
[默认=MOD11A1.005]
delta      从下载⽇期开始，往前下载多少天
[默认=10]
today    开始时间
[格式]YYYY-MM-DD
[默认=none] 默认为今天
enddate    结束⽇期（应该⽐today早，因为它是从today往前下载的）
[说明]如果delta与enddate都有设置，以enddate为主
jpg  是否下载jpg数据
convertmodis_gdal
参数说明
【hdfname】输⼊的hdf⽂件路径
【prefix】输出⽂件的名字前缀，程序会在后⾯⾃动加上⼀个波段名字
prefix = r"D:\mycode\GISandPython\2PyModis\tmp\out"
# 输出的⽂件名字为：
# D:\mycode\GISandPython\2PyModis\tmp\out_Image_Optical_Depth_Land_And_Ocean.tif # 其中，Image_Optical_Depth_Land_And_Ocean为第11波段的名字，为程序⾃动添加的【subset】处理的波段
# 【⽅法1】下标对应为1，则为要处理的波段
subset =[0for x in range(0,219)]#⽣成219个0的数组
subset[11]=1#处理第11个波段
subset[14]=1#处理第14个波段
【res】分辨率，单位即为坐标系的单位
【outformat】输出的格式（GDAL库中的格式）
1. GTiff：TIFF
【epsg】坐标系的EPSG编码，可
【wkt】坐标系的第⼆种设置⽅式：wkt格式的字符串（prj⽂件）
wkt=r"D:\WGS84.prj"# 投影⽂件
# 你可以在ArcGIS⾥创建⼀个对应坐标系的shp⽂件，⽣成的.prj即可⽤
【resampl】重采样⽅法
1. NEAREST_NEIGHBOR
【vrt】是否⽣成影像的头⽂件（默认为False：不⽣成）
⽰例
import pymodis
# 处理的波段
subset =[0for x in range(0,219)]
subset[11]=1#只处理第11波段
# 转换成TIFF格式
convertmodis = vertModisGDAL(
hdfname = r"D:\mycode\GISandPython\2PyModis\2001\MOD04_L2.A2001001.0110.061.2017220225602.hdf",    prefix = r"D:\mycode\GISandPython\2PyModis\tmp\4326",
subset= subset,
res=1,#设置的为WGS84坐标系，那么单位为°，则此为设置为1°
outformat='GTiff',
epsg=4326,#WGS84
resampl='NEAREST_NEIGHBOR'
)
convertmodis.run()
【程序运⾏结果】
【Projection情况】
【分辨率】1°
【坐标系】WGS84
相关链接
1.

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。