批量給MODIS拼接和重投影

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

2、是因?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 file name into

3、 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_TMP_%DAY%.hdf Result & del MOSA

4、IC_TMP_%DAY%.hdfdel *%DAY%.*.hdfset /a DAY= %DAY% + 16goto start4. 點(diǎn)擊MODISmosiac.bat，即可進(jìn)行拼接操作注意：該操作數(shù)據(jù)及BAT文件須放在一個文件夾（文件夾起名請用英文，MRT對中文路徑數(shù)據(jù)會打不開）下。代碼未考慮中間天數(shù)間斷的情況，比如對MODIS時間分辨率為16天的數(shù)據(jù)，在一個應(yīng)連續(xù)的等差數(shù)列中間有數(shù)據(jù)缺失，存在相差32天的情況，可能會報(bào)錯，請注意。同時，需要備份好你的輸入數(shù)據(jù)，因?yàn)槌绦蜻\(yùn)行中會自動刪除拼接好的數(shù)據(jù)。進(jìn)行到此處則得到了初步的拼接數(shù)據(jù)，下面進(jìn)行第二步B. 再進(jìn)行重投影，改格式，縮小范圍1. 新建

5、記事本，在記事本中輸入下面代碼，改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è)置你的保存路徑以及命名（最好以日期命名，到時候程序會自動按照你的規(guī)則全部按日期命名），選擇tif格式，選擇你要的投影類型，點(diǎn)擊Save Parameter File將投影文件保存。記住：這兩個文件必須和剛剛拼接后的數(shù)據(jù)放在一個文件夾中3. 雙擊modis.bat執(zhí)行批處理即可！操作如圖：總的來說，就是先進(jìn)行拼接得到一個hdf的臨時文件，然后在進(jìn)行重投影，轉(zhuǎn)格式，縮范圍的操作。整個過程你只需要安裝MRT即可進(jìn)行，不需要再安裝Cygwin模擬lunix環(huán)境，安