PIE SDK主成分变换

 1 /// <summary>
 2 /// 主成分正变换
 3 /// </summary>
 4 /// <param name="sender"></param>
 5 /// <param name="e"></param>
 6 private void toolStripButton1_Click(object sender, EventArgs e)
 7 {
 8     //1、参数设置
 9     PIE.CommonAlgo.ForwardPCA_Exchange_Info info = new ForwardPCA_Exchange_Info();
10     info.m_strInputFile = @"D:\data\02.测试数据\World\World.tif";
11     info.m_nOutDataType = 5;//float32
12     info.m_strOutputResultFile = @"D:\PCPT.tif";//输出文件
13     info.m_strOutputStatsFile = @"D:\PCPT.pcasta";//输出统计文件
14     info.m_nPCBands = 3;//输出的主成分波段数
15     info.m_strFileTypeCode = "GTiff";//文件格式类型        
16     info.m_bOutputLikeEnvi = true;//零均值处理
17     info.m_bPCBandsFromEigenvalus = false;//是否根据特征值排序PCA波段
18     info.m_bCovariance = false;//统计使用矩阵 true 协方差矩阵 false 使用相关系矩阵  前提是
19 bPCBandsFromEigenvalus为true，设置的才有效
20     //2、创建算法对象
21     ISystemAlgo algo = AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.TransformForwardPCAAlgo");
22     if (algo == null) return;
23     algo.Params = info;
24 
25     //3、执行算法
26     bool result = AlgoFactory.Instance().ExecuteAlgo(algo);
27     if (result != true) return;
28     ILayer layer = LayerFactory.CreateDefaultLayer(info.m_strOutputResultFile);
29     mapControlMain.ActiveView.FocusMap.AddLayer(layer);
30     mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);         
31 } 
32 
33 
34 /// <summary>
35 /// 主成分逆变换 （对主成分正变换的结果做逆变换，得到正变换之前的影像数据）
36 /// </summary>
37 /// <param name="sender"></param>
38 /// <param name="e"></param>
39 private void toolStripButton2_Click(object sender, EventArgs e)
40 {
41     //1、参数设置
42     PIE.CommonAlgo.InversePCA_Exchange_Info info = new InversePCA_Exchange_Info();
43     info.m_strInputPcaFile = @"D:\PCPT.tif";//输入文件 主成分正变换结果
44     info.m_strInputStatsFile = @"D:\PCPT.pcasta";//pac统计文件
45     info.m_strOutputResultFile = @"D:\InversePC.tif";//逆变换结果
46     info.m_strFileTypeCode = "GTiff";//输出结果
47     info.m_nOutDataType = 5;//float32
48 
49     //2、创建算法对象
50     ISystemAlgo algo = AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.TransformInversePCAAlgo");
51     if (algo == null) return;
52     algo.Params = info;
53 
54     //3、执行算法并加载结果图层
55     bool result = AlgoFactory.Instance().ExecuteAlgo(algo);
56     if (result != true) return;
57     ILayer layer = LayerFactory.CreateDefaultLayer(info.m_strOutputResultFile);
58     mapControlMain.ActiveView.FocusMap.AddLayer(layer);
59     mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll);      
60 }