CN116058874A - 子宫内膜蠕动波的测量方法及设备 - Google Patents

Abstract

本发明提供一种子宫内膜蠕动波的测量方法及设备。该方法包括：向被检者的子宫区域发射第一超声波，并接收子宫区域返回的第一超声回波，得到第一超声回波数据；根据第一超声回波数据确定感兴趣区域，感兴趣区域包括被检者子宫内膜上沿一方向的多个连续或者离散的点；向感兴趣区域发射第二超声波，并接收感兴趣区域返回的第二超声回波，得到第二超声回波数据；基于第二超声回波数据得到感兴趣区域的蠕动运动参数；根据感兴趣区域中的蠕动运动参数确定推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离；输出推动距离。对子宫内膜蠕动波在一方向上传播的距离进行量化与输出，有助于完善子宫内膜蠕动波的测量体系。

Description

子宫内膜蠕动波的测量方法及设备

技术领域

本发明实施例涉及医用超声技术领域，具体涉及一种子宫内膜蠕动波的测量方法及设备。

背景技术

子宫内膜蠕动波是指由子宫肌层收缩带动子宫内膜蠕动而产生的机械波。临床研究表明适当的子宫内膜蠕动可以辅助精子运输与胚胎着床，而过于激烈的子宫内膜蠕动则有可能降低受孕的成功率，因此可以通过子宫内膜蠕动波来对子宫内膜接受胚胎的能力即子宫内膜容受性进行评估。对子宫内膜容受性的正确评估在选择种植时机和评估妊娠率等方面具有重要的临床意义，因此对于子宫内膜蠕动波的测量也具有重要的临床意义。

目前对子宫内膜蠕动波的测量主要是通过对子宫内膜进行超声成像，依靠医生观察所采集的超声影像来对子宫内膜蠕动波的频次、方向和强度做出主观判断。通过测量得到的子宫内膜蠕动波的频次、方向和强度对子宫内膜容受性进行评估。临床研究发现，即使两个患者子宫内膜蠕动波的频次、方向和强度均相同，但是依然可能具有不同的临床表现，因此对于子宫内膜蠕动波的测量仍有待完善。

发明内容

本发明实施例提供一种子宫内膜蠕动波的测量方法及设备，用于弥补现有子宫内膜蠕动波的测量方法中存在的不足。

一个实施例中，提供一种子宫内膜蠕动波的测量方法，包括：

向被检者的子宫区域发射第一超声波，并接收子宫区域返回的第一超声回波，得到第一超声回波数据；

根据第一超声回波数据确定感兴趣区域，感兴趣区域包括被检者子宫内膜上沿一方向的多个连续或者离散的点；

向感兴趣区域发射第二超声波，并接收感兴趣区域返回的第二超声回波，得到第二超声回波数据；

基于第二超声回波数据得到感兴趣区域的蠕动运动参数；

根据感兴趣区域中的蠕动运动参数确定推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离；

输出推动距离。

一个实施例中，提供了一种子宫内膜蠕动波的测量方法，包括：

向被检者的子宫区域发射第一超声波，并接收所述子宫区域返回的第一超声回波，得到第一超声回波数据；

基于所述第一超声回波数据得到所述子宫区域的蠕动运动参数；

根据所述第一超声回波数据确定感兴趣区域，所述感兴趣区域包括所述被检者子宫内膜上沿一方向的多个连续或者离散的点；

根据所述感兴趣区域中的蠕动运动参数确定推动距离，所述推动距离用于度量子宫内膜蠕动波在所述方向上传播的距离；

输出所述推动距离。

一种实施例中，前述方法还包括：

获取被检者的子宫内膜长度；

确定相对推动距离，相对推动距离为推动距离与子宫内膜长度的比值；

输出相对推动距离。

一种实施例中，前述根据感兴趣区域中的蠕动运动参数确定推动距离，包括：

根据感兴趣区域中的蠕动运动参数生成时空分布图，时空分布图用于表示感兴趣区域中各点处的蠕动运动参数随时间的变化；

基于时空分布图确定推动距离。

一种实施例中，前述基于时空分布图确定推动距离，包括：

显示时空分布图；

根据用户在时空分布图上的点选操作或者对光标的调节操作确定子宫内膜蠕动波的起始位置点和终止位置点；

根据起始位置点和终止位置点确定推动距离。

一种实施例中，前述根据感兴趣区域中的蠕动运动参数确定推动距离，包括：

根据感兴趣区域中的蠕动运动参数生成感兴趣区域中各点的蠕动运动参数随时间变化的蠕动曲线；

根据各点的蠕动曲线的极值确定子宫内膜蠕动波的起始位置点和终止位置点；

根据起始位置点和终止位置点确定推动距离。

一种实施例中，前述根据所述感兴趣区域中的蠕动运动参数确定推动距离包括：

基于机器学习或深度学习算法根据所述感兴趣区域中的蠕动运动参数确定所述推动距离。

一种实施例中，前述根据第一超声回波数据确定感兴趣区域，包括：

根据第一超声回波数据生成子宫区域的超声图像；

根据用户在超声图像上的输入确定感兴趣区域，或者，从超声图像中识别出子宫内膜区域并且在子宫内膜区域中确定感兴趣区域。

一种实施例中，前述感兴趣区域包括被检者子宫内膜上沿宫颈至宫底方向的多个连续或者离散的点、沿宫底至宫颈方向的多个连续或者离散的点、沿宫颈至卵巢方向的多个连续或者离散的点和沿卵巢至宫底方向的多个连续或者离散的点中的至少一种。

一种实施例中，前述输出推动距离包括：

输出测量时间内各个子宫内膜蠕动波的推动距离、最大推动距离、最小推动距离和平均推动距离中的至少一项。

一种实施例中，前述输出相对推动距离包括：

输出测量时间内各个子宫内膜蠕动波的相对推动距离、最大相对推动距离、最小相对推动距离和平均相对推动距离中的至少一项。

一种实施例中，前述蠕动运动参数包括蠕动速度、蠕动加速度、组织位移和组织应变中的至少一项。

一种实施例中，提供了一种子宫内膜蠕动波的测量方法，包括：

获取被检者子宫区域的蠕动运动参数；

根据子宫区域的蠕动运动参数确定子宫内膜蠕动波在被检者子宫内膜一方向上的推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离。

一种实施例中，提供了一种超声成像设备，包括：

超声探头；

发射电路，用于按照设定模式将相应的发射序列输出至超声探头，以控制超声探头发射相应的超声波；

接收电路，用于接收超声探头输出的超声回波信号，输出超声回波数据；

显示器，用于输出可视化信息；

处理器，用于执行如上述任一项所述的子宫内膜蠕动波的测量方法。

一种实施例中，提供了一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机执行指令，计算机执行指令被处理器执行时用于实现如上述任一项所述的子宫内膜蠕动波的测量方法。

本发明实施例提供的子宫内膜蠕动波的测量方法及设备，通过向被检者的子宫区域发射第一超声波，并接收子宫区域返回的第一超声回波，得到第一超声回波数据；根据第一超声回波数据确定感兴趣区域，感兴趣区域包括被检者子宫内膜上沿一方向的多个连续或者离散的点；向感兴趣区域发射第二超声波，并接收感兴趣区域返回的第二超声回波，得到第二超声回波数据；基于第二超声回波数据得到感兴趣区域的蠕动运动参数；根据感兴趣区域中的蠕动运动参数确定推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离；输出推动距离。实现了对子宫内膜蠕动波在一方向上传播距离的量化与输出，为子宫内膜蠕动波的临床研究与诊断评估进一步提供了数值依据，有助于完善子宫内膜蠕动波的测量体系。

附图说明

图1为本发明提供的超声成像设备一实施例的结构框图；

图2为本发明提供的子宫内膜蠕动波的测量方法一实施例的流程图；

图3为本发明一实施例提供的显示界面示意图；

图4为本发明一实施例提供的时空分布图的示意图；

图5A和图5B为本发明一实施例提供的基于时空分布图确定推动距离的示意图；

图6为本发明又一实施例提供的基于时空分布图确定推动距离的示意图；

图7为本发明提供的子宫内膜蠕动波的测量方法又一实施例的流程图。

具体实施方式

下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中，很多细节描述是为了使得本申请能被更好的理解。然而，本领域技术人员可以毫不费力的认识到，其中部分特征在不同情况下是可以省略的，或者可以由其他元件、材料、方法所替代。在某些情况下，本申请相关的一些操作并没有在说明书中显示或者描述，这是为了避免本申请的核心部分被过多的描述所淹没，而对于本领域技术人员而言，详细描述这些相关操作并不是必要的，他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。

另外，说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时，方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此，说明书和附图中的各种顺序只是为了清楚描述某一个实施例，并不意味着是必须的顺序，除非另有说明其中某个顺序是必须遵循的。

本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。而本申请所说“连接”、“联接”，如无特别说明，均包括直接和间接连接(联接)。

目前对子宫内膜蠕动波的测量主要是通过对子宫内膜进行超声成像，依靠医生观察所采集的超声影像来对子宫内膜蠕动波的频次、方向和强度做出主观判断。通过测量得到的子宫内膜蠕动波的频次、方向和强度对子宫内膜容受性进行评估。

然而，申请人发现，即使两个患者子宫内膜蠕动波的频次、方向和强度均相同，但是依然可能具有不同的临床表现。比如，两个的频次、方向和强度均相同的蠕动波，临床表现上，可能一个蠕动波是从宫颈完整传递至宫底，而另一个蠕动波仅从宫颈传递至内膜中部。因此对于子宫内膜蠕动波的测量仍有待完善。本发明即提出了一种更完善的蠕动波测量方法和设备。

如图1所示，本发明提供的超声成像设备可以包括：超声探头20、发射/接收电路30(即发射电路310和接收电路320)、波束合成模块40、IQ解调模块50、存储器60、处理器70和人机交互装置。处理器70可以包括控制模块710和图像处理模块720。

超声探头20包括由阵列式排布的多个阵元组成的换能器(图中未示出)，多个阵元排列成一排构成线阵，或排布成二维矩阵构成面阵，多个阵元也可以构成凸阵列。阵元用于根据激励电信号发射超声波束，或将接收的超声波束变换为电信号。因此每个阵元可用于实现电脉冲信号和超声波束的相互转换，从而实现向人体组织的目标区域(例如本实施例中的子宫内膜区域)发射超声波、也可用于接收经组织反射回的超声波的回波。在进行超声检测时，可通过发射电路310和接收电路320控制哪些阵元用于发射超声波束，哪些阵元用于接收超声波束，或者控制阵元分时隙用于发射超声波束或接收超声波束的回波。参与超声波发射的阵元可以同时被电信号激励，从而同时发射超声波；或者参与超声波发射的阵元也可以被具有一定时间间隔的若干电信号激励，从而持续发射具有一定时间间隔的超声波。

本实施例中，用户通过移动超声探头20选择合适的位置和角度向子宫内膜区域10发射超声波并接收由子宫内膜区域10返回的超声波的回波，获得并输出该回波的电信号，回波的电信号是按以接收阵元为通道所形成的通道模拟电信号，其携带有幅度信息、频率信息和时间信息。

