前言
使用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校正》的相關內容。
示例數據