一种简单镜头成像的多通道PSF标定方法技术

本发明专利技术公开了一种简单镜头成像的多通道PSF标定方法，涉及图像复原技术领域，包括：提取彩色模糊图像的R、G、B通道图像；利用盲卷积图像复原算法分别估计R、G、B通道图像对应的PSF；基于每个通道的PSF对每个通道图像做去模糊处理；对每个通道单独做色散校正；将处理之后的R、G、B通道图像进行融合得到最终的清晰复原图像。本方法不同于传统方法对三个通道使用同一个PSF，分通道估计的PSF更加准确，去模糊效果更好，同时有效解决了色散在传统去模糊过程进一步加剧的问题，很好地提高最终图像复原质量。

A multichannel PSF calibration method for simple lens imaging

The invention discloses a multi-channel PSF a simple lens calibration method, which relates to the technical field of image restoration, including: R, G, B extraction channel color image blur image restoration algorithm; R, G and B respectively to estimate the channel image and the corresponding PSF image by blind deconvolution; each channel PSF for each channel to do image processing based on fuzzy; on each channel individually do dispersion correction; R, G, B channel image processing after fusion to get the final clear image. The method is different from the traditional method of three channels using the same PSF channel estimation of PSF is more accurate and fuzzy to better effect, and effectively solve the dispersion in the traditional fuzzy process to further aggravate the problem, greatly improve the final image quality.

【技术实现步骤摘要】
一种简单镜头成像的多通道PSF标定方法
本专利技术主要涉及到数字图像处理领域，特指一种简单镜头成像的多通道PSF标定方法。
技术介绍
目前，单反相机在人们的日常生活中发挥着越来越重要的作用。然而，为弥补单反镜头中镜片的几何畸变和像差，进一步提高成像质量，单反镜头的设计日益复杂，甚至包含数十个独立的光学器件。复杂镜头在提高成像质量的同时，也会增加镜头的体积和重量，导致镜头成本大大提高。近年来，随着计算摄影技术的发展，简单镜头结合后期图像复原算法正逐渐成为相机设计领域以及图像处理领域一个新的研究方向。简单镜头只包含一个镜片，受镜头像差和色散的影响，由简单镜头直接拍摄的图像是模糊的，而且边缘处存在色散，从而导致图像质量不高。需要首先标定得到简单镜头的点扩散函数(PointSpreadFunction，PSF)，PSF中包含简单镜头的像差等模糊信息，可用于后期的图像复原算法中。传统的PSF标定方法将彩色图像的三个通道作统一处理，即将彩色图像转换成灰度图像，通过标定得到一个PSF，并将这个PSF用于彩色图像三个通道的去模糊处理。这种方法主要存在以下两个问题：(1)受简单镜头中镜片球面形状的影响，由简单镜头得到的彩色图像的三个通道的PSF并不完全相同，如果统一处理会导致PSF估计不准；(2)色散效应在每个通道中都存在，而且在去模糊过程会进一步加剧，如果PSF估计不准，最终复原图像中的色散效应会加剧，从而影响图像复原质量。在中国专利申请号ZL201510222290.9中介绍了一种基于对称性的单透镜计算成像PSF快速标定方法，这种方法并没有考虑彩色图像R、G、B三个通道的差异性，对彩色图像估计的是一个统一PSF，这会影响PSF的估计精度，从而影响最终的图像复原效果。
技术实现思路
本专利技术为克服传统PSF标定方法中估计不准，色散效应在去模糊过程中加剧的问题，旨在提供一种提高简单镜头PSF标定精度后得到清晰复原图像的方法。一种简单镜头成像的多通道PSF标定方法，包括以下步骤：步骤一：提取原始模糊图像的R、G、B通道图像；步骤二：用盲卷积图像复原方法分别标定不同通道的PSF；其中采用的盲卷积图像复原算法为基于最大后验概率MAP的盲卷积图像复原算法，在最大后验概率模型下，盲卷积图像复原问题的统计学模型可以表述为：argmaxP(k,i|b)＝argmaxP(b|i,k)P(i)P(k)其中，k表示简单镜头的PSF；i表示清晰图像；b表示由简单镜头直接得到的模糊图像；P(k,i|b)表示在模糊图像b已知的条件下，与模糊图像b对应的PSF和清晰图像分别为k和i的概率；P(b|i,k)表示如果已知PSFk和清晰图像对应的模糊图像为b的概率；P(i)表示对清晰图像已知的先验概率；P(k)表示PSF的先验概率；考虑到对数与乘积的转换关系以及对数函数的单调性，对上述公式左右两端分别去负对数：-logP(k,i|b)＝-logP(b|k,i)-logP(i)-logP(k)则上述问题转化为能量最小化问题，左边的目标函数可以定义为：其中，表示数据拟合项；||▽(i)||p1表示图像先验；||k||p2表示PSF先验，λ1和λ2分别表示图像先验和PSF先验在目标函数中的权重。步骤三：根据所标定的PSF分别对每个通道进行图像复原；基于步骤二中得到的PSF，分别利用非盲卷积图像复原算法对R、G、B通道图像进行去模糊处理，得到对应的去模糊图像；所采用的非盲卷积图像复原目标函数可以表述为：其中，k表示简单镜头的PSF；i表示清晰图像；b表示由简单镜头直接得到的模糊图像；表示数据拟合项；||▽i||1表示图像先验，其中▽i表示图像的导数；λ表示PSF先验在目标函数中的权重。步骤四：对每个通道进行色散校正；步骤五：将处理后的R、G、B通道图像合并成最终复原图像。进一步的，所述步骤五的融合直接在matlab软件中完成。进一步的，所述步骤二中λ1＝0.45、λ2＝0.55。进一步的，所述步骤三中λ＝0.8。本专利技术充分考虑到简单镜头中镜片的球面特性，对每个通道单独估计PSF，并分通道进行色散校正，这种方法可以提高PSF的估计精度，而且很好地避免色散效应在去模糊过程进一步加剧，能提高最终复原图像的清晰度。本专利技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本专利技术的实践了解到。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本专利技术所述方法的流程示意图；图2为由简单镜头得到的模糊图像对应的R、G、B通道图像；图3为本专利技术R、G、B通道图像分别标定出的PSF；图4为本专利技术最终复原的清晰图像。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。一种简单镜头成像的多通道PSF标定方法，如图1所示，包括如下步骤：步骤一：提取原始模糊图像的R、G、B通道图像，图像如图2所示；步骤二：用盲卷积图像复原方法分别标定不同通道的PSF：即利用盲卷积图像复原方法，分别对模糊图像的R、G、B通道图像进行PSF标定，得到对应的R、G、B通道图像的PSF。其中采用的盲卷积图像复原算法为基于最大后验概率MAP的盲卷积图像复原算法，在最大后验概率模型下，盲卷积图像复原问题的统计学模型可以表述为：argmaxP(k,i|b)＝argmaxP(b|i,k)P(i)P(k)(1)其中，k表示简单镜头的PSF；i表示清晰图像；b表示由简单镜头直接得到的模糊图像；P(k,i|b)表示在模糊图像b已知的条件下，与模糊图像b对应的PSF和清晰图像分别为k和i的概率；P(b|i,k)表示如果已知PSFk和清晰图像i，对应的模糊图像为b的概率；P(i)表示对清晰图像已知的先验概率；P(k)表示PSF的先验概率；考虑到对数与乘积的转换关系以及对数函数的单调性，对公式(1)左右两端分别去负对数：-logP(k,i|b)＝-logP(b|k,i)-logP(i)-logP(k)则上述问题转化为能量最小化问题，左边的目标函数可以定义为：其中，表示数据拟合项；||▽(i)||p1表示图像先验；||k||p2表示PSF先验，λ1和λ2分别表示图像先验和PSF先验在目标函数中的权重。优选的，取λ1＝0.45，λ2＝0.55，因为是估计PSF，所以PSF的先验权重λ2大一些，这样使目标函数更加精准，在求解过程中更加趋于真实值，最终复原图像的清晰度更高。R、G、B三个通道所估计出的PSF如图3所示。步骤三：根据所标定的PSF分别对每个通道进行图像复原。即基于步骤二中得到的PSF，分别利用非盲卷积图像复原算法对R、G、B通道图像进行去模糊处理，得到对应的去模糊图像；所采用的非盲卷积图像复原目标函数可以表述为：其中，k表示简单镜头的PSF；i表示清晰图像；b表本文档来自技高网...

