CN112674799A - 超声弹性成像方法、电子设备及存储介质 - Google Patents

超声弹性成像方法、电子设备及存储介质 Download PDF

Info

Publication number
CN112674799A
CN112674799A CN202110006122.1A CN202110006122A CN112674799A CN 112674799 A CN112674799 A CN 112674799A CN 202110006122 A CN202110006122 A CN 202110006122A CN 112674799 A CN112674799 A CN 112674799A
Authority
CN
China
Prior art keywords
sliding window
echo data
ultrasonic echo
data
correction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110006122.1A
Other languages
English (en)
Other versions
CN112674799B (zh
Inventor
朱超超
付传卿
安力春
丁勇
王�琦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao Hisense Medical Equipment Co Ltd
Original Assignee
Qingdao Hisense Medical Equipment Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Qingdao Hisense Medical Equipment Co Ltd filed Critical Qingdao Hisense Medical Equipment Co Ltd
Priority to CN202110006122.1A priority Critical patent/CN112674799B/zh
Publication of CN112674799A publication Critical patent/CN112674799A/zh
Application granted granted Critical
Publication of CN112674799B publication Critical patent/CN112674799B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本申请公开了一种超声弹性成像方法、电子设备及存储介质,用于解决相关技术仅对深度方向进行应变估计导致超声弹性成像的结果的可信度不高的问题。本申请实施例中,首先采用二维窗口对施压前后的超声回波数据进行处理,通过利用超声回波数据的深度方向数据和横向方向数据进行散射子位置校正,然后基于散射子位置校正结果来进一步进行应变估计,是的应变估计结果中能够融合深度和横向方向的信息,故而使得最终的应变估计结果更为准确。

Description

超声弹性成像方法、电子设备及存储介质
技术领域
本申请涉及医用超声设备技术领域,特别涉及一种超声弹性成像方法、电子设备及存储介质。
背景技术
对于生命体的生命组织,当有病变时可能会引起组织内部病变导致组织内部的弹性不同。例如有肿块、结节的部位一般比较″硬″,而正常皮肤较有弹性柔软。
在超声系统中,由于正常组织和病变组织在弹性分布上有明显差异,可对组织施加作用力后让组织产生相应的形变,通过对组织的形变数据进行分析处理,可得到组织的弹性分布,从而依据组织的弹性分布来对组织进行诊断分析。
相关技术中,通常通过信号处理方法能够将组织的弹性分布信息给提取出来。其大致处理过程可包括获取超声回波数据后,得到位移数据,然后对位移数据进行应变估计得到组织的弹性分布信息。然而,相关技术中对组织的应变估计往往只考虑到了深度方向上的位移变化,并没有考虑到水平方向的位移变化,因此会产生误差,造成超声弹性成像的结果的可信度不高。
发明内容
本申请的目的是提供一种医用超声设备的扫描控制方法、电子设备及存储介质,以克服相关技术仅对深度方向进行应变估计导致超声弹性成像的结果的可信度不高的问题。
第一方面,本申请提供一种超声弹性成像方法,所述方法包括:
获取对目标组织施压前的第一超声回波数据和对所述目标组织施压后的第二超声回波数据;其中,所述第一超声波数据和所述第二超声波数据中均包括所述目标组织的深度方向以及横向方向,所述横向方向垂直于所述深度方向;
基于第一滑动窗口的尺寸对所述第一超声回波数据和所述第二超声回波数据进行采样,得到每个第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据;其中,所述第一滑动窗口的第一边平行于所述深度方向,且所述第一滑动窗口的第二边平行于所述横向方向;
针对每个所述第一滑动窗口位置分别执行:基于所述第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据,对所述第一滑动窗口位置的散射子位置进行校正,得到所述第一滑动窗口位置的散射子位置校正结果;
基于每个所述第一滑动窗口位置分别对应的散射子位置校正结果,对所述第一超声回波数据和所述第二超声回波数据的位移数据进行应变估计,得到所述目标组织的超声弹性成像结果。
在一些实施例中,每个所述第一滑动窗口位置内均包括多条接收线,且基于每条接收线上的所述第一超声回波数据和所述第二超声回波数据得到每条接收线各自的校正因子,所述基于所述第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据,对所述第一滑动窗口位置的散射子位置进行校正,得到所述第一滑动窗口位置的散射子位置校正结果,包括:
采用所述第一滑动窗口内的每条接收线各自的校正因子以及校正关系,对所述第一滑动窗口的散列子位置进行校正,得到所述散射子位置校正结果;
其中,在所述校正关系中,包括所述深度方向和所述横向方向上对不同接收线的校正因子的累加结果,以及所述深度方向和所述横向方向上对不同接收线的校正因子参考量的累加结果;且,所述校正因子参考量与所述校正因子和所述第一滑动窗口的指定边长均成正比;
所述散射子位置校正结果与所述校正因子的累加结果成反比、且与所述校正因子参考量的累加结果成正比。
在一些实施例中,针对每条接收线、所述校正因子与所述接收线上的所述第一超声回波数据和所述第二超声回波数据的信号指标成正比,所述信号指标包括信号幅度和相位中的至少一种。
在一些实施例中,所述信号指标为所述第一超声波数据和所述第二超声波数据的信号幅度的乘积,或者为所述第一超声波数据和所述第二超声波数据的所述相位之和。
在一些实施例中,所述校正关系表达式为:
Figure BDA0002883445970000031
其中,P表示所述散射子位置校正结果、winH表示所述第一滑动窗口的所述第一边的长度、winW表示所述第一滑动窗口的所述第二边的长度、f(m,n)表示位于(m,n)位置的接收线上的校正因子、b表示常量。
在一些实施例中,不同生物组织作为所述目标组织时,各生物组织分别对应有各自的所述第一滑动窗口的尺寸。
在一些实施例中,所述基于每个所述第一滑动窗口位置分别对应的散射子位置校正结果,对所述第一超声回波数据和所述第二超声回波数据的位移数据进行应变估计,包括:
基于第二滑动窗口的尺寸,对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算,得到每个所述第二滑动窗口位置分别对应的应变估计结果。
在一些实施例中,所述对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算,包括:
基于以下应变估计公式,对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算:
Figure BDA0002883445970000041
其中,Strain表示所述第二滑动窗口位置的应变估计结果、winW1表示所述第二滑动窗口的宽度、winH1表示所述第二滑动窗口的高度、Disp表示位移数据、(i,j)表示位置坐标。
第二方面,本申请另一实施例还提供了一种电子设备,包括至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行本申请实施例提供的任一超声弹性成像方法。
第三方面,本申请另一实施例还提供了一种计算机存储介质,其中,所述计算机存储介质存储有计算机程序,所述计算机程序用于使计算机执行本申请实施例中的任一超声弹性成像方法。
本申请实施例中首先采用二维窗口对施压前后的超声回波数据进行处理,通过利用超声回波数据的深度方向数据和横向方向数据进行散射子位置校正,然后基于散射子位置校正结果来进一步进行应变估计,是的应变估计结果中能够融合深度和横向方向的信息,故而使得最终的应变估计结果更为准确。
本申请的其它特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本申请而了解。本申请的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对本申请实施例中所需要使用的附图作简单地介绍,显而易见地,下面所介绍的附图仅仅是本申请的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为根据本申请一个实施例的应用场景示意图;
图2为根据本申请一个实施例的电子设备的结构示意图;
图3为根据本申请一个实施例的超声弹性成像方法的流程示意图;
图4和图5为根据本申请实施例提供的第一滑动窗口和超声回波数据的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。其中,所描述的实施例是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本申请保护的范围。
并且,在本申请实施例的描述中,除非另有说明,″/″表示或的意思,例如,A/B可以表示A或B;文本中的″和/或″仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况,另外,在本申请实施例的描述中,″多个″是指两个或多于两个。
以下,术语″第一″、″第二″仅用于描述目的,而不能理解为暗示或暗示相对重要性或者隐合指明所指示的技术特征的数量。由此,限定有″第一″、″第二″、的特征可以明示或者隐合地包括一个或者更多个该特征,在本申请实施例的描述中,除非另有说明,″多个″的含义是两个或两个以上。
相关技术中,对组织的应变估计往往只考虑到了深度方向上的位移变化,并没有考虑到水平方向的位移变化,因此会产生误差,造成弹性成像的结果的可信度不高,从而会影响诊断结果。
为了提高超声弹性成像的准确性,本发明提出了一种超声弹性成像方法与电子设备以解决该技术问题。
在本申请实施例中,考虑到组织的运动的复杂性,利用二维窗口来进行散射子位置校正,位置校正过程中参考了组织的深度变化和横向变化,故此能够得到更加准确和可靠的组织弹性分布信息。
下面结合附图,对本申请提供的医用超声设备的扫描控制方法进行说明。
如图1所示示出了本申请实施例提供的应用场景示意图。在该应用场景中,该应用环境中例如可以包括网络10、控制设备20、医用超声设备30以及数据库40。其中:
控制设备20用于对医用超声设备的扫描功能进行控制。控制设备可以独立于医用超声设备也可以设置在超声控制设备内部。
数据库中用于存储超声设备30的相关参数,也可以存储超声设备采集的超声图像以便于调用。
本申请实施例还提供一种电子设备,用于执行医用超声设备的超声弹性成像方法。本申请实施例中的电子设备可以至少包括至少一个处理器、以及至少一个存储器。其中,存储器存储有程序代码,当程序代码被处理器执行时,使得处理器执行本说明书根据本申请各种示例性实施方式的医用超声设备的超声弹性成像方法中的步骤。下面参照图2来描述根据本申请的这种实施方式的电子设备60。图2的电子设备60仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。
如图2,电子设备60以通用电子设备的形式表现。电子设备60的组件可以包括但不限于:上述至少一个处理单元21、上述至少一个存储单元22、连接不同系统组件(包括存储单元22和处理单元21)的总线23。
总线23表示几类总线结构中的一种或多种,包括存储器总线或者存储器控制器、外围总线、处理器或者使用多种总线结构中的任意总线结构的局域总线。
存储单元22可以包括易失性存储器形式的可读介质,例如随机存取存储器(RAM)221和/或高速缓存存储单元222,还可以进一步包括只读存储器(ROM)223。
存储单元22还可以包括具有一组(至少一个)程序模块224的程序/实用工具225,这样的程序模块224包括但不限于:操作系统、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。
电子设备60也可以与一个或多个外部设备24(例如键盘、指向设备、监控设备等)通信,还可与一个或者多个使得用户能与电子设备60交互的设备通信,和/或与使得该电子设备60能与一个或多个其它电子设备进行通信的任何设备(例如路由器、调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口25进行。并且,电子设备60还可以通过网络适配器26与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。如图2所示,网络适配器26通过总线23与电子设备60的其它模块通信。应当理解,尽管图中未示出,可以结合电子设备60使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理器、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
在一些可能的实施方式中,本申请提供的医用超声设备的超声弹性成像方法的各个方面还可以实现为一种程序产品的形式,其包括程序代码,当程序产品在计算机设备上运行时,程序代码用于使计算机设备执行本说明书描述的根据本申请各种示例性实施方式的医用超声设备的超声弹性成像方法中的步骤。
程序产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以是但不限于电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
为了便于理解本申请实施例提供的医用超声设备的超声弹性成像方法,下面结合附图,对此进行进一步说明。
如图3所示,为本申请实施例中的医用超声设备的超声弹性成像方法的流程示意图,包括以下步骤:
在步骤301中,获取对目标组织施压前的第一超声回波数据和对所述目标组织施压后的第二超声回波数据;其中,所述第一超声波数据和所述第二超声波数据中均包括所述目标组织的深度方向以及横向方向,所述横向方向垂直于所述深度方向;
其中,如图4所示,每条接收线(如图4中的一列方向)上具有多个点,接收线方向上表示了组织的深度方向,垂直于接收线方向上(即图4中的水平方向)即为横向方向的数据。
在步骤302中,基于第一滑动窗口的尺寸对所述第一超声回波数据和所述第二超声回波数据进行采样,得到每个第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据;其中,所述第一滑动窗口的第一边平行于所述深度方向,且所述第一滑动窗口的第二边平行于所述横向方向;
为了采用目标组织受压力影响下,在深度方向和横向方向的运动信息进行应变估计。本申请实施例中,第一滑动窗口如图4中的矩形框A和矩形框B。每第一滑动窗口位置内均包括多条接收线。由此,每个第一滑动窗口位置内均可以采集到想要的第一超声回波数据和第二超声回波数据。如图4所示A滑动窗口内包括多条接收线(其中,一个点表示一条接收线位置)。想要的因为超声回波数据包括施压前后的超声回波数据,故此每条接收线对应有目标组织被施压前后的超声回波数据。
在步骤303中,针对每个所述第一滑动窗口位置分别执行:基于所述第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据,对所述第一滑动窗口位置的散射子位置进行校正,得到所述第一滑动窗口位置的散射子位置校正结果;
为了能够准确的校正散射子位置,本申请实施例中,每条接收线上的所述第一超声回波数据和所述第二超声回波数据得到每条接收线各自的校正因子。该校正因子可用于对散射子位置进行校正。实施时,校正因子与接收线上的第一超声回波数据和第二超声回波数据的信号指标成正比,所述信号指标例如可包括信号幅度和相位中的至少一种。由此,校正因子可以根据超声回波信号的情况来确定,能够反映超声回波信号的分布,间接反映目标组织施压前后的运动情况。
在实施时,信号指标为第一超声波数据和第二超声波数据的信号幅度的乘积,或者为第一超声波数据和第二超声波数据的相位之和。由此,假若针对接收线L1,其施压前的超声回波信号为X1、施压后的超声回波信号为X2,则该接收线L1上的校正因子为超声回波信号为X1和超声回波信号为X2的信号幅度的乘积,或者为这两个信号的相位之和。
实施时,为了能够结合纵向和横向的信号准确的校正散射子位置,本申请实施例中,针对每个第一滑动窗口,采用该第一滑动窗口内的每条接收线各自的校正因子以及校正关系,对该第一滑动窗口的散列子位置进行校正,得到该第一滑动窗口的散射子位置校正结果。其中,在校正关系中,包括深度方向和横向方向上对不同接收线的校正因子的累加结果,以及深度方向和横向方向上对不同接收线的校正因子参考量的累加结果;且,校正因子参考量与校正因子和第一滑动窗口的指定边长均成正比;相应的,散射子位置校正结果与校正因子的累加结果成反比、且与校正因子参考量的累加结果成正比。
由此,能够结合第一滑动窗口的深度方向和横向方向的校正因子对散射子位置进行校正,从而实现结合信号的纵向和横向信息来校正散射子位置。以便于后续根据散射子位置进行应变估计时,能够融合横向运动信息。
一种可能的实施方式中,校正关系表达式为公式(1)所示:
Figure BDA0002883445970000091
其中,P表示散射子位置校正结果、winH表示第一滑动窗口的第一边的长度及第一滑动窗口的高度、winW表示第一滑动窗口的第二边的长度即第一滑动窗口的宽度、f(m,n)表示位于(m,n)位置的接收线上的校正因子、b表示常量,通常可取值为1。
实施时,第一滑动窗口可每间隔一定数量超声回波数据点进行移动,如图5所示,在横向方向上第一滑动窗口位置A与第一滑动窗口位置B具有重叠区域。在纵向方向上第一滑动窗口位置A与第一滑动窗口位置C具有重叠区域。一个第一滑动窗口位置可对应的的一个散射子位置校正结果。实施时,第一滑动窗口的尺寸可基于目标组织的类型确定。例如心脏的第一滑动窗口的尺寸与扁桃体腺体的第一滑动窗口的尺寸可不同。具体尺寸可基于不同组织部位的实验结果确定。由此,不同生物组织作为所述目标组织时,各生物组织分别对应有各自的所述第一滑动窗口的尺寸,由此可提高对不同生物组织进行超声弹性成像结果的准确性。
在得到每个第一滑动窗口位置对应的散射子校正位置之后,可以在步骤304中,基于每个所述第一滑动窗口位置分别对应的散射子位置校正结果,对所述第一超声回波数据和所述第二超声回波数据的位移数据进行应变估计,得到所述目标组织的超声弹性成像结果。
例如,实施时可以采用二维窗口的差分估计方法进行应变估计。二维窗口包括了深度方向和横向方向的超声回波信号,故此能够得到准确的应变估计结果。
例如,可以基于第二滑动窗口的尺寸,对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算,得到每个所述第二滑动窗口位置分别对应的应变估计结果。
一种差分运算的表达式如公式(2)所示:
差分运算:
Figure BDA0002883445970000101
其中,Strain表示所述第二滑动窗口位置的应变估计结果、winW1表示所述第二滑动窗口的宽度、winH1表示所述第二滑动窗口的高度、Disp表示位移数据、(i,j)表示位置坐标。
由此,基于横向和纵向两个方向通过差分运算能够得到融合了纵向深度信息和横向运动信息的应变估计结果。
例如,假设原始的超声回波数据量为P1*P2。假设第一滑动窗口的尺寸为m*n,则可得到(P1-m+1)*(P2-n+1)个数据,则对这些数据的施压前后的超声回波数据进行应变估计。例如,假若原始的超声回波数据量为100*100,第一滑动窗口的尺寸为5*5,则可以得到96*96个位置的散射子校正结果。对这些位置的超声回波数据进行应变估计,加上第二滑动窗口尺寸为7*7,则最终得到大小为90*90的超声弹性成像结果。
在一些可能的实施方式中,本申请提供的一种视频展示方法的各个方面还可以实现为一种程序产品的形式,其包括程序代码,当程序产品在计算机设备上运行时,程序代码用于使计算机设备执行本说明书上述描述的根据本申请各种示例性实施方式的一种视频展示方法中的步骤。
程序产品可以采用一个或多个可读介质的任意组合。可读介质可以是可读信号介质或者可读存储介质。可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。
本申请的实施方式的用于视频展示的程序产品可以采用便携式紧凑盘只读存储器(CD-ROM)并包括程序代码,并可以在电子设备上运行。然而,本申请的程序产品不限于此,在本文件中,可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了可读程序代码。这种传播的数据信号可以采用多种形式,包括一一但不限于——电磁信号、光信号或上述的任意合适的组合。可读信号介质还可以是可读存储介质以外的任何可读介质,该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于——无线、有线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言的任意组合来编写用于执行本申请操作的程序代码,程序设计语言包括面向对象的程序设计语言一诸如Java、C++等,还包括常规的过程式程序设计语言一诸如″C″语言或类似的程序设计语言。程序代码可以完全地在用户电子设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户电子设备上部分在远程电子设备上执行、或者完全在远程电子设备或服务端上执行。在涉及远程电子设备的情形中,远程电子设备可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)一连接到用户电子设备,或者,可以连接到外部电子设备(例如利用因特网服务提供商来通过因特网连接)。
应当注意,尽管在上文详细描述中提及了装置的若干单元或子单元,但是这种划分仅仅是示例性的并非强制性的。实际上,根据本申请的实施方式,上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。反之,上文描述的一个单元的特征和功能可以进一步划分为由多个单元来具体化。
此外,尽管在附图中以特定顺序描述了本申请方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程图像缩放设备的处理器以产生一个机器,使得通过计算机或其他可编程图像缩放设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程图像缩放设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程图像缩放设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本申请的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例做出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本申请范围的所有变更和修改。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。

Claims (10)

1.一种超声弹性成像方法,其特征在于,所述方法包括:
获取对目标组织施压前的第一超声回波数据和对所述目标组织施压后的第二超声回波数据;其中,所述第一超声波数据和所述第二超声波数据中均包括所述目标组织的深度方向以及横向方向,所述横向方向垂直于所述深度方向;
基于第一滑动窗口的尺寸对所述第一超声回波数据和所述第二超声回波数据进行采样,得到每个第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据;其中,所述第一滑动窗口的第一边平行于所述深度方向,且所述第一滑动窗口的第二边平行于所述横向方向;
针对每个所述第一滑动窗口位置分别执行:基于所述第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据,对所述第一滑动窗口位置的散射子位置进行校正,得到所述第一滑动窗口位置的散射子位置校正结果;
基于每个所述第一滑动窗口位置分别对应的散射子位置校正结果,对所述第一超声回波数据和所述第二超声回波数据的位移数据进行应变估计,得到所述目标组织的超声弹性成像结果。
2.根据权利要求1所述的方法,其特征在于,每个所述第一滑动窗口位置内均包括多条接收线,且基于每条接收线上的所述第一超声回波数据和所述第二超声回波数据得到每条接收线各自的校正因子,所述基于所述第一滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据,对所述第一滑动窗口位置的散射子位置进行校正,得到所述第一滑动窗口位置的散射子位置校正结果,包括:
采用所述第一滑动窗口内的每条接收线各自的校正因子以及校正关系,对所述第一滑动窗口的散列子位置进行校正,得到所述散射子位置校正结果;
其中,在所述校正关系中,包括所述深度方向和所述横向方向上对不同接收线的校正因子的累加结果,以及所述深度方向和所述横向方向上对不同接收线的校正因子参考量的累加结果;且,所述校正因子参考量与所述校正因子和所述第一滑动窗口的指定边长均成正比;
所述散射子位置校正结果与所述校正因子的累加结果成反比、且与所述校正因子参考量的累加结果成正比。
3.根据权利要求2所述的方法,其特征在于,针对每条接收线、所述校正因子与所述接收线上的所述第一超声回波数据和所述第二超声回波数据的信号指标成正比,所述信号指标包括信号幅度和相位中的至少一种。
4.根据权利要求3所述的方法,其特征在于,所述信号指标为所述第一超声波数据和所述第二超声波数据的信号幅度的乘积,或者为所述第一超声波数据和所述第二超声波数据的所述相位之和。
5.根据权利要求2-3中任一所述的方法,其特征在于,所述校正关系表达式为:
Figure FDA0002883445960000021
其中,P表示所述散射子位置校正结果、winH表示所述第一滑动窗口的所述第一边的长度、winW表示所述第一滑动窗口的所述第二边的长度、f(m,n)表示位于(m,n)位置的接收线上的校正因子、b表示常量。
6.根据权利要求1-4中任一所述的方法,其特征在于,不同生物组织作为所述目标组织时,各生物组织分别对应有各自的所述第一滑动窗口的尺寸。
7.根据权利要求1-4中任一所述的方法,其特征在于,所述基于每个所述第一滑动窗口位置分别对应的散射子位置校正结果,对所述第一超声回波数据和所述第二超声回波数据的位移数据进行应变估计,包括:
基于第二滑动窗口的尺寸,对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算,得到每个所述第二滑动窗口位置分别对应的应变估计结果。
8.根据权利要求7所述的方法,其特征在于,所述对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算,包括:
基于以下应变估计公式,对每个所述第二滑动窗口位置内的所述第一超声回波数据和所述第二超声回波数据的位移数据以及散射子位置校正结果进行差分运算:
Figure FDA0002883445960000031
其中,Strain表示所述第二滑动窗口位置的应变估计结果、winW1表示所述第二滑动窗口的宽度、winH1表示所述第二滑动窗口的高度、Disp表示位移数据、(i,j)表示位置坐标。
9.一种电子设备,其特征在于,包括至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行权利要求1-8中任一项所述的方法。
10.一种计算机存储介质,其特征在于,所述计算机存储介质存储有计算机程序,所述计算机程序用于使计算机执行权利要求1-8中任一项所述的方法。
CN202110006122.1A 2021-01-05 2021-01-05 超声弹性成像方法、电子设备及存储介质 Active CN112674799B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110006122.1A CN112674799B (zh) 2021-01-05 2021-01-05 超声弹性成像方法、电子设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110006122.1A CN112674799B (zh) 2021-01-05 2021-01-05 超声弹性成像方法、电子设备及存储介质

Publications (2)

Publication Number Publication Date
CN112674799A true CN112674799A (zh) 2021-04-20
CN112674799B CN112674799B (zh) 2022-11-25

Family

ID=75457170

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110006122.1A Active CN112674799B (zh) 2021-01-05 2021-01-05 超声弹性成像方法、电子设备及存储介质

Country Status (1)

Country Link
CN (1) CN112674799B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113628183A (zh) * 2021-08-06 2021-11-09 青岛海信医疗设备股份有限公司 一种超声检测对象的容积确定方法及超声设备

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6508768B1 (en) * 2000-11-22 2003-01-21 University Of Kansas Medical Center Ultrasonic elasticity imaging
US20070083116A1 (en) * 2005-09-21 2007-04-12 Fuji Photo Film Co., Ltd. Ultrasonic diagnostic apparatus
US20070167772A1 (en) * 2005-12-09 2007-07-19 Aloka Co., Ltd. Apparatus and method for optimized search for displacement estimation in elasticity imaging
US20100256494A1 (en) * 2007-11-16 2010-10-07 Takashi Azuma Ultrasonic imaging system
US20120065505A1 (en) * 2010-09-15 2012-03-15 Samsung Medison Co., Ltd. Ultrasound strain imaging based on lateral displacement compensation
US20120095333A1 (en) * 2010-10-13 2012-04-19 Jingfeng Jiang Coupled Axial and Lateral Displacement Estimation for Elasticity Imaging
US20120226158A1 (en) * 2011-03-04 2012-09-06 Greenleaf James F System and Method for Correcting Errors in Shear Wave Measurements Arising From Ultrasound Beam Geometry
CN102824194A (zh) * 2011-06-14 2012-12-19 深圳迈瑞生物医疗电子股份有限公司 一种弹性成像中的位移检测方法及装置
CN103735287A (zh) * 2013-12-05 2014-04-23 中国科学院苏州生物医学工程技术研究所 一种血管内超声弹性成像二维多级混合位移估计方法
US20140309531A1 (en) * 2012-10-01 2014-10-16 Olympus Medical Systems Corp. Ultrasonic observation apparatus, operation method of the same, and computer readable recording medium

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6508768B1 (en) * 2000-11-22 2003-01-21 University Of Kansas Medical Center Ultrasonic elasticity imaging
US20070083116A1 (en) * 2005-09-21 2007-04-12 Fuji Photo Film Co., Ltd. Ultrasonic diagnostic apparatus
US20070167772A1 (en) * 2005-12-09 2007-07-19 Aloka Co., Ltd. Apparatus and method for optimized search for displacement estimation in elasticity imaging
US20100256494A1 (en) * 2007-11-16 2010-10-07 Takashi Azuma Ultrasonic imaging system
US20120065505A1 (en) * 2010-09-15 2012-03-15 Samsung Medison Co., Ltd. Ultrasound strain imaging based on lateral displacement compensation
US20120095333A1 (en) * 2010-10-13 2012-04-19 Jingfeng Jiang Coupled Axial and Lateral Displacement Estimation for Elasticity Imaging
US20120226158A1 (en) * 2011-03-04 2012-09-06 Greenleaf James F System and Method for Correcting Errors in Shear Wave Measurements Arising From Ultrasound Beam Geometry
CN102824194A (zh) * 2011-06-14 2012-12-19 深圳迈瑞生物医疗电子股份有限公司 一种弹性成像中的位移检测方法及装置
US20140309531A1 (en) * 2012-10-01 2014-10-16 Olympus Medical Systems Corp. Ultrasonic observation apparatus, operation method of the same, and computer readable recording medium
CN103735287A (zh) * 2013-12-05 2014-04-23 中国科学院苏州生物医学工程技术研究所 一种血管内超声弹性成像二维多级混合位移估计方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113628183A (zh) * 2021-08-06 2021-11-09 青岛海信医疗设备股份有限公司 一种超声检测对象的容积确定方法及超声设备

Also Published As

Publication number Publication date
CN112674799B (zh) 2022-11-25

Similar Documents

Publication Publication Date Title
US7945117B2 (en) Methods and systems for registration of images
JP6985371B2 (ja) 検出精度を向上させるために関心領域の異なるビューからの複数の画像を用いるコンピュータ支援検出
US20190220977A1 (en) Cross-Domain Image Analysis and Cross-Domain Image Synthesis Using Deep Image-to-Image Networks and Adversarial Networks
US7565019B2 (en) Method of volume-panorama imaging processing
JP4917733B2 (ja) 尤度最大化を利用した画像位置合わせシステム及び方法
US8139710B2 (en) Systems, methods and apparatus to image objects
JP5432287B2 (ja) レポート機能を備えた画像システム及び動作方法
Jafari et al. Cardiac point-of-care to cart-based ultrasound translation using constrained CycleGAN
Virga et al. Use the force: deformation correction in robotic 3D ultrasound
CN112674799B (zh) 超声弹性成像方法、电子设备及存储介质
CN111223158B (zh) 心脏冠脉图像的伪影校正方法和可读存储介质
EP2428815A1 (en) 3D ultrasound system for extending view of image and method for operating the 3D ultrasound system
US20210035285A1 (en) Case-adaptive medical image quality assessment
CN112773402A (zh) 智能辅助引导方法、超声诊断设备及存储介质
US11484286B2 (en) Ultrasound evaluation of anatomical features
US20210342638A1 (en) X-ray image synthesis from ct images for training nodule detection systems
CN112967386A (zh) 生物力学建模方法、装置、电子设备及存储介质
US10881379B2 (en) Method of visualizing a sequence of ultrasound images, computer program product and ultrasound system
KZ Tehrani et al. Physically inspired constraint for unsupervised regularized ultrasound elastography
CN111803128A (zh) 一种乳腺组织弹性成像方法、装置、设备和介质
CN110477948B (zh) 弹性成像方法及装置、成像设备、存储介质
CN110930394B (zh) 测量肌肉肌纤维束线斜率和羽状角的方法及终端设备
US8326007B2 (en) Methods and apparatus for combined 4D presentation of quantitative regional measurements and morphology
CN113792736A (zh) 基于深度学习的医学图像处理方法、设备及介质
US10163529B2 (en) Display processing method and apparatus

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant