CN110811688B - 多角度平面波重复复合的超快超声多普勒血流估计方法 - Google Patents
多角度平面波重复复合的超快超声多普勒血流估计方法 Download PDFInfo
- Publication number
- CN110811688B CN110811688B CN201911213180.0A CN201911213180A CN110811688B CN 110811688 B CN110811688 B CN 110811688B CN 201911213180 A CN201911213180 A CN 201911213180A CN 110811688 B CN110811688 B CN 110811688B
- Authority
- CN
- China
- Prior art keywords
- blood flow
- beta
- signal
- frame
- signals
- 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
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/06—Measuring blood flow
-
- 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/488—Diagnostic techniques involving Doppler signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Hematology (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种多角度平面波重复复合的超快超声多普勒血流估计方法。所述方法包括:控制超声换能器循环发射偏转角度依次为{0,‑β,β,‑β}的平面波;接收经待扫描区域反射的平面波的回波信号;搜索回波信号中偏转角度依次为{0,‑β,β}和偏转角度依次为{β,‑β,0}的三帧连续信号,得到多个帧组;对每个帧组中的三帧连续信号进行叠加,得到复合射频信号;根据复合射频信号计算血流速度;对血流速度进行彩色编码生成二维血流速度图像。在相干复合中,最后一个射频回波信号被重复地用于下一次复合,与所有射频回波信号只复合一次的传统相干复合技术相比,脉冲重复频率提高了50%,从而提高了超快超声多普勒技术最大血流速度可检测范围。
Description
技术领域
本发明涉及超声成像技术领域,特别是涉及多角度平面波重复复合的超快超声多普勒血流估计方法。
背景技术
血管***中的血流速度分布是指在血管腔的径向切面上血流速度的空间分布。血液流动过程中,由于与血管壁之间的摩擦力和其自身的黏性,在径向上分层流动,每层流速不同,管腔中心层的速度最大,趋向血管壁的位置,血流速度逐渐降低。在健康的血管中,血流速度分布呈抛物线型,而在发生动脉粥样硬化的病变血管中,血流速度分布受斑块影响发生变化,不再呈抛物线型,甚至出现湍流和涡流。血流动力学以血流速度分布为基础,可根据血流速度分布计算剪切应力、速度剪切率和壁面剪切率等血流动力学参数。
目前,能够获取血流速度分布的超声技术主要有两类:图像测速技术和超声多普勒技术。其中,图像测速技术是在连续两帧B超图像中跟踪血流散斑以获取血流速度分布,但是图像测速技术的处理时间滞后于临床信号的采集,血流速度分布不能实时显示。超声多普勒技术以反射法测量多普勒频移为原理估计血流速度,具有处理速度快,测量精度高等优势,但是,根据奈奎斯特采样定理,当脉冲重复频率小于血流多普勒频移信号最大频率的两倍时,会出现速度混叠,造成测量失败。
超声多普勒技术包括脉冲多普勒双工扫描成像和彩色多普勒血流成像两类。脉冲多普勒双工扫描成像将B超成像和脉冲多普勒技术结合,在显示周围组织结构的同时显示血流多普勒频移的频谱,具有测距和测速的双重功能。由于牺牲了一部分脉冲重复频率用于B超成像,该技术能够检测的最大血流速度较小。而彩色多普勒血流成像在B超图像上同步显示二维彩色血流图像,用红蓝两种颜色区别速度方向,并使用两种颜色的亮度表示速度大小。与脉冲多普勒相比,彩色多普勒血流成像在空间上更直观地显示了血流速度分布。然而,为了保证彩色血流成像的帧频,该技术在每个扫描位置上的脉冲发射数量有限,对血流的观测时间较短,信噪比较低。综上所述,以上两种常用的超声多普勒技术皆受限于脉冲重复频率。
为了解决上述问题,新一代超快超声随之出现,主要发展方向是平面波成像技术,采用全孔径发射超声信号,通过一次发射即可得到整个成像区域的回波信号,其脉冲重复频率等于帧频,在5cm的成像深度范围内,高达每秒15000帧。但是,由于缺少发射焦点,回波信号信噪比低,成像质量差。
相干复合平面波算法的提出改进了平面波的成像质量。该算法通过依次延迟全孔径内所有阵元的激发时间,从多个对称偏转的发射角度获得同一成像区域的多帧平面波图像,并将多帧图像相干叠加得到复合图像。由于相干复合平面波算法平均了多帧图像,因此可以有效平滑噪声,提高信噪比。并且,由于发射角度不同、可以有效解决成像目标的边缘缺失等问题。
发射角度的设置包括极性偏转角序列和线性偏转角序列。在极性偏转角序列中,角度对称的两次扫描时间间隔更短,更利于抵消噪声保留有用信息。因此,基于极性偏转角序列相干复合获取的多帧图像,其信噪比优于基于线性偏转角序列的复合结果。然而,使用多帧图像生成一帧复合图像,相干复合平面波算法在提高信噪比的同时导致了脉冲重复频率降低,从而降低了超快超声多普勒技术最大血流速度的可检测范围。
发明内容
本发明的目的是提供一种多角度平面波重复复合的超快超声多普勒血流估计方法,将脉冲重复频率提高50%,从而提高超快超声多普勒技术最大血流速度可检测范围。
为实现上述目的,本发明提供了如下方案:
多角度平面波重复复合的超快超声多普勒血流估计方法,所述方法包括:
控制超声换能器向待扫描区域循环发射偏转角度依次为{0,-β,β,-β}的平面波;
接收经所述待扫描区域反射的平面波的回波信号;所述回波信号为偏转角度按照基本循环单元循环的信号;所述基本循环单元的角度排列依次为0、-β、β、-β;
搜索所述回波信号中偏转角度依次为{0,-β,β}和偏转角度依次为{β,-β,0}的三帧连续信号为一个帧组,得到多个帧组,多个所述帧组满足相邻的两个帧组中前一帧组的末尾帧信号为后一帧组的起始帧信号;
对每个帧组中的三帧连续信号进行叠加,得到复合射频信号;
根据所述复合射频信号计算血流速度;
对所述血流速度进行彩色编码生成二维血流速度图像。
可选的,控制超声换能器向待扫描区域发射偏转角为0的平面波时,设置超声换能器全孔径中的所有阵元同步激发;
控制超声换能器向待扫描区域发射偏转角为-β或β的平面波时,设置超声换能器全孔径中的所有阵元基于偏转角度计算的延迟时间依次激发。
可选的,在所述搜索所述回波信号中偏转角度依次为{0,-β,β}和偏转角度依次为{β,-β,0}的三帧连续信号为一个帧组,得到多个帧组之后,在所述对每个帧组中的三帧连续信号进行叠加,得到复合射频信号之前,还包括:
对每个帧组中的三帧连续信号进行局域配准。
可选的,所述对每个帧组中的三帧连续信号进行局域配准具体包括:
针对每个帧组的三帧连续信号都执行以下操作:
计算起始帧信号相对于中间帧信号的第一偏移量;
计算末尾帧信号相对于中间帧信号的第二偏移量;
根据所述第一偏移量和所述第二偏移量对所述每个帧组的三帧连续信号进行局域配准。
可选的,所述根据所述复合射频信号计算血流速度具体包括:
在连续的多帧所述复合射频信号中分别提取各帧复合射频信号的血流多普勒信号,得到血管腔内不同径向位置的血流多普勒信号;
对每个所述血流多普勒信号均进行正交解调,得到每个所述血流多普勒信号的正交分量和同相分量;
根据每个所述血流多普勒信号的正交分量和同相分量,分别计算得到每个血流多普勒信号的相移;
根据每个血流多普勒信号的相移分别计算得到每个血流多普勒信号的血流速度。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
控制超声换能器循环发射偏转角度依次为{0,-β,β,-β}的四角度平面波序列,在相干复合过程中使用偏转角为{0,-β,β}或{β,-β,0}的连续的三帧射频回波信号。即用于相干复合的三角度射频回波信号中,最后一个射频回波信号被重复地用于下一次复合,与所有射频回波信号只复合一次的传统相干复合技术相比,脉冲重复频率提高了50%,从而提高了超快超声多普勒技术最大血流速度可检测范围。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍。
图1为本发明提供的多角度平面波重复复合的超快超声多普勒血流估计方法的流程图;
图2为本发明实施例提供的扫描示意图;
图3为本发明实施例提供的四偏转角循环发射平面波时间序列示意图;
图4为本发明实施例提供的基于三偏转角平面波的复合射频信号时间序列示意图;
图5为本发明实施例提供的三偏转角平面波的修正及复合示意图;
图6为本发明实施例测量的径向血流速度分布图;
图7为本发明实施例提供的血流速度成像图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。
本发明的目的是提供一种多角度平面波重复复合的超快超声多普勒血流估计方法,将脉冲重复频率提高50%,从而提高超快超声多普勒技术最大血流速度可检测范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明提供的多角度平面波重复复合的超快超声多普勒血流估计方法流程图,如图1所示,该方法包括:
S101:控制超声换能器向待扫描区域循环发射偏转角度依次为{0,-β,β,-β}的平面波。
具体控制方式为:控制超声换能器向待扫描区域发射偏转角为0的平面波时,设置超声换能器全孔径中的所有阵元同步激发;控制超声换能器向待扫描区域发射偏转角为-β或β的平面波时,设置超声换能器全孔径中的所有阵元基于偏转角度计算的延迟时间依次激发。
本实施例中,设置循环发射平面波扫描血流的扫描范围为X×Z。
S102:接收经待扫描区域反射的平面波的回波信号;回波信号为偏转角度按照基本循环单元循环的信号;基本循环单元的角度排列依次为0、-β、β、-β。
以扫描区域内的点(x,z)为例,说明超声换能器接收到的平面波的回波信号具体控制过程为:
第n个阵元发射的平面波信号到达扫描区域内的点(x,z)的传输时间为该信号经点(x,z)反射后,由第n'个阵元接收的传输时间为对于点(x,z),信号总传输时间为t(x,z)=tgo(x,z)+tback(x,z)。
以全孔径接收经扫描区域内所有散射点反射的信号,该信号称为通道数据(Channel Data,CD)。
再基于超声延迟叠加法,对通道数据执行波束合成生成射频回波信号RF。以生成扫描区域内点(x,z)处的射频回波信号RF(x,z)为例,根据公式:
生成,其中,2a是波束合成中使用的合成孔径。
对每帧通道数据执行波束合成,获取多角度平面波的射频回波信号,其偏转角度依次为{…,0,-β,β,-β,0,-β,β,-β…}。本实施例中,最终得到的射频回波信号即为平面波回波信号。
S103:搜索回波信号中偏转角度依次为{0,-β,β}和偏转角度依次为{β,-β,0}的三帧连续信号为一个帧组,得到多个帧组,多个帧组满足相邻的两个帧组中前一帧组的末尾帧信号为后一帧组的起始帧信号。
S104:对每个帧组中的三帧连续信号进行局域配准。
针对每个帧组的三帧连续信号都执行以下操作:
计算起始帧信号相对于中间帧信号的第一偏移量。
计算末尾帧信号相对于中间帧信号的第二偏移量。
根据第一偏移量和第二偏移量对每个帧组的三帧连续信号进行局域配准。
本实施例中,局域配准的具体控制过程为:将三帧连续的偏转角度为{0,-β,β}或{β,-β,0}的射频回波信号{RF0,RF-β,RFβ}或{RFβ,RF-β,RF0},分别分割为H×K(横向分为H块,纵向分为K块)个局域射频回波信号。分别计算局域射频回波信号RF0(h,k)和RFβ(h,k)相对于RF-β(h,k)的偏移量,根据偏移量对RF0(h,k)、RFβ(h,k)和RF-β(h,k)进行配准,其中,1≤h≤H,1≤k≤K。
S105:对每个帧组中的三帧连续信号进行叠加,得到复合射频信号。
S106:根据复合射频信号计算血流速度。
在连续的多帧复合射频信号中分别提取各帧复合射频信号的血流多普勒信号,得到血管腔内不同径向位置的血流多普勒信号。根据公式R(k)=RFk(n,z)提取第k帧复合射频信号中不同径向位置的血流多普勒信号,RFk(n,z)表示第k帧复合射频信号中第n线第z个采样点,1≤k≤K。
对每个血流多普勒信号均进行正交解调,得到每个血流多普勒信号的正交分量和同相分量。根据公式R(k)=r(k)+j*I(k)分别对每个血流多普勒信号进行正交解调;其中r(k)为正交解调后的R(k)信号的同相分量,I(k)为正交解调后的R(k)信号的正交分量,j为虚数单位。
根据每个血流多普勒信号的正交分量和同相分量,分别计算得到每个血流多普勒信号的相移。采用自相关算法根据公式计算每个血流多普勒信号的相移;其中K表示连续的多帧复合射频信号为K帧,k为第k帧血流多普勒信号,k-1为第k-1帧血流多普勒信号,R(k)表示第k帧血流多普勒信号,I(k)为正交解调后的R(k)信号的正交分量,R(k-1)表示第k-1帧血流多普勒信号,I(k-1)为正交解调后的R(k-1)信号的正交分量。
根据每个血流多普勒信号的相移分别计算得到每个血流多普勒信号的血流速度。根据超声多普勒公式计算每个血流多普勒信号的血流速度;c表示声速,fp表示脉冲重复频率,f0表示超声换能器的中心频率,表示对应血流多普勒信号的相移。
S107:对血流速度进行彩色编码生成二维血流速度图像。
下面给出一种具体的实施例对本发明的技术方案进行举例说明。如图2所示,建立血流模型,设置超声换能器中心最大速度为1.5m/s,偏转倾斜角度为45°。设置超声换能器循环发射平面波,包括四个偏转角度为:{0°,-15°,15°,-15°}。如图3所示,根据设置的四偏转发射循环,循环发射平面波进行扫描。基于超声延迟叠加法,对超声换能器全孔径接收的每帧通道数据进行波束合成,以获取多角度平面波的射频回波时间序列,其偏转角度为{…,0°,-15°,15°,-15°,0°,-15°,15°,-15°,…}。
对三帧连续的偏转角度为{0°,-15°,15°}或{15°,-15°,0°}的射频回波信号进行局域分割获取局域射频回波信号。计算0°和15°局域射频回波信号相对于-15°局域射频回波信号的偏移,并根据偏移量对三个角度的局域射频回波信号进行配准。将配准后的三帧射频回波信号叠加,生成一帧复合射频信号,其中,每次复合的最后一帧射频回波信号被重复用于下一次复合(如图4、图5所示)。
在连续的K=10帧复合射频信号中,提取不同径向位置的血流多普勒信号R(k),并进行正交解调,公式为R(k)=r(k)+j*I(k),其中,1≤k≤10,r(k)为正交解调后的R(k)信号的同相分量,I(k)为正交解调后的R(k)信号的正交分量,j为虚数单位。
基于自相关算法计算血流多普勒信号的相移公式为其中K表示连续的多帧复合射频信号为K帧,k为第k帧血流多普勒信号,k-1为第k-1帧血流多普勒信号,R(k)表示第k帧血流多普勒信号,I(k)为正交解调后的R(k)信号的正交分量,R(k-1)表示第k-1帧血流多普勒信号,I(k-1)为正交解调后的R(k-1)信号的正交分量。
如图6所示,根据超声多普勒公式提取50个不同的径向位置的血流速度信息v,公式为:其中c表示声速,在人体组织中通常为1540m/s,fp=10000Hz表示脉冲重复频率(在超快超声中等于复合射频信号的帧频),f0=10MHz表示超声换能器的中心频率,表示对应血流多普勒信号的相移。
最后,基于成像范围内的所有血流速度信息,进行彩色编码,获取视野内的二维血流速度图像,结果如图7所示。
本发明中循环发射偏转角为{0,-β,β,-β}的四角度平面波序列,而在相干复合过程中仅使用偏转角为{0,-β,β}或{β,-β,0}的连续的三帧射频回波信号。每次复合使用的三个偏转角平面波都为极性偏转角序列。
本发明引入局域射频回波信号匹配技术,与整帧回波信号直接复合的传统复合方法不同,本发明在三角度回波信号复合之前先进行局域匹配,相对于整帧配准的全局配准方法,本发明对不同径向位置的血流射频信号进行局域配准的方法使血流的运动伪影最小化。
另外,假设未复合前的脉冲重复频率为PRFsingle,对于传统的三角度平面波相干复合方法,其脉冲重复频率为在本发明中,每次复合使用的三角度平面波,最后一帧被重复用于下一次复合,所以,三角度平面波重复复合后的脉冲重复频率为因此本发明的重复复合的方法使脉冲重复频率提高了50%。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (5)
1.多角度平面波重复复合的超快超声多普勒血流估计方法,其特征在于,所述方法包括:
控制超声换能器向待扫描区域循环发射偏转角度依次为{0,-β,β,-β}的平面波;
接收经所述待扫描区域反射的平面波的回波信号;所述回波信号为偏转角度按照基本循环单元循环的信号;所述基本循环单元的角度排列依次为0、-β、β、-β;
搜索所述回波信号中偏转角度依次为{0,-β,β}和偏转角度依次为{β,-β,0}的三帧连续信号为一个帧组,得到多个帧组,多个所述帧组满足相邻的两个帧组中前一帧组的末尾帧信号为后一帧组的起始帧信号;
对每个帧组中的三帧连续信号进行叠加,得到复合射频信号;
根据所述复合射频信号计算血流速度;
对所述血流速度进行彩色编码生成二维血流速度图像。
2.根据权利要求1所述的多角度平面波重复复合的超快超声多普勒血流估计方法,其特征在于,
控制超声换能器向待扫描区域发射偏转角为0的平面波时,设置超声换能器全孔径中的所有阵元同步激发;
控制超声换能器向待扫描区域发射偏转角为-β或β的平面波时,设置超声换能器全孔径中的所有阵元基于偏转角度计算的延迟时间依次激发。
3.根据权利要求1所述的多角度平面波重复复合的超快超声多普勒血流估计方法,其特征在于,
在搜索所述回波信号中偏转角度依次为{0,-β,β}和偏转角度依次为{β,-β,0}的三帧连续信号为一个帧组,得到多个帧组之后,在所述对每个帧组中的三帧连续信号进行叠加,得到复合射频信号之前,还包括:
对每个帧组中的三帧连续信号进行局域配准。
4.根据权利要求3所述的多角度平面波重复复合的超快超声多普勒血流估计方法,其特征在于,
所述对每个帧组中的三帧连续信号进行局域配准具体包括:
针对每个帧组的三帧连续信号都执行以下操作:
计算起始帧信号相对于中间帧信号的第一偏移量;
计算末尾帧信号相对于中间帧信号的第二偏移量;
根据所述第一偏移量和所述第二偏移量对所述每个帧组的三帧连续信号进行局域配准。
5.根据权利要求1所述的多角度平面波重复复合的超快超声多普勒血流估计方法,其特征在于,
所述根据所述复合射频信号计算血流速度具体包括:
在连续的多帧所述复合射频信号中分别提取各帧复合射频信号的血流多普勒信号,得到血管腔内不同径向位置的血流多普勒信号;
对每个所述血流多普勒信号均进行正交解调,得到每个所述血流多普勒信号的正交分量和同相分量;
根据每个所述血流多普勒信号的正交分量和同相分量,分别计算得到每个血流多普勒信号的相移;
根据每个血流多普勒信号的相移分别计算得到每个血流多普勒信号的血流速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911213180.0A CN110811688B (zh) | 2019-12-02 | 2019-12-02 | 多角度平面波重复复合的超快超声多普勒血流估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911213180.0A CN110811688B (zh) | 2019-12-02 | 2019-12-02 | 多角度平面波重复复合的超快超声多普勒血流估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110811688A CN110811688A (zh) | 2020-02-21 |
CN110811688B true CN110811688B (zh) | 2021-10-01 |
Family
ID=69542064
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911213180.0A Active CN110811688B (zh) | 2019-12-02 | 2019-12-02 | 多角度平面波重复复合的超快超声多普勒血流估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110811688B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021108953A1 (zh) * | 2019-12-02 | 2021-06-10 | 云南大学 | 多角度平面波重复复合的超快超声多普勒血流估计方法 |
CN111227867B (zh) * | 2020-03-10 | 2021-01-12 | 云南大学 | 一种超声多普勒血流成像方法及*** |
CN113768533B (zh) * | 2020-06-10 | 2024-05-14 | 无锡祥生医疗科技股份有限公司 | 超声显影装置和超声显影方法 |
CN112843497B (zh) * | 2021-01-05 | 2022-09-16 | 中国科学院上海高等研究院 | 一种基于射频偏转腔技术的质子束流扫描装置及扫描方法 |
CN113397591B (zh) * | 2021-06-16 | 2022-11-29 | 复旦大学 | 微纳马达实时成像与追踪方法、装置及微纳马达控制方法 |
CN115990034A (zh) * | 2021-10-19 | 2023-04-21 | 复旦大学 | 一种随机空间采样的超快超声血流成像方法及*** |
CN114391876A (zh) * | 2022-02-14 | 2022-04-26 | 逸超医疗科技(武汉)有限公司 | 多普勒频谱成像方法、装置、设备及存储介质 |
CN116172612B (zh) * | 2023-02-02 | 2023-12-15 | 逸超医疗科技(北京)有限公司 | 一种基于超快平面波复合超声微细血流成像方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101214159A (zh) * | 2007-12-29 | 2008-07-09 | 中国人民解放军第四军医大学 | 利用超声多普勒信息显示血流图像的方法及装置 |
CN102697524A (zh) * | 2012-05-04 | 2012-10-03 | 成都优途科技有限公司 | 全聚焦超声成像方法及其在血流成像中的运用 |
CN103829974A (zh) * | 2012-11-27 | 2014-06-04 | 飞依诺科技(苏州)有限公司 | 能量多普勒血流复合成像方法及*** |
CN104398271A (zh) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | 血管与斑块的三维力学及组织特性成像检测方法 |
CN105615919A (zh) * | 2014-11-21 | 2016-06-01 | �田睦 | 快速二维多普勒速度和方向成像 |
CN106580369A (zh) * | 2016-11-30 | 2017-04-26 | 珠海威泓医疗科技有限公司 | 一种多角度平面波相干彩色多普勒成像方法 |
CN108272439A (zh) * | 2017-01-06 | 2018-07-13 | 佳能株式会社 | 被检体信息处理装置和图像显示方法 |
CN108852415A (zh) * | 2018-05-07 | 2018-11-23 | 深圳市德力凯医疗设备股份有限公司 | 一种经颅三维脑血管复合成像方法及*** |
CN109157247A (zh) * | 2018-08-20 | 2019-01-08 | 云南大学 | 一种最优帧间隔追踪估计的超声血流速度场测量 |
CN109416400A (zh) * | 2016-06-23 | 2019-03-01 | 皇家飞利浦有限公司 | 具有大线性阵列的快速合成聚焦超声成像 |
EP3513735A1 (en) * | 2017-04-25 | 2019-07-24 | Sogang University Research Foundation | Device and method for generating ultrasound vector doppler image using plane wave synthesis |
CN110431444A (zh) * | 2017-01-16 | 2019-11-08 | 皇家飞利浦有限公司 | 用于超快成像的***和方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DK2940487T3 (en) * | 2011-10-19 | 2016-11-21 | Verasonics Inc | Estimation and display for vector Doppler imaging using flat wave transmissions |
CN110811686B (zh) * | 2015-06-05 | 2022-08-12 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体成像方法及超声流体成像*** |
CN106725599B (zh) * | 2016-12-30 | 2019-12-20 | 飞依诺科技(苏州)有限公司 | 一种基于平面波的成像方法及装置 |
-
2019
- 2019-12-02 CN CN201911213180.0A patent/CN110811688B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101214159A (zh) * | 2007-12-29 | 2008-07-09 | 中国人民解放军第四军医大学 | 利用超声多普勒信息显示血流图像的方法及装置 |
CN102697524A (zh) * | 2012-05-04 | 2012-10-03 | 成都优途科技有限公司 | 全聚焦超声成像方法及其在血流成像中的运用 |
CN103829974A (zh) * | 2012-11-27 | 2014-06-04 | 飞依诺科技(苏州)有限公司 | 能量多普勒血流复合成像方法及*** |
CN104398271A (zh) * | 2014-11-14 | 2015-03-11 | 西安交通大学 | 血管与斑块的三维力学及组织特性成像检测方法 |
CN105615919A (zh) * | 2014-11-21 | 2016-06-01 | �田睦 | 快速二维多普勒速度和方向成像 |
CN109416400A (zh) * | 2016-06-23 | 2019-03-01 | 皇家飞利浦有限公司 | 具有大线性阵列的快速合成聚焦超声成像 |
CN106580369A (zh) * | 2016-11-30 | 2017-04-26 | 珠海威泓医疗科技有限公司 | 一种多角度平面波相干彩色多普勒成像方法 |
CN108272439A (zh) * | 2017-01-06 | 2018-07-13 | 佳能株式会社 | 被检体信息处理装置和图像显示方法 |
CN110431444A (zh) * | 2017-01-16 | 2019-11-08 | 皇家飞利浦有限公司 | 用于超快成像的***和方法 |
EP3513735A1 (en) * | 2017-04-25 | 2019-07-24 | Sogang University Research Foundation | Device and method for generating ultrasound vector doppler image using plane wave synthesis |
CN108852415A (zh) * | 2018-05-07 | 2018-11-23 | 深圳市德力凯医疗设备股份有限公司 | 一种经颅三维脑血管复合成像方法及*** |
CN109157247A (zh) * | 2018-08-20 | 2019-01-08 | 云南大学 | 一种最优帧间隔追踪估计的超声血流速度场测量 |
Non-Patent Citations (3)
Title |
---|
"A Comparison Between Compounding Techniques Using Large Beam-Steered Plane Wave Imaging for Blood Vector Velocity Imaging in a Carotid Artery Model;Saris Anne E C M,等;《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL》;20161130;第63卷(第11期);第1758-1771页 * |
High frame-rate vector flow estimation using speckle tracking with recursive plane-wave compounding;Che-Chou Shen,等;《2010 IEEE International Ultrasonics Symposium, San Diego》;20101231;第1307-1310页 * |
High Frame-Rate Vector Flow Estimation Using Speckle Tracking with Recursive Plane-Wave Compounding;Che-Chou Shen,等;《2010 IEEE International Ultrasonics Symposium》;20101231;第1307-1310页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110811688A (zh) | 2020-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110811688B (zh) | 多角度平面波重复复合的超快超声多普勒血流估计方法 | |
US11944497B2 (en) | Ultrasonic blood flow imaging display method and ultrasonic imaging system | |
US20180192989A1 (en) | Methods and systems for ultrasonic imaging | |
US8539838B2 (en) | Extended field of view ultrasonic imaging with a two dimensional array probe | |
US8469887B2 (en) | Method and apparatus for flow parameter imaging | |
US20130197365A1 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic apparatus control method | |
US20070038091A1 (en) | Ultrasonic diagnostic apparatus and ultrasonic transmission method | |
US9474510B2 (en) | Ultrasound and system for forming an ultrasound image | |
US10799208B2 (en) | Compressional sound speed imaging using ultrasound | |
US8306296B2 (en) | Clutter signal filtering using eigenvectors in an ultrasound system | |
WO2021179573A1 (zh) | 一种超声多普勒血流成像方法及*** | |
US20140121519A1 (en) | Ultrasonic diagnostic apparatus and ultrasonic diagnostic apparatus control method | |
EP2609866B1 (en) | Providing motion mode image in ultrasound system | |
EP2610639A2 (en) | Estimating motion of particle based on vector doppler in ultrasound system | |
KR101652727B1 (ko) | 초음파 도플러 영상 장치 및 그 제어 방법 | |
EP2462872B1 (en) | Providing additional information corresponding to change of blood flow with a time in ultrasound system | |
CN108113702A (zh) | 一种超声波扫查方法和装置 | |
US7371219B2 (en) | Ultrasound diagnosis apparatus operable in doppler mode | |
Ricci et al. | Multi-line measurements of blood velocity vectors in real-time | |
US20080030581A1 (en) | Multizone Color Doppler Beam Transmission Method | |
WO2021108953A1 (zh) | 多角度平面波重复复合的超快超声多普勒血流估计方法 | |
US20220287685A1 (en) | Method and system for estimating motion from overlapping multiline acquisitions of successive ultrasound transmit events | |
US20220183655A1 (en) | Three dimensional volume flow quantification and measurement | |
Chen et al. | A 3D motion compensation method for high frame rate volumetric ultrasound imaging based on velocity vector estimation: a simulation study | |
US20240111046A1 (en) | Method and system for flow processing on channel data for application of nonlinear beamforming |
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 |