CN114668419A - 超声弹性成像方法、装置、设备及存储介质 - Google Patents

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

Info

Publication number
CN114668419A
CN114668419A CN202210234626.3A CN202210234626A CN114668419A CN 114668419 A CN114668419 A CN 114668419A CN 202210234626 A CN202210234626 A CN 202210234626A CN 114668419 A CN114668419 A CN 114668419A
Authority
CN
China
Prior art keywords
echo signal
sparse
displacement value
deformation displacement
deformation
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
CN202210234626.3A
Other languages
English (en)
Other versions
CN114668419B (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.)
Yichao Technology Wuhan Co ltd
Original Assignee
Yichao Medical Technology Wuhan 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 Yichao Medical Technology Wuhan Co ltd filed Critical Yichao Medical Technology Wuhan Co ltd
Priority to CN202210234626.3A priority Critical patent/CN114668419B/zh
Publication of CN114668419A publication Critical patent/CN114668419A/zh
Application granted granted Critical
Publication of CN114668419B publication Critical patent/CN114668419B/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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
    • 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
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties

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

本申请适用于超声成像领域，提供了一种超声弹性成像方法、装置、设备及存储介质，通过获取目标组织区域形变前与形变后的相邻稀疏信号，对第一稀疏回波信号和第二稀疏回波信号进行相位差估算，并根据目标相位差计算第一形变位移值，对第一稀疏回波信号和第二稀疏回波信号进行求模运算，确定线性图像数据，并根据线性图像数据计算第二形变位移值，对第一形变位移值和第二形变位移值进行集合运算，得到最终形变位移值，并根据最终形变位移值生成目标组织区域的目标应变弹性图像。本申请同时考虑横向位移和纵向位移，从而在最小化数据量的基础上提高应变力的计算精度，实现实时、高精度的超声弹性成像。

Description

