结合相位一致性、梯度幅值和结构显著性的医学图像质量评价方法
【技术领域】
[0001] 本发明属于医学图像评价领域，尤其是涉及一种医学图像质量评价方法。
【背景技术】
[0002] 优质的医学图像往往有助于提高医师诊断效率，并减少医疗事故的发生，但目前 医学图像质量评价研究要远滞后于自然图像质量评价研宄。这主要受限于硬件设备的约 束，公开数据源的数量有限，从事医学影像工作人员专业背景的制约。目前还没有一个可靠 性较高的方法，因此对医学图像质量评价的研究是一项必要工作。
[0003] 目前，基于统计数据的图像质量评狁方法和基于结构特征的图像质量评价方 法（SSIM)仍是使用最多的医学图像质量评价指标，其中基于统计数据的图像质量评价 方法主要有均方误差（Mean-square Error, MSE)和峰值信噪比（Peak Signal-noise Katies RSNR)。这些评价指标用于度量各种图像的压缩算法的效率，也用于度量其他各稗医 学图像处理技术对图像产生的影响。事实上，这些评价指标与人的主观评价结果的一致性 并不算高，在评价图像的诊断质量时也并不实用。
[0004] 基于统计数据的评价方法主要从数学统计的角度出发，将图像看成孤立点集合， 而忽略局域像素间的相关性。这就决定了该方法在处理结构性失真时不能给出准确的质量 评价结果。同时它将信号误差等同于质量的视觉感知，这也不符合HVS视觉特性。PSNR对 于噪声失真能做出较准确的判断,但对于其他类型失真也并不适用。ssm及其改进方法主 要基于人眼视觉系统高度自适应提取场景中的结构信息的特性来进行图像质量评估。SSM 方法从图像局部亮度、对比度和结构:Ξ：方面来度量参考图像与失真之间的相似性，得到图 像质量图谱后采用求均值的方法获得一个整体的质量分数。该类方法性能要优于早期方 法,但依旧达不到实际应周所期望的可靠性。
【发明内容】
[0005] 为了克服传统医学图像质量评价指标与人的主观评价结果的一致性不太高，在评 价诊断质量时实用性较差的不足，本发明提供了一种可靠性比较高、实用性良好的结合相 位一致性、梯度幅值和结构显著性的医学图像质量评价方法。
[0006] 本发明解决其技术问题所采用的技术方案是：
[0007] -种结合相位一致性、梯度幅值和结构显著性的医学图像质量评价方法，包括如 下步骤：
[0008] 1)使用Log-Gabor小波计算参考图像和失真图像的相位一致性图谱；
[0009] 2)使用Scharr算子计算参考图像和失真图像的梯度幅值图谱；
[0010] 3)结合1)得到的相位一致性图谱和2)得到的梯度幅值图谱得到图像局部质量图 谐；
[0011] 4)计算参考图像的基于人眼视觉特征的结构显著图谱；
[0012] 5)用所述结构显著性图谱作为聚合策略对3)中的图像质量图谱进行加权求和， 采用结构显著性模型SSM来体现人眼对图像结构信息的敏感度，用SSM(x)表示该权重图^ 图像质量分数]:QS(x)可表示为：
[0014] 其中，X表示图像在空域Ω上的每个像素；
[0015] 6)当图像质量分数高于预设的分数阈值时，判定医学图像质量较高。
[0016] 进一步，所述步骤1)中，在2D图像中，每个像素的相位一致性定义为：
[0021] 其中，m和η分别表不方向尺度个数，Aji(x)和Φ ji(x)分别表不在第i个方向第 j个对数尺度上的幅度和周部相位偏差，Π是在第i个方向上的估计噪声，W:{是第i个方 向上的权重函数^ Φ(χ，γ)是在（x，y)位置处的局部相位加权均值。
[0022] 再进一步，所述步骤2)中，图像梯度幅值使用卷积核来表示，使用的梯度算子为 Scharr算子，梯度幅值GM:的定义为：
[0024] 其中，Gx和Gy分别表示图像f(x)在水平和竖直方向上偏导数。
[0025] 更进一步，所述步骤3)中，假设Ref表示参考图像，Dst表示失真图像，PCJP PC2分别为参考图像Ref和失真图像Dst的相位一致性图谱，(;;和G 2为参考图像[tef和失真图 像Dst的梯度幅度图谱，PCS为两幅相位一致性图像的相似性图谱AMS为两幅梯度幅度图的 相似性图谱，Q为失真图像的局部质量图；用PCs 00表示相位一致性图的相似性，用GMs (X) 表>丁《梯度幅度图的相似性，PCs定义如1、：
[0027] 其Φ，(3;是一个小正数常量；
[0028] GMs定义如下：
[0030] 其中，〇2是一个小正数常量；
[0031] 局部质量图谱Q(x)表示为：
[0033] 所述步骤4)中，用Pb(x,y, Θ)来表示（x,y)位置Θ方向成为边缘的概率,Pb的 计算基于亮度，色彩，纹理三个通道的梯度方向信息计算得到；
[0034] 对于亮度通道，使羯被一条直径分割成两半的圆作为窗h通过调整直径的方向^ 计算出两个半圆所对应的直方图g和b用g和h计算（X，y)处两个直方图的距离：
[0036] 纹理和色彩通道也用同样的方式进行处理；
[0037] 然后引入三个尺度来检测精细结构和粗糙结构，最后把多尺度上的局部亮度，色 彩和纹理信息结合起来，全局边缘概率定义为：
[0039] 其中，s是尺度的索弓丨，i是特征通道的索引（亮度，色彩，纹理），Gi，σ (i，s) (x, y, )测量通道i中半径为σ (i，s)，圆心为（X，y)的圆被一条角度为Θ的直径分为两 个半圆之后，两个半圆的直方图之差。Pb(x，y，Θ)测量在（x，y)处位置Θ方向成为边缘的 概率。Di,s和γ的值通过使用BSDS训练图像的F-measure的梯度上升得到。结构显著 性模型SSM由上述全局边緣概率的方法实现。
[0040] 本发明的有益效果主要表现在：获得了与专业评估高度一致的结果，提高客观图 像质量评价方法的可靠性，实用性良好。
【附图说明】
[0041] 图1是结合相位一致性、梯度幅值和结构显著性的医学图像质量评价方法的算法 流程图。
【具体实施方式】
[0042] 下面结合附图对本发明作进一步描述。
[0043] 参照图i，一种结合相位一致性、梯度幅值和结构显著性的医学图像质量评价方 法，包括如下步骤：
[0044] 1)使周在计算机视觉和感知领域应用极为广泛的log-Gabor小波计算参考图像 和失真图像的相位一致性图谱；
[0045] 在2D图像中，每个像素的相位一致性可以定义为；
[0050] 其中m和η分别表示方向尺度个数，Aji (X)和Φ ji(x)分别表示在第i个方向第 j个对数尺度上的幅度和局部相位偏差。Π是在第i个方向上的估计噪声，Wi是第i个方 向上的权重函数。Φ(Χ，Υ)是在（X, y)位置处的局部相位加权均值。
[0051] 2)使用Scharr算子计算参考图像和失真图像的梯度幅值图谱；
[0052] 图像梯度幅值使用卷积核来表示，使用的梯度