发射电路310用于根据处理器70的控制模块710的控制产生发射序列，发射序列用于控制多个阵元中的部分或者全部向生物组织发射超声波，发射序列参数包括发射用的阵元位置、阵元数量和超声波束发射参数(例如幅度、频率、发射次数、发射间隔、发射角度、波型、聚焦位置等)。某些情况下，发射电路310还用于对发射的波束进行相位延迟，使不同的发射阵元按照不同的时间发射超声波，以便各发射超声波束能够在预定的感兴趣区域聚焦。不同的工作模式，例如B图像模式、C图像模式和D图像模式(多普勒模式)，发射序列参数可能不同，回波信号经接收电路320接收并经后续的模块和相应算法处理后，可生成反映组织解剖结构的B图像、反映组织解剖结构和血流信息的C图像以及反映多普勒频谱图像的D图像。

接收电路320用于从超声探头20接收超声回波的电信号，并对超声回波的电信号进行处理。接收电路320可以包括一个或多个放大器、模数转换器(ADC)等。放大器用于在适当增益补偿之后放大所接收到的超声回波的电信号，模数转换器用于对模拟回波信号按预定的时间间隔进行采样，从而转换成数字化的信号，数字化后的回波信号依然保留有幅度信息、频率信息和相位信息。接收电路320输出的数据可输出给波束合成模块40进行处理，或者，输出给存储器60进行存储。

波束合成模块40和接收电路320信号相连，用于对接收电路320输出的信号进行相应的延时和加权求和等波束合成处理，由于被测组织中的超声波接收点到接收阵元的距离不同，因此，不同接收阵元输出的同一接收点的通道数据具有延时差异，需要进行延时处理，将相位对齐，并将同一接收点的不同通道数据进行加权求和，得到波束合成后的超声图像数据，波束合成模块40输出的超声图像数据也称为射频数据(RF数据)。波束合成模块40将射频数据输出至IQ解调模块50。在有的实施例中，波束合成模块40也可以将射频数据输出至存储器60进行缓存或保存，或将射频数据直接输出至处理器70的图像处理模块720进行图像处理。

波束合成模块40可以采用硬件、固件或软件的方式执行上述功能，例如，波束合成模块40可以包括能够根据特定逻辑指令处理输入数据的中央控制器电路(CPU)、一个或多个微处理芯片或其他任何电子部件，当波束合成模块40采用软件方式实现时，其可以执行存储在有形和非暂态计算机可读介质(例如，存储器60)上的指令，以使用任何适当波束合成方法进行波束合成计算。

IQ解调模块50通过IQ解调去除信号载波，提取信号中包含的组织结构信息，并进行滤波去除噪声，此时获取的信号称为基带信号(IQ数据对)。IQ解调模块50将IQ数据对输出至处理器70的图像处理模块720进行图像处理。在有的实施例中，IQ解调模块50还将IQ数据对输出至存储器60进行缓存或保存，以便图像处理模块720从存储器60中读出数据进行后续的图像处理。

处理器70用于配置成能够根据特定逻辑指令处理输入数据的中央控制器电路(CPU)、一个或多个微处理器、图形控制器电路(GPU)或其他任何电子部件，其可以根据输入的指令或预定的指令对外围电子部件执行控制，或对存储器60执行数据读取和/或保存，也可以通过执行存储器60中的程序对输入数据进行处理，例如根据一个或多个工作模式对采集的超声数据执行一个或多个处理操作，处理操作包括但不限于调整或限定超声探头20发出的超声波的形式，生成各种图像帧以供后续人机交互装置的显示器80进行显示，或者调整或限定在显示器80上显示的内容和形式，或者调整在显示器80上显示的一个或多个图像显示设置(例如超声图像、界面组件、定位感兴趣区域)。

图像处理模块720用于对波束合成模块40输出的数据或IQ解调模块50输出的数据进行处理，以生成扫描范围内的信号强弱变化的灰度图像，该灰度图像反映组织内部的解剖结构，称为B图像。图像处理模块720可以将B图像输出至人机交互装置的显示器80进行显示。

人机交互装置用于进行人机交互，即接收用户的输入和输出可视化信息；其接收用户的输入可采用键盘、操作按钮、鼠标、轨迹球等，也可以采用与显示器集成在一起的触控屏；其输出可视化信息采用显示器80。

存储器60可以是有形且非暂态的计算机可读介质，例如可为闪存卡、固态存储器、硬盘等，用于存储数据或者程序，例如，存储器60可以用于存储所采集的超声数据或处理器70所生成的暂不立即显示的图像帧，或者存储器60可以存储图形用户界面、一个或多个默认图像显示设置、用于处理器、波束合成模块或IQ解码模块的编程指令。

需要说明的是，图1所示的结构仅为示意，还可以包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。图1中所示的各组件可以采用硬件和/或软件实现。图1所示的超声成像设备可以用于执行本发明任一实施例所提供的子宫内膜蠕动波的测量方法。

请参考图2，本发明一实施例提供的子宫内膜蠕动波的测量方法可以包括：

S201、向被检者的子宫区域发射第一超声波，并接收子宫区域返回的第一超声回波，得到第一超声回波数据。

本实施例中可以通过移动超声成像设备的超声探头以选择合适的位置和角度向被检者的子宫区域实时发射第一超声波，并接收子宫区域返回的第一超声回波，以便实时获取第一超声回波数据。本实施例中的第一超声回波数据可以是第一超声回波接收路径上任一环节得到的数据，例如可以是RF数据、IQ数据、图像数据等。

