1. Field of the Invention
The present invention relates to an object information acquiring apparatus and a control method thereof.
2. Description of the Related Art
A conventional ultrasound diagnostic apparatus used for medical image diagnoses employs an ultrasound probe including transducers having an ultrasound wave transmission/reception function. When an ultrasound beam formed from a synthesized wave of ultrasound waves is transmitted toward an object from the ultrasound probe, the ultrasound beam is reflected in an area of the object interior where acoustic impedance varies, or in other words a tissue boundary. By receiving an echo signal generated by the reflection and reconstructing an image on the basis of an intensity of the echo signal, a tissue condition in the object interior can be reproduced on screen as an ultrasound echo image.
Japanese Patent Application Publication No. 2009-28366 (Patent Literature 1: PTL 1) discloses a method of obtaining a three-dimensional ultrasound image of a wide area by performing a mechanical scanning operation using an ultrasound probe. More specifically, in this method, an ultrasound image is obtained while continuously moving a linear array probe in a direction (to be referred to hereafter as an elevation direction) which is orthogonal to and intersects an element array direction (to be referred to hereafter as a lateral direction). The linear array probe is capable of reconstructing a single tomographic slice image by performing electronic scanning using an ultrasound beam. Hence, by overlapping tomographic slice images created in respective positions in the elevation direction, a three-dimensional ultrasound image of an entire mechanical scanning area can be obtained. This three-dimensional image acquiring method is advantaged in terms of both speed and cost.
Japanese Patent Application Publication No. 2010-183979 (Patent Literature 2: PTL 2), meanwhile, discloses means for improving a resolution of ultrasound imaging using adaptive signal processing. A CAPON method, for example, is a type of adaptive signal processing using a spatial averaging method, which is employed in the radar field. The CAPON method serving as a type of adaptive signal processing may also be combined with a frequency domain interferometry (FDI) method. When adaptive signal processing is used, a frequency spectrum of a reception signal received during ultrasound imaging can be flattened with a high degree of precision, and as a result, an ultrasound image having a greatly improved spatial resolution in comparison with a conventional image can be obtained.
In the conventional example described in Japanese Patent Application Publication No. 2009-28366, however, on the single tomographic slice image reconstructed by performing electronic scanning with an ultrasound beam using the linear array probe, the image resolution in the elevation direction is much poorer than the image resolution in the lateral direction.
A first reason for this is that a pixel density in the elevation direction must be reduced to a certain extent. By reducing a scanning speed of the mechanical scan performed by the probe in order to increase a scanning pitch of the tomographic slice image, the pixel density in the elevation direction can be increased, but in this case, the duration of a physical load on an examinee increases. A second reason is that an effective aperture angle of the linear array probe in the elevation direction is smaller than the aperture angle in the array direction, and therefore a reconstruction resolution in the elevation direction is poor. This problem can be solved to a certain extent by using a two-dimensional array probe, but in this case, a required electrical circuit scale increases due to an increase in a number of transmission/reception elements, making practical application difficult in terms of cost.
The conventional example described in Japanese Patent Application Publication No. 2010-183979 describes means for improving the image resolution in the lateral direction in relation to a single tomographic slice image (a two-dimensional ultrasound image), but when this method is applied to the elevation direction, an increase in a calculation amount occurs. In other words, increases occur in the scale of a required signal processing circuit and an image memory, making practical application to an apparatus difficult in terms of cost. Moreover, when such a signal processing circuit is provided, a large increase in processing time may occur, making real time image display difficult.
The present invention has been designed in consideration of these problems, and an object thereof is to provide an object information acquiring apparatus with which image display speed and image resolution requirements can both be satisfied.
An object information acquiring apparatus according to the present invention is configured as described below.
More specifically, the present invention provides an object information acquiring apparatus comprising:
a probe in which a plurality of elements that receive acoustic waves propagating from an object and convert said acoustic waves into electric signals are arranged in at least a first direction;
a scanning unit that moves said probe in a second direction that intersects said first direction;
a generating unit that determines intensities of said acoustic waves in respective positions of an object interior using said plurality of electric signals, generates a plurality of first image data corresponding to tomographic images of said object in said second direction using a plurality of acoustic signals based on said intensities, and generates second image data using said plurality of acoustic signals; and
a display control unit into which said first image data and said second image data are input, and which displays on a display unit an image representing information relating to said object interior,
wherein said display control unit displays on said display unit a display based on said first image data, and
when said second image data are input from an identical position of said object, switches said display from said display based on said first image data to a display based on said second image data.
Further, the present invention provides an object information acquiring apparatus comprising:
a probe in which a plurality of elements that receive acoustic waves propagating from an object and convert said acoustic waves into electric signals are arranged in at least a first direction;
a scanning unit that moves said probe in a second direction that intersects said first direction;
a generating unit that determines intensities of said acoustic waves in respective positions of an object interior using said plurality of electric signals, generates a plurality of first image data corresponding to tomographic images of said object in said second direction using a plurality of acoustic signals based on said intensities, and generates second image data and third image data using said plurality of acoustic signals; and
a display control unit into which said first, second, and third image data are input, and which displays on a display unit an image representing information relating to said object interior,
wherein said display control unit displays on said display unit a display based on said second image data, and
when said third image data are input from an identical position of said object, switches said display from said display based on said second image data to a display based on said third image data.
Further, the present invention provides a control method for an object information acquiring apparatus having:
a probe in which a plurality of elements that receive acoustic waves propagating from an object and convert said acoustic waves into electric signals are arranged in at least a first direction;
a scanning unit that moves said probe in a second direction that intersects said first direction;
a generating unit that determines intensities of said acoustic waves in respective positions of an object interior from said plurality of electric signals, and generates image data using a plurality of acoustic signals based on said intensities; and
a display unit that displays an image of said object based on said image data,
the control method comprising the steps of:
generating by said generating unit a plurality of first image data corresponding to tomographic images of said object in said second direction;
displaying on said display unit a display based on said first image data;
generating by said generating unit second image data using said plurality of acoustic signals; and
switching by said display unit from said display based on said first image data to a display based on said second image data in an identical position of said object.
According to the present invention, it is possible to provide an object information acquiring apparatus with which both an image display speed and an image resolution can be improved.
Further features of the present invention will become apparent from the following description of exemplary embodiments with reference to the attached drawings.
Respective configurations of the present invention will be described in further detail below with reference to the drawings.
An object information acquiring apparatus according to the present invention uses a technique of obtaining object information in the form of image data by transmitting an acoustic wave to an object and receiving an acoustic wave (an echo signal) reflected in the object interior. An acoustic wave is a type of elastic wave, typically an ultrasound wave but also including elastic waves known as sound waves and ultrasound waves. A probe receives the acoustic wave propagating from the object interior. As described above, the information obtained from the object interior reflects differences in the acoustic impedance of tissue in the object interior.
An ultrasound diagnostic apparatus that performs a diagnosis on an object such as an organism will be described below as a representative example of the object information acquiring apparatus.
The ultrasound probe 4 is constituted by an oscillator in which electrodes are formed on either end of a piezoelectric ceramic, represented by PZT, or a piezoelectric material (a piezoelectric body) such as a polymer piezoelectric element, represented by PVDF, for example. Note that PZT is lead zirconium titanate, and PVDF is polyvinylidine difluoride. When the pulse-form or continuous wave-form transmission analog signal 100 is applied to the electrodes of the oscillator, the piezoelectric body expands and contracts. As a result of this expansion and contraction, pulse-form or continuous wave-form ultrasound waves are generated from the respective transducers, and by synthesizing these ultrasound waves, a transmission beam is formed. The respective transducers also expand and contract upon reception of propagating ultrasound waves, and generate electric signals as a result. These electric signals are output as the ultrasound wave reception analog signal 101. Here, elements employing different conversion methods may be used as the transducers. For example, the aforesaid oscillator may be used as the elements that transmit the ultrasound waves, while transducers employing a light detection method may be used as the elements that receive the ultrasound waves. A transducer employing a light detection method detects an ultrasound wave by converting the ultrasound wave into an optical signal, and is constituted by a Fabry-Perot resonator or a Fiber Bragg grating, for example.
The transmission/reception control unit 2 is controlled by software of the MPU 1 to control the transmission unit 3 and the reception unit 5 respectively on the basis of commands and information from an input operation unit. The transmission unit 3 is constituted by a pulsar drive circuit that supplies N channels of transducers constituting the ultrasound probe 4 with a number of transmission analog signals 100 corresponding to the N channels. The reception unit 5 first implements analog amplification processing on the weak reception analog signals 101 output from the N channels of transducers using a first stage LNA amplifier. Next, the reception unit 5 implements further analog amplification processing using a TGC (Time Gain Compensation) amplifier. Signals in an unnecessary frequency band are cut from the output of this amplifier using an AAF (Anti Aliasing Filter), whereupon A/D conversion processing is performed on each channel using a high-speed sampling (CLOCK) A/D converter. As a result, N channels of echo detection data 102 converted into reception digital signals are output.
An image generation unit 6 outputs two-dimensional image data 103 known as a B mode image by executing phase alignment processing, signal processing, and image generation on the input echo detection data 102. A DSC (Digital Scan Converter) 8 serves as display control unit for writing the input two-dimensional image data 103 (first image data, second image data, or the like) temporarily to an image storage unit 7 and outputting the two-dimensional image data 103 in the form of a video signal 104 in alignment with a timing of a horizontal synchronization frequency. A display unit 9 displays the B mode image upon input of the video signal 104.
Operations of the respective constituent elements will now be described in further detail.
The transmission/reception timings on the respective slice surfaces differ from each other. Here, a time from transmission to reception is calculated from a propagation distance and an acoustic velocity in order to adjust a reception time of the signals to be added in each reception element, whereupon reflection signals from identical P points are added together. In so doing, as shown in
Next, a point Q having an identical distance from a point S1 to the P point in the perpendicular direction in the slice surface SL#(n+1) will be considered. In this case, a triangle formed by the points S1, P, R1 and a triangle formed by the points S1, Q, R1 are clearly congruent, and therefore a time required to reach R1 from S1 via the P point is identical to a time required to reach R1 from S1 via the Q point. This relationship is identical in relation not only to the transmission/reception element in the R1 position, but also to the other reception elements in the same transmission/reception element group, and therefore, in positions of the slice surface SL#(n+1), identical added signals are obtained from a linear delay-and-sum result focusing on the P point and a delay-and-sum result focusing on the Q point. Hence, in the two-dimensional delay-and-sum processing performed in relation to the P point, linear delay-and-sum processing may be performed first on each slice surface to determine the delay-and-sum signals of the P point and the Q point, whereupon appropriate linear delay-and-sum processing is performed in the elevation direction to add together the delay-and-sum signals of the P point and the Q point.
More specifically, the image generation unit 6 according to the first embodiment, shown in
When the synthetic aperture method shown in
At the output timings according to the first embodiment, shown in
At the output timing according to this embodiment, shown in
Here, an outline of an operation performed during adaptive signal processing will be described. Adaptive signal processing is known in the field of radar as a method of estimating a target distance with a high degree of precision. Proc. Acoustics, Speech Signal Process., pp. 489-492 (March 2005) describes a method of improving resolution by employing adaptive signal processing when generating ultrasound echo image data. Further, a method in which both frequency domain interferometry (FDI) and adaptive signal processing are performed is known as a technique for improving spatial resolution in a depth direction. Conf Proc IEEE Eng Med Biol Soc. 2010; 1: 5298-5301 and Japanese Patent Application Publication No. 2010-183979 disclose results of an operation for forming an image of a layer structure of a blood vessel wall by performing an FDI method and a CAPON method, which is a type of adaptive signal processing, using electric signals output by a probe.
In adaptive signal processing, a processing parameter is varied adaptively in accordance with a reception signal. The CAPON method, which is a type of adaptive signal processing, is a method of processing a plurality of input signals such that in a condition where sensitivity to a focus position is fixed, power is minimized.
The FDI method is a method of analyzing reception signals (electric signals output from a probe) at each frequency, and estimating a reception power in a focus position using phase information relating to a plurality of frequency components. When a plurality of frequencies that are phase-aligned in a certain reference position are considered, a product of a distance from the reference position and a wave number is found to be proportionate to an amount of variation in the phase. In other words, when a certain focus distance is set and both the distance from the reference position to the focus distance and the frequency, or in other words the wave number, are known, it is possible to calculate a degree of phase variation. By applying the degree of phase variation to reception signals of respective frequencies and adding the reception signals together, the reception power at the focus distance can be estimated.
By combining the FDI method with adaptive signal processing, the reception power in the focus position can be estimated not based on a phase variation amount/weighting determined in advance in relation to the reception signals analyzed at each frequency component, but based on a phase variation amount/weighting calculated in accordance with the signals using adaptive signal processing.
In the present invention, the adaptive signal processing is not limited to the CAPON method, and a MUSIC method, an ESPRIT method, and so on may be used instead.
As described above, a typical ultrasound diagnostic apparatus forms an image by obtaining an envelope of a received waveform. When the FDI method and the CAPON method are applied in this case to improve the resolution further, it is envisaged that a plurality of reflection layers will exist in the FDI processing range. In an atmospheric observation radar, correlation between a plurality of reflection waves from the plurality of reflection layers can be suppressed by making an observation time sufficiently long, but during medical ultrasound imaging, the observation time of a single processing range is short, and therefore correlation between the plurality of reflection waves cannot be suppressed. A plurality of reflection waves from close reflection layers are therefore considered to have a high correlation.
It is known that when adaptive signal processing such as the CAPON method or the MUSIC method is applied as is to a plurality of reflection waveforms having a high correlation, unintended operations such as canceling out of a desired signal occur. By applying a frequency averaging method in this situation, operations of the FDI method and the CAPON method can be checked, and therefore a frequency averaging method is preferably used when FDI and CAPON are applied to medical ultrasound imaging.
According to Japanese Patent Application Publication No. 2010-183979, even when an observation subject having a different frequency characteristic exists, a calculation referral signal that takes the frequency characteristic of the subject into account can be generated by synthesizing reference signals. As a result, an improvement in spatial resolution in the depth direction can be achieved through adaptive signal processing.
Next, whitening processing is performed using the calculation referral signal, as shown in Equation 1 (S05). When the referral signal is g (t) and the Fourier transform implemented thereon is G (ω), whitening can be performed as shown in Equation (1), and as a result, a corrected correlation Hwhi (ω), which is a signal having a flattened frequency spectrum, can be calculated (S06). Note that η denotes noise power.
[Math. 1]
H
whi(ω)=H(ω)/(|G(ω)|2+η) (1)
Next, frequency domain interferometry and frequency averaging are applied to the flattened signal. More specifically, a correlation matrix R having i, j components is formed, as shown in Equation (2) (S07).
[Math. 2]
r
ij
=H
whi(ωi)Hwhi(ωj)H (2)
Next, a partial correlation matrix R′ is calculated using frequency averaging, as shown in Equation (3) (S08, S09). A depth direction power distribution P (r) is then estimated using the partial correlation matrix R′ thus determined (S10). Here, C is a constraint vector relative to a focus depth r, and kn is a wave number corresponding to an nth frequency. As a result of the processing described above, a plurality of depth direction power distributions corresponding to the plurality of calculation referral signals are calculated.
[Math. 3]
P(r)=1/(CHR′−1C)
C=[e
jk
r
, . . . , e
jk
r] (3)
Next, phases φ1 (f), φ2 (f) are determined by subjecting the two reference signals to Fourier transform (S21). An aopt for minimizing Σ(φ′(f)) 2 when φ′(f)=φ(f)−af is then searched for in relation to each reference signal. By determining the φ′(f) of the reference signals at aopt, the phase is flattened (S22).
Next, an amplitude and a phase are interpolated using a predetermined interpolation ratio (also referred to as an interpolation coefficient) α, as shown in Equation (4), whereby synthesized REF3 (f) is calculated (S23). Note that an interpolation ratio α is an arbitrary value that satisfies 0≦α≦1. Here, REF1 (f) and REF2 (f) are frequency components of f1 (t), f2(t) following power correction and phase correction, respectively.
[Math. 4]
REF1(f)=k′1exp(jφ′1(f))
REF2(f)=k′2exp(jφ′2(f))
REF3(f)={(1−α)k′1+αk′2(f)}exp(j((1−α)φ′1(f)+αφ′2(f))) (4)
Finally, a waveform of the calculation referral signal is determined by subjecting REF3 (f) to inverse Fourier transform (S24). Note that amplitude correction is preferably performed at this time to ensure that the signal power is fixed. The referral signal synthesis block 17 determines a plurality of calculation referral signals corresponding to a plurality of interpolation ratios α by varying the interpolation ratio α, and outputs the plurality of calculation referral signals to the adaptive calculation circuit 16. Values and variation steps of the interpolation ratios α may be set appropriately.
When the adaptive signal processing described above is performed, a process up to image generation is complicated. As a result, the processing time up to image generation is further increased in comparison with the delay-and-sum processing of the synthetic aperture method used in the first embodiment, and therefore the output timing shift increases further. In this embodiment, data obtained as a result of the adaptive signal processing correspond to the second image data.
At an output timing shown in
In a third embodiment, a case in which calculation circuits of two systems, namely the synthetic aperture processing and the adaptive signal processing described in the first and second embodiments, are provided simultaneously will be described. As regards the output timing of the calculation circuits for the two systems, the processing time of the calculation circuit used for adaptive signal processing is longer than that of the calculation circuit for the synthetic aperture processing, which is determined by the data amount of a tomographic slice image taken in advance, and therefore a time deviation occurs between image generation and output.
According to the embodiments described above, image generation can be performed in stages without impairing an image display speed in an ultrasound diagnostic apparatus that generates a three-dimensional ultrasound image using tomographic slice images obtained while continuously moving a linear array probe in an elevation direction. Simultaneously, an improvement in spatial resolution in the elevation direction can be achieved.
While the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.
This application claims the benefit of Japanese Patent Application No. 2011-201931, filed on Sep. 15, 2011, and, Japanese Patent Application No. 2012-175738, filed on Aug. 8, 2012, which are hereby incorporated by reference herein in their entirety.
Number | Date | Country | Kind |
---|---|---|---|
2011-201931 | Sep 2011 | JP | national |
2012-175738 | Aug 2012 | JP | national |