The present invention relates to an ultrasound imaging technique for imaging the inside of a test subject, through the use of ultrasound waves.
Ultrasound imaging is a technique for non-invasively creating an image of the inside of a test subject such as a human body, through the use of ultrasound waves (being a sound wave not intended for hearing, and generally a high-frequency sound wave having 20 kHz or higher). By way of example, a medical ultrasound imaging apparatus will be briefly explained. An ultrasound probe transmits the ultrasound waves to the inside of a patient body, and receives echo signals reflected from the inside of the patient. The received signals are subjected to signal processing in one or both of the ultrasound probe and the main unit of the ultrasound imaging apparatus, and thereafter transferred to a monitor, and then, an ultrasound image is displayed thereon. More specifically, for example, a transmit beamformer in the main unit of the ultrasound imaging apparatus generates signals of a transmission beam, allowing the signals to pass through the transmit-receive separation circuit (T/R), and thereafter transfers the signals to the ultrasound probe. The ultrasound probe sends out ultrasound waves. After receiving the echo signals from the internal body, the ultrasound probe transmits the signals to the main unit of the imaging apparatus. In the main unit of the imaging apparatus, the received signals pass through the transmit-receive separation circuit again, being subjected to beamforming process in the receive beamformer, and then those signals are transmitted to an image processor. The image processor executes various image processing through the use of various filters, a scan converter, and the like. Finally, the monitor displays an ultrasound image.
As described above, a general ultrasound diagnostic apparatus is made up of three techniques; transmit beamforming, receive beamforming, and a backend image processing. Particularly, since the beamformers for transmitting and receiving perform signal processing at an RF (high-frequency) level, algorithms and implementation architecture in the beamformers decide a basic image quality of the ultrasound image. Therefore, the beamformers serve as major parts of the apparatus.
The receive beamformer assigns a delay time to each received signal (received data) in plural elements that constitute the ultrasound probe, the delay time distributing an amount of delay in a concave form, in association with the relations between a focal position and the element positions, and after virtually achieving focus (focused) at a certain point in space, the received signal data items are summed up. This method is referred to as a beamforming according to a delay-and-sum method. In this delay-and-sum method, the received data items that are received by the plural elements in the ultrasound diagnostic apparatus are multiplied by a fixed weight vector stored in the diagnostic apparatus, then being weighted and summed up. This process is also performed in the transmit beamformer in a similar manner, not only in the receive beamformer.
On the other hand, as a basic problem of the ultrasound imaging apparatus, it is known that lateral resolution is subject to constraints. Since transmitting and receiving of the ultrasound waves are performed by an array having a finite opening size, there is an impact of diffraction on the edge of the opening. If an infinitely long array is prepared, there is a possibility that the resolution is enhanced infinitely in the same manner as in the depth direction. In actual, however, a physical restriction in designing the apparatus, i.e., the length of the array, has hampered the enhancement of the lateral resolution. In recent years, it is attempted that the aforementioned fixed weight vector used for delaying, upon the delay-and-sum by the beamformer, is made to vary adaptively for the time-series transmit-receive data items, one by one, thereby obtaining an ultrasound image of higher definition, and this attempt is coming to attention. Accordingly, there is a possibility that this brings a marked improvement in the lateral resolution, being one of essential problems in the beamforming technique.
Particularly in recent years, an improved lateral resolution technique is being reported which is achieved by applying an adaptive signal processing technique including the MVDR method (Minimum Variance Distortionless Response; Capon method) that has been developed in the field of mobile communication, to the beamformer of the received data (Non Patent Documents 1 to 6). The technique described in those adaptive methods is implemented by adaptively varying complex components of the weight vector used for the delay-and-sum, on the basis of the correlative matrix of the received data. In other words, conventionally, the weight vector has been a fixed value, but in the adaptive method, the weight vector is obtained according to an operation with the use of the received signal, for each sample point in the time direction of the received signal, and the received signal is multiplied by thus obtained weight vector.
In the adaptive signal processing, similar to the conventional delay-and-sum method, distortion (aberration) of the ultrasound wavefront presents a significant challenge, the distortion being caused by inhomogeneous sound-velocity distribution in a medium and scattering from micro scatterers. Also in the adaptive signal processing, the apparatus sets the focus of the receive beamformer, assuming that the sound-velocity is constant and the medium is homogeneous. Therefore, when distortion exists in sound-wave propagation, there is a problem that an image appears blurred, or an image is formed on a position different from the actual position. In the conventional delay-and-sum method, correction of the wavefront distortion (wave aberration) is a long standing problem, and aberration correction technique utilizing cross-correlation processing has been studied. Also in the adaptive beamformer, it is a challenge to significantly change the image quality of an ultrasound image caused by inhomogeneity of the medium.
Following prior arts are known as describing techniques of a beamformer relating to the present invention; the Patent documents 1 to 3 regarding delay time adjustment using plural beams, the Patent document 4 regarding combined delay outputs, and the Patent document 5 regarding a technique for correcting wavefront distortion.
It is known that in an ultrasound diagnostic apparatus using an adaptive beamformer, beams are highly directive. However, actually, the medium of a test subject is inhomogeneous, and thus ultrasound signals being transmitted are reflected on a focus and scattered around the focus, resulting in that information regarding the object on the focus also exists around the focus. Beams with high directivity from the adaptive beamformer are incapable of acquiring ultrasound signals that exist around the focus, and this may cause a problem that an image quality is deteriorated. Furthermore, since the medium of the test subject is inhomogeneous, there may occur wavefront distortion due to the sound-velocity inhomogeneity, or the like. The wavefront distortion may cause reflected waves from the medium surrounding the focus to mix into the received sound wave. The reflected waves from the surrounding medium may become a noise signal correlated with the information from the focus, and there is a problem that this correlative noise may deteriorate the image quality.
Therefore, there is required a technique which allows the focus information scattered around the focus to be reproduced, in compatible with reducing unnecessary correlative noise from the medium surrounding the focus. In other words, an adaptive beamforming technique is demanded, being able to reduce the number of the ultrasound signals failing to be picked up, after scattered around the focus, highly robust against the correlative noise caused by the wavefront distortion, and provided with an image-quality compensation method with relatively small processing load.
An object of the present invention is to provide an ultrasound imaging apparatus that is capable of compensating for the image-quality deterioration caused by inhomogeneity of the test subject medium.
The ultrasound imaging apparatus of the present invention is provided with plural elements configured to receive ultrasound signals from a test subject, a delay unit configured to delay each of the received signals on the plural elements, in association with a position of a predetermined receive focus, and generate post-delay received signals, a peripheral information operator configured to acquire information items as to plural points within a region including the receive focus and a peripheral region thereof, respectively from the post-delay received signals, a peripheral information combiner configured to combine the information items acquired as to the plural points respectively, and generate a beamforming output by using the information being combined, and an image processor configured to use the beamforming output to generate image data.
According to the present invention, it is possible to acquire ultrasound signals scattered around the focus, reduce of the number of the ultrasound signals failing to be picked up, achieve a highly robust property against the correlative noise caused by the wavefront distortion, and compensate for image-quality deterioration without significantly increasing the processing load.
a) is a perspective view illustrating a schematic configuration of the ultrasound imaging apparatus of the present embodiment, and
a) to
a) to
a) to
a) illustrates a mismatch between the signal spreading area by reflection at a focus (elliptic region 2205) and the point spread function 2204 based on the fixed steering angle, in association with the position of the ultrasound image,
a) is a graph showing the signal strength or the intensity distribution of the B-mode image in the twelfth embodiment, and
a) illustrates the ultrasound images (B-mode images) showing six point scatterers when no wavefront distortion occurs in the ultrasound beam, and
Preferred embodiments of the present invention will be explained.
The ultrasound imaging apparatus of the first embodiment is provided with plural elements configured to receive ultrasound signals from a test subject, a delay unit configured to delay each of the received signals on the plural elements, in association with a position of a predetermined receive focus and generate post-delay received signals, a peripheral information operator, a peripheral information combiner, and an image processor. The peripheral information operator acquires information items as to plural points within a region including the receive focus and the peripheral region of the receive focus, respectively from the post-delay received signals. The peripheral information combiner combines the information items respectively acquired as to the plural points, and uses the information being combined to generate a beamforming output. The image processor uses the beamforming output to generate image data. Specific explanations will be provided in the following.
Firstly, with reference to
As illustrated in
The ultrasound probe 101 is provided with plural elements (ultrasound transducers) 106 arranged in the form of an array. The transmit beamformer 104 generates signals for a transmission beam, and transfers the transmission beam to the ultrasound probe 101 via the transmit-receive separation circuit 107. The ultrasound probe 101 transmits ultrasound waves from the plural elements, directed to the inside of the test subject 100. The ultrasound probe 101 receives echo signals reflected inside the body. The received signals go through the transmit-receive separation circuit 107 again and subjected to beamforming operation, and the like, in the receive beamformer 108. The received signals after the beamforming operation is applied, are transferred to the image processor 109, and various image processes such as various filters, and a scan converter are executed on the received signals, and then, an ultrasound image is generated. The ultrasound image is transferred to the monitor 103, and displayed thereon.
In the receive beamformer 108 of the present embodiment, under the control of the controller 111, the active channel setter (not illustrated) sets an active channel in a part of the finite aperture of the ultrasound probe 101, so as to perform the receive beamforming process. In other words, the elements 106 within a partial range are assumed as the active channel 201, among the plural elements 106 constituting the ultrasound probe 101 which has received echoes associated with one transmit ultrasound beam, and the receive beamformer generates one line of image data (raster: beamforming output y(n)) in the ultrasound propagating direction, by using the received signals in the active channel 201. As illustrated in
The delay circuit 204 corresponds to a block provided as a previous stage of the peripheral information operator 205, configured to give a delay time to each of the received signals (received data) on the plural elements 106 constituting the active channel 201, in association with the element position, so as to perform processing to focus the received signals on a certain point virtually existing in a space. The delay time given to each of the elements 106 is prepared in advance, as a delay-time set, in association with the object positions in the space of the imaging target. The delay circuit 204 selects the delay-time set in response to the focal position being set, and gives different delay times to the respective received signals on the plural elements. With this configuration, it is possible to perform the focusing process covering the entire space of the imaging target.
The peripheral information operator 205 performs the operation to collect information items that are scattered around the receive focus due to the reflection from the receive focus. By way of example, the peripheral information operator 205 performs the adaptive signal processing, thereby obtaining adaptive weights as to two or more points, on the receive focus and in its periphery (peripheral region). With this configuration, it is possible to collect peripheral information with respect to the receive focus. It is to be noted that the periphery (peripheral region) corresponds to an area where signals are scattered due to the reflection of the signals on the receive focus, having a nearly elliptical shape (oval sphere in the three-dimensional case); the size being different depending on the conditions such as receive elements, a transmit frequency, a receive frequency, a ultrasound pulse waveform, the size and pitch of the element array, and a point spread function, roughly having a long axis spreading around ±5° or less in the active channel direction, assuming as the central axis, the axis connecting the receive focus and the center point of the active channel, and a short axis in a length equivalent to or a little shorter than the active channel direction, in the ultrasound propagating direction orthogonal to the active channel direction. On this occasion, the focus is positioned at an intersection of the short axis and the long axis. At the same time, the periphery (peripheral region) corresponds to an area where the wavefront distortion causes that the reflected waves from the medium surrounding the focus exist in the received waves in mixed manner, and the reflected waves from the surrounding medium may become obvious in the received signals in the form of noise signals, being correlated with the information from the focus.
The peripheral information combiner 206 combines (compounds) the peripheral information items collected by the peripheral information operator 205 and obtains a final beamforming output y(n). By way of example, the peripheral information combiner 206 combines the adaptive weights of at least two points obtained by the peripheral information operator 205, and this combined weight being obtained is assigned on the signals of the elements 106, after the delay process is performed by the delay circuit 204, and weighted signals are subjected to beamforming, thereby combining the peripheral information items.
This finally obtained beamforming output y(n) is sequentially transferred to the image processor 109, raster by raster, respectively in association with the active channels 202, 201, and 203, and an ultrasound image is generated on the basis of thus combined beamforming output.
As described above, in the first embodiment, the peripheral information operator 205 collects the information items from the surrounding of the receive focus, and even though the adaptive beamformer with high directivity is used, it is possible to collect the information items reflected on the focus, and scattered in the surrounding thereof. Furthermore, by using thus collected peripheral information, an asymmetry property of the correlative signals is utilized, and this cancels (performs decorrelation of) the correlative noise caused by the wavefront distortion. Therefore, even when the medium of the test subject is inhomogeneous, it is possible to prevent the image quality from being deteriorated, resisting the influence from the correlative noise caused by the wavefront distortion, thereby achieving a highly robust ultrasound imaging apparatus.
In the second embodiment, the peripheral information operator 205 in the ultrasound imaging apparatus of the first embodiment performs the adaptive beamforming, thereby obtaining an adaptive weight as information.
The peripheral information operator 205 obtains the adaptive weights as to plural points, by using steering vectors that are directional vectors connecting a predetermined position on a surface of the array made up of the plural elements 106 as described in the first embodiment, and the plural points.
Furthermore, the peripheral information operator 205 incorporates a matrix operator configured to generate a covariance matrix by using the post-delay received signals, and a weight vector operator configured to obtain adaptive weight vectors as to the plural points, from the steering vectors. The peripheral information combiner 206 incorporates a weight combiner configured to sum the adaptive weights as to plural points obtained by the peripheral information operator 205 and generate a combined weight, and an inner-product operator configured to multiply the post-delay received signals by the combined weight to generate a beamforming output.
It is further possible to arrange a fixed apodization multiplier between the peripheral information operator 205 and the weight combiner, the fixed apodization multiplier being configured to multiply each of the adaptive weights as to the plural points obtained in the peripheral information operator, by a fixed weight being predetermined.
The inner-product operator multiplies each of the post-delay received signals by the combined weight, and thereafter sums the post-delay received signals, so as to generate the beamforming output.
With reference to
In the second embodiment, the peripheral information operator 205 obtains adaptive weights as to two or more points on the receive focus and in the periphery thereof, thereby collecting the peripheral information of the receive focus. Specifically, upon calculating the adaptive weight, plural steering vectors “a” are used so as to obtain the adaptive weight also as to a point at which the direction from the elements in the active channel of the ultrasound probe is displaced from the receive focus (i.e., a point at which the direction of the steering vector is different from the focus direction). The peripheral information combiner 206 combines the adaptive weights of two or more points obtained by the peripheral information operator 205, and uses thus obtained combined weight to perform beamforming process on the received signals. The steering vectors are directional vectors connecting the center position of the active channel and the points as described above. Details of the steering vectors will be explained with reference to the formula (4) described below. Hereinafter, more specific explanations will be given.
As illustrated in
The matrix operator 300 receives inputs from the delay circuit 204 (post-delay received data x(n)), and generates the spatial covariance matrix R(n). After inputted in the adaptive weight operator 301, in each of the plural weight vector operators 3021, 3022, and 3023, the spatial covariance matrix calculates weight vectors w1(n) w2(n), and w3(n), as to the predetermined three points with different steering vector directions, out of the points including the receive focus and in the periphery thereof, and outputs the calculated weight vectors to the peripheral information combiner 206. Details of the steering vectors will be explained in the following.
The peripheral information combiner 206 is provided with the weight combiner 306 and the inner-product operator 307. The weight combiner 306 sums the plural weight vectors w1(n), w2(n), and w3(n) to combine the vectors, and generates a combined weight wsum(n). The inner-product operator 307 uses this combined weight wsum(n) to assign weight on the post-delay received data x(n) that is inputted via the bypass line 207, and then sums the data to obtain the beamforming output y(n). With this configuration, it is possible to obtain information in the form of one beamforming output, the information being acquired from the plural points within an identical raster and being different in the steering vector direction. Specifically the inner-product operator 307 incorporates a multiplier 3071 configured to multiply the post-delay received data x(n) by the combined weight wsum(n), and an adder 3072 configured to calculate the sum of the post-delay data after the multiplication.
If needed, the fixed apodization multiplier 305 may be added to the peripheral information combiner 206. The fixed apodization multiplier 305 multiplies the weight vectors w1(n), w2(n), and w3(n), by a predetermined fixed weight, so as assign the weight on the weight vectors.
The beamforming output y(n) obtained by the peripheral information combiner 206 is transferred to the image processor 109. The image processor 109 generates an ultrasound image on the basis of the beamforming result y(n).
In the following explanation, operations of each part will be explained in more detail, in the case where the received data is processed on each of the elements 106 in one active channel 201 in association with one-time transmission and receiving, so as to generate the beamforming data for one raster.
If the number of the elements constituting the active channel 201 (the number of channels) is K, it is possible to express K post-delay received data items at a certain snapshot time n, as the vector x(n) on the left-hand side of the following formula (1).
[Formula 1]
x(n)=[x1(n),x2(n), . . . ,xK(n)]T (1)
In the adaptive beamformer according to a conventional technique, this post-delay received data x(n) is inputted in the adaptive processor, so as to generate adaptive weight vector w(n)=[w1(n), w2(n), . . . , wK(n)]T being made up of weight values wK(n) for the respective K channels at the snapshot time n. Furthermore, this adaptive weight vectors w(n) and the post-delay received data x(n) are subjected to an inner-product operation in the inner-product operator, and the adaptive beamformer output y(n) at a certain snapshot time n is obtained. In the present invention, plural steering vectors “a” are used to obtain the adaptive weight vectors also for the periphery of the focal position.
As illustrated in
In the formula (2), the ensemble average number N may indicate as a uniform mean as shown on the rightmost side of the formula (2), assuming the total number of samples before and after ξ(n) of the target snapshot as N=2S+1 points. Another method for averaging in the time direction may be taken, the method multiplying each sample in the time direction by an arbitrary weight, such as a trapezoidal weight, and calculating the arithmetic mean. The spatial covariance matrix R(n) outputted from the matrix operator 300 is subsequently inputted into the adaptive weight operator 301.
The adaptive weight operator 301 that has received the spatial covariance matrix R(n) calculates weight vectors w(n) according to the MVDR method. In the present invention, plural steering vectors “a” are used to obtain the peripheral information of the focus. The steering vectors ap is expressed by the formula (4).
In the above formula (4), “p” represents the number of the steering vectors, and it is an integer satisfying 0<p<P+1, when the total number is assumed as P. As shown in the formula (4), the steering vectors ap is a directional vector having K vector elements (from 0 to (K−1)), the number being equal to the number of the active channels. The steering vector is expressed by a function of the receive frequency fp and the angle (θp, φp) formed by the normal vector direction on the surface of the elements 106 and the steering vector (hereinafter, referred to as “steering angle”). Here, “θp” represents an aperture angle from the normal vector, and “φp” represents a swivel angle from the array direction of the elements 106. When the array of the elements 106 in the ultrasound probe 101 is a one-dimensional (linear) array, the steering angle is expressed by the final form of the formula (4). Here, “λp” represents a wavelength of a sound wave associated with the frequency fp, and “d” represents a distance between the center of the elements 106 (pitch between the elements).
The weight vector w (n) obtained by the MVDR method, as to the direction of the steering vector ap as described above is calculated by the formula (5) in this example. Therefore, the weight vectors w(n) are computed as to the steering vectors ap being different respectively in the weight vector operators 3021 to 3023, thereby obtaining the adaptive weight vectors w1(n) to wP(n) corresponding to the number P of the steering vectors ap (i.e., the number of the weight vector operators 3021 to 3023).
In the formula (5), R(n) represents the spatial covariance matrix at a certain snap shot n in the time direction, being generated by the formula (2), and the superscript “−1” represents an inverse matrix.
It is to be noted that the number P of the steering vectors may be any value, as for as it is an integer at least two. Any method for selecting the direction of the steering vector is applicable. In the present embodiment, it is assumed that the focus is positioned on the central axis of the channel 201; i.e., if the number of the active channels is an even number, it is positioned on the normal line (the normal line is orthogonal to the active channel surface) passing through the middle point between the K/2th and the (K+2)/2th elements, out of the K elements 106 constituting the active channel 201. If the number of the active channels is an odd number, it is assumed that the focus is positioned on the normal line (the normal line is orthogonal to the active channel surface) passing through the center of the (K+1)/2th element, out of the K elements constituting the active channel 201. It is desirable that one of the steering vectors faces the focus direction. Another (P−1) steering vectors are configured in such a manner as facing any directions within the region surrounding the focus in the imaging space.
As described above, since the delay process is performed in the delay circuit 204, the received data vector x(n), for instance, at a certain clock time n in the linear scanning becomes data with aligned wavefronts in the direction of θ=0° being the direction of the normal vector. Therefore, in the formula (4), if it is assumed that θp=0°, the steering vectors ap are equal to [1, 1 . . . 1]T, and the adaptive weight vector w(n) in the focus direction is obtained. The adaptive weight vector w(n) in the focus direction corresponds to the adaptive weight vector w(n) obtained by the conventional MVDR method. In other words, if it is assumed that P=1 and θp=0° in the configuration of the present embodiment, it corresponds to the configuration of the adaptive beamformer using an ordinary MVDR method. The beamformer according to the present embodiment is different from the beamformer according to the ordinary MVDR method, in the point that various plural steering vectors ap are set, and adaptive weight vectors w1(n) to wP(n) are obtained for the respective vectors.
The present invention is not limited to the configuration that one of the plural steering vectors ap faces the focus direction (front direction) (θp=0°, ap=[1, 1 . . . 1]T), but all the steering vectors may face the directions different from the focus direction. In particular, as a direction vector constituting the plural steering vectors ap, the direction vector in the focus direction (θp=0°) is not used on purpose, and only the direction vector in the direction (θp<0°, 0°<θp) excluding the focus direction may be used. In this situation, the focus direction is not used on purpose, and this enhances the effect of the correlative noise canceling (decorrelation) more. Therefore, it is one of the preferred embodiments of the present invention.
The P adaptive weight vectors w1(n) to wP(n) outputted from the adaptive weight operator 301 are inputted into the peripheral information combiner 206. The peripheral information combiner 206 sums the P adaptive weight vectors w1(n) to wP(n) in the weight combiner 306, takes an arithmetic average, and calculate the combined weight vector wsum(n) as shown in the formula (6). By way of example, as shown in
It is to be noted that in the stage prior to obtaining the combined weight vector wsum(n), a fixed apodization multiplier 305 may be provided, so as to multiply each of the weight vectors wp(n) by a fixed apodization. By way of example, it is possible to multiply the weight vectors by the fixed apodization by having a distribution where a value of the adaptive weight vector wP(n) in the direction of θ=0° is made larger and values in another directions are made smaller, and this is implemented by the computation as shown in the formula (7).
By way of example, in the configuration of
Next, the combined weight vector wsum(n) is inputted into the inner-product operator 307 within the peripheral information combiner 206. The inner-product operator 307 incorporates the multiplier 3071 and the adder 3072, performs the inner-product operation between the combined weight vector wsum(n) and the received data vectors x(n) transferred from the delay circuit 204 via the bypass line 207, as shown in the formula (8), and then obtains the beamforming output y(n). Specifically, the multiplier 3071 obtains the product of the weight vector and each channel element of the post-delay received data vectors (1 to K). The adder 3072 calculates a total sum of the K products obtained by the multiplier 3071, and obtains a final output of the peripheral information combiner 206 (beamforming output y(n), scalar value).
[Formula 8]
y(n)=wsumH(n)×(n) (8)
This beamforming output y(n) from the inner-product operator 307 is transferred to the image processor 109. The processing above is repeated for the N samples constituting one raster, from the first sample n=1 to the final sample n=N, and then, N beamforming outputs y(n); y(1), y(2) . . . y(N), are transferred sequentially to the image processor 109. Consequently, the receive beamformer 108 completes the computations of the beamforming outputs y(n) for one raster in association with the active channel 201.
The beamforming outputs y(n) for one raster obtained by the formula (8), are acquired respectively for the active channels, shifting on the received array from the active channel 201 to the active channels 202, and 203, and transferred to the image processor 109, raster by raster. In the image processor 109, a scan converter associated with the scanning system of the ultrasound probe 101, arranges all the rasters and generates a two-dimensional image. In addition, backend image processing such as various filter processing, and computations by measurement applications are performed. Finally the monitor 103 displays an ultrasound image and a computational result of the measurement applications.
As described above, in the second embodiment, plural steering vectors ap are set, and the adaptive weight vectors w1(n) to wP(n) are obtained for the respective steering vectors, and combined weight vector wsum(n) obtained by combining those adaptive weight vectors is used to perform beamforming of the received data vectors x(n). Then, one beamforming output y(n) is obtained. With reference to
a) illustrates an example that data is received once from the active array 201, focusing on the receive focus 2201, and the receive beamforming is performed. When the receive beamforming is performed according to the conventional delay-and-sum method, signals in the region indicated by the point spread function 2202 around the receive focus 2201 may also be acquired by one-time receive beamforming. On the other hand, when adaptive beamforming is performed for this case, since the beam directivity is higher than the delay-and-sum method, the range of the signals being obtained by one-time receive beamforming is represented by the point spread function 2204, and it is smaller than the range represented by the point spread function 2202 of the delay-and-sum method. Here, it is conceivable that the signals of the object on the focus 2201 are reflected by the focus 2201 and distributed randomly in the periphery of the focus 2201 as indicated by the elliptic region 2205. This elliptic region 2205 corresponds to the periphery (peripheral region) of the first embodiment.
In the second embodiment, as described above, the adaptive weight vectors obtained by the adaptive beamformer, as to the respective plural points (here three points) having different steering angles are combined. This process is equivalent to allow the regions of the point spread function 2204 to overlap one another, the regions corresponding to the respective peripheries of the plural points as shown in
By using the weight vectors having different steering angles, it is possible to utilize an asymmetry property of noise on the right side and noise on the left side of the central axis, being correlated with each other. In other words, adaptive weights having different steering angles are combined, thereby canceling (performing decorrelation of) the correlative noise signals in the received signals, and therefore, it is possible to selectively eliminate highly correlative noise signals, without failing to pick up the information of the signals on the focus.
According to the second embodiment, the effects as described above bring about improvements against the disadvantages such as; failing to pick up signals due to the sharp beam directivity of the adaptive beamformer, and being vulnerable to influence of the wavefront distortion. Therefore, this provides robustness, being impervious to such influence.
In the second embodiment, as shown in
In the second embodiment, it is alternatively possible to perform spatial average operation using a subarray matrix, as one of other methods of the arithmetic algorithm in the aforementioned matrix operator 300. The subarray matrix R̂SUBl in the spatial average operation is expressed by a product of the subspace vectors ξ̂l(n) (Formula (10)) as shown in the formula (9). The subspace vector ξ̂l(n) is a vector obtained by extracting a partial substance (L elements) from the post-delay received data (that is represented by generalized complex signal vector ξ(n) here, but that may also be represented by a real vector x(n), instead) in association with K active channels. Therefore, the total number of the subspace vectors is equal to (K−L+1) (0<l (l is lower case of L)<(K−L+1)).
[Formula 9]
R
SUBl
̂(n)=ξl̂(n)ξl̂H(n) (9)
[Formula 10]
ξl̂(n)=[ξl(n),ξl+1(n), . . . ,ξl+L−1(n)]T (10)
Samples are displaced one by one, establishing association between the main diagonal component of this subarray matrix and the main diagonal component of the spatial covariance matrix R(n), resulting in a spatial averaging process of the (K−L+1) subarray matrixes, and the subarray spatial covariance matrix R̂(n) as shown in the formula (11) is obtained. When the adaptive weight operator 301 computes the subarray spatial covariance matrix R̂(n), it is used to substitute for R(n) in the above formula (5), and weight vectors wp (n) is computed as shown in the formula (12). It is to be noted that the output from the matrix operator 300 for this case is (L×L) in size, and the number of the elements constituting the weight vectors wp (n) becomes L.
As another alternative example of the spatial averaging method, a forward and backward spatial averaging method may be employed. In this situation, the backward subarray matrix R˜SUBl (n) may be obtained by a product of the backward subspace vector ξ˜l(n) as shown in the formula (13). The backward subspace vector is expressed by the formula (14). According to the computation similar to the forward spatial averaging, the backward subarray spatial matrix R˜(n) may be calculated as shown in the formula (15). Next, as shown in the formula (16), it is subjected to the arithmetic averaging together with the subarray spatial matrix R̂(n) for the case of the aforementioned forward spatial averaging, and finally the forward and backward subarray spatial covariance matrix RFB(n) may be obtained. Similar to the forward spatial averaging, this forward and backward subarray spatial covariance matrix RFB(n) is substituted for R(n) in the formula (5) in the computations in the adaptive weight operator 301, thereby computing the weight vectors wp (n) as shown in the formula (17). Also in this case, the output from the matrix operator 300 is (L×L) in size, and the number of elements constituting the weight vectors wp(n) is L.
Plural (P) weight vectors wp calculated in the formula (12) or in the formula (17) according to the spatial averaging method are transferred to the peripheral information combiner 206, similar to the case where the spatial averaging method is not employed. The weight combiner 306 sums the weight vectors, and calculates wsum (Formula (7), Formula (8)). Furthermore, the inner-product operator 307 performs the inner product calculation using the post-delay received data transferred via the bypass line 207, and eventually, the result is transferred to the image processor 109 as the beamforming data y(n).
In this example, the spatial averaging method is used under the condition that the number of the elements of the weight vectors wp is L. Therefore, in order to perform computations finally in the inner-product operator 307, a block is additionally required, so as to perform the computation that generates vectors g(n) having L components, out of (n) having K components (Formula (18)). For example, as shown in
The beamforming output y(n) using the spatial averaging method is expressed by the formula (19).
As described above, the matrix operator 300 performs the spatial averaging process using the subarray matrix, and produces an effect that restrains correlated noises being included in the received ultrasound signals. Therefore, the spatial averaging process using the subarray matrix is combined with the configuration for synthesizing plural adaptive weight vectors around the focus, thereby obtaining an ultrasound image that is much less influenced by noise.
In the second embodiment as described above, the MVDR is taken as an example for explaining the adaptive beamforming method. However, the algorithm applied in the adaptive weight operator 301 may be any other algorithm as far as it uses the spatial covariance matrix calculated in the matrix operator 300. Any other methods may be applicable, such as the MMSE method, APES method, Eigenspace-MV method (ESMV, EIBMV) using the spatial covariance matrix and its eigenvalues and eigenvectors, ESPRIT method, and MUSIC method.
The third embodiment is different from the second embodiment in the point that the peripheral information combiner 206 multiplies the post-delay received signals, by the adaptive weights respectively for the plural points obtained by the peripheral information operator 205, and plural inner-product operators for generating pre-synthesis beamforming outputs respectively for the plural adaptive weights, and an output combiner are provided. The output combiner sums and combines the pre-synthesis beamforming outputs respectively as to the plural adaptive weights, and generates a beamforming output to be used for generating image data.
It is further possible to arrange between the peripheral information operator 205 and plural inner-product operators, a fixed apodization multiplier configured to multiply the adaptive weights for the respective plural points obtained by the peripheral information operator 205, by fixed weights being predetermined respectively. It is alternatively possible to arrange between the plural inner-product operators and the output combiner, the fixed apodization multiplier configured to multiply the plural post-delay received signals for the respective adaptive weights, by the fixed weights being predetermined respectively. It is further alternatively possible to arrange between the plural multipliers and the plural adders constituting the plural inner-product operators, the fixed apodization multiplier configured to multiply the plural signals after multiplication, by fixed weights being predetermined respectively.
The peripheral information combiner 206 combines the information items that is obtained by the peripheral information operator 205 from the received signals in one active channel, and generates a final beamforming output that is used by the image processor to generate an image. By way of example, the peripheral information combiner 206 sums the pre-synthesis beamforming outputs for the respective plural adaptive weights, obtained by the peripheral information operator 205 from the received signals in one active channel, and generates a final beamforming output that is used by the image processor 109 with respect to each active channel, so as to generate an image.
With reference to
In the second embodiment as shown in
The post-delay received data x(n) is inputted in the multipliers 3071-1, 3071-2, and 3071-3, respectively via the bypass inputs 2071, 2072, and 2073. The output combiner 500 is arranged in the subsequent stage of the plural inner-product operators 307.
With this configuration, in the plural inner-product operators 307, the bypass input 2071 is subjected to the beamforming summation by using the weight vectors w1(n), the bypass input 2072 is subjected to the beamforming summation by using the weight vectors w2(n), and the bypass input 2073 is subjected to the beamforming summation by using the weight vectors w2(n). Then, plural pre-synthesis beamforming outputs y1(n), y2(n), and y3(n) respectively in association with the steering vectors are calculated (Formula (20) in the following). In the output combiner 500 in the subsequent stage, arithmetic averaging of the beamforming outputs y1(n), y2(n), and y3(n) respectively associated with the steering vectors is calculated, and the combined beamforming output ysum(n) is obtained as shown in the formula (21).
As illustrated in
As described above, the second embodiment is different from the third embodiment in the point that the weights are combined and thereafter the inner-product operation is performed; or the beamforming outputs after the inner-product operation are combined. On the other hand, the final output is identical since this is a linear operation, and y(n) in the formula (19) being the final beamforming output of the second embodiment becomes a value being equal to the final combined beamforming output ysum(n) of the third embodiment.
In the fourth embodiment, the delay unit (delay circuit 204) generates post-delay received signals respectively for different plural receive focuses. The peripheral information operator 205 and the peripheral information combiner 206 acquire information items as to the receive focus and plural points in the periphery of the receive focus, as to each of the plural receive focuses, and generates a beamforming output. Byway of example, the delay unit (delay circuit 204) is prepared more than one, and generates the post-delay received signal as to each of the receive focuses being different for the plural delay units, respectively. The peripheral information operator 205 and the peripheral information combiner 206 are arranged for each of the delay units, acquire information from the post-delay received signal generated by the delay unit for each receive focus, and generates the beamforming output.
On this occasion, the positions of the plural receive focuses of the received signals in the active channel at a certain point of time, may partially overlap the positions of the plural receive focuses for the active channel at a different point of time.
With reference to
In the second and the third embodiments, explanations are provided for the case that the focus is positioned on the central axis (the normal line passing through the active channel center, with respect to the active channel surface) for every received signal. In the fourth embodiment, as illustrated in
In the present embodiment, the post-delay data obtained as to each of the plural focuses 601, 602, 603 at a certain depth (at a certain snapshot time n) is subjected to arithmetic operation for collecting and combining the periphery information.
The signal processing operations of the fourth embodiment will be explained. The delay circuit 2042 forms a delay concave surface to obtain the focus 602 on the central axis 1600 of the active channel at a certain snapshot time n, the delay circuits 2041 and 2043 form the delay concave surfaces respectively for obtaining the focuses 601 and 603 on the axes displaced from the central axis 1600, and then the delay circuits generate the post-delay received data.
As shown in
As shown in
Similarly, as shown in
Those beamforming outputs y1sum(n), y2sum(n), and y3sum(n) are transferred to the image processor 109.
The processing above allows acquisition of ultrasound image data of the three focuses 601, 602, and 603, by using just one-time transmit-receive signals to and from the active channel 201, for instance. This processing indicates that when the direction of the focus 602 is assumed as the main beam direction, the image may be generated also in the sub-beam directions (the directions of the focuses 601 and 603). In other words, when the configuration of the fourth embodiment of the present invention is employed, i.e., to acquire signals of the object scattered around the focus, it is found that the configuration of the present invention is applicable to the technique such as the sub-beam processing and parallel beamforming. Therefore, there are obtained effects such as reduction of failing to pick up the signals and reduction of noise.
Furthermore, this process is also applicable to the plural active channels 201, 202, and 203. By way of example, the focus is set for the active channels 202 and 203 on the same point as the focus 602 for the active channel 201, in such a manner as viewing in the slant direction from the active channels 202 and 203, thereby newly acquiring the data regarding the focus 602. Therefore, as for the plural active channels 201, 202, and 203, data items acquired by viewing the identical point 602 from different directions are obtained respectively, and then, the obtained data items are combined and overlapped one another. This process is referred to as synthetic aperture (combined aperture) process. That is, by using the configuration of the fourth embodiment, the technique of the present invention for acquiring signals of the object, scattered around the focus, is applicable to the synthetic aperture process.
According to the present embodiment, an effect of the present invention, that is, reduction of the wavefront distortion, may be produced also in the off-axis beamforming technique focusing on an off-axis point, the parallel beamforming, and the synthetic aperture beamforming.
The internal configurations of the peripheral information operators 2051 to 2053 and the peripheral information combiners 2061 to 2063 of the present embodiment as shown in
In the fifth embodiment, a storing unit (temporary storage 800) is arranged between the plural inner-product operators of the third embodiment and the output combiner, the storage being configured to store the pre-synthesis beamforming outputs respectively for the plural adaptive weights generated by the respective inner-product operators. With reference to
In the fifth embodiment, there is provided the temporary storage 800 as a memory configured to store the beamforming outputs y1(n), y2(n), and y3(n) immediately before the output combiner 500. The temporary storage 800 stores the beamforming outputs y1(n), y2(n), and y3(n) obtained by the plural inner-product operators 307, every updating at the snapshot n, or every updating of the raster. With this configuration, the output combiner 500 is allowed to read and combine the beamforming outputs y1(n), y2(n), and y3(n) stored in the temporary storage 800, and it is possible to combine the beamforming outputs at various time samples within an identical raster, or combine the beamforming outputs across different rasters.
Furthermore, those beamforming outputs y1(n), y2(n), and y3(n) are results of the beamforming using the weight vectors w1(n), w2(n), and w3(n) obtained as to different steering vectors at every snapshot time n. Therefore, the beamforming outputs y1(n), y2(n), and y3(n) indicate beamforming outputs in the respective steering vector directions. This allows the output combiner 500 to select any combination of the beamforming outputs among different rasters and among different steering vector directions, and combine thus selected beamforming outputs. The beamforming outputs to be selected among the different steering vector directions may be predetermined, or set by the operator.
As illustrated in
The peripheral information combiner 206 provided with the temporary storage 800 of the present embodiment may substitute for the peripheral information combiners 2061, 2062, 2063 of the fourth embodiment. This configuration allows synthesis of the beamforming outputs of different time samples as to each of the plural focuses in an identical raster, or synthesis of the beamforming outputs across different rasters, and the like.
In the sixth embodiment, the output combiner 500 acquires and sums the pre-synthesis beamforming outputs different in the time direction, from the storing unit (temporary storage 800) and the plural inner-product operators, so as to generate a beamforming output.
With reference to
In the present embodiment, the temporary storage 800 is provided with a memory 900, and a bypass line 901 for the beamforming outputs (hereinafter, also referred to as “post-beamforming data”) y1(n) to yP(n). The memory 900 has a storage area (address) for storing (P×N) post-beamforming data items, that is, post-beamforming items y1(n) to yP(n) for every N snapshots in the time direction, in association with each of the P steering vectors. Specifically, as shown in
Operations of the temporary storage 800 of the present embodiment are described as the following. The temporary storage 800 receives the post-beamforming data items y1(n) to yP(n) from the plural inner-product operators 307, at a certain snapshot time n. Those post-beamforming data items y1(n) to yP(n) are transferred to both the memory 900 and the bypass line 901. In the memory 900, the post-beamforming data items y1(n) to yP(n) are stored in a predetermined write address (WA) in the memory areas 9001, 9002, and 9003, respectively. Simultaneously, the post-beamforming data items y1(n−1) to yP(n−1) . . . y1(n−i) to yP(n−i) at the sample times n−1 to n−i, stored in the read address (RA) are outputted. Those readout outputs at the sample times n−1 to n−i and the post-beamforming data items y1(n) to yP(n) passing through the bypass line at the current snapshot n are combined as an output, and this output is transferred from the output bus 902 to the output combiner 500.
In the configuration of the present embodiment, even at the sample time n, the beamforming outputs of the plural steering vectors 1 to P are outputted, combining the information items at the time (n−1) . . . (n−i), i.e., being obtained by going back to y1(n−1) to yP(n−1) . . . y1(n−i) to yP(n−i), with the current post-beamforming data items y1(n) to yP(n), and this allows the output combiner in the subsequent stage to combine all of those outputs together.
The information stored in the memory 900 is configured to be overwritten every time the raster is updated, so as to overwrite the information sequentially by the post-beamforming data in a new raster.
As shown in
It is to be noted that the read address (RA) may be determined in advance, or an output adjuster 503 may be provided and configured to designate the read address. In the latter case, the output adjuster 503 may transfer to the temporary storage 800, for instance, the information regarding which data is outputted, out of the current post-beamforming data items y1(n) to yP(n) and the previous data held in the memory, and a signal for controlling the output timing and the output data address.
In the seventh embodiment, the plural inner-product operators generate the pre-synthesis beamforming outputs respectively for the plural adaptive weights for each active channel, and the output combiner 500 acquires from the storing unit (temporary storage 800) and plural inner-product operators, the pre-synthesis beamforming outputs generated for different active channels and sums the outputs, so as to generate the beamforming output.
With reference to
The current post-beamforming data items y1(r, n), y2(r, n), and y3(r, n) at the snapshot time n are written into the write address WA as shown in the figure, and simultaneously those data items are also outputted to the output combiner 500 via the bypass line 901.
In
According to the present embodiment, the post-beamforming data items y1(r,1), y2(r,1) . . . yP(r,1) to y1(r,N), y2(r,N) . . . yP(r,N) of the current scanning raster (the rth raster) at the snapshot time 1 to N can be held in the memory 900-1, and further, the post-beamforming data items of the (r−1)th raster and the (r−2)th raster are also stored respectively in the memory 900-2 and in the memory 900-3. Therefore, the post-beamforming data items, relating to various rasters, snapshot times, and steering vectors, are respectively stored in the memories 900-1 to 900-3. Accordingly, at least one desired post-beamforming data item is selected from those post-beamforming data items, and outputted to the output combiner 500.
By way of example, as shown in
According to this process, post-beamforming data items acquired from plural rasters are combined with regard to an identical focus 1500, and compared to the case of using the post-beamforming data obtained from one steering vector in a single raster, it is possible to achieve higher resolution of the point image, and to reduce noise caused by combining plural rasters.
It is to be noted that the post-beamforming data items stores in the memories 900-1, 900-2, and 900-3 are overwritten every time the raster is updated, and the latest post-beamforming data items for the three rasters are constantly stored in the memories 900-1, 900-2, and 900-3, respectively.
Further in the seventh embodiment, the output adjuster 503 may be provided. In this case, the output adjuster 503 transfers to the temporary storage 800, a signal to control the data address to be outputted and its timing to the output combiner 500.
With reference to
Therefore, as shown
In other words, as seen from
As shown in
According to the present embodiment, though not clearly illustrated in
Also in the present embodiment, the output adjuster 503 may be provided, so as to control the data address to be outputted and its timing to the output combiner 500.
In the seventh and the eighth embodiments, information in the adjacent rasters is used, and thus at the clock time of the rth raster being currently scanning, the (r−1)th combined data is transferred to the image processor 109. Therefore, there is a certain delay for the timing of image processing.
As illustrated in
In the ninth embodiment, an explanation will be made as to the plural elements 106 that are in the form of a two-dimensional array, i.e., being arranged two-dimensionally. The peripheral information operator 205 generates pre-synthesis beamforming outputs as to plural points being set two-dimensionally. The storing unit (temporary storage) 800 stores those pre-synthesis beamforming outputs.
With reference to
The adaptive weight operator 301 varies the steering angle being calculated, in the two-dimensional direction where the elements 106 of the probe 101 are arranged, and obtains adaptive weight vectors as to P×Q steering angles. Within the Q temporary storages 800, the post-beamforming data items obtained as to the adaptive weight vectors for the respective steering vectors are stored at every snapshot time. With this configuration, it is possible to collect the peripheral information (post-beamforming data) three-dimensionally.
In the present embodiment, the steering vectors ap expressed by the aforementioned formula (4) are extended in the two-dimensional direction, and the present invention allows the steering vector direction to be selected, in the form of any combination of (θp, φq), particularly using 1≦p≦P and 1≦q≦Q in the two dimensional angles. The formula (24) expresses the steering vector a(p, q) for this case.
[Formula 24]
a
(p,q)=[exp{ψ1(θp,φq,f(p,q))},exp{ψ2(θp,φq,f(p,q))}, . . . ,exp{ψK(θp,φq,f(p,q))}] (24)
In the tenth embodiment, an explanation will be provided as to the adjuster 501 that adjusts at least either one of the number of the steering vectors and the directions thereof.
a) to
d) to
In addition, as shown in
As for a combination of the steering angle and the number of the steering vectors, the adjuster 501 stores data in advance in such a manner as storing plural types of combinations of the steering angle and the number of the vectors, and depending on the focal position (the position of the imaging target point), the imaging portion of the test subject 100, and the imaging sequence, the adjuster 501 selects a suitable combination of the steering angle and the number of the vectors, so as to set the selected combination in the angle adjuster 502 and in the output adjuster 503. It is alternatively possible that the operator selects via the console 110, a combination of the steering angle and the number of vectors, or the operator inputs in the adjuster 501 via the console 110, any steering angle and any number of the vectors.
In the eleventh embodiment, there is further provided a memory configured to store in advance, a distribution of the spread angles of the steering vectors directed to plural points, in association with the position of the receive focus within the image that is generated by the image processor. The adjuster uses the steering vectors having the spread angles read out from the memory, in association with the position of the receive focus, so as to obtain the adaptive weights. The spread angle of the steering vector stored in the memory is configured to be small for the receive focus being close to the send focus upon transmitting ultrasound wave signals to the test subject, whereas it is configured to be large for the receive focus being distant from the send focus. In addition, the spread angle of the steering vector is configured to be smaller at the edge part of the image, compared to the center part of the image.
With reference to
In general, the focus upon sending is set on at least one fixed position when one ultrasound image is taken, and the receive beamformer 108 varies the receive focus, thereby obtaining an image for each imaging target point within the imaging area. Therefore, the focus 2210 upon sending is determined at a position within the image. On the other hand, as explained with reference to
In the periphery of the first raster and the last raster in the scanning direction among the rasters upon receiving, that is, on the image edge parts, a part of the point spread function 2204 where signals are allowed to be acquired, being determined by the spread angles of the steering vectors, runs off the image as shown in
The aforementioned mismatch and extending off the image of the point spread function 2204 may result in that intensity unevenness shows up in a B-mode image. For example, as illustrated in
In the eleventh embodiment, in order to reduce the intensity unevenness in the depth direction, the steering spread angle is varied in accordance with the signal spreading region of the object, so that the point spread function 2204 upon beamforming matches the elliptic region 2205. In the image edge part, a steering angle is used which restrains spreading of the point spread function 2204 determined by the spread angles of the steering vectors, so that the point spread function may not expand unnecessarily out of the imaging target. By way of example, beamforming according to the following methods may be performed at the image edge part; a method of reducing the steering angle gradually, as the object becomes deeper, a method of preparing a fixed small steering angle, a method of using an output of normal adaptive beamforming (an output using the steering vectors of θp=0°) for the edge part, without performing synthesis of the steering angles, and the like.
In order to implement the configuration above, the adjuster 501 is provided in the memory 2215 as shown in
The adjuster 501 receives from the controller 111, information regarding the type of the probe 101 and the imaging sequence, and reads out from the memory 2215, data defining the distribution of the spread angles of the steering vectors in association with the information. According to the data being read out, the adjuster 501 sets the spread angles (steering angles) of the steering vectors in the angle adjuster 502 and in the output adjuster 503 in the receive beamformer 108, for each scanning direction and depth (time) direction of the raster.
Accordingly, it is possible to allow the point spread function 2204 upon beamforming to coincide with the elliptic region 2205 indicating the spread of the signals of the object. Further in the image edge part, the point spread function 2204 is controlled not to expand unnecessarily out of the imaging target. Therefore, as shown in
In the eleventh embodiment, the configuration of the receive beamformer 108 is not limited to the configuration of
In this example here, the case where the send focus is fixed to one or plural points is explained. If a dynamic focus for sending is used so as to change the send focus in the depth direction, there may occur unevenness in the elliptic region 2205 indicating the spread of signals of the object. Therefore, the tenth embodiment is applicable.
a) to
With reference to
Specifically, the beamforming of any of the first to the ninth embodiments is performed initially, with the use of the fixed steering angle as shown in
The intensity distribution or the signal strength distribution 230 being obtained is fitted to a predetermined function for converting the distribution of the intensity or the signal strength into the steering angle distribution, so as to obtain the steering angle distribution 231 as shown in
The function for converting the distribution of the intensity or signal strength into the steering angle distribution is obtained in advance by experiments, or computations (including a computer simulation), and the memory 2215, or the like, stores the function in the form of database like a table, or in the form of mathematical formula.
As described in the eleventh embodiment, by using the distribution of the spread angles (steering angles) of the steering vectors as stored in the memory 2215, the receive beamforming is performed as shown in
The obtained steering angle distribution 231 is compared with the steering angle distribution 232 used for generating the B-mode image of
With the configuration above, it is possible to perform modification with the feedback of the steering angle distribution 232, or the distribution of intensity or signal strength 230 of the image.
As explained in the present embodiment, the steering angle distribution may be obtained actively, or modified with the feedback operation. This process is applicable not only to the intensity unevenness when the adaptive beamformer is utilized, but also the intensity unevenness due to energy inhomogeneity of acoustic signals that occurs also in the beamforming according to a conventional DAS (Delay And Sum: delay addition process). Therefore, this may contribute to enhancement of an ultrasound B-mode image inherently.
A mode-switching part in the console may perform switching, in association with the various types of the probe, imaging conditions, and imaging sequences. For example, the switching part may switch to the mode that applies a set of various steering angles, being different depending on the probe types, such as a linear probe, a convex probe, a sector probe, a two-dimensional array probe, and a mechanical 3D imaging probe, so as to generate an image. Alternatively, for example, the switching part may switch to the mode that applies a set of various steering angles, being different depending on the imaging conditions, such as the send focus and receive focus, send/receive frequency, frame rate, parallel beamforming count, tissue harmonics, and contrast, and the imaging sequence, so as to generate an image.
It is further possible to configure such that the operator sets the arithmetic parameters respectively of the peripheral information operator 205 and the peripheral information combiner 206, through the use of the operating portion 1403 that is manipulated by the operator. On this occasion, the adjuster 501 as shown in
When the images 2301 and 2304 including the wavefront distortion are compared with the images 2300 and 2303 not including the wavefront distortion, it is found that even though the same beamforming method is applied, there exist belt-like high-intensity areas around the point scatterers, in the case where the wavefront distortion occurs, and it is difficult to distinguish the point scatterers.
When the image 2301 obtained by the conventional delay-and-sum method is compared with the image 2302 obtained by combining the steering vectors in the conventional delay-and-sum method, it is found that even though the steering vectors are combined in the delay-and-sum method, the belt-like high intensity areas around the scatterers are hardly improved, failing to remove the influence of the wavefront distortion.
On the other hand, when the image 2304 obtained according to the adaptive beamforming is compared with the image 2305 obtained by combining the steering vectors in the adaptive beamforming of the second or the third embodiment, obviously the belt-like high intensity area is decreased in the image 2305, and the six point scatterers are separated and recognizable. Accordingly, it is found that the influence of the wavefront distortion is removed by combining the steering vectors in the adaptive beamforming according to the second or the third embodiment.
As indicated by the dotted line 2401, the intensity peak position in the intensity distribution of the image 2301 according to the delay-and-sum method is displaced from the real position of the object. As indicated by the dashed-dotted line 2402, the intensity peak position in the intensity distribution of the image 2304 according to the conventional adaptive beamformer becomes closer to the real position of the object, but the magnitude of the intensity is lowered. As indicated by the solid line 2403, in the intensity distribution of the image 2305 according to the second or the third embodiment of the present invention, the signal strength is maintained to be equivalent to the intensity distribution by the delay-and-sum method (the dotted line 2401), the peak intensity is located on the real position of the object (the position zero on the horizontal axis), and therefore this ascertains the effect produced by compensating for the wavefront distortion.
As described above, according to the present invention, the adaptive beamformer is provided with the ability of avoiding the failure in picking up information due to sharp beam directivity. In addition, there is provided another ability of canceling and reducing the unnecessary correlative noise from the medium around the focus. Accordingly, the present invention provides the ultrasound imaging apparatus having a robust property against the wavefront distortion, caused by sound-velocity inhomogeneity within a living body, a distribution of scatterers, and body motion influence. The present invention further allows a point image made up of plural perspective angles to become high in resolution, without performing the focus calculation (delay calculation) one by one for each target point in the test subject 100. Therefore, according to the present invention, it is possible to provide the ultrasound imaging apparatus with the adaptive beamformer that needs relatively small processing loads for solving the problems above.
Here is a summary of the configurations of the present invention.
(1) The first embodiment relates to a configuration of the apparatus provided with the peripheral information operator and the peripheral information combiner. In other words, the ultrasound imaging apparatus is provided with plural elements configured to receive ultrasound signals from a test subject, a delay unit configured to delay each of the signals received by the plural elements in association with a predetermined position of the receive focus, and generate post-delay received signals, the peripheral information operator configured to acquire from the post-delay received signals, information items as to plural points on the receive focus and in the region surrounding the receive focus, the peripheral information combiner configured to combine the information items respectively acquired as to the plural points and generate a beamforming output by using the information items being combined, and an image processor configured to generate image data by using the beamforming output.
(2) The fourth embodiment relates to a configuration for performing synthetic aperture in the apparatus as described in the aforementioned item (1). In other words, the fourth embodiment relates to the ultrasound imaging apparatus where the delay unit generates the post-delay received signals respectively as to plural different receive focuses, and the peripheral information operator and the peripheral information combiner acquire the information items as to the plural points on the receive focus and in the region surrounding the receive focus, with respect to each of the plural receive focuses, and generate the beamforming output.
(3) The fourth embodiment also relates to a configuration that there are plural delay units in the apparatus as described in the aforementioned item (2). In other words, there are more than one delay unit, and each of the delay units generates the post-delay received signals for the receive focus that is different for each of the plural delay units. The peripheral information operator and the peripheral information combiner are provided for each of the delay units, acquire the information items from the post-delay received signals generated by the delay unit for the each of the receive focuses, and generate the beamforming output.
(4) The fourth embodiment also relates to a configuration that the synthetic aperture is performed in the apparatus as described in the aforementioned item (2) or in the item (3). In other words, there is further provided an active channel setter configured to set active channels sequentially to the plural elements, at the positions being different in time series, and transfer to the delay unit, the received signals of the elements included in the active channel. On this occasion, the positions of the plural receive focuses as to the received signals in the active channel at a certain point of time partially overlap the positions of the plural receive focuses as to the active channel at a different point of time.
(5) The second embodiment relates to a configuration that the adaptive beamforming is applied, in any of the apparatus in the aforementioned items (1) to (4). In other words, the peripheral information operator performs the adaptive beamforming, thereby obtaining the adaptive weight as the information item.
(6) The second embodiment relates to a configuration that the peripheral information operator in the apparatus as described in the aforementioned item (5) uses the steering vectors being directional vectors connecting a predetermined element among the plural elements and the plural points, so as to obtain the adaptive weights as to the plural points.
(7) The second embodiment relates to a configuration that in the apparatus as described in the aforementioned item (6), the covariance matrix is used to obtain the adaptive weight vectors. In other words, the peripheral information operator includes the matrix operator configured to use the post-delay received signals to generate the covariance matrix, and the weight vector operator configured to obtain the adaptive weight vectors as to the plural points, from the covariance matrix and the steering vectors.
(8) The second embodiment relates to a configuration that in the apparatus as described in any of the aforementioned items from (5) to (7), the adaptive weights are combined. In other words, the peripheral information combiner includes the weight combiner configured to sum the adaptive weights as to the plural points obtained by the peripheral information operator and generate the combined weight, and the inner-product operator configured to perform inner-product operation between the combined weight and the post-delay received signals, and generate the beamforming output.
(9) The second embodiment relates to a configuration that in the apparatus as described in the aforementioned item (8), the adaptive weight is multiplied by the fixed weight. That is, between the peripheral information operator and the weight combiner, there is arranged a fixed apodization multiplier configured to multiply the adaptive weights as to the plural points obtained by the peripheral information operator, respectively by the predetermined fixed weights.
(10) The second embodiment also relates to a configuration that in the apparatus as described in the aforementioned item (8), assigning weights on the post-delay received signals and summing the signals are performed. In other words, the inner-product operator multiplies each of the post-delay received signals by the combined weight, and thereafter, sums the post-delay received signals to generate the beamforming output.
(11) The third embodiment relates to a configuration that in the apparatus as described in any of the aforementioned items from (5) to (7), the received signals as to the plural points are subjected to beamforming using the respective adaptive weights, and then those signals are added. In other words, the peripheral information combiner includes the plural inner-product operators configured to perform inner-product operation between the adaptive weights as to the plural points obtained by the peripheral information operator and the post-delay received signals respectively, and generate the pre-synthesis beamforming outputs for the respective plural adaptive weights, and the output combiner configured to add and combine the pre-synthesis beamforming outputs as to the plural adaptive weights, and generate the beamforming output that is used for generating the image data.
(12) The third embodiment also relates to a configuration that in the apparatus as described in the aforementioned item (11), the plural adaptive weights are multiplied respectively by the fixed weights. That is, between the peripheral information operator and the plural inner-product operators, there is arranged the fixed apodization multiplier configured to multiply the adaptive weights as to the plural points obtained by the peripheral information operator, by predetermined fixed weights, respectively.
(13) The third embodiment also relates to a configuration that in the apparatus as described in the aforementioned item (11), the post-delay received signals after the beamforming are multiplied by the fixed weights, respectively. That is, between the plural inner-product operators and the output combiner, there is arranged the fixed apodization multiplier configured to multiply the post-delay received signals after the beamforming as to the plural points, by predetermined fixed weights, respectively.
(13-A) The third embodiment also relates to a configuration that in the apparatus as described in the aforementioned item (11), assigning weights on the post-delay received signals and adding the signals are performed. In other words, the plural inner-product operators are made up of a set of multipliers and plural adders, and the plural multipliers multiply the post-delay received signals by the adaptive weights as to the plural points respectively, and thereafter, in the plural adders, the post-delay received signals are added, thereby generating the pre-synthesis beamforming outputs.
(13-B) The third embodiment also relates to a configuration that in the apparatus as described in the aforementioned item (13-A), multiplication by the fixed weights is performed. That is, between the plural multipliers and the plural adders, there is arranged the fixed apodization multiplier configured to multiply the post-multiplication received signals through the plural multipliers, by the fixed weights being predetermined as to the plural points obtained by the peripheral information operator.
(14) The third embodiment also relates to a configuration that in the apparatus as described in any of the aforementioned items (1) to (13), the beamforming outputs are combined within an identical raster. In other words, there is further provided an active channel setter configured to set active channels sequentially to the plural elements, at the positions being different in time series, and transfer to the delay unit, the received signals of the elements included in the active channels. The peripheral information combiner combines the information items that the peripheral information operator obtains from the received signals in one of the active channels, and generates a final beamforming output that is used by the image processor.
(15) The third embodiment also relates to a configuration that in the apparatus as described in the aforementioned items (11), the beamforming outputs are combined within an identical raster. In other words, there is further provided the active channel setter configured to set active channels sequentially to the plural elements, at the positions being different in time series, and transfer to the delay unit, the received signals of the elements included in the active channels. The peripheral information combiner sums the pre-synthesis beamforming outputs respectively as to the plural adaptive weights that the peripheral information operator obtains from the received signals in one of the active channels, and generates a final beamforming output that is used by the image processor, with respect to each of the active channels.
(16) The fifth embodiment relates to a configuration that in the apparatus as described in the aforementioned item (11), synthesizing is performed by the use of a temporary storage. In other words, between the beamformer and the adder, there is provided the storing unit configured to store the pre-synthesis beamforming outputs generated by the beamformer, as to each of the plural adaptive weights.
(17) The sixth embodiment relates to a configuration that in the apparatus as described in the aforementioned item (16), synthesizing is performed using plural samples in the time direction. That is, the adder acquires from the storing unit and the beamformer, the pre-synthesis beamforming outputs different in the time direction, and sums those pre-synthesis beamforming outputs, so as to generate the beamforming output.
(18) The seventh embodiment and the eighth embodiment relate to a configuration that in the apparatus as described in the aforementioned item (16), synthesizing is performed between plural rasters. In other words, there is further provided the active channel setter configured to set active channels sequentially to the plural elements, at the positions being different in time series, and transfer to the delay unit, the received signals of the elements included in the active channels. The beamformer generates the pre-synthesis beamforming outputs respectively as to the plural adaptive weights, for each of the active channels, the adder acquires from the storing unit and the beamformer, the pre-synthesis beamforming outputs generated for the different active channels and sums the beamforming outputs, so as to generate the beamforming output.
(19) The ninth embodiment relates to a configuration that in the apparatus as described in the aforementioned item (16), the plural elements correspond to a two-dimensional array arranged two-dimensionally, and the peripheral information operator generates the pre-synthesis beamforming outputs as to the plural points set in the two-dimensional direction, and the storing unit stores the pre-synthesis beamforming outputs.
(20) The tenth embodiment relates to a configuration that in the apparatus as described in the aforementioned item (6), there is provided the adjuster configured to adjust at least either one of the number of the steering vectors and the directions of the steering vectors.
(21) The eleventh embodiment relates to a configuration that in the apparatus as described in the aforementioned item (20), setting of the steering angle distribution within an image is performed. In other words, there is provided the memory configured to store in advance a distribution of spread angles of the steering vectors directed to the plural points, in association with the positions of the receive focus within the image that is generated by the image processor, and the adjuster obtains the adaptive weights, by using the steering vectors with the spread angles being read from the memory, in association with the position of the receive focus.
(22) The eleventh embodiment relates to a relation between the steering angle distribution and the send focus in the apparatus as described in the aforementioned item (21). In other words, the spread angle of the steering vector stored in the memory is set to be small on the receive focus being close to the send focus upon sending the ultrasound signal to the test subject, and the spread angle is set to be large at the receive angle being distant from the send focus.
(23) The eleventh embodiment relates to a relation between the steering angle distribution and the image edge in the apparatus as described in the aforementioned item (21) or item (22). In other words, the spread angle of the steering vector stored in the memory is set to be smaller at the edge portion, compared to the center of the image.
(24) The twelfth embodiment relates to a configuration that in the apparatus as described in any of the aforementioned items (21), (22), and (23), the steering angle distribution within the image is computed from the B-mode image. In other words, the adjuster employs a predetermined function to convert the distribution of the B-mode image intensity or signal strength of the test subject imaged in advance, into the spread angle distribution of the steering vectors being obtained, and uses thus obtained spread angle distribution of the steering vectors.
(25) The twelfth embodiment relates to a configuration that in the apparatus as described in any of the aforementioned items (21), (22), and (23), the steering angle distribution within the image is corrected. In other words, the adjuster uses the predetermined function to convert the distribution of the B-mode image intensity or signal strength, obtained by setting the steering vectors with the use of the spread angle distribution of the steering vectors that is stored in the memory, into the distribution of the spread angle of the steering vectors. With the use of a difference between the spread angle distribution of the steering vectors within the memory, and the spread angle distribution of the steering vectors obtained using the function, the spread angle distribution, the image intensity, or the signal strength of the steering vectors within the memory is corrected.
(26) The apparatus as described above is further provided with the monitor for displaying the image and the operating portion, and at least either of the monitor and the operating portion is provided with an accepting portion for accepting from the operator, a setting of the steering vectors in the adjuster.
Number | Date | Country | Kind |
---|---|---|---|
2012-030991 | Feb 2012 | JP | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2013/051305 | 1/23/2013 | WO | 00 |