The present invention relates to multi-line artifact mitigation in ultrasound imaging during relative motion between ultrasound probe and an object being imaged.
Ultrasonic imaging systems are known, e.g. for producing real-time images of internal portions of the human body. An array of transducers is controlled to produce a transmit (TX) beam which propagates in a predetermined direction from the array. Reflected pressure pulses are received by the receive transducers which may be a sub-set or super-set of the transmit transducers. The reflected pressure pulses may be focused in a receive (RX) beam. Round-trip (RT) beams are, to a first approximation, the multiplication of the TX and RX beams. The collection of transducer compensating delays and signal summing circuitry for forming the transmit, and receive and round trip beams is referred to as a beamformer and is described, for example, in U.S. Pat. No. 4,140,022, which is incorporated herein by reference. The beamformer outputs a radio frequency (RF) signal representing amplitudes of received pressure pulses. A scan converter is disclosed for example in U.S. Pat. Nos. 4,468,747 and 4,471,449, the entire contents of which are incorporated herein by reference, for converting the RF signals output by the beamformer to information in X-Y coordinates used for display of an image on a monitor screen.
The number of lines of data sent to the scan converter is determined by the beam widths of the receive beams. Too few lines results in spatial aliasing which is exhibited as scintillation artifacts in the single lateral dimension for 2D scanning or in elevation and azimuth dimensions for 3D scanning. Scintillation artifacts result when the transducer array is shifted relative to the object. Detection and compression are non-linear operations which increase the lateral spatial frequency band widths. Accordingly, even if the beams going into the detector are not spatially aliased, it is possible that they exhibit spatial aliasing at the output of the detector.
U.S. Pat. Nos. 5,318,033 and 5,390,674 disclose a method for overcoming the problem of spatial aliasing by laterally upsampling using an interpolation filter for filtering the RF signals output by the beamformer. In this method, the TX, RX, and RT beams are collocated and the upsampling is performed on the RT beams.
According to receive multi-line imaging techniques more than one RX beam is acquired for each TX beam. Accordingly, there are more RT beams available for the detector, one for each TX/RX beam pair. The RX beams are displaced from the TX beam so that the RX beams straddle the TX beam.
In multi-line imaging, the location of each RT beam is displaced from both the constituent TX and RX beams. The RT beams are asymmetrical and the amplitudes of the RT beams are less than if the TX and RX beams are collocated. The displacements, asymmetries and amplitude losses associated with the RT beams cause jail-bar artifacts (alternating groupings or stripes aligned in the axial scan direction). Jail-bar artifacts are different from scintillation artifacts in that jail-bar artifacts occur even when there is no motion. TX focus is fixed and RX focus is dynamic. Therefore, the displacements, asymmetries, and amplitudes associated with the RT beams, and thus the jail-bar artifacts, are depth dependent. Jail-bar artifacts may be reduced by broadening or flattening the TX beams as described in U.S. Pat. Nos. 4,644,795 and 6,585,648 or by lateral filtering following the detector or compressor. However, these approaches tend to reduce lateral resolution.
Multi-line Artifact Mitigation (MAM), also referred to as Parallel Artifact Mitigation (PAM), is a technique for eliminating or at least reducing jail-bar artifacts while preserving spatial resolution and is described in U.S. Pat. No. 5,318,033. Various schemes exist, but a common element in all of the MAM schemes is that a filter is applied to received multi-line data prior to detection (the standard form of RF interpolation for scintillation reduction operates on collocated TX, RX, and RT data or, in other words the case of no multi-line). The two or more RT beams that are filtered typically arise from either different RX beam locations arising from a common TX beam or RX beams at the same location arising from different TX beams, i.e. common TX and common RX, respectively. MAM improves mutual similarities between all synthesized RT beams.
However, MAM assumes that the tissue is stationary with respect to the ultrasound probe. Excessive motion reintroduces jail-bar artifacts because the phases of RF data used in MAM varies from that assumed resulting variable amounts of destructive interference.
Excessive motion is defined as motion causing displacements of approximately ⅕ wavelength of the ultrasound signals during the period between successive transmit events used to synthesize the data. For 2D scanning, this period is typically about 200 μsec. At 3 MHz and a wavelength of 0.5 mm, the excessive motion is reached at an axial velocity of approximately 25 cm/sec. For 3D scanning there is typically a fast scan and slow scan dimension. The period between transmit events in the fast scan dimension is about the same as for 2D scanning. However, the period between transmit events in the slow scan dimension may be as high as 25 times larger, reducing the excessive motion axial velocity threshold to approximately 1 cm/sec.
In the example of
An object of the present invention is to eliminate or at least reduce jail-bar artifacts in ultrasound imaging caused by relative motion between the ultrasound probe and the subject being imaged.
The object of the present invention is met by a method of ultrasound imaging including scanning a patient or object using an ultrasound probe, monitoring for excessive relative motion between the object being imaged and the ultrasound probe, and implementing a jail-bar reduction process when excessive relative motion is detected.
Methods to detect the motion include image analysis, Doppler analysis, jail-bar detection, or use of a motion detector within the ultrasound probe. The relative motion may be caused by motion of the ultrasound probe, motion of object, or part of the object (i.e., beating heart or heart values), or a combination thereof. Excessive motion is generally defined as motion which causes jail-bar artifacts. The threshold of excessive motion will be different for different scanning modes. While the excessive speed is lower, and therefore more easily surpassed, in 3D imaging, the techniques of the present invention may be applied to 2D imaging of rapidly moving structures, such as a heart value.
When imaging a 3D volume with a 2D array which produces a single transmit beam per transmit firing, scanning may be effected similarly to a TV raster scan in which the transmit beam is scanned across one row rapidly (fast scan dimension). Once one row is completely scanned, the next row down is scanned, the vertical dimension being a slow scan dimension. The technique according to the present invention is particularly suitable for eliminating jail-bar artifacts in the “slow scan” dimension of the above-described TV raster scanning method.
Methods which may be used to reduce jail-bars includes reducing or turning off the MAM and implementing a further jail-bar reduction technique such as for example, spatial filtering, temporal filtering, dropping multi-line order, beam broadening, and normalization of average A-line amplitudes between lines.
Alternatively, the MAM may be maintained and the RF data may be pre-aligned to counter the misalignment caused by the relative motion between the ultrasound probe and the object being imaged.
Other objects and features of the present invention will become apparent from the following detailed description considered in conjunction with the accompanying drawings. It is to be understood, however, that the drawings are designed solely for purposes of illustration and not as a definition of the limits of the invention, for which reference should be made to the appended claims. It should be further understood that the drawings are not necessarily drawn to scale and that, unless otherwise indicated, they are merely intended to conceptually illustrate the structures and procedures described herein.
In the drawings:
The determination of relative motion may be accomplished by image analysis, Doppler analysis, jail-bar detection, a motion sensor arranged in or associated with the ultrasound probe 10 or any other known or hereafter developed apparatus or technique. Image analysis compares correlations between successive image data to determine whether something has moved. This is typically performed after the beamformer RF output data has been detected and log-compressed within the echo processor 36 with feedback path 110 as depicted in
Jail-bar detection analyzes the image or region of interest for jail-bars. This may include comparing brightness of type A synthesized lines to type B synthesized lines. Alternately, type A lines in a position P are compared to type B lines at position P at a later time. A disadvantage of the jail-bar detection approach is the need to periodically return to standard MAM in order to see if the motion has subsided.
The motion sensor may include a position, velocity, or acceleration detector arranged in the ultrasound probe. An example of this type of motion sensor is disclosed in U.S. Pat. No. 4,852,577. Other examples of motion detectors include magnetic position devices, such as the “Flock of Birds®” sensor manufactured by Ascension Technologies (Burlington, Vt.) or the “FASTRAK®” from Polhemus (Colchester, Vt.), and video imaging of markers on probes used to detect probe motion.
If excessive motion is detected, jail-bar mitigation is implemented, step 705. The jail-bar mitigation implemented at step 705 is maintained until the excessive motion is no longer detected, step 707.
The spatial filtering method applies a lateral low pass filter to the area in which excessive motion is detected. This equalizes the response of each beam but at the expense of blurring the image.
Temporal filtering with line interleaving flips the type A and type B line positions between frames (i.e., volumes in 3D imaging) and applies a heavy time average. An example of temporal filtering is disclosed in U.S. Pat. No. 5,980,458, the entire contents of which are incorporated herein by reference.
Dropping multi-line order includes dropping back from a high order multi-line such as 4× to a lower order multi-line such as 2× (i.e., 2 RX beams per TX beam). It is easier to control jail-bars in lower order multi-line with a post-detection lateral filter (i.e., a filter within echo processor 36 of
Reduction of the TX and RX apertures increases the beam sizes. This blurs the image but reduces certain types of jail-bars. Aperture sizes would be set in beamformer 30 of
Normalization determines the differences in average amplitudes of the various line types (i.e., A, B, C, D, etc.) and applies gains so the average brightnesses are equal between lines of different types. Normalization would be implemented in echo processor 36 and controlled by adaptive controller 300 via control path 210 in
For the sample sequence shown for lines l1 and l2 in
A number of samples from the two lines l1 and l2, taken from the input and the first tap of the delay line, respectively, are applied to a correlator 82. The correlator performs a cross-correlation of range aligned data samples of the two lines l1 and l2 to detect the condition of relative motion between the two lines. This cross correlation is performed in the conventional manner by sequentially shifting sample sequences from the two lines relative to each other, multiplying aligned samples after each shift, and summing the products to produce a correlation factor. The value and direction of shift for which the correlation factor is at a maximum indicates the amount and direction of motion that has occurred in the period between the acquisition of the two lines l1 and l2. The peak of the correlation factor is then used as the control input of a selector or multiplexer 84 to select the sample at the input of the selector which would, in the absence of motion, be range-aligned with the sample at the input of the delay line 80.
Thus, when sample Sb2 of line l2 is at the input of the delay line and there is no motion when the line samples are acquired, there would be a high degree of correlation between the two lines when no relative shift occurred, indicating no motion. The selector 84 would then select sample Sb1 to be used in calculating an interpolated value Xb with sample Sb2 at the output of summer 36. As shown in
But when there is motion away from the receiver (in a direction opposite that of the arrowheads in
In a similar manner motion toward the receiver would be detected by the correlator 82 and sample Sc1 would be selected to be used with sample Sb2 to interpolate value Xc, again at a range half increment due to motion.
The interpolated values may be at fractional range increments as the foregoing example illustrates. To bring the samples back into alignment with the same ranges as the received lines, the interpolated line values can be processed through an axial transversal filter with coefficients chosen to compute an interpolated value at each whole range increment along the interpolated line. This alignment would be useful if further line filtering were to be done using a multitap filter, for instance, or if the scan converter requires sample data points to be spatially organized in a uniform grid or pattern.
While the embodiment of
The above-embodiments describe a single threshold control in which the jail-bar artifact mitigation is either on or off. However, a graduated control may also be implemented in which one of a plurality of levels of jail-bar artifact mitigation is implemented depending the level of relative motion. Instead of a graduated control, a continuous control, i.e., sliding scale, approach may be also be implemented, wherein jail-bar artifact mitigation is initiated when excessive motion is detected and the degree of jail-bar mitigation is increased as the relative motion increases.
Thus, while there have been shown and described and pointed out fundamental novel features of the invention as applied to a preferred embodiment thereof, it will be understood that various omissions and substitutions and changes in the form and details of the devices illustrated, and in their operation, may be made by those skilled in the art without departing from the spirit of the invention. For example, it is expressly intended that all combinations of those elements and/or method steps which perform substantially the same function in substantially the same way to achieve the same results are within the scope of the invention. Moreover, it should be recognized that structures and/or elements and/or method steps shown and/or described in connection with any disclosed form or embodiment of the invention may be incorporated in any other disclosed or described or suggested form or embodiment as a general matter of design choice. It is the intention, therefore, to be limited only as indicated by the scope of the claims appended hereto.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB2006/050896 | 3/23/2006 | WO | 00 | 9/26/2007 |
Number | Date | Country | |
---|---|---|---|
60665757 | Mar 2005 | US |