CN101442937B - 具有运动补偿的超声合成发射聚焦 - Google Patents
具有运动补偿的超声合成发射聚焦 Download PDFInfo
- Publication number
- CN101442937B CN101442937B CN2007800167563A CN200780016756A CN101442937B CN 101442937 B CN101442937 B CN 101442937B CN 2007800167563 A CN2007800167563 A CN 2007800167563A CN 200780016756 A CN200780016756 A CN 200780016756A CN 101442937 B CN101442937 B CN 101442937B
- Authority
- CN
- China
- Prior art keywords
- thread
- launching
- line
- motion
- image field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
- A61B8/5276—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts due to motion
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/12—Diagnosis using ultrasonic, sonic or infrasonic waves in body cavities or body tracts, e.g. by using catheters
-
- 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
- A61B8/14—Echo-tomography
- A61B8/145—Echo-tomography characterised by scanning multiple planes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
-
- 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
-
- 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
- G01S7/52095—Details related to the ultrasound signal acquisition, e.g. scan sequences using multiline receive beamforming
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
一种超声诊断成像***通过为多线接收发射沿阵列间隔的多个波束而生成具有扩展聚焦范围的图像。多个发射波束的接收多线空间对准并与各自接收多线间的相位调整合并以防止非预期相位对消。合并的多线生成扩展发射焦点效应,以便使用合并多线生成的图像呈现出扩展聚焦范围。为防止运动伪影,可根据图像运动调整多线阶数。
Description
本发明涉及医学诊断超声***,具体而言涉及使用多线接收波束形成器扩展聚焦区域并具有降低运动影响的超声***。
合成聚焦原理一直是大量调查研究的主题,应用该原理将超声图像聚焦于像场中的每个点。例如,美国专利4,604,697(Luthra等人)描述了一种合成聚焦技术,在该技术中,从超声换能器阵列的每个元件中序列发射超声脉冲。由该阵列的所有元件接收从每次发射所接收到的回波信号并将其进行存储。在接收完全部信号之后,根据每个换能器元件相对于图像中每个点的位置以及超声信号到达每个点并从其返回的飞行时间的知识,可以在像场中每个点处形成聚焦回波信号。将合适的接收信号进行合并以形成图像中每个点的相干回波信号。选择用于形成每个点的所存储数据值为图像中每个点提供波束形成。虽然该方法能生成位于像场中每个点处的聚焦信号,但是其仍有若干缺陷。一个缺陷是必须存储来自整个像场的r.f.信号以进行处理。这需要进行大量的信息存储。第二个缺陷是需要进行大量处理以对图像中每个点的数据进行选择和加权,然后将合适的加权数据进行合并以计算图像数据点。第三个缺陷是由于单换能器元件发射的能量有限,因此该方法只对浅的穿透深度有效。
采用合成聚焦基本原理的一个具体应用是常规延迟及求和接收波束形成器,在该波束形成器中施加于来自每个接收元件的信号上的延迟等同于合成聚焦技术中的数据选择。常规波束形成器是这些原理的有限应用,因其发射聚焦于特定聚焦区域的发射波束并且动态聚焦仅沿着这一单发射波束的回波。从而需要多次发射以扫描整个像场。由此产生的功效是不需要存储所有发射到图像中每个点的数据;立即处理从一次发射接收到的数据以形成沿波束方向的相干回波信号。其局限性在于每个接收到的波束在发射时只聚焦于所选的聚焦区域。但是由于启动多个换能器元件发射波束,从而能够获得适当穿透,因此提高了较深处的信噪比。
美国专利6,231,511(Bae等人)以及Bae等人后续发表的题为“A Studyof Synthetic-Aperture Imaging with Virtual Source Elements in B-ModeUltrasound Imaging Systems”的文章(IEEE Trans.UFFC47卷2000年第6期1510页及之后各页)中提出合并标准聚焦波束形成器与合成聚焦的相关方面以提高常规聚焦区域之外的横向分辨率,从而实现在像场中所有点处的发射聚焦效果。该方法是基于标准发射波束焦点处的“虚源元件”假说提出的,该标准发射波束从该“虚源”向内和向外辐射能量。在标准发射聚焦波束被发射之后,由接收孔径的换能器元件接收能量并将其进行存储。在扫描完整个像场之后,根据由包绕着所述场中每个点的每个虚源场的各元件所接收到的信号计算该点处的回波信号。由于虚源模型是关于所发射焦点的沙漏形(hourglass-shaped)场,而根据接收到的多条扫描线的信号计算出在从焦点起的深度上除去的各点,因此将仅由一个波束成像出焦点处的像点。结果为示出了从发射焦点向内向外各点处横向分辨率得到改善的图像。但是如上所描述的基本合成孔径方法,必须存储大量数据用于处理来自每个接收孔径中每个元件的r.f.信号。另外,生成的图像被认为在焦点附近较暗,这是因为只有一个发射和接收对该像点及其分辨率有贡献,而多个发射和接收对除该发射焦点以外的各点有贡献。因此希望在图像中至少很大一部分上实现发射聚焦而不必存储大量r.f.数据。
根据本发明的原理,描述了一种诊断超声***和方法,其实现在有意义的场深上的发射聚焦而不必存储r.f.信号数据。发射超声波束,其声穿透多个扫描线位置的至少一部分,并且沿着多条扫描线(例如,平行的或通过时分复用的)并行处理接收波束。一系列此类发射生成与合用扫描线位置相关的多条扫描线。合并相关的扫描线数据以生成在有意义的场深上有效发射聚焦的回波数据。被合并的来自多个发射时间间隔的扫描线的数量是检测到的运动的函数。在下文描述的示例中,超声***用多线波束形成器接收同时发生的波束。在一些临床应用中,本发明的实现能够通过减少常要进行多区聚焦的需求来提高存在运动时超声成像的帧频。
在附图中:
图1a-1c示出了本发明的三波束示例的波束轮廓;
图2a-2d示出了本发明的六波束示例的波束轮廓;
图3a-3d示出了本发明的四波束示例,其中为图示清楚偏移前四个波束方向图(beam pattern);
图4a-4d示出了图3中四波束示例的延续,从而示出了接收波束的对准;
图5以方框图的形式示出了根据本发明原理构建的超声***;
图6a和6b示出了本发明的四波束示例的两个不同波束的延迟和加权特性;
图7以方框图的形式示出了一种超声***,该***适应性响应运动来调整本发明实现中合并的多线数量。
首先参照图1a-1c,重叠的波束轮廓示出了发射三个发射波束之后在每种情况下接收来自每个发射波束的三个波束。图1a示出了低于波束中心处强度峰值的稳定水平的发射波束轮廓10,所述波束由发射该波束的换能器阵列8发射并从该换能器阵列8中延伸出来。由设计者选择发射波束轮廓水平,其可以为3dB、6dB、20dB或低于波束中心处最大强度的其它水平。可以看到波束轮廓通过常规发射聚焦而在波束轮廓宽度最窄处的焦点12周围进行。换能器阵列8下面示出了波束20的正交视图,可以看到该波束包括中心瓣20a和位于主瓣20a两侧的旁瓣。发射波束在聚焦区域12处达到其最密集聚焦并由此发散。在其他实现中可以使用发散的发射波束。也可实现在相当深度处的聚焦。
发射具有包绕多条接收线14、16和18的宽度的发射波束10、20。通常从越小的发射孔径进行发射能产生越宽的波束。即:启动阵列8中较少数量的元件而该非阵列中全部数量的元件以发射波束。随后沿三个接收线位置14、16和18接收和聚焦发射回波。如下面所讨论的,以三种不同方式对接收孔径的换能器元件所接收到的回波进行延迟并求和,以便响应一个发射波束而形成位于不同线位置14、16和18处的多条线。在该示例中,接收线14沿发射波束10、20的中心进行接收,横向导引并聚焦接收线14和18以在中心线两侧进行接收。在该示例中,仅有外侧线14和18的近场和远场部分位于发射波束轮廓10之内。在这些区域,外侧线14和18接收到来自中心线位置两侧的发射能量,由此对中心线位置两侧像场中的目标进行取样,进而有效使用近场和远场中发射波束的横向扩展能量进行图像接收和分辨率。
在图1b中,已经通过将发射孔径向右移动一条接收线的间隔来发射第二波束。第二发射波束具有与第一发射波束相同的波束轮廓,并由波束轮廓曲线10′勾勒出。如在第一波束的情况下,响应第二发射在接收线位置16′、18′和22处同时接收三条接收线并使其形成波束。结果,接收线16′与来自第一发射的接收线16对准,接收线18′与来自第一发射的接收线18对准,而接收线22位于第二发射中心线18′的右侧。像第一组接收线一样,存储第二组接收多线16′、18′和22用于接下来的处理。
在图1c中,已经从再次向右移动一条接收线的中心孔径位置发射出第三波束。该发射波束由波束轮廓10″勾勒出,并在该发射之后同时接收三条接收线18″、22′和24。像前一次的接收线一样,这三条接收线全部或部分位于其发射波束的波束轮廓中,且具有与上述波束各线相同的间隔。结果,接收线18″与第二发射的接收线18′和第一发射的接收线18轴线对准,而接收线22′与第二发射的接收线22轴线对准。现在,位于接收线18、18′和18″路径中的目标由来自各自不同发射波束的三条接收线进行取样。按照下文所述的方式合并这些共同对准的波束以生成沿该线的图像数据线,所述波束与任意单条线的情况相比聚焦在更深的场深,从而产生扩展发射聚焦效应。由于合并来自三次波束发射的回波能量产生生成的图像数据,因此将在更深的场深上进行有效地聚焦。
发射和接收以这种方式持续穿过像场,直到扫描完整个像场。每当采集完给定线位置的最大接收线数量(本示例中为三)时,共同处理各接收线以生成该位置处的图像数据线。由于接收到的r.f.信号在接收时形成多条线的波束,从而不必存储来自任何发射的求和前r.f.数据,仅有限地需要存储在线位置处的之前各行直到采集完该位置处的所有接收线,这时可对它们进行处理并释放线存储器以用于存储后面各线。
图2a-2d给出了本发明的另一示例,其中发射波束轮廓30包括六个接收线位置处的全部或部分接收线,在图2a中以31-36标示。不同于第一个示例,在该示例中,发射波束中心没有接收线。取而代之,中心接收线33和34在波束中心的两侧间隔有接收线间距的一半。外侧接收线32和35在近场和远场处位于波束轮廓30之内,最外侧接收线31和36只有近场部分位于发射波束轮廓之内。可根据下面说明确定在图像数据形成中不使用近场行部分31和36。
如图2b中第二发射波束的轮廓30′所示,下一发射波束发射到第一发射波束右侧一个接收线距。在接收线位置32′、33′、34′、35′、36′和37处再次同时接收六条接收线并形成波束。可以看到这些接收线的前五个与第一发射波束的接收线32、33、34、35和36对准,从而在这些线位置中的每处提供用于处理的第二接收线。图2c示出了发射第三发射波束和接收该波束轮廓30″内的六条接收线之后的结果。现在有在线位置33″、34″、35″和36″处接收到的全部或部分三条接收线,在位置37′处接收到的全部或部分两条接收线,而在位置38处的部分一条接收线。如图2d所示的波束轮廓30″′,随着第四发射,有在线位置34″′、35″′和36″′处接收到的全部或部分四条接收线,在线位置37′处接收到的全部或部分三条接收线,在线位置38′处接收到的全部或部分两条接收线,而在位置39处有接收到的部分一条接收线。可以看到,当以这种方式持续扫描时,除了在末端处将如本示例中所示地接收到较少条线以外,在跨越孔径大部分的线位置处会接收到全部或部分六条接收线。当合并接收线数据形成图像数据时,每个接收线位置处的采样次数越多将产生越有效的发射聚焦。
可以使用数量更多的同时接收线,例如:八条、十二条或十六条间隔开的同时接收线,这样可以在发射时使用较少的数量F来声穿透较宽范围的接收线位置。
图3和图4示出了本发明的另一个示例,其使用来自每个发射波束的四条接收线。在这些附图中,逐次的波束和接收线组并不重叠,而是为了图示清楚而竖直对准。每个发射波束40-1、40-2、40-3等由指向下的虚线箭头示出,并且将来自每个发射波束的接收线在对应发射波束的两侧以实线箭头示出。图3a-3d示出了前四个发射-接收序列,且发射波束40-n向右移动用于每个相继发射时间间隔的一个接收线间隔。在该发射-接收序列的最后可以看到接收到的四条接收线都与接收线44对准,这四条接收线为:来自第二波束的接收线44-1、来自第三波束的接收线44-2以及来自第四波束的接收线44-3。合并这四条接收线以生成位于像场中这些接收线处的图像数据线。
图4a-4d示出了前述发射-接收序列的继续,其具有每个发射波束四条同时接收线的四个多发射-接收时间间隔。在这些附图中,接下来四个发射波束以40-5、40-6、40-7和40-8示出。如这些附图所示,现在有在线位置44、45、46、47和48处接收到的四条接收线。在这些位置中每个位置处接收到四条接收线之后,能够合并这四条线形成一个图像数据线,并删除所存储的接收线,以便可在相同存储位置存储随后各线。每当接收到包括四条对准线的组的另外第四线时,能够合并该组的四条线以在该位置处形成图像数据线,并形成供随后各线使用的存储器。序列以这种方式继续,接收每个发射波束跨像场的四条接收线,从而可在除像场横向末端之外全部区域上的每个线位置处合并四条接收线。
图5以方框图的形式示出了根据本发明原理构建的超声成像***。超声探头102包括由各换能器元件构成的换能器阵列104。由发射波束形成器106在相应的延迟时间启动所选择的换能器元件组,以便沿着该阵列从(一个或多个)期望起点开始发射波束,该波束在期望方向上聚焦于所选择聚焦区域。发射波束形成器通过发射/接收开关耦合到换能器元件,该发射/接收开关可包括保护接收器输入使之免于所施加的高电压发射脉冲的交叉点开关。向多线处理器110a-110n的输入施加由阵列104中每个换能器元件响应于每个发射波束所接收到的回波。每个多线处理器包括施加其自身延迟集合的接收波束形成器,以及如果必要,也包括变迹加权器,其用于对从各阵列元件接收到的回波进行加权以形成来自相同发射波束的经不同导引的接收波束。例如,在美国专利6,695,783(Henderson等人)和美国专利5,318,033(Savord)中可获得用于多线处理器110a-110n的合适的多线波束形成器。多线处理器110a-110n的输出耦合到线存储器112,该线存储器112至少存储接收到的多线直到已经采集到形成显示数据线所需的所有多线。向对应的各个乘法器116a-116n施加形成特定显示数据线所使用的多线组,以生成相应线位置的显示数据。如果必要,可以由变迹加权器114a-114n对来自每条线的回波数据进行加权。总之,这些加权器会根据每条线的往返脉冲响应对该条线进行加权。能够通过设振幅项(x,z)为像场中位置(x,z)处点通过发射波阵面的声透射振幅而获得合适的加权算法,方位角位置x=0对应于发射波束的中心轴。设X为与接收到的多线相对于发射波束轴线的方位角。为形成深度Z处的图像点而向该接收到的多线施加的权重为:
权重(X,z)=振幅(X,Z)
为了确定适当的延迟特性,设传播_时间(x,z)为发射波阵面到达位置(x,z)处的点所需要的传播时间,方位角x=0再次对应于发射波束的中心轴。设X为接收线相对于发射波束轴线的方位角。为形成深度Z处的图像点而向该接收到的多线所施加的延迟为:
延迟(X,Z)=传播_时间(X,Z)-传播_时间(0,Z)
其中,传播_时间(0,Z)为到达位于轴线外同样深度处的点的时间。
例如,函数振幅(X,Z)和传播_时间(X,Z)可通过模拟发射场获得。计算传播时间的合适方式是使用来自以若干频率的单频模拟的场相位延迟。通过对以若干频率的场振幅求平均可以计算出振幅。另外,可以向各加权器施加依赖深度的标准化。其用公因子乘以给定深度处的所有加权器。例如,可以对标准化进行选择以便斑点区域具有深度相关的均匀亮度。通过根据深度的函数改变各加权器,可以随深度动态改变孔径的大小和形状(变迹)。
不需要从模拟本***中所使用的确切发射特性获得振幅和传播时间。例如,设计者可选择使用不同的孔径大小或不同的变迹。
来自每条线的各回波由乘法器116a-116n进行加权,并由延迟线118a-118n进行延迟。总之,如上所示,这些延迟与发射波束中心到接收线位置的位置相关。使用延迟来补偿存在于具有不同发射-接收波束位置组合的多线间存在的相移变化,从而合并信号的相差不会引起信号对消。
应当领会到在数字***中,可通过在存储器中存储经加权的多线回波数据以及随后读出数据(这实现了必要的延迟)来实现延迟线。也可以使用不同长度的移位存储器和时钟信号来实现数字延迟,或者可以使用如前面提到的美国专利6,695,783中描述的内插波束形成器。延迟的信号由加法器120进行合并,且合成信号被耦合到图像处理器122。该图像处理器可执行扫描转换或其他处理来改进所显示的图像。在图像显示器124上显示生成的图像。
在图5的***中,延迟线118和加法器120实现了对从若干接收多线接收的各信号进行重聚焦,所述接收多线在给定方向上共同对准。重聚焦操作为从使用每个多线的不同发射波束位置引起的相差进行调整,从而防止合并信号中的非预期相位对消。加权器114与发射波束和多线位置的接近度相关地加权多线的贡献,从而将更高的权重给予具有更高信噪比的接收波束。这产生沿每条接收线的扩展的场深以及由于在每条接收线方向上多个取样的合并引起的穿透性提高(信噪比提高)。
在图4a-4d的示例中可以看到,每条接收到的多线只与其他共同对准的多线合并使用一次。这意味着本发明的实现不需要第二组加权器114、乘法器116和延迟器118,因为这些加权和延迟效应可由多线处理器110的加权器和延迟器实现,这使得该多线聚焦于合适的接收方向并解决了与相对于要进行合并的多线的相位和发射波束变化,所有这些都在一个过程步骤中执行。在接收到给定线位置的所有多线之后,将它们进行简单求和以提供显示线数据。或者,将给定线位置的每条接收到的多线存储在该线的线累加器中。该线位置上每条后续接收到的多线都添加到累加器的内容中直到在累加器中已经合并了该位置处全部补集的多线。累加器的内容然后向前发送给图像处理器,释放累加器以供另一线位置使用。
图6a和6b示出了可在合并给定显示线位置的各接收多线中使用的加权和延迟特性曲线的示例。图6a示出了接收到的多线的示例性加权和延迟特性曲线,该多线距发射波束的中心(例如图2b中接收到的多线33′)较远。示出的接收到的多线33′沿z方向从位于z=0处的换能器阵列表面延伸到位于图右侧的最大扫描深度。用虚线描绘出多线的中心,其中多线33′的该部分位于所选波束轮廓30′之外,而下面是其响应,设计者认为其可以接受。于是,加权特性曲线82用这样的权重对该多线加权,即该权重在波束轮廓之外时最小,而在多线将要用于重聚焦时为非零水平。在另外的实现中,加权特性曲线82在近场可下降为零。这是因为在施加必要延迟之后回波可能太过靠近阵列或在阵列之后,从而不容许精确接收。因此,在极近场中不会在多线合并中使用来自于横向远离多线的多线。
可以看到相位特性曲线84在发射波束焦点处跨越零相位调整,从而在近场中保持几乎恒定,而在远场中逐渐降低。
图6b示出了接收到的多线33的示例性加权和延迟特性曲线,该多线距其发射波束中心较近。由于更加靠近,如加权特性曲线86所示,在合并中该多线被施以较大的总体加权。增加聚焦区域附近的加权以平衡该区域中较远处多线(例如33′)加权的下降。由于横向较远多线的延迟特性曲线84是根据较靠近中心的多线33计算得到的,因此相位调整特性曲线88被描绘为直线。尽管在给定的实现中,设计者会选择使用紧邻极近场中线位置的发射波束的单线回波,但是在该示例中,使用多线33进行整个图像深度上的重聚焦。
本发明的实现可与各种接收函数一起使用。例如,实现可对来自经聚焦的子孔径的信号起作用。在另一个变更中,代替本身形成接收波束,可以使用有限数量的正交函数,例如傅立叶分量。然后可以在傅立叶空间合并不同的发射。合并的信号并不精确对应于接收到的图像线。由此通过合并傅立叶分量或子孔径形成接收波束。变迹改变时由波束得到的函数和若干波束的多种组合同样落入本发明的范围。在其它的实现中,可通过内插额外的中间多线来增加少量接收到的多线,然后应用增加数量的多线执行本发明的处理过程。
图7示出了图5的超声***的变更,其中关注运动的影响。图5所示的示例合并了从具有不同发射-接收时间间隔的数条共同对准多线所接收的r.f.信号。由于所合并的是r.f.信号而非检测到的信号包络,因此在接收待合并的多线期间图像中没有运动是十分重要的。无论像场中的物质的显著运动还是探头相对于像场的显著运动都会产生r.f.信号差值,该r.f.信号差值在合并来自不同多线的信号时会导致对消而非预期的信号增强。这一问题在具有其间为合并采集更多数量多线的更长采集时间间隔的情况下更为严重。图7的示例通过检测相对的图像运动以及相应于存在运动调整多线的阶数解决了这个问题。运动检测器130接收来自图像处理器122的连续图像帧,并比较图像内容找出指示运动的差别。一种适于执行这种操作的技术为使连续图像的一个或多个相同区域中图像像素值相关联的相关技术。其可以通过计算两幅图像数据间的零延迟标准化互相关来实现。如果相关因子较高(即,图像几乎完全相同),则使用更多数量的多线(这需要更多数量的发射-接收时间间隔)来提高场深。当相关因子较低(图像不同)时,其可能是由于像场中的运动或由探头引起,就合并较少数量的多线。适合的相关技术为美国专利6,126,598中描述的绝对差最少求和(MSAD)技术,该文献以引用方式并入本文。如该专利所描述的,为运动检测经MSAD处理的数据可为图像帧期间或图像帧之间发射的图像区域或参考线。也可使用多普勒技术检测运动,美国专利5,127,409(Daigle)对此进行了详细描述。
将运动检测的结果耦合到多线控制器132,该控制器根据运动调整多线的发射和接收,当存在运动时使用低位的多线。本示例中所作的调整为控制发射波束形成器,以便在存在运动时发射横向包绕较少多线的F数字更大的波束,使得所合并的多线跨越的发射时期更短。可以使用具有较大聚焦区域的较小发射孔径,其在较大范围的深度上更加均匀,从而需要较少发射波束来覆盖视场。接收时,多线控制器132控制所采用的多线处理器的数量,存在运动时接收的多线越少需要多线处理器的数量越少。多线控制器132也耦合到线存储器112,以便当存在运动时命令加权器114、延迟线118和加法器120对较少数量的多线进行重聚焦。这些调整的结果为需要更短的时间来采集待合并的多线,与较长采集时间间隔相比,较短的时期区域受运动的影响较小。因此,当运动检测器检测到运动影响较低的稳定像场时,多线控制器增加多线阶数。示例性调整范围为当像场稳定时使用16×多线(响应一个发射波束接收十六个多线以及合并来自十六个发射波束的多线),随着运动量的增加,降低多线阶数到8×、4×、2×和1×(单条线进行发射和接收)。例如,当在受检者身体上移动探头以搜索特定解剖结构或好的观看视角时,将使用标准单线成像。同样可以根据超声***的成像模式自动改变多线阶数。例如,当***由B模式成像变为多普勒成像时,可由多线控制器中断多线采集。多线的等级可在成像过程中由多线控制器自动设定,和/或可向用户提供手动控制来设定多线等级。
本领域的技术人员将很容易想到上述示例落入本发明的范围内的变更。例如,代替沿着连续的线位置进行发射的是,发射可跳过多线发射间的各线(多路复用)。发射可沿每隔一个线位置、每四个线位置或具有其他多线间隔的时间间隔进行,从而减少形成图像所需的发射事件的数量,并提高采集率。这也是降低运动伪影的方法。用较少数量的发射波束接收以较宽线间隔的时间间隔的多线是另一种解决运动问题的方法。该方法可用上面所提到的内插线进行扩增来增加所要合并的线的数量。该方法还会提高给定数量多线的场深的改善。所用的取样通常是发射孔径数量F的函数,其决定尼奎斯特(Nyquist)采样要求。不一定在图像上顺序地进行扫描;可以使用其它的发射序列,这将在存在运动时具有不同的响应。虽然上述示例具有变迹,但是没有必要对发射、接收或多线合并进行变迹。不必将待合并的多线全部精确地轴向对准。待合并的多线相互间可有横向偏移,从而在该多线合并时实现内插。本发明的原理也可通过执行在仰角尺寸和方位尺寸上的处理而应用于三维成像,也可应用于相控阵列(扇形扫描)成像以及线性阵列(正交扫描)成像并利用曲线阵列。
Claims (20)
1.一种用于对像场进行成像的超声诊断成像***,包括:
探头,其包括换能器元件阵列;
发射波束形成器,其耦合到所述换能器元件阵列,用于发射包绕多个横向间隔线位置的波束;
多线接收波束形成器,其耦合到所述换能器元件阵列,用于响应一个发射波束而在所述线位置处生成所述像场的多条接收线;
多线合并电路,其对由所述多线接收波束形成器响应多个发射波束而生成的与合用线位置相关的多条接收线的信号做出响应,从而用于以作为深度(Z)的函数而变化的相位特性曲线对接收线的信号进行延迟并合并所述多条接收线以生成发射重聚焦图像数据;
运动检测器,其用于检测所述像场中运动的存在;以及
显示器,其利用所述图像数据生成图像;
其中,由所述多线合并电路所合并的接收线的信号的数量是所述运动存在的函数。
2.如权利要求1所述的超声诊断成像***,其中,所述多线合并电路包括加法器,所述加法器用于合并来自所选数量的不同多线的数据。
3.如权利要求2所述的超声诊断成像***,其中,所述多线合并电路还包括对所述不同多线的数据做出响应的多个延迟器,其用于调整所述多线间的相位变化。
4.如权利要求3所述的超声诊断成像***,其中,所述多线合并电路还包括多个加权电路,其用于对多线数据进行加权。
5.如权利要求1所述的超声诊断成像***,其中,所述运动检测器还包括对来自不同图像的图像数据做出响应的相关器。
6.如权利要求5所述的超声诊断成像***,其中,所述相关器用于使用绝对差计算的求和生成相关因子。
7.如权利要求1所述的超声诊断成像***,其中,所述发射波束形成器进一步用于发射呈现如下波束宽度的波束,所述波束宽度为在所述像场中所检测到的运动的函数。
8.如权利要求1所述的超声诊断成像***,其中,所述多线接收波束形成器进一步用于生成一数量的接收线,所述接收线的数量为所述像场中所检测到的运动的函数。
9.一种用于在像场中存在运动的情况下生成具有扩展聚焦范围的超声图像的方法,包括:
从阵列换能器发射多个发射波束,每个发射波束集中于沿着所述换能器的阵列的不同位置,并且每个发射波束包绕多个横向间隔线位置,所述横向间隔线位置与另外波束的横向间隔线位置空间相关;
用所述阵列换能器接收回波信号;
并行处理响应一个发射波束所接收到的所述回波信号,以生成位于所述波束的所述横向间隔线位置处的多条接收线的回波信号;
针对额外的发射波束重复所述并行处理步骤;
检测所述像场中运动的存在;
向接收线的信号施加作为深度(Z)的函数而变化的相位特性曲线;
合并来自空间相关的不同发射波束的各接收线的回波信号以生成发射重聚焦图像数据,所述被合并的来自发射波束的各接收线的回波信号的数量与所述像场中运动的存在相关;以及
使用所述图像数据生成图像。
10.如权利要求9所述的方法,其中,所述运动的存在由超声***操作者通过观察图像进行检测,
其中,所述被合并的来自发射波束的各接收线的回波信号的数量通过调整超声***控制进行设定。
11.如权利要求9所述的方法,其中,检测所述运动存在的所述步骤还包括通过比较不同图像的图像数据来检测所述运动存在。
12.如权利要求11所述的方法,其中,所述比较步骤还包括通过计算绝对差求和使不同图像的图像数据相关。
13.如权利要求9所述的方法,其中,检测所述运动存在的所述步骤还包括通过多普勒技术检测所述运动存在。
14.如权利要求9所述的方法,其中,所述发射步骤还包括发射呈现如下波束宽度的多个发射波束,所述波束宽度为所述像场中所检测到的运动的函数。
15.如权利要求9所述的方法,其中,所述并行处理步骤还包括响应一个发射波束生成一数量的接收线,所述接收线的数量为所述像场中所检测到的运动的函数。
16.如权利要求9所述的方法,还包括存储响应不同发射波束所生成的一数量的接收线,所述接收线的数量为在所述像场中所检测到的运动的函数。
17.一种用于对像场进行成像的超声诊断成像***,包括:
探头,其包括换能器元件阵列;
发射波束形成器,其耦合到所述换能器元件阵列,用于发射包绕多个横向间隔线位置的波束;
多线接收波束形成器,其耦合到所述换能器元件阵列,用于响应一个发射波束在所述线位置处生成所述像场的多条接收线;
多线合并电路,其对由所述多线接收波束形成器响应多个发射波束而生成的与合用线位置相关的多条接收线的信号做出响应,从而用于以作为深度(Z)和所述发射波束中心的位置的函数而变化的相位特性曲线对接收线的信号进行延迟并合并所述多条接收线以生成发射重聚焦图像数据;
运动检测器,其用于检测所述像场中运动的存在;以及
显示器,其利用所述图像数据生成图像;
其中,由所述发射波束形成器所发射的所述波束为所述运动存在的函数。
18.如权利要求17所述的超声诊断成像***,其中,所述发射波束的横向尺寸为所述运动存在的函数。
19.如权利要求17所述的超声诊断成像***,其中,所述发射波束的横向间隔为所述运动存在的函数。
20.一种用于对像场进行成像的超声诊断成像***,包括:
探头,其包括换能器元件阵列;
发射波束形成器,其耦合到所述换能器元件阵列,用于发射包绕多个横向间隔线位置的波束;
多线接收波束形成器,其耦合到所述换能器元件阵列,用于响应一个发射波束在所述线位置处生成所述像场的多条接收线;
多线合并电路,其对由所述多线接收波束形成器响应多个发射波束而生成的与合用线位置相关的多条接收线的信号做出响应,从而用于以作为深度(Z)和所述发射波束中心的位置的函数而变化的相位特性曲线对接收线的信号进行延迟并合并所述多条接收线以生成发射重聚焦图像数据;
运动检测器,其用于检测所述像场中运动的存在;以及
显示器,其利用所述图像数据生成图像;
其中,由所述多线接收波束形成器所使用的所述线位置的间隔为所述运动存在的函数。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US74714906P | 2006-05-12 | 2006-05-12 | |
US60/747,149 | 2006-05-12 | ||
PCT/US2007/066752 WO2007133879A2 (en) | 2006-05-12 | 2007-04-17 | Ultrasonic synthetic transmit focusing with motion compensation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101442937A CN101442937A (zh) | 2009-05-27 |
CN101442937B true CN101442937B (zh) | 2012-07-18 |
Family
ID=38694586
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2007800167563A Active CN101442937B (zh) | 2006-05-12 | 2007-04-17 | 具有运动补偿的超声合成发射聚焦 |
Country Status (6)
Country | Link |
---|---|
US (2) | US9107564B2 (zh) |
EP (1) | EP2019625B1 (zh) |
JP (1) | JP4975098B2 (zh) |
KR (1) | KR20090016666A (zh) |
CN (1) | CN101442937B (zh) |
WO (1) | WO2007133879A2 (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5074963B2 (ja) * | 2008-03-13 | 2012-11-14 | 日立アロカメディカル株式会社 | 超音波診断装置 |
JP5419404B2 (ja) * | 2008-09-04 | 2014-02-19 | キヤノン株式会社 | 光音響装置 |
CN101664321B (zh) * | 2009-09-07 | 2011-06-15 | 无锡祥生科技有限公司 | 组织声速实时可调的超声诊断设备及其波束合成方法 |
EP2566394B1 (en) | 2010-05-03 | 2016-12-14 | Koninklijke Philips N.V. | Ultrasonic tracking of ultrasound transducer(s) aboard an interventional tool |
WO2013116783A1 (en) * | 2012-02-03 | 2013-08-08 | Los Alamos National Security, Llc | Windowed time-reversal music technique for super-resolution ultrasound imaging |
WO2013116807A1 (en) | 2012-02-03 | 2013-08-08 | Los Alamos National Security, Llc | Systems and methods for synthetic aperture ultrasound tomography |
JP6000196B2 (ja) * | 2012-09-27 | 2016-09-28 | 富士フイルム株式会社 | 超音波診断装置、音速決定方法およびプログラム |
KR20140056597A (ko) | 2012-10-30 | 2014-05-12 | 삼성전자주식회사 | 빔 포밍 계수 결정 방법, 빔 포밍 방법 및 초음파 이미징 장치 |
JP6047041B2 (ja) * | 2013-03-22 | 2016-12-21 | 富士フイルム株式会社 | 超音波診断装置、超音波診断装置の信号処理方法およびプログラム |
JP6331297B2 (ja) * | 2013-09-10 | 2018-05-30 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置、及び超音波測定方法 |
CN103536316B (zh) * | 2013-09-22 | 2015-03-04 | 华中科技大学 | 一种空时平滑相干因子类自适应超声成像方法 |
KR102164456B1 (ko) * | 2014-02-13 | 2020-10-12 | 삼성전자주식회사 | 초음파 측정 장치 및 초음파 측정 방법 |
JP5851549B2 (ja) * | 2014-04-21 | 2016-02-03 | 日立アロカメディカル株式会社 | 超音波診断装置 |
WO2015177658A1 (en) * | 2014-05-23 | 2015-11-26 | Koninklijke Philips N.V. | Motion gated-ultrasound thermometry using adaptive frame selection. |
US9763646B2 (en) * | 2014-06-12 | 2017-09-19 | General Electric Company | Method and systems for adjusting a pulse generated for ultrasound multi-line transmit |
WO2016098929A1 (ko) * | 2014-12-19 | 2016-06-23 | 삼성전자주식회사 | 초음파 영상 장치 및 그 제어방법 |
JP6406019B2 (ja) * | 2015-01-09 | 2018-10-17 | コニカミノルタ株式会社 | 超音波信号処理装置、及び超音波診断装置 |
JP2017046811A (ja) * | 2015-08-31 | 2017-03-09 | セイコーエプソン株式会社 | 超音波デバイス、超音波モジュール、及び超音波測定機 |
EP3475727B1 (en) | 2016-06-23 | 2021-11-03 | Koninklijke Philips N.V. | Rapid synthetic focus ultrasonic imaging with large linear arrays |
EP3263036A1 (en) | 2016-06-30 | 2018-01-03 | Esaote S.p.A. | Method and system for performing retrospective dynamic transmit focussing beamforming on ultrasound signals |
EP3507617B1 (en) * | 2016-09-02 | 2023-10-11 | Koninklijke Philips N.V. | Ultrasound probe with multiline digital microbeamformer |
US20200187910A1 (en) * | 2017-05-09 | 2020-06-18 | The University Of North Carolina At Chapel Hill | Adaptive multifocus beamforming ultrasound methods and systems for improved penetration and target sensitivity at high frame-rates |
CN111343925B (zh) | 2017-10-12 | 2023-10-13 | 皇家飞利浦有限公司 | 具有患者自适应剪切波生成的超声剪切波成像 |
CA3096219A1 (en) * | 2018-04-09 | 2019-10-17 | Butterfly Network, Inc. | Methods and apparatuses for offloading ultrasound data |
CN112166339B (zh) * | 2018-05-15 | 2024-07-05 | 皇家飞利浦有限公司 | 具有声速像差校正的合成透射聚焦超声*** |
CN112119327A (zh) * | 2018-05-15 | 2020-12-22 | 皇家飞利浦有限公司 | 具有声速像差映射的合成透射聚焦超声*** |
JP2024518474A (ja) | 2021-05-10 | 2024-05-01 | コーニンクレッカ フィリップス エヌ ヴェ | コヒーレントに合成された超音波画像生成、並びに関連するシステム、方法、及び装置 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4604697A (en) | 1983-08-05 | 1986-08-05 | Interspec, Inc. | Body imaging using vectorial addition of acoustic reflection to achieve effect of scanning beam continuously focused in range |
JPH078492A (ja) * | 1993-06-28 | 1995-01-13 | Toshiba Corp | 超音波診断装置 |
US5793701A (en) * | 1995-04-07 | 1998-08-11 | Acuson Corporation | Method and apparatus for coherent image formation |
JPH10118063A (ja) * | 1996-10-24 | 1998-05-12 | Aloka Co Ltd | 超音波診断装置 |
US5873830A (en) * | 1997-08-22 | 1999-02-23 | Acuson Corporation | Ultrasound imaging system and method for improving resolution and operation |
KR100280197B1 (ko) | 1997-11-10 | 2001-02-01 | 이민화 | 초음파영상화시스템의초음파신호집속방법및장치 |
US6126598A (en) * | 1998-10-01 | 2000-10-03 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging system with adaptive spatial compounding |
US6390980B1 (en) | 1998-12-07 | 2002-05-21 | Atl Ultrasound, Inc. | Spatial compounding with ultrasonic doppler signal information |
US6228031B1 (en) | 1999-02-17 | 2001-05-08 | Atl Ultrasound | High frame rate ultrasonic diagnostic imaging systems with motion artifact reduction |
US6432056B1 (en) * | 1999-10-08 | 2002-08-13 | Atl Ultrasound | Ultrasonic diagnostic imaging system with high frame rate synthetic transmit focus |
US6626834B2 (en) | 2001-01-25 | 2003-09-30 | Shane Dunne | Spiral scanner with electronic control |
US6508770B1 (en) * | 2001-03-08 | 2003-01-21 | Acuson Corporation | Aperture compounding for medical imaging |
US6482157B2 (en) * | 2001-03-30 | 2002-11-19 | Koninklijke Philips Electronics N.V. | Ultrasonic diagnostic imaging systems with blended multiline for 2D and 3D applications |
US6589177B1 (en) * | 2002-11-15 | 2003-07-08 | Koninklijke Philips Electronics N.V. | Method and apparatus for obtaining B-flow and B-mode data from multiline beams in an ultrasound imaging system |
ATE487149T1 (de) * | 2003-12-16 | 2010-11-15 | Koninkl Philips Electronics Nv | Diagnostisches ultraschall-bildgebungsverfahren und -system mit einer automatischen steuerung von auflösung und bildwiederholrate |
JP4625663B2 (ja) * | 2004-08-18 | 2011-02-02 | 株式会社東芝 | 超音波診断装置 |
US8317712B2 (en) * | 2006-05-12 | 2012-11-27 | Koninklijke Philips Electronics N.V. Eindhoven | Retrospective dynamic transmit focusing for spatial compounding |
EP2019626B1 (en) * | 2006-05-12 | 2015-09-02 | Koninklijke Philips Electronics, N.V. | Incoherent retrospective dynamic transmit focusing |
-
2007
- 2007-04-17 WO PCT/US2007/066752 patent/WO2007133879A2/en active Application Filing
- 2007-04-17 JP JP2009509922A patent/JP4975098B2/ja active Active
- 2007-04-17 CN CN2007800167563A patent/CN101442937B/zh active Active
- 2007-04-17 KR KR1020087027392A patent/KR20090016666A/ko not_active Application Discontinuation
- 2007-04-17 US US12/299,633 patent/US9107564B2/en active Active
- 2007-04-17 EP EP07760750.5A patent/EP2019625B1/en active Active
-
2015
- 2015-07-08 US US14/793,817 patent/US9345455B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
EP2019625B1 (en) | 2016-06-08 |
US9345455B2 (en) | 2016-05-24 |
CN101442937A (zh) | 2009-05-27 |
EP2019625A4 (en) | 2010-01-27 |
EP2019625A2 (en) | 2009-02-04 |
US9107564B2 (en) | 2015-08-18 |
US20090099456A1 (en) | 2009-04-16 |
WO2007133879A2 (en) | 2007-11-22 |
JP2009536854A (ja) | 2009-10-22 |
KR20090016666A (ko) | 2009-02-17 |
US20150305720A1 (en) | 2015-10-29 |
WO2007133879A3 (en) | 2008-02-21 |
JP4975098B2 (ja) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101442937B (zh) | 具有运动补偿的超声合成发射聚焦 | |
CN101442938B (zh) | 具有多线波束形成器的超声合成发射聚焦 | |
EP2019600B1 (en) | Retrospective dynamic transmit focusing for spatial compounding | |
JP2009536853A5 (zh) | ||
GB2099582A (en) | Ultrasonic image methods and apparatus | |
WO2009151142A1 (en) | Ultrasonic apparatus and control method therefor | |
JP7044723B2 (ja) | 大型線形アレイを備える高速合成集束超音波イメージング | |
US20080030581A1 (en) | Multizone Color Doppler Beam Transmission Method | |
EP1247490A1 (en) | Ultrasonic diagnosing apparatus | |
JP4632517B2 (ja) | 送受信波の整相方法及び超音波診断装置 | |
JP6396508B2 (ja) | 超音波トランスデューサアセンブリ | |
JPH0616784B2 (ja) | 超音波診断装置 |
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 |