S202、根据第一超声回波数据确定感兴趣区域，感兴趣区域包括被检者子宫内膜上沿一方向的多个连续或者离散的点。

为了实现对子宫内膜蠕动波的定向测量，本实施例中在获取第一超声回波数据之后，便可以据此确定测量方向。具体的，可以根据第一超声回波数据将被检者子宫内膜上沿一方向的多个连续或者离散的点确定为感兴趣区域，该方向便是子宫内膜蠕动波的测量方向。本实施例中的感兴趣区域可以是直线、曲线、折线等。

一种可选的实施方式中，根据第一超声回波数据确定感兴趣区域具体可以包括：根据第一超声回波数据生成子宫区域的超声图像；根据用户在超声图像上的输入确定感兴趣区域，或者，从超声图像中识别出子宫内膜区域并且在子宫内膜区域中确定感兴趣区域。

例如可以通过超声成像设备的图像处理模块对第一超声回波数据进行处理，以得到被检者子宫区域的超声图像，包括但不限于B型超声图像、C型超声图像等，然后基于超声图像确定感兴趣区域。基于超声图像确定感兴趣区域可以采用手动方式或者自动方式实现。

当采用手动方式确定感兴趣区域时，可以在超声成像设备的显示屏中显示超声图像，并响应于用户在超声图像上的输入确定感兴趣区域，例如用户可以在超声图像的子宫内膜区域沿着期望进行子宫内膜蠕动波测量的方向，通过鼠标、键盘、触控屏等外部输入设备描记出一条线，可以是直线、曲线、折线等，可以将用户输入的这条线上的点确定为感兴趣区域。

当采用自动方式确定感兴趣区域时，可以不对超声图像进行显示，而是直接从超声图像中自动识别出子宫内膜区域，然后按照预设规则在子宫内膜区域中自动确定出感兴趣区域。例如可以基于深度学习算法预先训练识别模型。可选的，当采用自动方式确定感兴趣区域之后，还可以在超声图像对感兴趣区域进行显示，并且可以根据用户的输入对感兴趣区域所在的位置进行调整。

子宫内膜蠕动波的传播方向多种多样，用户可以根据实际需要确定测量方向，例如可以沿着宫颈到宫底的方向、沿宫底至宫颈方向、沿宫颈至卵巢方向和沿卵巢至宫底方向中的至少一个进行测量，即感兴趣区域可以包括被检者子宫内膜上沿宫颈至宫底方向的多个连续或者离散的点、沿宫底至宫颈方向的多个连续或者离散的点、沿宫颈至卵巢方向的多个连续或者离散的点和沿卵巢至宫底方向的多个连续或者离散的点中的至少一种。当然，也可以是前述方向之外的其他方向。

需要说明的是，当用户需要对多个方向上的子宫内膜蠕动波进行测量时，可以确定多个感兴趣区域，每个感兴趣区域包括一个方向上多个连续或者离散的点，多个感兴趣区域可以部分相交。

S203、向感兴趣区域发射第二超声波，并接收感兴趣区域返回的第二超声回波，得到第二超声回波数据。

在确定感兴趣区域之后，便可以通过超声探头向感兴趣区域实时发射用于获取蠕动运动参数的第二超声波。需要说明的是，本实施例中用于获取感兴趣区域的第一超声波和用于获取蠕动运动参数的第二超声波可以共用相同的扫描序列，也可以采用完全不同的扫描序列。当扫描序列独立时，其发射接收参数(例如频率、聚焦方向、发射间隔、发射位置等)也可以独立设置。本实施例中的第二超声回波数据可以是第二超声回波接收路径上任一环节得到的数据，例如可以是RF数据、IQ数据、图像数据等。

S204、基于第二超声回波数据得到感兴趣区域的蠕动运动参数。

本实施例中可以基于第二超声回波数据生成感兴趣区域的超声图像，通过对感兴趣区域的超声图像进行相关处理，以得到感兴趣区域的蠕动运动参数。可以理解的是，当对空间中某个目标位置持续发射一段时间的超声波并接收返回的超声回波时，如果该目标位置是运动的，则不同时刻得到的超声回波将会发生变化。通过相关对比便可以检测出各时刻下超声回波的变化量或变化速度，基于此便可以确定目标位置的位移或者速度。因此可以通过向蠕动中的子宫内膜持续发射一段时间的第二超声波，然后采用相关对比获得感兴趣区域的蠕动运动参数，蠕动运动参数可以包括蠕动速度、蠕动加速度、组织位移和组织应变中的至少一项。

一种可选的实施方式中，可以通过位移检测方法得到感兴趣区域中的各个点的位移或速度；也可以通过位移检测方法得到感兴趣区域中的各个点的位移，再通过计算位移在时间上的梯度，得到各个点的速度。例如可以基于块匹配的方式，对于某一时刻某位置处的超声回波数据，在另一时刻的超声回波数据的不同位置进行搜索，找出与之互相关最大的位置，两者之间的距离即为这两个时刻下该位置处组织的位移，进一步结合两个时刻之间的时间差便可以得到蠕动速度、蠕动加速度等其他蠕动运动参数。又例如还可以采用基于超声多普勒效应的方式，以类似常规血流成像的原理检测出某位置处的组织在各个时刻的运动速度；或者还可以采用其他基于信号自相关或者互相关的方法来确定感兴趣区域的运动参数，本实施例对此不做限制。

S205、根据感兴趣区域中的蠕动运动参数确定推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离。

基于感兴趣区域中随时间变化的蠕动运动参数确定子宫内膜蠕动波在一方向上传播的距离，即推动距离。推动距离可以作为评价子宫内膜容受性的重要指标之一，有助于对子宫内膜蠕动波进行进一步的定量研究。

