图像恢复中的滤波方法技术

本发明专利技术公开了一种图像恢复中的滤波方法，包括以下步骤：1）对待处理图像进行分块预处理；2）建立图像在格林空间的状态空间模型；3）生成图像的固定滞后平滑器并对图像进行滤波处理。本发明专利技术将格林空间的滤波计算方法应用于图像处理中，在计算量上相比增广卡尔曼滤波算法大大降低，在未知噪声统计特性的情况下实现图像恢复，有效的提高了图像噪声滤波处理的适用性以及图像恢复的鲁棒性。

【技术实现步骤摘要】
图像恢复中的滤波方法
本专利技术属于计算机图像处理领域，具体涉及一种图像恢复中的滤波方法。其中，图像恢复是指通过计算机处理，对质量下降的图像加以重建或恢复的处理过程。
技术介绍
图像滤波是图像恢复中的一个重要步骤。现有技术中比较常用的是Hilbert空间的卡尔曼（Kalman）滤波。众所周知，卡尔曼滤波是一个不断地预测、修正的递推滤波过程，由于其在求解时不需要存储大量的观测数据，并且当得到新的观测数据时，可随时算得新的参数滤波值，便于实时的处理观测结果，因此被越来越多地应用于动态定位数据处理中。然而，在Hilbert空间使用Kalman滤波对系统进行滤波处理，我们需要事先知道系统的准确模型及噪声的统计特性，这在实际生产生活中都是很困难的。同时，卡尔曼滤波方法在初始迭代过程中，滤波器所得结果是与初值的选取密切相关的，一般认为该滤波方法在滤波的初始阶段误差较大、效果也差；只有当时间充分大时，滤波器的效果才与初值的选取无关，而滤波器趋于稳定状态。
技术实现思路
本专利技术所要解决的技术问题是：提出一种图像恢复中的滤波方法，解决卡尔曼滤波要求事先已知噪声的统计特性的问题，用来解决图像恢复问题。本专利技术的技术问题通过以下的技术方案予以解决：一种图像恢复中的滤波方法，包括以下步骤：1）对待处理图像进行分块预处理；2）建立图像在格林空间的状态空间模型；3）生成图像的固定滞后平滑器并对图像滤波。本专利技术与现有技术对比的有益效果是：本专利技术将格林空间的H∞滤波（即固定滞后平滑器）理论应用到图像恢复中，研究证实，H∞滤波在一定条件下就是格林空间卡尔曼滤波,因此可将Hilbert空间的滤波问题转化为格林空间的卡尔曼滤波问题，用来解决图像恢复问题。这种格林空间的滤波理论可以在没有事先知道噪声统计特性的情况下，对图像进行滤波并实现图像恢复。另外，采用H∞滤波器的滤波方法使得滤波误差始终维持在一个固定范围内，能够保证在整个图像范围内滤波效果都是较好的，包括初始阶段，有效的提高了滤波器对图像噪声的广泛适用性以及图像恢复的鲁棒性。【附图说明】图1是本专利技术实施例的技术方案流程框图。图2是本专利技术实施例对待处理图像分块的原理图。【具体实施方式】下面结合具体实施方式并对照附图对本专利技术做进一步详细说明。如图1所示，为本具体实施方式中的图像恢复中基于格林空间的固定滞后平滑器的滤波方法。利用格林空间中滤波器进行图像恢复的模块流程图，包括以下步骤：U1）：对待处理图像进行分块预处理；此步骤中，包含对待处理图像进行分块，以及确定分块后图像的每个块所代表的物理信息及状态向量的关系。具体地，对图像的分块预处理过程（如图2所示）可以表述为：每4个相邻的像素点按照如下规则整理成大小是4×1的列向量（X1,1(i),BX3,3(i)），然后选取这样的9个相邻向量组合得到一个36×1的状态向量X(i)。（上述数字4和9可以改为其它数，是数字n的平方即可，n≥2,但n最好不要太大超过7,如果太大会影响滤波效果。）X1,1(i)=[x1,1(i)x2,1(i)x1,2(i)x2,2(i)]TX2,1(i)=[x3,1(i)x4,1(i)x3,2(i)x4,2(i)]TX3,1(i)=[x5,1(i)x6,1(i)x5,2(i)x6,2(i)]TX1,2(i)=[x1,3(i)x2,3(i)x1,4(i)x2,4(i)]TX2,2(i)=[x3,3(i)x4,3(i)x3,4(i)x4,4(i)]TCX3,3(i)=[x5,5(i)x6,5(i)x5,6(i)x6,6(i)]T其中，列向量Xm,n(i)(1≤m,n≤3)是由第m,n个小方块内所包含的4个像素点的灰度值所组成的，而列向量X(i)是由9个小方块对应的列向量Xm,n(i)组成的。至此，依照该方法实现了对图像的分块，将大型的图像数据矩阵分成若干个小型矩阵，使得我们对图像的处理更加容易实现，也可在提高运算速度的同时减少计算量。U2）：建立图像在格林空间的状态空间模型；此步骤中，要建立图像在格林空间的状态空间模型。首先应依据U1）中对图像的分块结果以及每个向量所含的物理意义，列写系统模型的基本方程；其次再确定模型中各系数矩阵的具体表达式。根据格林空间系统模型的一般表达式，我们令图像具有如下形式的观测方程和量测方程：X(i+1)=FX(i)+Gu(i)(1)Y(i)=HX(i)+v(i)(2)其中，F,G分别具有如下形式：一般在实际图像中，通过图像的拍摄、传输等过程都存在对图像模糊化的现象。因此，矩阵H是参考点扩散函数（PSF）的特点来进行选取的。通常，矩阵H是由点扩散系数h组成的，而h满足下式：其中，r是模糊半径。在式(1)(2)所表示的系统中，X(i)表示每个块内的所有像素点灰度值，X(0)为初始状态，控制量u用来处理图像中所含的噪声、模糊等干扰因素，可选取它具有如下形式：u=[u1u2u3]T其中的v(i)表示能量有界的外部扰动，其统计特性未知。Y表示已知的观测量。根据以上的系统模型公式(1)(2)，我们可以设计相应的基于格林空间的平滑器来实现对图像的滤波并达到图像恢复的目的。此处，我们将图像的状态空间模型建立在格林空间内。其突出优点是，在格林空间内，允许噪声v(i)的内积为负定的（即v(i)vT(i)e0），这在普通的Hilberlt空间内是不被允许的，这样更便于理论上对滤波器的设计。另外，我们假设控制量u(i)和外部扰动v(i)与零均值的白噪声具有相同的统计特性，即有同时，关于初始状态X(0)的假设如下：初始状态X(0)独立于u(i)和v(i)，EX(0)=μ0，并且E[(X(0)-μ0)(X(0)-μ0)T]=P0其中，符号E为数学期望，T为转置号，δi,j=1(i=j),δi,j=0(i≠j)。此步骤中，还应确定滤波器（即基于格林空间的平滑器）存在的条件(因为我们设计的滤波器是在一定的性能指标下存在的，所以要确定这种滤波器存在需要满足的条件)。首先，初始条件X(0)记作X0，控制量u满足关系：u=u0+uu，这里u0是已知量，uu是未知量。其次，准则函数的应具有的形式：X0已知时满足：X0未知时满足：其中，表示控制量u的估计值，w表示噪声，γ是滤波器鲁棒性能指标。滤波器存在的条件是满足如下的Riccati方程：Pj=FPjFT+GGT-FPjHT(I+HPjHT)-1HPjFT且，P0=Π0。满足以上条件的图像才存在格林空间滤波器。U3）：设计图像的固定滞后平滑器对图像进行滤波处理，实现图像恢复。此步骤中，我们将对U2）中确定的系统方程进行格林空间的滤波器设计。涉及到的方法是基于经典Kalman滤波器和Mendel的输入白噪声估值器，设计相应的有限域固定滞后平滑器，并最终利用设计的平滑器对图像滤波以实现图像恢复。有限域固定滞后平滑器的问题可描述为：基于观测(y(t+N),y(t+N-1),B,y(1)),N>0,求状态x(t)的固定滞后平滑器，N固定。采用有限域固定滞后平滑器，根据后面块的图像信息估计前面块的图像信息。此处，我们可以选取系统（1）（2）的控制量u具有白噪声的形式。因为，系统是完全可观、可控的，且F是稳定矩阵，所以有相应的关于控制量u的稳态输入估值器，其形式如下：带初值且稳本文档来自技高网...

【技术保护点】
一种图像恢复中的滤波方法，其特征在于：包括以下步骤：1）对待处理图像进行分块预处理；2）建立图像在格林空间的状态空间模型；3）生成图像的固定滞后平滑器并对图像进行滤波处理，实现图像恢复。

【技术特征摘要】
1.一种图像恢复中的滤波方法，其特征在于：包括以下步骤：1)对待处理图像进行分块预处理；2)建立图像在格林空间的状态空间模型；3)生成图像的固定滞后平滑器并对图像进行滤波处理，实现图像恢复；所述步骤3)中图像的固定滞后平滑器具有如下形式的观测方程和量测方程：X(i+1)＝FX(i)+Gu(i)(1)Y(i)＝HX(i)+v(i)(2)其中，F,G分别具有如下形式：矩阵H是参考点扩散函数PSF的特点来进行选取的由点扩散系数h组成，而h满足下式：其中，r是模糊半径，x，y分别为选取的点的横坐标，纵坐标；X(i)表示每个块内的所有像素点灰度值，X(0)为初始状态，控制量u(i)用来处理...

【专利技术属性】
技术研发人员：王好谦，刘翠，张颖，
申请(专利权)人：清华大学深圳研究生院，
类型：发明
国别省市：