The present invention relates to an apparatus for reducing side lobes in ultrasonic images using a nonlinear filter, and more particularly, to an apparatus for reducing side lobes in ultrasonic images capable of improving ultrasonic image quality by removing a side lobe signal component using a nonlinear filter that estimates the side lobe signal using channel data resulting from applying a focusing delay to signals received by a plurality of receiving elements of an array transducer, subtracts the side lobe signal estimated from a receive focusing delayed channel data, and uses the magnitude of the estimated side lobe signal as the filter coefficients of a nonlinear side lobe reduction filter.
Generally, ultrasonic images are used in diagnosing lesions. In medical ultrasonic imaging ultrasonic signals are transmitted via a transducer and the magnitude of ultrasonic signals received after being reflected from an inside of a human body is converted to brightness.
Despite the advantages of safety and real-time imaging capability, the ultrasonic images have a low resolution problem compared to other medical images. To solve this problem, a method of using an array transducer to focus ultrasonic waves of short pulse widths in order to transmit and receive the ultrasonic waves is being applied in a general medical ultrasonic imaging system.
Taking a close look at ultrasonic field in the ultrasonic focusing system, we can see that a main lobe is formed with respect to the scan line direction (i.e., axial direction) of ultrasonic image and that side lobes are formed at both sides of the main lobe due to leakage of ultrasonic signals. When the echoes from the target in the main lobe direction are received, signals from the target in the side lobe directions are also received with the result that the signals of the reflector in the side lobe act as noise in ultrasonic images and lower the resolution of the ultrasonic images.
Consequently, various attempts for reducing side lobes in ultrasonic images are recently being made, and details thereof are disclosed in detail in [Document 1], [Document 2], etc. as below.
However, in the cases of [Document 1] and [Document 2] below, since a method of applying weighting values to each received channel data is used, there is a problem of having to perform excess computation to reduce side lobes, and this problem is aggravated further as the number of channels increases.
(Patent Document 1) [Document 1] Korean Unexamined Patent Application Publication No. 2009-0042152 (published on Apr. 29, 2009)
(Patent Document 2) [Document 2] Korean Registered Patent No. 971433 (announced on Jul. 14, 2010)
(Non-patent Document 1) None
The present invention has been devised to solve the above-mentioned problems of the prior art, and the present invention is directed to providing an apparatus for reducing adverse effects of side lobes in ultrasonic images using a nonlinear filter capable of filtering a side lobe signal component, in which a waveform of the side lobe signal is estimated after obtaining channel signals by focusing of signals received via an array transducer having a plurality of receiving elements, the estimated side lobe signal is subtracted in a process of summing the channel signals delayed for focusing, and the estimated side lobe signal is used as filter coefficients.
According to an embodiment of the present invention, an apparatus for reducing side lobes in ultrasonic images using a nonlinear filter includes an array transducer to receive an ultrasonic signal reflected from an imaging point from each of receiving elements to output the reflected ultrasonic signal as a channel signal of a corresponding receiving element, a focusing module to temporally align the channel signals of the receiving elements, a summation unit to add the temporally aligned channel signals and output the summed signal in order to form an ultrasonic image, a side lobe computation module to calculate a waveform of a side lobe signal generated due to a leakage of the ultrasonic signal, and a filter unit to filter the summed signal of the summation unit based on the magnitude of the side lobe signal computed by the side lobe computation module in order to improve ultrasonic image quality.
In addition, the filter unit may include a subtraction filter to perform a first filtering by subtracting the side lobe signal from the summed signal and a nonlinear filter to perform a second filtering whose input is fed from the signal firstly filtered by the subtraction filter.
In addition, the nonlinear filter may perform the second filtering using [Expression 1] below.
Here, Bpixel is a brightness value of an ultrasonic image pixel, B filtered is a brightness value of a filtered pixel, γ is a scale factor, and QFp is an image quality factor at each image pixel.
In addition, the image quality factor is the sum of the side lobe signal waveform calculated by the side lobe computation module in order to evaluate the ultrasonic image quality.
In addition, when calculating the spatial frequency of a sinusoidal signal waveform in the received channel data, the side lobe computation module uses zero-appending to extend the length of the received channel data.
As above, according to the present invention, an apparatus for reducing side lobes in ultrasonic images using a nonlinear filter is capable of improving an ultrasonic image quality by removing a side lobe component by filtering the side lobe component with a nonlinear filter that calculates a waveform of the side lobe signal in the received channel signals after applying focusing delay to the signals received via an array transducer having a plurality of receiving elements, subtracts the side lobe signal calculated using a subtraction filter in a process of summing up the channel signals delayed for focusing, and uses an output of the subtraction filter as an input to a second filter.
Hereinafter, an apparatus for reducing side lobes in ultrasonic images using a nonlinear filter according to an embodiment of the present invention will be described in detail with reference to the accompanying drawings.
As illustrated in
The array transducer 20 has a plurality of elements configured to transmit ultrasonic waves into a human body and receive signals reflected from a tissue of the human body.
The signals received from the human body tissue (i.e., an imaging point) arrive at each of the receiving elements at different times due to positions at which the receiving elements of the array transducer 20 are arranged. The focusing module 21 applies a time delay to each of the plurality of channel signals in which a difference in arrival time has occurred as above in order to perform focusing which temporally aligns the channel signals as if the channel signals had arrived at the same time.
In addition, the side lobe computation module 22 serves to compute a waveform of a side lobe signal using a spatial frequency characteristic of side lobe signal components included in the channel signals delayed for focusing by a method to be described below.
The filter unit 24 may improve ultrasonic image quality using a plurality of filters 25 and 26 that remove a side lobe signal component using the calculated waveform of the side lobe signal. A subtraction filter 25 performs a first filtering by subtracting the sum of the waveform of side lobe signal calculated by the side lobe computation module 22 from a summed signal of the summation unit 23, and a nonlinear filter 26 performs a second filtering on the signal firstly filtered by the subtraction filter 25.
As illustrated in
In this manner, when the reflected signals impinge on the transducer from a direction adjacent to and at random angles with the scan line direction of the transducer, the signals are incident on the receiving elements with different phases.
Consequently, the signals incident from the random incident angles are shown as signals having a specific frequency referred to as a spatial frequency when viewed from the transducer.
The spatial frequency will be described with reference to
A wavelength λ of a sinusoid observed at a position in the x-axis when an ultrasonic wave having a wavelength λ0 is incident at a predetermined angle θ can be expressed as [Equation 1].
The number of cycles of a sinusoid appearing in a channel signal received from a receiving element of a transducer whose size is D is defined as a cycle per aperture (CPA), and this can be written as [Equation 2] below. Here, the CPA refers to a spatial frequency of a signal periodically appearing across the aperture of a transducer array.
According to [Equation 2], the spatial frequency of the sinusoid appearing on the receiving element of the transducer whose size is D varies in accordance with the incident angle θ.
The received channel signal can be written as [Equation 3] below.
Here, Xk is a channel signal received by a kth receiving element at time n, and the sum of all of the channel signals is s(n). Signals simultaneously incident on a receive channel from various directions may be modeled as a sum of sinusoids having various spatial frequencies in accordance with an incident angle.
The discrete Fourier transform of [Equation 3] can be shown as [Equation 4] below.
[Equation 4] is modeled by summing sinusoids having frequencies which are integer multiples of the reciprocal of the channel length, and the signals summed in a receive focusing process can be expressed as [Equation 5] below using [Equation 3] and [Equation 4].
Consequently, sinusoidal components having frequencies which are integer multiples of the reciprocal of the channel length are removed from a finally focused signal and only a direct current (DC) component remains. The DC component of the channel signal becomes a signal that arrives at all receive channels with the same phase.
Referring to
Depending on the incident angle, the received channel signals may be modeled as a sum of sinusoids having various frequencies. Modeling the channel signals as a sum of sinusoids Em(n) having integer frequencies and sinusoids Om(n) having frequencies of (integer+0.5), we obtain [Equation 6] as follows.
EO(n) represents a DC component of a channel, k represents a channel number, m represents a null or a side lobe index, and n represents a sampling time.
As shown in [Equation 5], when the sinusoids Em(n) having the integer frequencies are summed, the sum becomes zero such that they are removed in a focusing process. However, when the sinusoids Om(n) having the frequencies of (integer+0.5) are summed in the focusing process, the positive and negative signals cannot be canceled out and remain as a side lobe.
A method for calculating a signal component having a frequency of (integer+0.5) from channel signals will be described.
The conventional discrete Fourier transform is not appropriate for accurately estimating a waveform of a side lobe signal since it calculates a sinusoidal signal having an integer spatial frequency. Consequently, the data length needs to be extended by appending an appropriate number of zeros in order to calculate a magnitude of a sinusoid having a frequency of (integer+0.5) using the orthogonality of sinusoids. This process is referred to as zero appending.
For example, in the case of a system having 64 channels (when there are 64 receiving elements of the transducer) an extended data length of the first side lobe is given by [Equation 7] below.
Here, round( ) is a function that makes the argument take an integer value by the process of rounding off. Since CPA=2 in the extended data when an amplitude of a first side lobe is calculated using the orthogonality principle, this can be written as [Equation 8] below.
Here, M=85. Ai(n) corresponds to a complex amplitude of a first side lobe. From [Equation 8] the channel data waveform of the first side lobe can be expressed as [Equation 9] below when CPA=1.5 is applied to the channel data in which the length N=64.
A waveform of a channel signal may also be calculated for other sinusoids in channel data having the frequency of (integer+0.5) using the same method after appending zeroes whose number corresponds to a half-wavelength of each spatial frequency. Since all side lobe signals have different frequency components, each of them should be separately calculated after appropriately extending to a different channel data length. Consequently, a case of calculating all the waveforms of side lobes up to a degree P and adding all of the channel data can be represented as follows:
sidelobep(n) is a sum of side lobe signals up to the degree P at a corresponding pixel, where the index P indicates the Pth side lobe, and this value is a clutter component that degrades ultrasonic image quality. Consequently, the ultrasonic image quality may be improved by removing the clutter component.
The side lobe computation module 22 may compute the magnitude of a side lobe signal by following the above computation process.
The sum of side lobe signals in the received channel signal calculated using [Equation 10] corresponds to a magnitude of a signal that degrades image quality. Thus, to evaluate the image quality, an image quality factor may be defined as QF, which can be expressed as follows:
QFp=sidelobep [Equation 11]
In [Equation 11], a time index n in the depth direction is removed for generalization.
The sum of side lobe signals calculated by the side lobe computation module 22 may be subtracted from the channel data in the focusing process to reduce the adverse effects of the side lobes. This operation can be expressed as:
B
filtered
=B
pixel
−QF
p [Equation 12]
The subtraction filter 25 of the filter unit 24 may perform a first filtering, which subtracts the sum of the side lobe signals calculated by applying [Equation 12].
Considering the fact that a wideband transmit pulse is usually used in an ultrasonic imaging system, the ultrasonic image quality cannot be sufficiently improved when we employ the subtraction process only.
To cope with the above problem, the filter unit 24 applies the nonlinear filter 26. That is, the nonlinear filter 26 may be applied as a filtering means capable of further increasing a filtering effect of removing the side lobe components from the ultrasonic images.
Hereinafter, a method of designing the nonlinear filter 26 will be described.
When it is expected that a poor image quality may result due to an increase of clutter in ultrasonic images, a nonlinear filter can be designed so as to decrease the brightness of a corresponding pixel using the following equation:
Here, Bpixel is a brightness value of an ultrasonic image pixel, and Bfiltered is a brightness value of a filtered pixel. The effect of the filter may be controlled using the scale factor γ and the image quality factor QFP of the degree P.
The nonlinear filter designed as in [Equation 13] suppress the magnitude of a signal coming from a direction of a side lobe at each image pixel.
The clutter may be further reduced when the subtraction filter 25 to which [Equation 12] is applied and the nonlinear filter 26 to which [Equation 13] is applied are combined. That is, when a nonlinear filter in which an output of [Equation 12] is input to [Equation 13] is designed, this can be represented as follows:
The nonlinear filter 26 whose input is the output of the subtraction filter 25 is applied according to [Equation 14]. The ultrasonic image quality may be further improved by more faithfully removing the side lobe components.
Hereinafter, a test result using computer simulations to verify an effect of removing side lobes from ultrasonic images according to the present invention will be described.
A wideband ultrasonic pulse was used, and a transmit pulse had a Gaussian shape with a duration of 5 cycles. Conditions for performing the computer simulation are shown in [Table 1] below.
A point spread function (PSF) of a wire target at a 35 mm depth in a phantom used for observing an effect of filtering side lobe signals was calculated. Here, an image had a 10-mm width and a 2-mm height, the image was log compressed over a dynamic range of 50 dB, and apodization was not applied to the image.
As shown in
In
The results of calculating contrast, contrast-to-noise ratio (CNR), and SNR in the regions marked with white windows in
Although the contrast significantly increases when the side lobes were removed as shown in
A wire target placed at a 38-mm depth in a background of random scatterers was imaged, and
The prior art is shown in the solid line G41, a case of subtracting side lobe components is shown in the dotted line G42, an application of the nonlinear filter of [Equation 13] is shown in the dash-dot line G43, and an application of the nonlinear filter of [Equation 14] is shown in the dashed line G44, where a scale factor γ=10 was used in each of the last three cases. Although the width of the main lobe increased as shown in the dotted line G42 when the side lobe signals were subtracted, it can be seen that both the main lobe width and the side lobes were decreased as can be seen in the dash-dot line G43 and the dashed line G44, which correspond to [Equation 13] and [Equation 14], respectively.
The prior art is shown in the solid line G51, a case of subtracting the side lobe components is shown in the dotted line G52, an application of the nonlinear filter of [Equation 13] is shown in the dash-dot line G53, and an application of the nonlinear filter of [Equation 14] is shown in the dashed line G54. A scale factor γ=10 was used in each of the last three cases.
Test conditions for obtaining channel data from wire targets in a water tank and a human body are presented in [Table 3] below.
The wire targets in the water tank were vertically placed at intervals of 10 mm over a depth of 10 mm to 70 mm.
In
The area indicated by an arrow in
A scale factor γ=0.01, 0.1, 1, and 10 was used in obtaining
The granular pattern of the images gets finer as the scale factor γ is increased. Consequently, when making an actual diagnosis, it is preferable that an ultrasonographer control the ultrasonic image quality by adjusting the scale factor γ while observing the image quality.
Although particular embodiments of the present invention have been illustrated and described above, it should be apparent to those of ordinary skill in the art that the technical spirit of the present invention is not limited to the accompanying drawings and the above description and that various modifications are possible within the scope not departing from the spirit of the present invention. Also, the modifications should be viewed as belonging to the claims of the present invention within the scope not contrary to the spirit of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
10-2015-0117275 | Aug 2015 | KR | national |