首页 > 专利 > 杭州电子科技大学 > 一种基于噪声水平函数的信号依赖噪声分段线性去噪方法专利详情

一种基于噪声水平函数的信号依赖噪声分段线性去噪方法   0    0

有效专利 查看PDF
专利申请流程有哪些步骤?
专利申请流程图
申请
申请号:指国家知识产权局受理一件专利申请时给予该专利申请的一个标示号码。唯一性原则。
申请日:提出专利申请之日。
2019-01-28
申请公布
申请公布指发明专利申请经初步审查合格后,自申请日(或优先权日)起18个月期满时的公布或根据申请人的请求提前进行的公布。
申请公布号:专利申请过程中,在尚未取得专利授权之前,国家专利局《专利公报》公开专利时的编号。
申请公布日:申请公开的日期,即在专利公报上予以公开的日期。
2019-06-11
授权
授权指对发明专利申请经实质审查没有发现驳回理由,授予发明专利权;或对实用新型或外观设计专利申请经初步审查没有发现驳回理由,授予实用新型专利权或外观设计专利权。
2020-08-04
预估到期
发明专利权的期限为二十年,实用新型专利权期限为十年,外观设计专利权期限为十五年,均自申请日起计算。专利届满后法律终止保护。
2039-01-28
基本信息
有效性 有效专利 专利类型 发明专利
申请号 CN201910080360.X 申请日 2019-01-28
公开/公告号 CN109767405B 公开/公告日 2020-08-04
授权日 2020-08-04 预估到期日 2039-01-28
申请年 2019年 公开/公告年 2020年
缴费截止日
分类号 G06T5/00 主分类号 G06T5/00
是否联合申请 独立申请 文献类型号 B
独权数量 1 从权数量 2
权利要求数量 3 非专利引证数量 0
引用专利数量 3 被引证专利数量 0
非专利引证
引用专利 CN105911153B、CN103578477B、CN102436646A 被引证专利
专利权维持 3 专利申请国编码 CN
专利事件 许可 事务标签 公开、实质审查、授权、实施许可
申请人信息
申请人 第一申请人
专利权人 杭州电子科技大学 当前专利权人 杭州电子科技大学
发明人 宋佳忠、张钰、李锦彧 第一发明人 宋佳忠
地址 浙江省杭州市下沙高教园区2号大街 邮编 310018
申请人数量 1 发明人数量 3
申请人所在省 浙江省 申请人所在市 浙江省杭州市
代理人信息
代理机构
专利代理机构是经省专利管理局审核,国家知识产权局批准设立,可以接受委托人的委托,在委托权限范围内以委托人的名义办理专利申请或其他专利事务的服务机构。
杭州君度专利代理事务所 代理人
专利代理师是代理他人进行专利申请和办理其他专利事务,取得一定资格的人。
黄前泽
摘要
本发明公开了一种基于噪声水平函数的信号依赖噪声分段线性去噪方法。现有技术中,对信号依赖噪声去噪的图像融合方法需要通过相机进行多次捕获,在实际应用中消耗了大量人力和计算,十分消耗资源和时间,因此实际用途并不能采取。本发明如下:一、求取被处理图像的噪声水平函数曲线并分段;二、划分最终匹配块。三、在最终匹配块中进一步划分搜索框。四、更新搜索框中心位置的像素。本发明能够精确地对信号依赖噪声被处理图像进行噪声估计,并清楚的展现被处理图像的噪声水平函数。此外,本发明解决了传统的对信号依赖图像去噪方法不能参照实际的随灰度水平变换的噪声进行去噪的问题。
  • 摘要附图
    一种基于噪声水平函数的信号依赖噪声分段线性去噪方法
