CN102170526B - 散焦模糊核计算及其散焦模糊图像清晰化处理方法 - Google Patents
散焦模糊核计算及其散焦模糊图像清晰化处理方法 Download PDFInfo
- Publication number
- CN102170526B CN102170526B CN201110069135A CN201110069135A CN102170526B CN 102170526 B CN102170526 B CN 102170526B CN 201110069135 A CN201110069135 A CN 201110069135A CN 201110069135 A CN201110069135 A CN 201110069135A CN 102170526 B CN102170526 B CN 102170526B
- Authority
- CN
- China
- Prior art keywords
- centerdot
- defocus
- image
- defocusing blurring
- picture
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种散焦模糊核计算及其散焦模糊图像清晰化处理方法,该方法预先在照相机或摄像机镜头中加载一个设计好的特定样式遮光板,更多地保留图像拍摄中的高频信息,然后通过拍摄同一场景的清晰图像和散焦模糊图像各一张,基于经典的图像退化模型建立这两幅图像间的关系,结合最优的散焦模糊核值具有对称和稀疏分布的特点,有效估计信噪比参数,最后设计一个演化算法计算该模糊图像所在位置的散焦模糊核。改造过的光圈使在这个位置拍摄的所有散焦模糊图像保留更多高频信息,利用所计算的散焦模糊核对它们进行解卷积能够得到更清晰的图像。
Description
技术领域
本发明属于数字图像恢复技术领域,具体涉及散焦模糊核计算及其散焦模糊图像清晰化处理方法。
背景技术
在许多实际环境中使用数码相机拍摄图像常常会得到模糊的图像,如所拍目标快速移动、所拍目标不在焦平面上、场景中大气扰动等。那些因为拍摄目标未处于成像***的焦平面上产生的图像模糊称为散焦模糊,如图1所示。散焦模糊图像清晰化技术在光学、天文学、计算机视觉及其计算机图形学等许多学科中都是一个重要研究方向,在环境检测和物证鉴定等多领域具有非常巨大的应用价值。
通常地,人们通过调节光圈大小控制景深,减小光圈意味着减少进光量,使得信噪比减小,相机的景深增加,因此为了提高信噪比,人们常常使用大光圈的镜头,从而也减小了景深。然而,普通的大光圈镜头会抑制散焦模糊图像的高频信息,从而丢弃了图像的许多重要细节信息,因此许多研究人员致力于设计合适的光圈尽可能多地保存拍摄中的高频图像信息,以改进对散焦模糊图像清晰化效果。一般地,散焦核同光圈形状和散焦程度密切相关,散焦核的形状由光圈形状决定,散焦程度和成像焦平面的距离有关,如图2a-2b所示。
迄今为止,国内外许多学者纷纷投入到模糊图像复原领域的研究中,他们大致从以下两个方面对模糊图像进行复原。如果具备足够的有关原始图像的先验知识,那么为原始图像建立一个数学模型并根据它对模糊图像进行拟合会更有效,此时,模糊图像的复原是一个检测问题。然而在实际生活中,我们往往缺乏原始图像的先验知识,这时就必须根据一些假设对图像的退化过程建立模型,这是一种估计方法。近些年来,人们在研究相机或摄像机成像的过程中发现,开放式光圈(如图2a所示)在图像退化模糊的情况下往往会大幅削减其高频信息,然而高频信息又能够直接反应图像的清晰程度,这便增加了人们事后对模糊图像恢复的难度,影响了图像的复原效果。
发明内容
本发明针对现有技术所存在的缺陷,而提供一种散焦模糊核计算及其散焦模糊图像清晰化处理方法。该方法能够实现对散焦模糊图像的清晰化处理。
为了达到上述目的,本发明采用如下技术方案:
散焦模糊核计算及其散焦模糊图像清晰化处理方法,该方法包括如下步骤:
(1)将用于产生对称散焦模糊核的遮光板***相机的镜头;
(2)利用装有经步骤(1)改造过的镜头的相机拍摄同一场景清晰图像一张,离相机不同距离的散焦模糊图像各一张;
(3)配准图像,估计图像中的信噪比并利用经典的图像退化模型建立清晰图像和模糊图像间的关系;
(4)基于清晰图像和散焦模糊图像设计演化算法计算不同距离的散焦模糊核;
(5)利用计算得到的散焦模糊核对对应散焦位置上的其他模糊图像解卷积,实现散焦模糊图像的清晰化处理。
在本发明的一实施例中,所述步骤(1)中遮光板具有对称形状,且产生的散焦模糊核值具有对称和稀疏分布的特点。
在本发明的一实施例中,所述步骤(2)中拍摄图像时需固定相机和拍摄目标,确定清晰图像的焦平面位置,然后以10厘米为间隔调节相机的成像焦平面远离清晰图像的焦平面位置。
在本发明的一实施例中,所述步骤(3)中对清晰图像进行一定的缩放以匹配对应的散焦模糊图像。
进一步的,所述步骤(3)中使用的经典图像退化模型:g=k*f+n,其中f,g,k分别是清晰图像,模糊图像和散焦模糊核,n是噪声且服从分布N(0,σ2)。
在本发明的一实施例中,所述步骤(4)中使用演化算法对散焦模糊核进行实数编码,计算参考的模糊图像:
其评估函数使用相关系数:
在本发明的一实施例中,所述步骤(5)利用步骤(4)所求的散焦模糊核和常用的约束最小平方滤波或者维纳滤波方法解卷积模糊图像,以此达到对散焦模糊图像的清晰化处理。
本发明根据散焦模糊图像的统计特征设计具有最佳聚焦深度鉴别能力的特殊样式遮光板;对图像的深度信息进行有效估计;基于所得的先验知识设计演化算法计算所需的卷积核;调节相关参数,利用解卷积算法实现对散焦模糊图像的清晰化处理。
本发明可广泛用于安防监控市场,既可以用于生产作为原有摄像头的替代品,也可以对原有摄像头进行升级,极大地节约了建设成本,而且使复原的图像达到更清晰的效果,如图4所示。
附图说明
以下结合附图和具体实施方式来进一步说明本发明。
图1为薄透镜成像示意图。
图2a为普通样式光圈的结构示意图。
图2b为普通样式光圈对应的散焦模糊核形状的示意图。
图3a为本发明中光圈的结构示意图。
图3b为本发明中光圈对应的散焦模糊核形状的示意图。
图4为散焦模糊图像复原效果示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
本发明为复原常见的散焦模糊图像,预先在照相机或摄像机镜头中加载一个设计好的特定样式光圈,更多地保留图像拍摄中的高频信息,然后设计一个演化算法计算该模糊图像所在位置的散焦模糊核,并利用所计算的散焦模糊核对它们进行解卷积能够得到更清晰的图像。
基于上述原理,本发明的具体实施如下:
第一步:将预先设计好的特殊样式遮光板***相机镜头。这些特殊样式的遮光板大都具有对称的形状,主要用于产生对称的散焦模糊核,如图3a-3b所示。对称的散焦模糊核在频域中具有更多的零频率,结合图像中的零频率信息,可以更多地保留高频信息,便于以后散焦模糊图像的清晰化处理。
第二步:固定装有改造过镜头的相机和拍摄目标,即固定相应的拍摄位置和焦距,拍摄同一目标清晰图像一张f,记下清晰图像成像焦平面的位置;沿着光轴方向水平移动拍摄目标远离焦平面,以10cm为间隔拍摄模糊图像一张g,并记录其位置信息。
第三步:利用SIFT特征配准清晰图像和模糊图像,估计图像中的信噪比并利用经典的图像退化模型建立清晰图像和模糊图像间的关系。由于模糊图像是沿着光轴方向移动目标产生,因此模糊图像大小可以认为清晰图像进行等比例缩放得到,如图1所示,其中缩放参数为:
λ和f0分别是薄透镜的f-数和焦距。
清晰图像和模糊图像间关系如下:
g=k*f+n
其中f,g,k分别是清晰图像,模糊图像和散焦模糊核,n是噪声且服从分布N(0,σ2)。信噪比为δ·(|Gx|2+|Gy|2)其中δ是调节参数,Gx和Gy分别是g沿x轴和y轴方向上梯度的福利叶变换。
同随机白噪声不同,自然图像通常具有独特的统计性质,即图像的梯度分布服从长尾分布:
其中h(z)=|z|2是一个长尾分布函数,x是清晰图像,i和j分别是行下标和列下标。
第四步:使用演化算法对散焦模糊核k进行实数编码。散焦模糊核值呈稀疏分布,即核的大多部分值为0,其灰度直方图在0点值最大,衰减快于高斯分布。在以后的演化算法中,用两个长尾分布的线性组合作为真实散焦模糊核值的分布:
其中k(i,j)是散焦核矩阵中第i行第j列的元素,w1,w2,β1和β2分别是模糊图像过程有待确定的参数。
随机生成m个序列
si=a11a12…a1na21a22…a2n…an1…ann
每一个序列si表示一个矩阵
基于矩阵A对称构造散焦核矩阵k:
基于两幅图像的相关系数设计适应值函数:
具体演化算法如下(计算散焦模糊核):
1,初始化当前代数gen=0;
2,根据公式(1)随机生成m个长度为L的序列作为初始种群;
3,for gen=1 to GEN do
4,基于每个序列a11...a1n,a21...a2n,...,an1...ann构造矩阵
5,基于矩阵A对称构造散焦核矩阵k:
6,基于公式(2)计算每个序列A的适应值;
7,从种群中选择适应值最大的前am个序列作为优秀种群S;
8,repeat
9,从S中选择两个序列a和b;
10,对a和b进行5次杂交变异,产生10个新个体;
11,把10个新个体中具有最优适应值的个体加入到优秀种群S中;
12,until|S|=m
13,end for。
第五步:利用所求的散焦模糊核k和常用的约束最小平方滤波或者维纳滤波方法解卷积对应位置的模糊图像g,以此达到对散焦模糊图像的清晰化处理。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (7)
1.散焦模糊核计算及其散焦模糊图像清晰化处理方法,其特征在于,所述方法包括如下步骤:
(1)将用于产生对称散焦模糊核的遮光板***相机的镜头;
(2)利用装有经步骤(1)改造过的镜头的相机拍摄同一场景清晰图像一张,离相机不同距离的散焦模糊图像各一张;
(3)配准图像,估计图像中的信噪比并利用经典的图像退化模型建立清晰图像和模糊图像间的关系;
(4)基于清晰图像和散焦模糊图像设计演化算法计算不同距离的散焦模糊核;
(5)利用计算得到的散焦模糊核对对应散焦位置上的其他模糊图像解卷积,实现散焦模糊图像的清晰化处理。
2.根据权利要求1所述的散焦模糊核计算及其散焦模糊图像清晰化处理方法,其特征在于,所述步骤(1)中遮光板具有对称形状,且产生的散焦模糊核值具有对称和稀疏分布的特点。
3.根据权利要求1所述的散焦模糊核计算及其散焦模糊图像清晰化处理方法,其特征在于,所述步骤(2)中拍摄图像时需固定相机和拍摄目标,确定清晰图像的焦平面位置,然后以10厘米为间隔调节相机的成像焦平面远离清晰图像的焦平面位置。
4.根据权利要求1所述的散焦模糊核计算及其散焦模糊图像清晰化处理方法,其特征在于,所述步骤(3)中对清晰图像进行一定的缩放以匹配对应的散焦模糊图像。
5.根据权利要求1所述的散焦模糊核计算及其散焦模糊图像清晰化处理方法,其特征在于,所述步骤(3)中使用的经典图像退化模型:g=k*f+n,其中f,g,k分别是清晰图像,模糊图像和散焦模糊核,n是噪声且服从分布N(0,σ2),信噪比为δ·(|Gx|2+|Gy|2),其中δ是调节参数,Gx和Gy分别是g沿x轴和y轴方向上梯度的福利叶变换。
7.根据权利要求1所述的散焦模糊核计算及其散焦模糊图像清晰化处理方法,其特征在于,所述步骤(5)利用步骤(4)所求的散焦模糊核和常用的约束最小平方滤波或者维纳滤波方法解卷积模糊图像,以此达到对散焦模糊图像的清晰化处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110069135A CN102170526B (zh) | 2011-03-22 | 2011-03-22 | 散焦模糊核计算及其散焦模糊图像清晰化处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110069135A CN102170526B (zh) | 2011-03-22 | 2011-03-22 | 散焦模糊核计算及其散焦模糊图像清晰化处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102170526A CN102170526A (zh) | 2011-08-31 |
CN102170526B true CN102170526B (zh) | 2012-09-05 |
Family
ID=44491487
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110069135A Expired - Fee Related CN102170526B (zh) | 2011-03-22 | 2011-03-22 | 散焦模糊核计算及其散焦模糊图像清晰化处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102170526B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102708555B (zh) * | 2012-03-27 | 2014-10-08 | 陕西师范大学 | 去除彩色图像散焦模糊的方法 |
CN102867297B (zh) * | 2012-08-31 | 2015-08-19 | 天津大学 | 一种低照度图像采集的数字化处理方法 |
CN104616257A (zh) * | 2015-01-26 | 2015-05-13 | 山东省计算中心(国家超级计算济南中心) | 一种模糊退化数字图像在司法中的复原取证方法 |
CN105118031B (zh) * | 2015-08-11 | 2017-11-03 | 中国科学院计算技术研究所 | 一种恢复深度信息的图像处理的方法 |
CN105574850A (zh) * | 2015-12-07 | 2016-05-11 | 天津大学 | 基于方差比特征的图像模糊区域检测与模糊核计算方法 |
CN106920229B (zh) * | 2017-01-22 | 2021-01-05 | 北京奇艺世纪科技有限公司 | 图像模糊区域自动检测方法及*** |
CN107734294A (zh) * | 2017-09-26 | 2018-02-23 | 中国科学院长春光学精密机械与物理研究所 | 监控图像复原***及方法 |
CN108254853B (zh) * | 2018-01-17 | 2023-08-11 | 宁波舜宇仪器有限公司 | 一种显微成像***及其实时对焦方法 |
CN110517197B (zh) * | 2019-08-16 | 2022-09-09 | 西北工业大学 | 基于参数模型的单幅图像去模糊方法 |
CN114663317A (zh) * | 2022-05-18 | 2022-06-24 | 武汉博观智能科技有限公司 | 违法行为的图像智能复原方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364302A (zh) * | 2008-09-28 | 2009-02-11 | 西安理工大学 | 一种散焦模糊图像的清晰化处理方法 |
CN101968881A (zh) * | 2010-10-27 | 2011-02-09 | 东南大学 | 一种运动模糊和散焦复合模糊的图像复原方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002023241A (ja) * | 2000-07-07 | 2002-01-23 | Tamron Co Ltd | フォーカス情報表示機能を有するカメラ |
-
2011
- 2011-03-22 CN CN201110069135A patent/CN102170526B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101364302A (zh) * | 2008-09-28 | 2009-02-11 | 西安理工大学 | 一种散焦模糊图像的清晰化处理方法 |
CN101968881A (zh) * | 2010-10-27 | 2011-02-09 | 东南大学 | 一种运动模糊和散焦复合模糊的图像复原方法 |
Non-Patent Citations (1)
Title |
---|
熊秉衡等.全息检测中的焦散线现象与干涉条纹的关系.《光学学报》.2004,第24卷(第9期),1219-1223. * |
Also Published As
Publication number | Publication date |
---|---|
CN102170526A (zh) | 2011-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102170526B (zh) | 散焦模糊核计算及其散焦模糊图像清晰化处理方法 | |
CN104091312B (zh) | 一种根据图像频谱信息提取模糊核先验的单透镜成像方法 | |
CN109410146A (zh) | 一种基于Bi-Skip-Net的图像去模糊算法 | |
CN112446380A (zh) | 图像处理方法和装置 | |
CN104217404A (zh) | 雾霾天视频图像清晰化处理方法及其装置 | |
CN103440653A (zh) | 双目视觉立体匹配方法 | |
CN106056629A (zh) | 通过运动物体检测和扩展去除鬼影的高动态范围成像方法 | |
CN105205779A (zh) | 一种基于图像变形的眼部图像处理方法、***及拍摄终端 | |
CN104899834A (zh) | 基于sift算法的模糊图像识别方法及装置 | |
CN114066831B (zh) | 一种基于两阶段训练的遥感图像镶嵌质量无参考评价方法 | |
CN110163818A (zh) | 一种用于海事无人机的低照度视频图像增强方法 | |
CN204719330U (zh) | 波前编码成像*** | |
CN107025660A (zh) | 一种确定双目动态视觉传感器图像视差的方法和装置 | |
CN106097383A (zh) | 一种针对遮挡问题的目标跟踪方法及设备 | |
CN103761712A (zh) | 基于自适应光学***点扩散函数重建的图像盲卷积方法 | |
CN107944354A (zh) | 一种基于深度学习的车辆检测方法 | |
CN105631890B (zh) | 基于图像梯度和相位一致性的失焦图片质量评价方法 | |
CN111260687B (zh) | 一种基于语义感知网络和相关滤波的航拍视频目标跟踪方法 | |
CN106027911A (zh) | 一种对地观测星载可见光传输型相机的在轨调焦方法 | |
CN116363168A (zh) | 一种基于超分辨网络的遥感视频目标跟踪方法及*** | |
CN110097530B (zh) | 基于超像素聚类与联合低秩表示的多聚焦图像融合方法 | |
Werth et al. | Silo: A machine learning dataset of synthetic ground-based observations of leo satellites | |
CN107220945A (zh) | 多重退化的极模糊图像的复原方法 | |
Sun et al. | Adaptive image dehazing and object tracking in UAV videos based on the template updating Siamese network | |
CN106023097A (zh) | 一种基于迭代法的流场图像预处理算法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120905 Termination date: 20180322 |
|
CF01 | Termination of patent right due to non-payment of annual fee |