Multi-stage digital ultrasound beamformer

FIELD OF THE INVENTION

0001This invention relates to diagnostic ultrasound imaging in general, and to digital beamforming in particular.

BACKGROUND OF THE INVENTION

DESCRIPTION OF THE RELATED ART

0002The importance of diagnostic ultrasound imaging is widely recognized, and has grown as imaging resolution and the range of available uses and features have steadily increased. Once an expensive luxury available only in the best-equipped hospitals, diagnostic ultrasound imaging is now a commonly and almost routinely offered procedure even in some individual physician's offices.

0003Diagnostic medical ultrasound machines are among the most sophisticated signal-processing systems in the civilian world and among the few that can press modern processors to their performance limits. The quality of an ultrasound image is consequently directly affected by many factors relating to both the hardware involved and the routines used to generate the necessary pattern of ultrasonic signals transmitted into the body and to receive and process their echo returns.

0004<figref idref="f0001">Figure 1</figref> illustrates the main components of a standard ultrasonic imaging system. The user enters various conventional scan parameters into an input/output (I/O) unit 100, which may comprise conventional hardware, system software, and devices as a keyboard, knobs, a mouse, and/or buttons. The input unit is connected to a processing system 200, which will typically be an electrically connected and cooperating group of processors 210 such as microprocessors and digital signal processors with necessary system software such as, for example, an operating system. The I/O unit 100 may therefore not need any processing capability of its own but rather simply act as a peripheral device for the main processing system 200.

0005As in known systems, the processing system 200 sets, adjusts, and monitors the operating parameters of a conventional transmission control circuit 220. This control circuit 220 generates and applies electrical control and driving signals to an ultrasonic probe, that is, transducer 300, which includes an array 310 of electroacoustic (usually piezoelectric, such as PZT) elements. As is well known in the art, the piezoelectric elements generate ultrasonic waves when electrical signals of the proper frequency are applied to them.

0006To create the ultrasonic signals that are transmitted into the body (an "interrogation region," or a "region of interest" 400), all or some subset of these elements are electrically excited at a given frequency and are individually phased and time-shifted relative to each element's position in the array in such a way that the in-phase portions of the signals form a coherent transmit (Tx) beam: The time delays are chosen such that the echo signals arriving from a desired angular direction add coherently, whereas those signals arriving from other directions do not add coherently and so tend to cancel each other out. The amplitudes of the excitation signals applied to the array elements are typically varied (apodized) to control the resulting beam width and to suppress side lobes. The aperture of the array may also be varied, that is, the "width" of the activated elements (not all elements are usually activated at once) may be changed, just as the aperture of a camera may be varied to affect focal depth and provide uniform lateral resolution.

0007The various known methods for activating the elements according to a particular signal profile to best image a given area are generally referred to as transmit beamforming. By changing the amplitude and phasing of the transmitted signals, the transmit beam can be not only focused at a particular depth, but also steered off-axis, that is, not perpendicular to the center element among those activated. In <figref idref="f0001">Figure 1</figref>, for example, the transmit beam has been steered at an angle to the face of the array 310. Out-of-phase portions of the individual transmit signals will tend to cancel each other. The transmit beam typically converges at a focal depth, beyond which it once again diverges. The transmit beam is steered in the azimuth/lateral direction LAT and the elevation direction EL, and is focused in the depth/axial direction AX so as to concentrate the ultrasonic energy of the beam onto desired points, such as a structure 410, within the interrogation region 400.

0008After any reflecting structure in the interrogation region is insonified, the transducer is typically switched by a switching circuit 230 from the transmit mode to a receive mode, such that the ultrasound reflected back onto the array face from any structures 410 within the imaged area will cause the elements to create corresponding electrical signals. The ultrasound sensed by any particular element, however, will be a combination of all the ultrasound energy reflected back from anywhere in the imaged region, regardless of which elements the energy originated from. In the general case, each element will receive some ultrasound energy created by all elements that were energized to form the transmit beam. In the receive (Rx) mode, each element will therefore generate a time-varying signal that in practice will differ from that of all other elements' received signals. Again, however, the in-phase portion of the reflected energy will dominate. The problem is how best to determine just what these coherent portions are, that is, which sensed values at the different elements represent the return signal from the desired focal point that, with proper time delay, corresponds to the portions of the return signals at the other array elements. This process is known as receive beamforming and the primary task of a receive beamformer is to determine and sum the coherent signals from targets received by all the active Rx array elements.

0009As is well understood, the piezoelectric array elements 310 convert the small mechanical vibrations caused by the echo signal into corresponding radio-frequency (RF) electrical signals. Amplification and other conventional signal conditioning are then applied to the return signals by a reception controller 240. This processing includes various procedures as needed to identify the echo signals that correspond to each scanned element in the interrogation region and form the receive beam. The reception controller 240, all or part of which is normally integrated into the processing system 200 itself, also converts the ultrasonic, radio-frequency (RF) return signals (typically on the order of a few to tens of megahertz) into lower frequency ranges for processing. This is well known in the art of ultrasonic imaging. The reception control circuitry includes a receive beamformer, aspects of which are described in greater detail below. Note that the transmissioin and reception controllers (or, for that matter, the transmit and receive beamformers) so not have to be implemented as separate components but may be combined in whole or in part with respect to both their hardware and software.

0010In conventional B-mode scanning, each point within the interrogation region is then represented as an intensity (brightness) value. The interrogation region can therefore be represented as a discretized pattern (matrix) of brightness or signal intensity values, which are stored as frame data in a memory 250. Other scan modes, such as for Doppler imaging, create and store corresponding values.

0011The interrogation region is normally not in the same shape as what the user wants to see displayed; even when it is, the digital acoustic intensity values that make up the frame data are normally not in a form suitable for driving a conventional gray-tone or color display directly. The acoustic intensity values for a selected sub-set (scan plane) are then typically applied to a conventional scan converter 260, which converts the digital acoustic values into display intensity values that are suitable for use in driving a display device 120, which may be part of the I/O unit 100. The display device 120 typically includes a screen (for example, LED or CRT) that is divided into an X-Y (or polar) matrix or pattern of picture elements or "pixels" that make up an image that the user can view and interpret.

