systems and methods for time-delay inverted beamforming are provided. One method includes sampling, substantially uniformly in time, a set of continuous-time element signals to form a set of sampled element signals and mapping, for each sampled element signal, one or more samples substantially equally spaced in time to one or more samples non-equally spaced in time corresponding to the respective contribution of each element signal to a continuous-time beamsum signal. The method also includes forming a beamsum signal sampled substantially uniformly in time from the non-uniformly spaced, mapped samples corresponding to each sampled element signal.
|
1. A beamforming method, for use in an ultrasound system, comprising:
sampling using a beamformer processor of the ultrasound system, uniformly in time, a set of continuous-time element signals to form a set of sampled element signals;
mapping using the beamformer processor of the ultrasound system, for each sampled element signal, samples of element signals equally spaced in time to output samples non-equally spaced in time corresponding to the respective contribution of each sampled element signal to a continuous-time beamsum signal to generate non-equally spaced, mapped samples based on beam direction and element position;
forming using the beamformer processor of the ultrasound system, a beamsum signal sampled uniformly in time from the non-equally spaced, mapped samples corresponding to each sampled element signal; and
generating an image based on the beamsum signal.
16. A non-transitory computer readable medium encoding one or more executable routines, which, when executed by a processor, cause the processor to perform acts of a beamformer method for use in an ultrasound system, comprising:
sampling using a beamformer processor of the ultrasound system, uniformly in time, a set of continuous-time element signals to form a set of sampled element signals;
mapping using the beamformer processor of the ultrasound system, for each sampled element signal, samples of element signals equally spaced in time to output samples non-equally spaced in time corresponding to the respective contribution of each sampled element signal to a continuous-time beamsum signal to generate non-equally spaced, mapped samples based on beam direction and element position;
forming using the beamformer processor of the ultrasound system, a beamsum signal sampled uniformly in time from the non-equally spaced, mapped samples corresponding to each sampled element signal; and
generating an image based on the beamsum signal.
5. An ultrasound system, comprising:
a transducer array comprising one or more transducers;
a transmitter configured to transmit electrical signals to the transducer array;
a receiver configured to receive electrical signals generated by the transducer array when the one or more transducers detect an ultrasonic echo;
transmitter/receiver switching circuitry coupled to the transducer array and configured to switch the transmitter and the receiver;
an analog-to-digital converter configured to receive analog data from the receiver corresponding to the ultrasonic echo and to convert the analog data to a digital element signal; and
a beamformer processor configured to receive and process the digital element signal to form a beamsum signal, wherein the beamformer processor is configured to sample, uniformly in time, a set of continuous-time element signals to form a set of sampled element signals, to map, for each sampled element signal, samples of element signals equally spaced in time to output samples non-equally spaced in time corresponding to the respective contribution of each sampled element signal to a continuous-time beamsum signal to generate non-equally spaced, mapped samples based on beam direction and element position, to form the beamsum signal sampled uniformly in time from the non-equally spaced, mapped samples corresponding to each sampled element signal, and to generate an image based on the beamsum signal.
2. The method of
3. The method of
4. The method of
6. The system of
7. The system of
8. The system of
9. The system of
10. The system of
11. The system of
12. The system of
13. The system of
14. The system of
15. The system of
17. The non-transitory computer readable medium of
18. The non-transitory computer readable medium of
19. The non-transitory computer readable medium of
20. The non-transitory computer readable medium of
|
The subject matter disclosed herein relates generally to imaging systems and methods and, more particularly, to systems and methods for performing inverted beamforming in ultrasonic imaging systems.
Medical diagnostic ultrasound is an imaging modality that employs ultrasound waves to probe the acoustic properties of the body of a patient and produce a corresponding image. Generation of sound wave pulses and detection of returning echoes is typically accomplished via a plurality of transducers located in the probe. Such transducers typically include electromechanical elements capable of converting electrical energy into mechanical energy for transmission and mechanical energy back into electrical energy for receiving purposes. Some ultrasound probes include up to thousands of transducers arranged as linear arrays or a two-dimensional matrix of elements.
A beamformer is typically associated with the ultrasound probe to combine the signals from each transducer element so as to enhance the response to received acoustic energy from a desired focus depth and direction. In certain ultrasound systems, both the inputs and output of the beamformer, the element signals, and the beamsum signal, are sampled uniformly in time. In certain systems, the beamformer focus function is implemented by calculating for each beamsum sample a corresponding time in each element signal that accounts for the propagation time for acoustic energy emitted from the transducer, reflected from scatterers at the desired focus point and returning to that element. In general, this propagation time used in the focus function will not correspond to one of uniformly spaced element samples that are available. In some traditional beamformers, the value at the focus time is obtained by interpolating the available uniformly spaced element samples to approximate the value that would have been obtained had the element signal been sampled at the focus time.
This interpolation can be performed in dedicated hardware, such as a custom application-specific integrated circuit (ASIC), or in a computer processor, such as a general purpose processor (CPU, or central processing unit) or a graphical processing unit (GPU). The first implementation may be called a hardware beamformer, and the second implementation may be called a software beamformer.
In some applications, it is useful to calculate more than one beamsum signal from the same set of element signals, i.e., from the data acquired during a single transmit event. This increases the maximum image display rate by reducing the number of transmit firings required to form an image over a desired two-dimensional scan area or three-dimensional volume. In a hardware beamformer, the element signal samples are typically processed sequentially as they are produced and then discarded. Forming multiple beamsums typically requires duplicating the interpolation hardware for each beamsum. This tends to be costly in terms of both power and expense, especially in volumetric imaging when many (e.g., approximately 16-64) receive beams are desired for each transmit event.
In a software beamformer, the element signal samples for one or more transmit firings are typically stored in memory temporarily, and the interpolation software accesses the required samples from the memory. In a software beamformer, the number of receive beams that can be calculated is limited primarily by the speed at which the processor can perform the calculations. In many applications, the software beamformer has significant advantages in power consumption and expense compared with the hardware beamformer.
The calculation speed is determined by both the raw speed of the processor, i.e., the rate at which math operations can be performed in the processor, and by the memory bandwidth, the rate at which data can be read and written to memory from the processor. The memory bandwidth may be the limiting factor in applications, such as software beamforming, in which only a small number of math operations needs to be performed on each of a very large number of samples. This is especially true when the memory in a specific application cannot be addressed sequentially, since the interface between processor and memory is typically highly optimized for sequential memory access. Unfortunately, conventional approaches to calculating multiple beamsums in a software beamformer require repeated, non-sequential access to data stored in memory, thus limiting the system efficiency.
In one embodiment, a method includes sampling, substantially uniformly in time, a set of continuous-time element signals to form a set of sampled element signals and mapping, for each sampled element signal, one or more samples substantially equally spaced in time to one or more samples non-equally spaced in time corresponding to the respective contribution of each element signal to a continuous-time beamsum signal. The method also includes forming a beamsum signal sampled substantially uniformly in time from the non-uniformly spaced, mapped samples corresponding to each sampled element signal.
In another embodiment, an ultrasound system includes a transducer array having one or more transducers, a transmitter adapted to transmit electrical signals to the transducer array, a receiver adapted to receive electrical signals generated by the transducer array when the one or more transducers detect an ultrasonic echo, transmitter/receiver switching circuitry coupled to the transducer array and adapted to switch the transmitter and the receiver, and an analog-to-digital converter adapted to receive analog data from the receiver corresponding to the ultrasonic echo and to convert the analog data to a digital element signal. The system also includes a beamformer adapted to receive and process the digital element signal to form a beamsum signal. The beamformer is adapted to sample, substantially uniformly in time, a set of continuous-time element signals to form a set of sampled element signals, to map, for each sampled element signal, one or more samples substantially equally spaced in time to one or more samples non-equally spaced in time corresponding to the respective contribution of each element signal to a continuous-time beamsum signal, and to form a beamsum signal sampled substantially uniformly in time from the non-uniformly spaced, mapped samples corresponding to each sampled element signal.
In another embodiment, a computer-readable medium encoding one or more executable routines, which, when executed by a processor, cause the processor to perform acts that include sampling, substantially uniformly in time, a set of continuous-time element signals to form a set of sampled element signals, mapping, for each sampled element signal, one or more samples substantially equally spaced in time to one or more samples non-equally spaced in time corresponding to the respective contribution of each element signal to a continuous-time beamsum signal, and forming a beamsum signal sampled substantially uniformly in time from the non-uniformly spaced, mapped samples corresponding to each sampled element signal.
These and other features, aspects, and advantages of the present invention will become better understood when the following detailed description is read with reference to the accompanying drawings in which like characters represent like parts throughout the drawings, wherein:
The time-delay function in beamforming is typically described and implemented as a mapping from uniformly spaced output sample times in the beamsum signal to non-uniformly spaced sample times in the element signals. Since the input samples that are available are uniformly spaced in time, the desired input values are obtained by interpolation. However, as described in more detail below, provided herein are systems and methods directed toward a novel approach to beamforming that inverts this mapping. More specifically, in the provided embodiments, the time-delay function is rewritten to map the uniformly spaced input sample times, which are available in the element signals, to non-uniformly spaced output sample times in the beamsum signal. In provided embodiments, the desired uniformly spaced output sample values are obtained by interpolation from a non-uniform to a uniform sampling time grid. In presently disclosed embodiments of this “inverted beamforming” method, it is only necessary to access each element sample once, regardless of the number of beamsum signals that are calculated from a given set of element samples. The foregoing feature may be advantageous, for example, in software beamforming, particularly when large numbers of beamsum signals are calculated using the same set of element samples, where the memory bandwidth may be a limiting factor. These and other features of the presently disclosed inverted beamforming methods and systems are described in more detail below.
Turning now to the drawings,
In the illustrated embodiment, a piezoelectric transducer array 18 is provided with a front face or acoustic lens structure 19 that is adapted to contact a subject 16, such that an ultrasonic scan may be performed to analyze internal features of the subject 16. Typically, the same transducer elements both generate and receive ultrasound energy in a pulse-echo mode, although different elements on the transducer may be used for these functions in some embodiments.
The transducer array 18 is connected via transmitter/receiver switching circuitry 22 to a transmitter 24 and a receiver 26. The transmitter/receiver switching circuitry 22 switches the electrical connections between the transducer array 18 and the transmitter 24 and receiver 26. In operation, the transmitter 24 is connected to the transducer array 18 when ultrasound energy has to be transmitted into the body of the subject 16, and the receiver circuitry 26 is connected to the transducer array 18 when the transducer array 18 receives the echo signals from the tissue layers of the subject 16.
That is, the illustrated transducer array 18 comprises a two-way transducer. When ultrasound waves are transmitted into a subject 16, the ultrasound waves are backscattered off the tissue and blood within the subject 16. The transducer elements of the transducer array 18 receive the backscattered waves at different times, depending on the distance into the tissue they return from, and the angle with respect to the surface of the transducer array 18 at which they return. The transducer elements are responsive to the backscattered waves and convert the ultrasound energy from the backscattered waves into electrical signals.
The electrical signals received by the transducer array 18 are routed through the transmitter/receiver switching circuitry 22 to the receiver 26. The receiver 26 amplifies the received signals after proper gain compensation, and an analog-to-digital converter (ADC) 20 converts these received analog signals from each transducer array element to digitized signals sampled uniformly in time, which are stored temporarily in system memory 33. The digitized signals correspond to the backscattered waves received by each transducer element at various times. After digitization, the signals still preserve the amplitude and phase information of the backscattered waves.
The system processor 32 retrieves the element signals from the system memory 33 and transfers them to beamformer memory 23. In some implementations, the system processor 32 may modify the element signals, such as converting them to baseband signals or compressing the signals. In other implementations, these additional processing steps may be perform in dedicated hardware before the signals are stored in the system memory 33.
A beamformer processor 21 reads the element signals from the beamformer memory 23, performs the beamforming steps as described in more detail below, produces the beamsummed signal, and writes the beamsummed signal to the beamformer memory 23. The system processor 32 reads the beamsummed signal from the beamformer memory 23 and processes the beamsummed signals from various transmit events to form an image, which may be displayed on a display 42 if desired in the given embodiment.
In the illustrated embodiment, the beamforming processing is performed in a processor (i.e., beamformer processor 21) separate from the system processor 32. However, it should be noted that in other embodiments, the beamforming processing may be performed in the system processor 32, or in any other suitable circuitry in the system 10.
In the illustrated embodiment, the system 10 also includes a user interface 44 that is in communicative coupling with the system processor 32 and the display unit 42 to enable a user to communicate with the processor 32, for example, to input one or more desired imaging or display parameters. In certain embodiments, the system processor 32 may also be coupled to a remote connectivity module 46 having a web server 48 and a remote connectivity interface 50 for coupling the ultrasound system to a network, via link 52. System processor 32 may be further coupled to an image database 54 to receive ultrasound image data. In turn, the image database 54 may be in communicative coupling with imaging workstation 56.
icts=2|R(i)|, (1)
where i is the sample number, an integer value, c is the sound propagation speed, and ts is the sampling time interval; the factor of two accounts for the two-way propagation of sound from the origin to the focus point R(i) and back to the origin. For an element at xe, these equally spaced output samples correspond to a set of unequally spaced sampling times (1/c)[|R(i)|+|R(i)−xe|], which is the time required for sound to propagate from the origin to the focus point R(i) and back to the element at xe. Since the available element samples are equally spaced in time, the desired element sample values are typically obtained by interpolation.
It should be noted that for simplicity, the beamsum is described herein as corresponding to a set of focus points which lie along a (straight) line oriented in a particular direction. However, it will be appreciated that more generally, with the proper choice of focus function, a beamformer may produce a beamsum which corresponds to a set of focus points which lie along a curve, and presently disclosed embodiments are not limited to those in which the beamsum corresponds to focus points which lie along a straight line.
The methods and techniques provided herein invert this procedure by instead mapping the equally spaced samples of a particular element signal to a set of unequally spaced output samples and performing the interpolation on the output:
icts=|R|+|R−xe|. (2)
Equation (2) can be solved for R≡|R|:
With the definitions:
β≡2R/cts (4)
g(|xe|)≡|xe|/cts
h(xe·{circumflex over (R)})≡xe·{circumflex over (R)}/cts={circumflex over (x)}e·{circumflex over (R)}g,
then
where β(i,g,h) is the non-integer sample number in the output corresponding to the integer sample number i in the signal for element xe and beam direction
The singularity in β at i=h is unphysical. From Equation (2), i c ts→|xe|/c ts=g for R→0, so that the minimum relevant value of i is g, which is non-negative. From Equation (4), we see that when h is non-negative, it is never larger than g. Thus, the denominator in Equation (5) will never vanish in practice. More precisely, the minimum value of i is ceil(g), the smallest integer not smaller than g. The special case xe=0 must also be avoided, but for this case, the desired element sampling times are just the uniformly spaced beamsum sample times.
These equations guide the presently disclosed embodiments of the inversion of the typical focusing equation for imposing beamforming time-delays. This inverted beamforming may be better understood by considering the schematics 62 and 64 in
In
In the inverted beamforming approach illustrated in the schematic 64 of
An embodiment of a method 100 suitable for performing this inverted beamforming approach is shown in
Embodiments of the presently disclosed inverted beamforming systems and methods provided herein may offer a variety of advantages over conventional beamforming approaches. For example, in conventional beamforming, whether performed in hardware or software, each element signal must be interpolated to implement the beamforming time-delay operation. However, in the presently disclosed inverted beamforming approach, the beamsum signal, not the element signals, is interpolated. In many ultrasound imagers, the number of beamsum signals calculated for a given transmit firing is much smaller than the number of element signals. Accordingly, in certain embodiments, by utilizing the disclosed inverted beamforming approach, fewer computations may be necessary, thus providing for reductions in power and monetary cost.
In addition, in conventional beamforming, the sampling frequency of the digital element signals is typically much greater than the Nyquist sampling frequency in order to reduce the cost and complexity of the interpolation hardware, or the computational power needed when the interpolation is performed in software. The Nyquist sampling frequency is the theoretical minimum sampling rate need to reconstruct a signal at any arbitrary time from its sampled values. The cost and complexity of an ADC typically increases with sampling rate, but generally, using a relatively high sampling rate minimizes the total system cost of a traditional beamformer. However, in the presently disclosed inverted beamforming approach, only the beamsum signal is interpolated, so the sampling rate of the large number of element signals may be reduced, and the complexity of the interpolating software of the small number of beamsum signals may be increased to achieve substantially the same performance as a traditional beamformer but at a lower system cost.
The implementation of a conventional software beamformer, as depicted in
However, in the inverted beamforming approach illustrated in the schematic 126 of
As described in detail above, the inverted beamforming approach produces, as an intermediate step, non-uniformly spaced contributions to the beamsum from each element signal. Typically, however, it is more convenient to have uniformly spaced beamsum samples for subsequent processing steps such as filtering and scan conversion. The techniques discussed above describe one method (i.e., interpolation) for converting the non-uniformly spaced samples to uniformly spaced samples. However, it should be noted that presently disclosed embodiments are not limited to the use of interpolation. Indeed, any of a variety of suitable methods for converting the non-uniformly spaced samples to uniformly spaced samples may be utilized. For example, in one embodiment described in more detail below, an oversampling method may be used instead of interpolation. However, again, any suitable method may be utilized in other embodiments.
In general, in one suitable oversampling embodiment, an output buffer for a beamsum signal with an increased or high sampling rate with respect to the input signal sampling rate is provided. In one embodiment, the sampling rate may be approximately twenty times the sampling rate that would normally be used for the element and beamsum signals. For example, when imaging with a center frequency of 5 MHz, a sampling rate of 50 MHz might traditionally be used for the element and beamsum signals, so the output buffer would use a sampling rate of 1 GHz. In one embodiment, the output buffer is initialized to zero and then each non-uniformly spaced element sample is added to the value in the buffer at the nearest position corresponding to the value β taking into account the higher sampling rate in the beamsum compared with the element signal. For example, if the oversampling factor p is twenty, then the sample would be added to the value in position Round (p β), where the Round( ) operation returns the nearest integer to its argument. Once all the time-delayed samples for all the elements for a given beam direction are added to the output buffer, the buffer is lowpass filtered and decimated to the desired output sampling rate.
In some embodiments, by utilizing the oversampling method shown in
This written description uses examples to disclose the invention, including the best mode, and also to enable any person skilled in the art to practice the invention, including making and using any devices or systems and performing any incorporated methods. The patentable scope of the invention is defined by the claims, and may include other examples that occur to those skilled in the art. Such other examples are intended to be within the scope of the claims if they have structural elements that do not differ from the literal language of the claims, or if they include equivalent structural elements with insubstantial differences from the literal languages of the claims.
Patent | Priority | Assignee | Title |
Patent | Priority | Assignee | Title |
5197037, | Oct 31 1991 | Koninklijke Philips Electronics N V | Method and apparatus for the simultaneous performance of the beam formation and scan conversion in a phased array system |
5844139, | Dec 30 1996 | General Electric Company | Method and apparatus for providing dynamically variable time delays for ultrasound beamformer |
6052427, | Oct 05 1998 | Arch Development Corp. | Multidimensional interpolation |
6747584, | Oct 10 2001 | National Research Council of Canada | Aggregate beamformer for use in a directional receiving array |
7737892, | Mar 21 2003 | Qinetiq Limited | Time delay beamformer and method of time delay beamforming |
7740583, | Jun 30 2004 | General Electric Company | Time delay estimation method and system for use in ultrasound imaging |
8241216, | Jun 06 2008 | Siemens Medical Solutions USA, Inc | Coherent image formation for dynamic transmit beamformation |
8241217, | Jun 29 1995 | TeraTech Corporation | Portable ultrasound imaging data |
20050228284, | |||
CN101571518, | |||
CN1714752, | |||
EP539700, | |||
EP1808122, |
Executed on | Assignor | Assignee | Conveyance | Frame | Reel | Doc |
Jan 16 2013 | General Electric Company | (assignment on the face of the patent) | / | |||
Jan 16 2013 | RIGBY, KENNETH WAYNE | General Electric Company | ASSIGNMENT OF ASSIGNORS INTEREST SEE DOCUMENT FOR DETAILS | 029643 | /0927 |
Date | Maintenance Fee Events |
Nov 21 2019 | M1551: Payment of Maintenance Fee, 4th Year, Large Entity. |
Nov 22 2023 | M1552: Payment of Maintenance Fee, 8th Year, Large Entity. |
Date | Maintenance Schedule |
Jun 28 2019 | 4 years fee payment window open |
Dec 28 2019 | 6 months grace period start (w surcharge) |
Jun 28 2020 | patent expiry (for year 4) |
Jun 28 2022 | 2 years to revive unintentionally abandoned end. (for year 4) |
Jun 28 2023 | 8 years fee payment window open |
Dec 28 2023 | 6 months grace period start (w surcharge) |
Jun 28 2024 | patent expiry (for year 8) |
Jun 28 2026 | 2 years to revive unintentionally abandoned end. (for year 8) |
Jun 28 2027 | 12 years fee payment window open |
Dec 28 2027 | 6 months grace period start (w surcharge) |
Jun 28 2028 | patent expiry (for year 12) |
Jun 28 2030 | 2 years to revive unintentionally abandoned end. (for year 12) |