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
示例代码
/// <summary> ///PCA融合算法测试,本算法实现了将PAN.tif和MS.tif进行最近邻域法PCA融合 /// </summary> public override void OnClick() { #region 1、参数设置 PIE.CommonAlgo.Pansharp_Exchange_Info info = new PIE.CommonAlgo.Pansharp_Exchange_Info(); info.PanFilePath = @"D:\Data\PAN.tif"; info.MssFilePath = @"D:\Data\MS.tif"; info.MULChannels = new List<int> { , , }; info.HighChannel = ; info.OutputFilePath = @"D:\Data\result8.tif"; info.FileTypeCode = "Gtiff"; info.ResampleMode = ; info.AlgoType = ; //PCA融合 //info.AlgoType = 1; //pansharp 融合 info.BMultiThread = true; PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.PansharpFuseAlgo"); if (algo == null) return; #endregion //2、算法执行 PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents; algo.Name = "PCA融合"; algo.Params = info; algo.Params = info; (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnExecuteCompleted +=PansharpFuseAlgo_OnExecuteCompleted; (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnProgressChanged += PansharpFuseAlgo_OnProgressChanged; PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar; statusBar.ShowProgress(, , ""); PIE.SystemAlgo.AlgoFactory.Instance().AsynExecuteAlgo(algo); } int PansharpFuseAlgo_OnProgressChanged(double complete, string msg, SystemAlgo.ISystemAlgo algo) { PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar; statusBar.UpdateProgress((int)complete, msg); return ; } void PansharpFuseAlgo_OnExecuteCompleted(SystemAlgo.ISystemAlgo algo) { PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar; statusBar.UpdateProgress(, ""); statusBar.HideProgress(); //3、结果显示 ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\result8.tif"); m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnExecuteCompleted -= PansharpFuseAlgo_OnExecuteCompleted; (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnProgressChanged -= PansharpFuseAlgo_OnProgressChanged; }

2.4.示例截图