0012<figref idref="f0002">Figure 2</figref> illustrates the main components of a receive beamformer 500. The illustrated beamformer 500 is simplified for purposes of illustration and the sizes and distances are not to scale. For example, even a one-dimensional linear array will usually have scores or even hundreds of array elements and not just the eight shown in <figref idref="f0002">Figure 2</figref>. Moreover, it is common not to activate all the ultrasound elements at once, but rather to activate predetermined sub-sets of them according to some pattern. In <figref idref="f0002">Figure 2</figref>, the transmit beam is assumed to have been focused on a focal point FP; the lines from the point 500 back to the array elements represent the acoustic return signals, that is, the echoes of the transmitted ultrasound.

0013In the illustrated case, the focal point FP, for purposes of easy illustration only, lies perpendicular to the face of the uppermost array element. Assuming that the ultrasound is propagating through tissue that is at least substantially acoustically homogeneous, ultrasound energy of the transmit beam that is reflected back from the focal point FP will reach the uppermost array element before reaching the others. The dashed line 520 illustrates the travel of coherent ultrasound back towards the array elements 310; that is, line 520 illustrates the coherent portion of reflection of a transmit beam.

0014As the wavy lines 530 illustrate, the reflected ultrasound will cause each array element to generate a corresponding analog electrical signal. The dots on the wavy lines represent the instantaneous values that correspond to the values of received data from a respective one of the array elements. Usually after amplification by a bank of amplifiers 540, the analog signals from the electroacoustic array elements 310 are converted into digital values by a bank of analog-to-digital converters 550, whose sampling rate is controlled by clocking circuitry 555.

0015A time series of the converted digital values is then usually buffered in first-in-first-out latches 560 (for speed, although they could be stored in other types of memory devices). In <figref idref="f0002">Figure 2</figref>, each stored digital value is represented as a small square on the signal path within the latches 560. The solid (filled in) squares within the latches illustrate values corresponding to an echo from the same target.

0016The digital values are then processed in a processing module (570), which will generally involve the main processors 210 (<figref idref="f0001">Figure 1</figref>) executing the computer-readable software code that defines not only the implemented beamforming routine, but also such procedures as may involve a dedicated processor part of the greater processor to form an output signal that represents the value to be assigned to the focal point, for example, a brightness value. The processing will usually include such procedures as time-gating, gain compensation, diffraction compensation, dynamic receive delay, dynamic receive aperture size control, and dynamic receive apodization.