需要特别强调的是，推动距离是子宫内膜蠕动波在一方向上传播的距离，而不是某一点的移动距离。下面以感兴趣区域包括沿宫颈到宫底方向上的多个点：S₁、S₂、S₃、S₄、S₅、S₆、S₇、S₈、S₉和S₁₀为例，来对推动距离进行说明。当子宫肌层收缩时带动点S₁移动d₁mm，点S₁带动点S₂移动d₂mm，点S₂带动点S₃移动d₃mm，点S₃带动点S₄移动d₄mm，点S₄带动点S₅移动d₅mm，点S₅带动点S₆移动d₆mm，可以理解的是，随着子宫内膜蠕动波的传播，其强度会逐渐减弱，即d₁>d₂>d₃>d₄>d₅>d₆，当点S₆的移动微弱到不足以带动点S₇移动时，传播终止。则此次子宫内膜蠕动波在宫颈到宫底方向上的推动距离为点S₁至点S₆之间的距离，而d₁是点S₁的移动距离。

假设在测量时间内从被检者的宫颈到宫底方向上检测到2次方向和强度均相同的子宫内膜蠕动波，其中一次从点S₁传播到点S₁₀，即是1次从宫颈完整传播到宫底的蠕动波；而另一次从点S₁传播到点S₅，即是1次从宫颈传播到子宫内膜中部但未到达宫底的蠕动波。若采用现有技术仅从方向和强度上来对子宫内膜蠕动波进行度量，这两次子宫内膜蠕动波将被认为是相同的。而采用本实施例提供的推动距离进行度量时，可知这两次子宫内膜蠕动波从宫颈到宫底方向上传播的距离并不相同，即便这两次子宫内膜蠕动波的方向和强度均相同，也应认为是两种不同的子宫内膜蠕动波，可能有着不同的临床意义，具有潜在研究价值。

S206、输出推动距离。

本实施例中在确定推动距离之后，还可以通过文字、图像、语音、表格等方式输出推动距离，供用户查看。可以获取测量时间内所有子宫内膜蠕动波的推动距离，最大推动距离为测量时间内所有推动距离的最大值，最小推动距离为测量时间内所有推动距离的最小值，平均推动距离为测量时间内所有推动距离的平均值，其反应了子宫内膜蠕动波在测量时间内在测量方向上传播的平均距离。测量时间可以根据观测需求由用户进行设置。推动距离可以在测量过程中进行实时输出，即检测到一个子宫内膜蠕动波，便实时输出其对应的推动距离，同时对最大推动距离、最小推动距离和平均推动距离进行实时更新；也可以在完成所有测量之后进行输出。一种可选的实施方式中，可以输出测量时间内各个子宫内膜蠕动波的推动距离、最大推动距离、最小推动距离和平均推动距离中的至少一项。

为便于用户将推动距离与解剖位置进行对应，可以将子宫内膜区域的超声图像与推动距离显示在同一界面上。请参考图3，图3左侧显示了子宫内膜区域的超声图像，在超声图像中采用折线标识出了感兴趣区域；右侧显示了推动距离，具体包括m次子宫内膜蠕动波的推动距离n₁ cm、n₂ cm、n₃ cm……n_m cm，以及m次子宫内膜蠕动波的平均推动距离n cm。

本实施例提供的子宫内膜蠕动波的测量方法，通过向被检者的子宫区域发射第一超声波，并接收子宫区域返回的第一超声回波，得到第一超声回波数据；根据第一超声回波数据确定感兴趣区域，感兴趣区域包括被检者子宫内膜上沿一方向的多个连续或者离散的点；向感兴趣区域发射第二超声波，并接收感兴趣区域返回的第二超声回波，得到第二超声回波数据；基于第二超声回波数据得到感兴趣区域的蠕动运动参数；根据感兴趣区域中的蠕动运动参数确定推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离；输出推动距离。实现了对子宫内膜蠕动波在一方向上传播距离的量化与输出，为子宫内膜蠕动波的临床研究与诊断评估进一步提供了数值依据，有助于完善子宫内膜蠕动波的测量体系。

本发明的一个实施例中，也可以不专门针对所确定的感兴趣区域发射第二超声波并基于第二超声回波数据获得感兴趣区域的蠕动运动参数，而是直接根据第一超声回波数据计算获得子宫区域的蠕动运动参数，然后再基于第一超声回波数据确定感兴趣区域，并根据该感兴趣区域内的已经计算出的蠕动运动参数确定该推动距离。本实施例中，可以包括下列步骤：

向被检者的子宫区域发射第一超声波，并接收该子宫区域返回的第一超声回波，得到第一超声回波数据；

基于该第一超声回波数据得到该子宫区域的蠕动运动参数；

根据该第一超声回波数据确定感兴趣区域，该感兴趣区域包括该被检者子宫内膜上沿一方向的多个连续或者离散的点；

根据该感兴趣区域中的蠕动运动参数确定推动距离，该推动距离用于度量子宫内膜蠕动波在该方向上传播的距离；

输出该推动距离。

本实施例中，各步骤的具体细节可以与前述实施例中类似，在此不再赘述。

在上述实施例的基础上，下面将进一步阐述如何根据感兴趣区域中的蠕动运动参数确定推动距离。一种可选的实施方式中，根据感兴趣区域中的蠕动运动参数确定推动距离可以包括：根据感兴趣区域中的蠕动运动参数生成时空分布图，时空分布图用于表示感兴趣区域中各点处的蠕动运动参数随时间的变化；基于时空分布图确定推动距离。

