技术领域本发明涉及一种图像处理技术，特别是一种基于退化图像序列的图像盲复原方法。
背景技术：
光学成像系统在成像时会受到诸多外界因素的影响，主要包括大气散射、吸收等造成的图像雾化，由于成像系统颤振造成的图像模糊等，若不能很好的消除这些负面因素，将严重影响成像系统的成像质量，损伤所得图像的细节信息，降低其应用价值。光学成像系统在成像时必须使成像器件获得足够的能量，因此光照条件成为其首要关注的因素。在光照条件不良的情况下，成像系统必须采用提高感光度或者延长曝光时间的手段以使成像器件获得足够的曝光量。然而前者会导致所得的图像含有大量噪声；后者则易于使成像系统受到外界因素的影响而发生颤振并造成图像模糊。为缓解成像系统的不稳定性，可考虑为成像系统增加稳像设备，然而系统的制造成本将迅速增长，经济性欠佳。因此，将所得的退化图像进行后期复原就成为一个经济而有效的选择。影响光学系统成像的另一个因素是大气环境，光子从地面景物发出后，要与大气分子或气溶胶粒子进行无数次碰撞，经吸收、反射和散射才能到达成像系统，这些因素会造成成像器件表面上光强分布的显著变化，这种变化首先表现为光能量衰减，可用Beer定律描述；其次表现为由于多次散射造成的光能量重新分布，可用Rayleigh散射和Mie散射理论进行分析，通过Monte-Carlo模拟方法可知其将在成像器件表面形成弥散斑，造成图像模糊。上述两种效应经叠加综合，将使所得的图像在视觉效果上产生对比度及分辨率下降，可称其为图像雾化。在成像条件不良的情况下，上述两种影响因素通常同时存在并产生作用，所得图像的退化是大气雾化与成像系统颤振共同作用的结果。然而，目前国内外对图像复原方面的研究仅局限于雾化或模糊单个影响因素，因此在模型建立、求解方法和处理效果上都存在局限性。在问题求解方面，带雾模糊遥感图像的复原将同时涉及去雾和盲去卷积，而盲去卷积又包含了点扩散函数估计和非盲去卷积两个层面的问题，具有非常典型的病态特征，即所得图像中轻微的测量误差也会在复原图像中得到放大，形成明显的噪声及振铃等负面效应，在基于单幅图像的去卷积问题中，由于已知信息量不足，这些负面效应将更加明显。正则化是缓解盲去卷积问题病态性的一种典型方法，它通过在原问题中引入约束条件，使所得的新问题趋于良性，达到限制负面效应的作用。此外，对于现代光学成像系统而言，实现连续成像模式并非难事，在面对成像系统不稳定及大气环境不佳的成像情况时，为了避免获取较暗且带有明显噪声的图像，可不提高成像器件感光度，而是选择延长曝光时间在连续成像模式下获取同一场景的带雾模糊图像序列，凭借其包含的丰富图像信息，采用合适的联合优化估计算法，将有望获得高质量的复原效果。
技术实现要素：
本发明的目的在于提供一种基于退化图像序列的图像盲复原方法。本发明能够在仅已知带雾模糊退化图像序列的情况下，去除图像雾化和模糊，提高图像对比度和清晰度，同时增加图像细节，复原效果好。本发明的过程包括以下步骤：步骤1，截取配准后的每一幅退化序列图像中的重叠区域获得待复原图像序列；步骤2，获取待复原图像序列中每幅图像的透过率和大气光值，基于该透过率和大气光值分别获取平均透过率和平均大气光值；步骤3，所获取的待复原图像序列中的每一幅图像中减去平均大气光干扰项，获得初步去雾结果的模糊图像序列；步骤4，用高斯概率模型对噪声序列发生的概率进行建模、用自然图像梯度先验概率模型对复原图像发生的概率进行建模、用基于l1范数的稀疏概率模型对模糊核序列发生的概率进行建模，对三个子模型相乘得到基于模糊图像序列的盲复原贝叶斯后验概率模型，并取负自然对数；步骤5，通过轮换迭代方法，获得复原模糊核序列；所述每次迭代中，先固定上次迭代所得的复原模糊核序列，对复原图像进行优化，然后固定优化所得的复原图像，对所有复原模糊核序列中的每一个模糊核依次进行优化；其中，模糊图像序列任意一幅图像和高斯型模糊核分别为初始化的复原图像和复原模糊核序列；步骤6，对所获得的最终复原模糊核序列和模糊图像序列进行去卷积后除以平均透过率得到最终复原图像。本发明与现有技术相比，具有以下优点：(1)本发明优点在于利用带雾模糊退化图像序列，在贝叶斯后验概率框架下，采用高斯概率模型、自然图像梯度先验概率模型和基于l1范数的稀疏概率模型分别对噪声序列、复原图像和复原模糊核序列发生的概率进行建模，经综合形成图像复原模型，具有很高的准确度；(2)为了对所得的模型进行求解，采用了一种轮换迭代优化方法，能够获得清晰度高、细节丰富的图像复原效果。下面结合说明书附图对本发明做进一步描述。附图说明图1为本发明实施例的流程示意图。图2～图4为本发明实施例的原始退化图像序列。图5～图7为本发明实施例将图2～图4经过配准后得到的重叠区域图像序列。图8～图10为本发明实施例将图5～图7减去大气干扰项后得到的图像序列。图11～图13为本发明实施例所得的退化序列图像的最终复原模糊核。图14为本发明实施例的最终复原图像。具体实施方式结合图1，本实施例基于退化图像序列的图像盲复原方法的实施步骤如下：步骤1，对退化序列图像进行配准，截取每一幅图像中的重叠区域得到待复原图像序列；步骤2，用暗原色先验方法，计算待复原图像序列中每幅图像的透过率以及大气光值，并分别计算所得透过率的平均值和所得大气光值的平均值，分别定义为平均透过率和平均大气光值；步骤3，将步骤1所得的待复原图像序列中的每一幅图像中减去平均大气光干扰项，得到序列图像的初步去雾结果，即为模糊图像序列；所述大气光干扰项用平均大气光值和平均透过率计算获得；步骤4，用高斯概率模型对噪声序列发生的概率进行建模、用自然图像梯度先验概率模型对复原图像发生的概率进行建模、用基于l1范数的稀疏概率模型对模糊核序列发生的概率进行建模；步骤5，将步骤4得到的三个子模型相乘得到基于模糊图像序列的盲复原贝叶斯后验概率模型，并取负自然对数得到待优化问题；步骤6，用步骤3得到的模糊图像序列中的任意一幅图像和高斯型模糊核分别对复原图像和复原模糊核序列进行初始化，并设置最大迭代次数；步骤7，在每次迭代中，固定上次迭代所得的复原模糊核序列，对复原图像进行优化；步骤8固定优化所得的复原图像，对所有复原模糊核序列中的每一个模糊核依次进行优化；步骤9，判断迭代次数是否小于最大迭代次数，若小于最大迭代次数，则重复执行步骤7和步骤8，否则判定得到最终复原模糊核序列；步骤10，利用步骤9所得的最终复原模糊核序列对步骤3所得的模糊图像序列进行去卷积操作；步骤11，对步骤10得到的结果除以步骤2得到的平均透过率，得到最终复原结果。如图2～图4所示，退化图像序列中每一幅图像的退化程度均不一样，且成像视场存在一定的位移。首先需要对序列图像进行配准，配准后截取每一幅图像中的重叠区域，如图5～图7所示。序列图像的成像时间间隔非常短，在序列图像的成像时间内，可以认为大气状态没有发生改变，即序列图像中的透过率和大气光值不变。因此在估计得到每幅图像的透过率值和大气光值后，用取平均的方法可以减小误差。本实施例采用暗原色先验的方法，估计图像的透过率值和大气光值。图像的退化，一方面是大气的干扰造成的，另一方面是由于成像平台颤振引起的。整个像质退化过程可以用公式表示为：yi=hi⊗(tif)+Ai(1-ti)]]>式中，y、t、h和A分别表示待复原图像序列、大气透过率序列、颤振模糊核序列和大气光值序列，i表示图像索引，I表示清晰图像。上式中每一幅图像透过率的计算步骤为：步骤2.1，根据计算原始退化图像的暗原色yidark，其中，y表示待复原图像序列，i表示序列图像的索引，x表示图像的像素索引，min操作符表示取最小值，c表示图像的r、g、b颜色通道，yic表示原始退化图像的c通道颜色分量，Ω(x)表示以像素点x为中心的方形区域，z为Ω(x)区域内的像素索引；步骤2.2，从所述暗原色yidark的所有点中找出最亮的1％个点，记这些像素点的位置集合为Qi。在原始退化图像yi中，找出集合Qi对应位置的像素中的最亮点，以该点的像素值作为大气光值Ai；步骤2.3，用暗原色先验的方法计算透过率的表达式为：(ti)x=1-ωminc∈{r,g,b