PIE SDK大氣校訂

1. 算法功能簡介

    大氣校訂的目的消除大氣對太陽和來自目標的輻射產生吸取和散射做用的 影響，從而得到目標反射率、輻射率、地表溫度等真實物理模型參數。大多數情 況下，大氣校訂同時也是反演地物真實反射率的過程。

    PIE 的大氣校訂模塊是基於 6S 大氣輻射傳輸模型。 6S 模型假定無雲大氣的 狀況下，考慮了水汽、 CO二、 O3 和 O2 的吸取、分子和睦溶膠的散射以及非均一 地面和雙向反射率的問題。光譜積分的步長爲 2.5nm，能夠模擬機載觀測、設置 目標高程、解釋 BRDF 做用和臨近效應，增長了兩種吸取氣體的計算（ CO、N2O）。

    PIE支持算法功能的執行，下面對大氣校訂算法功能進行介紹。

2. 算法功能實現說明

2.1. 實現步驟

第一步	算法參數設置
第二步	算法執行
第三步	結果顯示

2.2. 算法參數

算法名稱	大氣校訂
C#算法DLL	PIE.CommonAlgo.dll
C#算法名稱	PIE.CommonAlgo.AtmosphericCorrectionAlgo
參數結構體	DataProcess_AtmCor_Info
參數說明
DataType	­int	數據類型 1——DN值 2——輻亮度 3——表觀反射率
InputFile	String	輸入文件（多光譜數據） (*.tif;*.tiff;*.bmp;*.img;*.jpg;*.ldf)
InputXML	String	輸入文件的元數據文件
OutputSR	String	輸出文件路徑 (輸出類型*.tif;*.tiff; *.img)
AtmModel	int	大氣模式(默認是0，就是系統會根據影像的中心經緯度和成像時間，自動肯定一種大氣模式) 0——系統自動選擇大氣模式 1——熱帶大氣模式 2——中緯度夏季大氣模式 3——中緯度冬季大氣模式 4——副極地夏季大氣模式 5——副極地冬季大氣模式 6——美國1962大氣模式
AerosolType	int	氣溶膠類型（必需設置成1-5，選擇一種氣溶膠類型） 1——大陸型氣溶膠 2——海洋型氣溶膠 3——城市型氣溶膠 4——沙塵型氣溶膠 5——煤煙型氣溶膠 6——平流層型氣溶膠
InitialVIS	float	初始能見度 默認爲40.0KM
FileTypeCode	String	根據輸出類型得到文件編碼類型 .tif/.tiff——GTiff .img—————HFA 其餘—————ENVI
AeroRetrieval	Int	是否逐像元反演氣溶膠(0表明否，1表明是，默認是1)
SatelliteID	Int	傳感器類型

2.3. 示例代碼

項目路徑	百度雲盤地址下/PIE示例程序/10.算法調用/圖像預處理/  PIE.ImagePreprocessing. AtmosphericCorrectionAlgo
數據路徑	百度雲盤地址下/PIE示例數據/ 柵格數據/01.GF1/GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.xml
視頻路徑	百度雲盤地址下/PIE視頻教程/10.算法調用/圖像預處理/大氣校訂算法.avi
示例代碼
1 /// <summary> 2 ///大氣校訂算法測試,本算法實現了將GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff進行大氣校訂 3 ///大氣模式爲系統自動選擇大氣模式，氣溶膠類型爲大陸性氣溶膠，初始能見度爲40km，逐項元反演氣溶膠爲是 4 /// </summary> 5 public override void OnClick() 6 { 7 #region 一、參數設置 8 PIE.CommonAlgo.DataProcess_AtmCor_Info info = new PIE.CommonAlgo.DataProcess_AtmCor_Info(); 9 10 info.InputFile = @"D:\Data\GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff"; 11 info.InputXML = @"D:\Data\GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.xml"; 12 info.OutputSR = @"D:\Data\result2.tif"; 13 info.AtmModel = 0; 14 info.AerosolType = 1; 15 info.InitialVIS = 40; 16 info.AeroRetrieval = 1; 17 info.FileTypeCode = "Gtiff"; 18 info.DataType = 1; 19 20 PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.AtmosphericCorrectionAlgo"); 21 if (algo == null) return; 22 #endregion 23 24 //二、算法執行 25 PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents; 26 algo.Name = "大氣校訂"; 27 algo.Params = info; 28 bool result = PIE.SystemAlgo.AlgoFactory.Instance().ExecuteAlgo(algo); 29 30 //三、結果顯示 31 ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\result2.tif"); 32 m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); 33 } View Code

2.4. 示例截圖