CN102123668B - 使用未聚焦发送波束的高帧率定量多普勒流成像 - Google Patents
使用未聚焦发送波束的高帧率定量多普勒流成像 Download PDFInfo
- Publication number
- CN102123668B CN102123668B CN200980131855.5A CN200980131855A CN102123668B CN 102123668 B CN102123668 B CN 102123668B CN 200980131855 A CN200980131855 A CN 200980131855A CN 102123668 B CN102123668 B CN 102123668B
- Authority
- CN
- China
- Prior art keywords
- doppler
- image
- angle
- parameter
- medium
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
- G01S7/52028—Extracting wanted echo signals using digital techniques
-
- 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/13—Tomography
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52034—Data rate converters
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
- G01S15/8984—Measuring the velocity vector
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8995—Combining images from different aspect angles, e.g. spatial compounding
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52047—Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
提供了具有面向像素的处理的超声成像***,其中,通过以下步骤实现了产生多普勒速度图像的方法:将未聚焦声信号发送到基本整个场上方的介质内;响应于发射,在换能器阵列上接收散射的和反射的超声信号;处理所接收的超声信号,以提取信息来构造对应于所述介质中至少一个点的多普勒速度信号;以及根据所处理的多普勒速度信号,在显示设备上产生多普勒速度图像。描述了采集序列和信号处理算法,其提供了流体流参数的改进量化,包括血流区和组织区之间的改进辨别。描述了用于对于大视野的实时和后采集可视化的很高帧率频谱多普勒和向量多普勒采集模式。
Description
相关申请的交叉引用
本申请基于35U.S.C.§119(e)要求于2008年6月26日提交的第61/076,057号美国临时专利申请的权益,该临时申请的全部内容通过引用并入本文。
技术领域
本公开针对这样的***和方法,该***和方法捕获和处理超声数据并根据超声数据产生代表流体流的图像。
背景技术
超声成像已经发展成诊断多种疾病状态和状况的有效工具。超声设备的市场多年来已经显现稳步增长,这由图像质量方面的改进和区分各种组织的能力推动。不幸的是,对于超声***,仍存在设备花费对于重大采用来说过高的许多应用。示例是下面的应用领域,如乳腺癌检测、***成像、肌肉骨骼成像和介入放射学。在这些和其它领域中,超声成像的诊断效力依赖于用于各种组织类型的区分和识别的优秀空间和对比分辨率。这些性能能力仅在具有更广泛处理能力的更昂贵超声***上得到。
超声成像一直需要大量的信号和图像处理方法,尤其对于采用多达128个或更多换能元件的阵列***,每个换能元件具有独特的信号处理要求。过去的十年已经在除市场最底层的那些***以外的几乎所有***中显现出向数字信号处理的改进精度和灵活性的过渡。通过使用高度集成的数字电路,该过渡从长远角度看具有降低***成本的可能性。不幸的是,超声***的低制造量导致用于这些独特电路的大量开支和固定花费,因而向数字信号处理的过渡未显著降低***成本。
医学超声中的多普勒方法包括用于对血流进行成像和量化的多项相关技术。对于静止目标,从目标反射回换能器的脉冲的往返行程时间对于每次传输是相同的。相反,来自移动对象的连续的回声影像(echographic)返回将关于发送脉冲在不同时间到达,并且通过使这些回声互相关,可以估计对象速度。由于超声路径是定向的(沿着波束轴线),所以仅轴向运动产生多普勒信号。对于波束横向的流是不可检测的,因此传统多普勒方法中获得的速度幅度仅代表流速向量的轴向分量。为了估计流速向量的真实幅度,采用了向量多普勒方法。通常,这些方法依赖于多波束角数据来估计流向量和流速向量的方向。
已经发展了几种基于多普勒的方法来呈现血流的不同方面。通常,流场的“空间成像”用于定位血管(以测量血管尺寸)和观察流结构。“流成像”与回声影像成像结合用于“双功”模式,该“双功”模式将覆盖图中的两类图像与以灰阶呈现的回声影像振幅和以彩色绘制的流速相结合。由于流成像在采集时间和处理载荷两方面要求更高,所以在感兴趣区(ROI)内计算流场,该感兴趣区是较大回声影像图像的子集。
在选自ROI内的小得多的取样容积内,流速的详细量化是可能的。能够单独地采样和处理的最小容量由轴向长度(发送脉冲长度)和侧向波束宽度(成像面之内和之外)给出。任何方法的空间分辨率取决于取样容积的大小,还取决于对于那个位置的***灵敏度设置。
频谱多普勒方法报告流速的频谱以及在心动周期(cardiaccycle)期间流速的频谱如何变化,并且频谱多普勒方法通常将频谱通过图形呈现为频谱图并通过扬声器可听见地呈现频谱。此外,频谱多普勒方法对关于一系列发送获得的流速功率谱进行计算,并且通常将频谱通过图形呈现为频谱图并通过扬声器可听见地呈现频谱。对血液速度的全时变频谱的访问允许在取样区内准确地计算平均流速和峰值流速,并提供全部超声多普勒方法的流分布的最完整特征描述。
感兴趣区内速度场的彩色流多普勒(ColorFlowDoppler)成像是这样一种方法,该方法使用调色板呈现流并且通过使用暖(略带红色的)色调和冷(略带蓝色的)色调区分不同流向(通常朝向换能器或离开换能器),调色板通常将较高速度比较低速度渲染得更明亮。不对非常慢地移动的区域和静止区进行染色,并且使用“壁滤波器(wallfilter)”来设置最小截止速度。彩色流多普勒可以在感兴趣区内提供近似平均流速,但是由于保持合理的帧率所需要的短采集序列的原因,所以精度是有限的。
彩色流多普勒需要采集紧邻的一系列相同发送—接收事件或者“系综(ensemble)”,以通过多种方法检测和量化运动,本质上查找在信号的到达时间或相位方面的相关差异。脉冲重复频率(PRF)能够与由从换能器到图像的最大深度并再次返回的声音往返行程时间所允许的一样快,但是通常将其调整至在没有混淆的情况下显现峰值血液速度所允许的最小值。在ROI中对于每条多普勒扫描线通常使用具有8到16个之间的脉冲—回声事件的系综。发送波束焦点参数的选择通常导致比用于回声影像成像的那些参数宽2至3倍的多普勒扫描线。在每个波束方向上发送一整套脉冲的需要通常导致对于彩色流多普勒比对于回声影像成像更慢的帧率。来自慢帧率的伪影(artifact)在多普勒成像中可能经常比在灰阶回声影像中更明显,这是因为在心动周期的一部分期间,流中的显著变化可能发生,并且甚至轻微的探测运动可能在全部ROI上导致明显的流。
使用小的ROI可能提高帧率,但是可能限制流异常情况的评定。例如,使用10条多普勒线和12个脉冲的系综的彩色流ROI需要120个事件,这与全帧回声影像图像类似。
通常,高质量多普勒成像在技术上比回声影像成像在很大程度上更困难,这是因为来自血液的反向散射与组织相比十分弱。对产生整齐的且无伪影的彩色流图像的公知基本挑战包括:
·对高度可重复的发送脉冲的需求和采集硬件中很低的噪声和相位抖动。
·流信号经常具有与各种噪声源相同的数量级,但是平均值对帧率和其它运动伪影(motionartifact)具有不利影响。
·组织的散射振幅和血液的散射振幅之间大对比度导致辨别血管壁(强回声)和移动血液(弱回声)的困难,甚至在速度对比度高时。此外,血流速度在血管壁附近经常很慢,其经常与心动周期同步移动(脉动)。
·多普勒脉冲通常长于回声影像脉冲,并且必须小心来记录具有不同分辨率的流图像和回声图像。对于小的血管壁,这是特别有挑战性的,因为多普勒脉冲的取样容积可能大于血管壁直径。
发明内容
已经发展了多种方法来解决这些问题,多种方法既记载在技术著作中又记载在在先专利中。本文描述的实施方式基于申请人的在先元件级超声数据的基于像素的处理和使用未聚焦发送波束来提高帧率,该元件级超声数据的基于像素的处理是共同待审的第11/911,633号美国专利申请的主题。对于使用频谱多普勒处理的定量成像,描述了未聚焦发送(例如,平面波)的使用。所公开的实施方式描述了用于实时流和运动量化的新方法和用于通过使用所记录的高PRF数据的后处理来发展新成像模式的新方法。
依照一个实施方式,提供了产生多普勒速度图像的方法,其包括:将未聚焦声信号发送到基本整个场上方的介质内;响应于发射,在换能器阵列上接收散射的和反射的超声信号;处理所接收的超声信号,以提取信息来构造对应于所述介质中至少一个点的多普勒速度信号;以及根据所处理的多普勒速度信号,在显示设备上产生多普勒速度图像。
依照本公开的另一实施方式,提供了包括超声处理方法,其包括:产生未聚焦声信号;在多个接收元件处接收所述未聚焦的声信号的散射回声和反射回声并从散射回声和反射回声中获得多普勒速度回声信号;将给定像素映射到所述多普勒速度回声信号的区域内;对于所述给定像素,将所存储的多普勒速度回声信号的映射区组织成阵列;处理所述阵列以对于所述给定像素产生信号响应;以及使用所述信号响应来对于所述给定像素获得多普勒速度声信息。
依照前述实施方式的另一方面,所述方法包括产生一组给定像素的初始步骤,该组给定像素被选择成代表产生所述多普勒速度声信号的换能器的视野中的区域,其中,组中的每个给定像素对于所述多个接收元件具有已知的空间关系。
本公开的公开实施方式还指向超声成像方法和***,该超声成像方法和***在商业CPU上执行的软件中执行所有信号处理和图像形成。本方法中需要的唯一定制硬件是用于声脉冲的发送、数据采集和来自换能器的所接收信号的信号调节。作为重要的益处,新结构允许***动态范围方面的改进和新的采集模式,***动态范围方面的改进为在低成本扫描头设计中新换能器材料的利用做准备,新的采集模式提供显著新的诊断信息。
所公开的基于软件的超声***结构通过使设计基于商业计算机主板促进了来自计算机业的高容量低成本处理技术。虽然一些目前的超声***在它们的设计中包含计算机主板,但是计算机仅用于用户界面和某一***控制并且不参加任何实时处理任务。在所公开的结构中,计算机主板替换了几乎所有现有硬件,而不是补充它。使软件中的***基于通用平台以最低可能的***成本提供了灵活的高性能的成像***。该方法不需要定制集成电路,降低了***复杂度和上市时间。此外,随着计算机业实现CPU处理能力的进一步提升,它们可以被***容易地采用,以增强成像性能或者提供新的操作模式和信息提取。
依照面向像素的处理的一个实施方式,步骤包括:产生声信号;在多个接收元件处接收所述声信号的至少一个回声并根据所述声信号的至少一个回声获得回声信号;存储来自于每个所述多个接收元件的每个回声信号;将给定像素映射到所存储的回声信号的区域内;对于给定像素将所存储的回声信号的映射区组织成阵列;处理阵列以对于所述给定像素产生信号响应;以及对于所述给定像素使用所述信号响应来获得声信息。
依照前述实施方式的另一方面,提供了包括产生一组给定像素的初始步骤,该组给定像素被选择成代表产生所述声信号的换能器的视野中的区域,其中,阵列组中的均匀给定像素对于所述多个接收元件具有已知的空间关系。优选地所述方法还包括对于阵列中的给定像素根据所述声信息产生图像。
依照前述实施方式的另一方面,所述声信息可以用于下面中的一个或多个,包括但不限于:测量并显示空间数据;测量并显示时间数据;测量并显示血流数据;以及测量并显示对由声信号或声发送波导致的诱发机械位移做出响应的组织位移。
依照前述实施方式的另一方面,所述方法包括:产生多个声信号;从所述多个声信号中接收回声;在多个产生和接收期间合并所接收的回声,以增强从所接收的回声中获得的声信息。
依照前述实施方式的另一方面,合并和平均所存储的回声信号。而且,信号响应包括所存储的回声信号的平均。
依照前述实施方式的另一方面,所述方法包括以下步骤的多次循环的合并结果:产生声信号,接收回声并且从所接收的回声中获得回声信号以得到增强的声信息。
依照前述实施方式的另一方面,所述增强的声信息包括空间复合,该空间复合提高根据所述增强的声信息产生的最终图像的对比分辨率。
依照前述实施方式的另一方面,所合并的信号代表与移动组织或移动血细胞关联的多普勒信息。
依照前述实施方式的另一方面,以比处理所述阵列的速率高的速率进行回声信号的接收、获取和存储。
依照本公开的另一实施方式,提供了包括超声处理方法,其包括:产生声信号;在多个接收元件处接收所述声信号的至少一个回声并根据所述声信号的至少一个回声获得回声信号;存储来自于每个所述多个接收元件的每个回声信号;将给定三维像素映射到所存储的回声信号的区域内;对于所述给定三维像素将所存储的回声信号的映射区组织成阵列;处理阵列以对于所述给定三维像素产生信号响应;以及对于所述给定三维像素使用所述信号响应来获得三维声信息。
依照前述实施方式的另一方面,关于上述第一实施方式的所有方面可应用于本公开的第二实施方式。
依照本公开的另一实施方式,提供了处理声回声的方法,其包括:存储从多个接收元件中接收的声回声信号;将给定像素映射到所存储的回声信号的区域内;对于所述给定像素将所存储的回声信号的映射区组织成阵列;对所述阵列执行运算来对于所述给定像素产生信号响应;以及对于所述给定像素使用所述信号响应来获得声信息。
依照本公开的另一实施方式,提供了超声处理***,其包括:模块,该模块适合于产生声信号,在所述模块中的多个接收元件处接收所述声信号的至少一个回声以及根据所述声信号的至少一个回声获得多个回声信号;和用于处理的装置,其与所述模块通信并适于将给定像素映射到从所述模块中接收的所存储的回声信号的区域内,以对于所述给定像素将所存储的回声信号的映射区组织成阵列,对所述阵列执行运算来对于所述给定像素产生信号响应,以及对于所述给定像素使用所述信号响应来获得声信息。
依照前述实施方式的另一方面,处理装置适于最初产生一组给定像素,其中,对于所述模块中的接收元件,组中的每个给定像素具有已知的空间关系。理想地,将处理装置配置成在所述阵列中对于所述给定像素根据声信息产生图像。替换地或者与之结合,提供了用于显示图像的装置,其从处理装置接收信号响应,以在计算机显示屏上或者以打印的形式或者以本领域技术人员已知的其它方式产生图像。
依照本公开的另一实施方式,提供了超声处理***,其包括:模块,该模块适合于产生声信号,在所述模块中的多个接收元件处接收所述声信号的至少一个回声以及根据所述声信号的至少一个回声获得多个回声信号;和用于处理的装置,其与所述模块通信并适于将给定三维像素映射到从所述模块中接收的所存储的回声信号的区域内,以对于所述给定三维像素将所存储的回声信号的映射区组织成阵列,对所述阵列执行运算来对于所述给定三维像素产生信号响应,以及对于所述给定三维像素使用所述信号响应来获得声信息。
总之,对在商业可用计算平台上实现的基于软件的超声***结构的改变的益处包括:
—明显更低的硬件成本
—通过对于定制的集成电路(ASIC)避免延长的设计周期,较低的开发成本和较快的上市时间
—计算机技术上成本/性能提升的直接促进
—在商业环境和学术环境中,许多新处理方法的发展的灵活性
—对于成本敏感的应用领域,基于图像质量改进,提高的诊断能力
—在成本已经成为采用的障碍的专门应用中增加的超声应用。
附图说明
当结合下面附图时,因为从本公开的下面详细描述中本公开的前述和其它特征和优势变得更好理解,所以将更容易地体会到本公开的上述和其它特征和优势,其中:
图1示出了传统多普勒采集序列和图像空间;
图2示出了单角未聚焦发送多普勒采集序列和图像空间;
图3示出了被分解成轴向流分量和横向流分量的流速向量;
图4示出了传统超声***中用于交织回声影像、彩色流和频谱多普勒成像的采集序列;
图5示出了通过使用未聚焦平面波发送用于交织回波影像、彩色流和频谱多普勒成像的采集序列;
图6示出了作为多普勒频率或者相当地多普勒速度的函数的多普勒速度能量谱;
图7示出了用于本公开处理的***结构的高层图示;
图8是面向像素的处理的一个实施方式的基于软件结构的示意图示;
图9是依照面向像素的处理形成的插件模块的图;
图10是用于依照面向像素的处理形成的128元件线性阵列的采集数据的示意图示;以及
图11是面向像素的处理中使用的像素映射处理的图示。
具体实施方式
图1中,对于由N(8≤N≤16)个发送—接收事件的系综组成的单波束方向示出了传统多普勒采集序列和图像。通过使用延迟相加波束形成和形成图像所需的几个其它处理步骤,传统的超声***形成聚焦的发送波束并动态聚焦所接收的数据。在申请人的标题为“UltrasoundImagingSystemwithPixelOrientedProcessing(具有面向像素处理的超声成像***)”的第11/911,633号美国专利申请中,下面结合图7—图10描述了图像重构方法,该图像重构方法与传统的波束形成相比大大降低了处理负荷并允许使用多种非传统发送场。
例如,一个非传统发送场是平焦发送模式,在平焦发送模式中,使所有换能器元件同相发射以产生一部分平面波(对于线性阵列),该部分平面波可以用于通过单个脉冲声穿透整个视野,从而获得极快的帧率。一致的相位发送对于线性阵列产生平焦发送脉冲或平面波发送脉冲。
这在图1中示出,其中,换能器元件或头20包括被布置为线性阵列的多个单独的发送器/接收器22。从换能器头20侧面延伸的实心平行线示出了从换能器头20向在组织30中支撑的血管28发射的能量的二维额定边界24。额定边界代表线性阵列的成像空间。虚线32示出了又一二维边界,虚线32在其从侧面26向血管28的方向延伸时发散,代表用于成像和多普勒用途的有用能量的外边界。如图1中所示,血管28运送流体,流体由大体上指向图1右侧的流速向量箭头34指示。
图像中央的矩形框36限定了图1中的具体感兴趣区(ROI),感兴趣区ROI与血管28的中央部分重叠,在血管28的中央部分中,血管具有双路分支31。ROI框36具有与各个换能器元件20的端部对准的平行竖直线38。为了在ROI内产生血管28的图像,选择与ROI框36竖直对准的换能器元件20来发送波束。所选择的换能器元件限定发送子孔径40,当激活所选择的发送器22时,从发送子孔径40中形成了聚焦的发送波束42。
发送波束42的形状由从发送子孔径40中的每个发送器22中得到的合并波形限定,合并波形以额定扫描线44的纵轴线为中心。发送波束42的反向散射导致返回至换能器头20中的全部发送器/接收器22的反射波,换能器头20中的全部发射器/接收器22共同限定了接收孔径46,如图1中所示。
B型(B-mode)中的序列操作由图1底部的脉冲图连同N个用于彩色流成像的脉冲系综的图示出。由于传统的B型和彩色流成像系综是本领域技术人员公知的,所以这些在此将不进行详细描述。简单地说,脉冲图上示出的每个脉冲代表由脉冲重复间隔(PRI)分离的每条扫描线38上的脉冲重复。在彩色流成像中,对于每条扫描线38存在N个脉冲的系综。
与图1中示出的聚焦发送波束42对照,用一致相位激发的曲线阵列将产生具有弯曲波前的脉冲,这产生似乎从阵列的曲率中心发出的一部分球面波。这还可以由以下步骤完成:首先激活线性阵列的中央换能器,接着顺序地逐渐向外激活相邻的多对换能器至阵列端部。
还能够以直至由往返行程时间允许的最大值的帧率采集需要N个脉冲系综的多普勒帧,该帧率由PRF/N(通常,1kHz≤PRFmax≤12kHz)给出。在系综期间接收的N个数据组通过使用下面描述的基于像素的方法重构,并且在由平焦发送场声穿透的整个视野上通过使用用于(轴向)流速和能量估计的传统互相关多普勒方法进一步处理。生成的多普勒帧率很快,这允许在产生单条传统多普勒扫描线所需的时间内采集(并处理)全帧多普勒彩色流图像。这在快速变化的高速流的成像中是特别有用的。而且,声平面波的使用将在单次检测中探测流场,因而降低了由来自传统聚焦发送波束的广角声穿透导致的速度扩展。
改进的测量精度
借助未聚焦的发送多普勒成像方法,对于整个帧仅需单个多普勒系综采集。这允许使用比借助传统多发送波束方法可能的多普勒系综长得多的多普勒系综,而仍支持高帧率。传统多普勒血流成像方法对于全帧流图像使用多至128个系综,因而必须将最大系综长度限制至小于N=16个脉冲,以便防止显著影响帧率。借助仅使用单个系综的未聚焦发送多普勒方法,N的值可能比16大很多倍,而仍允许比传统方法高很多的采集帧率。较长的未聚焦发送系综允许在血液速度估计中提高的精度,因为多普勒频率估计(从中得到血液速度)中的不确定度约是系综总时间的倒数。
例如,具有5KHz的PRF的系综长度10将具有2毫秒的总采集时间,这导致500Hz的频率不确定度。在3MHz的典型换能器频率下,这将转换成约13厘米/秒(cm/sec)的血液速度不确定度,对于在血流量化方面的任何尝试的明显误差。此外,借助与10毫秒回声采集周期结合的典型64系综传统流图像,该示例的帧率将小于8帧每秒。借助未聚焦发送多普勒方法,可以使用100个脉冲的系综长度,这提供了20毫秒采集周期和仅1.3厘米/秒的血液速度不确定度。对于帧的回声成像部分再次假设10毫秒采集周期,帧率约是33帧每秒。
多角未聚焦发送多普勒
如在回声影像成像情况下一样,在所发送的场中缺少聚焦导致比用聚焦发送波束可获得的旁瓣干扰明显更大的旁瓣干扰(因而更差的侧向分辨率)。熟知的是,发送—接收事件的合并波束图由发送波束图和接收波束图的产物给出;并且由于平焦发送具有一致的样式,所以其不提供聚焦增益(因此没有对侧向分辨率的贡献)。通过对于几个不同平面波方向(即,对于关于换能器面以不同角发射的平面波)合并平焦系综,可以大大降低使用平面波发送场获得的旁瓣电平。(通过使线性阵列相位同步来倾斜平面波的角相当于移动曲线阵列的表面曲率中心并产生点源的合成阵列。)
由于平焦波前在整个深度上几乎保持平坦,所以在整个视野上旁瓣抑制几乎是均匀的。即使少至五个不同平面波角结合来在整个图像中提供好的旁瓣抑制和侧向分辨率,也能够在比在较小ROI内部对血流进行成像的传统多普勒彩色流图像帧通常需要的时间更少的时间内获得五个系综。
本公开的方法使用这样的算法,该算法采用为灰阶回声影像成像开发的多角方法来收集多角多普勒数据。由于在每个系综中存在N个脉冲并存在M个角,所以可以以两种方式收集数据:(a)在一个角下收集全部N个脉冲,然后改变角直至全部角完成为止,或者(b)对于M个角中的每一个收集一个脉冲,然后对于N个脉冲重复。基于最大预期流速做出选择,最大预期流速转而确定对最小PRF的需求,以避免混叠和其它伪影。
图2示出了本公开的单角未聚焦线性多普勒序列。这里,换能器头20具有被同时激活的所有单独换能器元件22,所有单独换能器元件22产生到组织30内的未聚焦或“平焦”波前48。在图2底部绘制的脉冲以B型示出了单个角或多个角并且对于色彩模式示出了系综,其中平焦的M比聚焦波束的M大得多。换能器头20从未聚焦波前48接收反向散射的信号,将模拟信号发送到处理器以进行存储和处理。
更具体地,一旦已收集数据,方法继续进行到使用至少两种可能方法中的一种或多种从ROI中重构图像。第一种对于每个脉冲重构多角数据,然后使用多普勒互相关来处理重构系综。该方法在多普勒图像中产生最佳侧向分辨率,但是不保留向量流信息,并且产生非常像聚焦发送波束在传统***中执行的谱展宽。可选地,在第二方法中,可以如下所述对于每个角进行估计并然后关于每个角向量合并多普勒速度和能量。两种处理方法的组合也是可行的。
向量多普勒和流检测
传统多普勒实现仅提供轴向流信息的基本原因是实际的:在一个方向上用发送波束(扫描线)覆盖ROI所需的时间已经处于临床应用的极限。在不由于改变流状况和非故意探测运动而降低帧率和引入伪影和误差的情况下,以不同发送角添加波束将是不可能的。
可以使用新的多角平焦多普勒方法来以高帧率获得向量流信息,因为使用不同的发送波束探测全部流场,每条发送波束在独特的方向上传播。通过利用轴向流幅度和波束轴线和流向量之间的角(多普勒角)之间的已知关系,可以合并数据,以单独地估计流方向和流速。此外,虽然在每个系综中噪声将显现为流,但是其可能与角完全不相关。因而,多角未聚焦线性多普勒将通过使用定向信息既提供良好侧向分辨率又提供噪声抑制。当然,还可以使用向量多普勒信息来提供流的向量显示(例如,流线图),并且最重要地,可以客观地测量绝对速度幅度(即,没有专门的角修正),只要流位于图像平面内。在该实施例中,图像平面是关于以直角与扫描线38相交的平面。在线性阵列情况下,图像平面将以直角与突出穿过每个发送器—接收器22中央并且垂直于每个发送器—接收器22表面的线相交。该平面代表正由超声***成像的空间中的二维区域。
图3示出了被分解成轴向流分量和横向流分量的流向量。竖直虚线44是超声发送的轴向,水平虚线50是横向流向。将血管28内的血流示出为拉长的箭头52。
必须从开始保存向量流信息,以便利用通过使用多角发送来辨别噪声和流的能力。可以沿着为多角采集而描绘的图案进行数据收集。可以将上面概括的处理步骤合并在矩阵表达中,使得新处理算法一起处理多普勒数据记录的N×M矩阵,提高侧向分辨率,辨别组织区和血流区,降低噪声并以快的帧率提供向量流信息。
流速向量的幅度和方向二者的信息帮助改进真实流和噪声之间的辨别以及慢流体流和可能也在运动(“壁运动”)中的组织之间的辨别。例如,慢流信号幅度经常接近***本底噪声,并且必须选择多普勒速度幅度阈值足够高于本底噪声,以防止流显示的污染。流向估计也是有噪声的;然而,真实流向是稳定的(在用于进行重复测量的小时间间隔上),而噪声方向是随机的,具有均值零,并且平均几个测量降低噪声并相干地合并流信号。向量信息还允许使用包括相邻取样容积在内的其它流相干滤波器来提高整个图像的SNR。
众所周知,动脉经常显现与心动周期一致的脉动。通过向量方向估计来改进壁运动和近壁流之间的辨别,因为壁运动对于血管壁轴线主要是横向的,而流通常是纵向的。因此,运动方向上的尖锐不连续可以用来扩增血管壁和腔之间辨别的其它方式。
彩色能量多普勒
传统互相关处理产生多普勒速度和多普勒能量的估计。后一量通常对于流更敏感,可以用来检测和绘制小血管。多角平焦采集和向量处理的噪声降低益处和侧向分辨率增强益处还延伸至彩色能量多普勒模式。
对于本文描述的实施方式,存在多项优势,包括但不限于:
(a)未聚焦发送多普勒血流成像在高帧率下提供全帧流图像。单个系综足以在整个图像空间上测量流,因而避免了ROI选择的过程和多普勒区大小和帧率之间的折中。
(b)用未聚焦发送多普勒方法可实现的较长系综长度在没有显著帧率降低的情况下提供了改进的血液速度测量精度。
(c)单角发送(使用线性换能器的平面波平焦)提供窄角激励并且降低***固有的谱展宽。
(d)多角未聚焦发送(如线性发送)允许整个图像上的高帧率向量多普勒测量。
(e)多角多普勒测量通过使用多个度量(例如,速度估计的方差、向量方向和振幅的角拟合的均方误差、或者多普勒能量的多角估计)允许流区和无流区之间改进的辨别。
将该方法扩展至使用被引导至在不同方向上传播的未聚焦发送平面波的组合并且在整个视野上提高彩色流图像的侧向分辨率。该方法能够适于使用被引导至在不同方向上传播的未聚焦发送平面波的组合来在整个视野上获得向量流方向和幅度,并且适于用标准换能器这样做。流必须在图像的平面内,以提供绝对流幅度。
通过使用上面引用的向量流信息和多普勒角和多普勒速度幅度之间的关系,该方法还能够适于使用被引导至在不同方向上传播的未聚焦发送平面波的组合来提高真实流和噪声之间的辨别。此外,这些方法能够适于曲线阵列(圆形波前,使表面曲率中心移位以创造点源的合成阵列)。
具有任何通常几何形状的换能器阵列(包括“相控阵”或“扇面阵列”、“球面阵列”、“2D阵列”)能够适于产生直线波前、环状波前、平面波前、球面波前或其它波前,这些波前产生依照本公开实施方式的“角分集”。
本公开的实施方式还延伸至基于像素的处理的软件实现,以包括多普勒和向量多普勒处理。本公开包括基于像素的多普勒和向量多普勒处理(例如,FPGA、ASIC)的硬件实现以及传统接收波束成形处理的硬件实现的使用,以包括并适应使用本文描述的平面波和其它未聚焦波束和采集序列和处理方法。前述内容还适用于彩色能量多普勒处理。
高帧率全场频谱多普勒
用于在单个图像点处量化流的传统频谱多普勒采集序列可以交织三个模式:(a)回声影像发送—接收序列60,(b)在包含该图像点的回声影像帧内对于跨越感兴趣区(ROI)的彩色线的彩色流系综62,和(c)使用图4中示出的单个聚焦发送波束的较长高PRF序列64。该三个模式经常被称为“三重模式”,三个采集模式60、62、64的该交织提供了目标点对于流的量化保持在期望位置上的信心,并且其还限制了由采集信息需要的大量时间导致的总体帧率。
目前使用中的超声***全部使用聚焦发送波束来进行三重模式扫描,如图1中所示。该方法提供了良好的灵敏度和空间分辨率,但是其花费相当多时间来扫描换能器的成像场内的感兴趣区。平面波发送场对三重扫描的应用允许对回声影像模式和彩色流模式两者的非常高帧率成像,因而最小化了为对准频谱多普勒而分配的采集时间的开销。而且,由于对准不需要与诊断扫描需要的图像质量一样好的图像质量,所以其甚至可能不需要使用多角平面波序列来提高图像分辨率和对比度。图4的三重模式成像实施例是通过使用图5中的未聚焦平面波发送呈现的。注意,频谱多普勒发送波束可能对于最高的旁瓣抑制被紧紧地聚焦在取样点上,或者可能使用未聚焦平面波。后者有具有单个多普勒角的优势(聚焦波束具有多普勒角的展开),但是导致稍大的取样容积。
在图4中示出的实施例中,对于回声影像成像过程60的一个帧,使用了128个发送—接收事件,对于彩色流帧(系综长度N=12)62使用了144个发送—接收事件,并且M=128的系综用于获得多普勒频谱64。该频谱被显示为一列像素(未示出),在该列像素中,像素亮度代表那个频率(速度)间隔内的多普勒能量。在该实施例中,分配给频谱多普勒的时间比率小于全部帧采集时间的三分之一(32%)。每一帧计算新频谱并将新频谱显示为在滚动显示中的另一竖直像素带。
在图5中,示出了用于通过使用未聚焦平面波发送对回声影像60方法、彩色流62方法和频谱多普勒64方法进行交织的采集序列。通过替换图4实施例中的聚焦发送,对于回声影像成像60的一个帧仅需要一个发送—接收事件,对于彩色流帧成像62,仅使用12个发送—接收事件的一个系综,并且128个用于获得多普勒频谱图像64。分配给频谱多普勒的时间比率是全部的90%。
多点频谱多普勒
在使用未聚焦发送波过程中固有的优势是能够立刻声穿透全部回声影像区,因而允许对图像空间中的任何点应用频谱多普勒处理。使用多栅取样的传统***还可以提供几个点处的频谱,但是这些通过适当的PRF和帧率的实际考虑被显示与沿着单个波束线的轴线上的点。未聚焦发送允许图像中任何地方的多个点处的流之间的量化比较。通过使用完整频谱或者对于在对采集PRF没有影响的情况下在任何数量的图像点处在心动周期期间跟踪的单个频谱参数(如峰值速度),可以进行这种比较。给定未聚焦平面波发送的宽范围,保持具有多个图像点处频谱参数的实时显示的快速PRF仅由数据处理和显示的速度限制。
在没有PRF或图像帧率的折中的情况下通过超声图像中的多个点处的频谱多普勒处理提供同时量化流信息的能力,提供了具有减少检查时间的复杂流异常的改进诊断。
聚焦的发送波束可以用于减少旁瓣。如果感兴趣的目标点位于单波束方向附近,则自适应算法可以自动形成聚焦发送波束,该聚焦发送波束不比声穿透目标点必需的波束宽。因而,可以调整发送波束以降低旁瓣杂乱,同时声穿透期望区域。可以提供用户界面控制来调整发送波束宽度,以实时评定宽度对图像质量的影响。如下面更全面讨论的,如果使用了宽波束(弱聚焦的)声穿透,则可以执行数据的后处理,以产生当前不可得到的新的流显示。
存储的高PRF数据的后处理
存储高PRF频谱数据的长序列允许后处理,以对声穿透区域中任一或全部点处的流进行量化。对于未聚焦平面波发送,对于未聚焦平面波发送,这种后处理将提供崭新的量化流和组织运动图像。可以对持续几个心动周期的高PRF数据记录进行后处理,以产生从频谱多普勒信息中得到的量化流参数的颜色叠加。与流参数基于系综数据的彩色多普勒图像相比,频谱多普勒处理提供高度准确的流信息。
在图6中,图解了从多普勒谱中获得的临床感兴趣的几个参数,并且可以对于每个图像帧计算这些参数中的每一个,以制作描绘不同流特征的影片。例如,峰值速度、对应于峰值多普勒能量的速度以及集成的频谱能量,均是这种新的频谱多普勒影片的可能参数选择。多普勒速度能量谱是多普勒频率或者相当的多普勒速度的函数。通过定位能量谱的模态获得峰值能量速度,并且类似地通过谱的第一矩给出平均速度。使用能量阈值电平限定峰值速度,低于该能量阈值电平,认为速度估计是不可靠的。类似地,谱宽度由速度之间的距离限定,对于该速度,能量已经降低至峰值能量电平以下的某一电平(例如,6dB)。壁滤波器截止用于从血流显示中消除组织运动。
可选择地,可以显示在完整心动周期期间代表流图的新的一类图像。例如,可以计算和绘制在一个或多个心动周期期间检测的最大流速,如最大谱宽度能够的(湍流的可能指示)。这种新模态的临床效用是未知的,但是有前景的,因为现在可以在整个流区域上产生目前仅在单点处可获得的信息。而且,以不通过检查声谱仪来预期的方式分析数据的能力可能对委派有远程地或者恰在检查发生之后审查数据的临床医生具有兴趣。本发明中预期了新检查和数据记录方案以及新后检查处理步骤的建立。
实现上述过程的硬件是独特的,因为其允许存储大量的非常高PRF的接收数据,持续几个心动周期。最大数据速率主要由通过PCIe总线的传输速率能力限制。当前***数据速率允许对于接收数据的64个通道高达1.5GB/s的传输,因而使高PRF多普勒超声数据能够连续流动。在用新数据覆盖连续流动数据之前连续流动数据可以被存储的最长时间主要由主计算机内存的容量限制。
可以实施执行多点频谱多普勒处理,提供了视野中各个空间点处的量化流信息,但是当前在信号处理带宽方面的限制可能限制提取的点的数量和/或信息的类型。存储所接收的数据并使用后处理来产生频谱多普勒信息克服了该处理限制并且能够产生如由用户选择的其它频谱参数的全新图像类型。
该方法的关键优势
通过使用未聚焦发送脉冲,在没有来自于必须发送多个波束的PRF损失的情况下对于图像空间中的每个点可获得频谱多普勒数据。因此,几个图像点可以通过使用临床准则来选择,并且对于相同的发送事件被处理和比较,从而使来自快速变化的流的采集伪影最小化。
未聚焦发送三重模式成像(回声影像、彩色流、频谱多普勒)以非常高的PRF提供了全帧图像和多普勒数据。单个彩色流系综足以在整个回声影像图像空间上使流成像,因而提供了更多用于频谱多普勒采集的时间。
单角发送(例如,使用线性换能器的平面波平焦)提供了窄角激励,并且降低了***固有的频谱展宽。
硬件能够存储大量持续许多心动周期的很高PRF的数据。未聚焦发送允许使用频谱多普勒方法处理数据,以获得图像空间中每个点处完整的频谱。然后通过使用刻画频谱特征的多个参数中的任何一个或者其随时间的变化,可以在每个点处量化流(和还可能的组织运动)。通过使用这些参数中的任一个来进行前所未有的流的量化和表示,可以产生图像。
未聚焦激励和长记录的高PRF数据的优势是在检查期间可以以声谱仪未预料的方式进行数据的回顾性分析,虽然通过使用实时显示可以保证数据质量。执行灵活回顾性分析的这种能力具有对于病人扫描的传统临床回顾和对于远程医疗的潜在应用。
还可以将本文描述的检测和量化流的所有方法应用于组织运动响应于心脏刺激或呼吸刺激或者响应于外部施加力的特征刻画。
图7是***层次框图,其表示用于实现本公开的过程的高层***结构70。应理解这仅是一个代表性的实施方式,并且所示出的结构70不是本公开的所有实施方式的必要条件。
结构70包括通过PCI快速总线(PCI-express)74连接到多通道收发器和数据采集***76的主计算机72。主计算机72具有用户界面和控制78以及显示屏80,两者连接到处理器82,处理器82利用基于像素的应用处理软件84。多通道收发器和数据采集***76硬件连接到超声换能器86,超声换能器86用于对于声介质90中的区域88进行成像。由于这些组件是容易商业获得的,所以本文中将不详细描述它们。
面向像素的处理
依照本公开的一个实施方式的基于软件的方法和***结构实现软件中的全部实时处理功能。所提出的结构示意性地示出在图8中。
基于软件的***中的唯一定制硬件组件是计算机的扩展总线的插件模块,该计算机包含脉冲发生和信号采集电路以及用于存储信号数据的一大块扩展内存。信号采集过程包括对在发送脉冲以后从换能器元件中的每一个返回的信号进行放大和数字化。通常,除了由换能器自身提供的自然带通滤波以外,信号在数字化之前的唯一滤波是用于A/D转换的低通抗混叠滤波。以与牵涉的频率一致的恒定频率对信号进行取样,并且以最少的处理将数字化的数据存储在内存中。信号采集的直接设计允许用现成组件将电路实现在相对少量的电路板面积内。
图9中示出了对插件模块的更详细查看。示出了多个采集通道,每一个由发送器、接收器、预放大器、A/D转换器和内存块组成。在接收期间,将换能器信号数字化并直接写往单独的内存块。内存块是双端口的,意味着在正从A/D转换器侧写入采集数据的同时,可以从计算机侧读取内存块。内存块显现为***CPU的正常扩展内存。应当注意,插件模块的容量不限于标准计算机扩展卡的正常容量,因为***优选地被装入定制的外壳中。同样,可以使用多个插件模块来适应大量换能器元件,每个模块处理换能器孔径的子集。
用于插件模块的组件(包括放大器、A/D转换器和关联的接口电路)和用于发送脉冲产生和信号采集的必需组件是容易商业获得的组件,本文中将不对它们进行详细描述。从接收的回声中获得的回声信号的RF数据存储所需的内存块实质上是与可商业获得的插件扩展内存卡中得到的电路相同的电路,外加了用于写入被数字化的信号数据的第二直接内存存取端口。(所接收的回声信号数据通常被称为RF数据,因为其包括由换能器产生的高频电振荡)。内存被映射到中央处理器的地址空间内并可以以与位于计算机主板上的其它CPU内存类似的方式被存取。内存的容量是使得对于多达256或更多个独立的发送/接收周期其能容纳单独的通道接收数据。由于身体中超声脉冲的往返行程的最大实际穿透深度约是500个波长,所以四倍于中心频率的典型取样频率将需要存储来自于单独换能器元件的多至4000个样本。对于16比特的取样精度和128个换能器通道,最大深度接收数据采集将对于每个发送/接收事件需要约一个兆字节(megabyte)的存储。为存储256个事件将因此需要256MB的存储,并且全部加起来,能够在几个插件卡上构造128个通道的***。
基于软件的超声***的另一方面是计算机主板和其关联组件。所提出设计的主板应当优选地支持多处理器CPU构造,以获得所需的处理功率。配备有电源、内存、硬盘存储器、DVD/CD-RW驱动器和监视器的完整多处理器计算机***是本领域技术人员公知的,其能够容易商业购买,因此将不更详细地进行描述。
基于软件的超声***必须真正地实现“高性能”,这意味着可与现有高端***相比的图像质量,以便对于医疗保健业提供显著益处。该层次的性能不能简单地通过将当前***的流通处理方法转换成软件实现来获得,因为在流通结构中简单添加对于一秒种的实时成像所需的全部处理操作给出了数量,该数量超出当前用几个通用处理器可获得的通常的每秒种操作数量。因此,需要新的处理方法,新的处理方法实现比流通方法高得多的效率。
在本发明的基于软件超声***结构的一个实施方式中,用于信号和图像处理的输入数据包括在一个或多个发送事件后从单独的换能器通道中获得的RF样本集。对于实施例,让我们考虑具有128个元件线性换能器阵列的典型2D成像扫描模式,如图10中所示。
在该情况中,“发送事件”将由来自多个换能器元件的定时脉冲组成,以产生多个声波,多个声波在介质中合并以形成从位于特定元件位置处的换能器上的原点向外发出的聚焦超声波束。多个发送事件(总共128个)产生横穿换能器面的宽度递增地顺序地发射的超声波束,因而查询了整个图像帧。对于这些发送波束中的每一个,从换能器中的128个接收器元件中的每一个中收集所接收的回声数据,并将所接收的回声数据组织成数据阵列,每一列代表由对应换能器元件接收的取样回声信号。因而,每个阵列具有对应于128个换能器元件的128列和在深度上对应于取得的样本数的多个行(在该情况下,我们将假设4096行导致4096个样本)。然后这些128个数据阵列构成足以产生一个完整图像帧的RF数据集。
值得注意的是,在流通结构中,上述的RF数据集甚至不存在(至少不是同时全部),因为当数据从换能器中流入时发生了波束和图像形成。换句话说,当发送事件后数据返回每个元件时,将它们处理和合并(被称为波束成形)以产生表示沿着单个波束(扫描线)的聚焦返回的单RF信号。将该RF信号处理(再次实时地)成回声振幅样本,回声振幅样本存储在存储阵列中。当已经处理了所有波束方向时,于是将回声振幅数据内插并格式编排成像素图像,以进行显示。由于所有处理是实时地发生的,所以处理电路必须能够“跟上”从换能器阵列流入的数据。
在本发明的基于软件的结构中,在处理之前存储所有输入数据。这将采集速率与处理速率脱开,若需要,这允许处理时间长于采集时间。这是高频扫描中的区别优势,在此情况下采集深度短,取样频率高。例如,10MHz扫描头可能具有约四厘米的可用成像深度。在此情况下,声音在组织中的速度指定了128个发送/接收事件中的每一个在52毫秒内采集并存储它们的数据,这是十分高的采集数据速率。在流通结构中,这些采集数据将在高处理速率下实时地被形成为扫描线。在本发明的基于软件的结构中,RF数据的存储允许处理占用长达显示的帧周期,显示的帧周期对于组织运动的实时可视化通常是33毫秒(30帧/秒)。对于128个像素列(与扫描线粗略类比),这将允许每列258毫秒的处理时间,而不是流通结构的52毫秒。对于典型扫描深度,该存储策略具有大大降低与流通结构相比的处理最大速率的效果。
输入数据的存储降低最大处理速率,但是不一定减少处理步骤的数量。为了实现这个目的,采取了超声数据处理的新方法。第一步骤是要认识到当***处于成像模式下时的最终目的是在输出显示器上产生图像。超声图像具有依赖于采集***的物理参数(如频率和阵列维数)的基本分辨率,并且可以被表示为具有对回声振幅或某一其它组织(声音)特性进行编码的像素值的矩形阵列。该矩形像素阵列的密度必须提供图像分辨率的足够空间取样。要认识到,显示图像不需要仅由像素的矩形阵列组成,而是可以由代表不同几何形状的任何任意的像素集组成。下一步骤是以该图像阵列中的像素之一开始并考虑RF数据集中的哪些样本点贡献该像素强度的计算,并且确定访问和处理它们的最有效方式。该方法是与当前流通结构使用的那个方法相比完全不同的方法,因为仅需要处理显示屏上贡献像素的信息。在本发明的方法中,显示图像上的小区域与大图像区域相比将花费更少的总处理时间,因为小区域包含更少的像素。相比之下,必须将流通处理方法设计成在不依赖图像区域大小的条件下处理最大数据流带宽。
在处理了恰当地表示超声图像所需的像素阵列后,可以将阵列以适当的大小提供给计算机显示屏,以供观看。计算机的图形处理器通常能够在不需要附加CPU处理的条件下执行由简单的缩放和内插组成的这个操作。
我们接下来考虑我们超声图像的单个像素的处理策略。在该讨论中,我们将假设,我们的目标是获得在像素关于换能器阵列的对应空间位置处的回声强度。可以类似地获得其它声音参数。我们的第一步是找到包含对回声强度计算做出贡献的样本的采集RF数据区。为了实现这个目的,对于图10的扫描方法,我们首先找到最接近于相交像素位置出现的采集扫描线,然后使用对应的单独元件数据阵列。图11对于超声图像中的示例像素示出了该映射过程。在图11中,所标示的像素映射到最接近的扫描采集线,该扫描采集线在该情况下是扫描线4,扫描线4的RF数据存在于第四单独元件RF数据阵列中(其表示从第四发送/接收事件中收集的数据)。多于一个RF数据阵列可以被选择为贡献像素信号,但是对于该实施例,我们将仅考虑单个数据阵列。
我们的下一步骤是标出单独元件阵列中的区域,该区域包含对像素的强度计算有贡献的样本。该映射过程相当复杂并依赖于几个因素。换能器元件每一个具有敏感区,该敏感区决定换能器元件将如何对从图像场中的特定点返回的信号做出响应。对于给定图像点,仅需要考虑具有预定阈值以上的敏感度的元件,因为如果敏感度过低,则元件将不对像素质量贡献有用信息。该敏感度阈值然后确定在所映射的区中包括的元件数据列的数量。如图5中所示,在换能器最右侧上的元件不包括在所映射的数据区内。
所映射的数据区的起始深度由每个单独换能器元件处的回声的到达时间决定。如图11中所示,后来及时捕获离图像点更远的元件的图像点信号,因此数据集在存储器中的起始点更深。最后,所映射的数据区需要的深度范围依赖于产生的发送脉冲的持续时间。较长的发送脉冲将对于较长的时间段激励图像点,产生遍布RF内存的更大深度跨度的回声信号。
幸运的是,可以对于给定像素方格对着手确定映射数据区的因素中的许多进行预计算,因为该方格关于实时图像序列的多个帧不变化。通过使用预计算的因子,能够快速且有效地确定给定像素的映射数据区,这在实时成像过程中节省了大量计算。
在选择出映射像素的RF数据之后,我们可以将其组织成矩阵RFPnm,如下所示:
记号“Pnm”指的是行n、列m中的图像像素。矩阵列是图11的竖线,其中假设每条竖线中样本的编号j是相同的。样本的编号j依赖于在捕获由发送脉冲产生的信号所需的时间内RF数据的范围。下标k是RF数据阵列中具有来自图像点的足够信号强度以参加强度计算的通道编号。
本说明书中引用的和/或申请数据页中列出的所有美国专利、美国专利申请公开、美国专利申请、外国专利、外国专利申请和非专利申请通过引用全部并入本文。
根据上述内容,将理解,虽然为了说明本文已经描述了公开的具体实施方式,但是可以在不背离本公开的精神和范围的情况下进行各种修改。例如,已经通过使用矩阵操作实现上述的产生像素或三维像素声音信息的处理操作,但是认识到可以使用标准数学运算或者甚至基于硬件的处理方法来实现一些或全部处理步骤。因此除了由权利要求限定的以外,本公开不是限制性的。
Claims (31)
1.一种产生高帧率定量超声多普勒图像的方法,包括:
在M个方向角下将多个多普勒系综发射到换能器的视野上方的介质内,所述多个多普勒系综中的每一个均包含N个未聚焦的声脉冲信号;
响应于每个声发射从所述换能器中的接收元件阵列中的多个接收元件中接收并存储散射的和反射的超声信号;
关于至少一个图像点处理所接收的超声信号以从与所述至少一个图像点对应的介质中的点中提取定量的运动信息;以及
在显示设备上产生从所述定量的运动信息中得到的多普勒图像,
其特征在于,根据预期流速选择如下两种方式中的一种对从所述多个多普勒系综的发射中获得的信号进行收集:
a)在一个角下收集每个系综的全部N个脉冲,并随后改变角直至完成在全部角下的收集;
b)对于M个角中的每一个,收集一个脉冲,并随后重复收集N个脉冲。
2.根据权利要求1所述的方法,其中,使用面向像素的处理方法实时地处理所接收和存储的散射和反射的超声信号的至少子集,所述面向像素的处理方法包括:
产生一组像素的初始步骤,该组像素被选择成表示接收所述声信号的所述换能器的视野中的范围,其中,该组像素中的每个给定像素对于所述多个接收元件具有已知的空间关系;
将给定图像像素映射到所存储的超声信号的区域内;
对于给定像素,将所存储的超声信号的映射区域组织成阵列;
通过在所述范围内的每个像素位置处使用矩阵处理,从所记录的回声信号的至少子集中重构所接收的散射超声信号的振幅和相位;
处理发送—接收事件的至少一个系综,其包括时间序列的声信号的至少一次发射和所散射和反射的超声信号的接收和存储,以得出所述范围中每个像素位置处的介质运动的估计多普勒参数;
产生在所述范围中每个图像像素位置处的一个或多个多普勒参数的图像,以用于产生显示。
3.根据权利要求2所述的方法,其中,所述参数涉及流体流或组织运动。
4.根据权利要求2所述的方法,其中,
多个未聚焦声发射信号以提供多个传播角到介质内的方式变化,并且处理根据所述多个传播角的多普勒参数估计,以计算每个图像像素位置处的速度向量。
5.根据权利要求2所述的方法,其中,
多个未聚焦声发射信号以提供多个传播角到介质内的方式变化,并且处理根据所述多个传播角的多普勒参数估计,以在介质中的真实速度与噪声之间辨别。
6.根据权利要求2所述的方法,其中,
多个未聚焦声发射信号以提供多个传播角到介质内的方式变化,并且处理根据所述多个传播角的多普勒参数估计,以计算每个图像像素位置处的多普勒能量。
7.根据权利要求2所述的方法,其中,处理包括使用多普勒频移频谱分析来得到每个图像像素位置处的介质运动的多普勒参数。
8.根据权利要求7所述的方法,其中,得到的介质运动的多普勒参数包括最大速度。
9.根据权利要求7所述的方法,其中,得到的介质运动的多普勒参数包括模、均值、方差、最大值和谱带中能量中的至少一个参数。
10.根据权利要求7所述的方法,其中,系综的处理包括:
在不降低多普勒帧的采集速率的情况下,实时地产生和显示至少一个图像像素位置的多普勒参数估计。
11.根据权利要求7所述的方法,其中,发射包括发射多个变化的未聚焦声信号,以提供到介质内的多个方向角的传播,并且处理包括使用来自于多个方向角的多普勒参数估计来计算每个图像像素位置处的速度向量。
12.根据权利要求7所述的方法,其中,发射包括发射多个变化的未聚焦声信号,以提供到介质内的多个方向角的传播,并且处理包括使用来自于多个方向角的多普勒参数估计来辨别由介质中真实运动产生的数据和由噪声产生的数据。
13.根据权利要求7所述的方法,其中,发射包括发射多个变化的未聚焦声信号,以提供到介质内的多个方向角的传播,并且处理包括使用来自于多个方向角的多普勒参数估计来计算频谱参数的绝对值,频谱参数包括每个图像像素位置处的频谱带中的模、均值、方差和最大能量。
14.根据权利要求2所述的方法,其中,产生图像包括产生图像或时间序列的图像,所述图像或时间序列的图像显示所述换能器阵列的整个视野上的得到的参数值。
15.根据权利要求2所述的方法,其中,产生图像包括产生所述图像的多个帧,所述图像的多个帧将介质的运动参数显示为影片序列的二维图像,其中,图像显示的速率能够与实时速率相同、慢于或者快于实时速率。
16.根据权利要求2所述的方法,其中,所述多普勒参数包括下面之中的至少一个:
对应于具有最大能量的多普勒移频的峰值速度;以及
作为湍流指示的最大谱宽度。
17.一种产生超声多普勒速度图像的***,包括:
模块,所述模块包括换能器,所述换能器将多个未聚焦声信号发射到所述换能器的整个视野的至少一部分上方的介质内并响应于每个声发射从所述换能器中的接收元件阵列中的多个接收元件中接收并存储散射和反射的超声信号;以及
处理器,所述处理器连接到所述模块,并被配置为关于至少一个图像点接收所储存的超声信号以从与所述至少一个图像点对应的介质中的点中提取定量运动信息并且在显示设备上产生从所述定量运动信息中得到的多普勒超声速度图像,
其特征在于,根据预期流速选择如下两种方式中的一种对多个未聚焦的超声信号进行收集:
a)在一个角下收集每个系综的全部N个脉冲,并随后改变角直至完成在全部角下的收集;
b)对于M个角中的每一个,收集一个脉冲,并随后重复收集N个脉冲。
18.根据权利要求17所述的***,其中,所述处理器被配置成通过使用面向像素的处理方法实时地处理所接收和存储的散射和反射的超声信号的至少子集,所述面向像素的处理方法包括:
产生一组像素的初始步骤,该组像素被选择成表示接收所述声信号的所述换能器的视野中的范围,其中,该组像素中的每个给定像素对于所述多个接收元件具有已知的空间关系;
将给定图像像素映射到所存储的超声信号的区域内;
对于给定像素,将所存储的超声信号的映射所述范围组织成阵列;
通过在所述范围内的每个像素位置处使用矩阵处理,从所接收的超声信号的至少子集中重构所接收的散射超声信号的振幅和相位;
处理发送—接收事件的至少一个系综,其包括时间序列的声信号的至少一次发射和响应于所述至少一次发射的所散射和反射的超声信号的接收和存储,以得出所述范围中每个像素位置处的介质运动的估计多普勒参数;
产生所述范围中每个图像像素位置处的一个或多个多普勒参数的图像,以用于产生显示。
19.根据权利要求18所述的***,其中,
多个未聚焦声发射信号以提供多个传播角到介质内的方式变化,并且处理根据所述多个传播角的多普勒参数估计,以计算每个图像像素位置处的速度向量。
20.根据权利要求18所述的***,其中,
多个未聚焦声发射信号以提供多个传播角到介质内的方式变化,并且处理根据所述多个传播角的多普勒参数估计,以在介质中的真实速度与噪声之间辨别。
21.根据权利要求18所述的***,其中,
多个未聚焦声发射信号以提供多个传播角到介质内的方式变化,并且处理根据所述多个传播角的多普勒参数估计,以计算每个图像像素位置处的多普勒能量。
22.根据权利要求18所述的***,其中,处理包括使用多普勒频移频谱分析来得到每个图像像素位置处的介质运动的多普勒参数。
23.根据权利要求22所述的***,其中,得到的介质运动的多普勒参数包括最大速度。
24.根据权利要求22所述的***,其中,得到的介质运动的多普勒参数包括模、均值、方差、最大值和谱带中能量中的至少一个参数。
25.根据权利要求22所述的***,其中,系综的处理包括:
在不降低多普勒帧的采集速率的情况下,在所述显示设备上实时地产生和显示至少一个图像像素位置的多普勒参数估计。
26.根据权利要求22所述的***,其中,发射包括从所述换能器中发射多个变化的未聚焦声信号,以提供从所述换能器到介质内的多个方向角的传播,并且处理包括在所述处理器中使用来自于多个方向角的多普勒参数估计来计算每个图像像素位置处的速度向量。
27.根据权利要求22所述的***,其中,发射包括发射多个变化的未聚焦声信号,以提供从所述换能器到介质内的多个方向角的传播,并且处理包括在所述处理器中使用来自于多个方向角的多普勒参数估计来辨别由介质中真实运动产生的数据和由噪声产生的数据。
28.根据权利要求22所述的***,其中,发射包括发射多个变化的未聚焦声信号,以提供从所述换能器到介质内的多个方向角的传播,并且处理包括在所述处理器中使用来自于多个方向角的多普勒参数估计来计算频谱参数的绝对值,频谱参数包括每个图像像素位置处的频谱带中的模、均值、方差和最大能量。
29.根据权利要求18所述的***,其中,产生图像包括在所述显示设备上产生图像或时间序列的图像,所述图像或时间序列的图像显示所述换能器阵列的整个视野上的得到的参数值。
30.根据权利要求18所述的***,其中,产生图像包括产生所述图像的多个帧,所述图像的多个帧将介质的运动参数在所述显示设备上显示为影片序列的二维图像,其中,图像显示的速率能够与实时速率相同、慢于或者快于实时速率。
31.根据权利要求18所述的***,其中,所述多普勒参数包括下面之中的至少一个:
对应于具有最大能量的多普勒移频的峰值速度;以及
作为湍流指示的最大谱宽度。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US7605708P | 2008-06-26 | 2008-06-26 | |
US61/076,057 | 2008-06-26 | ||
PCT/US2009/048440 WO2009158399A1 (en) | 2008-06-26 | 2009-06-24 | High frame rate quantitative doppler flow imaging using unfocused transmit beams |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102123668A CN102123668A (zh) | 2011-07-13 |
CN102123668B true CN102123668B (zh) | 2015-11-25 |
Family
ID=41082386
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980131855.5A Active CN102123668B (zh) | 2008-06-26 | 2009-06-24 | 使用未聚焦发送波束的高帧率定量多普勒流成像 |
Country Status (8)
Country | Link |
---|---|
US (1) | US10914826B2 (zh) |
EP (1) | EP2303131B1 (zh) |
JP (1) | JP5858783B2 (zh) |
KR (1) | KR101820422B1 (zh) |
CN (1) | CN102123668B (zh) |
CA (1) | CA2728998C (zh) |
DK (1) | DK2303131T3 (zh) |
WO (1) | WO2009158399A1 (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9220478B2 (en) | 2010-04-14 | 2015-12-29 | Maui Imaging, Inc. | Concave ultrasound transducers and 3D arrays |
US9265484B2 (en) | 2011-12-29 | 2016-02-23 | Maui Imaging, Inc. | M-mode ultrasound imaging of arbitrary paths |
US9282945B2 (en) | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
US9339256B2 (en) | 2007-10-01 | 2016-05-17 | Maui Imaging, Inc. | Determining material stiffness using multiple aperture ultrasound |
US9510806B2 (en) | 2013-03-13 | 2016-12-06 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
US9526475B2 (en) | 2006-09-14 | 2016-12-27 | Maui Imaging, Inc. | Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging |
US9572549B2 (en) | 2012-08-10 | 2017-02-21 | Maui Imaging, Inc. | Calibration of multiple aperture ultrasound probes |
Families Citing this family (75)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1229839A4 (en) | 1999-10-25 | 2005-12-07 | Therus Corp | USING FOCUSED ULTRASOUND FOR VASCULAR SEALING |
US6626855B1 (en) | 1999-11-26 | 2003-09-30 | Therus Corpoation | Controlled high efficiency lesion formation using high intensity ultrasound |
US20070213616A1 (en) | 2005-10-20 | 2007-09-13 | Thomas Anderson | Systems and methods for arteriotomy localization |
US8473239B2 (en) | 2009-04-14 | 2013-06-25 | Maui Imaging, Inc. | Multiple aperture ultrasound array alignment fixture |
KR101659910B1 (ko) | 2008-08-08 | 2016-09-27 | 마우이 이미징, 인코포레이티드 | 다중 구경 의료용 초음파를 통한 이미징 및 애드온 시스템의 동기화 |
US9513368B2 (en) * | 2010-06-30 | 2016-12-06 | General Electric Company | Method and system for ultrasound data processing |
WO2012051305A2 (en) | 2010-10-13 | 2012-04-19 | Mau Imaging, Inc. | Multiple aperture probe internal apparatus and cable assemblies |
US8903192B2 (en) * | 2010-10-14 | 2014-12-02 | Massachusetts Institute Of Technology | Noise reduction of imaging data |
JP5936857B2 (ja) * | 2011-01-27 | 2016-06-22 | 株式会社東芝 | 超音波診断装置及びその画像処理プログラム |
JP5905080B2 (ja) * | 2011-05-02 | 2016-04-20 | ヴェラゾニックス,インコーポレーテッド | オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成 |
RU2623301C2 (ru) | 2011-06-30 | 2017-06-23 | Конинклейке Филипс Н.В. | Способ и устройство для автоматизированной доплеровской оценки угла и скорости потока |
US9192359B2 (en) * | 2011-10-19 | 2015-11-24 | Verasonics, Inc. | Estimation and display for vector doppler imaging using plane wave transmissions |
TW201336478A (zh) * | 2011-12-01 | 2013-09-16 | Maui Imaging Inc | 使用以回音為基及多孔徑都卜勒超音波之移動偵測 |
US10231694B2 (en) | 2011-12-16 | 2019-03-19 | Koninklijke Philips N.V. | Automatic blood vessel identification by name |
WO2013088314A1 (en) * | 2011-12-16 | 2013-06-20 | Koninklijke Philips Electronics N.V. | Automated doppler pulse cycle selection |
KR101406806B1 (ko) * | 2011-12-27 | 2014-06-13 | 삼성메디슨 주식회사 | 초음파 영상을 제공하는 초음파 시스템 및 방법 |
KR101323330B1 (ko) * | 2011-12-28 | 2013-10-29 | 삼성메디슨 주식회사 | 결정 데이터에 기초하여 벡터 도플러 영상을 제공하는 초음파 시스템 및 방법 |
KR20130102913A (ko) * | 2012-03-08 | 2013-09-23 | 삼성메디슨 주식회사 | 조직의 이동 속도 및 방향 획득 방법 및 장치 |
JP6399999B2 (ja) | 2012-03-26 | 2018-10-03 | マウイ イマギング,インコーポレーテッド | 重み付け係数を適用することによって超音波画像の質を改善するためのシステム及び方法 |
KR101484958B1 (ko) * | 2012-03-28 | 2015-01-21 | 삼성메디슨 주식회사 | 초음파 시스템 및 초음파 영상 획득 방법 |
KR101433032B1 (ko) * | 2012-04-13 | 2014-08-21 | 서강대학교산학협력단 | 평면파를 이용한 기능성 혈류 영상 생성 방법 및 장치 |
CN102697525B (zh) * | 2012-05-04 | 2014-06-25 | 成都优途科技有限公司 | 全聚焦眼底彩超成像方法 |
CN102697524B (zh) * | 2012-05-04 | 2014-08-13 | 成都优途科技有限公司 | 全聚焦超声成像方法及其在血流成像中的运用 |
CN104582582B (zh) | 2012-08-21 | 2017-12-15 | 毛伊图像公司 | 超声成像***存储器架构 |
KR102070262B1 (ko) * | 2012-11-29 | 2020-03-02 | 삼성전자주식회사 | 초음파 프로브 장치 및 그의 제어 방법 |
US9025159B2 (en) | 2012-12-10 | 2015-05-05 | The Johns Hopkins University | Real-time 3D and 4D fourier domain doppler optical coherence tomography system |
DE112013007011T5 (de) * | 2013-04-30 | 2016-01-14 | Alpinion Medical Systems Co., Ltd. | Bildvergrösserungsverfahren und medizinische Ultraschallvorrichtung dafür |
CN105338909A (zh) * | 2013-04-30 | 2016-02-17 | 爱飞纽医疗机械贸易有限公司 | 影像放大方法和用于该影像放大方法的超声波医疗装置 |
EP3034004A4 (en) * | 2013-08-12 | 2017-05-17 | Samsung Electronics Co., Ltd. | Method for producing elastic image and ultrasonic diagnostic apparatus |
US9883848B2 (en) | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
CN106456118B (zh) * | 2013-11-19 | 2019-12-31 | 港大科桥有限公司 | 超声流体向量成像装置及其方法 |
WO2015180069A1 (zh) | 2014-05-28 | 2015-12-03 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像方法和*** |
JP5844430B1 (ja) * | 2014-06-27 | 2016-01-20 | 日立アロカメディカル株式会社 | 超音波診断装置 |
US10401493B2 (en) | 2014-08-18 | 2019-09-03 | Maui Imaging, Inc. | Network-based ultrasound imaging system |
US11026655B2 (en) * | 2014-09-26 | 2021-06-08 | Samsung Electronics Co., Ltd. | Ultrasound diagnostic apparatus and method of generating B-flow ultrasound image with single transmission and reception event |
KR102387708B1 (ko) * | 2015-01-30 | 2022-04-19 | 삼성메디슨 주식회사 | 향상된 hprf 도플러 영상을 위한 가이드를 제공하는 방법 및 초음파 시스템 |
US20160258904A1 (en) * | 2015-03-04 | 2016-09-08 | Bp Corporation North America Inc. | Process Control |
CA2981305A1 (en) * | 2015-04-01 | 2016-10-06 | Verasonics, Inc. | Method and system for coded excitation imaging by impulse response estimation and retrospective acquisition |
CN110013274B (zh) | 2015-04-29 | 2022-04-22 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流成像显示方法及超声成像*** |
CN104905815B (zh) * | 2015-05-22 | 2018-05-01 | 深圳先进技术研究院 | 一种线阵超声探头 |
EP3310261A4 (en) | 2015-06-19 | 2019-01-23 | Neural Analytics, Inc. | TRANSCRANIAL DOPPLER PROBE |
WO2017035838A1 (zh) * | 2015-09-06 | 2017-03-09 | 深圳迈瑞生物医疗电子股份有限公司 | 超声灰阶成像***及方法 |
US10588605B2 (en) * | 2015-10-27 | 2020-03-17 | General Electric Company | Methods and systems for segmenting a structure in medical images |
WO2017120388A1 (en) | 2016-01-05 | 2017-07-13 | Neural Analytics, Inc. | Systems and methods for determining clinical indications |
JP2019500155A (ja) | 2016-01-05 | 2019-01-10 | ニューラル アナリティクス、インコーポレイテッド | 一体型プローブ構造 |
US11589836B2 (en) | 2016-01-05 | 2023-02-28 | Novasignal Corp. | Systems and methods for detecting neurological conditions |
US10856846B2 (en) | 2016-01-27 | 2020-12-08 | Maui Imaging, Inc. | Ultrasound imaging with sparse array probes |
US11298110B2 (en) | 2016-02-26 | 2022-04-12 | Omar MANSOUR | Doppler measurement system and method |
CN107949331B (zh) * | 2016-06-30 | 2021-04-13 | 深圳迈瑞生物医疗电子股份有限公司 | 超声流体频谱多普勒成像方法和*** |
CA3032236C (en) | 2016-08-05 | 2024-05-28 | Norwegian University Of Science And Technology (Ntnu) | Ultrasound blood-flow monitoring |
US11717255B2 (en) | 2016-08-05 | 2023-08-08 | Cimon Medical As | Ultrasound blood-flow monitoring |
EP3512431A1 (en) * | 2016-09-16 | 2019-07-24 | Institut National de la Sante et de la Recherche Medicale (INSERM) | Method for imaging a sample with blood and associated devices |
US11246560B2 (en) | 2016-09-20 | 2022-02-15 | Samsung Medison Co., Ltd. | Ultrasound probe, ultrasound imaging apparatus, ultrasound imaging system, and method for controlling thereof |
KR102660559B1 (ko) | 2016-09-20 | 2024-04-26 | 삼성메디슨 주식회사 | 초음파 프로브, 초음파 영상 장치, 초음파 영상 시스템 및 그 제어방법 |
JP6745209B2 (ja) * | 2016-12-28 | 2020-08-26 | 株式会社日立製作所 | 超音波診断装置 |
JP7171583B2 (ja) * | 2017-01-16 | 2022-11-15 | コーニンクレッカ フィリップス エヌ ヴェ | 超高速撮像のためのシステム及び方法 |
EP4053594A3 (en) | 2017-07-06 | 2022-11-30 | Esaote S.p.A. | Ultrasound method and system for extracting signal components relating to spatial locations in a target region in the spatial and temporal frequency domain |
WO2019104468A1 (zh) * | 2017-11-28 | 2019-06-06 | 北京深迈瑞医疗电子技术研究院有限公司 | 一种造影成像方法以及超声成像设备 |
ES2951015T3 (es) * | 2018-02-07 | 2023-10-17 | Cimon Medical As | Monitorización del flujo sanguíneo por ultrasonidos |
WO2019155223A1 (en) * | 2018-02-07 | 2019-08-15 | Norwegian University Of Science And Technology (Ntnu) | Ultrasound blood-flow monitoring |
CN111787862B (zh) * | 2018-02-19 | 2023-09-29 | 国家医疗保健研究所 | 用于心脏的同时4d超快血流和组织多普勒成像并获得量化参数的方法和装置 |
US11771396B2 (en) * | 2018-03-01 | 2023-10-03 | Siemens Medical Solutions Usa, Inc. | Quantification of blood flow with ultrasound B-mode imaging |
US11506783B2 (en) * | 2018-03-09 | 2022-11-22 | Purdue Research Foundation | Method of processing an image |
EP3816670A4 (en) * | 2018-05-09 | 2022-02-09 | Furuno Electric Co., Ltd. | METEOROLOGICAL OBSERVATION DEVICE, METHOD AND PROGRAM |
EP3586759A1 (en) * | 2018-06-28 | 2020-01-01 | Koninklijke Philips N.V. | Methods and systems for performing color doppler ultrasound imaging |
GB201810711D0 (en) * | 2018-06-29 | 2018-08-15 | King S College London | Ultrasound Method and Apparatus |
US11647989B2 (en) * | 2018-09-11 | 2023-05-16 | Philips Image Guided Therapy Corporation | Devices, systems, and methods for multimodal ultrasound imaging |
KR102158177B1 (ko) | 2018-11-15 | 2020-09-22 | 서강대학교산학협력단 | 2D 샘플볼륨 기반의 스펙트럴 도플러 영상에서 Neural Network를 이용하여 움직임을 추적하는 기법 및 이를 이용한 초음파 의료 영상 장치 |
WO2020139775A1 (en) * | 2018-12-27 | 2020-07-02 | Exo Imaging, Inc. | Methods to maintain image quality in ultrasound imaging at reduced cost, size, and power |
CN110604591A (zh) * | 2018-12-29 | 2019-12-24 | 深圳迈瑞生物医疗电子股份有限公司 | 调整多普勒参数值的方法以及超声设备 |
CN109828029B (zh) * | 2019-03-28 | 2021-08-27 | 烟台中凯检测科技有限公司 | 一种基于原始数据的超声相控阵检测***和方法 |
CN111820946A (zh) * | 2020-06-23 | 2020-10-27 | 华中科技大学 | 一种用于多普勒超声探测的柔性测速装置及其应用 |
JP7472691B2 (ja) | 2020-07-09 | 2024-04-23 | コニカミノルタ株式会社 | 超音波診断装置、および、超音波信号処理方法 |
US20220199229A1 (en) * | 2020-12-03 | 2022-06-23 | Wavebase Inc. | Method and system for enhancing medical ultrasound imaging devices with computer vision, computer aided diagnostics, report generation and network communication in real-time and near real-time |
KR102636497B1 (ko) * | 2023-09-19 | 2024-02-14 | 주식회사 엣지케어 | 혈류 측정장치 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1650190A (zh) * | 2002-04-30 | 2005-08-03 | 皇家飞利浦电子股份有限公司 | 用于组织和血流成像的合成聚焦的超声诊断成像*** |
CN1981708A (zh) * | 2005-11-23 | 2007-06-20 | 美国西门子医疗解决公司 | 超声成像技术指导血栓治疗的造影剂增强超声治疗*** |
WO2008065570A1 (en) * | 2006-11-30 | 2008-06-05 | Koninklijke Philips Electronics, N.V. | Catheter with ultrasound transducer and variable focus lens used in aneurysm assessment |
CN101203183A (zh) * | 2005-04-14 | 2008-06-18 | 维拉声学公司 | 利用面向像素处理的超声成像*** |
Family Cites Families (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63117739A (ja) | 1986-11-06 | 1988-05-21 | アロカ株式会社 | 超音波ドプラ診断装置 |
JPH062134B2 (ja) * | 1989-09-08 | 1994-01-12 | 株式会社東芝 | 超音波診断装置 |
US5396890A (en) | 1993-09-30 | 1995-03-14 | Siemens Medical Systems, Inc. | Three-dimensional scan converter for ultrasound imaging |
US5365929A (en) * | 1993-10-04 | 1994-11-22 | Advanced Technology Laboratories, Inc. | Multiple sample volume spectral Doppler |
US5842473A (en) | 1993-11-29 | 1998-12-01 | Life Imaging Systems | Three-dimensional imaging system |
JP4056069B2 (ja) | 1994-12-27 | 2008-03-05 | 東芝医用システムエンジニアリング株式会社 | 超音波診断装置 |
US6248073B1 (en) | 1995-06-29 | 2001-06-19 | Teratech Corporation | Ultrasound scan conversion with spatial dithering |
JPH1033535A (ja) | 1996-07-30 | 1998-02-10 | Toshiba Corp | 超音波ドプラ診断装置および超音波ドプラ診断の方法 |
US5846200A (en) | 1996-11-08 | 1998-12-08 | Advanced Technology Laboratories, Inc. | Ultrasonic diagnostic imaging system for analysis of left ventricular function |
US6095980A (en) * | 1997-10-02 | 2000-08-01 | Sunnybrook Health Science Centre | Pulse inversion doppler ultrasonic diagnostic imaging |
US5910119A (en) * | 1998-05-12 | 1999-06-08 | Diasonics, Inc. | Ultrasonic color doppler velocity and direction imaging |
US6135956A (en) | 1998-10-01 | 2000-10-24 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging system with spatial compounding of resampled image data |
US6221020B1 (en) * | 1999-04-22 | 2001-04-24 | G.E. Medical Systems Global Technology Company, Llc | System and method for providing variable ultrasound analyses in a post-storage mode |
US6234968B1 (en) | 1999-06-15 | 2001-05-22 | Acuson Corporation | 3-D diagnostic medical ultrasound imaging using a 1-D array |
US6733455B2 (en) | 1999-08-20 | 2004-05-11 | Zonare Medical Systems, Inc. | System and method for adaptive clutter filtering in ultrasound color flow imaging |
US6773399B2 (en) | 2001-10-20 | 2004-08-10 | Zonare Medical Systems, Inc. | Block-switching in ultrasound imaging |
US6251073B1 (en) | 1999-08-20 | 2001-06-26 | Novasonics, Inc. | Miniaturized ultrasound apparatus and method |
US6685645B1 (en) | 2001-10-20 | 2004-02-03 | Zonare Medical Systems, Inc. | Broad-beam imaging |
US6468213B1 (en) | 1999-09-14 | 2002-10-22 | Ecton, Inc. | Medical diagnostic ultrasound system and method |
US6551246B1 (en) * | 2000-03-06 | 2003-04-22 | Acuson Corporation | Method and apparatus for forming medical ultrasound images |
US6390981B1 (en) * | 2000-05-23 | 2002-05-21 | Koninklijke Philips Electronics N.V. | Ultrasonic spatial compounding with curved array scanheads |
US6866631B2 (en) | 2001-05-31 | 2005-03-15 | Zonare Medical Systems, Inc. | System for phase inversion ultrasonic imaging |
US6618206B2 (en) | 2001-10-20 | 2003-09-09 | Zonare Medical Systems, Inc. | System and method for acoustic imaging at two focal lengths with a single lens |
US6663567B2 (en) | 2002-03-19 | 2003-12-16 | Zonare Medical Systems, Inc. | System and method for post-processing ultrasound color doppler imaging |
US6620103B1 (en) * | 2002-06-11 | 2003-09-16 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging system for low flow rate contrast agents |
US6866632B1 (en) | 2002-09-18 | 2005-03-15 | Zonare Medical Systems, Inc. | Adaptive receive aperture for ultrasound image reconstruction |
FR2848673B1 (fr) | 2002-12-17 | 2005-02-18 | Thales Ultrasonics Sas | Procede de conversion d'images pour presentation sous forme bitmap |
US6980419B2 (en) | 2003-03-12 | 2005-12-27 | Zonare Medical Systems, Inc. | Portable ultrasound unit and docking station |
JP4557573B2 (ja) * | 2003-03-14 | 2010-10-06 | 株式会社東芝 | 超音波診断装置及び超音波診断装置の作動方法 |
US20050049494A1 (en) | 2003-08-29 | 2005-03-03 | Arthur Gritzky | Method and apparatus for presenting multiple enhanced images |
US8469887B2 (en) * | 2003-12-19 | 2013-06-25 | General Electric Company | Method and apparatus for flow parameter imaging |
JP4660126B2 (ja) | 2004-06-18 | 2011-03-30 | 株式会社東芝 | 超音波血流イメージング装置 |
JP4625663B2 (ja) * | 2004-08-18 | 2011-02-02 | 株式会社東芝 | 超音波診断装置 |
JP4679141B2 (ja) | 2004-12-27 | 2011-04-27 | 株式会社東芝 | 超音波診断装置および超音波診断画像の表示方法 |
US8211024B2 (en) | 2005-06-06 | 2012-07-03 | Siemens Medical Solutions Usa, Inc. | Medical ultrasound pressure gradient measurement |
JP5342104B2 (ja) | 2006-10-30 | 2013-11-13 | 株式会社東芝 | 超音波診断装置及び画像処理プログラム |
US7984651B2 (en) | 2006-11-10 | 2011-07-26 | Penrith Corporation | Transducer array imaging system |
WO2008068709A1 (en) | 2006-12-04 | 2008-06-12 | Koninklijke Philips Electronics N.V. | Method and apparatus for multiline color flow and angio ultrasound imaging |
-
2009
- 2009-06-24 CN CN200980131855.5A patent/CN102123668B/zh active Active
- 2009-06-24 CA CA2728998A patent/CA2728998C/en active Active
- 2009-06-24 KR KR1020117002009A patent/KR101820422B1/ko active IP Right Grant
- 2009-06-24 EP EP09770932.3A patent/EP2303131B1/en active Active
- 2009-06-24 DK DK09770932T patent/DK2303131T3/en active
- 2009-06-24 US US12/490,780 patent/US10914826B2/en active Active
- 2009-06-24 WO PCT/US2009/048440 patent/WO2009158399A1/en active Application Filing
- 2009-06-24 JP JP2011516580A patent/JP5858783B2/ja active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1650190A (zh) * | 2002-04-30 | 2005-08-03 | 皇家飞利浦电子股份有限公司 | 用于组织和血流成像的合成聚焦的超声诊断成像*** |
CN101203183A (zh) * | 2005-04-14 | 2008-06-18 | 维拉声学公司 | 利用面向像素处理的超声成像*** |
CN1981708A (zh) * | 2005-11-23 | 2007-06-20 | 美国西门子医疗解决公司 | 超声成像技术指导血栓治疗的造影剂增强超声治疗*** |
WO2008065570A1 (en) * | 2006-11-30 | 2008-06-05 | Koninklijke Philips Electronics, N.V. | Catheter with ultrasound transducer and variable focus lens used in aneurysm assessment |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9526475B2 (en) | 2006-09-14 | 2016-12-27 | Maui Imaging, Inc. | Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging |
US9339256B2 (en) | 2007-10-01 | 2016-05-17 | Maui Imaging, Inc. | Determining material stiffness using multiple aperture ultrasound |
US9282945B2 (en) | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
US9220478B2 (en) | 2010-04-14 | 2015-12-29 | Maui Imaging, Inc. | Concave ultrasound transducers and 3D arrays |
US9247926B2 (en) | 2010-04-14 | 2016-02-02 | Maui Imaging, Inc. | Concave ultrasound transducers and 3D arrays |
US9265484B2 (en) | 2011-12-29 | 2016-02-23 | Maui Imaging, Inc. | M-mode ultrasound imaging of arbitrary paths |
US9572549B2 (en) | 2012-08-10 | 2017-02-21 | Maui Imaging, Inc. | Calibration of multiple aperture ultrasound probes |
US9510806B2 (en) | 2013-03-13 | 2016-12-06 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
Also Published As
Publication number | Publication date |
---|---|
EP2303131B1 (en) | 2015-04-22 |
EP2303131A1 (en) | 2011-04-06 |
CA2728998C (en) | 2021-11-16 |
JP2011526181A (ja) | 2011-10-06 |
KR101820422B1 (ko) | 2018-02-28 |
US10914826B2 (en) | 2021-02-09 |
CA2728998A1 (en) | 2009-12-30 |
US20090326379A1 (en) | 2009-12-31 |
DK2303131T3 (en) | 2015-05-04 |
JP5858783B2 (ja) | 2016-02-10 |
WO2009158399A1 (en) | 2009-12-30 |
CN102123668A (zh) | 2011-07-13 |
KR20110027806A (ko) | 2011-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102123668B (zh) | 使用未聚焦发送波束的高帧率定量多普勒流成像 | |
CN101203183B (zh) | 利用面向像素处理的超声成像*** | |
KR102055738B1 (ko) | 초음파 프로세싱 시스템에서 혈류 벡터 속도 이미지를 생성하는 방법 | |
US5840032A (en) | Method and apparatus for three-dimensional ultrasound imaging using transducer array having uniform elevation beamwidth | |
US9384530B2 (en) | Enhanced ultrasound image formation using qualified regions of overlapping transmit beams | |
Giangrossi | Development and real-time implementation of novel 2-D and 3-D imaging techniques on a research scanner |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |