基于自适应判决反馈均衡技术的盲图像复原方法
专利汇可以提供基于自适应判决反馈均衡技术的盲图像复原方法专利检索，专利查询，专利分析的服务。并且一种基于自适应判决反馈均衡技术的盲图像复原方法，其包括如下步骤：退化图像通过前向 滤波器 得到部分复原图像，同时将NL滤波器的 迭代 结果输入到反馈滤波器以得到退化图像的部分退化因素，然后对部分复原图像和部分退化因素求差产生图像的估计并进行投影，然后对图像的估计和NL滤波器的输出求差产生每个 像素 点的迭代误差，进而通过运算器和加法器得到盲图像复原的总误差，更新前向滤波器和反馈滤波器的参数并重复进行上述步骤的迭代循环，直到相邻两次迭代过程的盲图像复原总误差小于参考 阈值 后，终止迭代循环，输出图像。本 发明 的优点是： 算法 结构简单、收敛性好、计算复杂度低，在低 信噪比 情况下能够取得比较满意的退化图像复原效果。，下面是基于自适应判决反馈均衡技术的盲图像复原方法专利的具体信息内容。
1.一种基于自适应判决反馈均衡技术的盲图像复原方法，其特征在于包括如 下步骤：
第一步，区分退化图像的目标区域和背景区域；
第二步，退化图像通过前向滤波器卷积滤波得到部分复原图像，同时将NL滤 波器的迭代结果输入到反馈滤波器以得到退化图像的部分退化因素；
第三步，通过减法器对部分复原图像和部分退化因素求差产生图像的估计；
第四步，通过NL滤波器对图像的估计进行投影；
第五步，通过减法器对图像的估计和NL滤波器的输出求差产生每个像素点的 迭代误差；
第六步，通过运算器和加法器对每个像素点的迭代误差根据相应算法进行求和 以得到盲图像复原的总误差；
第七步，通过运算器更新前向滤波器和反馈滤波器的参数并重复进行步骤二至 步骤六以进行下一次的迭代循环，从而得到下一次的盲图像复原的总误差，比较判 断相邻两次迭代过程的盲图像复原总误差之差和参考阈值，如果相邻两次迭代过程 的盲图像复原总误差之差的绝对值小于参考阈值，则终止迭代循环，输出图像，否 则，重复进行步骤二至步骤七。
2.如权利要求1所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的步骤一为通过将退化图像的目标区域灰度值和背景区域灰度值分 别与设定灰度阈值进行比较，如果大于或等于设定灰度阈值，则为目标区域，如果 小于设定灰度阈值，则为背景区域。
3.如权利要求2所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的前向滤波器是一个N×N的矩阵：
[Uk(1，1)，...，Uk((Nx+1)/2，(Ny+1)/2)，...，Uk(Nx，Ny)]，
这里的N取大于3的奇数，其中，Uk(1，1)表示矩阵初始位置的滤波器参数第K 次迭代的结果，Uk((Nx+1)/2，(Ny+1)/2)表示矩阵中心位置的滤波器参数第K次迭代 的结果，Uk(Nx，Ny)表示矩阵末尾位置的滤波器参数第K次迭代的结果。
4.如权利要求3所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的反馈滤波器Wk(x，y)是一个N×N的矩阵：
[Wk(1，1)，...，Wk((Nx+1)/2，(Ny+1)/2)，...，Wk(Nx，Ny)]，N取大于3的奇数，
其中，Wk(1，1)表示矩阵初始位置的滤波器参数第K次迭代的结果，Wk((Nx+1)/2， (Ny+1)/2)表示矩阵中心位置的滤波器参数第K次迭代的结果，Wk(Nx，Ny)表示矩 阵末尾位置的滤波器参数第K次迭代的结果。
5.如权利要求4所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的NL滤波器是非线性滤波器。
6.如权利要求5所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的NL滤波器为：
其中Dsup为目标区域内部所有像素的集合，而 为目标区域外部所有像素 的集合，LB为模糊图像背景的灰度值。
7.如权利要求6所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的步骤四为对目标区域内小于零的像素灰度值投影为零，大于或等于 零的像素灰度值保持不变，在背景区域内的像素灰度值投影为LB。
8.如权利要求7所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的步骤六中的运算器采用最小均方算法。
9.如权利要求8所述的基于自适应判决反馈均衡技术的盲图像复原方法，其 特征在于所述的步骤七中的运算器采用共轭梯度法。
技术领域
本发明属于计算机图像处理中盲图像复原领域，尤其是涉及一种基于自适应判决 反馈均衡技术的盲图像复原方法。
背景技术
发明内容
自适应判决反馈均衡技术利用的是将判决后的信号作为后向抽头的输入信号，以 消除噪声对输出信号的影响。由于图像的退化可以看作是像素点之间的干扰，所以将 这种技术引入到二维退化图像复原中，可实现退化图像的复原。
本发明的基于自适应判决反馈均衡技术的盲图像复原方法主要包括如下步骤：
第一步，区分退化图像的目标区域和背景区域；
第二步，退化图像通过前向滤波器卷积滤波得到部分复原图像，同时将NL滤波 器的迭代结果输入到反馈滤波器以得到退化图像的部分退化因素；
第三步，通过减法器对部分复原图像和部分退化因素求差产生图像的估计；
第四步，通过NL滤波器对图像的估计进行投影；
第五步，通过减法器对图像的估计和NL滤波器的输出求差产生每个像素点的迭 代误差；
第六步，通过运算器和加法器对每个像素点的迭代误差根据相应算法进行求和以 得到盲图像复原的总误差；
第七步，通过运算器更新前向滤波器和反馈滤波器的参数并重复进行步骤二至步 骤六以进行下一次的迭代循环，从而得到下一次的盲图像复原的总误差，比较判断相 邻两次迭代过程的盲图像复原总误差之差和参考阈值，如果相邻两次迭代过程的盲图 像复原总误差之差的绝对值小于参考阈值，则终止迭代循环，输出图像，否则，重复 进行步骤二至步骤七。
本发明具有以下优点：(1)算法结构简单、收敛性好、计算复杂度低；(2)使用 一个前向滤波器和一个反馈滤波器模拟退化因子的倒数，避免了使用单个逆滤波器造 成的高频噪声放大问题，在低信噪比情况下能够取得比较满意的复原效果；(3)只需 要知道图像的目标区域信息而不需要退化函数的任何信息便可进行盲图像复原，因此 特别适合应用在退化因子未知时图像的复原，是一种很有发展前景的图像复原方法。
附图说明
下面参照附图和具体实施方式对本发明的基于自适应判决反馈均衡技术的盲图像 复原方法作进一步详细地说明。
图1是本发明的方法流程示意图；
图2是本发明的信噪比为15db时的图像复原效果示意图；
图3是本发明的信噪比为35db时的图像复原效果示意图。
具体实施方式
第一步，通过将退化图像的目标区域灰度值和背景区域灰度值分别与设定灰度阈 值进行比较，如果大于或等于设定灰度阈值，则为目标区域；如果小于设定灰度阈值， 则为背景区域，从而区分退化图像的目标区域和背景区域。
第二步，退化图像通过前向滤波器卷积滤波得到部分复原图像，所述的前向滤波 器是一个N×N的矩阵：[Uk(1，1)，...，Uk((Nx+1)/2，(Ny+1)/2)，...，Uk(Nx，Ny)]，这里的N 取大于3的奇数，其中，Uk(1，1)表示矩阵初始位置的滤波器参数第K次迭代的结果， Uk((Nx+1)/2，(Ny+1)/2)表示矩阵中心位置的滤波器参数第K次迭代的结果，Uk(Nx，Ny) 表示矩阵末尾位置的滤波器参数第K次迭代的结果。
这里，前向滤波器的中心位置参数在第1次迭代的初值设为1，矩阵中其余位置 参数都设为零，退化图像g(x，y)与Uk(x，y)卷积滤波后得到的部分复原图像pk(x，y)为：
pk(x，y)＝g(x，y)*Uk(x，y) g(x，y)为以(x，y)为中心位置的N×N像素矩阵。
将NL滤波器的迭代结果输入到反馈滤波器以得到退化图像的部分退化因素，(关 于NL滤波器的说明参照下面第四步的介绍。)
这里的反馈滤波器Wk(x，y)为一个N×N的矩阵：
[Wk(1，1)，...，Wk((Nx+1)/2，(Ny+1)/2)，...，Wk(Nx，Ny)]，N取大于3的奇数，
其中，Wk(1，1)表示矩阵初始位置的滤波器参数第K次迭代的结果，Wk((Nx+1)/2， (Ny+1)/2)表示矩阵中心位置的滤波器参数第K次迭代的结果，Wk(Nx，Ny)表示矩阵 末尾位置的滤波器参数第K次迭代的结果，这里，反馈滤波器初始值全设为零，部分 退化因素表示为：
qk(x，y)＝fNLk(x，y)*Wk(x，y)
这里的fNLk(x，y)为第k次迭代NL滤波器的输出结果。
第三步，通过减法器对部分复原图像和部分退化因素求差产生图像的估计；
第四步，通过NL滤波器对图像的估计进行投影，即将目标区域内小于零的像素 灰度值投影为零，大于或等于零的像素灰度值保持不变，在背景区域内的像素灰度值 投影为LB，所述的NL滤波器是一个非线性的滤波器，其定义如下：
其中Dsup为目标区域内部所有像素的集合，而 为目标区域外部所有像素的集 合，LB为背景区域的像素灰度值。
第五步，通过减法器对图像的估计和NL滤波器的输出求差产生每个像素点的迭 代误差，因此，上述投影过程产生的每个像素点第k次迭代的误差为：
ek(x，y)＝fNLk(x，y)-fk(x，y)，这里的fk(x，y)第k次迭代的图像估计。
第六步，通过运算器和加法器对每个像素点的迭代误差根据相应算法进行求和以 得到盲图像复原的总误差，在本实施例中，采用最小均方(LMS)算法，通过加法器 对每个像素点第k次迭代的误差进行求和得到的盲图像复原总误差表示为：
第七步，通过运算器更新前向滤波器和反馈滤波器的参数并重复进行步骤二至步 骤六以进行下一次的迭代循环，在本实施例中，采用共轭梯度法进行更新参数，然后 比较判断相邻两次迭代过程的盲图像复原总误差之差和参考阈值，如果相邻两次迭代 过程的盲图像复原总误差之差的绝对值小于参考阈值，则终止迭代循环，输出图像。
这里设相邻两次迭代过程的代价函数值分别为Jold、Jnew，令ΔJ＝|Jold-Jnew|， δ＝10-6为参考阈值，若ΔJ＜δ，则终止迭代循环，输出图像。
图2所示的是采用本发明方法对信噪比为15db的图像复原效果示意图，(a)为 128*256的二值图，背景灰度值为零；(b)为退化图像，信噪比为15db；(c)中白色 区域为图像的目标区域，这里阈值设定为60；(d)为使用本发明方法取得的复原图像。
图3所示的是采用本发明方法对信噪比为35db的图像复原效果示意图，(a)为 128*128的灰度图，背景灰度值为零；(b)为退化图像，信噪比为20db；(c)中白色 区域为图像的目标区域，这里阈值设定为12；(d)为使用本发明方法取得的复原图像。
以上所述的仅是本发明的优选实施方式。应当指出，对于本领域的普通技术人员 来说，在不脱离本发明原理的前提下，还可以作出若干变型和改进，这些变化也应视 为属于本发明的保护范围。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种图像篡改盲取证方法 | 2020-05-20 | 829 |
| 分布式无线传感器网络数据漂移盲校准方法 | 2020-05-15 | 838 |
| 基于深度姿态预估和多特征融合的行人衣着属性识别方法 | 2020-05-25 | 289 |
| 一种基于OFDM反向散射通信系统半盲接收机设计方法 | 2020-05-20 | 266 |
| 一种不需要参考发射波形的穿墙雷达墙体参数估计方法 | 2020-05-25 | 209 |
| 一种结合灰度熵和盲反卷积的采样率自适应分块压缩感知方法 | 2020-05-26 | 532 |
| 一种基于kernel-free图像复原的高速移动对象机器视觉尺寸检测方法 | 2020-05-12 | 563 |
| 一种利用预处理和后处理来提高光学显微系统分辨率的方法 | 2020-05-21 | 483 |
| 一种基于电流波数扫描的透镜三维轮廓测量装置及方法 | 2020-05-26 | 876 |
| 基于多流注意对抗网络的动态场景模糊图像盲复原方法 | 2020-05-11 | 986 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