法律状态
序号 法律状态公告日 法律状态 法律状态信息
1 2023-03-14 专利实施许可合同备案的生效 IPC(主分类): G06T 5/00 合同备案号: X2023980032721 专利申请号: 201910080360.X 申请日: 2019.01.28 让与人: 杭州电子科技大学 受让人: 杭州晟视科技有限公司 发明名称: 一种基于噪声水平函数的信号依赖噪声分段线性去噪方法 申请公布日: 2019.05.17 授权公告日: 2020.08.04 许可种类: 普通许可 备案日期: 2023.02.24
2 2020-08-04 授权
3 2019-06-11 实质审查的生效 IPC(主分类): G06T 5/00 专利申请号: 201910080360.X 申请日: 2019.01.28
4 2019-05-17 公开
权利要求
权利要求书是申请文件最核心的部分,是申请人向国家申请保护他的发明创造及划定保护范围的文件。
1.一种基于噪声水平函数的信号依赖噪声分段线性去噪方法,其特征在于:步骤一、以a1为步长遍历被处理图像,选取出M个初始匹配块,初始匹配块的长宽均为s1;
步骤二、分别计算M个初始匹配块各自的平均值avgk如式(1)所示,k=1,2,…,M;
式(1)中,eijk为第k个初始匹配块上第i行像素与第j列像素的交点的像素值;
步骤三、将步骤二求得的M个初始匹配块各自的平均值avgk进行归一化,得到M个归一化均值avg′k如式(2)所示;
式(2)中,avgmax为M个初始匹配块各自的平均值avgk中的最大值;avgmin为M个初始匹配块各自的平均值avgk中的最小值;
步骤四、计算所有匹配块gk的总协方差矩阵Cg如式(3)所示;
式(3)中,gk为第k个初始匹配块所对应的像素值矩阵; 为矩阵gk的转置矩阵;
步骤五、计算总协方差矩阵Cg的各个特征值和特征向量,选取总协方差矩阵Cg的各个特征值中最小的那个特征值对应的特征向量作为最小方差向量
分别计算M个初始匹配块的局部噪声标准差 如式(4)所示,k=1,2,…,M;
式(4)中, 为矩阵 的欧几里德范数;
建立横坐标为局部均值,纵坐标为局部标准差的平面直角坐标系;将M个初始匹配块对应的坐标点加入平面直角坐标系中;将平面直角坐标系横坐标0~1的区间均分为b1个区域;在b1个区域中各取纵坐标最小的那个点作为最佳点,得到b1个最佳点;通过最小二乘法将b1个最佳点拟合成噪声水平函数曲线;
步骤六、将噪声水平函数曲线上横坐标0~1的区间沿横坐标方向均分为b2段局部曲线;
并分别对b2段局部曲线取平均值,得到b2个局部平均值δi′,i=1,2,…,b2;
步骤七、以a2为步长遍历被处理图像,选取出N个最终匹配块;最终匹配块的长宽均为s2;
步骤八、分别计算N个最终匹配块各自的平均值I_avgk如式(5)所示,k=1,2,...,N;
式(5)中,e′ijk为第k个最终匹配块上第i行像素与第j列像素的交点的像素值;
步骤九、将步骤八求得的N个最终匹配块各自的平均值I_avgk进行归一化,得到N个归一化均值I_avgk′如式(6)所示;
式(6)中,I_avgmax为N个最终匹配块各自的平均值I_avgk中的最大值;I_avgmin为N个最终匹配块各自的平均值I_avgk中的最小值;
步骤十、k=1,2,…,N,依次执行步骤十一;
步骤十一、更新像素值,完成去噪;
11-
1.将1赋值给i、j;
11-
2.以第k个最终匹配块的第i行、第j列的像素点作为中心点,选取长宽均为If的图像作为第i·j个固定搜索框;选取与第i·j个固定搜索框同心,且长宽均为It的的图像作为第i·j个中心框;
11-
3.将1赋值给p、q;
11-
4.以第i·j个固定搜索框的第p行、第q列的像素点作为中心点,长宽均为It的图像作为第i·j个固定搜索框内第p·q个相似框;
11-
5.计算第i·j个固定搜索框内第p·q个相似框的权重值wijkpq如式(7)所示;
式(7)中,D为第i·j个固定搜索框内第p·q个相似框对应的矩阵与第i·j个中心框对应的矩阵求差所得的矩阵内所有元素的平均值;hk为第k个最终匹配块的去噪滤波系数值,其值等于
进入步骤11-6;
11-
6.若p11-
7.将被处理图像中第k个最终匹配块的第i·j个固定搜索框的中心点像素值γ′ijk更新为 的表达式如式(8)所示;
式(8)中,wijk,max为第k个最终匹配块内第i·j个固定搜索框所有权重值中的最大值;
γijkpq为第k个最终匹配块内第i·j个固定搜索框的第p·q个相似框的中心点像素值;
11-
8.若i
2.根据权利要求1所述的一种基于噪声水平函数的信号依赖噪声分段线性去噪方法,其特征在于:所述的噪声水平函数曲线为一条一元二次函数曲线。

