CN112790785B - 剪切波弹性成像方法及系统 - Google Patents

剪切波弹性成像方法及系统 Download PDF

Info

Publication number
CN112790785B
CN112790785B CN202011593839.2A CN202011593839A CN112790785B CN 112790785 B CN112790785 B CN 112790785B CN 202011593839 A CN202011593839 A CN 202011593839A CN 112790785 B CN112790785 B CN 112790785B
Authority
CN
China
Prior art keywords
sound pressure
array element
intensity
shear wave
wave elastography
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.)
Active
Application number
CN202011593839.2A
Other languages
English (en)
Other versions
CN112790785A (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.)
Chison Medical Technologies Co ltd
Original Assignee
Chison Medical Technologies 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 Chison Medical Technologies Co ltd filed Critical Chison Medical Technologies Co ltd
Priority to CN202011593839.2A priority Critical patent/CN112790785B/zh
Publication of CN112790785A publication Critical patent/CN112790785A/zh
Application granted granted Critical
Publication of CN112790785B publication Critical patent/CN112790785B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了一种剪切波弹性成像方法及系统，方法包括：获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；获取超声探头中的每个阵元在各个焦点位置处产生的声压h；获取声压权重w；根据所述各个焦点位置处的声压分量p、每个阵元在各个焦点位置产生的声压h及声压权重w，获取各个阵元的强度q；当各个阵元的强度q满足预设条件时，控制所述超声探头按照所述各个阵元的强度q进行剪切波弹性成像，本发明提供的方法及系统通过控制多个焦点位置处的超声强度分布均匀，在待测组织中生成质量好的剪切波弹性成像。

Description

