EP2820445A2 - Clutter suppression in ultrasonic imaging systems - Google Patents
Clutter suppression in ultrasonic imaging systemsInfo
- Publication number
- EP2820445A2 EP2820445A2 EP13707439.9A EP13707439A EP2820445A2 EP 2820445 A2 EP2820445 A2 EP 2820445A2 EP 13707439 A EP13707439 A EP 13707439A EP 2820445 A2 EP2820445 A2 EP 2820445A2
- Authority
- EP
- European Patent Office
- Prior art keywords
- local
- magnitude
- spatial
- compounding
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
Definitions
- the present invention relates generally to ultrasonic imaging systems, e.g., for medical imaging, and particularly to methods and systems for suppressing clutter effects in ultrasonic imaging systems.
- Ultrasonic medical imaging plays a crucial role in modern medicine, gradually becoming more and more important as new developments enter the market.
- One of the most common ultrasound imaging applications is echocardiography, or ultrasonic imaging of the cardiac system.
- Other widespread applications are obstetrics and gynecology, as well as abdominal imaging, to name a few.
- Ultrasonic imaging is also used in various other industries, e.g., for flaw detection during hardware manufacturing.
- Ultrasonic imaging systems typically produce relatively noisy images, making the analysis and/or diagnosis of these images a task for highly trained experts.
- One of the most problematic imaging artifacts is clutter, i.e., undesired information that appears in the imaging plane, obstructing data of interest.
- One of the main origins of clutter in ultrasonic imaging is effective imaging of objects outside the probe's mainlobe, also referred to as sidelobe clutter.
- sidelobe clutter For example, in echocardiography, the dominant reflectors outside the probe's mainlobe are typically the ribcage and the lungs.
- Another origin of clutter is multi-path reflections, also called reverberations.
- the geometry of the scanned tissue with respect to the probe, as well as the local reflective characteristics of the tissue causes a noticeable percentage of the transmitted energy to bounce back and forth in the tissue before reaching the probe.
- the signal measured for a specific range with respect to the probe may include contributions from other ranges, in addition to the desired range. If the signal emanating from other ranges is caused by highly reflective elements, it may have a significant effect on the image quality.
- a common method for enhancing the visibility of the desired ultrasonic image relative to the clutter, particularly in patients with low echogenicity (a common phenomenon among obese patients), is administering contrast agents.
- contrast agents enhance the ultrasonic backscatter from blood and aid in its differentiation from the surrounding tissue. This method is described, for example, by Krishna et al., in a paper entitled “Sub-harmonic Generation from Ultrasonic Contrast Agents,” Physics in Medicine and Biology, vol. 44, 1999, pages 681-694, which is incorporated herein by reference.
- An additional class of currently available methods for handling clutter is a family of clutter rejection algorithms, used in color-Doppler flow imaging. These methods estimate the flow velocity inside cardiac chambers or other blood vessels and suppress the effect of slow-moving objects; assuming that the blood flow velocity is significantly higher than the motion velocity of the surrounding tissue. These methods are described, for example, by Herment et al. in a paper entitled “Improved Estimation of Low Velocities in Color Doppler Imaging by Adapting the Mean Frequency Estimator to the Clutter Rejection Filter," IEEE Transactions on Biomedical Engineering, vol. 43, 1996, pages 919-927, which is incorporated herein by reference.
- Embodiments of the present invention provide methods and devices for reducing clutter effects in ultrasonic imaging systems.
- the present invention provides a method of ultrasound imaging, said method comprising: transmitting ultrasound radiation towards a target and receiving reflections of the ultrasound radiation from a region of the target in a main reflected signal and one or more auxiliary reflected signals, wherein each one of the reflected signals comprises an input dataset and is associated with a beam with a different and distinct beam pattern; compounding the input datasets from the main reflected signal and one or more auxiliary reflected signals, by the use of a compounding function, said compounding function using parameters derived from spatial analysis of the input datasets, said compounding function parameters being derived from at least one of:
- Fig. 1A is a schematic, pictorial illustration of an ultrasonic imaging system, in accordance with an embodiment of the present invention
- Fig. IB is a schematic, pictorial illustration of a probe used in an ultrasonic imaging system, in accordance with an embodiment of the present invention
- Fig. 2 is a schematic illustration of the beam pattern of two exemplary receive beams, defined in terms of gain as a function of the azimuth angle, in accordance with an embodiment of the present invention
- Fig. 3A is a schematic illustration of the expected magnitude ratio between the two receive beams of Fig. 2, defined for a point- like reflector as a function of the azimuth angle, wherein the azimuth scale is different than that in Fig. 2, providing a zoom on angles close to 0°, in accordance with an embodiment of the present invention
- Fig. 3B is a schematic illustration of the expected phase difference between the two receive beams of Fig. 2, defined for a point- like reflector as a function of the azimuth angle, wherein the azimuth scale matches that of Fig. 3A, in accordance with an embodiment of the present invention.
- the present invention relates to methods and systems for suppressing clutter effects in ultrasonic imaging systems.
- these methods and systems also result in enhancement of the ultrasound images' lateral resolution, i.e., the resolution along directions perpendicular to the centerline of an ultrasound beam.
- Fig. 1A is a schematic, pictorial illustration of an ultrasonic imaging system 20, in accordance with an embodiment of the present invention.
- System 20 comprises an ultrasound scanner 22 that scans organs of a patient using ultrasound radiation.
- a display unit 24 displays the scanned images.
- a probe 26, connected to scanner 22 by a cable 28, is typically held against the patient body in order to image a particular body structure, such as the heart (referred to as a "target").
- the probe may be adapted for insertion into the body, e.g., in transesophageal or transvaginal configurations.
- the probe transmits and receives ultrasound beams required for imaging.
- Scanner 22 comprises control and processing circuits for controlling probe 26 and processing the signals received by the probe.
- Fig. IB is a schematic, pictorial illustration of probe 26 used in imaging system 20, in accordance with an embodiment of the present invention.
- the probe comprises an array of transducers 30, e.g., piezoelectric transducers, which are configured to operate as a phased array.
- the transducers convert electrical signals produced by scanner 22 into a beam of ultrasound radiation transmitted into the patient's body.
- the transducers receive the ultrasound radiation reflected from different body tissues, and convert it into electrical signals sent to scanner 22 for processing.
- Probe 26 typically comprises several dozens and up to several hundreds of transducers 30 arranged in a horizontal linear array.
- the horizontal aperture of the array typically of the order of several centimeters, affects the minimal horizontal beam-width of the probe and the resulting horizontal angular resolution.
- the vertical beam-width may be adjusted by an acoustic lens.
- Some probes, referred to as "11 ⁇ 2 dimensional probes,” comprise several rows of transducers in the vertical dimension, providing a vertical sector-like beam pattern.
- Other probes comprise a complete two- dimensional (or multi-dimensional) array of transducers, enabling control over both horizontal and vertical directional beam patterns.
- the terms “horizontal” and “vertical” are used here solely for convenience, as the array may be positioned during imaging in any appropriate orientation relative to the patient' s body.
- a single array of transducers may generate different beam patterns, whose mainlobe may be pointed at different directions.
- the direction of the mainlobe is typically set by adjusting the relative phase and/or the relative time delay of the signals fed to the different transducers.
- the direction of the mainlobe is usually set by adjusting the phase shift and/or time delay introduced to the signal received by each transducer.
- apodization is a process of multiplying the output signal of each transducer by a multiplicative coefficient ("apodization coefficient"), before combining the outputs of all transducers to generate an overall array output signal.
- a reciprocal process is performed on transmission (in many cases a constant apodization is used on transmission).
- k be the transducer index (k should go over all transducers even if the transducer array comprises of more than one dimension), s k be the signal measured by transducer k, a k be the apodization coefficient of transducer k on reception, and (p k be the phase shift for transducer k on reception.
- Beamforming on reception, for obtaining the signal S at time t, may be performed using eq. (1):
- the signal reaching the different transducers may be sampled (by one or more analog-to-digital converters in scanner 22) at different beamforming stages, e.g., for each transducer separately, for different groups of transducers, or for all transducers combined. Eq. (1) and eq. (2) and combinations or variations thereof are applicable to all these cases.
- some systems perform sampling for all transducers combined (i.e., at the beam level), but use more than one set of apodization coefficients a k and/or more than one set of phase shifts (p k and/or more than one set of time-delays x k , a setting commonly referred to as multi-line acquisition, or MLA.
- MLA multi-line acquisition
- the beam pattern used on transmission is typically wider than those used on reception, so as to provide ultrasound energy to most or all of the volume covered by the different receive beams.
- Certain scanners 22 apply filtering matched to the transmitted waveform and/or down-conversion analogically, before sampling, whereas others apply matched filtering and/or down-conversion digitally.
- the matched filtering and/or down- conversion may be applied in any of the beamforming stages.
- Some scanners 22 sample the real data at one or more of the beamforming stages, whereas others use quadrature receivers, providing complex samples, including both a real in-phase component and an imaginary quadrature component.
- each receive beam is typically generated using a specific matched filter, corresponding to one of the transmission waveforms.
- a planar array is used instead of a phased array. In such cases, only a small set of predefined time-delays and/or phase-shifts is supported by the transducer array, wherein the beam formed by each such set is typically sampled separately. Further probe designs do not use transducer arrays, but instead allocate a single transducer to each transmit and/or receive beam.
- Probe 26 typically scans the region of interest, which may be one-dimensional (ID), two-dimensional (2D) or three-dimensional (3D), by one or more of the following:
- one or more auxiliary receive beams are used for suppressing clutter effects.
- the auxiliary receive beams may be generated using any configuration known in the art, including but not limited to those described hereinabove.
- Auxiliary receive beams may be allocated separately for each beam used normally for imaging ("main beam”), or for groups of main beams.
- the latter configuration may be useful, e.g., in MLA configurations, where multiple main beams are typically used concurrently on reception, corresponding to a single transmission beam, and one or more auxiliary beams may be defined.
- the auxiliary receive beams may be generated at the same time as the corresponding main receive beam or beams, or at different times. Combinations can also be considered, i.e., some of the auxiliary receive beams may be generated at the same time as the corresponding main receive beam or beams, and others may be generated at different times. Various sequences of relative timing between the different beams may be defined, e.g., alternating generation of main and auxiliary receive beams. If some or all of the auxiliary receive beams are generated at different times than the corresponding main receive beam or beams, separate transmission beams may be used for the main and auxiliary receive beams or groups thereof, using the same or different beam patterns.
- one main receive beam and one auxiliary receive beam are utilized, corresponding to a single transmission beam or two separate transmission beams.
- the auxiliary receive beam may have a different beam pattern than the main receive beam, e.g., it may be pointed at the same direction but its beam width may be wider (or narrower).
- multiple concurrent main receive beams e.g., 16 or 32
- one auxiliary receive beam is added, to be used for all concurrent receive beams.
- the beam patterns for all main receive beams are similar, but each main receive beam points at a slightly different direction.
- the auxiliary receive beam may be wider than some of them, e.g., its width may match that of the transmission beam.
- multiple main receive beams may be used with a single transmission beam, wherein one or more auxiliary receive beams may be added for each main receive beam or for each group of main receive beams.
- one or more receive beams may be treated as both main receive beams and auxiliary receive beams, in the context of different computations.
- compound information is provided by one or more main receive beams and one or more auxiliary receive beams so as to provide clutter suppressed outputs.
- the process may be applied in various processing phases of scanner 22, and to different types of data, e.g., analogically or digitally, to real or complex samples, before or after matched filtering and/or down- sampling.
- the clutter suppression process may be applied separately for different groups of transducers.
- "synthetic" main and auxiliary receive beams may be defined for computation purposes only, for each group of transducers separately, and the clutter suppressed outputs for the different transducer groups may be used for further stages of beamforming.
- the information for each main receive beam and/or auxiliary receive beam may be described as an array of real or complex measurements, each of which corresponding to a certain volume covered by the mainlobe (and sidelobes) of the applicable beam, between consecutive iso-time surfaces of the ultrasound wave within the medium (with respect to the probe 26), typically but not necessarily matching constant time intervals.
- Each such volume is commonly referred to as a volume pixel, or voxel.
- the samples are commonly referred to as range-gates, since the speed of sound does not change significantly while traversing soft tissues, so that iso-time surfaces can approximately be referred to as iso-range surfaces.
- a ID, 2D or 3D region of interest is scanned by scanner 22.
- the scanning may be performed by any method known in the art, including but not limited to the techniques mentioned hereinabove.
- different beams may have the same phase center but different beam directions.
- a polar coordinate system is typically used in 2D scanning
- a spherical coordinate system is typically used in 3D scanning.
- each voxel When using a polar coordinate system, the location of each voxel may be defined by the corresponding range-gate index and the angular direction of the beam with respect to the broadside of probe 26, wherein the probe's broadside is defined by a line perpendicular to the surface of probe 26 at its phase center, and wherein said angular direction may be defined in a Euclidian space by an azimuth angle, or by the u coordinate in sine-space.
- each voxel may be defined by the corresponding range-gate index and the angular direction of the beam with respect to the broadside of probe 26, wherein the angle direction may be defined either by the azimuth and elevation angles and/or by the (u,v) coordinates in sine-space.
- different beams may have different phase centers, wherein the beams may or may not be parallel to each other.
- a 2D or 3D Cartesian coordinate system may be used, wherein the location of each voxel is defined by its corresponding range-gate index and the location of the phase center.
- the dataset collected by one or more receive beams, which may or may not be concurrent, in each time-swath (“frame") may be organized in a ID, 2D or 3D array ("scanned data array"), using any voxel arrangement known in the art, wherein each index into the arrays relates to a different axis (e.g., in a polar coordinate system, a range-gate index and an azimuth index may be utilized), so that voxels which are adjacent to each other in one or more axes of the coordinate system also have similar indices in the corresponding axes.
- a axis e.g., in a polar coordinate system, a range-gate index and an azimuth index may be utilized
- the scanned data array may use a coordinate system which does not correspond to the scanning method, e.g., all receive beams may have the same phase center but different beam directions, but the scanned data array may utilize Cartesian rather than polar or spherical coordinates. This may be obtained by spatial interpolation and/or extrapolation, using any method known in the art, e.g., nearest neighbor interpolation, linear interpolation, spline or smoothing spline interpolation, and so forth.
- the scanned data array adheres to the following requirements, in which case the scanned data arrays are referred to as "common configuration scanned data arrays":
- All the data in a scanned data array corresponds to receive beams with the same beam pattern, except, perhaps, for translations and/or rotations.
- All the data in a scanned data array corresponds to receive beams using the same central frequency on reception.
- All the data in the scanned data array corresponds to receive beams used in conjunction with transmit beams using the same waveforms and the same transmit beam pattern, except, perhaps, for translations and/or rotations.
- the data in a scanned data array originates from receive beams which may be concurrent, divided into multiple groups each of which is concurrent (this is the most common case when using MLA configurations), or non-concurrent. In other embodiments, one or more of requirements (a)-(d) may not be met.
- the compounding of information provided by a main receive beam and one or more associated auxiliary receive beams is performed by computing a function of the local signal in two or more receive beams ("compounding function"), wherein the term local signal refers to a certain volume region, which is partially or wholly covered by more than one of the receive beams, e.g., a certain range-gate index of the different receive beams in cases where multiple receive beams have the same sampling frequency.
- the local signal in each such receive beam may be generated by setting the matched filters to fundamental imaging, i.e., using a reception frequency band centered at the center frequency of the corresponding transmission beam, or to harmonic imaging, i.e., using a reception frequency band centered at a frequency which is a natural number multiplied by the center frequency of the corresponding transmission beam (harmonies of the transmission frequencies).
- different receive beams may use different harmonies of the transmission frequencies, including but not limited to first harmonies (i.e., fundamental imaging).
- the output of the compounding function for each entry into the scanned data array and/or each frame may be used instead of the measured signal for the main receive beam in any further processing performed by ultrasound imaging systems, such as:
- Doppler processing e.g., for color-Doppler flow imaging, tissue-Doppler imaging, and/or pulse-Doppler studies. Note that Doppler processing typically entails analyzing data from two or more pulses for each applicable voxel.
- the compounding function may use parameters ("compounding function parameters") derived from at least one of:
- the local difference and/or ratio between spatial functions of the local magnitude and/or phase and/or complex signal in the different receive beams e.g., the ratio of local magnitude derivatives along a certain spatial axis; the ratio of local normalized magnitude derivatives along a certain spatial axis, wherein a local normalized derivative is defined as the local signal derivative divided by the local signal; the ratio of local magnitude spatial derivati es herein the spatial derivative is computed using a "Laplacian filter", i.e., the ratio of local magnitude normalized spatial derivatives, wherein the spatial derivative is computed using a "Laplacian filter”; the ratio of normalized magnitude spans, wherein a normalized signal span is defined as the difference between the maximal and minimal value, divided by the mean value, wherein an outlier rejection scheme may be applied prior to the statistical functions minimum, maximum and/or mean; and so forth.
- the ratio of local magnitude derivatives along a certain spatial axis e.g., the ratio of local normalized magnitude derivatives along a certain
- Temporal functions of the local phase difference and/or magnitude difference and/or magnitude ratio between different beams e.g., various temporal low-pass, band-pass or high-pass filters.
- Spatial-temporal functions of the local phase difference and/or magnitude difference and/or magnitude ratio between different beams e.g., spatial derivatives of the local phase difference and/or magnitude difference and/or magnitude ratio between different beams, applied after local temporal filtering; local difference and/or ratio between spatial functions of the magnitude and/or phase and/or complex signal in the different beams, applied after local temporal filtering; and so forth.
- spatial filtering may be applied to the scanned data array for a specific frame, whereas the temporal filtering may be applied to a specific voxel or group of adjacent voxels over multiple consecutive frames.
- non-clutter signals originate from angular directions close to that of the mainlobe of a main receive beam. For example, when using a main receive beam and an auxiliary receive, whose mainlobes point at the same spatial angle, one would expect that for reflectors at the center of the mainlobes, which by definition generate non-clutter signals:
- any local phase difference and/or magnitude difference and/or magnitude ratio may correspond to more than one angular direction of the signal source.
- a main receive beam and an auxiliary receive beam whose mainlobes point at the same spatial angle, wherein both the main receive beam and the auxiliary receive beam have a beam pattern which is symmetric about the center of their mainlobes (e.g., a conical beam pattern), at any given range with respect to the phase center of the transducer array
- the local phase difference and magnitude difference and magnitude ratio to be equal along concentric circular curves, whose center falls at the centerline of the mainlobes, wherein the curves lie in the corresponding range with respect to the phase center of the transducer array.
- Some deviations from this rule may be caused by physical effects within the medium, such as diffraction, which affect the actual beam pattern within the scanned region.
- the local phase difference By using more than one of the following parameters: the local phase difference, the local magnitude difference, and the local magnitude ratio, one may reduce the ambiguity in the information regarding the angular direction of the signal source. By utilizing information from multiple auxiliary receive beams for a given main receive beam, said ambiguity may be further reduced.
- the signal level as a function of spatial location may be defined, for example, by the local signal level, the local spatial derivatives and/or the local normalized spatial derivatives, wherein a local normalized spatial derivative is defined as the local spatial derivative divided by the local signal.
- the lateral resolution of various receive beams may be different due to differences in the corresponding beam patterns, and that the axial resolution of various receive beams may be different due to dissimilar waveforms on transmission and/or reception. Therefore, in some embodiments, before comparing the signal level and/or its magnitude as a function of spatial location between different receive beams, spatial filtering may be applied to one or more datasets, each of which corresponding to a specific type of beam pattern and/or waveform, in order to make sure that the lateral and/or axial resolution of the local signal in the compared beams would match.
- spatial filtering may be applied to one or more of the common configuration scanned data arrays before comparing the signal level and/or its magnitude as a function of spatial location between different receive beams. For example, when using narrow main receive beams and wide auxiliary receive beams, wherein for each main receive beam there is a corresponding auxiliary beam whose mainlobe points at the same direction, one may apply a spatial low-pass filter to the common configuration scanned data arrays acquired with the main receive beam. (c) Different target regions may have different dynamic models.
- the ribcage which often causes clutter reflections, moves at much slower rates than the heart muscle within the region of interest.
- Various temporal filters may be employed in order to estimate the signal levels for tissues with predefined dynamic properties.
- the clutter suppression processes described herein may be applied to specific frequency bands. Additionally or alternatively, one may apply the clutter suppression processes described herein separately to outputs of several temporal filters, and then compound the results in order to obtain the final dataset.
- the compounding function parameters may be derived and/or utilized in one or more of the following manners ("compounding function implementation manners”):
- Methods (b)-(e) are designed to reduce the time and/or memory complexity of the computations.
- one or more of the compounding function parameters may be defined as adhering to at least one of the following options:
- Such compounding function parameters are referred to herein as "generalized metrics.” They may be computed using any method known in the art. In some embodiments, e.g., in option (d), segmentation techniques may applied to at least one of the following (commonly referred to as the "segmentation technique parameters"): (a) the local received signal in one or more of the receive beams; and (b) the local value of one or more of the compounding function parameters.
- segmentation technique parameters commonly referred to as the "segmentation technique parameters”
- the local received signal in one or more of the receive beams and the local value of one or more of the compounding function parameters.
- radial search techniques see, for example, U.S.
- Patent 5,457,754 by Han et al, issued on Oct. 10, 1995, and title "Method for automatic contour extraction of a cardiac image" may be used for segmentation purposes. Such techniques may be combined with knowledge-based algorithms, aimed at performance enhancement, which may either be introduced during post-processing, or as a cost-function, incorporated with the initial boundary estimation.
- Another example for an applicable segmentation method is solving a constrained optimization problem, based on active contour models (see, for example, a paper by Mishra et al., entitled “A GA based approach for boundary detection of left ventricle with echocardiographic image sequences," Image and Vision Computing, vol. 21, 2003, pages 967-976, which is incorporated herein by reference).
- iterative region growing techniques may then be employed, wherein the predefined threshold for one or more of the segmentation technique parameters may be changed between consecutive iterations.
- the generalized metrics may be defined as real numbers whose values range from 0.0 to 1.0, wherein the value 0.0 is assigned to voxels that are substantially unaffected by clutter effects, and the value 1.0 is assigned to voxels in which substantially all of the measured signal emanates from clutter effects ("normalized generalized metrics").
- the normalized generalized metrics can be referred to as the estimated probability for the corresponding voxel to be substantially affected by clutter effects only; in option (c) immediately above, the normalized generalized metrics can be referred to as the estimated percentage of the received energy within the corresponding voxel that originates from clutter effects.
- Such definitions simplify signal calibration, especially when the compounding function is linear.
- signal calibration refers to making sure that the signal energy in voxels completely unaffected by clutter effects changes as little as possible due to the clutter suppression process, whereas the signal energy in voxels whose sole energy source is clutter would be suppressed as much as possible by the clutter suppression process.
- the generalized metric for a set of voxels in a set of frames may be derived using the following four-step process:
- compounding function parameter set (a) Select one or more compounding function parameters to be used ("compounding function parameter set").
- parameter #1 the local magnitude ratio between two different receive beams (a main and an auxiliary beam); parameter #2: the local phase difference between two different receive beams;
- parameter #1 the local ratio of spatial magnitude derivatives along a first spatial axis (e.g., the X axis) between two different receive beams
- parameter #2 the local ratio of spatial magnitude derivatives along a second spatial axis (e.g., the Y axis) between two different receive beams
- parameter #1 the local ratio of normalized spatial magnitude derivatives along a first spatial axis (e.g., the X axis) between two different receive beams
- parameter #2 the local ratio of normalized spatial magnitude derivatives along a second spatial axis (e.g., the Y axis) between two different receive beams
- a normalized local signal derivative is defined as the local signal derivative divided by the local signal
- parameter #1 the ratio of local magnitude spatial derivatives between two different receive beams, wherein the spatial derivative is computed using a "Laplacian filter”
- parameter #2 the local magnitude ratio between two different beams
- parameter #1 the ratio of local normalized magnitude spatial derivatives between two different receive beams, wherein the normalized spatial derivative is computed by dividing the local output of a "Laplacian filter" (applied to the local signal magnitude) by the local signal magnitude;
- parameter #2 the local magnitude ratio between two different receive beams;
- parameter #1 the ratio of the normalized magnitude spans between two different receive beams, wherein a normalized magnitude span is defined as the difference between the maximal and minimal magnitudes (with or without outlier rejection) within a local region, divided by the mean magnitude in said region, wherein the local region is defined so that its width is greater than its height
- parameter #2 the ratio of the normalized magnitude spans between two different beams, wherein the local region is defined so that its height is greater than its width.
- the results for each relevant voxel may be used as a generalized metric or a normalized generalized metric for that voxel. Multiple generalized metrics may thus be defined for each voxel.
- the compounding model may be defined using any method known in the art. Additionally or alternatively, it may be defined in one of more of the following ways:
- a predefined model (“predefined compounding model”), which may be theoretically determined and/or experimentally derived. Such a model may be based on the beam pattern of applicable transmit and/or receive beams. For example, the beam width in one or more axes of each such beam may be taken into account. It may be independent of spatial location and time, but it may also be dependent on spatial location and/or time. For example, one may use the fact that in the far-field, the nominal transmit beam width changes as a function of the distance from the probe, reaching a minimum at the focal length. Another example would be to use different models for different beam patterns, if more than one beam pattern type is in use.
- a simplified adaptive model (“simplified adaptive compounding model”), based on an assumption regarding the shape of the probability density function (PDF) for a voxel to be substantially affected (or substantially unaffected) by clutter effects as a function of the applicable compounding function parameter subset ("compounding model PDF").
- PDF probability density function
- compounding model PDF For example, a one-dimensional or multi-dimensional Gaussian function may be used, wherein the function's dimensions should correspond to the number of parameters in the compounding function parameter subset.
- the free parameters in the example of the Gaussian function - the vector of expectation values and the covariance matrix) may be determined by statistical analysis of the values of the compounding function parameter set over some or all of the relevant voxels in the appropriate receive beams.
- all relevant voxels in all frames may be used at once.
- the computations may be performed separately for each frame and/or each group of voxels.
- groups of voxels may be separated in one or more of the scanned data array axes (e.g., different range swaths and/or different angular sections with respect to the broadside of the probe).
- the computations may be performed for a subset of the frames and/or a subset of the voxels, wherein spatial and/or temporal interpolation and/or extrapolation may be used to estimate the compounding model PDF for any voxel.
- the interpolation and/or extrapolation may employ any method known in the art.
- an adaptive model (“adaptive compounding model”), wherein the compounding model PDF, i.e., the PDF for a voxel to be substantially affected (or substantially unaffected) by clutter effects as a function of the applicable compounding function parameter subset, is computed directly, without any prior assumptions regarding the shape of the PDF.
- the compounding model PDF may be computed using any method known in the art. For example, a one- or multi- dimensional histogram of the values of the applicable compounding function parameter subset for all relevant voxels may be calculated, wherein the histogram's dimensions should correspond to the number of parameters in the compounding function parameter subset. In order to adhere to the definition of a PDF, the histogram should be normalized such that the sum of all its elements would equal 1.0. If one assumes that:
- the primary source of the measured signal lies in a spatial angle close to the center of the mainlobe of the two-way beam pattern (this is because the beam pattern provides higher gain to this angular region, and most soft tissues have approximately the same reflective properties), and therefore voxels which are almost clutter free are also more ubiquitous than voxels which are strongly affected by clutter; and (ii)the compounding function parameter subset provides information regarding the local clutter levels and/or the angular direction of the primary source of the measured signal;
- the result may be used as an estimate of the PDF for a voxel to be substantially unaffected by clutter effects. If one goes over each histogram element, and replaces its value by 1.0 minus the original value, the result may be used as an estimate of the PDF for a voxel to be substantially affected by clutter effects.
- the histogram Before computing the histogram, one may apply outlier rejection in one or more of its axes, i.e., for one or more of the parameters in the compounding function parameter subset, thus limiting the histogram's dynamic range.
- the bins used in each histogram axis may be of a uniform width.
- the bins used in one or more of the histogram axes may be of adaptive width, wherein narrower bin widths are allotted to regions with higher PDF values.
- Voronoi diagrams may be used instead of histograms.
- the compounding model PDF may be computed using all relevant voxels in all frames at once.
- the compounding model PDF may be computed separately for each frame and/or each group of voxels.
- groups of voxels may be separated in one or more of the scanned data array axes (e.g., different range swaths and/or different angular sections with respect to the broadside of the probe).
- the compounding model PDF may be computed for a subset of the frames and/or a subset of the voxels, wherein spatial and/or temporal interpolation and/or extrapolation may be used to estimate the compounding model PDF for any voxel.
- the interpolation and/or extrapolation may utilize any method known in the art. For example, one may compute the compounding model PDF for the first frame out of every set of N f consecutive frames, and use that compounding model PDF for all N f consecutive frames.
- the predefined compounding model is expected to be most computationally efficient, but it cannot take into account medium- specific artifacts, e.g., local attenuation, scattering, refraction and diffraction effects, which may occur along the beam, changing its effective beam pattern.
- the adaptive compounding model is most accurate, yet it may be more computationally intensive. In order to reduce time complexity, one may use the compounding function implementation manners defined above.
- normalized generalized metrics are generalized metrics defined as real numbers whose value ranges from 0.0 to 1.0, wherein the value 0.0 is assigned to voxels that are substantially unaffected by clutter effects, and the value 1.0 is assigned to voxels in which substantially all of the measured signal emanates from clutter effects.
- normalized generalized metrics may be obtained instead of un-normalized generalized metrics by utilizing the probability for a voxel to be substantially affected by clutter effects only (as a function of the applicable compounding function parameter subset) instead of the compounding model PDF.
- the probability for a voxel to be substantially affected by clutter effects only (“compounding model probability function”) may be derived from the compounding model PDF as follows:
- ⁇ p n ⁇ be the compounding function parameter subset
- fpo F CiPn) be the compounding model PDF
- P( ⁇ p n ⁇ ) be the compounding model probability function
- P( ⁇ p n ⁇ ) for a given ⁇ p n ⁇ may be estimated by the sum over all values of fpo F CiPn)) m the model, which are lower than or equal to the value of f PDF for the specific ⁇ p n ⁇ .
- fpo F CiPn) is non " discrete
- an integral should be used instead of a sum. Note that a PDF is by definition normalized so that its sum (or integral) over the entire model is 1.0, so that the compounding model probability function ranges from 0.0 to 1.0, as expected.
- P( ⁇ p n ⁇ ) for a given ⁇ p n ⁇ may be estimated by the sum over all values of fpo F CiPn)) m me model, which are higher than or equal to the value of f PDF for the specific ⁇ p n ⁇ .
- fpo F CiPn) is non-discrete, an integral should be used instead of a sum.
- the compounding model PDF and/or the compounding model probability function may be utilized for estimating the local beam pattern within a given medium, and/or features of that beam pattern, such as the width of the mainlobe in one or more axes. For instance, if the transmit beam pattern is common to a main receive beam and its associated auxiliary receive beams, the receive beam pattern and/or features of the receive beam pattern may be estimated. Appropriate models should be defined for that purpose, which may be theoretically determined and/or experimentally derived ("beam pattern measurement models").
- the compounding model PDF and/or the compounding model probability function are expected to contain information regarding the elevation beam width on receive.
- a single transmit beam a single main receive beam and two auxiliary receive beams - one wider than the main receive beam in azimuth only, and the other in elevation only, one may use each auxiliary receive beam separately in order to deduce information regarding the beam width on receive along the applicable axis.
- the generalized metric and/or the normalized generalized metric may be utilized to estimate the angular direction of the primary energy source of the signal measured on receive for one or more voxels in one or more frames, using beam pattern measurement models. For example, when using a single transmit beam, a single main receive beam and a single auxiliary receive beam, wherein the azimuth beam width of all beams is the same, but the elevation beam width of the main receive beam is lower than that of the transmit beam and the auxiliary receive beam, the generalized metric and/or the normalized generalized metric are expected to be indicative of the elevation angle from which most of the received energy originates for each voxel or each frame.
- Various configurations of transmit beams, main receive beams and/or auxiliary receive beams may be utilized for estimating the spatial angle of the main energy source for each voxel on receive.
- the compounding function may be defined as a linear function of the local information for a main receive beam and one or more associated auxiliary receive beams. For example, if a single auxiliary receive beam is utilized, and a single normalized generalized metric m is calculated for each entry into the scanned data array p and/or each frame n separately (or for a subset of the entries into the scanned data array and/or for a subset of the frames), the compounding function may be defined as follows:
- the compounding function may be a linear combination of the local signal for the main receive beam, S main , and the auxiliary receive beam, S aux :
- the clutter suppression weight should therefore be computed in such a way that S out would be as close as possible to S free .
- this method may be able to extract relevant information from voxels in which clutter is the main source of the measured signal.
- a normalization factor may be introduces into eq. (6): ⁇ main ( ,n)-W( ,n)S aux ( ,n)
- S out for such voxels may the defined as equaling S main , so that the signal in such voxels remains unchanged.
- S out for such voxels may the set to the mean or median signal level of S out in the immediate vicinity of the current voxel, wherein the immediate vicinity may be defined in space and/or in time, and wherein voxels for which the value provided by eq. (7) has been replaced are excluded from the computation.
- Such techniques may also be used in other occurrences of "divide-by-zero" effects.
- the clutter suppression weight for each applicable entry into the scanned data array p and/or each applicable frame n may be estimated as follows:
- (a) Define a region surrounding (p, n), wherein the region may be spatial (i.e., including a set of near-by entries into the scanned data array) and/or temporal (i.e., including a set of consecutive frames, whose index is close to n) ("spatial-temporal region").
- the clutter suppression weight for (p, n) may be set to the local complex main-auxiliary ratio for to the voxel within the spatial-temporal region, whose generalized metric and/or normalized generalized metric values are most indicative of significant clutter effects.
- This exemplary function is based on the following assumptions: (i) the clutter signal contribution changes slowly enough in space and/or in time, so that one may consider it to be approximately constant within small spatial-temporal regions; and
- the voxel within the spatial-temporal region which is found to have the most substantial clutter effect may be treated as almost pure clutter, and therefore the local complex main- auxiliary ratio for that voxel is expected to significantly reduce the local clutter level.
- one may compute a linear combination of the outputs of a compounding function and the main receive beam signal, wherein the weights are based on the normalized generalized metric.
- a potential benefit of this technique is reducing the effects of the compounding function on clutter free voxels. For example, if we base our method on eq. (7), we obtain:
- two or more normalized generalized metrics may be defined and utilized by the compounding function.
- two normalized generalized metrics, ml and m2 we can use the following equation, wherein for each applicable entry into the scanned data array and/or for each applicable frame, the sum of all coefficients is, by definition, 1.0:
- S 1 and S 2 are the outputs of clutter suppression processes based on ml and m2 respectively, applied using any method herein above; and S 1 ⁇ 2 is the output of a clutter suppression process based on both ml and m2.
- S 1 ⁇ 2 may be computed as a simple local statistical function of S-L and S 2 , e.g., local minimum (between S-L and S 2 ), local average (between S-L and S 2 ), and so forth.
- two or more normalized generalized metrics may be defined and utilized in the compounding function, wherein at least one normalized generalized metric is used directly and at least one normalized generalized metric is used to compute a clutter suppression weight.
- ml and m2 wherein a clutter suppression weight W2 is computed using m2 the following equation may be used:
- a generalized metric and/or a normalized generalized metric for computing a clutter suppression weight and/or in a compounding function, it may undergo spatial filtering and/or temporal filtering. Similarly, before utilizing a clutter suppression weight in a compounding function, it may also undergo spatial filtering and/or temporal filtering.
- the filtering may be applied using the coordinate system of the scanned data array, wherein a time axis may be formed by combining the scanned data arrays collected in multiple frames, wherein consecutive frames are placed one alongside the other.
- any filter known in the art may be applied for that purpose.
- a linear low-pass filter, band-pass filter of high-pass filter may be applied in one or more of the spatial axes and/or the time axis.
- non-linear filters such as a median filter may be applied in one or more of the spatial axes and/or the time axis.
- Such filters are expected to reduce local irregularities in the generalized metrics, normalized generalized metrics and/or clutter suppression weights.
- non-linear filters applied in one or more of the spatial axes and/or the time axis, may be considered. For example, one may use the local minimum, local maximum, local mean plus local standard deviation, local mean minus local standard deviation, a certain percentile of the local signal, and so on.
- a local minimum is expected to reduce misdetections of voxels as clutter affected, since the presence of voxels for which the normalized generalized metric value is low (i.e., voxels determined to include low clutter levels, and therefore mostly relevant information) also reduces the normalized generalized metric value in their immediate vicinity.
- a local maximum is expected to increase the detection probability of clutter affected voxels, since the presence of voxels for which the normalized generalized metric value is high (i.e., voxels determined to include mostly clutter) also increases the normalized generalized metric value in their immediate vicinity.
- the output of the non-linear filter may be produced as follows: if the local median or the local mean or the local weighted mean is higher than a predefined value, use a certain statistical operator, and otherwise - another statistical operator. Examples for such statistical operators are the local minimum, local maximum, local mean, local mean plus local standard deviation, local mean minus local standard deviation, a certain percentile of the local signal and so forth.
- Such filters when applied to a normalized generalized metric, may allow one, for example, to increase the detection probability of clutter affected voxels in clutter affected regions and, at the same time, reduce the probability of misdetections in clutter free regions, assuming that clutter levels change relatively gradually in space and time.
- a transfer function may be applied to a generalized metric and/or a normalized generalized metric and/or a clutter suppression weight in various stages of the clutter suppression process. Any transfer function known in the art may be used, with various purposes in mind. Transfer functions applied to a normalized generalized metric, which by definition ranges from 0.0 to 1.0, should preferably: (a) yield values between 0.0 and 1.0; and (b) for each pair of values, wherein the second value is higher than the first one, the result of the transfer function for the second value should be higher than or equal to the result of the transfer function for the first value. [0098] For example, a transfer function may be applied to the normalized generalized metric m in eq. (3) and/or eq. (8) before using it for clutter suppression. Some exemplary transfer functions:
- a time axis may be formed by combining the scanned data arrays collected in multiple frames, wherein consecutive frames are placed one alongside the other.
- the spatial and/or temporal processing may comprise of various spatial and/or temporal filters known in the art.
- the local value of the output of a compounding function may be replaced by that of the corresponding voxel in the main receive beam if the local signal magnitude of the compounding function output is greater than the local signal magnitude of the main receive beam.
- Such a process may be considered, because the clutter suppression process is supposed to reduce clutter levels, without amplifying clutter free regions. For instance, high signal magnitudes may result from numerical effects such as "divide-by-zero.”
- (a) Define a block within the scanned data array, centered or approximately centered at the current voxel ("data block"). Depending on the dimensions of the scanned data array (which may include spatial and potentially also temporal axes), said block may be ID, 2D, 3D or four-dimensional (4D).
- the value of the current voxel within the compounding function output should be replaced by that of the corresponding voxel in the main receive beam if the ratio between the output of the previous step for the main receive beam and for the output of the compounding function is higher-than (or lower-than) than a predefined constant.
- the value of the current voxel within the compounding function output may be replaced by the mean signal, weighted mean signal or median signal within the data block of the compounding function output (either including or excluding the current voxel).
- Such a post process may be used, for instance, in order to take care of small dark regions within the compounding function output, which appear in bright continuous (and larger) regions of the image.
- a plurality of main receive beams may be used concurrently, e.g., in MLA configurations.
- one or more auxiliary receive beams may be defined, wherein each auxiliary receive beam is used in conjunction with one or more main receive beam.
- one possible clutter suppression method is to compute at least one of the following “compounding process functions":
- the common compounding model PDF may be transformed for each main- auxiliary beam pair before using it for generalized metric and/or normalized generalized metric computation.
- main- auxiliary beam pair data transformations one or more of which may be used when adding information to the common compounding model PDF (similar but inverse transformation may be employed when using the common compounding model PDF values for generalized metric and/or normalized generalized metric calculations):
- a constant may be added to one or more of the information's components, corresponding to one or more of the compounding function parameters referred to by the compounding model PDF.
- This transformation is referred to as translation of the compounding model PDF.
- two or more main receive beams and one auxiliary receive beam are used, wherein the beam patterns of the main receive beams are identical except for spatial rotation about the phase center of the transducer array, one may expect that the compounding model PDF calculated for each of the two or more main-auxiliary beam pairs would be similar, but translated in one or more of the compounding function parameters referred to by the compounding model PDF.
- either a linear or a non-linear transformation may be applied to one or more of the information's components, corresponding to one or more of the compounding function parameters referred to by the compounding model PDF.
- the linear or non-linear transformation may be applied to one or more compounding function parameters separately, but it may also be applied to two or more of the compounding function parameters combined.
- auxiliary receive beams may be used concurrently, either with a single main receive beam or with a plurality of main receive beams.
- One or more of the following techniques may be employed is such cases:
- One or more of the aforementioned clutter suppression schemes may be applied for each main-auxiliary beam pair separately. After that, for each applicable small spatial volume and/or for each applicable frame, a function may be applied to clutter suppression outputs for one or more of the associated main- auxiliary beam pairs ("multi-pair combining function"). Any function known in the art may be used, e.g., mean, maximum, minimum, median, a predefined percentile, weighted sum and so forth. For ordinal operators such as the median, the operator computation should be based on sorting the relevant values according to their magnitudes.
- the multi-pair combining function may also take into account one or more generalized metrics and/or normalized generalized metrics, as computed for one or more of the main- auxiliary beam pairs, e.g., a weighted sum may be calculated, wherein higher weights are assigned to values for main- auxiliary beam pairs in which the signal is considered to include lower clutter levels.
- clutter suppression output may then be referred to as a new (synthetic) main or auxiliary receive beam, to be used in conjunction with another auxiliary and/or main receive beam as inputs to one or more of the above described clutter suppression schemes. This process may be iteratively repeated several times.
- a compounding function may be applied to a main receive beams and a plurality of auxiliary receive beams, to produce a clutter suppressed output, wherein the compounding function employs one or more clutter suppression weights. For example, if we extend eq. (6) to relate to multiple auxiliary receive beams, we obtain:
- the weight vector W m in eq. (11), or similar coefficients in other clutter suppression equations referring to multiple auxiliary receive beam, may be estimated using any technique known in the art. For example, one may employ mean- square- error optimization, that is, for eq. (11), the weight vector W m for each applicable small spatial and/or temporal region may be set so as to minimize the expectation value of the output signal power, that is,
- Another example would be to compute the weight vector W m in eq. (11), or similar coefficients in other clutter suppression equations referring to multiple auxiliary receive beam, for each main- auxiliary beam pair separately, without taking the other main-auxiliary beam pairs into account.
- This methodology may be useful, for example, when the mainlobe of the main receive beam and the one or more auxiliary receive beams point at the same spatial direction, and for each main- auxiliary beam pair the beam width of the main and auxiliary receive beams are the same along a certain direction but different along another direction, wherein the directions along with the beams have different widths in different main- auxiliary beam pairs may or may not be perpendicular.
- different main-auxiliary beam pairs may be utilized to suppress clutter originating from different spatial angles with respect to the mainlobe of the main receive beam.
- a compounding function may be applied to a main receive beams and a plurality of auxiliary receive beams, to produce a clutter suppressed output, wherein the compounding function is linear, but does not employ clutter suppression weights, for instance, as described in eq. (3).
- the generalized metric and/or normalized generalized metric may be based on compounding function parameters which relate to more than one auxiliary receive beam associated with one or more main receive beams. For example, when using two auxiliary receive beams associated with a single main receive beam, one may use the following compounding function parameters: (i) the local relative magnitude between the main receive beam and the first auxiliary receive beam;
- a compounding function may be applied to a main receive beams and a plurality of auxiliary receive beams, to produce a clutter suppressed output, wherein:
- the compounding function employs one or more clutter suppression weights
- one or more of the generalized metrics and/or normalized generalized metrics may be based on compounding function parameters which relate to more than one auxiliary receive beam associated with one or more main receive beams.
- a compounding function may be applied to a main receive beams and a plurality of auxiliary receive beams, to produce a clutter suppressed output, wherein any compounding function may be employed.
- one or more receive beams may be employed as both main receive beams and as auxiliary receive beams, within the context of different main- auxiliary beam pairs.
- the phase centers of all associated auxiliary receive beams are identical to that of the corresponding main receive beam.
- there is at least one main- auxiliary beam pair for which the main and auxiliary receive beams have different phase centers (“non- concentric main-auxiliary beam pair").
- An exemplary use for such configurations is suppression of reverberation effects, in addition to sidelobe clutter effects.
- Reverberation artifacts measured for different receive phase centers are expected to be different, so that compounding functions may suppress both sidelobe clutter and reverberations.
- the scanned data array associated with the auxiliary receive beam may be translated and/or rotated with respect to the scanned data array associated with the corresponding main receive beam.
- spatial registration may be applied between the scanned data arrays associated with the main and auxiliary receive beams in non- concentric main- auxiliary beam pair ("main-auxiliary beam pair registration"). The registration may be performed for the entire cine-loop, for each frame separately, in each frame - for each entry into the scanned data array or group of adjacent entries, or for each entry into the scanned data array or group of adjacent entries in the entire cine-loop.
- the spatial registration may be performed using any technique known in the art.
- the spatial registration may be rigid, accounting for global translations and/or global rotation. Additionally or alternatively, the spatial registration may be non-rigid, also taking into account local deformations resulting from physical effects within the medium.
- some translation parameters such as the global translation and the global rotation, may be directly derived from the relative geometry of the centerlines of the main and auxiliary receive beams, which are known and controllable. Additionally or alternatively, some translation parameters may be initially estimated according to the relative geometry of the centerlines of the main and auxiliary receive beams, and then fine-tuned according to the respective acquired scanned data arrays.
- the spatial angle between the center of a voxel and the centerline of the main receive beam (“main spatial angle”) is typically different than the spatial angle between the center of said voxel and the centerline of the respective auxiliary receive beam (“auxiliary spatial angle”).
- auxiliary spatial angle the difference between the main spatial angle and the auxiliary spatial angle may be different for different voxels.
- the local relative magnitude and/or the local phase difference between the main and the auxiliary receive beam, as calculated for reflectors along the centerline of the main receive beam may change as a function of range (“parameter range variability").
- the compounding function may be adjusted in order to handle this parameter range variability.
- main-auxiliary beam pair data transformations as defined hereinabove, wherein the transformation for some voxels may be derived from the main spatial angle and auxiliary spatial angle for each voxel or group of voxels, wherein the main spatial angles and auxiliary spatial angles may be computed based on the geometric location of the center of the voxel with respect to the phase centers of the main and auxiliary receive beams, and/or from the outputs of the main-auxiliary beam pair registration process, if applied.
- U.S. patent 8,045,777 also referenced herein above, describes a clutter suppression method based on computing a linear combination of the signals in a main receive beam and in one or more auxiliary receive beams, wherein the coefficients of the linear combination ("temporal clutter suppression weights") are computed so as to suppress reflections due to clutter, wherein a reflection is determined as associated with clutter if the determined de-correlation time is above a specified threshold ("temporal clutter suppression process").
- temporary clutter suppression process When a single auxiliary receive beam is associated with each main receive beam, the following equation may be used:
- Ft is a temporal low-pass filter, applied for each applicable entry into the scanned data array and/or each applicable frame.
- This technique may be combined with the clutter suppression techniques hereof, using compounding functions which also depend on temporal clutter suppression weights.
- one or more generalized metrics and/or normalized generalized metrics may be used in addition to one or more temporal clutter suppression weights, as can be seen in the following exemplary equation:
- m is a normalized generalized metric, used so as to adjust the effect of the temporal clutter suppression process in each voxel based on the estimated local clutter contribution.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
Abstract
Description
Claims
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US13/408,136 US9451932B2 (en) | 2004-12-30 | 2012-02-29 | Clutter suppression in ultrasonic imaging systems |
PCT/IB2013/050392 WO2013128301A2 (en) | 2012-02-29 | 2013-01-16 | Clutter suppression in ultrasonic imaging systems |
Publications (1)
Publication Number | Publication Date |
---|---|
EP2820445A2 true EP2820445A2 (en) | 2015-01-07 |
Family
ID=47790277
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
EP13707439.9A Ceased EP2820445A2 (en) | 2012-02-29 | 2013-01-16 | Clutter suppression in ultrasonic imaging systems |
Country Status (3)
Country | Link |
---|---|
EP (1) | EP2820445A2 (en) |
CN (1) | CN104272134B (en) |
WO (1) | WO2013128301A2 (en) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9282945B2 (en) * | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
CN104105449B (en) | 2011-12-01 | 2018-07-17 | 毛伊图像公司 | Use the motion detection based on ping and porous doppler ultrasound |
EP2797515A4 (en) | 2011-12-29 | 2015-07-22 | Maui Imaging Inc | M-mode ultrasound imaging of arbitrary paths |
EP2816958B1 (en) | 2012-02-21 | 2020-03-25 | Maui Imaging, Inc. | Determining material stiffness using multiple aperture ultrasound |
KR102103137B1 (en) | 2012-03-26 | 2020-04-22 | 마우이 이미징, 인코포레이티드 | Systems and methods for improving ultrasound image quality by applying weighting factors |
KR102176193B1 (en) | 2012-08-10 | 2020-11-09 | 마우이 이미징, 인코포레이티드 | Calibration of Multiple Aperture Ultrasound Probes |
KR102176319B1 (en) | 2012-08-21 | 2020-11-09 | 마우이 이미징, 인코포레이티드 | Ultrasound imaging system memory architecture |
WO2014160291A1 (en) | 2013-03-13 | 2014-10-02 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
US9883848B2 (en) | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
JP6722656B2 (en) | 2014-08-18 | 2020-07-15 | マウイ イマギング,インコーポレーテッド | Network-based ultrasound imaging system |
CN107278271A (en) * | 2015-03-05 | 2017-10-20 | 克里斯塔维尔医学影像有限公司 | Clutter recognition in ultrasonic image-forming system |
WO2017132517A1 (en) | 2016-01-27 | 2017-08-03 | Maui Imaging, Inc. | Ultrasound imaging with sparse array probes |
CN106805997B (en) * | 2016-12-26 | 2020-08-07 | 乐普(北京)医疗器械股份有限公司 | Elastic imaging method and device |
JP7253496B2 (en) * | 2017-04-28 | 2023-04-06 | コーニンクレッカ フィリップス エヌ ヴェ | Power Doppler Imaging System and Method Including Improved Clutter Suppression |
WO2018206797A1 (en) | 2017-05-11 | 2018-11-15 | Koninklijke Philips N.V. | Methods and systems for controlling the generation of a compound ultrasound image |
CN110868908B (en) * | 2017-06-15 | 2022-04-08 | 富士胶片株式会社 | Medical image processing apparatus, endoscope system, and method for operating medical image processing apparatus |
WO2020061877A1 (en) * | 2018-09-27 | 2020-04-02 | 深圳迈瑞生物医疗电子股份有限公司 | Space compounding method and system, and computer-readable storage medium |
CN113640809B (en) * | 2021-08-23 | 2023-08-25 | 飞依诺科技股份有限公司 | Ultrasonic imaging noise reduction method and ultrasonic imaging equipment |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5457754A (en) | 1990-08-02 | 1995-10-10 | University Of Cincinnati | Method for automatic contour extraction of a cardiac image |
US6283919B1 (en) | 1996-11-26 | 2001-09-04 | Atl Ultrasound | Ultrasonic diagnostic imaging with blended tissue harmonic signals |
KR100232257B1 (en) * | 1997-09-04 | 1999-12-01 | 이민화 | Ultrasound color doppler imaging system for minimizing transient response to a clutter signal |
US6346124B1 (en) | 1998-08-25 | 2002-02-12 | University Of Florida | Autonomous boundary detection system for echocardiographic images |
BR0308375A (en) * | 2002-03-13 | 2005-01-11 | Raytheon Canada Ltd | Noise suppression system and method |
EP1842079A4 (en) * | 2004-12-30 | 2010-07-07 | Crystalview Medical Imaging Lt | Clutter suppression in ultrasonic imaging systems |
US8254654B2 (en) | 2007-10-31 | 2012-08-28 | University Of Southern California | Sidelobe suppression in ultrasound imaging using dual apodization with cross-correlation |
-
2013
- 2013-01-16 WO PCT/IB2013/050392 patent/WO2013128301A2/en active Application Filing
- 2013-01-16 CN CN201380022661.8A patent/CN104272134B/en not_active Expired - Fee Related
- 2013-01-16 EP EP13707439.9A patent/EP2820445A2/en not_active Ceased
Non-Patent Citations (2)
Title |
---|
None * |
See also references of WO2013128301A2 * |
Also Published As
Publication number | Publication date |
---|---|
WO2013128301A2 (en) | 2013-09-06 |
WO2013128301A3 (en) | 2013-11-07 |
CN104272134B (en) | 2017-05-31 |
CN104272134A (en) | 2015-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9451932B2 (en) | Clutter suppression in ultrasonic imaging systems | |
EP2820445A2 (en) | Clutter suppression in ultrasonic imaging systems | |
US20210378633A1 (en) | Method and apparatus to produce ultrasonic images using multiple apertures | |
US20130343627A1 (en) | Suppression of reverberations and/or clutter in ultrasonic imaging systems | |
US8045777B2 (en) | Clutter suppression in ultrasonic imaging systems | |
Szasz et al. | Beamforming through regularized inverse problems in ultrasound medical imaging | |
US8068647B2 (en) | Method and apparatus for real-time motion correction for ultrasound spatial compound imaging | |
US20130258805A1 (en) | Methods and systems for producing compounded ultrasound images | |
US9081097B2 (en) | Component frame enhancement for spatial compounding in ultrasound imaging | |
US10908269B2 (en) | Clutter suppression in ultrasonic imaging systems | |
US11087466B2 (en) | Methods and system for compound ultrasound image generation | |
US11096672B2 (en) | Clutter suppression in ultrasonic imaging systems | |
CN115869010A (en) | System and method for beamforming sound speed selection | |
Hossack | Influence of elevational motion on the degradation of 2D image frame matching | |
Khodadadi | Ultrasound elastography: Direct strain estimation | |
KR101610877B1 (en) | Module for Processing Ultrasonic Signal Based on Spatial Coherence and Method for Processing Ultrasonic Signal | |
Zhang | Speckle removal in medical ultrasound images by compounding and filtering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PUAI | Public reference made under article 153(3) epc to a published international application that has entered the european phase |
Free format text: ORIGINAL CODE: 0009012 |
|
17P | Request for examination filed |
Effective date: 20140602 |
|
AK | Designated contracting states |
Kind code of ref document: A2 Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR |
|
AX | Request for extension of the european patent |
Extension state: BA ME |
|
DAX | Request for extension of the european patent (deleted) | ||
17Q | First examination report despatched |
Effective date: 20150918 |
|
REG | Reference to a national code |
Ref country code: DE Ref legal event code: R003 |
|
STAA | Information on the status of an ep patent application or granted ep patent |
Free format text: STATUS: THE APPLICATION HAS BEEN REFUSED |
|
18R | Application refused |
Effective date: 20171129 |