使用GDAL對靜止衛星圓盤資料進行校正（以FY2為例子）
阿新 • • 發佈：2018-11-03
使用GDAL對靜止衛星圓盤資料進行校正（以FY2為例子）
前言
使用GDAL對靜止衛星資料的校正的方式與之前寫的校正MODIS資料一樣，只不過對於靜止衛星的圓盤資料而言，經緯度查詢是固定的，不會發生變化。
經緯度查詢表資料
對於FY2氣象衛星而言，經緯度查詢表文件（NOM_ITG_2288_2288(0E0N)_LE.zip）可以從網上下載。資料解壓之後，裡面有三個檔案，分別是：
| 序號 | 名稱 | 說明 |
|---|---|---|
| 1 | NOM_ITG_2288_2288(0E0N)_LE.dat | 經緯度查詢表資料 |
| 2 | NOM_ITG_2288_2288(0E0N)_LE.hdr | 經緯度查詢表資料標頭檔案 |
| 3 | Readme for FY latlon.txt | 資料說明檔案 |
在資料說明檔案中，詳細說明了經緯度查詢表資料的使用方法，同時針對FY2系列衛星的經度進行了說明，由於上述經緯度資料是以中心點為0度經度生成的，所以對於FY2的資料，需要在經度資料上加上衛星所在的經度。
| 序號 | 衛星名稱 | 星下點經度 |
|---|---|---|
| 1 | FY-2F | 112°E |
| 2 | FY-2E | 104.5°E |
| 3 | FY-2D | 86.5°E |
| 4 | FY-2C | 104.5°E |
影像校正
將上面的經緯度查詢表資料按照部落格的內容構造一個Geolocation資訊，然後呼叫gdalwarp就可以處理。
需要注意的是，需要事先對經緯度資料中的經度全部加上衛星的星下點經度，然後再構造GeoLocation資訊。也可以不用這麼處理，直接校正，校正結束之後，修改影象的仿射變換引數，先通過GDALDataset::GetGeotransform函式獲取仿射變換引數，然後給第一個值加上衛星的星下點經度之後，再使用GDALDataset::SetGeotransform函式更新到校正結果資料中即可。
使用GDAL工具校正
具體的處理方式參考《使用GDAL工具對FY3系列衛星資料進行校正》以及《使用GDAL對HDF資料進行校正》兩篇部落格。
自己編寫程式碼校正
具體處理步驟參考《使用GDAL對HDF資料進行geoloc校正》的相關內容。
示例資料
