ULTRASONIC MEASUREMENT METHOD, AND APPARATUS THEREFOR

Abstract: A method and apparatus for ultrasonic measurement transmit ultrasonic pulses having a plurality of different frequencies into an object under investigation and derive a mean value of attenuation coefficients indicative of attentua­ tion within the object, as well as the frequency dependence of the attenuation coefficient. If the frequency dependence of the reflection intensity from the object under investigation is constant or proportional to frequency within a range of measured frequencies, it is possible to obtain the attenuation coefficient and the coefficient indicative of its frequency dependence within a small area of the object by measuring the echo intensity of emitted ultrasonic pulses having three different frequencies. Furthermore, extending the above-­ mentioned processing over all amall areas of the entire object enables us to obtain a two-dimensional distribution image indicative of the attenuation coefficient and its frequency dependence.

Description: This invention relates to improvements in an ultrasonic measurement method for subjecting an object to an ultrasonic transmission and receiving reflected ultrasonic waves from the interior of the object to measure the acoustic characteristics of the object. More particularly, the invention relates to an ultrasonic measurement method and apparatus for providing information relating to attenuation that accompanies propagation of ultrasonic waves in the interior of an object. Ultrasonic measurement techniques find application widely in such fields as material testing, SONAR and medical diagnosis. In particular, ultrasound scanner systems for medical purposes have recently been developed. The principle of operation of an ultrasound scanner apparatus resides in use of a pulse-echo method and utilizes a phenomenon wherein an ultrasonic pulse transmitted into a living body, which is the object undergoing measurement, is reflected at a boundary where there is a difference in acoustic impedence. The reflected wave (echo) is received and processed to display a tomograph of the living body by a so-called B-mode method. The echo contains a variety of information such as the ultrasonic attenuation, acoustic impedence and propagation velocity of sound. With an apparatus of this kind, however, these various items of information cannot be separated from one another distinctly, so that the general practice is merely to display the amplitude of the echo. More specifically, the propagation velocity of sound is assumed to be constant and, with regard to attenuation ascribable to ultrasonic propagation, the value of echo amplitude arbitrarily corrected by a so-called STC (sensitivity time control) circuit is luminance-modulated, with the modulated amplitude being displayed on a cathode-ray tube as a tomograph. Accordingly, the tomograph obtained is a two-dimensional distribution at the acoustic impedence interface rendered qualitatively into a picture, so that the morphological information relating to the position and shape of the bioligical tissue inevitably forms the core of the information utilized. In other words, the state of the art is such that such biological tissue characteristics as degree of attenuation and prpagation velocity of sound are not measured. Several attempts at attaining attenuation information relating to biological tissue have been reported. As will be described below in further detail, an echo signal contains two types of information, namely attenuation due to propagation through biological tissue, and intensity of reflection at an interface or boundary where there is a difference in acoustic impedence. Both of these quantities are unknown. Therefore, isolating the affects of these two quantities is extremely difficult at the present time. If the reflected intensity is assumed to be independent of the frequency of the ultrasonic waves and ultrasonic waves having a plurality of different frequencies are transmitted and received with regard to the same portion of the object under measurement followed by measuring the sound pressure ratio of each frequency component of the echo, then it will be possible to eliminate the influence of the reflected intensity and derive an attenuation coefficient ascribable to propagation without the affects of the reflected intensity. The foregoing assumption holds in the case of an acoustic boundary having a sufficiently wide spread in comparison with the wavelength of the ultrasonic waves, e.g. in the case of a planar reflector. In actuality, however, a scatterer having a size approximately equivalent to or less than the wavelengths used can reside at the biological tissue. It is therefore difficult to consider that the foregoing assumption will always hold for the entirety of a bioligical tissue. In view of the foregoing circumstances, the applicants and one of the inventors of the present patent application has proposed, in Japanese Patent Application Laid-Open No.56-147082 and in JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 69 No 6, June 1981 pages 1838 to 1840, a method and apparatus for measuring attenuation information ascribable to ultrasonic propagation in a living body in which the influence of reflected intensity at an acoustic boundary is reduced. However, the ultrasonic attenuation information obtained according to the invention of this earlier application is not an attenuation coefficient but an equation of the first degree involving a mean value of ultrasonic attenuation coefficients over a prescribed interval. SUMMARY OF THE INVENTION An object of the present invention is to provide an ultrasonic measurement method and apparatus for measuring by approximation the attenuation coefficient of an object as well as the frequency dependence of the attenuation coefficient. According to the present invention, the foregoing object is attained by providing an ultrasonic measurement method and an ultrasonic measurement apparatus as defined in anyone of claims 1 to 4. In the method of the present invention, information relating to the attenuation coefficient of the object is derived from the mean value of attenuation coefficients measured on the basis of two different frequencies, and either the mean value of attenuation coefficients or a mean value of a function of attenuation coefficient and frequencies measured on the basis of three different frequencies inclusive of the first two frequencies is used. The information relating to the attenuation coefficient includes an attenuation coefficient of ultrasonic waves in the object, and the information relating to the attenuation coefficient includes the frequency dependence of the attenuation coefficient of ultrasonic waves in the object. The apparatus for practicing the foregoing ultrasonic measurement method includes means for transmitting ultrasonic pulses having a plurality of different frequencies into an object under investigation, means for detecting echos of the ultrasonic pulses reflected from within the object, means for displaying quantitative information obtained, and means for subjecting the detected echos, which will exhibit a plurality of frequencies, to information processing to obtain quantitative information relating to the ultrasonic characteristics of the object. This last-mentioned means obtains the mean value of attenuation coefficients having two different frequencies, and either the mean value of attenuation coefficients having three different frequencies inclusive of the first two frequencies or the mean value of the function of attenuation coefficient and frequencies which are said-three different frequencies inclusive of the first two frequencies, and derives information relating to the attenuation coefficient of the object on the basis of the mean values obtained. Other features and advantages of the present invention will be apparent from the following description taken in conjunction with the accompanying drawings, in which like reference characters designate the same or similar parts throughout the figures thereof. BRIEF DESCRIPTION OF THE DRAWINGS Figs. 1 and 2 are views useful in describing the fundamental principle of the present invention; Figs. 3, 4 and 5 are views useful in describing the calibration of an ultrasound field of ultrasonic probe used in an embodiment of the present invention; Fig. 6 is a block diagram illustrating an embodiment of an apparatus for practicing the ultrasonic measurement method of the present invention; Figs. 7 and 8 are waveform diagrams illustrating pulse waveforms useful in describing the operation of the embodiment shown in Fig. 6; Figs. 9(a) through 9(f) are flowcharts illustrating the operating steps of the apparatus embodied in Fig. 6; Fig. 10 is a view illustrating an example of an ultrasonic tomograph displayed by the apparatus embodied in Fig. 6; Fig. 11 (a) and (b) is a functional block diagram illustrating the functional construction of an arithmetic circuit in the apparatus embodied in Fig. 6; and Fig. 12 is a view useful in describing another display mode of the apparatus embodied in Fig. 6. DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of an ultrasonic measurement method according to the present invention will now be described in detail with reference to the accompanying drawings. The principle of the present invention will described in conjunction with Fig. 1. An ultrasonic pulse 10 transmitted into an object 4 by an ultrasonic probe 1 is reflected at each of two surfaces 2, 3 where acoustic characteristics are discontinuous, resulting in reflected pulses 30, 40, respectively, which are detected by the probe 1. If we let t represent the elapsed time from emission of the ultrasonic pulse to reception of the reflected pulse, the distance x from the ultrasonic probe 1 to a surface of ultrasonic characteristic discontinuity will be given by the following equation, where V is the prpagation velocity of sound: x = V·t/2 If we let I o (f) represent the intensity of the pulse 10 emitted at a frequency f, and let I(f, x) represent the intensity of a received reflected pulse, the following equation, indicative of an approximation, will hold: Taking the natural logarithm of both sides of Eq. (1) will result in the following: Here, α(f,x) is the attenuation coefficient ascribable to ultrasonic propagation within the object 4. In Eq. (1), f a(f,x) ·g(x) is the reflection intensity of the acoustic characteristic discontinuity surface with frequency dependence being taken into consideration. With regard to the value of a(f,x), a(f,x) = 0 will hold at an acoustic characteristic discontinuity surface sufficiently larger than the wavelength (λ= V/f ), whereas a(f,x) = 4 will hold at an acoustic characteristic discontinuity surface sufficiently smaller than the wavelength. Within a certain frequency range, therefore, a(f,x) will be constant, and the relation 0≦a(f,x)≦4 may be considered to hold for a living body. Further, x b(x) may be considered as being the effect of a weakening in the reflection intensity at the position of the ultrasonic probe owing to spread of the reflected wave. At a sufficiently broad acoustic characteristic discontinuity surface, b(x) = 0 will hold, whereas b(x) = -2 will prevail at a small acoustic characteristic continuity surface. In general, therefore, the relation -2≦b(x)≦0 will hold. In the foregoing considerations, the emitted ultrasonic wave 10 ideally is a fine pencil beam, and the effects of mutual interference among the waves reflected at a plurality of acoustic characteristic discontinuity surfaces located close to one another is neglected. The principal of carrying out information processing based on the echo intensity distribution to obtain quantitative information relating to the ultrasonic characteristics of the object 4 starts with Eq. (2). To facilitate the mathematical perspective, we will subject both sides of Eq. (2) to a first order partial differentiation with respect to x, assuming that the echo intensity I(f,x) is obtained continuously with regard to the distance x and frequency f. This will result in the following: Performing a first order partial differentiation with respect to the frequency f gives us: Multiplying both sides by f and arranging gives the following equation: Subjecting both sides of Eq. (3) to a second order partial differentiation with respect to nf, namely the log of the frequency, provides the following equation: If the ultrasonic attenuation coefficient of the object under investigation is proportional to the power β(x) of the frequency, α(f,x) = α o (x) f β(x) will hold, so that the left side of Eq. (5) will become as follows: The left side of Eq. (6) takes on the following form, noting that holds: Dividing Eq. (8) by Eq. (7) gives us: and dividing the square of Eq. (7) by Eq. (8) gives us: whereby the attenuation coefficient α(f,x) and the frequency dependence β(x) thereof can be found. That is, by using the following relation: derived from Eq. (9) and Eqs. (5), (6), as well as the relation α(f,x) = α o (x) f β(x) obtained from β(x) in Eq. (9) and α(f,x) in Eq. (10), α o (x) may be found from the equation α o (x) = α(f,x)/f β(x) . While the first term in the numerator and denominator of Eq. (11) is a measured quantity, the second term is not a measured quantity and is dependent upon the intensity of reflection at an acoustic characteristic discontinuity surface. This term is a source of error when obtaining β(x). If the coefficient a(f,x) indicating the frequency dependence of the reflection intensity is constant within the range of frequencies measured, that is, if holds, then a(f,x) will be a function solely of x (the assumption will be a(x) hereafter), so that Eqs. (5) and (6) may be written as follows, respectively: Therefore, Eq. (11) may be written as follows from Eqs. (12), (13): Likewise, Eq. (10) may be written as follows: It will be understood that these have fewer error terms in comparison with Eqs. (11) and (10). In other words, the error in Eq. (14) is solely in the denominator, and in Eq. (15) the error term is solely in the numerator. Furthermore, if the coefficient a(x) is constant and independent of the distance x, then, from Eq. (14) will not contain an error, as indicated by the following equation: In other words, if the coefficient a(f,x) indicating the frequency dependence of the reflection intensity is constant and independent of the measured range of frequency f and of the distance x, then the frequency dependence will be f a , and the frequency dependence β(x) of the attenuation coefficient α(f,x) can be calculated accurately by using Eq. (16). In this case, it is also possible to accurately calculate α(f,x), as shown by the following equation, based on Eq. (15): The case described above is one in which the echo signal intensity I(f,x) is obtained continuously with relation to the distance x and frequency f. In an actual measurement, however, the echo signal is obtained discretely based upon the given positions of the scatterers. Accordingly, some correction is necessary. More specifically, differences must be used in place of the differentials employed in Eqs. (5) and (6). As shown in Fig. 1, taking x₁ and x₂ as the positions at which the scatterers 2, 3 discretely reside, the difference in Eq. (2) between these two positions is obtained and a transformation is made, giving the following equation: To obtain an approximate expression corresponding to Eq. (4) by numerical differentiation, we take the difference in Eq. (18) at two frequencies f₁, f₃, multiply by and arrange to arrive at the following equation: To obtain an appoximate expression corresponding to Eq. (5), we multiply by (f₃ + f₁)/2 to arrive at the following equation: ,where From α(f,x) =α o (x)f β(x) we have: which gives us the following: Likewise, to obtain an approximate expression corresponding to Eq. (6) by numerical differentiation, we use data regarding three frequencies f₁, f₂, f₃ (f₁<f₂<f₃), take a difference quotient of second order with relation to the frequency logarithms ℓnf₁, ℓnf₂, ℓnf₃, and make a transformation to obtain the following equation: where and A = 2/{ℓn(f1/f2) ℓn(f2/f3)ℓn(f1/f3)} Therefore, we have the following: Dividing Eq. (23) by Eq. (21) gives the following equation: Thus, Eq. (24) corresponds to Eq. (9) and gives an approximate expression of a mean value of β(x) between x₁ and x₂. H₃/H₂ takes on the following form from Eqs. (20), (22): where Further, from Eqs. (21) and (23) we have: Thus, Eq. (26) corresponds to Eq. (10) and gives a mean value of α(f,x) between x₁ and x₂ and between f₁ and f₃. (H₂)²/H₃ takes on the following form from Eqs. (20), (22): The first term in the denominator of Eq. (25) and the first term in the numerator thereof are measured quantities, as in Eq. (11), but the second terms, namely D₂ and D₃, arise owing to frequency dependence ascribable to the intensity of reflection at the acoustic characteristic discontinuity surface, also as in Eq. (11). These represent error terms with respect to β(x). If the coefficient a(f,x) of the frequency dependence f a(f,x) of the reflection intensity is constant within the range of frequencies measured, namely if a(f₁,x₁) = a(f₂,x₁) = a(f₃,x₁) = a(x₁) and a(f₁,x₂) = a(f₂,x₂) = a(f₃,x₂)= a(x₂) hold, then a will be a function solely of x and D₃ will be zero. Accordingly, Eqs. (20), (22) may then be written as follows: Therefore, Eq. (25) takes on the following form from Eqs. (28) and (29): And Eq. (27) takes on the following form from Eqs. (28) and (29): It may be understood that these two equations have fewer error terms in comparison with Eqs. (25) and (27). More specifically, the error term in Eq. (30) is solely [a(x₂)-a(x₁)]·ℓn(f₃/f₁) in the denominator, and the error term in Eq. (31) is solely [a(x₂)-a(x₁)]·ℓn(f₃/f₁) in the numerator. If the coefficient a(x) of the frequency dependence f a(x) of the reflection intensity is constant and independent of the distance x, then, from a(x₂) = a(x₁) , Eq. (30) may be written as the following equation, in which the error is reduced even further: In other words, if the coefficient a(f,x) indicating the frequency dependence of the reflection intensity is constant and independent of the frequency f within the measured range and of the distance x, then a(f,x) = a will hold and, in accordance with Eq. (32), the mean value β(x) of the frequency dependence of the attenuation coefficient α(f,x) can be measured with even greater accuracy. In this case, in accordance with Eq. (31), it is also possible to measure α(f,x) with greater accuracy as shown by the following equation: In this case, α o (x) also is obtained with greater accuracy. Thus, as set forth above, ultrasonic pulses having frequencies of f₁, f₂, f₃ are projected into the object and α, α o , β can be calculated from the reflection intensity of these pulses. In Eqs. (20), (22), D₂ and D₃ are error terms which cannot be found from measured values. In Eq. (20), D₂ vanishes when a(f,x) is constant and independent of f₁, f₃, x₁, x₂. In Eq. (22), D₃ vanishes when a(f,x) is constant with regard to the frequencies f₁, f₂, f₃ even if it is dependent upon x₁, x₂. In cases where the foregoing conditions do not hold, however, the relative error can be evaluated in the following way. Specifically, since 0≦a(f,x)≦4 holds, we may write: Therefore, we have the following relation: This gives the relative error ascribable to D₂ in Eq. (20). Likewise, we have: This gives the relative error ascribable to D₃ in Eq. (22). In conclusion, the error terms D₂, D₃ can be neglected if, from Eqs. (34), (35), the value of the following is sufficiently large: As set forth above, attentuation coefficients α(f,x), α o (x) and their frequency dependence β(x) can be approximately measured by three different frequencies f₁, f₂, f₃. In particular, if the coefficient a(f,x) is constant and independent of the observed range of frequency f and of the distance x, then it is possible to measure α(f,x), α o (x) and (x) with greater accuracy. IN actuality, however, the coefficient a(f,x) indicative of frequency dependence of the reflection intensity is an unknown quantity. Therefore, the actual calculation of α(f,x), α o (x) and β(x) is performed from H₂ and H₃ in the following equations based on measured quantity I(f,x): In the following, we will discuss the case in which the ultrasonic beam are emitted into the object comprising the tissues of which coefficient a(f,x) is proportional to frequency f, namely a(f,x) = a o (x).f holds. We start with Eq.(2) again. To facilitate the mathematical perspective, we assume the echo intensity I(f,x) is obtained continuously with regard to the distance x and frequency f. Performing a first order partial differentiation with respect to f.ℓn f gives us: Subjecting both sides of Eq.(36) to a first order partial diffentiation with respect to x and arranging give us the following equation: ,where Subjecting both sides of Eq.(37) to partial differentiation with respect to fℓnf provides the following equation: ,where Assuming the following equation as described previously holds: α(f,x)=α o (x).f β(x) (39) We get the equation as follows: Thus, the left side of Eq.(38) will become as follows: Dividing Eq.(41) by Eq.(40) gives us; We get the following equation from Eq.(42) and Eq.(40). Thus, the attenuation coefficient α(f,x) and the frequency dependence β(x) of the coeffieient α(f,x) can be calculated as follows by using Eq.(42), Eq.(37) and Eq.(38): From Eq.(43), Eq.(37), (38), we have: From β(x) in Eq.(42) and α(f,x) in Eq.(43), we can calculate α o (x) by using the equation α(f,x)=α o (x).f β(x) . While the first term in the numerator and the denominator of Eq.(44) is a measured quantity, the terms on and after the second are derived from the intensity of reflection at an acoustic charasteristic discontinuity surface and are a source of error when obtaining β(x). If the coefficent a(f,x) is constant within the measured range of frequencies , that is, if a(f,x)=0 hold, then a(f,x) will be a function solely of x, so that Eqs. (37) and (38) may be written as follows, respectively: Therefore, Eq.(44) may be written from Eqs.(46) and (47) as follows: Likewise, we have the following equation from Eq.(45): It will be understood that these Eqs.(48),(49) have fewer error terms in comparison with Eqs.(43) and (44). If a(x) is constant and independent of distance x, Eq.(48) will not contain an error, as indicated by the following equation: In other words, if the coefficient a(f,x) indicating the frequency dependence of the reflection intensity is constant and independent of the measured range of frequency f and the distance x, then the frequency dependence will be f a ,and the frequency dependence β(x) of the attenuation coefficient α(f,x) can be calculated accurately by using Eq.(50). In this case, it is also possible to accurately calculated α(f,x) as shown by the following equation, based on Eq.(49): Next, we will discuss the case in which the coefficient a(f,x) is not constant but proportional to frequency f within the range of frequencies measured. In such case, we can assume that a(f,x)=a₀(x).f , thus obtaining the following Eqs.(37) and (38): Therefore, Eq.(44) will become from Eqs.(52) and (53) as follows: likewise, Eq.(43) will become: We will find out that Eqs.(54) and (55) have fewer error terms in comparison with Eqs.(44) and (45). In other words, the error term in Eq.(54) is solely in the denominator, and the error term in Eq.(55) is solely in both the denominator and the numerator. The case described above is one in which the echo signal intensity I(f,x) is obtained continuously with relation to the distance x and frequency f. In an actual measurement, however, echo signal is obtained discretely based upon the given positions of the scatterers. Accordingly, some correction is necessary. More specifically, differences must be used in place of the differentials employed in Eqs.(37) and (38). As shown in Fig. 1, taking x₁ and x₂ as the positions at which the scatterers 2, 3 discretely reside, the difference in Eq. (2) between these two positions is obtained and a transformation is made, giving the following equation as well as Eq.(18): The equation corresponding to Eq.(37) is as follows: ,where Likewise, the equation corresponding to Eq.(38) is as follows: ,where We can approximate the left side of Eq.(57) to the following equation, noting α(f,x)=α o (x)f β(x) and Eq.(40): Likewise, we approximate the left side of Eq.(59) to the following equation, noting Eq.(41): Multiplying both sides of Eq.(61) by and arranging gives us the following equations in stead of Eqs.(57) and (58): ,where Multiplying both sides of Eq.(62) by and arranging give us the following equations in stead of Eqs.(59) and (60). ,where The equation which gives β(x) and corresponds to Eq.(42) becomes as follows: Likewise, the equation which gives α(f,x) and corresponds to Eq.(43) becomes as follows: If the coefficient a(f,x) of the frequency dependence f a(f,x) of reflection intensity is constant within the range of frequencies measured, namely if a(f1 , x₁)=a(f₂,x₁)=a(f₃,x₁)=a(x₁) and a(f₁,x₂)=a(f₂,x₂)=a(f₃,x₂)=a(x₂) hold, then a(f,x) will be a function solely of x, thus Eqs.(63) and (64) will become as follows: What K₂' implies is the integral value, namely mean value, of the attenuation coefficient α(f,x) for two different frequencies f₁, f₃ between x₂ and x₁, and what K ' / 3 implies is the mean value of the function of three different frequencies f₁, f₂, f₃, and the attenuation coefficient for said three different frequencies between x₂ and x₁. Accordingly, β(x) takes on the following form Eqs.(65),(67) and (68) Likewise, α(f,x) takes on the following form from Eqs.(66), (67) and (68): We find out that these two equations,69) and (70) have fewer error terms in comparison with Eqs.(65) and (66) respectively. If a(x) is constant and independent of the distance x, namely a(x₁)=a(x₂), Eq.(69) may be written as the following equation, in which the error is reduced even further: ,where Likewise, we will obtain α(f,x) in the following form: We will again discuss the case in which a(f,x) is not constant but proportional to frequency f within the measured range of frequency as well as discussed in the case in which the echo signal is obtained continuously. If a(f,x)=a₀(x).f holds within the measured range of frequency, Eqs.(63) and (64) will become as follows: Accordingly, β(x) takes on the following form from Eqs.(65) and (73): α(f,x) will become as follows from Eqs.(66) and (74): We can find out that Eqs.(75) and (76) have fewer error terms in comparison with Eqs.(65) and (66) respectively. As set forth above, attenuation coefficient α(f,x), α₀(x) and their frequency dependence β(x) can be approximately measured by three different frequencies f₁,f₂, and f₃. In particular, if the coefficient a(f,x) is constant and independent of the measured range of frequency f and the distance x, or a(f,x) is proportional to frequency f within the range, then it is possible to measure α(f,x),α₀(x) and β(x) with great accuracy. In actuality, however, because the a(f,x) indicative of frequency dependence of reflection intensity is an unknown quantity, the calculation of α(f,x), α₀(x) and β(x) is performed from K ' / 2 and K ' / 3 in the following equations based on the measured quantity I(f,x); We have discussed the theoretical background of the present invention in the foregoing description, then we will describe the present invention more concretely. The ultrasonic beam emitted from the ultrasonic probe 1 and transmitted through a substance such as water which attenuates the beam very little develops a difference in sound field owing to the resonance frequency and the aperture size of the probe. The sound field has the approximate form shown in (a) of Fig. 2. The intensity of the sound field on the central axis of the ultrasonic probe 1 varies in accordance with the distance x from the probe, as illustrated in (b) of Fig. 2. It should be noted that the longitudinal axis in Fig. 2(b) indicates the intensity I x at the distance x normalized by maximum intensity I o . To make practical use of the apparatus, therefore, a correct measurement cannot be performed unless the apparatus is properly calibrated. To this end, a variation in sound pressure is measured beforehand using a standard medium, and the echo amplitude (sound pressure) from the object under investigation is divided by a standard sound pressure to achieve normalization. This eliminates the influence of the probe sound field characteristic and makes it possible to obtain a more universal value for the measured attenuation. An appropriate method of measurement using a standard medium will now be described with reference to Fig. 3. As shown, a perfect reflector 102 made of stainless steel or the like is submerged in a bath of degased water 100. The amplitude of an echo received from the perfect reflector 102 serves as a standard sound pressure. A standard sound pressure curve is obtained as shown in Fig. 4 by varying the distance between the ultrasonic probe 1 and the perfect reflector 102 and measuring the amplitude of the echo received at each distance. For an apparatus suited to practical use, the standard sound pressure curve would be measured in advance and stored beforehand within the apparatus. If it is desired to remeasure the standard sound pressure curve, the apparatus should be fabricated to include a step-like perfect reflector 102a as shown in Fig. 5. Specifically, if the ultrasonic probe 1 is arranged so as to be scanned horizontally in the direction of the arrows by a scanning section 8, the distance between the ultrasonic probe 1 and the reflective surface of the reflector 102a can be varied in a stepwise manner. Here the amplitude of the echo from each "step" of the perfect reflector 102a would be measured and the measured amplitudes stored in memory successively to obtain the standard sound pressure curve. Let this sound pressure curve be represented by J(f,x). In the illustrated embodiment, J(f,x) is stored in the internal memory of a control circuit in a manner to be described below. Therefore, calibrated echo intensity represented by I and(f,x), I (f,x)=I(f,x)/J(f,x) holds. Accordingly, when we calculate the attenuation coefficient etc, we have to substitute I and(f ,x) for I(f,x). Described next will be an embodiment of the present invention illustrated in toe block diagram of Fig. 6. The illustrated embodiment is for practicing the ultrasonic measurement method of the present invention and is composed of a transmitting system and a receiving system. The transmitting system includes a transmitting circuit 5 connected to the ultrasonic probe 1, which is set at the surface of an object under investigation, such as a living body 4, for the purpose of effecting a mutual conversion between an electric signal and ultrasonic waves. The receiving system includes a receiving circuit 7, which is also connected to the probe 1 and to the transmitting circuit 5, a logarithmic amplifier 18 which recieves the output of the receiving circuit 7, a detector 19 which receives the output of the logarithmic amplifier 18, and an STC circuit 20 which receives the output of the detector 19. The ultrasonic probe 1 is scanned by the scanning section 8 under the control of a control circuit 17. The control circuit 17, which is shown in Fig. 11(b), controls the entirety of the ultrasonic measurement apparatus of the present invention and includes a CPU 400 of a mini- or microcomputer, a ROM/RAM 401 in which the control program of the ultrasonic measurement apparatus is stored, a memory 402 which the CPU 400 uses as a working storage device for intermediate data, and an I/O port 403 serving as an input/output control interface between the CPU 400 and the scanning section 8, an arithmetic circuit 11 and other circuits. The ultrasonic echo received from the interior of the living body 4 is weakened by attenuation and spread of the ultrasonic beam. The purpose of the STC circuit is to vary the gain of the receiving system timewise during each sweep to apply a correction to, e.g., attenuation ascribable to distance. The transmitting circuit 5 supplies the probe 1 with a wide band drive pulse that undergoes sharp attenuation, as shown in Fig. 7. The probe 1 has a wide-band characteristic and preferably comprises a piezo-electric polymer [polyvinylidene fluoride (PVDF)] transducer, a composite transducer consisting of a piezo-electric polymer and inorganic matter, or a PZT transducer equipped with an acoustic matching layer. The result of using such a probe is that the probe will transmit a wide-band ultrasonic pulse of the kind shown in Fig. 8 into the interior of the object 4. It goes without saying that the frequencies f₁, f₂, f₃ preset for calculating α, α o and β fall within the frequency band of this wide-band ultrasonic pulse. An ultrasonic echo reflected or scattered at a surface (e.g., surface 2 or 3 in Fig. 1) within the object 4 having a discontinuous acoustic characteristic is received by the same probe 1 and applied thereby to the receiving circuit 7. Here the echo is converted into an amplified echo signal (A-mode signal) and then applied to an A/D converter 9 for conversion into a digital quantity stored in a memory 10. As performed in the prior art, a B-mode picture is obtained by scanning the probe 1 above the object 4 by means of the scanning section 8, and collecting the A-mode signals. Many scanning methods are available, including a mechanical sector scan, a compound scan, a linear electronic scan and a sector electronic scan. The details of these scanning methods do not have a direct bearing upon the present invention and will therefore not be described here. The echo signal obtained by such scanning enters the logarithmic amplifier 18 via the receiving circuit 7. The signal is logarithmically amplified by the amplifier 18, is detected by the detector 19 and is then subjected to an STC correction by STC circuit 20. The STC-corrected echo signal is stored in a memory 22 and is developed by a picture output amplifier circuit 14 into a B-mode picture for display on a display unit 15. The circuitry for producing the B-mode picture output is well-known and need not be described here. If the writing (entry) of data into the memory 22 is halted and the information within the memory 22 is output repeatedly (i.e., played back), a "frozen" B-mode picture is developed for display by the display unit 15. On the other hand, the method in which the B-mode picture can be displayed simultaneously with the scanning, is actually available in real-time ultrasound scanner apparatus using a linear electronic scan method. In this apparatus, we may utilize so-called "freezing function" to obtain desired tomograph. More particularly, an observer can make a freeze switch (not shown in fig. 6) function so that writing the data into the memory 22 may be prohibited and the data in memory 22 may be developed to be displayed in display unit 15 as a stationary picture continuously. Writing the data into memory 22 will restart by releasing the freeze switch, thus providing a dynamic picture to be displayed in display unit 15. As shown in Fig. 10, the frozen picture, indicated at 200, has a region of interest 202, in which α(f,x), α o (x) and β(x) are desired to be measured, designated thereon by a circuit 16 for setting the region of interest (ROI). This is performed by manipulating a control member (not shown) of the setting circuit 16 to designate the intervals of the area of interest by a frame 202, as shown in Fig. 10. This method is similar to a method of measuring the distance between two points by so-called caliper measurement. In accordance with a predetermined algorithm shown in Figs. 9(a) through 9(f), the A mode signal stored in the memory 10 is processed by the arithmetic circuit 11 to obtain β(x), α o (x) and α(f,x), the values whereof are then stored in a memory 12. The values of β(x), α o (x) and α(f,x) are selected from the memory 12 when required, the selected value is converted into an analog quantity by a D/A converter 13, and the resulting analog signal is applied to the picture output amplifier circuit 14 where the signal is amplified at delivered to the display unit 15 as a visible picture. The operation of the ultrasonic measurement apparatus of the present invention will now be described with reference to the flowcharts of Figs. 9(a) through 9(f). The operational flow is divided into two portions, namely a portion for processing up to the B-mode display and a portion for calculation of α(f,x), α o (x) and β(x). The program starts with processing for the B-mode display. The CPU 400 drives the transmitting circuit 5 via a control line 26 at a step 300. As mentioned above, the probe 1 is driven by the transmitting circuit 5 and projects ultrasonic waves into the object 4, assumed here to be a living body. The CPU 400 sets the scanning section 8 into operation via a control line 25 at a step 302. The scanning section 8 responds by sending a synchronizing signal to the probe 1 and CPU 400. As a result, the sweep of the probe 1 and the operation of the scanning section 8 and CPU 400 are synchronized. At the same time that the probe 1 starts to be swept in step 302, the CPU 400 executes a step 304 to instruct the A/D converter circuits 9, 21, via control lines 27, 24, respectively, to begin an analog-to-digital conversion operation. The CPU 400 simultaneously executes a step 306, in which the STC circuit 20 is placed in operation via a control line 28 in order to perform an STC correction. It goes without saying that the operations indicated by steps 300 through 306 are performed substantially simultaneously. Thus, one sweep of an A-mode signal for a B-mode display is stored in the memory 22, and an A-mode signal for calculating α(f,x), α o (x) and β(x), which is an object of the invention, is stored in the memory 10 (step 308). The foregoing sweeping operation is performed for a full set of scanning lines. In other words, it is determined at a step 312 whether a full scan has beem completed, with steps 300 through 312 being repeated when the decision in step 312 is NO. As previously described, when the freeze-switch functioned, the flow goes from step 312 to step 314 unconditionally. In this manner, A-mode signals for the entire scanning area are stored in the memories 10, 22. When a full scan is completed, i.e., when the answer is YES in step 312, the A-mode signals stored in the memory 22 are luminance-modulated by the picture output amplifier 14 and displayed on the display unit 15 (step 314). It should be noted that a control line 29 carries a synchronizing signal for the display unit 15. Described next will be the arithmetic circuit 11 for calculating α(f,x), α o (x) and β(x), the control circuitry therefore, and the calculation algorithm. Fig. 11(a), (b) is a functional block diagram showing the arithmetic circuit 11 and a portion of the control circuit 17. Illustrated are both the CPU 400, ROM/RAM 401, memory 402, I/O port 403 of the control circuit 17 and the various components of arthmetic circuit 11. The role played by control circuit 17 is as described above. Though the details are not specifically shown in Fig. 11, the elements constituting the arithmetic circuit 11 operate under the control of the CPU 400. The flow of data in Fig. 11 is indicated by the bold arrows, while the flow of control signals is shown by the slender arrows. The algorithm for processing executed by the arithmetic circuit 11 under the control of the CPU 400 will now be described in conjunction with the flowcharts of Figs. 9(b) through 9(f). The flowcharts are shown on the left side of Figs. 9(b) through 9(f), while the results of associated processing are indicated in the form of graphs or the like on the right side at positions corresponding to the relevant process steps. As described above, the CPU 400, at the steps 300 through 314 for scanning the probe 1 to produce the B-mode display, causes the A/D converter circuit 9 to convert the echo signal from the receiving circuit 7 into a digital signal, and stores the digital signal in the memory 10. When an entire scan is completed, therefore, echo signals for an entire zone are stored in the memory 10 in digital form. Meanwhile, the display unit 15 presents a B-mode display. As shown at step 316, the CPU 400 applies a pixel address of an area designated by the region of interest setting circuit 16 to a designated interval distance calculating unit 55 as a scanning line interval and time interval, as shown in Fig. 10. In other words, scanning lines k, ℓ and arrival times t i , t j of the ultrasonic waves are stored in the calculating unit 55. The program moves to step 318, where a standard sound pressure curve J(f,x) in the memory 402 of the CPU 400 is stored in a portion of the memory 10. As mentioned above, the standard sound pressure curve J(f,x) stored in the memory 402 may be remeasured through the method shown in Fig. 5 to update the contents of the memory. Next, the CPU sets an adding and averaging unit 50 into operation at a step 320. The adding and averaging unit 50 then reads, out of the memory 10, the A-mode signals in the interval designated by the values of k, ℓ, t i , t j stored in the designated interval distance calculating unit 55, adds and computes the mean of these signals, produces a signal 1000 indicative of the result and stores the signal in a memory 60 at a step 322. For the purpose of later processing, the A-mode signal 1000 stored in the memory 60 is also stored in a memory 61 at a step 324. The CPU 400 then starts a detection unit 51 at a step 326. The detection unit 51 responds by subjecting the A-mode signal stored in the memory 60 to detection processing. An A-mode signal 1001 resulting from the detection processing is stored in the memory 60 at a step 328. Next, the CPU starts a high-frequency component extracting unit 52 at a step 330. The extracting unit 52 sends the detected A-mode signal 1001 through a high-pass filter to produce an A-mode signal 1002 from which the low frequency component has been removed (step 332). Then, at a step 334, a portion of the A-mode signal 1002 below a predetermined level (L o ) is cut off to remove noise, thereby resulting in an A-mode signal 1003 which is to be stored in the memory 60. Next, the CPU 400 starts a maximum value detecting unit 53 at a step 336, the unit 53 responding by detecting the position (time) of the maximum value of the A-mode signal 1003. The maximum values detected are 2000, 2001, 2002 and 2003 (step 338). Thus, in the present case, the designated interval contains four points at which the acoustic characteristic is discontinuous. These points appear at times (positions) t o , t₁, t₂, t₃. These values are stored into memory 60. The CPU 400 then starts an echo intensity calculating unit 54 at a step 340. The echo intensity calculating unit 54 is adapted to find the maximum and minimum values from among the values at t o , t₁ t₂, t₃ in the time interval t i -t j stored in the designated interval distance calculating unit 55. In the example of Fig. 9(e), these values are located at t o and t₃ (step 342). Next, a Hamming window of a predetermined width is applied with regard to the A-mode signal corresponding to t o and t₃, a fast Fourier transform is performed to obtain I(f), and I(f₁,x₁), I(f₂,x₁), I(f₃,x₁), I(f₁,x₂), I(f₂,x₂), I(f₃,x₂) corresponding to f₁, f₂, f₃ are found (step 344). Next, the CPU 400 places a normalizing unit 56 into operation at a step 346. In order to calibrate the sound field characteristic, the normalizing unit 56 divides I(f,x), obtained in step 344, by the referential value J(f,x) stored in the memory 10, thus to obtain I(f,x) (steps 348, 350). The program then moves to a step 352, where the CPU starts an attenuation coefficient calculating unit 57. At step 354, using X₁=C₀t₀, x =C₀t₃(where C₀=1530 m/s is a mean value of propagation velocity in living body), the calculating unit 57 calculates the mean value H₂ or K₂' on the basis of Eqs. (28') or (73') according to the frequency dependence of the coefficient a(f,x). At step 356, the calculating unit 57 calculates the value H₃ or K₃' on the basis of Eqs.(29) or (74). From the values of both H₂ or K₂' and H₃ or K₃', the calculating unit 57 calculates β(x), α(f,x), according to Eqs.(32), (33) or (75), (76) at step 360, thus providing α₀(x) from the equation α o (x)=α(f,x)/f β(x) at step 362. As described above, we substitute I(f,x) for I(f,x) in the equations appeared at stop 354 thruough 356. Next, at a step 364, the calculated values of H₂, H₃, β(x), α(f,x) and α o (x) are stored in the memory 12. Thus, the sought values of α, β, and α o are obtained. In steps 368, 370, the values of either H₂, H₃, or K₂', K₃' in addition to β(x), α(f,x) and α o (x) stored in the memory 12 are delivered to the display unit 15 through the picture output amplifier 14. These values are displayed by the display unit 15 together with the B-mode picture read out of the memory 21. Though α(f,x) is expressed in Neper/cm units, a conversion into dB/cm units can be made within the arithmetic circuit 11 if necessary. Alternatively, it is also possible to display β(x). In addition, α o (x) can be displayed in dB/cm MHz units. It should be noted that with a second display mode presented by the inventive apparatus, rather than measuring the attenuation coefficient solely within the designated interval, H₂, H₃, K₂', K₃', β(x), α o (x) and α(f,x) per unit pixel (interval) with respect to the entire picture can be measured and the distribution thereof displayed as an image. This becomes feasible if an algorithm similar to that described above is applied with respect to each unit pixel Δ s of the image, as shown in Fig. 12. Though the distribution image naturally has poorer spacial resolution in comparison with the B-mode image, a fundamental difference is that attentuation coefficient-related information absolutely unmeasurable with the prior art can be obtained as the displayed information. In the illustrated embodiment, a wide-band probe is used as the probe 1. However, as described in the specification of Japanese Patent Application Laid-Open No. 56-147082, it is permissible to employ a probe having a plurality of different frequency bands. Further, to facilitate the description of the illustrated embodiment, the elements constituting the arithmetic unit 11 have been described as if they were constituted by hardware. Though hardware may of course be used, the functions of these elements can also be implemented by the software of the CPU 400. According to the present invention as described and illustrated above, the intensity of an echo from an object under investigation is measured at a plurality of different frequencies, whereby the attenuation coefficient of the object as well as the frequency dependence of the attenuation coefficient can be measured by approximation. It is also possible to obtain a distribution image of the values obtained. Therefore, unlike the conventional ultrasonic tomographic measurement method and apparatus with which only morphological information can be obtained, the present invention makes it possible to obtain quantitative information relating to attenuation within the object.

Claims: An ultrasonic measurement method for obtaining an attenuation coefficient α of an object under investigation and/or the frequency dependence β thereof, where α = αo(f)·f β(x) holds, comprising the steps of: (a) transmitting ultrasonic pulses having a plurality of different frequencies into the object under investigation, (b) detecting echos of the ultrasonic pulses reflected from an area of interest within the object, (c) subjecting the detected echos, which exhibit a plurality of frequencies to a digital processing to extract the echo intensities I(f,x) corresponding to three different frequencies f₁, f₂, f₃ within the plurality of transmitted frequencies for the distances x₁, x₂ of two points of said area of interest, characterized by the steps of: (d) calculating first and second mean values of the attenuation coefficients expressed by where A is expressed by (e) and deriving the attenuation coefficient α and/or the frequency dependence β, based on the following equations: An ultrasonic measurement method for obtaining an attenuation coefficient α of an object under investigation and/or the frequency dependence β thereof, where α = αo(x)·f β holds, comprising the steps of: (a) transmitting ultrasonic pulses having a plurality of different frequencies into the object under investigation, (b) detecting echos of the ultrasonic pulses reflected from an area of interest within the object, (c) subjecting the detected echos, which exhibit a plurality of frequencies to a digital processing to extract the echo intensities I(f,x) corresponding to three different frequencies f₁, f₂, f₃ within the plurality of transmitted frequencies for the distances x₁, x₂ of two points of said area of interest, characterized by the steps of: (d) calculating first and second mean values of attentuation coefficients expressed by where B is expressed by (e) and deriving the attenuation coefficient α and/or the frequency dependence β, based on the following equations: An ultrasonic measurement apparatus for practicing the ultrasonic method according to claim 1, comprising transmission means (5) for transmitting ultrasonic pulses having a plurality of different frequencies into an object (4) under investigation, detection means (7) for detecting echos of the ultrasonic pulses reflected from an area of interest within the object (4), processing means (11) for subjecting the detected echos which exhibit a plurality of frequencies to extract the echo intensities I(f,x) corresponding to three different frequencies f₁, f₂,f₃ within the plurality of transmitted frequencies for the distances x₁, x₂ of two points of said area of interest, characterized by means for calculating first and second mean values of the attenuation coefficient expressed by where A is expressed by and for deriving the attenuation coefficient α of the object and/or the frequency dependence β thereof, from said mean values of attenuation coefficients on the basis of the following equations: and display means (15) for displaying said quantitative information in addition to a tomograph. An ultrasonic measurement apparatus for practicing the ultrasonic method according to claim 2, comprising transmission means (5) for transmitting ultrasonic pulses having a plurality of different frequencies into an object (4) under investigation, detection means (7) for detecting echos of the ultrasonic pulses reflected from an area of interest within the object (4), processing means (11) for subjecting the detected echos which exhibit a plurality of frequencies to extract the echo intensities I(f,x) corresponding to three different frequencies f₁, f₂,f₃ within the plurality of transmitted frequencies for the distances x₁, x₂ of two points of said area of interest, characterized by means for calculating first and second mean values of the attenuation coefficient expressed by where B is expressed by and for deriving the attenuation coefficient α of the object and/or the frequency dependence β thereof, from said mean values of attenuation coefficients on the basis of the following equations: and display means (15) for displaying said quantitative information in addition to a tomograph.

Inventors: FUJII, TADASHI, HAYAKAWA, YOSHINORI
Applicants: HAYAKAWA, YOSHINORI, TERUMO KABUSHIKI KAISHA TRADING AS TERUMO CORPORATION
Application year: 1984
Main classification: G01S15
IPC codes: G01S7, G01N29, A61B8, G01S15