PIE SDK同态滤波

1.算法功能简介

同态滤波是减少低频增加高频，从而减少光照变化并锐化边缘或细节的图像滤波方法。

同态滤波的流程为：空间域图像→对数运算→傅里叶正变换→同态滤波――傅里叶逆变换→指数运算→同态滤波结果。不同空间分辨率的遥感图像，使用同态滤波的效果不同。如果图像中的光照可以认为是均匀的，那么，进行同态滤波产生的效果不大。但是，如果光照明显是不均匀的，那么同态滤波有助于表现出图像中暗处的细节。

PIE SDK支持算法功能的执行，下面对同态滤波算法功能进行介绍。

2.算法功能实现说明

2.1. 实现步骤

第一步	算法参数设置
第二步	算法执行
第三步	结果显示

2.2. 算法参数

算法名称	同态滤波
C#算法DLL	PIE.CommonAlgo.dll
C#算法名称	PIE.CommonAlgo.ImgProFiltHomoAlgo
参数结构体	StImageHomoInfo
参数说明
InputFilePath	String	输入文件 (*.tif;*.tiff; *.img)
OutputFilePath	String	输出文件路径 (*.tif;*.tiff; *.img)
XMLFile	String	XML文件路径 (*.xml)
LowGains	Float	低频增益：0.25
HighGains	Float	高频增益：2
HighPass	Int	滤波类型（巴特沃斯高通变换：0，高斯高通变换：1）
CutFreq	Float	截止频率：50
Constant	Float	常数或阶数：1
FuncName	String	功能名称
FileTypeCode	String	根据输出类型获得文件编码类型 .tif/.tiff——GTiff .img—————HFA 其他—————ENVI
LowBands	IList<Int>	输出影像的波段（至少选择一个波段，{ 0, 1, 2, 3 }）

2.3. 示例代码

项目路径	百度云盘地址下/PIE示例程序/10.算法调用/图像处理/ ImageProcessing. ImgProFiltHomoAlgo
数据路径	百度云盘地址下/PIE示例数据/栅格数据/04.World/World.tif
视频路径	百度云盘地址下/PIE视频教程/10.算法调用/图像处理/同态滤波算法.avi
示例代码
/// <summary> /// 同态滤波算法测试,本算法实现了将World.tif进行同态滤波，滤波类型为巴特沃斯高通变换,截止频率为50 /// 低频增益：0.25,高频增益：2,常数或阶数：1 /// </summary> public override void OnClick() { #region 1、参数设置 PIE.CommonAlgo.StImageHomoInfo info = new PIE.CommonAlgo.StImageHomoInfo(); info.InputFilePath = @"D:\Data\World.tif"; info.OutputFilePath = @"D:\Data\ip_result17.tif"; info.LowGains = 0.25F; info.HighGains = ; info.HighPass = ; info.CutFreq = ; info.Constant = ; info.FileTypeCode = "GTiff"; info.LowBands = new List<int> { , , }; PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.ImgProFiltHomoAlgo"); 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\ip_result17.tif"); m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); }

2.4. 示例截图