VB.NET + EmguCV实现PCA降维_vbemg-程序员宅基地

技术标签：主成分分析 VB.NET 图像处理相关 OpenCV PCA Emgu PCA降维 openCV PCA Emgu PCA

PCA。Principal Component Analysis。主成分分析。

顾名思义，就是要从多元事物中解析出主要的影响因素。揭示事物的本质，简化复杂的问题。

PCA的目的是将高维的数据通过线性变换投影到较低维空间。

投影后的值应该尽可能的分散，才能尽量多的保留原始的信息。在数学上可以用最大方差来表示。

投影的低维空间应该相互之间不存在相关性。在数学上可以用协方差来表示相关性。

维度内的方差与维度间的协方差都可以表示为内积的形式，通过数学推导，可以得出我们想要的方差最大和协方差矩阵为零可以等价为将协方差矩阵对角化。具体推导及原理可参考
http://blog.codinglabs.org/articles/pca-tutorial.html

算法步骤如下：

设有m条n维数据。想要将数据降维到k维

①将原始数据按列组成n行m列矩阵X

②将X的每一行（代表一个属性字段）进行零均值化，即减去这一行的均值

③求出协方差矩阵C = 1/m * X （
X转置）

④求出协方差矩阵的特征值及对应的特征向量

⑤将特征向量按对应特征值大小从上到下按行排列成矩阵，取前k行组成矩阵P

⑥Y = P X 即为降维到k维后的数据

VB.NET + EmguCV3.0的代码实现如下：

样本X是500幅29*29的图片

Dim n = 29 * 29
Dim m = 500
Dim X As Emgu.CV.Matrix(Of Double) = New Matrix(Of Double)(m, n)
Dim fileCollection As Collections.ObjectModel.ReadOnlyCollection(Of String) = FileIO.FileSystem.GetFiles(Application.StartupPath & "\img")
For i = 0 To fileCollection.Count - 1
	Dim img As Image(Of Gray, Byte) = New Image(Of Gray, Byte)(fileCollection(i))
	Dim j As Integer = 0
	For x1 = 0 To img.Width - 1
		For y1 = 0 To img.Height - 1
			X.Data(i, j) = img.Data(y1, x1, 0)
			j = j + 1
		Next
	Next
Next

计算每一维度上的均值，进行零均值化

'计算每一维度上的均值
Dim Xxmean(n - 1) As Double
For j = 0 To n - 1
	For i = 0 To m - 1
		Xxmean(j) = Xxmean(j) + X.Data(i, j)
	Next
	Xxmean(j) = Xxmean(j) / (m - 1)
Next
'每一维度上零均值化
For j = 0 To n - 1
	For i = 0 To 499
		X.Data(i, j) = X.Data(i, j) - Xxmean(j)
	Next
Next

计算协方差矩阵

'协方差矩阵C = 1/mXXt
Dim C As Emgu.CV.Matrix(Of Double) = New Matrix(Of Double)(n, n)
For i = 0 To n - 1
	For j = 0 To n - 1
		Dim sum As Double = 0
		For z = 0 To m - 1
			sum = sum + X.Data(z, i) * X.Data(z, j)
		Next
		C.Data(i, j) = sum / m
	Next
Next

计算对称矩阵的特征值和特征向量

因为协方差矩阵是对称阵，可以通过雅可比过关法来求特征值与特征向量。本文中并未实现该方法，而是通过EmguCv中的Eigen方法来求特征值与特征向量。

'特征向量，按行排列~
Dim eigenVector As Emgu.CV.Matrix(Of Double) = New Matrix(Of Double)(n, n)
'特征值输出矩阵，按降序排列。特征值与特征向量的排序是一致的
Dim eigenValues As Emgu.CV.Matrix(Of Double) = New Matrix(Of Double)(n, 1)
CvInvoke.Eigen(C, eigenValues, eigenVector)

特征值降序排列，取前k行，并与样本X相乘，求得降维后的结果