0017Because of its importance in ultrasonic imaging, there is a constant effort in the industry to improve receive beamformers. There are, indeed, hundreds of technical papers and even issued patents on the subject. These include the following issued U.S. patents: <ul id="ul0001" list-style="none" compact="compact"><li><u>U.S. Patent No.</u><u>Title</u></li><li><patcit id="pcit0001" dnum="US7199738B"><text>7,199,738</text></patcit>Sigma delta beamformer and method with reduced artifact</li><li><patcit id="pcit0002" dnum="US6937176B"><text>6,937,176</text></patcit> Ultrasonic signal acquisition in the digital beamformer</li><li><patcit id="pcit0003" dnum="US6716171B"><text>6,716,171</text></patcit> System and method for interfacing an ultrasound transducer with a computing device performing beamforming processing</li><li><patcit id="pcit0004" dnum="US6676602B"><text>6,676,602</text></patcit> Two dimensional array switching for beamforming in a volume</li><li><patcit id="pcit0005" dnum="US6635019B"><text>6,635,019</text></patcit> Scanhead assembly for ultrasonic imaging having an integral beamformer and demountable array</li><li><patcit id="pcit0006" dnum="US6602194B"><text>6,602,194</text></patcit> Dual beamformer ultrasound system for 2D and 3D imaging</li><li><patcit id="pcit0007" dnum="US6592524B"><text>6,592,524</text></patcit> Transmit beamformer delay architecture and method for diagnostic medical ultrasound</li><li><patcit id="pcit0008" dnum="US6565510B"><text>6,565,510</text></patcit> Method and apparatus for servicing remote ultrasound beamformer from central service facility</li><li><patcit id="pcit0009" dnum="US6500120B"><text>6,500,120</text></patcit> Beamforming system using analog random access memory</li><li><patcit id="pcit0010" dnum="US6491634B"><text>6,491,634</text></patcit> Sub-beamforming apparatus and method for a portable ultrasound imaging system</li><li><patcit id="pcit0011" dnum="US6447452B"><text>6,447,452</text></patcit> System and method for reduction of parallel beamforming artifacts</li><li><patcit id="pcit0012" dnum="US6436048B"><text>6,436,048</text></patcit> Ultrasonic diagnostic imaging system with scanhead elevation beamforming</li><li><patcit id="pcit0013" dnum="US6436044B"><text>6,436,044</text></patcit> System and method for adaptive beamformer apodization</li><li><patcit id="pcit0014" dnum="US6363033B"><text>6,363,033</text></patcit> Method and apparatus for transmit beamformer system</li><li><patcit id="pcit0015" dnum="US6245016B"><text>6,245,016</text></patcit> Ultrasound imaging system having post-beamformer signal processing using deconvolution algorithm</li><li><patcit id="pcit0016" dnum="US6172939B"><text>6,172,939</text></patcit> Method and apparatus for transmit beamformer system</li><li><patcit id="pcit0017" dnum="US6155980A"><text>6,155,980</text></patcit> Ultrasonic imaging system with beamforming using unipolar or bipolar coded excitation</li><li><patcit id="pcit0018" dnum="US6146328A"><text>6,146,328</text></patcit> Method and apparatus for ultrasonic beamforming using golay-coded excitation</li><li><patcit id="pcit0019" dnum="US6123671A"><text>6,123,671</text></patcit> Method and apparatus for distributed, agile calculation of beamforming time delays and apodization values</li><li><patcit id="pcit0020" dnum="US6113545A"><text>6,113,545</text></patcit> Ultrasonic beamforming with improved signal-to-noise ratio using orthogonal complementary sets</li><li><patcit id="pcit0021" dnum="US6110116A"><text>6,110,116</text></patcit> Method and apparatus for receive beamformer system</li><li><patcit id="pcit0022" dnum="US6108273A"><text>6,108,273</text></patcit> Transmit beamformer with frequency dependent focus</li><li><patcit id="pcit0023" dnum="US6104673A"><text>6,104,673</text></patcit> Method and apparatus for transmit beamformer system</li><li><patcit id="pcit0024" dnum="US6042547A"><text>6,042,547</text></patcit> Method and apparatus for receive beamformer system</li><li><patcit id="pcit0025" dnum="US6029116A"><text>6,029,116</text></patcit> Method and apparatus for a baseband processor of a receive beamformer system</li><li><patcit id="pcit0026" dnum="US5995450A"><text>5,995,450</text></patcit> Method and apparatus for transmit beamformer system</li><li><patcit id="pcit0027" dnum="US5984869A"><text>5,984,869</text></patcit> Method and apparatus for ultrasonic beamforming using golay-coded excitation</li><li><patcit id="pcit0028" dnum="US5928152A"><text>5,928,152</text></patcit> Method and apparatus for a baseband processor of a receive beamformer system</li><li><patcit id="pcit0029" dnum="US5921932A"><text>5,921,932</text></patcit> Method and apparatus for a baseband processor of a receive beamformer system</li><li><patcit id="pcit0030" dnum="US5882307A"><text>5,882,307</text></patcit> Method and apparatus for receive beamformer system</li><li><patcit id="pcit0031" dnum="US5856955A"><text>5,856,955</text></patcit> Method and apparatus for transmit beamformer system</li><li><patcit id="pcit0032" dnum="US5851187A"><text>5,851,187</text></patcit> Method and apparatus for ultrasonic beamforming with spatially encoded transmits</li><li><patcit id="pcit0033" dnum="US5844139A"><text>5,844,139</text></patcit> Method and apparatus for providing dynamically variable time delays for ultrasound beamformer</li><li><patcit id="pcit0034" dnum="US5827188A"><text>5,827,188</text></patcit> Method and apparatus for receive beamformer system</li><li><patcit id="pcit0035" dnum="US5797847A"><text>5,797,847</text></patcit> Method and apparatus for complex bandpass filtering and decimation in ultrasound beamformer</li><li><patcit id="pcit0036" dnum="US5720289A"><text>5,720,289</text></patcit> Method and apparatus for a geometric aberration transform in an adaptive focusing ultrasound beamformer system</li><li><patcit id="pcit0037" dnum="US5696737A"><text>5,696,737</text></patcit> Transmit beamformer with frequency dependent focus</li><li><patcit id="pcit0038" dnum="US5685308A"><text>5,685,308</text></patcit> Method and apparatus for receive beamformer system</li><li><patcit id="pcit0039" dnum="US5676147A"><text>5,676,147</text></patcit> Ultrasonic receive beamformer with phased sub-arrays</li><li><patcit id="pcit0040" dnum="US5675554A"><text>5,675,554</text></patcit> Method and apparatus for transmit beamformer</li><li><patcit id="pcit0041" dnum="US5653236A"><text>5,653,236</text></patcit> Apparatus for real-time distributed computation of beamforming delays in ultrasound imaging system</li><li><patcit id="pcit0042" dnum="US5617862A"><text>5,617,862</text></patcit> Method and apparatus for beamformer system with variable aperture</li><li><patcit id="pcit0043" dnum="US5608690A"><text>5,608,690</text></patcit> Transmit beamformer with frequency dependent focus</li><li><patcit id="pcit0044" dnum="US5581517A"><text>5,581,517</text></patcit> Method and apparatus for focus control of transmit and receive beamformer systems</li><li><patcit id="pcit0045" dnum="US5573001A"><text>5,573,001</text></patcit> Ultrasonic receive beamformer with phased sub-arrays</li><li><patcit id="pcit0046" dnum="US5570691A"><text>5,570,691</text></patcit> Method and apparatus for real-time, concurrent adaptive focusing in an ultrasound beamformer imaging system</li><li><patcit id="pcit0047" dnum="US5566675A"><text>5,566,675</text></patcit> Beamformer for phase aberration correction</li><li><patcit id="pcit0048" dnum="US5555534A"><text>5,555,534</text></patcit> Method and apparatus for doppler receive beamformer system</li><li><patcit id="pcit0049" dnum="US5551433A"><text>5,551,433</text></patcit> Method and apparatus for a geometric aberration transform in an adaptive focusing ultrasound beamformer system</li><li><patcit id="pcit0050" dnum="US5544128A"><text>5,544,128</text></patcit> Multi-beam digital beamforming method and apparatus</li><li><patcit id="pcit0051" dnum="US5517537A"><text>5,517,537</text></patcit> Integrated acoustic leak detection beamforming system</li><li><patcit id="pcit0052" dnum="US5388461A"><text>5,388,461</text></patcit> Beamforming time delay correction for a multi-element array ultrasonic scanner using beamsum-channel correlation</li><li><patcit id="pcit0053" dnum="US5388079A"><text>5,388,079</text></patcit> Partial beamforming</li><li><patcit id="pcit0054" dnum="US5369624A"><text>5,369,624</text></patcit> Digital beamformer having multi-phase parallel processing</li><li><patcit id="pcit0055" dnum="US4644795A"><text>4,644,795</text></patcit> High resolution multiline ultrasonic beamformer</li><li><patcit id="pcit0056" dnum="US4542653A"><text>4,542,653</text></patcit> Apparatus and method for beamforming in an ultrasonic transducer array</li></ul>

0018<patcit id="pcit0057" dnum="US20020143253A" dnum-type="L"><text>United States Published Patent Application 2002/0143253</text></patcit> discloses a method for multi-line acquisition in an ultrasound imaging system in which groups of beams are acquired and aligned in response to transmit pulses. After envelope detection, amplitude information is combined in a way so as to reduce multiline image artifacts.

0019A conventional digital ultrasound beamformer usually generates one transmit beam, but receives multiple receive beams back, since the single transmit beam may reflect back from different structures. Conventional single-stage beamformers are usually able to process these multiple beams by applying dynamic receive delay, dynamic receive aperture size control and dynamic receive apodization. Even so, prior art beamformers fail to take into consideration certain information that can be obtained from multiple receive beams over multiple transmits.

SUMMARY

0020The invention provides an ultrasound imaging system that comprises a main processing system including at least one processor and at least one memory; an ultrasound probe that includes an array of electroacoustic elements. The array may be one-dimensional (linear or non-linear, that is, curved) or have a higher dimension, such as two-dimensional. A beamforming system includes a transmit (Tx) beamformer and a receive (Rx) beamformer. The transmit beamformer is provided for activating at least one element of the array to generate a series of transmit beams that insonify a region of interest. The receive beamformer includes a first stage that is provided, for each transmit beam, for determining a plurality of primary receive beams. The receive beamformer also includes a second stage (which may be implemented with the first stage as an integrated beamformer, or separately) provided for computing at least one secondary receive beam as a temporally aligned combination of at least two of the primary receive beams.

