批量給MODIS拼接和重投影參考模板

1、A. 先進(jìn)行拼接處理1. 先新建一個(gè)DATA文件夾，并將HDF數(shù)據(jù)拷貝到該文件夾下2. 在DATA文件夾下再新建一個(gè)Result文件夾（用于存放拼接后的數(shù)據(jù)）3. 將以下代碼內(nèi)容粘貼到一個(gè)txt文件中，更改txt文件擴(kuò)展名為bat即可，假設(shè)你命名該bat文件名為MODISmosaic.bat；程序?qū)ёx：rem 表示開始的為注釋 ；MOSAICINPUT.TXT這是程序運(yùn)行時(shí)自動(dòng)生成的，不用管；MRTDATADIR為你的MRT安裝文件中data的路徑，“c:/MRT/bin/mrtmosaic.exe”改成你的mrtmosaic.exe的安裝路徑。"set /a DAY= %DAY%

2、+ 16 "則是因?yàn)檩斎霐?shù)據(jù)是16天間隔的，根據(jù)你的數(shù)據(jù)進(jìn)行相應(yīng)修改。rem Set the MRTDATADIR environmental var to the MRT data directory.set MRTDATADIR=C:MRTdataset /a DAY=2011097 rem *batch data start time*set /a DEADLINE=2011273 rem *batch data end time*:startif %DAY% leq %DEADLINE% (goto ORDER) else exit:ORDERrem *save the fi

3、le name into a notepad*dir *%DAY%.*.hdf/a/b/s > MOSAICINPUT.TXTrem *execute mosaic *rem Set the mrtmosaic.exe directory.c:/MRT/bin/mrtmosaic.exe -i MOSAICINPUT.TXT -s "1 0 0 0 0 0 0 0 0 0 0 0" -o MOSAIC_TMP_%DAY%.hdfrem *copy the result to a file and delete the input data*copy MOSAIC_TM

4、P_%DAY%.hdf Result & del MOSAIC_TMP_%DAY%.hdfdel *%DAY%.*.hdfset /a DAY= %DAY% + 16goto start4. 點(diǎn)擊MODISmosiac.bat，即可進(jìn)行拼接操作注意：該操作數(shù)據(jù)及BAT文件須放在一個(gè)文件夾（文件夾起名請(qǐng)用英文，MRT對(duì)中文路徑數(shù)據(jù)會(huì)打不開）下。代碼未考慮中間天數(shù)間斷的情況，比如對(duì)MODIS時(shí)間分辨率為16天的數(shù)據(jù)，在一個(gè)應(yīng)連續(xù)的等差數(shù)列中間有數(shù)據(jù)缺失，存在相差32天的情況，可能會(huì)報(bào)錯(cuò)，請(qǐng)注意。同時(shí)2 / 6，需要備份好你的輸入數(shù)據(jù)，因?yàn)槌绦蜻\(yùn)行中會(huì)自動(dòng)刪除拼接好的數(shù)據(jù)。進(jìn)行到此處則得到了

5、初步的拼接數(shù)據(jù)，下面進(jìn)行第二步B. 再進(jìn)行重投影，改格式，縮小范圍1. 新建記事本，在記事本中輸入下面代碼，改txt后綴為bat，定義modis.bat文件。rem Set the MRTDATADIR environmental var to the MRT data directory.set MRTDATADIR=C:MRTdatafor %i in (*.hdf) do C:MRTbinresample -p my.prm -i %i -o %iout.tif其中，MRTDATADIR為你的MRT安裝文件中data的路徑2. 利用MRT的圖形界面工具定義投影文件選擇一期剛剛拼接好的影像，在UL Corner和LR Corner可以輸入經(jīng)緯度來縮小你需要的區(qū)域，設(shè)置你的保存路徑以及命名（最好以日期命名，到時(shí)候程序會(huì)自動(dòng)按照你的規(guī)則全部按日期命名），選擇tif格式，選擇你要的投影類型，點(diǎn)擊Save Parameter File將投影文件保存。記?。哼@兩個(gè)文件必須和剛剛拼接后的數(shù)據(jù)放在一個(gè)文件夾中3. 雙擊modis.bat執(zhí)行批處理即可！操作如圖：總的來說，就是先進(jìn)行拼接得到一個(gè)hdf的臨時(shí)文件，然后在進(jìn)行重投影，轉(zhuǎn)格式，縮范圍的操作。整個(gè)過程你只需要安裝MRT即可進(jìn)行，不需要再安裝Cygwin模擬lunix環(huán)境