子宫内膜蠕动波的时空分布图用于表示感兴趣区域中各点处的蠕动运动参数随时间的变化。首先需要为时空分布图建立坐标系，该坐标系包括第一坐标轴和第二坐标轴，其中第一坐标轴表示时间，第二坐标轴表示空间位置；然后根据蠕动运动参数对应的时间和空间位置，将蠕动运动参数显示在坐标系中相应位置处。在时空分布图中可以采用不同的颜色或灰度表示蠕动运动参数的大小或方向。由于时空分布图包含有子宫内膜蠕动波的时间信息和空间位置信息，因而可以根据时空分布图确定子宫内膜蠕动波的推动距离。

为便于理解，图4示出了一种时空分布图的示意图。在图4所示的时空分布图中，横轴表示时间，纵轴表示空间位置。对于蠕动波，同一次蠕动波可以同时包括正向运动和反向运动，当然也可以只包括其中一种。图4中，深色平行四边形表示正向运动，浅色平行四边形表示同次蠕动波的反向运动。在时空分布图中，深色平行四边形和浅色平行四边形可以同时以成对的方式存在(此时其所代表的蠕动波同时包括正向运动和反向运动)，也可以单独存在(此时其所代表的蠕动波只包括一种运动)。

图4示出了测量时间内两次子宫内膜蠕动波，其中第一次的推动距离大于第二次的推动距离。可以理解的是，若确定了多个感兴趣区域，则可以分别生成与感兴趣区域相对应的多个时空分布图，并基于每个时空分布图得到对应的推动距离，以便进行对比分析。

基于时空分布图确定推动距离既可以采用手动方式进行确定，也可以采用自动方式进行确定。当采用手动方式基于时空分布图确定推动距离时，具体可以包括：显示时空分布图；根据用户在时空分布图上的点选操作或者对光标的调节操作确定子宫内膜蠕动波的起始位置点和终止位置点；根据起始位置点和终止位置点确定推动距离。为便于用户将时空分布图上各点与超声图像中的解剖位置进行对应，还可以将子宫内膜区域的超声图像与时空分布图显示在同一界面上。

根据用户在时空分布图上的点选操作确定推动距离可以参考图5A所示，用户可以在显示的时空分布图中通过外部输入设备进行点选操作，对子宫内膜蠕动波的起始位置点和终止位置点进行标记。如图5A所示，可以根据用户在时空分布图上的点选操作确定第1次子宫内膜蠕动波起始位置点的纵坐标为p₁，终止位置点的纵坐标为p₂，则可以确定第1次子宫内膜蠕动波的推动距离为|p₂-p₁|。对于第2次子宫内膜蠕动波可以采用相同的方法确定推动距离，此处不再赘述。

根据用户在时空分布图上对光标的调节操作确定推动距离可以参考图5B所示，在时空分布图上显示有可调节光标(如图5B中所示的两条虚线)，用户可以通过外部输入设备对光标进行调节操作，以使两个光标分别对应子宫内膜蠕动波的起始位置点和终止位置点，则子宫内膜蠕动波的推动距离为两个光标之间的距离△d。

可以理解的是，采用手动方式确定起始位置点和终止位置点并不限于上述的点选操作和对光标的调节操作，还可以采用如矩形、平行四边形等二维图形将二维时空分布图里某个子宫内膜蠕动波框住，或者针对曲线、三维图等的类似测量操作等。

当采用自动方式基于时空分布图确定推动距离时，具体可以包括：根据时空分布图生成感兴趣区域中各点的蠕动运动参数随时间变化的蠕动曲线；根据各点的蠕动曲线的极值确定子宫内膜蠕动波的起始位置点和终止位置点；根据起始位置点和终止位置点确定推动距离。例如可以通过对时空分布图进行纵向扫描，以获取蠕动曲线。每一条蠕动曲线都展示了感兴趣区域中一个点的蠕动运动参数随时间变化的轨迹。请参考图6所示，图6中左侧为时空分布图，右侧为蠕动曲线示意图，为便于显示图6中仅示出了p₁、p₂和p₃共3个位置点处的蠕动曲线。识别每一条蠕动曲线上的极值，为了排除干扰因素的影响，可以仅识别每一条蠕动曲线上大于预设阈值的极值，记录该极值所对应的空间位置及时间。然后根据子宫内膜蠕动波在空间和时间上的连续性，确定每次子宫内膜蠕动波的起始位置点与终止位置点。如图6所示，第1次子宫内膜蠕动波的起始位置点为p₁，终止位置点为p₂，则第1次子宫内膜蠕动波的推动距离为|p₂-p₁|；第2次子宫内膜蠕动波的起始位置点为p₁，终止位置点为p₃，则第2次子宫内膜蠕动波的推动距离为|p₃-p₁|。

一种可选的实施方式中，根据感兴趣区域中的蠕动运动参数确定推动距离还可以包括：根据感兴趣区域中的蠕动运动参数生成感兴趣区域中各点的蠕动运动参数随时间变化的蠕动曲线；根据各点的蠕动曲线的极值确定子宫内膜蠕动波的起始位置点和终止位置点；根据起始位置点和终止位置点确定推动距离。一方面可以先根据感兴趣区域中的蠕动运动参数生成时空分布图，然后通过对时空分布图进行纵向扫描得到感兴趣区域中各点的蠕动曲线；另一方面也可以根据感兴趣区域中蠕动运动参数所对应的时间和空间位置直接生成蠕动曲线。最后基于蠕动曲线的极值确定推动距离。

