Envi调用MODIS Reprojection Tool(MRT)对MODIS产品进行批处理拼

Envi调用MODIS Reprojection Tool（MRT）对MODIS产品进行批处理拼接、重投影、裁切

1 熟悉MRT

MODIS产品的类型不同，一景HDF格式的影像包含的波段也各不相同。MRT处理时需要选择处理波段，0表示不作处理，1表示处理，首先要确定影像的波段数。 1、拼接

Mrtmosaic.exe程序用来拼接影像。调用方式为：

mrtmosaic -i \0 0 0 0 1 1 1 0 0 0 0 \命令行帮助如下

Usage: mrtmosaic -i input_filenames_file -t -h -o output_filename -s spectral_subset \ -g filename for the log file

需要指定三个参数：

-I 指定输入文件，可以采取两种方式 1、1.hdf 2.hdf 3.hdf

2、把影像的完整路径保存到txt文件中，作为input参数。如-I MOD092008001.txt。

-s 指定需要处理的波段，同样可以采用两种方式

1、直接给出，如 –s “0 0 0 0 0 0 0 0 0 0 0 0”，注意英文引号； 2、指定一个txt路径，让程序读取；

-o 指定输出路径，一般直接给出

1、直接给出，如 –o g:\\tmp.hdf。注意直接存为HDF格式，便于后续处理。

2、重投影、裁切

Resample.exe用来重投影、裁切是MRT程序的核心。调用方式为： 命令行帮助如下：

Usage: resample -p parameter_file [options]

Options that override parameter file specifications: -i input_file_name -o output_file_name

-r resampling_type [NN BI CC NONE]

-t projection_type [AEA ER GEO HAM IGH ISIN LA LCC MERCAT MOL PS SIN TM UTM]

-j projection_parameter_list \ -s spectral_subset \

If using the -s switch, the SDSs should be represented as an array of 0s and 1s. A '1' specifies to process that SDS;

'0' specifies to skip that SDS. Unspecified SDSs will not be processed. If the -s switch is not specified, then all SDSs will be processed. -a spatial_subset_type [INPUT_LAT_LONG INPUT_LINE_SAMPLE OUTPUT_PROJ_COORDS]

-l spatial_subset \

-or- \ -or- \

NOTE: line/sample must be specified for the highest resolution of all SDSs specified to be processed in the product.

-u UTM_zone -x pixel_size

-g filename for the log file

可以只指定1个参数：

用 –p 读入prm参数文件，进行处理。例如resample -p \，prm文件如下： INPUT_FILENAME = g:\\tmp_%1.hdf输入文件

SPECTRAL_SUBSET = ( 1 1 1 1) 前面mrtmosaic拼接结果有4个波段

SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG 经纬度裁切方式

SPATIAL_SUBSET_UL_CORNER = ( 33.0 108.0 ) 左上纬经度 SPATIAL_SUBSET_LR_CORNER = ( 28.0 117.0 ) 右下纬经度

OUTPUT_FILENAME = F:\\MRT_out\\myd%1.tif 输出路径，不同波段自动区分

RESAMPLING_TYPE = NEAREST_NEIGHBOR 最近邻采样方法

OUTPUT_PROJECTION_TYPE = UTM 输出文件投影方式utm

OUTPUT_PROJECTION_PARAMETERS = (投影参数 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0

0.0 0.0 0.0 0.0 0.0 0.0 )

DATUM = WGS84 大地水准面

UTM_ZONE = 50投影分带带号

如果不会设置，可以先在图形界面里设置一次，把设置的结果保存下来。 打开刚刚保存的参数文件

去掉以#号打头的注释，文件显示如下： INPUT_FILENAME = G:\\TmpMosaic.hdf

SPECTRAL_SUBSET = ( 1 1 1 1 1 1 1 1 1 1 1 1 )

SPATIAL_SUBSET_TYPE = INPUT_LAT_LONG

SPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 ) SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )

OUTPUT_FILENAME = G:\\mod092001001.tif

RESAMPLING_TYPE = NEAREST_NEIGHBOR

OUTPUT_PROJECTION_TYPE = UTM

OUTPUT_PROJECTION_PARAMETERS = (

0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 )

DATUM = WGS84

UTM_ZONE = 50

如果研究区的经纬度范围不同，可以修改

SPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 ) SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )两行 如果需要不同的投影，如经纬度投影，可以修改 OUTPUT_PROJECTION_TYPE = UTM

OUTPUT_PROJECTION_PARAMETERS = ( 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 )

DATUM = WGS84

UTM_ZONE = 50

如果拼接时只选择了1个波段，可以修改

百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读，免费范文网，提供经典小说综合文库Envi调用MODIS Reprojection Tool(MRT)对MODIS产品进行批处理拼在线全文阅读。