首页 / 专利库 / 视听技术与设备 / 振铃效应 / 一种基于边缘图块锐化的模糊图像非盲复原方法

一种基于边缘图锐化的模糊图像非盲复原方法

阅读:924发布:2020-05-08

专利汇可以提供一种基于边缘图锐化的模糊图像非盲复原方法专利检索,专利查询,专利分析的服务。并且本 发明 为一种基于边缘图 块 锐化的模糊图像非盲复原方法,包括以下步骤:建立边缘图块 数据库 ;提取模糊图像G的边缘分割图像G1和平滑区域分割图像G2;对边缘分割图像G1进行网格划分,在划分好的网格区域中针对有 像素 值的边缘区域 位置 图块进行匹配锐化处理,获得锐化更新后的边缘区域图像G1′;对平滑区域分割图像G2进行复原,获得复原的平滑区域图像G2′;获得模糊图像G的复原图像G′。本发明的优点是:通过对边缘区域图像进行网格划分,对于每个划分的单元网格,搜索确定图块数据库中与其灰度相似度匹配系数最大的清晰图块,完成边缘区域的的图块锐化更新,实现边缘清晰度提升的同时,避免了复原 算法 迭代 引起的 振铃效应 。,下面是一种基于边缘图锐化的模糊图像非盲复原方法专利的具体信息内容。

1.一种基于边缘图锐化的模糊图像非盲复原方法,其特征是:包括以下步骤:
S1:建立边缘图块数据库
S2:提取模糊图像G的边缘分割图像G1和平滑区域分割图像G2;
S3:对边缘分割图像G1进行网格划分,在划分好的网格区域中针对有像素值的边缘区域位置图块进行匹配锐化处理,获得锐化更新后的边缘区域图像G1′;
S4:对平滑区域分割图像G2进行复原,获得复原的平滑区域图像G2′;
S5:获得模糊图像G的复原图像G′。
2.根据权利要求1所述的一种基于边缘图块锐化的模糊图像非盲复原方法,其特征是:
所述步骤S1中建立边缘图块数据库包括以下步骤:
S11:建立基础边缘种子图块;
S12:对每个基础边缘种子图块在平水平方向、竖直方向及对线方向进行平移获得拓展图块,平移像素为+1像素和-1像素;
S13:对每个基础边缘种子图块及拓展图块进行旋转。
3.根据权利要求1所述的一种基于边缘图块锐化的模糊图像非盲复原方法,其特征是:
所述步骤S2中提取模糊图像G的边缘分割图像G1和平滑区域分割图像G2包括以下步骤:
S21:计算模糊图像G在局部窗口内的所有像素梯度绝对值
其中,M表示像素大小为n×n的局部窗口,窗口中心位于像素(x,y),hw和vw分别为水平方向和竖直方向的梯度算子,hw=[1 -1],vw=[1 -1]T,Ntotal为局部窗口M的像素总数:
Ntotal=n×n;
S22:根据像素梯度绝对值F(x,y)得到模糊图像G的边缘区域
其中,V为边缘提取阈值,EdgeMapl为边缘提取区域;
S23:通过闭操作得到优化的连通边缘区域:
EdgeMap=EdgeMapl·E
其中,E为形态学结构元素矩阵,符号·为闭运算;
S24:获得边缘分割图像G1和平滑区域分割图像G2
G1=G.*EdgeMap
G2=G.*~EdgeMap
其中,~EdgeMap表示对二值化边缘提取区域模板的非操作。
4.根据权利要求1所述的一种基于边缘图块锐化的模糊图像非盲复原方法,其特征是:
所述步骤S3中对边缘分割图像G1进行网格划分时,划分的单元网格像素大小与边缘图块数据库中图块像素大小一致。
5.根据权利要求1或4所述的一种基于边缘图块锐化的模糊图像非盲复原方法,其特征是:所述步骤S3中匹配锐化处理,包括以下步骤:
S31:对单元网格图块X进行亮度归一化处理
其中,XN为灰度值归一化的单元网格图块X,Xmax和Xmin分别为网格单元图块X像素值的最大值和最小值,Dmax和Dmin数据图像像素值的最大值和最小值,对于8位的图像块Dmax=255,Dmin=0;
S32:选取边缘图块数据库中与图块XN的灰度相似度匹配系数最大的边缘图块,作为图块XN的替换图块XNsharp,
其中,ε(m,n)为灰度相似度匹配系数,Y(m,n)为边缘图块数据库对应位置下的边缘图块,M为图像块Y(m,n)和图块XN的各自像素总数,XNi表示图块XN第i个像素点的灰度值,Y(m,n)i表示边缘图块Y(m,n)第i个像素点的灰度值;
S33:对图块XNsharp灰度值去归一化,得到与单元网格图块X像素值均值和对比度分布一致的图块
其中,Xsharp为灰度值去归一化的XNsharp,XNsharpmax和XNsharpmin分别为XNsharp的最大值和最小值,Tmax和Tmin为单元网格图块X像素值的最大值和最小值;
S34:对边缘分割图像G1中各个有像素值的单元网格边缘区域位置图块完成匹配锐化更替,得到锐化更新后的边缘区域图像G1′。
6.根据权利要求1所述的一种基于边缘图块锐化的模糊图像非盲复原方法,其特征是:
所述S4中采用RL非盲复原方法对平滑区域分割图像G2进行复原
其中,L为RL复原的迭代次数,H*为H的伴随矩阵,H为模糊图像G的点扩散函数能量分布的模糊核。
7.根据权利要求1所述的一种基于边缘图块锐化的模糊图像非盲复原方法,其特征是:
所述步骤S5中复原图像G’通过以下公式获得
G’=G′1+G2′。

