Envi 调用 MODIS Reprojection Tool 裁切
(MRT )对MODIS 产品进行批处理拼接、 重投影、
1熟悉MRT
MODIS产品的类型不同,一景 HDF格式的影像包含的波段也各不相同。 需要选择处理波段,0表示不作处理,1表示处理,首先要确定影像的波段数。
MRT处理时
1、拼接 Mrtmosaic.exe
程序用来拼接影像。调用方式为:
mrtmosaic -i \\命令行帮助如下
Usage: mrtmosaic -i in put_file names_file -t -h -o output_file name
-s spectral_subset \
需要指定三个参数:
-I指定输入文件,可以采取两种方式
1、1.hdf 2.hdf 3.hdf 2、把影像的完整路径保存到 txt文件中,作为in put参数。如-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_zo ne -x pixel_size
-g file name for the log file
可以只指定1个参数:
用-p读入prm 参数文件,进行处理。例如 INPUT_FILENAME = g:\\tmp_%1.hdf SPECTRAL_SUBSET = ( 1 1 1 1) resample -p \
,prm文件如下: 输入文件 前面 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