[发明专利]中值滤波方法及装置有效
申请号: | 201810336653.5 | 申请日: | 2018-04-16 |
公开(公告)号: | CN108573478B | 公开(公告)日: | 2019-12-20 |
发明(设计)人: | 王行;李骊;盛赞;周晓军;李朔;杨淼 | 申请(专利权)人: | 北京华捷艾米科技有限公司 |
主分类号: | G06T5/00 | 分类号: | G06T5/00 |
代理公司: | 32200 南京经纬专利商标代理有限公司 | 代理人: | 李明 |
地址: | 100080 北京市*** | 国省代码: | 北京;11 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 像素 中值滤波 中位 邻域像素 图像 加权 衰减因子 相似度 求和 滤波 边缘保持 边缘模糊 扩散效应 优化图像 中心像素 传统的 去噪 噪声 查找 联合 | ||
1.一种中值滤波方法,其特征在于,所述方法包括:
计算第一图像中邻域像素与中心像素的第一类像素值相似度和距离衰减因子值;
根据所述第一类像素值相似度和所述距离衰减因子值计算所述邻域像素的第一类权值的中位值;
查找所述第一类权值的中位值对应的第二类像素值的中位值;
将第二图像中邻域像素的第二类像素值中与第二类像素值的中位值接近的部分像素的第二像素值加权求和,得到加权求和值;
所述根据所述第一类像素值相似度和所述距离衰减因子值计算所述邻域像素的第一类权值的中位值,包括:
根据以下公式计算所述邻域像素的第一类权值:w(p,q)=s(p,q)r(p,q),q∈N(p);其中,s(p,q)代表第一类像素值相似度,r(p,q)代表距离衰减因子值;
根据以下公式计算所述第一类权值的中位值:
其中,α是常数系数,且0<α<1,N(p)代表p像素的邻域像素。
2.根据权利要求1所述的方法,其特征在于,采用以下公式计算第一图像中的邻域像素与中心像素的第一类像素值相似度:
s(p,q)=f(|I(p)-I(q)|),其中,f随|I(p)-I(q)|单调递减;
其中,p表示当前计算的像素位置;q代表p像素的邻域像素位置;I(p)代表p像素的第一类像素值;I(q)代表q像素的第一类像素值;
和/或,根据以下公式计算所述距离衰减因子:r(p,q)=g(|p-q|),其中,g随|p-q|单调递减。
3.根据权利要求1所述的方法,其特征在于,采用以下任一方式查找第一类权值的中位值对应的第二类像素值的中位值:二分法、排序法或者累计直方图法。
4.根据权利要求1至3任一项所述的方法,其特征在于,所述查找所述第一类权值的中位值对应的第二类像素值的中位值之前,所述方法包括:
根据以下公式计算邻域像素的第一类权值的直方图:
其中,δ(·)为狄拉克delta函数,h(p,n)是p像素的邻域像素的第一类权值直方图,w(p,q)代表p像素邻域内q像素的第一类权值,D(q)代表q像素的第二类像素值,N(p)代表p像素的邻域像素;
根据以下公式计算邻域像素的第一类权值累计直方图:
其中,c(p,m)表示p像素的邻域像素的第一类权值累计直方图。
5.根据权利要求4所述的方法,其特征在于,所述查找所述第一类权值的中位值对应的第二类像素值的中位值,包括:根据以下公式查找所述第一类权值的中位值对应的第二类像素值的中位值:
6.根据权利要求5所述的方法,其特征在于,所述将第二图像中邻域像素的第二类像素值中与第二类像素值的中位值接近的部分像素的第二像素值加权求和,包括:
根据以下公式将邻域像素中第二类像素值与第二类像素值的中位值接近的部分像素的第二类像素值加权求和:
其中,N(p)代表p像素的邻域像素,u(·)是单位阶跃函数,v(p,·)表示邻域各像素的第二类权值,T为第二类像素门限值,K(p)为归一化因子,其计算公式为:
7.根据权利要求6所述的方法,其特征在于,所述第二类权值的计算公式为v(p,q)=h(|d(p)-D(q)|),其中,h随|d(p)-D(q)|单调递减。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于北京华捷艾米科技有限公司,未经北京华捷艾米科技有限公司许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/201810336653.5/1.html,转载请声明来源钻瓜专利网。