PIE SDK ISODATA分類

 

1.算法功能簡介

    ISODATA（IterativeSelf-OrganizingDataAnalysisTechniqueAlgorithm）即迭代式自組織數據分析技術， 其大體原理是首先計算數據空間中均勻分佈的類均值， 而後用最小距離規則將剩餘的像元進行迭代聚合；每次迭代都從新計算均值，且 根據所得的新均值，對像元進行再分類；這一處理過程持續到每一類的像元數變化少於所選的像元變化閥值或者達到了迭代的最大次數。

    ISODATA 算法經過設置初始參數而引入人機對話環節，並使用歸併和分裂 等機制，當兩類聚中心小於某個閥值時，將它們合併爲一類。當某類的標準差大 於某一閥值時或其樣本數目超過某一閥值時，將其分裂爲兩類，在某類樣本數目 小於某一閥值時，將其取消。這樣根據初始類聚中心和設定的類別數目等參數迭代，最終獲得一個比較理想的分類結果。 ISODATA 算法是一種經常使用的聚類分析方法，是一種非監督學習方法。

    PIE SDK支持算法功能的執行，下面對ISODATA分類算法功能進行介紹。

2.算法功能實現說明

2.1. 實現步驟

第一步	算法參數設置
第二步	算法執行
第三步	結果顯示

2.2算法參數

算法名稱	ISODATA分類
C#算法DLL	PIE.CommonAlgo.dll
C#算法名稱	PIE.CommonAlgo.ISODataClassificationAlgo
參數結構體	ISODataClassification_Exchange_Info
參數說明
InputFilePath	String	輸入文件 (*.tif;*.tiff; *.img)
OutputFilePath	String	輸出文件路徑 (*.tif;*.tiff; *.img)
ProspClassNum	Int	預期類數（2~255的正整數，默認爲8）
InitClassNum	Int	初始類數（2~255的正整數，默認爲5）
MinSam	Int	最少像元數（大於0的正整數，默認爲5）
MaxMerge	Int	最大合併對數（大於0的正整數，，默認爲1）
MaxLoop	Int	最大迭代次數（大於0的正整數，默認爲5）
Dev	double	最大標準差（正數，默認爲9.8）
MinDis	double	最小中心距離（正數，默認爲6.4）
FuncName	String	功能名稱
FileTypeCode	String	根據輸出類型得到文件編碼類型 .tif/.tiff——GTiff .img—————HFA 其餘—————ENVI
ListBands	IList<Int>	輸入影像的波段（至少選擇一個波段，{ 0, 1, 2, 3 }）

2.3. 示例代碼

項目路徑	百度雲盤地址下/PIE示例程序/10.算法調用/圖像處理/ ImageProcessing. ISODataClassificationAlgo
數據路徑	百度雲盤地址下/PIE示例數據/柵格數據/01.GF1/GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff
視頻路徑	百度雲盤地址下/PIE視頻教程/10.算法調用/圖像處理/ISODATA分類算法.avi
示例代碼
1 /// <summary> 2 ///ISODATA分類算法測試,本算法實現了將GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff進行ISODATA分類 3 ///預期類數爲8，初始類數爲5，最少像元數爲5，最大合併對數爲1，最大迭代次數爲5，最大標準差爲9.8，最小中心距離爲6.4 4 /// </summary> 5 public override void OnClick() 6 { 7 #region 一、參數設置 8 PIE.CommonAlgo.ISODataClassification_Exchange_Info info = new PIE.CommonAlgo.ISODataClassification_Exchange_Info(); 9 10 info.InputFilePath = @"D:\Data\GF1_PMS1_E116.5_N39.4_20131127_L1A0000117600-MSS1.tiff"; 11 info.OutputFilePath = @"D:\Data\ip_result1.tif"; 12 info.ProspClassNum = 8; 13 info.InitClassNum = 5; 14 info.MinSam = 5; 15 info.MaxLoop=5; 16 info.MaxMerge = 1; 17 info.Dev = 9.8; 18 info.MinDis = 6.4; 19 info.LowBands = new List<int> { 0, 1, 2, 3 }; 20 21 info.FileTypeCode = "Gtiff"; 22 23 PIE.SystemAlgo.ISystemAlgo algo = PIE.SystemAlgo.AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.ISODataClassificationAlgo"); 24 if (algo == null) return; 25 #endregion 26 27 //二、算法執行 28 PIE.SystemAlgo.ISystemAlgoEvents algoEvents = algo as PIE.SystemAlgo.ISystemAlgoEvents; 29 algo.Name = "ISODATA分類"; 30 algo.Params = info; 31 bool result = PIE.SystemAlgo.AlgoFactory.Instance().ExecuteAlgo(algo); 32 33 //三、結果顯示 34 ILayer layer = PIE.Carto.LayerFactory.CreateDefaultLayer(@"D:\Data\ip_result1.tif"); 35 m_HookHelper.ActiveView.FocusMap.AddLayer(layer); m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); 36 } View Code

2.4. 示例截圖