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 |
示例代码 |
|
1 /// <summary> 2 /// 同态滤波算法测试,本算法实现了将World.tif进行同态滤波,滤波类型为巴特沃斯高通变换,截止频率为50 3 /// 低频增益:0.25,高频增益:2,常数或阶数:1 4 /// </summary> 5 public override void OnClick() 6 { 7 #region 1、参数设置 8 PIE.CommonAlgo.StImageHomoInfo info = new PIE.CommonAlgo.StImageHomoInfo(); 9 10 info.InputFilePath = @"D:DataWorld.tif"; 11 info.OutputFilePath = @"D:Dataip_result17.tif"; 12 info.LowGains = 0.25F; 13 info.HighGains = 2; 14 info.HighPass = 0; 15 info.CutFreq = 50; 16 info.Constant = 1; 17 info.FileTypeCode = "GTiff"; 18 info.LowBands = new List<int> { 0, 1, 2 }; 19 20 PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.ImgProFiltHomoAlgo"); 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:Dataip_result17.tif"); 32 m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); 33 } |