Ultrasound image pickup apparatus
11160533 · 2021-11-02
Assignee
Inventors
Cpc classification
G01S7/52046
PHYSICS
International classification
A61B8/00
HUMAN NECESSITIES
Abstract
A transmitting beamformer performs convergence transmission that forms a transmission focus of an ultrasonic beam in a subject. A receiving beamformer comprises a virtual sound source method-based delay amount calculation part that obtains delay amount of a received signal with regarding the transmission focus as a virtual sound source, and a correction operation part that corrects the delay amount obtained by the virtual sound source method-based delay amount calculation part depending on position of imaging point. Delay amounts can be thereby obtained with good accuracy for imaging points in a wide area.
Claims
1. An ultrasonic imaging apparatus, comprising: an ultrasonic element array in which a plurality of ultrasonic elements are arranged along a predetermined direction; a transmitting beamformer that forms an ultrasonic beam to be transmitted into a subject by the ultrasonic element array; a computer coupled to the ultrasonic element array and the transmitting beamformer; and a controller coupled to the computer that generates image data using results outputted by a receiving beamformer, wherein the transmitting beamformer performs convergence transmission that forms a transmission focus of the ultrasonic beam in the subject, and wherein the computer is programmed to: perform phasing of a plurality of received signals obtained by receiving the ultrasonic waves reflected in the subject with the ultrasonic element array by delaying received signals, obtain delay amounts of the received signals based on the transmission focus as a virtual sound source, and correct the obtained delay amounts based on positions of an imaging point, wherein, the imaging point is outside a transmitted sonic wave region bound by two virtual lines projecting from different ultrasonic elements of the ultrasound element array and through the transmission focus, wherein the computer is further programmed to: correct a first delay amount (D.sub.1) obtained for the imaging point using a second delay amount (D.sub.2) obtained for a point located on one of the two virtual lines of the transmitted sonic wave region, wherein the point located on one of the two virtual lines of the transmitted sonic wave region is a projection of the imaging point in a direction perpendicular to a sound axis of the ultrasonic beam.
2. The ultrasonic imaging apparatus according to claim 1, wherein computer is further programed to: obtain the second delay amount (D.sub.2) by an operation using the first delay amount (D.sub.1), and an angle θ formed by a line connecting the imaging point and the transmission focus, and a sound axis of the ultrasonic beam.
3. The ultrasonic imaging apparatus according to claim 1, wherein the computer is further programed to: weight the first delay amount (D.sub.1) and the second delay amount (D.sub.2) with weight values and to add the first delay amount (D.sub.1) and the second delay amount (D.sub.2).
4. The ultrasonic imaging apparatus according to claim 3, wherein the weight values used for the weighting vary depending on an angle θ formed by a virtual line connecting the imaging point and the transmission focus, and a sound axis of the ultrasonic beam.
5. The ultrasonic imaging apparatus according to claim 3, wherein the computer is programmed to receive input of the weight values.
6. The ultrasonic imaging apparatus according to claim 1, wherein the computer is programmed to: store phased outputs, and determine a phased output obtained from a received signal of the ultrasonic beam of a certain transmission, and synthesizes it with a phased output obtained from a received signal of the ultrasonic beam of another transmission.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
MODES FOR CARRYING OUT THE INVENTION
(18) Hereafter, embodiments of the present invention will be explained with reference to the drawings. However, the present invention is not limited to the following embodiments.
First Embodiment
(19) The ultrasonic imaging apparatus of the first embodiment of the present invention comprises, as shown in
(20) The transmitting beamformer 602 performs convergence transmission that forms the transmission focus 203 of the ultrasonic beam 104 in the subject. The receiving beamformer 603 comprises a virtual sound source method-based delay amount calculation part 609 that obtains delay amount of a received signal with regarding the transmission focus 203 as a virtual sound source, and a correction operation part 610 that corrects the delay amount obtained by the virtual sound source method-based delay amount calculation part 609 depending on position of imaging point.
(21) For example, when an imaging point locates outside the region between transmitted sonic wave ends 206, which are two lines connecting end ultrasonic elements 600a and 600b at the both ends among ultrasonic elements 600 that transmit the ultrasonic beam 104 into the subject, and the transmission focus 203, respectively (for example, imaging point 802 shown in
(22) Specifically, as shown in
(23) The point 803 is a point locating on the line of the transmitted sonic wave end 206 or between two of the transmitted sonic wave ends 206. The point 803 is preferably a point obtained by projecting the imaging point 802 on the nearest transmitted sonic wave end 206, because such a point reduces calculation complexity. The imaging point 802 is projected on the transmitted sonic wave end 206 by, for example, moving the imaging point 802 along the direction perpendicular to the sound axis 702 of the ultrasonic beam 104.
(24) The correction operation part 610 obtains a corrected delay amount D by, for example, weighting the aforementioned delay amount (D.sub.1) and the delay amount (D.sub.2), and adding them. For obtaining weight values for the weighting, there is used, for example, a function of the angle θ formed by the line connecting the imaging point 802 and the transmission focus 203, and the sound axis 702 of the ultrasonic beam 104 as a variable (
(25) Hereafter, still more specific explanation will be made. As shown in
(26) The virtual sound source method-based delay amount calculation part 609 obtains delay amount (D.sub.1) of a received signal by the virtual sound source method according to the distance of a desired imaging point from the transmission focus 203 with supposing that the virtual sound source locates at the transmission focus 203. Since the calculation method of the delay amount according to the virtual sound source method is a widely known method, detailed explanation thereof is omitted in this specification. But as explained with reference to
(27) However, if the actual wave face of transmitted sonic waves is obtained by simulation, the wave face 700 of transmitted sonic waves is substantially perpendicular to the transmission sound axis 702 around the transmission focus 203, as shown in
(28) Therefore, in this embodiment, the correction operation part 610 obtains a delay amount (D) by correcting the delay amount (D.sub.1) obtained by the virtual sound source method-based delay amount calculation part 609 for a desired imaging point. The delay amount (D) for an imaging point in the neighborhood of the transmission focus 203 may be thereby obtained with good accuracy. Moreover, even if the imaging point locates outside the transmitted sonic wave end 206, the delay amount (D) may be obtained with good accuracy.
(29) The correction operation performed by the correction operation part 610 will be explained with reference to
(30) The delay amount (D.sub.2) obtained for the point 803 by the virtual sound source method may also be obtained an operation performed by the virtual sound source method-based delay amount calculation part 609, or it may also be obtained by calculation from the value of the delay amount (D.sub.1) obtained for the imaging point 802. By calculating the delay amount (D.sub.2) from the delay amount (D.sub.1), calculation complexity imposed on the receiving beamformer 603 may be reduced. Specifically, when the point 803 is a point on the transmitted sonic wave end 206, the delay amount (D.sub.2) is obtained from the value of the delay amount (D.sub.1) in accordance with the equation (1).
(31) (Several 1)
D.sub.2=(D.sub.1|cos θ|)/cos θ.sub.1 (1)
θ.sub.1 is the angle formed by the transmitted sonic wave end 206 and the sound axis 702, and θ is the angle formed by a line connecting the imaging point and the transmission focus 203, and the sound axis 702.
(32) The correction operation part 610 obtains a corrected delay amount D of the imaging point 802 by weighting the delay amount D.sub.1 obtained for the imaging point 802 by the virtual sound source method, and the delay amount (D.sub.2) obtained for the point 803 by the virtual sound source method, and adding them, as shown by the equation (2).
(33)
(34) In the equation (2), a is a weight value, and is obtained in accordance with the following equations (3).
(35)
(36) As shown by the equation (3), the weight value a is a function of the variable sin θ, and such a weight value as drawn in the graph of
(37) When the absolute value of sin θ is 1, the imaging point 802c is at a position on the same horizontal level as that of the transmission focus 203. In such a case, a.sub.0 is set as a (a.sub.0 is a constant larger than 1) as shown in
(38) When the absolute value of sin θ is not smaller than the absolute value of sin θ.sub.1 and smaller than 1, a is set to be a value larger than 1 and smaller than a.sub.0 depending on the value of sin θ as shown in
(39) Thus, for the imaging points 802a and 802e locating within the region between the transmitted sonic wave ends 206 among the imaging points locating on the scanning line 901 as shown in
(40) Operation of the whole receiving beamformer 603 mentioned above is explained with reference to
(41) First, in the step 1901 shown in
(42) Then, in the step 1902, the correction operation part 610 obtains the delay amount (D.sub.2) for a point obtained by projecting each imaging point on the transmitted sonic wave end 206 from the delay amount (D.sub.1) for each imaging point obtained in the step 1901 in accordance with the equation (1) mentioned above. As θ.sub.1 contained in the equation (1), a value obtained by an operation from depth of the transmission focus 203 received from the transmitting beamformer 602 and the positions of the both ends of the driven ultrasonic elements 600 is used. Alternatively, a value obtained beforehand for each depth of the transmission focus 203 may also be used as θ.sub.1. As θ, a value obtained by an operation from the positions of the imaging point on the scanning line 901 and the transmission focus 203 is used. Alternatively, a value obtained beforehand for every combination of each imaging point and the transmission focus 203 may be used as θ. Furthermore, the correction operation part 610 substitutes such values as mentioned above for θ.sub.1 and θ in the equation (3) to obtain weight value a for each imaging point. The obtained values of a, D.sub.2, and the delay amount D.sub.1 obtained in the step 1901 are substituted for those of the equation (2) to calculate the corrected delay amount (D).
(43) The process advances to the step 1903, and the correction operation part 610 receives signals received by each ultrasonic element 600 via the transmission/reception separation circuit (T/R) 604, and performs phasing of them by delaying them by the corrected delay amount (D), and adding them. By performing this operation for each imaging point on the scanning line 901, phasing of the received signals is carried out for every imaging point on the scanning line 901 to generate an image of one raster (phased output), and it is delivered to the image processing part 605.
(44) The image processing part 605 performs processings for putting the phased outputs (rasters) of a plurality of scanning lines 901 in order etc. to generate an ultrasonogram, and displays it on the image display part 607.
(45) As described above, according to the present invention, there may be obtained a delay amount (D) by correcting the delay amount (D.sub.1) obtained by the virtual sound source method. Since the delay amount (D) may be obtained with good accuracy irrespective of whether the imaging point locates within the main lobe or not, sufficiently accurate phased output (raster) is obtained even for the scanning line 901 at a position remote from the transmission focus 203. Therefore, the area in which the scanning line 901 can be set is wide, and a plurality of phased outputs (rasters) for a wide area are generated by one time of transmission. A highly precise image is thereby generated with a small number of times of transmission.
(46) Although
(47) Although the weight value a is set on the basis of the angle θ in the aforementioned embodiment, there may also be employed a configuration that the weight value a is set on the basis of the distance from the transmission focus 203 to the imaging point 802, or the like, or on the basis of both the angle θ and the distance.
(48) Further, although the method of moving the imaging point 802 along the direction perpendicular to the sound axis is explained as the projection method for obtaining the point 803, for which the delay amount D.sub.2 is obtained, the method is not limited to this method. It is of course also possible to obtain a point 803 of at a position on a shape further closer to the actual wave face, and obtain the delay amount D.sub.2 for such a point 803.
Second Embodiment
(49) The ultrasonic imaging apparatus of the second embodiment of the present invention will be explained with reference to
(50) As shown in
(51) Thus, the correction operation part 610 may read a weight value from the weight table 1100, and used it in the step 1902 shown in
(52) The configurations and operation of the apparatus of this embodiment other than those explained above are the same as those of the first embodiment, and therefore explanations thereof are omitted.
Third Embodiment
(53) The ultrasonic imaging apparatus of the third embodiment of the present invention will be explained with reference to
(54) Although the ultrasonic imaging apparatus of the third embodiment has the same configurations as those of the ultrasonic imaging apparatus of the second embodiment (
(55) The weight changing module 1200 receives each value for weight value a, or the values for the whole table from the operator through the console 608, and substitutes them for values or table stored in the weight table 1100. The weight values a may be thereby changed to appropriate values depending on characteristics of a subject to be imaged, and therefore a delay amount (D) suitable for the subject may be obtained with the correction operation part 610. Since the other configurations and operation are the same as those of the second embodiment, explanations thereof are omitted.
(56) Further, the weight changing module 1200 may also be provided in the receiving beamformer 603 that is not provided with the weight table 1100. In this case, the weight changing module 1200 changes value of a parameter to be used for the calculation of the weight value a by the correction operation part 610 to a value received via the console 608. For example, it changes value of the constant a.sub.0 to be used for the operation of the equation (3) to a value received via the console 608. Since the value of a.sub.0 may be thereby changed depending on subject to a value suitable for the subject, a delay amount (D) suitable for the subject can be obtained by the correction operation part 610.
(57) Since other configurations and operation are the same as those of the first and second embodiments, explanations thereof are omitted.
Fourth Embodiment
(58) The ultrasonic imaging apparatus of the fourth embodiment of the present invention will be explained with reference to
(59) The receiving beamformer 603 of the ultrasonic imaging apparatus of the fourth embodiment performs beamforming according to the aperture synthesis method. As shown in
(60) In the receiving beamformer 603, an inter-transmission synthesis part 1300 is disposed, and in the main body 50, a frame memory 1302 is also disposed besides the beam memory 1301. The other configurations are the same as those of the third embodiment.
(61) The receiving beamformer 603 sets a plurality of scanning lines for each transmission, and carries out phasing addition of the received signals for a plurality of imaging points on each scanning line by using the delay amounts (D) obtained by the virtual sound source method-based delay amount calculation part 609 and the correction operation part 610. Phased outputs (rasters) are thereby obtained. A plurality of the obtained phasing outputs (rasters) are sent to the beam memory 1301 via the inter-transmission synthesis part 1300 and stored therein. The inter-transmission synthesis part 1300 reads out a plurality of phased results for a specific (the same) imaging point from a plurality of phased outputs (rasters) stored in the beam memory 1301, and synthesizes them (aperture synthesis). For example, as shown in
(62) The aperture synthesis image obtained from the synthesis processing is stored in the frame memory 1302, sent to the image processing part 605, and displayed on the image display part 607. The image processing part 605 displays the image obtained by the aperture synthesis on the image display part 607.
(63) Since the other configurations and operation such as those of the weight changing module 1200 are the same as those of the third embodiment and the first embodiment, explanations thereof are omitted.
(64) Since the receiving beamformer 603 of the present invention correct the delay amount (D.sub.1) obtained by the virtual sound source method to obtain the corrected delay amount (D), a phased output (raster) is obtained with good accuracy even for an imaging point at a position remote from the transmission focus 203 (scanning line 901), irrespective of whether the imaging point is inside the main beam or not. Therefore, a plurality of rasters are obtained for a wide area with one time of transmission. Accordingly, by storing them in the beam memory 1301, and carrying out the aperture synthesis with rasters obtained by another transmission, rasters of higher accuracy are generated, and used to generate an image.
Fifth Embodiment
(65) The ultrasonic imaging apparatus of the fifth embodiment of the present invention will be explained with reference to
(66) As shown in
(67) The directional angle information memory part 1400 receives transmission aperture P and transmission frequency f among the set parameters from the transmitting beamformer 602, obtains transmission directional angle θ.sub.a by calculation in accordance with the equation (4), and memorizes it.
[Equation 4]
Sin θ.sub.a=v/(f.Math.P) (4)
The symbol v represents acoustic velocity in a subject, and it is supposed that oscillators of the ultrasonic elements 600 have a rectangular shape.
(68) The transmission interval information memory part 1401 receives information on transmission interval of the transmission beam (ultrasonic beam) 104 in the direction along the probe array 100 from the transmitting beamformer 602, and memorizes it. The phasing area control part 1402 receives the transmission directional angle θ.sub.a and the transmission interval 503b from the directional angle information memory part 1400 and the transmission interval information memory part 1401, respectively, determines a phasing area 105b that defines area of imaging points on the basis of them, and specifies it for the virtual sound source method-based delay amount calculation part 609. The virtual sound source method-based delay amount calculation part 609 calculates delay amounts for imaging points in the specified phasing area 105b.
(69) The operation of the phasing area control part 1402 will be further explained with reference to
(70) According to this embodiment, as shown in
[Equation 5]
θ.sub.b=b.Math.θ.sub.a (5)
(71) The coefficient b is set by the phasing area control part 1402 with reference to an equation or table defined beforehand so that the minimum width 502b of the phasing area 105b becomes larger than the transmission interval 503b read from the transmission interval information memory part 1401, as shown in
(72) Since the receiving beamformer 603 of the present invention obtains the corrected delay amount (D) by correcting the delay amount (D.sub.1) obtained by the virtual sound source method, it obtains the delay amount (D) with good accuracy even for an imaging point at a position remote from the transmission focus 203. Therefore, phased output (raster) is obtained with good accuracy by using the corrected delay amount (D).
(73) By extending the phasing area 105b as shown in
(74) Since a plurality of rasters in a wide area may be obtained by one time of transmission by extending the phasing area 105b, by storing these in the beam memory 1301, and performing the aperture synthesis with them together with rasters obtained by other transmissions, an ultrasonogram of high precision can be generated for an area 501 (
(75) If the overlapping areas of the phasing areas 105b of a plurality of times of transmissions are small, block noises are generated in an image obtained after the aperture synthesis. Therefore, it is preferable to set the coefficient b so that the phasing areas 105b of at least three or more times of transmissions overlap with each other.
(76) Since the other configurations and operation are the same as those of the fourth embodiment, explanations thereof are omitted.
Sixth Embodiment
(77) The ultrasonic imaging apparatus of the sixth embodiment of the present invention will be explained with reference to
(78) The operator manually sets the value of the coefficient b of the equation (5) explained for the fifth embodiment, or finely tuning the coefficient b set by the phasing area control part 1402, by using the console 608.
(79) The value of the coefficient b is thereby appropriately adjusted according to scattering state of sonic waves, which varies depending on subjects, and therefore an appropriate phasing area 105b can be set. Accordingly, even by high frame rate imaging, an ultrasonogram of higher precision can be generated.
DESCRIPTION OF NUMERICAL NOTATIONS
(80) 100 Probe array 104 Transmission beam (ultrasonic beam) 105a and 105b Phasing area 202 Transmission caliber center point 203 Transmission focus 206 Transmitted sonic wave end 501 Phasing area of aperture synthesis for a plurality of times of transmission 503a and 503b Transmission interval 504 Omission region where phasing areas do not overlap 609 Virtual sound source method-based delay amount calculation part 610 Correction operation part 700 Transmitted sonic wave face at transmission focus depth 901 Virtual scanning line 1100 Weight table