The invention relates generally to ultrasound imaging using multiple ultrasound transmit beams, and more particularly to isolating ultrasound transmit beams and reducing cross-transmit beam interference in a multi-beam system using a Doppler method.
Diagnostic Ultrasound is one of the most versatile, least expensive, and widely used diagnostic imaging modalities in use today. With the advent of three-dimensional ultrasound and Doppler Tissue Imaging (DTI), much effort has been invested in increasing the frame rate in ultrasound imaging. One particular method involves receive multi-line beam processing where numerous ultrasound receive beams are calculated for each transmit beam or event.
A problem with this method is that to receive energy along a given scan line direction, ultrasound transmit energy needs to be supplied along that line of sight. To solve this problem, there are basically two approaches.
The first approach involves widening or “fattening” the transmit beam so that it encompasses a larger area or volume. This technique suffers from decreased resolution (both detailed and contrast) and from decreased sensitivity.
The second approach involves transmitting or “firing” multiple focused and compact transmit beams into the human body simultaneously. The problem with this method is cross-transmit beam interference (i.e., a form of cross-talk), that is, energy from one transmit beam contaminates the receive beams clustered along another transmit beam, and vice versa.
Several solutions have been presented to solve this problem of cross-transmit beam interference. Some of these solutions include aggressive nulling of the receive beamform to exclude energy from other transmit beams, coded excitation, spatial diversity, that is, placing the transmit beams as far apart as possible, and frequency diversity. For example, U.S. Pat. No. 6,179,780 describes various methods for overcoming the problem of cross talk, including using a receive beam synthesizer, using coded transmissions, using non-uniform scanning sequences, and using different transmit center frequencies. To the inventor's knowledge, these methods have not, as of yet, been employed commercially.
The present invention provides a solution to cross-transmit beam interference in a multi-beam system by providing a novel method of isolating the energy from the desired transmit beam, and the means for mitigating the energy and susceptibility to the “other” transmit beam(s).
The inventive method for isolating ultrasound transmit beams and reducing cross-transmit beam interference in a multi-beam system comprises the steps of performing a first transmit event by simultaneously transmitting at least two of ultrasound beams at disjoint spatial locations, each of the transmitted ultrasound beams generating an echo return; generating a sequence of transmit events; applying a phase factor to each of the transmitted ultrasound beams in each transmit event; in each successive transmit event, modulating the phase factor by a unique amount for each of the transmitted ultrasound beams; and, linearly combining the echo returns from two or more transmit events by constructively adding energy from a desired transmitted ultrasound beam and destructively interfering energy from the remaining transmitted ultrasound beams
The foregoing and other objects, aspects, features, and advantages of the invention will become more apparent from the following description and from the claims.
The invention is further described in the detailed description that follows, by reference to the noted drawings by way of non-limiting illustrative embodiments of the invention. As should be understood, however, the invention is not limited to the precise arrangements and instrumentalities shown. In the drawings, like reference characters generally refer to the same parts throughout the different views. Also, the drawings are not necessarily to scale, emphasis instead generally being placed upon illustrating the principles of the invention.
Reference will now be made in detail to the preferred embodiments of the present invention. While the invention will be described in conjunction with the preferred embodiments, it will be understood that they are not intended to limit the invention to these embodiments. On the contrary, the invention is intended to cover alternatives, modifications and equivalents, which may be included within the spirit and scope of the invention as defined by the appended claims. Furthermore, in the following detailed description of the present invention, numerous specific details are set forth in order to provide a thorough understanding of the present invention. However, it will be recognized by one of ordinary skill in the art that the present invention may be practiced without these specific details. In other instances, well known methods, procedures, components, and circuits have not been described in detail as not to unnecessarily obscure aspects of the present invention.
Referring to
Using dynamic receive beamforming, receive beams 108 and 114 are acquired or received by the ultrasound transmitter/receiver 102. However, the receiver 102 also receives a beam or signal 116, which is a reflected component of return beam or signal 114. The signal 116 contaminates the return beam or signal 108. Likewise, the receiver 102 also receives a beam or signal 110, which is a reflected component of return signal 108. The signal 110 contaminates the return signal 114. This cross-contamination of return signals 108 and 114 is referred to as cross-transmit beam interference, and degrades the contrast resolution of the ultrasound image.
In order to remove the contaminating signals 114, and 116 from return signals 108, and 110, respectively, a two coefficient finite image response (FIR) is applied to each of the return signals 108, 110, 114, 116, according to the equations A and B shown below.
((B3+N1)+(B4+(−N2)))/2=((B3+B4)/2)+((N1−N2)/2)=average of B3 and B4 Equation A:
((B1+N3)−(((−B2)+N4))/2=((B1−B2)/2)−((N3−N4)/2)=average of B1 and B2 Equation B:
Where B1, B2, B3, B4 are transmit beams and N1, N2, N3, N4 are nodes.
In a simple embodiment such as that shown in
Thus, in this simple embodiment, we will only have round trip reconstructed beams at odd degree values (as corresponding to the example table above). Focusing only on the constructive interference of the “Good” or non-contaminated energy produces the following equations:
RT
−43=(+R43X−44−−R−43X−42)/2
RT
−41=(−R−41X−42−+R−41X−40)/2
RT
−39=(+R39X−40−−R39X−38)/2
Where:
RT
1=(+R1X0++R1X2)/2
RT
3=(+R3X2++R3X4)/2
RT
5=(+R5X4++R5X6)/2
Note that the desired energy associated with RT−43 has the transmit toggling in polarity every other transmit beam (+, −, +, −). Hence the “minus” sign in its equation. Conversely, the sign to coherently add the energy for RT1 is associated with transmit beams that are always of the same polarity. Hence coherent summation requires the receive beams be “summed”.
The above equations are an oversimplification of what really happens, because the negative degree roundtrip beams, e.g. RT−43, are also susceptible to “Bad” or contaminated energy from the positive degree transmit events, and vice versa. The following equation includes the effect of the “Bad” energy.
RT
−43={(+R−43X−44+BADR−43X1)−(−R−43X−42+BADR−43X3)}/2
Rearranging the terms in this equation yields:
RT
−43={(+R−43X−44+BADR−43X1)−(−R−43X−42+BADR−43X3)}/2
The desired “Good” energy in the first half of the equation adds coherently, whereas the “Bad” energy from the 2nd half of the equation is appropriately destructed. This is easy to see for the other “negative” degreed angles.
The technique illustrated with the above equations will also work for the positive degreed roundtrip angles as shown below.
RT
1={(+R1X0+BADR1X−44)+(+R1X2−BADR1X−42)}/2
Rearranging the terms in this equation yields:
RT
1={(+R1X0+R1X2)+(BADR1X−44−BADR1X−42)}/2
Again one can see that the bad energy from the opposite side transmit beams are appropriately canceled out.
In more advanced and preferred embodiments, there will be numerous receive beams for each transmit event, and in the simple positive/negative polarity case, the span of the receive beams will overlap each other by fifty percent.
In the embodiment shown in
In a preferred embodiment, we would have eight or more receive beams per transmit beam, and the overlap would be seventy-five percent or greater. This is illustrated in
The equation defining how to combine the receive beams for group 250 is:
RT
250
=a*X
1
R
7
−b*X
2
R
5
+c*X
3
R
3
−d*X
4
R
1
There are some constraints on the coefficients that should improve the performance and achieve the desired results.
Constrain #1: The sum of the coefficients should equal one:
a+b+c+d=1
This causes the average energy in the multiple receive beams to have unity gain.
Constrain #2: The coefficients should interpolate to a location between the X2 and X3 transmit beams, and should be specifically located closer to X2 (as is graphically indicated in
1*a+2*b+3*c+4*d=2.25
Note that the 1, 2, 3, 4 correspond to the spatial locations of the transmit beams X1, X2, X3, and X4, and the value 2.25 corresponds to the desired location of the interpolated output.
Constrain #3: The coefficients need to cancel out the energy from non-toggling polarity transmit beams from group 260 in
a−b+c−d=0
One solution that meets the above constraints is:
a=0.025
b=0.60
c=0.475
d=−0.10
For the group of receive lines defined by 255 (to right of group 250), the coefficients can be swapped, yielding:
RT
255
=d*X
1
R
7
−c*X
2
R
5
+b*X
3
R
3
−a*X
4
R
1
Note that swapping the coefficients will modify Constraint #2, such that the resultant output beam will interpolate to “2.75” (still between X2 and X3, but now closer to X3).
Likewise, these coefficients can be applied to the groups 260 and 265 (to right of 260):
RT
260
=a*X
101
R
7
+b*X
102
R
5
+c*X
103
R
3
+d*X
104
R
1
RT
265
=d*X
101
R
7
+c*X
102
R
5
+b*X
103
R
3
+a*X
104
R
1
Note the difference in the sign of the coefficients.
As will be obvious to one skilled in the art, the round trip beams defined by RT250, RT255, RT260, and RT265 will be accurately located and will reject leakage energy from the “other” group of transmit beams.
A further embodiment of this invention is its use in conjunction with U.S. Provisional Patent Application No. 60/747,148, titled “ULTRASONIC SYNTHETIC TRANSMIT FOCUSING WITH A MULTILINE BEAMFORMER”, incorporated herein by reference. In this case, one can describe the RT260 round-trip beam as follows:
RT
260(t)=a*X1R7(t−d1)+b*X2R5(t−d2)+c*X3R3(t−d3)+d*X4R1(t−d4)
In this equation, “t” refers to the time during which the ultrasound echoes are coming from increasing depths in the body, and the delays, d1, d2, d3, d4, are calculated to retrospectively beamform the transmit beam as defined in the above provisional patent application. By applying Constrain #3 (a−b+c−d=0) to the above RT260(t) equation, one can achieve the benefits of improved transmit focusing and mitigation of energy from the undesired transmit beams.
Referring to
Referring to
Next, the transmitter 202 transmits a third beam at a positive five-degree angle and transmitter 204 transmits a simultaneous third beam at a negative forty-one degree angle (Step 412). Receiver 220 receives a third return signal and a third reflected signal, and receiver 222 also receives a third return signal and a third reflected signal (Step 414). The data processing unit again executes the signal averaging algorithm to determine the average of return signal 208b and the third return signal, and the average of the return signal 214b and the other third return signal (Step 412). This sequence of steps repeats until the desired tissue area (not shown) has been scanned.
The aforementioned embodiments all involve two simultaneous transmit beams, such that one sequence of beams maintains a normal polarity, while the second set of transmit beams toggle in polarity. An aspect of this invention is to support more than two transmit beams, such that for any given transmit beam sequence, energy from all other transmits is mitigated. The following example will demonstrate a simultaneous four beam sequence. Four simultaneous transmit beams 510 can be co-planar for scanning a 2D image, as is illustrated in
In one embodiment, assume there are four beam sequences, referred to as Xa, Xb, Xc, and Xd. Each beam will cover different portions of the scanned regions. Furthermore, each beam will proceed through four different transmit waveforms.
For Xa, shown in
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)
Note that “t” refers to time, “n” refers to the transmit event, “f” refers to the nominal transmit frequency (e.g. 5.0 MHz), and “w(t)” refers to a time windowing function. For the example in
Also, for Xb, shown in
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)
However, for Xc (and Xd), one needs yet another sequence that can be uniquely distinguished. In this case, one can advance (or retard) the “phase” of the transmit waveform. Xc, which uses an advancing phase term, as shown in
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)
And for Xd, which uses a retarding phase term, as seen in
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)
For purposes of illustrating this particular embodiment, each of the four transmit beam sequences will simultaneously receive four beams per transmit, as is shown in
RT
XA@2.5
=a*X
A1
R
4
+b*X
A2
R
3
+c*X
A3
R
2
+d*X
A4
R
1
Since there are more concurrent transmit beams than in prior embodiments, there will be some additional constraints on the selection of the a,b,c,d coefficients.
Solving for a,b,c,d yields the very simple result:
a=b=c=d=0.25
For someone skilled in the art, it would be a simple matter to come up with similar sets of coefficients for the other transmits: Xb, Xc, and Xd.
Hence it is a further inventive aspect of this invention to allow simultaneous transmit beams to be used with matrix transducers using micro-beamformers. One inventive element replicates the micro-beamformer electronics, one for each simultaneous transmit beam. For example, in the case where two beams are simultaneously transmitted, there will be two micro-beamformers per patch (per group of elements). Each micro-beamformer will produce a distinct transmit wave field, will be combined with the transmit wave fields from the other micro-beamformers associated with a single patch, will be amplified, and will drive the patch elements to send sound waves into the body (see
It is implicit in all of the aforementioned embodiments that they were designed for use in a “fundamental” mode of black-and-white gray scale imaging. Fundamental mode is where the transmit frequency is the same as the receive frequency. There is another mode of operation, referred to as Tissue Harmonic Imaging (THI), which is quite common in the current clinical practice of diagnostic ultrasound. In THI, harmonic frequencies are generated during the transmission and propagation of the transmit waveforms. These harmonics (often the second harmonic) are then selectively isolated on receive using bandpass filters. For example, the transmit waveform might be centered at 2.5 MHz, and the receive filters are set to 5.0 MHz to selectively receive the desired second harmonic.
In THI, in order to reject cross-beam contamination from simultaneous transmit as described by this invention, one needs to control the transmit in such a way that the desired phase relationship is observed on receive. For example, in the 2× multi-beam transmit embodiment, it is desired that the first sequence of beams have common receive phase, whereas the second set of transmits have the polarity of the receive signal toggle by 180 degrees every other transmit. In order to achieve this 180 degree toggling on the receive harmonic, the transmits for this sequence need to toggle between 0 and 90 degrees. In other words, the transmit sequence would toggle between a windowed cosine burst and a windowed sine burst. In the 4× multi-beam transmit embodiment, the various transmit sequences would need to be advanced (or retarded) by 45 degrees to achieve the desired 90 degree shift on receive (for the second harmonic).
As would be known to someone skilled in the art, the transmit phase shift would be approximately 1/H of the desired phase shift observed on receive, where “H” refers to the receive harmonic. Also known to one skilled in the art is that this phase relationship is not always exact, and may need to be finely adjusted based upon empirical measurements.
In a preferred embodiment, the data processing unit can be an FPGA (field programmable gate array), or an ASIC (application specific integrated circuit). The processing can also be performed using DSPs (digital signal processing units) or other computational units. In a preferred embodiment, two transmitter/receivers are used with one of the transmission beams switching between zero and 180 degree phases. In other embodiments, three or more ultrasound transmitters are used with the transmission beams transmitting at 0, 90, 180, and 270 degree phases. In still another embodiment, one beam would always be in phase (zero degrees), one beam would advance at +90 increments, one beam would advance at −90 increments, and one beam would toggle between 0 and 180 degrees.
Variations, modifications, and other implementations of what is described herein may occur to those of ordinary skill in the art without departing from the spirit and scope of the invention. Accordingly, the invention is not to be defined only by the preceding illustrative description.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB07/55235 | 12/19/2007 | WO | 00 | 6/18/2009 |
Number | Date | Country | |
---|---|---|---|
60871014 | Dec 2006 | US |