【摘要】 本发明提供一种运动模糊和散焦复合模糊的图像复原方法，该方法能够对运动模糊和散焦同时出现的图像进行参数估计和图像复原，它包含如下步骤：(1)建立高斯白噪声模板，将降质图像与白噪声模板卷积，达到去除噪声的目的；(2)通过图像能量谱估计图像的模糊主方向和副方向(3)计算图像的主方向导数矩阵和副方向导数矩阵；(4)分别对主方向导数矩阵和副方向导数矩阵实施自相关运算和方向累加运算；(5)根据主方向导数自相关的累加曲线和副方向导数自相关的累加曲线，估计主方向模糊长度和副方向模糊长度；(6)根据获得的主方向模糊长度和副方向模糊长度，建立复合模糊模型；(7)利用维纳滤波对降质图像进行复原。 【专利类型】发明申请 【申请人】东南大学 【申请人类型】学校 【申请人地址】210096 江苏省南京市四牌楼2号 【申请人地区】中国 【申请人城市】南京市 【申请人区县】玄武区 【申请号】CN201010522739.0 【申请日】2010-10-27 【申请年份】2010 【公开公告号】CN101968881A 【公开公告日】2011-02-09 【公开公告年份】2011 【授权公告号】CN101968881B 【授权公告日】2012-07-18 【授权公告年份】2012.0 【IPC分类号】G06T5/00 【发明人】路小波; 李楠 【主权项内容】1.一种运动模糊和散焦复合模糊的图像复原方法，其特征在于，包括如下步骤： 步骤1降质图像去噪 利用低通滤波对降质图像g(x，y)进行去噪处理，其中x和y分别为行坐标和列坐标，且都为大于0的整数，得到去噪后的降质图像f′(x，y)，其算法如下： 均值为零的高斯分布p(x，y)的函数形式pz(i1，i2)可以表示为： 其中下标z表示均值为零的高斯分布的函数形式，i1和i2分别为横坐标和纵坐标，且都为实数，σ为标准差，e为自然对数函数的底数，π为圆周率，设成像系统的加性噪声n(x，y)服从均值为零的高斯分布p(x，y)，这里σ取0.5，建立一个3×3的高斯白噪声模板n′(x，y)： 将降质图像g(x，y)与高斯白噪声模板n′(x，y)进行卷积，得到M1×M2大小的去噪后的降质图像f′(x，y)： f′(x，y)＝g(x，y)*n′(x，y) 其中*为卷积操作， 步骤2模糊方向识别 利用公知的傅立叶变换来估计去噪后的降质图像f′(x，y)的主方向θ和副方向θ′，假设此时主方向θ的角度为α，-90°≤α＜90°，副方向θ′的角度为β，0°≤β＜180°，旋转图像至主方向θ的角度为0°，具体方法如下：设散焦与运动模糊同时发生时运动模糊的方向为主方向θ，与主模糊方向垂直的方向为副方向θ′，计算步骤1得到的去噪后的降质图像f′(x，y)的傅立叶变换F(ω1，ω2)： ω1和ω2为频率变量，F(ω1，ω2)为复数，j为虚数单位，去噪后的降质图像的傅立叶变换F(ω1，ω2)可以表示为： F(ω1，ω2)＝R(ω1，ω2)+jI(ω1，ω2) 其中R(ω1，ω2)为实部，I(ω1，ω2)为虚部， 计算待估计降质图像f′(x，y)的能量谱E(ω1，ω2)： E(ω1，ω2)＝R(ω1，ω2)2+I(ω1，ω2)2 能量谱E(ω1，ω2)中各元素为实数，将能量谱E(ω1，ω2)显示出来，在图中心位置会出现一个类似于椭圆的白色亮斑，亮斑周围会有连续的亮环，在所有穿过亮斑的弦中，最长的弦lmax所处的方向，取0°到180°之间的角做为β的值，与此垂直的方向角为α，α＝β-90°，此时主方向θ的角度为α，副方向θ′的角度为β，以垂直于去噪后的降质图像f′(x，y)，且过去噪后的降质图像f′(x，y)中心的直线为旋转轴，如果0°≤α＜90°，以α角顺时针旋转去噪后的降质图像f′(x，y)，如果-90°≤α＜0°，以α+90°角顺时针旋转去噪后的降质图像f′(x，y)，降质图像f′(x，y)旋转后，主方向θ的角度为 0°，副方向θ′的角度为90°，得到M3×M4大小的旋转后的去噪后的降质图像f″(x，y)； 步骤3图像导数计算 对步骤2中得到的旋转后的去噪后的降质图像f″(x，y)，计算主方向的导数Da(x，y)和副方向的导数Db(x，y)： Da(x，y)＝f″(x，y)*La(x，y) Db(x，y)＝f″(x，y)*Lb(x，y) 其中下标a和b分别表示主方向和副方向，La(x，y)和Lb(x，y)分别为主方向和副方向的导数卷积模板，此时主方向θ的角度为0°，副方向θ′的角度为90°，主方向卷积模板La(x，y)＝[-1，1]，副方向卷积模板 此时主方向的导数Da(x，y)和副方向的导数Db(x，y)大小都为M3×M4； 步骤4图像导数的自相关计算 步骤4.1利用公知的自相关运算，对步骤3中得到的主方向导数Da(x，y)和副方向导数Db(x，y)，分别计算主方向导数Da(x，y)的自相关Ca(x，y)和副方向导数Db(x，y)的自相关Cb(x，y)： 其中下标a和b分别表示主方向和副方向， 表示相关操作，此时主方向导数的自相关Ca(x，y)和副方向导数的自相关Cb(x，y)的大小都为(2×M3-1)×(2×M4-1)； 步骤4.2对主方向导数的自相关Ca(x，y)按照与主方向θ垂直的方向进行累加得ACa(x，y)，对副方向导数的自相关Cb(x，y)按照与副方向θ′垂直的方向进行累加得ACb(x，y)，累加方法如下： 其中m和n分别为行坐标和列坐标，且都为大于0的整数，主方向导数自相关的累加ACa(x，y)的大小为1×(2×M4-1)，副方向导数的自相关的累加ACb(x，y)的大小为(2×M3-1)×1； 步骤5模糊参数识别 设主方向模糊长度为d，副方向模糊半径为r，对步骤4中得到的主方向导数自相关的累加ACa(x，y)和副方向导数的自相关的累加ACb(x，y)分别绘制曲线图，在主方向导数自相关的累加ACa(x，y)曲线图中将出现一个正的最高峰kap，在正最高峰两侧，出现负的左侧最低峰kap1和负的右侧最低峰kap2，下标ap，ap1和ap2分别表示主方向导数自相关的累加ACa(x，y)曲线图中正的最高峰，负的左侧最低峰和负的右侧最低峰，两负最低峰kap1和kap2之间的距离s1的值为2d+4r，在ACa(x，y)曲线图中将出现一个正的最高峰kbp，在正最高峰两侧，出现负的左侧最低峰kbp1和负的右侧最低峰kbp2，下标bp，bp1和bp2分别表示副方向导数自相关的累加ACb(x，y)曲线图中正的最高峰，负的左侧最低峰和负的右侧最低峰，两负最低峰kbp1和kbp2之间的距离s2的值为4r，计算主方向模糊长度d和副方向模糊半径r： 步骤6建立复合模糊模型 根据步骤5得到的主方向模糊长度d和副方向模糊半径r，通过对标准运动模糊模型hq(x，y)和标准散焦模型hp(x，y)进行线性运算和均一化处理来建立复合模糊模型h(x，y)，下标q和p分别表示运动模糊和散焦，建立复合模糊模型h(x，y)的方法如下： 步骤6.1对标准运动模糊模型hq(x，y)与标准散焦模型hp(x，y)进行卷积操作，得到复合模糊模型的理论模型hv(x，y)，下标v表示复合模糊模型的理论模型，具体方法如下： 标准散焦模型hp(x，y)的函数形式hpt(i1，i2)为： 其中i1，i2为横坐标和纵坐标，且都为实数，下标pt表示标准散焦模型的数学模型，r为副方向模糊半径且以r为散焦模糊半径，r为大于零的实数，标准散焦模型hp(x，y)是标准散焦模型的函数形式hpt(i1，i2)的矩阵形式，标准散焦模型hp(x，y)中的任意元素为大于或等于零的实数， 假设运动模糊方向水平，标准运动模糊模型hp(x，y)的函数形式hpt(i1，i2)为： 其中i1，i2为横坐标和纵坐标，且都为实数，下标qt表示标准运动模糊模型的数学模型，d为主方向模糊长度且以d为运动模糊长度，d为大于零的实数，标准散焦模型hq(x，y)是标准运动模糊模型的函数模型hqt(i1，i2)的矩阵形式，标准运动模糊模型hq(x，y)中的任意元素为大于或等于零的实数， 对标准散焦模型hp(x，y)和标准运动模糊模型hq(x，y)进行卷积，得到复合模糊模型的理论模型hv(x，y)，将复合模糊模型的理论模型hv(x，y)表示为函数形式hvt(i1，i2)： 其中下标vt表示复合模糊模型的理论模型的函数形式， 步骤6.2复合模糊模型的理论模型hv(x，y)确定后，在能量守恒的前提下，对复合模糊模型的理论模型hv(x，y)进行均一化处理，得到最终的复合模糊模型h(x，y)，方法如下： 步骤6.2.1 计算复合模糊模型的理论模型hv(x，y)中大于零的点的个数w， 步骤6.2.2 将复合模糊模型的理论模型hv(x，y)中大于零的点所对应的值替换为 得到最终的复合模糊模型h(x，y)，将复合模糊模型h(x，y)表示成函数形式ht(i1，i2)： 其中下标t表示复合模糊模型的函数形式， 步骤7图像复原 根据步骤6得到的复合模糊模型h(x，y)，利用公知的维纳滤波进行复原，得到得到复原图像 【当前权利人】东南大学 【当前专利权人地址】江苏省南京市四牌楼2号 【统一社会信用代码】12100000466006770Q 【引证次数】4.0 【被引证次数】28 【他引次数】4.0 【被自引次数】4.0 【被他引次数】24.0 【家族引证次数】4.0 【家族被引证次数】28

未经允许不得转载：http://www.zhongzhencnc.com/1785514908.html