WO2012002420A1 - 超音波診断装置及び超音波診断方法 - Google Patents
超音波診断装置及び超音波診断方法 Download PDFInfo
- Publication number
- WO2012002420A1 WO2012002420A1 PCT/JP2011/064890 JP2011064890W WO2012002420A1 WO 2012002420 A1 WO2012002420 A1 WO 2012002420A1 JP 2011064890 W JP2011064890 W JP 2011064890W WO 2012002420 A1 WO2012002420 A1 WO 2012002420A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- variation
- interest
- ultrasonic
- region
- attenuation
- 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/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0858—Detecting organic movements or changes, e.g. tumours, cysts, swellings involving measuring tissue layers, e.g. skin, interfaces
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
- A61B6/46—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient
- A61B6/467—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient characterised by special input means
- A61B6/469—Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with special arrangements for interfacing with the operator or the patient characterised by special input means for selecting a region of interest [ROI]
-
- 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/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
Definitions
- the present invention relates to an ultrasonic diagnostic apparatus and an ultrasonic diagnostic method, and in particular, uses ultrasonic waves to capture and display an ultrasonic image of a diagnosis region of a subject,
- the present invention relates to an ultrasonic diagnostic apparatus and an ultrasonic diagnostic method for diagnosing tissue characteristics of a subject based on a change in attenuation.
- an ultrasonic diagnostic apparatus that captures and displays an ultrasonic image of a diagnostic part of a subject using ultrasonic waves is known, and measures acoustic characteristics such as ultrasonic sound velocity, attenuation, and scattering.
- acoustic characteristics such as ultrasonic sound velocity, attenuation, and scattering
- Patent Document 1 irradiation is performed while changing the irradiation angle of the transmission vibrator, and the elapsed time during reception while changing the incident angle of the reception vibrator is measured and stored in a memory.
- the transmission / reception angle of the transducer is changed on the sound velocity memory, and the corresponding sound ray path is set for each, and the required time and destination Compare the elapsed time stored in the memory to obtain the error data, correct the sound speed distribution stored in the sound speed memory to minimize the error data, and obtain the sound speed by the corrected sound speed distribution
- a sound speed measurement method is disclosed.
- Patent Document 2 discloses a method for extracting and drawing information associated with propagation attenuation in a living body by taking a spatial difference between detection signals that have been divided into bands near the depth of interest.
- Patent Documents 3 to 5 Other conventional techniques related to the present invention include Patent Documents 3 to 5 and Non-Patent Documents 1 to 3.
- the conventional sound velocity or attenuation measurement method is macroscopic, and it is difficult for such a macro acoustic property measurement method to detect changes in the micro acoustic property structure depending on the type of lesion. It has been reported.
- Non-Patent Document 1 Hiroyuki Hachiya “Acoustic characteristics of biological tissue and characteristics of ultrasonic tomographic images” (Medical Imaging Technology vol.21 No.2 March 2003), when liver becomes cirrhosis and cirrhosis progresses Necrotic tissues are joined together, and in order to repair, the surrounding tissues fibrosis and form nodules, and the liver lobules are replaced with regenerative nodules.
- FIG. 12 shows an example of scatterer arrangement.
- FIG. 12 (A) shows a normal liver, and each hepatic lobule structure has a random size of about 1.0 to 1.5 mm, but as shown in FIG. 12 (B).
- Non-Patent Document 2 Koichi Akamatsu “tissue characterization by sound velocity measurement” (clinician vol.12 no.11 1986), there is no significant difference between the macro sound velocity values of normal liver and cirrhosis. There is no report.
- the macro sound speed and attenuation measurement methods that have been proposed in the past have a problem that the micro structural changes of the sound speed and attenuation may not be captured.
- the present invention has been made in view of such circumstances, and an ultrasonic diagnostic apparatus capable of capturing a micro structural change in sound velocity / attenuation due to a disease change, which is difficult with conventional absolute value measurement of sound velocity / attenuation, and An object is to provide an ultrasonic diagnostic method.
- the present invention provides an ultrasonic wave including a plurality of ultrasonic transducers that transmit ultrasonic waves to a subject, receive ultrasonic waves reflected by the subject, and output ultrasonic detection signals.
- an ultrasonic diagnostic apparatus including an acoustic probe, a target region setting unit that sets a target region in a subject, and a variation measuring unit that measures ultrasonic speed variation or attenuation variation of the ultrasonic wave in the target region.
- the ultrasonic diagnostic apparatus includes a focus point setting unit that sets at least one focus point in the focus region, and the variation measurement unit is configured to detect each of the ultrasonic transducers of the reflected wave from the focus point. The sound speed variation is measured based on the change in the reception time of the received signal.
- the variation measuring unit measures the attenuation variation based on a change in amplitude or frequency of a reception signal of each ultrasonic transducer of a reflected wave from the point of interest.
- the variation measuring unit measures the attenuation variation based on a change in attenuation of each point of interest in the region of interest.
- the variation measuring unit measures the sound speed variation by comparing the RF signal obtained as a result of matching and addition at a constant sound speed or an image obtained therefrom for each set sound speed.
- the ultrasonic diagnostic apparatus further includes a transmission focus instruction unit that causes the ultrasonic probe to transmit and focus an ultrasonic wave in a region of interest.
- the transmission focus instruction unit determines the effective range of each transmission focus in advance.
- the present invention relates to a region-of-interest setting step in which the ultrasonic diagnostic apparatus sets a region of interest in a subject, and sets at least one point of interest in the region of interest.
- An ultrasonic diagnostic method is provided that performs a property diagnosis step of diagnosing tissue properties based on an index.
- the present invention provides an ultrasonic diagnostic apparatus in which a focused area setting step for setting a focused area in a subject and at least one focused point in the focused area are set.
- a step a step of measuring an attenuation variation based on a change in amplitude or frequency of a received signal of each ultrasonic transducer of a reflected wave from a point of interest, an index calculating step for calculating a variation index from the attenuation variation, and a calculated variation
- An ultrasonic diagnostic method is provided that performs a property diagnosis step of diagnosing tissue properties based on an index.
- the ultrasonic diagnostic method further executes a transmission focus instruction step for transmitting and focusing ultrasonic waves to the region of interest with respect to the ultrasonic transducer.
- the transmission focus instruction step determines the effective range of each transmission focus in advance.
- the present invention it is possible to capture a micro structural change of sound speed / attenuation due to a disease change, which is difficult to measure with the conventional absolute value of sound speed / attenuation, and represents sound speed variation or attenuation variation.
- the tissue property can be diagnosed using the variation index.
- FIG. 1 is a block diagram showing a schematic configuration of an ultrasonic diagnostic apparatus according to an embodiment of the present invention.
- (A) to (C) are explanatory views schematically showing processing for measuring sonic variation and attenuation variation. It is explanatory drawing which shows typically a mode that two types of media exist in the path
- FIG. 1 is a block diagram showing a schematic configuration of an ultrasonic diagnostic apparatus according to an embodiment of the present invention.
- the ultrasonic diagnostic apparatus 10 transmits an ultrasonic beam from an ultrasonic probe 300 to a subject OBJ and is reflected by the subject OBJ (ultrasonic wave).
- This is a device that receives (echo) and creates and displays an ultrasonic image from the detection signal of the ultrasonic echo.
- CPU (Central Processing Unit) 100 controls each block of the ultrasonic diagnostic apparatus 10 in accordance with an operation input from the operation input unit 200.
- the operation input unit 200 is an input device that receives an operation input from an operator, and includes an operation console 202 and a pointing device 204.
- the console 202 is a display mode between a keyboard that accepts input of character information (for example, patient information), a mode that displays an amplitude image (B-mode image) alone, and a mode that displays a determination result of local sound velocity values.
- the pointing device 204 is a device that receives an input for designating an area on the screen of the display unit 104, and is, for example, a trackball or a mouse. Note that a touch panel can be used as the pointing device 204.
- the storage unit 102 is a storage device that stores a control program for causing the CPU 100 to control each block of the ultrasonic diagnostic apparatus 10, and is, for example, a hard disk or a semiconductor memory.
- the display unit 104 is, for example, a CRT (Cathode Ray Tube) display or a liquid crystal display, and displays an ultrasonic image (moving image and still image) and various setting screens.
- CTR Cathode Ray Tube
- ultrasonic image moving image and still image
- the ultrasonic probe 300 is a probe used in contact with the subject OBJ, and includes a plurality of ultrasonic transducers 302 constituting a one-dimensional or two-dimensional transducer array.
- the ultrasonic transducer 302 transmits an ultrasonic beam to the subject OBJ based on the drive signal applied from the transmission circuit 402, receives an ultrasonic echo reflected from the subject OBJ, and outputs a detection signal.
- the ultrasonic transducer 302 includes a vibrator configured by forming electrodes on both ends of a piezoelectric material (piezoelectric body).
- a piezoelectric body constituting such a vibrator
- Examples of the piezoelectric body constituting such a vibrator include a piezoelectric ceramic such as PZT (lead zirconate titanate) and a polymer such as PVDF (polyvinylidene difluoride).
- a piezoelectric element can be used.
- the piezoelectric body expands and contracts, and ultrasonic waves are generated in each vibrator by the expansion and contraction of the piezoelectric body.
- a pulsed electric signal is sent to the electrode of the vibrator
- a pulsed ultrasonic wave is generated
- a continuous wave electric signal is sent to the electrode of the vibrator
- a continuous wave ultrasonic wave is generated.
- the ultrasonic waves generated in the respective vibrators are combined to form an ultrasonic beam.
- the piezoelectric body of each vibrator expands and contracts to generate an electric signal.
- the electrical signal generated in each transducer is output to the receiving circuit 404 as an ultrasonic detection signal.
- the ultrasonic transducer 302 it is possible to use a plurality of types of elements having different ultrasonic conversion methods.
- a transducer constituted by the piezoelectric body may be used as an element that transmits ultrasonic waves
- an optical transducer of an optical detection type may be used as an element that receives ultrasonic waves.
- the light detection type ultrasonic transducer converts an ultrasonic signal into an optical signal for detection, and is, for example, a Fabry-Perot resonator or a fiber Bragg grating.
- the live mode is a mode for displaying, analyzing, and measuring an ultrasonic image (moving image) obtained by transmitting and receiving ultrasonic waves by bringing the ultrasonic probe 300 into contact with the subject OBJ.
- the CPU 100 When the ultrasound probe 300 is brought into contact with the subject OBJ and ultrasound diagnosis is started by an instruction input from the operation input unit 200, the CPU 100 outputs a control signal to the transmission / reception unit 400, and the ultrasound Transmission of the beam to the subject OBJ and reception of ultrasonic echoes from the subject OBJ are started.
- the CPU 100 sets the transmission direction of the ultrasonic beam and the reception direction of the ultrasonic echo for each ultrasonic transducer 302.
- the CPU 100 selects a transmission delay pattern according to the transmission direction of the ultrasonic beam, and selects a reception delay pattern according to the reception direction of the ultrasonic echo.
- the transmission delay pattern is pattern data of a delay time given to a drive signal in order to form an ultrasonic beam in a desired direction by ultrasonic waves transmitted from a plurality of ultrasonic transducers 302, and a reception delay pattern. Is pattern data of delay time received by the plurality of ultrasonic transducers 302.
- the transmission delay pattern and the reception delay pattern are stored in the storage unit 102 in advance.
- the CPU 100 selects a transmission delay pattern and a reception delay pattern from those stored in the storage unit 102, and outputs a control signal to the transmission / reception unit 400 according to the selected transmission delay pattern and reception delay pattern to transmit / receive ultrasonic waves. Take control.
- the transmission circuit 402 generates a drive signal in accordance with a control signal from the CPU 100 and applies the drive signal to the ultrasonic transducer 302. At this time, the transmission circuit 402 delays the drive signal applied to each ultrasonic transducer 302 based on the transmission delay pattern selected by the CPU 100. Here, the transmission circuit 402 adjusts (delays) the timing at which the drive signal is applied to each ultrasonic transducer 302 so that the ultrasonic waves transmitted from the plurality of ultrasonic transducers 302 form an ultrasonic beam. Perform focus. Note that the timing of applying the drive signal may be adjusted so that the ultrasonic waves transmitted from the plurality of ultrasonic transducers 302 reach the entire imaging region of the subject OBJ.
- the receiving circuit 404 receives and amplifies the ultrasonic detection signal output from each ultrasonic transducer 302. As described above, since the distances between the ultrasonic transducers 302 and the ultrasonic wave reflection sources in the subject OBJ are different, the time for the reflected wave to reach the ultrasonic transducers 302 is different.
- the reception circuit 404 includes a delay circuit, and the difference in arrival time of the reflected wave (delay time) according to the reception delay pattern set based on the sound speed selected by the CPU 100 (hereinafter referred to as assumed sound speed) or the sound speed distribution. Each detection signal is delayed by an amount corresponding to. Next, the reception circuit 404 performs reception focus processing by matching and adding detection signals given delay times.
- the addition circuit of the reception circuit 404 adds them.
- the phases of the ultrasonic detection signals from the other ultrasonic reflection sources cancel each other.
- the received signal from the ultrasonic reflection source XROI becomes the largest, and the focus is achieved.
- a sound ray signal hereinafter referred to as an RF signal
- the A / D converter 406 converts the analog RF signal output from the receiving circuit 404 into a digital RF signal (hereinafter referred to as RF data).
- the RF data includes phase information of the received wave (carrier wave).
- the RF data output from the A / D converter 406 is input to the signal processing unit 502 and the cine memory 602, respectively.
- the cine memory 602 sequentially stores the RF data input from the A / D converter 406.
- the cine memory 602 stores information related to the frame rate input from the CPU 100 (for example, parameters indicating the depth of the reflection position of the ultrasonic wave, the density of the scanning line, and the visual field width) in association with the RF data.
- the signal processing unit 502 corrects the attenuation by the distance according to the depth of the reflection position of the ultrasonic wave by STC (Sensitivity Time gain Control) on the RF data, and then performs the envelope detection processing to obtain the B mode.
- Image data image data representing the amplitude of ultrasonic echoes by the brightness (luminance) of a point
- STC Sensitivity Time gain Control
- the B-mode image data generated by the signal processing unit 502 is obtained by a scanning method different from a normal television signal scanning method. Therefore, a DSC (Digital Scan Converter) 504 converts (raster conversion) the B-mode image data into normal image data (for example, image data of a television signal scanning system (NTSC system)).
- the image processing unit 506 performs various necessary image processing (for example, gradation processing) on the image data input from the DSC 504.
- the image memory 508 stores the image data input from the image processing unit 506.
- the D / A converter 510 converts the image data read from the image memory 508 into an analog image signal and outputs the analog image signal to the display unit 104. Thereby, an ultrasonic image (moving image) photographed by the ultrasonic probe 300 is displayed on the display unit 104.
- the detection signal subjected to the reception focus process in the reception circuit 404 is an RF signal, but the detection signal not subjected to the reception focus process may be an RF signal.
- a plurality of ultrasonic detection signals output from the plurality of ultrasonic transducers 302 are amplified by the reception circuit 404, and the amplified detection signals, that is, RF signals are A / D converted by the A / D converter 406.
- the RF data is generated by the conversion.
- the RF data is supplied to the signal processing unit 502 and stored in the cine memory 602.
- the reception focus process is performed digitally in the signal processing unit 502.
- the cine memory playback mode is a mode for displaying, analyzing, and measuring an ultrasonic diagnostic image based on RF data stored in the cine memory 602.
- the CPU 100 switches the operation mode of the ultrasonic diagnostic apparatus 10 to the cine memory playback mode.
- the CPU 100 instructs the cine memory reproduction unit 604 to reproduce the RF data designated by the operation input from the operator.
- the cine memory reproduction unit 604 reads RF data from the cine memory 602 based on a command from the CPU 100 and transmits the RF data to the signal processing unit 502 of the image signal generation unit 500.
- the RF data transmitted from the cine memory 602 is subjected to predetermined processing (processing similar to that in the live mode) in the signal processing unit 502, DSC 504, and image processing unit 506, and converted into image data.
- the data is output to the display unit 104 via the D / A converter 510. Accordingly, an ultrasonic image (moving image or still image) based on the RF data stored in the cine memory 602 is displayed on the display unit 104.
- the freeze button on the console 202 When the freeze button on the console 202 is pressed while an ultrasonic image (moving image) is displayed in the live mode or the cine memory playback mode, the ultrasonic image displayed when the freeze button is pressed is displayed on the display unit 104. A still image is displayed. Thereby, the operator can display and observe a still image of the region of interest (ROI: Region : of Interest).
- ROI Region : of Interest
- the analysis / measurement designated by the operation input from the operator is performed.
- the data analysis measurement unit 106 acquires RF data before image processing is performed from the A / D converter 406 or the cine memory 602, and uses the RF data.
- Operator-specified analysis / measurement for example, tissue strain analysis (hardness diagnosis), blood flow measurement, tissue motion measurement, or IMT (Intima-Media Thickness) value measurement )I do.
- the data analysis measurement unit 106 performs processing for measuring a local sound speed value and calculating an index representing variation in sound speed or attenuation, which will be described in detail later.
- the analysis / measurement result by the data analysis measurement unit 106 is output to the DSC 504 of the image signal generation unit 500.
- the DSC 504 inserts the analysis / measurement result into the image data of the ultrasonic image by the data analysis measurement unit 106 and outputs the result to the display unit 104. Thereby, the ultrasonic image and the analysis / measurement result are displayed on the display unit 104.
- a mode for displaying the B mode image alone When the display mode switching button is pressed, a mode for displaying the B mode image alone, a mode for displaying the determination result of the sound speed / attenuation variation on the B mode image (for example, depending on the sound speed / attenuation variation).
- the display mode is switched between a mode in which color-coded or luminance is changed, or a display in which dots having the same sound speed / attenuation variation are connected by lines), and a mode in which the B-mode image and the sound speed / attenuation variation determination result are displayed side by side.
- the operator can find a lesion, for example, by observing the determination result of the sound speed / attenuation variation.
- a region of interest is set, the sound speed variation or attenuation variation of the region of interest is measured, and the tissue property is diagnosed.
- the present invention forms a pseudo point reflection by applying transmission focus, obtains the time difference from the reception time approximated at a constant sound speed from the reception data of each element, and determines the sound speed variation from the variation.
- attenuation (scattering, absorption) variation is measured from variation in amplitude and frequency variation approximated with constant attenuation. The measured variation is intended to be useful for diagnosis of tissue properties.
- FIG. 2 schematically shows a process for measuring sonic variation and attenuation variation.
- FIG. 2 (A) shows that the variation in sound velocity is measured by obtaining the variation from the reception time approximated at a constant sound velocity.
- FIG. 2B shows that attenuation variation is measured from variation in amplitude approximated by constant attenuation.
- FIG. 2C shows the measurement of the variation in attenuation from the variation in the center frequency approximated with constant attenuation.
- a pseudo point reflection is formed by applying the transmission focus, and the sound speed variation and attenuation variation are measured from the received data of each element.
- FIG. 2A shows a wavefront (reception time) actually measured in each element by a solid line, and a wavefront approximated assuming that the medium of the subject OBJ is constant in sound speed, by a broken line. It shows the variation in the reception time due to the sonic variation.
- FIG. 2B shows the amplitude after logarithmic compression actually measured in each element by a solid line, and shows the amplitude after logarithmic approximation approximated on the assumption that the medium of the subject OBJ is constant in attenuation.
- the variation in amplitude after logarithmic compression due to variation in attenuation (absorption and scattering) at each azimuth position is shown.
- FIG. 2C shows the center frequency actually measured in each element by a solid line, the center frequency approximated assuming that the medium of the subject OBJ is constant attenuation, and the attenuation at each azimuth position. (Absorption and scattering) The dispersion of the center frequency due to the dispersion is shown.
- the reception time, amplitude, and center frequency of the pseudo point reflection reception wave formed by setting the transmission focus have variations from the reception time, amplitude, and center frequency assuming constant sound velocity and constant attenuation. Yes.
- the reflected wave reception time from the grid point X, the amplitude, the reception time assuming the constant sound velocity / attenuation of the center frequency, the amplitude, the variation from the center frequency and the spatial frequency, the different sound speeds and It is possible to acquire information on the variation in the mixing ratio of the medium having attenuation and the spatial frequency.
- FIG. 3 shows the medium 1 and the medium 2 gathered to one side so that the mixing ratio can be clearly seen.
- L is the total length of the path from the sound source to the element (path length)
- ⁇ 1 and ⁇ 2 are the average mixing ratio of the medium 1 and the medium 2 that do not depend on the propagation path, and ⁇ depends on the path. This represents the change in the mixing ratio.
- the mixing ratio of the medium 1 and the medium 2 in the path shown in FIG. 3 is shifted from the average mixing ratio ⁇ 1: ⁇ 2 by ⁇ and becomes ( ⁇ 1 + ⁇ ) :( ⁇ 2 ⁇ ),
- the medium 1 exists as a length L * ( ⁇ 1 + ⁇ )
- the medium 2 exists as a length L * ( ⁇ 2 ⁇ ).
- the reception time t at which the element receives the ultrasonic wave emitted from the sound source in FIG. 3 is given by the following equation.
- the variation index according to the above formula (1) is an index that well represents the degree of variation. Become.
- the degree of variation is an index that well represents the degree of disease change.
- the reception time, the reception time that does not include the path variation, and the path length are unknown.
- the reception time (element reception time) can be obtained by using a known phase aberration analysis method (for example, see Patent Document 3: Japanese Patent Laid-Open No. 6-105841).
- a certain signal is used as a reference signal to detect these phase differences, the phase difference detection results of adjacent elements are compared, and the difference is defined as D.
- a sound source is assumed at a depth L from the element surface in a target composed of a variety of different sound velocity media (mixed media).
- the sound velocity (average sound velocity) and the depth obtained by assuming that the medium to the sound source is uniform are obtained from the ultrasonic element reception signals emitted from the sound source shown in FIG.
- the reception time T (X) at the element located at a distance X from directly below the sound source can be given by the following equation.
- T (X) ⁇ (L2 + X2) * ((1 / v1) * ⁇ 1 + (1 / v2) * ⁇ 2 + (1 / v3) * ⁇ 3 +7)
- ⁇ (A) represents the square root of A
- ⁇ n and vn represent the mixing ratio of the medium n and the sound speed.
- the average sound speed and each path length can be obtained by looking at each element reception time in total. Even if the change ⁇ of the mixing ratio by each path is included, it is considered that the influence is small if each element signal is viewed in total.
- Patent Document 4 Japanese Patent Application Laid-Open No. 2007-7045
- This is a method in which an average sound speed (and depth) is assumed and a value that maximizes the sharpness and contrast of the image of the sound source is obtained.
- the average reception time may be obtained by least squares fitting, and the average sound speed (and depth) corresponding thereto may be obtained.
- a variation index is obtained when another medium different from the mixed medium exists in front of each element. Specifically, the change in each path length and reception time is obtained by removing the influence of another medium.
- a local region is set so that the vicinity of the boundary between the mixed medium and another medium is the lower surface, and a plurality of grid points are set on the lower surface of the local region. Set. If the average sound velocity in this region is obtained, the sound source depth L 'in the local region is obtained, and from this, the respective path lengths to the respective lattices separated by the distance X' are obtained.
- the ambient sound velocity at each lattice on the lower surface of the local point is obtained.
- the environmental sound speed value is a sound speed value at which the contrast and sharpness of the image are the highest, and can be obtained by the above-described image analysis method.
- the waveform of the virtual received wave WX when the point of interest is the reflection point is calculated.
- the initial value of the assumed sound speed in the local area is set, the assumed sound speed is changed by one step, the received wave at each lattice point on the lower surface of the local area is calculated from the environmental sound speed, and the received wave is a delay determined by the assumed sound speed.
- a virtually synthesized virtual received wave WSUM is calculated.
- an error between the virtual received wave WX and the virtual synthesized received wave WSUM is calculated.
- the error between the virtual received wave WX and the virtual combined received wave WSUM can be calculated by a method of cross-correlating each other, a method of multiplying the virtual received wave WX by a delay obtained from the virtual combined received wave WSUM, or conversely, a virtual It is calculated by a method of performing phase matching addition by multiplying the combined received wave WSUM by the delay obtained from the virtual received wave WX.
- the time at which the ultrasonic wave propagated at the environmental sound velocity value at the point of interest arrives at each element may be set as the delay with the point of interest as the reflection point.
- an equiphase line is extracted from the phase difference of the combined received wave between adjacent elements, and the equal phase line is used as a delay, or the combined elements are simply combined.
- the phase difference at the maximum (peak) position of the received wave may be used as the delay.
- the cross-correlation peak position of the combined received wave from each element may be set as a delay.
- the error at the time of phase matching addition is obtained by a method of setting the peak to peak of the waveform after the matching addition or a method of setting the maximum amplitude after the envelope detection.
- the local sound speed value in the local region is determined. That is, the assumed sound speed value at which the difference between the virtual received wave WX and the virtual synthesized received wave WSUM is minimized is determined as the local sound speed value in the local region.
- the reception time (received wave) at the point of interest in the region of interest and each lattice point on the lower surface of the local region is obtained in advance, and the combined received waves are superimposed by superimposing the received waves of each lattice with a delay determined by the assumed sound velocity in the region of interest.
- the reception time at the point of interest and each lattice point on the lower surface of the local region can be obtained using the above-described image analysis method and phase aberration analysis method.
- the reception time (received wave) of each lattice point on the lower surface of the local region is obtained by the image analysis and phase aberration analysis methods, and then assumed in the region of interest.
- the minimum sum is set as the delay of the element, an image is generated based on the delay, and the image is analyzed, for example, the local sound speed is determined from the condition that the sharpness is maximum. good.
- the reception time (received wave) of each lattice point on the lower surface of the local region is obtained by a method of image analysis and phase aberration analysis.
- a signal that is matched and added with each delay is set as the received signal of each virtual element, and an image is generated from the received signal of each virtual element based on the assumed sound speed of the region of interest, and the image is analyzed.
- the local sound speed may be determined from the condition that the sharpness is maximized.
- the method of removing the lattice points and the lower surface of the local region is not particularly limited to a plane, and can be set as an arbitrary curved surface on the element side (near side) from the point of interest. For example, it may be set on a boundary surface of a tissue or a lesion.
- the change from the reception time determined by the average sound speed is obtained from the reception time of the signal received by each element, and a low frequency cut process is performed on it.
- a change due to another medium is removed, and each element position is converted into each lattice position on the lower surface of the local region.
- the conversion from each element position to each lattice position on the lower surface of the local area is based on the propagation of ultrasonic waves from the point of interest to each element from the local point of interest in the local area and the ambient sound speed at each lattice point or each element reception time. Since the path is obtained, it is possible to propagate in the reverse direction from each element position along this propagation path.
- the attenuation variation can be obtained as follows by the same method as the sound speed.
- Attenuation There are three types of attenuation: diffusion attenuation due to the spread of sound waves, absorption attenuation due to sound waves being absorbed into the medium and converted to heat, and scattering attenuation due to scattering of living tissue.
- the attenuation is given by exp ( ⁇ x).
- ⁇ represents an attenuation coefficient
- x represents a propagation distance.
- Non-Patent Document 3 “Ultrasonic Handbook” (Maruzen, 1999), it can be assumed that in biological tissue, the attenuation coefficient ⁇ is approximately proportional to the frequency in the frequency range of the MHz band. Due to the attenuation proportional to the frequency, the center frequency of the Gaussian pulse is shifted in proportion to the propagation distance. Using this, attenuation can be obtained from the center frequency shift.
- the attenuation variation index is considered in the same manner as in the case of sound velocity.
- a (x) A (0) ⁇ L * ( ⁇ 1 * ⁇ 1 + ⁇ 2 * ⁇ 2) ⁇ L * ⁇ * ( ⁇ 1 ⁇ 2)
- F (x) F (0) ⁇ L * ( ⁇ 1 * ⁇ 1 + ⁇ 2 * ⁇ 2) ⁇ L * ⁇ * ( ⁇ 1 ⁇ 2) (4)
- ⁇ is an attenuation coefficient including a term depending on frequency
- ⁇ is a constant determined by ⁇ and a band of a pulse wave (assuming a Gaussian pulse).
- the path length L is obtained from the average sound speed by the above formula (3), for example. Further, in order to obtain the change, an amplitude or center frequency that does not depend on the path is necessary, and it is obtained by fitting to the amplitude or center frequency of each element reception signal. At this time, the fitting may be performed assuming an average ⁇ or ⁇ based on the depth of the sound source obtained by Expression (3).
- the method for obtaining the variation index when there is a non-uniform layer is the above equation (6) from each path length in the local region of FIG. 5 and the average attenuation of the amplitude or center frequency of each grating on the lower surface of the local region. What is necessary is just to obtain
- FIG. 6 and 7 collectively show how to find the sonic variation and the attenuation variation in each step.
- FIG. 6 shows a case where there is no other medium
- FIG. 7 shows a case where there is another medium.
- each path length is derived in step S100 of FIG. If the depth of the point of interest (sound source) can be derived, each path length connecting the point of interest and each element can be derived.
- the average sound speed and depth to the point of interest can be obtained using a known image analysis technique. Specifically, it can be obtained as a value that maximizes the sharpness and contrast characteristics of the image of the point of interest.
- the amplitude of each element reception signal is obtained, and the value obtained by fitting the amplitude obtained by assuming the average attenuation and the depth up to the point of interest in the above equation (5) to obtain a value that minimizes the error.
- the method to adopt can be considered.
- the amplitude is affected not only by attenuation but also by diffusion, transmission focus, nonlinear characteristics, etc., it is difficult to apply the above formula (5).
- the difference between the two frequencies after logarithmic compression is proportional to the attenuation coefficient by using the method using the amplitude ratio of the two frequencies disclosed in Japanese Patent Publication No. 3-24868. Can be obtained from
- the center frequency of the received signal of each element is obtained, and an error is obtained by fitting the center frequency obtained by assuming the average attenuation and the depth up to the point of interest in Equation (5). It is conceivable to adopt a value that minimizes. At this time, if the center frequency at the time of transmission is known, fitting can be performed with higher accuracy.
- step S110 in FIG. 6 the reception time, amplitude, and center frequency of each element are derived.
- the reception time As a method of obtaining the reception time, as described above, it can be obtained using a known phase aberration analysis method.
- the phase aberration analysis may be performed starting from a “reception time not including path variation” obtained in advance together with the average sound speed by image analysis of the point of interest.
- a method of obtaining the amplitude there is a method of obtaining the value at the reception time described above after performing envelope detection on each element reception signal and converting it to amplitude information.
- a peak value in a predetermined range may be acquired based on the “reception time not including path variations”.
- a predetermined range based on the reception time described above is obtained from each element reception signal, and the following formula is obtained after frequency conversion: ⁇ f * P (f) df / ⁇ P (f) df There is a method to find the center of gravity.
- f is a frequency
- P (f) is a spectral density at f.
- the center frequency may be a frequency at which the spectral density peaks, or the center of the half width. Or you may make it obtain
- step S120 of FIG. 6 the reception time, amplitude, and center frequency not including the path variation of each path are derived.
- the reception time can be obtained if the average sound speed is obtained in step S100 described above.
- a curve that minimizes the error may be fitted to each element reception signal obtained in step S110.
- the reception time since the path length connecting the point of interest and each element is obtained in step S100, the reception time may be calculated assuming the average sound speed, and the average sound speed (reception time) that minimizes the error may be obtained.
- the amplitude is obtained by fitting a curve that minimizes the error to the amplitude of each element reception signal obtained in step S110.
- the amplitude is calculated assuming the average attenuation of Equation (5), and the average attenuation (amplitude) that minimizes the error is obtained. May be.
- the center frequency of each element is obtained at the same time.
- a curve that minimizes the error may be fitted to the center frequency of each element reception signal obtained in step S110.
- the center frequency may be calculated assuming average attenuation, and the average attenuation (center frequency) that minimizes the error may be obtained.
- fitting can be performed with high accuracy.
- step S130 of FIG. 6 a change in reception time, amplitude, and center frequency of each path is derived. This can be obtained by subtracting the value obtained in step S120 from the reception time, amplitude, and center frequency obtained in step S110.
- step S140 of FIG. 6 an index is derived by dividing the change by the path length, and then a variation index is derived.
- the standard deviation or the maximum value of the index of each route may be used as the variation index.
- the flowchart in FIG. 7 is substantially the same as the flowchart in FIG. 6 described above.
- “each path” becomes “each path in the local region” and “each element reception time, amplitude, The difference is that the “center frequency” is a “value on the lower surface of the local region”.
- each path length in the local region is derived (pseudo point reflection depth is derived).
- a local region is set such that each point of interest (sound source) is on the upper surface and a boundary near another medium is on the lower surface, and each path length in the local region is obtained.
- the depth of the point of interest in the local region is obtained.
- a method of obtaining together with the average sound speed in the local region is preferably used. As described above, there are various methods for obtaining the average sound speed (local sound speed) in the local region.
- step S210 of FIG. 7 values on the lower surface of the local region of each element reception time, amplitude, and center frequency are derived.
- the reception time (or average sound speed) of each grating on the lower surface of the local region is obtained by using a known image analysis and phase aberration analysis method, and this is used as a delay.
- the reception time (received wave) at the point of interest is obtained by image analysis and phase aberration analysis.
- a signal obtained by matching and adding the reception wave at the point of interest with each delay is set as a reception signal of each virtual element.
- the local reception time at the point of interest is obtained by performing the phase aberration analysis of the reception signal of each virtual element.
- each lattice is regarded as a virtual element, and the latest time among the times obtained by subtracting the delay from each element reception time at the point of interest is adopted as the local reception time of each virtual element.
- the received wave of each grating on the lower surface of the local area is considered to be the same, a representative received wave is determined, and the received wave at the point of interest is inverted by the received wave representing the lower grating of the local area.
- the local reception time at the point of interest is obtained by performing deconvolution.
- the deconvolution process can be performed on each element reception signal or on its frequency space.
- the local reception time may be searched so that the error from the reception wave) is minimized.
- search algorithms For example, a quasi-Newton method may be used.
- the center frequency As a method of obtaining the center frequency, the local reception time or average sound speed of the local area and the reception time or average sound speed of each lattice on the lower surface of the local area are obtained in advance. Then, the propagation path of the point of interest ⁇ each lattice ⁇ each element is obtained. The center frequency at the time of transmission is assumed to be known.
- the center frequency at the time of transmission is unknown, if a uniform attenuation is assumed in all paths of the lattice ⁇ each element, the attenuation coefficient can be obtained and the shift amount can be obtained (however, the center frequency at the time of transmission has better accuracy. ).
- the attenuation of each lattice ⁇ each element is obtained in advance from the center frequency shift amount, and the propagation path of each lattice ⁇ each element is obtained.
- the attenuation of each grating ⁇ each element is corrected from the amplitude of each element, and the amplitude of each grating is obtained.
- step S220 in FIG. 7 the reception time, amplitude, and center frequency not including the path variation in the local region are derived.
- the reception time it is obtained from the average sound speed and the path length of the local area obtained in step S200. Alternatively, it may be obtained by fitting a curve that minimizes the error at the reception time of each grid obtained in step S210. At this time, since the path length connecting the point of interest and each grid is obtained in step S200, the reception time may be calculated assuming the average sound speed, and the average sound speed (reception time) that minimizes the error may be obtained.
- the amplitude As a method for obtaining the amplitude, it is possible to obtain the amplitude by fitting a curve with a minimum error to the amplitude of each lattice obtained in step S210. At this time, since the path length connecting the point of interest and each grid is obtained in step S200, the amplitude is calculated assuming the average attenuation of equation (5), and the average sound speed (amplitude) that minimizes the error is obtained. Also good.
- a curve having a minimum error is fitted to the center frequency of each lattice obtained in step S210.
- each path length is obtained in step S200, it is possible to calculate the center frequency assuming the average attenuation and obtain the average attenuation (center frequency) that minimizes the error.
- steps S230 and S240 is the same as the processing in steps S130 and S140 in FIG.
- the lower surface of the local region set when there is another medium does not necessarily have to be near the boundary as shown in FIG. 5, or may be a curved surface instead of a flat surface.
- the transmission focus is set not only on the local area but also on the lower surface of the local area.
- the amount to be standardized includes the depth and the amount close to the depth in addition to the path length, and the reception time / frequency shift amount of the middle element (or lattice). Normalization by these is not necessary when evaluating the variation with the depth of the point of interest (region of interest) (the depth excluding the other medium if there is another medium) being constant.
- the sound speed / attenuation is an amount dependent only on ⁇ separated from the depth, so the sound speed / attenuation may be used as an index (in this case, the ratio of the depth to the aperture is Seems to be unrelated).
- the sound speed variation is obtained based on the sharpness of an image (B-mode image) obtained by performing envelope detection processing on an RF signal or RF signal (RF data) as a result of matching and adding at a constant sound speed.
- B-mode image an image obtained by performing envelope detection processing on an RF signal or RF signal (RF data) as a result of matching and adding at a constant sound speed.
- FIG. 8 shows the case where the sound speed of the target tissue is uniform
- FIG. 9 shows the case where the sound speed of the target tissue is non-uniform.
- the sound speed at both ends for example, the difference between the sharpness at 1400 m / sec and 1650 m / sec, is used as the variation index with respect to the maximum value of the sharpness. Also, it is standardized by the maximum value in consideration of luminance dependency.
- the maximum value of sharpness may be simply set, or the half value width may be set.
- variation indicators there are the following as variation indicators. That is, there is a sound speed / attenuation range in which the absolute value of the difference from each sound speed / attenuation approximate curve and the mean value of the squares are the minimum predetermined ratio. In addition, there are sound velocity / attenuation ranges on both sides in contact with the measured reception time, amplitude, and center frequency. In addition, in the case of the sound speed, there are a sound speed range in which the focus index of the matched and added image has a maximum predetermined ratio, a standard deviation of sound speed / attenuation obtained by dividing the opening into small openings, and the like.
- the average sound velocity and average attenuation itself may vary depending on the position of the point of interest. Therefore, the standard deviation of the variation in the region of interest of the average sound speed / average attenuation at each point of interest may be used as the variation index.
- the variation index described above is an index based on the variation in reception time, amplitude, and center frequency, but it is considered that the spatial frequency varies. Specifically, since it is considered that the frequency with respect to the azimuth position corresponding to the change of the reception time / amplitude and the center frequency shown in FIG. 2 changes, it may be a variation index based on this.
- the change in reception time / amplitude and center frequency can be obtained by the flowchart of FIG. 6 when there is no other medium, and by the flowchart of FIG. 7 when there is another medium.
- the magnitude of the change in the reception time, amplitude, and center frequency increases with the depth of the point of interest. There is no need. However, since there is a slight difference in how to increase the amount of change for each route, the amount of change may be corrected by normalizing with the length of each route. That is, either the change obtained by the flowchart of FIG. 6 or FIG. 7 or the change normalized by the path length may be used as an index.
- the center frequency or band of the frequency distribution with respect to the azimuth position of the index obtained from the above or a variable based on them is obtained as a variation index.
- the uniform and small liver lobule is replaced with a non-uniform and large nodule, so that the center frequency moves to the low frequency side or the band is widened. Can be diagnosed.
- the center frequency is obtained by ⁇ f * P (f) df / ⁇ P (f) df.
- f represents the frequency
- P (f) represents the amplitude of the frequency f.
- the center frequency may be a frequency having a maximum amplitude, a center frequency in a band having a maximum amplitude, or a frequency at which the integral value of P (f) becomes a half value.
- f0 represents the center frequency. This may remain dispersed.
- the band may be a band having a predetermined ratio of the amplitude of the center frequency and the maximum amplitude, or the integral value of P (f) around the center frequency and the frequency of the maximum amplitude is a predetermined ratio of the total integral value. It is good also as a zone
- the frequency distribution distortion may be obtained as a variation index. This is obtained from the third moment ⁇ (f ⁇ f0) 3 * P (f) df / ⁇ P (f) df of the frequency distribution.
- the method for obtaining the variation index using the received time / amplitude / center frequency change and the path length standardized change as an index has been described.
- the received time / amplitude / center frequency may be directly used as an index instead of the change.
- the component of “reception time / amplitude / center frequency not including path variation” is included in the extremely low frequency of the frequency distribution, the extremely low frequency component may be excluded when calculating the variation index.
- the variation index may be obtained based on the spatial frequency of variation due to the position of the target point of average sound speed and average attenuation.
- a two-dimensional frequency distribution of average sound speed and average attenuation in the region of interest is obtained, and obtained from the center frequency, band, and distortion.
- FIG. 10 is a flowchart showing the flow of the entire process for obtaining a variation index of sound speed or attenuation.
- step S300 in FIG. 10 a region of interest is set.
- step S310 a predetermined number of transmission focal points are set along each line in the region of interest, the corresponding transmission focus is performed, and each element is received.
- the effective range of each transmission focus may be determined in advance as follows.
- a transmission focus No (n) is set, and a specified line width is added / subtracted to a predetermined specified line No to obtain a line No (m), and each of the focus No (n) and the line No (m) Read the element reception signal.
- the set sound speed No (k) is set, the reception focus of the set sound speed No (k) is performed on the reception signal of the line No (m) of the transmission focus No (n), and the index or the image is stored. .
- This process is repeated while changing the set sound speed No.
- the line No is changed, and the set sound speed No (k) is changed and transmitted to the new line No as described above. Implement focus.
- the ambient sound speed (average sound speed) at each depth is calculated from the indicators or images of all lines for each set sound speed, the standard deviation in the depth direction of the environmental sound speed is calculated, and the minimum point is determined as the actual depth of focus.
- the effective range of the transmission focus No (n) is obtained. Then, the transmission focus No is changed, and the effective area for the next transmission focus No is obtained in the same manner as above.
- transmission focus may be performed corresponding to each point of interest.
- step S320 a predetermined number (i0) of points of interest are set in the region of interest. There may be one point of interest, that is, the predetermined number i0 may be one.
- step S340 each element reception data of the transmission focus corresponding to the i-th point of interest is selected, and a variation index representing the sound speed variation or attenuation variation of the i-th point of interest is obtained from the corresponding data.
- the method of obtaining the variation index is obtained by the method shown in the flowchart of FIG. 6 described above when there is no other medium.
- it is obtained by the method shown in the flowchart of FIG.
- a local region is set separately from the region of interest, and each element reception data of the transmission focus corresponding to each lattice point on the lower surface of the local region is also used.
- step S350 of FIG. 10 the point of interest number i is incremented by 1 (1 is added to i), and in the next step S360, has i exceeded the set number of points of interest (predetermined number i0)? Judge whether or not.
- step S340 the process returns to step S340, and the process for obtaining the variation index representing the sonic variation or attenuation variation of the point of interest i is repeated.
- step S370 the sum of the variation indices at each point of interest i is taken and used as the variation index in the region of interest.
- the variation index of the region of interest is the sum of the variation indexes of the respective points of interest.
- An index may be used.
- the dispersion index at each point of interest in the shallow region within the region of interest is compared with the variation index in the deep region.
- the difference index of the variation index in the depth direction at each point of interest, or the difference value from the variation index at the point of interest in the shallow area, or the normalized value is summed for each line, and the variation index of the region of interest It is good.
- the sum of the variation indexes in the deep region may be taken as the variation index of the region of interest.
- the variation in the time, amplitude, and frequency of the received signal occurs in the propagation process from the pseudo point reflection, and the variation is sought.
- variations in time, amplitude, and frequency of these received signals occur not only in the propagation process from the pseudo point reflection, but also due to ambient interference when forming the pseudo point reflection.
- the transmission focus is not concentrated at one point due to sound speed, attenuation, or scattering unevenness, or if the surrounding scattering is larger than the point of interest due to scattering unevenness, interference from ambient scattering is received, and as a result, reception Variations in the time, amplitude, and frequency of signals occur.
- FIG. 11 shows an example in which the received signals vary.
- the horizontal axis represents the element position
- the vertical axis represents the reception time
- the gray scale represents the amplitude of the received signal.
- FIG. 11A shows a received signal having a clean wavefront, and the wavefront having a parabolic shape indicates that the signal arrived first in the middle element. Further, FIG. 11B shows a reception signal that varies.
- the present invention is intended for measurement including the influence of ambient scattering caused by sound speed, attenuation, and scattering unevenness in pseudo point reflection formation by transmission focus.
- the histogram shape feature amount may be used as a variation index such as skewness or kurtosis, or an indicator for each point of interest.
- the standard deviation of the distribution in the region of interest, the histogram shape feature amount, or the texture feature amount by the co-occurrence matrix for example, uniformity, contrast, correlation, entropy, etc. may be used as the variation index.
- the histogram feature amount, the texture feature amount, and the like can be similarly used as the variation index.
- a variation index may be obtained from a plurality of feature quantities using, for example, a multiple regression equation.
- a large amount of data on specific lesions and variations in sound speed or attenuation is accumulated, and based on this data, the correspondence between the variation index value and the state of tissue properties is statistically determined to diagnose tissue properties
- a threshold value for the variation index is set in advance.
- a variation index is obtained by the method described above, and the tissue characteristic is diagnosed by comparing it with a preset threshold value. This facilitates diagnosis of tissue properties.
Abstract
Description
=L*(1/v1)*ρ1+L*(1/v2)*ρ2
+L*Δρ*((1/v1)-(1/v2))
このうち、経路に依存しない(経路バラツキを含まない)受信時刻
L*(1/v1)*ρ1+L*(1/v2)*ρ2
を引けば、経路に依存した受信時刻の変化分が、
L*Δρ*((1/v1)-(1/v2))
で与えられることがわかる。
しかし、経路によって経路長Lにおける混合割合の変化分Δρは異なるため、上記式(1)で表される指標は、経路によって異なる。
1/ 平均音速 =((1/v1)*ρ1+(1/v2)*ρ2)
に分離することができる。
ここで、記号√(A)は、Aの平方根を表し、ρn、vnは媒質nの混合割合及び音速を表す。またここでは混合割合の変化分Δρは含まないとする。
深さ=L ・・・・・・(3)
上記式(3)の平均音速は、前述した平均音速であり、また各経路長は、深さLと素子位置Xから求めることができる。
F(x)=F(0)-L*(β1*ρ1+β2*ρ2)-L*Δρ*(β1-β2)・・・・・・(4)
ただし、振幅に関しては、拡散、送信フォーカス、非線形特性、指向性などの影響を無視している。ここで、αは、周波数に依存する項を含む減衰係数、βは、αとパルス波(ガウシアンパルスを仮定)の帯域によって決まる定数である。
F(x)=F(0)-L*(β1*ρ1+β2*ρ2)・・・・・・(5)
上記式(4)から式(5)を引き、求めた振幅の変化分L*Δρ*(α1-α2)または中心周波数の変化分L*Δρ*(β1-β2)を経路長Lで割れば、次の式(6)が得られる。
これにより経路長に依存しない指標が得られる。
-(送信経路長)*(α1*ρ1+α2*ρ2)
を加えればよいだけである。
∫f*P(f)df/∫P(f)df
によって重心を求める方法がある。
この値を、(中心周波数[各素子の受信信号]-中心周波数[送信時])から引いた値が、格子→各素子の伝播パスの減衰による中心周波数シフト量を表している。
=∫f2*P(f)df/∫P(f)df-f02の平方根をとることによって求められる。ここで、f0は中心周波数を表す。これは分散のままでもよい。また、帯域として他に中心周波数の振幅や最大振幅の所定比となる帯域としてもよいし、中心周波数や振幅最大の周波数を中心としたP(f)の積分値が全積分値の所定比となる帯域としてもよい。
Claims (11)
- 超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して超音波検出信号を出力する複数の超音波トランスデューサを含む超音波探触子と、
前記被検体内の着目領域を設定する着目領域設定部と、
前記着目領域における超音波の音速バラツキまたは減衰バラツキを計測するバラツキ計測部と、
を備える超音波診断装置。 - 請求項1に記載の超音波診断装置であって、さらに、前記着目領域内に少なくとも1つ以上の着目点を設定する着目点設定部を有し、前記バラツキ計測部は、前記着目点からの反射波の前記各超音波トランスデューサの受信信号の受信時刻の変化に基づいて前記音速バラツキを計測する超音波診断装置。
- 前記バラツキ計測部は、前記着目点からの反射波の前記各超音波トランスデューサの受信信号の振幅または周波数の変化に基づいて前記減衰バラツキを計測する請求項2に記載の超音波診断装置。
- 前記バラツキ計測部は、前記着目領域内の前記各着目点の減衰の変化に基づいて前記減衰バラツキを計測する請求項2に記載の超音波診断装置。
- 前記バラツキ計測部は、一定音速で整合加算した結果のRF信号又はこれから得られる画像を設定音速毎に比較して、音速バラツキを計測する請求項1に記載の超音波診断装置。
- 請求項1から5のいずれか1項に記載の超音波診断装置であって、さらに、前記超音波探触子に対して前記着目領域に前記超音波を送信フォーカスさせる送信フォーカス指示部を備える超音波診断装置。
- 前記送信フォーカス指示部は、事前に各送信フォーカスの有効域を判定する請求項6に記載の超音波診断装置。
- 超音波診断装置が、
被検体内の着目領域を設定する着目領域設定ステップと、
前記着目領域内に少なくとも1つ以上の着目点を設定するステップと、
前記着目点からの反射波の各超音波トランスデューサの受信信号の受信時刻の変化に基づいて音速バラツキを計測するステップと、
前記音速バラツキからバラツキ指標を算出する指標算出ステップと、
前記算出したバラツキ指標に基づいて組織性状を診断する性状診断ステップと、
を実行する超音波診断方法。 - 超音波診断装置が、
被検体内の着目領域を設定する着目領域設定ステップと、
前記着目領域内に少なくとも1つ以上の着目点を設定するステップと、
前記着目点からの反射波の各超音波トランスデューサの受信信号の振幅または周波数の変化に基づいて前記減衰バラツキを計測するステップと、
前記減衰バラツキからバラツキ指標を算出する指標算出ステップと、
前記算出したバラツキ指標に基づいて組織性状を診断する性状診断ステップと、
をを実行する超音波診断方法。 - 請求項8または9に記載の超音波診断方法であって、さらに、前記超音波トランスデューサに対して前記着目領域に前記超音波を送信フォーカスさせる送信フォーカス指示ステップを実行する超音波診断方法。
- 前記送信フォーカス指示ステップは、事前に各送信フォーカスの有効域を判定する請求項10に記載の超音波診断方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201180032787.4A CN102970936B (zh) | 2010-06-30 | 2011-06-29 | 超声诊断装置和超声诊断方法 |
JP2012522652A JP5798117B2 (ja) | 2010-06-30 | 2011-06-29 | 超音波診断装置及び超音波診断装置の作動方法 |
EP11800882.0A EP2589341A1 (en) | 2010-06-30 | 2011-06-29 | Ultrasound diagnostic device and ultrasound diagnostic method |
US13/730,231 US10874375B2 (en) | 2010-06-30 | 2012-12-28 | Ultrasound diagnostic device and ultrasound diagnostic method |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010149209 | 2010-06-30 | ||
JP2010-149209 | 2010-06-30 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/730,231 Continuation US10874375B2 (en) | 2010-06-30 | 2012-12-28 | Ultrasound diagnostic device and ultrasound diagnostic method |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012002420A1 true WO2012002420A1 (ja) | 2012-01-05 |
Family
ID=45402125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2011/064890 WO2012002420A1 (ja) | 2010-06-30 | 2011-06-29 | 超音波診断装置及び超音波診断方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US10874375B2 (ja) |
EP (1) | EP2589341A1 (ja) |
JP (1) | JP5798117B2 (ja) |
CN (1) | CN102970936B (ja) |
WO (1) | WO2012002420A1 (ja) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013176260A1 (ja) * | 2012-05-25 | 2013-11-28 | 富士フイルム株式会社 | 超音波診断装置、超音波診断方法、及び超音波診断プログラム |
WO2013176259A1 (ja) * | 2012-05-25 | 2013-11-28 | 富士フイルム株式会社 | 超音波診断装置、超音波診断方法、及び超音波診断プログラム |
JP2017205482A (ja) * | 2016-05-17 | 2017-11-24 | 東芝メディカルシステムズ株式会社 | 支援装置及び支援方法 |
US11266377B2 (en) | 2016-05-17 | 2022-03-08 | Canon Medical Systems Corporation | Support apparatus and support method |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014200417A1 (en) * | 2013-06-10 | 2014-12-18 | Medscienta Ab | Method and system for determining a property of a non-homogeneous material |
KR20150118732A (ko) * | 2014-04-15 | 2015-10-23 | 삼성전자주식회사 | 초음파 장치 및 그 제어 방법 |
WO2017056566A1 (ja) * | 2015-09-29 | 2017-04-06 | 富士フイルム株式会社 | 音速算出システムおよび音速算出方法 |
CN108366784A (zh) * | 2015-11-30 | 2018-08-03 | 奥林巴斯株式会社 | 超声波观测装置、超声波观测装置的工作方法以及超声波观测装置的工作程序 |
WO2017110756A1 (ja) * | 2015-12-24 | 2017-06-29 | オリンパス株式会社 | 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム |
DE102016125016B4 (de) * | 2016-10-20 | 2022-03-10 | OSRAM Opto Semiconductors Gesellschaft mit beschränkter Haftung | Verfahren zur ultraschallmikroskopischen Vermessung von Halbleiterproben, Computerprogramm zur ultraschallmikroskopischen Vermessung von Halbleiterproben, Computerprogrammprodukt und Ultraschallmikroskop |
EP3698722A4 (en) * | 2017-10-17 | 2020-12-23 | FUJIFILM Corporation | SOUND WAVE DIAGNOSTIC DEVICE AND METHOD OF CONTROLLING A SOUND WAVE DIAGNOSTIC DEVICE |
JP6930668B2 (ja) * | 2018-10-22 | 2021-09-01 | 株式会社Lily MedTech | 超音波診断システム |
CN109738517B (zh) * | 2018-12-21 | 2021-07-02 | 河海大学 | 一种轻型节电扫射式声波探测器及其探测方法 |
JP7435322B2 (ja) * | 2020-07-03 | 2024-02-21 | コニカミノルタ株式会社 | 超音波探傷装置、超音波探傷方法、および、プログラム |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60220051A (ja) * | 1984-04-13 | 1985-11-02 | 株式会社日立製作所 | 超音波診断装置 |
JPS6125536A (ja) * | 1984-07-17 | 1986-02-04 | アロカ株式会社 | 超音波診断装置 |
JPS6187538A (ja) * | 1984-10-05 | 1986-05-02 | 株式会社日立製作所 | 超音波診断装置 |
JPH0373136A (ja) * | 1989-08-11 | 1991-03-28 | Fujitsu Ltd | 超音波診断装置 |
JPH0435653A (ja) * | 1990-05-31 | 1992-02-06 | Fujitsu Ltd | 超音波診断装置 |
JPH05111484A (ja) * | 1991-10-24 | 1993-05-07 | Matsushita Electric Ind Co Ltd | 医用超音波画像評価装置 |
JPH0670929A (ja) * | 1992-08-25 | 1994-03-15 | Fuji Electric Co Ltd | 超音波診断装置 |
JP2003339698A (ja) * | 2002-05-23 | 2003-12-02 | Hitachi Medical Corp | 超音波診断装置 |
JP2004073697A (ja) * | 2002-08-22 | 2004-03-11 | Toshiba Corp | 超音波照射装置 |
WO2009060751A1 (ja) * | 2007-11-09 | 2009-05-14 | Hitachi Medical Corporation | 超音波診断装置とその作動方法及び超音波画像診断プログラム |
Family Cites Families (38)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS59208457A (ja) * | 1983-05-12 | 1984-11-26 | Hitachi Ltd | 超音波送受信方法とその装置 |
US4574635A (en) * | 1984-11-02 | 1986-03-11 | North American Philips Corporation | Monitoring of frequency shift of ultrasound pulses in tissue |
JPS61154651A (ja) * | 1984-12-28 | 1986-07-14 | 株式会社東芝 | 超音波診断装置 |
US4830015A (en) | 1986-09-16 | 1989-05-16 | Kabushiki Kaisha Toshiba | Method and system for measuring an ultrasound tissue characterization |
JPH0324868A (ja) | 1989-06-22 | 1991-02-01 | Nippon Chemicon Corp | 垂直偏向回路 |
JPH06105841A (ja) | 1991-06-04 | 1994-04-19 | Yokogawa Medical Syst Ltd | 超音波診断装置 |
JPH0595946A (ja) | 1991-10-09 | 1993-04-20 | Yokogawa Medical Syst Ltd | 超音波の音速測定法及び音速測定手段を備えた超音波診断装置 |
JPH0751270A (ja) | 1993-08-13 | 1995-02-28 | Hitachi Medical Corp | 超音波診断装置 |
US5465723A (en) * | 1994-06-02 | 1995-11-14 | Vingmed Sound A/S | Method and apparatus for ultrasound imaging |
US5531117A (en) * | 1994-09-26 | 1996-07-02 | General Electric Company | Closed loop maximum likelihood phase aberration correction in phased-array imaging systems |
US5810731A (en) * | 1995-11-13 | 1998-09-22 | Artann Laboratories | Method and apparatus for elasticity imaging using remotely induced shear wave |
JP2001170046A (ja) | 1999-12-17 | 2001-06-26 | Olympus Optical Co Ltd | 生体組織性状診断装置 |
JP2001238884A (ja) | 2000-02-29 | 2001-09-04 | Toshiba Corp | 超音波診断装置及び超音波による組織性状の定量解析方法 |
JP4244300B2 (ja) | 2003-03-24 | 2009-03-25 | 富士フイルム株式会社 | 超音波送受信装置 |
US7611465B2 (en) * | 2003-07-15 | 2009-11-03 | Board Of Regents, The University Of Texas System | Rapid and accurate detection of bone quality using ultrasound critical angle reflectometry |
JP4713112B2 (ja) | 2003-09-30 | 2011-06-29 | 富士フイルム株式会社 | 超音波診断装置 |
US20050070796A1 (en) * | 2003-09-30 | 2005-03-31 | Fuji Photo Film Co., Ltd. | Ultrasonic diagnosing apparatus |
JP4594836B2 (ja) | 2004-09-29 | 2010-12-08 | 富士フイルム株式会社 | 超音波撮像装置 |
JP4575738B2 (ja) | 2004-09-29 | 2010-11-04 | 富士フイルム株式会社 | 超音波画像境界抽出方法及び超音波画像境界抽出装置、並びに、超音波撮像装置 |
JP4585326B2 (ja) * | 2005-02-08 | 2010-11-24 | 富士フイルム株式会社 | 超音波撮像装置及び超音波撮像方法 |
JP2006255014A (ja) | 2005-03-15 | 2006-09-28 | Fuji Photo Film Co Ltd | 超音波診断装置 |
JP4713921B2 (ja) | 2005-04-15 | 2011-06-29 | 株式会社東芝 | 超音波診断装置及び超音波信号処理プログラム |
JP4817728B2 (ja) | 2005-06-29 | 2011-11-16 | 株式会社東芝 | 超音波診断装置 |
US20070083110A1 (en) * | 2005-10-09 | 2007-04-12 | Sonowise, Inc. | Programmable phase velocity in an ultrasonic imaging system |
FR2899336B1 (fr) * | 2006-03-29 | 2008-07-04 | Super Sonic Imagine | Procede et dispositif pour l'imagerie d'un milieu viscoelastique |
US8870771B2 (en) | 2007-05-04 | 2014-10-28 | Barbara Ann Karmanos Cancer Institute | Method and apparatus for categorizing breast density and assessing cancer risk utilizing acoustic parameters |
US8545407B2 (en) | 2007-05-16 | 2013-10-01 | Super Sonic Imagine | Method and device for measuring a mean value of visco-elasticity of a region of interest |
JP5371344B2 (ja) * | 2007-09-18 | 2013-12-18 | 富士フイルム株式会社 | 超音波診断方法及び装置 |
JP5448328B2 (ja) * | 2007-10-30 | 2014-03-19 | 株式会社東芝 | 超音波診断装置及び画像データ生成装置 |
JP5637653B2 (ja) * | 2008-08-18 | 2014-12-10 | 株式会社東芝 | 医用画像処理装置、超音波診断装置、及び医用画像処理プログラム |
JP2010051553A (ja) | 2008-08-28 | 2010-03-11 | Konica Minolta Medical & Graphic Inc | 超音波診断装置および方法 |
JP2010099452A (ja) | 2008-09-25 | 2010-05-06 | Fujifilm Corp | 超音波診断装置及び超音波診断方法 |
JP5470900B2 (ja) | 2009-02-23 | 2014-04-16 | コニカミノルタ株式会社 | 超音波診断装置 |
US10368843B2 (en) * | 2009-11-25 | 2019-08-06 | Koninklijke Philips N.V. | Ultrasonic shear wave imaging with focused scanline beamforming |
US9351707B2 (en) | 2010-04-05 | 2016-05-31 | Hitachi Aloka Medical, Ltd. | Methods and apparatus to determine shear wave propagation property |
JP5467922B2 (ja) * | 2010-04-30 | 2014-04-09 | 日立アロカメディカル株式会社 | 超音波診断装置 |
CN102970937B (zh) | 2010-06-30 | 2015-07-08 | 富士胶片株式会社 | 超声诊断装置和超声诊断方法 |
JP5389118B2 (ja) | 2011-08-19 | 2014-01-15 | 株式会社東芝 | 超音波診断装置 |
-
2011
- 2011-06-29 EP EP11800882.0A patent/EP2589341A1/en not_active Withdrawn
- 2011-06-29 WO PCT/JP2011/064890 patent/WO2012002420A1/ja active Application Filing
- 2011-06-29 CN CN201180032787.4A patent/CN102970936B/zh active Active
- 2011-06-29 JP JP2012522652A patent/JP5798117B2/ja active Active
-
2012
- 2012-12-28 US US13/730,231 patent/US10874375B2/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS60220051A (ja) * | 1984-04-13 | 1985-11-02 | 株式会社日立製作所 | 超音波診断装置 |
JPS6125536A (ja) * | 1984-07-17 | 1986-02-04 | アロカ株式会社 | 超音波診断装置 |
JPS6187538A (ja) * | 1984-10-05 | 1986-05-02 | 株式会社日立製作所 | 超音波診断装置 |
JPH0373136A (ja) * | 1989-08-11 | 1991-03-28 | Fujitsu Ltd | 超音波診断装置 |
JPH0435653A (ja) * | 1990-05-31 | 1992-02-06 | Fujitsu Ltd | 超音波診断装置 |
JPH05111484A (ja) * | 1991-10-24 | 1993-05-07 | Matsushita Electric Ind Co Ltd | 医用超音波画像評価装置 |
JPH0670929A (ja) * | 1992-08-25 | 1994-03-15 | Fuji Electric Co Ltd | 超音波診断装置 |
JP2003339698A (ja) * | 2002-05-23 | 2003-12-02 | Hitachi Medical Corp | 超音波診断装置 |
JP2004073697A (ja) * | 2002-08-22 | 2004-03-11 | Toshiba Corp | 超音波照射装置 |
WO2009060751A1 (ja) * | 2007-11-09 | 2009-05-14 | Hitachi Medical Corporation | 超音波診断装置とその作動方法及び超音波画像診断プログラム |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013176260A1 (ja) * | 2012-05-25 | 2013-11-28 | 富士フイルム株式会社 | 超音波診断装置、超音波診断方法、及び超音波診断プログラム |
WO2013176259A1 (ja) * | 2012-05-25 | 2013-11-28 | 富士フイルム株式会社 | 超音波診断装置、超音波診断方法、及び超音波診断プログラム |
JP2013244218A (ja) * | 2012-05-25 | 2013-12-09 | Fujifilm Corp | 超音波診断装置、超音波診断方法、及び超音波診断プログラム |
JP2013244219A (ja) * | 2012-05-25 | 2013-12-09 | Fujifilm Corp | 超音波診断装置、超音波診断方法、及び超音波診断プログラム |
US20150057543A1 (en) * | 2012-05-25 | 2015-02-26 | Fijifilm Corporation | Ultrasound diagnostic device, ultrasound diagnostic method and ultrasound diagnostic program storage medium |
US10595817B2 (en) | 2012-05-25 | 2020-03-24 | Fujifilm Corporation | Ultrasound diagnostic device, ultrasound diagnostic method and ultrasound diagnostic program storage medium |
US10631822B2 (en) | 2012-05-25 | 2020-04-28 | Fijifilm Corporation | Ultrasound diagnostic device, ultrasound diagnostic method and ultrasound diagnostic program storage medium |
JP2017205482A (ja) * | 2016-05-17 | 2017-11-24 | 東芝メディカルシステムズ株式会社 | 支援装置及び支援方法 |
JP7019301B2 (ja) | 2016-05-17 | 2022-02-15 | キヤノンメディカルシステムズ株式会社 | 支援装置 |
US11266377B2 (en) | 2016-05-17 | 2022-03-08 | Canon Medical Systems Corporation | Support apparatus and support method |
Also Published As
Publication number | Publication date |
---|---|
JP5798117B2 (ja) | 2015-10-21 |
US10874375B2 (en) | 2020-12-29 |
EP2589341A1 (en) | 2013-05-08 |
US20130123628A1 (en) | 2013-05-16 |
JPWO2012002420A1 (ja) | 2013-08-29 |
CN102970936B (zh) | 2015-07-29 |
CN102970936A (zh) | 2013-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5798117B2 (ja) | 超音波診断装置及び超音波診断装置の作動方法 | |
JP5808325B2 (ja) | 超音波診断装置及び超音波診断装置の作動方法 | |
JP6462340B2 (ja) | 超音波診断装置、画像処理装置及び画像処理方法 | |
EP1974672B9 (en) | Ultrasonic imaging apparatus and ultrasonic velocity optimization method | |
JP5230106B2 (ja) | 超音波診断装置、imt計測方法及びimt計測プログラム | |
JP5925438B2 (ja) | 超音波診断装置 | |
US20120259225A1 (en) | Ultrasound diagnostic apparatus and ultrasound image producing method | |
JP6008581B2 (ja) | 超音波診断装置、超音波診断装置の制御方法、及び超音波診断プログラム | |
JP2017104526A (ja) | 解析装置 | |
JP2012157387A (ja) | 超音波診断装置及び画像生成制御プログラム | |
JP5623160B2 (ja) | 超音波診断装置及びその作動方法 | |
JP6008580B2 (ja) | 超音波診断装置、超音波診断装置の制御方法、及び超音波診断プログラム | |
JP2006095000A (ja) | 超音波撮像装置 | |
US8708910B2 (en) | Ultrasonic diagnosis apparatus and ultrasonic diagnosis method | |
US20120203111A1 (en) | Ultrasonic diagnostic apparatus, ultrasonic image processing apparatus, and ultrasonic image acquisition method | |
JP2013244162A (ja) | 超音波診断装置 | |
JP2006223737A (ja) | 超音波診断装置 | |
JP2012249929A (ja) | 超音波画像生成方法 | |
JP2013244132A (ja) | 超音波診断装置および方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
WWE | Wipo information: entry into national phase |
Ref document number: 201180032787.4 Country of ref document: CN |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 11800882 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012522652 Country of ref document: JP |
|
REEP | Request for entry into the european phase |
Ref document number: 2011800882 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2011800882 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |