PIESDK距离分类和最大似然分类

  1 using PIE.Carto;
  2 using PIE.CommonAlgo;
  3 using PIE.Framework;
  4 using PIE.SystemAlgo;
  5 using System;
  6 using System.Collections.Generic;
  7 using System.Linq;
  8 using System.Text;
  9 
 10 namespace DistanceClassify
 11 {
 12     //距离分类
 13     class DistanceClassifyCommand : DesktopCommand
 14     {
 15         /// 

 16         /// 监督分类参数
 17         /// 
 18         private PIE.CommonAlgo.SupervisedClassification_Exchange_Info m_DataInfo = null;
 19 
 20         public DistanceClassifyCommand()
 21         {
 22             this.Name = "DistanceClassifyCommand";
 23             this.Caption = "DistanceClassifyCommand";
 24         }
 25 
 26         /// 

 27         /// 单击方法
 28         /// 
 29         public override void OnClick()
 30         {
 31             int selBandNums = 3;//选择分类数据的波段
 32             m_DataInfo = new PIE.CommonAlgo.SupervisedClassification_Exchange_Info();
 33             m_DataInfo.InputFilePath = @"D:\data\China1\World.tif";
 34             m_DataInfo.OutputFilePath = @"D:\data\temp\WorldDisClassify.tif";
 35             m_DataInfo.CStart = 0;//行列起止值
 36             m_DataInfo.CEnd = 2046;
 37             m_DataInfo.RStart = 0;
 38             m_DataInfo.REnd = 1022;
 39             m_DataInfo.SelBandNums = selBandNums;//波段的个数
 40             m_DataInfo.ClassifierType = 1;//分类类型 0 代表最大似然;1代表距离分类
 41 
 42             //多光谱波段的的集合
 43             List<int> listBandIndex = new List<int>();
 44             for (int i = 0; i )
 45             {
 46                 listBandIndex.Add(i + 1);
 47             }
 48             m_DataInfo.SelBandIndexs = listBandIndex;
 49 
 50             #region ROIStatistics ROI统计算法
 51             ROIStatistics_Exchange_Info roiDataInfo = new PIE.CommonAlgo.ROIStatistics_Exchange_Info();
 52             roiDataInfo.ClassifierType = m_DataInfo.ClassifierType;
 53             roiDataInfo.currentFileName = m_DataInfo.InputFilePath;
 54             roiDataInfo.selBandIndex = m_DataInfo.SelBandIndexs.ToArray();
 55             roiDataInfo.selBandNums = m_DataInfo.SelBandNums;
 56 
 57             IList listLayer = m_HookHelper.ActiveView.FocusMap.GetAllLayer();
 58             foreach (ILayer item in listLayer)
 59             {
 60                 if (item.Name == "roi_layer")//得到利用分类工具生成的roi文件
 61                 {
 62                     roiDataInfo.pROILayer = item as IGraphicsLayer;
 63                 }
 64             }
 65 //创建并执行roi统计算法
 66             ISystemAlgo algo = AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.ROIStatisticsAlgo");
 67             algo.Params = roiDataInfo;
 68             bool result = AlgoFactory.Instance().ExecuteAlgo(algo);//同步执行roi统计算法
 69             #endregion
 70 
 71             #region DistanceClassify算法
 72 //执行完成后得到roi统计的参数会发生变化
 73             roiDataInfo = (ROIStatistics_Exchange_Info)algo.Params;
 74            int roiNum = roiDataInfo.vecROIName.Count;
 75 
 76             m_DataInfo.ROINums = roiNum; //ROI个数
 77             m_DataInfo.ListRoiNames = roiDataInfo.vecROIName;//roi名称列表
 78             m_DataInfo.ListRoiColors = roiDataInfo.vecROIColor;//roi颜色集合
 79             m_DataInfo.ROIMean = new List<double>();
 80             m_DataInfo.ROICof = new List<double>();
 81             for (int i = 0; i //roi均值集合
 82             {
 83                 for (int j = 0; j )
 84                 {
 85                     m_DataInfo.ROIMean.Add(roiDataInfo.vecMean[i * selBandNums + j]);
 86                 }
 87                 for (int k = 0; k )
 88                 {
 89                     m_DataInfo.ROICof.Add(roiDataInfo.vecCof[i * selBandNums * selBandNums + k]);
 90                 }
 91             }
 92             ISystemAlgo distanceAlgo = AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.DistanceClassificationAlgo");//最大似然法就将DistanceClassificationAlgo替换为MLClassificationAlgo
 93             distanceAlgo.Params = m_DataInfo;
 94             ISystemAlgoEvents systemEvents = distanceAlgo as ISystemAlgoEvents;
 95             systemEvents.OnExecuteCompleted += systemEvents_OnExecuteCompleted;
 96             result = AlgoFactory.Instance().ExecuteAlgo(distanceAlgo);
 97             #endregion
 98         }
 99 
100         /// 

101         /// 算法执行完成事件
102         /// 
103         /// 
104         void systemEvents_OnExecuteCompleted(ISystemAlgo algo)
105         {
106             ISystemAlgoEvents systemEvents = algo as ISystemAlgoEvents;
107             systemEvents.OnExecuteCompleted -= systemEvents_OnExecuteCompleted;
108 ILayer layer = LayerFactory.CreateDefaultLayer(m_DataInfo.OutputFilePath);
109             if (layer == null)
110             {
111                 System.Windows.Forms.MessageBox.Show("分类后图层为空");
112                 return;
113             }
114             m_HookHelper.ActiveView.FocusMap.AddLayer(layer);
115             m_HookHelper.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);
116         }
117     }
118 }