除了上述方法外，还可以基于机器学习或深度学习等人工智能方法根据感兴趣区域中的蠕动运动参数自动确定推动距离。具体地，机器学习算法可以通过分析子宫内膜的形态变化规律或者通过观察得到子宫内膜蠕动波的起始位置点和终止位置，从而得到子宫内膜蠕动波的推动距离。

不同的被检者其子宫内膜长度往往不同，为便于对不同被检者的子宫内膜蠕动波进行比较，本发明还提供一种对子宫内膜蠕动的相对推动距离进行测量的方法。如图7所示，本实施例提供的子宫内膜蠕动波的测量方法在上述任一实施例的基础上还可以包括：

S701、获取被检者的子宫内膜长度。

本实施例中被检者的子宫内膜长度可以是预先测量得到并存储于存储设备中的；也可是由用户在基于第一超声回波数据所生成的子宫区域的超声图像上进行手工描记所确定的；还可以是通过人工智能方法对基于第一超声回波数据所得到的超声图像进行处理得到的。本实施例对于获取被检者的子宫内膜长度的具体实施方式不做限制。

S702、确定相对推动距离，相对推动距离为推动距离与子宫内膜长度的比值。

将子宫内膜蠕动波的推动距离与子宫内膜长度的比值确定为该子宫内膜蠕动波的相对推动距离，以消除由子宫内膜长度不同而带来的差异。

S703、输出相对推动距离。

本实施例中在确定相对推动距离之后，还可以通过文字、图像、语音、表格等方式输出相对推动距离，供用户查看。具体的，可以输出测量时间内各个子宫内膜蠕动波的相对推动距离、最大相对推动距离、最小相对推动距离和平均相对推动距离中的至少一项。例如可以在图3所示显示界面中增加对相对推动距离的显示。

本实施例提供的子宫内膜蠕动波的测量方法，在上述任一实施例的基础上，进一步的将推动距离与子宫内膜长度的比值确定为相对推动距离，并输出相对推动距离，有效地消除了由子宫内膜长度不同而带来的差异，便于在不同被检者之间对子宫内膜蠕动波进行比较。

本发明实施例还提供一种对子宫内膜蠕动波进行离线测量的方法，可以基于存储器中所存储的超声回波数据对子宫内膜蠕动波进行离线测量。例如医护人员在完成对被检者的扫查后，可以将得到的超声回波数据和/或根据超声回波数据所确定的蠕动运动参数存储在存储设备中，或者，发送给远程的其他医护人员，以进行远程会诊。具体可以包括：获取被检者子宫区域的蠕动运动参数；根据子宫区域的蠕动运动参数确定子宫内膜蠕动波在被检者子宫内膜一方向上的推动距离，推动距离用于度量子宫内膜蠕动波在一方向上传播的距离。可以从存储设备中直接获取被检者子宫区域的蠕动运动参数，也可以从存储设备中获取子宫区域的超声回波数据，然后根据超声回波数据确定被检者子宫区域的蠕动运动参数。而在得到蠕动运动参数之后，便可以根据上述任一实施例提供的方法来确定推动距离。

本文参照了各种示范实施例进行说明。然而，本领域的技术人员将认识到，在不脱离本文范围的情况下，可以对示范性实施例做出改变和修正。例如，各种操作步骤以及用于执行操作步骤的组件，可以根据特定的应用或考虑与系统的操作相关联的任何数量的成本函数以不同的方式实现(例如一个或多个步骤可以被删除、修改或结合到其他步骤中)。

另外，如本领域技术人员所理解的，本文的原理可以反映在计算机可读存储介质上的计算机程序产品中，该可读存储介质预装有计算机可读程序代码。任何有形的、非暂时性的计算机可读存储介质皆可被使用，包括磁存储设备(硬盘、软盘等)、光学存储设备(CD-ROM、DVD、Blu Ray盘等)、闪存和/或诸如此类。这些计算机程序指令可被加载到通用计算机、专用计算机或其他可编程数据处理设备上以形成机器，使得这些在计算机上或其他可编程数据处理装置上执行的指令可以生成实现指定的功能的装置。这些计算机程序指令也可以存储在计算机可读存储器中，该计算机可读存储器可以指示计算机或其他可编程数据处理设备以特定的方式运行，这样存储在计算机可读存储器中的指令就可以形成一件制造品，包括实现指定功能的实现装置。计算机程序指令也可以加载到计算机或其他可编程数据处理设备上，从而在计算机或其他可编程设备上执行一系列操作步骤以产生一个计算机实现的进程，使得在计算机或其他可编程设备上执行的指令可以提供用于实现指定功能的步骤。

虽然在各种实施例中已经示出了本文的原理，但是许多特别适用于特定环境和操作要求的结构、布置、比例、元件、材料和部件的修改可以在不脱离本披露的原则和范围内使用。以上修改和其他改变或修正将被包含在本文的范围之内。

前述具体说明已参照各种实施例进行了描述。然而，本领域技术人员将认识到，可以在不脱离本披露的范围的情况下进行各种修正和改变。因此，对于本披露的考虑将是说明性的而非限制性的意义上的，并且所有这些修改都将被包含在其范围内。同样，有关于各种实施例的优点、其他优点和问题的解决方案已如上所述。然而，益处、优点、问题的解决方案以及任何能产生这些的要素，或使其变得更明确的解决方案都不应被解释为关键的、必需的或必要的。本文中所用的术语“包括”和其任何其他变体，皆属于非排他性包含，这样包括要素列表的过程、方法、文章或设备不仅包括这些要素，还包括未明确列出的或不属于该过程、方法、系统、文章或设备的其他要素。此外，本文中所使用的术语“耦合”和其任何其他变体都是指物理连接、电连接、磁连接、光连接、通信连接、功能连接和/或任何其他连接。