说明书全文

一种基于边缘图锐化的模糊图像非盲复原方法

技术领域

[0001] 本发明涉及图像处理领域,尤其涉及一种基于边缘图块锐化的模糊图像非盲复原方法。

背景技术

[0002] 图像模糊问题是数字成像采集设备经常会遇到的问题。由于成像光学系统存在离焦、在的低照度拍摄时的手持成像设备需要长曝光时间,或者拍摄过程中运动场景和成像传感器发生相对运动,都会导致获取的图像出现模糊退化,造成图像分辨率和探测能的下降。需要引入图像复原算法来改善图像质量
[0003] 传统的非盲复原方法,即使在点扩散函数事先知道的情况下,需要利用图像先验信息,对图像复原约束能量方程建模,开展图像复原求解。这种算法往往会造成图像噪声的放大,同时,在图像边缘区域也容易出现严重的振铃条纹,影响复原图像质量。而图像边缘信息,是人眼视觉评价最着重关注的区域,属于中高频信息,对边缘区域的清晰度提升和改善,对于图像主观视觉感受有着最直接的影响。

发明内容

[0004] 本发明主要解决了现有非盲复原方法会放大图像噪声、在图像边缘出现严重振铃条纹的问题,提供了一种能够有效提升图像边缘区域分辨率,同时避免了振铃等伪信息出现的基于边缘图块锐化的模糊图像非盲复原方法。
[0005] 本发明解决其技术问题所采用的技术方案是,一种基于边缘图块锐化的模糊图像非盲复原方法,包括以下步骤:
[0006] S1:建立边缘图块数据库
[0007] S2:提取模糊图像G的边缘分割图像G1和平滑区域分割图像G2;
[0008] S3:对边缘分割图像G1进行网格划分,在划分好的网格区域中针对有像素值的边缘区域位置图块进行匹配锐化处理,获得锐化更新后的边缘区域图像G1′;
[0009] S4:对平滑区域分割图像G2进行复原,获得复原的平滑区域图像G2′;
[0010] S5:获得模糊图像G的复原图像G′。
[0011] 利用局部窗口梯度绝对值和形态学闭操作,有效提取了图像中的主要边缘区域,将图像区分为边缘区域图像和平滑区域图像,为后续的处理提供了基础。针对图像边缘区域,通过建立边缘图块数据库采用图块匹配锐化更替的方法,来提升其清晰度和空间分辨率,通过对边缘区域图像进行网格划分,对于每个划分的单元网格,搜索确定图块数据库中与其灰度相似度匹配系数最大的清晰图块,完成边缘区域的图块锐化更新,实现边缘清晰度提升的同时,避免了复原算法迭代引起的振铃效应。针对平滑区域,结合已知的点扩散函数,引入RL非盲复原算法,开展清晰度复原提升。平滑区域边缘纹理信息较少,人眼视觉敏感程度较低,设置合理的迭代参数,可以很好地改善该区域清晰度,最终得到主观视觉感受较好的整幅复原结果。
[0012] 作为上述方案的一种优选方案,所述步骤S1中建立边缘图块数据库包括以下步骤:
[0013] S11:建立基础边缘种子图块;
[0014] S12:对每个基础边缘种子图块在平水平方向、竖直方向及对线方向进行平移获得拓展图块,平移像素为+1像素和-1像素;
[0015] S13:对每个基础边缘种子图块及拓展图块进行旋转。
[0016] 作为上述方案的一种优选方案,所述步骤S2中提取模糊图像G的边缘分割图像G1和平滑区域分割图像G2包括以下步骤:
[0017] S21:计算模糊图像G在局部窗口内的所有像素梯度绝对值
[0018]
[0019] 其中,M表示像素大小为n×n的局部窗口,窗口中心位于像素(x,y),hw和vw分别为水平方向和竖直方向的梯度算子,hw=[1 -1],vw=[1 -1]T,Ntotal为局部窗口M的像素总数:Ntotal=n×n;
[0020] S22:根据像素梯度绝对值F(x,y)得到模糊图像G的边缘区域
[0021]
[0022] 其中,V为边缘提取阈值,EdgeMapl为边缘提取区域;
[0023] S23:通过闭操作得到优化的连通边缘区域:
[0024] EdgeMap=EdgeMapl·E
[0025] 其中,E为形态学结构元素矩阵,符号·为闭运算;
[0026] S24:获得边缘分割图像G1和平滑区域分割图像G2
[0027] G1=G.*EdgeMap
[0028] G2=G.*~EdgeMap
[0029] 其中,~EdgeMap表示对二值化边缘提取区域模板的非操作。
[0030] 作为上述方案的一种优选方案,所述步骤S3中对边缘分割图像G1进行网格划分时,划分的单元网格像素大小与边缘图块数据库中图块像素大小一致。
[0031] 作为上述方案的一种优选方案,所述步骤S3中匹配锐化处理,包括以下步骤:
[0032] S31:对单元网格图块X进行亮度归一化处理
[0033]
[0034] 其中,XN为灰度值归一化的单元网格图块X,Xmax和Xmin分别为网格单元图块X像素值的最大值和最小值,Dmax和Dmin数据图像像素值的最大值和最小值,对于8位的图像块Dmax=255,Dmin=0;
[0035] S32:选取边缘图块数据库中与图块XN的灰度相似度匹配系数最大的边缘图块,作为图块XN的替换图块XNsharp,
[0036]
[0037] 其中,ε(m,n)为灰度相似度匹配系数,Y(m,n)为边缘图块数据库对应位置下的边缘图块,M为图像块Y(m,n)和图块XN的各自像素总数,XNi表示图块XN第i个像素点的灰度值,Y(m,n)i表示边缘图块Y(m,n)第i个像素点的灰度值;
[0038] S33:对图块XNsharp灰度值去归一化,得到与单元网格图块X像素值均值和对比度分布一致的图块
[0039]
[0040] 其中,Xsharp为灰度值去归一化的XNsharp,XNsharpmax和XNsharpmin分别为XNsharp的最大值和最小值,Tmax和Tmin为单元网格图块X像素值的最大值和最小值;
[0041] S34:对边缘分割图像G1中各个有像素值的单元网格边缘区域位置图块完成匹配锐化更替,得到锐化更新后的边缘区域图像G1′。
[0042] 作为上述方案的一种优选方案,所述S4中采用RL非盲复原方法对平滑区域分割图像G2进行复原
[0043]
[0044] 其中,L为RL复原的迭代次数,H*为H的伴随矩阵,H为模糊图像G的点扩散函数能量分布的模糊核。
[0045] 作为上述方案的一种优选方案,所述步骤S5中复原图像G’通过以下公式获得[0046] G’=G′1+G2′。
[0047] 本发明的优点是:通过对边缘区域图像进行网格划分,对于每个划分的单元网格,搜索确定图块数据库中与其灰度相似度匹配系数最大的清晰图块,完成边缘区域的的图块锐化更新,实现边缘清晰度提升的同时,避免了复原算法迭代引起的振铃效应。附图说明
[0048] 图1为本发明的一种流程示意图。
[0049] 图2为本发明中建立边缘图块数据库的一种流程示意图。
[0050] 图3为本发明中基础边缘种子图块的一种示意图。
[0051] 图4为本发明中提取边缘分割图像G1和平滑区域分割图像G2的一种流程示意图。
[0052] 图5为本发明中匹配锐化处理的一种流程示意图。

