一种基于超声回波射频信号的肺超声特殊征象B线识别分类方法

一种基于超声回波射频信号的肺超声特殊征象b线识别分类方法
技术领域
1.本发明属于超声检测及成像领域，具体涉及一种基于超声回波射频信号的肺超声特殊征象b线识别分类方法。


背景技术：

2.随着现代超声医学的不断发展进步，曾经被视作超声盲区的肺脏超声检查也由于其无辐射、可床旁诊断等优点成为临床上检测肺脏疾病、监测重症病人生理指标的重要手段。b线是超声检测、诊断肺部疾病的一种重要超声伪影，在超声影像上表现为从胸膜线产生并延伸至屏幕底部的离散、激光样垂直高回声伪影，并与肺滑动同步移动。
3.肺超声影像中b线数量的增加与肺水肿、肺纤维化、肺炎等肺病变之间显著的相关性已经成为临床上对患者进行床旁无创、快速诊断的重要依据，医生可以按照肺超声检查国际推荐标准，快速评估患者的病情，并做出相应的诊断。
4.然而，在相同的评估标准之下，不同的医生对同一患者的评估、甚至同一医生在不同时间对同一患者的评估都可能会产生一定偏差，这与评估者的临床经验、评估时的生理状态等因素有关。因此，能够辅助医生对肺超声进行评估的算法可以大大提高医生诊断的稳定性和准确性。已经有一些研究人员提出的基于超声图像的深度学习算法在辅助评估中的到了较好的准确性和敏感性，但该方法计算复杂，并且是基于超声图像进行计算的，只利用了幅度和空间信息。


技术实现要素：

