US9895137B2 - Ultrasonic diagnosis device and transmission/reception method - Google Patents
Ultrasonic diagnosis device and transmission/reception method Download PDFInfo
- Publication number
- US9895137B2 US9895137B2 US14/413,174 US201414413174A US9895137B2 US 9895137 B2 US9895137 B2 US 9895137B2 US 201414413174 A US201414413174 A US 201414413174A US 9895137 B2 US9895137 B2 US 9895137B2
- Authority
- US
- United States
- Prior art keywords
- detection
- displacement
- detection position
- shear wave
- ultrasonic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0833—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
- A61B8/085—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
-
- 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/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- 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/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- 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/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Definitions
- the present invention relates to an ultrasonic diagnostic apparatus, and more particularly, to an ultrasonic technology of transmitting and receiving an ultrasonic wave to detect a shear wave in a subject.
- an elastography technology As a method of diagnosing breast cancer, cirrhosis, vascular disease and the like, there has been known a method (an elastography technology) of diagnosing hardness in a subject such as a living body from an ultrasonic echo signal, instead of palpation of a doctor.
- an engaged person touches and presses an ultrasonic probe to a surface of the subject, thereby generating a displacement in a tissue in the subject.
- an ultrasonic beam for displacement detection is irradiated to acquire an echo signal
- a displacement in the compression direction is estimated on the basis of the acquired echo signal and a distortion, which is a spatial differentiating amount of the displacement, is calculated.
- a value relating to the hardness, for example, a Young's modulus is calculated from the distortion and stress.
- NPL 1 there is an elastography technology of using a focused beam of ultrasonic waves to apply a radiation pressure into a subject, displacing a target tissue while suppressing an influence of an interposed layer, irradiating an ultrasonic beam for displacement detection to acquire an echo signal, and diagnosing hardness of the subject such as a living body on the basis of the acquired echo signal.
- a target to be imaged is limited to an organ that can be easily pressed from a body surface.
- a slip surface which is an interposed layer, exists between the body surface and the liver, it is difficult to perform the pressing causing a sufficient displacement in the liver.
- the tissue is displaced by the ultrasonic waves, it is possible to reduce an influence of the interposed layer such as the slip surface, which reduces the dependency on the procedure.
- the pressing or radiation pressure is applied to generate a shear wave in the subject
- the displacement of the tissue in the subject which is accompanied by propagation of the shear wave
- is detected by the ultrasonic wave for displacement detection and a velocity of the shear wave and the like are estimated to estimate the hardness of the subject and the like.
- the estimation precision of the shear wave velocity is increased as the number of detection positions of the shear wave by the ultrasonic wave for displacement detection increases. Therefore, a method (a parallel tracking method) of repeating an operation of irradiating an ultrasonic beam pulse for displacement detection in parallel in a time-division manner to a plurality of detection positions so as to detect the displacement at the plurality of detection positions has been known.
- a time resolution of the shear wave measurement at the entire measurement positions is limited by PRT (Pulse Repetition Time) of the ultrasonic beam for displacement detection. That is, as the number of the detection positions, to which the ultrasonic beam for displacement detection is irradiated in parallel in the time-division manner, is increased, a period (a time interval) for which the ultrasonic beam for displacement detection is irradiated to the same detection position is prolonged and the time resolution per one position is lowered.
- PRT Pulse Repetition Time
- NPL 1 also discloses a method of enlarging a beam diameter of the ultrasonic beam for detecting the shear wave, irradiating the ultrasonic wave to one place and generating reception waveforms at a plurality of positions from a reception waveform of the irradiated ultrasonic wave, so as to detect the shear waves at the plurality of detection positions at the same time without lowering the time resolution by using the parallel tracking technology. In this case, it is possible to maintain the time resolution. However, since the beam diameter is enlarged, the spatial resolution is lowered.
- An object of the present invention is to solve the above problems occurring in the elastography technology and to provide an ultrasonic diagnostic apparatus capable of measuring a shear wave velocity with high time resolution and spatial resolution, and a method of transmitting and receiving an ultrasonic beam.
- an ultrasonic diagnostic apparatus including a ultrasonic probe; a displacement generation unit configured to displace an inside of a subject and to thus generate a shear wave; a detection unit configured to sequentially transmit a first ultrasonic beam to a plurality of detection positions of the subject from the ultrasonic probe, to receive reflection signals and to detect the displacement at the detection positions by using the received reflection signals, thereby detecting the shear wave, and a control unit configured to calculate a velocity of the shear wave on the basis of an output of the detection unit, thereby outputting hardness information of the subject.
- a method of transmitting and receiving an ultrasonic beam includes displacing an inside of a subject; transmitting a first ultrasonic beam to one of a plurality of detection positions of the subject from a ultrasonic probe and receiving a reflection signal, and switching the detection position on the basis of the received reflection signal and detecting a shear wave velocity caused due to the displacement at the plurality of detection positions, thereby calculating hardness information of the subject.
- the present invention it is possible to measure the shear wave with the high time resolution and spatial resolution.
- FIG. 1 is a block diagram showing an example of a system configuration of an ultrasonic diagnostic apparatus according to a first illustrative embodiment.
- FIG. 2 illustrates a state where a focused beam of ultrasonic waves is converged on a focus F and a shear wave is generated by an ultrasonic probe in the first illustrative embodiment.
- FIG. 3 illustrates a beam forming of an ultrasonic wave in the first illustrative embodiment.
- FIG. 4 illustrates a sequence of transmitting and receiving the ultrasonic beam in the first illustrative embodiment.
- FIG. 5A is a graph showing waveforms of time-variation of displacement resulting from propagation of a shear wave at each detection position in the first illustrative embodiment.
- FIG. 5B is a graph showing a relation between time at which the displacement of the shear wave peaks and the detection position in the first illustrative embodiment.
- FIG. 6 illustrates a sequence of detection timings of shear waves at detection positions X 1 , X 2 , X 3 in the related art.
- FIG. 7 illustrates a sequence of detection timings of the shear waves at the detection positions X 1 , X 2 , X 3 in the first illustrative embodiment.
- FIG. 8 is a flowchart showing a shear wave detection operation in first and second illustrative embodiments.
- FIG. 9 illustrates an example of a display screen of a hardness measurement result on a display unit 5 in the first illustrative embodiment.
- FIG. 10 illustrates a sequence of detection timings of the shear waves at the detection positions X 1 , X 2 , X 3 in an application embodiment of the first illustrative embodiment.
- FIG. 11A is a graph showing waveforms of time-variation of a particle velocity of the shear wave at each detection position in the first illustrative embodiment.
- FIG. 11B is a graph showing another example of waveforms of time-variation of the particle velocity of the shear wave at each detection position in the first illustrative embodiment.
- FIG. 12A illustrates a relation between a temporal waveform of the displacement of the shear wave and a threshold TH when a damping ratio is equal to or smaller than 1 in the first illustrative embodiment.
- FIG. 12B illustrates a relation between a temporal waveform of the displacement of the shear wave and a threshold TH, ⁇ TH when the damping ratio is larger than 1 in the first illustrative embodiment.
- FIG. 13 illustrates a sequence of detection timings of the shear waves at the detection positions X 1 , X 2 , X 3 in a second illustrative embodiment.
- FIG. 14A illustrates an example of the correlation function of the shear wave in the second illustrative embodiment.
- FIG. 14B illustrates an example of a correlation function of the shear wave in the second illustrative embodiment.
- FIG. 15 illustrates a propagation direction of the shear wave when there is a boundary in a third illustrative embodiment.
- FIG. 16 illustrates a sequence of detection timings of the shear waves at the detection positions X 1 , X 2 , X 3 in the third illustrative embodiment.
- FIG. 17 is a flowchart showing a shear wave detection operation in the third illustrative embodiments.
- FIG. 18 illustrates an example of a display screen of a measurement result on the display unit 5 in the third illustrative embodiment.
- FIG. 19A illustrates an example (a case where there is no boundary) of a subtrahend characteristic of the displacement (amplitude) of the shear wave in the third illustrative embodiment.
- FIG. 19B illustrates another example (a case where there is a boundary) of the subtrahend characteristic of the amplitude of the shear wave in the third illustrative embodiment.
- FIG. 20 is a block diagram showing an example of a system configuration of an ultrasonic diagnostic apparatus configured to generate a shear wave by using a mechanical driving source according to a fourth illustrative embodiment.
- the ultrasonic diagnostic apparatus of the illustrative embodiments is configured to generate a shear wave in a subject tissue, to detect displacement of the subject tissue caused by propagation of the shear wave at a plurality of detection positions and to estimate a shear wave velocity, for example.
- the shear wave velocity is calculated and estimated by a linear approximation of the detection position and arrival time of the shear wave.
- the arrival time of the shear wave is estimated by detecting a peak value of a temporal waveform of the shear wave, for example.
- FIG. 1 shows an overall configuration of an example of the ultrasonic diagnostic apparatus according to the first to third illustrative embodiments.
- the ultrasonic diagnostic apparatus is provided with an ultrasonic probe 1 configured to transmit and receive an ultrasonic beam to and from a subject such as a living body, a transmission and reception changeover switch 2 connected to the ultrasonic probe 1 , a displacement generation unit 10 configured to generate a displacement in the subject, a detection unit 20 configured to detect the displacement generated in the subject and a control unit 3 configured to control the displacement generation unit 10 and the detection unit 20 .
- the ultrasonic diagnostic apparatus is further provided with a digital scan converter 4 to which an output of the detection unit 20 is input, a display unit 5 configured to display the output, and an input unit 6 connected to the control unit 3 and configured to perform a variety of inputs, which will be described later.
- the displacement generation unit 10 includes a displacement generation transmission wave waveform generation unit 11 , a focus position setting unit 12 , a displacement generation transmission wave beam generation unit 13 and a beam frequency setting unit 14 , as described later.
- the detection unit 20 includes a displacement detection transmission wave waveform generation unit 21 , a displacement detection transmission wave beam generation unit 22 , a displacement detection reception wave beam calculation unit 23 and a wave detection unit 25 .
- an ultrasonic beam for displacement detection may be also referred to as a first ultrasonic beam and an ultrasonic beam for displacement generation may be also referred to as a second ultrasonic beam.
- the control unit 3 includes a waveform analysis unit 26 configured to receive an output of the wave detection unit 25 and to analyze a waveform thereof and a hardness information calculation unit 28 configured to receive an output of the waveform analysis unit 26 and to calculate hardness information of the subject. Also, the control unit 3 includes a central control unit 29 configured to control the waveform analysis unit 26 , the hardness information calculation unit 28 , the displacement generation unit 10 and the detection unit 20 .
- the ultrasonic probe 1 is connected to the displacement generation transmission wave beam generation unit 13 , the displacement detection transmission wave beam generation unit 22 and the displacement detection reception wave beam calculation unit 23 through the transmission and reception changeover switch 2 .
- an outline arrow indicates above noted connection relations between the ultrasonic probe 1 , the displacement generation transmission wave beam generation unit 13 , the displacement detection transmission wave beam generation unit 22 , the displacement detection reception wave beam calculation unit 23 and the transmission and reception changeover switch 2 .
- the displacement generation transmission wave beam generation unit 13 in the displacement generation unit 10 is configured to generate a transmission wave signal for each element of the ultrasonic probe 1 having delay time or weighting applied thereto so that an ultrasonic beam is focused on a position set by the focus position setting unit 12 , by using the waveform generated in the displacement generation transmission wave waveform generation unit 11 , under control of the central control unit 29 .
- the transmission wave signal from the displacement generation transmission wave beam generation unit 13 is converted into an ultrasonic signal at the ultrasonic probe 1 , which is then irradiated to the subject such as the living body.
- the second ultrasonic beam which is the ultrasonic beam for displacement generation, is irradiated.
- the irradiation starting time and ending time of the ultrasonic beam for displacement generation are set at the beam frequency setting unit 14 .
- a beam frequency indicates a repetition frequency of the irradiation of the ultrasonic beam for displacement generation.
- an ultrasonic beam for displacement detection which is the first ultrasonic beam for detecting a displacement of a tissue in the subject, is irradiated from the ultrasonic probe 1 .
- the displacement detection transmission wave beam generation unit 22 is configured to generate a transmission wave signal for each element of the ultrasonic probe 1 having delay time or weighting applied thereto so that an ultrasonic beam for displacement detection is focused on a desired position, by using the waveform generated in the displacement detection transmission wave waveform generation unit 21 , under control of the central control unit 29 .
- An echo signal which is a reflection signal reflected in the subject and returned to the ultrasonic probe 1 , is converted into an electric signal at the ultrasonic probe 1 , which is then sent to the displacement detection reception wave beam calculation unit 23 via the transmission and reception changeover switch 2 .
- An output of the displacement detection reception wave beam calculation unit 23 is output to the waveform analysis unit 26 in the control unit 3 , as a signal indicating a change of the subject tissue resulting from propagation of the shear wave, after signal processing in the wave detection unit 25 such as envelope detection, log compression, bandpass filter processing, gain control and the like.
- the waveform analysis unit 26 is configured to calculate a displacement of the subject tissue caused due to the shear wave and to analyze the shear wave, based on the input signal. Regarding the analysis of the shear wave, a calculation of detecting a peak value of a temporal waveform of the shear wave and calculating peak time is performed, for example.
- An output of the waveform analysis unit 26 is sent to the displacement detection transmission wave beam generation unit 22 via the central control unit 29 and is used to set a detection position and detection time of the shear wave. That is, in the ultrasonic diagnostic apparatus of this illustrative embodiment, the control is performed to switch a detection position of a plurality of detection positions, to which the ultrasonic beam for displacement detection (the first ultrasonic beam) is irradiated, based on the output of the detection unit 20 . Also, the detection time at the plurality of detection positions is controlled, based on the output of the detection unit 20 .
- the detection position and an analysis result are sent to the hardness information calculation unit 28 , as an output signal of the waveform analysis unit 26 .
- the analysis result is a signal indicating time at which the displacement of the subject tissue caused due to the shear wave peaks, for example.
- the hardness information calculation unit 28 is configured to calculate and output a shear wave velocity, a Young's modulus and the like indicating hardness information of the subject tissue, from the input signal.
- the signal indicating the hardness information, which is output from the hardness information calculation unit 28 is converted into an image signal at the digital scan converter 4 , which is then displayed on the display unit 5 , as a numerical value or image indicating the hardness information.
- a general computer including a central processing unit (CPU), a memory configured to store therein a program and data, and the like executes the program with the CPU, so that the central control unit 29 , the waveform analysis unit 26 , the hardness information calculation unit 28 and the like, which are parts of the block shown in FIG. 1 , can be implemented.
- the program is stored in the memory and the like of the computer configuration.
- at least one of the waveform analysis unit 26 and the hardness information calculation unit 28 can be configured by dedicated hardware.
- a first illustrative embodiment is an illustrative embodiment of the ultrasonic diagnostic apparatus capable of switching the detection time at a plurality of detection positions for measuring the shear wave in the subject to thus calculate a velocity of the shear wave, thereby estimating the hardness information of the subject.
- the velocity of the shear wave is calculated as the hardness information of the subject.
- the present invention is not limited thereto. That is, the other hardness information such as Young's modulus may be also used in the present invention.
- the ultrasonic diagnostic apparatus of this illustrative embodiment is provided with the ultrasonic probe 1 , the displacement generation unit 10 configured to displace an inside of a subject, the detection unit 20 configured to sequentially transmit a first ultrasonic beam to a plurality of detection positions of the subject from the ultrasonic probe 1 , to receive the reflection signals, to detect the displacement of the subject tissue at the detection positions by using the received reflection signals and to thus detect propagation of a shear wave, and the control unit 3 configured to calculate a shear wave velocity on the basis of an output of the detection unit 20 and to output the shear wave velocity as the hardness information of the subject.
- a method of transmitting and receiving an ultrasonic beam of this illustrative embodiment is an illustrative embodiment of a method of transmitting and receiving an ultrasonic beam including displacing an inside of a subject, transmitting a first ultrasonic beam to one of a plurality of detection positions of the subject from an ultrasonic probe to thus receive an echo signal (a reflection signal), and switching the detection position and detecting a shear wave velocity caused due to the displacement at the plurality of detection positions to thus calculate hardness information of the subject on the basis of the received reflection signal.
- the displacement generation unit 10 is configured to emit the ultrasonic beam for displacement generation, which is the second ultrasonic beam, and to thus displace the subject tissue, thereby generating the shear wave, before irradiating the first ultrasonic beam (ultrasonic beam for displacement detection).
- the detection position and detection time at which the first ultrasonic beam is irradiated so as to detect the shear wave by the detection unit 20 are changed by switching transmission time of the first ultrasonic beam to each of the plurality of detection positions.
- the linear array-type ultrasonic probe 1 is contacted to a body surface of a living body, which is the subject, and the second ultrasonic beam (the ultrasonic beam for displacement generation) irradiated from the ultrasonic probe 1 is focused on a desired cross-sectional plane in the living body on the basis of the output of the displacement generation transmission wave beam generation unit 13 of the displacement generation unit 10 .
- a propagation direction of the ultrasonic beam for displacement generation in the desired cross-sectional plane is perpendicular to the body surface.
- an outline arrow indicates a propagation direction of the shear wave and a solid arrow indicates a direction of the displacement.
- a beamforming of the ultrasonic beam for displacement generation is implemented by calculating a distance between a focus to which the ultrasonic beam is converged and each element 100 of the ultrasonic probe 1 and applying time, which is calculated by dividing the distance difference between the elements 100 by a sound velocity of a target, to the transmission wave signal for each element 100 , as the delay time.
- a radiation pressure is generated in correspondence to absorption or scattering of the ultrasonic wave associated with the propagation.
- the radiation pressure is the maximum at the focus and a displacement is generated in a living tissue in the focus area.
- FIG. 3 shows a case where the focus position is changed by controlling delay time of driving pulses for the elements 100 of the ultrasonic probe 1 .
- a focus position of the ultrasonic beam for displacement detection can be also changed by controlling the delay time of the transmission wave time to be supplied to each element 100 . Thereby, it is possible to change the detection position to which the ultrasonic beam for displacement detection is irradiated.
- FIG. 4 pictorially shows irradiation and reception sequences of a transmission wave beam 7 for displacement generation, which is the second ultrasonic beam, a transmission wave beam 8 for displacement detection, which is the first ultrasonic beam, and a reception wave beam 9 for displacement detection of the ultrasonic diagnostic apparatus of this illustrative embodiment.
- a transmission wave beam 7 for displacement generation which is the second ultrasonic beam
- a transmission wave beam 8 for displacement detection which is the first ultrasonic beam
- a reception wave beam 9 for displacement detection of the ultrasonic diagnostic apparatus of this illustrative embodiment.
- 1 in a vertical direction is referred to as ON and 0 is referred to as OFF.
- the transmission wave beam 7 for displacement generation and the transmission wave beam 8 for displacement detection are irradiated to the focus F and the detection position, respectively, when they are ON.
- the ON state of the reception wave beam 9 for displacement detection means that the connection of the displacement detection transmission wave beam generation unit 22 and the ultrasonic probe 1 is disconnected and the displacement detection reception wave beam calculation unit 23 is connected to the ultrasonic probe 1 by the transmission and reception changeover switch 2 , so that a reception wave signal, which is the echo signal (the reflection signal), is acquired and a phasing addition calculation is performed for beamforming.
- the switching of the ON and OFF states is performed by controlling the transmission and reception changeover switch 2 with a voltage amplitude value, for example. As shown in FIG.
- the ON state is first made in order of the transmission wave beam 8 for displacement detection and the reception wave beam 9 for displacement detection and a reference signal, which is used for the calculation of detecting the displacement of the subject tissue caused due to the shear wave and the shear wave velocity, is obtained.
- the transmission wave beam 7 for displacement generation which is the second ultrasonic beam, is irradiated to the focus F of FIG. 2 , so that a shear wave is generated in the subject tissue.
- the irradiation (ON) and non-irradiation (OFF) of the transmission wave beam 7 for displacement generation are repeated with a predetermined frequency PRF p .
- a value of the frequency PRF p is set in the beam frequency setting unit 14 of the displacement generation unit 10 .
- the measurement can be performed while controlling the ON and OFF frequency PRF p of the transmission wave beam 7 for displacement generation.
- a transmission frequency (a frequency of the transmission wave beam, not the repetition frequency PRF p ) of the transmission wave beam 7 for displacement generation is set in the vicinity of a central frequency at which the high-sensitive transmission and reception are realized by the ultrasonic probe 1 .
- the number of irradiation (ON) times of the transmission wave beam 7 for displacement generation is three times, for example.
- the number of irradiation times is not limited to three times.
- the transmission wave beam 7 for displacement generation is irradiated, the transmission wave beam 8 for displacement detection and the reception wave beam 9 for displacement detection become ON in order so as to detect the displacement of the subject tissue caused due to the shear wave at the detection position.
- a signal obtained by the transmission and reception beam for displacement detection is subject to the signal processing such as the bandpass filter and the calculation processing in the wave detection unit 25 of the detection unit 20 , and a signal equivalent to the frequency PRF p is extracted and output as a signal relating to the shear wave.
- the calculation of the displacement detection at the detection position is performed using the obtained reference signal and a signal obtained by the transmission and reception wave beams for displacement detection in the waveform analysis unit 26 .
- a correlation calculation, a phase detection and the like which are the well known technologies, are used.
- the transmission and reception wave beams for displacement detection are repeatedly made to be ON with a repetition frequency PRF d , so that a temporal waveform of the displacement caused due to the shear wave is detected.
- the frequency PRF d is set to a sampling interval, for example, ( 1/10) ⁇ T that can be sufficiently measured for a period T of the shear wave.
- the frequency PRF d is set in the displacement detection transmission wave waveform generation unit 21 of the detection unit 20 .
- the detection position of the displacement is a plurality of positions x 1 , x 2 , x 3 , which are located along the propagation direction of the shear wave parallel with the subject surface and are spaced at an equal interval.
- the hardness information calculation unit 28 of the control unit 3 is configured to calculate a gradient of an approximation straight line 50 , which is obtained by a linear approximation of the peak time and the detection position of the shear wave, from a graph of the peak time and the detection position, thereby calculating a shear wave velocity, which is a parameter indicating the hardness on the basis of output of the waveform analysis unit 26 .
- the shear wave velocity that is calculated using a peak value or zero cross value is a group velocity. That is, in this illustrative embodiment, the control unit 3 is configured to calculate a group velocity of the shear wave on the basis of the peak value or zero cross value of the shear wave, and to calculate the hardness information of the subject on the basis of the group velocity. As described later, a phase velocity and a viscosity of the shear wave are calculated on the basis of a spectrum value of the shear wave and the hardness information may be calculated on the basis of the phase velocity and the viscosity. An operator can set (select) the group velocity of the shear wave, and the phase velocity and the viscosity of the shear wave, as the hardness information to be obtained, by the input unit 6 .
- a conventional technology (parallel tracking technology) of detecting the displacement caused due to the shear wave at the plurality of detection positions x 1 , x 2 , x 3 is described with reference to FIG. 6 .
- this technology an operation of irradiating the ultrasonic beam for displacement detection in parallel in a time-division manner to the plurality of detection positions x 1 , x 2 , x 3 is repeated for an entire measurement time period, so that the displacement at the plurality of detection positions is detected at the same time for the entire measurement time period. That is, the operation of irradiating the ultrasonic beam for displacement detection in order of the respective detection positions x 1 , x 2 , x 3 , x 1 , x 2 , x 3 .
- a vertical axis in FIG. 6 indicates a detection position and displacement at the detection position, and a horizontal axis indicates the measurement time.
- a dotted line indicates a detection line (irradiation timing of the ultrasonic for displacement detection) of the shear wave in a time direction.
- curves indicating the shear waves W 1 , W 2 , W 3 indicate waveforms (displacement) of the shear waves arriving at the detection positions x 1 , x 2 , x 3 . Therefore, it can be seen that a displacement amount of the shear wave is detected at a position at which the detection line and the curves of the shear waves W 1 , W 2 , W 3 intersect as shown in FIG. 6 , for the entire measurement time period. As shown in FIG. 6 , the displacement of the shear wave is detected at the respective detection positions x 1 , x 2 , x 3 at an equal time interval PRT′.
- a time resolution of the shear wave displacement detection is described.
- the time resolution of the shear wave displacement detection is limited by a pulse repetition time (PRT) of the ultrasonic beam for displacement detection. That is, since it is not possible to detect the displacement of the shear wave with a period less than the PRT of the ultrasonic beam for displacement detection, the minimum time resolution is the PRT.
- the PRT of the ultrasonic beam for displacement detection is 1/PRFd, and when the detection position of the shear wave is only one, it is possible to irradiate the ultrasonic beam for displacement detection to the corresponding detection position with the PRT, so that the minimum time resolution is the same as the PRT.
- FIG. 4 the PRT of the ultrasonic beam for displacement detection is 1/PRFd
- the minimum time resolution PRT′ at each detection position is PRT* 3 .
- a detection sequence of the shear wave at the plurality of detection positions in this illustrative embodiment is described with reference to FIG. 7 .
- a measurement time period is allotted to each detection position by using a fact that there is a time difference in arrival time of the shear wave at the plurality of detection positions, and the shear wave is detected at only one detection position for one measurement time period.
- the switching of the detection position at which the shear wave is measured is controlled on the basis of the output of the waveform analysis unit 26 .
- the ultrasonic beam for displacement detection which is the first ultrasonic beam
- the ultrasonic beam for displacement detection is transmitted to only the detection position x 1 with the minimum period (PRT) during which the ultrasonic beam can be irradiated, so that the displacement of the shear wave W 1 is detected at the detection position x 1 (the measurement time period of the detection position x 1 ).
- the detection is not performed at the other detection positions x 2 , x 3 .
- the detection at the detection position x 1 is over, and at the same time, the ultrasonic beam for displacement detection, which is the first ultrasonic beam, is transmitted to the detection position x 2 with the minimum period (PRT) during which the ultrasonic beam can be irradiated, so that the detection of the shear wave W 2 at the detection position x 2 is performed (the measurement time period of the detection position x 2 ).
- PRT minimum period
- the ultrasonic beam for displacement detection which is the first ultrasonic beam
- the ultrasonic beam for displacement detection is transmitted to the detection position x 3 with the minimum period (PRT) during which the ultrasonic beam can be irradiated, so that the detection of the shear wave W 3 at the detection position x 3 is performed (the measurement time period of the detection position x 3 ).
- the peak value of the displacement of the shear wave is detected by the waveform analysis unit 26 of the control unit 3 , for example, the CPU by using the well-known technology such as HPF (High Pass Filter).
- HPF High Pass Filter
- the fact that there is a time difference in arrival time of the shear wave at the plurality of detection positions is used and thus the ultrasonic beam for displacement detection is continuously irradiated to only one detection position until the peak of the displacement of the shear wave at the one detection position is detected. Therefore, the time resolution PRT′ of the displacement detection of the shear wave is the same as the PRT.
- the irradiation position of the ultrasonic beam for displacement detection is moved to a next detection position at timing at which the peak is detected, so that it is possible to detect the peak of the displacement of the shear wave at each detection position.
- the measurement starts in step S 00 .
- an operator sets a hardness measuring range (ROI: Region of Interest) while seeing a cross-sectional image displayed on the display unit 5 , for example, a B mode image or strain image, and pushes a start switch in a GUI (Graphic User Interface) on a display screen of the display unit 5 or a start switch of the input unit 6 .
- ROI Hardness measuring range
- the range may be set through the display screen of the display unit 5 or the input unit 6 such as a track ball, a mouse, a keyboard, a touch panel and the like.
- the control unit 3 may read out a position corresponding to a measurement part such as the liver, the mammary gland and the like from the memory (not shown), which is a storage device, thereby setting the range.
- the start switch is arranged on a panel, which is the input unit 6 of the ultrasonic diagnostic apparatus.
- the detection positions x 1 , x 2 , x 3 . . . xn are set within the ROI set by the operator.
- the specific position of the detection position such as an interval of the detection positions can be automatically set by the control unit 3 or can be set through the input unit 6 by the operator.
- step S 10 the detection of the shear wave W 1 at the detection position x 1 shown in FIG. 7 starts.
- the displacement detection transmission wave beam generation unit 22 irradiates the ultrasonic beam for displacement detection, which is the first ultrasonic beam, to the detection position x 1 at a time interval of the PRT, the displacement detection reception wave beam calculation unit 23 receives a reflected wave from the subject, and the waveform analysis unit 26 calculates the displacement and detects a peak via the wave detection unit 25 .
- step S 12 of FIG. 8 the beam for displacement detection is continuously transmitted and received at the detection position x 1 until a peak value of the shear wave W 1 is detected at the detection position x 1 .
- a signal for setting a next displacement detection position to the detection position x 2 is output from the waveform analysis unit 26 to the displacement detection transmission wave beam generation unit 22 directly or via the central control unit 29 .
- the displacement detection transmission wave beam generation unit 22 of the detection unit 20 can switch the detection time at which the ultrasonic beam for displacement detection is transmitted to each of the detection positions x 1 , x 2 , x 3 . . . . In other words, it is possible to dynamically switch the detection position at which the displacement is detected, depending on the output of the detection unit 20 . Also, the first detection position x 1 and the time T 1 at which a peak value is made, which are shown in FIG. 7 , are recorded in the memory (not shown) by the central control unit 29 .
- the detection unit 20 can set the time at which the ultrasonic beam for displacement detection is transmitted to the second detection position x 2 , i.e., the time at which the detection position is dynamically switched to the second detection position x 2 , depending on an analysis result of the shear wave at the first detection position.
- step S 20 like step S 10 , the detection of the shear wave at the detection position x 2 starts. Also, in step S 22 , the same processing as step S 12 is executed for the shear wave W 2 detected at the detection position x 2 .
- the detection position x 2 and the time T 2 at which the peak value is made are recorded in the memory (not shown) by the control unit 3 .
- the shear wave velocity in the hardness information calculation unit 28 at least two detection positions and peak time are required, as can be clearly seen from FIG. 5B .
- the shear wave velocity is calculated in step S 50 .
- the shear wave velocity is calculated in the hardness information calculation unit 28 by using the peak time and detection position obtained from the waveform analysis unit 26 .
- a value calculated at the hardness information calculation unit 28 is the hardness information such as a shear viscosity, a Young's modulus, a shear wave elastic coefficient and the like, in addition to the shear wave velocity.
- step S 24 the detection position by the ultrasonic beam for displacement detection, which is transmitted from the displacement detection transmission wave beam generation unit 22 , is dynamically switched, so that the peak detection of the shear wave W 3 at the detection position x 3 starts.
- step S 40 when a peak cannot be detected at the detection position x 3 or the detection position x 3 reaches a limit of a measurable position, the shear wave velocity is estimated in step S 50 and the measurement is over in step S 60 .
- the description that a peak cannot be detected means that the phase or displacement amount obtained by the correlation calculation in the waveform analysis unit 26 is below a noise caused due to an influence of the pulse wave and the like.
- the limit of a measurable position is a position of the beam for displacement detection generated in the ultrasonic probe 1 or a pre-designated position. The designation of the position is made by an input of the operator through the input unit 6 or by the readout of the central control unit 29 of the control unit 3 from the memory (not shown).
- FIG. 9 illustrates an example of a display screen of a hardness measurement result measured in the hardness information calculation unit 28 of the apparatus of this illustrative embodiment and output to the display unit 5 .
- the display unit 5 displays the hardness information.
- the display unit 5 can display a degree of reliability of the hardness information and can display, as the degree of reliability, a standard deviation, a correlation function of waveforms of two or more shear waves, and the like. Further, it is possible to display the pulse repetition time (PRT) or pulse repetition frequency of the first ultrasonic beam.
- PRT pulse repetition time
- FIG. 9 illustrates a peak position 51 and a wave surface 52 of the shear wave when the shear wave is detected with the linear array-type probe.
- a value of the shear wave velocity estimated within the ROI, a value of the standard deviation indicating a degree of reliability of the value, and the like are displayed together.
- a graph showing a relation between the peak position and the detection position is displayed at the right of the display screen of FIG. 9 . From FIG. 9 , the operator can check with eyes whether the measurement is correctly performed.
- an effective PRT which is the time resolution of the shear wave measurement
- a PRF which is an inverse of the PRT
- a color mapping may be displayed by displaying a difference in the shear wave velocity with a color value on the basis of a plurality of measurement results.
- a color scaling function may be added to the configuration of the ultrasonic diagnostic apparatus shown in FIG. 1 .
- the shear wave velocity which is the hardness information of the subject calculated and output by the control unit 3
- the standard deviation which indicates the degree of reliability of the shear wave velocity
- the PRT the physical speed of the shear wave
- the PRF the information (ON, OFF) indicating whether the parallel tracking method is applied or not
- the wave surface 52 of the shear wave, the peak position 51 of the shear wave, and the like are displayed on the display screen of the display unit 5 of this illustrative embodiment, which can improve the convenience of the user.
- the control unit 3 can set the interval of the detection positions.
- the operator can set the interval of the detection positions through the input unit 6 .
- the interval of the detection positions is 1 ⁇ 2 ⁇ to ⁇ when a wavelength of the shear wave is denoted with ⁇ . The wavelength is determined from the shear wave velocity assumed at a measurement target part and the irradiation time information of the focused beam for displacement generation.
- the shear wave velocity in the liver is about 1 to 5 m/s.
- the wavelength can be calculated by multiplying the shear wave velocity and the period and is about 1 to 5 mm. Therefore, the control unit 3 sets an interval of the deposition positions of the shear wave within a range of 0.5 mm to 5 mm.
- the interval of the detection positions may be input by the operator through the display unit 5 capable of functioning as an input unit or the input unit 6 . That is, the display unit 5 and the input unit 6 shown in FIG. 1 can set the interval of the detection positions. In this case, the operator can manually flexibly change the interval of the detection positions, based on the measurement result. Alternatively, the interval of the detection positions can be automatically set. In this case, the setting of the interval of the detection positions is performed by a processing unit, which is the waveform analysis unit 26 of the control unit 3 . The interval of the detection positions may be changed to measure the hardness, and an average value, a variance value and the like of the hardness for the changed interval of the detection positions may be displayed on the display unit 5 .
- time at which the first ultrasonic beam is transmitted to a second detection position is set based on the detection result of the shear wave, which is the output of the detection unit 20 at a first detection position.
- the peak is detected at the detection position x 1 , and at the same time, the detection starts at the detection position x 2 .
- peak time T 2 at the detection position x 2 is estimated on the basis of the detection position x 1 , the peak time T 1 and irradiation time information of the focused beam for displacement generation.
- FIGS. 11A and 11B illustrate the particle velocity (Velocity) of the shear wave shown in FIG. 5A , i.e., the displacement of the shear wave/PRT.
- the particle velocity is different from the propagation velocity of the shear wave, which is a value to be measured.
- the particle velocity three parameters of a positive peak, a negative peak and a zero cross value exist per one detection position.
- one or more parameters may be used to estimate the group velocity of the shear wave, as the shear wave velocity.
- FIG. 11A illustrates time T 11 , T 12 , T 13 , T 14 at which positive peaks are made in particle velocity waveforms W 11 , W 12 , W 13 , W 14 of the shear waves.
- FIG. 11B illustrates time T 22 , T 23 , T 24 at which zero cross values are made.
- a merit of using the particle velocity is that when the positive peak is detected, the detectable numbers are increased, as compared to the configuration where the peak is detected from the displacement. Also, since the waveform of the shear wave is steeper, it is possible to detect the peak robustly against a variation due to a body motion and the like.
- the particle velocity and the particle acceleration may be detected using a filter configured to detect a gradient or outline, for example, gradient filter processing, Laplacian filter processing and the like, which are the well-known technologies as regards the signal processing of the temporal waveform of the shear wave.
- a filter configured to detect a gradient or outline, for example, gradient filter processing, Laplacian filter processing and the like, which are the well-known technologies as regards the signal processing of the temporal waveform of the shear wave.
- the zero cross value and the like may be used in addition to the peak value, as the value detected in the measurement flow of FIG. 8 .
- the ultrasonic diagnostic apparatus of the first illustrative embodiment it is possible to measure the shear wave with the high time resolution and spatial resolution. Thereby, it is possible to estimate the shear wave velocity and to evaluate the hardness with high precision. Also, since it is possible to suppress the time, which is required to detect the shear wave at each measurement position, to the minimum, it is possible to reduce the calculation cost. Further, since it is possible to determine the next detection position or time from the waveform of the shear wave observed at any one measurement position, it is possible to maximize the PRT per one place.
- the second illustrative embodiment is an illustrative embodiment of the ultrasonic diagnostic apparatus capable of evaluating the frequency dependency of the shear wave.
- the configuration of the ultrasonic diagnostic apparatus of the second illustrative embodiment is the same as the apparatus configuration of the first illustrative embodiment shown in FIG. 1 .
- the hardness information is measured using only the peak position and peak time of the shear wave.
- the shear wave velocity obtained in this case is a value relating to the group velocity.
- the frequency described here means a frequency of the shear wave.
- the shear wave velocity has the frequency dependency.
- an entire waveform of the shear wave is detected and spectrum-analyzed to calculate spectrum values.
- the FFT analysis which is the well-known technology, and the like are used.
- the detection unit 20 can calculate a phase velocity and a viscosity of the shear wave on the basis of the spectrum values of the shear wave and estimate the hardness information of the subject on the basis of the phase velocity and the viscosity.
- FIGS. 12A and 12B illustrate time waveforms of the shear wave resulting from a difference in the damping ratio.
- FIG. 12A illustrates a time response waveform when the damping ratio is equal to or smaller than 1
- FIG. 12B illustrates a time response waveform when the damping ratio is larger than 1.
- FIG. 13 illustrates a detection sequence of the shear wave according to the second illustrative embodiment.
- the detection sequence can be also implemented by the program processing of the CPU configuring the waveform analysis unit 26 , the hardness information calculation unit 28 and the central control unit 29 .
- the displacement measurement is continuously performed until the measurement of the entire waveform is over, even after time T 31 , T 32 , T 33 at which the peak values of the shear waves W 31 , W 32 , W 33 are detected.
- the measurement position is changed after the measurement of the entire waveform is over. As shown in FIGS.
- the ending time of the entire waveform is set to time at which an absolute value of the displacement becomes smaller than an absolute value
- the threshold a detectable minimum displacement amount, i.e., about 0.05 ⁇ m is used, for example.
- the displacement of the shear wave becomes larger than the threshold TH at time Tth( 1 ), increases over time, becomes maximum and then decreases, and becomes smaller than the threshold TH at time Tth( 2 ).
- time Tth( 2 ) is set as the measurement ending time.
- the damping ratio is larger than 1, and the displacement repeatedly becomes negative and positive over time.
- the time becoming a candidate of the measurement ending is time Tth( 2 ), Tth( 4 ) and Tth( 6 ).
- the earliest time of the time candidates is used as the ending time.
- the latest time Tth( 6 ) is used.
- the ending time can be determined in the same manner.
- the measurement target is the particle velocity, the particle acceleration, the phase or the phase velocity
- the threshold to be set is changed in conformity to the parameter to be measured.
- time zones (TD 2 , TD 3 ) during which the measurement is performed at the same time at two detection positions. That is, there is a time zone during which the measurement time periods at the two detection positions overlap.
- the parallel tracking method (refer to NPL 1) that is the well known technology is applied only for the time zones (TD 2 , TD 3 ) during which the measurement is performed at the same time at two detection positions.
- PRT′ is PRT* 2 .
- the parallel tracking method is applied at two detection positions only for the specific measurement time zones (TD 2 , TD 3 ), as shown in FIG. 13 .
- the interval of the detection positions may be manually changed by the operator through the input unit 6 or may be automatically set by the apparatus.
- a measurement flow of the shear wave of the second illustrative embodiment is the same as the measurement flow of the first illustrative embodiment shown in FIG. 8 .
- the time at which the measurement is over with the respective detection values is when the absolute value of the displacement amount is below the threshold TH.
- the shear wave is spectrum-analyzed to estimate the frequency characteristic of the shear wave velocity and the viscosity parameter such as the shear viscosity.
- a measurement method of measuring the group velocity with the method of the first illustrative embodiment at first time and measuring the viscosity parameter at second time may be performed.
- the value displayed on the display unit 5 in step S 60 is the viscosity parameter such as the shear viscosity, in addition to the values described in the first illustrative embodiment.
- a graph showing a relation between the shear wave velocity and the frequency may be displayed.
- ‘Parallel: ON’ may be displayed on the display unit 5 in the display example of FIG. 9 , and the operator can visually check the detection method of the shear wave.
- the degree of reliability may be determined from a correlation function of the waveforms of two or more shear waves detected at different positions, in addition to the standard deviation described in the first illustrative embodiment with reference to FIG. 9 .
- the waveforms of the shear waves at the different detection positions are similar to each other.
- the correlation function of the waveforms of the two or more shear waves which is a value indicating a degree of similarity of the waveforms of two or more shear waves, may be used.
- the correlation function is described with reference to FIGS. 14A and 14B .
- the correlation function can be calculated using temporal waveforms of two different shear waves or waveforms of spectrum values thereof.
- FIG. 14A illustrates a graph of a correlation function, which is calculated using temporal waveforms W 11 , W 12 of the shear waves, and a maximum value A and a variance value ⁇ of the correlation function.
- FIG. 14B illustrates a graph of a correlation function, which is calculated using frequency spectra W 11 ′, W 12 ′ of the displacement of the shear waves, and a maximum value A′ and a variance value ⁇ ′ of the correlation function.
- the maximum value or variance value or both values are displayed on the display unit 5 , so that the operator can evaluate the degree of reliability of the measurement result.
- the waveform that is used for calculation of the correlation function the waveforms of the particle velocity, the phase rotation amount, the phase, the acceleration and the like of the shear wave may be also used, in addition to the temporal waveform of the shear wave.
- the ultrasonic diagnostic apparatus of the above-described second illustrative embodiment it is possible to accomplish the effects capable of evaluating the frequency dependency of the shear wave and also evaluating the degree of reliability of the measurement result, in addition to the effects of the apparatus of the first illustrative embodiment.
- the third illustrative embodiment is an illustrative embodiment of the ultrasonic diagnostic apparatus configured to detect hardness of the subject, considering a reflected wave of the shear wave, too. That is, in the third illustrative embodiment, the detection unit 20 of the ultrasonic diagnostic apparatus is an embodiment of the configuration of detecting a traveling wave or reflected wave of the shear wave. In the meantime, although the configuration of the ultrasonic diagnostic apparatus of the third illustrative embodiment is the same as the apparatus configuration of the first illustrative embodiment shown in FIG. 1 , the detection unit 20 can detect a traveling wave or reflected wave of the shear wave.
- the first and second illustrative embodiments only one direction is considered as regards the propagation direction of the shear wave and the reflection of the shear wave is not considered.
- a measurement method relating to a case where the shear wave is reflected is described.
- a case where only the peak value is detected is described.
- a case where the viscosity is measured can be also considered in the same manner.
- FIG. 15 illustrates a vertical incidence, i.e., a case where a propagation direction of the shear wave is perpendicular to the boundary of the media.
- the shear wave before passing through a boundary 53 is referred to as an incident wave.
- the density of the media 1 , 2 is about 1000 kg/m 3 when they are assumed as soft parts of the living body.
- the shear wave is incident to the boundary.
- the incident wave is divided into a transmitted wave and a reflected wave. Due to the shear wave velocities c 1 , c 2 in the media 1 , 2 , an amplitude and a phase of the reflected wave of the shear wave are different.
- FIG. 16 illustrates a detection sequence of the reflected wave of the shear wave when c 1 >>c 2 .
- the boundary exists between the detection positions x 2 , x 3 .
- a peak is detected at the detection position x 2 and the peak detection of the shear wave starts at the detection position x 3 starts.
- time T 43 at which the shear wave peaks at the detection position x 3 is predicted from the detection positions x 1 , x 2 and peak time T 41 , T 42 at the respective positions.
- the shear wave velocity at the detection position x 3 is slower than the shear wave velocity at the detection position x 2 .
- ⁇ T is calculated by the central control unit 29 by reading a value of the shear wave velocity from the storage medium (not shown), depending on the measurement target.
- the detection position is changed to a position, for example, a position xr located in an opposite direction to the propagation direction.
- a detection flow of the shear wave according to the third illustrative embodiment is described with reference to FIG. 17 .
- the sequence can be also implemented by the program processing of the CPU configuring the waveform analysis unit 26 , the hardness information calculation unit 28 and the central control unit 29 .
- steps S 42 , S 44 , S 46 , S 48 , S 52 , S 62 different from the detection flow of FIG. 8 described in the first illustrative embodiment are described.
- a position corresponding to r is output to the displacement detection transmission wave beam generation unit 22 directly from the waveform analysis unit 26 or through the central control unit 29 .
- the position xr x 2 ⁇ 2 ⁇ (x 2 ⁇ x 1 ).
- an optimal value is set depending on the detection position in step S 30 .
- step S 46 the detection of the reflected wave of the shear wave starts at the detection position xr.
- step S 48 a peak value of the shear wave is detected at the detection position xr, and peak time T 44 is output to the hardness information calculation unit 28 .
- step S 52 a shear wave velocity, a coefficient of reflection and the like are calculated using the detection position, the peak value and the peak time.
- a reflection position is also estimated at the same time.
- the reflection position can be estimated using the detection positions x 2 , xr and time T 42 , T 44 at which the peak values are made at the corresponding positions.
- the calculation is performed in the hardness information calculation unit 28 .
- a measurement result is displayed.
- FIG. 18 illustrates an example of the display screen of the display unit 5 in the third illustrative embodiment.
- the coefficient of reflection is displayed on the display unit 5 , in addition to the shear wave velocity.
- a wave surface 55 of the reflected wave of the shear wave and a boundary position 54 are displayed.
- a graph showing positions of the traveling wave and the reflected wave may be displayed, as the peak position relative to the detection position.
- the operator can check whether the measurement is correctly performed while seeing the display unit 5 . Also, the operator can specify a site of lesion and observe a symptom from the information of the boundary position, the coefficient of reflection and the like.
- an output value of the hardness information calculation unit 28 is the coefficient of reflection of 0, for example.
- the measurement position of the reflected wave is one place.
- the number of the measurement positions may be two or more places.
- the entire wave surface may be measured to calculate the frequency characteristic and the viscosity parameter.
- a detection method of the reflected wave is described when c 1 ⁇ c 2 .
- the threshold since there is a transmitted wave as well as the reflected wave, it is not possible to evaluate whether the reflection is made, by using the threshold. In this case, it is determined whether there is a reflected wave from the damping characteristic of the displacement at the detection positions x 1 , x 2 , x 3 .
- the damping characteristic of the shear wave is described with reference to FIGS. 19A, 19B .
- the maximum values P 1 , P 2 , P 3 of the displacement of the shear waves at the respective detection positions and time at which the maximum values are made are expressed by a damping function.
- the maximum displacement P is expressed by B*exp( ⁇ t). ⁇ is a damping integer and t is time.
- c 1 ⁇ c 2 is described with reference to FIG. 19B .
- the shear wave W 43 is a displacement waveform after passing through the boundary
- the maximum displacement P 3 ′ is a value smaller than a value when there is no boundary.
- the shear wave W 43 is regarded as a transmitted wave.
- the threshold 10% is used.
- the damping function may be a function showing a general damping of a wave, such as Hankel function, in addition to the exponential function.
- a sign of the detected peak of the shear wave is the same as those of the peaks detected at the detection positions x 1 , x 2 .
- the operator can select whether or not to detect the reflected wave by using the input unit 6 .
- the parallel tracking method described in the second illustrative embodiment it is possible to detect both the transmitted wave and the reflected wave, irrespective of the coefficient of reflection.
- this illustrative embodiment can be also applied to a general case where the propagation direction of the shear wave is not perpendicular to the boundary. According to the ultrasonic diagnostic apparatus of the above-described third illustrative embodiment, it is possible to accomplish the effect capable of evaluating the reflected wave of the shear wave, in addition to the effects of the apparatus of the first illustrative embodiment.
- the displacement generation unit 10 is configured to emit the second ultrasonic beam, thereby displacing the subject. That is, the radiation pressure is generated by the focused beam of ultrasonic waves, so that the shear wave is generated.
- the well known methods such as a mechanical driving (a DC motor, a vibration pump and the like), a manual pressing, a pressing by an electric pulse, and the like may be also used, in addition to the radiation pressure.
- a mechanical driving configured to displace the subject by a mechanical vibration is used to touch and press the ultrasonic probe to a surface of the subject such as the living body.
- a mechanically-driven displacement generation unit is an illustrative embodiment of the ultrasonic diagnostic apparatus configured to displace the subject by the mechanical driving.
- FIG. 20 illustrates a configuration example of the ultrasonic diagnostic apparatus configured to generate the shear wave by using the mechanical driving for measurement according to the fourth illustrative embodiment. Here, only differences from the FIG. 1 are described.
- the displacement is generated from a mechanically-driven displacement generation unit 30 , instead of the displacement generation unit 10 .
- a waveform used for the mechanical driving is generated by a displacement generation transmission wave waveform generation unit 31 and is then output to a displacement generation driving unit 33 .
- a driving frequency is set in a driving frequency setting unit 34 .
- the displacement generation driving unit 33 is configured to drive an actuator (not shown) such as a motor. The actuator is touched to a body surface of the living body, thereby generating a displacement and a shear wave in the living body.
- a merit of using the mechanical driving is that since it is possible to ignore a temperature increase in the living body due to the mechanical driving, there are no limitations on the measurement time and the number of measurement times. Also, when the actuator is operated by a sinusoidal wave, it is possible to measure a viscosity parameter of a specific frequency component. However, the spatial resolution is lower, as compare to the configuration of the first to third illustrative embodiments where the radiation pressure is used. On the other hand, when the radiation pressure is used, the measurement should be performed at the measurement time and measurement interval at which the temperature increase in the living body does not exceed a limit value. However, the spatial resolution is higher because the focused beam can be irradiated to any position.
- the present invention is not limited to the above-described illustrative embodiments and includes a variety of modified embodiments.
- the above-described illustrative embodiments have been described in detail for better understandings and it is not construed to include all the described configurations.
- a part of the configuration of any illustrative embodiment can be replaced with the configuration of the other illustrative embodiments, and the configuration of the other illustrative embodiments can be added to the configuration of any illustrative embodiment.
- a part of the configuration of each illustrative embodiment can be added/deleted/replaced with the other configurations.
- the above-described illustrative embodiments relate to the ON/OFF control of the transmission and reception, the processing method of the received signal, and the feedback of the transmission and reception control based on the results thereof. That is, when implementing the illustrative embodiments, a convex-type or sector-type probe may be used instead of the linear array-type probe, a two dimensional probe may be used and a type of a device to be used is not particularly limited inasmuch as it can transmit and receive the ultrasonic wave.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Vascular Medicine (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physiology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
-
- 1: ultrasonic probe
- 2: transmission and reception changeover switch
- 3: control unit
- 4: digital scan converter
- 5: display unit
- 6: input unit
- 7: transmission wave beam for displacement generation
- 8: transmission wave beam for displacement detection
- 9: reception wave beam for displacement detection
- 10: displacement generation unit
- 11: displacement generation transmission wave waveform generation unit
- 12: focus position setting unit
- 13: displacement generation transmission wave beam generation unit
- 14: beam frequency setting unit
- 20: detection unit
- 21: displacement detection transmission wave waveform generation unit
- 22: displacement detection transmission wave beam generation unit
- 23: displacement detection reception wave beam calculation unit
- 25: wave detection unit
- 26: waveform analysis unit
- 28: hardness information calculation unit
- 29: central control unit
- 30: mechanically-driven displacement generation unit
- 31: displacement generation transmission wave waveform generation unit
- 33: displacement generation driving unit
- 34: driving frequency setting unit
- 50: approximation straight line
- 51: peak position
- 52: wave surface of shear wave
- 53: boundary
- 54: boundary position
- 55: wave surface of reflected wave of shear wave
- 100: element of ultrasonic probe
Claims (20)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2013042991 | 2013-03-05 | ||
JP2013-042991 | 2013-03-05 | ||
PCT/JP2014/051950 WO2014136502A1 (en) | 2013-03-05 | 2014-01-29 | Ultrasonic diagnosis device and transmission/reception method |
Publications (2)
Publication Number | Publication Date |
---|---|
US20150173718A1 US20150173718A1 (en) | 2015-06-25 |
US9895137B2 true US9895137B2 (en) | 2018-02-20 |
Family
ID=51491028
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/413,174 Active US9895137B2 (en) | 2013-03-05 | 2014-01-29 | Ultrasonic diagnosis device and transmission/reception method |
Country Status (5)
Country | Link |
---|---|
US (1) | US9895137B2 (en) |
EP (1) | EP2965694A4 (en) |
JP (1) | JP5753633B2 (en) |
CN (1) | CN104302233B (en) |
WO (1) | WO2014136502A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10327738B2 (en) | 2015-03-31 | 2019-06-25 | Samsung Medison Co., Ltd. | Ultrasound imaging apparatus and method of processing ultrasound image thereof |
US10485515B2 (en) | 2015-01-28 | 2019-11-26 | Konica Minolta, Inc. | Ultrasound diagnostic device and method for controlling ultrasound diagnostic device |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015029651A1 (en) * | 2013-08-26 | 2015-03-05 | 日立アロカメディカル株式会社 | Diagnostic ultrasound apparatus and elasticity evaluation method |
JP5730978B2 (en) * | 2013-11-08 | 2015-06-10 | 日立アロカメディカル株式会社 | Ultrasonic diagnostic apparatus and method |
JP6475456B2 (en) * | 2014-09-26 | 2019-02-27 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ultrasonic diagnostic equipment |
CA2965755A1 (en) | 2014-10-28 | 2016-05-06 | Super Sonic Imagine | Imaging methods and apparatuses for performing shear wave elastography imaging |
JP6224572B2 (en) * | 2014-11-26 | 2017-11-01 | ジーイー・メディカル・システムズ・グローバル・テクノロジー・カンパニー・エルエルシー | Ultrasonic diagnostic apparatus and program |
KR101868019B1 (en) * | 2015-03-31 | 2018-06-18 | 삼성전자주식회사 | Ultrasound Imaging Apparatus and Method for processing a ultrasound image thereof |
JP6390516B2 (en) * | 2015-05-27 | 2018-09-19 | コニカミノルタ株式会社 | Ultrasonic diagnostic apparatus and control method of ultrasonic diagnostic apparatus |
CN105212968B (en) * | 2015-10-29 | 2019-01-04 | 无锡海斯凯尔医学技术有限公司 | Elastomeric check method and apparatus |
JP6987496B2 (en) * | 2015-12-04 | 2022-01-05 | キヤノンメディカルシステムズ株式会社 | Analyst |
US11071524B2 (en) | 2015-12-04 | 2021-07-27 | Canon Medical Systems Corporation | Analyzing apparatus |
US11006928B2 (en) * | 2016-02-10 | 2021-05-18 | Siemens Medical Solutions Usa, Inc. | Sound speed imaging using shear waves |
WO2017169332A1 (en) * | 2016-03-28 | 2017-10-05 | 株式会社日立製作所 | Ultrasonic imaging device and ultrasonic reception signal processing method |
JP6843159B2 (en) * | 2016-06-10 | 2021-03-17 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | Use of reflected shear waves to monitor damage growth in thermal ablation |
JP6731321B2 (en) * | 2016-09-06 | 2020-07-29 | 株式会社日立製作所 | Ultrasonic diagnostic equipment |
US20190257931A1 (en) * | 2016-09-15 | 2019-08-22 | National Institute Of Maritime, Port And Aviation Technology | Ultrasonic three-dimensional measurement apparatus |
KR101931748B1 (en) * | 2017-03-28 | 2019-03-13 | 삼성메디슨 주식회사 | Ultrasound diagnostic apparatus and method for operating the same |
WO2018226688A1 (en) * | 2017-06-05 | 2018-12-13 | Mayo Foundation For Medical Education And Research | Estimating phase velocity dispersion in ultrasound elastography using a multiple signal classification |
CN107440740B (en) * | 2017-07-21 | 2021-06-25 | 无锡海斯凯尔医学技术有限公司 | Method and device for determining viscoelasticity of medium |
CN107505232B (en) * | 2017-07-21 | 2019-09-03 | 无锡海斯凯尔医学技术有限公司 | Motion information acquisition methods and device |
US11766244B2 (en) * | 2018-05-03 | 2023-09-26 | Koninklijke Philips N.V. | Shear wave amplitude reconstruction for tissue elasticity monitoring and display |
CN109077754B (en) * | 2018-07-06 | 2020-04-14 | 深圳大学 | Method and equipment for measuring tissue mechanical characteristic parameters |
CN110545732B (en) * | 2018-09-25 | 2022-06-03 | 深圳迈瑞生物医疗电子股份有限公司 | Method and device for elastic imaging and storage medium |
JP7236312B2 (en) * | 2019-04-04 | 2023-03-09 | 富士フイルムヘルスケア株式会社 | Ultrasound diagnostic device, signal processing device, and program |
CN110840488B (en) * | 2019-11-27 | 2021-09-28 | 深圳开立生物医疗科技股份有限公司 | Imaging method, system and device based on shear wave |
KR102617857B1 (en) | 2020-03-12 | 2023-12-22 | 지멘스 메디컬 솔루션즈 유에스에이, 인크. | Liver disease activity estimation with ultrasound medical imaging |
CN116528769A (en) * | 2021-01-04 | 2023-08-01 | 深圳迈瑞生物医疗电子股份有限公司 | Viscoelastic measurement method and ultrasonic imaging system |
US20220354464A1 (en) * | 2021-04-23 | 2022-11-10 | University Of Rochester | Reverberant shear wave gradients for ultrasound elastography |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06319736A (en) | 1993-05-17 | 1994-11-22 | Toshiba Corp | Ultrasonic doppler diagnostic apparatus |
US5999486A (en) * | 1998-07-23 | 1999-12-07 | Colorado School Of Mines | Method for fracture detection using multicomponent seismic data |
US6023977A (en) * | 1997-08-01 | 2000-02-15 | Acuson Corporation | Ultrasonic imaging aberration correction system and method |
US20050004463A1 (en) * | 2003-04-09 | 2005-01-06 | Shigao Chen | Method and apparatus for shear property characterization from resonance induced by oscillatory radiation force |
US20060173319A1 (en) * | 2005-01-21 | 2006-08-03 | Chikayoshi Sumi | Clinical apparatuses |
WO2007083745A1 (en) | 2006-01-20 | 2007-07-26 | Hitachi Medical Corporation | Elastic image display method and elastic image display |
US20090216119A1 (en) | 2008-02-27 | 2009-08-27 | Liexiang Fan | Sparse tissue property measurements in medical ultrasound imaging |
US20100016718A1 (en) | 2008-07-16 | 2010-01-21 | Siemens Medical Solutions Usa, Inc. | Shear Wave Imaging |
US20100069751A1 (en) | 2008-09-18 | 2010-03-18 | General Electric Company | Systems and methods for detecting regions of altered stiffness |
US20100280373A1 (en) * | 2009-05-04 | 2010-11-04 | Liexiang Fan | Feedback in medical ultrasound imaging for high intensity focused ultrasound |
WO2011004661A1 (en) * | 2009-07-07 | 2011-01-13 | 株式会社 日立メディコ | Ultrasonic diagnosis apparatus and ultrasonic measurement method |
US20110028838A1 (en) | 2009-07-31 | 2011-02-03 | Super Sonic Imagine | Method and Apparatus for Measuring a Physical Parameter in Mammal Soft Tissues by Propagating Shear Waves |
WO2011027644A1 (en) * | 2009-09-04 | 2011-03-10 | 株式会社 日立メディコ | Ultrasonic diagnostic device |
US20110066030A1 (en) * | 2009-09-11 | 2011-03-17 | Mr Holdings (Hk) Ltd. | Systems and methods for shear wave field formation |
WO2011064688A1 (en) | 2009-11-25 | 2011-06-03 | Koninklijke Philips Electronics N.V. | Ultrasonic shear wave imaging with focused scanline beamforming |
US20110184287A1 (en) * | 2009-12-10 | 2011-07-28 | Mcaleavey Stephen A | Methods And Systems For Spatially Modulated Ultrasound Radiation Force Imaging |
US20110216311A1 (en) * | 2010-03-02 | 2011-09-08 | Li-Cor, Inc. | Method and apparatus for locking a laser with a resonant cavity |
JP2012100997A (en) | 2010-11-12 | 2012-05-31 | Hitachi Medical Corp | Ultrasonograph and operation method therefor |
WO2012077579A1 (en) | 2010-12-08 | 2012-06-14 | 株式会社日立メディコ | Ultrasound diagnosis apparatus |
US20120215101A1 (en) * | 2011-02-17 | 2012-08-23 | Siemens Medical Solutions Usa, Inc | Viscoelasticity Measurement Using Amplitude-Phase Modulated Ultrasound Wave |
US20120232387A1 (en) | 2011-03-08 | 2012-09-13 | Fujifilm Corporation | Ultrasound diagnostic apparatus and method of determining elasticity index reliability |
US20140221833A1 (en) * | 2011-04-08 | 2014-08-07 | Canon Kabushiki Kaisha | Subject information acquisition apparatus |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2913875B1 (en) * | 2007-03-21 | 2009-08-07 | Echosens Sa | DEVICE FOR MEASURING VISCOELASTIC PROPERTIES OF BIOLOGICAL TISSUES AND METHOD USING THE DEVICE |
CN101784234B (en) * | 2007-05-16 | 2014-07-09 | 超声成像公司 | Method and device for measuring a mean value of visco-elasticity of a region of interest |
CN102078205A (en) * | 2011-03-04 | 2011-06-01 | 深圳市一体医疗科技股份有限公司 | Displacement estimating method for measuring elasticity of viscoelastic medium and application method |
-
2014
- 2014-01-29 CN CN201480000976.7A patent/CN104302233B/en not_active Expired - Fee Related
- 2014-01-29 JP JP2014545027A patent/JP5753633B2/en not_active Expired - Fee Related
- 2014-01-29 EP EP14759767.8A patent/EP2965694A4/en not_active Withdrawn
- 2014-01-29 WO PCT/JP2014/051950 patent/WO2014136502A1/en active Application Filing
- 2014-01-29 US US14/413,174 patent/US9895137B2/en active Active
Patent Citations (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH06319736A (en) | 1993-05-17 | 1994-11-22 | Toshiba Corp | Ultrasonic doppler diagnostic apparatus |
US6023977A (en) * | 1997-08-01 | 2000-02-15 | Acuson Corporation | Ultrasonic imaging aberration correction system and method |
US5999486A (en) * | 1998-07-23 | 1999-12-07 | Colorado School Of Mines | Method for fracture detection using multicomponent seismic data |
US20050004463A1 (en) * | 2003-04-09 | 2005-01-06 | Shigao Chen | Method and apparatus for shear property characterization from resonance induced by oscillatory radiation force |
US20060173319A1 (en) * | 2005-01-21 | 2006-08-03 | Chikayoshi Sumi | Clinical apparatuses |
WO2007083745A1 (en) | 2006-01-20 | 2007-07-26 | Hitachi Medical Corporation | Elastic image display method and elastic image display |
US20100220901A1 (en) | 2006-01-20 | 2010-09-02 | Takeshi Matsumura | Elastic Image Display Method and Elastic Image Display Device |
US20090216119A1 (en) | 2008-02-27 | 2009-08-27 | Liexiang Fan | Sparse tissue property measurements in medical ultrasound imaging |
US20100016718A1 (en) | 2008-07-16 | 2010-01-21 | Siemens Medical Solutions Usa, Inc. | Shear Wave Imaging |
US20100069751A1 (en) | 2008-09-18 | 2010-03-18 | General Electric Company | Systems and methods for detecting regions of altered stiffness |
US20100280373A1 (en) * | 2009-05-04 | 2010-11-04 | Liexiang Fan | Feedback in medical ultrasound imaging for high intensity focused ultrasound |
WO2011004661A1 (en) * | 2009-07-07 | 2011-01-13 | 株式会社 日立メディコ | Ultrasonic diagnosis apparatus and ultrasonic measurement method |
US20120123263A1 (en) * | 2009-07-07 | 2012-05-17 | Hitachi Medical Corporation | Ultrasonic diagnostic apparatus and ultrasonic measurement method |
US20110028838A1 (en) | 2009-07-31 | 2011-02-03 | Super Sonic Imagine | Method and Apparatus for Measuring a Physical Parameter in Mammal Soft Tissues by Propagating Shear Waves |
JP2013500752A (en) | 2009-07-31 | 2013-01-10 | スーパー ソニック イマジン | Method and apparatus for measuring physical parameters in mammalian soft tissue by propagating shear waves |
WO2011027644A1 (en) * | 2009-09-04 | 2011-03-10 | 株式会社 日立メディコ | Ultrasonic diagnostic device |
US20120136250A1 (en) * | 2009-09-04 | 2012-05-31 | Hitachi Medical Corporation | Ultrasonic diagnostic device |
US20110066030A1 (en) * | 2009-09-11 | 2011-03-17 | Mr Holdings (Hk) Ltd. | Systems and methods for shear wave field formation |
WO2011064688A1 (en) | 2009-11-25 | 2011-06-03 | Koninklijke Philips Electronics N.V. | Ultrasonic shear wave imaging with focused scanline beamforming |
US20130131511A1 (en) | 2009-11-25 | 2013-05-23 | Koninklijke Philips Electronics N.V. | Ultrasonic shear wave imaging with focused scanline beamforming |
US20110184287A1 (en) * | 2009-12-10 | 2011-07-28 | Mcaleavey Stephen A | Methods And Systems For Spatially Modulated Ultrasound Radiation Force Imaging |
US20110216311A1 (en) * | 2010-03-02 | 2011-09-08 | Li-Cor, Inc. | Method and apparatus for locking a laser with a resonant cavity |
JP2012100997A (en) | 2010-11-12 | 2012-05-31 | Hitachi Medical Corp | Ultrasonograph and operation method therefor |
WO2012077579A1 (en) | 2010-12-08 | 2012-06-14 | 株式会社日立メディコ | Ultrasound diagnosis apparatus |
US20130289402A1 (en) | 2010-12-08 | 2013-10-31 | Hitachi Medical Corporation | Ultrasound diagnosis apparatus |
US9439620B2 (en) | 2010-12-08 | 2016-09-13 | Hitachi, Ltd. | Ultrasound diagnosis apparatus |
US20120215101A1 (en) * | 2011-02-17 | 2012-08-23 | Siemens Medical Solutions Usa, Inc | Viscoelasticity Measurement Using Amplitude-Phase Modulated Ultrasound Wave |
JP2012170823A (en) | 2011-02-17 | 2012-09-10 | Siemens Medical Solutions Usa Inc | Viscoelasticity measurement using amplitude and phase-modulated ultrasonic wave |
US20120232387A1 (en) | 2011-03-08 | 2012-09-13 | Fujifilm Corporation | Ultrasound diagnostic apparatus and method of determining elasticity index reliability |
JP2012183261A (en) | 2011-03-08 | 2012-09-27 | Fujifilm Corp | Ultrasonic diagnostic apparatus and method for determining elasticity index reliability |
US8852102B2 (en) | 2011-03-08 | 2014-10-07 | Fujifilm Corporation | Ultrasound diagnostic apparatus and method of determining elasticity index reliability |
US20140221833A1 (en) * | 2011-04-08 | 2014-08-07 | Canon Kabushiki Kaisha | Subject information acquisition apparatus |
Non-Patent Citations (5)
Title |
---|
Extended European Search Report (EESR), dated Oct. 14, 2016. |
International Search Report from International Application No. PCT/JP14/051950 dated Mar. 11, 2014. |
Jeremy J. Dahl et al., "A Parallel Tracking Method for Acoustic Radiation Force Impulse Imaging", IEEE transaction on ultrasonics, ferroelectrics, and frequency control, vol. 54, No. 2, 2007, pp. 301-312. |
Shigao Chen, et al "Quantifying Elasticity and Viscosity from Measurement of Shear Wave Speed Dispersion", The Journal of the Acoustical Society of America, American Institute of Physics for the Acoustical Society of America, New York, NY vol. 115, No. 6, Jun. 1, 2004, pp. 2781-2785. |
Thomas Deffieux, et al "On the Effects of Reflected Waves in Transient Shear Wave Elastography", IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, IEEE, US, vol. 58, No. 10, Oct. 1, 2011, pp. 2032-2035. |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10485515B2 (en) | 2015-01-28 | 2019-11-26 | Konica Minolta, Inc. | Ultrasound diagnostic device and method for controlling ultrasound diagnostic device |
US10327738B2 (en) | 2015-03-31 | 2019-06-25 | Samsung Medison Co., Ltd. | Ultrasound imaging apparatus and method of processing ultrasound image thereof |
Also Published As
Publication number | Publication date |
---|---|
EP2965694A1 (en) | 2016-01-13 |
EP2965694A4 (en) | 2016-11-16 |
US20150173718A1 (en) | 2015-06-25 |
WO2014136502A1 (en) | 2014-09-12 |
JP5753633B2 (en) | 2015-07-22 |
JPWO2014136502A1 (en) | 2017-02-09 |
CN104302233B (en) | 2016-10-12 |
CN104302233A (en) | 2015-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9895137B2 (en) | Ultrasonic diagnosis device and transmission/reception method | |
JP5730978B2 (en) | Ultrasonic diagnostic apparatus and method | |
US9439620B2 (en) | Ultrasound diagnosis apparatus | |
US8494791B2 (en) | Methods and systems for improved correlation of shear displacement waveforms | |
CN104510499B (en) | Shear wave detection in medical ultrasound imaging | |
CN104434216B (en) | Shear wave estimation from analytical data | |
EP2926739B1 (en) | Acquisition control for elasticity ultrasound imaging | |
US9239373B2 (en) | Adaptive image optimization in induced wave ultrasound imaging | |
US10959704B2 (en) | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing method | |
KR102051293B1 (en) | Classification preprocessing in medical ultrasound shear wave imaging | |
US11154277B2 (en) | Tissue viscoelastic estimation from shear velocity in ultrasound medical imaging | |
KR101922522B1 (en) | Sound speed imaging using shear waves | |
RU2641067C1 (en) | Ultrasonic system and method for measurement using shear wave | |
US20240245391A1 (en) | Shear wave imaging based on ultrasound with increased pulse repetition interval | |
JP2006115937A (en) | Ultrasonic diagnostic apparatus | |
US11446007B2 (en) | Ultrasonic shearwave imaging with patient-adaptive shearwave generation | |
EP2853918B1 (en) | Shear wave detection in medical ultrasound imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: HITACHI ALOKA MEDICAL, LTD., JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TABARU, MARIE;YOSHIKAWA, HIDEKI;ASAMI, REI;AND OTHERS;SIGNING DATES FROM 20140930 TO 20141026;REEL/FRAME:034646/0226 |
|
AS | Assignment |
Owner name: HITACHI, LTD., JAPAN Free format text: MERGER AND CHANGE OF NAME;ASSIGNORS:ALOKA CO., LTD.;HITACHI ALOKA MEDICAL, LTD.;REEL/FRAME:040515/0750 Effective date: 20110401 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |
|
AS | Assignment |
Owner name: FUJIFILM HEALTHCARE CORPORATION, JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:HITACHI, LTD.;REEL/FRAME:058026/0559 Effective date: 20211013 |
|
AS | Assignment |
Owner name: FUJIFILM HEALTHCARE CORPORATION, JAPAN Free format text: CORRECTIVE ASSIGNMENT TO CORRECT THE THE PROPERTY AND APPLICATION NUMBERS PREVIOUSLY RECORDED AT REEL: 058026 FRAME: 0559. ASSIGNOR(S) HEREBY CONFIRMS THE ASSIGNMENT;ASSIGNOR:HITACHI LTD.;REEL/FRAME:058917/0853 Effective date: 20211013 |