推荐Envi调用MODISReprojectionToolMRT对MODIS产品进行批处理拼接

1、(MRT )对 MODISEnvi 调用 MODIS Reprojection Tool产品进行批处理拼接、重投影、裁切1熟悉MRTMODIS产品的类型不同，一景HDF格式的影像包含的波段也各不相同。 MRT 处理时需要选择处理波段，0表示不作处理，1表示处理，首先要确定影像的波 段数。1、拼接Mrtmosaic.exe 程序用来拼接影像。调用方式为:mrtmosaic -i "g:n%1.txt" -s " 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 1 1 1 0 0 0 0 " -o

2、"g:temp.hdf"命令行帮助如下Usage: mrtmosaic -i in put_file names_file -t -h -o output_file name-s spectral_subset "b1 b2 . bN"-g file name for the log file需要指定三个参数：-I 指定输入文件，可以采取两种方式1、1.hdf 2.hdf 3.hdf2、把影像的完整路径保存到 txt 文件中，作为 input 参数。如 -IMOD092008001.txt 。-s 指定需要处理的波段，同样可以采用两种方式1、直接给出，女

3、口 - s “ 0 0 0 0 0 0 0 0 0 0 0 0”，注意英文引号;2、指定一个 txt 路径，让程序读取；-o 指定输出路径，一般直接给出1、直接给出，如 P g:tmp.hdf 。注意直接存为HDF格式，便于后续处理。2、重投影、裁切Resample.exe 用来重投影、裁切是 MRT 程序的核心。调用方式为： 命令行帮助如下：Usage: resample -p parameter_file optionsOptions that override parameter file specifications:-i input_file_name-o output_file_n

4、ame-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 "p1 p2 . p15"-s spectral_subset "b1 b2 . bN"If using the -s switch, the SDSs should be represented as anarray of 0s and 1s. A '1' spec

5、ifies 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 "ULlat ULlong LRlat LRlong"

6、;-or- "ULline ULsample LRline LRsample (0-based)"-or- "ULprojx ULprojy LRprojx LRprojy"NOTE: line/sample must be specified for thehighestresolutio nto be processed in the product.of all SDSs specified-u UTM_zo ne-x pixel_size-g file name for the log file可以只指定1个参数：用-p读入prm参数文件，进行处

7、理。例如 resample -p"g:prrmMOD092008001.prm"， 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

8、 = F:MRT_outmyd%1.tif 输出路径，不同波段自动区分 RESAMPLING_TYPE = NEAREST_NEIGHBOR 最近邻采样方法 OUTPUT_PROJECTION_TYPE = UTM 输出文件投影方式 utm OUTPUT_PROJECTION_PARAMETERS =(投影参数0.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.0 )DATUM = WGS84 大地水准面UTM_ZONE = 50 投影分带带号如果不会设置，可以先在图形界面里设置一次，把设置的结果保存下来。打开刚刚保存的参数文件去掉以#号

9、打头的注释，文件显示如下：INPUT_FILENAME = G:TmpMosaic.hdfSPECTRAL_SUBSET = ( 1 1 1 1 1 1 1 1 1 1 1 1 )SPATIAL_SUBSET_TYPE = INPUT_LAT_LONGSPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 )SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )OUTPUT_FILENAME = G:mod092001001.tifRESAMPLING_TYPE = NEAREST_NEIGHBOROUTPUT_PROJECTION_TYPE

10、 = UTMOUTPUT_PROJECTION_PARAMETERS =(0.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.0 )DATUM = WGS84UTM_ZONE = 50如果研究区的经纬度范围不同，可以修改SPATIAL_SUBSET_UL_CORNER = ( 32.0 107.0 )SPATIAL_SUBSET_LR_CORNER = ( 29.0 116.0 )两行如果需要不同的投影，如经纬度投影，可以修改OUTPUT_PROJECTION_TYPE = UTMOUTPUT_PROJECTION_PARAMETERS

11、=(0.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.00.0 0.0 0.0 )DATUM = WGS84UTM_ZONE = 50如果拼接时只选择了 1个波段，可以修改SPECTRAL_SUBSET = ( 1 )当然这句INPUT_FILENAME = G:TmpMosaic.hdf批处理是必须改掉的了关于像素分辨率，一般留空，也就是默认不改变影像的像素大小2关于DOS批处理裁切MOD产品单次处理Dir g:mod11*A%1*.hdf /s /b >> g:n%1.txtf:mrtmodisbi n'mrtmosaic -i &q

12、uot;g:n%1.txt" -s " 0 0 0 0 0 0 0 0 0 0 0 0 0 00 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 0 " -o "g:temp_%1.hdf"f:mrtmodisbi n'resample -p "g:prm%1.prm"del g:prm%1.prm /q /fdel g:n%1.txt /q /fdel g:temp_%1.hdf /q /f批量处理or /f %i in (e:lis.txt) do 单次处理文件名%i如果双击运行，用%i号

13、表示变量。在命令行中%i表示变量其中e:lis.txt是年份及天数2008001200800920080172008025200803320080412008049G:n%1.txt在运行时为g:n2008001.txt ，其中包括2008年第1天的多景MODIS产品影像名。Dir g:mod11*A%1*.hdf /s /b >> g:n%1.txt用来列出g盘（数据存储盘）以mod11打头的指定日期的的hdf格式文件，即MODIS影像。这里要说明一下，MODIS产品命名的规则MYD11A2.A2002185.h27v05.005.2007221162948.hdf从左至右，前7

14、位表示产品类型MYD11指下午星（MYD）第11号产品（landsurface temperature ） ，A2表示处理级别。A2002185表示过境时间是2002年第185天，h27v05是MODIS产品在全球的轨道行列号，中国一般是h27v05、h27v06、h28v05、h28v06 ,005 表示 HDF5 文件格式，早期 MODIS 采用 HDF4格式存储，所以有时也可见到004的产品。2007等一串数字表示影像处理的年月曰时分秒，hdf表示存储格式为EOS-HDF 了解了 MODIS产品命名规则之后，可以用 DOS dir命令来挑选符合条件的 MODIS产品路径了。>>

15、; g:n%1.txt表示把结果转存为g:n%1.txt%1表示接收的第一个参数，即20080013Envi生成时间列表2008001200800920080172008025这样的列表文件可以用在 Envi中定义一个pro来生成：pro daylist,startdate,e nddate,result_txt_path;example daylist,2008001,2008013,'example.txt'JJa=(startdate mod 1000) /8 *8+(startdate /1000)*1000+1b=(e nddate mod 1000) /8 +1)*8+(e nddate /1000)*1000+1help,a,bstep=8ope nw,l un ,result_txt_path,/get_l unres=a;for i=a,b,step do beg inprin tf,lu n,stri ng(r