pyModis开发（一）convertmodis_gdal module模块
pyModis开发（一）convertmodis_gdal module模块
pyModis是一个专门针对MODIS数据的开源python开发包。该开发包主要实现MODIS数据的批量下载和投影转换（正弦投影转其他投影等）、拼接、裁切、数据抽取等处理功能。对于MODIS数据的处理，MRT和Gdal是两个常用的工具和开发库。pyModis结合了两者，提供了多种处理方式。pyMODIS开发包库包含：downmodis module 、parsemodis module 、convertmodis module、convertmodis_gdal module 、qualitymodis module、optparse_required module 和
optparse_gui module 七个主要模块。
convertmodis_gdal module模块
利用python gdal库来转换MODISHDF文件。可以创建GeoTiff文件（或其他gdal支持的文件格式）或者多个带的HDF拼接文件。模块包括三个类和三个函数：
Classes:
- file_info
- createMosaicGDAL
- convertModisGDAL
Functions:
- getResampling()
- raster_copy()
- raster_copy_with_nodata()
convertModisGDAL类
convertModisGDAL(hdfname, prefix, subset, res, outformat=“GTiff”,
epsg=None, wkt=None, resampl=‘NEAREST_NEIGHBOR’, vrt=False)
利用GDAL库将MODIS数据从hdf格式转换成GDAL格式，可以实现hdf的投影转换、子数据集（波段数据）抽取、格式转换、重采样等功能。
参数
- hdfname (str) – 输入数据的文件路径列表，list。
- prefix (str) – 输出数据的前缀
- subset (str) – 子数据集或者波段数据，可以是"(0 1 1)"类似的字符串（支持的字符串格式还包括：‘011’,'0,1,1’等等），也可以是0,1,1的列表，0表示不需要处理，1表示处理。如果设为None，则表示对整个数据集里所有子集进行处理。
- res (int) – 输出分辨率，如果设置为0，则采用源数据的分辨率。
- outformat (str) – 输出格式字符串, GDAL支持的格式都可以，默认“GTiff“。
- epsg (int) – 输出文件投影的EPSG编码。
- wkt (str) – 输出文件投影的WKT字符串。
- resampl (str) – 重采样类型，默认为’NEAREST_NEIGHBOR’,支持的重采样类型还包括：‘AVERAGE’, ‘BILINEAR’, ‘CUBIC’, ‘CUBIC_SPLINE’, ‘LANCZOS’,‘MODE’。
- vrt (bool) – 输入文件是否为通过createMosaicGDAL创建的GDAL VRT 文件，True为使用VRT文件，False为不使用VRT。注意，使用vrt时，只处理第一个子数据集（波段数据）。当需要处理VRT中的多个波段或子数据集时，该操作无效。
主要函数
run(quiet=False)
处理主函数，对已选图层（波段）子集进行重投影变换。如果vrt为True，程序会直接调用run_vrt_separated()，且只会处理vrt的第一个子数据集（波段数据）。
run_vrt_separated()
针对VRT文件的处理函数，对createMosaicGDAL创建的VRT文件进行重投影变换，对应的createMosaicGDAL类的函数write_vrt的参数separated设置为True。
用法示例
对MODIS数据进行多幅拼接
from pymodis import convertmodis_gdal
mosaic_g = createMosaicGDAL(self.__hdfs, self.__subset)
if mosaic_g is None:
return False
else:
mosaic_name_list = []
out_name_list = []
for k in list(mosaic_g.file_infos.keys()):
mosaic_names = "{pref}_{band}.vrt".format(pref=self.__mosaic_name, band=k)
out_names = "{pref}_{band}".format(pref=self.__outname, band=k)
mosaic_name_list.append(mosaic_names)
out_name_list.append(out_names)
if mosaic_g.write_vrt(self.__mosaic_name):
self.__mosaic_files = mosaic_name_list
self.__out_names = out_name_list
return True
else:
return False
createMosaicGDAL类
createMosaicGDAL(hdfnames, subset, outformat=‘HDF4Image’)
利用GDAL拼接MODIS数据从hdf格式到GDAL支持的文件格式。
参数：
- hdfnames (list) – 待拼接的文件路径列表
- subset (str) – 待拼接的子集或波段数据，可以是"(0 1 1)"类似的字符串（支持的字符串格式还包括：‘011’,'0,1,1’等等），也可以是0,1,1的列表，0表示不需要处理，1表示处理。如果设为None，则表示对整个数据集里所有子集进行处理。
- outformat (str) – 输出格式字符串,支持的格式包括：HDF4Image、GTiff、HFA，不支持VRT文件输出。
主要函数
run(output, quiet=False)
处理主函数，创建拼接处理过程。
参数
output (st) - 输出文件的名称。
write_mosaic_xml(prefix)
输出MODIS拼接的XML元数据文件。
参数
prefix(str) - 包含元数据的XML文件的名称前缀
write_vrt(output, separate=True, quiet=False)
输出VRT文件
参数
- output(str) - 输出文件的前缀
- separate(bool) - 是否为各个波段各自输出一个VRT文件，True表示为每个波段输出一个VRT文件，False表示输出到一个文件。
file_info类
GDAL文件信息类,该类拷贝自gdal_merge.py
参数
filename(str) - 待读取的文件名称列表
返回值
1表示成功，0表示文件无法打开
主要函数
copy_into(t_fh, s_band=1, t_band=1, nodata_arg=None)
将此文件映像复制到目标文件中。
这个方法将计算file_info对象文件和目标gdal的重叠区域。对象，并复制用于公共窗口区域的图像数据。假设这些文件在一个相同的投影中。不作检查或纠正。但是，如果目标文件的分辨率不同，或者图像像素类型不同，则会执行适当的重采样和转换(使用常规的GDAL提升/降级规则)。
参数
- t_fh – gdal.Dataset object for the file into which some or all of this file may be copied.
- s_band –
- t_band –
- nodata_arg –
返回值
1代表成功（或者不需要拷贝），0表示失败。
init_from_name(filename)
初始化
getResampling(res)函数
getResampling(res)
Return the GDAL resampling method
Parameters: res (str) – the string of resampling method
raster_copy函数
raster_copy(s_fh, s_xoff, s_yoff, s_xsize, s_ysize, s_band_n, t_fh, t_xoff, t_yoff, t_xsize, t_ysize, t_band_n, nodata=None)[source]¶
Copy a band of raster into the output file.
Function copied from gdal_merge.py
raster_copy_with_nodata函数
raster_copy_with_nodata(s_fh, s_xoff, s_yoff, s_xsize, s_ysize, s_band_n, t_fh, t_xoff, t_yoff, t_xsize, t_ysize, t_band_n, nodata)[source]
Copy a band of raster into the output file with nodata values.
Function copied from gdal_merge.py