基于稀疏表示的非一致图像盲复原方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN105957025A(43)申请公布日2016.09.21(21)申请号201610254970.3(22)申请日2016.04.21(71)申请人天津大学地址300072天津市南开区卫津路92号(72)发明人杨爱萍王南梁斌何宇清魏宝强(74)专利代理机构天津市北洋有限责任专利代理事务所12201代理人李丽萍(51)Int.Cl.G06T5/00(2006.01)权利要求书2页说明书7页附图3页(54)发明名称基于稀疏表示的非一致图像盲复原方法(57)摘要本发明公开了一种基于稀疏表示的非一致图像盲复原方法，包括：根据相机三维抖动模型，并结合自然图像的过完备字典表示构建非一致图像模糊退化模型；输入待复原的模糊图像和过完备字典，求得初始稀疏系数，并对参数进行初始化设置；将自然图像的过完备字典表示，以及模糊核和稀疏系数的稀疏性作为模型的正则约束，并采用交替迭代方法将非一致盲复原模型的求解转化为多个简单的子问题，从而实现对模糊图像y的盲复原。本发明提出的算法对自然条件下获取的模糊图像有更好的复原效果，复原图像细节明显、不失真、噪声低，具有更好的视觉效果，且求解方法有可扩展性。CN105957025ACN105957025A权利要求书1/2页1.一种基于稀疏表示的非一致图像盲复原方法，包括以下步骤：步骤一、构建图像盲复原模型，包括：1-1)根据相机三维抖动模型，构建式(1)所示的非一致图像模糊退化模型：式(1)中，y为模糊图像，x为原始清晰图像，ε为加性噪声，表示清晰图像在偏移角度为θk时的偏移图像；偏移角度为θk时，清晰图像像素点j处的像素值xj在Cijk的作用下对应到模糊图像的像素点i处；wk是偏移角度为θk时的权值，称之为模糊核；当模糊核或清晰图像已知时，式(1)所示模型表示为如下两个线性模型：y＝∑kCijkwkx+ε(2)y＝∑jCijkxjw+ε(3)过完备字典对图像的稀疏表示如式(4)所示：式(4)中，Ri是提取图像第i个图像块的矩阵操作符，D是过完备字典，αi是第i个图像块在过完备字典表示下的稀疏表示系数；1-2)根据上述构建的非一致模糊退化模型，并结合自然图像的过完备字典表示，构建式(5)所示的图像复原模型：式(5)中，第一项是由非一致模糊模型构造的数据保真项；第二项和第三项分别是过完备字典表示和稀疏系数αi的正则化约束项；第四项是模糊核约束；其中，模糊核w具有稀疏特性，满足归一化和非负性，采用l1范数对该模糊核w进行约束；η(ηi)、λ(λi)和γ(γi)为平衡参数，η(ηi)＝0.05、λ(λi)＝0.1、γ(γi)的取值范围为4.0～5.5；原始清晰图像x和模糊图像y均为过完备字典D排列的向量形式；步骤二、输入待复原的模糊图像y和过完备字典D，其中，过完备字典D是采用块邻域梯度字典学习方法对Berkeleysegmentationdataset数据库中的200幅图像随机提取20000个8×8像素点图像块进行自适应学习得到的；步骤三、采用交替迭代方法对待复原的模糊图像y进行盲复原，包括：3-1)首次迭代时，稀疏系数α由min||α||0s.t.得到，并设原始清晰图像x为输入的模糊图像y；3-2)根据稀疏系数α和原始清晰图像x估计模糊核w，将上述式(5)简化为：并采用迭代收缩软阈值算法估计模糊核w；2CN105957025A权利要求书2/2页3-3)根据模糊核w和稀疏系数α估计原始清晰图像x，将模糊图像y分解为多个不同的重叠块，对每个重叠块进行估计原始清晰图像x的操作，上述式(5)简化为：式(7)中，A＝∑kCijkwk，利用快速傅里叶变换算法估计原始清晰图像x；3-4)根据原始清晰图像x和模糊核w估计稀疏系数α，上述式(5)简化为：式(8)中，第二部分为l1范数的加权优化形式，权值ωi采用稀疏系数反函数形式；由于过完备字典方法中各个图像块之间是相互独立的，因此，将式(8)简化为式(9)，从而将原始清晰图像x稀疏系数的求解问题转化为各个图像块的稀疏系数的求解问题，采用YALL1算法估计稀疏系数α；3-5)当迭代次数小于2时，返回步骤3-2)；否则，执行步骤3-6)；3-6)计算相邻两次迭代求得的模糊核w之间的差值平方差，若该差值平方差的值小于10-4，则停止迭代，输出的图像即为原始清晰图像x，否则返回步骤3-2)。3CN105957025A说明书1/7页基于稀疏表示的非一致图像盲复原方法技术领域[0001]本发明涉及一种计算机图像处理方法，尤其是涉及一种图像复原方法。背景技术[0002]成像设备在采集图像时，人为抖动或设备固有的机械抖动会使采集到的图像出现整体模糊、目标物体边界不明确、细节信息丧失等问题，给地形勘测、目标识别、自动导航等计算