The present invention relates to an ultrasound imaging technology for capturing an image inside a test object by using ultrasound waves.
Ultrasound imaging technology is a technology for non-invasively imaging the inside of a test object including a human body by using ultrasound waves (a sound wave not intended to be heard, sound waves at a high frequency equal to or greater than 20 kHz, generally). As an example, a medical ultrasound imaging apparatus will be described simply. An ultrasound probe transmits an ultrasound wave into the body of a patient and receives an echo signal reflected from the inside of the patient's body. The reception signal is subjected to signal processing in any one of or both the ultrasound probe and an ultrasound imaging apparatus main body. Thereafter, the reception signal is transferred to an image display unit, and an ultrasound image is displayed. More specifically, for example, a signal of a transmission beam is generated by a transmission beamformer in the ultrasound imaging apparatus main body, and then, the signal passes through a transmission/reception separation circuit (T/R), thereby being sent to the ultrasound probe. The ultrasound probe transmits ultrasound waves. After receiving an echo signal from the inside of a human body, the ultrasound probe transmits a signal to the imaging apparatus main body. In the imaging apparatus main body, the reception signal passes through the transmission/reception separation circuit again and is subjected to phasing processing by a reception beamformer, thereby being transmitted to an image processing unit. In the image processing unit, various types of image processing such as diverse filtering and scan converting are executed. Eventually, an ultrasound image is displayed by the image display unit.
In this manner, a general ultrasound diagnostic apparatus is configured by three technologies such as transmission beamforming, reception beamforming, and the back-end image processing. Particularly, since the beamformer performs signal processing at the RF (high frequency) level during transmission and reception, algorithm of the beamformer or mounting architecture decides the basic image quality of an ultrasound image. Therefore, the beamformer is a key unit in the apparatus.
The reception beamformer applies a delay time in which a delay amount is distributed in a concave surface-type in accordance with a relationship between a focal position and an element position with respect to each reception signal (reception data) of multiple elements configuring the ultrasound probe. The reception beamformer sets a focal point to be focused on a certain point in a virtual space and performs adding of reception signal data. This method is called phasing performed by a delay adding method. In the delay adding method, reception data received by the multiple elements of the ultrasound diagnostic apparatus is multiplied by a fixed weight vector which is saved in a diagnostic apparatus. Then, adding is performed after performing apodization. The method is applied to not only the reception beamformer but is applied similarly to the transmission beamformer.
Meanwhile, it is understood that there is a limitation on azimuthal direction resolution as a basic problem of the ultrasound imaging apparatus. Since transmission and reception of an ultrasound wave are performed by an array which has a finite open caliber, there is an occurrence of an influence to diffraction at the edge of the opening portion. It is possible to infinitely improve the resolution in accordance with an improved depth direction if an infinitely long array can be prepared. However, realistically, due to the physical limitation on apparatus design such as the length of a transmission/reception array, improvement of azimuthal direction resolution has been restricted. Recently, an attempt attracting attention thereto has been made. That is, the aforementioned fixed weight vector which is used for performing delaying during delay adding of the beamformer is adaptively varied with respect to every single item of transmission/reception data in chronological order, thereby acquiring a more accurate ultrasound image. Accordingly, it is possible to particularly improve azimuthal direction resolution which has been an essential problem in the beamforming technology.
Particularly in recent years, reports on a technology for improving azimuthal direction resolution have begun to be made. In the technology, technologies of adaptive signal processing including the MVDR method (minimum variance distortionless response; the capon method) that have been developed in the mobile communication field are applied to a beamformer of reception data. The aforementioned adaptive techniques are realized by adaptively varying a complex component in the weight vector for performing delay adding, based on a correlation matrix of reception data. In other words, in the related art, the weight vector used to be a fixed value, but in the adaptive technique, the weight vector is obtained through computation by using a reception signal for each sample point in a temporal direction of the reception signal, and the obtained weight vector is multiplied by the reception signal.
In adaptive signal processing, similar to the delay adding method in the related art, there is a problem of distortion in ultrasound wave-front resulting from distribution of heterogenous sound velocities in a medium, and a randomly scattered micro body. In adaptive signal processing, the focal point of the reception beamformer set in the apparatus is set while assuming that the medium is uniform in sound velocity and is homogeneous (homogeneity). Therefore, when there is a presence of distortion in sound wave propagation, the distortion may lead to blurring of an image or image formation at a place different from the actual position. In the delay adding method in the related art, correction of wave-front distortion has been a problem for a long time. Accordingly, an aberration correction technology utilizing cross-correlation processing is examined, and similarly in the adaptive beamformer, significantly varying the image quality of an ultrasound image due to heterogeneity inside the medium has become a problem.
As the related art of the beamformer, for example, PTL 1 discloses a technology of a delay adding beamformer which utilizes a fixed null direction.
According to the technology in PTL 1, in order to decrease an artifact due to a side lobe, a main lobe is oriented in a predetermined direction of the target from which information is intended to be acquired so as to generate a first reception signal. Then, the null direction is oriented in the predetermined direction, thereby generating a second reception signal. The second reception signal is an artifact (noise) signal which negligibly includes information regarding a predetermined direction of the target. Therefore, the second reception signal is subtracted from the first reception signal. Thus, a useful signal component in an ultrasound image can be prevented from lacking, and the artifact can be decreased in the first reception signal (particularly, in paragraph 0060 in PTL 1).
PTL 1: JP-A-2010-158374
Since an actual test object medium is heterogenous, a transmitted ultrasound signal is reflected by a focal point and scatters over the periphery of the focal point. Accordingly, target information at the focal point also exists on the periphery of the focal point. Since an adaptive beamformer has sharp directivity, the ultrasound signal existing on the periphery of the focal point cannot be acquired, thereby leading to a problem of deterioration in the image quality of an ultrasound image.
In addition, due to the heterogenous test object medium, there is an occurrence of wave-front distortion resulting from non-uniformity of sound velocity and the like, and the wave-front distortion causes a problem that a reflected wave from the medium on the periphery of the focal point may be mixed into a reception sound wave. Since the reflected wave from the medium on the periphery thereof is a noise signal (correlative noise) which is correlated to a signal from the focal point, it is difficult to be eliminated through a general noise decreasing technology. The correlative noise causes deterioration in the image quality of an ultrasound image.
An object of the present invention is to provide an ultrasound imaging apparatus which can compensate deterioration in the image quality resulting from heterogeneity of a test object medium.
A reception beamformer of an ultrasound imaging apparatus according to the present invention includes a phasing synthesizing unit which performs synthesizing after performing phasing processing of the signal received by the ultrasound element array for each of two or more steering directions and a steering direction instruction unit which instructs the phasing synthesizing unit to perform phasing in the two or more steering directions. The two or more steering directions include at least two directions other than a direction of a reception focal point.
According to the present invention, correlative noise can be decreased by synthesizing a signal which is subjected to phasing processing for two or more steering directions including two directions other than a direction of a reception focal point direction. Therefore, it is possible to realize compensation of deterioration in the image quality of an ultrasound image caused by correlative noise resulting from wave-front distortion and improvement of an S/N ratio.
a) a perspective view illustrating a schematic configuration of an ultrasound imaging apparatus in the first embodiment, and
a) and 23(b) are diagrams illustrating image contrast and null positions which are acquired from different test objects of the present invention.
An ultrasound diagnostic apparatus in an embodiment of the present invention will be described.
The ultrasound diagnostic apparatus in a first embodiment will be described with reference to
As long as each of the two directions 21 and 22 is a direction other than the direction 20 of the reception focal point 10, it is possible to acquire a certain effect of compensating deterioration in the image quality resulting from heterogeneity of a test object medium. In addition, as illustrated in
Hereinafter, an ultrasound imaging apparatus of the first embodiment will be described more specifically.
More descriptions will be given regarding the overall configuration of the ultrasound imaging apparatus with reference to
As shown in
The transmission beamformer 104 generates a signal for a transmission beam. The transmission beam signal is transferred to the ultrasound probe 106 via the transmission/reception separation circuit 107. The ultrasound probe 106 transmits an ultrasound wave from the ultrasound element 105 of the ultrasound element array 101 toward the inside of a body of a test object 100. An echo signal which is reflected inside the body is received by the ultrasound element array 101 of the ultrasound probe 106. The reception signal passes through the transmission/reception separation circuit 107 again, thereby being subjected to phasing computation processing and the like in the reception beamformer 108.
As described above, the reception beamformer 108 includes the phasing synthesizing unit 113 and the steering direction instruction unit 112. The steering direction instruction unit 112 instructs the phasing synthesizing unit 113 to perform phasing in the two or more steering directions. As shown in
The phasing synthesizing unit 113 includes a phasing processing unit 204 and a synthesizing unit 206. The phasing processing unit 204 performs adding after performing phasing of each of the reception signals output by the multiple ultrasound elements 105 of the ultrasound element array 101 with respect to a point 11 of a steering direction 21. Then, phasing-adding is also performed with respect to a point 12 of the steering direction 22 in a similar manner. When there is another steering direction, phasing-adding is performed in a similar manner. Phasing processing of the phasing processing unit 204 may be performed through any processing method. For example, delay adding processing or adaptive phasing processing can be adopted.
The synthesizing unit 206 receives a result of phasing-adding from the phasing processing unit 204 for each of the multiple steering directions 21 and 22, and the like. Then, the synthesizing unit 206 synthesizes the result thereof by performing adding processing and the like, thereby outputting a phasing signal y(n).
An output (the phasing signal y(n)) of the synthesizing unit 206 is transferred to the image processing unit 109, and various types of image processing such as diverse filtering and scan converting are executed, thereby generating an ultrasound image. The ultrasound image is transferred to the image display unit 103 and is displayed therein.
In this manner, according to the present invention, by synthesizing a signal which is subjected to phasing processing, with the two directions 21 and 22 forming predetermined angles θL and θR with respect to the direction 20 of the reception focal point 10 along the ultrasound element array 106 in the array direction, it is possible to acquire the phasing signal y(n) in which noise (correlative noise Nc) indicating a correlationship is decreased, with respect to a signal S from the reception focal point 10 (the point to be received, the reception focal point). Accordingly, the phasing signal y(n) having a significant S/N ratio can be acquired. Thus, the image quality of an ultrasound image can be improved.
Followings are additional descriptions. Since the test object medium is heterogenous, a transmitted ultrasound signal is reflected by the reception focal point 10, and since the reflected ultrasound signal scatters over points on the periphery of the focal point, target information at the reception focal point 10 also exists on the periphery of the reception focal point 10. In addition, wave-front distortion occurs due to non-uniformity of sound velocity and the like resulting from heterogeneity of the test object medium. The wave-front distortion causes a reflected wave from the medium on the periphery of the reception focal point 10 to be mixed into a reception sound wave. The reflected wave from the medium on the periphery thereof is a noise signal (correlative noise) which is correlated to a signal from the reception focal point 10. Therefore, there is a need to reproduce the information of the focal point scattered over the periphery of the focal point and to decrease unnecessary correlative noise from the medium on the periphery of the focal point.
Wave-front distortion occurs not only through a propagation process of a reflected wave but also occurs through a propagation process of a transmission sound wave from when the sound wave is transmitted until the sound wave reaches the reception focal point 10, thereby obviously causing an increase of the aforementioned unnecessary correlative noise.
Therefore, according to the present invention, synthesizing is performed for a signal that is subjected to phasing processing for the two directions 21 and 22 each of which forms predetermined angles θL and θR with respect to the direction 20 of the reception focal point 10 along the ultrasound element array 106 in the array direction. For example, according to the present invention as shown in
In addition, as shown in
Here, since the angle θL and the angle θR are oriented in directions different from each other in a transverse direction, S_left and S_right, and Nc_left and Nc_right are not only different from each other in a positional phase and a temporal phase of a signal but also have components of directionality reverse to each other. In other words, S_left and S_right, and Nc_left and Nc_right are asymmetrical to each other in both time and space. S_left and S_right, and Nc_left and Nc_right are signals which cancel out each other when being added. Therefore, as (S_left+Nc_left+Nu) and (S_right+Nc_right+Nu) are added by the synthesizing unit 206 of the phasing synthesizing unit 113, the component of the correlative noise Nc can be cancelled. In addition, particularly, when the angles θL and θR are set to the null angle θnull, the signals of S_left and S_right included in the output signal of the phasing processing unit 204 occupy a small proportion in the overall signal. Therefore, when synthesizing is performed for a signal which is subjected to phasing processing for the two directions 21 and 22 on the right and left each of which forms predetermined angles θL and θR, an effect of decreasing the correlative noise Nc becomes greater than an effect of decreasing the signals S_left and S_right indicating the information of the reception focal point 10. Therefore, as a result, a proportion contributing to improvement of SNR can be increased.
Moreover, as illustrated in
Meanwhile, the non-correlative noise Nu has no directionality and becomes (1/√n) times by being added n times. Accordingly, the non-correlative noise Nu is uniformly decreased by adding the outputs of the phasing processing unit 204 in the multiple steering directions.
As described above, compared to the original output signal (S+Nc+Nu) of the phasing processing unit 204 in the main lobe direction 20, the signal (S′+Nc′+Nu′)=(⅓) ×{(S+N+R)+(S_left+Nc_left+Nu)+(S_right+Nc_right +Nu) which is obtained by synthesizing the output of the phasing processing unit 204 in the main lobe direction 20 in addition to the directions of null beams 21 and 22 on the right and left increases in SNR (Expression (1)). Here, the outputs of the phasing processing unit 204 for three directions are added while being ⅓ each in weight. However apodization can be performed in weights different from one another.
A phasing processing method for the phasing synthesizing unit 113 may be any method. However, by adopting adaptive phasing processing, the component of the signal S indicating the information of the reception focal point 10 can have increased directivity. In other words, while the desired signal component S in the synthesis result (⅓)×(S+S_left+S_right) of the desired signal component is negligibly decreased, an effect of decreasing the correlative noise Nc can be enhanced further, thereby being preferable.
In this manner, according to the first embodiment, by synthesizing the phasing processing results of the steering directions for at least two directions other than the direction of the reception focal point, and utilizing asymmetrical postulate of correlative noise in time and space, it is possible to decrease correlative noise resulting from wave-front distortion. Moreover, since information of the periphery of the reception focal point 10 can be collected, information which is reflected by the focal point and scatters over the periphery can also be collected. Accordingly, even though the medium of the test object is heterogenous, deterioration in the image quality can be prevented. Thus, it is possible to acquire the ultrasound imaging apparatus which is unlikely to be affected by correlative noise due to wave-front distortion, has high robustness, and has a relatively small processing load.
With reference to
According to the second embodiment, the phasing processing unit 204 inside the phasing synthesizing unit 113 of the first embodiment includes delay adding units 204-1 to 204-P in the number P in order to perform phasing processing through delay adding processing. The number P is the maximum number of the steering directions which can be designated by the steering direction instruction unit 112. In addition, the steering direction instruction unit 112 includes a null angle calculation unit 501 which calculates the null angle based on a condition of the ultrasound element array 101 and an irradiation condition of an ultrasound wave, and the directions of the null angles obtained by performing calculation are designated to the phasing processing unit 204 as the steering direction. The directions of the null angles are set to include at least two angles other than the direction of the reception focal point 10. The delay adding units 204-1 to 204-P of the phasing processing unit 204 are subjected to phasing through delay adding for each of the designated steering directions.
In addition, according to the second embodiment, an active channel setting unit 111a inside the control unit 111 sets active channels 201, 202, 203, and the like at a portion of the finite caliber of the ultrasound element array 101. The reception beamformer 108 performs reception beamforming for a reception signal of the active channel. In other words, elements 105 in the number K within a predetermined range among the ultrasound elements 105 configuring the ultrasound element array 101 which has received a reception echo with respect to one ultrasound transmission beam is set to be the active channel 201, and the reception beamformer 108 generates one item of image data (a raster: a phasing output y(n)) in the ultrasound propagation direction by using the reception signal of the active channel 201. As shown in
The null angle calculation unit 501 calculates the null angle for the set active channel 201 and the like.
Hereinafter, descriptions will be given more specifically with reference to the flow in
First, the control unit 111 generates a control signal which includes information indicating a condition of the probe, an irradiation condition of an ultrasound wave, and a selection condition of the null angle corresponding to the current imaging condition, thereby outputting the control signal to the null angle calculation unit 501 (Step 61).
The null angle calculation unit 501 receives an output of the control unit 111 (Step 62), and the null angle is calculated by using the following Expression (3) and the like while having the information included in the control signal as a parameter (Step 63). Expression (3) calculates the null angle only for transmission. However, the null angle can be calculated for reception as well by performing calculation in a similar manner. In addition, it is possible to calculate the superimposed transmission/reception null angle.
In calculation of the null angle, the primary null angle can be confirmed based on a frequency of an ultrasound wave, an imaging condition such as a frequency waveform, the number of the elements (the number of the active channel elements) of the ultrasound element 105 included in the currently used active channels 201 to 203, the caliber for transmission and reception, the channel intervals (the intervals in the ultrasound element 105), and characteristics of a frequency. In other words, the null angle can be calculated by using the control signal transmitted from the control unit 111.
For example, when the ultrasound element 105 is a square piston-type ultrasound element, directivity of an ultrasound beam transceived by the active channel 201 of the ultrasound element array 101 can be calculated through the following Expression (2). D(θ) is intensity of an ultrasound beam at a steering angle θ.
In Expression (2), θ represents the steering angle, a represents the caliber of the active channels (201 to 203), and k represents the number of waves. The value of k is obtained by dividing a sound velocity c by a frequency f of an ultrasound wave, thereby being represented by k=c/f.
angle θnull having a value of zero in Expression (2) is the null angle in this case. The null angle can be analytically obtained through Expression (3).
As seen in Expression (3), there are multiple null angles θnull instead of being one. As shown in
The effect is maximized when n=±1. Therefore, it may be configured that only a processing unit of n=±1 included and two null angles of n=±1 must be used without performing null selection.
The null angle calculation unit 501 outputs the directions of the null angles which are selected from the null angles obtained by performing calculation in accordance with the selection condition of the null angle, with respect to the phasing processing unit 204 as the steering directions. As a form of the information designating the steering direction, for example, the information may be physical angle information such as θ=1×10° and −10°. In addition, a control signal for activating (turn ON) the delay adding units 204-1 to 204-P of the phasing processing unit 204 and a signal which shares the aforementioned angle information may be directly transferred to the delay adding units 204-1 to 204-n. The aforementioned method exhibits a practically preferable embodiment. For example, θL=10° is designated for the delay adding unit 204-1, and delay adding processing is instructed. Then, θR=−10° is designated for the delay adding unit 204-n, and delay adding processing is instructed.
The delay adding units 204-1 to 204-P of the phasing processing unit 204 generates a phasing signal in a case where the directions of the null angles designated from the null angle calculation unit 501 is set to be the steering direction, thereby outputting the generated signal (Step 65). In this case, it is desirable that a phasing signal is generated in the front direction (θ=0°: the direction of the reception focal point 20) as well. In other words, when the null angles are θ=10° and −10°, in the delay adding units 204-1, 204-2, and 204-3, delay adding processing is performed for three directions, that is, the front directions of θ=0°, θ=10°, and θ=−10°, and the results thereof are output to the synthesizing unit 206.
The synthesizing unit 206 receives inputs of the phasing signals for the null direction and the front direction from the delay adding units 204-1 to 204-P (Step 66). The synthesizing unit 206 performs synthesizing by adding the inputs, thereby outputting the added result (Step 67). When performing adding, for example, simple adding may be performed for each of three directions, that is, the front directions of θ=0°, θ=10°, and θ=−10° (the adding result is preferably multiplied by ⅓). Otherwise, for each steering direction, adding may be performed after multiplying the phasing signal by a predetermined weight. For example, it is possible to increase the proportion of the phasing signal for the front direction. In contrast, the synthesis result can be adjusted by multiplying the phasing signal for the null direction by a significant weight.
The synthesis signal computed by the synthesizing unit 206 is output to the image processing unit 109. The image processing unit 109 receives the synthesis signal (Step 68), and performs processing similar to that for the phasing result in the related art, thereby generating one item of image data (a raster). Then, while being gradually shifted in position, each of the active channels 201, 202, and 203 acquires the raster. The rasters are arrayed so as to be processed into an ultrasound image. The ultrasound image is displayed by the image display unit 103. In this manner, the flow of the second embodiment according to the present invention is completed.
In the above-described flow in
In addition, in the above described configuration, calculation is performed for the null angle every time a change is made for the condition indicated by the control signal from the control unit 111. However, without being limited to the configuration, it is possible to arrange a memory 502 in the null angle calculation unit 501, and to store a table (the null angle LUT (LUT: Look Up Table)) in which the calculation results of the null angles in a set correspond to conditions (a condition of the probe, an irradiation condition of an ultrasound wave, and a selection condition of the null angle) indicated by the control signal, in advance. In this manner, the null angle calculation unit 501 can read out the null angles corresponding to the conditions indicated by the control signal from the null angle LUT so as to set the null angles for the delay adding units 204-1 to 204-P. Therefore, the computation amount can be decreased.
Regarding the null angle LUT, the values prepared by obtaining the null angles in advance for all the conditions that can be set (a condition of the probe, an irradiation condition of an ultrasound wave, and a selection condition of the null angle) can be stored in the memory 502. However, without being limited thereto, the null angle LUT for only a portion of the conditions may be prepared and stored in the memory 502. In a case of a condition which is not included in the null angle LUT, the null angle calculation unit 501 may be configured to obtain the null angle by performing calculation.
In addition, it is also possible to have a configuration in which the null angle obtained by the null angle calculation unit 501 by performing calculation is caused to correspond to the condition indicated by the control signal at the moment and the result thereof is stored in the memory 502, thereby successively generating the null angle LUT.
As illustrated in
As described above, in the second embodiment, by acquiring the phasing results to be obtained through delay adding processing of the steering directions for at least the two directions other than the direction of the reception focal point, synthesizing the acquired results, and utilizing asymmetrical postulate of correlative noise in time and space, it is possible to decrease correlative noise resulting from wave-front distortion. Moreover, information which is reflected by the focal point and scatters over the periphery can also be collected. Accordingly, even though the medium of the test object is heterogenous, deterioration in the image quality can be prevented. Thus, it is possible to acquire the ultrasound imaging apparatus which is unlikely to be affected by correlative noise due to wave-front distortion, has high robustness, and has a relatively small processing load.
Other configurations of the ultrasound imaging apparatus in the second embodiment are similar to those in the first embodiment, thereby omitting the descriptions.
As a third embodiment, the ultrasound imaging apparatus will be described with reference to
As shown in
The output signal of a delay adding unit 204 in the present invention is the same as multi-look (multi-directional) reception data in known aperture synthesis processing. Accordingly, by applying the present invention when performing synthesizing of multi-look reception data among multiple transmissions in the known aperture synthesis processing, phasing synthesis processing for the two or more steering directions of the present invention can be realized through aperture synthesis processing.
Specifically, items of data in the number P which are generated for each ultrasound wave transmission by the delay adding units 204-1 to 204-P are successively saved in the frame memory unit 701 in
The data subjected to adding processing by the frame adding unit 702 is data similar to what is synthesized by the synthesizing unit 206 in the second embodiment (the data after being subjected to delay adding in multiple directions). In other words, the two or more steering directions from the steering direction instruction unit 112 are input to the frame adding unit 702 as synthesis angle information. The frame adding unit 702 reads out data regarding the designated two or more steering directions from the frame memory 701, thereby performing adding processing. Accordingly, aperture synthesis processing and data adding (adding of the two or more steering directions) in the null direction can be performed at the same time.
The ultrasound imaging apparatus in the fourth embodiment will be described with reference to
In the fourth embodiment, the steering direction instruction unit 112 includes a null angle detection unit 207 and detects the null angle by using the phasing signal of the actual reception signal of the ultrasound element array 101. In order to enable the detection, the phasing synthesizing unit 113 includes the delay adding units 204-1 to 204-Q in the number Q in the phasing processing unit 204, and a memory unit 205 which stores the phasing outputs of the delay adding units 204-1 to 204-Q in the number Q.
Moreover, the null angle detection unit 207 includes the null angle calculation unit 501 described in the second embodiment and is configured to perform scanning of the null angle on the periphery of the null angle obtained by performing calculation so as to detect the null angle.
Other configurations are similar to those in the second embodiment, thereby omitting the descriptions.
Hereinafter, an operation of the reception beamformer 108 in the fourth embodiment will be described with reference to the flow chart in
The second embodiment is configured to obtain the null angle from a condition of the probe, an irradiation condition of an ultrasound wave, and the like by performing calculation. Therefore, the null angle is decided depending on the condition of the apparatus side such as setting of the apparatus or an irradiation condition of an ultrasound wave. However, directivity of the actual reception signal is affected by heterogenous sound wave propagation in a living body so that directivity of a reception beam varies for each sample of the reception signal. Therefore, in the fourth embodiment, the null angle detection unit 207 detects the null angle for each sample of the reception signal based on the actual reception signal.
The control unit 111 generates and outputs control signals which indicate a condition of the probe, an irradiation condition of an ultrasound wave, and a selection condition of the null angle in accordance with the imaging condition at the moment, similarly to the second embodiment (Step 61). The null angle calculation unit 501 receives the control signals and obtains the null angle by performing calculation (Steps 62 and 63). The details for Steps 61 to 63 are the same as those described in the second embodiment. As the null angle calculation unit 501 obtains the null angle by performing calculation, the null angle detection unit 207 requests the memory unit 205 for data of the phasing signal.
Meanwhile, the delay adding units 204-1 to 204-Q of the phasing synthesizing unit perform delay adding processing for each of the steering directions in the number Q which are different from one another by a predetermined angle, within the preset angle range from θmin to θmax, thereby outputting the phasing signals (Step 92). The memory unit 205 stores the phasing signal for each steering angle within the angle range from θmin to θmax (Step 93). Physically, for example, the angle range from θmin to θmax is the steering angle from −60° to 60°.
In accordance with the request from the null angle detection unit 207 in Step 91, the memory unit 205 outputs the phasing signal for the steering directions in the number Q within the angle range from θmin to θmax to the null angle detection unit 207 (Step 94).
The null angle detection unit 207 receives the phasing signal for the steering directions in the number Q within the angle range from θmin to θmax, the phasing signal within a predetermined angle range θ1 to θ2 is scanned out of the angle range from θmin to θmax, and extracts the minimum (minimal) phasing outputs as the null angle θnull_scan, thereby outputting the extracted results (Steps 95 and 96). This steps (θ1 and θ2) are repeated as many the number of the sets.
Here, the angle range from θ1 to θ2 is an angle range having the preset extent including the null angle which is obtained by the null angle calculation unit 501 through calculation. The angle range from θ1 to θ2 is the angle range which is assumed to have the actual null angle between the angles thereof. In other words, the set of (θ1 and θ2) is the angle range having the preset extent which becomes θ1<θnull <θ2 with respect to θnull obtained by the null angle calculation unit 501. In addition, the number of sets of (θ1 and θ2) is the same number as the number of the null angles θnull which has decided depending on a null selection condition included in the control signal. Regarding this, it is the same as that described in the second embodiment.
By performing the processing, the null angle detection unit 207 detects the actual null angle θnull_scan in the same number as the number based on a selection condition of the null angle.
The null angle detection unit 207 instructs the memory unit 205 regarding the detected null angle θnull_scan. The memory unit 205 transmits the null angle information used in the delay adding unit 204. The transmission form of the detected null angle θnull_scan is the same as that in the second embodiment. For example, the null angle is transmitted as physical angle information such as θ=10° and −10°.
The memory unit 205 selects the phasing signal for the steering directions of the null angle θnull_scan and the front direction (θ=0°) which are transmitted from the null angle detection unit 207, thereby outputting the selected signal to the synthesizing unit 206. In other words, when the null angle is θ=10° and −10°, the memory unit 205 outputs the phasing signals for three directions, that is, the front directions of θ=0°, θ=10°, and θ=−10°. In this case, the direction of θ=0° is not necessarily included.
The synthesizing unit 206 receives the phasing signal of the steering direction of the null angle θnull_scan and the front direction (θ=0° from the memory unit 205 (Step 97) and performs synthesizing by performing adding, thereby outputting the result to the image processing unit 109 (Step 98). When performing adding, simple adding may be performed, or adding may be performed after multiplying each phasing signal by a predetermined weight.
The image processing unit 109 receives the phasing signal which is output by the synthesizing unit 206 and is synthesized. By using the received signal, similar to that in the second embodiment, an ultrasound image is generated, and the image is displayed by the image display unit 103 (Step 99). In this manner, the flow of processing in the fourth embodiment of the present invention is completed.
Regarding one sample of the reception signal, if outputting of the synthesis signal performed by the synthesizing unit 206 is completed in Step 98, the procedure returns to each of Steps 91 and 92, thereby performing processing for the reception signal of the succeeding sample.
As described above, in the fourth embodiment, since the null angle having the directivity of the actual reception signal can be detected, it is possible to set the null angle in accordance with variation of the null angle caused by an influence of heterogenous sound wave propagation in a living body. Accordingly, even though the medium of the test object is heterogenous, deterioration in the image quality can be prevented, and it is unlikely to be affected by correlative noise resulting from wave-front distortion. Thus, it is possible to acquire the ultrasound imaging apparatus having high robustness.
In the present embodiment, the null angle calculation unit may obtain the null angle by performing calculation, and the null angle detection unit 207 may scan the null angle only for a predetermined angle range θ1 to θ2 including the null angle obtained by performing calculation. Accordingly, a computation amount of the null angle detection unit 207 can be decreased. However, the present invention is not limited to the aforementioned configuration. The null angle detection unit 207 need not include the null angle calculation unit 501 and may scan the null angle for the entire angle range.
In addition, the null angle calculation unit 501 may be configured to include the memory 502 as described in the second embodiment and obtain the null angle with reference to the null angle LUT in the memory 502.
The ultrasound imaging apparatus in a fifth embodiment will be described with reference to
The ultrasound imaging apparatus in the fifth embodiment is similar to that in the second embodiment. However, the point where the phasing processing unit 204 performs adaptive phasing processing is different from that in the second embodiment. In other words, as illustrated in
As shown in
In this manner, in the fifth embodiment, since the reception beamformer 108 performs adaptive phasing processing, it is possible to acquire an effect in which the correlative noise Nu can be more effectively decreased while the desired signal component S is negligibly decreased.
Hereinafter, adaptive phasing processing will be described specifically. When the number of the ultrasound elements 105 (the number of channels) configuring the active channel 201 is k, post-delaying reception data of the delay circuit 511 at a certain snapshot time n can be represented by a vector x(n) on the left side in the following Expression (4), by using an output xk(n) of the ultrasound element 105.
[Expression 4]
x(n)=[x1(n),x2(n), . . . ,xL(n)]T (4)
The post-delaying reception data x(n) is input to each of the adaptive processing units 512-1 to 512-P. The adaptive processing units 512-1 to 512-P generate an adaptive weight vector wp(n)=[wp
After synthesizing the adaptive weight vectors w1(n), w2(n) to wp(n), the adaptive synthesizing unit 516 performs inner product computation (apodization computation) with the post-delaying reception data x(n), thereby acquiring the adaptive beamformer output y(n) at a certain snapshot time n.
As shown in
In the matrix computation unit 300, the post-delaying reception data vector x(n) obtains a spatial covariance matrix R(n) through the following Expression (5). When obtaining the matrix R(n), the vector x(n) may be used as the actual signal as it is, or a vector which is converted into data in the complex number by executing Hilbert transformation, baseband modulation, or the like. Here, as a more general form, descriptions will be given by exemplifying a case where the vector x(n) is converted into a complex data ξ(n) so as to obtain the spatial covariance matrix R(n) as shown in Expression (5). The matrix R(n) in Expression (5) is obtained by taking an ensemble average of the product between the complex vector ξ(n) which is represented by Expression (6) and the (complex) conjugate transposition vector ξH(n) thereof.
In Expression (5), an ensemble average number N can be uniformly averaged as shown on the far right side in Expression (5), as the points in total N=2S+1 including the samples in the number S in front and rear of the vector ξ(n) in the target snapshot. Moreover, as the method of averaging the temporal direction, in addition thereto, it is possible to adopt a method in which after each sample in the temporal direction is multiplied by an arbitrary weight such as a trapezoidal weight, adding average is obtained. The spatial covariance matrix R(n) output from the matrix computation unit 300 is successively input to the adaptive beam steering unit 301.
In the adaptive beam steering unit 301 which receives the spatial covariance matrix R(n), the weight vector wp(n) is calculated by using the MVDR method. A steering vector ap designated by the null angle calculation unit 501 is represented by Expression (7).
In Expression (7) shown above, the reference sign p represents the number of the steering vectors and is an integer satisfying 0<p<P+1 as the total number P. As shown in Expression (7), the steering vector ap is a directional vector which has the vector element (0 to (K−1)) in the number K equal to the number of the active channels. The steering vector ap is represented by a function of a reception frequency fp, and an angle (hereinafter, referred to as a steering angle) (θp, φp) formed by the normal line vector direction on the surface of the ultrasound element 105 and the steering vector.
The reference sign θp represents an open angle from the normal-line vector, and the reference sign φp represents a swiveling angle from the array direction of the ultrasound element 105. When the arrangement of the ultrasound elements 105 of the ultrasound probe 106 is a one-dimensional straight-line (linear) array, the steering angle is represented by the last expression in Expression (7). Here, the reference sign λp represents a wavelength of a sound wave corresponding to the frequency fp, and the reference sign d represents the interval (the pitch of the elements) between the element centers of the ultrasound element 105.
Regarding the direction of the above-described steering vector ap, the weight vector wp(n) obtained by the MVDR method is calculated through Expression (8) herein. Accordingly, by computing the weight vector wp(n) for each of the steering vectors ap different from one another by the adaptive processing units 512-1 to 512-P, it is possible to acquire the adaptive weight vector w1(n) to wp(n) as many the number of P of the steering vector ap.
In Expression (8), the reference sign R(n) represents a spatial covariance matrix at a certain snapshot n in the temporal direction generated through Expression (5), and the superscript −1 represents an inverse matrix.
The reception focal point 10 is on the central axis of the active channel 201, that is, when the number of the active channels is an even number, the reception focal point 10 is positioned on the normal line (the normal line orthogonal to the active channel surface) passing through the center point between the elements on the K/2nd and (K+2)/2nd among the ultrasound elements 105 in the number K configuring the active channels 201. In addition, when the number of the active channels is an odd number, the reception focal point 10 is positioned on the normal line (the normal line orthogonal to the active channel surface) passing through the center of the ultrasound element 105 on (K+1)/2nd among the ultrasound elements 105 in the number K configuring the active channels 201.
Since delaying processing is executed in a delay circuit 204, for example, the reception data vector x(n) at a certain time n in linear scanning becomes data in which wave-fronts are aligned in a direction of θ=0° that is a direction of the normal line vector. Accordingly, in Expression (7), when θp=0°, the steering vector ap is ap=[1, 1, and so forth 1]T, thereby obtaining the adaptive weight vector w(n) regarding the focal point direction. The adaptive weight vector w(n) regarding the focal point direction coincides with the adaptive weight vector w(n) obtained by the MVDR method in the related art.
The adaptive weight vectors w1(n) to wp(n) in the number P output from the adaptive beam steering unit 301 are input to the adaptive synthesizing unit 516. As shown in
In the previous stage where the synthesis weight vector wsum(n) is obtained, a fixed apodization multiplying unit 305 may be arranged so as perform multiplication for the fixed apodization with respect to each of the weight vectors wp(n). For example, it is possible to perform multiplication for a fixed apodization bp applied with a distribution in which the value of the adaptive weight vector wp(n) having the direction of θ=0° is increased and the values of other directions are decreased. This is realized by performing computation through Expression (10).
For example, in the configuration in
Subsequently, the synthesis weight vector wsum(n) is input to an apodization computation unit 307 inside the adaptive synthesizing unit 516. The apodization computation unit 307 includes a multiplying unit 3071 and performs an inner product computation as shown in Expression (11) between the synthesis weight vector wsum(n) and the post-delaying reception data vector x(n) which is transmitted through a bypass line 207 from the delay circuit 204, thereby acquiring the phasing output y(n). Specifically, the multiplying unit 3071 calculates the product between the channel elements (1 to K) of each of the weight vectors and the post-delaying reception data vectors. An adding unit 3072 calculates the sum total of the products for the elements in the number k acquired by the multiplying unit 3071 of the apodization computation unit 307, as the final output (the phasing output y(n), a scalar value) of the adaptive synthesizing unit 516.
[Expression 11]
y(n)=wsumH(n)×(n) (11)
The phasing output y(n) of the apodization computation unit 307 is output to the image processing unit 109. The above-described processing is repeated from the first sample n=1 to the last sample n=N among the N samples configuring one raster. Then, the apodization computation unit 307 outputs the sequential phasing outputs Y(1), Y(2) to Y(N) to the sequential image processing unit 109. The phasing output y(n) of one raster which can acquired through Expression (11) is acquired individually while shifting the active channel 201 to the active channels 202 and 203 on the ultrasound element array 101. In the image processing unit 109, all the rasters are arrayed by a scan converter adapted for the scanning method of the ultrasound probe 106, thereby generating a two-dimensional image. In addition, various types of back-end image processing such as diverse filtering and computation of measurement application are performed. Eventually, the image display unit 103 displays an ultrasound image or a computation result of the measurement application.
In this manner, in the fifth embodiment, by using adaptive phasing processing, while the desired signal component S is prevented from decreasing, it is possible to acquire an effect in which the correlative noise Nu can be more effectively decreased.
In the fifth embodiment, as another technique for computing the algorithm in the above-described matrix computation unit 300, it is possible to perform spatial averaging computation by using a sub-array matrix. A sub-array matrix R̂SUB1 in spatial averaging computation is represented by the product of the subspace vector ξ̂1(n) (Expression (13)) as shown in Expression (12). The subspace vector ξ̂1(n) is a vector which is obtained by partially eliminating the component (corresponding to the element L) out of the post-delaying reception data (represented by the general complex signal vector ξ(n) in this case, but the actual signal vector x(n) may be adopted in place thereof) with respect to the active channels in the number K. Accordingly, the total number of the subspace vector becomes the number K−L+1(0 <1(small letter 1)<K−L+1).
[Expression 12]
R̂
SUB1(n)=ξ̂l(n)ξ̂lH(n) (12)
[Expression 13]
ξ̂l(n)=[ξl(n), ξl+1(n), . . . , ξl+L−1(n)]T (13)
When the main diagonal component of the sub-array matrix is caused to coincide with the main diagonal component of the spatial covariance matrix R(n) so as to be shifted by one sample, it becomes spatial averaging processing for the sub-array matrix in the number K−L+1, thereby acquiring the sub-array spatial covariance matrix R̂(n) in Expression (14). When computing the sub-array spatial covariance matrix R̂ (n) by the adaptive beam steering unit 301, the sub-array spatial covariance matrix R̂(n) can replace R(n) in the above-described Expression (8), and thus, it is possible to compute the weight vector wp(n) as shown in Expression (15). In this case, the output of the matrix computation unit 300 has a size of L×L, and the number of elements configuring the weight vector wp(n) becomes L.
In addition, as another example of the spatial averaging method, a frontward-rearward spatial averaging method may be adopted. In this case, a rearward sub-array matrix R to SUB1(n) can be obtained by the product of a rearward the subspace vector ξto 1(n) as shown in Expression (16). The rearward subspace vector is represented by Expression (17). In addition, as shown in Expression (18), by performing computation similarly to the frontward spatial averaging, the rearward sub-array spatial matrix R to (n) can be calculated as shown in Expression (18). Subsequently, as shown in Expression (19), by performing arithmetic averaging with the sub-array spatial matrix R̂(n) obtained in the above-described frontward spatial averaging, it is possible to eventually obtain a frontward/rearward sub-array spatial covariance matrix RFB(n). Similar to frontward spatial averaging, by replacing R(n) in Expression (8) with the frontward/rearward sub-array spatial covariance matrix RFB(n) when computing the adaptive beam steering unit 301, the weight vector wp(n) can be computed as shown in Expression (20). In this case as well, the output of the matrix computation unit 300 has a size of L×L, and the number of elements configuring the weight vector wp(n) becomes L.
The weight vector wp in multiple numbers (the number P) calculated through Expression (15) or Expression (20) by using the spatial averaging method is transmitted to the adaptive synthesizing unit 516, similarly to a case where the spatial averaging method is not used.
Here, when the spatial averaging method is used, the number of elements of the weight vector wp is L. Accordingly, in order to eventually perform computation in an inner product computation unit 307, a block is additionally required so as to perform computation (Expression (21)) for producing from ξ(n) having the components in the number K to the vector g(n) having the components in the number L. For example, as illustrated in
The phasing output y(n) when the spatial averaging method is used can be represented by Expression (22).
In this manner, in the matrix computation unit 300, by performing spatial averaging processing of the sub-array matrix, it is possible to acquire an effect of preventing correlative noise included in an ultrasound wave reception signal. Accordingly, by combining spatial averaging processing of the sub-array matrix in the configuration in which multiple adaptive weight vectors on the periphery of the focal point is synthesized, it is possible to acquire an ultrasound image which is negligibly affected by noise.
As described above, in the fifth embodiment, MVDR is exemplified as an example of a technique for adaptive beam forming. However, the algorithm applied to the adaptive beam steering unit 301 is acceptable as long as the algorithm used the spatial covariance matrix calculated by the matrix computation unit 300. Moreover, any types of other techniques such as the MMSE method; the APES method; the eigen space-MV method (ESMV, EIBMV) utilizing the spatial covariance matrix, the eigen value, and the eigen vector thereof; the ESPRIT method, and the MUSIC method, can be used.
The ultrasound imaging apparatus in a sixth embodiment will be described with reference to
Specifically, as illustrated in
Similarly to the fourth embodiment, the adaptive processing units 612-1 to 612-Q obtain the phasing outputs for the steering directions in the number Q within the angle range from θmin to θmax. As illustrated in
Operations of the memory unit 205 and the null angle detection unit 207 are similar to those in the fourth embodiment, thereby omitting the descriptions.
In the sixth embodiment, since the phasing signals obtained through adaptive phasing processing for the steering direction are synthesized, it is possible to enhance the effect to decrease correlative noise resulting from wave-front distortion. In addition, since the null angle is detected by using the phasing signal obtained through adaptive phasing processing, there is an effect of achieving high detection accuracy of the null angle as well.
The ultrasound imaging apparatus in a seventh embodiment will be described with reference to
In the present embodiment, the null angle is detected by the phasing signal obtained through delay adding processing, and adaptive phasing processing is performed only for the steering direction designated by the null angle detection unit 207. Therefore, it is possible to synthesize the phasing signals obtained through adaptive phasing processing for the steering direction with a computation amount less than that in the sixth embodiment. Thus, it is possible to enhance the effect of decreasing correlative noise resulting from wave-front distortion with a small computation amount.
The adaptive phasing unit 612 and the synthesizing unit 206 in
The ultrasound imaging apparatus in an eighth embodiment will be described with reference to
The eighth embodiment is configured to perform aperture synthesizing in the reception beamformer 108 of the fifth embodiment in
With reference to
As shown in
The frame memory unit 701 successively saves the phasing signals in the number R which are subjected to adaptive phasing processing and are generated by the phasing processing units 204 in the number R for each of the ultrasound wave transmissions. The frame adding unit 702 receives the two or more steering directions from the steering direction instruction unit 112 as synthesis angle information, and reads out the data from the frame memory 701 regarding the designated two or more steering directions, thereby performing adding processing. Accordingly, aperture synthesis processing and data adding in the null direction (adding of the two or more steering directions) can be performed at the same time.
In the eighth embodiment, since the phasing signals obtained through adaptive phasing processing for the steering direction are synthesized, compared to the third embodiment, the effect of decreasing correlative noise resulting from wave-front distortion is significant. In addition, since the null angle can be calculated by using the phasing signal obtained through adaptive phasing processing, there is an effect of achieving high detection accuracy of the null angle as well.
The ultrasound imaging apparatus in a ninth embodiment will be described with reference to
The ninth embodiment provides another embodiment of the adaptive synthesizing unit different from those in the fifth to eighth embodiments. The inside of the adaptive synthesizing unit in the fifth to eighth embodiments may have the configuration as that in the fifth embodiment illustrated in
In the fifth embodiment illustrated in
According to such a configuration, in the inner product computation unit 307, the bypass input 2071 is subjected to phasing-adding by using the weight vector w1(n). The bypass input 2072 is subjected to phasing-adding by using the weight vector w2(n), and the bypass input 2073 is subjected to phasing-adding by using the weight vector w3(n), thereby calculating the pre-synthesizing multiple phasing outputs y1(n), y2(n), and y3(n) corresponding to each of the steering vectors (Expression (23) below). In the following stage thereof, the output synthesizing unit 500 obtains the arithmetic average of the phasing outputs y1(n), y2(n), and y3(n) corresponding to each of the steering vectors, thereby acquiring the synthesizing-phasing output ysum(n) as shown in Expression (24).
As illustrated in
As described above, the fifth embodiment and the eighth embodiment are different from each other in the embodiment whether inner product computation is performed after the weights are synthesized or the phasing output is synthesized after inner product computation. Meanwhile, the final outputs are identical to each other due to linear computation. Therefore, the values of y(n) in Expression (22) which is the last phasing output is the fifth embodiment and the last synthesizing-phasing output ysum(n) of the eighth embodiment are the same.
In addition, the mode switching unit of the console may be formed so as to be able to switch the mode in accordance with switching in different probes, an imaging condition, and an imaging sequence. For example, it is possible to adopt a switching unit which can switch the mode so as to be able to generate an image by applying a set of the steering directions different from each other, depending on a probe such as a linear probe, a convex probe, a sector probe, a two-dimensional array probe, and a mechanical type 3D imaging probe. In addition, for example, it is possible to adopt a switching unit which can generate an image by applying a set of the steering directions different from each other depending on an imaging condition and an imaging sequence such as transmission/reception focus, a transmission/reception frequency, a frame rate, the number of parallel beamformers, tissue harmonics, and imaging contrast.
In addition, it is possible to provide a configuration in which a practitioner can set parameters of imaging conditions by using the handy operation unit 1103.
As it is clearly shown in
In
According to
In
According to
Lastly, with reference to
As seen in the line profiles 1601 and 1602 in
100 test object
101 ultrasound element array
102 apparatus main body
103 image display unit
104 transmission beamformer
106 ultrasound probe
107 transmission/reception separation circuit (T/R)
108 reception beamformer
109 image processing unit
110 console
111 control unit
204 phasing processing unit
205 memory unit
206 peripheral information synthesizing unit
207 null angle detection unit
301 adaptive beam steering unit
501 null angle calculation unit
701 channel memory
702 frame adding unit
1001, 1002 knob portion
Number | Date | Country | Kind |
---|---|---|---|
2013-003863 | Jan 2013 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2014/050344 | 1/10/2014 | WO | 00 |