超声弹性成像方法、装置、设备及存储介质
技术领域
本申请涉及超声成像领域，尤其涉及一种超声弹性成像方法、装置、设备及存储介质。
背景技术
在医疗彩超设备系统中，超声弹性成像是利用声波检测组织的硬度属性，将超声探头置于人体表面，通过有规律的施加压力，使得人体组织在外加压力的情况发生形变。当组织受到压力时，组织内将产生一个沿着探头的纵向压力方向的应变，若组织内部弹性模量分布不均匀，组织内的应变分布也会所有差异。弹性模量较大，则应变较小，形变也较小，即组织移动较小；弹性模量较小，则应变较大，形变也较大，即组织移动较大。
通过外部施加压力产生硬度图的方法为准静态弹性成像方法，其容易受到外部施压者的主观意愿影响，并且只考虑沿着纵向压力方向计算组织硬度信息，而实际组织在发生形变时会产生横向的位移信息，从而导致最终硬度测量不准。
发明内容
本申请提供了一种超声弹性成像方法、装置、设备及存储介质，以解决当前超声弹性成像结果的准确度低的技术问题。
为了解决上述技术问题，第一方面，本申请提供一种超声弹性成像方法，包括：
获取目标组织区域形变前与形变后的相邻稀疏信号，相邻稀疏信号包括第一稀疏回波信号和与第一稀疏回波信号相邻的第二稀疏回波信号，第一稀疏回波信号为目标组织区域形变前的稀疏信号，第二稀疏回波信号为目标组织区域形变后的稀疏信号；
对第一稀疏回波信号和第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值；
对第一稀疏回波信号和第二稀疏回波信号进行求模运算，确定第一稀疏回波信号和第二稀疏回波信号的线性图像数据，并根据线性图像数据计算第二形变位移值；
对第一形变位移值和第二形变位移值进行集合运算，得到目标组织区域的最终形变位移值，并根据最终形变位移值生成目标组织区域的目标应变弹性图像。
本申请通过获取目标组织区域形变前与形变后的相邻稀疏信号作为超声成像的运算数据，以最小化数据量；再进行两个阶段的实时应变力计算，即对第一稀疏回波信号和第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值，以及对第一稀疏回波信号和第二稀疏回波信号进行求模运算，确定第一稀疏回波信号和第二稀疏回波信号的线性图像数据，并根据线性图像数据计算第二形变位移值，以提高超声成像的实时性以及同时考虑横向位移和纵向位移，从而在最小化数据量的基础上提高应变力的计算精度；最后对第一形变位移值和第二形变位移值进行集合运算，得到目标组织区域的最终形变位移值，并根据最终形变位移值生成目标组织区域的目标应变弹性图像，从而实现实时、高精度的超声弹性成像。
作为优选，获取目标组织区域形变前与形变后的相邻稀疏信号，包括：
获取目标组织区域形变前与形变后的多帧原始回波信号；
基于预设时间间隔，提取多帧原始回波信号中相邻的两帧原始回波信号；
对相邻的两帧原始回波信号进行下采样，得到第一稀疏回波信号和第二稀疏回波信号。
本申请通过信号提取和下采样，以最小化运算数据的数据量，提高后续运算的运算效率和实时性。
作为优选，对第一稀疏回波信号和第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值，包括：
基于预设的零相位迭代算法，对第一稀疏回波信号和第二稀疏回波信号在同一位置纵向点上进行相位差估算，得到第一稀疏回波信号与第二稀疏回波信号之间的第一相位差；
基于零相位迭代算法，对第一稀疏回波信号和第二稀疏回波信号在与纵向点相邻的多个横向点进行相位差估算，得到多个横向点对应的第二相位差；
对多个第二相位差进行均值运算，得到第一形变位移值。
本申请通过纵向相位差估算和横向相位差估算，以同时考虑组织的纵向应力形变和横向应力形变，提高形变位移的精度。
作为优选，零相位迭代算法的计算公式为：
Figure BDA0003534475720000031
其中，Dpre(i+u,p)为上一帧信号相对于当前帧信号的形变位移值，Z1pre(i+u+k,p+h)为第一稀疏回波信号，Z1pos(i+u+k,p+h)为第二稀疏回波信号，i为横向方向的信号线号，p为纵向方向的信号点号，u为横向内的计算次数，j为虚数因子，K为Mask窗口的横向尺寸，H为Mask窗口的纵向尺寸，w0为角频率，conj为共轭复数。
作为优选，对第一稀疏回波信号和第二稀疏回波信号进行求模运算，确定第一稀疏回波信号和第二稀疏回波信号的线性图像数据，并根据线性图像数据计算第二形变位移值，包括：
对第一稀疏回波信号和第二稀疏回波信号进行求模运算，得到第一稀疏回波信号的第一高位宽模值和第二稀疏回波信号的第二高位宽模值；
基于预设的数据动态范围，对第一高位宽模值和第二高位宽模值进行压缩到低位宽，得到第一稀疏回波信号对应的第一线性图像数据和第二稀疏回波信号对应的第二线性图像数据；
根据第一线性图像数据和第二线性图像数据，计算第二形变位移值。
本申请通过求模运算和数据动态范围压缩，以考虑相邻稀疏信号作为复数信号的线性图像位移，为第二形变位移值的计算提供数据源并进一步提高数据的计算效率，同时提高后续成像的应变弹性图像的显示精度。
作为优选，根据第一线性图像数据和第二线性图像数据，计算第二形变位移值，包括：
基于预设光流法，对第一线性图像数据和第二线性图像数据进行运算，得到第二形变位移值。
作为优选，对第一形变位移值和第二形变位移值进行集合运算，得到目标组织区域的最终形变位移值，并根据最终形变位移值生成目标组织区域的目标应变弹性图像，包括：
对第一形变位移值和第二形变位移值进行矩阵求和运算，得到最终形变位移值；
根据最终形变位移值，生成初始应变弹性图像；
对初始应变弹性图像进行纵向差分运算，并对纵向差分运算后的初始应变弹性图像进行归一化，映射得到目标应变弹性图像。
第二方面，本申请还提供一种超声弹性成像装置，包括：
获取模块，用于获取目标组织区域形变前与形变后的相邻稀疏信号，相邻稀疏信号包括第一稀疏回波信号和与第一稀疏回波信号相邻的第二稀疏回波信号，第一稀疏回波信号为目标组织区域形变前的稀疏信号，第二稀疏回波信号为目标组织区域形变后的稀疏信号；
第一计算模块，用于对第一稀疏回波信号和第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值；
第二计算模块，用于对第一稀疏回波信号和第二稀疏回波信号进行求模运算，确定第一稀疏回波信号和第二稀疏回波信号的线性图像数据，并根据线性图像数据计算第二形变位移值；
成像模块，用于对第一形变位移值和第二形变位移值进行集合运算，得到目标组织区域的最终形变位移值，并根据最终形变位移值生成目标组织区域的应变弹性图像。
第三方面，本申请还提供一种计算机设备，包括处理器和存储器，存储器用于存储计算机程序，计算机程序被处理器执行时实现如第一方面的超声弹性成像方法。
第四方面，本申请还提供一种计算机可读存储介质，其存储有计算机程序，计算机程序被处理器执行时实现如第一方面的超声弹性成像方法。
需要说明的是，上述第二方面至第四方面的有益效果请参见上述第一方面的相关描述，在此不再赘述。
附图说明
图1为本申请实施例示出的超声弹性成像方法的流程示意图；
图2为本申请实施例示出的原始回波信号的示意图；
图3为本申请实施例示出的相邻原始回波信号的提取示意图；
图4为本申请实施例示出的下采样效果示意图；
图5为本申请实施例示出的第一形变位移值的计算示意图；
图6为本申请实施例示出的装置的结构示意图；
图7为本申请实施例示出的计算机设备的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
如相关技术记载，通过外部施加压力产生硬度图的方法为准静态弹性成像方法，其容易受到外部施压者的主观意愿影响，并且只考虑沿着纵向压力方向计算组织硬度信息，而实际组织在发生形变时会产生横向的位移信息，从而导致最终硬度测量不准。
为此，本申请实施例提供一种超声弹性成像方法，通过获取目标组织区域形变前与形变后的相邻稀疏信号作为超声成像的运算数据，以最小化数据量；再进行两个阶段的实时应变力计算，即对第一稀疏回波信号和第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值，以及对第一稀疏回波信号和第二稀疏回波信号进行求模运算，确定第一稀疏回波信号和第二稀疏回波信号的线性图像数据，并根据线性图像数据计算第二形变位移值，以提高超声成像的实时性以及同时考虑横向位移和纵向位移，从而在最小化数据量的基础上提高应变力的计算精度；最后对第一形变位移值和第二形变位移值进行集合运算，得到目标组织区域的最终形变位移值，并根据最终形变位移值生成目标组织区域的目标应变弹性图像，从而实现实时、高精度的超声弹性成像。
请参照图1，图1为本申请实施例提供的一种超声弹性成像方法的流程示意图。本申请实施例的超声弹性成像方法可应用于计算机设备，该计算机设备与超声设备通信连接，以控制超声设备的工作状态等参数。计算机设备包括但不限于笔记本电脑、平板电脑、桌上型计算机、物理服务器和云服务器等设备。如图1所示，本实施例的超声弹性成像方法包括步骤S101至步骤S104，详述如下：
步骤S101，获取目标组织区域形变前与形变后的相邻稀疏信号，所述相邻稀疏信号包括所述第一稀疏回波信号和与所述第一稀疏回波信号相邻的第二稀疏回波信号，所述第一稀疏回波信号为所述目标组织区域形变前的稀疏信号，所述第二稀疏回波信号为所述目标组织区域形变后的稀疏信号。
在本步骤中，相邻稀疏信号是对相邻两帧原始回波信号(一帧为形变前的回波信号，另一帧为形变后的回波信号)进行下采样得到的两帧稀疏信号，其中原始回波信号是超声设备的探头扫描人体组织时发射超声信号后接收到的回波信号，探头可以是线阵、凸阵、相控阵等探头。
步骤S102，对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值。
在本步骤中，纵向点的相位差估算包括纵向点在纵向方向的相位差估算和在横向方向的相位差估算，相位差估算为前一帧信号相对于当前帧信号的形变位移估算。可选地，可利用零相位迭代算法实现相位差估算。
步骤S103，对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值。
在本步骤中，第一稀疏回波信号和第二稀疏回波信号为复数信号，通过求模运算确定线性图像数据，以考虑图像显示的形变位移值。
步骤S104，对所述第一形变位移值和所述第二形变位移值进行集合运算，得到所述目标组织区域的最终形变位移值，并根据所述最终形变位移值生成所述目标组织区域的目标应变弹性图像。
在本步骤中，每两帧相邻稀疏信号得到一个第一形变位移值和一个第二形变位移值，则多帧相邻稀疏信号对应有多个第一形变位移值和多个第二形变位移值，多个第一形变位移值组成第一矩阵，多个第二形变位移值组成第二矩阵，则集合运算为对第一矩阵和第二矩阵进行求和。
在一实施例中，在图1所示实施例的基础上，上述步骤S101，包括：
获取所述目标组织区域形变前与形变后的多帧原始回波信号；
基于预设时间间隔，提取多帧所述原始回波信号中相邻的两帧原始回波信号；
对相邻的两帧所述原始回波信号进行下采样，得到所述第一稀疏回波信号和所述第二稀疏回波信号。
在本实施例中，如图2所示的原始回波信号，获取目标组织区域内的多帧原始回波信号，记为Z(i，p，t)＝I(i，p，t)+Q(i，p，t)×j，其中i为横向方向的信号线号，p为纵向方向的信号点号，t为时间方向帧号，I表示复数信号的实部，Q表示复数信号的虚部，多帧原始回波信号包括组织形变前的原始信号帧和组织形变后的原始信号帧。可以理解的是，形变前的原始回波信号为未施加外力而得到的回波信号，形变后的原始回波信号为外力直接作用于人体组织得到的回波信号。
可选地，如图3A所示的交叉扫描方式，通过以预设时间间隔按照交叉扫描方式提取形变前和形变后的相邻两帧原始回波信号，分别记为Zpre(i，p)＝Ipre(i，p)+Qpre(i，p)×j和Zpos(i，p)＝Ipos(i，p)+Qpos(i，p)×j。假设用户按压探头的频率为f，形变前和形变后的时间间隔可进行调节，设定为n，每条扫描线耗费的时间为PRT，则时间间隔计算公式为：
Figure BDA0003534475720000081
可选地，若时间间隔n代表一帧中扫描的线数，则n值可覆盖目标组织区域的全部线数，如图3B的常规扫描方式，将线号34至36组成一帧形变前的原始回波信号，将线号为35至37组成一帧形变后的原始回波信号；n值也可覆盖部分线数，如图3A的交叉扫描方式，将线号30至31生成前一帧信号，再重复一次30至31生成后一帧信号，剩余的信号线遵从交叉扫描方式，即从32至33生成前一帧信号，再重复一次32至33生成后一帧信号。
可选地，如图4所示的稀疏信号提取示意图，41为原始回波信号，42为稀疏信号，下采样的的范围大小可以根据目标组织区域的大小进行动态配置。可选地，图4所示的横向线不进行下采样，纵向点的下采样率设为M，下采样为低通滤波，低通滤波采用窗函数设计的FIR滤波器，计算公式为：
Figure BDA0003534475720000082
Figure BDA0003534475720000083
Figure BDA0003534475720000084
其中，L为滤波器的阶数，fc为滤波器的归一化截止频率，Z1pre(i,p)为第一稀疏回波信号，Zpre(i,M×p+k)为形变前的原始回波信号，Z1pos(i,p)为第二稀疏回波信号，Zpos(i,M×p+k)为形变后的原始回波信号。
在一实施例中，在图1所示实施例的基础上，上述步骤S102，包括：
基于预设的零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在同一位置纵向点上进行相位差估算，得到所述第一稀疏回波信号与所述第二稀疏回波信号之间的第一相位差；
基于所述零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在与所述纵向点相邻的多个横向点进行相位差估算，得到多个横向点对应的第二相位差；
对多个所述第二相位差进行均值运算，得到所述第一形变位移值。
在本实施例中，如图5所示，采用零相位迭代算法在相邻两帧信号的同一位置纵向点上进行相位差计算，相位差设定为D(i,p)。可选地，所述零相位迭代算法的计算公式为：
Figure BDA0003534475720000091
其中，Dpre(i+u,p)为上一帧信号相对于当前帧信号的形变位移值，Z1pre(i+u+k,p+h)为所述第一稀疏回波信号，Z1pos(i+u+k,p+h)为所述第二稀疏回波信号，i为横向方向的信号线号，p为纵向方向的信号点号(即纵向点)，u为横向内的计算次数，j为虚数因子，K为Mask窗口的横向尺寸，H为Mask窗口的纵向尺寸，w0为角频率，conj为共轭复数。
在计算完当前纵向点的相位差后，再计算与当前纵向点在同一横向位置相邻多点的相位差。最后，将所有横向位置的相位差进行求平均，得到第一位移值：
Figure BDA0003534475720000092
在一实施例中，在图1所示实施例的基础上，上述步骤S103，包括：
对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，得到所述第一稀疏回波信号的第一高位宽模值和所述第二稀疏回波信号的第二高位宽模值；
基于预设的数据动态范围，对所述第一高位宽模值和所述第二高位宽模值进行压缩到低位宽，得到所述第一稀疏回波信号对应的第一线性图像数据和所述第二稀疏回波信号对应的第二线性图像数据；
根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值。
在本实施例中，对第一稀疏回波信号Z1pre(i，p)和第二稀疏回波信号Z1pos(i，p)进行求模运算，计算公式为：
Figure BDA0003534475720000101
Figure BDA0003534475720000102
再对高位宽模值Modulepre和Modulepos进行数据动态范围压缩到低位宽，得到相邻两帧相邻稀疏信号的线性图像数据，计算公式为：
Imgpre(i,p)＝int{Dr×20log10[Modulepre(i,p)]+0.5}∈[0,255]；
Imgpos(i,p)＝int{Dr×20log10[Modulepre(i,p)]+0.5}∈[0,255]；
其中Dr为斜率参数，用于调整图像的显示范围；log10为以10为底的对数运算。
可选地，所述根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值，包括：
基于预设光流法，对所述第一线性图像数据和所述第二线性图像数据进行运算，得到所述第二形变位移值。
在可选实施例中，采用机器视觉中的光流法，根据线性图像数据Imgpre(i,p)和Imgpos(i,p)，计算第二位移值D2pre(i,p)。光流法计算公式为：
D2pre(i,p)＝[A(i,p)TA(i,p)]-1A(i,p)Tb(i×p)；
其中：A(i,p)＝{GradientX[Imgpre(i,p)],GradientY[Imgpre(i,p)]}；
b(i×p)＝{Imgpre(i,p)-Imgpos(i,p)}；
GradientX表示i方向的梯度函数；GradientY表示p方向的梯度函数。
在一实施例中，在图1所示实施例的基础上，上述步骤S104，包括：
对所述第一形变位移值和所述第二形变位移值进行矩阵求和运算，得到所述最终形变位移值；
根据所述最终形变位移值，生成初始应变弹性图像；
对所述初始应变弹性图像进行纵向差分运算，并对纵向差分运算后的初始应变弹性图像进行归一化，映射得到所述目标应变弹性图像。
在本实施例中，集合运算为矩阵求和运算，最终形变位移值为：
Dpre(i,p)＝D1pre(i,p)+D2pre(i,p)；
根据最终位移值Dpre(i,p)可计算初始应变弹性图像，对初始应变弹性图像作纵向差分运算，对差分后的图像Stresspre(i,p)进行图像归一化并映射成目标应变弹性图像Stress(i,p)。可选地，采用SG－I滤波进行纵向差分运算，计算公式为：
Figure BDA0003534475720000111
Figure BDA0003534475720000112
其中(2M+1)为滤波器的长度，max[]为求最大值函数。
为了执行上述方法实施例对应的超声弹性成像方法，以实现相应的功能和技术效果。参见图6，图6示出了本申请实施例提供的一种超声弹性成像装置的结构框图。为了便于说明，仅示出了与本实施例相关的部分，本申请实施例提供的超声弹性成像装置，包括：
获取模块601，用于获取目标组织区域形变前与形变后的相邻稀疏信号，所述相邻稀疏信号包括所述第一稀疏回波信号和与所述第一稀疏回波信号相邻的第二稀疏回波信号，所述第一稀疏回波信号为所述目标组织区域形变前的稀疏信号，所述第二稀疏回波信号为所述目标组织区域形变后的稀疏信号；
第一计算模块602，用于对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值；
第二计算模块603，用于对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值；
成像模块604，用于对所述第一形变位移值和所述第二形变位移值进行集合运算，得到所述目标组织区域的最终形变位移值，并根据所述最终形变位移值生成所述目标组织区域的应变弹性图像。
作为优选，所述获取模块601，包括：
获取单元，用于获取所述目标组织区域形变前与形变后的多帧原始回波信号；
提取单元，用于基于预设时间间隔，提取多帧所述原始回波信号中相邻的两帧原始回波信号；
下采样单元，用于对相邻的两帧所述原始回波信号进行下采样，得到所述第一稀疏回波信号和所述第二稀疏回波信号。
作为优选，所述第一计算模块602，包括：
第一估算单元，用于基于预设的零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在同一位置纵向点上进行相位差估算，得到所述第一稀疏回波信号与所述第二稀疏回波信号之间的第一相位差；
第二估算单元，用于基于所述零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在与所述纵向点相邻的多个横向点进行相位差估算，得到多个横向点对应的第二相位差；
第一运算单元，用于对多个所述第二相位差进行均值运算，得到所述第一形变位移值。
作为优选，所述零相位迭代算法的计算公式为：
Figure BDA0003534475720000121
其中，Dpre(i+u,p)为上一帧信号相对于当前帧信号的形变位移值，Z1pre(i+u+k,p+h)为所述第一稀疏回波信号，Z1pos(i+u+k,p+h)为所述第二稀疏回波信号，i为横向方向的信号线号，p为纵向方向的信号点号，u为横向内的计算次数，j为虚数因子，K为Mask窗口的横向尺寸，H为Mask窗口的纵向尺寸，w0为角频率，conj为共轭复数。
作为优选，所述第二计算模块603，包括：
第二运算单元，用于对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，得到所述第一稀疏回波信号的第一高位宽模值和所述第二稀疏回波信号的第二高位宽模值；
压缩单元，用于基于预设的数据动态范围，对所述第一高位宽模值和所述第二高位宽模值进行压缩到低位宽，得到所述第一稀疏回波信号对应的第一线性图像数据和所述第二稀疏回波信号对应的第二线性图像数据；
第三计算单元，用于根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值。
作为优选，所述第三计算单元具体用于：
基于预设光流法，对所述第一线性图像数据和所述第二线性图像数据进行运算，得到所述第二形变位移值。
作为优选，所述成像模块604，包括：
第三运算单元，用于对所述第一形变位移值和所述第二形变位移值进行矩阵求和运算，得到所述最终形变位移值；
生成单元，用于根据所述最终形变位移值，生成初始应变弹性图像；
映射单元，用于对所述初始应变弹性图像进行纵向差分运算，并对纵向差分运算后的初始应变弹性图像进行归一化，映射得到所述目标应变弹性图像。
上述的超声弹性成像装置可实施上述方法实施例的超声弹性成像方法。上述方法实施例中的可选项也适用于本实施例，这里不再详述。本申请实施例的其余内容可参照上述方法实施例的内容，在本实施例中，不再进行赘述。
图7为本申请一实施例提供的计算机设备的结构示意图。如图7所示，该实施例的计算机设备7包括：至少一个处理器70(图7中仅示出一个)处理器、存储器71以及存储在所述存储器71中并可在所述至少一个处理器70上运行的计算机程序72，所述处理器70执行所述计算机程序72时实现上述任意方法实施例中的步骤。
所述计算机设备7可以是智能手机、平板电脑、桌上型计算机和云端服务器等计算设备。该计算机设备可包括但不仅限于处理器70、存储器71。本领域技术人员可以理解，图7仅仅是计算机设备7的举例，并不构成对计算机设备7的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如还可以包括输入输出设备、网络接入设备等。
所称处理器70可以是中央处理单元(Central Processing Unit，CPU)，该处理器70还可以是其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现成可编程门阵列(Field－Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
所述存储器71在一些实施例中可以是所述计算机设备7的内部存储单元，例如计算机设备7的硬盘或内存。所述存储器71在另一些实施例中也可以是所述计算机设备7的外部存储设备，例如所述计算机设备7上配备的插接式硬盘，智能存储卡(Smart Media Card，SMC)，安全数字(Secure Digital，SD)卡，闪存卡(Flash Card)等。进一步地，所述存储器71还可以既包括所述计算机设备7的内部存储单元也包括外部存储设备。所述存储器71用于存储操作系统、应用程序、引导装载程序(BootLoader)、数据以及其他程序等，例如所述计算机程序的程序代码等。所述存储器71还可以用于暂时地存储已经输出或者将要输出的数据。
另外，本申请实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述任意方法实施例中的步骤。
本申请实施例提供了一种计算机程序产品，当计算机程序产品在计算机设备上运行时，使得计算机设备执行时实现上述各个方法实施例中的步骤。
在本申请所提供的几个实施例中，可以理解的是，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意的是，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read－Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述的具体实施例，对本申请的目的、技术方案和有益效果进行了进一步的详细说明，应当理解，以上所述仅为本申请的具体实施例而已，并不用于限定本申请的保护范围。特别指出，对于本领域技术人员来说，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。

Claims (10)

1.一种超声弹性成像方法，其特征在于，包括：
获取目标组织区域形变前与形变后的相邻稀疏信号，所述相邻稀疏信号包括第一稀疏回波信号和与所述第一稀疏回波信号相邻的第二稀疏回波信号，所述第一稀疏回波信号为所述目标组织区域形变前的稀疏信号，所述第二稀疏回波信号为所述目标组织区域形变后的稀疏信号；
对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值；
对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值；
对所述第一形变位移值和所述第二形变位移值进行集合运算，得到所述目标组织区域的最终形变位移值，并根据所述最终形变位移值生成所述目标组织区域的目标应变弹性图像。
2.如权利要求1所述的超声弹性成像方法，其特征在于，所述获取目标组织区域形变前与形变后的相邻稀疏信号，包括：
获取所述目标组织区域形变前与形变后的多帧原始回波信号；
基于预设时间间隔，提取多帧所述原始回波信号中相邻的两帧原始回波信号；
对相邻的两帧所述原始回波信号进行下采样，得到所述第一稀疏回波信号和所述第二稀疏回波信号。
3.如权利要求1所述的超声弹性成像方法，其特征在于，所述对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值，包括：
基于预设的零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在同一位置纵向点上进行相位差估算，得到所述第一稀疏回波信号与所述第二稀疏回波信号之间的第一相位差；
基于所述零相位迭代算法，对所述第一稀疏回波信号和所述第二稀疏回波信号在与所述纵向点相邻的多个横向点进行相位差估算，得到多个横向点对应的第二相位差；
对多个所述第二相位差进行均值运算，得到所述第一形变位移值。
4.如权利要求3所述的超声弹性成像方法，其特征在于，所述零相位迭代算法的计算公式为：
Figure FDA0003534475710000021
其中，Dpre(i+u,p)为上一帧信号相对于当前帧信号的形变位移值，Z1pre(i+u+k,p+h)为所述第一稀疏回波信号，Z1pos(i+u+k,p+h)为所述第二稀疏回波信号，i为横向方向的信号线号，p为纵向方向的信号点号，u为横向内的计算次数，j为虚数因子，K为Mask窗口的横向尺寸，H为Mask窗口的纵向尺寸，w0为角频率，conj为共轭复数。
5.如权利要求1所述的超声弹性成像方法，其特征在于，所述对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值，包括：
对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，得到所述第一稀疏回波信号的第一高位宽模值和所述第二稀疏回波信号的第二高位宽模值；
基于预设的数据动态范围，对所述第一高位宽模值和所述第二高位宽模值进行压缩到低位宽，得到所述第一稀疏回波信号对应的第一线性图像数据和所述第二稀疏回波信号对应的第二线性图像数据；
根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值。
6.如权利要求5所述的超声弹性成像方法，其特征在于，所述根据所述第一线性图像数据和所述第二线性图像数据，计算所述第二形变位移值，包括：
基于预设光流法，对所述第一线性图像数据和所述第二线性图像数据进行运算，得到所述第二形变位移值。
7.如权利要求1所述的超声弹性成像方法，其特征在于，所述对所述第一形变位移值和所述第二形变位移值进行集合运算，得到所述目标组织区域的最终形变位移值，并根据所述最终形变位移值生成所述目标组织区域的目标应变弹性图像，包括：
对所述第一形变位移值和所述第二形变位移值进行矩阵求和运算，得到所述最终形变位移值；
根据所述最终形变位移值，生成初始应变弹性图像；
对所述初始应变弹性图像进行纵向差分运算，并对纵向差分运算后的初始应变弹性图像进行归一化，映射得到所述目标应变弹性图像。
8.一种超声弹性成像装置，其特征在于，包括：
获取模块，用于获取目标组织区域形变前与形变后的相邻稀疏信号，所述相邻稀疏信号包括所述第一稀疏回波信号和与所述第一稀疏回波信号相邻的第二稀疏回波信号，所述第一稀疏回波信号为所述目标组织区域形变前的稀疏信号，所述第二稀疏回波信号为所述目标组织区域形变后的稀疏信号；
第一计算模块，用于对所述第一稀疏回波信号和所述第二稀疏回波信号进行纵向点的相位差估算，并根据估算得到的目标相位差计算第一形变位移值；
第二计算模块，用于对所述第一稀疏回波信号和所述第二稀疏回波信号进行求模运算，确定所述第一稀疏回波信号和所述第二稀疏回波信号的线性图像数据，并根据所述线性图像数据计算第二形变位移值；
成像模块，用于对所述第一形变位移值和所述第二形变位移值进行集合运算，得到所述目标组织区域的最终形变位移值，并根据所述最终形变位移值生成所述目标组织区域的应变弹性图像。
9.一种计算机设备，其特征在于，包括处理器和存储器，所述存储器用于存储计算机程序，所述计算机程序被所述处理器执行时实现如权利要求1至7任一项所述的超声弹性成像方法。
10.一种计算机可读存储介质，其特征在于，其存储有计算机程序，所述计算机程序被处理器执行时实现如权利要求1至7任一项所述的超声弹性成像方法。
CN202210234626.3A 2022-03-07 2022-03-07 超声弹性成像方法、装置、设备及存储介质 Active CN114668419B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210234626.3A CN114668419B (zh) 2022-03-07 2022-03-07 超声弹性成像方法、装置、设备及存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210234626.3A CN114668419B (zh) 2022-03-07 2022-03-07 超声弹性成像方法、装置、设备及存储介质

Publications (2)

Publication Number Publication Date
CN114668419A true CN114668419A (zh) 2022-06-28
CN114668419B CN114668419B (zh) 2023-05-23

Family

ID=82072798

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210234626.3A Active CN114668419B (zh) 2022-03-07 2022-03-07 超声弹性成像方法、装置、设备及存储介质

Country Status (1)

Country Link
CN (1) CN114668419B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116712101A (zh) * 2023-08-04 2023-09-08 深圳市联影高端医疗装备创新研究院 超声图像生成方法、装置、计算机设备和存储介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040147836A1 (en) * 2003-01-27 2004-07-29 Lagendijk Jan Jacob Wouter Device for image guided automated insertion of an elongated hollow needle to a desired location in an animal body for effecting radiation therapy of cancerous tissue
CN101014086A (zh) * 2007-01-31 2007-08-08 天津大学 采用运动检测和自适应加权滤波的去隔行装置
CN104605888A (zh) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 超声信号的时延估计方法及系统
CN107198545A (zh) * 2017-06-06 2017-09-26 苏州国科昂卓医疗科技有限公司 生物组织的弹性位移及应变估计方法、装置
US20200405269A1 (en) * 2018-02-27 2020-12-31 Koninklijke Philips N.V. Ultrasound system with a neural network for producing images from undersampled ultrasound data
CN112904349A (zh) * 2021-01-28 2021-06-04 清华大学 一种超声弹性成像方法
CN112998750A (zh) * 2021-02-22 2021-06-22 深圳华声医疗技术股份有限公司 超声图像合成方法、装置、超声设备及存储介质

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040147836A1 (en) * 2003-01-27 2004-07-29 Lagendijk Jan Jacob Wouter Device for image guided automated insertion of an elongated hollow needle to a desired location in an animal body for effecting radiation therapy of cancerous tissue
CN101014086A (zh) * 2007-01-31 2007-08-08 天津大学 采用运动检测和自适应加权滤波的去隔行装置
CN104605888A (zh) * 2015-01-22 2015-05-13 飞依诺科技(苏州)有限公司 超声信号的时延估计方法及系统
CN107198545A (zh) * 2017-06-06 2017-09-26 苏州国科昂卓医疗科技有限公司 生物组织的弹性位移及应变估计方法、装置
US20200405269A1 (en) * 2018-02-27 2020-12-31 Koninklijke Philips N.V. Ultrasound system with a neural network for producing images from undersampled ultrasound data
CN112904349A (zh) * 2021-01-28 2021-06-04 清华大学 一种超声弹性成像方法
CN112998750A (zh) * 2021-02-22 2021-06-22 深圳华声医疗技术股份有限公司 超声图像合成方法、装置、超声设备及存储介质

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116712101A (zh) * 2023-08-04 2023-09-08 深圳市联影高端医疗装备创新研究院 超声图像生成方法、装置、计算机设备和存储介质
CN116712101B (zh) * 2023-08-04 2023-11-21 深圳市联影高端医疗装备创新研究院 超声图像生成方法、装置、计算机设备和存储介质

Also Published As

Publication number Publication date
CN114668419B (zh) 2023-05-23

Similar Documents

Publication Publication Date Title
Tay et al. Ultrasound despeckling for contrast enhancement
US6322505B1 (en) Medical diagnostic ultrasound system and method for post processing
US8111947B2 (en) Image processing apparatus and method which match two images based on a shift vector
JP4217741B2 (ja) ２つの画像の位置合わせパラメータ判定方法、装置、プログラム
US20180220997A1 (en) System and method for accelerated clutter filtering in ultrasound blood flow imaging using randomized ultrasound data
Vedula et al. Towards CT-quality ultrasound imaging using deep learning
US8403850B2 (en) Rapid two/three-dimensional sector strain imaging
US20230386048A1 (en) Ultrasound image segmentation method and apparatus, terminal device, and storage medium
CN114668419B (zh) 超声弹性成像方法、装置、设备及存储介质
Kuzmin et al. Multi-frame elastography using a handheld force-controlled ultrasound probe
JP6207972B2 (ja) 超音波診断装置、画像処理装置及び画像処理プログラム
JP2015532848A (ja) 画像組織（ｔｉｓｓｕｅ）動作推定
CN111507912B (zh) 乳腺图像增强方法、装置及计算机可读存储介质
CN110766661B (zh) 磁共振成像的水脂分离方法、磁共振成像方法和设备
CN108961161B (zh) 一种图像数据处理方法、装置和计算机存储介质
CN116797683A (zh) 一种数字减影血管造影控制系统及控制方法
WO2012069833A1 (en) Process and apparatus for data registration
Khodadadi et al. Edge-preserving ultrasonic strain imaging with uniform precision
Biswas et al. An FPGA-based architecture of DSC–SRI units specially for motion blind ultrasound systems
CN110930394B (zh) 测量肌肉肌纤维束线斜率和羽状角的方法及终端设备
Ko et al. Accurate depth estimation of skin surface using a light‐field camera toward dynamic haptic palpation
CN111369425A (zh) 图像处理方法、装置、电子设备和计算机可读介质
JP3176652B2 (ja) 画像歪補正装置、及び画像間演算装置、画像動き算出装置
CN116596797A (zh) 图像处理方法和装置
CN117575922A (zh) 一种基于数字化治疗的牙齿与牙冠位姿确定方法及系统

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20221018

Address after: 100000303-24, Floor 3, Building 4, Yard 9, Yike Road, Life Science Park, Changping District, Beijing (cluster registration)

Applicant after: Yichao Medical Technology (Beijing) Co.,Ltd.

Address before: 430206 No. 412, east of floor 4, building B1, phase II, Wuhan hi tech medical device Park, No. 818, Gaoxin Avenue, Donghu New Technology Development Zone, Wuhan, Hubei Province (Wuhan area of free trade zone)

Applicant before: Yichao medical technology (Wuhan) Co.,Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20230425

Address after: 430000 West facing 701-715, Building B1, 7th Floor, Phase II, Block B, Wuhan High tech Medical Equipment Park, No. 818 Gaoxin Avenue, Donghu New Technology Development Zone, Wuhan, Hubei Province (Wuhan area of the Free Trade Zone)

Applicant after: Yichao Technology (Wuhan) Co.,Ltd.

Address before: 100000303-24, Floor 3, Building 4, Yard 9, Yike Road, Life Science Park, Changping District, Beijing (cluster registration)

Applicant before: Yichao Medical Technology (Beijing) Co.,Ltd.

GR01 Patent grant
GR01 Patent grant