剪切波弹性成像方法及系统
技术领域
本发明涉及超声扫查技术领域，具体涉及一种剪切波弹性成像方法及系统。
背景技术
最近的研究表明，剪切波弹性成像可以显著缩短进行慢性肝炎患者肝纤维化的肝活检，用于评估肝脏弹性的剪切波弹性成像技术基于内部机械推动产生的剪切波，这种内部推动是由超声系统所产生，并直接进入肝脏的发射波束，其脉冲长度通常需要100μs的持续时间才能在组织中产生如图1所示的剪切波，剪切波是一种机械波，可以从焦点沿横向方向传播产生。
由于剪切波在沿横向传播时会衰减，因此一个焦点只能覆盖一定宽度的区域，例如大约一厘米；如果ROI(目标区域)的宽度为几厘米宽，则将重复生成焦点以将焦点集中在ROI中的多个位置，每个位置都具有跟踪数据采集，因此为了生成剪切波弹性成像，通常需要在ROI中不同时间上依次生成多个单焦点，以覆盖整个横向区域。另一种方法是通过同时产生多个焦点而进行一次快速跟踪数据采集这种方法，将有助于减少由于一系列单焦点而可能产生的运动问题。
然而由于声场干扰用梳状焦点法产生的多个焦点可能并不相似，现有技术采用孔径分割方法，产生超声波来克服这种焦点强度不均匀的分布，也无法达到良好的效果。例如，图2显示了通过将整个孔径分为三个子孔径来产生三个焦点，每个子孔负责一个焦点，因为3个子孔径之间会相互干扰，因此即使每个子孔径都集中在一个焦点位置，这三个位置的最终焦点强度可能会有所不同，图3示出了将64个通道分成三组子孔径时在三个焦点位置的最终强度分布的模拟结果，在位置-9mm，0mm，9mm处的三个焦点，它们的深度z等于15mm，可以清楚地看到，中央的焦点几乎是两侧的两个焦点的两倍，即使发射子孔径，即每个单独焦点的有效阵元数量相同，产生的最终焦点的强度通常也不相同。
发明内容
因此，本发明要解决的技术问题在于克服现有技术中在剪切波弹性成像时，超声探头产生多个焦点位置处的强度分布不均的问题，从而提供一种剪切波弹性成像方法及系统。
为达到上述目的，本发明提供如下技术方案：
第一方面，本发明实施例提供一种剪切波弹性成像方法，包括如下步骤：
获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；
获取所述超声探头中的每个阵元在各个焦点位置处产生的声压h；
获取声压权重w；
根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q；
当各个阵元的强度q满足预设条件时，控制所述超声探头按照所述各个阵元的强度q进行剪切波弹性成像。
在一实施例中，所述根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q的步骤，包括：
构建由各个焦点位置处的声压分量p组成的第一声压矩阵P、每个阵元在各个焦点位置处产生的声压h组成的第二声压矩阵H，以及由各个声压权重w组成的声压权重矩阵W，则由各个阵元的强度q组成的源强度矩阵Q通过以下公式获取：
Q＝WH*T(HWH*T)-1P。
在一实施例中，当源强度q满足不满足预设条件时，更新声压权重w，利用更新后的声压权重w更新所述源强度q，直至满足所述预设条件。
在一实施例中，所述预设条件包括：源强度矩阵Q的标准偏差值小于预设阈值；
当源强度矩阵Q的标准偏差值大于预设阈值时，更新权重声压权重矩阵W，利用更新后的声压权重矩阵W更新所述源强度矩阵Q，直至当源强度矩阵Q的标准偏差值小于预设阈值或迭代次数大于预设迭代截止次数值时，停止迭代。
在一实施例中，各个焦点位置处的声压分量p，通过以下公式获取：
Figure BDA0002869385540000021
其中，
Figure BDA0002869385540000022
ρ0是组织的密度，c是组织中声音的速度，k是波数，λ是波长，qn是阵元n的源强度，L是阵元的宽度，r是从阵元中心到场点的距离，θmn是从阵元n到场点m的角度，pm是超声第m个焦点的声压，包括幅度和相位。
在一实施例中，每个阵元以单位声强度在各个焦点位置产生的声压h，通过以下公式获取：
Figure BDA0002869385540000023
其中，hmn是第n个阵元以单位声强度在第m个位置产生的声压。
在一实施例中，所述声压权重w，通过以下公式表示：
Figure BDA0002869385540000024
其中，v是迭代次数，N是阵元的数量。
在一实施例中，所述源强度矩阵Q的标准偏差值通过以下公式计算：
Figure BDA0002869385540000031
Figure BDA0002869385540000032
第二方面，本发明实施例提供一种剪切波弹性成像系统，包括：
焦点声压获取模块，用于获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；
阵元声压获取模块，用于获取所述超声探头中的每个阵元在各个焦点位置处产生的声压h；
声压权重获取模块，用于获取声压权重w；
阵元的强度获取模块，用于根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q；
阵元的强度更新模块，用于当个阵元的强度q满足预设条件时，所述超声探头按照当前的各个阵元的强度q发射超声波进行剪切波弹性成像。
第三方面，本发明实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行本发明实施例第一方面的剪切波弹性成像方法。
第四方面，本发明实施例提供一种计算机设备，包括：存储器和处理器，所述存储器和所述处理器之间互相通信连接，所述存储器存储有计算机指令，所述处理器通过执行所述计算机指令，从而执行本发明实施例第一方面的剪切波弹性成像方法。
本发明技术方案，具有如下优点：
1、本发明提供的剪切波弹性成像方法及系统，首先获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；获取所述超声探头中的每个阵元在各个焦点位置处产生的声压h；获取声压权重w；根据各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q；当各个阵元的强度q满足预设条件时，控制所述超声探头按照各个阵元的强度q进行剪切波弹性成像，本发明通过控制多个焦点位置处的强度分布均匀，在待测组织中生成质量好的剪切波弹性成像。
2.本发明提供的剪切波弹性成像方法及系统，当源强度矩阵Q的标准偏差值大于预设阈值时，更新权重声压权重矩阵W，利用更新后的声压权重矩阵W更新所述源强度矩阵Q，直至当源强度矩阵Q的标准偏差值小于预设阈值或迭代次数大于预设迭代截止次数值时停止迭代，通过这种迭代加权的方式控制每个阵元的强度相同，使得最终生成的多个焦点位置处的强度分布均匀。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为剪切波弹性生成过程的示意图；
图2为现有技术中提供的采用孔径分割方法产生的焦点力示意图；
图3为采用孔径分割方法产生的焦点强度分布示意图；
图4为本发明实施例中提供的剪切波弹性成像方法的一个示例的流程图；
图5为一个阵元在远场中的声压幅度的坐标图；
图6为本发明实施例提供的迭代加权处理方式的迭代流程图；
图7为本发明实施例提供的通过迭代加权处理方式得到的焦点强度分布示意图；
图8为本发明实施例中提供的剪切波弹性成像系统的一个具体示例的模块组成图；
图9为本发明实施例提供的计算机设备一个具体示例的组成图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
实施例1
声波在组织中的传播比剪切波传播的快得多(约1000倍)，因此可以完全跟踪剪切波在组织中沿横向方向的传播，通过测量感兴趣区域的剪切波速度，可以提供组织的二维定量弹性图。为此需要非常高的运动检测帧频率，因为剪切波速度通常约为每秒几米，为了获得剪切波弹性成像，超声系统需要将跟踪信号发送到组织并接收反向散射的回波，可以传输动态激励以在体内产生剪切波，检测剪切波传播速度的方法是使用相同的成像换能器将跟踪脉冲重复发送到人体，并接收反射信号以监视组织位移。
在实际应用中，超声设备的超声探头阵列发射的超声波可以产生聚集于不同位置的单个推动力，可以同时产生多个推动力，即会产生一个焦点或多个焦点，本发明实施例主要针对的是产生多个焦点时的情况，提供一种剪切波弹性成像方法，如图4所示，该方法包括如下步骤：
步骤S1：获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个。
实际应用中根据期望剪切波弹性成像的目标区域大小，并以此来确定超声设备所需的焦点数、焦点位置和声压分布，由于单个焦点的扫描范围有限，通常需要同时产生多个焦点而进行一次快速跟踪数据。
在一具体实施例中，期望的剪切波弹性成像为二维剪切波弹性成像，获取超声弹性检测设备的超声探头阵列中N个阵元同时参与产生的M(M不小于2)个焦点的声压，获取超声探头发射的超声波在各个焦点位置处的声压分量p。
超声探头阵列中一个阵元在远场中的声压幅度如图5所示，声压幅度通过以下公式表示：
Figure BDA0002869385540000051
其中，
Figure BDA0002869385540000052
ρ0是组织的密度，c是组织中声音的速度，k是波数，q是阵元的源强度，L是阵元的宽度，ω是角频率，r是从阵元中心到场点(被测组织中的位置点)的距离，θ是从阵元到场点的角度。
假设有N个阵元同时参与产生M焦点，通常N>>M，对于第m∈(1,M)个焦点的声压通过pm以下公式表示：
Figure BDA0002869385540000053
其中，
Figure BDA0002869385540000061
ρ0是组织的密度，c是组织中声音的速度，k是波数，λ是波长，qn是阵元n∈(1,N)的源强度，L是阵元的宽度，r是从阵元中心到场点的距离，θmn是从阵元n到场点的角度，pm是超声第m个焦点的声压，包括幅度和相位。
因此，由各个焦点声压分量p组成的第一声压矩阵P表示为：
Figure BDA0002869385540000062
步骤S2：获取所述超声探头中的每个阵元在各个焦点位置处产生的声压H。
本发明实施例中，第n个阵元在第m个焦点产生的声压h，通过以下公式计算，通过以下公式获取：
Figure BDA0002869385540000063
由每个阵元在各个焦点位置处产生的声压h组成的第二声压矩阵H表示为：
Figure BDA0002869385540000064
步骤S3：获取声压权重w。
本发明实施例中的声压权重w，通过以下公式表示：
Figure BDA0002869385540000065
由各个声压权重w组成的声压权重矩阵W。
步骤S4：根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q。
构建由每个阵元的声强度q构成的源强度矩阵Q：
Figure BDA0002869385540000066
通过第一声压矩阵P、第二声压矩阵H以及声压权重矩阵W，获取源强度矩阵Q，具体通过以下公式获得：
Q＝wH*T(HWH*T)-1P (8)
步骤S5：当各个阵元的强度q满足预设条件时，控制所述超声探头按照所述各个阵元的强度q进行剪切波弹性成像。
由于计算出的源强度矩阵Q可能对每个阵元都不具有相同的强度，因此本发明实施例通过迭代加权的处理方式以使阵元的强度均匀分布。本发明的迭代停止条件，包括：当源强度矩阵Q的标准偏差值大于预设阈值时，更新权重声压权重矩阵W，利用更新后的声压权重矩阵W更新所述源强度矩阵Q，直至当源强度矩阵Q的标准偏差值小于预设阈值σ或迭代次数大于预设迭代截止次数值V时，停止迭代，整个迭代流程如图6所示。其中源强度矩阵Q的标准偏差值，通过以下公式计算：
Figure BDA0002869385540000071
Figure BDA0002869385540000072
在一具体实施例中，在迭代之前设置初始的第一声压矩阵P和初始的声压权重矩阵W均为单位矩阵，即W＝I，P＝[1…1]T，仅作为举例，不以此为限，在实际应用中根据具体需求做合理设置初始值。在迭代停止之后得到最终源强度矩阵Q中的元素为每个阵元的源强度值，控制每个阵元按照最终源强度矩阵Q确定的最终值发射信号，相位按照计算得到的相位发射。
基于图3中示例相同的超声系统设备配置，采用的本发明实施提供的方法后，如图7所示，在三个焦点位置的最终强度分布的模拟结果，可以清楚地看到，在中间焦点位置处的强度与在其他两个焦点位置处的强度非常相似，达到了焦点强度均匀分布的效果，从而使得最终的生成剪切波弹性成像质量更好。
实施例2
本发明实施例提供一种剪切波弹性成像系统，如图8所示，包括：
焦点声压获取模块1，用于获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；此模块执行实施例1中的步骤S1所描述的方法，在此不再赘述。
阵元声压获取模块2，用于获取所述超声探头中的每个阵元在各个焦点位置处产生的声压h；此模块执行实施例1中的步骤S2所描述的方法，在此不再赘述。
声压权重获取模块3，用于获取声压权重w；模块执行实施例1中的步骤S3所描述的方法，在此不再赘述。
阵元的强度获取模块4，用于根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q；此模块执行实施例1中的步骤S4所描述的方法，在此不再赘述。
阵元的强度q更新模块5，用于当个阵元的强度q满足预设条件时，所述超声探头按照当前的各个阵元的强度q发射超声波进行剪切波弹性成像。此模块执行实施例1中的步骤S5所描述的方法，在此不再赘述。
本发明实施例提供的剪切波弹性成像系统，通过迭代加权的处理方式，将达到预设迭代停止条件时的源强度矩阵，作为最终源强度发射超声信号，进行使得多个焦点的声压强度分布均匀，超声探头阵列根据最终阵元强度，在待测组织中发出超声波生成质量好的剪切波弹性成像。
实施例3
本发明实施例提供一种计算机设备，如图9所示，该设备可以包括处理器51和存储器52，其中处理器51和存储器52可以通过总线或者其他方式连接，图9以通过总线连接为例。
存储器52作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块，如本发明实施例中的对应的程序指令/模块。处理器51通过运行存储在存储器52中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及数据处理，即实现上述方法实施例1中的剪切波弹性成像方法。
存储器52可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器51所创建的数据等。此外，存储器52可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器52可选包括相对于处理器51远程设置的存储器，这些远程存储器可以通过网络连接至处理器51。上述网络的实例包括但不限于互联网、企业内部网、企业内网、移动通信网及其组合。
一个或者多个模块存储在存储器52中，当被处理器51执行时，执行实施例1中的剪切波弹性成像方法。
上述计算机设备具体细节可以对应参阅实施例1中对应的相关描述和效果进行理解，此处不再赘述。
本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)、随机存储记忆体(Random Access Memory，RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive，缩写：HDD)或固态硬盘(Solid-StateDrive，SSD)等；存储介质还可以包括上述种类的存储器的组合。
显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (11)