0021According to a related method of operation, at least one element of an electroacoustic array of elements is activated to generate a series of transmit beams that insonify a region of interest. For at least one of the transmit beams, a plurality of primary receive beams is then determined. A secondary beamformer stage then computes at least one secondary receive beam as a temporally aligned combination of at least two of the primary receive beams. An image may then be generated as a function of at least one of the secondary receive beams.

0022The secondary beamformer stage may temporally align the at least two primary receive beams that correspond to the ultrasonic return signals from different ones of the transmit beams but that are received on a common receive beam line (a line of ultrasound propagation that passes through a focal point and a corresponding beam origin on the array) and then compute each secondary receive beam as a function of the at least two temporally aligned primary receive beams. One example of a suitable function is the possibly weighted, time-delayed sum of the at least two primary receive beams. One example of temporal alignment comprises computing a time delay of at least one of the primary receive beams relative to another one of the primary receive beams.

0023One example of a method for computing the time delay comprises: simulating conditions corresponding to anticipated transmission conditions of a first prime beam; determining a first instant transmit pulse pressure field (700) for a transmit pulse on the first prime beam for the conditions at a specified time, which corresponds to a specified depth; determining the centers of the transmit pulse along at least one radius direction from a first beam origin (A) of the first prime beam; determining a first transit time (R0) between the first transmit origin (A) and an intersection (B) of the instant transmit pulse pressure field at the specified depth and the first prime beam; determining a second transit time (R1) between a second beam origin (C) and an intersection (D) of a second pulse pressure field for a second primary receive beam; determining a time differential as a predetermined function of the first and second transit times; and in computing at least one of the secondary receive beams, applying to the second primary receive beam a delay computed as a predetermined function of the time differential.

BRIEF DESCRIPTION OF THE DRAWINGS

0024<figref idref="f0001">Figure 1</figref> illustrates the main components of an ultrasound imaging system as found in the prior art.

0025<figref idref="f0002">Figure 2</figref> illustrates the operation of a transmit and primary receive beamformer as in the prior art.

0026<figref idref="f0003">Figure 3</figref> illustrates the distribution of transmit beams in a one-dimensional array and in particular a numbering convention used to illustrate the principles of the invention.

0027<figref idref="f0004 f0005">Figures 4a-4j</figref> illustrate one example of how a secondary series of primary receive beams may be formed from a series of transmit beams.

0028<figref idref="f0006">Figures 5</figref> and <figref idref="f0007">6</figref> are two different representations of how primary receive beams from different transmits may align so as to form the basis for creating a secondary receive beam.

0029<figref idref="f0008">Figure 7</figref> illustrates the calculation of delay values used to form secondary receive beams.

0030<figref idref="f0009">Figure 8</figref> illustrates the components of an ultrasound imaging system that includes two-stage beamforming.

0031<figref idref="f0010">Figures 9a-9c</figref> illustrate two-stage beamforming for a two-dimensional array.

DETAILED DESCRIPTION

0032The multi-stage beamformer and its novel methods of operation according to the invention are described primarily with reference to a one-dimensional array (which may be linear or curved, or a one-dimensional subset of the elements of a higher-dimensional array, etc.), but the principles of the invention are not restricted to this and indeed a 2-D example is illustrated. The invention may be used with any array having any number of elements, as long as there are enough elements to allow for receive element/beam "overlap" as explained below.

0033To better understand the principles of the invention, assume that the invention is to be implemented with an N-element linear array, as is illustrated in <figref idref="f0003">Figure 3</figref>. Transmit (Tx) beams may originate either at the center of an element or between two elements. For the purposes of illustration only, let the Tx beam with an origin between array elements (m-1) and m be numbered (2m-2); the beam with its origin in the center of element m be numbered (2m-1) and the beam with its origin between array elements m and (m+1) be numbered 2m. The beams with origins at the left and right edges of the array (viewed as in <figref idref="f0003">Figure 3</figref>) are then labeled 0 and 2N, respectively. Thus, the transmit beam with its origin substantially between elements 4 and 5 is beam 8; the transmit beam with its origin substantially at the center of element 2 is Tx3; and so on. For an N-element linear array there will therefore in general be (2N+1) possible beam origins, for example 257 beam origins for a 128-element array.

0034Note that "origin" does not mean that the transmitted beam must be steered perpendicular to the array face from that point; thus, many transmit beams may originate from the same point, but may be steered to different angles and depths. Also, that a beam has its origin at a particular point does not necessarily define which array elements are used to generate the beam. For example, Tx5 could be generated by energizing array element 3 alone, or by activating elements 2-4, or just 2 and 4, or elements 1-5, or elements 1, 3 and 5, etc., all of which represent a symmetric combination of elements that has the enter of element 3 as the center of the combination. Even this symmetry is not necessary, however: It is also possible, using known methods, to calculate delays such that a particular Tx beam has an origin not at the center of the set of energized elements, that is, from an asymmetric pattern of elements. For example, Tx8 could also be determined when elements 1-15 are used to form the transmit beam.

0035How a particular transmit beam is generated will in general depend on such considerations as a desired aperture and apodization. This invention does not require any particular aperture, apodization, etc. as long as Tx mean origins can be determined. Another way to think of a beam origin is that it is the location at which the system computes the delays for all elements in the aperture. A beam is the line from the beam origin and through the focal point FP. This is true for both transmit and receive beams, although the direction of sound propagation will of course be reversed (to as opposed to <i>from</i> the focal point). An important calculation in both transmit and receive beamforming is the calculation of delays, and the beam origin is thus the point at which the delay is zero, meaning that other delays are measured relative to it. A delay length may then be defined as the difference between the length between an element to the focal point and the beam origin to the focal point.

0036The numbering convention illustrated in <figref idref="f0003">Figure 3</figref> may also be used to number receive beams. Note that it is not necessary for the same elements or even the same number of elements to be used to receive echo signals as are used to generate the transmit beam. This will again be a question of which aperture, etc., is applied. As one example, elements 1-4 could be energized to form Tx4, but then the signals received at elements 3-5 could be used to form a receive beam with zero or "origin" delay on line 7; thus, the system could transmit on Tx4 but "listen" on Rx7.

