CN102348415A - 用于获取和处理部分三维超声数据的系统和方法 - Google Patents

用于获取和处理部分三维超声数据的系统和方法 Download PDF

Info

Publication number
CN102348415A
CN102348415A CN2010800115310A CN201080011531A CN102348415A CN 102348415 A CN102348415 A CN 102348415A CN 2010800115310 A CN2010800115310 A CN 2010800115310A CN 201080011531 A CN201080011531 A CN 201080011531A CN 102348415 A CN102348415 A CN 102348415A
Authority
CN
China
Prior art keywords
data
ultrasound data
frame
partial
ultrasound
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.)
Pending
Application number
CN2010800115310A
Other languages
English (en)
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.)
Individual
Original Assignee
Individual
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
Priority claimed from US12/625,875 external-priority patent/US20100138191A1/en
Priority claimed from US12/625,885 external-priority patent/US20100185085A1/en
Application filed by Individual filed Critical Individual
Publication of CN102348415A publication Critical patent/CN102348415A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8993Three dimensional imaging systems
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/483Diagnostic techniques involving the acquisition of a 3D volume of data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52025Details of receivers for pulse systems
    • G01S7/52026Extracting wanted echo signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • G01S7/5209Details related to the ultrasound signal acquisition, e.g. scan sequences using multibeam transmission
    • G01S7/52093Details related to the ultrasound signal acquisition, e.g. scan sequences using multibeam transmission using coded signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52085Details related to the ultrasound signal acquisition, e.g. scan sequences
    • G01S7/52095Details related to the ultrasound signal acquisition, e.g. scan sequences using multiline receive beamforming
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8959Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52034Data rate converters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/5205Means for monitoring or calibrating

Landscapes

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

Abstract

用于获取和处理三维超声数据的方法包括获取部分三维超声数据。部分三维超声数据包括通过收集超声靶平面和收集至少一个超声偏置平面来收集的部分三维超声数据帧。该方法还包括处理部分三维超声数据。

Description

