WO2015030076A1 - 超音波画像撮像装置及び超音波画像表示方法 - Google Patents

超音波画像撮像装置及び超音波画像表示方法 Download PDF

Info

Publication number
WO2015030076A1
WO2015030076A1 PCT/JP2014/072481 JP2014072481W WO2015030076A1 WO 2015030076 A1 WO2015030076 A1 WO 2015030076A1 JP 2014072481 W JP2014072481 W JP 2014072481W WO 2015030076 A1 WO2015030076 A1 WO 2015030076A1
Authority
WO
WIPO (PCT)
Prior art keywords
boundary
condition
stress
ultrasonic image
subject
Prior art date
Application number
PCT/JP2014/072481
Other languages
English (en)
French (fr)
Inventor
松村　剛
Original Assignee
日立アロカメディカル株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 日立アロカメディカル株式会社 filed Critical 日立アロカメディカル株式会社
Publication of WO2015030076A1 publication Critical patent/WO2015030076A1/ja

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48—Diagnostic techniques
    • A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/899—Combination of imaging systems with ancillary equipment
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52023—Details of receivers
    • G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52053—Display arrangements
    • G01S7/52057—Cathode ray tube displays
    • G01S7/52071—Multicolour displays; using colour coding; Optimising colour or information content in displays, e.g. parametric imaging

