超声多普勒包络确定方法、装置、计算机设备和存储介质与流程

1.本发明涉及医学信号处理技术领域，特别是涉及一种超声多普勒包络确定方法、装置、计算机设备、存储介质和计算机程序产品。


背景技术：

2.近年来，随着超声多普勒技术的迅速发展，超声多普勒能够无损伤地检测人体血管中的血流状况，从而被广泛用来诊断血管疾病。大量研究表明，血管疾病会从多普勒信号最大频率曲线的变化上体现出来。因此，能否获得准确的超声多普勒血流信号的正反向最大频率曲线(即频谱包络)对于血管疾病的诊断起着重要的作用。传统的频谱包络是由经验丰富的超声医师绘制，不仅效率低下，且需要大量人工干预。
3.相关技术中，多普勒频谱的包络提取技术主要包括两种：一种是利用图像处理技术，比如图像的边缘检测等来求得多普勒频谱的包络，不过这种基于图像处理的方式存在计算量大、复杂度高以及实时性不高的问题。另一种是利用超声多普勒血流信号的功率谱密度积分曲线的频率特性，来进行最大频率的估计，而这种基于数字频率估计的方法虽然计算量较小，计算效率高，但是却也存在一定的局限性，如果存在反向血流速度时，则无法正常工作。


技术实现要素：

4.基于此，有必要针对上述技术问题，提供一种能够同时确定正反向波形的超声多普勒包络确定方法、装置、计算机设备、存储介质和计算机程序产品。
5.本技术提供了一种超声多普勒包络确定方法。所述方法包括：
6.获取超声图谱中每一列信号对应的功率谱密度信号；
7.对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点；
8.连接所述超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，得到所述超声图谱对应的波形包络线。
9.进一步地，所述对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点，包括：基于所述功率谱密度信号对应的第一频率区间对所述功率谱密度信号进行第一积分处理，得到对应的第一积分曲线；根据所述第一积分曲线确定所述功率谱密度信号的最大能量点；根据所述最大能量点确定所述功率谱密度信号的两个最小能量点；基于所述两个最小能量点对应的第二频率区间对相应区间的功率谱密度信号进行第二积分处理，得到对应的第二积分曲线；根据所述第二积分曲线确定所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
10.进一步地，所述根据所述第一积分曲线确定所述功率谱密度信号的最大能量点，包括：获取所述第一积分曲线的首尾两个端点，通过第一参考线连接所述两个端点，所述第一参考线为直线；获取所述第一参考线与所述第一积分曲线之间的第一交叉点，将所述第
一交叉点作为所述功率谱密度信号的最大能量点。
11.进一步地，所述根据所述最大能量点确定所述功率谱密度信号的两个最小能量点，包括：根据所述最大能量点将所述功率谱密度信号划分为两个信号区间；分别获取每个信号区间内纵坐标最小的点，作为所述功率谱密度信号的最小能量点。
12.进一步地，所述根据所述第二积分曲线确定所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点，包括：获取所述第二积分曲线的首部第一端点和尾部第二端点，通过第二参考线连接所述第一端点和第二端点，所述第二参考线为直线；确定所述第二参考线与所述第二积分曲线之间的第二交叉点、正向频率点和反向频率点，所述正向频率点为所述第二积分曲线到所述第二参考线的正向最大距离对应的频率点，所述反向频率点为所述第二积分曲线到所述第二参考线的反向最大距离对应的频率点；获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值，当所述第一平均值差值大于或等于预设阈值时，将所述正向频率点确定为正向最大频率点；获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值，当所述第二平均值差值大于或等于预设阈值时，将所述反向频率点确定为反向最大频率点。
13.进一步地，所述方法还包括：当所述第一平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述正向频率点，并返回执行获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值的步骤；当所述第二平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述反向频率点，并返回执行获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值的步骤。
14.本技术还提供了一种超声多普勒包络确定装置，所述装置包括：
15.功率谱密度信号获取模块，用于获取超声图谱中每一列信号对应的功率谱密度信号；
16.频率点确定模块，用于对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点；
17.包络线获取模块，用于连接所述超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，得到所述超声图谱对应的包络线。
18.本技术还提供了一种计算机设备。所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以上实施例中任意一项所述的方法中的步骤。
19.本技术还提供了一种计算机可读存储介质。所述计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以上实施例中任意一项所述的方法中的步骤。
20.本技术还提供了一种计算机程序产品。所述计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以上实施例中任意一项所述的方法中的步骤。
21.上述超声多普勒包络确定方法、装置、计算机设备、存储介质和计算机程序产品，通过获取超声图谱中每一列信号对应的功率谱密度信号，对功率谱密度信号进行积分处理，以获取功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点，进而连接超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，从而得到超声图
谱对应的波形包络线。本实施例通过将超声信号转换为功率谱密度信号，并基于功率谱密度信号同时确定超声图谱对应的正向包络线和反向包络线，从而实现得到由正向包络线和反向包络线组成的波形包络线，且具有计算过程简单、计算效率高等优势。
附图说明
22.图1为本技术中超声多普勒包络确定方法的流程示意图。
23.图2为本技术中信号转换示意图。
24.图3为本技术中对功率谱密度信号进行积分处理步骤的流程示意图。
25.图4为本技术中第一积分曲线的示意图。
26.图5为本技术中确定功率谱密度信号中最大能量点和最小能量点的示意图。
27.图6为本技术中第二积分曲线的示意图。
28.图7为本技术中波形包络线的示意图。
29.图8为本技术中超声多普勒包络确定装置的结构框图。
30.图9为本技术中计算机设备的内部结构图。
具体实施方式
31.为使本发明的技术方案和有益效果能够更加明显易懂，下面通过列举具体实施例的方式进行详细说明。其中，附图不一定是按比例绘制的，局部特征可以被放大或缩小，以更加清楚的显示局部特征的细节；除非另有定义，本文所使用的技术和科学术语与本技术所属的技术领域中的技术和科学术语的含义相同。
32.本技术提供了一种超声多普勒包络确定方法，本实施例以该方法应用于终端进行举例说明，可以理解的是，该方法也可以应用于服务器，还可以应用于包括终端和服务器的系统，并通过终端和服务器的交互实现。具体地，如图1所示，上述方法包括以下步骤：
33.步骤102，获取超声图谱中每一列信号对应的功率谱密度信号。
34.其中，超声图谱是通过超声多普勒探头作为传感器采集的超声信号所形成的图谱，超声图谱中的每一列信号包括若干个信号点。对于每一个信号点，基于对应的频率和能量都能在功率谱中唯一地表示出来。
35.在本实施例中，如图2所示，对于超声图谱中任一列信号n对应的若干个信号点，根据该列信号中每个信号点分别对应的频率和能量可以将其转换为通过功率谱表示的功率谱密度信号s(n)，其中，功率谱的横坐标为频率，纵坐标为能量。
36.步骤104，对功率谱密度信号进行积分处理，获取功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
37.其中，正向最大频率点表示对应列信号的最大速度点，反向最大频率点则表示对应列信号的最小速度点。在本实施例中，通过对超声图谱中每一列信号n分别对应的功率谱密度信号s(n)进行积分处理，从而得到各个列信号中的正向最大频率点和反向最大频率点。
38.步骤106，连接超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，得到超声图谱对应的波形包络线。
39.具体地，通过将超声图谱中每一列信号分别对应的正向最大频率点依次连接，从
而得到对应的正向包络线，通过将超声图谱中每一列信号分别对应的反向最大频率点依次连接，从而得到对应的反向包络线，由正向包络线和反向包络线则组成的波形则为超声图谱对应的波形包络线。
40.上述实施例中，通过获取超声图谱中每一列信号对应的功率谱密度信号，对功率谱密度信号进行积分处理，以获取功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点，进而连接超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，从而得到超声图谱对应的波形包络线。本实施例通过将超声信号转换为功率谱密度信号，并基于功率谱密度信号同时确定超声图谱对应的正向包络线和反向包络线，从而实现得到由正向包络线和反向包络线组成的波形包络线，且具有计算过程简单、计算效率高等优势。
41.在一个实施例中，如图3所示，上述对功率谱密度信号进行积分处理，获取功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点的步骤，具体可以包括：
42.步骤302，基于功率谱密度信号对应的第一频率区间对功率谱密度信号进行第一积分处理，得到对应的第一积分曲线。
43.其中，第一频率区间是指功率谱密度信号所对应的频率区间。第一积分处理实质是对超声图谱中列信号的灰度从低频到高频的累加过程。第一积分曲线是对功率谱密度信号进行第一积分处理后得到的积分结果。
44.在本实施例，如图4所示，通过对图2中的功率谱密度信号s(n)进行第一积分处理，从而得到对应的第一积分曲线p(n)。具体地，第一积分曲线p(n)是对超声图谱中第n列信号对应的功率谱密度信号s(n)积分的离散数据点曲线，即对功率谱密度信号s(n)随频率的增加进行积分得到的离散数据点曲线。可以理解的是，对于超声图谱中不同的列信号n，其对应的功率谱密度信号s(n)可能存在不同，则进行积分处理后得到的第一积分曲线p(n)也存在不同。
45.步骤304，根据第一积分曲线确定功率谱密度信号的最大能量点。
46.由于功率谱的纵坐标为能量，最大能量点则为图2所示的功率谱中纵坐标最大的点。又由于第一积分曲线p(n)是对图2所示的功率谱密度信号s(n)进行积分处理后得到的，因此，在本实施例中，通过图4所示的第一积分曲线p(n)可以较为方便地确定功率谱密度信号s(n)的最大能量点。
47.具体地，如图4所示，首先确定第一积分曲线p(n)的首尾两个端点，然后通过第一参考线k1连接该两个端点，其中，第一参考线k1为直线。进而获取第一参考线k1与第一积分曲线p(n)之间的第一交叉点m1，将该第一交叉点m1作为功率谱密度信号的最大能量点。
48.步骤306，根据最大能量点确定功率谱密度信号的两个最小能量点。
49.由于最大能量点为功率谱中纵坐标最大的点，则最小能量点为功率谱中纵坐标最小的点。
50.在本实施例中，基于上述步骤确定的最大能量点也即第一交叉点m1将功率谱密度信号s(n)划分为两个信号区间，并分别获取每个信号区间内纵坐标最小的点作为功率谱密度信号的最小能量点。
51.具体地，如图5所示，通过第一交叉点m1在横坐标方向上的位置将功率谱密度信号s(n)进行区域划分，得到s(1)到s(m1)的信号区间以及s(m1)到s(n)的信号区间(其中，s(1)
为坐标原点对应的信号点，s(m1)为m1信号点对应的横坐标，s(n)为远离原点的信号点n对应的横坐标)。并搜索s(1)到s(m1)的信号区间中纵坐标最小的点，得到lowest_p，以及搜索s(m1)到s(n)的信号区间中纵坐标最小的点，得到lowest_r，从而得到两个最小能量点lowest_p和lowest_r。
52.步骤308，基于两个最小能量点对应的第二频率区间对相应区间的功率谱密度信号进行第二积分处理，得到对应的第二积分曲线。
53.在本实施例中，在通过上述步骤确定两个最小能量点后，则以该两个最小能量点分别对应的横坐标所组成的区间为第二频率区间，并对该区间所对应的功率谱密度信号进行第二积分处理，从而得到对应的第二积分曲线。
54.具体地，如图5所示，若最小能量点lowest_p至lowest_r的区间所对应的功率谱密度信号为s(p_r)(s(p_r)是从s(n)中截取的从lowest_p至lowest_r的信号区间对应的信号),则以lowest_p至lowest_r所组成的第二频率区间作为积分区间对s(p_r)进行第二积分处理，从而得到如图6所示的第二积分曲线p(m)。
55.步骤310，根据第二积分曲线确定功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
56.具体地，如图6所示，首先通过获取第二积分曲线p(m)的首部第一端点和尾部第二端点，然后通过第二参考线k2连接第一端点和第二端点，其中，第二参考线k2为直线。第一端点为第二积分曲线p(m)中离原点最近的端点，第二端点则为第二积分曲线p(m)中远离原点的端点。
57.进而确定第二参考线k2与第二积分曲线p(m)之间的第二交叉点m2、正向频率点locate_p和反向频率点locate_r。其中，正向频率点locate_p为第二积分曲线p(m)到第二参考线k2的正向最大距离对应的频率点，反向频率点locate_r为第二积分曲线p(m)到第二参考线k2的反向最大距离对应的频率点。
58.获取第一端点至第二交叉点m2之间以正向频率点locate_p分割的两边信号对应的第一平均值差值，即获取第一端点至正向频率点locate_p之间所有信号点分别对应的频率的第一平均值，并获取正向频率点locate_p至第二交叉点m2之间所有信号点分别对应的频率的第二平均值，并计算第一平均值与第二平均值之间的差值，该差值即为第一平均值差值。当第一平均值差值大于或等于预设阈值时，将该正向频率点locate_p确定为正向最大频率点。其中，预设阈值的范围可以是0至255之间的任意值。
59.同理，获取第二交叉点m2至第二端点之间以反向频率点locate_r分割的两边的第二平均值差值，当第二平均值差值大于或等于预设阈值时，将该反向频率点locate_r确定为反向最大频率点。
60.而当第一平均值差值小于预设阈值时，则向第二交叉点的方向移动正向频率点，即确定locate_(p 1)为正向频率点，并返回执行获取第一端点至第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值的步骤，直到第一平均值差值大于或等于预设阈值时，则对应的正向频率点可以确定为正向最大频率点。
61.同理，当第二平均值差值小于预设阈值时，向第二交叉点的方向移动反向频率点，即确定locate_(r-1)为反向频率点，并返回执行获取第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值的步骤，直到第二平均值差值大于或等于预设
阈值时，则对应的反向频率点可以确定为反向最大频率点。
62.通过对超声图谱中每一列信号对应的功率谱密度信号分别进行如图3所示的处理，即可得到每一列信号对应的正向最大频率点和反向最大频率点。通过将超声图谱中每一列信号分别对应的正向最大频率点依次连接，从而得到对应的正向包络线，通过将超声图谱中每一列信号分别对应的反向最大频率点依次连接，从而得到对应的反向包络线，由正向包络线和反向包络线组成的波形则为超声图谱对应的波形包络线(如图7所示)。
63.上述实施例中，通过对功率谱密度信号进行第一积分处理，从而确定最大能量点和最小能量点，以实现对信号进行滤波去噪的目的，进而对滤波去噪后的信号进行第二积分处理，从而能够较为准确的确定超声图谱中对应列信号的正向最大频率点和反向最大频率点，通过预设阈值的方式还可以灵活调节波形的紧贴程度。
64.应该理解的是，虽然如上所述的各实施例所涉及的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，如上所述的各实施例所涉及的流程图中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
65.基于同样的发明构思，本技术实施例还提供了一种用于实现上述所涉及的超声多普勒包络确定方法的超声多普勒包络确定装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似，故下面所提供的一个或多个超声多普勒包络确定装置实施例中的具体限定可以参见上文中对于超声多普勒包络确定方法的限定，在此不再赘述。
66.在一个实施例中，如图8所示，提供了一种超声多普勒包络确定装置，包括：功率谱密度信号获取模块802、频率点确定模块804和包络线获取模块806，其中：
67.功率谱密度信号获取模块802，用于获取超声图谱中每一列信号对应的功率谱密度信号；
68.频率点确定模块804，用于对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点；
69.包络线获取模块806，用于连接所述超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，得到所述超声图谱对应的包络线。
70.在一个实施例中，频率点确定模块还可以包括：第一积分处理单元，用于基于所述功率谱密度信号对应的第一频率区间对所述功率谱密度信号进行第一积分处理，得到对应的第一积分曲线；最大能量点确定单元，用于根据所述第一积分曲线确定所述功率谱密度信号的最大能量点；最小能量点确定单元，用于根据所述最大能量点确定所述功率谱密度信号的两个最小能量点；第二积分处理单元，用于基于所述两个最小能量点对应的第二频率区间对相应区间的功率谱密度信号进行第二积分处理，得到对应的第二积分曲线；频率点确定单元，用于根据所述第二积分曲线确定所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
71.在一个实施例中，最大能量点确定单元具体还用于：获取所述第一积分曲线的首尾两个端点，通过第一参考线连接所述两个端点，所述第一参考线为直线；获取所述第一参
考线与所述第一积分曲线之间的第一交叉点，将所述第一交叉点作为所述功率谱密度信号的最大能量点。
72.在一个实施例中，最小能量点确定单元具体还用于：根据所述最大能量点将所述功率谱密度信号划分为两个信号区间；分别获取每个信号区间内纵坐标最小的点，作为所述功率谱密度信号的最小能量点。
73.在一个实施例中，频率点确定单元具体还用于：获取所述第二积分曲线的首部第一端点和尾部第二端点，通过第二参考线连接所述第一端点和第二端点，所述第二参考线为直线；确定所述第二参考线与所述第二积分曲线之间的第二交叉点、正向频率点和反向频率点，所述正向频率点为所述第二积分曲线到所述第二参考线的正向最大距离对应的频率点，所述反向频率点为所述第二积分曲线到所述第二参考线的反向最大距离对应的频率点；获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边的第一平均值差值，当所述第一平均值差值大于或等于预设阈值时，将所述正向频率点确定为正向最大频率点；获取所述第二交叉点至第二端点之间以反向频率点分割的两边的第二平均值差值，当所述第二平均值差值大于或等于预设阈值时，将所述反向频率点确定为反向最大频率点。
74.在一个实施例中，频率点确定单元具体还用于：当所述第一平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述正向频率点，并返回执行获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边的第一平均值差值的步骤；当所述第二平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述反向频率点，并返回执行获取所述第二交叉点至第二端点之间以反向频率点分割的两边的第二平均值差值的步骤。
75.上述超声多普勒包络确定装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
76.在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图9所示。该计算机设备包括通过系统总线连接的处理器、存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的通信接口用于与外部的终端进行有线或无线方式的通信，无线方式可通过wifi、移动蜂窝网络、nfc(近场通信)或其他技术实现。该计算机程序被处理器执行时以实现一种超声多普勒包络确定方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
77.本领域技术人员可以理解，图9中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
78.在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
79.获取超声图谱中每一列信号对应的功率谱密度信号；
80.对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点；
81.连接所述超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，得到所述超声图谱对应的波形包络线。
82.在一个实施例中，处理器执行计算机程序时还实现以下步骤：基于所述功率谱密度信号对应的第一频率区间对所述功率谱密度信号进行第一积分处理，得到对应的第一积分曲线；根据所述第一积分曲线确定所述功率谱密度信号的最大能量点；根据所述最大能量点确定所述功率谱密度信号的两个最小能量点；基于所述两个最小能量点对应的第二频率区间对相应区间的功率谱密度信号进行第二积分处理，得到对应的第二积分曲线；根据所述第二积分曲线确定所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
83.在一个实施例中，处理器执行计算机程序时还实现以下步骤：获取所述第一积分曲线的首尾两个端点，通过第一参考线连接所述两个端点，所述第一参考线为直线；获取所述第一参考线与所述第一积分曲线之间的第一交叉点，将所述第一交叉点作为所述功率谱密度信号的最大能量点。
84.在一个实施例中，处理器执行计算机程序时还实现以下步骤：根据所述最大能量点将所述功率谱密度信号划分为两个信号区间；分别获取每个信号区间内纵坐标最小的点，作为所述功率谱密度信号的最小能量点。
85.在一个实施例中，处理器执行计算机程序时还实现以下步骤：获取所述第二积分曲线的首部第一端点和尾部第二端点，通过第二参考线连接所述第一端点和第二端点，所述第二参考线为直线；确定所述第二参考线与所述第二积分曲线之间的第二交叉点、正向频率点和反向频率点，所述正向频率点为所述第二积分曲线到所述第二参考线的正向最大距离对应的频率点，所述反向频率点为所述第二积分曲线到所述第二参考线的反向最大距离对应的频率点；获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值，当所述第一平均值差值大于或等于预设阈值时，将所述正向频率点确定为正向最大频率点；获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值，当所述第二平均值差值大于或等于预设阈值时，将所述反向频率点确定为反向最大频率点。
86.在一个实施例中，处理器执行计算机程序时还实现以下步骤：当所述第一平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述正向频率点，并返回执行获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值的步骤；当所述第二平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述反向频率点，并返回执行获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值的步骤。
87.在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
88.获取超声图谱中每一列信号对应的功率谱密度信号；
89.对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中
的正向最大频率点和反向最大频率点；
90.连接所述超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，得到所述超声图谱对应的波形包络线。
91.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：基于所述功率谱密度信号对应的第一频率区间对所述功率谱密度信号进行第一积分处理，得到对应的第一积分曲线；根据所述第一积分曲线确定所述功率谱密度信号的最大能量点；根据所述最大能量点确定所述功率谱密度信号的两个最小能量点；基于所述两个最小能量点对应的第二频率区间对相应区间的功率谱密度信号进行第二积分处理，得到对应的第二积分曲线；根据所述第二积分曲线确定所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
92.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取所述第一积分曲线的首尾两个端点，通过第一参考线连接所述两个端点，所述第一参考线为直线；获取所述第一参考线与所述第一积分曲线之间的第一交叉点，将所述第一交叉点作为所述功率谱密度信号的最大能量点。
93.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据所述最大能量点将所述功率谱密度信号划分为两个信号区间；分别获取每个信号区间内纵坐标最小的点，作为所述功率谱密度信号的最小能量点。
94.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取所述第二积分曲线的首部第一端点和尾部第二端点，通过第二参考线连接所述第一端点和第二端点，所述第二参考线为直线；确定所述第二参考线与所述第二积分曲线之间的第二交叉点、正向频率点和反向频率点，所述正向频率点为所述第二积分曲线到所述第二参考线的正向最大距离对应的频率点，所述反向频率点为所述第二积分曲线到所述第二参考线的反向最大距离对应的频率点；获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值，当所述第一平均值差值大于或等于预设阈值时，将所述正向频率点确定为正向最大频率点；获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值，当所述第二平均值差值大于或等于预设阈值时，将所述反向频率点确定为反向最大频率点。
95.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：当所述第一平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述正向频率点，并返回执行获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值的步骤；当所述第二平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述反向频率点，并返回执行获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值的步骤。
96.在一个实施例中，提供了一种计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：
97.获取超声图谱中每一列信号对应的功率谱密度信号；
98.对所述功率谱密度信号进行积分处理，获取所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点；
99.连接所述超声图谱中每一列信号分别对应的正向最大频率点和反向最大频率点，
得到所述超声图谱对应的波形包络线。
100.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：基于所述功率谱密度信号对应的第一频率区间对所述功率谱密度信号进行第一积分处理，得到对应的第一积分曲线；根据所述第一积分曲线确定所述功率谱密度信号的最大能量点；根据所述最大能量点确定所述功率谱密度信号的两个最小能量点；基于所述两个最小能量点对应的第二频率区间对相应区间的功率谱密度信号进行第二积分处理，得到对应的第二积分曲线；根据所述第二积分曲线确定所述功率谱密度信号对应的列信号中的正向最大频率点和反向最大频率点。
101.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取所述第一积分曲线的首尾两个端点，通过第一参考线连接所述两个端点，所述第一参考线为直线；获取所述第一参考线与所述第一积分曲线之间的第一交叉点，将所述第一交叉点作为所述功率谱密度信号的最大能量点。
102.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：根据所述最大能量点将所述功率谱密度信号划分为两个信号区间；分别获取每个信号区间内纵坐标最小的点，作为所述功率谱密度信号的最小能量点。
103.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取所述第二积分曲线的首部第一端点和尾部第二端点，通过第二参考线连接所述第一端点和第二端点，所述第二参考线为直线；确定所述第二参考线与所述第二积分曲线之间的第二交叉点、正向频率点和反向频率点，所述正向频率点为所述第二积分曲线到所述第二参考线的正向最大距离对应的频率点，所述反向频率点为所述第二积分曲线到所述第二参考线的反向最大距离对应的频率点；获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值，当所述第一平均值差值大于或等于预设阈值时，将所述正向频率点确定为正向最大频率点；获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值，当所述第二平均值差值大于或等于预设阈值时，将所述反向频率点确定为反向最大频率点。
104.在一个实施例中，计算机程序被处理器执行时还实现以下步骤：当所述第一平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述正向频率点，并返回执行获取所述第一端点至所述第二交叉点之间以正向频率点分割的两边信号对应的第一平均值差值的步骤；当所述第二平均值差值小于预设阈值时，向所述第二交叉点的方向移动所述反向频率点，并返回执行获取所述第二交叉点至第二端点之间以反向频率点分割的两边信号对应的第二平均值差值的步骤。
105.本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(read-only memory，rom)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(reram)、磁变存储器(magnetoresistive random access memory，mram)、铁电存储器(ferroelectric random access memory，fram)、相变存储器(phase change memory，
pcm)、石墨烯存储器等。易失性存储器可包括随机存取存储器(random access memory，ram)或外部高速缓冲存储器等。作为说明而非局限，ram可以是多种形式，比如静态随机存取存储器(static random access memory，sram)或动态随机存取存储器(dynamic random access memory，dram)等。本技术所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。本技术所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。
106.应当理解，以上实施例均为示例性的，不用于包含权利要求所包含的所有可能的实施方式。在不脱离本公开的范围的情况下，还可以在以上实施例的基础上做出各种变形和改变。同样的，也可以对以上实施例的各个技术特征进行任意组合，以形成可能没有被明确描述的本发明的另外的实施例。因此，上述实施例仅表达了本发明的几种实施方式，不对本发明专利的保护范围进行限制。