WO2012105152A1 - 超音波診断装置,および方法 - Google Patents
超音波診断装置,および方法 Download PDFInfo
- Publication number
- WO2012105152A1 WO2012105152A1 PCT/JP2011/080448 JP2011080448W WO2012105152A1 WO 2012105152 A1 WO2012105152 A1 WO 2012105152A1 JP 2011080448 W JP2011080448 W JP 2011080448W WO 2012105152 A1 WO2012105152 A1 WO 2012105152A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- elastic modulus
- displacement
- region
- ultrasonic
- diagnostic apparatus
- Prior art date
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/42—Details of probe positioning or probe attachment to the patient
-
- 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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
-
- 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/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5238—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
- A61B8/5246—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image combining images from the same or different imaging techniques, e.g. color Doppler and B-mode
Definitions
- the present invention relates to an ultrasonic diagnostic apparatus for diagnosing a subject by ultrasonic transmission / reception, and more particularly to an ultrasonic diagnostic technique for detecting a difference in hardness inside the subject.
- a method of diagnosing the hardness (distortion, elastic modulus, etc.) inside the subject from ultrasound echo signals Technology.
- the worker presses the ultrasonic probe against the surface of the subject and compresses it to cause displacement in the tissue inside the living body (hereinafter referred to as the conventional method).
- the displacement in the compression direction is estimated from echo signals before and after compression of the living tissue by compression, and the distortion, which is the spatial differential amount of the displacement, is obtained and imaged.
- This method is extremely effective for organs (such as the mammary gland) that exist as an imaging target where compression from the body surface is easy, but may not always be effective for all imaging targets. .
- organs such as the mammary gland
- a slip surface exists as an intervening layer between the body surface and the liver, it is difficult to perform compression that causes sufficient displacement.
- a boundary condition is required.
- the boundary condition becomes complicated and it is difficult to estimate the elastic modulus. Therefore, there is a technique for diagnosing hardness by using a focused beam as a displacement-generating transmission beam, applying radiation pressure inside the subject, and suppressing the influence of the intervening layer to displace the target tissue.
- ARFI Acoustic Radiation Force Impulse Imaging
- the amount of tissue displacement generated in the direction in which the focused beam travels is imaged, or the shear wave is estimated from the estimation of the propagation speed of the shear wave generated in the direction perpendicular to the direction in which the focused beam travels with tissue displacement at the focal point.
- An elastic coefficient such as an average shear modulus in the propagation region (hereinafter referred to as elastic modulus) is measured and imaged.
- elastic modulus an average shear modulus in the propagation region
- the elastic modulus is measured at a plurality of positions in the imaging field, and the dense strain distribution (relative value of hardness) and the sparse elastic modulus (absolute value of hardness) obtained by the conventional method are used.
- the term “absolute value” includes a signed value, that is, a distance from 0, and an absolute value as a synonym for a relative value.
- an absolute value is mainly used as a synonym for a relative value. The word is used.
- the amount of attenuation of shear waves is larger than that of longitudinal waves.
- the area where the elastic modulus can be measured when a focused beam is irradiated to one focal point is limited to the range of shear wave propagation distance (about 5 to 10 mm). .
- the elastic modulus is measured in a region wider than this range, it is necessary to measure the elastic modulus by generating shear waves at a plurality of locations in the desired region.
- the irradiation time of the focused beam (several hundred ⁇ s to 1 ms) is several hundred to several thousand compared to the irradiation time of the ultrasonic beam used for the B-mode image. It is twice as long, and the temperature rise of biological tissue and ultrasound probe is proportional to the irradiation time. Therefore, the safety related to the temperature rise is lower than when measuring a B-mode image.
- the elastic modulus measurement time interval should be sufficiently large (1 to 2 seconds). If it is increased, the frame rate (number of screen updates per unit time) decreases.
- Patent Document 2 the position at which the elastic modulus is measured (the position where shear waves are generated).
- the selection method If a plurality of hardnesses are included in the elastic modulus measurement area per location, the measured elastic modulus becomes an average value, and the measurement accuracy decreases.
- An object of the present invention is to solve the above-mentioned problems, and to select an optimal part for measuring the elastic modulus of a subject, and then measure the elastic modulus, thereby making it possible to increase the accuracy of the elastic modulus measurement.
- a diagnostic device and an ultrasonic display method are provided.
- an ultrasonic diagnostic apparatus using an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal, the first displacement detection beam And a displacement calculator that calculates distortion information in the region 1 based on an echo signal from the received subject, and a displacement that radiates a focused beam into the subject to displace the tissue in the subject.
- the displacement of the shear wave generated by the focused beam is detected on the basis of the echo signal from the generated object and the received second displacement detection beam.
- an ultrasonic diagnostic apparatus having an elastic modulus calculation unit for detecting the elastic modulus of the image and a display unit for displaying a strain image based on the strain information and the elastic modulus.
- an ultrasonic diagnostic apparatus using an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal, the first displacement detection A distortion calculation unit that calculates distortion information in the region 1 based on the echo signal from the subject radiated and received, and a region 2 included in the region 1 is selected based on the distortion information.
- an ultrasonic diagnostic apparatus including an elastic modulus calculation unit for detecting a displacement of a shear wave generated by a focused beam and detecting an elastic modulus in the region 2 is provided.
- an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal from the subject is used, and an image is generated based on the received echo signal.
- the first displacement detection beam is emitted to receive the echo signal from the subject, and the distortion information in the first region is calculated and calculated.
- a distortion image based on the distortion information is displayed on the display unit, a focused beam is emitted into the subject to displace the tissue of the subject, and a second displacement detection beam is emitted to receive an echo signal from the subject.
- the displacement of the shear wave generated by the focused beam is detected, the elastic modulus in the second region included in the first region is calculated based on the displacement of the shear wave, and the calculated elastic modulus is displayed on the display unit.
- An ultrasonic display method is provided.
- the present invention it is possible to provide a high-accuracy ultrasonic diagnostic apparatus and method using a hybrid system that synthesizes and displays a strain image based on strain information and an elastic modulus obtained by shear wave generation. Further, by providing a single elastic modulus measurement site due to the generation of shear waves, a highly safe ultrasonic diagnostic apparatus and ultrasonic display method can be provided.
- FIG. 1 is a diagram illustrating an overall system configuration of an ultrasonic diagnostic apparatus according to Embodiment 1.
- FIG. It is a figure explaining the display method of the distortion image based on distortion information based on Example 1.
- FIG. It is a figure which shows the process flowchart of a hybrid system based on Example 1.
- FIG. It is a figure explaining the selection method of the elastic modulus measurement position based on Example 1.
- FIG. It is a figure explaining the beam forming of the ultrasonic wave based on Example 1.
- FIG. FIG. 3 is a diagram illustrating a measurement diagram using an ultrasonic probe according to the first embodiment. It is a figure explaining the estimation method of the propagation speed of a shear wave based on Example 1.
- FIG. 3 is a diagram illustrating an overall system configuration of an ultrasonic diagnostic apparatus according to a second embodiment.
- 6 is a diagram illustrating a measurement diagram using an ultrasonic probe according to Example 2.
- FIG. It is a figure explaining the sequence of the ultrasonic transmission / reception based on Example 2.
- FIG. It is a figure which shows the shear wave penetration based on Example 2.
- FIG. It is a figure explaining the displacement direction of a shear wave and the propagation direction of a shear wave based on Example 2.
- FIG. 10 is a diagram illustrating a processing flowchart of a hybrid system according to a second embodiment. It is a figure explaining the selected elastic modulus measurement position based on Example 4. FIG. It is a figure explaining the magnitude
- FIG. 1 is a diagram illustrating an overall configuration of an ultrasonic diagnostic system according to a first modification of the first embodiment.
- FIG. 1 is a diagram illustrating an overall configuration of an ultrasonic diagnostic system according to a second modification of the first embodiment.
- each function program executed by the processing unit of the computer is expressed by a function, means, unit, or the like.
- a program for calculating the elastic modulus is referred to as an elastic modulus calculating function, an elastic modulus calculating means, an elastic modulus calculating unit, and the like.
- a method of combining a strain image based on strain information and an elastic modulus measured by shear wave generation is referred to as a hybrid method.
- Embodiment 1 is an ultrasonic diagnostic apparatus using an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal, and radiates and receives a first displacement detection beam. Based on an echo signal from the sample, a distortion calculation unit 24 that calculates distortion information in the region 1, a displacement generation unit 10 that radiates a focused beam into the subject and displaces tissue in the subject, and a second The displacement detection beam is emitted, and the displacement of the shear wave generated by the focused beam is detected based on the echo signal from the received object, and the elastic modulus in the region 2 included in the region 1 is detected.
- the present invention relates to an ultrasonic diagnostic apparatus having an elastic modulus calculation unit 34 and a display unit 7 for displaying a strain image based on strain information and an elastic modulus.
- the present embodiment is the above-described ultrasonic diagnostic apparatus, which includes a measurement position selection unit 40 that selects at least one elastic modulus detection position that detects an elastic modulus based on strain information, and irradiates a focused beam.
- An ultrasonic diagnostic apparatus in which at least one focal position to be determined is determined from at least one elastic modulus detection position selected by the measurement position selection unit 40, and the measurement position selection unit 40 operates
- An ultrasonic diagnostic apparatus that enables a person to select an elastic modulus detection position while viewing an image displayed on the display unit 7 or the ultrasonic diagnostic apparatus described above, wherein the display unit 7 displays a distortion image to be displayed.
- the present invention relates to an ultrasonic diagnostic apparatus that displays a color scale indicating the elastic modulus in the display range.
- the present embodiment is an ultrasonic diagnostic apparatus using an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal, and radiates a first displacement detection beam to receive the wave.
- a distortion calculation unit 24 that calculates distortion information in the region 1 based on the echo signal from the subject
- a measurement position selection unit 40 that selects the region 2 included in the region 1 based on the distortion information
- a focused beam is emitted into the subject to displace the tissue in the subject, and a second displacement detection beam is emitted. Based on the echo signal from the received subject, the focused beam is emitted.
- the present invention relates to an ultrasonic diagnostic apparatus including an elastic modulus calculation unit 34 that detects a displacement of a shear wave caused by the above and detects an elastic modulus in the region 2.
- this embodiment is the above-described ultrasonic diagnostic apparatus, and the measurement position selection unit 40 selects the standard deviation of the strain distribution in the region 1 or the maximum value and the minimum value when selecting the region 2.
- the elastic modulus calculator 34 uses the elastic modulus in the region 2 and the strain information in the region 2.
- the measurement position selection unit 40 is an ultrasonic diagnostic apparatus that selects a region 2 by extracting a contour of a distortion distribution in the region 1 by image processing when selecting the region 2, or the ultrasonic diagnostic device described above.
- the ultrasonic diagnostic apparatus further includes a display unit 7 that displays a distortion image based on information and an elastic modulus, and the measurement position selection unit 40 can select the region 2 based on the distortion image displayed on the display unit 7 by the operator. It is about.
- this embodiment uses an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal from the subject, and displays an image on a display unit based on the received echo signal.
- a sound wave display method that emits a first displacement detection beam, receives an echo signal from a subject, calculates distortion information in the first region, and displays a distortion image based on the calculated distortion information.
- a focused beam is emitted into the subject to displace the tissue of the subject, a second displacement detection beam is emitted to receive an echo signal from the subject, and a shearing caused by the focused beam
- the present invention relates to an ultrasonic display method for detecting a wave displacement, calculating an elastic modulus in a second region included in the first region based on the shear wave displacement, and displaying the calculated elastic modulus on a display unit 7. .
- this embodiment relates to the above-described ultrasonic diagnostic method, which relates to the ultrasonic diagnostic method for determining the focal position of the focused beam selected based on the distortion information. Further, the present embodiment is the above-described ultrasonic diagnostic method, and when selecting the second region for emitting the second displacement detection beam, the distortion information is obtained based on the distortion image displayed on the display unit.
- the present invention relates to an ultrasonic diagnostic method for selecting from uniform locations. Furthermore, the present embodiment relates to the ultrasonic diagnostic method described above, which displays a scale indicating the elastic modulus in the display range of the strain image displayed on the display unit.
- FIG. 1 shows a specific example of the overall system configuration of the ultrasonic diagnostic apparatus according to the first embodiment.
- An ultrasonic probe 1 used in contact with a skin of a subject is formed having an ultrasonic transmission / reception surface on which a plurality of transducers for transmitting and receiving ultrasonic waves to and from the subject are arranged.
- the transmission / reception selector switch 2 is connected.
- the central control unit 3 is a control unit for controlling the ultrasonic diagnostic apparatus.
- the displacement generation unit 10 that causes displacement in the subject
- the transmission / reception changeover switch 2 the first ultrasonic transmission / reception unit 20, the strain calculation unit 24, the second ultrasonic transmission / reception unit 30, and the elastic modulus that is the elastic modulus are calculated.
- the elastic modulus calculation unit 34, the measurement position selection unit 40, and the color scale setting unit 50 are configured to be controlled by the central control unit 3.
- the ultrasonic probe 1 is connected to a displacement-generating transmission beam generation unit 13, a first ultrasonic transmission / reception unit 20, and a second ultrasonic transmission / reception unit 30 via a transmission / reception changeover switch 2.
- the transmission / reception wave changeover switch 2 is connected to the ultrasonic probe 1, the displacement generation transmission beam generation unit 13, the first ultrasonic transmission / reception unit 20, and the second ultrasonic transmission / reception unit 30 via the central control unit 3. It is controlled to disconnect or disconnect.
- the displacement-generating transmission beam generator 13 generates a focused beam that is radiated into the subject to displace the tissue in the subject.
- the first ultrasonic transmission / reception unit 20 gives a delay time and a weight to the transmission signal of the element of the ultrasonic probe 1 using the waveform generated by the displacement detection transmission waveform generation unit (not shown).
- the central control unit 3 controls the ultrasonic beam for detecting displacement to be focused on a desired position of the subject (not shown).
- the echo signal reflected in the subject and returned to the probe is converted into an electrical signal by the ultrasonic probe 1 and sent to the first ultrasonic transmission / reception unit 20.
- the first ultrasonic transmission / reception unit 20 includes a signal processing circuit that performs phasing addition of echo signals and performs envelope detection, log compression, bandpass filter, gain control, and the like.
- the output signal from the first ultrasonic transmission / reception unit 20 is input to the monochrome DSC (digital scan converter) 5 and the displacement calculation unit 22.
- the monochrome DSC 5 tomographic image (B-mode image) information representing the luminance composed of monochrome is formed.
- the displacement calculation unit 22 calculates the displacement of each part by image correlation using adjacent tomographic images of two frames.
- the displacement information output from the displacement calculation unit 22 is input to the distortion calculation unit 24, and the distortion of each part is calculated based on the spatial differentiation of the displacement. This distortion information is input to the color DSC 4 and is subjected to hue modulation according to the value of the distortion information.
- An image 41 (hereinafter referred to as a distortion image) 41 hue-modulated by the color DSC 4 based on the distortion information is transmitted to the synthesizing unit 6 and superimposed on the B-mode image 45 on the display unit 7 as a distortion image as shown in FIG. It is displayed.
- the distortion information calculated by the distortion calculation unit 24 is transmitted to the color scale setting unit 50 via the measurement position selection unit 40 and the central control unit 3.
- the color scale setting unit 50 creates a color scale 43 corresponding to the distorted image 41 based on the maximum and minimum distortion values.
- the color scale 43 is displayed on the display unit 7 adjacent to the B-mode image 45 and the distorted image 41.
- hard (hard) indicating that the distortion is small and soft (soft) indicating that the distortion is large are displayed.
- the distorted image 41 is an image that is color-tone-modulated based on the distortion information, but may be an image that is monochrome-modulated based on the distortion information.
- the distortion information calculated by the distortion calculation unit 24 is transmitted to the measurement position selection unit 40.
- the measurement position selection unit 40 performs signal processing based on the two-dimensional strain information and determines the position in the tomographic plane where the elastic modulus is measured.
- the position information determined by the measurement position selection unit 40 is input to the focal position setting unit 12 of the displacement generation unit 10 via the central control unit 3.
- the displacement generator 10 will be described.
- the displacement-generating transmission beam generating unit 13 for generating the focused beam uses the waveform generated by the displacement-generating transmission waveform generating unit 11 to transmit the transmission signals of the respective elements 100 of the ultrasonic probe 1.
- the delay time and weight are given to the central position so that the ultrasonic beam is focused at the focus position determined based on the position set by the focus position setting unit 12, that is, the position determined by the measurement position selection unit 40. It is controlled by the control unit 3.
- the electrical signal from the displacement-generating transmission beam generating unit 13 is converted into an ultrasonic signal by the ultrasonic probe 1 via the transmission / reception change-over switch 2, and is directed toward the subject (not shown) for generating the displacement. Is irradiated.
- the second ultrasonic transmission / reception unit 30 gives a delay time and a weight to the transmission signal of the element of the ultrasonic probe 1 using the waveform generated by the displacement detection transmission waveform generation unit (not shown), Control is performed by the central control unit 3 so that the ultrasonic beam for displacement detection is focused on a desired position of the subject (not shown).
- the echo signal reflected in the subject and returned to the probe is converted into an electrical signal by the ultrasonic probe 1 and sent to the second ultrasonic transmission / reception unit 30.
- the second ultrasonic transmission / reception unit 30 includes a signal processing circuit that performs phasing addition of echo signals and performs envelope detection, log compression, bandpass filter, gain control, and the like.
- the output signal from the second ultrasonic transmission / reception unit 30 is input to the shear wave displacement calculation unit 32.
- the shear wave displacement calculation unit 32 calculates the displacement of each part by correlation calculation.
- the displacement information output from the shear wave displacement calculation unit 32 is input to the elastic modulus calculation unit 34, and values relating to the hardness such as the shear wave propagation velocity and the elastic modulus are calculated.
- the value relating to the hardness is transmitted to the synthesizing unit 6 and is displayed on the same screen as the B-mode image and the distortion image displayed on the display unit 7.
- the central control unit 3 the displacement calculation unit 22, the strain calculation unit 24, the shear wave displacement calculation unit 32, the elastic modulus calculation unit 34, the measurement position selection unit 40, etc., which are part of the block shown in FIG. This can be realized by executing a program in a central processing unit (CPU).
- CPU central processing unit
- the first ultrasonic transmission / reception unit 20 and the second ultrasonic transmission / reception unit 30 have been described separately as shown in FIG. The explanation of the function is given priority. From the viewpoint of mounting the apparatus, the first ultrasonic transmission / reception unit 20 and the second ultrasonic transmission / reception unit 30 may be configured by one ultrasonic transmission / reception unit.
- step S00 the diagnosis of hardness by the hybrid method is started.
- the start signal is input via an input device (not shown).
- a B-mode image or a distorted image is displayed on the display unit 7.
- distortion information in a desired region ROI (Region of Interest) input by an input device (not shown) is acquired (step S02).
- the ROI for acquiring distortion information is called ROI_s or the first region (region 1).
- ⁇ Information on this distortion is acquired in the usual way. That is, the ultrasonic probe 1 is pressed against the body surface of the subject, and ultrasonic waves are repeatedly transmitted and received while pressing the body surface to detect displacement and distortion.
- processing such as delay time and weighting of ultrasonic transmission and conversion from an echo signal to an electric signal is performed by the first ultrasonic transmitting / receiving unit 20, displacement calculation is performed by a displacement calculating unit 22, and distortion information is calculated. Is performed by the distortion calculation unit 24.
- Step S04 the B-mode image from the monochrome DSC 5 and the distortion image from the color DSC 4 are combined by the combining unit 6 and displayed on the display unit 7 as shown in FIG. (Step S04).
- the measurement position selection unit 40 determines an optimal position for performing the elastic modulus measurement (step S06).
- the region where the elastic modulus is measured is ROI_e or the second region (region 2), and ROI_e is included in the region ROI_s or the first region (region 1) where the strain is displayed. Is a smaller area.
- An optimal elastic modulus measurement region is, for example, ROI_e (1) in FIG. 2, and is a region having the same strain value in the ROI_e, that is, the same hardness.
- a region unfavorable for elastic modulus measurement is a case where two or more strain amounts are mixed in ROI_e as in ROI_e (2).
- the elastic modulus is measured in an ROI in which two or more strains are mixed, the averaged elastic modulus is measured, so that the measurement accuracy is lowered. For this reason, it is necessary to set a region where the strain amount is the same as much as possible (that is, the strain amount is uniform) as the elastic modulus measurement position.
- the optimal elastic modulus measurement position is determined by applying signal processing to strain information and strain images.
- the hybrid type signal processing in this embodiment will be described.
- a case will be described in which the linear array type ultrasonic probe 1 is brought into contact with the body surface of the subject and the transmission beam for generating displacement is focused on the target tomographic plane in the body. Further, a case will be described in which the propagation direction of the displacement-generating transmission beam is perpendicular to the body surface in the desired tomographic plane.
- a rectangular area (kernel K) for performing a signal processing method is created.
- the size of the kernel K is preferably the same as the size of the region ROI_e where the elastic modulus is measured, that is, the second region (region 2).
- the width in the depth direction is a value represented by -6 dB width of the beam shape in the depth direction.
- the length in the azimuth direction of ROI_e is determined by the propagation distance in the shear wave propagation direction.
- the propagation distance of the shear wave refers to the maximum distance from the focal position at which the displacement of the shear wave can be detected in the direction along the propagation direction of the shear wave.
- the length of ROI_e in the depth direction is 10 mm
- the length of ROI_e in the azimuth direction is 5 mm.
- the length ly of the kernel K in the depth direction is desirably 10 mm
- the length lx of the ROI_e in the azimuth direction is desirably 5 mm.
- the length of each side of ROI_e and kernel K is not limited to 10 mm in the depth direction and 5 mm in the azimuth direction.
- the size of the kernel K is read from a storage medium (not shown) and used for signal processing in the measurement position selection unit 40.
- the size of the kernel K may be input by an operator via an input medium (not shown).
- the kernel K42 moves up and down and right and left in the distorted image 41.
- the position of the kernel K42 is assumed to be P (x, y).
- P (x, y) is the center position of the kernel K42
- the x-axis is the azimuth direction
- the y-axis is the depth direction.
- the optimal elastic modulus measurement position is determined from the value calculated using the strain information in the kernel. For example, the standard deviation S (x, y) of the distortion in the kernel K is calculated for the kernel position P (x, y).
- the optimum elastic modulus measurement position is determined as the position where the standard deviation is the smallest.
- Max-min the difference between the maximum value and the minimum value of the strain amount in the kernel, Max-min (x, y), is calculated, and the optimum elastic modulus measurement position is determined as the position where Max-min is the smallest. Is done.
- a region having the same distortion amount may be determined using a known calculation formula.
- image processing is performed on a distorted image such as the distorted image 41 using a well-known two-dimensional filter (hereinafter referred to as filter G) 44 used in image processing.
- filter G a well-known two-dimensional filter
- the size of the filter G44 is preferably the same as the size of ROI_e for measuring the elastic modulus.
- the length ly in the depth direction is 10 mm
- the length lx in the azimuth direction is 5 mm.
- the length of each side of the filter G is not limited to 10 mm in the depth direction and 5 mm in the azimuth direction.
- the size of an appropriate filter G with respect to a part such as a mammary gland, a prostate, a blood vessel, and a liver, a carrier frequency of an ultrasonic beam for generating a shear wave, an F value, and a focal length is determined from a storage medium (not shown). It is read out and used for signal processing in the measurement position selection unit 40.
- the size of the filter G may be input by an operator via an input medium (not shown).
- the filter G44 is a well-known filter in the image processing technology, for example, and a Laplacian filter that is a filter for extracting the contour of the image is used. Distortion in the distorted image 41 is extracted by image processing using a Laplacian filter. When the outline of the distortion distribution is extracted, the distortion information in the distortion image 41 is divided into a plurality of areas. As a result of the image processing, for example, as shown in FIG. 4B, a plurality of regions R (1), region R (2), region R ( 3) is obtained.
- the number of regions is three, but the number of regions is an arbitrary integer of 2 or more.
- the region R (1) and the region R (2) are suitable elastic modulus measurement regions, and the region R (3) is not an optimal elastic modulus measurement region.
- the optimum region R (n) for elastic modulus measurement is output to the central control unit 3.
- the optimum region R (n) is automatically selected by the measurement position selection unit 40 from the appropriate regions having the largest area.
- the filter G a well-known filter can be applied in addition to a Laplacian filter.
- step S08 in FIG. 3 the focal position determined based on the optimal elastic modulus measurement position and region is irradiated with a focused beam, which is a displacement-generating transmission beam, to generate a shear wave, and the shear wave propagation speed.
- the elastic modulus is calculated by estimating.
- the optimum elastic modulus measurement position determined by the measurement position selection unit 40 is output to the central control unit 3, and further, the focal point F of the focal point F is transmitted via the central control unit 3.
- the position is input to the focal position setting unit 12.
- the optimal elastic modulus measurement position is determined as P (x1, y1) in FIG. 4A, for example, the position of the focal point F is F (x1- ⁇ x, y1). ).
- ⁇ x is a value corresponding to half of the width in the azimuth direction of ROI_e, that is, a value corresponding to half of the propagation distance of the shear wave.
- the displacement generating transmission beam is irradiated to the left end of the kernel K, and the displacement of the shear wave propagating in the + x direction is detected.
- the displacement direction of the shear wave is a direction perpendicular to the x axis, that is, the y direction.
- the central control unit 3 When converting to the position of the focal point F with respect to the optimum elastic modulus measurement position, first, the central control unit 3 performs the frequency and F value of the carrier of the ultrasonic beam for generating shear waves, parts such as the breast, prostate, blood vessels, and liver.
- the optimum conversion method based on the above is read out from a recording medium or the like (not shown), and then the conversion operation is performed in the central control unit 3 or the central processing unit (not shown) that functions as the processing unit. ing.
- the beam generation of the displacement generation ultrasonic wave is performed by the displacement generation transmission beam generation unit 13.
- the distance between the focal point and the position of each element 100 of the ultrasonic probe 1 is obtained, and the distance difference between the elements is divided by the sound velocity of the object. This is realized by transmitting the delay time calculated for each element.
- the position of the focal point F can be changed by controlling the delay time.
- the generation of this radiation pressure generates a shear wave in a direction parallel to the subject surface starting from the focal point.
- the raster to be used for detecting the displacement of the shear wave propagation severe ⁇ m to several tens ⁇ m
- the sampling point on the raster are determined.
- the displacement detection beam reception PRF Pulse Repetition : Frequency
- the PRF is set so as to satisfy the Nyquist theorem for the expected shear wave frequency. For example, when the raster is in the same direction as the shear wave displacement, the PRF is set to be twice or more the shear wave frequency.
- the processing such as the delay time and weighting of the ultrasonic transmission and the conversion from the echo signal to the electric signal is performed by the second ultrasonic transmission / reception unit 30, the displacement calculation shear wave. This is performed by the displacement calculator 32.
- the ultrasonic transmission / reception for detecting the displacement of the shear wave is performed by the raster used for the displacement detection described above.
- ultrasonic wave transmission / reception is performed once before irradiation of a displacement-generating transmission beam, and a reference signal for displacement calculation is acquired.
- a plurality of ultrasonic signals are acquired by performing ultrasonic transmission / reception a plurality of times during the time from when the displacement generating transmission beam is irradiated until the shear wave finishes propagating in the ROI_e.
- a correlation calculation is performed between the plurality of ultrasonic signals after irradiation of the displacement-generating transmission beam and the reference signal, and the displacement is calculated.
- the displacement information calculated at a plurality of times is converted into the time waveform information of the shear wave displacement.
- the time waveform for each position x (n) there is a time t (n) at which the displacement of the shear wave becomes maximum. From the relational expression between x (n) and t (n), the shear wave propagation velocity c is estimated. For example, as shown in FIG.
- ⁇ is the density of the tissue for which the elastic modulus was measured.
- an elastic modulus image is displayed in step S10 of FIG.
- the stress is calculated from the average value ⁇ ′ of the strain in the ROI_e measured for the elastic modulus and the elastic modulus E ′ calculated in step S08. Also, the elastic modulus distribution corresponding to the strain distribution is calculated from the stress and the strain distribution obtained by the conventional method.
- the calculation of the elastic modulus distribution in ROI_s corresponding to the strain distribution in ROI_s which is the first region (region 1) is performed by the color scale setting unit 50.
- the color scale setting unit 50 converts the distortion color scale into an elastic modulus color scale.
- a B-mode image 45, an elastic modulus image 46, and a color scale 43 indicating elastic modulus information are displayed on the display unit 7.
- the display range of the elastic modulus image 46 is the same as the display range of the distortion image 41 described above, and is displayed in color.
- the minimum value, the maximum value 47, etc. of the elastic modulus (absolute value) represented by the color scale 43 are displayed.
- the details of the elastic modulus determination method and the color scale minimum and maximum value calculation method using FIG. 8 will be described in the case of two effective elastic modulus measurement ROI_e. First, the measurement results ⁇ 1 and ⁇ 2 of the strain distribution are obtained from ROI1 and ROI2.
- the most probable ⁇ can be obtained by a method such as least square fitting even when there are three or more places. At this time, ⁇ is assumed to be uniform in the space. However, if there are three ROI_e, it is possible to devise a higher accuracy such as monotonically decreasing ⁇ in the depth direction.
- the elastic modulus is displayed as a color scale. However, if this method is used, in addition to the elastic modulus, other parameters related to mechanical elasticity such as the shear wave velocity and Poisson's ratio. Can also be displayed.
- the strain information is the absolute value of the hardness (such as the elastic modulus) in all the regions in the first region ROI_s where the image including the strain information that is the relative value of the hardness is displayed. ) Is converted into information, and an image including information on the absolute value of hardness is displayed in ROI_s.
- the elastic modulus measurement performed to display the elastic modulus image there is only one elastic modulus measurement location, so that the temperature rise of the living tissue near the focal point and the ultrasound probe 1 is reduced. Is possible.
- the elastic modulus information can be imaged in a region wider than the elastic modulus measurement range representing the absolute value of hardness. Furthermore, since the modulus of elasticity is measured by selecting a location with uniform hardness, it is possible to measure the modulus of elasticity with high accuracy by measuring the modulus of elasticity at one location.
- step S12 it is determined whether or not to end the diagnosis by the hybrid method.
- the hybrid end signal is input via an input device (not shown). If an end signal is input at the time of determining the end of diagnosis, the hybrid method is ended in step S14. After completion of the hybrid method, a distortion image according to the conventional method is displayed superimposed on the B mode, or only the B mode image is displayed.
- step S02 when selecting the elastic modulus measurement position, it is possible to select a portion having a uniform hardness again, so that the accuracy of measurement is improved and the robustness against body movement is improved.
- step S08 when the process returns to step S08, it is possible to reduce the signal processing time and the calculation cost.
- an ultrasonic beam for generating a shear wave is irradiated to the same place each time, safety is lower than when returning to step S02.
- the display of the strain image and the display of the elastic modulus image are alternately repeated.
- the time interval for performing the elastic modulus measurement is required to be 1 to 2 seconds or more in order to suppress the temperature rise in the living body.
- the frame rate of the elastic modulus image becomes 1 or less.
- the frame rate of the distorted image acquired in steps S02 and S04 is normally about 10. In order not to lose the real-time property of the ultrasonic image, it is desirable to update the strain image from the time when the strain image is displayed in step S04 until the time when the elastic modulus image is displayed in step S10.
- the elastic modulus image may be updated once after the distortion image is updated ten times.
- the elastic modulus measurement position selected in step S06 may be shifted in the imaging section while the strain image is being updated.
- the elastic modulus measurement position is, for example, the Motion Correction (body motion correction) method (H. Yoshikawa, et.al., Japan Journal of Applied Physics, Vol. 45, No. 5B, p. 4754, 2006) can be corrected at any time.
- the distribution of the standard deviation S of the strain amount and the difference Max-min between the maximum value and the minimum value of the strain amount as a result of signal processing. May be displayed on the display unit 7 in black and white or color, and the operator may look at the displayed image and determine the elastic modulus measurement position using an input device or the like (not shown).
- the elastic modulus measurement position can be determined by avoiding blood vessels and the like.
- the operator may determine the measurement position using an input device (not shown) by looking at the distortion image displayed in step S04 without performing signal processing. The determined position information is input to the central control unit 3.
- the distortion amount, standard deviation, black and white of Max-min, or color gradation displayed on the display unit 7 may be controlled to 3 to 256.
- the gradation is changed using an input device or the like (not shown). By setting the gradation to 3, it becomes easy to select a portion having a uniform amount of distortion. If the gradation is 256, it is possible to select a precise measurement site.
- the measurement position selection unit 40 can be omitted.
- the shape of the kernel K and the filter G a circular shape, an elliptical shape, a square shape, and other geometric shapes can be applied in addition to the rectangular shape described above.
- the measurement position selecting unit 40 does not display the place where the uniformity of distortion in the ROI is the maximum as the optimum part, but the distortion uniformity in the ROI exceeds a certain threshold. All parts are displayed as elastic modulus measurement position recommended parts. From the recommended range, the elastic modulus measurement position is determined based on the information input by the operator using the input unit 61.
- the threshold may be set so that the uniformity of distortion in the ROI is multiplied by a coefficient such as 0.8 for the maximum uniformity. Further, a method based on statistical information in which the uniformity is calculated for each place in the image and the threshold value is determined from the histogram may be used.
- a standard other than uniformity can be added to the selection of the elastic modulus measurement position.
- a standard other than uniformity can be added to the selection of the elastic modulus measurement position.
- FIGS. it is a method for more precisely determining a site where it is better to avoid irradiating ultrasonic waves for generating degeneration such as bones and blood vessels.
- the focus of the displacement generating beam is set at the position shown in FIG.
- the shape of the displacement-generating transmission beam can be calculated from the aperture diameter and frequency.
- this calculation is calculated by the displacement generation beam propagation path estimation unit 62, which is further added to the system configuration of FIG.
- the position input by the operator from the input unit 61 where it is better to avoid irradiation of bones, and the position exposed to intense ultrasonic waves are compared on the image.
- the overlap between the portion that avoids irradiation and the displacement generating beam exceeds a certain threshold value, it is determined that the focal position of the displacement generating beam is not desirable.
- this operation is sequentially performed on the entire image, a site suitable for elastic modulus measurement and a site not suitable can be separated.
- the displacement generation transmission beam shape may be estimated by using a general attenuation rate or a value estimated from an echo signal in the first ultrasonic transmission / reception unit 20.
- the measurement target is the inside of the subject, but a material with a known elastic modulus such as a coupler (not made of a material such as a polymer gel) is sandwiched between the ultrasound probe and the subject.
- ROI_s is set so as to include a coupler, and at least one of ROI_e is also set in the coupler.
- the uniformity includes a substantially uniform case that can be treated as substantially uniform.
- the second embodiment is an ultrasonic diagnostic apparatus using an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal, and radiates and receives a first displacement detection beam.
- a distortion calculation unit 24 that calculates distortion information in the region 1 based on an echo signal from the sample; a displacement generation unit that emits a focused beam into the subject to displace the tissue in the subject; Elasticity for detecting the elastic modulus in the region 2 included in the region 1 by detecting the displacement of the shear wave generated by the focused beam on the basis of the echo signal from the subject that radiated the displacement detection beam and received.
- a rate calculating unit a display unit for displaying a strain image based on strain information and an elastic modulus, and a measurement position selecting unit for selecting at least one elastic modulus detection position for detecting the elastic modulus based on the strain information.
- the measurement position selection unit 40 Fabric selected one place modulus detection position becomes uniform, two different focus positions for irradiating the focused beam, an ultrasonic diagnostic apparatus which is determined from the selected elastic modulus detection position.
- the present embodiment is the above-described ultrasonic diagnostic apparatus, in which the displacement generation unit has a beam time setting unit 14 for setting the transmission time of the focused beam, and the beam time setting unit 14 includes two different focal points.
- the present invention relates to an ultrasonic diagnostic apparatus for setting a transmission time while changing the ON / OFF switching period of a focused beam irradiating a position to be the same, and changing the ON / OFF switching period to be a chirp signal of the period or phase. .
- FIG. 10 shows a specific example of the overall configuration of the ultrasonic diagnostic system for carrying out the second embodiment. 1 is different from the system configuration of the first embodiment in FIG. 1 in that a beam time setting unit 14 is added to the displacement generation unit 10 and a hardness spectrum calculation unit 35 is installed instead of the elastic modulus calculation unit. It is.
- the output of the shear wave displacement calculation unit 32 is input to the hardness spectrum calculation unit 35, and a value related to hardness is calculated.
- the hardness spectrum calculator 35 can be realized by executing a program in a central processing unit (CPU) functioning as a processing unit.
- CPU central processing unit
- a method of transmitting a focused beam for generating displacement by the burst chirp method will be described with reference to FIGS.
- the two displacement-generating transmission beams are controlled so as to cause displacement alternately at the focal point F1 and the focal point F2 of the subject tissue shown in FIG.
- ON / OFF of the irradiation of the transmission beam for generating the displacement to each focal point is controlled by the central control unit 3, and the ON / OFF switching time is set by the newly added beam time setting unit 14.
- Fig. 12 shows the sequence of the displacement-generating transmission beam and the displacement-detecting transmission / reception beam.
- the sequence of the displacement-generating transmission beam at the time is shown.
- this generation ultrasonic beam transmission method is called a burst chirp method.
- ON / OFF of irradiation of the displacement detection transmission / reception beam is performed by the central control unit 3 controlling the second ultrasonic transmission / reception unit 30.
- This reference signal shear wave displacement calculation unit 32 is used for calculation of shear wave displacement.
- the displacement-generating transmission beam to the focal point F1 is always ON when 0 ⁇ t ⁇ T1.
- the displacement-generating transmission beam at the focal point F1 is turned off.
- the displacement-generating transmission beam at the focal point F2 is turned on, displacement occurs at the focal point F2, and a shear wave propagates.
- the displacement-generating transmission beam to the focal point F1 is OFF, and the displacement-generating transmission beam to the focal point F2 is ON when T1 ⁇ t ⁇ T1 + T1.
- the switching period of the two displacement-generating transmission beams is T1.
- ⁇ Change the ON / OFF switching period Tm of the displacement-generating transmission beam, and repeat the irradiation of the displacement-generating transmission beam to the two focal points and the transmission / reception of the displacement-detecting transmission / reception beam.
- m represents a cycle in which the focus F1 and the focus F2 are turned on for the mth time
- m 1, 2, 3,.
- the magnitude of the sound intensity of the burst signal to each focal point may be the same or different.
- shear waves generated at the focal points F1 and F2 interfere with each other while propagating, and cancel or amplify.
- heat is generated simultaneously with displacement at each focal point.
- the distance between the two focal points is d.
- the degree of interference increases.
- the temperature rise E between the focal points becomes larger than the temperature at the focal point due to heat conduction, and safety is lowered.
- the optimum value of d is a value at which the maximum value of the temperature rise is equal to the maximum value of the temperature rise at each focal point, and wave interference occurs.
- the optimum value d depends on the depth of focus, the irradiation time of the transmission beam for generating the displacement, the frequency, the diagnostic part (which affects the sound speed, ultrasonic absorption, thermal conductivity, etc. of the living body). Since the thermal conductivity of a living body is about 0.6 W / m / K, when the irradiation time is a few milliseconds, the range of temperature rise near the focal point is about the same as the width of the displacement-generating transmission beam. Yes, a distance d equal to or greater than the width of the displacement-generating transmission beam is required.
- the width of the displacement-generating transmission beam is approximately equal to the beam width of the focusing transducer, and is the diameter of the region (circle) where the ultrasonic energy density at the focal point is zero for the first time.
- FIG. 15 shows the relationship between the distance d and the temperature rise E.
- the temperature rise has a maximum value at the focal position, and is zero at a location farther than the beam width d. Therefore, as shown in the middle and upper parts of the figure, when the distance d is equal to or larger than the beam width, the maximum value of the temperature rise at each focal point and the maximum value of the total temperature rise obtained by adding the temperature rises. E_max is the same. Further, as shown in the lower part of the figure, when the distance d is smaller than the beam width, the maximum value E_max of the overall temperature rise is ⁇ E higher than the maximum value of the temperature rise at each focal position indicated by the dotted line. Only increase. Therefore, it can be seen that if the distance d is equal to or greater than the beam width, the safety with respect to the temperature rise is maintained.
- the setting of the optimum value d for measurement will be described.
- the F value is 1, and the carrier frequency is 2 MHz, the beam width is 1.8 mm.
- the F value is calculated as depth of focus / aperture diameter, and the beam width is calculated as (2.44 * F value * carrier signal wavelength).
- the wavelength ⁇ of the shear wave is about 0.2 mm.
- the maximum value of the propagation distance of the shear wave that can be detected by experimental data is about 6 mm.
- the value of d is set such that 10 ⁇ ⁇ d ⁇ 30 ⁇ using the beam width, maximum propagation distance, and shear wave wavelength.
- the focal depth is 2 cm
- the F value is 1
- the carrier frequency is 7 MHz
- the beam width is calculated as 0.5 mm.
- the wavelength ⁇ of the shear wave is about 0.1 mm.
- the maximum value of the propagation distance of the shear wave that can be detected by experimental data is about 3 mm. Therefore, d is set so that 5 ⁇ ⁇ d ⁇ 30 ⁇ .
- the central control unit 3 reads this value d from a memory (not shown) and sets it in the focus position setting unit 12. Further, a value related to the switching period is determined from the set value of d and the predicted value of the propagation speed of the shear wave.
- the switching period when the shear waves generated at the two focal points interfere to increase the amplitude is obtained, and the elastic modulus is measured.
- a condition for amplifying the interference wave will be described.
- the condition in which the interference wave is amplified and the absolute value of the displacement amount reaches the peak value (maximum value) is that the distance d between the two focal points is (n + 1/2) times the wavelength ⁇ .
- the curve 103 and Equation 1 can be used.
- the switching frequency fm at this time is represented by fM (n).
- the shear wave propagation velocity is a value inherent to the tissue properties.
- Equation 1 the value of the switching period Tm serving as a peak value
- Equation 2 Equation 2
- TM (n) d / c * (2 / (2n + 1))
- the shear wave propagation velocity c depends on the hardness of the tissue, and the harder the c, the larger the c. Therefore, from the TM (n) value, the shear wave propagation velocity, the elastic modulus, etc. It becomes possible to estimate the hardness of the tissue. It is desirable to control Tm, which is the ON / OFF switching cycle of the displacement-generating transmission beam, in the range of several tens of Hz to several kHz. In addition, if the hardness is estimated from TM (n) for a large n, the total irradiation time of the focused beam for generating displacement can be shortened, and the temperature rise can be suppressed.
- the central control unit 3 turns on / off the irradiation of the displacement generating transmission beam to each focal point so that the total irradiation time of the displacement generating transmission beam to each focal point is 1 ms or less. It is desirable to be controlled in
- this technology is characterized in that the elastic modulus is measured not by the period of the carrier signal of the displacement-generating transmission beam but by ON and OFF time control. Therefore, by increasing the carrier frequency, the beam width is narrow, and imaging can be performed with high spatial resolution.
- the hardness spectrum calculation unit 35 spectrum analysis is performed on the output signal from the shear wave displacement calculation unit 32 to calculate fM having the maximum amplitude value and TM corresponding thereto, the shear wave propagation speed, Calculate the values for hardness such as elastic modulus and shear modulus.
- FIG. 16 shows a signal waveform of a shear wave generated from one focal point and an interference waveform of the shear wave generated from two focal points.
- the minimum amplitude when the displacement generating transmission beam is irradiated to one focal point and the displacement of the shear wave can be observed is 1.
- FIG. 16A when the amplitude of the shear wave generated by the displacement-generating transmission beam to each focal point is 1, and when the switching period is TM, FIG.
- the amplitude of the interference wave is larger than the amplitude before the interference (ideally twice), and the generation efficiency of the displacement with respect to the transmission, that is, the transmission sensitivity is increased.
- the displacement-generating transmission wave there is also a method of improving safety by providing a time during which the beam is turned on by the last n% (n is a positive real number) of the time Tm when the beam is turned on. In this case, it should be noted that the switching cycle Tm does not change even when the time for turning OFF is provided.
- Steps S00 to S04 and Steps S10 to S14 are the same as the processing flow of FIG. Steps S06 and S08 will be described only for processing different from the first embodiment.
- step S06 the positions of the focus F1 and the focus F2 are set.
- the position of the two focal points for example, the central point of the two focal points (here, the central point of the straight line connecting the two focal points of each set) is set as POI (Point of Interest), and the distance between the two focal points is Is set.
- the distance between the two focal points is set to a value smaller than the distance at which the two shear waves interfere with each other and larger than the width of the displacement-generating transmission beam irradiated to each focal point.
- the POI coordinate is determined as POI (x1, y1), and the distance is determined as lx.
- the Information on the POI and distance is set in the focal position setting unit 12.
- the coordinates of the two focal points are F (x1 ⁇ x, y1) and F (x1 + ⁇ x, y1).
- ⁇ x 1x / 2.
- the value of n in Equations 1 and 2 and the optimal observation point are determined from the expected shear wave velocity.
- the observation point is determined based on the maximum point of the absolute value of the displacement of the shear wave or a plurality of positions including the maximum point within the shear wave propagation distance.
- the observation point is read and set from a storage medium (not shown).
- the raster used for detecting the amplitude (several ⁇ m to several tens of ⁇ m) of the shear wave propagation on the observation point and the sampling point on the raster are determined.
- the displacement detection beam reception PRF (Pulse Repetition : Frequency) is set so as to satisfy the Nyquist theorem for the expected shear wave frequency. For example, when the raster is in the same direction as the shear wave displacement, the PRF is set to be twice or more the shear wave frequency.
- the determined n and observation position may be displayed on the screen.
- an initial value Tstart, a final value Tend, and an interval ⁇ T of the switching cycle of the two-focus displacement-generating transmission beam irradiation are set.
- the initial value Tstart and the final value Tend are set to values at which a peak is obtained in a range that satisfies Equations 1 and 2 with respect to the measurement site and the distance d between the two focal points.
- the initial value Tstart, the final value Tend, and the interval ⁇ T are read from a storage medium (not shown) in the central control unit 2 according to the measurement site, depth, distance between the focal points, and the like. 14 is set.
- a signal is transmitted and a shear wave is generated at two focal points.
- step S24 an ultrasonic signal for detecting the shear wave is received on the raster for detecting the amplitude of the shear wave propagation.
- Displacement detection at each measurement point may be always detected from when the burst chirp signal is turned off until a shear wave arrives and passes through all observation points, or a shear wave arrives. It is also possible to obtain the time until passing through in advance from the shear wave velocity estimated as the distance between the focal point and the observation point and detect only that time. Since the latter can increase the PRF, it is possible to detect displacement with high accuracy.
- the second ultrasonic transmission / reception unit 30 extracts a signal corresponding to fm from the received wave signal by signal processing such as a bandpass filter, and then the shear wave displacement calculation unit 32 performs a known correlation calculation or the like. Calculate the displacement of the shear wave.
- the correlation calculation is a calculation performed using an echo signal for each time received by the reference signal and the displacement detection beam. By this calculation, the time waveform of the shear wave amplitude at each observation point is obtained.
- step S26 it is determined whether or not the immediately preceding switching cycle Tm is Tend. If it is not Tend, the process returns to step S22 to transmit a burst chirp signal at the next Tm + 1 switching cycle.
- step S22 when the reference signal is acquired again, the robustness of the correlation calculation due to the shift of the focus position during measurement is large. For transmission of a burst chirp signal after the switching period of Tm + 1, if the reference signal is not acquired and the correlation calculation is performed using the reference signal acquired first, the measurement time can be shortened. it can.
- step S08 in FIG. 17 the hardness spectrum calculation unit 35 calculates values relating to the hardness such as the shear wave propagation velocity, the elastic modulus, and the shear elastic modulus. If the diagnosis is not terminated in step S12 of FIG. 17, the process returns to step S02, step S08, step S20, or step S22.
- ⁇ TM is derived from Equation 2 as Equation 3.
- Equation 3 d / c * ( ⁇ 4 / ((2n + 1) * (2n + 3))) --- (Equation 3)
- the control is performed so that the transmission beam for generating the displacement to the two focal points is alternately turned ON / OFF, it is also possible to cause the displacement at the same time by controlling the ON / OFF to be turned ON / OFF at the same time.
- the condition that the interference wave is amplified and has a peak value is when d is (n + 1) times the wavelength ⁇ , and therefore, Equation 4 is obtained corresponding to Equation 1.
- next switching cycle T (m + 1) instead of changing to the next switching cycle T (m + 1) by one ON / OFF control for each switching cycle Tm, the next switching cycle is set after repeating ON / OFF several times for the same Tm. You may make it do. By repeating ON / OFF several times for the same Tm, measurement with higher sensitivity becomes possible.
- Tm is changed from large to small, it is also possible to change from small to large, and ⁇ Tm is changed based on an arbitrary function other than a fixed value such as a geometric series. Is also possible.
- a minimum value may be used.
- the irradiation of the displacement-generating transmission beam and the transmission / reception of the displacement detection beam are performed once.
- the spectrum analysis is performed to calculate the displacement for a plurality of fm.
- the system configuration of the third embodiment (not shown) is obtained by adding the beam time setting unit 14 and the hardness spectrum calculation unit 35 of FIG. 10 to the system configuration of FIG.
- the output signal from the shear wave displacement calculation unit 32 is input to both the elastic modulus calculation unit 34 and the hardness spectrum calculation unit 35 installed in parallel.
- the central control unit 3 is connected to both the beam time setting unit 14 and the hardness spectrum calculation unit 35. Either one of the elastic modulus calculation unit 34 and the hardness spectrum calculation unit 35 arranged in parallel is selected and used by the operator.
- the elastic modulus calculation unit 34 and the hardness spectrum calculation unit 35 may be collectively referred to as an elastic modulus calculation unit.
- the central control unit 3 stores an unillustrated elastic modulus measurement method according to the depth of the elastic modulus measurement position, the carrier frequency of the ultrasonic beam for shear wave generation, and the measurement site. It can be read from the medium.
- a focused ultrasound that irradiates a portion where a safer measurement is necessary (near a blood vessel or the like) or the number of channels that the ultrasonic probe 1 can be used is small.
- the burst chirp method can be applied when the acoustic intensity of the beam is smaller than a certain threshold.
- the depth from the body surface of the elastic modulus measurement position is a certain value If it is larger, or if a value represented by a function of these values and the carrier frequency of the ultrasonic beam for generating a shear wave is larger than a certain threshold value, the burst chirp method is applied.
- the fourth embodiment is an ultrasonic diagnostic apparatus using an ultrasonic probe that transmits an ultrasonic beam to a subject and receives an echo signal.
- the ultrasonic diagnostic device emits a first displacement detection beam and receives the received wave.
- a distortion calculation unit that calculates distortion information in the region 1
- a displacement generation unit 10 that radiates a focused beam into the subject and displaces the tissue in the subject
- a second Elasticity for detecting the elastic modulus in the region 2 included in the region 1 by detecting the displacement of the shear wave generated by the focused beam on the basis of the echo signal from the subject that radiated the displacement detection beam and received.
- a rate calculating unit a display unit for displaying a strain image based on strain information and an elastic modulus, and a measurement position selecting unit for selecting at least one elastic modulus detection position for detecting the elastic modulus based on the strain information.
- the measurement position selection unit 40 Ultrasonic diagnosis in which a plurality of elastic modulus detection positions at which the cloth is uniform are selected, and at least one focal position where the focused beam is irradiated is determined from the elastic modulus detection positions selected by the measurement position selection unit 40 Relates to the device.
- the output from the measurement position selection unit 40 is output.
- a specific example of the configuration of an ultrasonic diagnostic apparatus in which there are a plurality of optimum positions for performing the measured elastic modulus will be described.
- the position P output from the measurement position selection unit 40 is P (x1, y1), P (x2, y2), and P (x3, y3) will be described.
- the number of positions P output from the measurement position selection unit 40 is not limited to three and is an integer of 2 or more.
- the measurement position selection unit 40 calculates the distances L1, L2, and L3 between the positions P, and the distance information is output to the central control unit 3 together with the position information.
- the central control unit 3 outputs the focal position F for the three positions P to the focal position setting unit 12. At this time, the order of the focal positions output from the central control unit 3 is controlled by the magnitude of the distance between the focal points. In the example of FIG. 18, L3> L2> L1.
- the central control unit 3 determines the order of the three focus positions F to be output so that the distance between the two focus positions F to be output at the preceding and following times becomes large to small.
- the elastic modulus is measured for each of the three converging positions, and the elastic modulus image is updated each time the elastic modulus is measured.
- the processing flow of the hybrid system in the fourth embodiment will be described using the example of FIG. Here, the case where it applies to Example 1 is demonstrated. Only parts different from the first embodiment will be described.
- the focal points F determined for the positions P (x1, y1), P (x2, y2), P (x3, y3) are F (x1- ⁇ x, y1), F (x2- ⁇ x, y2), Let F (x3- ⁇ x, y3).
- step S06 in FIG. 3 the above-described three positions are selected, and the distance is calculated.
- step S08 in FIG. 3 the central control unit 3 outputs the focal point F (x1- ⁇ x, y1) to the focal position setting unit 12, and the displacement-generating transmission beam is output to the focal point F (x1- ⁇ x, y1).
- the elastic modulus is measured by irradiating the focused beam.
- step S10 if the measurement is not completed in step S12, the process returns to step S08, and then F (x3- ⁇ x, y3) is transferred from the central control unit 3 to the focal position setting unit 12.
- the elastic modulus is measured by irradiating the focal point F (x3- ⁇ x, y3) with the focused beam.
- step S10 After the elastic modulus image is displayed in step S10, if the measurement is not completed in step S12, the process returns to step S08, and then F (x2- ⁇ x, y2) is transferred from the central control unit 3 to the focal position setting unit 12.
- the elastic modulus is measured by irradiating the focused beam onto the focal point F (x2- ⁇ x, y2).
- step S10 After the elastic modulus image is displayed in step S10, if the measurement is not completed in step S12, the strain distribution is measured by returning to step S02, and a plurality of elastic modulus measurement positions are selected again in step S06. Is done. At this time, instead of selecting a plurality of elastic modulus measurement positions again, the three positions selected at the previous time may be used.
- step S06 by selecting a plurality of positions where the modulus of elasticity measurement is performed in step S06, calculation time and measurement time can be shortened.
- the second and subsequent positions are selected at positions away from the focal position irradiated at the previous time. It is possible to minimize the influence of the temperature rise caused when the focused beam corresponding to the position is irradiated on the temperature rise of the other focus positions. Moreover, since the focused beam is not continuously irradiated to the same place, it is possible to reduce the local temperature rise of heat.
- Example 5 another elastic modulus measurement position selection method in Examples 1 to 3 will be described.
- the description will be made on the premise of the system configuration of the first embodiment.
- a case where the size of the kernel K and the filter G is the same as that of the region ROI_e where the elastic modulus is measured will be described.
- the size of the kernel K or the filter G is reduced by half in the y direction, and signal processing for selecting the elastic modulus measurement position is performed.
- the + y direction is the propagation direction of the focused beam that is the displacement-generating transmission beam, and the + x direction shear wave propagation direction.
- the length of ROI_e in the y direction is determined by the width of the displacement generating beam in the depth direction.
- the kernel K or the filter G whose size in the y direction is reduced to half is called a kernel K ′ or a filter G ′.
- the length in the depth direction of the kernel K or the filter G is ly and the length in the azimuth direction is lx, the length in the depth direction of the kernel K ′ or the filter G ′ is ly / 2, the azimuth direction
- the length of is lx.
- the length in the depth direction is desirably a value that is 1 ⁇ 2 of the width in the depth direction of the displacement-generating transmission beam.
- the region R ′ (n) (n is a positive integer) having a region larger than the filter G ′, and the barycentric position P ′ (x) of R ′ (n) , Y) are all derived by image processing.
- the threshold value is read from a storage medium (not shown) or the like by the central control unit 3 and automatically set, or is inputted manually by a user via a recording medium (not shown).
- the measurement position selection unit 40 selects the kernel K ′ or Two positions P ′ (x, y) where the filter G ′ is continuous are searched. For example, when one position is P1 '(x', y ') and the other position is P2' (x ', y' + ly / 2), two positions P '(x, y).
- the x coordinate of F ′ is x′ ⁇ lx / 2
- the y coordinate is y ′ + ly / 4.
- it is a value corresponding to a half of the propagation distance of lx / 2 shear wave.
- y ′ + ly / 4 is a value corresponding to the midpoint of the y coordinates of the positions P1 ′ and P2 ′, that is, a half of the length in the depth direction of the region ROI_e where the elastic modulus measurement is performed.
- Information on the focal point F ′ is input to the focal position setting unit 12, and the focal point F ′ is irradiated with a transmission beam for generating displacement. Further, the displacement of the shear wave propagating from the focal point F ′ in the x direction is detected.
- the elastic modulus calculation unit 34 calculates the elastic modulus in ROI_e
- two elastic moduli are calculated: an elastic modulus at a location shallower than the focal point F ′ and an elastic modulus at a location deeper than the focal point F ′.
- Example 5 when measuring the elastic modulus per place, as shown in FIG. 20, the relationship between two strains ⁇ ′ and ⁇ ′′ and the two elastic moduli E ′ and E ′′ corresponding to these strains. Can be investigated. Therefore, two stresses are calculated from the two sets of strain and elastic modulus.
- the color scale setting unit 50 converts the distortion color scale into the elastic modulus color scale.
- the estimation of the elastic modulus related to the strain between ⁇ ′ and ⁇ ′′ can be obtained by interpolation from E ′ and E ′′.
- the desired distortion is not between ⁇ ′ and ⁇ ′′, it can be stopped by extrapolation.
- linear interpolation such as interpolation and extrapolation from two stresses, it is possible to improve the accuracy of conversion.
- the estimation accuracy of elastic modulus is improved in a wide range when ⁇ ′ and ⁇ ′′ are as different values as possible.
- the fitting process is performed by the color scale setting unit 50 or the central control unit 3, and the calculation of the fitting is performed by a program process of a central processing unit (not shown). If there are no two kernels K ′ or filters G ′ that are continuous in the depth direction, the size of the kernel or filter, or the threshold value is automatically or manually changed. Further, when two kernels K ′ or a filter G ′ that are continuous in the depth direction do not exist, signal processing is performed again using the kernel K and the filter G, and the elastic modulus measurement position is selected.
- the elastic modulus may be measured by the methods of Example 2 and Example 3.
- the kernel can be reduced to 1 / m (m is a positive integer).
- the elastic modulus and strain pair calculated by the elastic modulus calculator 34 is multiplied by m, so that the accuracy of the stress calculated by averaging is improved. Is possible.
- the positions of the two focal points are, for example, F ′ (x′ ⁇ lx / 2, y ′ + ly / 4), F ′′ (x ′) in FIG. + Lx / 2, y ′ + ly / 4).
- the size of the kernel K and the filter G may be reduced by half in the x direction, that is, in the shear wave propagation direction.
- the length in the azimuth direction of the reduced kernel K ′ and filter G ′ is preferably a value that is 1 ⁇ 2 of the propagation distance of the shear wave.
- each embodiment examples include living bodies such as liver, mammary gland, blood vessel, and prostate.
- Each embodiment is a technique for synthesizing an elastic modulus measurement using a shear wave and a distorted color image, and one of the features is that the acquisition of the distorted color image is performed prior to the elastic modulus measurement.
- strain color images By acquiring strain color images first, it is possible to select a location suitable for elastic modulus measurement, and greatly improve the accuracy of elastic modulus measurement and the accuracy of composite images based on elastic modulus measurement. .
- each embodiment when generating a shear wave, instead of an ultrasonic focused beam, mechanical drive (DC motor, vibration pump, etc.), manual compression, compression by electric pulse, movement of living tissue such as heart and blood vessels A known method such as the above may be used. Further, a two-dimensional probe may be used instead of the linear array type probe.
- Each element of the ultrasonic probe 1 may be changed to, for example, a ceramic, a piezoelectric element made of a polymer, a vibrator using the electrostatic force of silicon, or the like.
- the displacement may be calculated using a known calculation method such as cross-correlation calculation, least square sum, Doppler method in addition to correlation calculation.
- the configurations, functions, processing units, etc. of the above-described embodiments may be partly or entirely configured as a dedicated hardware configuration, a software configuration, or a configuration in which they are shared.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Animal Behavior & Ethology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
圧縮方向の変位が推定され,変位の空間微分量である歪みを求め,画像化する。この方法では,撮像対象として,体表からの圧迫が容易なところに存在する臓器(例えば乳腺など)には極めて有効な方法であるが,必ずしも全ての撮像対象に対して有効では無い場合もある。例えば,体表と肝臓の間に介在層として,すべり面が存在するため,十分な変位を生じさせるような圧迫が困難である。また,従来方式で,弾性率を推定しようとする場合,境界条件が必要になるが,すべり面が存在すると境界条件が複雑となり,弾性率の推定が困難である。
そこで,変位生成用送波ビームとしての集束ビームを用い,被検体内部に放射圧を印加し,介在層の影響を抑えて対象組織を変位させて,硬さの診断をする技術がある。例えば,特許文献1に記載されるARFI(Acoustic Radiation Force Impulse)Imagingがある。この技術では,集束ビームの進む方向に生じる組織の変位量を画像化したり,焦点での組織変位に伴って集束ビームの進む方向とは垂直な方向に生じるずり波の伝搬速度の推定からずり波伝搬領域内の平均的なずり弾性率といった弾性係数(以下,弾性率と呼ぶ)を測定し画像化したりする。ずり波生成による弾性率計測技術を用いると,上記のすべり面などの介在層の影響を低減する効果以外にも,超音波によって組織を変位させるため,手技依存性が少ない診断が期待される。
測定位置選定部40は,領域2を選定する際に,画像処理で領域1内の歪み分布の輪郭を抽出することによって選定する超音波診断装置,或いは上記の超音波診断装置であって,歪み情報に基づく歪み画像と弾性率を表示する表示部7を更に備え,測定位置選定部40は,操作者が表示部7に表示された歪み画像に基づき領域2を選定可能である超音波診断装置に関するものである。
3)が得られる。ここでは,説明を簡単にするため,領域の数を3つとしているが,領域の数は2以上の任意の整数である。領域複数の領域の中から,領域R(n)(n=1,2,3,・・・)が弾性率の計算領域ROI_e,すなわち,フィルタGより大きい領域を持つ領域R(n)が最適な弾性率測定の領域として選定される。
この硬さスペクトル計算部35は,処理部として機能する中央処理部(Central Processing Unit:CPU)におけるプログラム実行により実現することができる。
k*d =(2πfM(n)/c)*d = 2π(n+1/2)‐‐‐(数式1)
ただし,kは波数(=2π/λ),cはずり波の伝搬速度,nは0もしくは正の整数(n=0,1,2,・・・)である。ずり波の伝搬速度は,組織性状に固有の値である。
TM(n)= d/c*(2/(2n+1)) ‐‐‐(数式2)
例えば,n=1, d=2[mm]であるとき,c=1[m/s]ではTM(1)=1.3[ms](fM(1)=750[Hz])となり,c = 5[m/s]では,TM(1)= 0.3[ms](fM(1)=3.8[kHz])となる。先に記述したように,ずり波の伝搬速度cは組織の硬さに依存し,硬いほどcは大きくなるので,TM(n)の値から,ずり波の伝搬速度,および,弾性率等の組織の硬さを推定することが可能となる。変位生成用送波ビームのON/OFF切替え周期であるTmは,数十Hz~数kHzの範囲で制御するのが望ましい。また,大きなnに対するTM(n)から硬さを推定した方が,変位生成用の集束ビームの合計照射時間を短縮することができ,温度上昇を抑制することができる。温度上昇を抑制するため,各焦点への変位生成用送波ビームの合計照射時間が1ms以下となるように,各焦点への変位生成用送波ビームの照射のON/OFFが中央制御部3において制御されることが望ましい。
14に設定される。
ΔTM=d/c*(-4/((2n+1)*(2n+3)))‐‐‐(数式3)
また,2焦点への変位生成用送波ビームのON/OFFが交互になるように制御したが,同時にON/OFFとなるように制御して同時刻に変位を生じさせることも可能である。この場合は,干渉波が増幅してピーク値を持つ条件はdが波長λの(n+1)倍のときであるので,数式1に対応して,数式4のようになる。
k*d= (2πf/c)*d = 2π(n+1) ‐‐‐(数式4)
高精度な弾性率の測定方法として,1回目に,切替え周期の間隔ΔTを粗く設定して測定し,TMを求め,次の測定で,TMの付近の切替え周期でΔTをより細かい設定にしてより詳細なTMの値を求めるようにしてもよい。
2…送受切替スイッチ,
3…中央制御部,
4…カラーDSC,
5…白黒DSC,
6…合成部,
7…表示部,
10…変位生成部,
11…変位生成用送波波形生成部,
12…焦点位置設定部,
13…変位生成用送波ビーム生成部,
14…ビーム時間設定部,
20…第一超音波送受信部,
22…変位演算部,
24…歪み演算部,
30…第二超音波送受信部,
32…ずり波変位演算部,
34…弾性率演算部,
35…硬さスペクトル計算部,
41…歪み画像,
42…カーネルK,
43…弾性率カラースケール,
44…フィルタG,
45…Bモード像,
46…弾性率画像,
47…弾性率(絶対値,最大-最小),
40…測定位置選定部,
50…カラースケール設定部,
60…弾性率測定位置演算部,
61…入力部,
62…変位生成用ビーム伝搬経路推定部,
100…超音波探触子1の各素子。
Claims (15)
- 被検体に超音波ビームを送信し,エコー信号を受信する超音波探触子を用いる超音波診断装置であって,
第一の変位検出用ビームを放射し,受波した前記被検体からのエコー信号に基づき,領域1内の歪み情報を演算する歪み演算部と,
前記被検体内に集束ビームを放射して,前記被検体内の組織を変位させる変位生成部と,第二の変位検出用ビームを放射し,受波した前記被検体からのエコー信号に基づき,前記集束ビームにより生じるずり波の変位を検出して,前記領域1の内部に含まれる領域2内の弾性率を検出する弾性率演算部と,
前記歪み情報に基づく歪み画像と前記弾性率を表示する表示部を有する,
ことを特徴とする超音波診断装置。 - 請求項1に記載の超音波診断装置であって,
前記歪み情報に基づき,前記弾性率を検出する,少なくとも1つの弾性率検出位置を選定する測定位置選定部を更に有し,
前記集束ビームを照射する少なくとも1つの焦点位置が,前記測定位置選定部によって選定された少なくとも1つの前記弾性率検出位置から決定される,
ことを特徴とする超音波診断装置。 - 請求項2に記載の超音波診断装置であって,
前記測定位置選定部は,歪みの分布が均一となる弾性率検出位置を1箇所選定し,
前記集束ビームを照射する異なる2つの前記焦点位置が,選定された前記弾性率検出位置から決定される,
ことを特徴とする超音波診断装置。 - 請求項3に記載の超音波診断装置であって,
前記変位生成部は前記集束ビームの送波時間を設定する送波ビーム時間設定部を有し,
前記送波ビーム時間設定部は,異なる2つの前記焦点位置に照射する前記集束ビームのON/OFF切替え周期が同じで,かつ,前記ON/OFF切替え周期が,周期または位相のチャープ信号となるように変化させながら前記送波時間を設定する,
ことを特徴とする超音波診断装置。 - 請求項2に記載の超音波診断装置であって,
前記測定位置選定部は,歪みの分布が均一となる弾性率検出位置を複数箇所選定し,
前記集束ビームを照射する少なくとも1つの前記焦点位置が,前記測定位置選定部において選定された複数箇所の前記弾性率検出位置から決定される,
ことを特徴とする超音波診断装置。 - 請求項2に記載の超音波診断装置であって,
前記測定位置選定部は,操作者が前記表示部に表示された画像を見ながら前記弾性率検出位置を選択することが可能である,
ことを特徴とする超音波診断装置。 - 請求項1に記載の超音波診断装置であって,
前記表示部は,表示される前記歪み画像の表示範囲における前記弾性率を示すカラースケールを表示する,
ことを特徴とする超音波診断装置。 - 被検体に超音波ビームを送信し,エコー信号を受信する超音波探触子を用いる超音波診断装置であって,
第一の変位検出用ビームを放射し,受波した前記被検体からのエコー信号に基づき,領域1内の歪み情報を演算する歪み演算部と,
前記歪み情報に基づき,前記領域1の内部に含まれる領域2を選定する測定位置選定部と,
前記被検体内に集束ビームを放射して,前記被検体内の組織を変位させる変位生成部と,第二の変位検出用ビームを放射し,受波した前記被検体からのエコー信号に基づき,前記集束ビームにより生じるずり波の変位を検出して,前記領域2内の弾性率を検出する弾性率演算部を備える,
ことを特徴とする超音波診断装置。 - 請求項8に記載の超音波診断装置であって,
前記測定位置選定部は,前記領域2を選定する際に,前記領域1内の歪みの分布の標準偏差,もしくは,最大値と最小値の差がある閾値より小さくなる位置を求める,
ことを特徴とする超音波診断装置。 - 請求項8に記載の超音波診断装置であって,
前記弾性率演算部は,
前記領域2内の弾性率,および前記領域2内の前記歪み情報を用いて応力を計算し,前記領域1内の前記歪み情報,および,前記応力から,前記領域1の前記弾性率を演算する,ことを特徴とする超音波診断装置。 - 請求項8に記載の超音波診断装置であって,
前記歪み情報に基づく歪み画像と前記弾性率を表示する表示部を更に備え,
前記測定位置選定部は,操作者が前記表示部に表示された前記歪み画像に基づき前記領域2を選定可能である,
ことを特徴とする超音波診断装置。 - 被検体に超音波ビームを送信し,被検体からのエコー信号を受信する超音波探触子を用い,受信した当該エコー信号に基づき画像を表示部に表示する超音波表示方法であって,
第一の変位検出用ビームを放射して被検体からのエコー信号を受波し,第一の領域内の歪み情報を演算し,
演算した前記歪み情報に基づく歪み画像を前記表示部に表示し,
前記被検体内に集束ビームを放射して前記被検体の組織を変位させ,
前記第二の変位検出用ビームを放射して前記被検体からのエコー信号を受信し,前記集束ビームにより生ずるずり波の変位を検出し,
前記ずり波の変位に基づき,前記第一の領域に含まれる第二の領域内の弾性率を演算し,演算した前記弾性率を前記表示部に表示する,
ことを特徴とする超音波表示方法。 - 請求項12に記載の超音波表示方法であって,
前記歪み情報に基づき選定された前記集束ビームの焦点位置を決定する,
ことを特徴とする超音波表示方法。 - 請求項12に記載の超音波表示方法であって,
前記第二の変位検出用ビームを放射する前記第二の領域を選定する際,前記表示部に表示された前記歪み画像に基づき,前記歪み情報が均一な場所から選定する,
ことを特徴とする超音波表示方法。 - 請求項12に記載の超音波表示方法であって,
前記表示部に,表示される前記歪み画像の表示範囲における前記弾性率を示すスケールを表示する,
ことを特徴とする超音波表示方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP11857846.7A EP2671511B1 (en) | 2011-02-04 | 2011-12-28 | Ultrasound diagnostic apparatus and method |
US13/982,801 US20130317361A1 (en) | 2011-02-04 | 2011-12-28 | Ultrasound diagnostic apparatus and method |
CN201180066724.0A CN103347450B (zh) | 2011-02-04 | 2011-12-28 | 超声波诊断装置及方法 |
JP2012555710A JP5619191B2 (ja) | 2011-02-04 | 2011-12-28 | 超音波診断装置,および方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011-023152 | 2011-02-04 | ||
JP2011023152 | 2011-02-04 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012105152A1 true WO2012105152A1 (ja) | 2012-08-09 |
Family
ID=46602382
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/080448 WO2012105152A1 (ja) | 2011-02-04 | 2011-12-28 | 超音波診断装置,および方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US20130317361A1 (ja) |
EP (1) | EP2671511B1 (ja) |
JP (1) | JP5619191B2 (ja) |
CN (1) | CN103347450B (ja) |
WO (1) | WO2012105152A1 (ja) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454640A (zh) * | 2013-09-13 | 2013-12-18 | 北京东方惠尔图像技术有限公司 | 一种超声波成像装置和成像方法 |
CN103462643A (zh) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | 一种剪切波速度测量方法、装置和系统 |
WO2014207668A1 (en) * | 2013-06-26 | 2014-12-31 | Koninklijke Philips N.V. | Elastography measurement system and method |
JP2015009040A (ja) * | 2013-07-01 | 2015-01-19 | 株式会社東芝 | 超音波診断装置及び超音波イメージングプログラム |
JP2015058193A (ja) * | 2013-09-19 | 2015-03-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
JP2015084909A (ja) * | 2013-10-30 | 2015-05-07 | 株式会社東芝 | 超音波画像診断装置 |
EP2881041A1 (en) | 2013-11-08 | 2015-06-10 | Hitachi Aloka Medical, Ltd. | Apparatus and method for ultrasonic diagnosis |
CN104706377A (zh) * | 2013-12-13 | 2015-06-17 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
JP2015202373A (ja) * | 2014-04-16 | 2015-11-16 | 株式会社東芝 | 超音波診断装置及び制御プログラム |
EP2939598A4 (en) * | 2012-12-25 | 2016-08-17 | Hitachi Aloka Medical Ltd | ULTRASONIC DIAGNOSTIC DEVICE AND ELASTICITY EVALUATION METHOD |
JP2018516135A (ja) * | 2015-06-04 | 2018-06-21 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | がん悪性度マップにより拡張された精密診断及び治療に対するシステム及び方法 |
KR20190075773A (ko) * | 2017-12-21 | 2019-07-01 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그 표시 방법 |
US11995807B2 (en) | 2019-05-17 | 2024-05-28 | Supersonic Imagine | Method for ultrasound determination of a corrected image of a medium, and device for implementing this method |
Families Citing this family (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101580584B1 (ko) * | 2013-11-28 | 2015-12-28 | 삼성전자주식회사 | 탄성 영상 내에 종양을 표시하는 방법 및 이를 위한 초음파 장치 |
US20150148676A1 (en) * | 2013-11-28 | 2015-05-28 | Samsung Electronics Co., Ltd. | Method and ultrasound apparatus for marking tumor on ultrasound elastography image |
KR20150070859A (ko) * | 2013-12-17 | 2015-06-25 | 삼성전자주식회사 | 전단파를 이용하여 관심 영역에 대한 탄성 정보를 획득하는 방법 및 장치. |
CN103700111A (zh) * | 2013-12-30 | 2014-04-02 | 深圳市开立科技有限公司 | 一种绘制图像标尺的方法和设备 |
CN103720489B (zh) * | 2013-12-30 | 2015-10-28 | 中国科学院深圳先进技术研究院 | 病变组织生长监测方法和系统 |
FR3017041B1 (fr) * | 2014-01-31 | 2016-03-04 | Centre Nat Rech Scient | Procede et dispositif ultrasonore de caracterisation des milieux mous anisotropes, et ensemble de sonde ultrasonore pour un tel dispositif de caracterisation |
CN106461766B (zh) * | 2014-05-30 | 2021-01-26 | 皇家飞利浦有限公司 | 来自多个声学窗口的同步相控阵列数据采集 |
KR101643622B1 (ko) * | 2014-09-25 | 2016-07-29 | 삼성전자주식회사 | 초음파 영상 처리 방법 및 이를 위한 초음파 영상 장치 |
JP6006769B2 (ja) * | 2014-10-16 | 2016-10-12 | 株式会社日立製作所 | 超音波診断装置 |
KR102406937B1 (ko) * | 2014-11-07 | 2022-06-10 | 삼성메디슨 주식회사 | 초음파 영상장치 및 그 제어 방법 |
CN107405133B (zh) | 2014-12-24 | 2021-06-08 | 声科影像有限公司 | 适用于各向异性介质成像的剪切波弹性成像方法及装置 |
JP6419945B2 (ja) * | 2015-03-23 | 2018-11-07 | 富士フイルム株式会社 | 音響波画像生成装置およびその制御方法 |
TWI580960B (zh) * | 2015-05-08 | 2017-05-01 | 國立臺灣大學 | 彈性分布影像生成系統 |
US10660604B2 (en) | 2015-07-13 | 2020-05-26 | Otonexus Medical Technologies, Inc. | Apparatus and method for characterization of acute otitis media |
RU2641067C1 (ru) * | 2015-08-03 | 2018-01-15 | Конинклейке Филипс Н.В. | Ультразвуковая система и способ измерения с использованием сдвиговой волны |
CN105212968B (zh) * | 2015-10-29 | 2019-01-04 | 无锡海斯凯尔医学技术有限公司 | 弹性检测方法和设备 |
JP6987496B2 (ja) * | 2015-12-04 | 2022-01-05 | キヤノンメディカルシステムズ株式会社 | 解析装置 |
KR102669618B1 (ko) * | 2016-04-12 | 2024-05-27 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 그 제어 방법 |
US10675001B2 (en) * | 2016-06-04 | 2020-06-09 | Otonexus Medical Technologies, Inc. | Apparatus and method for characterization of a ductile membrane, surface, and sub-surface properties |
KR20180054360A (ko) * | 2016-11-15 | 2018-05-24 | 삼성메디슨 주식회사 | 초음파 진단 장치 및 초음파 진단 장치 제어 방법 |
US10646202B2 (en) * | 2017-01-23 | 2020-05-12 | Siemens Medical Solutions Usa, Inc. | Sheer speed imaging using coherence |
KR102695456B1 (ko) | 2017-01-25 | 2024-08-14 | 삼성메디슨 주식회사 | 대상체에 관한 횡파 탄성 데이터를 표시하는 초음파 진단 장치 그 동작 방법 |
WO2019075697A1 (zh) | 2017-10-19 | 2019-04-25 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性测量装置及弹性对比测量方法 |
JP6996035B2 (ja) * | 2017-11-02 | 2022-01-17 | 富士フイルムヘルスケア株式会社 | 超音波診断装置、および、生体組織の物性評価方法 |
WO2019196033A1 (zh) * | 2018-04-11 | 2019-10-17 | 深圳迈瑞生物医疗电子股份有限公司 | 超声弹性成像方法和系统 |
CN110536644B (zh) * | 2018-04-28 | 2023-08-11 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声弹性测量装置及方法 |
US11452503B2 (en) * | 2018-05-18 | 2022-09-27 | Siemens Medical Solutions Usa, Inc. | Shear wave imaging based on ultrasound with increased pulse repetition frequency |
EP3818944A4 (en) * | 2018-07-02 | 2021-09-01 | FUJIFILM Corporation | ACOUSTIC DIAGNOSTIC DEVICE AND METHOD OF CONTROLLING AN ACOUSTIC DIAGNOSTIC DEVICE |
WO2020037674A1 (zh) * | 2018-08-24 | 2020-02-27 | 深圳迈瑞生物医疗电子股份有限公司 | 一种超声成像装置和方法、超声弹性检测装置和方法 |
WO2020047805A1 (zh) * | 2018-09-06 | 2020-03-12 | 深圳迈瑞生物医疗电子股份有限公司 | 一种剪切波弹性成像方法和装置、及计算机存储介质 |
CN109919918A (zh) * | 2019-02-21 | 2019-06-21 | 清华大学 | 弹性成像的控制方法及装置、计算机设备及可读存储介质 |
CN112438751A (zh) * | 2019-09-04 | 2021-03-05 | 通用电气精准医疗有限责任公司 | 用于剪切波弹性成像的方法和系统以及存储有对应程序的介质 |
CN113545806A (zh) * | 2020-04-26 | 2021-10-26 | 深圳迈瑞生物医疗电子股份有限公司 | 前列腺弹性成像方法和超声弹性成像系统 |
CN116642410B (zh) * | 2023-07-26 | 2023-10-10 | 西安交通大学 | 一种非接触式cfrp结构损伤监测系统及方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003210460A (ja) * | 2002-01-18 | 2003-07-29 | Chikayoshi Sumi | ずり弾性率計測装置および治療装置 |
JP2007105400A (ja) * | 2005-10-17 | 2007-04-26 | Toshiba Corp | 超音波診断装置及び画像処理装置 |
JP2009201989A (ja) * | 2008-02-27 | 2009-09-10 | Siemens Medical Solutions Usa Inc | ずり情報を決定するための方法および粘弾性組織プロパティを決定するためのシステム |
US20100016718A1 (en) * | 2008-07-16 | 2010-01-21 | Siemens Medical Solutions Usa, Inc. | Shear Wave Imaging |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4601413B2 (ja) * | 2004-12-21 | 2010-12-22 | 株式会社日立メディコ | 超音波診断装置 |
US7632231B2 (en) * | 2006-03-22 | 2009-12-15 | Wisconsin Alumni Research Foundation | Ultrasonic strain imaging device and method providing parallel displacement processing |
KR101014558B1 (ko) * | 2008-07-16 | 2011-02-16 | 주식회사 메디슨 | 스트레인 정규화를 통한 탄성영상 형성 방법 및 그를 위한초음파 시스템 |
US9078592B2 (en) * | 2008-10-27 | 2015-07-14 | Wisconsin Alumni Research Foundation | Ultrasonic strain imaging device with selectable cost-function |
US8992426B2 (en) * | 2009-05-04 | 2015-03-31 | Siemens Medical Solutions Usa, Inc. | Feedback in medical ultrasound imaging for high intensity focused ultrasound |
US20110060222A1 (en) * | 2009-09-10 | 2011-03-10 | Board Of Regents Of The University Of Texas System | Method and apparatus for characterizing inhomgeneities using axial shear strain elastography |
US8500639B2 (en) * | 2009-09-11 | 2013-08-06 | Mr Holdings (Hk) Limited | Systems and methods for shear wave field formation |
CN101912278A (zh) * | 2010-08-12 | 2010-12-15 | 陈庆武 | 超声动态弹性成像探头及方法 |
-
2011
- 2011-12-28 WO PCT/JP2011/080448 patent/WO2012105152A1/ja active Application Filing
- 2011-12-28 US US13/982,801 patent/US20130317361A1/en not_active Abandoned
- 2011-12-28 EP EP11857846.7A patent/EP2671511B1/en not_active Not-in-force
- 2011-12-28 JP JP2012555710A patent/JP5619191B2/ja active Active
- 2011-12-28 CN CN201180066724.0A patent/CN103347450B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003210460A (ja) * | 2002-01-18 | 2003-07-29 | Chikayoshi Sumi | ずり弾性率計測装置および治療装置 |
JP2007105400A (ja) * | 2005-10-17 | 2007-04-26 | Toshiba Corp | 超音波診断装置及び画像処理装置 |
JP2009201989A (ja) * | 2008-02-27 | 2009-09-10 | Siemens Medical Solutions Usa Inc | ずり情報を決定するための方法および粘弾性組織プロパティを決定するためのシステム |
US20100016718A1 (en) * | 2008-07-16 | 2010-01-21 | Siemens Medical Solutions Usa, Inc. | Shear Wave Imaging |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2939598A4 (en) * | 2012-12-25 | 2016-08-17 | Hitachi Aloka Medical Ltd | ULTRASONIC DIAGNOSTIC DEVICE AND ELASTICITY EVALUATION METHOD |
CN105491959B (zh) * | 2013-06-26 | 2019-06-25 | 皇家飞利浦有限公司 | 弹性成像测量系统和方法 |
RU2667617C2 (ru) * | 2013-06-26 | 2018-09-21 | Конинклейке Филипс Н.В. | Система и способ эластографических измерений |
JP6994494B2 (ja) | 2013-06-26 | 2022-01-14 | コーニンクレッカ フィリップス エヌ ヴェ | エラストグラフィ測定システム及びその方法 |
JP2020072937A (ja) * | 2013-06-26 | 2020-05-14 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | エラストグラフィ測定システム及びその方法 |
US10251627B2 (en) * | 2013-06-26 | 2019-04-09 | Koninklijke Philips N. V. | Elastography measurement system and method |
CN105491959A (zh) * | 2013-06-26 | 2016-04-13 | 皇家飞利浦有限公司 | 弹性成像测量系统和方法 |
US20160143621A1 (en) * | 2013-06-26 | 2016-05-26 | Koninklijke Philips N. V. | Elastography measurement system and method |
WO2014207668A1 (en) * | 2013-06-26 | 2014-12-31 | Koninklijke Philips N.V. | Elastography measurement system and method |
JP2015009040A (ja) * | 2013-07-01 | 2015-01-19 | 株式会社東芝 | 超音波診断装置及び超音波イメージングプログラム |
CN103454640A (zh) * | 2013-09-13 | 2013-12-18 | 北京东方惠尔图像技术有限公司 | 一种超声波成像装置和成像方法 |
JP2015058193A (ja) * | 2013-09-19 | 2015-03-30 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | 超音波診断装置 |
CN103462643A (zh) * | 2013-09-29 | 2013-12-25 | 深圳市开立科技有限公司 | 一种剪切波速度测量方法、装置和系统 |
JP2015084909A (ja) * | 2013-10-30 | 2015-05-07 | 株式会社東芝 | 超音波画像診断装置 |
EP2881041A1 (en) | 2013-11-08 | 2015-06-10 | Hitachi Aloka Medical, Ltd. | Apparatus and method for ultrasonic diagnosis |
CN104706377A (zh) * | 2013-12-13 | 2015-06-17 | 株式会社东芝 | 超声波诊断装置、图像处理装置以及图像处理方法 |
JP2015202373A (ja) * | 2014-04-16 | 2015-11-16 | 株式会社東芝 | 超音波診断装置及び制御プログラム |
US11337674B2 (en) | 2014-04-16 | 2022-05-24 | Canon Medical Systems Corporation | Ultrasonography apparatus and control method |
JP2018516135A (ja) * | 2015-06-04 | 2018-06-21 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | がん悪性度マップにより拡張された精密診断及び治療に対するシステム及び方法 |
KR20190075773A (ko) * | 2017-12-21 | 2019-07-01 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그 표시 방법 |
KR102618496B1 (ko) | 2017-12-21 | 2023-12-27 | 삼성메디슨 주식회사 | 초음파 영상 장치 및 그 표시 방법 |
US11995807B2 (en) | 2019-05-17 | 2024-05-28 | Supersonic Imagine | Method for ultrasound determination of a corrected image of a medium, and device for implementing this method |
Also Published As
Publication number | Publication date |
---|---|
JPWO2012105152A1 (ja) | 2014-07-03 |
EP2671511A4 (en) | 2015-11-18 |
CN103347450B (zh) | 2015-07-08 |
CN103347450A (zh) | 2013-10-09 |
US20130317361A1 (en) | 2013-11-28 |
JP5619191B2 (ja) | 2014-11-05 |
EP2671511B1 (en) | 2017-11-08 |
EP2671511A1 (en) | 2013-12-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5619191B2 (ja) | 超音波診断装置,および方法 | |
JP5926193B2 (ja) | 超音波診断装置 | |
EP2881041B1 (en) | Apparatus and method for ultrasonic diagnosis | |
US11717270B2 (en) | Method, apparatus and system for imaging in ultrasonic scanning | |
EP2926739B1 (en) | Acquisition control for elasticity ultrasound imaging | |
JP5456047B2 (ja) | 超音波診断装置 | |
Hasegawa et al. | High-frame-rate echocardiography using diverging transmit beams and parallel receive beamforming | |
Jensen | Medical ultrasound imaging | |
JP5646290B2 (ja) | 超音波診断装置及びその作動方法 | |
JP5511201B2 (ja) | 信号処理装置、超音波装置、信号処理装置の制御方法及び超音波装置の制御方法 | |
JP5798117B2 (ja) | 超音波診断装置及び超音波診断装置の作動方法 | |
KR102011545B1 (ko) | 초음파 진단 장치, 화상 처리 장치 및 화상 처리 방법 | |
JP2013523324A (ja) | 超音波撮像の方法および装置 | |
JP2006122683A (ja) | 超音波マルチライン画像歪みを削減するアパーチャ・シェージング推定手法 | |
JP2013005876A (ja) | 超音波診断装置 | |
JP2015512273A (ja) | 超音波イメージングのための方法及び装置 | |
JP2015213673A (ja) | 超音波診断装置 | |
WO2018082458A1 (zh) | 一种定量剪切波弹性成像方法及系统 | |
WO2007032134A1 (ja) | 超音波診断装置 | |
EP2157442B1 (en) | Formation of an elastic image in an ultrasound system | |
JP7387249B2 (ja) | 超音波診断装置、医用画像処理装置及び医用画像処理プログラム | |
KR102220822B1 (ko) | Arfi 이미징을 위한 교정 | |
US10080550B2 (en) | Ultrasonic apparatus and control method for the same | |
JP2008253663A (ja) | 超音波診断装置およびその制御処理プログラム | |
JP2023006198A (ja) | 超音波診断装置及びプログラム |
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: 11857846 Country of ref document: EP Kind code of ref document: A1 |
|
ENP | Entry into the national phase |
Ref document number: 2012555710 Country of ref document: JP Kind code of ref document: A |
|
REEP | Request for entry into the european phase |
Ref document number: 2011857846 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011857846 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13982801 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |