Ultrasound diagnostic equipment for characterising tissue by analysis of backscatter

Abstract: Ultrasound diagnostic equipment comprises a tissue power unit (5) for calculating tissue scattering power of tissues of an object of measurement in an organ, a blood power unit (6) for calculating scattering power of blood in the organ and estimating blood scattering power at a portion of the tissues of the object of measurement by using the calculated scattering power of the blood, an integrated backscatter unit (7) for normalizing the tissue scattering power by using the estimated blood scattering power and calculating an integrated backscatter (IB) by the normalized tissue scattering power, and a display unit (9) for displaying the integrated backscatter (IB). Therefore, the integrated backscatter (IB) of the tissues can be quantitatively estimated without being influenced by attenuation. Furthermore, this estimated quantitative integrated backscatter (IB) of the tissues is displayed on the display unit accompanied with a B-mode image, an electrocardiogram, and the like, so that diagnostic accuracy can be increased.

Description: BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to ultrasound diagnostic equipment, more particularly, to ultrasound diagnostic equipment for displaying scattering coefficients b and n of tissues of an object to be measured in an organ by using scattering power spectrum of blood. In diagnostic techniques using ultrasound waves, it has become strong that needs to improve diagnostic accuracy by quantitative diagnosis of tissues. Now, acoustic characteristics of tissues, especially, attenuation characteristics and scattering characteristics, are used to indicate the characteristics of the tissues. The present invention relates to an equipment for calculating coefficients b and n quantitatively and exactly and for displaying them, when defining scattering frequency characteristics of the acoustic characteristics to the following equation. S ( f ) = b f n 2. Description of the Related Art Conventionally, when diagnosing a heart by using ultrasound waves and detecting scattering power of a cardiac wall by synchronising with a heartbeat, it is known that a change in ultrasound scattering power (IB) of a cardiac wall portion of a myocardial infarction is different from that of a normal cardiac wall portion. Futhermore, conventionally when diagnosing a heart by using ultrasound waves and calculating scattering coefficients, it is known that the scattering coefficients b and n of a cardiac wall portion of a myocardial infarction are different from those of a normal cardiac wall portion. However, when practically diagnosing a heart from a body surface by using ultrasound waves, the ultrasound waves are attenuated from the body surface to a measurement portion. In addition, the degree of attenuation from the body surface to each measurement portion is different from the others, and thus an absolute value of the obtained result cannot exactly indicate a characteristic of the portion. In order to solve the above problem, a method to measure the scattering power, which has high accuracy and does not depend on the attenuation (attenuation characteristics), is required. Furthermore, in order to solve the above problem, a method to measure the scattering coefficients, which has high accuracy and does not depend on the attenuation, are also required. EP-A-0 349 321, which falls under Article 54(3) EPC, discloses how to calculate the blood backscatter attenuation at a distance d by measuring the blood backscatter power at a point located at said distance d. This document discloses the embodiments as represented in figs.1 to 7 and in the corresponding passages of the description of this patent specification. These embodiments referred to as the first aspect are therefore comprised in the prior art under Article 54(3) EPC. It is an object of the present invention to provide ultrasound diagnostic equipment for calculating scattering coefficients b and n of tissues of an object of measurement in an organ exactly by normalizing scattering power spectrum of the tissues of the object of measurement by using scattering power spectrum of blood which has definite ultrasound scattering characteristics, and for displaying them. These objects are achieved by the features of claim 1. The tissue power spectrum unit may comprise a power spectrum unit, supplied with the received signals, for calculating scattering power spectrum of the tissues in the organ, and a tissue spatial average unit, connected to the power spectrum unit, for deriving a spatial average of the scattering power spectra of the tissues in the organ obtained by each scanning line. The blood power spectrum unit may calculate an average of the blood scattering power spectra for each duration of the blood flow speed's exceeding threshold levels or for each duration of the blood scattering power's exceeding a threshold level, and estimating blood scattering power spectrum at the portion of the tissues of the object of measurement by using the calculated average scattering power spectrum of the blood. The blood power spectrum unit may comprise a low pass filter having f -n characteristics, supplied with the received signals, a doppler detector, connected to the low pass filter, for calculating a blood flow speed, a blood flow power and scattering power spectrum of the blood, a duration unit, connected to the doppler detector, for detecting a duration of the blood flow speed's exceeding threshold levels or a duration of the blood flow power's exceeding a threshold level, a correction unit, connected to the doppler detector, for correcting a variable factor of the blood scattering power spectrum by an individual difference of volume percentage of blood corpuscle in each human, or for correcting the scattering power spectrum of the blood against the blood flow speed, an external designation unit, for externally designating a duration where the scattering power spectrum of the blood is calculated, an average unit, connected to the duration unit, the correction unit and the external designation unit, for averaging the scattering power spectra of the blood in a duration detected by the duration unit or designated by the external designation unit, an interpolation unit, connected to the average unit, for interpolating optional frequency components of each average scattering power spectrum of blood by an m-th order interpolation, an attenuation correction unit, connected to the interpolation unit, for correcting blood power spectrum in accordance with attenuation characteristics from the position of the blood to the position of the tissues, and a blood spatial average unit, connected to the attenuation correction unit and the scattering coefficient unit, for deriving a spatial average of the scattering power spectra of the blood obtained by each scanning line. The doppler detector may include a quadrature detector, an MTI-filter, a flow unit, and a doppler power spectrum unit. The MTI-filter may be applied to a series of data obtained by repeatedly transmitting ultrasound pulses to the same direction, and a complex Fourier transformation is carried out for data series of output from the MTI-filter for each transmitting, and the blood power spectrum is calculated by using results of the complex Fourier transformation. Each frequency component of the blood power spectrum and the tissue power spectrum may be modified so that the number of samples on the frequency axis become the same for each power spectrum. The scattering coefficient unit may comprise an ROI-setting unit, supplied with the received signals, for designating an ROI (region of interest) or a marker, and for activating a function of moving the ROI or the marker in accordance with a wall-movement, a normalizing unit, connected to the tissue power spectrum unit and the blood power spectrum unit, for calculating scattering coefficients b and n, and an angle correction unit, connected to the normalizing unit, for correcting the scattering coefficients b and n according to an angle between a cardiac muscle direction and an ultrasound beam. The ultrasound diagnostic equipment may further comprise an image unit, supplied with the received signals, for generating a B-mode image, a color flow image and/or an M-mode image, so that the display unit can display both the scattering coefficients b and n of the tissues and the B-mode image, the color flow image and/or the M-mode image. The scattering coefficient unit may obtain an oblique line with a linear approximation by applying a method of least squares to the normalized tissue scattering power spectrum, and the scattering coefficient n is determined by a slope of the oblique line and the scattering coefficient b is determined by a crossing point of a Y-axis and the oblique line, where an X-axis is established as logalithmic frequency and the Y-axis is established as decibel expression of the normalized power. The display unit may display the scattering coefficients b and n of the tissues as a time varying graph. The display unit may display the scattering coefficients b and n of the tissues accompanied with an electrocardiogram. The display unit may display a marker or a color at a duration of a blood flow speed's exceeding threshold levels or a duration of a blood flow power's exceeding a threshold level. The display unit may display an one-dimentional scattering coefficients b and n profile along a designated direction. The ultrasound diagnostic equipment may be used for diagnosing a myocardial infarction of a heart in the organ. According to an aspect of the present invention, there is also provided ultrasound diagnostic equipment for calculating and displaying scattering coefficients b and n of an organ from received signals, characterized in that the ultrasound diagnostic equipment comprises a probe, contacted to a body surface over the organ, for radiating an ultrasound beam to an optional portion using electrical pulses and receiving ultrasound waves scattered from the optional portion and converting the received ultrasound waves to received electrical signals, a transmitting circuit, for generating electrical pulses in accordance with a predetermined timing, a transmitting amplifier, connected between the probe and the transmitting circuit, for amplifying the electrical pulses and driving the probe, a receiving amplifier, connected to the probe, for amplifying electrical signals received by the probe, a tissue power spectrum unit, connected to the receiving amplifier, for calculating scattering power spectrum of a tissue region of an object of measurement in the organ, a blood power spectrum unit, connected to the receiving amplifier, for calculating scattering power spectrum of a blood region of the object of measurement in the organ, and for estimating blood scattering power spectrum at a portion of the tissue region by using the calculated scattering power spectrum of the blood, a scattering coefficient unit, connected to the tissue power spectrum unit and the blood power spectrum unit, for normalizing the tissue scattering power spectrum calculated by the tissue power spectrum unit by using the blood scattering power spectrum estimated by the blood power spectrum unit, and for calculating scattering coefficients b and n, an image unit, connected to the receiving amplifier, for generating a B-mode image and/or an M-mode image, and a display unit, connected to the scattering coefficient unit and the image unit, for displaying the scattering coefficients b and n calculated by the scattering coefficient unit and the B-mode image and/or M-mode image generated by the image unit. BRIEF DESCRIPTION OF THE DRAWINGS The present invention will be more clearly understood from the description of the preferred embodiments as set forth below with reference to the accompanying drawings, wherein: Figure 1 is a block diagram showing a construction of one embodiment according to a first aspect of the present invention; Figures 2A to 2C are diagrams for explaining a measuring method of IB; Figures 3A and 3B are explanatory diagrams of a blood scattering wave; Figure 4 is a block diagram showing an example of the first aspect of the present invention; Figures 5A to 5E are explanatory diagrams of a blood reference duration and a blood reference portion; Figures 6A to 6E are diagrams showing an example for displaying an IB; Figures 7A to 7D are diagrams showing another example for displaying an IB; Figure 8 is a block diagram showing a construction of one embodiment according to a second aspect of the present invention; Figures 9A and 9B are diagrams for explaining the concept of the second aspect of the present invention; Figures 10A and 10B are explanatory diagrams of estimating a scattering strength; Figure 11 is a block diagram showing an example of the second aspect of the present invention; Figures 12A to 12E are diagrams for explaining a comparison method between blood scattering power spectrum and tissue scattering power spectrum; Figures 13A to 13E are explanatory diagrams of a blood reference duration and a blood reference portion; Figures 14A to 14F are diagrams showing an example for displaying scattering coefficients; and Figures 15A to 15 D are diagrams showing another example for displaying scattering coefficients. DESCRIPTION OF THE PREFERRED EMBODIMENTS Below examples of ultrasound diagnostic equipment according to the present invention will be explained with reference to the drawings. First, an embodiment of ultrasound diagnostic equipment according to a first aspect will be explained with reference to Figs. 1 to 7D. Figure 1 is a block diagram showing a construction of one embodiment according to a first aspect. Ultrasound diagnostic equipment of the first aspect of the present invention is used for calculating and displaying an integrated backscatter IB of an organ from received ultrasound signals. The ultrasound diagnostic equipment is, for example, used for diagnosing a myocardial infarction of a heart in an organ. In Fig. 1, scattering power of tissues in an organ (tissue scattering power) at the tissues of an object of measurement (a tissue region B) is calculated by the tissue power unit 5 with using the received ultrasound signal reflected by the tissue region B. The scattering power of blood in an organ (blood scattering power) at a portion of blood of an object of measurement (a blood region A) is calculated by the blood power unit 6 in accordance with ultrasound signals reflected by the blood region A. The scattering power of blood at the tissue region B is estimated by the blood power unit 6 with using the calculated blood scattering power at the blood region A and the attenuation between the blood region A and the tissue region B. Note, the blood has definite ultrasound scattering characteristics, and attenuation characteristics. The tissue scattering power (scattering power of tissues), which is calculated by the tissue power unit 5, is normalized (divided) by the integrated backscatter unit 7 with using blood scattering power at the tissue region B estimated by the blood power unit 6, and an IB is determined by the normalized tissue scattering power. In Fig. 1, a probe 1, which is contacted to a body surface over an organ, is used for radiating an ultrasound beam 11 to a blood region A and a tissue region B, and for receiving ultrasound waves scattered from the blood region A and the tissue region B. A transmitting circuit 3 is used for generating electrical pulses in accordance with a predetermined timing. A transmitting amplifier 2, which is connected between the probe 1 and the transmitting circuit 3, is used for amplifying the electrical pulses and driving the probe 1. A receiving amplifier 4, which is connected to the probe 1, is used for amplifying electrical signals received by the probe 1. A tissue power unit 5, which is connected to the receiving amplifier 4, is used for calculating scattering power of the tissue region B of an object of measurement in the organ. A blood power unit 6, which is connected to the receiving amplifier 4, is used for calculating scattering power of a blood region A of the object of measurement in the organ, and for estimating blood scattering power at the tissue region B by using the calculated scattering power of the blood. An integrated backscatter unit 7, which is connected to the tissue power unit 5 and the blood power unit 6, is used for normalizing the tissue scattering power calculated by the tissue power unit 5 by using the blood scattering power estimated by the blood power unit 6, and for calculating an integrated backscatter IB. An image unit 8, which is connected to the receiving amplifier 4, is used for generating a B-mode image, a color flow image (CFI), an M-mode image, and the like. A display unit 9, which is connected to the integrated backscatter unit 7 and the image unit 8, is used for displaying the IB calculated by the integrated backscatter unit 7 and the B-mode image, the color flow image, the M-mode image, and the like, which are generated by the image unit 8. A tissue 10 is a tissue (tissues) of an object of measurement of the present embodiment. The tissue region B located in the tissue 10 is a region for calculating the IB. Note, the blood region A is a region of flowing blood whose scattering and attenuation characteristics are already known. Figures 2A to 2C are diagrams for explaining a measuring method of IB (Integrated Backscatter). As shown in Fig. 2A, an ultrasound beam 11 is radiated from the ultrasound probe 1, scattering power of the blood region A at a depth Z1 where blood exists is calculated, and blood scattering power at a depth Z3 is estimated by the calculated scattering power of the blood region A and the attenuation between the region A and the region B. Further, scattering power of the tissue region B at the depth Z3 is calculated. Note, as shown in Fig. 2A, a wall of dividing the blood and the tissues is located at a depth Z2. Next, a process of calculating an IB of the tissues at the depth Z3 by the scattering power of the tissue region B at the depth Z3 and the blood scattering power at the depth Z3 estimated by the scattering power of the blood at the depth Z1 will be explained with reference to equations. Assuming the frequency characteristics of transmitted ultrasound pulse signals to be I(f), the scattering power of tissues to be P T (Z), the scattering power of blood to be P B (Z), the characteristics including transmitting-receiving characteristics of an ultrasound probe and sound field characteristics to be F(f, Z), the round trip attenuation characteristics from the probe 1 to a depth Z to be A(f, Z), and the scattering characteristics of blood to be b B f 4 , the blood scattering power at a depth Z1 is indicated by the following equation (1). P B (Z1) = -∞ +∞ I(f) F(f, Z1) A(f, Z1) b B f 4 df The blood scattering power at the depth Z3 is approximately assumed to give the following equation from the scattering power of the equation (1). P B (Z3) = P B (Z1) -∞ +∞ G(f, Z3; Z1) A B (f, Z2; Z1) A T (f, Z3; Z2) df Where, G(f, Z3; Z1) indicates a correction term of a sound field change between the depth Z1 and Z3, A B (f, Z2; Z1) indicates attenuation characteristics of blood between the depth Z1 and Z2, and A T (f, Z3; Z2) indicates attenuation characteristics of the tissues between the depth Z2 and Z3. Therefore, a relationship indicated by the following equation is produced. A(f, Z3) = A(f, Z1) A B (f, Z2; Z1) A T (f, Z3; Z2) G(f, Z3; Z1) = F(f, Z3)/F(f, Z1) Further, the scattering power of the tissues at the depth Z3 is indicated by the following equation (3). P T (Z3) = -∞ +∞ I(f) F(f, Z3) A(f, Z3) b T f n df Next, normalizing the equation (3) by using the equation (2), the following equation (4) is obtained. P(Z3) = P T (Z3)/P B (Z3) Note, the equation (4) indicates an integrated backscatter IB of the tissues determined on the basis of the blood scattering power, in the present embodiment. Figure 2B is an explanatory graph chart of an IB corresponding to the equation (4) displayed by the display unit 9. The IB of the present embodiment is obtained by using scattering power of blood, which has a determined value with no individual difference as a standard, and thus the obtained IB has an absolute value. Therefore, the IB of the present embodiment gives two values; an absolute value and a rate of change, compared with the conventional IB's giving only rate of change. This means the IB of the present embodiment gives more information for diagnosis. Namely, the display unit 9 can display the IB of the tissue region B as a graph showing changes over time. Figure 2C is an explanatory electrocardiogram which can be incorporated to the graph of IB shown in Fig. 2B by the display unit 9, in order to enhance the timing information. Figures 3A and 3B are explanatory diagrams of a blood scattering wave. Figure 3A indicates scattering waves (amplitude of received signal) corresponding to a parenchyma (tissue) and a blood shown in Fig. 3B. As understood from the scattering waves, the scattering power of the blood is much smaller than the scattering power of the tissues. Further, large reflected waves may be frequently returned from the wall and/or the tissues and may overlap on the waves reflected from the blood. Therefore, it is difficult to calculate the scattering power of the blood with high accuracy due to the influence of reflections from the wall and/or the tissues. A blood scattering power can be calculated by using the doppler power of the blood flow obtained by a doppler detector 62 in Fig. 4, which will be described below, with decreasing the influence of reflections generated from the wall and/or the tissues. Figure 4 is a block diagram showing an example of the present invention. In Fig. 4, a received signal 41 indicates a signal received by a probe 1. The tissue power unit 5 is, for example, constituted by a power unit 51 and a tissue spatial average unit 52. The power unit 51, which is supplied with the received signals, is used for calculating scattering power of the tissue region B. The tissue spatial average unit 52, which is connected to the power unit 51, is used for deriving a spatial average of the scattering powers of the tissue region B obtained by each scanning line. The blood power unit 6 is, for example, constituted by a low pass filter 61, a doppler detector 62, a duration unit 63, a correction unit 64, an external designation unit 65, an average unit 66, an interpolation unit 67, an attenuation correction unit 68, and a blood spatial average unit 69. The low pass filter 61 has f -n characteristics, and is supplied with the received signals. Assuming that scattering strength (or a differential scattering cross section) is b B , it is already known that the scattering characteristics of blood are indicated by b B f 4 , and thus a signal having a constant scattering characteristics against frequency can be obtained by passing the signal through a filter having f -4 characteristics. Note, the n can be changed to another value. The doppler detector 62 is constituted by a quadrature detector 621, an MTI-filter 624, a flow unit 622, and a doppler power unit 623, and a blood flow speed and a blood flow power are calculated thereby. The duration unit 63 is connected to the doppler detector 62, or connected to the flow unit 622 and the doppler power unit 623. The duration unit 63 is used for detecting a duration of the blood flow speed's exceeding threshold levels or a duration of the blood flow power's exceeding a threshold level (with reference to Figs. 5C and 5D). The correction unit 64 is connected to the doppler detector 62, or connected to the doppler power unit 623. The correction unit 64 is used for correcting a variable factor of the blood scattering power by an individual difference of volume percentage of blood corpuscle in each human, or for correcting the scattering power of the blood against the blood flow speed. The external designation unit 65 is used for externally designating a duration where the blood scattering power is calculated. The average unit 66, which is connected to the duration unit 63, the correction unit 64 and the external designation unit 65, is used for averaging the scattering power of the blood in a duration detected by the duration unit 63 or designated by the external designation unit 65. The interpolation unit 67, which is connected to the average unit 66, is used for interpolating each average scattering power of the blood by an m-th order interpolation (with reference to Fig. 5E), when scattering power cannot be sufficiently obtained for the reason that the blood flow speed is too slow. Note, in Fig. 5E, the interpolation is a first-order interpolation. The attenuation correction unit 68, which is connected to the interpolation unit 67, is used for correcting blood power in accordance with attenuation from the blood region A to the tissue region B. The blood spatial average unit 69, which is connected to the attenuation correction unit 68 and the integrated backscatter unit 7 (an ROI-setting unit 71), is used for deriving a spatial average of the scattering powers of the blood obtained by each scanning line. The integrated backscatter unit 7 is, for example, constituted by an ROI-setting unit 71, a normalizing unit 72, and an angle correction unit 73. The ROI-setting unit 71, which is supplied with the received signals, is used for designating an ROI (region of interest) or a marker, and for activating a function of moving the ROI or the marker in accordance with a wall-movement. The normalizing unit 72, which is connected to the tissue power unit 5 (the tissue spatial average unit 52) and the blood power unit 6 (the blood spatial average unit 69), is used for calculating an integrated backscatter IB according to, for example, the equation (4). The angle correction unit 73, which is connected to the normalizing unit 72, is used for correcting the integrated backscatter IB according to an angle between a cardiac muscle direction and the ultrasound beam 11. More detailed explanation of this correction operation is given later with reference to Fig. 7. A reference numeral 81 denotes information of an image, an electrocardiogram, and the like from the ultrasound diagnostic equipment. A reference numeral 91 denotes a display control unit for controlling a display 92, which is connected to the integrated backscatter unit 7 (the normalizing unit 72) and is supplied with the information 81. Figures 5A to 5E are explanatory diagrams of a blood reference duration, and a blood reference portion and a tissue reference portion. Figure 5A indicates a B-mode image display area, calculable regions (B1, B2) of blood power, and detectable regions (T1, T2, T3) of tissues. Figure 5B indicates an electrocardiogram. Figure 5C indicates blood scattering power and durations of the blood power exceeding a specific threshold level with added markers. Figure 5D indicates-blood flow speed and durations of the blood flow speed (velocity) exceeding specific threshold levels with added markers. Namely, the display unit 9 can display a marker or a color at a duration of a blood flow speed's exceeding threshold levels or a duration of a blood flow power's exceeding a threshold level. Figure 5E indicates the scattering power interpolated by a first-order interpolation, which is shown by a dot-line in Fig. 5E. The integrated backscatter IB of the tissues for all durations can be obtained by using the above interpolation. Figures 6A to 6E are diagrams showing an example of displaying an IB. Figure 6A indicates calculation portions of an IB of tissues T1 and T2, and blood reference regions B1 and B2 by using a display area of a B-mode image. Figure 6B indicates an electrocardiogram. Figure 6C indicates examples of time varying graph of the IB of the tissues T1 and T2. Figure 6D indicates a marker in a blood reference time phase. Figure 6E indicates an example of a color coded image of an IB superimposed on an M-mode image. Note, Figs. 6A to 6E can be voluntarily indicated on the display unit in any combination. Figures 7A to 7D are diagrams showing an example of displaying an IB. Figure 7A indicates an IB, which is indicated by a color indication, on a B-mode image with a well-known color flow image (CFI). Furthermore, an ROI for indicating a region, where the blood power is calculated, is also indicated. In this diagram, the color flow image and the IB are simultaniously indicated, but it can be constituted so that the color flow image and the IB are separately indicated by different displays, or only the IB is indicated. Figure 7B indicates a marker which is externally established at the state of freezing of the B-mode image or the color flow image to indicate a cardiac muscle direction. Note, it is well known that the integrated backscatter IB varies depending on an angle between a scanning direction and the musle. Figure 7C indicates an example of the factor to correct the angle dependency. Figure 7D indicates an explanatory example of one-dimensional profile of an integrated backscatter IB of the tissues with reference to blood. Namely, the display unit 9 can display a spatial IB of the tissues as a profile. According to the first aspect as in the above descriptions, blood scattering power at a portion of tissues of an object of measurement is estimated by using a scattering power of blood which has the definite ultrasound scattering characteristics, scattering power of the tissues of the object of measurement is normalized (divided) by the estimated scattering power of the blood, and an integrated backscatter IB is obtained and displayed, so that the IB of the tissues can be quantitatively estimated without being influenced by the attenuation. Furthermore, this estimated quantitative IB of the tissues is displayed on the display unit accompanied with a B-mode image, an electrocardiogram, and the like, so that diagnostic accuracy can be increased. Next, an embodiment of ultrasound diagnostic equipment according to an aspect of the present invention will be explained with reference to Figs. 8 to 15D. Figure 8 is a block diagram showing a construction of one embodiment according to a second aspect of the present invention. Ultrasound diagnostic equipment of the second aspect of the present invention is used for calculating and displaying scattering coefficients b and n of an organ from received ultrasound signals. The ultrasound diagnostic equipment is, for example, used for diagnosing a myocardial infarction of a heart in an organ. Note, the scattering coefficients b and n are defined in the following equation. S ( f ) = b f n In Fig. 8, a tissue power spectrum unit 5a is used for calculating scattering power spectrum of tissues of an object of measurement in an organ. A blood power spectrum unit 6a is used for calculating scattering power spectrum of blood in the organ and estimating blood scattering power spectrum at a portion of the tissues of the object of measurement by using the calculated scattering power spectrum of the blood and the attenuation characteristics between a blood region A and a tissue region B. A scattering coefficient unit 7a is used for normalizing the tissue scattering power by using the estimated blood scattering power spectrum, and for calculating scattering coefficients b and n by the normalized tissue scattering power spectrum. As shown in Fig. 8, scattering power spectrum of tissues in an organ (tissue scattering power spectrum) at the tissues of an object of measurement (a tissue region B) is calculated by the tissue power spectrum unit 5a by using the received ultrasound signal reflected by the tissue region B. The scattering power spectrum of blood in an organ (blood scattering power spectrum) at a portion of blood of an object of measurement (a blood region A) is calculated by the blood power spectrum unit 6a by using ultrasound signals reflected by the blood region A. The scattering power of blood at the tissue region B is estimated by the blood power spectrum unit 6a by using the calculated blood scattering power spectrum at the blood region A and the attenuation characteristics between the blood region A and the tissue region B. Note, the blood has definite ultrasound scattering and attenuation characteristics. The tissue scattering power spectrum (scattering power spectrum of tissues), which is calculated by the tissue power spectrum unit 5a, is normalized (divided) by the scattering coefficient unit 7a with using blood scattering power spectrum at the tissue region B estimated by the blood power spectrum unit 6a, and scattering coefficients b and n are determined by the normalized tissue scattering power spectrum. In Fig. 8, the same constituents as in Fig. 1 will have the same name and reference numeral as in Fig. 1, and the explanation thereof will be omitted. A tissue power spectrum unit 5a, which is connected to the receiving amplifier 4, is used for calculating scattering power spectrum of the tissue region B of an object of measurement in the organ. A blood power spectrum unit 6a, which is connected to the receiving amplifier 4, is used for calculating scattering power spectrum of the blood region A of the object of measurement in the organ, and for estimating blood scattering power spectrum at the tissue region B by using the calculated scattering power spectrum of the blood. A scattering coefficient unit 7a, which is connected to the tissue power spectrum unit 5a and the blood power spectrum unit 6a, is used for normalizing the tissue scattering power spectrum calculated by the tissue power spectrum unit 5a by using the blood scattering power spectrum estimated by the blood power spectrum unit 6a, and for calculating scattering coefficients b and n. Note, the scattering coefficient n is, for example, determined by a slope of the normalized spectrum, and the the scattering coefficient b is, for example, determined by a crossing point of a Y-axis and the normalized spectrum (with reference to Fig. 9B). Note, an X-axis is established as logalithmic frequency and the Y-axis is established as decibel expression of the normalized power as shown in Fig. 9B. A display unit 9, which is connected to the scattering coefficient unit 7a and the image unit 8, is used for displaying the scattering coefficients b and n calculated by the scattering coefficient unit 7a and the B-mode image, the color flow image (CFI), the M-mode image, and the like, which are generated by the image unit 8. A tissue 10 is a tissue (tissues) of an object of measurement of the present embodiment. The tissue region B located in the tissue 10 is a region for calculating the scattering coefficients. Note, the blood region A is a region of flowing blood whose scattering and attenuation characteristics are already known. Figures 9A and 9B are diagrams for explaining the concept of the aspect of the present invention. As shown in Fig. 9A, an ultrasound beam 11 is radiated from the probe 1, the scattering power spectrum of the blood region A at a depth Z1 where blood exists is calculated, and the blood scattering power spectrum at a depth Z3 is estimated by the calculated scattering power spectrum of the blood region A and the attenuation characteristics between the blood region A and the tissue region B. Further, the scattering power spectrum of the tissue region B at the depth Z3 is calculated. Note, as shown in Fig. 9A, a wall dividing the blood and the tissues is located at a depth Z2. Next, a process of calculating scattering coefficients b and n of the tissues at the depth Z3 by the scattering power spectrum of the tissue region B at the depth Z3 and the blood scattering power spectrum at the depth Z3 estimated by the scattering power spectrum of the blood at the depth Z1, will be explained with reference to equations. Assuming the frequency characteristics of transmitted ultrasound pulse signals to be I(f), the scattering power spectrum of tissues to be S T (f, Z), the scattering power spectrum of blood to be S B (f, Z), the characteristics including transmitting-receiving characteristics of a probe and sound field characteristics to be F(f, Z), the round trip attenuation characteristics from the probe 1 to a depth Z to be A(f, Z), and the scattering characteristics of blood to be b B f 4 , the blood scattering power spectrum at a depth Z1 is indicated by the following equation (5). S B (f, Z1) = I(f) F(f, Z1) A(f, Z1) b B f 4 The scattering power spectrum at the depth Z3 is assumed to give the following equation (6) from the scattering power spectrum of the equation (5). S B (f, Z3) = G(f, Z3; Z1) A B (f, Z2; Z1) A T (f, Z3; Z2) S B (f, Z1) Where, G(f, Z3; Z1) indicates a correction term of a sound field change between the depth Z1 and Z3, A B (f, Z2; Z1) indicates attenuation characteristics of blood between the depth Z1 and Z2, and A T (f, Z3; Z2) indicates attenuation characteristics of the tissues between the depth Z2 and Z3. Note, A T (f, Z3; Z2) is unknown, so we assume it to specific characteristics. Further, the scattering power spectrum of the tissues at the depth Z3 is indicated by the following equation (7). S T (f, Z3) = I(f) F(f, Z3) A(f, Z3) b T f n Where, A(f, Z3) = A(f, Z1) A B (f, Z2; Z1) A T (f, Z3; Z2) G(f, Z3; Z1) = F(f, Z3)/F(f, Z1) Therefore, normalizing the equation (7) by using the equation (6), the following equation (8) is obtained. S(f, Z3) = S T (f, Z3)/(S B (f, Z3) · f -4 ) = b T /b B · f n The equation (8) is indicated by the following decibel expression equation (9). 10 log 10 S(f, Z3) = 10 log 10 (b T /b B ) + 10 log 10 f · n A graph indicating the equation (9) is shown in Fig. 9B. Note, an oblique line is obtained with a linear approximation by applying a method of least squares to the normalized scattering power spectrum 10 log 10 S(f, Z3), and a scattering coefficient n is, for example, determined by a slope of the oblique line and a scattering strength b (= b T /b B ) is determined by a crossing point of a Y-axis and the oblique line as shown in Fig. 9B. Namely, the scattering coefficient unit 7a can obtain an oblique line with a linear approximation by applying a method of least squares to the normalized tissue scattering power spectrum, and the scattering coefficient n is determined by a slope of the oblique line and the scattering coefficient b is determined by a crossing point of a Y-axis and the oblique line. Note, an X-axis is established as logalithmic frequency and the Y-axis is established as decibel expression of the normalized power as descrived in above with reference to Fig. 9B. Therefore, according to the present embodiment, scattering power spectrum of a blood region A at a depth Z1 is calculated, the blood scattering power spectrum at a depth Z3 is estimated, the tissue scattering power spectrum of a tissue region B at the depth Z3 is calculated, the tissue scattering power spectrum of the tissue region B is normalized using the estimated blood scattering power spectrum, a graph of the equation (9) is made, and then a scattering coefficient n is determined by a slope of the oblique line and a scattering coefficient b is determined by a crossing point of a Y-axis. Figures 10A and 10B are explanatory diagrams of estimating scattering coefficients (scattering strength) from the normalized scattering power spectrum 10 log 10 S(f, Z3). Figure 10A indicates an approach, wherein an oblique line is obtained with a linear approximation by applying a method of least squares to a normalized scattering power spectrum between an effective band (f 1 - f 2 ), and a scattering strength b at a center frequency f c of the effective band (which is a central portion between f 1 and f 2 on a logalithmic frequency axis). According to this approach, an estimated value b can be obtained with only a small fluctuation. Figure 10B indicates another approach, wherein an average power between f 1 and f c and an average power between f c and f 2 are calculated, and a scattering strength b is determined by the following equation (10). Scattering strength b at f c = [(average power between f 1 - f c ) + (average power between f c - f 2 )] / 2 Figures 3A and 3B, which are used for explaining an embodiment of the first aspect of the present invention, are an explanatory diagrams of a blood scattering wave. Figure 3A indicates scattering waves (amplitude of received signal) corresponding to a parenchyma (tissue) and a blood shown in Fig. 3B. As understood from the scattering waves, the scattering power of the blood is much smaller than the scattering power of the tissues. Further, large reflected waves may be frequently returned from the wall and/or the tissues and may overlap on the waves reflected from the blood. Therefore, it is difficult to calculate the scattering power of the blood with high accuracy due to the influence of reflections from the wall and/or the tissues. These multiple reflections can be removed from the scattering power spectrum of the blood by a doppler detector 62a in Fig. 11, which will be described below. Figure 11 is a block diagram showing an example of the second aspect of the present invention. In Fig. 11, a received signal 41 indicates a signal received by a probe 1. The tissue power spectrum unit 5a is, for example, constituted by a power spectrum unit 51a and a tissue spatial average unit 52a. The power spectrum unit 51a, which is supplied with the received signals, is used for calculating scattering power spectrum of the tissue region B. The tissue spatial average unit 52a, which is connected to the power spectrum unit 51a, is used for deriving a spatial average of the scattering power spectrum of the tissue region B obtained by each scanning line. The blood power spectrum unit 6a is, for example, constituted by a low pass filter 61, a doppler detector 62a, a duration unit 63, a correction unit 64a, an external designation unit 65, an average unit 66a, an interpolation unit 67a, an attenuation correction unit 68a, and a blood spatial average unit 69a. The low pass filter 61 is the same as in Fig. 4. The doppler detector 62a is constituted by a quadrature detector 621, an MTI-filter 624, a flow unit 622, and a doppler power spectrum unit 623a and the like, and a speed, scattering power spectrum and scattering power of the blood are calculated thereby. The duration unit 63 is connected to the doppler detector 62a, or connected to the flow unit 622 and the doppler power spectrum unit 623a. The duration unit 63 is used for detecting a duration of the blood flow speed's exceeding threshold levels or a duration of the blood flow power's exceeding a threshold level (with reference to Figs. 13C and 13D). The correction unit 64a is connected to the doppler detector 62a, or connected to the doppler power spectrum unit 623a. The correction unit 64a is used for correcting a variable factor of the blood scattering power spectrum by an individual difference of volume percentage of blood corpuscle in each human, or for correcting the scattering power spectrum of the blood against the blood flow speed. The external designation unit 65 is used for externally designating a duration where the blood scattering power spectrum is calculated. The average unit 66a, which is connected to the duration unit 63, the correction unit 64a and the external designation unit 65, is used for averaging the scattering power spectrum of the blood in a duration detected by the duration unit 63 or designated by the external designation unit 65. The interpolation unit 67a, which is connected to the average unit 66a, is used for interpolating scattering power spectrum of blood by an m-th order interpolation (with reference to Fig. 13E), when scattering power spectrum cannot be sufficiently obtained for the reason that the blood flow speed is too slow. Note, in Fig. 13E, the interpolation is a first-order interpolation. The attenuation correction unit 68a, which is connected to the interpolation unit 67a, is used for correcting blood power spectrum in accordance with attenuation characteristics between the blood region A and the tissue region B. The blood spatial average unit 69a, which is connected to the attenuation correction unit 68a and the scattering coefficient unit 7a (an ROI-setting unit 71), is used for deriving a spatial average of the scattering power spectra of the blood obtained by each scanning line. The scattering coefficient unit 7a is, for example, constituted by an ROI-setting unit 71, a normalizing unit 72a, and an angle correction unit 73. The ROI-setting unit 71, which is supplied with the received signals, is used for designating an ROI or a marker, and for activating a function of moving the ROI or the marker in accordance with a wall-movement. The normalizing unit 72a, which is connected to the tissue power spectrum unit 5a (the tissue spatial average unit 52a) and the blood power spectrum unit 6a (the blood spatial average unit 69a), is used for calculating scattering coefficients b and n. The angle correction unit 73, which is connected to the normalizing unit 72a, is used for correcting the scattering coefficients b and n according to an angle between a cardiac muscle direction and the ultrasound beam 11. A reference numeral 81 denotes information of an image, an electrocardiogram, and the like from the ultrasound diagnostic equipment. A reference numeral 91 denotes a display control unit for controlling a display 92, which is connected to the scattering coefficient unit 7a (the normalizing unit 72a) and is supplied with the information 81. Figures 12A to 12E are diagrams for explaining a method to compare a blood scattering power spectrum with tissue scattering power spectrum. Figure 12A indicates real components and imaginary components after a quadrature detection. Namely, Fig. 12A indicates two output waveforms as examples of output waveforms in a quadrature detector 621. When estimating a blood flow speed and the like by using the quadrature detector 621, K ultrasound pulses are transmitted in the same scanning line direction. A reference I denotes an imaginary component, and a reference R denotes a real component. Below, a calculation of the scattering power spectrum of blood flow using the quadrature detector will be explained. First, in order to remove influence from a slow wall-signal or a reflected signal of a parencyma, a well-known MTI-filter is applied to a series of data sampled at the same depth. Namely, the MTI-filter is applied to a series of data (R 1j , R 2j , ··· ··· R kj ), (I 1j , I 2j , ··· ··· (I kj ); (j = 1, 2, ··· ··· M), which are obtained by repeatedly transmitting ultrasound pulses and sampling the output of the quadrature detector at the timing corresponding to the same depth, and then the power spectrum is calculated by carrying out a complex Fourier transformation on the data series R ij , I ij , (where, j = 1 to M). Figure 12B indicates an average power spectrum which is related to i as described above. Note, the power spectrum of the tissues can be obtained by carrying out a Fourier transformation to the reflected signal. Figure 12C indicates estimated power spectrum of the tissues for data acquisition duration T with a sampling frequency f s . Figure 12D indicates blood power spectrum obtained by shifting the power spectrum of the blood as shown in Fig. 12B by a value corresponding to a mixing frequency. Figure 12E indicates blood power spectrum obtained by averaging power spectrum of the tisseus of Fig. 12C by N sx pieces on a frequency axis, where, N sx = (f k /M)/(f s /T) Note, scattering coefficients can be obtained by comparing Figs. 12D and 12E. Figures 13A to 13E are explanatory diagrams of a blood reference duration, a blood reference portion and a tissue reference portion. Figure 13A indicates a B-mode image display area, calculable regions (B1, B2) of blood power spectrum, and detectable regions (T1, T2, T3) of tissues. Figure 13B indicates an electrocardiogram. Figure 13C indicates blood scattering power and durations of the blood power exceeding a specific threshold level with added markers. Figure 13D indicates blood flow speed and durations of the blood flow speed (velocity) exceeding specific threshold levels with added markers. Namely, the display unit 9 can display a marker or a color at a duration of a blood flow speed's exceeding threshold levels or a duration of a blood flow power's exceeding a threshold level. Figure 13E indicates a specific frequency component of the scattering power spectrum interpolated by a first-order interpolation, which is shown by a dot-line in Fig. 13E. The scattering coefficients of the tissues in all durations can be obtained by using the above interpolated power spectra. Figures 14A to 14F are diagrams showing an example for displaying scattering coefficients. Figure 14A indicates calculation portions of scattering coefficients of tissues T1 and T2, and blood reference regions B1 and B2 by using a display area of a B-mode image. Figure 14B indicates an electrocardiogram. Figure 14C indicates an example of displaying a changing state with time of the scattering coefficient n of the tissues T1 and T2. Figure 14D indicates an example of displaying a changing state with time of the scattering coefficient b of the tissues T1 and T2. Figure 14E indicates a marker in a blood reference time phase. Figure 14F indicates an example of a color coded image of scattering coefficients on an M-mode image. Note, Figs. 14A to 14F can be voluntarily indicated on the display unit in any combination. Figures 15A to 15D are diagrams showing an example of displaying scattering coefficients. Figure 15A indicates scattering coefficients b and n, which are indicated by a color indication, on a B-mode image with a color flow image (CFI). Furthermore, an ROI for indicating a region of reference blood is also indicated. In this diagram, the color flow image and the scattering coefficients b and n are simultaniously indicated, but it can be constituted that the color flow image and the scattering coefficients b and n are separately indicated by different displays, or the scattering coefficients b and n only are indicated. Figure 15B indicates a marker which is externally established at the state of freezing of the B-mode image or the CFI to indicate a cardiac muscle direction. Note, it is well known that the scattering coefficients varY depending on an angle between a scanning direction and the musle. Figure 15C indicates an example of the factor to correct the angle dependency. Figure 15D indicates an explanatory example of one-dimensional profile of profiles scattering coefficients b and n of the tissues with reference to blood. Namely, the display unit 9 can display spatial scattering coefficients b and n of the tissues as a profile. According to an aspect of the present invention as in the above descriptions, blood scattering power spectrum at a portion of tissues of an object of measurement is estimated by using scattering power spectrum of blood which has a definite ultrasound scattering characteristics, the scattering power spectrum of the tissues of the object of measurement is normalized (divided) by the estimated scattering power spectrum, and scattering coefficients b and n are obtained and displayed, so that the scattering coefficients b and n of the tissues can be exactly estimated without being influenced by the attenuation. Furthermore, these estimated scattering coefficients b and n of the tissues are displayed on the display unit accompanied with a B-mode image, an electrocardiogram, and the like, so that diagnostic accuracy can be increased.

