Subject information acquisition apparatus

0Associated
Cases 
0Associated
Defendants 
0Accused
Products 
0Forward
Citations 
0
Petitions 
0
Assignments
First Claim
1. A subject information acquisition apparatus for transmitting, through a reference layer, an elastic wave to a subject, receiving, through the reference layer, the elastic wave reflected from inside the subject, and acquiring information inside the subject, the subject information acquisition apparatus comprising:
 a conversion device configured to receive the elastic wave and convert the elastic wave into an electrical signal; and
a memory and processor which function as each unit comprising;
a measuring unit configured to measure strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer interposed between the conversion device and the subject by using the electrical signal; and
a calculation unit configured to calculate a viscosity coefficient (TO) of the subject by using a viscosity coefficient and an elastic modulus of the reference layer, the strain of the subject, and the strain of the reference layer,whereinthe pressure applied to the subject and the reference layer is periodic pressure application including releasing the pressure, andthe calculation unit calculates the viscosity coefficient (η
) of the subject by using an area of a hysteresis loop drawn with strain values of the subject and strain values of the reference layer which are generated by the periodic pressure.
0 Assignments
0 Petitions
Accused Products
Abstract
Viscoelastic characteristics in a subject are imaged by a simple method. A viscoelasticity measurement reference layer whose elastic modulus and viscosity coefficient are known is included between an ultrasonic wave transmitting/receiving probe and the subject and distributions of elastic modulus and viscosity coefficient inside the subject are calculated from a change over time of strain generated in the viscoelasticity measurement reference layer and the subject according to a pressure applied to the subject which changes over time and known elastic modulus and viscosity coefficient of the viscoelasticity measurement reference layer.
5 Citations
No References
REFERENCE DEFORMABLE BODY, ULTRASONIC DIAGNOSTIC APPARATUS, AND ULTRASONIC DIAGNOSTIC METHOD  
Patent #
US 20110040185A1
Filed 04/14/2009

Current Assignee
Hitachi America Limited

Sponsoring Entity
Hitachi Medical Corporation

METHOD AND APPARATUS FOR DETERMINING VISCOELASTIC PARAMETERS IN TISSUE  
Patent #
US 20100160778A1
Filed 11/03/2009

Current Assignee
University of British Columbia

Sponsoring Entity
University of British Columbia

METHOD AND SYSTEM FOR AUTOMATIC DETECTION OF OBJECTS IN AN IMAGE  
Patent #
US 20080267499A1
Filed 04/30/2007

Current Assignee
General Electric Company

Sponsoring Entity
General Electric Company

SUBJECT INFORMATION ACQUISITION APPARATUS  
Patent #
US 20140221833A1
Filed 04/05/2012

Current Assignee
Canon Kabushiki Kaisha

Sponsoring Entity
Canon Kabushiki Kaisha

Ultrasound diagnostic apparatus  
Patent #
US 9,289,186 B2
Filed 03/18/2013

Current Assignee
Fujifilm Corporation

Sponsoring Entity
Fujifilm Corporation

24 Claims
 1. A subject information acquisition apparatus for transmitting, through a reference layer, an elastic wave to a subject, receiving, through the reference layer, the elastic wave reflected from inside the subject, and acquiring information inside the subject, the subject information acquisition apparatus comprising:
a conversion device configured to receive the elastic wave and convert the elastic wave into an electrical signal; and a memory and processor which function as each unit comprising; a measuring unit configured to measure strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer interposed between the conversion device and the subject by using the electrical signal; and a calculation unit configured to calculate a viscosity coefficient (TO) of the subject by using a viscosity coefficient and an elastic modulus of the reference layer, the strain of the subject, and the strain of the reference layer, wherein the pressure applied to the subject and the reference layer is periodic pressure application including releasing the pressure, and the calculation unit calculates the viscosity coefficient (η
) of the subject by using an area of a hysteresis loop drawn with strain values of the subject and strain values of the reference layer which are generated by the periodic pressure. View Dependent Claims (2, 3, 4, 5)
 6. A subject information acquisition apparatus for transmitting, through a reference layer, an elastic wave to a subject, receiving, through the reference layer, the elastic wave reflected from inside the subject, and acquiring information inside the subject, the subject information acquisition apparatus comprising:
a conversion device configured to receive the elastic wave and convert the elastic wave into an electrical signal; and a memory and processor which function as each unit comprising; a measuring unit configured to measure strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer interposed between the conversion device and the subject by using the electrical signal, the reference layer having a viscosity coefficient and an elastic modulus; and a calculation unit configured to calculate a viscosity coefficient (TO) of the subject by using the viscosity coefficient and the elastic modulus of the reference layer, the strain of the subject, and the strain of the reference layer, wherein the calculation unit calculates the viscosity coefficient (η
) of the subject by using a delay time between the change over time of the strain value of the subject and the change over time of the strain value of the reference layer which are generated by the application of the pressure. View Dependent Claims (7, 8)
 9. A subject information acquisition apparatus for transmitting, through a reference layer, an elastic wave to a subject, receiving, through the reference layer, the elastic wave reflected from inside the subject, and acquiring information inside the subject, the subject information acquisition apparatus comprising:
a conversion device configured to receive the elastic wave and convert the elastic wave into an electrical signal; and a memory and processor which function as each unit comprising; a measuring unit configured to measure strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer interposed between the conversion device and the subject by using the electrical signal; and a calculation unit configured to calculate a stress applied to the subject by using a viscosity coefficient of the reference layer, a Young'"'"'s modulus of the reference layer, and the strain of the reference layer, according to the following equation;  View Dependent Claims (10, 11, 12)
 13. A subject information acquisition method for acquiring information on inside of a subject by using echo signals from the subject, the method comprising:
transmitting, through a reference layer, an elastic wave to the subject; receiving, through the reference layer, the elastic wave reflected in the subject; measuring strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer; and calculating a viscosity coefficient (η
) of the subject by using a viscosity coefficient and an elastic modulus of the reference layer, the strain of the subject, and the strain of the reference layer,wherein the pressure applied to the subject and the reference layer is periodic pressure application including releasing the pressure, and calculating the viscosity coefficient (η
) of the subject is performed by using an area of a hysteresis loop drawn with strain values of the subject and strain values of the reference layer which are generated by the periodic pressure. View Dependent Claims (14, 15, 16, 17)
 18. A subject information acquisition method for acquiring information on inside of a subject by using echo signals from the subject, the method comprising:
transmitting, through a reference layer, an elastic wave to the subject; receiving, through the reference layer, the elastic wave reflected in the subject; measuring strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer; and calculating a viscosity coefficient (η
) of the subject by using a viscosity coefficient and an elastic modulus of the reference layer, the strain of the subject, and the strain of the reference layer,wherein calculating the viscosity coefficient (η
) of the subject is performed by using a delay time between the change over time of the strain value of the subject and the change over time of the strain value of the reference layer which are generated by the application of the pressure. View Dependent Claims (19, 20)
 21. A subject information acquisition method for acquiring information on inside of a subject by using echo signals from the subject, the method comprising:
transmitting, through a reference layer, an elastic wave to the subject; receiving, through the reference layer, the elastic wave reflected in the subject; measuring strain of the subject and strain of the reference layer which are generated when a pressure is applied to the subject and the reference layer; and calculating a stress applied to the subject by using a viscosity coefficient of the reference layer, a Young'"'"'s modulus of the reference layer, and the strain of the reference layer, according to the following equation;  View Dependent Claims (22, 23, 24)
1 Specification
This application is a Continuation of copending U.S. patent application Ser. No. 14/110,125 filed Mar. 6, 2014, which is a U.S. national stage application of International Patent Application No. PCT/JP2012/059952, filed Apr. 5, 2012, which claims foreign priority benefit of Japanese Patent Application No. 2011086512, filed Apr. 8, 2011. All of the abovenamed patent applications are hereby incorporated by reference herein in their entirety.
The present invention relates to a subject information acquisition apparatus, in particular to a subject information acquisition apparatus for imaging viscoelastic characteristics of a subject.
An ultrasonic diagnostic apparatus is known as a subject information acquisition apparatus. The ultrasonic diagnostic apparatus noninvasively images information inside the subject, so that the ultrasonic diagnostic apparatus is widely used in the medical field.
In the ultrasonic diagnostic apparatus, an ultrasonic crosssectional image in which a structure in a biological body is imaged from reflection echo generated by differences of reflectivity of tissues and a Doppler ultrasound image in which a blood flow velocity or the like is imaged by using an ultrasonic Doppler effect generated by a blood flow are conventionally used.
Further, recently, hardness of tissue is measured by an ultrasonic wave and the hardness is used to diagnose the tissue. This is because the hardness of tissue deeply related to a pathological state. For example, it is known that a sclerosing cancer such as breast cancer and thyroid cancer is harder than a normal tissue and a benign tumor.
In recent years, it is reported that viscosity is different between cancer and a benign tumor, so that not only tissue hardness measurement, but also viscoelastic characteristics evaluation including viscosity characteristics evaluation is required for tissue diagnosis. Therefore, it is required to measure distribution of viscoelastic characteristics in a biological body by using ultrasonic wave, form the distribution into an image, combine the image with a conventional image of tissue structure, and use the combined images to diagnose cancer or the like in a tissue.
As an ultrasonic diagnostic apparatus for calculating the viscosity of the subject, PTL 1 discloses an apparatus that uses an ultrasonic probe in which a pressure sensor is provided on a surface of a transducer transmitting and receiving an ultrasonic wave, measures strain distribution by ultrasonic wave, measures pressure distribution applied to the subject by the pressure sensor, and calculates values of elasticity and viscosity from the strain distribution and the pressure distribution.
As a method for measuring pressure, PTL 2 discloses an ultrasonic diagnostic apparatus in which a pressure measuring deformation portion whose elastic modulus is known is sandwiched between an ultrasonic probe and a biological tissue to measure pressure. The ultrasonic diagnostic apparatus measures deformation of the pressure measuring deformation portion by ultrasonic wave, calculates pressure (stress) applied to the pressure measuring deformation portion from a relationship between elastic modulus and strain, calculates elastic modulus distribution in the subject from the stress and strain distribution in the subject, and displays the elastic modulus distribution.
However, in the configuration described in PTL 1, the structure of the ultrasonic probe is complicated. Further, generally when strain distribution in a tissue is measured by using ultrasonic wave, applied pressure is very small, so that it is difficult for a normal pressure sensor to measure the pressure. Furthermore, if the pressure sensor is inserted immediately below the ultrasonic probe, efficiency and sensitivity of transmitting and receiving ultrasonic wave degrade, so that there is a problem that measurement of deep portion cannot be performed.
In PTL 2, although the configuration of the apparatus is simple, only the elasticity distribution in the subject is displayed and a method for measuring the viscosity distribution in the subject is not disclosed. Therefore, an apparatus which has a simple configuration and which can calculate viscosity of the subject is required.
PTL 1 Japanese Patent LaidOpen No. 2007222605
PTL 2 Japanese Patent LaidOpen No. 200566041
The present invention that can solve the above problem is a subject information acquisition apparatus for transmitting an elastic wave to a subject, receiving the elastic wave reflected from inside the subject, and acquiring information inside the subject. The subject information acquisition apparatus includes a conversion device that receives the elastic wave and converts the elastic wave into an electrical signal, a reference layer which is provided between the conversion device and the subject and whose viscosity coefficient (ηc) and elastic modulus (Ec) are known, a measuring unit that measures strain (ε) of the subject and strain (εc) of the reference layer which are generated when a pressure is applied to the subject and the reference layer by using the electrical signal, and a calculation unit that calculates a viscosity coefficient (η) of the subject by using a viscosity coefficient (ηc) and an elastic modulus (Ec) of the reference layer, the strain (ε) of the subject, and the strain (εc) of the reference layer.
The present invention can provide a subject information acquisition apparatus, such as an ultrasonic diagnostic apparatus, which can measure highlysensitive and highresolution strain distribution, calculate distributions of Young'"'"'s modulus and viscous coefficient from measurement of stress distribution based on the strain distribution, and image distribution of viscoelastic characteristics in a biological body even though having a simple apparatus configuration.
Hereinafter, an embodiment of the present embodiment will be described in detail with reference to the drawings.
In
In the embodiment of the present invention, the probe 1 periodically applies a pressure to the subject by manual operation or mechanical driving, and at the same time, the probe 1 transmits and receives an ultrasonic wave as an elastic wave to and from the subject through the viscoelasticity measurement reference layer 3 whose elastic modulus and viscosity coefficient are known, so that a strain distribution in the viscoelasticity measurement reference layer 3 and a strain distribution in the subject 4 are measured.
When a tissue has no viscosity, an amplitude value of strain periodically generated in the tissue of the subject 4 when a pressure is periodically applied to the subject 4 is represented by the following equation.
Here, E is an elastic constant (Young'"'"'s modulus) of the tissue and σ_{0 }is an amplitude of a stress (pressure) which is periodically applied. Therefore, when the tissue has no viscosity, the elasticity of the tissue can be obtained from the amplitude value of the strain. Conversely, when the viscoelasticity measurement reference layer 3 has no viscosity, the stress can be obtained from the amount of strain.
On the other hand, when a tissue has a viscosity, for example, if a periodic pressure having a triangular waveform with respect to time is applied, the amplitude value of the periodic strain generated in the tissue is represented by the following equation.
Here, τ is a ratio of the viscosity coefficient and the Young'"'"'s modulus of the tissue and T is a period of the pressure which is periodically applied. As described above, the value obtained from the measurement of the strain is a ratio of the Young'"'"'s modulus and the viscosity coefficient, so that the Young'"'"'s modulus and the viscosity coefficient cannot be obtained separately. Therefore, it is necessary to calculate the viscoelasticity of the subject 4 by using measured values of strain and stress while measuring variation of the strain and stress changing over time in time series by periodically pressing the subject 4. At this time, periodic amplitude of the strain of the viscoelasticity measurement reference layer 3 is represented by the following equation due to the viscosity in the layer.
It is understood that both (Equation 2) and (Equation 3) approximately approach (Equation 1) when the period is very large (T→∞), that is, when the measurement is performed by spending sufficient time. (Therefore, if the period of the pressure is set to be sufficiently long, at least the Young'"'"'s modulus can be measured by using (Equation 1)). However, for normal tissue and the viscoelasticity measurement reference layer 3, the ratios τ and τ_{c }of the Young'"'"'s modulus and the viscosity coefficient are a time from several seconds to several tens of seconds. If a period sufficiently larger than the time is used, the measurement period becomes a long time such as several minutes, so that it is difficult to use (Equation 1) for actual measurement. (Further, in this case, the viscosity coefficient cannot be measured.)
Next, a case in which the viscoelasticity of the subject is temporarily calculated by using a technique described in PTL 2 and a calculation of the viscoelasticity of the subject 4 on the basis of the present embodiment will be described while comparing both calculations.
On the other hand, as in the present embodiment, when the Young'"'"'s modulus Ec and the ratio τ_{c }of the Young'"'"'s modulus and the viscosity coefficient of the viscoelasticity measurement reference layer 3 are known, the stress amplitude σ_{0 }is obtained from the amplitude ε_{c }of periodic strain of the viscoelasticity measurement reference layer 3 by using (Equation 3). Specifically, the strains of the viscoelasticity measurement reference layer 3 are measured over the positions of the viscoelasticity measurement reference layer 3. The strain amplitude ε_{c }is obtained by measuring the maximum and minimum values of the strains and the stress amplitude σ_{0 }is obtained from (Equation 3).
When more general pressure is applied, the change over time of the stress applied to the tissue is represented by the following equation.
Here, σ(t) and ε_{c}(t) are changes over time of the stress and the strain of the viscoelasticity measurement reference layer 3. ε_{c }and η_{c }are the known Young'"'"'s modulus and viscosity coefficient of the viscoelasticity measurement reference layer 3.
Therefore, it is possible to obtain the change of the stress from the strain of the viscoelasticity measurement reference layer 3 measured over at least one cycle of the pressure by using (Equation 4) and obtain the stress amplitude σ_{0 }from the amplitude of the change of the stress.
In an example of the present embodiment, the strain distributions in the subject 4 and the viscoelasticity measurement reference layer 3 are measured in time series. Thereby, the changes over time of an amount of strain of a specific tissue in the subject 4 and an amount of strain of a region that presses the tissue in the viscoelasticity measurement reference layer 3 are compared and a time shift (delay time) between both amounts of strain is calculated. At this time, the delay time between both amounts of strain can be represented by the following equation.
Here, τ is a ratio of the viscosity coefficient and the Young'"'"'s modulus of the tissue in the subject 4 and τ_{c }is a ratio of the viscosity coefficient and the Young'"'"'s modulus of the viscoelasticity measurement reference layer 3. When the viscosity coefficient and the Young'"'"'s modulus of the viscoelasticity measurement reference layer 3 are known, the ratio τ of the viscosity coefficient and the Young'"'"'s modulus of the tissue are obtained from (Equation 5) and the Young'"'"'s modulus E of the subject 4 is calculated from the stress amplitude σ_{0 }obtained in the above description and the measured amplitude of strain of the subject 4 by using (Equation 2). Further, the viscosity coefficient of the subject 4 is calculated from the Young'"'"'s modulus E of the subject 4 and the ratio τ of the viscosity coefficient and the Young'"'"'s modulus.
Further, another method for obtaining the ratio T of the viscosity coefficient and the Young'"'"'s modulus of the tissue in the subject 4 will be described with reference to
When the strain and the stress are measured while a periodic pressure is applied and the values of the strain and the stress are plotted on a graph, a hysteresis loop as shown in
When (Equation 2) and (Equation 3) are applied the above equation, the following equation is established.
Here, T is a period of the periodic pressure, E and E are respectively the amplitude of strain of the viscoelasticity measurement reference layer 3 and the amplitude of strain of the subject 4 which are obtained by measurement, and τ_{c }is the ratio of the known viscosity coefficient and Young'"'"'s modulus of the viscoelasticity measurement reference layer 3. Therefore, the area S1 of the hysteresis loop of the lefthand side and values of each variable in the righthand side can be measured by measuring strain except for τ. Thus, the value of τ can be obtained by reversely solving this equation.
In particular, if the period T is set to a short time of about several seconds and periodical pressure is applied, the following approximate equation can be used.
Further, a hysteresis parameter described in (Reference Document 1) Nitta, Shiina, Ueno, “Hysteresis Parameter Imaging of Soft Tissue under QuasiStatic Deformation”, 2003 IEEE International Ultrasonics Symposium Proceedings, pp. 16061609, 2003 may be used.
The value of τ can be obtained by using the hysteresis parameter HP in the same manner as described above. When the hysteresis parameter is used, calculation is standardized by the reference area, so that there is an advantage that there is not so much effect of error in the measurement. Further, when the value of τ is obtained by reversely solving the equation, there is an advantage that multiple value function is difficult to be generated in general and process is simple.
In the description here, equations related to a case in which a periodical pressure having a triangular waveform with respect to time is applied are shown for simplicity of the description. However, the pressing method may be a method of periodically repeating pressing and releasing and the equations used in the above calculation according to a pattern of change over time may be used for an analysis formula based on the model shown in
As a measuring method of another embodiment of the present invention, it is possible to apply pressure in a steplike manner instead of periodically applying pressure and measure the strain distributions in the subject 4 and the viscoelasticity measurement reference layer 3 in time series. Pressure is applied instantaneously and thereafter the strain distributions in the subject 4 and the viscoelasticity measurement reference layer 3 are measured in time series. Even when pressure is applied instantaneously, the effects of the viscosity of the viscoelasticity measurement reference layer 3 are removed after sufficient time has elapsed. However, it is not practical because the measurement time is too long. Therefore, changes over time of the strain value are measured in a short time interval in the same manner as in the case in which pressure is applied periodically.
At this time, in the model shown in
The stress amplitude σ_{0 }is obtained from the measured strain of the viscoelasticity measurement reference layer 3 and the known Young'"'"'s modulus and viscosity coefficient by using (Equation 10). Further, the ratio τ of the viscosity coefficient and the Young'"'"'s modulus is obtained from the measured change of the strain of the subject 4. Specifically, the logarithm of the time derivative of the strain of the subject 4 is calculated by the equation described below.
Therefore, the reciprocal number of the ratio τ of the viscosity coefficient and the Young'"'"'s modulus, which is an inclination of the logarithms of the time derivative of the strain of the subject 4 at a plurality of different time points, from the logarithms of the time derivative. When the reciprocal number of the ratio τ is applied to (Equation 11), the Young'"'"'s modulus of the subject 4 is obtained, and the viscosity coefficient of the subject 4 is calculated from the Young'"'"'s modulus and the ratio τ.
A method for applying the measurement (calculation) of the Young'"'"'s modulus and the viscosity coefficient to viscoelasticity distribution measurement (calculation) of the subject 4 will be described with reference to
On the basis of the principle described above, the viscoelasticity measurement reference layer 3 whose viscosity coefficient and elastic modulus are known is provided between the conversion element array 2 that is a conversion device which receives an elastic wave reflected by the subject 4 and converts the elastic wave into an electrical signal and the subject 4, and the strain of the subject 4 and the strain of the reference layer 3 which are generated when a pressure is applied to the subject 4 and the reference layer 3 are measured on the basis of the elastic wave reflected by the subject 4. It is possible to calculate the viscosity coefficient of the subject by using the measured strain of the reference layer 3, the measured strain of the subject 4, and the known viscosity coefficient and elastic modulus of the reference layer 3.
As described above, the pressure applied to the subject 4 and the reference layer 3 is a periodic pressure including a release of the pressure, so that the viscoelasticity distribution calculation unit 8 that is a calculation unit can calculate the viscosity coefficient (1) of the subject 4 by using the change over time of the strain value of the subject 4 and the change over time of the strain value of the reference layer 3 which are generated by the periodic pressure.
Specifically, for example, as described in the above equation 6 to 8, the viscoelasticity distribution calculation unit 8 that is a calculation unit can calculate the viscosity coefficient (η) of the subject by using the area of the hysteresis loop drawn with the strain values of the subject 4 and the strain values of the reference layer 3 generated by the periodic pressure, or as described in the above equations 5, using the delay time between the change over time of the strain value of the subject 4 and the change over time of the strain value of the reference layer 3.
As described above, the pressure applied to the subject 4 and the reference layer 3 is instantaneously increased to a certain pressure and thereafter the certain pressure is held for a certain time period, so that the viscoelasticity distribution calculation unit 8 that is a calculation unit can calculate the viscosity coefficient (η) of the subject 4 by using the strain value of the subject 4 and the strain value of the reference layer 3 after the pressure applied to the subject 4 and the reference layer 3 is instantaneously increased to the certain pressure.
Although, in the above embodiment, the Young'"'"'s modulus is used as elasticity characteristics and the viscosity coefficient is used as viscosity characteristics, an elastic coefficient such as a stiffness coefficient and a pressure elastic coefficient and a viscous coefficient such as a viscosity and a viscosity coefficient can be used. Further, here, a calculation method of the viscosity characteristics using a storage elastic modulus as the elasticity characteristics and a loss elastic modulus as the viscosity characteristics will be described. When a stress σ(t) is applied to a material having viscoelastic characteristics, the phase of the strain ε(t) is delayed from the phase of the strain.
σ(t)=σ_{0 }sin ωt
ε(t)=ε_{0 }sin(ωt−δ)
δ represents a delay of the phase. At this time, a storage elastic modulus G′, a loss elastic modulus G″, and an absolute dynamic elastic modulus G can be represented by the equation below.
Here, when the amplitude of strain of the viscoelasticity measurement reference layer is ε_{c0}, the amplitude of strain of the subject is ε_{t0}, and the phase difference between the strain of the viscoelasticity measurement reference layer and the strain of the subject is δ_{tc}, these three values can be measured by an ultrasonic wave. A storage elastic modulus of the viscoelasticity measurement reference layer G_{c}′, a loss elastic modulus of the viscoelasticity measurement reference layer G_{c}″, an absolute dynamic elastic modulus G_{c}, and a phase of strain of the viscoelasticity measurement reference layer δ_{c }are separately measured, so that these values become known parameters. When the values that can be measured by an ultrasonic wave and the known parameters are used, a phase of strain of the subject δ_{t}, a storage elastic modulus of the subject G_{t}′, a loss elastic modulus of the subject G_{t}″ can be calculated by the equation below.
In other words, the viscosity characteristics (loss elastic modulus) of the subject can be calculated by using the elasticity characteristics (storage elastic modulus) of the viscoelasticity measurement reference layer, the viscosity characteristics (loss elastic modulus) of the viscoelasticity measurement reference layer, and the strains of the viscoelasticity measurement reference layer and the subject. The storage elastic modulus and the loss elastic modulus calculated as described above are values that can be measured by a normal viscoelastic characteristic measuring device. For example, these values can be directly compared with mechanical measurement results of viscoelasticity measured using an extracted tissue.
Next, an example of the present invention will be described along with the constituent elements of the above embodiment on the basis of an example 1.
Hereinafter, the example 1 of the present invention will be described with reference to the drawings.
[Specific Configuration Example of Viscoelasticity Measurement Reference Layer]
The viscoelasticity measurement reference layer 3 is formed of a material which propagates but does not absorb or scatter an ultrasonic wave as an elastic wave and which has shape retention property, flexibility, and appropriate elasticity. The propagation property of ultrasonic wave is important because image is formed by transmitting and receiving an ultrasonic wave through the viscoelasticity measurement reference layer 3. If the ultrasonic wave propagation is attenuated by absorption or scattering, efficiency of the transmission and reception degrades. Therefore, the viscoelasticity measurement reference layer 3 is formed of a uniform material which is transparent with respect to ultrasonic wave and which has a small attenuation constant and a propagation speed (sound speed) near the average sound speed in the subject 4. If the acoustic impedance of the viscoelasticity measurement reference layer 3 is largely different from that of the subject 4, reflection occurs on the boundary surface and efficiency of the transmission and reception for creating an image degrades, so that the acoustic impedance of the viscoelasticity measurement reference layer 3 is not so much different from that of the subject 4. However, if the acoustic impedance of the viscoelasticity measurement reference layer 3 is the same as that of the subject 4, it is difficult to detect the boundary surface described later, so that the acoustic impedance of the viscoelasticity measurement reference layer 3 is different from that of the subject 4 by an appropriate amount of acoustic impedance. The viscoelasticity measurement reference layer 3 has shape retention property, flexibility, and appropriate elasticity to measure the distribution of the stress (pressure) by using the amount of deformation of the viscoelasticity measurement reference layer 3. As a material having the above characteristics, aqueous gel such as polyvinyl alcohol, polyurethane, and rubber material can be used. However, any material having the above characteristics can be used. However, a material is used whose elastic coefficient such as Young'"'"'s modulus, stiffness coefficient, and pressure elastic coefficient and viscous coefficient such as viscosity and viscosity coefficient are known. The thickness of the viscoelasticity measurement reference layer 3 should be thin in the viewpoint of the propagation of the ultrasonic wave. However, the viscoelasticity measurement reference layer 3 should be appropriately deformed to measure the distribution of the stress (pressure). Further, it is necessary for the viscoelasticity measurement reference layer 3 to be put in an area where the strain distribution is detected. Therefore, the thickness of the viscoelasticity measurement reference layer 3 is 0.1 mm to 50 mm, and preferably 1 mm to 10 mm.
[Creating Reception Beam Signal and Displaying BMode Image]
The probe 1 mechanically or electronically performs beam scan and transmits and receives an ultrasonic wave to and from the subject 4. The conversion element array 2 is a conversion device in which conversion elements which are driven by a drive waveform from the transmitting unit 5 and generate an ultrasonic wave are arranged. The transmitting unit 5 generates a drive waveform to each conversion element of the conversion element array 2 and adjusts the drive timing of each conversion element on the basis of control from the control unit 11. A synthesized ultrasonic wave obtained by synthesizing ultrasonic waves generated from each conversion element by adjusting the drive timing of the conversion elements forms an ultrasonic wave transmission beam that converges to a predetermined point. In other words, the control unit 11 controls the transmitting unit 5, so that the ultrasonic wave transmission beam having a convergence point at a desired depth and a directivity in a desired direction is transmitted from the probe 1. The control unit 11 performs sector scan in which the direction of the ultrasonic wave transmission beam is swung and beam scan is performed. Or, the conversion elements that are driven on the conversion element array 2 are limited to form a transmission opening and the transmission opening is moved, so that linear scan is performed in which the ultrasonic wave transmission beams scan in substantially parallel.
The receiving unit 6 includes a receiving circuit and a phasing addition circuit. The receiving circuit amplifies received signals that are received by each conversion element of the conversion element array 2 and converts the received signals into a plurality of time series digital signals by performing AD conversion. The phasing addition circuit is to form an ultrasonic wave reception beam. The phasing addition circuit provides and adds (phasing of) a delay time controlled by the control unit 11 to the digitalized received signals to form a reception beam signal. The reception beam signal is a signal in which received timing of the signals received by the conversion elements are adjusted and the signals are summed up and it is possible to form an ultrasonic wave reception beam having directivity with respect to the receiver sensitivity and one or a plurality of convergence points. The control unit 11 performs beam scan of the ultrasonic wave reception beam in synchronization with the ultrasonic wave transmission beam. As a result, a reception beam signal is obtained by transmitting an ultrasonic wave transmission beam into the subject 4 from the probe 1 and receiving an ultrasonic reflection echo signal reflected from an internal tissue of the subject 4 in an ultrasonic wave reception beam. By the beam scan of the control unit 11, a plurality of reception beam signals corresponding to ultrasonic reflection echo signals with respect to a specific point or a specific direction in the subject 4 are obtained. The reception beam signals are arranged corresponding to the beam scan, so that a crosssectional image corresponding to the intensity of the ultrasonic reflection echo signals in the subject 4 is obtained. The control unit 11 repeats the beam scan, so that crosssectional images at different times of day are obtained.
A plurality of reception beam signals generated by the receiving unit 6 are inputted into the imaging unit 9. The imaging unit 9 forms a gray scale Bmode crosssectional image reflecting the crosssectional structure in the subject 4 from the plurality of reception beam signals. The imaging unit 9 arranges the plurality of reception beam signals according to the beam scan and performs signal processing, such as gain correction, filter processing, envelope demodulation, and log compression, on the reception beam signals. Further, the imaging unit 9 may include a digital scan conversion circuit for converting the signals into a display signal to the display unit 10 by performing image processing such as unsharp mask and image filter for displaying image after the signal processing and a DA conversion circuit for analog video signal. The imaging unit 9 causes the display unit 10 to display a Bmode crosssectional image according to display control from the control unit 11.
The display unit 10 is a display apparatus such as CRT and LCD. The display unit 10 displays an image based on the display signal from the imaging unit 9.
[Specific Configuration Example of a Strain Detecting Unit which is a Measuring Unit for Measuring Strain]
The plurality of reception beam signals generated by the receiving unit 6 are also inputted into the strain distribution detecting unit 7. The strain distribution detecting unit 7 that is a measuring unit for measuring strain detects strain distributions in the subject 4 and the viscoelasticity measurement reference layer 3 on the basis of the plurality of reception beam signals.
The strain distribution can be detected by performing correlation calculation between reception beam signals at the same position obtained in different beam scan periods in the same manner as in the known colorflow Doppler method and tissue tracking method.
The strain distribution detecting unit 7 has a memory circuit for storing the plurality of reception beam signals and a correlation calculation circuit. The plurality of reception beam signals obtained when the control unit 11 repeats the beam scan are stored in the memory circuit. A set of reception beam signals obtained by one beam scan forms frame reception beam signal data corresponding to a crosssectional image at a certain time point. Different frame reception beam signal data corresponds to a beam scan at a different time point. (Each frame reception signal data includes reception beam signals from a plurality of different positions.) Reception beam signals at a corresponding position between different frame reception beam signal data are ultrasonic reflection echo signals from the same position between beam scans at different time points. Therefore, a displacement at the corresponding position generated between beam scans at different time points can be measured by performing correlation calculation of these reception beam signals. The strain distribution is calculated from the displacement distribution in taking the difference along beam depth direction.
The correlation calculation circuit has a Hilbert transform filter at the input portion thereof. The correlation calculation circuit converts two reception beam signals of beam scans at different time points corresponding to the same position into an analysis signal and performs complex correlation calculation of these signals. Or, by using an IQ signal obtained by quadrature detection of the reception beam signal, an instantaneous phase delay in time series between two reception beam signals of beam scans at different time points corresponding to the same position is calculated, and displacement may be obtained by converting the instantaneous phase delay into a sound speed. Further, the accuracy of the strain measurement may be improved by using a combined autocorrelation method or the like known in the related technical field.
As a specific configuration of the memory circuit and the correlation calculation circuit, the configuration described below may be used.
A first example includes a memory unit that stores a plurality of frame reception beam signal data and a plurality of correlation calculation circuits corresponding to beam positions used in one beam scan. Two frame reception beam signal data corresponding to beam scans at different time points are selected from a plurality of frame reception beam signal data stored in the memory circuit and reception beam signals corresponding to the two frame reception beam signal data are inputted into correlation calculation circuits corresponding to each beam position. A displacement amount distribution in the beam depth direction at a corresponding beam position can be obtained by performing correlation calculation in the correlation calculation circuit. In this configuration, beam scan is repeatedly performed and correlation calculation is performed by selecting frame reception beam signal data from a plurality of obtained frame reception beam signal data, so that it is possible to measure a displacement between beam scans at any two time points. The amount of strain can be measured by calculating a difference between the strains in the depth direction.
Here, a plurality of correlation circuits that perform onedimensional correlation calculation for each beam position are used. However, it is possible to use a correlation circuit that performs twodimensional correlation calculation for frame reception beam signal data at different time points of beam scan.
Another example includes a plurality of memory circuits and a plurality of correlation calculation circuits corresponding to each beam position used in beam scan. Each of the plurality of memory circuits includes two FIFO type memories and a TMP memory for primarily storing a correlation calculation result. While the beam scan is repeatedly performed, reception beam signals are first stored in first FIFO type memories for each beam position. An output of the first FIFO type memory is inputted into a second FIFO type memory and the correlation calculation circuit. Further, an output of the second FIFO type memory is inputted into the correlation calculation circuit. Each of the two FIFO type memories has a memory capacity for storing a reception beam signal obtained at a beam position corresponding to each beam scan. Therefore, at a certain time point in repetitive beam scans, two reception beam signals at beam positions corresponding to two temporally consecutive beam scans are stored. The two reception beam signals are calculated by the correlation calculation circuit and the displacement amount distribution in the beam depth direction at a corresponding beam position is obtained. The displacement amount distribution obtained here is generated between the two temporally consecutive beam scans in the repetitive beam scans. To obtain displacement generated between a plurality of beam scans, the calculation result of the correlation calculation circuit is accumulatively added. Therefore, the calculation result of the correlation calculation is added to data in the TMP memory and the addition result is stored in the TMP memory. Further, the strain distribution is calculated by calculating differences of the displacement distribution in the depth direction. This configuration is suited to implement the strain distribution detecting unit 7 by pipeline type processing suitable to real time processing. In particular, in this configuration, the capacity of the memory circuit can be smaller than that of the example described above.
In both examples described above, the strain distribution in the beam depth direction generated between different beam scans is obtained from a plurality of correlation calculation circuits corresponding to each beam position. The strain distributions are arranged in each beam position, so that a strain distribution of twodimensional crosssectional image generated between different beam scans is obtained. In addition to the two examples described above, it is possible to detect a strain distribution of twodimensional crosssection using frame reception beam signal data obtained by beam scans at different time points by using a method known in the related art.
In the present example, in particular, an area including not only the subject 4, but also the viscoelasticity measurement reference layer 3 is used to measure the strain distribution. More specifically, the area is set so that the amount of displacement in a boundary portion between the subject 4 and the viscoelasticity measurement reference layer 3 is measured. Thereby, the strain distribution of the viscoelasticity measurement reference layer 3 is calculated. To detect the amount of displacement in the viscoelasticity measurement reference layer 3, a reception beam signal is cut out so that the reception beam signal includes a reflection echo of ultrasonic wave from the boundary portion and a series of correlation calculations described above may be performed. However, the calculation circuit as described below can be separately provided.
Generally, inside of the viscoelasticity measurement reference layer 3 is substantially transparent with respect to ultrasonic waves, so that a reflection echo signal of ultrasonic wave is very small and the reception beam signal is also very small. Therefore, the correlation calculation result of the reception beam signal in this area includes an error and the result may be unstable. On the other hand, the boundary portion between the subject 4 and the viscoelasticity measurement reference layer 3 locally generates a large reflection echo signal due to a difference of acoustic impedance between the subject 4 and the viscoelasticity measurement reference layer 3, so that the reception beam signal and the correlation calculation result thereof can be accurately detected. The viscoelasticity measurement reference layer 3 is formed of a material uniform in the thickness direction and has a strain distribution constant in the thickness direction, so that a strain distribution in the layer is obtained by using the strain values obtained from the displacement in the boundary portion.
Therefore, a threshold circuit is provided at an input portion of the correlation calculation circuit and a reception beam signal whose intensity is smaller than a certain value is set to 0, so that the correlation calculation is performed by masking a front portion and a portion corresponding to the thickness of the viscoelasticity measurement reference layer 3 of each reception beam signal. Thus, instability as described above can be eliminated.
Further, a portion in which the amplitude intensity of the reception beam signal is smaller than a certain value is determined to be inside of the viscoelasticity measurement reference layer 3, so that the boundary portion of the reception beam signal is directly detected and displacement in the boundary portion between beam scans whose boundary portion position is different may be obtained from a difference between the frame reception beam signal data of the beam scans. The strain distribution of the viscoelasticity measurement reference layer 3 may be calculated from displacement distribution of the boundary portion between the subject 4 and the viscoelasticity measurement reference layer 3, in taking the difference along the depth direction.
As described above, the amount of strain of the viscoelasticity measurement reference layer 3 is calculated from a front portion of each reception beam signal and the strain distribution of the subject 4 along the depth direction is calculated from a portion following the front portion. The strain distributions along the depth direction calculated from the reception beam signals are arranged at corresponding beam positions, so that a crosssection strain distribution of the subject 4 and a strain distribution of the viscoelasticity measurement reference layer 3 along the boundary between the viscoelasticity measurement reference layer 3 and the subject 4 are measured. These strain distributions are inputted from the strain distribution detecting unit 7 to the viscoelasticity distribution calculation unit 8, which is a unit for calculating the viscosity coefficient.
[Viscoelasticity Distribution Calculation Unit which is a Unit for Calculating the Viscosity Coefficient]
The viscoelasticity distribution calculation unit 8, which is a unit for calculating the viscosity coefficient, calculates the elastic modulus distribution and the viscosity coefficient distribution in the subject 4 from the strain distributions of the viscoelasticity measurement reference layer 3 and the subject 4 and the elastic coefficient and the viscous coefficient of the viscoelasticity measurement reference layer 3 according to the principle described above.
The viscoelasticity distribution calculation unit 8 may be formed by using a microprocessor, a memory, a control circuit that controls the microprocessor and the memory, and a bus circuit that transmits data, or may be implemented as processing software on a programmable processing circuit such as a generalpurpose PC and an FPGA.
Therefore, the process in the viscoelasticity distribution calculation unit 8 will be described below. When linear scan of an ultrasonic beam is performed to measure strain, each position of the ultrasonic beams corresponding to each reception beam signal substantially corresponds to a direction in which pressure is applied. Therefore, the amount of strain of the viscoelasticity measurement reference layer 3 obtained for each reception beam signal and the strain distribution of the subject 4 along the depth direction correspond to strains in the same column aligned in the vertical direction described in
For each frame reception beam signal data obtained by one beam scan, a crosssectional strain distribution data set at the frame time of the beam scan is obtained. The crosssectional strain distribution data set includes the amount of strain of the viscoelasticity measurement reference layer 3 and a string of strain values along the depth direction of the subject 4 obtained from each reception beam signal at the time corresponding to one beam scan. The data set of the amount of strain of the viscoelasticity measurement reference layer 3 and the strain values along the depth direction of the subject 4 for each reception beam signal is referred to as a beam strain distribution data set. The crosssectional strain distribution data set includes a plurality of beam strain distribution data sets corresponding to each reception beam signal at each time. In the crosssectional strain distribution data set, each beam strain distribution data set corresponds to an ultrasonic beam position in a beam scan. A plurality of crosssectional strain distribution data sets corresponding to scans at different times can be obtained by repeating beam scan.
Crosssectional strain distribution data sets at each time can be obtained by repeating beam scan while periodically applying pressure.
The area S1 of viscoelasticity measurement reference layer strain−subject strain hysteresis loop is calculated from a plurality of crosssectional strain distribution data sets obtained when one cycle or a plurality of cycles of application of pressure is completed. The maximum strain ε_{cmax }and the minimum strain ε_{cmin }of the viscoelasticity measurement reference layer 3 are calculated for each ultrasonic beam position from a corresponding beam strain distribution data set. (A beam strain distribution data set corresponding to a desired ultrasonic beam position is extracted from the crosssectional strain distribution data sets of each time point obtained during one cycle of application of pressure, and a maximum value and a minimum value of the strain of the viscoelasticity measurement reference layer 3 in the data set are obtained. The above operation is performed for each ultrasonic beam position. An average value of the strain values of the viscoelasticity measurement reference layer 3 in one cycle of application of pressure may be calculated by using the crosssectional strain distribution data sets obtained during a plurality of cycles of application of pressure, and the maximum value and the minimum value in one cycle may be obtained. As the maximum value and the minimum value, values in the beam strain distribution data set, that is, values obtained at a time point when the data is obtained at each scan, may be used. Or, values may be interpolated between these values with respect to time, and an interpolated value of a time point between time points when data is actually obtained may be estimated as the maximum value or the minimum value.) For each ultrasonic beam position, viscoelasticity measurement reference layer strain amplitude ε_{c }is calculated from a difference between the maximum strain ε_{cmax }and the minimum strain ε_{cmin}. Next, the area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop is calculated for each point in the depth direction in each ultrasonic beam position. First, an outward area is obtained by, from a strain value of the subject and a strain value of the viscoelasticity measurement reference layer at the time point along from a time when the minimum strain ε_{cmin }is given to a time when the maximum strain ε_{cmax }is given. (The strain value of the subject is given for a point at each depth in the beam strain distribution data and the strain value of the viscoelasticity measurement reference layer is given for each beam strain distribution data to which the point belongs. Therefore, a pair of the strain value of the subject and the strain value of the viscoelasticity measurement reference layer at each time for a point at each depth of the beam position is obtained. For example, trapezoidal integration approximation or Simpson integral approximation and the area S1 are obtained from the pair of these values of each time point.) Next, similarly, a return area is obtained along from the time when the maximum strain ε_{cmax }is given to the time when the minimum strain ε_{cmin }is given and a difference between the outward area and the return area is calculated, so that the area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop is obtained. At the same time, a maximum strain ε_{max }and a minimum strain ε_{min }of the subject are calculated at each point and viscoelasticity measurement reference layer strain amplitude E is calculated from the difference between them. Then, the elastic modulus and the viscosity coefficient at each point are calculated by reversely solving (Equation 8).
Accuracy of data for calculating the hysteresis loop area can be improved by repeating a plurality of cycles of pressure application. Specifically, the area is calculated by using all pairs of the viscoelasticity measurement reference layer strain and the subject strain which are obtained at each point and each time. In this case, the pairs obtained in a process in which the viscoelasticity measurement reference layer strain increases as time advances are used for calculating the outward area and the pairs obtained in a process in which the viscoelasticity measurement reference layer strain decreases as time advances are used for calculating the return area. Regarding the maximum strain and the minimum strain, similarly, the maximum value and the minimum value of the strains measured across the entire measurement time may be used. The area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop may be calculated for each pressure cycle and an average value of the area S1 may be used.
Further, a part of the amounts of strain of the viscoelasticity measurement reference layer 3 or an average strain value in the beam strain distribution data set is calculated and a pressure cycle time T may be calculated by measuring change over time of the average strain value. In other words, the pressure cycle time T can be obtained by calculating a time difference between time points of the maximum value and the minimum value of the average strain value changing over time in a process in which periodic pressure is applied. When a plurality of cycles of pressure is applied, the pressure cycle time T may be obtained by calculating an average value of the pressure cycle times. Although a part or an average of the amounts of strain of the viscoelasticity measurement reference layer 3 is used as the average strain value, an amount of strain in a specific position in the subject 4 or an average of the amounts of strain in the subject 4 may be used.
Distributions of the elastic modulus and the viscosity coefficient in a crosssection of the subject are calculated by performing the above calculation for all points in the depth direction in the ultrasonic beam position and for all the ultrasonic beam positions in the beam scan.
The calculated distribution data of elastic modulus and viscosity coefficient is inputted into the imaging unit 9. The imaging unit 9 forms a viscoelasticity crosssectional image reflecting the viscoelasticity in the subject 4 from a plurality of distribution data of elastic modulus and viscosity coefficient. In the same manner as forming the Bmode crosssectional image, the imaging unit 9 performs signal processing such as gain correction and filter processing and image processing such as unsharp mask and image filter to form the viscoelasticity crosssectional image. Further, the imaging unit 9 causes the display unit 10 to display the viscoelasticity crosssectional image along with the Bmode crosssectional image by style such as switching, simultaneous parallel, and superimposed synthetic according to display control from the control unit 11. Further, the imaging unit 9 may generate a strain crosssectional image on the basis of the strain distribution data in the subject 4 generated by the strain distribution detecting unit 7 and cause the display unit 10 to display the viscoelasticity crosssectional image along with the Bmode crosssectional image by style such as switching, simultaneous parallel, and superimposed synthetic according to display control from the control unit 11. The elastic modulus and the viscosity coefficient may be assigned to any one of hue, chroma, and brightness respectively and displayed. By such a display, it is possible to intuitively observe the two physical properties such as the elastic modulus and the viscosity coefficient. Further, when the elastic modulus and the viscosity coefficient are superimposed on the Bmode crosssectional image, a degree of transparency which reflects one or both of the elastic modulus and the viscosity coefficient may be set and superimposed. By such a superimposition, it is possible to efficiently observe an area having a noticeable elastic modulus or viscosity coefficient.
[Viscoelasticity Distribution Measurement Operation and Display of Viscoelasticity CrossSectional Image]
Next, the process performed by the control unit 11 when the viscoelasticity distribution measurement is performed will be described. As described above, the control unit 11 performs control of ultrasonic wave transmission/reception scan and causes the display unit 10 to display the Bmode crosssectional image. At this time, the control unit 11 may cause the strain distribution detecting unit 7 to operate and display the strain crosssectional image on the display unit 10 along with the Bmode crosssectional image or instead of the Bmode crosssectional depending on an input.
When a user disposes the probe 1 on a predetermined portion of the subject 4 while checking the Bmode crosssectional image on the display unit 10 and turns on the control switch 12, the control unit 11 starts the process of the viscoelasticity distribution calculation unit 8 and the viscoelasticity distribution measurement operation starts. While the user applies periodical pressure to the subject 4 via the probe 1, the viscoelasticity distribution calculation unit 8 operates and generates a crosssectional strain distribution data set for each time point. When a plurality of cycles of pressure is applied, the user inputs a cycle period signal from the control switch 12 at the start of a cycle and the end of a cycle. The control unit 11 calculates the pressure cycle T from the inputs of the cycle period signals and inputs the pressure cycle T into the viscoelasticity distribution calculation unit 8. When the periodical pressure application is completed, by turning off the control switch 12, generation of the crosssectional strain distribution data set in the viscoelasticity distribution calculation unit 8 ends, the area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop is calculated, and the Distributions of the elastic modulus and the viscosity coefficient are calculated. When the calculation process is completed, the control unit 11 causes the imaging unit 9 to operate and generate a viscoelasticity crosssectional image and causes the display unit 10 to display the viscoelasticity crosssectional image according to a display style inputted in advance.
Although, here, start and end of the pressure cycle are inputted from the control switch 12, it is not necessary to input the cycle period signal if the pressure cycle T is obtained from a change over time of one of strains in the obtained crosssectional strain distribution data set in the viscoelasticity distribution calculation unit 8 as described above.
Further, although the area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop is calculated after all the crosssectional strain distribution data sets are obtained, it is possible to specify the number of data acquisition pressure cycles and start calculation of the area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop, the elastic modulus, and the viscosity coefficient after the specified number of pressure cycles are completed. The area of the hysteresis loop, the elastic modulus, and the viscosity coefficient are calculated while the crosssectional strain distribution data set to be used is updated, so that the distributions of the elastic modulus and the viscosity coefficient can be calculated over time, and at the same time, the viscoelasticity crosssectional image can be updated for each pressure cycle.
Although, in the above description, the elastic modulus and the viscosity coefficient are calculated by using the area S1 of the viscoelasticity measurement reference layer strain−subject strain hysteresis loop, the elastic modulus and the viscosity coefficient may be calculated by using the hysteresis parameter.
It is possible to calculate the viscoelasticity measurement reference layer strain amplitude ε_{c }and calculate the stress amplitude σ_{0 }at a beam position corresponding to each beam strain distribution data by using (Equation 3). Thereafter, a delay time between the change of the viscoelasticity measurement reference layer strain and the strain change in the subject at each point may be calculated from changes over time of the viscoelasticity measurement reference layer strain at the beam position and the strains in the subject arranged in the depth direction in the beam position, and the elastic modulus and the viscosity coefficient may be calculated by using (Equation 5). (The delay time may be obtained from a difference of times when the strain value of subject and the strain value of viscoelasticity measurement reference layer reach the maximum value and/or the minimum value respectively. Here, the strain value of subject is obtained for each point in the subject, and a strain value of viscoelasticity measurement reference layer obtained with respect to a beam position to which the point belongs is used.) Or, the delay time may be obtained by performing a crosscorrelation calculation on timeseries data of the strain of the viscoelasticity measurement reference layer and the strain of the subject which change for each point.
Although, in the present example, a user applies periodic pressure by using a handheld type probe, it is possible to automatically apply periodic pressure by attaching an excitation unit such as a micro motor to the handheld type probe or providing a separate pressing unit for pressing the probe and driving the excitation unit or the separate pressing unit from the control unit 11. In this configuration, it is possible to accurately control a change over time of pressure application, such as the pressure cycle τ which is periodically applied, so that the calculation accuracy of the elastic modulus and the viscosity coefficient can be improved.
Furthermore, although, in the above description, the viscoelasticity measurement reference layer 3 is included in the probe 1, the viscoelasticity measurement reference layer 3 can be formed as an attachable and detachable coupler separate from the probe 1. The shape and the viscoelasticity of the coupler can be selected according to a region to be diagnosed. In this case, in the control unit 11, the elastic modulus and the viscosity coefficient can be inputted or selected according to the type of the coupler, and the viscoelasticity distribution calculation unit 8 may calculate the distribution of the elastic modulus and the viscosity coefficient of the subject by using the selected values.
In the present example, pressure is not periodically applied but applied in a steplike manner. The difference from the example 1 is mainly the process of the viscoelasticity distribution calculation unit 8, so that the process of the viscoelasticity distribution calculation unit 8 will be mainly described and only the difference will described.
When the control switch 12 of the probe 1 is turned on, the control unit 11 controls the viscoelasticity distribution calculation unit 8 so that the crosssectional strain distribution data set is stored corresponding to each time point of the beam scan. The probe 1 is instantaneously pressed in a steplike manner to a certain pressure level by a manual operation or a pressing unit and the pressure is held until data acquisition is completed. The viscoelasticity distribution calculation unit 8 uses (Equation 10) and a change over time of strain of the viscoelasticity measurement reference layer 3 of each beam strain distribution data set obtained at different time points thereafter and calculates the stress amplitude σ_{0 }at a corresponding beam position (When Equation 10 is used, the stress amplitude σ_{0 }may be obtained by applying the strain value of the viscoelasticity measurement reference layer 3 obtained after a specific time interval from the steplike pressure application, a known Young'"'"'s modulus of the viscoelasticity measurement reference layer 3, and a value of τ_{c }obtained from the Young'"'"'s modulus and a known viscosity coefficient to (Equation 10)). However, it is desired to improve measurement accuracy by using time intervals of a plurality of points and fitting to (Equation 10).
By using subject strain values in each beam strain distribution data set in a plurality of crosssectional strain distribution data sets obtained at different time points, a change over time of subject strain at each position in the depth direction at corresponding beam positions is extracted. The ratio τ of the Young'"'"'s modulus and the viscosity coefficient at each point in the subject is calculated from the change over time at each point by using (Equation 12). First, to calculate a differential of the subject strain, an approximate time differential value of the strain is calculated by dividing an increment of the subject strain for each specific time interval from the steplike pressure application corresponding to each crosssectional strain distribution data set by the time interval, and 1/τ is obtained from the inclination of the change over time of the approximate time differential value of the strain according to (Equation 12). Also at this time, it is desired to obtain 1/τ from the approximate differential value of strain at a plurality of time points by fitting. At the same time, the value of the following expression is obtained from an intercept of a change over time of the approximate time differential value of the strain.
Then, the Young'"'"'s modulus E is obtained from T obtained in the above calculation and the stress amplitude σ_{0}.
Distributions of the elastic modulus and the viscosity coefficient in a crosssection of the subject are calculated by performing the above calculation for all points in the depth direction in the ultrasonic beam position and for all the ultrasonic beam positions in the beam scan.
The calculated distributions of the elastic modulus and the viscosity coefficient are inputted into the imaging unit 9 in the same manner as in the example 1 and the viscoelasticity crosssectional image is displayed on the display unit 10.
Although, here, the approximate time differential value is obtained by dividing an increment of the subject strain by the time interval, a highorder central difference, forward difference, or backward difference may be used.
In the present example, the Young'"'"'s modulus and the viscosity coefficient of the subject are obtained from a change over time of strains in the subject and the viscoelasticity measurement reference layer against the steplike pressure application. In this case, regarding necessary time series time intervals, a sufficiently accurate calculation can be performed by using data at several to several tens of time points including a case in which the data is used for the fitting described above. Thereby, the necessary number of crosssectional strain distribution data sets for each time point is several to several tens, and the frame reception beam signal data used to generate the crosssectional strain distribution data sets includes several to several tens of frames. The beam scan for generating these data is completed within one to several seconds, so that there is an advantage that data used to generate viscoelasticity distribution data is acquired in a short time period. Also there is an advantage that measurement data can be acquired in a short time period without performing pressure application and pressure release, which are performed in the example 1.
While the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments. The scope of the following claims is to be accorded the broadest interpretation so as to encompass all such modifications and equivalent structures and functions.