以上应用了具体个例对本发明进行阐述，只是用于帮助理解本发明，并不用以限制本发明。对于本发明所属技术领域的技术人员，依据本发明的思想，还可以做出若干简单推演、变形或替换。

Claims (15)

1.一种子宫内膜蠕动波的测量方法，其特征在于，包括：

向被检者的子宫区域发射第一超声波，并接收所述子宫区域返回的第一超声回波，得到第一超声回波数据；

根据所述第一超声回波数据确定感兴趣区域，所述感兴趣区域包括所述被检者子宫内膜上沿一方向的多个连续或者离散的点；

向所述感兴趣区域发射第二超声波，并接收所述感兴趣区域返回的第二超声回波，得到第二超声回波数据；

基于所述第二超声回波数据得到所述感兴趣区域的蠕动运动参数；

根据所述感兴趣区域中的蠕动运动参数确定推动距离，所述推动距离用于度量子宫内膜蠕动波在所述方向上传播的距离；

输出所述推动距离。

2.一种子宫内膜蠕动波的测量方法，其特征在于，包括：

向被检者的子宫区域发射第一超声波，并接收所述子宫区域返回的第一超声回波，得到第一超声回波数据；

基于所述第一超声回波数据得到所述子宫区域的蠕动运动参数；

根据所述第一超声回波数据确定感兴趣区域，所述感兴趣区域包括所述被检者子宫内膜上沿一方向的多个连续或者离散的点；

根据所述感兴趣区域中的蠕动运动参数确定推动距离，所述推动距离用于度量子宫内膜蠕动波在所述方向上传播的距离；

输出所述推动距离。

3.如权利要求1或者2所述的方法，其特征在于，所述方法还包括：

获取所述被检者的子宫内膜长度；

确定相对推动距离，所述相对推动距离为所述推动距离与所述子宫内膜长度的比值；

输出所述相对推动距离。

4.如权利要求1或者2所述的方法，其特征在于，所述根据所述感兴趣区域中的蠕动运动参数确定推动距离，包括：

根据所述感兴趣区域中的蠕动运动参数生成时空分布图，所述时空分布图用于表示所述感兴趣区域中各点处的蠕动运动参数随时间的变化；

基于所述时空分布图确定所述推动距离。

5.如权利要求4所述的方法，其特征在于，所述基于所述时空分布图确定所述推动距离，包括：

显示所述时空分布图；

根据用户在所述时空分布图上的点选操作或者对光标的调节操作确定子宫内膜蠕动波的起始位置点和终止位置点；

根据所述起始位置点和所述终止位置点确定所述推动距离。

6.如权利要求1或者2所述的方法，其特征在于，所述根据所述感兴趣区域中的蠕动运动参数确定推动距离，包括：

根据所述感兴趣区域中的蠕动运动参数生成所述感兴趣区域中各点的蠕动运动参数随时间变化的蠕动曲线；

根据各点的蠕动曲线的极值确定子宫内膜蠕动波的起始位置点和终止位置点；

根据所述起始位置点和所述终止位置点确定所述推动距离。

7.如权利要求1或者2所述的方法，其特征在于，所述根据所述感兴趣区域中的蠕动运动参数确定推动距离包括：

基于机器学习或深度学习算法根据所述感兴趣区域中的蠕动运动参数确定所述推动距离。

8.如权利要求1或者2所述的方法，其特征在于，所述根据所述第一超声回波数据确定感兴趣区域，包括：

根据所述第一超声回波数据生成所述子宫区域的超声图像；

根据用户在所述超声图像上的输入确定所述感兴趣区域，或者，从所述超声图像中识别出子宫内膜区域并且在所述子宫内膜区域中确定所述感兴趣区域。

9.如权利要求1或者2所述的方法，其特征在于，所述感兴趣区域包括所述被检者子宫内膜上沿宫颈至宫底方向的多个连续或者离散的点、沿宫底至宫颈方向的多个连续或者离散的点、沿宫颈至卵巢方向的多个连续或者离散的点和沿卵巢至宫底方向的多个连续或者离散的点中的至少一种。

10.如权利要求1或者2所述的方法，其特征在于，所述输出所述推动距离，包括：

输出测量时间内各个子宫内膜蠕动波的推动距离、最大推动距离、最小推动距离和平均推动距离中的至少一项。

11.如权利要求3所述的方法，其特征在于，所述输出所述相对推动距离，包括：

输出测量时间内各个子宫内膜蠕动波的相对推动距离、最大相对推动距离、最小相对推动距离和平均相对推动距离中的至少一项。

12.如权利要求1或者2所述的方法，其特征在于，所述蠕动运动参数包括蠕动速度、蠕动加速度、组织位移和组织应变中的至少一项。

13.一种子宫内膜蠕动波的测量方法，其特征在于，包括：

获取被检者子宫区域的蠕动运动参数；

根据所述子宫区域的蠕动运动参数确定子宫内膜蠕动波在所述被检者子宫内膜一方向上的推动距离，所述推动距离用于度量子宫内膜蠕动波在所述方向上传播的距离。

14.一种超声成像设备，其特征在于，包括：

超声探头；

发射电路，用于按照设定模式将相应的发射序列输出至所述超声探头，以控制所述超声探头发射相应的超声波；

接收电路，用于接收所述超声探头输出的超声回波信号，输出超声回波数据；

显示器，用于输出可视化信息；

处理器，用于执行如权利要求1-13任一项所述的子宫内膜蠕动波的测量方法。

15.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机执行指令，所述计算机执行指令被处理器执行时用于实现如权利要求1-13任一项所述的子宫内膜蠕动波的测量方法。

Images