PIE SDK PCA融合

1.算法功能简介

    PCA 融合分三步实现，首先将多光谱数据进行主成分变换，然后用高分辨单波段替换第一主成分波段，最后进行主成份逆变换得到融合图像。

    PIE支持算法功能的执行，下面对PCA融合算法功能进行介绍。

2.算法功能实现说明

2.1. 实现步骤

第一步

算法参数设置

第二步

算法执行

第三步

结果显示

2.2. 算法参数

算法名称

PCA融合

C#算法DLL

PIE.CommonAlgo.dll

C#算法名称

PIE.CommonAlgo.PansharpFuseAlgo

参数结构体

Pansharp_Exchange_Info

参数说明

PanFilePath

String

输入高分辨率数据路径

MssFilePath

String

输入多光谱数据路径

MULChannels

IList<int>

输入的多光谱数据波段集合{0,1,2,3}

HighChannel

Int

融合选择的高分辨率数据波段

OutputFilePath

String

输出文件

(*.tif;*.tiff; *.img)

FileTypeCode

String

根据输出类型获得文件编码类型

.tif/.tiff——GTiff

.img—————HFA

其他—————ENVI

ResampleMode

Int

重采样类型

(最近邻点插值0/双线性插值1/三次卷积法插值2)

AlgoType

Int

融合类型

0代表pca融合 、1代表pansharp 融合

FuncName

String

功能名称

BMultiThread

Bool

是否多线程(true)

TRUE 是； FALSE 否

2.3. 示例代码

项目路径

百度云盘地址下/PIE示例程序/ 10.算法调用/图像预处理/PIE.ImagePreprocessing. PansharpFuseAlgo

       数据路径 

百度云盘地址下/PIE示例数据/栅格数据/01.GF1/PAN.tif ;MS.tif

视频路径

百度云盘地址下/PIE视频教程/10.算法调用/图像预处理/ PCA融合算法.avi

示例代码

1 1          /// <summary>
2 2         ///PCA融合算法测试,本算法实现了将PAN.tif和MS.tif进行最近邻域法PCA融合
3 3         /// </summary>
4 4         public override void OnClick()
5 5         {
6 6             #region 1、参数设置
7 7             PIE.CommonAlgo.Pansharp_Exchange_Info info = new PIE.CommonAlgo.Pansharp_Exchange_Info();
8 8 
9 9             info.PanFilePath = @"D:\Data\PAN.tif";
1010             info.MssFilePath = @"D:\Data\MS.tif";
1111             info.MULChannels = new List<int> { 0, 1, 2 };
1212             info.HighChannel = 0;
1313             info.OutputFilePath = @"D:\Data\result8.tif";
1414             info.FileTypeCode = "Gtiff";
1515             info.ResampleMode = 0;
1616             info.AlgoType = 0; //PCA融合
1717             //info.AlgoType = 1; //pansharp 融合
1818             info.BMultiThread = true;
1919 
2020             PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.PansharpFuseAlgo");
2121             if (algo == null) return;
2222             #endregion
2323 
2424             //2、算法执行
2525             PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents;
2626             algo.Name = "PCA融合";
2727             algo.Params = info;
2828             algo.Params = info;
2929             (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnExecuteCompleted +=PansharpFuseAlgo_OnExecuteCompleted;
3030             (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnProgressChanged += PansharpFuseAlgo_OnProgressChanged;
3131 
3232             PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar;
3333             statusBar.ShowProgress(0, 100, "");
3434             PIE.SystemAlgo.AlgoFactory.Instance().AsynExecuteAlgo(algo);
3535         }
3636         int PansharpFuseAlgo_OnProgressChanged(double complete, string msg, SystemAlgo.ISystemAlgo algo)
3737         {
3838             PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar;
3939             statusBar.UpdateProgress((int)complete, msg);
4040             return 0;
4141         }
4242 
4343         void PansharpFuseAlgo_OnExecuteCompleted(SystemAlgo.ISystemAlgo algo)
4444         {
4545             PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar;
4646             statusBar.UpdateProgress(100, "");
4747             statusBar.HideProgress();
4848             //3、结果显示
4949             ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\result8.tif");
5050             m_HookHelper.ActiveView.FocusMap.AddLayer(layer);            m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);          m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);
5151             (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnExecuteCompleted -= PansharpFuseAlgo_OnExecuteCompleted;
5252             (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnProgressChanged -= PansharpFuseAlgo_OnProgressChanged;
5353         }

View Code

2.4.示例截图