1.一种剪切波弹性成像方法，其特征在于，包括如下步骤：
获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；
获取所述超声探头中的每个阵元在各个焦点位置处产生的声压h；
获取声压权重w；
根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q；
当各个阵元的强度q满足预设条件时，控制所述超声探头按照所述各个阵元的强度q进行剪切波弹性成像。
2.根据权利要求1所述的剪切波弹性成像方法，其特征在于，所述根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q的步骤，包括：
构建由各个焦点位置处的声压分量p组成的第一声压矩阵P、每个阵元在各个焦点位置处产生的声压h组成的第二声压矩阵H，以及由各个声压权重w组成的声压权重矩阵W，则由各个阵元的强度q组成的源强度矩阵Q通过以下公式获取：
Q＝WH*T(HWH*T)-1P。
3.根据权利要求2所述的剪切波弹性成像方法，其特征在于，当源强度q满足不满足预设条件时，更新声压权重w，利用更新后的声压权重w更新所述源强度q，直至满足所述预设条件。
4.根据权利要求3所述的剪切波弹性成像方法，其特征在于，所述预设条件包括：源强度矩阵Q的标准偏差值小于预设阈值；
当源强度矩阵Q的标准偏差值大于预设阈值时，更新权重声压权重矩阵W，利用更新后的声压权重矩阵W更新所述源强度矩阵Q，直至当源强度矩阵Q的标准偏差值小于预设阈值或迭代次数大于预设迭代截止次数值时，停止迭代。
5.根据权利要求1所述的剪切波弹性成像方法，其特征在于，各个焦点位置处的声压分量p，通过以下公式获取：
Figure FDA0002869385530000011
其中，
Figure FDA0002869385530000012
ρ0是组织的密度，c是组织中声音的速度，k是波数，λ是波长，qn是阵元n的源强度，L是阵元的宽度，r是从阵元中心到场点的距离，θmn是从阵元n到场点m的角度，pm是超声第m个焦点的声压，包括幅度和相位。
6.根据权利要求5所述的剪切波弹性成像方法，其特征在于，每个阵元以单位声强度在各个焦点位置产生的声压h，通过以下公式获取：
Figure FDA0002869385530000013
其中，hmn是第n个阵元以单位声强度在第m个位置产生的声压。
7.根据权利要求5所述的剪切波弹性成像方法，其特征在于，所述声压权重w，通过以下公式表示：
Figure FDA0002869385530000021
其中，v是迭代次数，N是阵元的数量。
8.根据权利要求4所述的剪切波弹性成像方法，其特征在于，所述源强度矩阵Q的标准偏差值通过以下公式计算：
Figure FDA0002869385530000022
Figure FDA0002869385530000023
9.一种剪切波弹性成像系统，其特征在于，包括：
焦点声压获取模块，用于获取超声探头发射的超声波在各个焦点位置处的声压分量p，其中根据期望剪切波弹性成像的目标区域大小来决定所需要的焦点数、焦点位置和声压分布，所述焦点包括至少两个；
阵元声压获取模块，用于获取所述超声探头中的每个阵元在各个焦点位置处产生的声压h；
声压权重获取模块，用于获取声压权重w；
阵元的强度获取模块，用于根据所述各个焦点位置处的声压分量p、所述每个阵元在各个焦点位置产生的声压h及所述声压权重w，获取各个阵元的强度q；
阵元的强度更新模块，用于当个阵元的强度q满足预设条件时，所述超声探头按照当前的各个阵元的强度q发射超声波进行剪切波弹性成像。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机指令，所述计算机指令用于使所述计算机执行如权利要求1-8任一项所述的剪切波弹性成像方法。
11.一种计算机设备，其特征在于，包括：存储器和处理器，所述存储器和所述处理器之间互相通信连接，所述存储器存储有计算机指令，所述处理器通过执行所述计算机指令，从而执行如权利要求1-8任一项所述的剪切波弹性成像方法。
CN202011593839.2A 2020-12-29 2020-12-29 剪切波弹性成像方法及系统 Active CN112790785B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011593839.2A CN112790785B (zh) 2020-12-29 2020-12-29 剪切波弹性成像方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011593839.2A CN112790785B (zh) 2020-12-29 2020-12-29 剪切波弹性成像方法及系统

