一种蠕动波的参数测量方法及超声测量系统与流程

1.本发明涉及蠕动波测量领域，具体而言涉及一种蠕动波的参数测量方法及超声测量系统。


背景技术：

2.子宫内膜容受性，指的是子宫内膜处于一种允许囊胚定位、黏附、侵入并使内膜间质发生改变从而导致胚胎着床的状态，即子宫内膜接受受精卵的能力。对子宫内膜容受性的正确评价在选择种植时机、评估妊娠率等方面具有重要临床意义，是当今生殖评估标准体系中重要的一环。
3.子宫内膜蠕动波，指的是子宫肌层收缩带动子宫内膜蠕动而产生的机械波。蠕动波的频次、方向、强度等随月经周期变化而变化，进而辅助精子运输与胚胎着床，是评判子宫内膜容受性的重要指标之一。同时，子宫蠕动规律受子宫疾病的影响，因此对蠕动波的研究也有辅助诊断子宫病变的潜在价值。
4.然而目前没有辅助蠕动波临床诊断的客观手段，只能依靠医生肉眼重复观察所采集的超声影片作出的主观判断，耗时长、友好度低，且操作者重复性差，不同医生判断的蠕动频次、方向不一定一致。因此，客观辅助工具的缺乏是对蠕动波进行进一步研究与应用的主要阻碍。


技术实现要素：

