PIE SDK大气校正
1. 算法功能简介
大气校正的目的消除大气对太阳和来自目标的辐射产生吸收和散射作用的 影响,从而获得目标反射率、辐射率、地表温度等真实物理模型参数。大多数情 况下,大气校正同时也是反演地物真实反射率的过程。
PIE 的大气校正模块是基于 6S 大气辐射传输模型。 6S 模型假定无云大气的 情况下,考虑了水汽、 CO2、 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 1、参数设置 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 //2、算法执行 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 //3、结果显示 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 } |
2.4. 示例截图