使用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校正》的相關內容。

示例數據

校正前的圓盤數據
校正後的圓盤數據
校正後放大數據

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章