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 /// <summary> 2 ///PCA融合算法測試,本算法實現了將PAN.tif和MS.tif進行最近鄰域法PCA融合 3 /// </summary> 4 public override void OnClick() 5 { 6 #region 一、參數設置 7 PIE.CommonAlgo.Pansharp_Exchange_Info info = new PIE.CommonAlgo.Pansharp_Exchange_Info(); 8 9 info.PanFilePath = @"D:\Data\PAN.tif"; 10 info.MssFilePath = @"D:\Data\MS.tif"; 11 info.MULChannels = new List<int> { 0, 1, 2 }; 12 info.HighChannel = 0; 13 info.OutputFilePath = @"D:\Data\result8.tif"; 14 info.FileTypeCode = "Gtiff"; 15 info.ResampleMode = 0; 16 info.AlgoType = 0; //PCA融合 17 //info.AlgoType = 1; //pansharp 融合 18 info.BMultiThread = true; 19 20 PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.PansharpFuseAlgo"); 21 if (algo == null) return; 22 #endregion 23 24 //二、算法執行 25 PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents; 26 algo.Name = "PCA融合"; 27 algo.Params = info; 28 algo.Params = info; 29 (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnExecuteCompleted +=PansharpFuseAlgo_OnExecuteCompleted; 30 (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnProgressChanged += PansharpFuseAlgo_OnProgressChanged; 31 32 PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar; 33 statusBar.ShowProgress(0, 100, ""); 34 PIE.SystemAlgo.AlgoFactory.Instance().AsynExecuteAlgo(algo); 35 } 36 int PansharpFuseAlgo_OnProgressChanged(double complete, string msg, SystemAlgo.ISystemAlgo algo) 37 { 38 PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar; 39 statusBar.UpdateProgress((int)complete, msg); 40 return 0; 41 } 42 43 void PansharpFuseAlgo_OnExecuteCompleted(SystemAlgo.ISystemAlgo algo) 44 { 45 PIE.AxControls.IStatusBar statusBar = m_Application.StatusBar; 46 statusBar.UpdateProgress(100, ""); 47 statusBar.HideProgress(); 48 //三、結果顯示 49 ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\result8.tif"); 50 m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); 51 (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnExecuteCompleted -= PansharpFuseAlgo_OnExecuteCompleted; 52 (algo as PIE.SystemAlgo.ISystemAlgoEvents).OnProgressChanged -= PansharpFuseAlgo_OnProgressChanged; 53 } View Code

2.4.示例截圖