基於GLT方法的風雲三號C（FY-3C）衛星影像幾何校正

軟體開發 · 發表 2018-10-24 21:43:23

摘要：第一步、資料開啟風雲三號C 資料以 HDF5 格式儲存，包括一個數據檔案和一個幾何定位檔案。（1） File > Open As >Scientific Formats >HDF5 ，選擇風雲 3 的資料檔案： FY3...

第一步、資料開啟

風雲三號C 資料以 HDF5 格式儲存，包括一個數據檔案和一個幾何定位檔案。

（1） File > Open As >Scientific Formats >HDF5 ，選擇風雲 3 的資料檔案： FY3C_VIRRX_GBAL_L1_20170101_0235_1000M_MS.HDF 檔案，進入 Select HDF5 Datasets 面板。

圖1 HDF5 Datasets 面板

（2）檔案中包含很多資訊，選擇影象資料 EV_RefSB ，點選中間的載入箭頭，將所選資料載入到右邊的列表中，單擊右下角 Open Rasters ，開啟影象。

圖2 FY-3 EV_RefSB 資料（假彩色顯示）

（3）同樣的方法，點選 File > Open As >Scientific Formats >HDF5 ，選擇 FY-3 的幾何定位檔案： FY3C_VIRRX_GBAL_L1_20170101_0235_GEOXX_MS.HDF 檔案，進入 Select HDF5 Datasets 面板相同方法，同時將定位檔案（ Latitude 和 Longitude ）載入並在 ENVI 中開啟。

圖3 通過 HDF Datasets 載入 FY-3 經緯度資料

圖4 FY-3 經緯度資料

（4）部分 FY-3 幾何定位資料，即經緯度資料存在無效值，如無效值為 -999.900024 表示。對於存在無效值的經緯度資料，在構建 GLT 時會出現異常，因此需要對經緯度資料的無效值進行修復。

注：對於無無效值的經緯度資料，跳過此步驟。

圖5 FY-3 經緯度資料無效值檢視

無效值修復方法為：

1) 在 ToolBox 中，選擇 /Terrain/Replace Bad Values

2) 在彈出的對話方塊中選擇 FY3C_VIRRX_GBAL_L1_20170101_0235_GEOXX_MS.HDF 檔案，點選 OK 。

3) 點選 OK 後彈出引數面板， Max 填入 -200 ，設定輸出檔名稱，其它引數保持預設，點選 OK 執行。

注：min不填自動讀取影象最小值，Max大於經緯度最小-180即可

圖6 FY-3 經緯度資料無效值修復

第二步、幾何定位

（1）在 Toolbox 中，選擇 /Geometric Correction/Build GLT 。

（2）在彈出的對話方塊中， Input X Geometry Band 選擇修復後的經度 longitude 波段， Input Y Geometry 選擇修復後的緯度 latitude 波段。在接下來彈出的對話方塊中填寫輸出 GLT 檔案的投影資訊，這裡選擇 Geographic Lat/Lon ， WGS-84 ，如下圖所示。