3.根据权利要求1所述的一种基于噪声水平函数的信号依赖噪声分段线性去噪方法,其特征在于:s1及s2均为奇数。
说明书

技术领域

[0001] 本发明属于图像去噪领域,具体涉及一种基于噪声水平函数的信号依赖噪声分段线性去噪方法。

背景技术

[0002] 图像已成为人类科技生活中息息相关的重要一部分,然而图像被实际相机采集后,一定会有噪声降低图像质量。来自实际相机的噪声被更好的建模成信号依赖噪声。大部分针对去噪方法都只是单纯地对信号无关噪声或称为高斯噪声在做去噪,其中对信号依赖噪声的去噪十分少。到目前为止,传统的对信号依赖噪声方法仍存在着一些缺点和不足。对信号依赖噪声去噪的图像融合方法需要通过相机进行多次捕获,在实际应用中消耗了大量人力和计算,十分消耗资源和时间,因此实际用途并不能采取。还有通过估计图像信号依赖噪声的噪声方差,来计算出全局的统一的噪声方差,根据该噪声方差进行高斯方法滤波,然而由于图像噪声是信号依赖噪声,是随图像灰度水平变化的,单独使用全局的噪声方差并不能更好得对不同灰度水平的图像噪声去噪,会导致去噪效果不佳,峰值信噪比低。因此提出一个随灰度水平变化的自适应局部的对信号依赖噪声去噪方法具有重要意义。

发明内容

[0003] 本发明所要解决的技术问题是提供一种基于噪声水平函数的信号依赖噪声分段线性去噪方法,以解决传统方法中图像信号依赖噪声对图像造成的问题,实现更好的去噪可视化效果。
[0004] 本发明的具体步骤如下:
[0005] 步骤一、以a1为步长遍历被处理图像,选取出M个初始匹配块,初始匹配块的长宽均为s1。
[0006] 步骤二、分别计算M个初始匹配块各自的平均值avgk如式(1)所示,k=1,2,…,M。
[0007]
[0008] 式(1)中,eijk为第k个初始匹配块上第i行像素与第j列像素的交点的像素值。
[0009] 步骤三、将步骤二求得的M个初始匹配块各自的平均值avgk进行归一化,得到M个归一化均值avgk′如式(2)所示。
[0010]
[0011] 式(2)中,avgmax为M个初始匹配块各自的平均值avgk中的最大值;avgmin为M个初始匹配块各自的平均值avgk中的最小值。
[0012] 步骤四、计算所有匹配块gk的总协方差矩阵Cg如式(3)所示。
[0013]
[0014] 式(3)中,gk为第k个初始匹配块所对应的像素值矩阵; 为矩阵gk的转置矩阵。
[0015] 步骤五、计算总协方差矩阵Cg的各个特征值和特征向量,选取总协方差矩阵Cg的各个特征值中最小的那个特征值对应的特征向量作为最小方差向量
[0016] 分别计算M个初始匹配块的局部噪声标准差 如式(4)所示,k=1,2,…,M。
[0017]
[0018] 式(4)中, 为矩阵 的欧几里德范数。
[0019] 建立横坐标为局部均值,纵坐标为局部标准差的平面直角坐标系。将M个初始匹配块对应的坐标点加入平面直角坐标系中。将平面直角坐标系横坐标0~1的区间均分为b1个区域。在b1个区域中各取纵坐标最小的那个点作为最佳点,得到b1个最佳点。通过最小二乘法将b1个最佳点拟合成噪声水平函数曲线。
[0020] 步骤六、将噪声水平函数曲线上横坐标0~1的区间沿横坐标方向均分为b2段局部曲线;并分别对b2段局部曲线取平均值,得到b2个局部平均值δ′i,i=1,2,…,b2。
[0021] 步骤七、以a2为步长遍历被处理图像,选取出N个最终匹配块。最终匹配块的长宽均为s2。
[0022] 步骤八、分别计算N个最终匹配块各自的平均值I_avgk如式(5)所示,k=1,2,…,N。
[0023]
[0024] 式(5)中,e′ijk为第k个最终匹配块上第i行像素与第j列像素的交点的像素值。
[0025] 步骤九、将步骤八求得的N个最终匹配块各自的平均值I_avgk进行归一化,得到N个归一化均值I_avgk′如式(6)所示。
[0026]
[0027] 式(6)中,I_avgmax为N个最终匹配块各自的平均值I_avgk中的最大值;I_avgmin为N个最终匹配块各自的平均值I_avgk中的最小值。
[0028] 步骤十、k=1,2,…,N,依次执行步骤十一。
[0029] 步骤十一、更新像素值,完成去噪。
[0030] 11-1.将1赋值给i、j。
[0031] 11-2.以第k个最终匹配块的第i行、第j列的像素点作为中心点,选取长宽均为If的图像作为第i·j个固定搜索框;选取与第i·j个固定搜索框同心,且长宽均为It的的图像作为第i·j个中心框。
[0032] 11-3.将1赋值给p、q。
[0033] 11-4.以第i·j个固定搜索框的第p行、第q列的像素点作为中心点,长宽均为It的图像作为第i·j个固定搜索框内第p·q个相似框。
[0034] 11-5.计算第i·j个固定搜索框内第p·q个相似框的权重值wijkpq如式(7)所示。
[0035]
[0036] 式(7)中,D为第i·j个固定搜索框内第p·q个相似框对应的矩阵与第i·j个中心框对应的矩阵求差所得的矩阵内所有元素的平均值。hk为第k个最终匹配块的去噪滤波系数值,其值等于
[0037] 进入步骤11-6。
[0038] 11-6.若p
[0039] 11-7.将被处理图像中第k个最终匹配块的第i·j个固定搜索框的中心点像素值γ′ijk更新为 的表达式如式(8)所示。
[0040]
[0041] 式(8)中,wijk,max为第k个最终匹配块内第i·j个固定搜索框所有权重值中的最大值;γijkpq为第k个最终匹配块内第i·j个固定搜索框的第p·q个相似框的中心点像素值。
[0042] 11-8.若i
[0043] 进一步地,所述的噪声水平函数曲线为一条一元二次函数曲线。
[0044] 进一步地,s1及s2均为奇数。
[0045] 本发明具有的有益效果是:
[0046] 本发明能够精确地对信号依赖噪声被处理图像进行噪声估计,并清楚的展现被处理图像的噪声水平函数。此外,本发明解决了传统的对信号依赖图像去噪方法不能参照实际的随灰度水平变换的噪声进行滤波的问题,本发明参照噪声水平函数,自适应的对被处理图像不同匹配块进行不同的去噪系数调整进行去噪,从而实现了更好的图像效果,极大地提高对信号依赖噪声的滤波去噪。