Claims: Ultrasound diagnostic equipment for calculating and displaying scattering coefficients b and n of an organ from received signals when defining scattering frequency characteristics S(f) = b·f n , said ultrasound diagnostic equipment comprising: a tissue power spectrum means (5a), supplied with said received signals, for calculating scattering power spectrum of tissues of an object of measurement in said organ; a blood power spectrum means (6a), supplied with said received signals, for calculating scattering power spectrum of blood in said organ and estimating blood scattering power spectrum at the portion of the tissues of said object of measurement by using said calculated scattering power spectrum of the blood; a scattering coefficient means (7a), connected to said tissue power spectrum means (5a) and said blood power spectrum means (6a), for normalizing tissue scattering power spectrum by using said estimated blood scattering power spectrum, and for calculating scattering coefficients b and n of the tissues of the object of measurement in said organ by said normalized tissue scattering power spectrum; and a display means (9), connected to said scattering coefficient means (7a), for displaying said scattering coefficients b and n. Ultrasound diagnostic equipment as claimed in claim 1, wherein a normalizing unit (72a) is connected for calculating scattering coefficients b and n; and an angle correction unit (73) is connected to said normalizing unit (72a), for correcting said scattering coefficients b and n according to an angle between a cardiac muscle direction and an ultrasound beam (11; Fig. 15C). Ultrasound diagnostic equipment as claimed claim 1 or 2, wherein said tissue power means (5) or said tissue power spectrum means (5a) comprises: a power unit (51), supplied with said received signals, for calculating scattering power of the tissue in said organ, or a power spectrum unit (51a), supplied with said received signals, for calculating scattering power spectrum of the tissue in said organ; and a tissue spatial average unit (52), connected to said power unit (51), for deriving a spatial average of the scattering powers of the tissues in said organ obtained by each scanning line. Ultrasound diagnostic equipment as claimed in any of claims 1 to 3, wherein said blood power means (6) calculates an average of said blood scattering powers for each duration of the blood flow speed's exceeding threshold levels or for each duration of the blood scattering power's exceeding a threshold level, and estimating blood scattering power at the portion of the tissues of said object of measurement by using said calculated average scattering power of the blood. Ultrasound diagnostic equipment as claimed in any of claims 1 to 4, wherein said blood power means (6) comprises (Fig. 4): a low pass filter (61) having f -n characteristics, supplied with said received signals; a doppler detector (62), connected to said low pass filter (61), for calculating a blood flow speed and a blood flow power; a duration unit (63), connected to said doppler detector (62), for detecting a duration of the blood flow speed's exceeding threshold levels or a duration of the blood flow power's exceeding a threshold level (Figs. 5C, 5D); a correction unit (64), connected to said doppler detector (62), for correcting a variable factor of the blood scattering power by an individual difference of volume percentage of blood corpuscle in each human, or for correcting the scattering power of the blood against the blood flow speed; an external designation unit (65), for externally designating a duration where the blood scattering power is calculated; an average unit (66), connected to said duration unit (63), said correction unit (64) and said external designation unit (65), for averaging the scattering powers of said blood in a duration detected by said duration unit (63) or designated by said external designation unit (65); an interpolation unit (67), connected to said average unit (66), for interpolating each average scattering power of said blood by an m-th order interpolation (Fig. 5E); an attenuation correction unit (68), connected to said interpolation unit (67), for correcting blood power in accordance with attenuation from the position of said blood to the position of said tissues; and a blood spatial average unit (69), connected to said attenuation correction unit (68) and said integrated backscatter means (7), for deriving a spatial average of the scattering powers of said blood obtained by each scanning line. Ultrasound diagnostic equipment as claimed in claim 5, wherein said doppler detector (62) includes a quadrature detector (621), an MTI-filter (624), a flow unit (622), and a doppler power unit (623). Ultrasound diagnostic equipment as claimed in any of claims 1 to 6, wherein said integrated backscatter means (7) comprises: an ROI-setting unit (71), supplied with said received signals, for designating an ROI or a marker, and for activating a function of moving the ROI or the marker in accordance with a wall-movement; a normalizing unit (72), connected to said tissue power means (5) and said blood power means (6), for calculating an integrated backscatter (IB); and an angle correction unit (73), connected to said normalizing unit (72), for correcting said integrated backscatter (IB) according to an angle between a cardiac muscle direction and an ultrasound beam (11; Fig. 7C). Ultrasound diagnostic equipment as claimed in any of claims 1 to 7, wherein said ultrasound diagnostic equipment further comprises an image means (8), supplied with said received signals, for generating a B-mode image, a color flow image (CFI) and/or an M-mode image, so that said display means (9) displays both said integrated backscatter (IB) of the tissues and said B-mode image, said color flow image (CFI) and/or said M-mode image (Figs. 1, 6E, 7A). Ultrasound diagnostic equipment according to any of claims 1 to 8, wherein said display means (9) displays said integrated backscatter (IB) of the tissues as a time varying graph (Fig. 2B). Ultrasound diagnostic equipment according to any of claims 1 to 9, wherein said display means (9) displays said integrated backscatter (IB) of the tissues accompanied with an electrocardiogram (Figs. 2B, 2C). Ultrasound diagnostic equipment according to any of claims 1 to 9, wherein said display means (9) displays a marker or a color at a duration of a blood flow speed's exceeding threshold levels or a duration of a blood flow power's exceeding a threshold level (Fig. 6D). Ultrasound diagnostic equipment according to any of claims 1 to 9, wherein said display means (9) displays an one-dimensional integrated backscatter (IB) profile along a designated direction (Figs. 7B, 7D). Ultrasound diagnostic equipment according to any of claims 1 to 12, wherein said ultrasound diagnostic equipment is used for diagnosing a myocardial infarction of a heart in said organ. Ultrasound diagnostic equipment as claimed in any of claims 3 to 13, wherein said MTI-filter (624) is applied to a series of data obtained by repeatedly transmitting ultrasound pulses to the same direction, and a complex Fourier transformation is carried out for data series of output from said MTI-filter (624) for each transmitting, and said blood power spectrum is calculated by using results of said complex Fourier transformation. Ultrasound diagnostic equipment as claimed in any of claims 1 to 14, wherein each frequency component of said blood power spectrum and said tissue power spectrum are modified so that the number of samples on the frequency axis become the same for each power spectrum. Ultrasound diagnostic equipment as claimed in any of claims 1 to 15, wherein said scattering coefficient means (7a) obtains an oblique line with a linear approximation by applying a method of least squares to the normalized tissue scattering power spectrum, and said scattering coefficient n is determined by a slope of said oblique line and said scattering coefficient b is determined by a crossing point of a Y-axis and said oblique line, where an X-axis is established as logarithmic frequency and said Y-axis is established as decibel expression of the normalized power (Fig. 9B). Ultrasound diagnostic equipment as claimed in any of claims 1 to 16, wherein said tissue power spectrum unit (5a) comprises: a power spectrum unit (51a), supplied with said received signals, for calculating scattering power spectrum of said tissue region (B).

Inventors: MURAKAMI KEIICHI, SHIBA AKIRA, SHIMURA TAKAKI, YAMADA ISAMU, MURAKAMI, KEIICHI, SHIBA, AKIRA, SHIMURA, TAKAKI, YAMADA, ISAMU, SHIBA, AKIRA, MONAKU MANSION 307
Applicants: FUJITSU LTD, FUJITSU LIMITED
Application year: 1990
Main classification: A61B8
IPC codes: A61B8, G01S7