0037<figref idref="f0004 f0005">Figures 4a-4j</figref> illustrate how different primary receive beams are formed from two different transmit beams Tx8 and Tx12. Merely by way of illustration, the transmit beams are shown as being generated by four adjacent elements: elements 3-6 for Tx8 and elements 5-8 for Tx12. Similarly, it is assumed merely by way of example that the signals from three adjacent elements are used in the receive mode. For the sake of simplicity, it is also assumed that the transmit and receive lines are at the center of their respective groups of active elements. As explained above, these assumptions reflect just one of several possible choices. The elements used in the transmit are shown with a bold outline. The elements used for reception are indicated by a bracket above them and also an identification of their origin number below. Thus, in <figref idref="f0004">Figure 4a</figref>, Tx8 is generated by elements 3-6, and reception is on elements 5-7, to form Rx11.

0038In the invention, the beam that is formed from the echo signals from a single transmit is referred to as the primary receive beam. The notation PB(i, j) is used here to indicate the primary receive beam that results from a transmit beam whose origin is i and that is received by any chosen group of elements whose receive beam origin is at position j, that is, reception is on line j.

0039In <figref idref="f0004">Figures 4a-4e</figref>, the transmit beam is in each case Tx8, but the elements used to form five parallel receive beams from the same transmit of Tx8 shift from elements 1-3 (Rx3) to form PB(8,3) illustrated in <figref idref="f0004">Figure 4b</figref>, to elements 2-4 (Rx5) to form PB(8,5) illustrated in <figref idref="f0004">Figure 4c</figref>, to elements 3-5 (Rx7) to form PB(8,7) illustrated in <figref idref="f0004">Figure 4d</figref>, to elements 4-6 (Rx9) to form PB(8,9) illustrated in <figref idref="f0004">Figure 4e</figref>, to elements 5-7 (Rx11) to form PB(8,11) illustrated in <figref idref="f0004">Figure 4a</figref>. Thus, the shift in receive mode is one element (two beam numbers) per receive beam.

0040In <figref idref="f0005">Figures 4f-4j</figref>, the transmit beam is Tx12. Relative to what is illustrated in <figref idref="f0004">Figures 4a-4e</figref>, the transmit beam origin has thus been shifted to the right by four beam numbers (two elements). Note that, even for the same target, because of different Tx beams, the echo from this same target will be located at different depths of the prime beams, for example, on lines PB(8,5) and PB(12,5). The second-stage beamformer according to this invention (described in greater detail below) eliminates or at least substantially reduces this difference. In <figref idref="f0005">Figure 4f</figref>, the elements used in the receive mode for Tx12 are elements 2-4 (Rx5) to form PB(12,5); in <figref idref="f0005">Figure 4g</figref>, receive elements 3-5 (Rx7) form PB(12,7); in <figref idref="f0005">Figure 4h</figref>, receive elements 4-6 (Rx9) form PB(12,9); in <figref idref="f0005">Figure 4i</figref>, receive elements 5-7 (Rx11) form PB(12,11); and in <figref idref="f0005">Figure 4j</figref>, receive elements 5-7 (Rx13) form PB(12,13). Thus, each subsequent primary receive beam in <figref idref="f0005">Figures 4f-4j</figref> is created by shifting two beam numbers (one element) to the right relative to the previous primary receive beam.

0041<figref idref="f0006">Figure 5</figref> illustrates one of many possible examples of a portion of an ultrasound scan. In <figref idref="f0006">Figure 5</figref>, the numbered circles represent primary beams. In this illustrated case, eight primary beams are computed (in any known manner) for a Tx beam, then the aperture is "shifted" rightward (as illustrated in <figref idref="f0006">Figure 5</figref>) by four beam numbers, a new transmit beam is formed and new receive beams are computed, then the aperture is shifted four more beam numbers to the right, and so on. <figref idref="f0006">Figure 5</figref> illustrates the activation pattern in receive mode for five transmit beams. Again just by way of example, five transmit beams Tx8, Tx12, Tx16, Tx20, and Tx24 are assumed. In <figref idref="f0006">Figure 5</figref>, the darkened circles represent the primary receive beams that are computed for each respective transmit. Thus, the first transmit beam Tx8 results in primary receive beams PB(8,1), PB(8,3), PB(8,5), PB(8,7), PB(8,9), PB(8,11), PB(8,13), and PB(8,15). In table form, the primary receive beams for each of the five transmit beams is: <tables id="tabl0001" num="0001"><table frame="all"><title><b>Table 1</b></title><tgroup cols="2"><colspec colnum="1" colname="col1" colwidth="24mm" /><colspec colnum="2" colname="col2" colwidth="75mm" /><thead><row><entry align="center"><b>TX beam no.</b></entry><entry align="center"><b>Receive beam numbers</b></entry></row></thead><tbody><row><entry align="center" valign="middle">8</entry><entry valign="middle">1,3,5,7,9,11,13,15</entry></row><row><entry align="center" valign="middle">12</entry><entry valign="middle">  5, 7, 9, 11, 13, 15, 17, 19</entry></row><row><entry align="center" valign="middle">16</entry><entry valign="middle">     9,11,13,15,17,19,21,23</entry></row><row><entry align="center" valign="middle">20</entry><entry valign="middle">        13,15,17,19,21,23,25,27</entry></row><row><entry align="center" valign="middle">24</entry><entry valign="middle">           17,19,21,23,25,27,29,31</entry></row><row><entry align="center" valign="middle">···</entry><entry valign="middle">...</entry></row></tbody></tgroup></table></tables>

0042This is also illustrated in <figref idref="f0007">Figure 6</figref>, which shows and labels the various receive beams for each transmit beam, where the transmit beams are separated in time by an amount Δt, which will be a function of the chosen frame rate. As <figref idref="f0006">Figures 5</figref>, <figref idref="f0007">Figure 6</figref>, and Table 1 show, several of the receive beams include information relating to the same insonified structure, but from different transmit beams. Thus, in this simple example, there are four separate values for primary receive beams 13, 15, 17 and 19. For primary receive beam 13, for example, there are the values PB(8,13), PB(12,13), PB(16,13) and PB(16,13), and for receive beam 19 there are values PB(12,19), PB(16,19), PB(20,19) and PB(24,19).