5.在发明内容部分中引入了一系列简化形式的概念，这将在具体实施方式部分中进一步详细说明。本发明的发明内容部分并不意味着要试图限定出所要求保护的技术方案的关键特征和必要技术特征，更不意味着试图确定所要求保护的技术方案的保护范围。
6.针对现有技术的不足，本发明实施例第一方面提供了一种蠕动波的参数测量方法，所述方法包括：
7.向被测对象的子宫内膜发射第一超声波，接收所述子宫内膜返回的超声回波，以得到第一超声回波信号；
8.对所述第一超声回波信号进行处理，以得到所述子宫内膜中目标区域内随时间变化的蠕动参数；
9.基于所述目标区域内随时间变化的蠕动参数得到蠕动波在所述目标区域内传递的传递时间；
10.输出所述传递时间。
11.在一个实施例中，所述基于所述目标区域内随时间变化的蠕动参数得到蠕动波在所述目标区域内传递的传递时间，包括：
12.根据所述目标区域内不同位置处随时间变化的蠕动参数生成蠕动波时空分布图，所述蠕动波时空分布图表示所述蠕动参数随时间和空间的变化；
13.基于所述蠕动波时空分布图确定所述传递时间。
14.在一个实施例中，所述方法还包括：
15.向所述被测对象的子宫内膜发射第二超声波；
16.接收所述子宫内膜返回的第二超声回波，以得到第二超声回波信号；
17.对所述第二超声回波信号进行处理，以得到所述子宫内膜的超声图像；
18.根据所述超声图像确定所述目标区域。
19.在一个实施例中，所述基于所述目标区域内随时间变化的蠕动参数得到蠕动波在所述目标区域内传递的传递时间，包括：
20.显示所述蠕动波时空分布图；
21.获得对所述蠕动波时空分布图上蠕动波传递到所述目标区域内不同位置处的时间点进行的标注；
22.根据所述标注对应的时间点确定蠕动波在所述不同位置之间传递的传递时间。
23.在一个实施例中，所述获得对所述蠕动波时空分布图上蠕动波传递到所述目标区域内不同位置处的时间点进行的标注，包括：
24.接收对所述蠕动波时空分布图上蠕动波传递到所述不同位置处的时间点进行的点选操作，根据所述点选操作确定所述标注的位置；
25.或者，在所述蠕动波时空分布图上显示可调节光标，接收对所述可调节光标进行的调节操作，根据所述调节操作确定所述标注的位置。
26.在一个实施例中，所述基于所述蠕动波时空分布图确定所述传递时间，包括：
27.基于所述蠕动波时空分布图分别得到所述目标区域内至少两个位置处蠕动参数随时间变化的至少两条曲线；
28.提取所述至少两条曲线上蠕动波同一波动段的对应时间点，根据所述对应时间点之间的时间间隔确定蠕动波在所述至少两个位置之间传递的传递时间。
29.在一个实施例中，所述对应时间点包括所述至少两条曲线上同一波动段的峰值、同一波动段的起始点或同一波动段的终点所对应的时间点。
30.在一个实施例中，所述基于所述蠕动波时空分布图确定所述传递时间，包括：
31.基于所述蠕动波时空分布图得到预定位置处蠕动参数随时间变化的蠕动曲线；
32.提取所述蠕动曲线上相邻波动段上的对应时间点，根据所述对应时间点之间的间隔时间确定所述传递时间。
33.在一个实施例中，所述基于所述目标区域内随时间变化的蠕动参数得到蠕动波在所述目标区域内传递的传递时间，包括：
34.基于机器学习算法根据所述蠕动参数自动分析蠕动波的传递时间。
35.在一个实施例中，所述方法还包括：
36.获得多次蠕动波在所述目标区域内传递的多个传递时间；
37.确定多个所述传递时间的平均传递时间，输出所述传递时间还包括输出所述平均传递时间。
38.在一个实施例中，所述蠕动参数包括以下至少一项：蠕动速度、组织位移、组织应变。
39.在一个实施例中，所述方法还包括：基于所述目标区域内随时间变化的蠕动参数得到预定时间内所述目标区域内的蠕动次数。
40.本发明实施例第二方面提供了一种超声测量系统，所述超声测量系统包括超声探头、发射电路、接收电路、存储器、处理器和显示器，所述存储器上存储有由所述处理器运行的计算机程序，所述计算机程序在被所述处理器运行时执行上述蠕动波的参数测量方法的步骤。
41.本技术实施例的蠕动波的参数测量方法及超声测量系统将蠕动波的传递时间作为一项新的蠕动波相关参数进行量化并输出，为用户提供了蠕动波的客观测量工具。
附图说明
42.为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
43.图1示出根据本发明一个实施例的超声测量系统的示意性框图；
44.图2示出根据本发明一个实施例的蠕动波的参数测量方法的示意性流程图；
45.图3示出根据本发明一个实施例的基于蠕动波时空分布图手动测量传递时间的示意图；
46.图4示出根据本发明另一个实施例的基于蠕动波时空分布图手动测量传递时间的示意图；
47.图5示出根据本发明一个实施例的基于蠕动波时空分布图自动测量传递时间的示意图；
48.图6示出根据本发明另一实施例的基于蠕动波时空分布图自动测量传递时间的示意图；
49.图7示出根据本发明一个实施例的显示界面的示意图。
具体实施方式
50.为了使得本发明的目的、技术方案和优点更为明显，下面将参照附图详细描述根据本发明的示例实施例。显然，所描述的实施例仅仅是本发明的一部分实施例，而不是本发明的全部实施例，应理解，本发明不受这里描述的示例实施例的限制。基于本发明中描述的本发明实施例，本领域技术人员在没有付出创造性劳动的情况下所得到的所有其它实施例都应落入本发明的保护范围之内。
51.在下文的描述中，给出了大量具体的细节以便提供对本发明更为彻底的理解。然而，对于本领域技术人员而言显而易见的是，本发明可以无需一个或多个这些细节而得以实施。在其他的例子中，为了避免与本发明发生混淆，对于本领域公知的一些技术特征未进行描述。
52.应当理解的是，本发明能够以不同形式实施，而不应当解释为局限于这里提出的实施例。相反地，提供这些实施例将使公开彻底和完全，并且将本发明的范围完全地传递给本领域技术人员。
53.在此使用的术语的目的仅在于描述具体实施例并且不作为本发明的限制。在此使用时，单数形式的“一”、“一个”和“所述/该”也意图包括复数形式，除非上下文清楚指出另
外的方式。还应明白术语“组成”和/或“包括”，当在该说明书中使用时，确定所述特征、整数、步骤、操作、元件和/或部件的存在，但不排除一个或更多其它的特征、整数、步骤、操作、元件、部件和/或组的存在或添加。在此使用时，术语“和/或”包括相关所列项目的任何及所有组合。
54.为了彻底理解本发明，将在下列的描述中提出详细的结构，以便阐释本发明提出的技术方案。本发明的可选实施例详细描述如下，然而除了这些详细描述外，本发明还可以具有其他实施方式。
55.下面，首先参考图1描述根据本技术一个实施例的超声测量系统，该超声测量系统可以用于实现本技术实施例的蠕动波的参数测量方法。图1示出了根据本技术实施例的超声测量系统100的示意性结构框图。
56.如图1所示，超声测量系统100包括超声探头110、发射电路112、接收电路114、处理器116和显示器118。进一步地，超声测量系统还可以包括发射/接收选择开关120和波束合成电路122，发射电路112和接收电路114可以通过发射/接收选择开关120与超声探头110连接。
57.超声探头110包括多个换能器阵元，多个换能器阵元可以排列成一排构成线阵，或排布成二维矩阵构成面阵，多个换能器阵元也可以构成凸阵列。换能器用于根据激励电信号发射超声波，或将接收的超声波转换为电信号，因此每个阵元可用于实现电脉冲信号和超声波的相互转换，从而实现向被测对象的目标区域的组织发射超声波、也可用于接收经组织反射回的超声波回波。在进行超声检测时，可通过发射序列和接收序列控制哪些换能器阵元用于发射超声波，哪些换能器阵元用于接收超声波，或者控制换能器阵元分时隙用于发射超声波或接收超声波的回波。参与超声波发射的换能器阵元可以同时被电信号激励，从而同时发射超声波；或者，参与超声波束发射的换能器阵元也可以被具有一定时间间隔的若干电信号激励，从而持续发射具有一定时间间隔的超声波。
58.在超声成像过程中，发射电路112将经过延迟聚焦的发射脉冲通过发射/接收选择开关120发送到超声探头110。超声探头110受发射脉冲的激励而向被测对象的目标区域的组织发射超声波束，经一定延时后接收从目标区域的组织反射回来的带有组织信息的超声回波，并将此超声回波重新转换为电信号。接收电路114接收超声探头110转换生成的电信号，获得超声回波信号，并将这些超声回波信号送入波束合成电路122，波束合成电路122对超声回波数据进行聚焦延时、加权和通道求和等处理，然后送入处理器116。处理器116对超声回波信号进行信号检测、信号增强、数据转换、对数压缩等处理形成超声图像。处理器116得到的超声图像可以在显示器118上显示，也可以存储于存储器124中。
59.可选地，处理器116可以实现为软件、硬件、固件或其任意组合，并且可以使用单个或多个专用集成电路(application specific integrated circuit，asic)、单个或多个通用集成电路、单个或多个微处理器、单个或多个可编程逻辑器件、或者前述电路和/或器件的任意组合、或者其他适合的电路或器件。并且，处理器116可以控制所述超声测量系统100中的其它组件以执行本说明书中的各个实施例中的方法的相应步骤。
60.显示器118与处理器116连接，显示器118可以为触摸显示屏、液晶显示屏等；或者，显示器118可以为独立于超声测量系统100之外的液晶显示器、电视机等独立显示器；或者，显示器118可以是智能手机、平板电脑等电子设备的显示屏，等等。其中，显示器118的数量
可以为一个或多个。例如，显示器118可以包括主屏和触摸屏，主屏主要用于显示超声图像，触摸屏主要用于人机交互。
61.显示器118可以显示处理器116得到的超声图像。此外，显示器118在显示超声图像的同时还可以提供给用户进行人机交互的图形界面，在图形界面上设置一个或多个被控对象，提供给用户利用人机交互装置输入操作指令来控制这些被控对象，从而执行相应的控制操作。例如，在图形界面上显示图标，利用人机交互装置可以对该图标进行操作，用来执行特定的功能，例如在超声图像上绘制出感兴趣区域框等。
62.可选地，超声测量系统100还可以包括显示器118之外的其他人机交互装置，其与处理器116连接，例如，处理器116可以通过外部输入/输出端口与人机交互装置连接，外部输入/输出端口可以是无线通信模块，也可以是有线通信模块，或者两者的组合。外部输入/输出端口也可基于usb、如can等总线协议、和/或有线网络协议等来实现。
63.其中，人机交互装置可以包括输入设备，用于检测用户的输入信息，该输入信息例如可以是对超声波发射/接收时序的控制指令，可以是在超声图像上绘制出点、线或框等的操作输入指令，或者还可以包括其他指令类型。输入设备可以包括键盘、鼠标、滚轮、轨迹球、移动式输入设备(例如带触摸显示屏的移动设备、手机等等)、多功能旋钮等等其中之一或者多个的结合。人机交互装置还可以包括诸如打印机之类的输出设备。
64.超声测量系统100还可以包括存储器124，用于存储处理器执行的指令、存储接收到的超声回波、存储超声图像，等等。存储器可以为闪存卡、固态存储器、硬盘等。其可以为易失性存储器和/或非易失性存储器，为可移除存储器和/或不可移除存储器等。
65.应理解，图1所示的超声测量系统100所包括的部件只是示意性的，其可以包括更多或更少的部件。本技术对此不限定。
66.下面，将参考图2描述根据本技术一个实施例的蠕动波的参数测量方法。图2是本技术实施例的蠕动波的参数测量方法200的一个示意性流程图。
67.如图2所示，本技术实施例的蠕动波的参数测量方法200包括以下步骤：
68.在步骤s210，向被测对象的子宫内膜发射第一超声波，接收所述子宫内膜返回的超声回波，以得到第一超声回波信号；
69.在步骤s220，对所述第一超声回波信号进行处理，以得到所述子宫内膜中目标区域内随时间变化的蠕动参数；
70.在步骤s230，基于所述目标区域内随时间变化的蠕动参数得到蠕动波在所述目标区域内传递的传递时间；
71.在步骤s240，输出所述传递时间。
72.研究表明，除了蠕动波的频次、方向、强度会影响妊娠率以外，每一次蠕动波的传递时间(即每次蠕动波的持续时间)也能够作为评价子宫内膜容受性的重要指标。本技术实施例的蠕动波的参数测量方法200将蠕动波在目标区域内传递的传递时间作为一项新的蠕动波相关参数量化并输出，为用户提供了蠕动波的客观测量工具，有助于对蠕动波进行进一步的定量的临床研究与诊断评估，为日后完善蠕动波评判体系打下基础。
73.参照图1，在步骤s210中，发射电路112将发射脉冲通过发射/接收选择开关120发送到超声探头110，以激励超声探头110向被测对象的子宫内膜发射第一超声波。第一超声波可以针对整个子宫内膜，也可以只针对子宫内膜中的目标区域。目标区域可以是通过超
声图像(包括但不限于b型超声图像、c型超声图像)确定的，可以是通过蠕动波图像确定的，也可以是直接根据预设的扫描策略确定的。经一定延时后，在步骤s220中，接收电路114通过发射/接收选择开关120控制超声探头接收第一超声波的回波，得到第一超声回波信号，并将该第一超声回波信号发送给波束合成电路122，波束合成电路122对第一超声回波信号进行聚焦延时、加权和通道求和等处理，然后将该波束合成后的第一超声回波信号送入处理器116进行处理，以得到蠕动参数。
74.在本技术实施例中，处理器116可以对波束合成处理后任一环节的第一超声回波信号进行处理以得到蠕动参数，也可以基于第一超声回波信号生成超声图像，并对超声图像进行相关处理，以得到蠕动参数。上文中的第一超声波为蠕动波的检测序列，其可以与超声图像的发射扫描序列共用，也可以采取完全不同扫描序列。当扫描序列独立时，其发射接收参数(例如频率、聚焦方向、发射间隔、发射位置等)也可以独立设置，例如可以使相邻两帧回波的扫描间隔更短，从而获得时间上更高精度的检测。
75.对于空间中的某个目标位置，当持续一段时间发射超声波，并接收超声回波时，如果该目标位置是运动的，则不同时刻所得到的超声回波将发生变化，基于相关对比可以检测出各时刻下超声回波的变化量或变化速度，即位移检测或者速度检测。基于此，对于蠕动中的子宫内膜，可以采用相关对比获得其蠕动参数，蠕动参数可以包括以下至少一项：蠕动速度、组织位移、组织应变。作为一种实现方式，可以通过位移检测方法得到感兴趣区域中的各个点的位移或速度；作为另一种实现方式，可以通过位移检测方法得到感兴趣区域中的各个点的位移，再通过计算位移在时间上的梯度，得到各个点的速度。
76.基于位移检测确定蠕动参数的具体方法可以有多种。例如，可以采用基于块匹配的方式，对于某一时刻的某位置处的超声回波信号，在另一个时刻的超声回波信号的不同位置进行搜索，找出与之互相关最大的位置，其与原位置的差异即作为这两个时刻下的位移量，进一步结合两个时刻之间的时间差可以得到蠕动速度、蠕动加速度等蠕动参数。又例如，可以采用基于超声多普勒效应的方式，以类似常规血流成像的原理检测出某位置处的组织在各个时刻的运动速度。或者，也可以采用其他基于信号自相关或者互相关的位移检测方法，本技术实施例对此不做限制。
77.在一个实施例中，当根据超声图像确定目标区域时，则该蠕动波的参数测量方法200还包括：控制超声探头向被测对象的子宫内膜发射第二超声波；接收被测对象的子宫内膜返回的第二超声回波，以得到第二超声回波信号；对第二超声回波信号进行处理，以得到子宫内膜的超声图像，包括但不限于b型超声图像或c型超声图像等。之后，在超声图像中确定感兴趣区域，感兴趣区域与子宫内膜的目标区域相对应。在一些实施例中，可以根据第二回波信号获得子宫内膜的多帧超声图像；在获得子宫内膜的多帧超声图像后，可以根据需要，在多帧超声图像中的第一帧超声图像上确定感兴趣区域，也可以在中间帧的超声图像上确定感兴趣区域，或者在最后一帧超声图像上确定感兴趣区域。示例性地，感兴趣区域可以是点、线、框等，具体可以包括直线、曲线、离散的点、连续的点或任意形状的框，本技术实施例对感兴趣区域的形状不做限制。
78.示例性地，在超声图像中确定感兴趣区域可以实现为自动或手动方式。当采用手动方式确定感兴趣区域时，可以显示超声图像，并响应于用户在超声图像上的选择操作确定感兴趣区域，例如，用户可以通过鼠标等输入设备在超声图像中选择感兴趣区域。
79.当自动确定感兴趣区域时，则超声测量系统可以自动识别超声图像中的子宫内膜区域，并在识别到的子宫内膜区域内按预设规则自动选择感兴趣区域。例如，可以在子宫内膜区域内部选择一个或多个兴趣点以作为感兴趣区域；或者，可以在子宫内膜区域的某一方向上选择线段作为感兴趣区域，例如所选的线段可以是从宫颈端延伸到宫底端的线段，等等。当然，上述自动选择感兴趣区域的方式仅作为示例，本技术实施例也可基于其他预设条件自动选择感兴趣区域，本技术实施例对此不作限制。示例性地，当根据超声图像自动识别感兴趣区域后，还可以在超声图像上显示感兴趣区域所在的位置，并且可以根据用户输入对感兴趣区域所在的位置进行调整。
80.在步骤s220中，当基于超声图像得到蠕动参数时，可以针对超声图像中的每个像素点得到对应的蠕动参数，也可以只针对超声图像中子宫内膜区域的各像素点得到对应的蠕动参数。之后，可以从中提取感兴趣区域对应的蠕动参数。或者，也可以只针对超声图像中感兴趣区域对应的像素点得到对应的蠕动参数。对于多帧超声图像，可以得到每一帧超声图像对应的时刻上各像素点的蠕动参数，从而得到目标区域内各位置处随时间变化的蠕动参数。
81.在步骤s230，基于目标区域内随时间变化的蠕动参数得到蠕动波在目标区域内传递的传递时间，该传递时间能够作为评价子宫内膜容受性的重要指标，有助于对蠕动波进行进一步的定量研究与评估。示例性地，当目标区域为点时，该传递时间可以是两点之间的传递时间；当目标区域为线段时，该传递时间可以是从线段一端传递到另一端的传递时间；当目标区域为框时，该传递时间可以是从框的一边传递到另一边的传递时间。
82.在一个实施例中，可以根据目标区域内不同位置处随时间变化的蠕动参数生成蠕动波时空分布图，并基于蠕动波时空分布图确定蠕动波在目标区域内传递的传递时间。
83.具体地，蠕动波时空分布图表示蠕动参数随时间和空间的变化。生成蠕动波时空分布图可以包括：建立蠕动波时空分布图坐标系，其中，该坐标系包括第一坐标轴和第二坐标轴，第一坐标轴表示时间，第二坐标轴表示空间位置；根据蠕动参数对应的时间和空间位置，将蠕动参数显示在时空分布图坐标系中。示例性地，在蠕动波时空分布图中可以以不同的颜色或灰度表示蠕动参数的大小或方向。由于蠕动波时空分布图包含有蠕动波的时间信息和空间信息，因而根据蠕动波时空分布图可以进一步得到蠕动波在目标区域内传递的传递时间。
84.为了便于理解，图3示出了一种示例性的蠕动波时空分布图，该蠕动波时空分布图所对应的感兴趣区域为线形。在图3所示的蠕动波时空分布图中，横轴表示时间，纵轴表示空间位置，图中深色平行四边形表示正向的蠕动速度，浅色平行四边形表示反向的速度蠕动，可以理解的是，图3的蠕动波时空分布图整体上表示蠕动波从下方对应位置传递到上方对应位置，即在t1时刻蠕动波传递到目标区域的一端，t2时刻蠕动波传递到目标区域的另一端，t2-t1即为蠕动波在目标区域内传递的传递时间。
85.示例性地，若在超声图像中确定了多个感兴趣区域，则可以分别生成与感兴趣区域对应的多个蠕动波时空分布图，并基于每个蠕动波时空分布图得到对应的传递时间，以进行对比分析。例如，子宫内膜上下边界的蠕动表现可能是不同的，因而可以分别获得子宫内膜上下边界的蠕动波时空分布图，并基于各蠕动波时空分布图分别确定传递时间，从而更直观地观察各位置处蠕动波的差异性。
86.根据蠕动波时空分布图确定蠕动波传递时间的方式可以实现为系统自动确定或根据接收到的用户输入确定，即手动确定。其中，手动确定方式具体包括：显示该蠕动波时空分布图；获得对蠕动波时空分布图上蠕动波传递到目标区域内不同位置处的时间点进行的标注；根据该标注对应的时间点确定蠕动波在不同位置之间传递的传递时间。
87.在一个示例中，获得对蠕动波时空分布图上蠕动波传递到目标区域内不同位置处的时间点进行的标注，包括：接收对蠕动波时空分布图上蠕动波传递到不同位置处的时间点进行的点选操作，根据该点选操作确定所述标注的位置。继续参照图3，图3所示的蠕动波时空分布图的纵轴底部对应目标区域的起始位置，顶部对应于目标区域的结束位置，若接收到用户在蠕动波时空分布图上对蠕动波传递到起始位置的开始时间点t1和传递到结束位置的结束时间点t2进行的标注，则可以确定蠕动波在目标区域内传递的传递时间t＝t2-t1。
88.在另一个示例中，可以在蠕动波时空分布图上显示可调节光标，接收对可调节光标进行的调节操作，并根据接收到的调节操作确定用户选定的标注所在的位置。参见图4，图4左侧显示了可调节光标的初始位置，用户可以调节可调节光标，对其进行平移、宽度调节等操作，使其分别对应蠕动波开始点t1和结束时间点t2，即根据接收到的用户输入可以确定蠕动波的传递时间t＝
△
t＝t2-t1。
89.当然，获得对蠕动波时空分布图上蠕动波传递到目标区域内不同位置处的时间点进行的标注的具体方式不限于以上两种，例如用户也可以在蠕动波时空分布图上划线、绘制选框等等，只要使系统能够确定用户对蠕动波时空分布图进行的标注即可。
90.在其他实施例中，超声测量系统也可以根据蠕动波时空分布图自动测量传递时间，以简化用户操作。
91.具体地，在一个示例中，基于所述蠕动波时空分布图确定传递时间包括：基于蠕动波时空分布图分别得到目标区域内至少两个位置处蠕动参数随时间变化的至少两条曲线；提取至少两条曲线上蠕动波同一波动段的对应时间点，根据对应时间点之间的时间间隔确定蠕动波在至少两个位置之间传递的传递时间。其中，对应时间点可以是至少两条曲线上同一波动段波峰、波谷、同一波动段的起始点、同一波动段的终点、同一波动段与坐标轴的交点等特征点。示例性地，上述目标区域内至少两个位置至少包括目标区域的两端位置，则可以得到蠕动波进入目标区域到离开目标区域的传递时间。例如，若目标区域为线段，则目标区域内至少两个位置至少包括线段的两个端点位置。
92.如图5所示，在得到图5左侧的蠕动波时空分布图以后，对其纵向逐行扫描，则每一行得到的蠕动曲线为空间中对应位置处随时间变化的蠕动曲线。提取每一个蠕动曲线上的同一个波动段的对应时间点，例如，第一个波峰的最早时间点为t1，最晚时间点为t2，则蠕动波第一个波动段的传递时间为t＝t2-t1。
93.在另一个实施例中，基于蠕动波时空分布图确定传递时间，包括：基于蠕动波时空分布图得到预定位置处蠕动参数随时间变化的蠕动曲线；提取蠕动曲线上相邻波动段上的对应时间点，根据对应时间点之间的间隔时间确定蠕动波的传递时间。其中，可以提取蠕动波时空分布图中预设空间位置处随时间变化的蠕动参数，绘制该位置处蠕动参数随时间变化的蠕动曲线，并根据同一条蠕动曲线得到蠕动波的传递时间。
94.具体地，当子宫内膜连续蠕动时，后一次蠕动开始时，前一次蠕动到达目标区域末
端。例如，若目标区域为自宫颈端到宫底端的连线时，若子宫内膜连续蠕动，则当第一个波动段传递到宫底端时，宫颈端的第二个波动段开始。因此，可以根据相邻波动段上的对应时间点之间的间隔时间确定蠕动波的传递时间。
95.示例性地，参见图6，当子宫内膜连续蠕动时，子宫内膜上的点的蠕动(例如速度、位移、应变等)曲线为类波动式曲线，曲线上相邻两个波峰之间的时间差即为两次蠕动之间的时间间隔(即图6中的t1、t2和t3)，同时也是子宫内膜连续蠕动时的蠕动波传递时间。
96.在本技术实施例中，除了根据蠕动波时空分布图确定蠕动波的传递时间以外，也可以采用其他方式确定蠕动波在目标区域内传递的传递时间。例如，可以基于机器学习算法根据蠕动参数自动分析蠕动波的传递时间。具体地，机器学习算法可以通过分析子宫内膜的形态变化规律、或通过观察得到蠕动开始时刻与同一次蠕动的结束时刻等，从而得到蠕动波的传递时间。
97.在一些实施例中，还可以获得多次蠕动波在目标区域内传递的多个传递时间，并确定多个传递时间的平均传递时间以及输出该平均传递时间，平均传递时间能更好地反应蠕动波传递时间的一般水平。例如，参照图6，可以确定预设时间内每相邻两个波动段上对应时间点之间的时间间隔(即图6中的t1、t2和t3)，计算其平均值，从而得到多个波动段的平均时间间隔，即连续蠕动时每次蠕动波的平均传递时间。当然，除了基于预定位置处蠕动参数随时间变化的蠕动曲线确定蠕动波的传递时间以外，采用其他方式确定蠕动波的传递时间时，也可以确定多次蠕动波的平均传递时间，并输出该平均传递时间。
98.在一个实施例中，还可以基于目标区域内随时间变化的蠕动参数得到预定时间内目标区域内的蠕动次数。例如，可以根据蠕动参数得到蠕动波时空分布图，并基于蠕动波时空分布图得到预定时间内目标区域内的蠕动次数。示例性地，可以由用户对蠕动波传递到同一空间位置处的时间点进行标注，预定时间内用户标注的个数即为预定时间内目标区域内的蠕动次数。或者，也可以由超声测量系统对蠕动波时空分布图进行逐列扫描以得到蠕动曲线，并根据蠕动曲线自动分析得到预定时间内目标区域内的蠕动次数，例如，可以将蠕动曲线上波峰或波谷等出现的次数确定为蠕动次数。
99.在一个实施例中，可以将蠕动波的传递时间、平均传递时间、蠕动次数等参数与超声图像和蠕动波时空分布图显示在同一显示界面上，有助于使用户更好地识别定位上述参数和蠕动波时空分布图所对应的解剖位置。图7示出了一种示例性的显示界面，该显示界面左上方显示的是子宫内膜的超声图像，在超声图像中显示有折线形的感兴趣区域；左下方显示的是蠕动波时空分布图，该蠕动波时空分布图中显示有用户对蠕动波开始时间点t1和结束时间点t2进行的标注；显示界面右侧显示的是基于蠕动波时空分布图得到的参数，具体包括m次蠕动波的传递时间n1s、n2s、n3s
……
nms，m次蠕动波的平均传递时间ns，以及蠕动方向、蠕动频谱、最大幅度、平均幅度等蠕动波相关参数。
100.本技术实施例的蠕动波的参数测量方法将蠕动波的传递时间作为一项新的蠕动波相关参数进行量化并输出，为用户提供了蠕动波的客观测量工具。
101.现在重新参照图1，本技术实施例还提供了一种超声测量系统100，超声测量系统100可以用于实现上述蠕动波的参数测量方法200。超声测量系统100可以包括超声探头110、发射电路112、接收电路114、处理器116、显示器118以及存储器124等部件，各个部件的相关描述可以参照上文。以下仅对超声测量系统100的主要功能进行描述，而省略以上已经
描述过的细节内容。
102.其中，发射电路112用于激励超声探头110向被测对象的子宫内膜发射第一超声波；接收电路114用于控制超声探头110接收被测对象的子宫内膜返回的超声回波，以获得第一超声回波信号；处理器114用于对第一超声回波信号进行处理，以得到被测对象子宫内膜中目标区域内随时间变化的蠕动参数；以及基于目标区域内随时间变化的蠕动参数得到蠕动波在目标区域内传递的传递时间；处理器114还用于控制输出设备输出所述传递时间，例如，可以在显示器118上显示该传递时间。
103.超声测量系统100和由超声测量系统100实现的蠕动波的参数测量方法200的其他具体细节可以参照上文的相关描述，在此不做赘述。
104.另外，本发明实施例还提供了一种计算机存储介质，其上存储有计算机程序。在计算机可读存储介质上可以存储一个或多个计算机程序指令，处理器可以运行存储装置存储的程序指令，以实现本文的本发明实施例中(由处理器实现)的功能以及/或者其它期望的功能，例如以执行根据本发明实施例的蠕动波的参数测量方法的相应步骤，在计算机可读存储介质中还可以存储各种应用程序和各种数据，例如应用程序使用和/或产生的各种数据等。
105.例如，计算机存储介质例如可以包括存储卡、平板电脑的存储部件、个人计算机的硬盘、只读存储器(rom)、可擦除可编程只读存储器(eprom)、便携式紧致盘只读存储器(cd-rom)、usb存储器、或者上述存储介质的任意组合。
106.综上，本技术实施例的蠕动波的参数测量方法及超声测量系统将蠕动波的传递时间作为一项新的蠕动波相关参数进行量化并输出，为用户提供了蠕动波的客观测量工具。
107.尽管这里已经参考附图描述了示例实施例，应理解上述示例实施例仅仅是示例性的，并且不意图将本发明的范围限制于此。本领域普通技术人员可以在其中进行各种改变和修改，而不偏离本发明的范围和精神。所有这些改变和修改意在被包括在所附权利要求所要求的本发明的范围之内。
108.本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
109.在本技术所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。例如，以上所描述的设备实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个设备，或一些特征可以忽略，或不执行。
110.在此处所提供的说明书中，说明了大量具体细节。然而，能够理解，本发明的实施例可以在没有这些具体细节的情况下实践。在一些实例中，并未详细示出公知的方法、结构和技术，以便不模糊对本说明书的理解。
111.类似地，应当理解，为了精简本发明并帮助理解各个发明方面中的一个或多个，在对本发明的示例性实施例的描述中，本发明的各个特征有时被一起分组到单个实施例、图、或者对其的描述中。然而，并不应将该本发明的方法解释成反映如下意图：即所要求保护的
本发明要求比在每个权利要求中所明确记载的特征更多的特征。更确切地说，如相应的权利要求书所反映的那样，其发明点在于可以用少于某个公开的单个实施例的所有特征的特征来解决相应的技术问题。因此，遵循具体实施方式的权利要求书由此明确地并入该具体实施方式，其中每个权利要求本身都作为本发明的单独实施例。
112.本领域的技术人员可以理解，除了特征之间相互排斥之外，可以采用任何组合对本说明书(包括伴随的权利要求、摘要和附图)中公开的所有特征以及如此公开的任何方法或者设备的所有过程或单元进行组合。除非另外明确陈述，本说明书(包括伴随的权利要求、摘要和附图)中公开的每个特征可以由提供相同、等同或相似目的替代特征来代替。
113.此外，本领域的技术人员能够理解，尽管在此所述的一些实施例包括其它实施例中所包括的某些特征而不是其它特征，但是不同实施例的特征的组合意味着处于本发明的范围之内并且形成不同的实施例。例如，在权利要求书中，所要求保护的实施例的任意之一都可以以任意的组合方式来使用。
114.本发明的各个部件实施例可以以硬件实现，或者以在一个或者多个处理器上运行的软件模块实现，或者以它们的组合实现。本领域的技术人员应当理解，可以在实践中使用微处理器或者数字信号处理器(dsp)来实现根据本发明实施例的一些模块的一些或者全部功能。本发明还可以实现为用于执行这里所描述的方法的一部分或者全部的装置程序(例如，计算机程序和计算机程序产品)。这样的实现本发明的程序可以存储在计算机可读介质上，或者可以具有一个或者多个信号的形式。这样的信号可以从因特网网站上下载得到，或者在载体信号上提供，或者以任何其他形式提供。
115.应该注意的是上述实施例对本发明进行说明而不是对本发明进行限制，并且本领域技术人员在不脱离所附权利要求的范围的情况下可设计出替换实施例。在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。本发明可以借助于包括有若干不同元件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的单元权利要求中，这些装置中的若干个可以是通过同一个硬件项来具体体现。单词第一、第二、以及第三等的使用不表示任何顺序。可将这些单词解释为名称。