[发明专利]基于四元数小波变换幅值与相位的图像清晰度评价方法

说明书

技术领域

本发明涉及一种清晰度的评价方法，特别是基于四元数小波变换幅值与相位的图像清晰度评价方法。

背景技术

随着智能手机、数码相机等数字成像技术的不断发展，产生的图像越来越多。在图像处理领域，如何筛选出更好的图像并指导成像系统的控制，使系统始终工作在最佳状态，相应的图像清晰度评价已然成为最重要的一环。

图像清晰度评价方法分为主观图像清晰度评价方法和客观图像清晰度评价方法。由于图像清晰与否最终都是由人的主观感受来做出判断，因此主观图像清晰度评价方法是最好的评价方法。但每个人的认知水平不同，也没有公认的评价标准，故客观图像清晰度评价方法更受关注。

到目前为止，图像的客观清晰度评价方法可以分为基于空域的和基于频域的评价方法。一般基于空域的评价方法是通过清晰度评价函数对图像本身进行计算，得到清晰度的值。灰度图像的方差是空域最简单的一种图像清晰度评价方法，类似的还有求灰度图像梯度、特征值和拉普拉斯能量和的方法。这些评价方法计算简单，效率高，实时性好，但它们不能准确区分图像清晰度的微小变化。因此，研究人员提出了基于频域的图像清晰度评价方法。该类方法利用多尺度变换工具先把图像转换到频域，然后再提取其特征以得到清晰度评价函数。由于多种变换工具各有其优势，研究者们提出了基于小波、轮廓波和非下采样轮廓波等的清晰度评价方法。这些变换工具可以通过其频域系数将图像清晰度的微小变化放大。相较基于空域的评价方法，基于频域的评价方法对图像清晰度的评价更准确，但是对不同图像进行评价时不够全面。因为由四元数小波变换分解系数得到的幅值与相位信息具有很好的纹理表示能力，所以其为图像清晰度评价提供了一种新的途径。

发明内容

针对上述现有技术中存在的缺陷或不足，本发明的目的在于，提出一种基于四元数小波变换幅值与相位的图像清晰度评价方法，从而达到评价效果更准确，适用范围更广的目的。

实现本发明的技术方案是，基于四元数小波变换幅值与相位的图像清晰度评价方法，包括如下步骤：

1)准备待进行评价的灰度图像；

2)对灰度图像进行四元数小波变换(QWT)，得到变换系数；

3)由变换系数根据四元数代数计算得到各子带(LL、LH、HL、HH子带)对应的幅值与相位信息，其中低频LL子带水平方向、垂直方向、对角方向的(θ，φ，ψ)相位分别记为PhaseL_hor、PhaseL_ver、PhaseL_dia，高频LH子带水平方向的θ相位记为PhaseLH_hor，高频HL子带垂直方向的φ相位记为PhaseHL_ver，高频HH子带对角方向的ψ相位记为PhaseHH_dia；

4)由步骤3)得到的LL子带幅值计算其水平、垂直和正对角方向的梯度，分别记为GradL_hor、GradL_ver、GradL_dia；

5)将步骤3)和步骤4)得到的梯度值与对应方向的相位值进行处理之后得到清晰度评价指标。

所述步骤1)具体为准备大小为N×N的方形灰度图像X。

所述步骤2)具体为对X进行QWT处理，每个子带都将其表示为x＝a+bi+cj+dk的形式。其中，i、j、k为四元数的三个虚数单位并且两两正交，a、b、c、d则是X对应的QWT变换系数。

所述步骤3)具体为根据计算公式对各子带进行计算得到对应的幅值与相位，并记上相应的标记。该计算公式为：

其中，m为幅值，分别为对应子带水平方向、垂直方向、对角方向的相位。

所述步骤4)具体为由步骤3)得到的LL子带幅值m根据计算公式进行计算得到m在水平、垂直和正对角方向的梯度，分别记为GradL_hor、GradL_ver、GradL_dia。该计算公式为：

其中，f(x,y)指m位于(x,y)处的幅值。

所述步骤5)具体为由步骤3)和步骤4)计算得到的梯度和相位根据计算公式计算得到一组清晰度评价值。该计算公式为：

其中，x和y代表当前像素的位置。