0043Defined in terms of beam numbers, if array elements corresponding to a span of n beam numbers are used in the receive aperture, and the aperture is shifted m beams on each transmit, then, except near the edges, there will be n/m primary beams for each beam location. The values of n and m will typically be design choices representing a trade-off between resolution and speed. Note that shifting can be defined in terms of either beam numbers or element numbers and that the numbering in either case is simply a matter of choice of convention.

0044The formation of primary beams is referred to here as the first or primary beamforming stage, which is labeled Stage I in <figref idref="f0007">Figure 6</figref>. This invention provides for secondary beamforming that uses the information that conventional beamformers ignore. This secondary beamforming stage provides dynamic beamforming of the primary beams on the same line (as illustrated, having the same beam number).

0045The invention computes (using executable code stored on or in any suitable medium such as in the memory 250) <i>secondary</i> receive beams SB(p) as a function of the primary receive beams on line p. One suitable function is the weighted time-delayed sum of the primary receive beams PB(c,p), summed over all c (transmit beam numbers) for which there is a p (receive beam number) value. Expressed more abstractly: <maths id="math0001"><math display="block"><mi>SB</mi><mfenced><mi mathvariant="normal">p</mi></mfenced><mo mathvariant="normal">=</mo><mi mathvariant="normal">f</mi><mfenced><mi>PB</mi><mfenced><mi mathvariant="normal">c</mi><mo>⁢</mo><mi mathvariant="normal">p</mi></mfenced><mo mathvariant="normal">,</mo><mi mathvariant="normal">t</mi></mfenced><mo mathvariant="normal">=</mo><mi mathvariant="normal">Σ</mi><mo>⁢</mo><mfenced open="[" close="]"><msub><mi mathvariant="normal">w</mi><mi mathvariant="normal">c</mi></msub><mo mathvariant="normal">⋅</mo><mi>PB</mi><mo>⁢</mo><mfenced><mi mathvariant="normal">c</mi><mo mathvariant="normal">,</mo><mi mathvariant="normal">p</mi><mo mathvariant="normal">,</mo><mi mathvariant="normal">t</mi><mo mathvariant="normal">-</mo><msub><mi mathvariant="normal">δ</mi><mi mathvariant="normal">c</mi></msub></mfenced></mfenced></math><img file="EP2053418B1_D0001.tif" /></maths> where the sum is taken over c (all or a selected sub-set) for which there is a p value. Here, PB(c, p, t-δ_c) is the primary receive beam PB(c,p) but that is time-delayed by an amount δ_c and weighted (possibly simply by unity, meaning no differential weighting at all) by the value w_c. In other words, the values of at least some of the primary receive beams received for different transmits but on the same beam are combined. Note that the summation is over the whole length of the beams, that is, over all time t in a chosen interval. Moreover, w_c and δ_c may also be time- (depth-) dependent, and could have different values for different beams.

0046Thus, returning to the example above, SB(13) may be computed as the time-delayed, weighted sum of PB(8,13), PB(12,13), PB(16,13) and PB(16,13). In practice, what this means is that each secondary receive beam includes energy (information) from several transmit beams. Proper choices of aperture, shifting, transmit timing (all of which may be determined from case to case using known operating methods), and application of the invention should therefore lead to a synthesized, secondary receive beam that has less noise and is therefore clearer than any one of the primary receive beams processed in isolation.

0047As another example, assume as before that there are eight parallel Rx beams with a relative spacing of one element (two beam positions), but in this case assume that the shift per transmit is only two beam positions (one element) instead of four (two elements). Table 2 illustrates the correspondence between transmit and receive beam numbers: <tables id="tabl0002" num="0002"><table frame="all"><title><b>Table 2</b></title><tgroup cols="2"><colspec colnum="1" colname="col1" colwidth="23mm" /><colspec colnum="2" colname="col2" colwidth="72mm" /><thead><row><entry align="center" valign="middle"><b>Tx beam no.</b></entry><entry align="center" valign="middle"><b>Receive beam numbers</b></entry></row></thead><tbody><row><entry align="center" valign="middle">8</entry><entry valign="middle">1,3,5,7,9,11,13,15</entry></row><row><entry align="center" valign="middle">10</entry><entry valign="middle"> 3,5,7,9,11,13,15,17</entry></row><row><entry align="center" valign="middle">12</entry><entry valign="middle">  5,7,9,11,13,15,17,19</entry></row><row><entry align="center" valign="middle">14</entry><entry valign="middle">   7,9,11,13,15,17,19,21</entry></row><row><entry align="center" valign="middle">16</entry><entry valign="middle">    9,11,13,15,17,19,21,23</entry></row><row><entry align="center" valign="middle">18</entry><entry valign="middle">     11,13,15,17,19,21,23,25</entry></row><row><entry align="center" valign="middle">20</entry><entry valign="middle">      13,15,17,19,21,23,25,27</entry></row><row><entry align="center" valign="middle">22</entry><entry valign="middle">       15,17,19,21,23,25,27,29</entry></row><row><entry align="center" valign="middle">24</entry><entry valign="middle">        17,19,21,23,25,27,29,31</entry></row><row><entry align="center" valign="middle">26</entry><entry valign="middle">         19, 21, 23, 25, 27, 29, 31, 33</entry></row><row><entry align="center" valign="middle">···</entry><entry valign="middle">...</entry></row></tbody></tgroup></table></tables>

0048It can be seen in Table 2 that from beam number 15 and up, there are eight primary beams for each beam location. To form the final, that is, secondary beam number 15, SB(15), is formed from the primary beams PB(8,15), PB(10,15), PB(12,15), PB(14,15), PB(16,15), PB(18,15),PB(20,15) and PB(22,15). All these eight primary beams will, as before, be time-delayed and weight-summed to align to an assumed beam whose transmit origin and receive origin are at the beam origin of beam 15.

0049In Table 2, observe also that for beam number 9, there are 5 primary beams for each beam location. To form the final beam number SB(9), the primary beams PB(8,9), PB(10,9), PB(12,9), PB(14,9) and PB(16,9) may be time-delayed and weight-summed. In other words, it is possible to have a different number of terms in the sum for different secondary beams; this will typically occur close to the edges of the array. It would also be possible simply not to perform secondary beamforming near the array edges if this is preferred over the additional programming burden. Once the secondary beams are computed, they may be passed on for further processing and display just as primary beams are in conventional beamformers.