Publications (2)

Publication Number Publication Date
CN112790785A CN112790785A (zh) 2021-05-14
CN112790785B true CN112790785B (zh) 2022-08-05

Family

ID=75805555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011593839.2A Active CN112790785B (zh) 2020-12-29 2020-12-29 剪切波弹性成像方法及系统

Country Status (1)

Country Link
CN (1) CN112790785B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022141052A1 (zh) * 2020-12-29 2022-07-07 无锡祥生医疗科技股份有限公司 剪切波弹性成像方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6135971A (en) * 1995-11-09 2000-10-24 Brigham And Women's Hospital Apparatus for deposition of ultrasound energy in body tissue
CN101108268A (zh) * 2007-08-23 2008-01-23 上海交通大学 相控阵聚焦超声多模式热场形成方法
CN109069115A (zh) * 2017-06-06 2018-12-21 深圳迈瑞生物医疗电子股份有限公司 一种在超声扫描中成像的方法、装置及系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6135971A (en) * 1995-11-09 2000-10-24 Brigham And Women's Hospital Apparatus for deposition of ultrasound energy in body tissue
CN101108268A (zh) * 2007-08-23 2008-01-23 上海交通大学 相控阵聚焦超声多模式热场形成方法
CN109069115A (zh) * 2017-06-06 2018-12-21 深圳迈瑞生物医疗电子股份有限公司 一种在超声扫描中成像的方法、装置及系统

