1、Envi調(diào)用MODISReprojectionTool(MRT)對(duì)MODIS產(chǎn)品進(jìn)行批處理拼接、重投影、裁切1熟悉MODIS產(chǎn)品的類型不同，一景HDF格式的影像包含的波段也各不相同。MRT處理時(shí)需要選擇處理波段，0表示不作處理，1表示處理，首先要確定影像的波段數(shù)。1、拼接Mrtmosaic.exe程序用來拼接影像。調(diào)用方式為mrtmosaic-ig:n%1.txt-s0000000000000000000000000001110000-og:temp.hdf命令行幫助如下Usage:mrtmosaic-iinput_filenames_file-t-h-ooutput_filename-ssp

2、ectral_subsetb1b2.bN-gfilenameforthelogfile需要指定三個(gè)參數(shù)：-I指定輸入文件，可以采取兩種方式1、1.hdf2.hdf3.hdf2、把影像的完整路徑保存到txt文件中，作為input參數(shù)。如-IMOD092008001.txt。-s指定需要處理的波段，同樣可以采用兩種方式1、直接給出，如S000000000000注意英文引號(hào);2、指定一個(gè)txt路徑，讓程序讀??；-o指定輸出路徑，一般直接給出1、直接給出，如og:tmp.hdf。注意直接存為HDF格式，便于后續(xù)處理。2、重投影、裁切Resample.exe用來重投影、裁切是MRT程序的核心。調(diào)用方式為

3、:命令行幫助如下：Usage:resample-pparameter_fileoptionsOptionsthatoverrideparameterfilespecifications:-iinput_file_name-ooutput_file_name-rresampling_typeNNBICCNONE-tprojection_typeAEAERGEOHAMIGHISINLALCCMERCATMOLPSSINTMUTM-jprojection_parameter_listp1p2.p15-sspectral_subsetb1b2.bNIfusingthe-sswitch,theSDSss

4、houldberepresentedasanarrayof0sand1s.A1specifiestoprocessthatSDS;0specifiestoskipthatSDS.UnspecifiedSDSswillnotbeprocessed.Ifthe-sswitchisnotspecified,thenallSDSswillbeprocessed.-aspatial_subset_typeINPUT_LAT_LONGINPUT_LINE_SAMPLEOUTPUT_PROJ_COORDS-lspatial_subsetULlatULlongLRlatLRlong-or-ULlineULsa

5、mpleLRlineLRsample(0-based)-or-ULprojxULprojyLRprojxLRprojyNOTE:line/samplemustbespecifiedforthehighestresolutionofallSDSsspecifiedtobeprocessedintheproduct.-uUTM_zone-xpixel_size-gfilenameforthelogfile可以只指定1個(gè)參數(shù)：用p讀入prm參數(shù)文件，進(jìn)行處理。例如resample-pg:prrmMOD092008001.prm，prm文件如下：INPUT_FILENAME=g:tmp_%1.hdf輸

6、入文件SPECTRAL_SUBSET=(1111)前面mrtmosaic拼接結(jié)果有4個(gè)波段SPATIAL_SUBSET_TYPE=INPUT_LAT_LONG經(jīng)緯度裁切方式SPATIAL_SUBSET_UL_CORNER=(33.0108.0)左上緯經(jīng)度SPATIAL_SUBSET_LR_CORNER=(28.0117.0)右下緯經(jīng)度OUTPUT_FILENAME=F:MRT_outmyd%1.tif輸出路徑，不同波段自動(dòng)區(qū)分RESAMPLING_TYPE=NEAREST_NEIGHBOR最近鄰采樣方法OUTPUT_PROJECTION_TYPE=UTM輸出文件投影方式utmOUTPUT_PR

7、OJECTION_PARAMETERS=(投影參數(shù)0.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0)DATUM=WGS84大地水準(zhǔn)面UTM_ZONE=50投影分帶帶號(hào)如果不會(huì)設(shè)置，可以先在圖形界面里設(shè)置一次，把設(shè)置的結(jié)果保存下來打開剛剛保存的參數(shù)文件去掉以#號(hào)打頭的注釋，文件顯示如下：INPUT_FILENAME=G:TmpMosaic.hdfSPECTRAL_SUBSET=(111111111111)SPATIAL_SUBSET_TYPE=INPUT_LAT_LONG(32.0107.0)(29.0116.0)SPATIAL_SUBSET_UL_CO

8、RNER=SPATIAL_SUBSET_LR_CORNER=OUTPUT_FILENAME=G:mod092001001.tifRESAMPLING_TYPE=NEAREST_NEIGHBOROUTPUT_PROJECTION_TYPE=UTMOUTPUT_PROJECTION_PARAMETERS=(0.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0)DATUM=WGS84UTM_ZONE=50如果研究區(qū)的經(jīng)緯度范圍不同，可以修改(32.0107.0)SPATIAL_SUBSET_UL_CORNER=SPATIAL_SUBSET_LR_CORNER=(

9、29.0116.0)兩行如果需要不同的投影，如經(jīng)緯度投影，可以修改OUTPUT_PROJECTION_TYPE=UTMOUTPUT_PROJECTION_PARAMETERS=(0.00.00.00.00.00.00.00.00.00.00.00.00.00.00.0)DATUM=WGS84UTM_ZONE=50如果拼接時(shí)只選擇了1個(gè)波段，可以修改SPECTRAL_SUBSET=(1)當(dāng)然這句INPUT_FILENAME=G:TmpMosaic.hdf批處理是必須改掉的了關(guān)于像素分辨率，一般留空，也就是默認(rèn)不改變影像的像素大小。2關(guān)于DOS批處理裁切MOD產(chǎn)品單次處理Dirg:mod11*A%

10、1*.hdf/s/bg:n%1.txtf:mrtmodisbinmrtmosaic-ig:n%1.txt-s0000000000000000000000000001110000-og:temp_%1.hdff:mrtmodisbinresample-pg:prm%1.prmdelg:prm%1.prm/q/fdelg:n%1.txt/q/fdelg:temp_%1.hdf/q/f批量處理for/f%iin(e:lis.txt)do單次處理文件名%i如果雙擊運(yùn)行，用%i號(hào)表示變量。在命令行中i表示變量。其中e:lis.txt是年份及天數(shù)G:n%1.txt在運(yùn)行時(shí)為g:n2008001.txt,其

11、中包括2008年第1天的多景MODIS產(chǎn)品影像名。Dirg:mod11*A%1*.hdf/s/bg:n%1.txt用來列出g盤(數(shù)據(jù)存儲(chǔ)盤)以mod11打頭的指定日期的的hdf格式文件，即MODIS影像。這里要說明一下，MODIS產(chǎn)品命名的規(guī)則MYD11A2.A2002185.h27v05.005.2007221162948.hdf從左至右，前7位表示產(chǎn)品類型MYD11指下午星(MYD)第11號(hào)產(chǎn)品(landsurfacetemperature)，A2表示處理級(jí)別。A2002185表示過境時(shí)間是2002年第185天，h27v05是MODIS產(chǎn)品在全球的軌道行列號(hào)，中國一般是h27v05、h27v06、h28v05、h28v06,005表示HDF5文件格式，早期MODIS采用HDF4格式存儲(chǔ)，所以有時(shí)也可見到004的產(chǎn)品。2007等一串?dāng)?shù)字表示影像處理的年月日時(shí)分秒，hdf表示存儲(chǔ)格式為EOS-HDF。了解了MODIS產(chǎn)品命名規(guī)則之后，可以用DOSdir命令來挑選符合條