天天看點

預處理——HJ衛星資料

(遙感資料打開後,為了呈現更好的顯示結果,都進行2%的拉伸,下面不再贅述)

(一)圖像裁剪

1、打開2016年的HJ1A-CCD2衛星資料→File→Save as→Spatial Subset→選擇要裁剪的東平湖區域→OK

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

2、同樣的方法對Landsat8資料進行相同區域的裁剪用于幾何配準的baseimage,注意在選擇裁剪區域時點選下方的Subset By File檔案選擇上一步驟中裁剪好的HJ遙感影像資料。

預處理——HJ衛星資料
預處理——HJ衛星資料

(二)幾何配準

使用portal工具檢視HJ資料和Landsat8資料,發現湖的邊界位置錯開不重合,需要進行幾何配準。

預處理——HJ衛星資料

1、找到/Geometric Correction/Registration/Image Registration Workflow工具→base image選擇裁剪好的landsat8資料、warp image選擇裁剪好的HJ資料→next→切換到Advanced面闆,把Matching Band in Warp Image改為Band4。

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

2、針對類似地形複雜的區域,手動選幾個種子點能大大提高自動找點的準确度。

Seed tie points→start editing→在影像上均勻的選點,每選完一個點右擊accept as individual,顯示視窗會自動顯示另一幅影像資料,然後選擇對應的相同地物後右擊accept as individual,進行下一個點的配準。

預處理——HJ衛星資料
預處理——HJ衛星資料

2、選完20個種子點後點選stop editing→next

預處理——HJ衛星資料
預處理——HJ衛星資料

3、系統自動完成選點30個→show table檢視誤差大小,删掉誤差較大的點→最後保留了28個點誤差大小為0.846844

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

4、warping→warping method:folynomial、resampling:cubic convolution、output pixel size from:warp image→next→選擇儲存位置

預處理——HJ衛星資料
預處理——HJ衛星資料

下圖為集合配準後的結果,檢視兩幅資料,可以看到湖的邊界很好的重合。

預處理——HJ衛星資料

(三)大氣校正

1、建波譜響應函數:display→new plot window→import→ASCII…→選擇提前下載下傳好的HJ1ACCD2資料→data units:nanometers→OK→4個波段的波普響應能力圖像就繪制出來了

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

2、點選右側的show→對四個波段進行進行重命名,依次為b1、b2、b3、b4→export→輸出儲存成波普庫檔案spectral library

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

3、convert interleave→選擇幾何校正後的HJ資料→OK→BIL、Yes→OK→是

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

4、flaash atmospheric correction→input radiance image→選擇上個步驟中改為BIL的HJ資料→use single…、Single scale factor改為10→OK

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

5、設定輸入輸出檔案→輸入傳感器高度:650KM、區域的平均海拔0.05KM、分辨率30m→輸入資料擷取時間,在XML檔案中檢視→輸入大氣模型參數,查詢表格可知此資料應選MLS→氣溶膠模型:rural 不需要反演選擇none(詳見下圖闆面設定的内容)

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

6、multispectral setting→GUI→filter function file→open→spectral library→選擇之前儲存的.sli檔案→OK→OK

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

7、advanced setting→tile size改為100→OK→apply

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

可檢視大氣校正後的地物波譜曲線

預處理——HJ衛星資料

(四)裁剪湖泊區域

1、spectral profile→/Feature Extraction/Rule Based Feature Extraction Workflow→輸入進行過大氣校正後的資料→custom bands:normalized difference選擇3、4波段→next→輸入東平湖的矢量邊界.shp

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

2、subset data from rois→選擇要進行裁剪的資料→選擇矢量資料、Yes、選擇存儲位置→OK

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料

3、将黑色的背景去掉

/Raster Management/Edit ENVI Header→選擇要進行處理的資料→OK→edit attributes→data ignore value→data ignore value輸入0→OK→OK

預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料
預處理——HJ衛星資料