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

示例代码

         /// <summary>
///大气校正算法测试,本算法实现了将GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff进行大气校正
///大气模式为系统自动选择大气模式,气溶胶类型为大陆性气溶胶,初始能见度为40km,逐项元反演气溶胶为是
/// </summary>
public override void OnClick()
{
#region 1、参数设置
PIE.CommonAlgo.DataProcess_AtmCor_Info info = new PIE.CommonAlgo.DataProcess_AtmCor_Info(); info.InputFile = @"D:\Data\GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff";
info.InputXML = @"D:\Data\GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.xml";
info.OutputSR = @"D:\Data\result2.tif";
info.AtmModel = ;
info.AerosolType = ;
info.InitialVIS = ;
info.AeroRetrieval = ;
info.FileTypeCode = "Gtiff";
info.DataType = ; PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.AtmosphericCorrectionAlgo");
if (algo == null) return;
#endregion //2、算法执行
PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents;
algo.Name = "大气校正";
algo.Params = info;
bool result = PIE.SystemAlgo.AlgoFactory.Instance().ExecuteAlgo(algo); //3、结果显示
ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\result2.tif");
m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);
}

2.4. 示例截图

PIE SDK大气校正-LMLPHP

05-28 07:37