Definitions

  • the present invention relates to an ultrasound probe and an ultrasound imaging apparatus, and more particularly to an ultrasound imaging apparatus and an ultrasound image display method for displaying an ultrasound image.
  • the ultrasonic imaging apparatus transmits an ultrasonic wave to an object from an ultrasonic probe in contact with the surface of the object, receives an ultrasonic reflected wave from the inside of the object, and is a reception signal thereof Based on the reflected echo signal, an ultrasonic image (such as a tomographic image) relating to the biological information of the part of the subject is displayed.
  • an ultrasonic image such as a tomographic image
  • the ultrasonic diagnostic apparatuses described in Patent Literatures 1 and 2 are elastic from elastic information based on displacement and strain at each measurement point inside the subject and stress at each measurement point that occur when a compression force is applied to the subject. An image is generated (see Patent Document 1 and Patent Document 2). Then, the examiner distinguishes lesions such as normal tissues, cancer cells, and tumors by observing the elastic image while adjusting the compression force.
  • the stress distribution is estimated by obtaining the stress distribution by assuming the boundary conditions of the subject such as the shape of the subject and the restraint state.
  • the present invention takes into account changes in the internal stress (or stress distribution) of a subject due to the conditions of the subject (shape, boundary conditions, etc.) by considering the conditions of the subject (shape, boundary conditions, etc.).
  • An object of the present invention is to provide an ultrasonic imaging apparatus capable of measuring stress distribution with high accuracy by setting the conditions (shape, boundary conditions, etc.) of the subject with high accuracy.
  • An ultrasonic imaging apparatus includes a pressure measurement unit that measures a pressure applied to a subject, a condition setting unit that sets a condition relating to a measurement region based on an ultrasonic image of the subject, and the pressure And a stress information calculation unit that calculates the stress information of the measurement site under the condition.
  • the present invention provides an ultrasonic imaging apparatus and an ultrasonic image display method capable of measuring a stress distribution with high accuracy by setting conditions (shape, boundary conditions, etc.) of a subject with high accuracy. .
  • FEM finite element method
  • FIG. 1 is a block diagram showing an ultrasonic imaging apparatus according to the present embodiment.
  • the ultrasonic imaging apparatus obtains a tomographic image of a diagnostic region of a subject using ultrasonic waves and displays an elastic image representing the hardness or softness of a living tissue.
  • a probe (ultrasonic probe) 2 used in contact with a subject 1 has a plurality of transducers that transmit and receive ultrasonic waves to and from the subject 1.
  • the probe 2 is driven by ultrasonic pulses supplied from the transmission circuit 3.
  • the transmission / reception control circuit 4 controls the transmission timing of the ultrasonic pulses that drive the plurality of transducers of the probe 2 and performs a process of forming an ultrasonic beam toward the focal point set in the subject 1. .
  • the transmission / reception control circuit 4 performs a process of electronically scanning the ultrasonic beam in the arrangement direction of the transducers of the probe 2.
  • the probe 2 receives a reflected echo signal generated from the inside of the subject 1 and outputs it to the receiving circuit 5.
  • the reception circuit 5 takes in the reflected echo signal and performs reception processing such as amplification.
  • the reflected echo signal processed by the reception circuit 5 is input to the phasing addition circuit 6.
  • the phasing and adding circuit 6 amplifies the reflected echo signal by adding the phases of the reflected echo signals received by the plurality of transducers together.
  • the reflected echo signal phased and added by the phase adjusting and adding circuit 6 is input to the signal processing unit 7 and subjected to signal processing such as gain correction, log compression, detection, contour enhancement, and filter processing.
  • the reflected echo signal processed by the signal processing unit 7 is guided to the monochrome scan converter 8 and converted into two-dimensional tomographic image data (digital data) corresponding to the scanning surface of the ultrasonic beam.
  • the signal processing unit 7 and the black and white scan converter 8 constitute a tomographic image (B-mode image) image forming unit.
  • the tomographic image data output from the monochrome scan converter 8 is supplied to the display unit 10 via the switching addition unit 9.
  • the display unit 10 displays a tomographic image (B mode image) that is an ultrasonic image.
  • the reflected echo signal output from the phasing addition circuit 6 is guided to the frame data acquisition unit 11.
  • the frame data acquisition unit 11 acquires a reflected echo signal corresponding to the scanning plane (tomographic plane) of the ultrasonic beam as frame data, and stores a plurality of frame data in a storage unit such as a memory.
  • the displacement / strain calculation unit 12 sequentially fetches a plurality of frame data stored in the frame data acquisition unit 11 at different acquisition times, and obtains the displacement or distortion of the measurement point on the tomographic plane based on the acquired pair of frame data.
  • the obtained displacement or strain is output to the elasticity information calculation unit 13 as displacement frame data or strain frame data (displacement / strain frame data).
  • the pressure measurement unit 17 measures the pressure applied to the surface of the subject 1 from the probe 2 and outputs the measured pressure to the stress calculation unit 18 as pressure data.
  • the method for obtaining the pressure of the probe 2 can be selected as appropriate, and there is a method for obtaining the surface pressure of the subject 1 using a pressure sensor, an acoustic coupler, or the like.
  • the stress calculation unit 18 calculates the stress or stress distribution inside the subject 1 as stress frame data based on the pressure data, and outputs the calculated stress frame data to the elasticity information calculation unit 13.
  • the elastic information calculation unit 13 Based on the displacement / strain frame data and the stress frame data, the elastic information calculation unit 13 obtains elastic information such as the elastic modulus (eg, Young's modulus) of the living tissue at each measurement point as elastic frame data, and the elastic frame data is elastic. The information is output to the information processing unit 14.
  • the elastic modulus eg, Young's modulus
  • the elastic information processing unit 14 subjects the elastic frame data input from the elastic information calculation unit 13 to smoothing processing and contrast optimization processing within the coordinate plane of the elastic frame data, so that a time axis between a plurality of elastic frame data is obtained.
  • Various image processing is performed by performing direction smoothing processing and the like, and the elastic frame data subjected to the image processing is sent to the color scan converter 15.
  • the color scan converter 15 takes in elastic frame data (elastic frame data indicating elastic characteristics such as elastic modulus) processed by the elastic information processing unit 14 and assigns a color code for each pixel according to a color map corresponding to the elastic modulus. Thus, a color elastic image is generated.
  • elastic frame data elastic frame data indicating elastic characteristics such as elastic modulus
  • the color elasticity image (ultrasonic image) generated by the color scan converter 15 is displayed on the display unit 10 via the switching addition unit 9. Further, the switching addition unit 9 inputs the black and white tomographic image output from the black and white scan converter 8 and the color elastic image output from the color scan converter 15 and switches both images to display one of them. A display function, a function of causing one of the two images to be translucent, adding and synthesizing and superimposing them on the display unit 10, and a function of displaying both images side by side on the display unit 10. Further, the image data output from the switching addition unit 9 is stored in the cine memory 20 under the control of the device control interface unit 19. The image data stored in the cine memory 20 is displayed on the display unit 10 under the control of the device control interface unit 19.
  • the stress calculation unit 18 acquires the shape of the subject 1 and the boundary conditions (restraint state, etc.) of the subject 1, and calculates the stress applied to each measurement point in consideration of those conditions. Ask.
  • FIG. 2 is a block diagram showing an example of the configuration of the stress calculation unit 18 according to the present embodiment.
  • the stress calculation unit 18 includes a condition setting unit 180, a condition adjustment unit 181, a model selection unit 182, and a stress information calculation unit 184.
  • the condition setting unit 180 sets conditions for the measurement site.
  • the condition adjustment unit 181 adjusts the set condition.
  • the model selection unit 182 selects from the database 183 a model for simulating the stress at the measurement site.
  • the stress information calculation unit 184 calculates stress information (stress or stress distribution) based on the model, and generates stress frame data.
  • the probe 2 presses the subject 1 and scans the subject 1 with an ultrasonic beam.
  • the receiving circuit 5 continuously receives the reflected echo signal from the scanning plane.
  • the signal processing unit 7 and the monochrome scan converter 8 form a tomographic image, and the display unit 10 displays the tomographic image via the switching addition unit 9.
  • the frame data acquisition unit 11 captures the reflected echo signal, repeatedly acquires the frame data in synchronization with the frame rate, and stores it in a built-in storage unit (such as a frame memory) in chronological order. Then, the frame data acquisition unit 11 selects a plurality of pairs of frame data continuously with a pair of frame data having different acquisition times as one unit, and outputs the frame data to the displacement / distortion calculation unit 12.
  • the displacement / distortion calculation unit 12 performs one-dimensional or two-dimensional correlation processing on the selected pair of frame data, measures the displacement of a plurality of measurement points on the scanning plane, and moves the displacement frame data or the distortion frame data (displacement / distortion frame data). Distortion frame data).
  • a displacement vector detection method for example, a block matching method or a gradient method described in JP-A-5-317313 is known.
  • the block matching method an image is divided into blocks of N ⁇ N pixels, for example, and a block closest to the focused block (for example, a block at a measurement point) in the current frame is searched from the previous frame.
  • the displacement of the measurement point is obtained based on the block position, and displacement frame data is generated. Further, the displacement is calculated by calculating the autocorrelation in the same region of the paired RF signal frame data, and the displacement frame data is generated.
  • the distortion of each measurement point is obtained based on the displacement frame data.
  • the distortion calculation is calculated, for example, by spatially differentiating the displacement, and distortion frame data is generated.
  • the pressure measurement unit 17 is applied from the probe 2 to the surface (compression surface) of the subject 1 by a pressure sensor 16 or an acoustic coupler provided between the ultrasonic transmission / reception surface of the probe 2 and the subject 1.
  • the measured pressure is measured, and the measured pressure is output to the stress calculation unit 18 as pressure data.
  • the stress calculation unit 18 calculates the stress or stress distribution inside the subject 1 as stress frame data based on the pressure data.
  • the elasticity information calculation unit 13 obtains elasticity information of the measurement point using the displacement / strain frame data and the stress frame data. For example, the elasticity information calculation unit 13 calculates the elastic modulus (for example, Young's modulus E) of each measurement point on the scanning surface from the stress at each measurement point and the distortion obtained by the elasticity information calculation unit 13, and the elasticity frame Data is output to the elastic information processing unit 14.
  • the elastic information processing unit 14 performs processing such as smoothing processing on the elastic frame data and outputs the processed data to the color scan converter 15.
  • the color scan converter 15 generates a color elasticity image based on the elasticity information.
  • a color corresponding to the elastic modulus is assigned to each pixel with a gradation of color tone by 256 gradations.
  • a black and white scan converter may be used instead of the color scan converter 15.
  • the black and white scan converter generates a black and white elastic image based on the elasticity information.
  • the elasticity information can be visually recognized by increasing the brightness of areas with a high elastic modulus and decreasing the brightness of areas with a low elastic modulus, for example, lesions such as normal tissues, cancer cells, and tumors. Can be distinguished.
  • the stress calculation unit 18 calculates the stress or stress distribution inside the subject 1 as stress frame data based on the pressure data, and the elasticity information calculation unit 13 calculates the elastic modulus (for example, Elasticity information such as Young's modulus E) is calculated.
  • the stress calculation unit 18 obtains a stress distribution assuming the shape of the measurement site and the boundary conditions (such as restraint state), which are the conditions (shape, boundary conditions, etc.) of the subject 1, and the elasticity information calculation unit 13 performs each measurement.
  • the elastic modulus E of the point is calculated.
  • the measurement site is an isotropic elastic body and using three equations, a balance equation, a strain-displacement relationship equation, and a stress-strain relationship equation (Hook's law), as shown in equation (1)
  • the elastic modulus (eg, Young's modulus E) at each measurement point of the tissue is obtained from the strain ⁇ z in the z-axis direction and the stresses ⁇ x, ⁇ y, ⁇ z in the x, y, and z-axis directions.
  • the stress distribution is a distribution of stresses ⁇ x, ⁇ y, ⁇ z in space (coordinates (x, y, z)).
  • the stress ⁇ or the stress distribution ⁇ (x, y, z) can also be obtained by using a structural analysis simulation by a finite element method (FEM).
  • FEM finite element method
  • the compression direction is the z-axis direction.
  • FIG. 3A in order to obtain the stress distribution inside the subject 1, the shape of the measurement site and the boundary conditions (restraint state, etc.) are assumed.
  • FIGS. 3B and 3C are diagrams showing the results of analyzing the stress distribution when a pressure of 1 [kPa] is applied to the compression surface in the z-axis direction.
  • FIGS. 3B and 3C show changes in the stress distributions ⁇ z and ⁇ in depth when the thickness of the subject 1 is changed to 10 [mm], 30 [mm], and 50 [mm].
  • the stress distributions ⁇ and ⁇ z are estimated by setting the measurement site conditions as described above.
  • the stress distributions ⁇ and ⁇ z vary greatly depending on the set conditions (such as the shape of the subject 1 and boundary conditions). That is, the estimation results of the stress distributions ⁇ and ⁇ z depend on the set conditions.
  • the estimation results of the stress distributions ⁇ and ⁇ z depend on the shape (thickness) of the set measurement site.
  • strain (fixed) in xyz direction which is a boundary condition (constraint condition) of the bottom surface is changed to “unconstrained”
  • the stress distributions ⁇ and ⁇ z largely change particularly near the bottom surface.
  • the ultrasonic imaging apparatus according to the present embodiment can measure the stress or the stress distribution with high accuracy by accurately setting the shape of the measurement site and the boundary condition.
  • FIG. 4 is a view showing a tomographic image (B-mode image) and an elasticity image generated by applying pressure to a breast part including a mammary gland tissue.
  • B-mode image a tomographic image
  • elasticity image generated by applying pressure to a breast part including a mammary gland tissue.
  • the condition setting unit 180 sets conditions for the measurement part.
  • FIG. 5 is a diagram showing a state in which the measurement site condition is set based on the tomographic image (B-mode image) 40.
  • the condition setting unit 180 sets the conditions (shape, boundary conditions, etc.) of the measurement site based on the tomographic image (B mode image) 40 displayed on the display unit 10.
  • the distance (thickness) of the soft tissue portion from fat to the pectoralis major muscle) that is easily deformed by compression is about 30 mm
  • the hard tissue (rib) that is difficult to deform by compression is It is a bottom surface (support surface) that supports soft tissue, and is restrained (fixed) in the xyz direction.
  • the condition setting unit 180 refers to a model (stereoscopic model) described later, such as a boundary (boundary line or boundary surface), a boundary condition, a shape (including size) of a measurement site, and a compression position.
  • a model stereographic model
  • a boundary line 50 is set at the boundary between the great pectoral muscle that is a soft tissue and the rib that is a hard tissue, and the distance from the contact surface 41 to the boundary line (rib) 50 ( Approx.
  • the boundary condition between the upper surface (contact surface 41), the bottom surface (the surface in the y direction passing through the boundary line 50) and the side surface (the surface in the z direction or the compression direction sandwiched between the bottom surface and the upper surface) are set as “compression surface (contact surface 41) ⁇ no constraint, bottom surface ⁇ constraint (fixed) in xyz direction, side surface (z direction or surface in the compression direction) ⁇ unconstrained”.
  • compression surface (contact surface 41) ⁇ no constraint, bottom surface ⁇ constraint (fixed) in xyz direction, side surface (z direction or surface in the compression direction) ⁇ unconstrained”.
  • the boundary line 50 (or the boundary surface) is set between the soft tissue and the hard tissue is that the hard tissue has a small deformation due to the compression and generally has a small movement due to the compression.
  • Boundary conditions are set according to the characteristics of each surface (upper surface, bottom surface, and side surface).
  • the boundary conditions of the bottom surface are “bottom ⁇ constrained (fixed) in xy direction, unconstrained in z direction (free)”, “bottom ⁇ unconstrained in xy direction (free), constrained in z direction (fixed)”, etc. It is set according to the characteristics of the bottom surface.
  • the characteristics of each surface may be visually identified by changing the lines that define each surface and the display of the surface.
  • the characteristics of the boundary line 50 or boundary surface
  • the characteristics of the boundary surface change according to the characteristics of the bottom surface.
  • the boundary condition is “restrained in the bottom surface ⁇ xyz direction”.
  • the color of the boundary line 50 may change to red.
  • characters representing boundary conditions for example, “restraint (fixed) in bottom surface ⁇ xyz direction”
  • symbols, figures, and the like may be attached according to the characteristics of each surface.
  • FIG. 6 is a diagram showing a state in which a boundary line 50 is set between the thyroid gland (soft tissue) and the trachea (fixed tissue) based on the tomographic image (B-mode image) 40.
  • the reason why the boundary line 50 (or the boundary surface) is set between the soft tissue and the fixed tissue is that the fixed tissue hardly moves due to the compression.
  • Soft tissue can be deformed up to fixed tissue and does not deform beyond fixed tissue. That is, the boundary line 50 (or boundary surface) between the thyroid gland (soft tissue) and the trachea (fixed tissue) is fixed in the x-axis direction.
  • the condition setting unit 180 sets “restraint (fixed) in side surface (z direction or compression direction surface) ⁇ x axis direction” as the boundary condition.
  • the condition setting unit 180 can set a curve as the boundary line 50.
  • the condition setting unit 180 sets the boundary line 50 along the curved surface of the trachea.
  • the boundary line 50 (or boundary surface) between the thyroid gland (soft tissue) and the trachea (fixed tissue) is fixed in the normal direction of the boundary line 50.
  • the condition setting unit 180 sets “restraint (fixed) in side surface (curved surface) ⁇ normal direction” as the boundary condition.
  • condition setting unit 180 can set the boundary line 50 on the upper surface (contact surface 41).
  • the boundary condition of the upper surface (contact surface 41) varies depending on whether the shape of the probe 2 is linear or convex. Further, the boundary condition of the upper surface (contact surface 41) changes depending on the curvature of the contact surface of the probe 2. Accordingly, the boundary line 50 is set according to the boundary condition (shape) of the upper surface (contact surface 41).
  • the boundary line 50 can be input as an elastic image instead of the tomographic image (B-mode image) 40.
  • condition set by the condition adjustment unit 181 may be adjusted.
  • the condition adjustment unit 181 adjusts (changes, moves, rotates) boundary conditions such as the boundary line 50 on the ultrasonic image (tomographic image or elastic image) in accordance with an input from the device control interface unit 19 such as a trackball or a mouse. Zoom in and out).
  • the condition adjustment unit 181 may adjust the angle of the boundary line 50 set horizontally and the distance from the contact surface 41.
  • condition adjustment part 181 may adjust the shape of the boundary line 50 so that the boundary line 50 may follow the curved surface of a trachea.
  • the model selection unit 182 is a model corresponding to the conditions (boundary conditions and the like) set by the condition setting unit 180 or adjusted by the condition adjustment unit 181, and a model for simulating stress at the measurement site is stored in the database 183. select. That is, the model selection unit 182 reads from the database 183 a model (stress distribution template) that matches the set measurement site conditions.
  • the database 183 stores a plurality of models according to the characteristics of the measurement site.
  • the database 183 has various shapes (such as a cylindrical shape and a quadrangular prism shape), sizes, and positions (such as a bottom surface, a side surface, and a central portion) of a support (such as a hard tissue or a fixed tissue). , And a model that combines the size, position, and direction of the compression surface.
  • a cylindrical model shown in FIG. 7A is suitable.
  • the models shown in FIGS. 7C and 7D are suitable.
  • a plurality of models are stored in the database 183 in units of 1 mm.
  • the database 183 selects a cylindrical model (for example, the model of FIG. 7A) that meets the set condition, and selects stress characteristic information (for example, stress distribution ⁇ Tx) of the selected model.
  • the stress distributions ⁇ Tx (i, j), ⁇ Ty (i, j), and ⁇ Tz (i, j) are stress distributions in the x, y, and z axis directions simulated in advance in each model, and the coordinates (i: x axis) (Direction coordinates, j: coordinates in the z-axis direction).
  • the stress distribution ⁇ T (i, j) is expressed by Expression (4).
  • the stress calculation unit 18 may calculate the stress information by a structural analysis simulation of a finite element method (FEM) using a model for simulating the stress at the measurement site.
  • FEM finite element method
  • FIG. 3 shows the stress distribution ⁇ in the center line
  • the stress distribution ⁇ T at the coordinates (i, j) can be calculated by simulating the stress distribution of all lines by the finite element method (FEM).
  • the stress distribution of the model stress distribution template
  • FEM finite element method
  • the stress distribution ⁇ Tx, ⁇ Ty, ⁇ Tz, ⁇ T of the model (stress distribution template) is a stress distribution generated when a pressure of 1 [Pa] is uniformly applied to the compression surface, and is a standardized stress distribution. Therefore, the stress distribution when the pressure P [Pa] is applied to the compression surface is calculated by multiplying the stress distributions ⁇ Tx, ⁇ Ty, ⁇ Tz, ⁇ T of the model (stress distribution template) by P.
  • the stress information calculation unit 184 calculates stress information based on the model, and generates stress frame data.
  • the stress distribution is calculated in advance based on the pressure distribution P (i, j) [Pa] on the compression surface, It may be stored in the database 183.
  • the stress calculation unit 18 calculates the stress or stress distribution inside the subject 1 as stress frame data based on the pressure data P, and outputs the calculated stress frame data to the elasticity information calculation unit 13.
  • the elasticity information calculation unit 13 inputs the displacement frame data D (i, j), the strain frame data ⁇ z (i, j), and the stress frame data ⁇ z (i, j), ⁇ (i, j), and the equation ( 7), the elastic frame data E (i, j) of the elastic information (such as Young's modulus E) is obtained, and the elastic frame data E (i, j) is output to the elastic information processing unit 14.
  • the output elastic frame data E (i, j) is displayed on the display unit 10 as a color elastic image (ultrasonic image) via the elastic information processing unit 14, the color scan converter 15, and the switching addition unit 9. .
  • the ultrasonic imaging apparatus takes into account the subject caused by the condition of the subject (shape, boundary condition, etc.) by considering the condition of the subject (shape, boundary condition, etc.). It is possible to cope with changes in the stress (or stress distribution) inside the specimen, and the stress distribution can be measured with high precision by setting the conditions (shape, boundary conditions, etc.) of the specimen with high precision.
  • the ultrasonic image display method the pressure applied to the subject is measured, the condition relating to the measurement region is set based on the ultrasonic image of the subject, and the pressure is applied based on the pressure. The stress information of the measurement site under conditions is calculated.
  • the model selection unit 182 selects a model (stress distribution template) from the database 183.
  • the model selection unit 182 may generate a model (stress distribution template) that matches the conditions of the measurement site by selecting a plurality of models and using linear interpolation or the like. For example, when a cylindrical model having a thickness of 25 [mm] is generated, the stress ⁇ (10 mm) of the cylindrical model having a thickness of 10 [mm] and a cylindrical shape are obtained as shown in FIG.
  • the model selection unit 182 may generate a model adapted to the conditions of the measurement site by linearly complementing the stress ⁇ (30 mm) of the model having a thickness of 30 [mm] with a weight of 1/4: 3/4. .
  • ⁇ z (25 mm) ⁇ z (10 mm) ⁇ 1/4 + ⁇ z (30 mm) ⁇ 3/4 (8)
  • the model selection unit 182 can generate a model (stress distribution template) that matches the conditions of the measurement site by combining a plurality of models.
  • the condition setting unit 180 sets the boundary line 50 (or boundary surface) as the condition of the measurement region in order to refer to the model (stereoscopic model).
  • the condition setting unit 180 may set a plurality of boundary lines (or boundary surfaces).
  • FIG. 9 is a diagram illustrating a state in which a plurality of boundary lines are set based on the tomographic image (B-mode image) 40. As shown in FIG. 9, in addition to the boundary line 50 at the boundary between the great pectoral muscle that is a soft tissue and the rib that is a hard tissue, the boundary line 51, at each boundary also in the soft tissue (fat, mammary gland, and great pectoral muscle) 52 is set.
  • a boundary line 53 may be set at the boundary of the specific part (target part).
  • a boundary line 50, 51, 52, and 53 at the respective boundaries, more accurate stress information or elastic information. Can be obtained. It is also possible to set boundary conditions for a plurality of boundary lines 50, 51, 52, and 53, respectively.
  • the model selection unit 182 selects a plurality of models from the database 183 based on the set boundary lines 50, 51, 52, and 53, and combines the plurality of models so that the model ( Stress distribution template) can be generated.
  • the elastic information calculation unit 13 calculates an elastic modulus (for example, Young's modulus E), but may calculate a strain ratio as elastic information. In this case, the elastic information calculation unit 13 generates strain ratio frame data SR (i, j) instead of the elastic frame data E (i, j) of the elastic modulus (Young's modulus E), and the display unit 10 The strain ratio image is displayed as an elastic image.
  • an elastic modulus for example, Young's modulus E
  • the elastic information calculation unit 13 generates strain ratio frame data SR (i, j) instead of the elastic frame data E (i, j) of the elastic modulus (Young's modulus E), and the display unit 10
  • the strain ratio image is displayed as an elastic image.
  • the coupler distortion in the acoustic coupler in the line is “ ⁇ 0 (i)”
  • the displacement / distortion calculation unit 12 Assuming that the measurement point distortion of the measurement site calculated by the above equation is “ ⁇ (i, j)”, the distortion ratio SR (i, j) in the same line is expressed by Expression (9).
  • the stress information calculation unit 184 calculates the stress attenuation rate ⁇ a using the stress distribution ⁇ T (i, j) of the model (stress distribution template).
  • the attenuation rate ⁇ a of the stress ⁇ T at the depth (j: coordinate in the z-axis direction) of each line (i: coordinate in the x-axis direction) is an expression for the stress ⁇ T (i, 0) at the depth 0 (zero). It is represented by (10).
  • ⁇ a (i, j) ⁇ T (i, 0) / ⁇ T (i, j) (10)
  • the elastic information calculation unit 13 calculates a corrected strain ratio SR ′ (i, j) from which the influence of stress attenuation is removed, based on the attenuation rate ⁇ a, as shown in Expression (11).
  • a value obtained by gradationizing and displaying the value of the corrected strain ratio SR '(i, j) is a strain ratio image (elastic image).
  • the elastic information processing unit 14 performs image processing on the elastic frame data having the corrected strain ratio SR ′, and the color scan converter 15 takes in the elastic frame data (elastic frame data indicating the elastic characteristics of the corrected strain ratio SR ′). Then, according to a color map corresponding to the distortion ratio SR ′, a color elasticity code is generated for each pixel to generate a color elastic image, and the display unit 10 displays the distortion ratio image.
  • a strain ratio image (elastic image) is generated based on the corrected strain ratio SR ′ from which the influence of stress attenuation has been removed, and the elasticity information of the measurement site (region of interest of the subject 1) is evaluated with high accuracy. It becomes possible.
  • the distortion of the acoustic coupler is obtained based on the distortion frame data obtained by the displacement / distortion calculation unit 12.
  • the distortion of the acoustic coupler may be obtained based on a change in the thickness of the acoustic coupler obtained from the coordinates of the boundary between the acoustic coupler and the subject 1.
  • the condition setting unit 180 may set conditions (at least one of a boundary and a boundary condition) of a measurement site based on the luminance of the ultrasonic image.
  • the soft tissue When the measurement site is soft tissue, the soft tissue is supported by hard tissue such as bone.
  • hard tissue such as bone.
  • the condition setting unit 180 sets the boundary between the soft tissue and the hard tissue as a boundary line (or boundary surface).
  • Ultrasonic waves have a characteristic that they are reflected more strongly at boundaries where the acoustic impedance changes greatly. Since the acoustic impedance of soft tissue and hard tissue (such as bone) is greatly different, the acoustic impedance changes greatly at the boundary between the soft tissue and the hard tissue. Accordingly, the boundary between the soft tissue and the hard tissue is displayed with high luminance on the tomographic image (B-mode image), and the echo is attenuated at a site deep in the boundary (inside the hard tissue), so that the inside of the hard tissue has low luminance. .
  • the condition setting unit 180 detects boundary coordinates between the soft tissue and the hard tissue based on the luminance (including luminance change) of the ultrasonic image, and sets the boundary between the soft tissue and the hard tissue as a boundary line (or boundary surface). For example, the condition setting unit 180 detects boundary coordinates of a portion (in the hard tissue) having a luminance lower than a predetermined threshold, and sets the boundary coordinates (or boundary surface) based on the boundary coordinates. In addition, the condition setting unit 180 sets the distance between the set boundary line (or boundary surface) and the contact surface 41 as the shape (thickness) of the measurement site.
  • the condition setting unit 180 may set conditions (shape, boundary conditions, etc.) of the measurement part based on at least one of displacement and distortion of the measurement part.
  • the condition setting unit 180 detects the boundary coordinates between the soft tissue and the hard tissue based on the displacement or distortion (including change of the displacement or distortion) of the ultrasonic image, and sets the boundary between the soft tissue and the hard tissue as a boundary line (or boundary surface). ).
  • the condition setting unit 180 detects boundary coordinates of a portion (inside hard tissue) having a displacement or distortion smaller than a predetermined threshold, and sets the boundary coordinates (or boundary surface) based on the boundary coordinates.
  • the condition setting unit 180 detects boundary coordinates of a portion (inside the hard tissue) having a strain change rate in the z-axis direction that is larger than a predetermined threshold, and uses the boundary coordinates (or boundary surface) as the boundary coordinates.
  • Set The condition setting unit 180 sets the distance between the set boundary line (or boundary surface) and the contact surface 41 as the shape (thickness) of the measurement site.
  • condition setting unit 180 may set the boundary condition based on at least one of displacement and distortion of the measurement site of the ultrasonic image.
  • the condition setting unit 180 may set the boundary condition of the measurement region based on at least one of displacement and distortion in the boundary direction (boundary line direction or boundary surface direction) of the measurement region.
  • the condition setting unit 180 performs “boundary is constrained in xyz direction (fixed)” or “boundary is constrained in x direction (fixed)”. If the boundary has a displacement or distortion in the boundary direction larger than a predetermined threshold value, “boundary is unconstrained in xyz direction (free)” or “boundary is unconstrained in x direction (free)” is set.
  • the x-axis direction is larger than the longitudinal strain ⁇ z (i, jb) in the z-axis direction.
  • the condition setting unit 180 determines that “the boundary is constrained (fixed) in the xyz direction” or “the boundary is x Set “Restrain (fixed) by direction”.
  • the condition setting unit 180 evaluates the displacement ratio D1 or the difference D2 to determine whether the plurality of places moved in the same direction or the plurality of places moved in the reverse direction as the interlinkage between the plurality of places across the boundary. It is possible to evaluate how much the difference in the movement of the parts is. That is, the condition setting unit 180 can evaluate the displacement near the boundary (for example, above and below the boundary or inside and outside the hard tissue) by evaluating the displacement ratio D1 or the difference D2.
  • the condition setting unit 180 determines that “the boundary is unconstrained in the xyz direction (free)” or “the boundary is Set “No constraint (free) in x direction”.
  • the condition setting unit 180 determines that “the boundary is constrained (fixed) in the xyz direction”. Or, “boundary is restricted in x direction (fixed)” is set.
  • the x-axis direction displacement Dx (a1 ⁇ i ⁇ a2, b1 ⁇ j ⁇ b2) at the point may be obtained, and the boundary condition may be set based on the statistical variation (dispersion, standard deviation, etc.) of the displacement. That is, the condition setting unit 180 may obtain a plurality of displacements across the boundary and set the boundary condition based on the displacements across the boundary.
  • the boundary condition straint state
  • condition setting unit 180 determines that “the boundary is in the xyz direction Is set to "No constraint (free)" or "Boundary is not constraint in the x direction (free)".
  • condition setting unit 180 Is “constrained (fixed) in xyz direction” or “boundary is constrained (fixed) in x direction”.
  • condition setting unit 180 sets the boundary condition of the measurement part based on at least one of the ratio, difference, and variation of the displacements at a plurality of positions sandwiching the boundary of the measurement part.
  • FIG. 11 is a diagram showing that the display unit 10 displays the selected model.
  • a model stress distribution template
  • FEM finite element method
  • the display unit 10 displays the model 60 selected from the database 183 or the model 60 used for the finite element method (FEM).
  • the display unit 10 may be constructed by designing the model 60 as an icon and displaying the model together with the tomographic image (B-mode image) or the elasticity image.
  • the model selection unit 182 may select a model from the database 183 based on a body mark representing the positional relationship between the subject 1 and the probe 2.
  • the body mark may be a three-dimensional solid body mark created from the volume data of the measurement site. Since the body mark schematically represents the measurement site, the type and shape of the measurement site can be grasped from the body mark.
  • the body mark is stored in the database in association with the ultrasonic image of the measurement site.
  • a body mark corresponding to the measurement site of the subject 1 is associated with the model and stored in the database. For example, a body mark for the mammary gland is associated with the model for the mammary gland and stored in the database in the breast, and a body mark for the thyroid gland is associated with the model for the thyroid gland and stored in the database.
  • the model selection unit 182 selects the shape of the model from the body mark corresponding to the measurement site, and selects the arrangement of the compression surface (such as the size, position, and direction of the compression surface) from the arrangement of the probe 2.
  • the model selection unit 182 may select a model from the database 183 based on volume data corresponding to the positional relationship between the subject 1 and the probe 2.
  • the position information of the probe 2 is detected in real time by a technique such as RVS (Real-time Virtual Sonography), and a CT image and an MRI image corresponding to the cross-sectional position of the tomographic image observed in the ultrasonic image Is constructed from volume data as a three-dimensional solid body mark.
  • the model selection unit 182 selects the shape of the model corresponding to the measurement site based on the body mark constructed from the volume data of the subject 1 corresponding to the cross-sectional position of the ultrasound image, and arranges the probe 2. Select the placement of the compression surface (such as the size, position and direction of the compression surface).
  • the condition setting unit 180 tracks (tracks) the boundary coordinates (boundary line or boundary surface coordinates) of the measurement site, and the stress information calculation unit 184 calculates real-time stress information based on the tracked boundary coordinates. Generate stress frame data in real time. As described above, the condition setting unit 180 sets a boundary between the soft tissue and the hard tissue at the measurement site. Since the elasticity information is obtained by compressing the subject 1, the distance from the contact surface where the probe 2 contacts the subject 1 to the boundary of the measurement site becomes shorter or longer depending on the amount of compression. And even during the measurement of stress information. The condition setting unit 180 detects boundary coordinates, tracks and follows the boundary coordinates at an arbitrary time, and the stress information calculation unit 184 constructs stress frame data corresponding to the boundary coordinates in real time.
  • the distance from the contact surface of the probe 2 to the boundary greatly changes due to increase or decrease in the amount of compression. It is effective to build the stress frame data in real time.
  • the stress information calculation unit 184 may calculate the stress information by calculating its own weight based on the model.
  • the stress information calculation unit 184 calculates the stress information by calculating the own weight (initial stress) of the measurement site based on the model having a thickness of 30 [mm].
  • the stress information calculation unit 184 can also calculate the stress information by calculating the own weight (initial stress) of the measurement site at an arbitrary depth. Thus, by considering the initial stress due to the weight of the measurement part in the calculation of the stress information at the measurement point, it is possible to measure a more accurate stress distribution.
  • FIG. 12 is a diagram showing that the display unit 10 displays elasticity information and stress information.
  • an ROI region of interest
  • the display unit 10 displays the elastic information and stress information in the ROI in the information display area 61.
  • the display unit 10 displays in the information display area 61 the relationship between the elastic information and the stress information in the ROI at each time (each frame) acquired by changing the strength of compression over time.
  • the display unit 10 displays the relationship between the elastic information and the stress in the information display area 61, the elastic characteristics of the tissue in the ROI can be analyzed in detail.
  • the ROI is set for the tomographic image (B-mode image) 40, but the ROI may be set for the elastic image.
  • the approximate function is not limited to a linear function, and may be any function that appropriately evaluates the characteristics of the plotted values.
  • FIG. 13 is a diagram showing that the display unit 10 displays sound speed information and stress information in the ROI.
  • dynamic elastography SSI, ARFI, etc.
  • the display unit 10 displays the elasticity information (including the shear wave velocity) and the stress information of the tissue in the ROI in the information display area 61.
  • the information display area 61 shows the relationship between the elastic information (including the shear wave velocity) in the ROI and the stress information at each time (each frame) acquired by the display unit 10 by changing the strength of the compression with time. Is displayed.
  • the display unit 10 displays the relationship between the elastic information and the stress in the information display area 61, the elastic characteristics of the tissue in the ROI can be analyzed in detail.
  • the condition setting unit 180 may set the boundary of the set ROI as the condition.
  • the present invention is capable of measuring stress distribution with high accuracy by setting the conditions (shape, boundary conditions, etc.) of a subject with high accuracy, and an ultrasonic imaging apparatus for displaying an ultrasonic image. It is useful as an ultrasonic image display method.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biomedical Technology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Acoustics & Sound (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

　本発明は、被検体の条件（形状や境界条件など）の設定を高精度で行うことにより、応力分布を高精度に計測することができる超音波画像撮像装置及び超音波画像表示方法を提供する。　本発明に係る超音波画像撮像装置は、被検体に印加された圧力を計測する圧力計測部と、前記被検体の超音波画像に基づいて計測部位に関する条件を設定する条件設定部と、前記圧力に基づいて前記条件における前記計測部位の応力情報を演算する応力情報演算部とを備える。

Description

超音波画像撮像装置及び超音波画像表示方法
　本発明は、超音波探触子及び超音波画像撮像装置に関し、特に、超音波画像を表示するための超音波画像撮像装置及び超音波画像表示方法に関する。
　超音波画像撮像装置は、被検体の表面に当接された超音波探触子から被検体に超音波を送信し、被検体の内部から超音波の反射波を受信し、その受信信号である反射エコー信号に基づいて被検体の部位の生体情報に関する超音波画像（断層画像など）を表示する。
　特許文献１，２に記載されている超音波診断装置は、被検体に圧迫力を与えた際に生じる被検体内部の各計測点における変位や歪みと各計測点における応力に基づく弾性情報から弾性画像を生成する（特許文献１及び特許文献２参照）。そして、検者は、圧迫力を調整しながら弾性画像を観察することにより、正常組織、ガン細胞、及び腫瘍などの病変を鑑別する。
　一方、組織の弾性情報を高精度に取得するためには、被検体に圧迫力を加えて内部に生じた応力（又は、応力分布）を高精度に計測することが重要となる。しかしながら、通常の超音波診断装置は応力（又は、応力分布）を直接計測することができないため、応力（又は、応力分布）の様々な推定手法が考案されている。例えば、有限要素法（ＦＥＭ）は、被検体の形状や拘束状態などの被検体の境界条件を仮定して応力分布を求めることにより、応力分布を推定する。
特開平５－３１７３１３号公報 特開２０００－６０８５３号公報
　しかしながら、特許文献１，２に記載された超音波診断装置は、被検体の条件（形状や境界条件など）を考慮していないので、被検体の条件（形状や境界条件など）に起因する被検体の内部の応力（又は、応力分布）の変化に対応できない。また、有限要素法（ＦＥＭ）では、被検体の形状や境界状態の設定が実際の形状や境界条件と異なる場合、応力分布の推定精度が低下し、結果的に信頼性の低い弾性画像が生成されることになる。
　本発明は、被検体の条件（形状や境界条件など）を考慮することにより、被検体の条件（形状や境界条件など）に起因する被検体の内部の応力（又は、応力分布）の変化に対応でき、被検体の条件（形状や境界条件など）の設定を高精度で行うことにより、応力分布を高精度に計測することができる超音波画像撮像装置を提供することを目的とする。
　本発明に係る超音波画像撮像装置は、被検体に印加された圧力を計測する圧力計測部と、前記被検体の超音波画像に基づいて計測部位に関する条件を設定する条件設定部と、前記圧力に基づいて前記条件における前記計測部位の応力情報を演算する応力情報演算部とを備える。
　本発明は、被検体の条件（形状や境界条件など）の設定を高精度で行うことにより、応力分布を高精度に計測することができる超音波画像撮像装置及び超音波画像表示方法を提供する。
本実施の形態に係る超音波画像撮像装置を示すブロック図である。 本実施の形態に係る応力演算部の構成の一例を示すブロック図である。 有限要素法（ＦＥＭ）による構造解析シミュレーションを用いることにより、応力σ又は応力分布σ（ｘ，ｙ，ｚ）を求めることを示す図である。 乳腺組織を含む乳房部位に圧迫を印加して生成された断層画像（Ｂモード画像）及び弾性画像を示した図である。 断層画像（Ｂモード画像）に基づいて、計測部位の条件が設定された状態を示した図である。 断層画像（Ｂモード画像）に基づいて、甲状腺（軟組織）と気管（固定組織）の間に境界線が設定された状態を示した図である。 データベースに格納されるモデルの一例を示した図である。 線形補完により計測部位の条件に適合したモデルを生成することを示した図である。 断層画像（Ｂモード画像）に基づいて、複数の境界線が設定された状態を示した図である。 カプラ歪みと計測点歪みに基づいて、同一ラインにおける歪み比を算出することを示した図である。 表示部が選択されたモデルを表示することを示した図である。 表示部が弾性情報と応力情報を表示することを示した図である。 表示部がＲＯＩ内の音速情報と応力情報を表示することを示した図である。
　以下、本発明の実施の形態に係る超音波画像撮像装置について、図面を用いて説明する。図１は、本実施の形態に係る超音波画像撮像装置を示すブロック図である。超音波画像撮像装置は、超音波を利用して被検体の診断部位について断層画像を得るとともに生体組織の硬さ又は軟らかさを表す弾性画像を表示するものである。
　図１に示すように、被検体１に当接して用いられる探触子（超音波探触子）２は、被検体１との間で超音波を送信及び受信する複数の振動子を有する。探触子２は、送信回路３から供給される超音波パルスにより駆動される。送受信制御回路４は、探触子２の複数の振動子を駆動する超音波パルスの送信タイミングを制御して、被検体１内に設定される焦点に向けて超音波ビームを形成する処理を行う。また、送受信制御回路４は、探触子２の振動子の配列方向に電子的に超音波ビームを走査する処理を行う。
　探触子２は、被検体１内から発生する反射エコー信号を受信して受信回路５に出力する。受信回路５は、送受信制御回路４から入力されるタイミング信号に従って、反射エコー信号を取り込んで増幅などの受信処理を行う。受信回路５により処理された反射エコー信号は整相加算回路６に入力される。整相加算回路６は、複数の振動子により受信された反射エコー信号の位相を合わせて加算することにより反射エコー信号を増幅する。整相加算回路６において整相加算された反射エコー信号は、信号処理部７に入力され、ゲイン補正、ログ圧縮、検波、輪郭強調、及びフィルタ処理などの信号処理が施される。
　信号処理部７により処理された反射エコー信号は、白黒スキャンコンバータ８に導かれて、超音波ビームの走査面に対応した２次元の断層画像データ（ディジタルデータ）に変換される。これらの信号処理部７と白黒スキャンコンバータ８によって断層画像（Ｂモード画像）の画像構成部が構成される。白黒スキャンコンバータ８から出力される断層画像データは、切替加算部９を介して表示部１０に供給される。表示部１０は、超音波画像である断層画像（Ｂモード画像）を表示する。
　また、整相加算回路６から出力される反射エコー信号は、フレームデータ取得部１１に導かれる。フレームデータ取得部１１は、超音波ビームの走査面（断層面）に対応する反射エコー信号をフレームデータとして取得し、複数のフレームデータをメモリなどの記憶部に格納する。変位・歪み演算部１２は、フレームデータ取得部１１に格納されている取得時刻が異なる複数のフレームデータを順次取り込み、取り込んだペアのフレームデータに基づいて断層面における計測点の変位又は歪みを求め、求められた変位又は歪みを変位フレームデータ又は歪みフレームデータ（変位・歪みフレームデータ）として弾性情報演算部１３に出力する。
　圧力計測部１７は、探触子２から被検体１の表面に印加された圧力を計測し、計測された圧力を圧力データとして応力演算部１８に出力する。ここで、探触子２の圧力を求める方法は適宜選択可能であり、圧力センサーや音響カプラなどを用いて被検体１の表面圧力を求める方法などがある。
　応力演算部１８は、圧力データに基づいて被検体１の内部の応力又は応力分布を応力フレームデータとして演算し、演算された応力フレームデータを弾性情報演算部１３に出力する。
　弾性情報演算部１３は、変位・歪みフレームデータと応力フレームデータに基づいて各計測点の生体組織の弾性率（例えば、ヤング率）などの弾性情報を弾性フレームデータとして求め、弾性フレームデータを弾性情報処理部１４に出力する。
　弾性情報処理部１４は、弾性情報演算部１３から入力される弾性フレームデータに対して、弾性フレームデータの座標平面内におけるスムージング処理やコントラスト最適化処理を施し、複数の弾性フレームデータ間における時間軸方向のスムージング処理などを施すことにより、様々な画像処理を行い、画像処理された弾性フレームデータをカラースキャンコンバータ１５に送出する。
　カラースキャンコンバータ１５は、弾性情報処理部１４により処理された弾性フレームデータ（弾性率などの弾性特性を示す弾性フレームデータ）を取り込み、弾性率に対応するカラーマップに従って、画素ごとに色調コードを付与してカラー弾性画像を生成する。
　カラースキャンコンバータ１５により生成されたカラー弾性画像（超音波画像）は、切替加算部９を介して表示部１０に表示される。また、切替加算部９は、白黒スキャンコンバータ８から出力される白黒の断層画像と、カラースキャンコンバータ１５から出力されるカラー弾性画像とを入力し、両画像を切り替えていずれか一方を表示部１０に表示させる機能と、両画像の一方を半透明にして加算合成して表示部１０に重畳表示させる機能と、両画像を並べて表示部１０に表示させる機能とを有する。また、切替加算部９から出力される画像データは、装置制御インターフェイス部１９の制御に従って、シネメモリ２０に格納される。シネメモリ２０に格納された画像データは、装置制御インターフェイス部１９の制御に従って、表示部１０に表示される。
　本実施形態の特徴に係る応力演算部１８は、被検体１の形状や被検体１の境界条件（拘束状態など）を取得し、それらの条件を考慮して各計測点に印加される応力を求める。
　図２は、本実施の形態に係る応力演算部１８の構成の一例を示すブロック図である。図２に示すように、応力演算部１８は、条件設定部１８０、条件調整部１８１、モデル選択部１８２、及び応力情報演算部１８４を備える。条件設定部１８０は、計測部位の条件を設定する。条件調整部１８１は、設定された条件を調整する。モデル選択部１８２は、計測部位における応力をシミュレーションするためのモデルをデータベース１８３から選択する。応力情報演算部１８４は、モデルに基づいて応力情報（応力又は応力分布）を演算し、応力フレームデータを生成する。
　次に、本実施形態の基本的な動作について説明する。まず、探触子２が被検体１を圧迫して、被検体１に超音波ビームを走査する。受信回路５が走査面からの反射エコー信号を連続的に受信する。整相加算回路６から出力される反射エコー信号に基づいて、信号処理部７及び白黒スキャンコンバータ８が断層画像を構成し、切替加算部９を介して表示部１０が断層画像を表示する。
　フレームデータ取得部１１は、反射エコー信号を取り込んでフレームレートに同期させてフレームデータを繰り返し取得し、内蔵された記憶部（フレームメモリなど）に時系列順に記憶する。そして、フレームデータ取得部１１は、取得時刻が異なるペアのフレームデータを１単位として、複数ペアのフレームデータを連続的に選択して変位・歪み演算部１２に出力する。変位・歪み演算部１２は、選択された１つのペアのフレームデータを１次元もしくは２次元相関処理し、走査面における複数の計測点の変位を計測して変位フレームデータ又は歪みフレームデータ（変位・歪みフレームデータ）を生成する。
　変位ベクトルの検出法として、例えば特開平５－３１７３１３号公報に記載されているブロックマッチング法又はグラジェント法が知られている。ブロックマッチング法では、画像を例えばＮ×Ｎ画素からなるブロックに分け、現フレーム中の着目されているブロック（例えば、計測点のブロック）に最も近似しているブロックを前フレームから探索し、これらのブロック位置に基づいて計測点の変位が求められ、変位フレームデータが生成される。また、ペアのＲＦ信号フレームデータの同一領域における自己相関を計算して変位が算出され、変位フレームデータが生成される。
　変位フレームデータに基づいて、各計測点の歪みが求められる。歪み演算は、例えば変位を空間微分することによって計算され、歪みフレームデータが生成される。
　圧力計測部１７は、探触子２の超音波送受信面と被検体１との間に設けられた圧力センサー１６や音響カプラにより、探触子２から被検体１の表面（圧迫面）に印加された圧力を計測して、計測された圧力を圧力データとして応力演算部１８に出力する。応力演算部１８は、圧力データに基づいて被検体１の内部の応力又は応力分布を応力フレームデータとして演算する。
　弾性情報演算部１３は、変位・歪みフレームデータ、及び応力フレームデータを用いて、計測点の弾性情報を求める。例えば、弾性情報演算部１３は、各計測点における応力と、弾性情報演算部１３で求められた歪みから走査面上の各計測点の弾性率（例えば、ヤング率Ｅ）を演算し、弾性フレームデータを弾性情報処理部１４に出力する。弾性情報処理部１４は、弾性フレームデータにスムージング処理などの処理を施して、カラースキャンコンバータ１５に出力する。カラースキャンコンバータ１５は、弾性情報に基づいてカラー弾性画像を生成する。例えば、カラー弾性画像では、弾性率に応じた色が２５６階調化による色調のグラデーションで画素ごとに割り当てられる。なお、カラースキャンコンバータ１５に代えて、白黒スキャンコンバータが用いられてもよい。この場合、白黒スキャンコンバータは弾性情報に基づいて白黒弾性画像を生成する。白黒弾性画像では、弾性率が大きい領域は輝度を明るくし、弾性率が小さい領域は輝度を暗くすることにより弾性情報を視覚的に認識でき、例えば正常組織、ガン細胞、及び腫瘍などの病変を鑑別することができる。
　上記のように、応力演算部１８は、圧力データに基づいて被検体１の内部の応力又は応力分布を応力フレームデータとして演算し、弾性情報演算部１３は、各計測点の弾性率（例えば、ヤング率Ｅ）などの弾性情報を算出する。応力演算部１８は、被検体１の条件（形状や境界条件など）である計測部位の形状や境界条件（拘束状態など）を仮定して応力分布を求め、弾性情報演算部１３は、各計測点の弾性率Ｅを演算する。例えば、計測部位を等方性弾性体と仮定し、つりあい方程式、歪み－変位関係式、及び応力－歪み関係式（フックの法則）の３つの式を用いると、式（１）に示すように、組織の各計測点の弾性率（例えば、ヤング率Ｅ）は、ｚ軸方向の歪みεｚと、ｘ，ｙ，ｚ軸方向の応力σｘ，σｙ，σｚから求められる。なお、応力分布は、空間（座標（ｘ，ｙ，ｚ））における応力σｘ，σｙ，σｚの分布である。
Ｅ＝（σｚ－ν（σｘ＋σｙ）)／εｚ     （ν：ポアソン比）      ・・・・・（１）
　式（２）に示すように、（σｚ－ν（σｘ＋σｙ）)をσと定義すると、弾性率（ヤング率Ｅ）は式（３）となる。
σ≡（σｚ－ν（σｘ＋σｙ）)   ・・・・・（２）
Ｅ＝σ／εｚ    ・・・・・（３）
　また、図３に示すように、有限要素法（ＦＥＭ）による構造解析シミュレーションを用いることにより、応力σ又は応力分布σ（ｘ，ｙ，ｚ）を求めることもできる。ここで、圧迫方向はｚ軸方向である。図３（ａ）に示すように、被検体１内部の応力分布を求めるために、計測部位の形状や境界条件（拘束状態など）が仮定される。
＜計測部位の条件＞
形状：円柱形（直径：φ１００［ｍｍ］，厚さ：１０、３０、及び５０[ｍｍ]）
硬さ（弾性率）：ヤング率Ｅ＝１．０[ｋＰａ]
ポアソン比：ν＝０．４９５
境界条件：圧迫面→拘束なし，底面→ｘｙｚ方向で拘束（固定），側面→拘束なし
　図３（ｂ）及び（ｃ）は、z軸方向に１[ｋＰａ]の圧力を圧迫面に印加ときの応力分布を解析した結果を示した図である。図３（ｂ）及び（ｃ）は、被検体１の厚さを１０[ｍｍ]、３０[ｍｍ]、５０[ｍｍ]と変化させたときの応力分布σｚ，σの深度における変化を示している。組織の弾性情報を高精度に取得するためには、被検体１に圧迫力を加えて内部に生じた応力分布σ，σzを高精度に計測することが重要となる。通常の超音波診断装置は応力（又は、応力分布）を直接計測することができないため、応力分布σ，σzは上記のような計測部位の条件を設定することにより推定されるが、推定される応力分布σ，σzは設定される条件（被検体１の形状や境界条件など）に応じて大きく変化する。つまり、応力分布σ，σzの推定結果は設定される条件に依存している。
　例えば、図３（ｂ）及び（ｃ）に示すように、応力分布σ，σzの推定結果は設定される計測部位の形状（厚さ）に依存している。また、底面の境界条件（拘束条件）である“ｘｙｚ方向で拘束（固定）”を“拘束なし”に変更すると、特に底面の近傍において応力分布σ，σzが大きく変化する。本実施の形態に係る超音波画像撮像装置は、計測部位の形状や境界条件を正確に設定することにより、応力又は応力分布を高精度に計測することができる。
　次に、本実施の形態に係る超音波画像撮像装置が計測部位の形状や境界条件を設定することについて説明する。図４は、乳腺組織を含む乳房部位に圧迫を印加して生成された断層画像（Ｂモード画像）及び弾性画像を示した図である。図４に示すように、探触子２が被検体１に接触している接触面４１から最深部の肋骨まで、脂肪、乳腺、腫瘍、大胸筋、及び肋骨が断層画像（Ｂモード画像）４０に撮像されている。本実施の形態では、生体組織の各部位における計測点の弾性率（ヤング率Ｅ）が弾性情報として求められる。
　条件設定部１８０が、計測部位の条件を設定する。図５は、断層画像（Ｂモード画像）４０に基づいて、計測部位の条件が設定された状態を示した図である。図５に示すように、条件設定部１８０は、表示部１０に表示された断層画像（Ｂモード画像）４０に基づいて、計測部位の条件（形状や境界条件など）を設定する。例えば、図５に示すように、圧迫により変形しやすい軟組織部（脂肪から大胸筋まで）の距離（厚さ）は約３０［ｍｍ］であり、圧迫により変形しにくい硬組織（肋骨）は軟組織を支持する底面（支持面）であり、ｘｙｚ方向で拘束（固定）されている。
　この場合、条件設定部１８０は、後述するモデル（立体モデル）を参照するために、境界（境界線又は境界面）、境界条件、計測部位の形状（大きさを含む）、及び圧迫位置などの計測部位の条件を設定する。例えば、断層画像（Ｂモード画像）４０に基づいて、軟組織である大胸筋と硬組織である肋骨の境界に境界線５０が設定され、接触面４１から境界線（肋骨）５０までの距離（約３０［ｍｍ］）が計測され、上面（接触面４１）と底面（境界線５０を通るｙ方向の面）と側面（底面と上面に挟まれたｚ方向又は圧迫方向の面）の境界条件として“圧迫面（接触面４１）→拘束なし，底面→ｘｙｚ方向で拘束（固定），側面（ｚ方向又は圧迫方向の面）→拘束なし”が設定される。ここで、軟組織と硬組織との間に境界線５０（又は、境界面）を設定したのは、硬組織は圧迫による変形が小さく、一般的に圧迫による移動が小さいからである。
　境界条件は、各面（上面、底面、及び側面）の特性に応じて設定される。例えば、底面の境界条件は、“底面→ｘｙ方向で拘束（固定），ｚ方向で拘束なし（自由）”や“底面→ｘｙ方向で拘束なし（自由），ｚ方向で拘束（固定）”など、底面の特性に応じて設定される。各面（上面、底面、及び側面）の特性に応じて、各面を規定する線や面の表示が変化することにより、各面の特性が視覚的に識別されてもよい。例えば、底面の特性に応じて、境界線５０（又は、境界面）の特徴（線又は点線の色、形状、及び太さなど）が変化し、例えば、境界条件が“底面→ｘｙｚ方向で拘束（固定）”である場合は境界線５０（又は、境界面）の色が赤く変化してもよい。また、各面の特性に応じて、境界条件を表す文字（例えば、“底面→ｘｙｚ方向で拘束（固定）”）や記号や図形などが付されてもよい。
　また、条件設定部１８０は、底面の代わりに、側面に境界線５０を設定することも可能である。図６は、断層画像（Ｂモード画像）４０に基づいて、甲状腺（軟組織）と気管（固定組織）の間に境界線５０が設定された状態を示した図である。ここで、軟組織と固定組織との間に境界線５０（又は、境界面）を設定したのは、固定組織は圧迫による移動が小さいからである。軟組織は、固定組織までは変形可能であり、固定組織を超えて変形しない。つまり、甲状腺（軟組織）と気管（固定組織）の間の境界線５０（又は、境界面）は、ｘ軸方向に固定されている。この場合、条件設定部１８０は、境界条件として“側面（ｚ方向又は圧迫方向の面）→ｘ軸方向で拘束（固定）”を設定する。
　また、図６に示すように、条件設定部１８０は、曲線を境界線５０として設定することも可能である。例えば、条件設定部１８０は、気管の曲面に沿って境界線５０を設定する。甲状腺（軟組織）と気管（固定組織）の間の境界線５０（又は、境界面）は、境界線５０の法線方向に固定されている。この場合、条件設定部１８０は、境界条件として“側面（曲面）→法線方向で拘束（固定）”を設定する。
　また、条件設定部１８０は、上面（接触面４１）に境界線５０を設定することも可能である。探触子２の形状がリニアかコンベックスかにより上面（接触面４１）の境界条件が変化する。また、探触子２の接触面の曲率により上面（接触面４１）の境界条件が変化する。したがって、上面（接触面４１）の境界条件（形状）などに応じて、境界線５０が設定される。
　境界線５０は、断層画像（Ｂモード画像）４０の代わりに、弾性画像で入力するようにすることも可能である。
　また、条件調整部１８１が設定された条件を調整してもよい。条件調整部１８１は、トラックボールやマウスなどの装置制御インターフェイス部１９からの入力に従って、超音波画像（断層画像又は弾性画像）上で境界線５０などの境界条件を調整（変更、移動、回転、拡大、及び縮小など）する。例えば、図５に示すように、条件調整部１８１は、水平に設定された境界線５０の角度や接触面４１からの距離を調整してもよい。また、図６に示すように、条件調整部１８１は、境界線５０が気管の曲面に沿うように境界線５０の形状を調整してもよい。
　モデル選択部１８２は、条件設定部１８０により設定され又は条件調整部１８１により調整された条件（境界条件など）に対応するモデルであって、計測部位における応力をシミュレーションするためのモデルをデータベース１８３から選択する。つまり、モデル選択部１８２は、設定された計測部位の条件に適合したモデル（応力分布テンプレート）をデータベース１８３から読み込む。
　データベース１８３は、計測部位の特徴に応じて複数のモデルを格納している。例えば、図７に示すように、データベース１８３は、様々な形状（円柱形や四角柱形など）、大きさ、支持体（硬組織や固定組織など）の位置（底面や側面や中央部など）、及び圧迫面の大きさや位置や方向などを組み合わせたモデルを格納する。計測部位が乳房である場合は、図７（ａ）に示す円柱形のモデルが適している。また、計測部位が足（大腿直筋など）や腕（回内筋など）である場合は、図７（ｃ）や（ｄ）のモデルが適している。モデルの大きさは、１［ｍｍ］単位で複数のモデルがデータベース１８３に格納されている。
　例えば、図５に示すように、条件設定部１８０により設定された計測部位の条件が“形状：円柱形（直径：φ１００［ｍｍ］，厚さ：３０[ｍｍ]）”、“硬さ（弾性率）：ヤング率Ｅ＝１．０[ｋＰａ]”、“ポアソン比：ν＝０．４９５”、及び“境界条件：圧迫面→拘束なし，底面→ｘｙｚ方向で拘束（固定），側面→拘束なし”である場合、データベース１８３は、設定された条件に適合する円柱形のモデル（例えば、図７（ａ）のモデル）が選択され、選択されたモデルの応力特性情報（例えば、応力分布σＴｘ（ｉ，ｊ），σＴｙ（ｉ，ｊ），σＴｚ（ｉ，ｊ）又は応力分布σＴ（ｉ，ｊ））をデータベース１８３から読み出す。応力分布σＴｘ（ｉ，ｊ），σＴｙ（ｉ，ｊ），σＴｚ（ｉ，ｊ）は、各モデルにおいて予めシミュレーションされたｘ，ｙ，ｚ軸方向の応力分布であり、座標（ｉ：ｘ軸方向の座標，ｊ：ｚ軸方向の座標）における応力分布である。応力分布σＴ（ｉ，ｊ）は式（４）で表される。
σＴ≡（σＴｚ－ν（σＴｘ＋σＴｙ）)   ・・・・・（４）
　応力演算部１８は、計測部位における応力をシミュレーションするためのモデルを用いて、有限要素法（ＦＥＭ）の構造解析シミュレーションにより、応力情報を演算してもよい。図３では、中央ラインにおける応力分布σを示したが、全ラインの応力分布を有限要素法（ＦＥＭ）でシミュレーションすることにより、座標（ｉ，ｊ）における応力分布σＴを算出することができる。また、有限要素法（ＦＥＭ）の代わりに、他の解析手法を用いることにより、モデル（応力分布テンプレート）の応力分布が予め算出され、データベース１８３に格納されてもよい。
　モデル（応力分布テンプレート）の応力分布σＴｘ，σＴｙ，σＴｚ，σＴは、圧迫面に均一に１［Ｐａ］の圧力が加えられた場合に生じる応力分布であり、規格化された応力分布である。したがって、圧迫面にＰ［Ｐａ］の圧力が加えられた場合の応力分布は、モデル（応力分布テンプレート）の応力分布σＴｘ，σＴｙ，σＴｚ，σＴをＰ倍することにより算出される。
　つまり、式（５）又は式（６）に示すように、応力情報演算部１８４がモデルに基づいて応力情報を演算し、応力フレームデータを生成する。
σｚ（ｉ，ｊ）＝Ｐ×σＴｚ（ｉ，ｊ）    ・・・・・（５）
σ（ｉ，ｊ）＝Ｐ×σＴ（ｉ，ｊ）        ・・・・・（６）
　また、圧迫面に均一に１［Ｐａ］の圧力が加えられた場合に生じる応力分布の代わりに、圧迫面における圧力分布Ｐ（ｉ，ｊ）［Ｐａ］に基づいて応力分布が予め算出され、データベース１８３に格納されてもよい。
　応力演算部１８は、圧力データＰに基づいて被検体１の内部の応力又は応力分布を応力フレームデータとして演算し、演算された応力フレームデータを弾性情報演算部１３に出力する。
　弾性情報演算部１３は、変位フレームデータＤ（ｉ，ｊ）、歪みフレームデータεｚ（ｉ，ｊ）、及び応力フレームデータσｚ（ｉ，ｊ），σ（ｉ，ｊ）を入力し、式（７）に示すように、弾性情報（ヤング率Ｅなど）の弾性フレームデータＥ（ｉ，ｊ）を求め、弾性フレームデータＥ（ｉ，ｊ）を弾性情報処理部１４に出力する。
Ｅ（ｉ，ｊ）＝σ（ｉ，ｊ）／εｚ（ｉ，ｊ）      ・・・・・（７）
　出力された弾性フレームデータＥ（ｉ，ｊ）は、弾性情報処理部１４、カラースキャンコンバータ１５、及び切替加算部９を介して、カラー弾性画像（超音波画像）として表示部１０に表示される。
　以上のように、本実施の形態に係る超音波画像撮像装置は、被検体の条件（形状や境界条件など）を考慮することにより、被検体の条件（形状や境界条件など）に起因する被検体の内部の応力（又は、応力分布）の変化に対応でき、被検体の条件（形状や境界条件など）の設定を高精度で行うことにより、応力分布を高精度に計測することができる。なお、本実施の形態に係る超音波画像表示方法は、被検体に印加された圧力を計測し、前記被検体の超音波画像に基づいて計測部位に関する条件を設定し、前記圧力に基づいて前記条件における前記計測部位の応力情報を演算する。
　以上、本発明にかかる実施の形態について説明したが、本発明はこれらに限定されるものではなく、請求項に記載された範囲内において変更・変形することが可能である。
　図７では、モデル選択部１８２は、モデル（応力分布テンプレート）をデータベース１８３から選択する。この場合、モデル選択部１８２は、複数のモデルを選択して、線形補完などを用いることにより、計測部位の条件に適合したモデル（応力分布テンプレート）を生成してもよい。例えば、円柱形の厚さ２５［ｍｍ］のモデルを生成する場合、図８及び式（８）に示すように、円柱形の厚さ１０［ｍｍ］のモデルの応力σ（１０ｍｍ）と円柱形の厚さ３０［ｍｍ］のモデルの応力σ（３０ｍｍ）を１／４：３／４の重み付けで線形補完して、モデル選択部１８２は計測部位の条件に適合したモデルを生成してもよい。
σｚ（２５ｍｍ）＝σｚ（１０ｍｍ）×１／４＋σｚ（３０ｍｍ）×３／４
　　　　　　　　　　　　　　　　　　　　　　　　　・・・・・（８）
　このように、モデル選択部１８２は、複数のモデルを組み合わせることにより、計測部位の条件に適合したモデル（応力分布テンプレート）を生成することができる。
　また、条件設定部１８０は、モデル（立体モデル）を参照するために、境界線５０（又は、境界面）を計測部位の条件として設定する。この場合、条件設定部１８０は、複数の境界線（又は、境界面）を設定してもよい。図９は、断層画像（Ｂモード画像）４０に基づいて、複数の境界線が設定された状態を示した図である。図９に示すように、軟組織である大胸筋と硬組織である肋骨の境界に境界線５０の他に、軟組織（脂肪、乳腺、及び大胸筋）においてもそれぞれの境界に境界線５１，５２が設定される。また、腫瘍などの特定部位（ターゲット部位）がある場合は、特定部位（ターゲット部位）の境界に境界線５３が設定されてもよい。軟組織であっても、脂肪、乳腺、及び大胸筋ではそれぞれ応力特性が異なるので、それぞれの境界に境界線５０，５１，５２，５３を設定することにより、より高精度な応力情報又は弾性情報を取得することができる。なお、複数の境界線５０，５１，５２，５３において境界条件をそれぞれ設定することも可能である。モデル選択部１８２は、設定された複数の境界線５０，５１，５２，５３に基づいて複数のモデルをデータベース１８３から選択し、複数のモデルを組み合わせることにより、計測部位の条件に適合したモデル（応力分布テンプレート）を生成することができる。
　弾性情報演算部１３は、弾性率（例えば、ヤング率Ｅ）を演算するが、歪み比を弾性情報として演算してもよい。この場合、弾性情報演算部１３は、弾性率（ヤング率Ｅ）の弾性フレームデータＥ（ｉ，ｊ）の代わりに、歪み比フレームデータＳＲ（ｉ，ｊ）を生成し、表示部１０は、歪み比画像を弾性画像として表示する。
　例えば、２次元のフレームで考えた場合、図１０に示すように、ライン（ｉ：ｘ軸方向の座標）における音響カプラ内のカプラ歪みを“ε０（ｉ）”とし、変位・歪み演算部１２により演算された計測部位の計測点歪みを“ε（ｉ，ｊ）”とすると、同一ラインにおける歪み比ＳＲ（ｉ，ｊ）は式（９）で表される。
ＳＲ（ｉ，ｊ）＝ε０（ｉ）／ε（ｉ，ｊ）        ・・・・・（９）
　ただし、応力の減衰の影響により、深度方向（ｚ軸方向）の深度が大きくなるにつれて、計測部位の歪み比ＳＲ（ｉ，ｊ）は小さくなる。
　そこで、応力情報演算部１８４は、モデル（応力分布テンプレート）の応力分布σＴ（ｉ，ｊ）を用いて、応力の減衰率σａを演算する。例えば、各ライン（ｉ：ｘ軸方向の座標）の深度（ｊ：ｚ軸方向の座標）における応力σＴの減衰率σａは、深度０（ゼロ）の応力σＴ（ｉ，０）に対して式（１０）で表される。
σａ（ｉ，ｊ）＝σＴ（ｉ，０）／σＴ（ｉ，ｊ）  ・・・・・（１０）
　弾性情報演算部１３は、式（１１）に示すように、減衰率σａに基づいて、応力の減衰の影響を除去した補正歪み比ＳＲ’（ｉ，ｊ）を算出する。
ＳＲ‘（ｉ，ｊ）＝（ε０（ｉ）／ε（ｉ，ｊ））／σａ（ｉ，ｊ）
　　　　　　　　　　　　　　　　　　　　　　　　・・・・・（１１）
　補正歪み比ＳＲ’（ｉ，ｊ）の値を階調化して表示したものが歪み比画像（弾性画像）となる。弾性情報処理部１４は、補正歪み比ＳＲ’の弾性フレームデータに対して画像処理を行い、カラースキャンコンバータ１５は、弾性フレームデータ（補正歪み比ＳＲ’の弾性特性を示す弾性フレームデータ）を取り込み、歪み比ＳＲ’に対応するカラーマップに従って、画素ごとに色調コードを付与してカラー弾性画像を生成し、表示部１０が歪み比画像を表示する。
　このように、応力の減衰の影響を除去した補正歪み比ＳＲ’に基づいて歪み比画像（弾性画像）が生成され、計測部位（被検体１の関心部位）の弾性情報を高精度に評価することが可能となる。
　なお、音響カプラの歪みは、変位・歪み演算部１２により求められた歪みフレームデータに基づいて求められる。また、音響カプラの歪みは、音響カプラと被検体１の境界の座標により求められた音響カプラの厚さの変化に基づいて求められてもよい。
　条件設定部１８０は、超音波画像の輝度に基づいて計測部位の条件（境界及び境界条件の少なくとも１つ）を設定してもよい。
　計測部位が軟組織である場合、軟組織は骨などの硬組織により支持されている。例えば、被検体１の対象が乳房（軟組織）であれば肋骨が乳房を支持し、被検体１の対象が大腿直筋（軟組織）であれば大腿骨が支持している。したがって、条件設定部１８０は、軟組織と硬組織の境界を境界線（又は、境界面）として設定する。
　超音波は、音響インピーダンスが大きく変化する境界でより強く反射するという特性を有する。軟組織と硬組織（骨など）は音響インピーダンスが大きく異なるため、軟組織と硬組織の境界では音響インピーダンスが大きく変化する。したがって、断層画像（Ｂモード画像）上で軟組織と硬組織の境界は高輝度で表示され、境界から深度が深い部位（硬組織内部）ではエコーが減弱するので、硬組織内部は低輝度となる。条件設定部１８０は、超音波画像の輝度（輝度変化を含む）に基づいて軟組織と硬組織の境界座標を検出し、軟組織と硬組織の境界を境界線（又は、境界面）として設定する。例えば、条件設定部１８０は、所定の閾値より低い輝度を有する部分（硬組織内部）の境界座標を検出し、境界座標に基づいて境界線（又は、境界面）として設定する。また、条件設定部１８０は、設定された境界線（又は、境界面）と接触面４１との距離を計測部位の形状（厚さ）として設定する。
　条件設定部１８０は、計測部位の変位及び歪みの少なくとも１つに基づいて計測部位の条件（形状や境界条件など）を設定してもよい。
　軟組織と硬組織は硬さが異なるため、軟組織と硬組織の境界では変位又は歪み（例えば、ｚ軸方向の歪みεｚ（ｉ，ｊ））が変化する。例えば、骨は軟組織と比較して顕著に硬いので、骨と軟組織の境界で歪みεｚ（ｉ，ｊ）の大きさが急激に減少する。条件設定部１８０は、超音波画像の変位又は歪み（変位又は歪みの変化を含む）に基づいて軟組織と硬組織の境界座標を検出し、軟組織と硬組織の境界を境界線（又は、境界面）として設定する。例えば、条件設定部１８０は、所定の閾値より小さい変位又は歪みを有する部分（硬組織内部）の境界座標を検出し、境界座標に基づいて境界線（又は、境界面）として設定する。また、条件設定部１８０は、所定の閾値より大きいｚ軸方向の歪みの変化率を有する部分（硬組織内部）の境界座標を検出し、境界座標に基づいて境界線（又は、境界面）として設定する。条件設定部１８０は、設定された境界線（又は、境界面）と接触面４１との距離を計測部位の形状（厚さ）として設定する。
　また、条件設定部１８０は、超音波画像の計測部位の変位及び歪みの少なくとも１つに基づいて境界条件を設定してもよい。条件設定部１８０は、計測部位の境界方向（境界線方向又は境界面方向）の変位及び歪みの少なくとも１つに基づいて計測部位の境界条件を設定してもよい。
　例えば、条件設定部１８０は、境界が所定の閾値より小さい境界方向の変位又は歪みを有する場合は、“境界がｘｙｚ方向で拘束（固定）”又は“境界がｘ方向で拘束（固定）”を設定し、境界が所定の閾値より大きい境界方向の変位又は歪みを有する場合は、“境界がｘｙｚ方向で拘束なし（自由）”又は“境界がｘ方向で拘束なし（自由）”を設定する。
　また、境界線の方向がｘ軸方向である場合、条件設定部１８０は、境界座標（ｊ＝ｊｂ）において、ｘ軸方向の歪みεｘ（ｉ，ｊｂ）を評価することにより、境界座標における境界条件（拘束状態）を設定する。この場合、条件設定部１８０は、式（１２）に示すように、ｚ軸方向の歪みεｚ（ｉ，ｊｂ）とｘ軸方向の歪みεｘ（ｉ，ｊｂ）の比率（ポアソン比ν）を求め、比率（ポアソン比ν）に基づいて境界条件を設定してもよい。
ポアソン比ν（ｉ，ｊｂ）＝－εｘ（ｉ，ｊｂ）／εｚ（ｉ，ｊｂ）
　　　　　　　　　　　　　　　　　　　　　　　　・・・・・（１２）
　境界座標（ｊ＝ｊｂ）において、ポアソン比ν（ｉ，ｊｂ）が十分に大きい場合（所定の閾値より大きい場合）は、ｚ軸方向の縦歪みεｚ（ｉ，ｊｂ）に比べてｘ軸方向の横歪み－εｘ（ｉ，ｊｂ）が十分に大きいことを意味し、条件設定部１８０は、“境界がｘｙｚ方向で拘束なし（自由）”又は“境界がｘ方向で拘束なし（自由）”を設定する。一方、ポアソン比ν（ｉ，ｊｂ）が０（ゼロ）に近い場合（所定の閾値より小さい場合）は、条件設定部１８０は、“境界がｘｙｚ方向で拘束（固定）”又は“境界がｘ方向で拘束（固定）”を設定する。
　また、境界線の方向がｘ軸方向である場合、条件設定部１８０は、境界座標（ｊ＝ｊｂ）において、ｘ軸方向の変位Ｄｘ（ｉ，ｊｂ）を評価することにより、境界座標における境界条件（拘束状態）を設定する。この場合、条件設定部１８０は、境界を挟む複数個所の変位を求め、境界を挟む複数個所の変位に基づいて境界条件を設定してもよい。例えば、式（１３）及び式（１４）に示すように、条件設定部１８０は、硬組織（骨など）と軟組織の境界座標（ｊ＝ｊｂ）の上下（硬組織の内外）の２個所（ｉ，ｊｂ－１）及び（ｉ，ｊｂ＋１）の計測点におけるｘ軸方向の変位Ｄｘ（ｉ，ｊｂ－１）及びＤｘ（ｉ,ｊｂ＋１）を求め、変位の比Ｄ１又は差Ｄ２に基づいて境界条件を設定してもよい。
Ｄ１＝Ｄｘ（ｉ，ｊｂ－１）／Ｄｘ（ｉ，ｊｂ＋１）・・・・・（１３）
Ｄ２＝Ｄｘ（ｉ，ｊｂ－１）－Ｄｘ（ｉ，ｊｂ＋１）・・・・・（１４）
　条件設定部１８０は、変位の比Ｄ１又は差Ｄ２を評価することにより、境界を挟む複数箇所の連動性として、複数箇所が同方向に動いたのか、複数箇所が逆方向に動いたのか、複数箇所の動きの大きさの違いはどのくらいなのかを評価することができる。つまり、条件設定部１８０は、変位の比Ｄ１又は差Ｄ２を評価することにより、境界の近傍（例えば、境界の上下又は硬組織の内外）でのズリを評価することができる。例えば、ｘ軸方向のズリが大きい場合（変位の比Ｄ１又は差Ｄ２が所定の閾値より大きい場合）は、条件設定部１８０は、“境界がｘｙｚ方向で拘束なし（自由）”又は“境界がｘ方向で拘束なし（自由）”を設定する。一方、ｘ軸方向のズリが０（ゼロ）に近い場合（変位の比Ｄ１又は差Ｄ２が所定の閾値より小さい場合）は、条件設定部１８０は、“境界がｘｙｚ方向で拘束（固定）”又は“境界がｘ方向で拘束（固定）”を設定する。
　また、条件設定部１８０は、硬組織（骨など）と軟組織の境界座標（ｊ＝ｊｂ）の上下（硬組織の内外）の複数個所（ａ１≦ｉ≦ａ２，ｂ１≦ｊ≦ｂ２）の計測点におけるｘ軸方向の変位Ｄｘ（ａ１≦ｉ≦ａ２，ｂ１≦ｊ≦ｂ２）を求め、変位の統計的バラツキ（分散や標準偏差など）に基づいて境界条件を設定してもよい。つまり、条件設定部１８０は、境界を挟む複数個所の変位を求め、境界を挟む複数個所の変位に基づいて境界条件を設定してもよい。
　条件設定部１８０は、硬組織（骨など）と軟組織の境界座標（ｊ＝ｊｂ）において、ｘ軸方向の変位Ｄｘ（ｉ，ｊｂ）の近傍（ａ１≦ｉ≦ａ２，ｂ１≦ｊｂ≦ｂ２）の局所的なバラツキを評価することにより、境界座標（ｊ＝ｊｂ）の近傍における境界条件（拘束状態）を識別することができる。例えば、条件設定部１８０は、硬組織と軟組織の境界座標（ｊ＝ｊｂ）の計測点において、境界座標（ｊ＝ｊｂ）を中心とした局所領域として、５×５の局所領域（ｉ－２≦ｉ≦ｉ＋２，ｊｂ－２≦ｊｂ≦ｊｂ＋２）を設定し、その２５個の計測点におけるｘ軸方向の変位Ｄｘ（ｉ－２≦ｉ≦ｉ＋２，ｊｂ－２≦ｊｂ≦ｊｂ＋２）の標準偏差を求める。ｘ軸方向の変位Ｄｘ（ｉ－２≦ｉ≦ｉ＋２，ｊｂ－２≦ｊｂ≦ｊｂ＋２）の標準偏差が大きい場合（所定の閾値より大きい場合）は、条件設定部１８０は、“境界がｘｙｚ方向で拘束なし（自由）”又は“境界がｘ方向で拘束なし（自由）”を設定する。一方、変位Ｄｘ（ｉ－２≦ｉ≦ｉ＋２，ｊｂ－２≦ｊｂ≦ｊｂ＋２）の標準偏差が０（ゼロ）に近い場合（所定の閾値より小さい場合）は、条件設定部１８０は、“境界がｘｙｚ方向で拘束（固定）”又は“境界がｘ方向で拘束（固定）”を設定する。
　このように、条件設定部１８０は、計測部位の境界を挟む複数個所の変位の比、差、及びバラツキの少なくとも１つに基づいて計測部位の境界条件を設定する。
　図１１は、表示部１０が選択されたモデルを表示することを示した図である。モデル（応力分布テンプレート）を選択する場合は、計測部位における応力をシミュレーションするためのモデルがデータベース１８３に予め格納されている。また、有限要素法（ＦＥＭ）を用いる場合にも、計測部位における応力をシミュレーションするためのモデルが必要となる。表示部１０は、データベース１８３から選択されたモデル６０又は有限要素法（ＦＥＭ）に用いられるモデル６０を表示する。表示部１０は、モデル６０をアイコンとして図案化して構築し、断層画像（Ｂモード画像）や弾性画像とともにモデルを表示してもよい。
　モデル選択部１８２は、被検体１と探触子２の位置関係を表すボディマークに基づいてモデルをデータベース１８３から選択してもよい。ボディマークは、計測部位のボリュームデータから作成された３次元立体ボディマークであってもよい。ボディマークは計測部位を模式的に表したものであるので、ボディマークから計測部位の種類や形状などを把握することができる。また、ボディマークは、計測部位の超音波画像に関連付けてデータベースに格納されている。被検体１の計測部位に対応するボディマークがモデルに関連付けられてデータベースに格納されている。例えば、乳房には乳腺用のボディマークが乳腺用のモデルに関連付けられてデータベースに格納されており、甲状腺には甲状腺用のボディマークが甲状腺用のモデルに関連付けられてデータベースに格納されている。探触子２により取得される断層画像の断面位置を記録するために、ボディマーク上に探触子２の配置（位置や角度など）が記録される。モデル選択部１８２は、計測部位に対応するボディマークからモデルの形状を選択し、探触子２の配置から圧迫面の配置（圧迫面の大きさや位置や方向など）を選択する。
　モデル選択部１８２は、被検体１と探触子２の位置関係に対応するボリュームデータに基づいてモデルをデータベース１８３から選択してもよい。近年、ＲＶＳ（Ｒｅａｌ－ｔｉｍｅ　Ｖｉｒｔｕａｌ　Ｓｏｎｏｇｒａｐｈｙ）などの技術により、探触子２の位置情報がリアルタイムで検出されて、超音波画像で観察されている断層画像の断面位置に対応するＣＴ画像やＭＲＩ画像が３次元立体ボディマークとしてボリュームデータから構築される。モデル選択部１８２は、超音波画像の断面位置に対応して被検体１のボリュームデータから構築されるボディマークに基づいて、計測部位に対応するモデルの形状を選択し、探触子２の配置から圧迫面の配置（圧迫面の大きさや位置や方向など）を選択する。
　条件設定部１８０は、計測部位の境界座標（境界線又は境界面の座標）を追跡（トラッキング）し、応力情報演算部１８４は、追跡されている境界座標に基づいてリアルタイムの応力情報を演算し、応力フレームデータをリアルタイムで生成する。上記のように、条件設定部１８０は、計測部位の軟組織と硬組織との間に境界を設定する。被検体１を圧迫することにより弾性情報が所得されるため、探触子２が被検体１に接触する接触面から計測部位の境界までの距離は、圧迫量に応じて短くなったり長くなったりと、応力情報の計測中でも変化する。条件設定部１８０は、境界座標を検出し、任意の時刻における境界座標をトラッキングして追従し、応力情報演算部１８４は、境界座標に応じた応力フレームデータをリアルタイムで構築する。例えば、乳房などの非常に軟らかい組織では、圧迫量の増減によって探触子２の接触面から境界（肋骨などの硬組織）までの距離が大きく変化するため、トラッキングされている境界座標に応じた応力フレームデータをリアルタイムで構築することは効果的である。
　応力情報演算部１８４は、モデルに基づいて計測部位の自重を演算して、応力情報を演算してもよい。計測部位の厚さが大きいほど、圧迫面に近い浅部と圧迫面から遠い深部の間では計測部位の自重による初期応力の影響が異なり、深部ほど計測部位の自重による初期応力の大きさが大きくなる。被検体１の比重を１ｇ／ｃｍ３とすると、例えば厚さ３０［ｍｍ］の計測部位の底面（ｘｚ面）には、ｚ軸方向の初期応力（σｚ０＝３ｇｆ／ｃｍ２）が印加されている。よって、応力情報演算部１８４は、厚さ３０［ｍｍ］のモデルに基づいて計測部位の自重（初期応力）を演算して、応力情報を演算する。また、応力情報演算部１８４は、任意の深度の計測部位の自重（初期応力）を演算して、応力情報を演算することもできる。このように、計測部位の自重による初期応力を計測点における応力情報の演算に考慮することにより、より高精度な応力分布を計測することができる。
　図１２は、表示部１０が弾性情報と応力情報を表示することを示した図である。図１２に示すように、断層画像４０の画像領域にＲＯＩ（関心領域）が設定され、表示部１０は、ＲＯＩ内の弾性情報及び応力情報を情報表示領域６１に表示する。図１２では、表示部１０が、圧迫の強さを経時的に変化させて取得した各時刻（各フレーム）におけるＲＯＩ内の弾性情報と応力情報の関係を情報表示領域６１に表示している。表示部１０が弾性情報と応力の関係を情報表示領域６１に表示することにより、ＲＯＩ内の組織の弾性特性を詳細に解析することが可能となる。図１２では、断層画像（Ｂモード画像）４０にＲＯＩが設定されているが、弾性画像にＲＯＩが設定されてもよい。
　情報表示領域６１にプロットされた値は、近似関数（例えば、１次関数ｙ＝ａｘ＋ｂ）で近似されて、表示部１０は、近似線を表す値（傾き：ａ及びｙ切片：ｂなど）を表示してもよい。また、近似関数は、１次関数に限られず、プロットされた値の特性を適切に評価する任意の関数であってもよい。
　図１３は、表示部１０がＲＯＩ内の音速情報と応力情報を表示することを示した図である。図１３に示すように、放射圧により対象組織を圧迫し、組織を伝搬する横波（ｓｈｅａｒ　ｗａｖｅ：せん断弾性波）の音速などを弾性情報として構築及び表示する動的エラストグラフィ（ＳＳＩやＡＲＦＩなど）を用いる場合に、表示部１０は、ＲＯＩ内の組織の弾性情報（横波音速を含む）及び応力情報を情報表示領域６１に表示する。図１３では、表示部１０が、圧迫の強さを経時的に変化させて取得した各時刻（各フレーム）におけるＲＯＩ内の弾性情報（横波音速を含む）と応力情報の関係を情報表示領域６１に表示している。表示部１０が弾性情報と応力の関係を情報表示領域６１に表示することにより、ＲＯＩ内の組織の弾性特性を詳細に解析することが可能となる。
　条件設定部１８０は、設定されたＲＯＩの境界を前記条件として設定してもよい。
　本発明は、被検体の条件（形状や境界条件など）の設定を高精度で行うことにより、応力分布を高精度に計測することができ、超音波画像を表示するための超音波画像撮像装置及び超音波画像表示方法などとして有用である。
２　探触子
３　送信回路
４　送受信制御回路
５　受信回路
６　整相加算回路
７　信号処理部
８　白黒スキャンコンバータ
９　切替加算部
１０　表示部
１１　フレームデータ取得部
１２　変位・歪み演算部
１３　弾性情報演算部
１４　弾性情報処理部
１５　カラースキャンコンバータ
１６　圧力センサー
１７　圧力計測部
１８　応力演算部
１９　装置制御インターフェイス部
２０　シネメモリ
６１　情報表示領域
１８０　条件設定部
１８１　条件調整部
１８２　モデル選択部
１８３　データベース
１８４　応力情報演算部

Claims (15)

  1. 　被検体に印加された圧力を計測する圧力計測部と、
    　前記被検体の超音波画像に基づいて計測部位に関する条件を設定する条件設定部と、
    　前記圧力に基づいて前記条件における前記計測部位の応力情報を演算する応力情報演算部と
    　を備えることを特徴とする超音波画像撮像装置。
  2. 　前記条件設定部は、前記計測部位の軟組織と硬組織との間に境界及び境界条件の少なくとも１つを設定することを特徴とする請求項１に記載の超音波画像撮像装置。
  3. 　前記条件設定部は、前記超音波画像の輝度に基づいて前記条件の境界及び境界条件の少なくとも１つを設定することを特徴とする請求項１に記載の超音波画像撮像装置。
  4. 　前記条件設定部は、前記計測部位の変位及び歪みの少なくとも１つに基づいて前記条件の境界及び境界条件の少なくとも１つを設定することを特徴とする請求項１に記載の超音波画像撮像装置。
  5. 　前記条件設定部は、前記計測部位の境界方向の変位及び歪みの少なくとも１つに基づいて前記計測部位の境界条件を設定することを特徴とする請求項１に記載の超音波画像撮像装置。
  6. 　前記条件設定部は、前記計測部位の境界を挟む複数個所の変位に基づいて前記計測部位の境界条件を設定することを特徴とする請求項１に記載の超音波画像撮像装置。
  7. 　前記条件設定部は、前記複数個所の変位の比、差、及びバラツキの少なくとも１つに基づいて前記計測部位の境界条件を設定することを特徴とする請求項６に記載の超音波画像撮像装置。
  8. 　前記計測部位における応力をシミュレーションするためのモデルであって、前記条件に対応するモデルをデータベースから選択するモデル選択部を備え、
    　前記応力情報演算部は、前記モデルに基づいて前記計測部位の応力情報を演算することを特徴とする請求項１に記載の超音波画像撮像装置。
  9. 　前記モデル選択部は、前記被検体と探触子の位置関係を表すボディマークに基づいて前記モデルを前記データベースから選択することを特徴とする請求項８に記載の超音波画像撮像装置。
  10. 　前記モデル選択部は、前記超音波画像の断面位置に対応して前記被検体のボリュームデータから構築される前記ボディマークに基づいて前記モデルを前記データベースから選択することを特徴とする請求項９に記載の超音波画像撮像装置。
  11. 　前記条件設定部は、前記計測部位の境界を追跡し、
    　前記応力情報演算部は、追跡されている境界に基づいてリアルタイムの応力情報を演算することを特徴とする請求項１に記載の超音波画像撮像装置。
  12. 　前記応力情報演算部は、前記モデルに基づいて前記計測部位の自重を演算して、前記応力情報を演算することを特徴とする請求項８に記載の超音波画像撮像装置。
  13. 　前記超音波画像を表示する表示部を備え、
    　前記表示部は、前記超音波画像に設定されたＲＯＩ内の弾性情報と応力情報を表示することを特徴とする請求項１に記載の超音波画像撮像装置。
  14. 　前記条件設定部は、前記超音波画像に設定されたＲＯＩの境界を前記条件として設定することを特徴とする請求項１に記載の超音波画像撮像装置。
  15. 　被検体に印加された圧力を計測し、
    　前記被検体の超音波画像に基づいて計測部位に関する条件を設定し、
    　前記圧力に基づいて前記条件における前記計測部位の応力情報を演算することを特徴とする超音波画像表示方法。
PCT/JP2014/072481 2013-08-29 2014-08-27 超音波画像撮像装置及び超音波画像表示方法 WO2015030076A1 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013177664A JP5890358B2 (ja) 2013-08-29 2013-08-29 超音波画像撮像装置及び超音波画像表示方法
JP2013-177664 2013-08-29

Publications (1)

Publication Number Publication Date
WO2015030076A1 true WO2015030076A1 (ja) 2015-03-05

Family

ID=52586626

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/072481 WO2015030076A1 (ja) 2013-08-29 2014-08-27 超音波画像撮像装置及び超音波画像表示方法

Country Status (2)

Country Link
JP (1) JP5890358B2 (ja)
WO (1) WO2015030076A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015043911A (ja) * 2013-08-29 2015-03-12 日立アロカメディカル株式会社 超音波画像撮像装置及び超音波画像表示方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6791624B2 (ja) * 2015-10-30 2020-11-25 ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー 超音波診断装置及びその制御プログラム
JP7078487B2 (ja) 2018-08-02 2022-05-31 富士フイルムヘルスケア株式会社 超音波診断装置及び超音波画像処理方法
JP7294996B2 (ja) * 2019-11-28 2023-06-20 富士フイルムヘルスケア株式会社 超音波診断装置及び表示方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005120358A1 (ja) * 2004-06-09 2005-12-22 Hitachi Medical Corporation 弾性画像表示方法及び超音波診断装置
JP4221555B2 (ja) * 2002-07-31 2009-02-12 毅 椎名 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法
JP4721893B2 (ja) * 2005-12-15 2011-07-13 パナソニック株式会社 超音波診断装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5890358B2 (ja) * 2013-08-29 2016-03-22 日立アロカメディカル株式会社 超音波画像撮像装置及び超音波画像表示方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4221555B2 (ja) * 2002-07-31 2009-02-12 毅 椎名 超音波診断システム、歪み分布表示方法及び弾性係数分布表示方法
WO2005120358A1 (ja) * 2004-06-09 2005-12-22 Hitachi Medical Corporation 弾性画像表示方法及び超音波診断装置
JP4721893B2 (ja) * 2005-12-15 2011-07-13 パナソニック株式会社 超音波診断装置

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015043911A (ja) * 2013-08-29 2015-03-12 日立アロカメディカル株式会社 超音波画像撮像装置及び超音波画像表示方法

Also Published As

Publication number Publication date
JP5890358B2 (ja) 2016-03-22
JP2015043911A (ja) 2015-03-12

Similar Documents

Publication Publication Date Title
JP4966578B2 (ja) 弾性画像生成方法及び超音波診断装置
US11715202B2 (en) Analyzing apparatus and analyzing method
JP5304986B2 (ja) 超音波診断装置
JP5485508B2 (ja) 改良された軟部組織の超音波歪測定のための方法及び装置
JP5689073B2 (ja) 超音波診断装置、及び３次元弾性比算出方法
JP5411699B2 (ja) 超音波撮像装置
JP4843432B2 (ja) 超音波診断装置
JP5726081B2 (ja) 超音波診断装置及び弾性画像の分類プログラム
JP7490840B2 (ja) 解析装置及び解析プログラム
JP5113322B2 (ja) 超音波診断装置
JP5890358B2 (ja) 超音波画像撮像装置及び超音波画像表示方法
JP2018511447A (ja) 組織形態及び弾性情報の処理方法、並びに弾性検出装置
JP2016112285A (ja) 超音波診断装置
JP2015136449A (ja) 超音波診断装置及びビームフォーミング方法
JP5075830B2 (ja) 超音波診断装置
US11250564B2 (en) Methods and systems for automatic measurement of strains and strain-ratio calculation for sonoelastography
WO2019075697A1 (zh) 一种超声弹性测量装置及弹性对比测量方法
JP5680703B2 (ja) 超音波診断装置
JP5638641B2 (ja) 超音波診断装置
JP5128149B2 (ja) 超音波診断装置
JP6230801B2 (ja) 超音波撮像装置及び超音波画像表示方法
JP2017192575A (ja) 超音波診断装置
JP2012055742A (ja) 超音波診断装置
JP2013183982A (ja) 超音波診断装置及び弾性画像生成方法
KR20140036977A (ko) 초음파 진단 장치 및 그 제어 프로그램

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14839010

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 14839010

Country of ref document: EP

Kind code of ref document: A1