【技术保护点】
一种简单镜头成像的多通道PSF标定方法，其特征在于，包括以下步骤：步骤一：提取原始模糊图像的R、G、B通道图像；步骤二：用盲卷积图像复原方法分别标定不同通道的PSF；其中采用的盲卷积图像复原算法为基于最大后验概率MAP的盲卷积图像复原算法，在最大后验概率模型下，盲卷积图像复原问题的统计学模型可以表述为：arg max P(k,i|b)＝arg max P(b|i,k)P(i)P(k)其中，k表示简单镜头的PSF；i表示清晰图像；b表示由简单镜头直接得到的模糊图像；P(k,i|b)表示在模糊图像b已知的条件下，与模糊图像b对应的PSF和清晰图像分别为k和i的概率；P(b|i,k)表示如果已知PSF k和清晰图像对应的模糊图像为b的概率；P(i)表示对清晰图像已知的先验概率；P(k)表示PSF的先验概率；考虑到对数与乘积的转换关系以及对数函数的单调性，对上述公式左右两端分别去负对数：‑log P(k,i|b)＝‑log P(b|k,i)‑log P(i)‑log P(k)则上述问题转化为能量最小化问题，左边的目标函数可以定义为：

【技术特征摘要】
1.一种简单镜头成像的多通道PSF标定方法，其特征在于，包括以下步骤：步骤一：提取原始模糊图像的R、G、B通道图像；步骤二：用盲卷积图像复原方法分别标定不同通道的PSF；其中采用的盲卷积图像复原算法为基于最大后验概率MAP的盲卷积图像复原算法，在最大后验概率模型下，盲卷积图像复原问题的统计学模型可以表述为：argmaxP(k,i|b)＝argmaxP(b|i,k)P(i)P(k)其中，k表示简单镜头的PSF；i表示清晰图像；b表示由简单镜头直接得到的模糊图像；P(k,i|b)表示在模糊图像b已知的条件下，与模糊图像b对应的PSF和清晰图像分别为k和i的概率；P(b|i,k)表示如果已知PSFk和清晰图像对应的模糊图像为b的概率；P(i)表示对清晰图像已知的先验概率；P(k)表示PSF的先验概率；考虑到对数与乘积的转换关系以及对数函数的单调性，对上述公式左右两端分别去负对数：-logP(k,i|b)＝-logP(b|k,i)-logP(i)-logP(k)则上述问题转化为能量最小化问题，左边的目标函数可以定义为：其中，...

【专利技术属性】
技术研发人员：不公告发明人，
申请(专利权)人：长沙全度影像科技有限公司，
类型：发明
国别省市：湖南,43