具体实施方式

[0053] 下面通过实施例,并结合附图,对本发明的技术方案作进一步的说明。
[0054] 实施例:
[0055] 本实施例一种基于边缘图块锐化的模糊图像非盲复原方法,如图1所示,包括以下步骤:
[0056] S1:建立边缘图块数据库,如图2所示,包括以下步骤:
[0057] S11:建立基础边缘种子图块,本实施例中共建立12种像素大小为5×5的基础边缘种子图块,如图3所示,从左到右从上到下分别为单像素阶跃边缘种子图像块:
[0058]
[0059] 两个单像素阶跃边缘种子图像块:
[0060]
[0061] 双像素阶跃边缘种子图像块:
[0062]
[0063] 拐角边缘种子图像块:
[0064]
[0065] 环形边缘种子图像块:
[0066]
[0067] 十字边缘种子图像块:
[0068]
[0069] 单像素阶跃边缘对比度翻转种子图像块:
[0070]
[0071] 两个单像素阶跃边缘对比度翻转种子图像块:
[0072]
[0073] 双像素阶跃边缘对比度翻转种子图像块:
[0074]
[0075] 拐角边缘对比度翻转种子图像块:
[0076]
[0077] 环形边缘对比度翻转种子图像块:
[0078]
[0079] 十字边缘对比度翻转种子图像块:
[0080]
[0081] S12:对每个基础边缘种子图块在水平水平方向、竖直方向及对角线方向进行平移获得8个拓展图块,平移像素为+1像素和-1像素;
[0082] S13:对每个基础边缘种子图块及拓展图块进行顺时针旋转,以每个图块的中心为原点旋转一周,及旋转角度步长为3°,即每个图块旋转120次,构建的边缘图块数据库中共包含(8+1)×12×120=12960个清晰的边缘图像子块。
[0083] S2:提取模糊图像G的边缘分割图像G1和平滑区域分割图像G2,如图4所示,包括以下步骤:
[0084] S21:计算模糊图像G在局部窗口内的所有像素梯度绝对值
[0085]
[0086] 其中,M表示像素大小为5×5的局部窗口,窗口中心位于像素(x,y),hw和vw分别为水平方向和竖直方向的梯度算子,hw=[1 -1],vw=[1 -1]T,Ntotal为局部窗口M的像素总数:Ntotal=5×5;
[0087] S22:根据像素梯度绝对值F(x,y)得到模糊图像G的边缘区域
[0088]
[0089] 其中,V为边缘提取阈值,V=[0.001,0.01],EdgeMapl为边缘提取区域;
[0090] S23:通过闭操作得到优化的连通边缘区域:
[0091] EdgeMap=EdgeMapl·E
[0092] 其中,E为形态学结构元素矩阵, 符号·为闭运算,即利用E对EdgeMapl先膨胀再腐蚀得到结果;
[0093] S24:获得边缘分割图像G1和平滑区域分割图像G2
[0094] G1=G.*EdgeMap
[0095] G2=G.*~EdgeMap
[0096] 其中,~EdgeMap表示对二值化边缘提取区域模板的非操作。;
[0097] S3:对边缘分割图像G1进行网格划分,划分的单元网格像素大小与边缘图块数据库中图块像素大小一致,即将边缘分割图像G1划分为单元网格为5×5像素大小的网格,在划分好的网格区域中针对有像素值的边缘区域位置图块进行匹配锐化处理,获得锐化更新后的边缘区域图像G1′,如图5所示,包括以下步骤:
[0098] S31:为消除真实图像中亮度和对比度差异的干扰,对单元网格图块X进行亮度归一化处理
[0099]
[0100] 其中,XN为灰度值归一化的单元网格图块X,Xmax和Xmin分别为网格单元图块X像素值的最大值和最小值,Dmax和Dmin数据图像像素值的最大值和最小值,对于8位的图像块Dmax=255,Dmin=0;
[0101] S32:选取边缘图块数据库中与图块XN的灰度相似度匹配系数最大的边缘图块,作为图块XN的替换图块XNsharp,
[0102]
[0103] 其中,ε(m,n)为灰度相似度匹配系数,Y(m,n)为边缘图块数据库对应位置下的边缘图块,M为图像块Y(m,n)和图块XN的各自像素总数,即M=5×5=25,XNi表示图块XN第i个像素点的灰度值,Y(m,n)i表示边缘图块Y(m,n)第i个像素点的灰度值;
[0104] S33:对图块XNsharp灰度值去归一化,得到与单元网格图块X像素值均值和对比度分布一致的图块
[0105]
[0106] 其中,Xsharp为灰度值去归一化的XNsharp,XNsharpmax和XNsharpmin分别为XNsharp的最大值和最小值,Tmax和Tmin为单元网格图块X像素值的最大值和最小值;
[0107] S34:对边缘分割图像G1中各个有像素值的单元网格边缘区域位置图块完成匹配锐化更替,得到锐化更新后的边缘区域图像G1′。
[0108] S4:采用RL非盲复原方法对平滑区域分割图像G2进行复原,获得复原的平滑区域图像G2′
[0109]
[0110] 其中,L为RL复原的迭代次数,本实施例中L设置为30~40,H*为H的伴随矩阵,H为模糊图像G的点扩散函数能量分布的模糊核;
[0111] S5:获得模糊图像G的复原图像G′,复原图像G′由G1′和G2′组合得到。
[0112] 本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