CN101568304A - 多波束发射隔离 - Google Patents
多波束发射隔离 Download PDFInfo
- Publication number
- CN101568304A CN101568304A CNA2007800469277A CN200780046927A CN101568304A CN 101568304 A CN101568304 A CN 101568304A CN A2007800469277 A CNA2007800469277 A CN A2007800469277A CN 200780046927 A CN200780046927 A CN 200780046927A CN 101568304 A CN101568304 A CN 101568304A
- Authority
- CN
- China
- Prior art keywords
- emission
- ultrasonic
- wave
- transmit events
- beams
- 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
Links
Images
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/52077—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 with means for elimination of unwanted signals, e.g. noise or interference
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/07—Analysing solids by measuring propagation velocity or propagation time of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/46—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by spectral analysis, e.g. Fourier analysis or wavelet analysis
-
- 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/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/8927—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array using simultaneously or sequentially two or more subarrays or subapertures
-
- 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/52079—Constructional features
- G01S7/5208—Constructional features with integration of processing functions inside probe or scanhead
-
- 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/5209—Details related to the ultrasound signal acquisition, e.g. scan sequences using multibeam transmission
-
- 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/5209—Details related to the ultrasound signal acquisition, e.g. scan sequences using multibeam transmission
- G01S7/52093—Details related to the ultrasound signal acquisition, e.g. scan sequences using multibeam transmission using coded signals
-
- 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
-
- 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/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/8925—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being a two-dimensional transducer configuration, i.e. matrix or orthogonal linear arrays
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Signal Processing (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Medical Informatics (AREA)
- Mathematical Physics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Radiology & Medical Imaging (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Transducers For Ultrasonic Waves (AREA)
Abstract
一种用于在多波束系统中隔离超声发射波束和减小交叉发射波束干涉的方法包括发射在第一和第二正角的第一超声波束和发射在第一和第二负角的第二超声波束。该方法还涉及接收第一、第二、第三和第四复合信号,其中,每个复合信号包括返回信号和反射分量。该方法还包括将有限脉冲响应滤波器应用于第一和第三复合信号以及第二和第四复合信号以获得第一与第二复合信号的平均值和第二与第四复合信号的平均值并且去除反射分量。
Description
本发明大体涉及使用多个超声发射波束的超声成像,尤其涉及在使用多普勒方法的多波束系统中进行隔离超声发射波束和减小交叉发射波束干涉。
诊断超声是目前使用的最通用的、最便宜的和最广泛使用的诊断成像模式之一。伴随着三维超声和多普勒组织成像(DTI)的出现,人们在增加超声成像中的帧速方面付出了很多努力。一个特殊方法涉及接收多线波束处理,其中,针对每个发射波束或事件计算许多超声接收波束。
该方法存在的问题是为了沿着给定扫描线方向接收能量,超声发射能量需要沿着该视线被提供。为了解决该问题,基本上有两种方法。
第一方法涉及加宽或“展平”发射波束,使得其包围较大的面积或体积。该技术遭遇到分辨率减小(在细节和对比度两方面)和灵敏度减小的问题。
第二方法涉及同时将多个聚焦和密集发射波束发射或“击发”到人体中。该方法存在的问题是交叉发射波束干涉(即,一种串扰形式),也就是说,来自一个发射波束的能量污染沿另一发射波束聚集的接收波束,反之亦然。
提出了几种方案以解决该交叉发射波束干涉的问题。这些方案中的一些包括对接收波形主动置零以排除来自其他发射波束的能量、编码激励、空间分集,也就是说,尽可能远地放置发射波束和频率分集。例如,美国专利6,179,780描述了用于克服串扰问题的各种方法,包括使用接收波束合成器、使用编码传输、使用非均匀扫描序列以及使用不同发射中心频率。据申请人所知,这些方法到目前为止还未在商业上被利用。
本发明通过提供隔离来自预期发射波束的能量的新颖方法和用于减小能量和对“其他”(一个或多个)发射波束的易感性的装置而提供多波束系统中交叉发射波束干涉的解决方案。
用于在多波束系统中隔离超声发射波束和减小交叉发射波束干涉的本发明的方法包括以下步骤:通过在不相交的空间位置同时发射至少两个超声波束执行第一发射事件,每个发射超声波束生成回波返回;生成发射事件的序列;将相位因子应用于每个发射事件中的每个发射超声波束;在每个连续发射事件中,针对每个发射超声波束以唯一量调制相位因子;以及通过相长地增加来自预期发射超声波束的能量和相消地干涉来自其余发射超声波束的能量线性地组合来自两个或以上发射事件的回波返回。
根据以下描述以及根据权利要求,本发明的前述和其他目标、方面、特征和优势将更显而易见。
经由本发明的非限定性说明性实施例通过参考附图在以下详细描述中进一步描述了本发明。然而应当理解,本发明并不限于所示的精确布置和手段。在图中,在所有不同视图中相似的附图标记通常表示相同的部分。而且,附图不一定按比例绘制,而是通常强调突出显示本发明的原理。
图1是根据本发明的一个实施例的被定位成扫描人体组织的超声波束发射器的说明性示意图;
图2A是根据本发明的一个实施例的接收和发射波束的说明性示意图;
图2B是根据本发明的另一实施例的接收和发射波束的说明性示意图;
图2C是根据本发明的另一实施例的接收和发射波束的说明性示意图;
图3是根据本发明的一个实施例的超声发射事件、角和极性的说明性图表;
图4是根据本发明的一个实施例的用于在多波束系统中隔离发射超声波束和减小交叉发射波束干涉的方法的说明性流程图;
图5A是用于扫描2D图像的共面的四个同时发射波束的说明性示意图;
图5B是用于扫描体积的不共面的四个同时发射波束的说明性示意图;
图6A示出了当发射波形相同时的发射波形序列;
图6B示出了当每隔一次发射极性切换(polarity toggle)时的发射波形序列;
图6C示出了当发射波形使用超前相位项时的发射波形序列;
图6D显示了当发射波形使用滞后相位项时的发射波形序列;
图7是根据本发明的另一实施例的接收和发射波束的说明性示意图;
图8A是将声波发送到身体的不同发射波场的说明性示意图;以及
图8B是从身体返回的各面元的回波的累加的说明性示意图。
现在将详细参考本发明的优选实施例。尽管将结合优选实施例描述本发明,应当理解其并非旨在将本发明限制于这些实施例。相反地,本发明旨在涵盖可以被包括在由从属权利要求限定的本发明的精神和范围内的替换选择、修改和等同物。此外,在本发明的以下详细描述中,阐述许多特定细节以提供对本发明的彻底理解。然而,本领域的普通技术人员应当认识到本发明可以在没有这些特定细节的情况下被实施。在其他情况下,为了不必要地使本发明的方面晦涩,详细地描述了公知方法、程序、部件和电路。
参考图1,在简单的实施例中,对于每个扫描帧或扫描体积,利用两个同时超声发射波束;在其他实施例中,更多的同时超声发射波束被利用并且将在下面被论述。图1显示了超声发射器/接收器102,并且粗实箭头106、112对应于被定位成扫描人体组织的两个同时发射波束。围绕带有箭头106、112的这些粗线的实线104、120、122、124示出了大约6dB的能量波束宽度,其有效地限定对应于该轴向深度的发射波束的宽度(分辨率)。使用动态接收波束形成,由使用虚线的箭头示出的四个同时接收波束108、110、114、116被采集。图1包含针对每个发射波束106、112的两个接收波束114、116、108、110。多个同时发射事件被激发以在整个2D图像上扫描,或者在体积的情况下,在体积的侧向和高程尺度上扫描。超声发射器102产生在正四十五度角的一个超声波束106,和在负四十五度角的另一超声波束112。
使用动态接收波束形成,接收波束108和114由超声发射器/接收器102采集或接收。然而,接收器102也接收作为返回波束或信号114的反射分量的波束或信号116。信号116污染返回波束或信号108。类似地,接收器102也接收作为返回信号108的反射分量的波束或信号110。信号110污染返回信号114。返回信号108和114的这一交叉污染被称为交叉发射波束干涉,并且降低超声图像的对比度分辨率。
为了分别从返回信号108和110去除污染信号114和116,根据下面所示的方程A和B将两个系数有限图像响应(FIR)应用于每个返回信号108、110、114、116。
方程A:((B3+N1)+(B4+(-N2)))/2=((B3+B4)/2)+((N1-N2)/2)=B3和B4的平均值。
方程B:((B1+N3)-(((-B2)+N4))/2=((B1--B2)/2)-((N3-N4)/2)=B1和B2的平均值。
其中,B1、B2、B3、B4是发射波束而N1、N2、N3、N4是节点。
在诸如图1中所示的简单实施例中,我们可以假设每个发射波束有两个接收波束或信号,并且我们可以进一步假设当发射波束序列在视场上出现时,接收波束将重叠。下面的简表示出了简单实施例序列。
发射事件 | 发射A | 发射A极性 | 接收X | 接收Y | 发射B | 发射B极性 | 接收U | 接收V |
1 | -44度 | + | -45度 | -43度 | 0度 | + | -1度 | 1度 |
2 | -42 | - | -43 | -41 | 2 | + | 3 | 1 |
3 | -40 | + | -41 | -39 | 4 | + | 5 | 3 |
4 | -38 | - | -39 | -37 | 6 | + | 7 | 5 |
图2A对应于该简表,示出了本发明的简单实施例。图2A示出了对应于发射波束150、160的向下实箭头,和对应于接收波束位置165、168的向上虚箭头。假设在左边的发射事件150在极性上切换,而在右边的发射事件160保持相同极性。
因此,在该简单实施例中,我们将仅仅具有在奇度值的往返重建波束(其对应于上面的例表)。仅仅关注“良好”或无污染能量的相长干涉产生以下方程:
RT-43=(+R-43X-44--R-43X-42)/2
RT-41=(-R-41X-42-+R-41X-40)/2
RT-39=(+R-39X-40--R-39X-38)/2
其中:
RT-43是在-43度的往返波束位置。
R-43X-44是与@-44度的发射波束关联的在@-43度的接收波束。
并且,同时解出与发射波束“B”关联的往返:
RT1=(+R1X0++R1X2)/2
RT3=(+R3X2++R3X4)/2
RT5=(+R5X4++R5X6)/2
应当注意,与RT-43关联的预期能量每隔一次发射波束在极性上进行发射切换(+,-,+,-)。因此其方程中出现“减”号。相反地,相干地增加RT1的能量的符号与总是具有相同极性的发射波束关联。因此相干累加需要接收波束被“累加”。
以上方程是实际情况的过度简化,原因是例如RT-43的负度数往返波束也易于受到来自正度数发射事件的“有害”或污染能量,反之亦然。以下方程包括“有害(BAD)”能量的影响。
RT-43={(+R-43X-44+BADR-43X1)-(-R-43X-42+BADR-43X3)}/2
重排该方程中的项得到:
RT-43={(+R-43X-44+R-43X1)+(BADR-43X1-BADR-43X3)}/2
方程的第一半部中的预期“良好”能量相干地增加,而来自方程的第二半部的“有害”能量适当地被破坏。对于其他“负”度数角这容易看到。
随着以上方程示出的技术也将适用于如下所示的正度数往返角。
RT1={(+R1X0+BADR1X-44)+(+R1X2-BADR1X-42)}/2
重排该方程中的项得到:
RT1={(+R1X0+R1X2)+(BADR1X-44-BADR1X-42)}/2
再次可以看到来自相对侧发射波束的有害能量适当地被消除。
在更高级和优选的实施例中,将有针对每个发射事件的许多接收波束,并且在简单的正/负极性情况下,接收波束的跨度将彼此重叠百分之五十。图2B示出了每个发射波束有四个接收波束,其中,接收波束的跨度彼此重叠百分之五十。在图2B中,与图2A中相同,向下实箭头对应于发射波束210、220,并且向上虚箭头对应于接收波束位置230、240。与简单实施例相同,假设在左边的发射事件210在极性上切换,而在右边的发射事件240保持相同极性。
在图2B所示的实施例中,交叉波束抑制被减小,这是由于向正确往返波束位置“插值”需要使用诸如1/4,3/4的系数,这导致“良好”能量的正确放置,但是“有害”能量仅仅被减小6dB(减小1/2)。
在优选实施例中,我们将让每个发射波束具有八个或以上接收波束,并且重叠将为百分之七十五或以上。这在图2C中示出。圆形区域250、260示出了如何从对应于四个不同发射事件212的相同成角接收波束重建往返波束。由于往返波束将具有与它关联的四个不同系数,即,四分接口插值滤波器(four tap interpolation filter),抑制来自其他发射波束的“有害”能量的能力将被提高。
定义如何组合组250的接收波束的方程为:
RT250=a*X1R7-b*X2R5+c*X3R3-d*X4R1
对系数的一些约束条件将提高性能并获得预期结果。
约束条件#1:系数的总和应当等于一:
a+b+c+d=1
这导致多个接收波束中的平均能量具有单位增益。
约束条件#2:系数应当被插值到X2与X3发射波束之间的位置,并且应当特别地被定位成更靠近X2(其在图2C中被示意性地示出)。以方程形式对此进行描述得到:
1*a+2*b+3*c+4*d=2.25
应当注意1、2、3、4对应于发射波束X1、X2、X3和X4的空间位置,而值2.25对应于插值输出的预期位置。
约束条件#3:系数需要消除来自图2C中的组260的非切换极性发射波束的能量。这可以通过切换系数的极性并且保证它们总和为零而实现:
a-b+c-d=0
满足以上约束条件的一个解为:
a=0.025
b=0.60
c=0.475
d=-0.10
对于由255限定的接收线的组(在组250的右边),系数可以被交换,得到:
RT255=d*X1R7-c*X2R5+b*X3R3-a*X4R1
应当注意,交换系数将修改约束条件#2,使得结果输出波束将被插值到“2.75”(仍然在X2与X3之间,但是现在更靠近X3)。
类似地,这些系数可以被应用于组260和265(在260的右边):
RT260=a*X101R7+b*X102R5+c*X103R3+d*X104R1
RT265=d*X101R7+c*X102R5+b*X103R3+a*X104R1
应当注意系数的符号的差异。
本领域的技术人员将显而易见,由RT250、RT255、RT260和RT265定义的往返波束将被精确地定位并且将抑制来自发射波束的“其他”组的泄漏能量。
本发明的另一实施例是它与在此引用作为参考的名称为“ULTRASONICSYNTHETIC TRANSMIT FOCUSING WITH A MULTILINE BEAMFORMER”的美国临时专利申请No.60/747,148结合使用。在该情况下,可以如下描述RT260往返波束:
RT260(t)=a*X1R7(t-d1)+b*X2R5(t-d2)+c*X3R3(t-d3)+d*X4R1(t-d4)
在该方程中,“t”表示超声回波来自身体中的增加深度的持续时间,并且延迟d1、d2、d3、d4被计算以回溯地波束形成如以上临时专利申请中定义的发射波束。通过将约束条件#3(a-b+c-d=0)应用到以上RT260(t)方程,可以获得的好处是提高发射聚焦和减小来自非期望发射波束的能量。
参考图3,在一个实施例中,示出了包括发射角302、304和发射信号的极性306、308的超声发射事件301(发射信号的实例)的表。对于发射器204,发射角302以+2度的增量从-45度递增到-1度,并且发射信号的极性304保持为负(即,同相)。对于发射器202,发射角306以2度的增量从+1度递增到+45度,并且发射信号的极性308从正切换到负(即,反相180度),使得每隔一次信号发射与前一信号发射反相180度。
参考图3和4,针对每个发射器202和204重复先前描述的方法。例如,发射器/接收器202发射在正一度角的波束206a并且发射器204同时发射在负四十五度角的波束212a(步骤402)。接收器220接收返回信号208a和反射信号216a并且接收器222接收返回信号214a和反射信号210a(步骤404)。发射器202接着发射在负三度角的波束206b并且发射器204同时发射在负四十三度角的波束212b(步骤406)。接收器220接收返回信号208b和反射信号216b并且接收器222接收返回信号214b和反射信号210b(步骤408)。诸如计算机的数据处理单元执行信号平均算法以确定返回信号208a与208b和返回信号214a与214b的平均值(步骤410)。
接着,发射器202发射在正五度角的第三波束并且发射器204同时发射在负四十一度角的第三波束(步骤412)。接收器220接收第三返回信号和第三反射信号,并且接收器222也接收第三返回信号和第三反射信号(步骤414)。数据处理单元再次执行信号平均算法以确定返回信号208b与第三返回信号的平均值,和返回信号214b与另一第三返回信号的平均值(步骤412)。该步骤序列重复直到预期组织区域(未显示)被扫描。
前述实施例全部涉及两个同时发射波束,使得波束的一个序列保持正常极性,而第二组发射波束在极性上切换。本发明的一个方面是支持两个以上发射波束,使得对于任何指定发射波束序列,来自所有其他发射的能量被减小。以下例子将示范同时四个波束序列。四个同时发射波束510可以共面以用于扫描2D图像,如图5A中所示,或者它们可以不共面520,以用于扫描体积,如图5B中所示。为了发射不共面发射波束,使用元件530的2D矩阵换能器,如图5B中所示。应当注意以下例子同时适用于共面和不共面的情况。“有害”能量的抑制在时域中发生,因此交叉污染发射波束在空间中位于什么地方并不重要。
在一个实施例中,假设有四个波束序列,被称为Xa、Xb、Xc和Xd。每个波束将覆盖扫描区域的不同部分。此外,每个波束将行进通过四个不同发射波形。
对于图6A中所示的Xa,发射波形将是相同的。这些可以被表达为:
Xa(t,n=1)=cos(2*pi*f*t)*w(t)
Xa(t,n=2)=cos(2*pi*f*t)*w(t)
Xa(t,n=3)=cos(2*pi*f*t)*w(t)
Xa(t,n=4)=cos(2*pi*f*t)*w(t)
应当注意,“t”表示时间,“n”表示发射事件,“f”表示额定发射频率(例如5.0MHz),以及“w(t)”表示时间窗口函数。对于图6a、6b、6c和6d中的例子,w(t)可以是矩形窗口函数,其仅仅在-0.4到+0.4微秒时是开启的(=1)。在5MHz,这将导致仅仅具有四个周期的发射波形。假设w(t)对于所有发射序列(Xa、Xb、Xc和Xd)是相同的。此外,假设该四波形序列将重复,使得第五波形将使用波形#1:Xa(t,n=5)=Xa(t,n=1)。
同样,对于图6B中所示的Xb,发射波形将使用先前的方法,其中,每隔一次发射极性切换。这可以表达为:
Xb(t,n=1)=+cos(2*pi*f*t)*w(t)
Xb(t,n=2)=-cos(2*pi*f*t)*w(t)
Xb(t,n=3)=+cos(2*pi*f*t)*w(t)
Xb(t,n=4)=-cos(2*pi*f*t)*w(t)
然而,对于Xc(和Xd),需要可以被唯一区分的又一序列。在该情况下,可以超前(或滞后)发射波形的“相位”。如图6C中所示的使用超前相位项的Xc可以表达为:
Xc(t,n=1)=+cos(2*pi*f*t)*w(t)
Xc(t,n=2)=+sin(2*pi*f*t)*w(t)
Xc(t,n=3)=-cos(2*pi*f*t)*w(t)
Xc(t,n=4)=-sin(2*pi*f*t)*w(t)
并且对于如图6D中可见的使用滞后相位项的Xd,表达式为:
Xd(t,n=1)=+cos(2*pi*f*t)*w(t)
Xd(t,n=2)=-sin(2*pi*f*t)*w(t)
Xd(t,n=3)=-cos(2*pi*f*t)*w(t)
Xd(t,n=4)=+sin(2*pi*f*t)*w(t)
为了示出该特定实施例,四个发射波束序列中的每一个每次发射将同时接收四个波束,如图7中所示。以下方程对应于发射Xa的接收线的包围组:
RTXA@2.5=a*XA1R4+b*XA2R3+c*XA3R2+d*XA4R1
由于有比现有实施例中更多的同时发射波束,对a、b、c、d系数的选择将有一些附加约束条件。
约束条件1:a+b+c+d=1累加来自Xa的相关能量
约束条件2:a-b+c-d=0抑制来自Xb的能量
约束条件3:a+jb-c-jd=0抑制来自Xc的能量
约束条件4:a-jb-c+jd=0抑制来自Xd的能量
应当注意,“j”表示-1平方根的虚部,并且对应于与发射Xc和Xd相关的90度相移。
求解a、b、c、d得到非常简单的结果:
a=b=c=d=0.25
对于本领域的计算人员来说,得到其他发射Xb、Xc和Xd的系数的类似集合将是简单的。
图5B示出了使用2D矩阵换能器530使用四个同时发射波束扫描体积。在矩阵换能器上,期望使用全采样孔(所有元件电活性)提高图像质量和灵敏度。这与仅仅连接少量元件的稀疏阵列形成对比。可以通过使用定位于矩阵换能器的外壳中的微波束形成器获得全采样阵列。参见在此处引用作为参考的美国专利5,997,479和6,126,602。每个微波束形成器将适当地波束形成被称为面元的元件的小子集。如本领域的技术人员目前已知的,微波束形成器的使用将与同时发射波束和与本发明不相容。这是因为元件的每个面元或组在发射和接收时被限制到单一转向角。并且,本发明中的暗示是使用可以空间分离并且非共同定位的多个发射。
因此,本发明的另一个创造性方面是允许同时发射波束与使用微波束形成器的矩阵换能器一起使用。一个创造性元件重复微波束形成器电子装置,每个同时发射波束用一个。例如,在两个波束同时被发射的情况下,每个面元(每个元件组)将有两个微波束形成器。每个微波束形成器将产生不同发射波场,将与来自与单一面元相关联的其他微波束形成器的发射波场组合,将被放大,并且将驱动面元元件将声波发送到身体中(参见图8A)。另外,在接收时,共用面元元件将把返回声波转换成电信号,将被放大,并且将被发送到N个不同的微波束形成器。每个波束形成器然后将对沿着与发射期间使用的方向关联的方向的返回斑块回波进行延迟和累加(参见图8b)。在一般情况下,“N”个同时发射波束探视方向将需要“N”个微波束形成器。
在所有前述实施例中暗示了它们被设计用于黑白灰度级成像的“基本”模式中。基本模式是发射频率与接收频率相同。有另一操作模式,被称为组织谐波成像(THI),其在诊断超声的当前临床实践中是很常见的。在THI中,在发射波形的发射和传播期间生成谐波频率。这些谐波(常常是二次谐波)然后在接收时使用带通滤波器选择性地隔离。例如,发射波形可以以2.5MHz为中心,而接收滤波器被设置为5.0MHz以选择性地接收预期二次谐波。
在THI中,为了如本发明中所述抑制来自同时发射的交叉波束污染,需要以这样的方式控制发射使得在接收时观察到预期相位关系。例如,在2x多波束发射实施例中,期望波束的第一序列具有共同的接收相位,而第二组波束每隔一次发射使接收信号的极性切换180度。为了在接收谐波上获得该180度切换,该序列的发射需要在0和90度之间切换。换句话说,发射序列将在窗口余弦脉冲串和窗口正弦脉冲串之间切换。在4x多波束发射实施例中,各种发射序列将需要超前(或滞后)45度以在接收时获得预期90度相移(对于二次谐波)。
本领域的技术人员将会知道,发射相移将大约为在接收时观察到的预期相移的1/H,其中,“H”表示接收谐波。本领域技术人员也知道该相位关系并不总是正好的,而是需要基于经验测量精细地调整。
在优选实施例中,数据处理单元可以是FPGA(场可编程门阵列),或ASIC(特定用途集成电路)。处理也可以使用DSP(数字信号处理单元)或其他计算单元执行。在优选实施例中,两个发射器/接收器与在零和180度相位之间切换的发射波束之一一起使用。在其他实施例中,三个或以上超声发射器用于在0、90、180、和270度相位发射的发射波束。在又一实施例中,一个波束将总是同相(零度),一个波束将超前+90增量,一个波束将超前-90增量,且一个波束将在0和180度之间切换。
本领域的普通技术人员可以想到在此处描述的变化、修改和其他实现方式而不脱离本发明的精神和范围。因此,本发明并不仅仅由前面的示例性描述限定。
Claims (21)
1、一种用于在多波束系统中隔离超声发射波束和减小交叉发射波束干涉的方法,所述方法包括:
通过在不相交的空间位置同时发射一定数量的超声波束执行第一发射事件,所述数量为至少两个,所述至少两个发射超声波束的每一个生成回波返回;
生成发射事件随着时间的序列;
将相位因子应用于每个发射事件中的所述至少两个发射超声波束的每一个;以及
在每个连续发射事件中,针对每个发射超声波束以唯一量调制所述相位因子,其中,通过相长地增加来自预期发射超声波束的能量和相消地干涉来自其余发射超声波束的能量对来自两个或更多发射事件的所述回波返回进行组合。
2、根据权利要求1所述的方法,其中,所述多个发射超声波束等于每个发射事件有两种发射超声波束,并且其中,针对所述发射超声波束之一将所述相位因子简化为{+1+1+1...},并且针对其他发射超声波束将所述相位因子简化为{+1-1+1-1...}。
3、根据权利要求1所述的方法,其中,所述不相交的空间位置以与相控阵扇区和弯曲线阵式换能器中的一个相关联的度数被限定。
4、根据权利要求1所述的方法,其中,所述不相交的空间位置以与线性换能器相关联的侧向距离相偏移。
5、根据权利要求1所述的方法,其中,所述不相交的空间位置对应于不同发射焦深。
6、根据权利要求1所述的方法,其中,所述连续发射事件顺序地扫描2D图像和3D体积中的一个。
7、根据权利要求1所述的方法,其中,在累加节点形成接收波束之后所述至少两个发射超声波束被隔离。
8、根据权利要求7所述的方法,在所述形成接收波束期间还使用并行处理以针对所述至少两个发射超声波束的每一个产生一个或多个接收波束。
9、根据权利要求8所述的方法,其中,所述接收波束的每一个具有用于组合来自连续发射事件的能量的系数的唯一集合,其中,来自所述预期发射超声波束的能量被相长地增加,而来自其他非期望发射超声波束的能量被相消地干涉。
10、根据权利要求1所述的方法,其中,多波束系统包含利用微波束形成电子装置的超声换能器。
11、根据权利要求10所述的方法,其中,所述微波束形成电子装置波束形成至少一个面元,并且每个面元的组内处理器针对N个空间不相交的发射超声波束的每一个重复N次。
12、根据权利要求1所述的方法,其中,可以使用关于发射和接收的至少一个的时间延迟对相位因子调制进行近似。
13、根据权利要求1所述的方法,其中,使用组织谐波成像调制所述相位因子。
14、根据权利要求13所述的方法,其中,所述组织谐波成像包括至少两个谐波分量,并且应用于所述发射波束的相位因子调制量基本上被减半,其中,在接收期间在2xRF的观测相位因子通过与所述至少两个谐波的第二次谐波关联的非线性波传播被有效地加倍。
15、根据权利要求13所述的方法,其中,对于发射波形的第M次谐波分量,应用于所述发射波束的相位因子调制量基本上被减半,其中,在接收期间在第M次接收谐波M x Fxmit F观察到的观测相位因子通过与所述组织谐波成像的第二次谐波关联的非线性波传播被有效地加倍。
16、一种用于在超声成像中允许更快帧速的方法,所述方法包括:
使用具有一个或多个微波束形成器的矩阵阵列超声换能器同时发射多个超声波束,其中,矩阵换能器包括超声元件的2D阵列,所述超声元件包含换能器外壳中的电子装置以执行波束形成的某个方面,所述换能器外壳中的电子装置支持在不相交的空间位置中被波束形成的独立的和分离的同时发射的超声波束。
17、根据权利要求16所述的方法,还包括:
生成发射事件随着时间的序列,每个发射事件包括在不相交的空间位置的同时发射的多个超声波束,所述发射的超声波束的每一个生成回波返回;
将相位因子应用于每个发射事件中的所述发射的超声波束的每一个;以及
在每个连续发射事件中,针对每个发射的超声波束以唯一量调制相位因子,其中,通过相长地增加来自预期发射超声波束的能量和相消地干涉来自其余发射超声波束的能量对来自两个或以上发射事件的回波返回进行组合。
18、根据权利要求16所述的方法,其中,具有同时发射的波束的至少两个微波束形成器包括面元。
19、根据权利要求16所述的方法,其中,每个微波束形成器产生不同发射波场,并且所述面元中的所述微波束形成器的所述不同发射波场是可组合的。
20、根据权利要求16所述的方法,其中,每个面元的组内处理器针对空间不相交的发射超声波束的每一个重复N次。
21、根据权利要求16所述的方法,其中,可以使用关于发射和接收的至少一个的时间延迟对所述相位因子调制进行近似。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US87101406P | 2006-12-20 | 2006-12-20 | |
US60/871,014 | 2006-12-20 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101568304A true CN101568304A (zh) | 2009-10-28 |
Family
ID=39536811
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2007800469277A Pending CN101568304A (zh) | 2006-12-20 | 2007-12-19 | 多波束发射隔离 |
Country Status (6)
Country | Link |
---|---|
US (1) | US20100016725A1 (zh) |
EP (1) | EP2120721A2 (zh) |
JP (1) | JP2010514322A (zh) |
KR (1) | KR20090094447A (zh) |
CN (1) | CN101568304A (zh) |
WO (1) | WO2008075302A2 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106999149A (zh) * | 2014-11-25 | 2017-08-01 | 皇家飞利浦有限公司 | 多传感器超声探头及相关方法 |
Families Citing this family (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008051639A2 (en) | 2006-10-25 | 2008-05-02 | Maui Imaging, Inc. | Method and apparatus to produce ultrasonic images using multiple apertures |
US9282945B2 (en) | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
US10226234B2 (en) | 2011-12-01 | 2019-03-12 | Maui Imaging, Inc. | Motion detection using ping-based and multiple aperture doppler ultrasound |
US9788813B2 (en) | 2010-10-13 | 2017-10-17 | Maui Imaging, Inc. | Multiple aperture probe internal apparatus and cable assemblies |
US8235906B2 (en) * | 2007-12-31 | 2012-08-07 | The Brigham And Women's Hospital, Inc. | System and method for accelerated focused ultrasound imaging |
US8602993B2 (en) * | 2008-08-08 | 2013-12-10 | Maui Imaging, Inc. | Imaging with multiple aperture medical ultrasound and synchronization of add-on systems |
EP2419023A4 (en) * | 2009-04-14 | 2013-01-16 | Maui Imaging Inc | UNIVERSAL MEDICAL MULTIAPERTUR ULTRASOUND HEAD |
WO2010120907A2 (en) | 2009-04-14 | 2010-10-21 | Maui Imaging, Inc. | Multiple aperture ultrasound array alignment fixture |
US8500007B2 (en) * | 2009-10-02 | 2013-08-06 | Giftcodes.Com, Llc | System and method for merchant interaction with and tracking of the secondary gift card marketplace |
JP6274724B2 (ja) * | 2010-02-18 | 2018-02-07 | マウイ イマギング,インコーポレーテッド | 多開口超音波撮像を用いた点音源送信及び音速補正 |
US9668714B2 (en) | 2010-04-14 | 2017-06-06 | Maui Imaging, Inc. | Systems and methods for improving ultrasound image quality by applying weighting factors |
JP6092109B2 (ja) | 2010-10-13 | 2017-03-08 | マウイ イマギング,インコーポレーテッド | 凹面超音波トランスデューサ及び3dアレイ |
EP2638860B1 (en) * | 2010-11-09 | 2020-02-26 | Konica Minolta, Inc. | Beam-forming method, ultrasonic diagnosis device, program, and integrated circuit |
US9265484B2 (en) | 2011-12-29 | 2016-02-23 | Maui Imaging, Inc. | M-mode ultrasound imaging of arbitrary paths |
CN104135937B (zh) | 2012-02-21 | 2017-03-29 | 毛伊图像公司 | 使用多孔超声确定材料刚度 |
CN104620128B (zh) | 2012-08-10 | 2017-06-23 | 毛伊图像公司 | 多孔径超声探头的校准 |
CN104582582B (zh) | 2012-08-21 | 2017-12-15 | 毛伊图像公司 | 超声成像系统存储器架构 |
US9510806B2 (en) | 2013-03-13 | 2016-12-06 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
US9883848B2 (en) | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
US9763646B2 (en) | 2014-06-12 | 2017-09-19 | General Electric Company | Method and systems for adjusting a pulse generated for ultrasound multi-line transmit |
CN106794007B (zh) | 2014-08-18 | 2021-03-09 | 毛伊图像公司 | 基于网络的超声成像系统 |
KR101683518B1 (ko) * | 2015-07-22 | 2016-12-07 | 기아자동차 주식회사 | 비접촉 내구 진단 장치 및 방법 |
US10987085B2 (en) | 2015-12-10 | 2021-04-27 | 1929803 Ontario Corp | Systems and methods for automated fluid response measurement |
EP3386398B1 (en) | 2015-12-10 | 2024-01-10 | 1929803 Ontario Corp. D/b/a Ke2 Technologies | Systems for automated fluid response measurement |
CN113729764A (zh) | 2016-01-27 | 2021-12-03 | 毛伊图像公司 | 具有稀疏阵列探测器的超声成像 |
WO2019169508A1 (en) | 2018-03-09 | 2019-09-12 | 1929803 Ontario Corp. D/B/A Flosonics Medical | Dynamically controllable patient fluid control device |
US11109831B2 (en) | 2018-07-17 | 2021-09-07 | 1929803 Ontario Corp, (o/a FloSonics Medical) | Ultrasound patch for detecting fluid flow |
WO2021216723A1 (en) | 2020-04-22 | 2021-10-28 | Bfly Operations, Inc. | Methods and apparatuses for beamforming in ultrasound systems |
US11937976B2 (en) | 2020-07-06 | 2024-03-26 | 1929803 Ontario Corp | Ultrasound patch with integrated flexible transducer assembly |
WO2024076975A1 (en) * | 2022-10-03 | 2024-04-11 | Mayo Foundation For Medical Education And Research | Fast ultrasound imaging method for motion detection using comb detection beams |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5997479A (en) * | 1998-05-28 | 1999-12-07 | Hewlett-Packard Company | Phased array acoustic systems with intra-group processors |
US6179780B1 (en) * | 1999-08-06 | 2001-01-30 | Acuson Corporation | Method and apparatus for medical diagnostic ultrasound real-time 3-D transmitting and imaging |
KR100393370B1 (ko) * | 2001-04-25 | 2003-07-31 | 주식회사 메디슨 | 직교 골레이 코드를 이용하는 초음파 영상 형성 방법 및장치 |
US20050148874A1 (en) * | 2003-12-19 | 2005-07-07 | Brock-Fisher George A. | Ultrasonic imaging aberration correction with microbeamforming |
-
2007
- 2007-12-19 US US12/519,902 patent/US20100016725A1/en not_active Abandoned
- 2007-12-19 JP JP2009542370A patent/JP2010514322A/ja not_active Withdrawn
- 2007-12-19 KR KR1020097012686A patent/KR20090094447A/ko not_active Application Discontinuation
- 2007-12-19 EP EP07859460A patent/EP2120721A2/en not_active Withdrawn
- 2007-12-19 WO PCT/IB2007/055235 patent/WO2008075302A2/en active Application Filing
- 2007-12-19 CN CNA2007800469277A patent/CN101568304A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106999149A (zh) * | 2014-11-25 | 2017-08-01 | 皇家飞利浦有限公司 | 多传感器超声探头及相关方法 |
CN106999149B (zh) * | 2014-11-25 | 2020-09-29 | 皇家飞利浦有限公司 | 多传感器超声探头及相关方法 |
Also Published As
Publication number | Publication date |
---|---|
US20100016725A1 (en) | 2010-01-21 |
WO2008075302A2 (en) | 2008-06-26 |
KR20090094447A (ko) | 2009-09-07 |
JP2010514322A (ja) | 2010-04-30 |
EP2120721A2 (en) | 2009-11-25 |
WO2008075302A3 (en) | 2008-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101568304A (zh) | 多波束发射隔离 | |
CN110063749B (zh) | 超声波测定装置、超声波图像装置及超声波测定方法 | |
US9198636B2 (en) | Continuous transmit focusing method and apparatus for ultrasound imaging system | |
EP2053420B1 (en) | Method of removing an effect of side lobes in forming an ultrasound synthetic image by motion estimation and compensation | |
US9645118B2 (en) | Ultrasonic imaging apparatus and ultrasonic imaging method | |
US7393324B2 (en) | Ultrasonic transmitting and receiving apparatus | |
US8245577B2 (en) | Pulse period jitter for artifact detection or reduction in ultrasound imaging | |
US5355888A (en) | High resolution phased array echo imager | |
US8038620B2 (en) | Fresnel zone imaging system and method | |
WO2008051639A3 (en) | Method and apparatus to produce ultrasonic images using multiple apertures | |
JP5702326B2 (ja) | 超音波プローブおよびそれを備える超音波診断装置 | |
EP1315981A2 (en) | Ultrasonic nonlinear imaging at fundamental frequencies | |
US6585648B1 (en) | System, method and machine readable program for performing ultrasonic fat beam transmission and multiline receive imaging | |
Nomura et al. | Feasibility of low-frequency ultrasound imaging using pulse compressed parametric ultrasound | |
US4552021A (en) | Electro-sound transducer eliminating acoustic multi-reflection, and ultrasonic diagnostic apparatus applying it | |
KR102245671B1 (ko) | 음향 방사력-기반 초음파 이미징에서의 적응식 클러터 필터링 | |
KR20160090059A (ko) | 초음파 프로브, 이를 포함하는 초음파 진단장치 및 초음파 신호생성방법 | |
US20110282204A1 (en) | Ultrasonic diagnostic apparatus | |
JP6396508B2 (ja) | 超音波トランスデューサアセンブリ | |
CN102908167A (zh) | 医用超声诊断仪的云成像处理系统及其方法 | |
US11914066B1 (en) | Multiplexed phased array multibeam sonar | |
Nikolov et al. | K-space model of motion artifacts in synthetic transmit ultrasound imaging | |
Kirkebø | Beamforming for Imaging: A Brief Overview | |
Tong et al. | Multi-transmit beam forming for fast cardiac imaging: quantitative analysis of the cross-talk between MLT beams | |
Singh et al. | A Modified Sum Focussing Approach to Ultrasonic Speckle Reduction |
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 |
Open date: 20091028 |