实施方案

[0047] 以下对本发明作进一步说明。
[0048] 一种基于噪声水平函数的信号依赖噪声分段线性去噪方法的具体步骤如下:
[0049] 步骤一、输入含有信号依赖噪声的被处理图像。被处理图像的分辨率为l·z;以a1为步长遍历被处理图像,a1=3。选取出M个初始匹配块,初始匹配块的长宽均为s1,s1为奇数(本实施例中取值为7), 第k个初始匹配块的中心点为被处理图像的第 行第 列相交处的那个像素点。k
=1,2,…,M。 表示k除以 所得的余数; 表示k除以
后向下取整所得值。
[0050] 步骤二、分别计算M个初始匹配块各自的平均值avgk如式(1)所示,k=1,2,…,M。
[0051]
[0052] 式(1)中,eijk为第k个初始匹配块上第i行像素与第j列像素的交点的像素值。
[0053] 步骤三、将步骤二求得的M个初始匹配块各自的平均值avgk进行归一化,得到M个归一化均值avg′k如式(2)所示。
[0054]
[0055] 式(2)中,avgmax为M个初始匹配块各自的平均值avgk中的最大值;avgmin为M个初始匹配块各自的平均值avgk中的最小值。
[0056] 步骤四、利用主成分分析法来估计局部噪声方差,计算所有匹配块gk的总协方差矩阵Cg如式(3)所示。
[0057]
[0058] 式(3)中,gk为第k个初始匹配块所对应的像素值矩阵(该矩阵的任意一个元素的大小均为第k个初始匹配块对应位置的像素值); 为矩阵gk的转置矩阵。
[0059] 步骤五、计算总协方差矩阵Cg的各个特征值和特征向量;由于在图像噪声处理中认为噪声的方差值较小,故选取总协方差矩阵Cg的各个特征值中最小的那个特征值对应的特征向量作为最小方差向量
[0060] 分别计算M个初始匹配块的局部噪声标准差 如式(4)所示,k=1,2,…,M。
[0061]
[0062] 式(4)中, 为矩阵 的欧几里德范数(即 最大的特征值开根号)。
[0063] 由于求得的局部噪声标准差包括了噪声标准差和一定的图像本身纹理标准差,即δn为我们需求的噪声标准差,δI为图像纹理标准差。
[0064] 建立横坐标为局部均值,纵坐标为局部标准差的平面直角坐标系。将M个初始匹配块对应的坐标点加入平面直角坐标系中,得到含M个离散点的离散坐标点图。将平面直角坐标系横坐标0~1的区间均分为b1个区域(即为第i个区域为横坐标[(i-1)b,ib))。在b1个区域中各取纵坐标最小的那个点作为最佳点,得到b1个最佳点。通过最小二乘法将b1个最佳点拟合成噪声水平函数曲线。噪声水平函数曲线为一条一元二次函数曲线(即f=c1·x2+c2·x+c3)。
[0065] 步骤六、对步骤六所得的噪声水平函数线进行分段线性处理。具体为:将噪声水平函数曲线上横坐标0~1的区间沿横坐标方向均分为b2段局部曲线(每段局部曲线对应横坐标的长度为1/b2);并分别对b2段局部曲线取平均值,得到b2个局部平均值δi′,i=1,2,…,b2。
[0066] 步骤七、以a2为步长遍历完被处理图像,a2=5,选取出N个最终匹配块。最终匹配块的长宽均为s2,s2为奇数(本实施例中取值为31), 第k个最终匹匹配块的中心点为被处理图像的第 行第
列相交处的那个像素点。k=1,2,…,N。
[0067] 步骤八、分别计算N个最终匹配块各自的平均值I_avgk如式(5)所示,k=1,2,…,N。
[0068]
[0069] 式(5)中,ei′jk为第k个最终匹配块上第i行像素与第j列像素的交点的像素值。
[0070] 步骤九、将步骤八求得的N个最终匹配块各自的平均值I_avgk进行归一化,得到N个归一化均值I_avgk′如式(6)所示。
[0071]
[0072] 式(6)中,I_avgmax为N个最终匹配块各自的平均值I_avgk中的最大值;I_avgmin为N个最终匹配块各自的平均值I_avgk中的最小值。
[0073] 步骤十、k=1,2,…,N,依次执行步骤十一。
[0074] 步骤十一、更新像素值,完成去噪。
[0075] 11-1.将1赋值给i、j。
[0076] 11-2.以第k个最终匹配块的第i行、第j列的像素点作为中心点,选取长宽均为If的图像作为第i·j个固定搜索框,If=21;选取与第i·j个固定搜索框同心,且长宽均为It的的图像作为第i·j个中心框,It=7。
[0077] 11-3.将1赋值给p、q。
[0078] 11-4.以第i·j个固定搜索框的第p行、第q列的像素点作为中心点,长宽均为It的图像作为第i·j个固定搜索框内第p·q个相似框。
[0079] 11-5.计算第i·j个固定搜索框内第p·q个相似框的权重值wijkpq如式(7)所示。
[0080]
[0081] 式(7)中,D为第i·j个固定搜索框内第p·q个相似框对应的矩阵与第i·j个中心框对应的矩阵求差所得的矩阵内所有元素的平均值。hk为第k个最终匹配块的去噪滤波系数值,其值等于 (即步骤六求得的第 段局部曲线的平均值)。 为e的 次方。
[0082] 进入步骤11-6。
[0083] 11-6.若p
[0084] 11-7.将被处理图像中第k个最终匹配块的第i·j个固定搜索框的中心点像素值γi′jk更新为 的表达式如式(8)所示。
[0085]
[0086] 式(8)中,wijk,max为第k个最终匹配块内第i·j个固定搜索框所有权重值中的最大值;γijkpq为第k个最终匹配块内第i·j个固定搜索框的第p·q个相似框的中心点像素值。
[0087] 11-8.若i
版权所有:盲专网 ©2023 zlpt.xyz  蜀ICP备2023003576号