'k为降维后的维度
Dim k As Integer = 99
'取前k行组成P
Dim P As Emgu.CV.Matrix(Of Double) = New Matrix(Of Double)(n, k)
For i = 0 To k - 1
	For j = 0 To n - 1
		P.Data(j, i) = eigenVector.Data(j, i)
	Next
Next
'Y = PX即为降维到k维后的输出
Dim Y As Emgu.CV.Matrix(Of Double) = New Matrix(Of Double)(k, m)
For i = 0 To k - 1
	For j = 0 To m - 1
		Dim sum As Double = 0.0
		For x1 = 0 To n - 1
			sum = sum + P.Data(x1, i) * X.Data(j, x1)
		Next
		Y.Data(i, j) = sum
	Next
Next

http://blog.csdn.net/u013162930/article/details/77688165

本文链接：https://blog.csdn.net/u013162930/article/details/77688165

原作者删帖不实内容删帖广告或垃圾文章投诉

智能推荐

使用 arm-linux-androideabi-addr2line 工具定位 libunity.so 崩溃问题-程序员宅基地

文章浏览阅读710次，点赞13次，收藏7次。它的名称有点不同 - aarch64-linux-android-addr2line。尽管该实用程序的名称不包含单词arm，但它成功地解密了arm64-v8下的堆栈跟踪，并通常感知arm64-v8的字符。是 NDK 自带的调试工具，可以用来分析 so 崩溃时输出的的内存地址。之后就是通过 cmd 进入到这个路径。找到了 64 位所需的实用程序。_arm-linux-androideabi

javaweb-邮件发送_javaweb发送邮件-程序员宅基地

javaweb-邮件发送摘要: 本文介绍了邮件传输协议(SMTP和POP3)以及电子邮件的发送和接收过程。还讨论了纯文本文件、带图片和附件的邮件发送方法，以及通过servlet方式注册邮箱和使用springboot框架发送邮件的实现。

element-ui table 设置表格滚动条位置_element table 滚动条位置-程序员宅基地

文章浏览阅读4.3k次，点赞6次，收藏11次。在切换不同页面时（被 keep-alive 缓存的组件间切换），页面中的element-ui table的滚动条位置没有停留在原来的位置。目前需要切换不同的页面返回来后，滚动条保持在原来的位置。_element table 滚动条位置

前端开发经验总结_属性值[session.getattribute("strpath")]引用["],在值内使用时必须-程序员宅基地

文章浏览阅读2.6k次。我设置nowrap和不设置nowrap效果都一样。就是表格随着文字自动伸展，为什么？ →回答问题：TD元素noWrap属性的行为与TD元素的width属性有关。 ◆如果未设置TD宽度，则noWrap属性是起作用的。◆如果设置了TD宽度，则noWrap属性是不起作用的。 http://www.blueidea.com/tech/web/2003/943.as_属性值[session.getattribute("strpath")]引用["],在值内使用时必须被转义。

JS如何把Object对象的数据输出到控制台中_前端怎么通过控制台查看字段取值-程序员宅基地

文章浏览阅读3.4k次，点赞2次，收藏3次。前端时间在编写程序时遇到这样一个问题，即在前端页面通过一些js框架自带的异步请求返回的数据类型为Object数据类型，笔者根据网上查阅到的资料，找到以下这种简单的方式：//把Object类型转为Json数据格式，再通过console命令在控制台中打印出来console.log("xhr的值为："+JSON.st..._前端怎么通过控制台查看字段取值

8.cc.Button组件使用详解_cc button.start-程序员宅基地

文章浏览阅读556次。1. cc.Button添加按钮的方法 2种方式 (1)直接创建带Button组件的节点; (2) 先创建节点，再添加组件;按钮组件, 按钮是游戏中最常用的组件, 点击然后响应事件;按钮的过渡效果: 过渡: 普通状态, 鼠标滑动到物体上, 按下状态, 禁用状态 (1)没有过渡，只有响应事件; (2)颜色过渡, 过渡效果中使用颜色; (3)精灵..._cc button.start