PIE SDK主成分变换

1.算法功能简介

　  主成分变换(Principal Component Analysis，PCA)又称K-L（Karhunen-Loeve）变换或霍特林（Hotelling）变换，是基于变量之间的相关关系，在尽量不丢失信息前提下的一种线性变换的方法，主要用于数据压缩和信息增强。

　　主成分正变换，一般意义的K-L变换就是指正变换，该过程通过对图像进行统计，在波段协方差矩阵的基础上计算特征值，构造主成分。根据主成分与特征值的关系，可以选择少数的主成分作为输出结果。

　　主成分逆变换，如果在正变换中选择的主成分数目与波段/变量数目相同，那么逆变换结果将完全等同于原始影像。如果选择的主成分数目少于波段数逆变换结果相当于压抑了图像中的噪声。受选择的主成分数目的影响，逆变换结果图像的各个“波段”与原始图像波段可能会有较大的差异而不再具有原始图像波段的物理意义。

　　PIESDK提供了正变换和逆变换的算法，只需要设置对应的参数条件就可以执行，下面介绍下两种算法使用方法。

下面的示例代码需要安装DevExpress三方界面库

2.算法功能实现说明

2.1. 实现步骤

第一步	算法参数设置
第二步	算法执行
第三步	结果显示

2.2. 算法参数

算法名称	主成分正变换
C#算法DLL	PIE.CommonAlgo.dll
C#算法名称	PIE.CommonAlgo.TransformForwardPCAAlgo
参数结构体	ForwardPCA_Exchange_Info
参数说明
m_strInputFile	String	输入文件
m_strOutputResultFile	String	输出影像路径
m_strOutputStatsFile	String	输出统计文件
m_strFileTypeCode	String	输出文件类型
m_nPCBands	int	输出波段数量
m_nOutDataType	int	输出文件类型 0、字节型（8位）； 1、无符号整形（16位）； 2、整形（16位）； 3、无符号长整形（32位）； 4、长整形（32位）； 5、浮点型（32位）； 6、双精度浮点型（64位）
m_eigenvalues	IList<string>	特征值 -返回
m_bPCBandsFromEigenvalus	bool	根据特征值排序选择PCA波段
m_bOutputLikeEnvi	bool	零均值处理
m_bCovariance	bool	统计使用矩阵 -true- 使用协方差矩阵 -false- 使用相关系数矩阵
m_accumulate_contribute	IList<string>	百分比 -返回

算法名称	主成分逆变换
C#算法DLL	PIE.CommonAlgo.dll
C#算法名称	PIE.CommonAlgo.TransformInversePCAAlgo
参数结构体	InversePCA_Exchange_Info
参数说明
m_m_nOutDataType	int	输出文件字节类型 0、字节型（8位）； 1、无符号整形（16位）； 2、整形（16位）； 3、无符号长整形（32位）； 4、长整形（32位）； 5、浮点型（32位）； 6、双精度浮点型（64位）
m_strFileTypeCode	String	输出文件格式
m_strInputPcaFile	String	输入PCA结果文件
m_strInputStatsFile	String	输入PCA结果统计文件
m_strOutputResultFile	String	输出文件路径

2.3. 示例代码

项目路径	百度云盘地址下/PIE示例程序/10.算法调用/图像处理/ImageTransform
数据路径	百度云盘地址下/PIE示例数据/栅格数据/04.World/World.tif
视频路径	百度云盘地址下/PIE视频教程/10.算法调用/图像处理/主成分变换.avi
示例代码
/// <summary> /// 主成分正变换 /// </summary> /// <param name="sender"></param> /// <param name="e"></param> private void toolStripButton1_Click(object sender, EventArgs e) { //1、参数设置 PIE.CommonAlgo.ForwardPCA_Exchange_Info info = new ForwardPCA_Exchange_Info(); info.m_strInputFile = @"D:\data\02.测试数据\World\World.tif"; info.m_nOutDataType = ;//float32 info.m_strOutputResultFile = @"D:\PCPT.tif";//输出文件 info.m_strOutputStatsFile = @"D:\PCPT.pcasta";//输出统计文件 info.m_nPCBands = ;//输出的主成分波段数 info.m_strFileTypeCode = "GTiff";//文件格式类型 info.m_bOutputLikeEnvi = true;//零均值处理 info.m_bPCBandsFromEigenvalus = false;//是否根据特征值排序PCA波段（如果为true，m_nPCBands 为0个波段） info.m_bCovariance = true;//统计使用矩阵 true 协方差矩阵 false 使用相关系矩阵 前提是 bPCBandsFromEigenvalus为true，设置的才有效 //2、创建算法对象 ISystemAlgo algo = AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.TransformForwardPCAAlgo"); if (algo == null) return; algo.Params = info; //3、执行算法 bool result = AlgoFactory.Instance().ExecuteAlgo(algo); if (result) { if (info.m_nPCBands == && info.m_bPCBandsFromEigenvalus) { PIE.CommonAlgo.ForwardPCA_Exchange_Info info1 = algo.Params as PIE.CommonAlgo.ForwardPCA_Exchange_Info; PIE.Plugin.FrmPCABandSelect frm = new PIE.Plugin.FrmPCABandSelect(); frm.Init(info1.m_eigenvalues, info1.m_accumulate_contribute); if (frm.ShowDialog() != DialogResult.OK) return; info1.m_nPCBands = frm.nOutBandCount;//选择输出的波段号 algo.Params = info1; result = AlgoFactory.Instance().ExecuteAlgo(algo); } } if (result == false) return; mapControlMain.AddLayerFromFile(info.m_strOutputResultFile, ); mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); } /// <summary> /// 主成分逆变换 （对主成分正变换的结果做逆变换，得到正变换之前的影像数据） /// </summary> /// <param name="sender"></param> /// <param name="e"></param> private void toolStripButton2_Click(object sender, EventArgs e) { //1、参数设置 PIE.CommonAlgo.InversePCA_Exchange_Info info = new InversePCA_Exchange_Info(); info.m_strInputPcaFile = @"D:\PCPT.tif";//输入文件 主成分正变换结果 info.m_strInputStatsFile = @"D:\PCPT.pcasta";//pac统计文件 info.m_strOutputResultFile = @"D:\InversePC.tif";//逆变换结果 info.m_strFileTypeCode = "GTiff";//输出结果 info.m_nOutDataType = ;//float32 //2、创建算法对象 ISystemAlgo algo = AlgoFactory.Instance().CreateAlgo("PIE.CommonAlgo.dll", "PIE.CommonAlgo.TransformInversePCAAlgo"); if (algo == null) return; algo.Params = info; //3、执行算法并加载结果图层 bool result = AlgoFactory.Instance().ExecuteAlgo(algo); if (result != true) return; ILayer layer = LayerFactory.CreateDefaultLayer(info.m_strOutputResultFile); mapControlMain.ActiveView.FocusMap.AddLayer(layer); mapControlMain.ActiveView.PartialRefresh(ViewDrawPhaseType.ViewAll); }

2.4. 示例截图

图一、主成分正变换

图二：主成分逆变换