本发明涉及一种超声图像定量评估系统，包括：图像预处理模块，用于对获取的关节超声图像进行预处理；图像分割模块，用于对关节超声图像进行图像分割；定量分析指标提取模块，用于对关节超声图像进行定量分析指标提取；参数联合分析模块，用于对关节超声图像进行多参数联合分析。本发明通过在超声图像上自动分割出软骨前后壁和内部区域，并对分割出的感兴趣区域进行定量分析，提取出与超声回波强度、声衰减、软骨厚度和散射子统计学分布等相关的定量分析指标，再利用多变量分析的方法，对多种指标进行综合分析，从而实现对关节软骨病变程度的无创及定量评估。
图像预处理模块，用于对获取的关节超声图像进行预处理，包括对获取的关节超声图像进行插值、归一化和衰减补偿的预处理，具体为
1)对关节超声图像进行深度方向的插值，以使不同深度采集到的关节超声图像具有相同的像素大小；
2)用各帧最大强度值对每帧关节超声图像作归一化，消除增益对超声回声强度的影响；同时，使用同一个超声系统，在相同的成像配置下，采集均匀仿体的关节超声图像，用仿体信号对软骨相应深度的信号作归一化，用于校正深度对超声回声强度的影响；
3)使用衰减补偿方法对关节超声图像中软骨上层覆盖的软组织引起的声衰减进行补偿；
定量分析指标提取模块，用于对关节超声图像进行定量分析指标提取；
参数联合分析模块，用于对关节超声图像进行多参数联合分析。
2.根据权利要求1所述的超声图像定量评估系统，其特征在于，对关节超声图像进行图像分割包括：从所述超声图像中分割出软骨前壁、软骨后壁和软骨内部所包含的感兴趣区域。
3.根据权利要求2所述的超声图像定量评估系统，其特征在于，对所述关节超声图像进行定量分析指标提取包括：根据关节软骨病变的B超表现，从预处理和图像分割后得到的感兴趣区域中提取出与超声回声强度、声衰减、软骨厚度和散射子统计学分布相关的定量分析指标。
4.根据权利要求3所述的超声图像定量评估系统，其特征在于，对关节超声图像进行多参数联合分析包括：对定量分析指标提取模块提取出的各定量分析指标，使用联合分析或机器学习的方法进行综合分析。
5.根据权利要求2所述的超声图像定量评估系统，其特征在于，使用手动或者基于区域、边界、数学形态学、小波理论或者神经网络的图像分割算法，从关节超声图像中分割出软骨前壁、软骨后壁和软骨内部所包含的感兴趣区域。
6.根据权利要求3所述的超声图像定量评估系统，其特征在于，从预处理和图像分割后得到的感兴趣区域中提取出与超声回声强度、声衰减、软骨厚度和散射子统计学分布相关的定量分析指标包括：
7.根据权利要求6所述的超声图像定量评估系统，其特征在于，在进行所述步骤1)时，具体包括以下步骤：①计算软骨前壁和软骨后壁的超声回声强度平均值和标准差，用以反映软骨前壁和软骨后壁的平均回声强度和均匀度，其中计算公式如下：式中，Raw和SVRaw分别表示软骨前壁的平均回声强度和均匀度；Rpw和SVRpw分别表示软骨后壁的平均回声强度和均匀度； 和 分别表示软骨前壁和软骨后壁第i位置处的超声回声强度；m表示软骨沿着横向的长度分量；
②计算软骨内部的超声回声强度平均值和标准差，用以反映软骨内部的超声回声强度和均匀度，其中计算公式如下：式中，Rcartilage和SVRcartilage分别表示软骨内部的超声回声强度和均匀度； 表示软骨内部第i位置处的超声回声强度；
③计算软骨前壁和软骨后壁的超声回波强度与软骨内部的超声回声强度的比值，用以反映软骨前壁和软骨后壁与软骨内部的回声对比度，其计算公式如下：式中，contrastaw表示软骨前壁与软骨内部的回声对比度；contrastpw表示软骨后壁与软骨内部的回声对比度。
8.根据权利要求6所述的超声图像定量评估系统，其特征在于，在进行所述步骤2)时，具体包括以下内容：根据关节超声图像，获得软骨上边界的位置信息，对软骨上边界各位置到超声探头表面的垂直距离进行高通滤波获得滤波后的距离信号，并求取滤波后的距离信号的标准差用以反映软骨表面的粗糙度，其中粗糙度的计算公式如下：式中，URI表示软骨表面的粗糙度；di表示第i位置处滤波后的距离信号；表示滤波后距离信号的平均值；m表示软骨沿着横向的长度分量。
9.根据权利要求6所述的超声图像定量评估系统，其特征在于，在进行所述步骤3)时，具体包括以下步骤：①从软骨前壁和软骨后壁下方各选出长度为n的感兴趣区域，并将感兴趣区域等分成若干份，对感兴趣区域内部各位置处的超声回声强度和对应深度进行线性拟合，用拟合直线的斜率表征回声幅度衰减系数，从而获得软骨前壁和软骨后壁所包含的感兴趣区域各段范围及整个范围内的回声幅度衰减系数，用于反映与软骨病变相关的声衰减特性的改变；
②同理，获得软骨内部所包含的感兴趣区域各段范围及整个范围内的回声幅度衰减系数。
10.根据权利要求6所述的超声图像定量评估系统，其特征在于，在进行所述步骤4)时，具体包括以下内容：计算不同位置处，软骨前壁和软骨后壁的垂直距离，所有垂直距离的平均值、标准值和最小值可分别反映软骨厚度的平均值、厚度均匀性和最小厚度。
11.根据权利要求6所述的超声图像定量评估系统，其特征在于，在进行所述步骤5)时，具体包括以下内容：根据散射子统计学分布模型，利用超声包络信号拟合出软骨内部的散射子分布参数，其中散射子统计学分布模型所采用的Nakagami分布模型如下：式中，r为软骨内部的局部超声回声强度，且r≥0；n为形状参数，且n＞0；Ω为尺度参数，且Ω＞0；Γ(n)表示关于形状参数n的Gamma函数。
一种超声图像定量评估系统
技术领域
[0001] 本发明涉及一种超声图像定量评估系统，具体是关于一种可以用于评估关节软骨病变的超声图像定量评估系统。
背景技术
[0002] 多种病因导致的骨关节炎是老年人最常见的慢性关节退行性疾病。骨关节炎的早期识别和准确的关节退化分级有助于对病情进行及时和适当的干预，从而控制疾病进展，显著改善预后。临床上，常规X线摄影(conventional radiography,CR)常用于评价骨关节结构改变，但其无法显示软组织病变。磁共振成像(magnetic resonance imaging,MRI)可用于检测关节软骨、滑膜、半月板和软组织病变，但其价格昂贵，扫描时间长，且对软骨早期病变不敏感。而超声诊断具有实时、便携、易于使用、价格低廉等优势，并且可以显示骨与软骨的结构改变，以及关节内的炎性改变，因此适用于大样本人群中骨关节炎的筛查。但是，目前临床超声诊断存在主观性强，无法给出定量结果等问题。
[0003] 基于金标准关节镜检查，关节病变程度可以分为三度：I度下，探针可触及关节软骨软化，可观察到少量的表面纤维化，闭合性的软骨分离及泡状改变；II度下，可观察到少量的关节软骨纤维束样改变，呈蟹肉状外观；III度下，可观察到软骨坏死脱落，软骨下骨外露并出现“象牙化”。但是关节镜检查是一种有创的检查技术，不适合骨关节炎的常规筛查。
[0004] 目前，关于使用定量超声的方法来评价软骨退化的研究主要应用于离体实验或者在体关节镜检查中。这些研究通过使用高频探头来提高轴向分辨率，通过从超声射频信号中提取与软骨表面声反射特性和粗糙度等相关的定量指标，来表征病变引起的软骨表面微结构的变化，或者对超声射频信号进行功率谱分析，提取出与软骨背向散射特性、声衰减等相关的指标。然而对于临床超声检查，由于探头频率与成像深度之间的折衷，离体实验或者关节镜检查中使用的高频探头无法探及关节软骨所在深度，因此不适用于关节软骨的临床检查。并且，目前大部分临床使用的超声机器无法提供超声射频数据。因此，离体实验和关节镜检查中常用的反映软骨表面微结构变化的指标和根据超声射频数据频谱分析提取得到的指标无法应用于临床超声检查中。
[0005] 因此，亟待为临床提供一种无创/定量的关节软骨病变评估方法，以更好的应用于大范围人群的病变筛查和分级诊断。
发明内容
[0006] 针对上述问题，本发明的目的是提供一种超声图像定量评估系统，该系统通过对超声图像进行定量分析，提取对软骨病变敏感的指标，并利用多参数联合分析对多种指标综合分析，对关节软骨病变程度进行无创及定量的评估。
[0007] 为实现上述目的，本发明采取以下技术方案：一种超声图像定量评估方法，包括：
[0008] 一、对获取的关节超声图像进行预处理的步骤；二、对关节超声图像进行图像分割的步骤；三、对关节超声图像进行定量分析指标提取的步骤；四、对关节超声图像进行多参数联合分析的步骤。
[0009] 所述的超声图像定量评估方法，优选地，所述对获取的关节超声图像进行预处理的步骤，包括：对获取的关节超声图像进行插值、归一化和衰减补偿的预处理。
[0010] 所述的超声图像定量评估方法，优选地，所述对关节超声图像进行图像分割的步骤，包括：从所述超声图像中分割出软骨前壁、软骨后壁和软骨内部所包含的感兴趣区域。
[0011] 所述的超声图像定量评估方法，优选地，所述对所述关节超声图像进行定量分析指标提取的步骤，包括：根据关节软骨病变的B超表现，从预处理和图像分割后得到的感兴趣区域中提取出与超声回声强度、声衰减、软骨厚度和散射子统计学分布相关的定量分析指标。
[0012] 所述的超声图像定量评估方法，优选地，所述对关节超声图像进行多参数联合分析的步骤，包括：对步骤三中提取出的各定量分析指标，使用联合分析或机器学习的方法进行综合分析，实现软骨病变的无创定量分级，最终得到软骨病变的分级结果。
[0013] 所述的超声图像定量评估方法，优选地，所述对获取的关节超声图像进行插值、归一化和衰减补偿的预处理，包括以下步骤：1)对关节超声图像进行深度方向的插值，以使不同深度采集到的关节超声图像具有相同的像素大小；2)用各帧最大强度值对每帧关节超声图像作归一化，消除增益对超声回声强度的影响；同时，使用同一个超声系统，在相同的成像配置下，采集均匀仿体的关节超声图像，用仿体信号对软骨相应深度的信号作归一化，用于校正深度对超声回声强度的影响；3)使用衰减补偿方法对关节超声图像中软骨上层覆盖的软组织引起的声衰减进行补偿。
[0014] 所述的超声图像定量评估方法，优选地，使用手动或者基于区域、边界、数学形态学、小波理论或者神经网络的图像分割算法，从关节超声图像中分割出软骨前壁、软骨后壁和软骨内部所包含的感兴趣区域。
[0015] 所述的超声图像定量评估方法，优选地，所述从预处理和图像分割后得到的感兴趣区域中提取出与超声回声强度、声衰减、软骨厚度和散射子统计学分布相关的定量分析指标，包括以下步骤：1)超声回声强度计算的步骤；2)软骨表面粗糙度计算的步骤；3)声衰减速度计算的步骤；4)软骨厚度计算的步骤；5)软骨内部散射子统计学分布计算的步骤。
[0016] 所述的超声图像定量评估方法，优选地，在进行所述步骤1)时，具体包括以下步骤：
[0017] ①计算软骨前壁和软骨后壁的超声回声强度平均值和标准差，用以反映软骨前壁和软骨后壁的平均回声强度和均匀度，其中计算公式如下：
[0018]
[0019]
[0020]
[0021]
[0022] 式中，Raw和SVRaw分别表示软骨前壁的平均回声强度和均匀度；Rpw和SVRpw分别表示软骨后壁的平均回声强度和均匀度； 和 分别表示软骨前壁和软骨后壁第i位置处的超声回声强度；m表示软骨沿着横向的长度分量；
[0023] ②计算软骨内部的超声回声强度平均值和标准差，用以反映软骨内部的超声回声强度和均匀度，其中计算公式如下：
[0024]
[0025]
[0026] 式中，Rcartilage和SVRcartilage分别表示软骨内部的超声回声强度和均匀度；表示软骨内部第i位置处的超声回声强度；
[0027] ③计算软骨前壁和软骨后壁的超声回波强度与软骨内部的超声回声强度的比值，用以反映软骨前壁和软骨后壁与软骨内部的回声对比度，其计算公式如下：
[0028]
[0029]
[0030] 式中，contrastaw表示软骨前壁与软骨内部的回声对比度；contrastpw表示软骨后壁与软骨内部的回声对比度。
[0031] 所述的超声图像定量评估方法，优选地，在进行所述步骤2)时，具体包括以下内容：
[0032] 根据关节超声图像，获得软骨上边界的位置信息，对软骨上边界各位置到超声探头表面的垂直距离进行高通滤波获得滤波后的距离信号，并求取滤波后的距离信号的标准差用以反映软骨表面的粗糙度，其中粗糙度的计算公式如下：
[0033]
[0034] 式中，URI表示软骨表面的粗糙度；di表示第i位置处滤波后的距离信号；表示滤波后距离信号的平均值；m表示软骨沿着横向的长度分量。
[0035] 所述的超声图像定量评估方法，优选地，在进行所述步骤3)时，具体包括以下步骤：
[0036] ①从软骨前壁和软骨后壁下方各选出长度为n的感兴趣区域,并将感兴趣区域等分成若干份，对感兴趣区域内部各位置处的超声回声强度和对应深度进行线性拟合，用拟合直线的斜率表征回声幅度衰减系数，从而获得软骨前壁和软骨后壁所包含的感兴趣区域各段范围及整个范围内的回声幅度衰减系数，用于反映与软骨病变相关的声衰减特性的改变；
[0037] ②同理，获得软骨内部所包含的感兴趣区域各段范围及整个范围内的回声幅度衰减系数。
[0038] 所述的超声图像定量评估方法，优选地，在进行所述步骤4)时，具体包括以下内容：
[0039] 计算不同位置处，软骨前壁和软骨后壁的垂直距离，所有垂直距离的平均值、标准值和最小值可分别反映软骨厚度的平均值、厚度均匀性和最小厚度。
[0040] 所述的超声图像定量评估方法，优选地，在进行所述步骤5)时，具体包括以下内容：
[0041] 根据散射子统计学分布模型，利用超声包络信号拟合出软骨内部的散射子分布参数，其中散射子统计学分布模型所采用的Nakagami分布模型如下：
[0042]
[0043] 式中，r为软骨内部的局部超声回声强度，且r≥0；n为形状参数，且n>0；Ω为尺度参数，且Ω>0；Γ(n)表示关于形状参数n的Gamma函数。
[0044] 本发明由于采取以上技术方案，其具有以下优点：本发明通过在超声图像上自动分割出软骨前后壁和内部区域，并对分割出的感兴趣区域进行定量分析，提取出与超声回波强度、声衰减、软骨厚度和散射子统计学分布等相关的定量分析指标，再利用多变量分析的方法，对多种指标进行综合分析，从而实现对关节软骨病变程度的无创及定量评估。本发明提出的方法同样适合于其它在体组织或器官(如颈动脉斑块、甲状腺、肝脏、肾脏等)病变程度的无创定量评估。
附图说明
[0045] 图1为本发明一实施例提供的方法的流程示意图。
具体实施方式
[0046] 为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0047] 如图1所示，本实施例提供的超声图像定量评估方法，其包括以下步骤：
[0048] 一、对获取的超声图像进行预处理的步骤：
[0049] 在本实施例中，获取的超声图像是按照统一标准操作流程(例如outcome measures in rheumatology，OMERACT准则)通过超声系统采集到标准格式(如DICOM格式)的关节超声图像，然后通过对该超声图像进行插值、归一化和衰减补偿等预处理，尽量减少不同检查者之间，由于成像深度、增益、软骨上层组织声衰减等因素造成的对回波信号强度的影响，具体为：
[0050] 1)对超声图像进行深度方向的插值，以使不同深度采集到的超声图像具有相同的像素大小；
[0051] 2)用各帧最大强度值对每帧超声图像作归一化，消除增益对超声回声强度的影响；同时，使用同一个超声系统，在相同的成像配置下，采集均匀仿体的超声图像，用仿体信号对软骨相应深度的信号作归一化，用于校正深度对超声回声强度的影响；
[0052] 3)使用适当的衰减补偿方法(例如最佳功率谱移位估计法)对超声图像中软骨上层覆盖的软组织引起的声衰减进行补偿。
[0053] 二、对超声图像进行图像分割的步骤：
[0054] 在本实施例中，使用手动或者基于区域、边界、数学形态学、小波理论或者神经网络的图像分割算法，从超声图像中分割出软骨前壁、软骨后壁和软骨内部所包含的感兴趣区域(region of interest，ROI)。
[0055] 三、对超声图像进行定量分析指标提取的步骤：
[0056] 由于临床超声探头频率和采集系统的限制，依赖于高频探头和射频数据分析的定量分析方法无法应用于临床采集得到的图像分析中。因此，在本实施例中，根据关节软骨病变的B超表现，从预处理和图像分割后得到的ROI中提取出与超声回声强度、声衰减、软骨厚度和散射子统计学分布等相关的定量分析指标，具体为：
[0057] 1)超声回声强度计算：
[0058] ①计算软骨前壁和软骨后壁的超声回声强度平均值和标准差，用以反映软骨前壁和软骨后壁的平均回声强度和均匀度，其中计算公式如下：
[0059]
[0060]
[0061]
[0062]
[0063] 式中，Raw和SVRaw分别表示软骨前壁的平均回声强度和均匀度；Rpw和SVRpw分别表示软骨后壁的平均回声强度和均匀度； 和 分别表示软骨前壁和软骨后壁第i位置处的超声回声强度；m表示软骨沿着横向的长度分量。
[0064] ②计算软骨内部的超声回声强度平均值和标准差，用以反映软骨内部的超声回声强度和均匀度，其中计算公式如下：
[0065]
[0066]
[0067] 式中，Rcartilage和SVRcartilage分别表示软骨内部的超声回声强度和均匀度；表示软骨内部第i位置处的超声回声强度。
[0068] ③计算软骨前壁和软骨后壁的超声回波强度与软骨内部的超声回声强度的比值，用以反映软骨前壁和软骨后壁与软骨内部的回声对比度，其计算公式如下：
[0069]
[0070]
[0071] 式中，contrastaw表示软骨前壁与软骨内部的回声对比度；contrastpw表示软骨后壁与软骨内部的回声对比度。
[0072] 2)软骨表面粗糙度计算：
[0073] 根据超声图像，获得软骨上边界的位置信息，对软骨上边界各位置到超声探头表面的垂直距离进行高通滤波获得滤波后的距离信号，并求取滤波后的距离信号的标准差用以反映软骨上壁(即软骨表面)的粗糙度，其中粗糙度的计算公式如下：
[0074]
[0075] 式中，URI表示软骨表面的粗糙度；di表示第i位置处滤波后的距离信号；表示滤波后距离信号的平均值。
[0076] 3)声衰减速度计算：
[0077] ①从软骨前壁和软骨后壁下方各选出长度为n(本实施例中n＝0.1mm～5mm)的ROI,并将ROI等分成若干份(例如四份)，对ROI内部各位置处的超声回声强度和对应深度进行线性拟合，用拟合直线的斜率表征回声幅度衰减系数，从而获得软骨前壁和软骨后壁所包含的ROI各段范围及整个范围内的回声幅度衰减系数，用于反映与软骨病变相关的声衰减特性的改变。
[0078] ②同理，获得软骨内部所包含的ROI各段范围及整个范围内的回声幅度衰减系数。
[0079] 4)软骨厚度计算：
[0080] 通过现有方法计算不同位置处，软骨前壁和软骨后壁的垂直距离，所有垂直距离的平均值、标准值和最小值可分别反映软骨厚度的平均值、厚度均匀性和最小厚度。
[0081] 5)软骨内部散射子统计学分布计算：
[0082] 根据散射子统计学分布模型(本实施例中采用Nakagami分布模型)，利用超声包络信号拟合出软骨内部的散射子分布参数(在Nakagami分布模型中，拟合参数为n和Ω)，其中Nakagami分布模型如下：
[0083]
[0084] 式中，r为软骨内部的局部超声回声强度，且r≥0；n为形状参数，且n>0；Ω为尺度参数，且Ω>0；Γ(n)表示关于形状参数n的Gamma函数。
[0085] 四、对超声图像进行多参数联合分析的步骤：
[0086] 对步骤三中提取出的各定量分析指标，使用联合分析或机器学习的方法进行综合分析，实现软骨病变的无创定量分级，最终得到软骨病变的分级结果(即I、II、III度软骨病变)。
[0087] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
