SUPPRESSION OF REVERBERATIONS AND/OR CLUTTER IN ULTRASONIC IMAGING SYSTEMS
First Claim
1. A method for clutter suppression in ultrasonic imaging, said method comprising:
- transmitting an ultrasonic radiation towards a target medium via a probe;
receiving reflections of the ultrasonic radiation from said target medium in a reflected signal via a scanner, wherein the reflected signal is spatially arranged in a scanned data array, which may be one-, two-, or three-dimensional, so that each entry into the scanned data array corresponds to a pixel or a volume pixel (either pixel or volume pixel being collectively a “
voxel”
), and wherein the reflected signal may also be divided into frames, each of which corresponding to a specific timeframe is a cine-loop;
the method including the following steps;
step 110—
computing one or more self-similarity measures between two or more voxels or groups of voxels within a cine-loop or within a processed subset of the cine-loop, so as to assess their self-similarity;
step 120—
for at least one of;
(i) each voxel;
(ii) each group of adjacent voxels within the cine-loop or the processed subset of the cine-loop, and (iii) each group of voxels which are determined to be affected by clutter, based on one or more criteria, at least one of which relates to the self-similarity measures computed in step 110, computing one or more clutter parameters, at least one of which also depends on the self-similarity measures computed in step 110; and
step 130—
for at least one of;
(i) each voxel;
(ii) each group of adjacent voxels within the cine-loop or the processed subset of the cine-loop; and
(iii) each group of voxels which are determined to be clutter affected voxels, based on one or more criteria, at least one of which relates to the self-similarity measures computed in step 110,applying clutter suppression using the corresponding suppression parameters.
0 Assignments
0 Petitions
Accused Products
Abstract
There is provided a method for reverberation and/or clutter suppression in ultrasonic imaging, comprising: computing similarity measures between voxels within a cine-loop or subset of the cine-loop, to assess their spatial and/or temporal self-similarity; for at least one of: (i) each voxel; (ii) each group of adjacent voxels within the cine-loop or subset of the cine-loop; and (iii) each group of voxels which are determined to be affected by reverberations and/or clutter, based on one or more criteria; computing one or more reverberation and/or clutter parameters; and for at least one of: (i) each voxel; (ii) each group of adjacent voxels within the cine-loop or subset of the cine-loop; and (iii) each group of voxels which are determined to be reverberation and/or clutter affected voxels, based on one or more criteria, applying reverberation and/or clutter suppression using the corresponding reverberation and/or clutter suppression parameters.
15 Citations
33 Claims
-
1. A method for clutter suppression in ultrasonic imaging, said method comprising:
-
transmitting an ultrasonic radiation towards a target medium via a probe; receiving reflections of the ultrasonic radiation from said target medium in a reflected signal via a scanner, wherein the reflected signal is spatially arranged in a scanned data array, which may be one-, two-, or three-dimensional, so that each entry into the scanned data array corresponds to a pixel or a volume pixel (either pixel or volume pixel being collectively a “
voxel”
), and wherein the reflected signal may also be divided into frames, each of which corresponding to a specific timeframe is a cine-loop;the method including the following steps; step 110—
computing one or more self-similarity measures between two or more voxels or groups of voxels within a cine-loop or within a processed subset of the cine-loop, so as to assess their self-similarity;step 120—
for at least one of;
(i) each voxel;
(ii) each group of adjacent voxels within the cine-loop or the processed subset of the cine-loop, and (iii) each group of voxels which are determined to be affected by clutter, based on one or more criteria, at least one of which relates to the self-similarity measures computed in step 110, computing one or more clutter parameters, at least one of which also depends on the self-similarity measures computed in step 110; andstep 130—
for at least one of;
(i) each voxel;
(ii) each group of adjacent voxels within the cine-loop or the processed subset of the cine-loop; and
(iii) each group of voxels which are determined to be clutter affected voxels, based on one or more criteria, at least one of which relates to the self-similarity measures computed in step 110,applying clutter suppression using the corresponding suppression parameters. - View Dependent Claims (2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33)
wherein the term “
local”
refers to a certain voxel or to a certain voxel and one or more adjacent voxels in the scanned data array.
-
-
13. The method of claim 11, wherein the identification of reverberation and/or clutter affected voxels is based on comparing the one or more measures of temporal variability computed in step 110 to one or more corresponding thresholds (“
- identification thresholds”
), and when more than one measure of temporal variability is used, the identification of reverberation and/or clutter affected voxels is performed by applying one or more logic criteria to the results of comparing the measures of temporal variability to the corresponding identification thresholds, by applying an AND and/or an OR and/or a XOR and/or a NOT operator between the results.
- identification thresholds”
-
14. The method of claim 1, wherein step 120 includes the identification of reverberation and/or clutter affected voxels, wherein the identification of reverberation and/or clutter affected voxels is performed for each cine-loop and/or each frame and/or each spatial region within the cine-loop and/or one or more spatial regions within each frame.
-
15. The method of claim 13, wherein each of the identification thresholds is one of:
-
(a) Predefined, either as a global threshold or as a threshold which depends on the index of the entry into the scanned data array and/or on the frame index;
or(b) Adaptively determined for each cine-loop and/or each frame and/or each spatial region.
-
-
16. The method of claim 15, wherein the adaptive determination of the identification thresholds is performed employing the following method, wherein there is an assumption that the values of the temporal variability measure values are divided into two separate populations, one of which corresponds to reverberation and/or clutter affected voxels and the other to voxels substantially unaffected by reverberation and/or clutter, said method comprising the following steps:
-
(a) Select the set of voxels for which the identification threshold would be computed (the “
identification threshold voxel set”
);(b) Produce a list of the values of a temporal variability measure corresponding to the identification threshold voxel set, and sort this list in either ascending or descending order, to obtain the “
sorted temporal variability measure list”
;(c) For each element of the sorted temporal variability measure list; (i) Compute the mean value, m1, for all elements whose index into the sorted temporal variability measure list is lower than (alternatively, whose index is lower than or equal to) the current element'"'"'s index, and the mean value, m2, for all elements whose index is higher than (alternatively, whose index is higher than or equal to) the current element'"'"'s index; (ii) Compute the sum of the squared differences, S1, between the value of each element whose index is lower (alternatively, whose index is lower than or equal to) than the current element'"'"'s index and the value of m1; (iii) Compute the sum of the squared differences, S2, between the value of each element whose index is higher than (alternatively, whose index is higher than or equal to) the current element'"'"'s index and the value of m2; and (d) Set the identification threshold to the value of the temporal variability measure corresponding to the element of the sorted temporal variability measure list for which the value of S1+S2 is minimal.
-
-
17. The method of claim 14, wherein the identification of reverberation and/or clutter affected voxels is based on comparing the one or more measures of temporal variability computed in step 110 to one or more corresponding thresholds (“
- identification thresholds”
), wherein the identification thresholds are adaptively determined for each cine-loop and/or each frame and/or each spatial region, wherein there is an assumption that the values of the temporal variability measure values are divided into two separate populations, one of which corresponds to reverberation and/or clutter affected voxels and the other to voxels substantially unaffected by reverberation and/or clutter, and wherein the adaptive determination of the identification thresholds is performed employing the following steps;(a) Select the set of voxels for which the identification threshold would be computed (the “
identification threshold voxel set”
);(b) Produce a list of the values of a temporal variability measure corresponding to the identification threshold voxel set, and sort this list in either ascending or descending order, to obtain the “
sorted temporal variability measure list”
;(c) For each element of the sorted temporal variability measure list; (i) Compute the mean value, m1, for all elements whose index into the sorted temporal variability measure list is lower than (alternatively, whose index is lower than or equal to) the current element'"'"'s index, and the mean value, m2, for all elements whose index is higher than (alternatively, whose index is higher than or equal to) the current element'"'"'s index; (ii) Compute the sum of the squared differences, S1, between the value of each element whose index is lower (alternatively, whose index is lower than or equal to) than the current element'"'"'s index and the value of m1; (iii) Compute the sum of the squared differences, S2, between the value of each element whose index is higher than (alternatively, whose index is higher than or equal to) the current element'"'"'s index and the value of m2; and (d) Set the identification threshold to the value of the temporal variability measure corresponding to the element of the sorted temporal variability measure list for which the value of S1+S2 is minimal.
- identification thresholds”
-
18. The method of claim 1, wherein step 130 includes applying a reverberation and/or clutter suppression operator to reverberation and/or clutter affected voxels, as determined by step 120, wherein, at least one of the following suppression operators is employed:
-
(a) Set the signal value in reverberation and/or clutter affected voxels to a certain predefined constant; (b) Multiply the signal corresponding to reverberation and/or clutter affected voxels by a predefined constant, preferably between 0 and 1; (c) Subtract a certain predefined constant from the signal corresponding to reverberation and/or clutter affected voxels; (d) Apply a temporal high-pass or a temporal band-pass filter to reverberation and/or clutter affected voxels, so as to suppress the contribution of low temporal frequencies, the lower cut-off frequency of the filters being set so as to attenuate or to almost nullify low-frequency content; and (e) Replace the signal value in reverberation and/or clutter affected voxels by a function of the signal levels in their immediate spatial and/or temporal vicinity, said function being;
the mean;
weighted mean;
or median value of the signal level in the immediate spatial and/or temporal vicinity.
-
-
19. The method of claim 1, where applicable reverberation and/or clutter suppression parameters are one or more of the following:
-
(a) A temporal variability measure or a function of two or more temporal variability measures; (b) A function of (a), defined so that its values range from 0 to 1, receiving a certain constant value for voxels which are substantially unaffected by reverberation and/or clutter and another constant for voxels which are strongly affected by reverberation and/or clutter; (c) A Sigmoid function of (a);
-
-
20. The method of claim 1, wherein one or more of the following reverberation and/or clutter suppression operators are used in step 130 per processed voxel:
-
(a) Multiply the signal value by one or more reverberation and/or clutter suppression parameters; (b) Multiply the signal value by a linear function of the one or more reverberation and/or clutter suppression parameters; (c) Add to the signal value a linear function of the one or more reverberation and/or clutter suppression parameters; and (d) Apply a temporal high-pass filter or a temporal band-pass filter to the signal value, wherein the filter parameters depend on one or more reverberation and/or clutter suppression parameters.
-
-
21. The method of claim 1, wherein computing the one or more reverberation and/or clutter suppression parameters in step 120 includes detecting one or more ghost voxels or groups of voxels resulting from reverberation and/or clutter artifacts (the “
- ghost patterns”
) out of two or more similar voxels or groups of voxels (the “
similar patterns”
), said reverberation and/or clutter suppression parameters then being set so as to suppress ghost patterns without affecting the remaining similar patterns (the “
true patterns”
).
- ghost patterns”
-
22. The method of claim 21, wherein at least one of the following parameters is used to detect ghost patterns out of similar patterns (“
- ghost pattern parameters”
);(a) Mean signal magnitude and/or energy within each pattern and/or a subset of the voxels within each pattern; (b) Parameters derived from the spatial frequency distribution within each pattern and/or a subset of the voxels within each pattern; (c) Parameters relating to the information content within each pattern and/or a subset of the voxels within each pattern; and (d) Parameters relating to the distribution of the signal and/or signal magnitude and/or signal energy within each pattern and/or a subset of the voxels within each pattern.
- ghost pattern parameters”
-
23. The method of claim 21, wherein, in step 120, the detection of one or more ghost patterns out of two or more similar patterns also employs criteria based on whether one or more of the similar patterns is a ghost of one or more of the other similar patterns given one or more detected acoustic interfaces according to ghost image estimation by ray-tracing.
-
24. The method of claim 23, further comprising artifact sources search, and given the location of an artifact source, performing at least one of the following:
-
(a) For each artifact source, employing ghost image estimation by ray-tracing to assess the potential location of one or more ghosts of that artifact source which result from reverberations (the “
artifact source ghost targets”
);
or(b) For each scan line, estimating the contribution of the artifact sources to each range gate which results from sidelobe clutter, based on the sidelobe pattern of probe 26 for that scan line.
-
-
25. The method of claim 24, wherein the results of the artifact sources search are employed in step 110, for selecting the processed subset of the cine-loop, said processed subset of the cine-loop including, for each frame, at least one of:
- (i) one or more artifact sources; and
(ii) one or more artifact source ghost targets.
- (i) one or more artifact sources; and
-
26. The method of claim 24, wherein the artifact sources are selected by detecting continuous regions whose signal energy is relatively high.
-
27. The method of claim 26, wherein detecting continuous regions includes at least one of:
-
(a) Applying a non-linear filter to one or more frames of the cine-loop, said non-linear filter producing high values for areas where both the mean signal energy is relatively high and the standard deviation of the signal energy is relatively low; and
/or(b) Applying an energy threshold to the signal within one or more frames to detect high energy peaks, and then employing region growing methods to each such high energy peak to produce the artifact sources.
-
-
28. The method of claim 24, wherein the detection of artifact sources and/or of acoustic interfaces is performed by an edge detection and/or segmentation process.
-
29. The method of claim 1, further comprising:
-
tracking one or more patterns between two or more consecutive frames (“
pattern tracking”
), using a spatial registration method, said pattern tracking being utilized in at least one of the following steps;(a) In step 110, for selecting the processed subset of the cine-loop, wherein, once the processed subset of the cine-loop has been defined for a given frame (the “
subset reference frame”
), the processed subset of the cine-loop in one or more of the following frames is determined by pattern tracking for each voxel or group of voxels within the processed subset of the cine-loop for the subset reference frame;(b) In step 120, the detection of one or more ghost patterns out of two or more similar patterns also employs criteria based on the assumption that one of the similar patterns (“
similar pattern G”
) is considered more likely to be a ghost of another of the similar patterns (“
similar pattern O”
) if the relative motion of the two patterns over consecutive frames follow certain criteria, such as one or more of the following criteria;(i) The spatial distance traversed (between two or more consecutive frames) by the center of mass of each of similar pattern G and similar pattern O is approximately the same; (ii) The distance traversed (between two or more consecutive frames) along one or more specific axes by the center of mass of each of similar pattern G and similar pattern O is approximately the same; (iii) The angular rotation (between two or more consecutive frames) of similar pattern G and similar pattern O, with or without taking into account mirror reversal, is approximately the same; (iv) The magnitude of the angular rotation (between two or more consecutive frames) of similar pattern G and similar pattern O, with or without taking into account mirror reversal, is approximately the same; (v) The magnitude of the angular rotation (between two or more consecutive frames) of similar pattern G and similar pattern O, with or without taking into account mirror reversal, is approximately the same, but the angular rotations are in opposite directions; and (vi) The angular rotation (between two or more consecutive frames) with respect to a specific axis of similar pattern G and similar pattern O, with or without taking into account mirror reversal, is approximately the same.
-
-
30. The method of claim 1, wherein step 130 further comprises applying a reverberation and/or clutter suppression operator to reverberation and/or clutter affected voxels, as determined by step 120, wherein at least one of the following operators is employed:
-
(a) Set the signal value in reverberation and/or clutter affected voxels to a certain predefined constant; (b) Multiply the signal corresponding to reverberation and/or clutter affected voxels by a predefined constant, preferably between 0 and 1; (c) Subtract a certain predefined constant from the signal corresponding to reverberation and/or clutter affected voxels; (d) Replace the signal value in reverberation and/or clutter affected voxels by a function of the signal levels in their immediate spatial and/or temporal vicinity, said function may be, for example, the mean, weighted mean or median value of the signal level in the immediate spatial and/or temporal vicinity; and (e) For each group of spatially and/or temporally adjacent reverberation and/or clutter affected voxels (“
clutter affected voxel group”
), compute at least one of the following inter-voxel group parameters;(i) The ratio (“
voxel group ratio”
) between the value of a certain statistic of the signal and/or the signal magnitude and/or the signal energy and/or the signal videointensity within the group and within the corresponding group of voxels within the corresponding true pattern (“
true pattern voxel group”
), wherein the statistic may be, for example, the mean, median, maximum, certain predefined percentile, maximum likelihood;(ii) The relative angular rotation between the true pattern voxel group and the reverberation and/or clutter affected voxel group (“
voxel group angular rotation”
), using any registration technique known in the art, with or without taking into account mirror reversal; and(iii) The point spread function (PSF) that would approximately produce the clutter affected voxel group from the true pattern voxel group (“
voxel group PSF”
), wherein the PSF may be estimated after correcting for the voxel group ratio and/or applying mirror reversal and/or rotating the clutter affected voxel group to match the true pattern voxel group (or vise versa).
-
-
31. The method of claim 30, wherein after computing at least one of the inter-voxel group parameters, applying these parameters to the true pattern voxel group, by multiplying the true pattern voxel group by the voxel group ratio, and/or rotating the true pattern voxel group by the voxel group angular rotation, with or without mirror reversal, and/or applying the voxel group PSF to the true pattern voxel group, and subtracting the result multiplied by a certain constant, from the clutter affected voxel group.
-
32. The method of claim 1, wherein each of the one or more reverberation and/or clutter suppression parameters is at least one of:
-
(a) A similarity measure or a function or two or more similarity measures; (b) A parameter indicative of the probability for a voxel or a group of voxels to be reverberation and/or clutter affected voxels. The values of such parameters should increase with at least one of the following; (i) The value of one or more similarity measures (whose value increases with greater similarity) between the pattern to which the current voxel or group of voxels belong (the “
current pattern”
) and another pattern within the cine-loop increases;(ii) The value of one or more similarity measures (whose value decreases with greater similarity) between the pattern to which the current voxel or group of voxels belong (the “
current pattern”
) and another pattern within the cine-loop decreases; and(iii) Out of the group of patterns similar to the current pattern, including the current pattern itself, the current pattern is most likely to be a ghost pattern, based on ghost pattern parameters and/or on ghost image estimation by ray-tracing; (c) A function of (a) and/or of (b), defined so that its values would range from 0 to 1, receiving a certain constant for voxels which are substantially unaffected by reverberation and/or clutter and another constant for voxels which are strongly affected by reverberation and/or clutter; and (d) The result of applying a spatial and/or a temporal low-pass filter, using any method known in the art, to (a) or (b) or (c).
-
-
33. The method of claim 1, wherein one or more of the following reverberation and/or clutter suppression operators are used in step 130 per processed voxel:
-
(a) Multiply the signal value by one or more reverberation and/or clutter suppression parameters; (b) Multiply the signal value by a linear function of the one or more reverberation and/or clutter suppression parameters; and (c) Add to the signal value a linear function of the one or more reverberation and/or clutter suppression parameters.
-
Specification