5.本发明的目的在于克服上述现有技术的缺点，提供一种基于超声回波射频信号的肺超声特殊征象b线识别分类方法，以解决现有技术中通过超声图像对评估结果有偏差的问题。
6.为达到上述目的，本发明采用以下技术方案予以实现：
7.一种基于超声回波射频号的肺超声特殊征象b线识别分类方法，包括以下步骤：
8.步骤1，从多帧超声图像提取出含有b线的肺超声回波信号的射频数据，通过脉冲逆转方法提取出每一个射频数据的二次谐波数据，将多帧二次谐波数据组成三维矩阵；三维矩阵中每一列的二次谐波数据为一个波束；
9.步骤2，通过滑动窗提取三维矩阵中每一个波束的信号包络能量和信息熵，分别设定信号包络能量的阈值一和信息熵的阈值二，将滑动窗提取出的信息熵高于阈值二且包络能量小于阈值一的射频信号视为噪声，去除视为噪声的射频信号，得到去除噪声后的三维矩阵；
10.步骤3，从去除噪声后的三维矩阵中，提取出每一个波束的信息熵和排列熵；
11.步骤4，标记出每帧超声图像的b线和非b线对应的波束，b线的波束标记值为1，非b线的波束标记值为0；
12.步骤5，随机选择若干个波束对应的信息熵、排列熵以及该波束的b线标记值；将选择出的信息熵、排列熵和b线标记组成训练样本集x；
13.步骤6，将训练样本集x输入至非线性svm，通过训练样本集x训练非线性svm，得到判断b线的分类模型；所述非线性svm中通过核函数rbf得到分类模型；
14.步骤7，通过分类模型识别超声图像中的b线。
15.本发明的进一步改进在于：
16.优选的，步骤1的具体过程为：
17.1)从含有b线的肺超声回波信号射频数据中，提取并分离互为反向脉冲的两组回波信号；
18.2)将两组回波信号按对应位置相加得到脉冲逆转射频数据；
19.3)读取脉冲逆转射频数据头文件中的中心频率；
20.4)将脉冲逆转射频数据输入巴特沃斯带通滤波器，获得射频数据的二次谐波数据，将多帧的二次谐波数据形成三维矩阵。
21.优选的，所述巴特沃斯带通滤波器的频带宽度为步骤3)中心频率的40％～80％，通过频带宽度确定巴特沃斯带通滤波器的通带频率和阻带频率。
22.优选的，步骤2中，信号包络能量的计算公式为：
[0023][0024][0025][0026]
x(t)＝a(t)cos(w0t θ(t))(14)
[0027][0028][0029]
式(1)中，x(t)为原始射频信号，为射频信号经希尔伯特变换后的信号；式(3)中为令原信号为实部，希尔伯特变换后的信号作为虚部构造的解析信号；将式(4)带入式(3)得到的解析信号的绝对值为信号包络；式(6)中e为包络能量。
[0030]
优选的，步骤2中，信息熵的计算公式为：
[0031][0032]
式(7)中，h表示单个波束射频信号的信息熵，n表示单波束的射频信号的点数，f表示射频信号的概率密度，x为单波束对应的射频信号数据。
[0033]
优选的，步骤3中，提取每一个波束排列熵的过程为：
[0034]
1)设定延迟时间为1，嵌入维度为4，逐个提取单波束的重构矩阵；
[0035]
2)将重构矩阵每一行按照升序排列，得到新升序序列；
[0036]
3)提取新升序序列对于原重构矩阵的位置索引所构成的符号序列及其出现的次
数；
[0037]
4)每种符号序列出现次数除以序列总数作为该序列出现的概率；
[0038]
5)计算排列熵。
[0039]
优选的，所述排列熵的计算公式为：
[0040][0041]
式(11)中h
pe
表示单波束射频信号的排列熵。
[0042]
优选的，步骤5中，所述核函数rbf核函数为
[0043]
k(x,y)＝exp(-γ||x-y||2)
ꢀꢀꢀ
(19)。
[0044]
优选的，步骤5中，训练非线性svm的过程为：
[0045]
1)将训练样本集x作为原始空间；
[0046]
2)在原始空间中引入影射函数，使得训练样本集x中的函数k(ei,pi)满足下式：
[0047]
k(ei,pj)＝φ(ei)
·
φ(pj)
ꢀꢀꢀ
(20)
[0048]
式(13)中，k表示核函数，φ表示核函数的映射函数。
[0049]
3)在原始空间中输入式(13)，在特征空间中映射成新向量的内积；
[0050]
4)计算非线性svm在特征空间的最大分隔超平面，得到分类模型。
[0051]
优选的，步骤7中，识别出超声图像中的b线后进行标记。
[0052]
与现有技术相比，本发明具有以下有益效果：
[0053]
本发明提出了一种基于超声回波射频信号的肺超声特殊征象b线识别方法，该方法具体包括以下步骤：步骤一，输入在b模式下含有b线的肺超声回波信号射频数据，利用脉冲逆转技术提取二次谐波，采集多帧二维数据生成二次谐波数据三维矩阵；步骤二，利用滑动窗提取信号包络能量和信息熵作为去噪参数，根据实际情况调整阈值，去除噪声；步骤三，逐个波束提取对应射频数据的信息熵和排列熵作为特征参数；步骤四，分别标记b线和非b线位置对应波束，b线标记为1，非b线标记为0；步骤五，通过随机选择的方法选择一定数量的射频数据样本所对应特征参数作为训练样本集x；步骤六，将x输入训练非线性svm(支持向量机)，得到分类模型h；步骤七，用训练好的分类模型h进行b线检测；步骤八，在超声图像上标记出b线位置。本发明提出利用超声原始射频信号对肺超声b线进行识别，不仅可以利用信号幅度和空间信息，而且可以利用信号频率和相位等信息，不仅提高了信息的利用率，增加了识别的可靠性，而且避免不同超声仪器图像参数设置的影响，能够得到普适性更高的检测算法。本发明的优点在于还包括：
[0054]
(1)具有较高的b线识别准确率，能够辅助临床医生进行b线识别；
[0055]
(2)由于是基于超声回波射频信号的识别，受仪器参数影响小，具有很好的普适性；
[0056]
(3)本发明基于超声射频数据，相比于基于超声图像的识别方法，可以保留更多的组织信息，提高识别的可靠性；
[0057]
(4)本发明基于射频数据，相比于其他识别方法，射频数据不会受到超声仪器图像参数设置的影响，可以适用于各种仪器设备。
[0058]
(5)本发明提取多个射频数据特征，包括信息熵和排列熵特征，作为识别是否为b线的依据，可以减少误判的情况，提高b线识别的准确率。
[0059]
(6)本发明使用非线性svm，对于复杂的射频信号特征有很好的适配性，可以实现很好的识别准确率。
[0060]
进一步的，支持向量机(support vector machine，svm)是一种监督学习方式的二分类模型，其决策边界是对学习样本求解的最大边距超平面，是一种具有稀疏性和稳健性的分类器，可以对线性可分的样本进行较准确的二元分类。对于线性不可分的学习样本，非线性支持向量机利用核函数，通过非线性映射算法将二位线性不可分的样本映射到高维特征空间中，使得样本点再高维空间线性可分，采用与线性支持向量机相同的间隔最大化原则学习得到非线性支持向量机模型。一般来说升维意味着更高的计算复杂度，但由于核函数的展开定理，无需得到非线性映射的显性表达式，而是在高维特征空间中建立线性学习机，所以与线性模型相比，非线性模型并不显著增加计算复杂度，是一种对非线性数据分类效果很好，应用很广的常用分类模型。本发明基于超声原始射频数据，结合svm技术，能够得到很好的b线识别结果。
附图说明
[0061]
图1是基于超声回波射频信号的肺超声特殊征象b线识别分类方法总流程；
[0062]
图2是用训练集训练非线性svm后，对测试集b线分类的准确率折线图；
[0063]
图3是肺部超声b线位置识别及标注典型结果图；
[0064]
其中，(a)是肺部超声图像；(b)是b线位置识别及标注后超声图像。
具体实施方式
[0065]
下面结合附图对本发明做进一步详细描述：
[0066]
本发明公开了一种基于超声回波射频信号的肺超声特殊征象b线识别分类方法，参见图1，该识别分类方法包括以下步骤：
[0067]
步骤一，输入在b模式下含有b线的肺超声回波信号射频数据，利用脉冲逆转技术提取二次谐波，采集多帧二维数据生成二次谐波数据三维矩阵；具体的包括以下步骤：
[0068]
1)对于含有b线的肺超声回波信号射频数据，分别提取并分离每一帧射频数据中互为反向脉冲的两组回波信号；
[0069]
2)将1)中提取的两组回波信号按对应位置相加得到脉冲逆转射频数据；
[0070]
3)读取射频数据头文件中的中心频率；
[0071]
4)生成巴特沃斯带通滤波器，根据两倍中心频率即二次谐波频率选择合适的频带宽度，具体的，频带宽度为中心频率的40～80％，基于频带宽度设置合适的通带频率和阻带频率；
[0072]
5)将脉冲逆转射频数据输入巴特沃斯带通滤波器，获得二次谐波数据，将多帧的二次谐波数据形成三维矩阵。
[0073]
步骤二，利用滑动窗提取信号包络能量和信息熵作为去噪参数，根据实际情况调整阈值，去除噪声。具体包括以下步骤：
[0074]
针对的三维矩阵中的每一列二次谐波数据，根据二次谐波数据尺寸设置合适的滑动窗窗长和步长，依据各波束数据点数确定滑动窗总数，波束为三维矩阵中一列二次谐波数据，取该波束的信号包络能量和信息熵作为该段二次谐波数据的去噪参数，通过该去噪
参数对该波束去除噪声。
[0075]
具体的，其中信号包络能量的计算方法按照公式(1)-(6)进行：
[0076][0077][0078][0079]
x(t)＝a(t)cos(w0t θ(t))
ꢀꢀꢀ
(24)
[0080][0081][0082]
式(1)中，x(t)为原始射频信号，为射频信号经希尔伯特变换后的信号。式(3)中为令原信号为实部，希尔伯特变换后的信号作为虚部构造的解析信号。将式(4)带入式(3)得到的解析信号的绝对值即为信号包络。式(6)中e为包络能量。
[0083]
信息熵的计算方法按照公式(7)进行：
[0084][0085]
式(7)中，h表示单个波束射频信号的信息熵，n表示单波束的射频信号的点数，f表示射频信号的概率密度，x为单波束对应的射频信号数据。
[0086]
其中，逐个波束提取对应射频数据的信息熵时计算概率密度的方法，具体按照以下步骤进行：
[0087]
1)逐个计算单波束射频数据的范围，平均分为1000段；
[0088]
2)计算每段内数据个数，除以单波束射频数据总量，作为该范围段内数据的概率密度。
[0089][0090]
提取单波束各个滑动窗的平均包络能量与信息熵后，去除射频数据噪声的方法，按照以下方法进行：根据超声仪器设定的中心频率、声能等影响射频信号特征的参数调整合适的阈值，分别设定信息熵的判定标准为阈值一，满足包络能量的判定标准为阈值二，利用噪声信息熵数值大、包络能量数值小的特征对二次谐波数据进行去噪处理。将满足信息熵数值高于阈值二、包络能量小于阈值一的滑动窗段的射频信号视为噪声，滤除该段信号。
[0091]
步骤三，对去噪后的波束逐个提取对应射频数据的信息熵和排列熵作为特征参数。逐个波束提取对应射频数据信息熵的方法，按照公式(7)进行。
[0092]
所述步骤三中，逐个波束提取对应射频数据的排列熵的方法，具体按照以下步骤进行：
[0093]
1)设定延迟时间为1，嵌入维度为4，逐个提取单波束的重构矩阵；
[0094][0095]
式(8)中，y表示单波束射频信号关于时间的重构矩阵，x表示单波束对应射频信号数据，t表示延迟时间，m表示嵌入维数，n表示单波束的射频信号的点数且k＝n-(m-1)t。
[0096]
2)将重构矩阵每一行按照升序排列，得到新的升序序列；
[0097]
s(l)＝{j1,j2,
…
,jm},l＝1,2,
…
,k.且k≤m！
ꢀꢀꢀ
(29)
[0098]
式(9)中s表示将y矩阵中各重构分量按照升序排列所得到的向量中各元素在y矩阵的元素位置列索引构成的符号序列，l表示符号序列s对应的重构分量在矩阵y中的位置，j表示升序向量中各元素在y矩阵的元素位置列索引。
[0099]
3)提取新升序序列对于原重构矩阵的位置索引所构成的符号序列及其出现的次数；
[0100]
4)每种符号序列出现次数除以序列总数作为该序列出现的概率；
[0101]
a＝{p1,p1,
…
,pk}
ꢀꢀꢀ
(30)
[0102]
式(10)中a表示不同符号序列出现概率的集合，pk表示重构分量对应符号序列出现的概率。
[0103]
5)根据式(11)计算排列熵。
[0104][0105]
式(11)中h
pe
表示单波束射频信号的排列熵。
[0106]
步骤四，分别标记b线和非b线位置对应波束，b线标记为1，非b线标记为0；
[0107]
所述步骤四中标记b线和非b线位置对应波束，在超声图像中按照从胸膜线产生并延伸至屏幕底部的离散、激光样垂直高回声的伪影，b线特征，
[0108]
表现对所有波束对应射频信号进行标记，b线标为1，非b线标为0。
[0109]
步骤五，通过随机选择的方法选择一定数量的射频数据样本所对应特征参数作为训练样本集x；训练样本集x中包括三列，分别为信息熵、排列熵，以及是否为b线的标记。
[0110]
步骤六，将x输入训练非线性svm，得到分类模型h；
[0111]
步骤六所述非线性svm，设置c-svc的惩罚系数为1000，使用rbf核函数：
[0112]
k(x,y)＝exp(-γ||x-y||2)
ꢀꢀꢀ
(32)
[0113]
式(12)中k表示rbf核函数，x、y表示原始空间样本，γ表示超参数，||x-y||表示向量范数。
[0114]
步骤六所述训练非线性svm模型，按照以下步骤进行：
[0115]
1)将计算得到的信息熵e和排列熵p数据作为原始空间(样本集x)；
[0116]
2)引入一个映射函数φ，使得对所有ei∈e,pj∈p，函数k(ei,pi)满足条件：
[0117]
k(ei,pj)＝φ(ei)
·
φ(pj)
ꢀꢀꢀ
(33)
[0118]
式(13)中k表示核函数，φ表示核函数的映射函数。
[0119]
3)原始空间中样本输入式(13)中后，生成了升维的特征空间，获得特征空间中映
射成的新向量的内积。
[0120]
4)依照线性svm计算方法，求出非线性svm在特征空间的最大分隔超平面，得到非线性svm分类模型h。
[0121]
步骤七，用训练好的分类模型h进行b线检测；
[0122]
步骤八，在超声图像上标记出b线位置。
[0123]
步骤八所述在超声图像上标记，具体按照以下步骤进行：
[0124]
1)输入射频信号数据，提取识别结果；
[0125]
2)依照识别结果以及射频信号与超声图像位置的对应关系确定超声图像中b线的位置；
[0126]
3)利用半透明彩色条带标记出超声图像中b线位置。
[0127]
如图1所示，本发明公开一种基于超声回波射频信号的肺超声特殊征象b线识别分类方法，该方法具体包括以下步骤：
[0128]
进一步的，所述的一种基于超声回波射频信号的肺超声特殊征象b线识别分类方法，所述步骤一中，输入在b模式下含有b线的肺超声回波信号射频数据，利用提取二次谐波的方法，按照以下方法进行：
[0129]
1)分别提取每一帧射频数据中互为反向脉冲的两组回波信号；
[0130]
2)将1)中提取的两组回波信号按对应位置相加得到脉冲逆转射频数据；
[0131]
3)读取射频数据头文件中的中心频率；
[0132]
4)生成巴特沃斯带通滤波器，根据两倍中心频率即二次谐波频率选择合适的频带宽度，设置合适的通带频率和阻带频率；
[0133]
5)将脉冲逆转射频数据输入滤波器，获得二次谐波数据。
[0134]
为了验证本发明的有效性，进行了分类实验。
[0135]
实验采用了100帧在b模式脉冲逆转下含有b线的射频数据，每一帧数据规格为512
×
2336二维矩阵，将数据经过脉冲逆转相加、去噪、滤波处理后提取射频数据二次谐波信息熵和排列熵作为输入样本，将其中70帧作为训练集，30帧作为测试集，将训练集输入非线性svm训练获得分类模型h，利用分类模型h对测试集数据进行b线识别检测，最后根据识别结果在超声图像中标注b线位置。30帧测试集数据的分类实验准确率如表1所示：
[0136]
表1
[0137][0138]
图2是利用十折交叉验证方法检验分类模型性能的折线图，即将所有的特征参数随机分为十份，九份用于训练，一份用于识别，所得的具体准确率如表2所示：
[0139]
表2
[0140]
[0141]
由表1可以看出，该分类模型h对测试集的识别结果准确率均在90％以上，最高准确率为98.047％，平均准确率为95.377％。由表2可以看出，本模型识别准确率较好，平均准确率为94.631％。可以看出，本发明对b线识别的准确率具有较高的水平，有很好的临床参考价值。
[0142]
图3分别是对原始超声图像b线位置识别及标注前后对比图。可以看到在原图的基础上，本发明使用蓝色半透明条带清楚的标记出了b线位置，可以准确、稳定、快速地辅助临床医生进行诊断。
[0143]
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。