0050Known theoretical and experimental methods may be used to determine the weights w_c, assuming that weighing is even implemented. For example, one might choose weights to deemphasize more off-axis primary beams. As another example, one might choose to deemphasize primary beams occurring relatively farther back in time, especially where the image region is itself moving, such as for cardiac imaging, as opposed to more static objects such as abdominal structures. These considerations should also be taken into account in determining the degree of shifting of active array elements - again, the more primary beams are accumulated, the better the secondary beam is likely to be, but the longer this will take, and faster scanning is usually preferred for dynamic regions such as the heart. In practice, the number of elements activated at a time and the degree of shift will be part of a trade-off: The greater the resolution desired, the longer it will take for processing.

0051Part of the computation of the secondary receive beams involves a determination of the dynamic time delays δ_c used to properly align the primary beams temporally in the SB(p) calculations to make sure they really represent the same beam line. In conventional transmit beamformer delay computations, the transmit time is usually computed as a function of geometry -- the travel time of the Tx pulse from a Tx element to the Tx focal point, which is equal to the linear distance from the Tx element to the Tx focal point divided by the speed of sound. This assumes, however, that the acoustic wave pulse front is spherical, with the Tx element at the center; alternatively, these systems assume the Tx element to be a point source.

0052In the second stage of the two-stage beamformer according to the invention, the computations are based not upon signals from or to individual elements, but rather upon the prime transmit beams. In these delay computations, one should therefore also consider the transmit time which is the travel time of the Tx pulse from the Tx origin of the prime transmit beam to the Tx focal point. However, this travel time is not necessarily directly proportional to the linear distance from the Tx origin of the transmit beam to the Tx focal point divided by the speed of sound, because the wave front of the prime transmit beam is no longer necessarily spherical. One special case is where the transmit aperture of all the prime transmit beams is a single probe element; in this special case, geometrical calculations alone will typically suffice to compute the transmit delays. In general cases, however, the transmit apertures of prime transmit beams are a group of probe elements with a delay profile and an apodization; such an aperture cannot be viewed as a point source.

0053One way to compute the Tx delays δ_c of primary beams is through simulation of the instant Tx pulse pressure field of the transmit beams given knowledge of which type of body region (cardiac, abdominal, fetal, obstructing fatty tissue, etc.) the computations are to be used for, that is, what medium the ultrasound will pass though. The main steps of one example of such simulation include:

00541. Set simulation conditions (either computationally or by using physical models) to match the actual Tx conditions of a primary beam. The conditions should include the probe geometry, the Tx pulse, Tx aperture size, Tx apodization, probe element directivity, tissue attenuation, the Tx focal location and others.

00552. Set time (depth) to look, for example, at 50 µsec after the Tx; compute (or measure) the instant Tx pulse pressure field of this set of Tx conditions at the time specified.

00563. Compute the centers of the pulse along the radius directions from the Tx origin. Draw (literally or computationally) a curve along the pulse centers on the 2-D Tx pressure field in 2-D image cases; or draw a surface along the pulse centers on the 3-D Tx pressure field in 3-D image cases.

00574. Find the intersection of the pulse center curve and a transmit beam. Let R₀ be the distance (in time) from the Tx origin to the intersection.

00585. Find the intersection of the pulse center curve and a primary beam. Let R₁ be the distance (in time) from the beam origin of the prime beam to the intersection. The primary beam delay for the second stage BF at this depth for this Tx beam is equal to the difference between the R₀ and R₁.

00596. Compute the delay for all primary beams. This yields a delay profile at one depth.

00607. Repeat steps 2-6 for every depth along the whole field to determine the dynamic Tx delay profiles for primary beams from one prime transmit beam. It maybe possible to reduce the computation by interpolation between computed depths.

00618. Repeat steps 1-7 for all prime transmit beams to determine the Tx delay profiles for all primary beams from all prime transmit beams.

0062As one simplified example, see <figref idref="f0008">Figure 7</figref>. Assume a transmit beam Tx12 and a receive beam originating at beam origin 7, giving PB(12,7). Curve 700 represents the Tx pulse center 49 µs after transmitting. In <figref idref="f0008">Figure 7</figref>, Tx12 beam origin is labeled A and the point of intersection of the curve 700 and Tx12 is labeled B. Let R₀ be the distance from A to B. In this example, R₀=49 µs. The beam origin of PB(12,7) is labeled C and the intersection of curve 700 and PB(12,7) is labeled D. Let R₁ be the distance from C to D, and assume by way of example that R₁=50 µs. Curve 701 illustrates an assumed pulse center 50 µsec from C. Now if the system were to use a beam having an origin at location 7 to see a target located at D, it would fire a transmit beam from C, the pulse would travel 50 µsec to D and will then reflect back to C for reception. This would correspond to a primary receive beam PB(7,7), whose round-trip time would thus be 50µs+50µs =100µs. Using the two-stage beamformer, however, the system could fire Tx12, whose pulse originates at A, travels 49 µs and reaches the target at D. The echo from D will then take 50 µs from D to C for reception of the beam at location 7. In this case, the total round-trip time will be 49µs+50µs =99µs. In the second stage of the beamformer, the system may therefore apply a 1 µs delay to PB(12, 7) at that depth to correct the Tx beam offset to generate a beam that appears like PB(7,7)..

0063The two-stage beamformer according to the invention can be programmed flexibly in many ways to display various characteristics. In general, if the array and transmit controller are chosen and programmed so as to be able to form Np parallel beams from one transmit, and if Nt transmits are used to form an image, and if k primary beams are used to form one final, secondary beam, then the total number of final, secondary beams will be Nt*Np/k, since each transmit will give Np/k secondary beams.

0064If the frame rate is more important in a given application, one might, for example, configure the system such that k=2 to have Np/2 final beams per transmit. Thus, if Np=16, there will be eight final, secondary beams from one transmit, with few parallel artifacts due the second stage Tx beamformer of two primary beams. It would even be possible to set k=1 to form 16 beams per transmit (assuming by way of example that is the chosen aperture) to maximize frame rate, although at the cost of some parallel artifacts.

0065If image quality is of great importance, with Np=16 one could set k=16 to have a single final beam per transmit. Each secondary beam will then be formed dynamically from 16 primary beams; therefore, the final beams will have good Tx focusing along the depth.

0066For deep penetration, again with Np=16 as an example, one could also set k=16, yielding a single final, secondary beam per transmit. Each final beam may then be dynamically delay-summed from 16 primary beams, such that each final beam will comprise the acoustic power of 16 prime transmit beams. This should increase penetration roughly 12 dB for a 2.5MHz probe, which translates into a significant increase of depth of around 4.8 cm.