Also Published As

Publication number Publication date
CN112790785A (zh) 2021-05-14

Similar Documents

Publication Publication Date Title
JP6257096B2 (ja) 位相コヒーレンスによるグレーティングローブ抑圧のためのビーム形成法
US8313436B2 (en) Methods and apparatus for ultrasound imaging
EP3581961A1 (en) Method and apparatus for ultrasound imaging with improved beamforming
Nikolov et al. Virtual ultrasound sources in high-resolution ultrasound imaging
Rasmussen et al. 3-D ultrasound imaging performance of a row-column addressed 2-D array transducer: A measurement study
US8001843B2 (en) Circular arc wide beam transmission method and apparatus for ultrasonic imaging
JP2016019849A (ja) 超音波合成イメージングの装置と方法
JP2002523203A (ja) 高調波イメージングを用いて疎なアレイのサイドローブ性能を向上させる方法及び装置
CN103251429A (zh) 超声波成像装置
JP6014643B2 (ja) 超音波診断装置
US9883851B2 (en) System and method for shear wave generation with steered ultrasound push beams
Londhe et al. Superharmonic imaging for medical ultrasound: A review
CN106940883B (zh) 基于超声系统点扩散函数仿真和压缩感知的超声成像方法
CN112790785B (zh) 剪切波弹性成像方法及系统
US20170276775A1 (en) Acoustic imaging method and device
CN111580112B (zh) 基于平面波的水下声呐传感器阵列成像方法
JPWO2016132924A1 (ja) 超音波撮像装置および超音波信号の処理方法
CN112998745A (zh) 一种用于超声成像的发射波束形成方法、系统及诊断设备
WO2022141052A1 (zh) 剪切波弹性成像方法及系统
Peretz et al. Deep learning applied to beamforming in synthetic aperture ultrasound
Orlowska et al. The effect of different coherence-based beamforming techniques on the accuracy of high frame rate speckle tracking echocardiography
KR101551469B1 (ko) 평면파를 이용한 초음파 송수신집속 방법 및 시스템
JP7211150B2 (ja) 超音波診断装置、超音波画像生成方法及びプログラム
Boni et al. Prototype 3D real-time imaging system based on a sparse PZT spiral array
Kortbek et al. P2b-1 synthetic aperture focusing applied to imaging using a rotating single element transducer

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