KDE核密度估计(使用Ep函数)
2019-10-13 本文已影响0人
一技破万法
随着计算机视觉、模式识别、机器人技术和自动化技术的发展，视觉图像系统在各个领域都得到了广泛的应用，而应用的前提首先是对环境中的特定目标做出准确检测。本文介绍了基于Epanechnikov核函数的核密度估计方法进行初步的前景检测。
什么是核密度估计
核密度估计就是估计概率密度函数，从几个样本估计其服从的分布，即求出其概率密度函数，这样就可以求任意区间处的概率了。所以核密度就是一个从具体样本到普遍概率密度的过程，然后再用普遍指导具体的问题。
核密度估计函数的公式如下表示：
核密度估计公式
其中Kh为核函数，Ｈ为d*d大小的正定带宽矩阵。当简化带宽矩阵为单参数ｈ时，则公式如下：
核密度简化带宽公式
对于多维变量的核函数可以通过每个维度下单变量和函数的乘积来表示。
使用Ep函数
在ｉ时刻时，我们将位于某坐标(x,y)处的像素点与数据集中所有图像在该坐标点的像素点带入密度函数的核密度估计中去，可得如下公式：
核密度估计公式
其中Ｋ为选择的核函数，满足性质Ｋ(t)>=0和K对时间t的积分为１．
由于图像帧一般是多通道彩色图像，所以样本点都是多维向量，将d维颜色通道的核函数用每个色彩通道的核函数乘积(上文提到)表示，如下图公式所示：
核密度估计公式
单维度Epanechnikov函数表示如下：
单维度Ep函数
其中Ｃｄ表示ｄ维超球体体积，对于三维来说Ｃ３就是半径为１的球体体积。此时像素输入Ｘｔ，样本点为Ｘｉ，基于Ep函数的三通道彩色图像像素值的概率密度估计函数为：
概率密度估计函数
我们指定Th为一个小鱼１的阈值与计算得到的估计值进行比较决策，低于此阈值的像素点为褐色表示前景，反之为白色表示背景。当所有像素点都通过决策之后，便生成了一副二值检测图像。
比较决策
最后附上效果图：
原图像
三维KDE估计图
二值化图像
本程序是基于python 3.5.2,使用的库有numpy|math|matplotlib|opencv|mpl_toolkits等
如果有不明白的可以关注作者微信公众号并私信问题：
一技破万法
