modis數(shù)據(jù)預(yù)處理

《modis數(shù)據(jù)預(yù)處理》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。

1、MODIS數(shù)據(jù)預(yù)處理1.波段設(shè)置Modis影像有三種打開方式,一般我們用打開外部文件的方式打開科學(xué)數(shù)據(jù)集,因為需要數(shù)據(jù)集中的一些輔助信息(主要是太陽幾何,衛(wèi)星幾何).但是這樣打開之后顯示的波段從1開始的,而數(shù)據(jù)集中對應(yīng)的modis通道并不是這個順序.通過菜單欄中的basic_tools->preprocessing->data_specificutilities->viewHDFdatasetattributes可以打開數(shù)據(jù)集里每個要素的屬性表,在里面選中需要的HDF文件中的數(shù)據(jù)集，就會打開其屬性表，波段內(nèi)容如下Attri

2、bute3-5:"band_names""8,9,10,11,12,13lo,13hi,14lo,14hi,15,16,17,18,19,26"對應(yīng)打開的HDF文件里1KM輻亮度文件的波段數(shù),一共16個波段.其中13/14波段比較特殊,都有hi和lo兩組數(shù)據(jù),它們是傳感器高敏感度和低敏感度兩種狀態(tài)下獲取到的DN值,分別對應(yīng)于較暗地物和較亮地物,使用哪個文件根據(jù)需要而定.但是在太湖湖區(qū),13/14波段大部分區(qū)域效果都不太好.值會很大,出現(xiàn)溢出.可能是由于太湖的高渾濁度.1.幾何校正幾何校正有三種方法:1)用envi自帶模塊進

3、行幾何校正,通過菜單欄中的Map->GeoreferencesMODIS選中envi中已經(jīng)打開的需要校正的數(shù)據(jù)集，輸入研究區(qū)的地理位置，如下圖左，投影用UTM，基準(zhǔn)面用WGS-84，區(qū)域根據(jù)經(jīng)緯度確定。輸入完成，envi會自動校正，并執(zhí)行去蝴蝶結(jié)效應(yīng)算法，有點是能對我們需要的那些波段進行校正。缺點也很明顯。如下圖右，校正結(jié)束的圖像會失去原始圖像四個角的信息，這樣就無法和GLT校正的圖像很好的匹配起來，不利于一些后續(xù)的處理。2)用GLT，即是查找表法對圖像進行幾何校正Map->GeoreferencefrominputGeo

4、metry->buidGLT用來建立查找表。在彈出的對話框中選擇查找表的XY信息，其中X對應(yīng)圖像經(jīng)度信息，Y對應(yīng)緯度信息。然后只需要規(guī)定投影、基準(zhǔn)面和區(qū)位信息，就可以生成一個查找表文件。這個查找表文件的實質(zhì)也是兩幅圖像，分別在每個像元上保存著經(jīng)緯度值，并且像元位置是拉伸到我們規(guī)定的輸出投影上面去了，而且是逐像元的拉伸。那么剩下的矯正工作就只是把想要矯正的信息和查找表一一匹配起來，因此速度也很快。Map->GeoreferencefrominputGeometry->GeoreferencefromGLT就是上面所說的，把查

5、找表運用到需要矯正的每一個像元上。這里只需要選擇想要矯正的影響就好了。有一點需要注意：這里查找表的分辨率比輻亮度信息要低，因此如果直接用這個查找表對原始的輻亮度圖像進行矯正，只能得到原始圖像一小部分的信息。因為矯正時是按照像元行列號一一對應(yīng)查詢的。解決方案有兩個：可以把原始圖像分辨率降低，重采樣成和查找表一樣大小的圖像，再用查找表進行幾何校正，這樣運行速度快但是損失原始圖像信息；也可以將查找表重采樣成和原始的輻亮度圖像一樣大小的圖像，再用這個新的查找表對圖像進行矯正。這樣就不會損失原始圖像信息，但是計算速度會大大降低。視應(yīng)

6、用選擇。1)IDL批處理;forward_functionenvi_proj_createPROModis_gef_batchenvi,/restore_base_save_files;恢復(fù)ENVIsav文件envi_batch_init,log_file=’batch.txt’;開始批處理模式inpath=DIALOG_PICKFILE(/DIRECTORY,$TITLE="selectMODISfilespath")CD,inpathfilename=FILE_SEARCH('*.HDF');print,resultn

7、=N_ELEMENTS(filename)outpath=DIALOG_PICKFILE(/DIRECTORY,$TITLE="selectMODISoutpath")FORi=0,n-1DOBEGINin_name=inpath+filename[i]out_name='ReGeo'+filename[i];設(shè)置校正方法;0=RadianceEmissivity,1=ReflectanceEmissivitycalib_method=1;設(shè)置輸出方法;0=Standard,1=Projected,2=Standard

8、andProjectedout_method=1;設(shè)置輸出投影output_projection=envi_proj_create(/geographic);在輸出時設(shè)置去除蝴蝶效應(yīng)convert_modis_data,in_file=in_name,out_path=outpath,$out_root