用于获取和处理部分三维超声数据的系统和方法
相关申请的交叉引用
本申请是以下在先申请号的继续申请:于2009年11月25日提交的且名称为“Method and System for Acquiring and Transforming Ultrasound Data”的美国专利序列号12/625,875以及于2009年11月25日提交的且名称为“Dynamic Ultrasound Processing Using Object Motion Calculation”的美国专利序列号12/625,885,二者都通过引用被全部并入。
本申请还要求于2009年1月19日提交的且名称为“DynamicUltrasound Acquisition and Processing Using Object Motion Calculation”的美国临时申请序列号61/145,710以及于2009年2月17提交的且名称为“System and Methodfor Tissue Motion Measurement Using 3D Ultrasound”的美国临时申请序列号61/153,250的利益,二者都通过引用被全部并入。
本申请涉及以下申请:(1)于2007年7月20日提交的且名称为“Methodof Tracking Speckle Displacement Between Two Images”的美国专利序列号11/781,212,(2)于2007年7月20日提交的且名称为“Method ofModifyingData Acquisition Parameters of an Ultrasound Device”的美国专利序列号11/781,217,(3)于2007年7月20日提交的且名称为“Method ofProcessingSpatial-Temporal Data Processing”的美国专利序列号11/781,223,以及(4)于2009年9月23日提交的且名称为“System and Method for Flexible RateProcessing of Ultrasound Data”的美国专利序列号12/565,662,这些申请都通过引用被全部并入。
有关联邦资助的研究或开发的声明
本发明由来自美国国家心肺和血液研究所(the National Heart,Lung,and Blood Institute)的资金(#5R44HL071379)资助,并且美国政府因此可以拥有本发明的某些权利。
技术领域
本发明大体上涉及医学超声领域,并且更具体地涉及用于在超声数据获取和处理领域中获取和处理三维超声的新的和有用的方法和系统。
附图的简要说明
图1是本发明的优选的实施方案的示意图;
图2A和2B是优选的实施方案的方法的变化形式的示意图;
图3是优选的实施方案的方法的变化形式的流程图,该方法包括使用编码传输信号的快速获取的变化形式;
图4和5是用于优选的快速获取方法的编码传输信号的图形表示;
图6是优选的实施方案的方法的变化形式的流程图,该方法包括使用局部子集获取的快速获取过程的变化形式;
图7是用于优选的快速获取方法的局部子集获取的图形表示;
图8是优选的实施方案的方法的变化形式的流程图,该方法包括帧选择的变化形式;
图9A和9B是帧选择的图形表示;
图10A和10B是包括多阶段散斑跟踪的优选方法的变化形式的流程图;
图11是用于距离估计的多阶段散斑跟踪的图形表示;
图12是优选的动态获取方法的示意图;
图13是优选的动态获取方法的详细示意图;
图14A和14B是优选的动态处理方法的示意图;
图15A-15C是优选的动态处理方法的变化形式的详细示意图;以及
图16是本发明的优选的实施方案的示意图。
优选实施方案的描述
以下的对本发明的优选实施方案的描述并不是用来将本发明限制于这些优选的实施方案,而更确切地是使本领域的任何技术人员能够实现和使用本发明。
1.用于获取和处理部分三维超声的方法
如图1中所示的,优选实施方案的用于获取和处理部分三维超声的方法包括获取部分三维超声数据(S110)(其优选地包括扫描靶平面的子步骤S112和扫描至少一个偏置平面的子步骤S114)以及处理与部分三维超声数据相关的超声数据(S190)。该方法起作用来获取数据的部分三维体积,由于部分三维数据的减小的体积尺寸,数据的部分三维体积实质上比普通三维数据更容易处理。此外,该方法优选地包括从所收集的超声数据计算目标运动(S150)。数据的部分三维体积优选地实现普通三维超声的三维运动跟踪益处,但是在二维平面中被测量。作为另一附加部分,优选的方法可以包括基于目标运动来修改系统参数(S170),如图2A中所示的。参数可以包括数据生成参数S171(即动态获取)和/或处理参数S181(即动态处理),如图2B中所示的。多种另外的可选形式可以被应用于该方法,例如多阶段散斑跟踪、使用编码传输信号的数据的快速获取、使用帧子集获取的数据的快速获取、帧选择、和/或可以用于部分三维数据的任何合适的过程,如图2B中所示的。优选实施方案的变化形式还可以按任何合适的顺序、组合或排列使用。
包括获取部分三维超声数据的步骤S110起作用来生成数据的部分三维体积。部分三维超声数据集优选地主要由部分三维超声数据帧(即图像)组成。三维超声数据帧优选地限定所扫描的体积。步骤S110优选地包括扫描靶平面的子步骤S112和扫描至少一个偏置平面的子步骤S114。优选地,与靶平面和偏置平面相关联的数据被组合以形成部分三维超声数据帧。此外,多个偏置平面可以被获取以形成更详细的三维数据。可选地,任何合适的方法可以被用于获取部分三维体积。时间部分三维超声数据优选地被获取以测量运动。两个或更多个部分三维数据帧优选地被用于测量帧之间的运动。步骤S110优选地包括收集数据和准备数据的子步骤。收集数据的步骤起作用来例如从超声换能器或存储原始超声数据的设备收集原始超声数据。原始超声数据可以由原始超声数据的实数的或复数的、经解调的或频移的(例如基带数据)或任何合适的表示来表示。准备数据起作用来进行预先处理,以将原始数据转换为合适的形式,例如亮度模式(B模式)、运动模式(M模式)、多普勒或超声数据的任何其他合适的形式。此外,准备数据优选地包括从靶平面和偏置平面的扫描形成部分三维超声帧。所获取的数据可以可选地作为原始超声数据留下,或所获取的数据可以可选地以准备好的数据格式从外部设备收集。此外,预波束形成的或后波束形成的数据可以被获取。所获取的数据优选地来自超声设备,但是可以可选地是对运动敏感的任何合适的数据获取系统。所获取的数据可以可选地由中间设备例如数据存储单元(例如硬盘驱动器)、数据缓冲器或任何合适的设备提供。所获取的部分三维超声还可以作为处理数据和控制数据被输出。处理数据优选地是将在步骤S190中被处理的数据。控制数据可以在步骤S150中用于运动计算并且用于系统参数修改。处理数据和控制数据优选地以相同的格式,但是可以可选地以上文描述的变化的形式。
包括扫描靶平面的子步骤S112起作用来获取所关注的材料(组织)的数据图像。靶平面的扫描优选地被超声换能器执行,但是可以使用任何合适的设备。数据图像优选地是沿着靶平面(超声波束被检查的平面)搜集的二维图像,或可选地,一维数据、三维数据或任何合适的数据可以被获取。
包括扫描偏置平面的子步骤S114起作用来获取平行于靶平面并且从靶平面偏置的材料的数据图像。偏置平面优选地实质上平行于靶平面并且被定位在靶平面的前方或后方,优选地分隔开预定的距离。偏置平面的扫描也以与靶平面实质上相似的方法被执行,但是可选地,可以使用不同的超声换能器、波束形状、平面的方位和/或图像类型。
包括计算目标运动的步骤S150起作用来分析所获取的数据以探测组织运动、探针运动和/或影响所获取的数据的任何其他运动。目标运动优选地包括影响所获取的数据的任何运动,例如组织运动、组织变形、探针运动和/或任何合适的运动。所测量的运动可以是组织速度、位移、加速度、应变、应变率或探针、组织运动或组织变形的任何合适的特征的测量。目标运动优选地使用原始超声数据来计算,但是可以可选地使用超声数据的任何合适的形式。在不同的时间获取的至少两个数据帧(例如数据图像或体积)优选地用于计算一维、二维或三维运动。散斑跟踪优选地被使用,但是可选地,可以使用多普勒处理、块匹配、交叉相关处理、横向波束调制(lateral beam modulation)和/或任何合适的方法。此外运动测量可以使用组织运动的模型来改进和改善。此外运动测量可以使用组织运动的模型来改进和改善。目标运动(或运动数据)优选地在步骤S170中的处理参数的修改中被用作参数输入,但是可以可选地或额外地被直接地用在处理步骤S190中。
散斑跟踪是通过跟踪作为从被扫描的目标的超声干涉和反射的结果的超声散斑的核心(部分)的位置来实现的运动跟踪方法。超声散斑的图样在小的运动上非常相似,这允许跟踪散斑核心随时间推移在搜索窗(或区域)内的运动。搜索窗优选地是核心被预期在其内找到的窗,假定有正常的组织运动。优选地,搜索窗还取决于超声数据的帧率。较小的搜索窗可以按较快的帧率使用,假定有相同的组织速度。核心的尺寸影响运动测量的分辨率。例如,较小的核心将导致较高的分辨率。来自散斑跟踪的运动可以使用各种算法例如绝对差和(SAD)或标准化交叉相关来计算。
包括处理部分三维超声数据的步骤S190起作用来转换所获取的数据以用于超声成像、分析或任何其他合适的目的。处理的步骤优选地帮助图像特征的探测、测量和/或可视化。在超声数据的处理完成之后,方法优选地继续进行以输出经处理的数据(即所转换的数据)(S198)。所输出的数据可以被用于任何合适的操作,例如被存储、显示、传递到另一个设备或任何合适的用途。处理的步骤可以是任何合适的处理任务,例如空间或时间滤波(例如用于多普勒成像和彩色血流成像的壁滤波)、求和、加权、排序、分类、再采样或其他过程,并且可以为了任何合适的应用而被设计。优选地,步骤S190使用在步骤S110中获取的部分三维超声数据并且还可以使用任何在步骤S170中修改的参数,如下文描述的。作为实例,目标运动数据(在步骤S150中被计算的)可以被用于自动地识别或区分开目标特征,例如血液和组织。根据情况,速度、应变或应变率计算或任何合适的计算可以在步骤S190中被优化以仅以所关注的目标特征为目标。例如,应变计算可以忽略与血液相关联的超声数据,作为提高组织变形测量的精确度的方式。处理数据可以是原始超声数据(例如RF数据)或其他合适的形式的数据,例如被转换为合适的形式(即预处理)的原始数据。处理优选地在数据正在被获取时实时地对超声数据执行,但是可以可选地离线地或远程地对所保存或缓冲的数据执行。如图14B中所示的,部分三维超声数据的处理优选地包括以下子步骤:形成超声图像(S192),再采样超声图像(S194),并且进行时间处理(S196)。S190的处理步骤可以优选地以任何合适的顺序进行,并且子步骤S192、S194和S196可以全部或部分地以任何合适的组合进行。
包括形成超声图像的步骤S192起作用来从在步骤S110中获取的部分三维超声数据输出超声图像。来自步骤S110的部分三维超声数据优选地被转换为用于处理操作的格式。本步骤是可选的,并且例如在处理步骤基于原始超声数据的情况下不是必需的。超声图像优选地是超声数据或来源于超声信号的数据——包括原始超声数据(即射频(RF)数据图像)、B模式图像(来自原始超声数据的幅值或包络探测的图像)、彩色多普勒图像、功率多普勒图像、组织运动图像(例如速度和位移)、组织变形图像(例如应变和应变率)或任何合适的图像——的任何空间表示。
包括超声图像的再采样的步骤S194起作用来将基于运动数据的处理参数应用于超声数据的处理。再采样优选地是在空间上聚焦的,时间处理在步骤S196中发生,但是步骤S194和步骤S196可以可选地在实质上同一个步骤中实现。超声图像改善可以使用运动数据作为用于图像处理操作的滤波器来作出。例如,运动数据可以被用于识别具有高组织速度的区域,并且应用图像校正(锐化或聚焦)来解释图像中的由运动导致的失真。另外地或可选地,超声图像的再采样可以包括使用帧之间的空间转换的测量在空间上映射数据,以将数据映射到公共栅格。在空间上映射数据优选地包括通过将图像帧适应性地转换至公共空间基准帧来移动和另外地扭曲图像。这优选地与步骤S196的时间处理配合地使用以实现运动补偿的帧平均。
包括进行时间处理的步骤S196起作用来应用连续的超声数据图像的基于时间的处理。时间处理优选地描述帧到帧(即时间序列)处理。此外,进行时间处理的步骤可以根据被目标运动计算控制的参数来执行。时间处理可以包括时间积分、加权求和(有限脉冲响应(FIR)滤波)以及帧群成员使用之前的时间处理输出的加权求和(无限脉冲响应(HR)滤波)。帧平均的简单的方法由具有对于每个帧的恒定的加权的FIR滤波器描述。帧平均或持续性可以被用于抑制噪声。帧平均一般在假定没有运动的情况下进行。时间处理可以另外地利用在步骤S194中进行的数据的空间映射,以增强帧平均。例如,使用以每秒20帧(即50ms帧内时间)获取数据的系统以及具有100ms的目标稳定时间(即下层目标可以被考虑为是恒定的时间)的目标,仅两个帧可以被平均或处理而没有图像质量降级。使用帧之间的空间转换的测量,数据可以在时间处理之前被映射至公共栅格,以补偿目标运动,提供较大的时间处理窗并且最终从信噪比增加提供改进的图像质量。在本实施例中,假定在探针和目标运动被除去时目标稳定时间增加10倍(至1秒)。现在,20帧可以被平均而没有降级,将信噪比提高了大于3的倍数(假定白噪声)。
2.使用数据的快速获取——编码传输信号的变异方法
如图3中所示的,优选实施方案的方法还可以用于数据的快速获取。数据的快速获取的技术可以通过多种变化形式实施。优选实施方案的编码传输信号变化形式包括以下的另外的步骤:多路复用与至少一个传输波束信号多路复用的第一传输波束信号(S122),传输被多路复用的传输波束信号(S124),接收相应于传输波束信号的至少一个接收波束(S126),以及将所接收的波束解复用成它们各自的信号(S128)。快速获取的方法优选地被应用于由上文描述的方法收集的部分三维数据,但是快速获取的方法可以可选地被应用于全三维数据或任何合适的数据。优选实施方案的该变化形式起作用来并行化获取以产生较快的帧率,但是可以可选地用于任何合适的目的。快速获取步骤优选地是步骤S110的子步骤并且在步骤S112和/或S114中使用。然而,快速获取步骤可以可选地被使用来代替扫描靶平面和扫描偏置平面以获取数据的部分三维体积。
包括多路复用与至少一个传输波束信号多路复用的第一传输波束信号的步骤S122起作用来多路复用所述传输波束。该步骤还可以优选地起作用来允许多个传输波束被同时传输。优选地,传输波束信号使用正交或实质上正交的代码调制。然而,传输波束信号可以使用任何合适的调制技术来多路复用。优选地,每个传输波束的脉冲被编码以唯一地识别它。
包括传输被多路复用的传输波束信号的步骤S124起作用来传输来自超声系统的作为传输信号的被多路复用的波束。被多路复用的传输波束信号优选地以相似于规则的所传输的波束的方式传输,但是可选地多个超声换能器可以每个传输被多路复用的传输波束信号的一部分,或信号可以用任何适当的方式传输。
包括接收相应于每个传输波束信号的至少一个接收波束的步骤S126起作用来探测在被多路复用的传输波束的所传输的超声脉冲传播时产生的超声回声。如图4中所示的,本发明的优选实施方案的这些技术通过同时收集多个区域中的信号增加了对于基于超声的组织跟踪的数据获取速率。在信号接收期间,所有接收波束优选被同时收集。可选地,接收波束可以被顺序地收集。
包括多路分离所接收的波束的步骤S128起作用来分离被多路复用的所接收的波束。来自多个接收波束的信号的处理优选地使用编码方案并行地进行。所接收的波束信号优选地使用对传输代码特定的滤波器被多路传输、解码、解调、滤波或“分类”为它们各自的信号。解码滤波器优选地仅作用于它们各自的信号,丢弃其他信号,如图5中所示的。为了提高图像质量,代码优选地是正交的或几乎正交的。
3.使用数据的快速获取——帧子集获取的变异方法
如图6中所示的,作为数据的快速获取的技术的另一种另外的或可选的变化形式,优选的方法包括以高速率收集全帧的局部子集(S132),在步骤S150中计算局部子集的目标运动,以及组合局部子集的目标运动信息(即跟踪结果)来以较低的速率形成全帧图像。该帧子集获取变化形式起作用来实现精确的组织(散斑)跟踪所必需的高帧率。如在图7中示例的,全帧的两个区域A和B被获取。波束组A和B被用于收集这些帧子集。波束的每个组以精确的组织跟踪所需要的速率被收集。图像的其他区域优选地以相似的方式被收集。这些技术有时被用于血液的彩色流成像,其也需要高局部帧率以测量高速血流。根据每个波束的获取时间(例如图像深度)、组中的波束的数量以及局部帧率、来自多个组的波束可以被顺序地收集。例如,收集方案可以是:从组1收集波束1、从组2收集波束1、从组1收集波束2、从组2收集波束2等。如上文指出的,帧子集获取的方法和编码传输信号的方法可以被组合。优选地,全帧的每个子集(部分)被获取,并且然后局部跟踪结果被组合来以较低的速率形成全帧图像。
4.使用帧选择的变异方法
如图8中所示的,优选实施方案的方法还可以与帧选择一起使用。帧选择的步骤优选地包括以下子步骤:在步骤S110期间以数据获取速率捕获超声数据,设定帧集间数据速率(S142),选择帧以形成多个帧集(S146),以及在步骤S190期间以所控制的数据速率处理来自存储器的数据。本发明的优选的方法还可以包括设定帧集内数据速率的步骤S144。帧选择的步骤起作用来允许高帧率数据速率(获取数据速率)根据第二数据速率(帧集间数据速率)被显示或处理。因为处理部分三维超声数据可以包括处理器密集的操作,所以帧选择优选地允许实时处理发生,同时维持高帧率数据,如图9A和9B中所示的。帧集优选地是以处理操作所必需的速率的帧的选择,并且帧集优选地根据帧集间数据速率被间隔开,使得显示或其他操作(具有不同的帧率要求)可以被充分执行。此外,处理优选地在原始超声数据或未被处理的超声数据上发生,但是可以可选地在被预处理的超声数据上发生。详细的分析、另外的处理、缓慢运动回放、快速运动回放和/或其他操作可以对超声数据执行,假定超声数据被存储在存储器中,同时仍然提供实时显示。虽然优选的方法集中于超声散斑跟踪,但是,在期望处理从获取速率或动态处理去耦合的情况下,它也可以被应用于其他超声成像模式。在一个实施例中,执行以每秒100帧需要数据并且以每秒30帧显示输出的处理任务,处理要求可以被减小至小于完全处理要求的三分之一而不牺牲结果的质量。
在步骤S110期间,部分三维超声数据优选地以高至足以实现散斑跟踪的速率捕获。数据获取速率优选地确定在所收集的超声帧之间的时间,如由图9B中的t1指示的。例如,与心脏扩张和收缩相关联的大形变率(即~2Hz的峰应变率)的精确的散斑跟踪需要优选地大于每秒100帧的帧率。该帧率是大于以每秒30帧的实时可视化所需要的帧率的约3倍。在大多数的情况下,精确的散斑跟踪所需要的帧率大于实时可视化速率所需要的帧率。大量组织运动的特征确定可视化速率,与超声与组织散射的相互作用相反,这确定散斑跟踪速率(也被称为帧集内速率)。数据获取速率可以根据技术限制或数据处理要求被设定为任何合适的速率。最大可视化速率被人类视觉感知限制为约每秒30帧。然而,更低的可视化速率可以是合适的,如由组织运动的细节(例如组织加速度)确定的。
包括设定帧集间数据速率的步骤S142起作用来根据预先规定的速率从所获取的数据选择(或采样)包括帧集的帧。帧集间数据速率被定义为所处理的帧集之间的时间,如由图9B中的t2指示的。在设定帧集间数据速率时,步骤S142优选地包括从所获取的部分三维超声数据选择帧以形成多个帧集(S146)。步骤S146起作用来形成用于处理的帧集。帧集优选地根据帧集间数据速率以及帧集的任何合适的参数被间隔开。帧集间数据速率优选地被设定为期望的输出数据速率,例如显示速率。帧集间数据速率小于或等于数据获取速率。帧集间数据速率优选地是数据获取速率的整数倍,但是以其他方式优选地独立于数据获取速率。获取速率设定帧集间采样的最大速率。另外地或可选地,帧集的参数可以根据处理步骤S190的需要或任何合适的需要来设定。参数优选地是帧集间数据速率,但是可以可选地包括帧集内数据速率、帧的数量、帧集的数量、帧或帧集的定时(例如非线性间隔)、触发事件(从其他生理事件)、数据压缩、数据质量和/或帧集的任何合适的参数。在一个变化形式中,帧集间数据速率在获取期间(例如S171的一部分)优选地根据生理运动被动态地调整,以更好地跟踪组织的相对运动(即在大组织运动和加速度的帧集之间的较短的时间,以及在小组织运动的帧集之间的较长的时间)。在图9B所示的实施例中,帧集速率(或输出产品速率)是获取速率的四分之一(1/4)。
作为步骤S190的一部分,部分三维超声数据以所控制的数据速率从存储器被处理。可选地或另外地,部分超声数据的以所控制的数据速率的处理可以在目标运动的计算(S150)期间例如对散斑跟踪发生。处理优选地对帧的帧集单独地执行。帧集优选地根据帧集间数据速率被顺序地处理。所控制的数据速率优选地被理解为包括任何所设定的数据速率,控制向处理器传送的数据速率,例如以帧集间数据速率处理帧集,以帧集内数据速率处理帧集的帧,以及可选地,以产品数据速率输出数据。散斑跟踪优选地对两个或更多个帧的帧集执行。散斑跟踪优选地至少以对运动测量或可视化足够的速率(例如每秒30帧集)来处理帧集,但是更高的或更低的帧率可以可选地被用于其他应用和要求。例如,机器视觉算法可能需要较高的可视化数据速率。较低的可视化数据速率可以用于长期监测或事件探测。可选地,可以进行任何合适的处理操作,例如插值法。处理操作优选地需要比最终的期望输出数据速率高的帧率。数据优选地在数据的处理之后以产品速率被输出。产品速率优选地等于帧集间数据速率,但是可以可选地不同于帧集间数据速率,取决于处理操作。
优选的方法还包括设定帧集内数据速率(S144),其起作用来调整帧集内的帧之间的时间,如由图9B中的t3指示的。帧集的帧之间的时间被获取速率限制。然而,虽然帧集优选地包括一对被顺序地获取的帧,但是帧集可以可选地包括以数据获取速率获取的一对被非顺序地获取的帧(即以数据获取速率获取的每隔一个帧)。获取速率设定帧集内采样的最大速率。然而,可以优选地根据生理运动使用可变的帧集内数据速率,以优化散斑跟踪性能(即具有迅速地变化的散斑的帧之间的较短的时间和具有缓慢地变化的散斑的帧之间的较长的时间)。可变的帧集内数据速率优选地在获取参数的修改(S171)期间被设定。帧集内采样数据速率优选地是数据获取速率的倍数,但是以其他方式优选地独立于数据获取速率。也在图9B所示的实施例中,帧集是一对被顺序地获取的帧,并且因此帧集的帧之间的时间是所获取的帧之间的时间,并且帧集内速率被确定为是数据获取速率。
5.使用多阶段散斑跟踪的变异方法
此外,优选的实施方案的方法可以用于多阶段散斑跟踪,如图10A和10B中所示的。在优选实施方案的多阶段散斑跟踪变化形式中,计算目标运动的步骤S150包括跟踪第一图像和第二图像之间的散斑位移。该变化形式的步骤S150优选地包括以下子步骤:计算至少一个第一阶段位移估计(S152)以及使用第一阶段位移估计来计算至少一个第二阶段位移(S154)。步骤S150和步骤S150的子步骤优选地被应用于在上文描述的方法中收集的部分三维数据,但是步骤S150和步骤S150的子步骤可以可选地被应用于全三维数据或任何合适的数据。多阶段散斑跟踪起作用来减少用于图像交叉相关的计算或其他合适的运动计算。如图10B中所示的,过程分辨率位移估计(course resolution displacement estimate)优选地被用作第一阶段位移估计,并且更精细的分辨率位移估计优选地被用作第二阶段位移估计。如图11中所示的,多阶段散斑跟踪的多分辨率变化允许来自低分辨率图像的距离估计以引导高分辨率位移估计。与没有初始低分辨率估计的单个精细位移估计比较,这优选地减少目标运动计算的计算。
包括计算至少一个第一阶段位移估计的步骤S152起作用来计算较低精确度的和/或较低分辨率的位移估计。优选地,第一阶段位移估计是来自超声成像的粗(低分辨率和/或精确度)位移估计。粗位移优选地通过交叉相关至少两个数据图像来计算,并且交叉相关函数的峰值优选地被用作粗位移估计。此外,数据图像的分辨率可以在估计过程之前被减小。然而,可以使用任何计算位移估计的方法,例如较不精确但是在计算上较廉价的位移算法。优选地,至少一个第一阶段位移估计被传送至步骤S154。至少一个第一阶段位移估计可以可选地被传送至相继的第一阶段估计阶段以进行另外的第一阶段位移估计。每个相继的阶段估计阶段优选地具有比之前的估计阶段更精确和/或更精细的分辨率的结果(例如对于过程位移估计的更精细的分辨率)。在过程分辨率估计的情况下,每个粗估计阶段可以将数据图像分辨率初始地减小至优选地比之前的阶段更精细的分辨率。作为另一种附加形式,过程位移估计可以被向上取样以匹配后面的估计阶段的分辨率。在将第一阶段估计传送至步骤S154之前可以可选地使用任何合适的数量的第一阶段估计。
包括使用第一阶段位移估计来计算至少一个第二位移的步骤S154起作用来使用第一阶段位移估计来计算更精确和/或更精细分辨率的位移。第一位移估计优选地被用作搜索偏置,以引导至少一个更精细的位移估计,与仅使用高精度和/或高分辨率阶段的处理相比提高了计算效率。来自步骤S152的第一阶段位移估计优选地确定交叉相关的原始图像的区域。优选地,第二阶段位移估计是使用步骤S152的粗分辨率位移估计的精细分辨率位移估计。精细分辨率位移优选地是交叉相关函数的峰值的位置。更优选地,精细分辨率位移处理优选地以整数像素精确度提供横向和轴向运动的估计。第二阶段位移可以可选地使用任何合适的方法来计算,例如使用第一阶段位移估计作为起始点来降低计算要求的更精确的(以及一般在计算上花费更多的)位移计算。
优选实施方案的变化形式的另外的子步骤包括计算子象素位移估计步骤S156,步骤S156起作用来进一步提高位移估计的精确度。优选地,对于子象素位移处理,仅需要相关函数的局部搜索区域。子象素位移计算优选地由拟合来自S154的相关函数以估计相关函数峰值的位置(即子象素滞后)的参量模型实现,或者如果复杂的图像帧被用作输入则通过交叉相关函数阶段的零交叉来实现。然而,子象素位移计算可以由任何合适的方法或设备实现。
6.使用动态获取的变异方法
如图12中所示的,优选实施方案的方法还可以被用于数据的动态获取,作为修改系统参数(S170)的可能的变化形式。优选实施方案的动态获取变化形式包括基于目标运动来修改数据生成的参数的步骤S171。该变化形式起作用来通过基于目标运动调整数据生成过程而实时地对改进的超声数据输出优化超声数据获取。所计算的目标运动被包括在数据获取系统的反馈环路中,以优化数据获取过程。
包括修改数据生成的参数的步骤S171起作用来改变用于处理的超声数据的收集和/或组织。修改数据生成的参数优选地改变数据获取的输入和/或输出。步骤S171可以包括多个子步骤。如图13中所示的,设备的收集超声数据的操作可以被改变,如在步骤S172中的,和/或所获取的数据可以在处理之前被改变,如在步骤S176和S178中的。
包括调整超声获取设备的操作的步骤S172起作用来基于目标运动数据调整超声获取设备的设置。超声数据获取设备的控制输入优选地根据使用目标运动计算的参数来改变。数据获取的可能被修改的参数优选地包括发送和接收波束位置、波束形状、超声脉冲波形,频率、点火速度和/或超声设备的任何合适的参数。此外,超声设备的修改可以包括修改靶平面的扫描和/或偏置平面的扫描。此外,偏置距离、偏置平面的数量或部分三维超声数据获取的任何合适的参数可以被修改。步骤S172可以另外地或可选地修改获取超声数据的任何变化形式——例如使用编码传输信号的快速数据获取、使用子集获取的快速数据获取、帧选择、多阶段获取和/或任何合适的变化形式——的参数。作为可能的修改的实例,之前的跟踪结果可以指示在图像的一部分中的图像或运动中的很少的运动或没有运动。帧率、局部帧率或获取速率可以被减小至较低的数据速率或与图像的其他区域权衡获取速率。作为另一个实例,波束间隔可以被自动地调整以匹配组织位移,潜在地提高数据质量(即测量的相关性)。
另外地或可选地,如图13中所示的,优选实施方案的方法可以包括修改数据形成的参数的步骤S176和形成数据的步骤S178。另外的步骤S176和S178起作用来将图像(数据)形成阶段从其他处理阶段去耦合。图像形成优选地定义超声数据的时间和空间采样。步骤S176和S178优选地作为步骤S171的一部分被执行,并且可以在修改或不修改超声获取设备的参数(S172)或有或没有方法100的任何其他可选的步骤的情况下被执行。
包括修改数据形成的参数的步骤S176起作用来使用所计算的目标运动来改变数据形成的参数。数据形成的参数优选地包括图像数据点的时间和/或空间采样、接收波束形成参数,例如孔径变迹和元件数据滤波,或数据形成过程的任何合适的方面。
包括形成数据的步骤S178起作用来组织用于超声处理的图像数据。基于目标运动的参数优选地在数据形成过程中使用。数据形成(或图像形成)阶段优选地定义从所获取的或所准备的超声数据生成的图像数据的时间和空间采样。所形成的数据优选地是超声图像。超声图像优选地是超声数据或来源于超声信号的数据——包括原始超声数据(即射频(RF)数据图像)、B模式图像(来自原始超声数据的幅值或包络探测的图像)、彩色多普勒图像、功率多普勒图像、组织运动图像(例如速度和位移)、组织变形图像(例如应变和应变率)或任何合适的图像——的任何空间表示。例如,使用孔径数据(即预波束形成的元件数据),样本可以沿着连续的波束形成以产生相似于传统的波束形成的数据。
7.使用动态处理的变异方法
此外,优选实施方案的方法可以与数据的动态处理一起使用,作为修改系统参数(S170)的可能的变化形式,如图14A中所示的。包括修改处理参数的步骤S181起作用来利用目标运动计算以增强或改进数据处理。滤波器或信号处理操作的系数或控制参数优选地根据与步骤S150中计算的目标运动相关的参数输入来调整。更优选地,所计算的目标运动被用作参数输入以修改处理参数。参数输入可以另外地或可选地包括其他信息,例如在下文更详细地讨论的数据质量度量。步骤S181可以包括变化形式,取决于数据处理应用。例如,数据处理可以包括使用散斑跟踪的组织运动计算。在这种情况下,窗优选地被增加尺寸,并且对于在静态组织的区域中的散斑跟踪的情况,搜索区域减小。相反地,数据窗优选地被减小尺寸,并且对于在移动的或变形的组织的区域中的散斑跟踪,搜索区域增大。运动控制的数据处理的另一个实例是图像帧配准(image frame registration)。在这种情况下,运动估计可以被用于再采样和排列(align)B模式或原始数据样本,用于改进的滤波、平均或任何合适的信号处理。图像再采样系数优选地被调整以提供帧配准。作为另一个实例,参数输入可以确定用于例如在再采样超声图像时处理超声数据的系数,或可选地,新的坐标系。所修改的处理参数可以另外地在以下应用中使用:各种算法——包括彩色血流(二维多普勒)、B模式、M模式和图像扫描转换——的空间和时间采样;用于彩色血流和多普勒处理的壁滤波;时间和空间滤波器编程(例如滤波器响应截止);散斑跟踪窗尺寸、搜索尺寸、空间和时间采样;设置散斑抑制算法的参数;和/或任何合适的应用。
作为另外的变化形式,如图15A、15B和15C中所示的,步骤S181可以与包括计算数据质量度量(DQM)(S160)的优选的实施方案的变化形式一起使用。步骤S160优选地起作用来通过确定反映数据的质量的值来帮助数据处理的优化。DQM优选地涉及数据是有效的保证的水平。数据质量度量优选地对于每个样本、图像区域的样本的子集和/或形成DQM地图的每个像素来计算。DQM优选地从与组织速度、位移、应变和/或应变率,或更具体地,峰值相关性、组织位移的时间和空间变化(例如导数和方差)、以及相关幅值的空间和时间变化有关的计算获取。
数据质量度量(DQM)优选地从步骤150的散斑跟踪方法的参数计算,并且更优选地是上文描述的DQI。使用标准化交叉相关执行的散斑跟踪产生被称为DQI的可以被用作DQM的量。标准化交叉相关优选地通过在目标的变形之前和之后获取超声射频(RF)图像或信号来执行。然后,图像的图像区域或窗使用交叉相关函数在两次获取之间被跟踪。交叉相关函数测量两个区域之间的相似性,作为区域之间的位移的函数。相关函数的峰值幅值相应于最大化信号匹配的位移。该峰值是DQI。DQI优选地在0.0至1.0的标度上表示,其中0.0表示低质量数据,而1.0表示高质量数据。然而,可以使用任何合适的标度。与组织相关联的数据的DQI趋于具有比含有血液或噪声的区域中的数据高的值。如下文描述的,这种信息可以在超声数据的处理中用于分割和信号识别。DQM优选地在步骤S181中被用作参数输入以修改处理参数。DQM可以被单独地使用以修改处理参数(图15A),DQM可以与所计算的目标运动配合地使用以修改处理参数(图15B),和/或DQM和运动信息可以被用于修改第一处理参数和第二处理参数(图15C)。
包括修改处理参数的步骤S181的变化形式优选地利用目标运动计算和/或DQM以增强或改进数据处理。滤波器或信号处理操作的系数或控制参数优选地根据与在步骤S150中测量的目标运动和/或步骤S160的DQM相关的参数输入来调整。处理参数的修改可以直接地基于DQM(图15A)和/或所计算的目标运动(图14A和14B)。处理参数的修改可以如在图15B中配合地或如在图15C中同时地(例如单独地但是并行)可选地基于处理参数的组合。
DQM的使用优选地使多种方式能够控制数据的处理。例如,测量例如B模式、速度、应变和应变率可以基于DQM来加权或分类(过滤)。DQM可以优选地被用于多种解释。DQM可以被解释为数据的质量的量化评定。具有不足够高的质量的数据可以从超声数据过滤。作为实例,对组织的一部分的超声导出的速度测量可能遭受噪声。在将速度测量过滤成仅包括具有高于0.9的DQI的测量之后,噪声水平被减少并且测量改进了。DQM可以可选地被解释为组织标识符。如上文提到的,DQI可以被用于区分开目标——特别是血液和组织——的类型。因此,DQI可以在处理超声数据时用于分割和信号或区域识别。作为一个应用的实例,DQM或更具体地DQI可以被用于确定血液与心脏壁边界,并且可以被用于自动地识别解剖结构或特征。处理操作还可以通过基于所识别的特征(例如组织或血液)选择性地执行处理任务来优化。例如,当计算组织的应变率时,具有血液的区域(如由低DQI指示的)可以在计算过程期间被忽略。此外,较高帧率和较高分辨率的图像要求更高的处理能力。根据组织类型使用DQM来分割超声数据或图像,组织特定的处理操作可以被用于减少对在计算上花费大的过程的处理要求。在本变化形式中,对所关注的数据执行在计算上花费大的过程。较不关注的数据可以接收不同的过程或较低分辨率的过程以减少计算成本。
8.用于获取和处理部分三维超声的系统
如图16中所示的,在超声系统中的优选的三维(3D)运动跟踪系统包括部分三维超声获取系统210、运动测量单元220和超声处理器240。系统起作用来获取数据的部分三维体积,数据的部分三维体积实质上由于与全体积三维数据相比的减小的体积尺寸而更容易处理。功能是在二维平面中产生三维运动测量。
部分三维超声获取系统210起作用来收集组织数据的部分三维体积。部分三维体积是具有比其他维度实质上小的尺寸和/或分辨率的一个维度的体积(例如三维体积的板或切片)。部分三维超声系统优选地包括超声换能器212以及数据获取设备214,超声换能器212扫描靶平面和至少一个偏置平面。优选地,从靶平面和偏置平面收集的数据每个是二维(2D)数据图像。靶平面和偏置平面优选地被组合以形成部分三维体积。在不同的时间获取至少两个体积使组织运动能够在三维中被测量。多个超声换能器可以被用于获取靶平面和偏置平面。可选地,超声数据的任何合适数量的平面、换能器的布置和/或波束形状可以被用于收集组织数据的部分三维体积。数据获取设备214优选地处理部分三维超声数据的数据组织。此外,部分三维超声获取系统210可以被设计为实现上文描述的过程,例如使用编码传输信号的快速获取、使用帧子集获取的快速数据获取、帧选择和/或超声获取的任何合适的过程。
优选的实施方案的超声换能器212起作用来从靶平面和偏置平面获取超声数据。超声换能器212优选地相似于通常用于一维或二维超声感测的超声设备,并且主超声换能器212优选地传输和探测超声波束。然而,超声换能器212可以是任何合适的设备。所传输的波束优选地使数据能够从材料(组织)收集,所传输的波束通过该材料传播。脉冲和波束的特征被波束形成器控制。靶平面优选地是二维数据图像并且优选地是由超声波束检查的区域。所获取的数据优选地是原始超声数据。原始超声数据可以具有多种表示,例如原始超声数据的实数的或复数的、经解调的或频移的(例如基带数据)或任何合适的形式。原始超声数据可以被准备来形成亮度模式(B模式)、运动模式(M模式)、多普勒或超声数据的任何合适的被准备的形式。
优选的实施方案的靶平面优选地是所关注的平面的二维超声数据。靶平面优选地由超声换能器扫描,但是可以可选地由专用设备、多个换能器或任何合适的设备获取。
优选的实施方案的偏置平面优选地与靶平面相同,除了下文指出的以外。偏置平面优选地平行于靶平面,但是偏置了任何合适的距离。距离优选地与目标运动(例如在偏置方向上的预期组织运动或探针运动)的期望幅值相同或相似。此外,任何合适的数量的偏置平面可以被获取。
优选的实施方案的数据获取设备214起作用来将超声数据组织为三维体积数据。数据获取设备214优选地处理下列任务:将数据传递到外部设备、存储数据、缓冲数据和/或任何合适的数据任务。数据获取设备优选地将数据保持为原始数据形式(未经处理的),但是数据获取可以可选地执行任何合适的预处理操作。
优选的实施方案的运动测量单元220起作用来分析数据的部分三维体积以探测目标运动。目标运动优选地包括组织运动、探针运动和/或影响所获取的数据的任何合适的运动。目标运动优选地使用原始超声数据来计算。在不同的时间获取的至少两组数据优选地被用于计算一维、二维或三维运动。散斑跟踪优选地被使用,但是可选地,可以使用多普勒处理、交叉相关处理、横向波束调制和/或任何合适的方法。运动测量还可以使用目标运动模型(例如参数拟合、空间滤波等)被改进和改善。运动测量单元220还可以计算数据质量度量(DQM),数据质量度量可以由超声数据处理器或系统的任何合适的部分用作输入变量。
此外,优选的实施方案的系统包括系统参数修改器230。系统参数修改器230优选地使用由运动测量单元生成的目标运动信息,用于调整整个系统的方面。更优选地,系统参数修改器修改部分三维超声获取系统的参数或超声数据处理器的参数。此外,运动测量单元的DQM可以被用于确定系统参数修改器的操作。
优选的实施方案的超声数据处理器240起作用来将超声数据变换为另一种形式的数据。超声数据处理器还可以使用由系统参数修改器确定的处理参数。
可选的实施方案优选地在存储计算机可读指令的计算机可读介质中实现上文的方法。指令优选地由用于获取和处理部分三维超声数据的计算机可执行部件执行。计算机可读介质可以被存储在任何合适的计算机可读介质例如RAM、ROM、闪存、EEPROM、光学设备(CD或DVD)、硬盘驱动器、软盘驱动器或任何合适的设备上。计算机可执行部件优选地是处理器,但是指令可以可选地或另外地由任何合适的专用硬件设备执行。如上文描述的超声获取设备还可以与计算机可执行部件配合地使用。
如本领域的技术人员将从之前的详细描述以及从附图和权利要求认识到的,可以对本发明的优选实施方案作出修改和变化而不偏离在下文的权利要求中限定的本发明的范围。

Claims (21)

1.一种用于获取和处理三维超声数据的方法,包括:
获取包括部分三维超声数据帧的部分三维超声数据,其中部分三维超声数据帧通过以下操作来收集:
收集超声靶平面;以及
收集至少一个超声偏置平面;以及
处理所述部分三维超声数据。
2.根据权利要求1所述的方法,其中所述偏置平面实质上平行于所述靶平面并且从所述靶平面位移一段设定的距离,并且所述靶平面和所述至少一个偏置平面配合地组合以形成三维体积图像。
3.根据权利要求1所述的方法,其中获取包括部分三维超声数据帧的部分三维超声数据还包括:
将第一传输波束信号与第二传输波束信号多路复用;
传输被多路复用的传输波束信号;
接收相应于所述第一传输波束信号的至少一个接收波束和相应于所述第二传输波束信号的至少一个接收波束;以及
解复用所接收的波束。
4.根据权利要求3所述的方法,还包括使用实质上的正交的代码来调制所述传输波束信号。
5.根据权利要求1所述的方法,其中获取包括部分三维超声数据帧的超声数据还包括:
以高速率收集全超声数据帧的局部子集;
从所收集的超声数据计算所述局部子集的目标运动;以及
组合所述局部子集的目标运动信息来以较低的速率形成全帧图像。
6.根据权利要求1所述的方法,还包括:
设定帧集间数据速率;
设定帧集内数据速率;
从所获取的部分三维超声数据选择帧来以帧集间数据速率和帧集内数据速率形成多个帧集;以及
其中处理超声数据是对所述帧集进行的。
7.根据权利要求1所述的方法,还包括从所获取的部分三维超声数据计算目标运动。
8.根据权利要求7所述的方法,其中计算目标运动还包括:
从第一超声数据帧和第二超声数据帧计算至少一个第一阶段位移估计;以及
从所述第一超声数据帧、所述第二超声数据帧和所述第一阶段位移估计计算至少一个第二阶段位移估计。
9.根据权利要求8所述的方法,其中第一阶段位移估计是比所述第二阶段位移估计更低分辨率的位移估计。
10.根据权利要求8所述的方法,其中第一阶段位移估计是比所述第二阶段位移估计更低精确度的估计。
11.根据权利要求7所述的方法,还包括基于所计算的目标运动来修改系统参数。
12.根据权利要求11所述的方法,其中修改系统参数包括基于目标运动来修改数据生成的参数。
13.根据权利要求12所述的方法,其中修改数据生成的参数包括调整正在获取所述部分三维超声数据的超声获取设备的操作。
14.根据权利要求12所述的方法,其中修改数据生成的参数包括修改数据形成的参数以及在处理所述超声数据之前形成所述部分三维超声数据。
15.根据权利要求11所述的方法,其中修改系统参数包括基于目标运动来修改处理参数。
16.根据权利要求15所述的方法,还包括计算数据质量度量;其中对处理参数的修改此外基于所述数据质量度量。
17.根据权利要求16所述的方法,其中处理部分三维超声数据包括形成超声图像,再采样所述超声图像以及进行时间处理。
18.根据权利要求12所述的方法,其中修改系统参数此外包括基于目标运动来修改处理参数。
19.根据权利要求18所述的方法,还包括:
通过执行快速获取数据的技术的一种技术来获取部分三维超声数据;
设定帧集间数据速率;
从所获取的部分三维超声数据选择帧来以帧集间数据速率形成多个帧集;
其中处理超声数据是对所述帧集进行的;以及
其中计算目标运动包括:
从第一超声数据帧和第二超声数据帧计算至少一个第一阶段位移估计;以及
从所述第一超声数据帧、所述第二超声数据帧和所述第一阶段位移估计计算至少一个第二阶段位移估计。
20.一种用于获取和处理三维超声数据的系统,包括:
部分三维获取系统,其收集靶数据平面和偏置数据平面,以形成部分三维超声数据帧;
运动测量单元;以及
超声处理器。
21.根据权利要求20所述的系统,还包括系统参数修改器,所述系统参数修改器使用所述运动测量单元的输出来调整所述系统的设置。
CN2010800115310A 2009-01-19 2010-01-15 用于获取和处理部分三维超声数据的系统和方法 Pending CN102348415A (zh)

Applications Claiming Priority (9)

Application Number Priority Date Filing Date Title
US14571009P 2009-01-19 2009-01-19
US61/145,710 2009-01-19
US15325009P 2009-02-17 2009-02-17
US61/153,250 2009-02-17
US12/625,885 2009-11-25
US12/625,875 US20100138191A1 (en) 2006-07-20 2009-11-25 Method and system for acquiring and transforming ultrasound data
US12/625,885 US20100185085A1 (en) 2009-01-19 2009-11-25 Dynamic ultrasound processing using object motion calculation
US12/625,875 2009-11-25
PCT/US2010/021279 WO2010083468A1 (en) 2009-01-19 2010-01-15 System and method for acquiring and processing partial 3d ultrasound data

Publications (1)

Publication Number Publication Date
CN102348415A true CN102348415A (zh) 2012-02-08

Family

ID=42340113

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010800115310A Pending CN102348415A (zh) 2009-01-19 2010-01-15 用于获取和处理部分三维超声数据的系统和方法

Country Status (4)

Country Link
US (1) US20100185093A1 (zh)
EP (1) EP2387360A4 (zh)
CN (1) CN102348415A (zh)
WO (1) WO2010083468A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103251430A (zh) * 2012-02-16 2013-08-21 美国西门子医疗解决公司 超声波切变波成像中的相关信息的可视化

Families Citing this family (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100138191A1 (en) * 2006-07-20 2010-06-03 James Hamilton Method and system for acquiring and transforming ultrasound data
US20080021319A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of modifying data acquisition parameters of an ultrasound device
US20080021945A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of processing spatial-temporal data processing
US20100185085A1 (en) * 2009-01-19 2010-07-22 James Hamilton Dynamic ultrasound processing using object motion calculation
US9275471B2 (en) 2007-07-20 2016-03-01 Ultrasound Medical Devices, Inc. Method for ultrasound motion tracking via synthetic speckle patterns
US20100086187A1 (en) * 2008-09-23 2010-04-08 James Hamilton System and method for flexible rate processing of ultrasound data
JP5665040B2 (ja) * 2009-09-10 2015-02-04 学校法人上智学院 変位計測方法及び装置、並びに、超音波診断装置
EP2623033B1 (en) * 2010-09-30 2017-01-11 Konica Minolta, Inc. Ultrasound diagnostic apparatus
JP5944913B2 (ja) * 2010-10-28 2016-07-05 ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. 超音波画像内の不均一回転歪みを低減するためのコンピュータ可読媒体及び同コンピュータ可読媒体を含むシステム
JP6109498B2 (ja) * 2011-07-05 2017-04-05 東芝メディカルシステムズ株式会社 超音波診断装置及び超音波診断装置制御プログラム
JP5848539B2 (ja) * 2011-07-26 2016-01-27 日立アロカメディカル株式会社 超音波データ処理装置
TWI446897B (zh) * 2011-08-19 2014-08-01 Ind Tech Res Inst 超音波影像對齊裝置及其方法
KR101894391B1 (ko) * 2011-10-05 2018-09-04 삼성전자주식회사 진단영상 생성장치, 의료영상시스템 및 빔포밍 수행방법
US9392995B2 (en) * 2012-07-25 2016-07-19 General Electric Company Ultrasound imaging system and method
EP2808760B1 (en) * 2013-05-29 2023-08-16 Dassault Systèmes Body posture tracking
US10034657B2 (en) 2013-07-26 2018-07-31 Siemens Medical Solutions Usa, Inc. Motion artifact suppression for three-dimensional parametric ultrasound imaging
WO2015058154A2 (en) * 2013-10-20 2015-04-23 Oahu Group, Llc Method and system for determining object motion
EP3282951A4 (en) * 2015-04-01 2019-01-16 Verasonics, Inc. METHOD AND SYSTEM FOR IMAGING WITH CODED AERATION BY IMPULSE RESPONSE ESTIMATION AND RETROSPECTIVE DETECTION
JP7105062B2 (ja) * 2017-12-21 2022-07-22 株式会社ソニー・インタラクティブエンタテインメント 画像処理装置、コンテンツ処理装置、コンテンツ処理システム、および画像処理方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6201900B1 (en) * 1996-02-29 2001-03-13 Acuson Corporation Multiple ultrasound image registration system, method and transducer
US20040006273A1 (en) * 2002-05-11 2004-01-08 Medison Co., Ltd. Three-dimensional ultrasound imaging method and apparatus using lateral distance correlation function
CN1989910A (zh) * 2005-12-27 2007-07-04 株式会社东芝 超声波图像处理装置以及超声波图像处理装置的控制程序
CN101101277A (zh) * 2007-08-10 2008-01-09 华南理工大学 一种高清晰焊缝超声成像无损检测方法及检测系统
US20080021945A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of processing spatial-temporal data processing
US20080019609A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of tracking speckle displacement between two images
US20080021319A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of modifying data acquisition parameters of an ultrasound device

Family Cites Families (63)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5675554A (en) * 1994-08-05 1997-10-07 Acuson Corporation Method and apparatus for transmit beamformer
US5685308A (en) * 1994-08-05 1997-11-11 Acuson Corporation Method and apparatus for receive beamformer system
US5503153A (en) * 1995-06-30 1996-04-02 Siemens Medical Systems, Inc. Noise suppression method utilizing motion compensation for ultrasound images
GB9518094D0 (en) * 1995-09-05 1995-11-08 Cardionics Ltd Heart monitoring apparatus
DE69721803T2 (de) * 1996-11-07 2003-11-20 Tomtec Imaging Systems Gmbh Verfahren und anlage zur rekonstruktion eines ultraschallbildes
US5919137A (en) * 1996-12-04 1999-07-06 Acuson Corporation Ultrasonic diagnostic imaging system with programmable acoustic signal processor
US5800356A (en) * 1997-05-29 1998-09-01 Advanced Technology Laboratories, Inc. Ultrasonic diagnostic imaging system with doppler assisted tracking of tissue motion
US5876342A (en) * 1997-06-30 1999-03-02 Siemens Medical Systems, Inc. System and method for 3-D ultrasound imaging and motion estimation
US5873830A (en) * 1997-08-22 1999-02-23 Acuson Corporation Ultrasound imaging system and method for improving resolution and operation
US6099471A (en) * 1997-10-07 2000-08-08 General Electric Company Method and apparatus for real-time calculation and display of strain in ultrasound imaging
US6074348A (en) * 1998-03-31 2000-06-13 General Electric Company Method and apparatus for enhanced flow imaging in B-mode ultrasound
US5934288A (en) * 1998-04-23 1999-08-10 General Electric Company Method and apparatus for displaying 3D ultrasound data using three modes of operation
US6066095A (en) * 1998-05-13 2000-05-23 Duke University Ultrasound methods, systems, and computer program products for determining movement of biological tissues
US6270459B1 (en) * 1998-05-26 2001-08-07 The Board Of Regents Of The University Of Texas System Method for estimating and imaging of transverse displacements, transverse strains and strain ratios
DE19824108A1 (de) * 1998-05-29 1999-12-02 Andreas Pesavento Ein System zur schnellen Berechnung von Dehnungsbildern aus hochfrequenten Ultraschall-Echosignalen
US6056691A (en) * 1998-06-24 2000-05-02 Ecton, Inc. System for collecting ultrasound imaging data at an adjustable collection image frame rate
US6162174A (en) * 1998-09-16 2000-12-19 Siemens Medical Systems, Inc. Method for compensating for object movement in ultrasound images
US6142946A (en) * 1998-11-20 2000-11-07 Atl Ultrasound, Inc. Ultrasonic diagnostic imaging system with cordless scanheads
US6213947B1 (en) * 1999-03-31 2001-04-10 Acuson Corporation Medical diagnostic ultrasonic imaging system using coded transmit pulses
US6352507B1 (en) * 1999-08-23 2002-03-05 G.E. Vingmed Ultrasound As Method and apparatus for providing real-time calculation and display of tissue deformation in ultrasound imaging
US6312381B1 (en) * 1999-09-14 2001-11-06 Acuson Corporation Medical diagnostic ultrasound system and method
US6512838B1 (en) * 1999-09-22 2003-01-28 Canesta, Inc. Methods for enhancing performance and data acquired from three-dimensional image systems
US6443894B1 (en) * 1999-09-29 2002-09-03 Acuson Corporation Medical diagnostic ultrasound system and method for mapping surface data for three dimensional imaging
US6210333B1 (en) * 1999-10-12 2001-04-03 Acuson Corporation Medical diagnostic ultrasound system and method for automated triggered intervals
US6282963B1 (en) * 1999-10-12 2001-09-04 General Electric Company Numerical optimization of ultrasound beam path
US6350238B1 (en) * 1999-11-02 2002-02-26 Ge Medical Systems Global Technology Company, Llc Real-time display of ultrasound in slow motion
US6447450B1 (en) * 1999-11-02 2002-09-10 Ge Medical Systems Global Technology Company, Llc ECG gated ultrasonic image compounding
US6277075B1 (en) * 1999-11-26 2001-08-21 Ge Medical Systems Global Technology Company, Llc Method and apparatus for visualization of motion in ultrasound flow imaging using continuous data acquisition
US6527717B1 (en) * 2000-03-10 2003-03-04 Acuson Corporation Tissue motion analysis medical diagnostic ultrasound system and method
US6346079B1 (en) * 2000-05-25 2002-02-12 General Electric Company Method and apparatus for adaptive frame-rate adjustment in ultrasound imaging system
US6318179B1 (en) * 2000-06-20 2001-11-20 Ge Medical Systems Global Technology Company, Llc Ultrasound based quantitative motion measurement using speckle size estimation
ATE369794T1 (de) * 2000-11-15 2007-09-15 Aloka Co Ltd Utraschalldiagnosegerät
US6447454B1 (en) * 2000-12-07 2002-09-10 Koninklijke Philips Electronics N.V. Acquisition, analysis and display of ultrasonic diagnostic cardiac images
US6447453B1 (en) * 2000-12-07 2002-09-10 Koninklijke Philips Electronics N.V. Analysis of cardiac performance using ultrasonic diagnostic images
US6537221B2 (en) * 2000-12-07 2003-03-25 Koninklijke Philips Electronics, N.V. Strain rate analysis in ultrasonic diagnostic images
US6666823B2 (en) * 2001-04-04 2003-12-23 Siemens Medical Solutions Usa, Inc. Beam combination method and system
US6605042B2 (en) * 2001-08-10 2003-08-12 Ge Medical Systems Global Technology Company, Llc Method and apparatus for rotation registration of extended field of view ultrasound images
US6537217B1 (en) * 2001-08-24 2003-03-25 Ge Medical Systems Global Technology Company, Llc Method and apparatus for improved spatial and temporal resolution in ultrasound imaging
US6638221B2 (en) * 2001-09-21 2003-10-28 Kabushiki Kaisha Toshiba Ultrasound diagnostic apparatus, and image processing method
US6676603B2 (en) * 2001-11-09 2004-01-13 Kretztechnik Ag Method and apparatus for beam compounding
US6776759B2 (en) * 2002-02-27 2004-08-17 Ge Medical Systems Global Technology Company, Llc Method and apparatus for high strain rate rejection filtering
US7314446B2 (en) * 2002-07-22 2008-01-01 Ep Medsystems, Inc. Method and apparatus for time gating of medical images
US6994673B2 (en) * 2003-01-16 2006-02-07 Ge Ultrasound Israel, Ltd Method and apparatus for quantitative myocardial assessment
US7558402B2 (en) * 2003-03-07 2009-07-07 Siemens Medical Solutions Usa, Inc. System and method for tracking a global shape of an object in motion
US7131947B2 (en) * 2003-05-08 2006-11-07 Koninklijke Philips Electronics N.V. Volumetric ultrasonic image segment acquisition with ECG display
US7033320B2 (en) * 2003-08-05 2006-04-25 Siemens Medical Solutions Usa, Inc. Extended volume ultrasound data acquisition
US7536043B2 (en) * 2003-08-18 2009-05-19 Siemens Medical Solutions Usa, Inc. Flow representation method and system for medical imaging
US20050096543A1 (en) * 2003-11-03 2005-05-05 Jackson John I. Motion tracking for medical imaging
US20070276236A1 (en) * 2003-12-16 2007-11-29 Koninklijke Philips Electronics N.V. Ultrasonic diagnostic imaging system with automatic control of penetration, resolution and frame rate
US7088850B2 (en) * 2004-04-15 2006-08-08 Edda Technology, Inc. Spatial-temporal lesion detection, segmentation, and diagnostic information extraction system and method
US20050288589A1 (en) * 2004-06-25 2005-12-29 Siemens Medical Solutions Usa, Inc. Surface model parametric ultrasound imaging
US7366278B2 (en) * 2004-06-30 2008-04-29 Accuray, Inc. DRR generation using a non-linear attenuation model
US7983456B2 (en) * 2005-09-23 2011-07-19 Siemens Medical Solutions Usa, Inc. Speckle adaptive medical image processing
US8191359B2 (en) * 2006-04-13 2012-06-05 The Regents Of The University Of California Motion estimation using hidden markov model processing in MRI and other applications
US20100138191A1 (en) * 2006-07-20 2010-06-03 James Hamilton Method and system for acquiring and transforming ultrasound data
JP5148094B2 (ja) * 2006-09-27 2013-02-20 株式会社東芝 超音波診断装置、医用画像処理装置及びプログラム
US20080125657A1 (en) * 2006-09-27 2008-05-29 Chomas James E Automated contrast agent augmented ultrasound therapy for thrombus treatment
US20080114251A1 (en) * 2006-11-10 2008-05-15 Penrith Corporation Transducer array imaging system
US20080214934A1 (en) * 2007-03-02 2008-09-04 Siemens Medical Solutions Usa, Inc. Inter-frame processing for contrast agent enhanced medical diagnostic ultrasound imaging
US20100185085A1 (en) * 2009-01-19 2010-07-22 James Hamilton Dynamic ultrasound processing using object motion calculation
US20100086187A1 (en) * 2008-09-23 2010-04-08 James Hamilton System and method for flexible rate processing of ultrasound data
US9275471B2 (en) * 2007-07-20 2016-03-01 Ultrasound Medical Devices, Inc. Method for ultrasound motion tracking via synthetic speckle patterns
US20100081937A1 (en) * 2008-09-23 2010-04-01 James Hamilton System and method for processing a real-time ultrasound signal within a time window

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6201900B1 (en) * 1996-02-29 2001-03-13 Acuson Corporation Multiple ultrasound image registration system, method and transducer
US20040006273A1 (en) * 2002-05-11 2004-01-08 Medison Co., Ltd. Three-dimensional ultrasound imaging method and apparatus using lateral distance correlation function
CN1989910A (zh) * 2005-12-27 2007-07-04 株式会社东芝 超声波图像处理装置以及超声波图像处理装置的控制程序
US20080021945A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of processing spatial-temporal data processing
US20080019609A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of tracking speckle displacement between two images
US20080021319A1 (en) * 2006-07-20 2008-01-24 James Hamilton Method of modifying data acquisition parameters of an ultrasound device
CN101101277A (zh) * 2007-08-10 2008-01-09 华南理工大学 一种高清晰焊缝超声成像无损检测方法及检测系统

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103251430A (zh) * 2012-02-16 2013-08-21 美国西门子医疗解决公司 超声波切变波成像中的相关信息的可视化

Also Published As

Publication number Publication date
US20100185093A1 (en) 2010-07-22
WO2010083468A1 (en) 2010-07-22
EP2387360A4 (en) 2014-02-26
EP2387360A1 (en) 2011-11-23

Similar Documents

Publication Publication Date Title
CN102348415A (zh) 用于获取和处理部分三维超声数据的系统和方法
CN102348416A (zh) 使用目标运动计算的动态超声处理
US7925068B2 (en) Method and apparatus for forming a guide image for an ultrasound image scanner
US9275471B2 (en) Method for ultrasound motion tracking via synthetic speckle patterns
CN111432733A (zh) 用于确定超声探头的运动的设备和方法
US20080306384A1 (en) Apparatus and method for computing 3D ultrasound elasticity images
US20140364735A1 (en) Ultrasound waveform tomography with tv regularization
US6494837B2 (en) System and method for three-dimensional ultrasound imaging using a steerable probe
JP2006091011A (ja) マルチビュー視差表示
JPH08506656A (ja) 多次元信号処理及び表示
US10453193B2 (en) Methods and system for shading a two-dimensional ultrasound image
CN113543718B (zh) 用于确定超声探头的包括前后方向性的运动的装置和方法
US9332958B2 (en) Ultrasonic diagnostic apparatus and method of displaying ultrasonic image
CN107072639B (zh) 使用切变波进行测量的超声系统和方法
US9247922B2 (en) Ultrasonic diagnostic apparatus and ultrasonic image display method
CN101258525A (zh) 用于对心脏右心室进行可靠的3d评价的超声系统及其方法
CN104268412A (zh) 一种角道集射线层析偏移速度分析方法及装置
US8675450B2 (en) Displacement estimating method and displacement estimating apparatus
KR20150118734A (ko) 초음파 영상 장치 및 그 제어방법
CN103906473A (zh) 超声波成像装置、超声波成像方法以及超声波成像用程序
CN110811674A (zh) 超声波诊断装置、程序以及超声波诊断装置的动作方法
Almekkawy et al. Two-dimensional speckle tracking using parabolic polynomial expansion with riesz transform
CN103747740A (zh) 超声波数据处理设备
Tang et al. Three dimensional height information reconstruction based on mobile active sonar detection
CN103717145B (zh) 超声图像处理装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20120208