0067<figref idref="f0009">Figure 8</figref>, in which some components of <figref idref="f0001">Figure 1</figref> have been removed for clarity, illustrates how the Rx controller 840 includes both the first-stage and the second-stage beamformers 842, 844. Two common components of the Tx controller 820 are also illustrated, namely aperture and apodization controllers 822, 824.

0068It is not necessary to implement the two-stage beamformer of the invention as two hardware stages. It may, for example, prove practical to combine the second stage dynamic Tx delays into the first stage dynamic Rx delays. However, the first stage dynamic Rx delays can be computed as in the prior art as a function of geometry; this may be implemented using hardware calculators to reduce the need to transfer more control data. Computation of the delays δ_c for the second stage, however, is not as simple, as discussed above. These values may instead be pre-computed, pre-stored and loaded during live imaging from a table 860 stored either in the main system memory or in some other storage device; moreover, the values may be computed in real time, or off-line, or though a combination of both.

0069The computation of the SB(p) values may be accomplished in the second stage beamformer using a memory/accumulator, with a memory unit (either dedicated, for speed, or a portion of other memory such as the system memory) whose size is at least as big as needed to contain all the values defining the required number of beams, and at least equal to the maximum number of parallel beams per transmit. Once the SB(p) values are computed in the second receive beamformer stage, then these values may be passed on for subsequent conventional processing to form the basis of the display of the region of interest.

0070The principles of the invention may also be applied to higher dimensional arrays. <figref idref="f0010 f0011">Figure 9</figref>, for example, illustrates a 8x16 element portion of a 2-D array 900, that is, with i = 8 "rows" and j = 16 "columns" of elements. Let E(i = a,b; j = x,y) represent the portion of the array comprising elements in rows a to b and columns x to y.

0071As in the one-dimensional case, beam <i>numbers</i> may also be defined to label the different possible origins of the transmit beams in the 2-D array. One possible numbering scheme for transmit beam origins for a 3-by-3 section (the upper left section, viewed as in <figref idref="f0010">Figure 9a</figref>) of the 2-D array is illustrated in <figref idref="f0011">Figure 9b</figref>, in which circles indicate beam origins and adjacent two-dimensional labels in parentheses indicate the beam numbers. Thus, beam number (3,3) (and thus the corresponding beam Tx(3,3) originating at the indicated point) is in the center of element E(2,2), beam number (4,5) is at the middle of the edge between E(2,3) and E(3,3), and so on. Other numbering schemes might, for example, not include any origin locations other than corners and centers of elements; the numbering scheme is simply a convention used here to illustrate the principles of the invention, although some numbering scheme will normally be adopted to make programming the software controlling the transmit and receive beamformers easier.

0072As mentioned above, partial or entire rows or columns of elements can be activated in a 2-D array such that they will in effect form 1-D arrays. This is not necessary either in general for diagnostic ultrasound imaging or in particular for implementing this invention; rather, the principles of the invention may be practiced in higher dimensions than just 1-D. By way of example, assume that a first transmit beam Tx(4,12) is formed from the shaded 4x4 subset of elements E(i =1,4; j = 5,8), with a beam origin corresponding to beam number (4,12), which is shown as a small open circle. Assume by way of simplicity that a scan proceeds with a receive aperture chosen such that the first of the primary receive beams (indicated by solid dots) corresponds to the beam origin with beam number (4,4), then the aperture is shifted one element (two beam lines) in the j direction (to the right, viewed as in <figref idref="f0010">Figures 9a</figref> and <figref idref="f0011">9b</figref>), a new primary beam is computed for beam number (4,6), then the receive aperture is again shifted one element to the right, another primary beam is computed for beam number (4,8), then the receive aperture is again shifted one element to the right, and another primary beam is computed for beam number (4,10).

0073PB[(ti, tj); (ra, rb)] may then be used to represent the primary receive beam that results from a transmit beam having the origin indicated by beam number (ti, tj) as received on the beam line with the number (numbered with the same scheme) (ra, rb). Thus, in the example given in the previous paragraph, for transmit beam Tx(4,12), primary beams PB[(4,12); (4,4)], PB[(4,12); (4,6)], PB[(4,12); (4,8)], and PB[(4,12); (4,10)] are computed in any conventional manner.

0074Now the transmit aperture may be shifted, for example, two elements (four beam numbers) to generate a transmit beam Tx(4,16) from the origin at beam number (4,16), indicated in <figref idref="f0010">Figure 9c</figref> by the open circle at the center of the shaded aperture region elements E(i =1,4; j = 7,10) Using the same pattern for receive beams as before (indicated again by solid dots), four primary receive beams can be computed for PB[(4,16); (4,8)], PB[(4,16); (4,10)], PB[(4,16); (4,12)], and PB[(4,16); (4,14)].

0075In this simple example, primary receive beams PB[(4,12); (4,8)] (from the first transmit) and PB[(4,16); (4,8)] (from the second transmit) both are primary receive beams on the same beam line, namely (4,8). The same applies to PB[(4,12); (4,10)] and PB[(4,16); (4,10)] on beam line (4,10). Each of these pairs may then be combined into a single secondary receive beam SB(4,10) as a weighted (again, even by unity), time-delayed sum. These "overlapping" receive beams are marked in <figref idref="f0010">Figure 9c</figref> by being surrounded by a dashed rectangle.

0076It is not necessary for the shifting of either transmit or receive apertures to be in either the azimuthal or elevational direction alone. Even if shifting of one or the other, or both, contains a "diagonal" component, as long as more than one primary receive beam is formed for the same beam line, then these beams may be combined into a secondary beam. The degree of shifting on successive transmits and receives will in general determine the degree of "overlap" and thus how many terms (primary receive beams) are combined to form a single, corresponding secondary receive beam.

0077This invention provides a multi-stage beamformer (with features relating to both transmit and receive beamforming) that addresses several of the weaknesses of the prior art: It can candle a large number of parallel receive beams from a single transmit. In the second stage, it enables dynamic transmit focusing, dynamic transmit aperture size control and dynamic transmit apodization. Depending on some of the choices of parameters described above, such as aperture, degree of overlap, etc., this beamformer may also allow for an increased the frame rate. It also may improve the transmit beamformer quality and more fully use the acoustic available power.