WO2013073514A1 - 超音波診断装置及び方法 - Google Patents
超音波診断装置及び方法 Download PDFInfo
- Publication number
- WO2013073514A1 WO2013073514A1 PCT/JP2012/079336 JP2012079336W WO2013073514A1 WO 2013073514 A1 WO2013073514 A1 WO 2013073514A1 JP 2012079336 W JP2012079336 W JP 2012079336W WO 2013073514 A1 WO2013073514 A1 WO 2013073514A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- lattice point
- sound speed
- ultrasonic
- scanning line
- calculating
- 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
-
- 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
Definitions
- the present invention relates to an ultrasonic diagnostic apparatus and method, and in particular, determines the presence or absence of refraction of an ultrasonic scanning line, and based on the determination result, the sound speed (hereinafter referred to as “local sound speed”) in a part (diagnostic site) within a subject. It is related with the technique which calculates accurately.
- the reflection point X1 ROI is obtained.
- the distance L and the speed V can be obtained uniquely.
- the sound speed in the subject when the sound speed in the subject is constant, the sound speed V can be obtained, but when the internal sound speed is not constant as in the subject OBJ2 shown in FIG. In the above method, it is difficult to obtain the distance L to the reflection point (region) X2 ROI and the sound speeds V and V ′.
- Patent Document 1 proposes an ultrasonic diagnostic method capable of accurately calculating the local sound speed even when the sound speed in the subject is not constant.
- an ultrasonic scanning line is emitted from an ultrasonic probe to a subject at a predetermined interval, and a focused scan among received signals obtained by receiving ultrasonic waves reflected by the subject.
- the ambient sound velocity (reception time) which is the average sound velocity of the region from the upper lattice point to the ultrasonic probe, is calculated based on the reception signal of reflection at the lattice point (upper lattice point) set in the region of interest on the line.
- the ultrasonic probe is performed from each lower lattice point.
- the ambient sound speed (reception time) which is the average sound speed of the region up to the child, is calculated, while the assumed sound speed in the region of interest is assumed, and the propagation time from the upper lattice point to each lower lattice point is calculated.
- the incident angle of the ultrasonic wave incident on each lower lattice point from the upper lattice point the assumed sound velocity of the region of interest and the environmental sound velocity calculated in relation to the reflection at the lower lattice point
- the emission angle of the ultrasonic wave emitted from each lower lattice point is calculated, and the position of the ultrasonic probe element where the ultrasonic wave emitted from the lower lattice point with the calculated emission angle is incident and incident on the element.
- the ultrasonic reception time at the position of the element of the ultrasonic probe is calculated by adding the two propagation times, and the ultrasonic reception of the reflection at the calculated reception time and the upper lattice point is calculated.
- the assumed sound speed is corrected so that the error from the reception time at the position of the child element is minimized, and the corrected assumed sound speed is determined as the local sound speed in the region of interest.
- Patent Document 2 discloses a method for obtaining a sound speed distribution by sequentially determining a local sound speed from a shallow area based on focus accuracy in an ultrasonic image.
- an ultrasonic transmission transducer and a reception transducer are arranged at a predetermined distance from each other, and the received vibration from the transmission transducer is changed while changing the transmission and reception angles of these transducers.
- Measure the propagation time of the ultrasonic wave to the child find the error between this measured propagation time and the propagation time of the ultrasonic wave from the transmitting transducer to the receiving transducer calculated separately based on the virtual sound velocity distribution, A method for obtaining the sound speed distribution in the subject by correcting the virtual sound speed distribution so that this error is minimized is disclosed.
- each grid point is defined by the scanning line position and the reception time, and the spatial position is unknown.
- the spatial position of each scanning line is known, and the spatial position is given by approximating the same reception time to the same depth.
- the direction of each scanning line changes due to refraction at the shallower abdominal wall.
- Patent Document 3 The method described in Patent Document 3 has the following problems.
- a dedicated device configuration is required to transmit and receive at a desired angle.
- the present invention has been made in view of such circumstances, and can extract non-refractive scanning lines among sequentially scanned ultrasonic scanning lines, and use only the extracted non-refractive scanning lines. It is an object of the present invention to provide an ultrasonic diagnostic apparatus and method that can accurately calculate the sound speed in a region of interest.
- the accuracy can be improved by using only scanning lines without refraction.
- an ultrasonic diagnostic apparatus transmits an ultrasonic wave to a subject, receives an ultrasonic wave reflected by the subject, and outputs an ultrasonic detection signal
- An ultrasonic probe including a plurality of elements, reception signal acquisition means for acquiring reception signals from a plurality of reflection points based on the ultrasonic detection signal, and a scanning line refracted based on the acquired reception signal
- determining means for determining whether or not it is.
- a scanning line having no refraction among the plurality of scanning lines is extracted based on a reception signal of each reflection on the plurality of scanning lines.
- the received signal is a signal received by a plurality of elements.
- the reception signal acquisition unit acquires reception signals from a plurality of reflection points having different depths, and the determination unit determines the reception signal based on the reception signal. It has a local sound speed calculating means for calculating each local sound speed at a plurality of reflection points, and determines whether or not the scanning line is refracted based on a change in the depth direction of the calculated local sound speed.
- the local sound speed in each region where a plurality of reflection points exist is calculated based on reception signals from a plurality of reflection points having different depths, and the depths of these calculated local sound speeds are calculated.
- a scanning line without refraction is determined based on the change in the vertical direction.
- the determination unit calculates an inclination of the calculated local sound velocity with respect to the depth direction for each of the plurality of scanning lines, and the calculated inclination approximates zero.
- a scanning line that falls within a predetermined threshold is determined as a scanning line without refraction.
- the ultrasonic diagnostic apparatus includes a calculation unit that calculates an average value of local sound velocities calculated for each lattice point on the scanning line determined by the determination unit, and the calculated average The value is the local sound velocity of the region of interest including the reflection point on the scanning line without refraction.
- the reception signal acquisition unit acquires reception signals from a plurality of reflection points having different depths, and the determination unit determines the reception signal based on the reception signal. It has an environmental sound speed calculation means for calculating each environmental sound speed at a plurality of reflection points, and determines whether or not the scanning line is refracted based on a change in the depth direction of the calculated environmental sound speed.
- the reception signal acquisition unit acquires reception signals from a plurality of different scanning line positions, and the determination unit determines the plurality of reception signals based on the reception signals. It has an environmental sound speed calculation means for calculating each environmental sound speed at the scanning line position, and determines whether or not the scanning line is refracted based on a change in the scanning direction of the calculated environmental sound speed.
- the environmental sound speed which is the average sound speed in the region from the reflection point to the ultrasonic probe corresponding to the same reception time on the scanning lines without these refractions, is constant. become. Therefore, in still another aspect of the present invention, received signals from a plurality of different scanning line positions are acquired, and the ambient sound speed at each scanning line position is calculated based on these received signals. A scanning line group without refraction is determined based on the fluctuation in the scanning direction of the scanning line of the calculated environmental sound speed.
- the determination unit is configured to select a continuous scanning line group in which a variation in the scanning direction of the scanning line of the calculated environmental sound speed is within a preset threshold value. It is determined as a scanning line group without refraction.
- the determination unit may calculate the environmental sound speed calculated by the environmental sound speed calculation unit for each predetermined number of scanning line groups among the plurality of scanning lines.
- a standard deviation, variance, or a difference between the maximum value and the minimum value is calculated, and a scanning line group without refraction is determined based on the calculation result.
- the scanning line group is determined as a scanning line group without refraction.
- the reception signal acquisition unit includes lattice points corresponding to reflection points on a plurality of scanning lines, and upper lattice points set in a desired region of interest.
- a reception signal of reflection of a lower lattice point set between the upper lattice point and the ultrasonic probe is acquired, and the local sound speed calculation means is based on the reception signal of reflection at the lower lattice point
- An environmental sound speed calculating means for calculating an environmental sound speed that is an average sound speed in a region from the lower lattice point to the ultrasonic probe, and assuming an assumed sound speed in the region of interest, from the upper lattice point to the lower lattice point.
- an upper lattice point is set in a region of interest, and a lower lattice point is set between the upper lattice point and the ultrasonic probe, so that the sound speed (assumed sound velocity) of the region of interest is set.
- the first propagation time from the upper lattice point to the lower lattice point is calculated.
- the emission angle (refracted) of the ultrasonic wave incident on the lower lattice point at a predetermined incident angle from the upper lattice point is calculated according to Snell's law.
- Snell's law is a law that expresses that there is a fixed relationship between the propagation velocity of each sound wave in two media and the incident and exit angles at the interface between the two media. It is also called. Since the refraction angle at the lower lattice point can be obtained in this way, the position of the element of the ultrasonic probe where the ultrasonic wave is incident from the lower lattice point and the second propagation time until it is incident on the element. And can be calculated.
- the reception time at the element position of the ultrasonic probe obtained by adding the first propagation time and the second propagation time, and the ultrasonic probe of the reflection at the upper lattice point The assumed sound speed when the error from the actual reception time at the position of the child element is minimized is determined as the local sound speed in the region of interest.
- the reception signal acquisition unit includes lattice points corresponding to reflection points on a plurality of scanning lines, and upper lattice points set in a desired region of interest.
- a reception signal of reflection of a lower lattice point set between the upper lattice point and the ultrasonic probe is acquired, and the local sound speed calculation means is configured to reflect the reflection at the upper lattice point and the lower lattice point.
- an environmental sound speed calculating means for calculating an environmental sound speed that is an average sound speed in an area from each lattice point to the ultrasonic probe, and a first reception when the upper lattice point is a reflection point
- a first calculation means for calculating a wave based on the environmental sound speed calculated corresponding to the upper lattice point, and an assumed sound speed in the region of interest, and a propagation time from the upper lattice point to the lower lattice point
- Means for calculating and a second received wave from the lower grid point
- a second calculating means for calculating based on the environmental sound velocity calculated corresponding to the child point and the calculated propagation time; the first received wave calculated by the first calculating means; and the second calculating means.
- Local sound speed determination means for determining, as the local sound speed in the region of interest, the hypothetical sound speed with which the error from the second received wave calculated by the above is minimized.
- an upper lattice point is set in a region of interest, and a lower lattice point is set between the upper lattice point and the ultrasonic probe, so that the sound speed (assumed sound velocity) of the region of interest is set.
- the ambient sound velocity which is the average sound velocity in the region from each lattice point to the ultrasonic probe, is calculated based on the reception signals reflected at the upper lattice point and the lower lattice point. Then, the first received wave when the upper grid point is used as a reflection point is calculated based on the calculated environmental sound speed corresponding to the upper grid point.
- the propagation time from the upper lattice point to each lower lattice point is calculated, and based on this propagation time and the ambient sound velocity calculated corresponding to the lower lattice point, A second received wave from the lattice point is calculated. Then, the assumed sound speed when the calculated error between the first received wave and the second received wave is minimized is determined as the local sound speed in the region of interest. This utilizes the fact that the received wave from the upper lattice point and the received wave from the lower lattice point coincide with each other by Huygens' principle.
- the received signal acquisition means may detect an upper refraction point corresponding to a reflection point on the scan line when the determination means determines a scan line without refraction.
- a reception signal of reflection at a lower lattice point is acquired, and the local sound speed calculation means calculates a local sound speed in the region of interest based on the acquired reception signal.
- the received signal acquisition unit is a lattice point corresponding to a reflection point on the scan line when the determination unit determines a scan line having no refraction.
- a reception signal of reflection of an upper lattice point set in a desired region of interest and a lower lattice point set between the upper lattice point and the ultrasonic probe, and at the lower lattice point Based on the received reflection signal, the environmental sound speed calculating means for calculating the environmental sound speed, which is the average sound speed of the area from the lower lattice point to the ultrasonic probe, and the assumed sound speed in the region of interest, Means for calculating a first propagation time from the upper lattice point to the lower lattice point; an incident angle of an ultrasonic wave incident on the lower lattice point from the upper lattice point according to Snell's law; an assumed sound velocity of the region of interest; The ring of the region between the lower lattice point
- the received signal acquisition unit is a lattice point corresponding to a reflection point on the scan line when the determination unit determines a scan line having no refraction.
- a first calculation means for calculating the first received wave based on the environmental sound speed calculated corresponding to the upper lattice point, and an assumed sound speed in the region of interest; Means for calculating the propagation time to the lattice point, and from the lower lattice point A second calculating means for calculating a second received wave based on the environmental sound
- the received signal acquisition unit is a lattice point corresponding to a reflection point on the scan line when the determination unit determines a scan line having no refraction.
- the determination means Propagation time calculating means for calculating the propagation time of ultrasonic waves between the upper and lower lattice points on the scanned line without refraction, based on the received signal acquired by the received signal acquiring means, and the determining means Based on the positions of the upper and lower lattice points on the scanning line without refraction determined by the above and the propagation time of the ultrasonic wave between the upper and lower lattice points calculated by the propagation time calculating means Wearing between the upper and lower grid points Further comprising a local sound velocity calculation means for calculating the local sound velocity, the
- the local sound speed is accurately calculated by obtaining a reception signal of reflection only from lattice points on a scanning line having no refraction, and obtaining a position and propagation time between these lattice points. I can do it.
- the propagation time calculation means calculates a first propagation time from an upper lattice point on the target scanning line to each element of the ultrasonic probe.
- a time difference on the element that maximizes the time difference between the calculated first propagation time and the second propagation time on each element of the second propagation time calculating means and the ultrasonic probe is calculated.
- means for calculating the propagation time of ultrasonic waves from the upper lattice point to the lower lattice point is calculated.
- a first propagation time from the upper lattice point to each element of the ultrasonic probe and a first propagation time from the lower lattice point to each element of the ultrasonic probe can be easily calculated based on the propagation time of 2.
- the ultrasonic diagnostic apparatus further includes a steer angle adjustment unit that adjusts a steer angle of a scanning line transmitted and received from the ultrasonic probe, and the determination unit includes the steer angle adjustment unit.
- the steering angle is adjusted, it is determined whether or not the scanning line is refracted based on the acquired received signal.
- the scanning line is incident so as to be substantially orthogonal to the boundary surface of different media, the refraction of the scanning line becomes small. Accordingly, it is possible to transmit a scanning line in which the refraction of the scanning line is reduced by adjusting the steering angle of the scanning line.
- the ultrasonic diagnostic apparatus further includes display means for displaying a scanning line having no refraction determined by the determination means.
- An ultrasonic diagnostic method transmits an ultrasonic wave from an ultrasonic probe including a plurality of elements to a subject and receives an ultrasonic wave reflected by the subject.
- the present invention it is possible to easily determine a scanning line without refraction based on a received signal of each element of an ultrasonic probe that can be acquired by a normal apparatus configuration, and only a scanning line without refraction.
- the effect that the sound speed (local sound speed) in the region of interest can be calculated with high accuracy, and that the displacement detection accuracy in the azimuth direction can be improved in functions such as transverse wave speed measurement and lateral speckle tracking. effective.
- FIG. 1 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus according to the present invention.
- the figure used to explain the steer angle of the ultrasonic beam emitted from the ultrasonic probe The flowchart which shows 1st Embodiment of the process sequence which calculates the local sound speed in the attention area
- the figure used in order to demonstrate the modification in the case of calculating a local sound speed by the 1st calculation method or the 2nd calculation method A graph showing changes in local sound speed with respect to the depth direction of biological phantoms with different media
- FIG. 1 is a block diagram showing an embodiment of an ultrasonic diagnostic apparatus according to the present invention.
- the ultrasonic diagnostic apparatus 10 shown in FIG. 1 transmits an ultrasonic beam from the ultrasonic probe 300 to the subject OBJ, receives an ultrasonic beam (ultrasonic echo) reflected by the subject OBJ, and This is an apparatus for creating and displaying an ultrasonic image from a detection signal of a sound echo.
- a CPU (Central Processing Unit) 100 controls each block of the ultrasonic diagnostic apparatus 10 according to 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 has 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 speed.
- Display mode switching button for switching, freeze button for instructing switching between live mode and freeze mode, cine memory playback button for instructing cine memory playback, analysis for instructing analysis / measurement of ultrasonic images Includes a measurement button.
- 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 extracts a control program for controlling the control of each block of the ultrasonic diagnostic apparatus 10 by the CPU 100, a parameter, and a program for calculating a local sound speed by extracting a scanning line without refraction according to the present invention.
- a control program for controlling the control of each block of the ultrasonic diagnostic apparatus 10 by the CPU 100 a parameter, and a program for calculating a local sound speed by extracting a scanning line without refraction according to the present invention.
- a hard disk or a semiconductor memory 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), a scanning line direction, a local sound velocity map, and various setting screens according to the present invention. indicate.
- CTR Cathode Ray Tube
- ultrasonic image moving image and still image
- scanning line direction a scanning line direction
- local sound velocity map a local sound velocity map
- various setting screens according to the present invention. indicate.
- the ultrasonic probe 300 is a probe used in contact with the subject OBJ, and includes a plurality of elements 302 constituting a one-dimensional or two-dimensional ultrasonic transducer array.
- the plurality of elements 302 transmit an ultrasonic beam to the subject OBJ based on the drive signal applied from the transmission circuit 402, receive an ultrasonic echo reflected from the subject OBJ, and output a detection signal.
- Each element 302 of the ultrasonic probe 300 includes a vibrator formed by forming electrodes on both ends of a piezoelectric material (piezoelectric body).
- piezoelectric body constituting the vibrator
- the piezoelectric body constituting the vibrator include piezoelectric ceramics such as PZT (lead zirconate titanate) and polymer piezoelectric elements such as PVDF (polyvinylidene difluoride).
- PZT lead zirconate titanate
- PVDF polyvinylidene difluoride
- 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 element 302 of the ultrasonic probe 300 it is also 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 element 302.
- the CPU 100 selects a transmission delay pattern according to the transmission direction of the ultrasonic beam and also 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 the drive signal in order to form an ultrasonic beam in a desired direction by ultrasonic waves transmitted from the plurality of elements 302
- the reception delay pattern is The pattern data of the delay time given to the detection signal in order to extract the ultrasonic echo from the desired direction by the ultrasonic waves received by the plurality of elements 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 element 302.
- the transmission circuit 402 has delay circuits ⁇ 1 to ⁇ N for each element 302 as shown in FIG. 2, and delays the drive signal applied to each element 302 based on the transmission delay pattern selected by the CPU 100.
- the transmission circuit 402 adjusts (delays) the timing at which the drive signal is applied to each element 302 so that the ultrasonic waves transmitted from the plurality of elements 302 form an ultrasonic beam, as shown in FIG.
- the timing for applying the drive signal to each element 302 is adjusted (delayed) so as to adjust the direction of the ultrasonic beam (steer angle ⁇ ).
- the timing of applying the drive signal may be adjusted so that the ultrasonic waves transmitted from the plurality of elements 302 at a time reach the entire imaging region of the subject OBJ.
- the receiving circuit 404 receives and amplifies an ultrasonic detection signal output from each element 302 of the ultrasonic probe 300. As described above, since the distance between each element 302 and the ultrasonic wave reflection source in the subject OBJ is different, the time for the reflected wave to reach each element 302 is different.
- the reception circuit 404 includes a delay circuit, and a difference (delay in arrival time) of the reflected wave according to a reception delay pattern set based on a sound speed selected by the CPU 100 (hereinafter referred to as “assumed sound speed”) or a sound speed distribution. Each detection signal is delayed by an amount corresponding to (time).
- the reception circuit 404 performs reception focus processing by matching and adding detection signals given delay times.
- the arrival time of the ultrasonic detection signal from the other ultrasonic reflection source is different.
- the phases of the ultrasonic detection signals from other ultrasonic reflection sources cancel each other.
- the received signal from the ultrasonic wave reflection source X ROI becomes the largest and is focused.
- RF signal a sound ray 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 due to 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 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 is generated.
- 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, television signal scan system (NTSC system image data)).
- 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 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 elements 302 are amplified in the reception circuit 404, and the amplified detection signals, that is, RF signals are A / D converted in the A / D converter 406.
- RF data is generated.
- 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 according to 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 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 by the data analysis / measurement unit 106 into the image data of the ultrasonic image and outputs it 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, a mode for displaying the determination result of the local sound speed superimposed on the B mode image (for example, color coding or luminance according to the local sound speed).
- the display mode is switched between a mode in which a B-mode image and a local sound speed value determination result image are displayed side by side, or a display in which the local sound speed is equalized by a line.
- the operator can find a lesion, for example, by observing the determination result of the local sound speed.
- a B-mode image obtained by performing at least one of the transmission focus process and the reception focus process may be displayed on the display unit 104 based on the determination result of the local sound speed.
- the ultrasonic diagnostic apparatus 10 uses an ultrasonic beam (hereinafter referred to as “scan line”) based on reception signals at each element 302 of the ultrasonic probe 300. )) Is extracted.
- scan line an ultrasonic beam
- a method for extracting a scanning line without refraction will be described later.
- the display unit 104 can display a scanning line without refraction.
- FIG. 3 is a flowchart showing a first embodiment of a processing procedure for calculating the local sound velocity in the region of interest of the subject.
- a region of interest of the subject is set (step S1).
- This region of interest may be set by an operator using a pointing device on a still image of an ultrasonic image displayed on the display unit 104, or may be automatically set at a predetermined position and a predetermined size by a control program.
- binarization processing may be performed on the ultrasonic image
- labeling processing may be performed in which the same numbers are assigned to pixels in which white portions (or black portions) are continuous, and the settings may be automatically performed in the order of the labeled numbers.
- lattice points including the upper lattice point and the lower lattice point
- the ambient sound speed at each lattice point is obtained (step S2).
- the position of each grid point is defined by the scanning line position and the reception time. That is, as shown in FIG. 4, reflection points with different depths on the scanning lines 1, 2,..., N, which are emitted from the ultrasound probe 300 to the region of interest of the subject OBJ, are latticed. Set as a point.
- each scanning lines 1, ..., reception time on n are the same reflection point, likewise the upper grid points B 1 , B 2, B 3, ..., B n and upper grid points C 1, C 2, C 3, ..., C n , ... are also reflected with the same reception time on each scanning line 1, 2, ..., n. Is a point.
- the lower lattice point A and the upper lattice points B and C are illustrated as lattice points having the same depth, but in actuality, between each lattice point and the ultrasonic probe 300. Since the sound speed of the area is not uniform, it becomes a reflection point of different depth in the space, and each scanning line 1, 2, ..., n that is linearly scanned is also refracted due to the difference in the sound speed of the propagation area of the scanning line, All scan lines are not necessarily parallel.
- the range and number of each grid point are determined in advance.
- the range of the lattice point used for the local sound speed calculation is wide, the error of the local sound speed value becomes large, and if it is narrow, the error with the virtual received wave becomes large. Therefore, the range of the lattice point is determined based on these factors.
- the interval between the lattice points in the x direction is determined by the balance between the resolution and the processing time.
- the interval between the lattice points in the x direction is, for example, 1 mm to 1 cm.
- the distance between the lower lattice point and the upper lattice point in the y direction is narrow, the error in error calculation becomes large, and if it is wide, the error in local sound speed becomes large.
- the interval between the lattice points in the y direction is determined based on the setting of the image resolution of the ultrasonic image, and is 1 cm as an example.
- the lattice point X1 ROI is measured.
- the distance L and the speed V can be obtained uniquely.
- the environmental sound speed at a certain grid point is the average sound speed in the region from the grid point to the ultrasonic probe, and is the sound speed at which the contrast and sharpness of the image are the highest. Therefore, as a method for determining the environmental sound speed, for example, a method (for example, Japanese Patent Laid-Open No. 8-317926) for determining from the contrast of the image, the spatial frequency in the scanning direction, dispersion, and the like can be applied. Further, it is preferable to focus so as to form transmission focal points at fine depth intervals in the region of interest so that the environmental sound speed can be obtained with sufficient accuracy.
- the reception time of each element reflected from the lattice point can be obtained from the environmental sound velocity obtained in this way. That is, for each element reception signal reflected from the lattice point, a delay is determined assuming a certain sound speed, and when the contrast and sharpness of an image generated using the delay are the highest, the delay is This means that the reception time of each element is closest, so that it can be regarded as the reception time of each element with its sound speed (environmental sound speed), that is, with a delay. Instead of the environmental sound speed, the reception time of each element may be obtained by a method such as phase aberration analysis and used thereafter.
- the propagation times ⁇ T1, ⁇ T2, and ⁇ T3 from the lattice point (upper lattice point) B on the scanning line of interest to the lattice points (lower lattice points) A1, A2, and A3 on the peripheral scanning line are Can be calculated.
- ⁇ Tn max (TBi-TAni)
- ⁇ Tn Propagation time from lattice point B to lattice An
- TBi Propagation time from lattice point B to element i of the ultrasonic probe (reception time of lattice point B reflection at element i)
- TAni propagation time from the lattice point An to the element i of the ultrasonic probe (reception time of the lattice point An reflection at the element i) It is.
- TBi and TAni indicate the propagation time of one way of the propagation path.
- the propagation time of each element of the ultrasonic probe or the propagation time from the environmental sound speed (the reception time or the shortest time at the element on the target scanning line). It is obtained by subtracting half of the reception time.
- FIG. 5 shows a curve indicating a reception time at each element i of the reflection ultrasonic probe at the lattice point B, and each element i of the reflection ultrasonic probe at the reflection of the lattice points A1, A2, and A3.
- a curve indicating the reception time is shown.
- the reception wave from the lattice point B (curve indicating the reception time) and the reception wave from the lattice points A1, A2, A3,... are propagated only from the lattice point B to each lattice point. This coincides with a virtual composite received wave (curve envelope indicating each reception time) which is virtually combined with delay.
- ⁇ Tn calculated by the above [Expression 2] is a grid point from the grid point B necessary for the received wave from the grid point B and the virtual synthesized received wave from the grid points A1, A2, A3.
- the propagation time to An is shown.
- the position of the lower lattice points A1, A2, A3, instead of obtaining each propagation time ⁇ Tn independently as in the formula [2], the position of the lower lattice points A1, A2, A3,.
- Each propagation time ⁇ Tn can be obtained based on the positions of the lower lattice points A1, A2, A3.
- the positions of the upper lattice point B and the lower lattice points A1, A2, A3, Since the propagation time from the upper lattice point to the lower lattice point is given in the scanning line of interest, the depth of the upper lattice point and the lower lattice point is assumed by assuming the local sound speed between the upper lattice point and the lower lattice point. A distance in the vertical direction is given, so that the propagation time from the upper lattice point B to each of the lower lattice points A1, A2, A3. Each propagation time is compared with each propagation time ⁇ Tn obtained above, and the assumed local sound speed when the error is minimized is determined as the true sound speed (local sound speed).
- the local sound velocity can be calculated for each scanning line.
- received signals from the upper grid points C, D,... Having different depth directions are used (see FIG. 4), and the local points at the grid points having different depth directions are used in the same manner as described above. Calculate the speed of sound.
- FIG. 6 is a diagram schematically showing a method of calculating the local sound speed by the refraction model calculation disclosed in Patent Document 1.
- the direction parallel to the element surface S2 on which each element 302 of the ultrasonic probe 300 is arranged is defined as the X direction
- the direction perpendicular to the X direction is defined as the Y direction.
- the upper grid point representing the region of interest ROI in the region A in the subject OBJJ is set as B ROI
- the lower grid points are set as A 1, A 2,.
- the spatial coordinates of these lattice points are given on the assumption that each scanning line is not refracted.
- a region between the boundary surface S1 connecting the lower lattice points A1, A2,..., An,... And the upper lattice point B ROI in the subject OBJ is defined as a region A, and the boundary surface S1 and the ultrasonic probe 300 are A region between the element surface S2 is a region B. It is assumed that the sound speeds in the regions A and B are constant.
- the reception time at each element 302 is obtained by tracing the sound ray refracted at the boundary surface between the regions A and B according to Snell's law.
- the sound ray passing through the lower lattice point X ′ Can be expressed by the following equation according to Snell's law.
- the propagation time from the upper lattice point B ROI to the lower lattice points A1, A2, can be calculated by assuming the assumed sound velocity V A because the distance between the respective lattice points can be obtained.
- the calculated reception time and the measured reception time are calculated.
- the assumed sound speed when the error is minimized is determined as the true sound speed (local sound speed) in the region of interest.
- the local sound velocity can be calculated for each scanning line.
- received signals from the upper grid points C, D,... Having different depth directions are used (see FIG. 4), and the local points at the grid points having different depth directions are used in the same manner as described above. Calculate the speed of sound.
- FIG. 7 is a diagram schematically showing a method of calculating the local sound speed using the Huygens principle disclosed in Patent Document 1.
- the propagation point T and the delay time ⁇ T of the received waves (W A1 , W A2 ,...) From the lower lattice points A1, A2 ,.
- the local sound speed at the lattice point B ROI is obtained from the positional relationship between the lattice points A1, A2 ,.
- the Huygens principle indicates that the received wave W X from the upper grid point B ROI matches the received wave W SUM virtually combined with the received waves from the lower grid points A1, A2 ,. Use.
- the spatial coordinates of the upper lattice point B ROI and the lower lattice points A1, A2,... are given on the assumption that each scanning line is not refracted.
- an error between the received wave W X and the combined received wave W SUM is calculated.
- Error between the received wave W X resultant received wave W SUM a method, a method of phase matching addition is multiplied by the delay obtained from the resultant received wave W SUM to the receiving wave W X, or synthetic reception reversed cross-correlating with each other It is calculated by a method of multiplying the wave W SUM by the delay obtained from the received wave W X and adding the phase matching.
- the lattice point B ROI is used as a reflection point, and the time at which the ultrasonic wave propagated at the sound velocity V arrives at each element may be set as the delay.
- an equiphase line is extracted from the phase difference of the combined reception wave between adjacent elements, and the equal phase line is used as a delay, or simply a combination of each element.
- 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 value of the amplitude after the envelope detection.
- the error between the received wave W X and the synthesized received wave W SUM varies depending on the assumed sound speed V A. Then, the assumed sound speed when the error is minimum (maximum during phase matching addition) is determined as the true sound speed (local sound speed) in the region of interest.
- the local sound velocity can be calculated for each scanning line.
- received signals from the upper grid points C, D,... Having different depth directions are used (see FIG. 4), and the local points at the grid points having different depth directions are used in the same manner as described above. Calculate the speed of sound.
- the local sound velocities at the upper lattice points B, C, D,... Calculated by the first to third calculation methods are the local sounds in the region between the upper lattice points and the lower lattice points A1, A2,. It corresponds to the speed of sound, that is, the lower grid points are shared, and there is an overlap in each area. However, the lower grid points may be set separately for the upper grid points so that the areas do not overlap. At this time, a grid point common to each upper grid point is separately set at a reference depth with respect to each upper grid point, and the local sound velocity at each upper grid point is matched so that the environmental sound speeds at this grid point match. The speed of sound may be calculated.
- FIGS. 9A to 9C are graphs showing the results of measuring the local sound speed at the depth of the living body phantom on the premise that the scanning line of interest and the surrounding scanning lines are not refracted, each being shallower than the measurement point. Different sound velocity media are installed in different regions with different boundary shapes. In these graphs, it is the influence of noise that the amplitude fluctuates greatly.
- the average value of the local sound speeds on the scanning line determined as having no refraction as described above is calculated, and the calculated average value is set as the local sound speed in the region of interest.
- the local sound speed is calculated by applying the local sound speed calculation method disclosed in Patent Document 1 again. It may be.
- the environmental sound speed of each scanning line without refraction may be averaged, the local sound speed may be obtained based on the depth profile, the received signal or the focus index is averaged, and the environmental sound speed is obtained based thereon, Local sound speed may be obtained.
- FIG. 10 is a flowchart showing a second embodiment of the processing procedure for calculating the local sound velocity in the region of interest of the subject.
- the same step number is attached
- steps S3 'and S4' are different from the processes in steps S3 and S4 in the first embodiment.
- step S3 ' the fluctuation of the environmental sound speed in the scanning direction in a predetermined scanning line width is calculated based on the environmental sound speed for each scanning line obtained in step S2. At this time, a plurality of environmental sound speeds obtained in the depth direction of the scanning line may be averaged, and the average environmental sound speed fluctuation may be calculated.
- FIG. 11 is a graph showing an example of a change in environmental sound speed with respect to the scanning line direction of the scanning line.
- step S3 ' the standard deviation, variance, or difference between the maximum value and the minimum value of continuous environmental sound speeds corresponding to a preset number of scanning lines is calculated as continuous environmental sound speed fluctuation information.
- the fluctuation of the environmental sound speed in the scanning line direction at the lattice point of the region of interest may be included.
- step S4 ' a scanning line group in which the fluctuation of the environmental sound speed obtained in step S3' is equal to or less than a preset threshold is extracted, and this scanning line group is determined as a scanning line group without refraction.
- a scanning line without refraction is extracted in this way, only the received signals from the upper and lower lattice points on the scanning line without refraction are used, and the refraction model calculation disclosed in Patent Document 1
- the local sound speed is calculated using the principle or the like.
- the steering angle ⁇ can be adjusted by delaying the drive signal applied from the transmission circuit 402 to the element i of the ultrasonic probe 300 by the delay time ⁇ i shown in the following equation.
- reception focus is performed so as to form a focus at each depth in the direction of the steer angle ⁇ , and the ambient sound speed at each lattice point is obtained.
- V sound velocity (eg, known sound velocity in subcutaneous fat, etc.)
- the method of the first embodiment or the second embodiment is applied, and the scanning line group without refraction or the scanning line with low refraction is obtained.
- a group can be determined.
- the region of interest can also be applied to objects with non-uniform sound speed by making the region of interest small.
- the local sound velocity may be obtained independently in each region, the result of the region close to the ultrasonic probe (shallow region) may be used.
- the scan line group without refraction may be determined including the determination result without refraction of the scan line in the shallow region.
- each scanning line indicating the extracted scanning line direction without refraction may be displayed on the display unit 104.
- the scanning line without refraction is not limited to the case where the scanning line is not refracted at all, but includes a scanning line with low refraction.
- the change in the scanning line direction is 1 depending on the accuracy of the required local sound velocity.
- SYMBOLS 10 Ultrasound diagnostic apparatus, 100 ... Central processing unit (CPU), 102 ... Storage part, 104 ... Display part, 106 ... Data analysis measurement part, 200 ... Operation input part, 202 ... Console, 204 ... Pointing device, 300 DESCRIPTION OF SYMBOLS ... Ultrasonic probe, 302 ... Ultrasonic transducer (element), 400 ... Transmission / reception part, 402 ... Transmission circuit, 404 ... Reception circuit, 500 ... Image signal generation part, 502 ... Signal processing part, 506 ... Image processing part, 508: Image memory, 510: D / A converter, 600: Playback unit
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
複数の素子を含む超音波探触子から超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信し、深さの異なる複数の反射点からの受信信号を取得し、前記複数の反射点が存在する各領域における局所音速を算出する。これらの算出した局所音速の深さ方向に対する変化(傾き)を算出し、前記算出した傾きによりゼロ近傍の所定の閾値以内となる走査線を、屈折のない走査線として判定する。
Description
本発明は超音波診断装置及び方法に係り、特に超音波走査線の屈折の有無を判定し、その判定結果に基づいて被検体内の一部(診断部位)における音速(以下、「局所音速」という)を精度よく算出する技術に関する。
音速が一定の媒質からなる被検体OBJ1内の音速Vは、下記のようにして算出することができる。図12(a)に示すように、被検体OBJ内の反射点(領域)X1ROIから超音波探触子300Aまでの距離をLとすると、反射点X1ROIで超音波が反射されてから反射点X1ROIの直下の超音波トランスデューサ(素子)302A0で受信されるまでの経過時間Tは、T=L/Vである。素子302A0からX方向(素子302Aの配列方向)に距離X離れた位置にある素子302Aiで受信されるまで経過時間をT+ΔTとすると、素子302A0と302Aiとの間の遅延時間ΔTは、下記の[数1]式により表される。
従って、超音波が送波されて反射点X1ROIで時間T後に反射された後、各素子により受信されるまでの経過時間[2T,2T+ΔT]を測定することにより、反射点X1ROIまでの距離Lと速度Vを一意に求めることができる。
上記のように、被検体内の音速が一定の場合には、音速Vを求めることが可能であるが、図12(b)に示す被検体OBJ2のように、内部の音速が一定でない場合には、上記の方法では、反射点(領域)X2ROIまでの距離L及び音速V,V´を求めることは困難である。
これに対し、被検体内の音速が一定でない場合であっても、局所音速を精度よく算出することができる超音波診断方法が提案されている(特許文献1)。
この超音波診断方法は、超音波探触子から超音波走査線を所定の間隔で被検体に出射し、被検体によって反射される超音波を受信して得た受信信号のうち、着目する走査線上の着目領域に設定された格子点(上格子点)での反射の受信信号に基づいて上格子点から超音波探触子までの領域の平均音速である環境音速(受信時刻)を算出するとともに、前記上格子点と前記超音波探触子との間に設定された各走査線上の格子点(下格子点)での反射の受信信号に基づいて、各下格子点から超音波探触子までの領域の平均音速である環境音速(受信時刻)を算出し、一方、前記着目領域における仮定音速を仮定し、上格子点から各下格子点までの伝播時間を算出する。
また、スネルの法則により前記上格子点から各下格子点に入射する超音波の入射角と、前記着目領域の仮定音速と前記下格子点での反射に関連して算出した環境音速とに基づいて各下格子点から出射する超音波の出射角を算出し、前記下格子点から前記算出した出射角で出射する超音波が入射する前記超音波探触子の素子の位置と該素子に入射するまでの伝播時間とを算出する。
そして、超音波探触子の素子の位置における超音波の受信時刻を、前記2つの伝播時間を加算することにより算出し、この算出した受信時刻と前記上格子点での反射の超音波探触子の素子の位置における受信時刻との誤差が最小になるように前記仮定音速を修正し、その修正した仮定音速を着目領域における局所音速として判定するようにしている。
また、特許文献2には、超音波画像におけるフォーカス精度に基づき、浅い領域から局所的な音速を順次決定することにより音速分布を求める方法が開示されている。
特許文献4には、超音波送波振動子と受波振動子とを所定距離離間して配置し、これらの振動子の送波及び受波角度を変えながら、送波振動子から受波振動子までの超音波の伝播時間を計測し、この計測した伝播時間と、別途仮想音速分布に基づいて計算した送波振動子から受波振動子までの超音波の伝播時間との誤差を求め、この誤差が最小となるように仮想音速分布を修正することにより被検体内の音速分布を求める方法が開示されている。
特許文献1に記載の超音波診断方法は、局所音速の測定の前提として、着目領域に設定される上格子点に対し、浅い領域に設定される各下格子点の相対的な空間位置が既知でなければならない。
しかしながら、各格子点は走査線位置と受信時刻とによりその位置が定義され、空間位置は未知である。これに対し、特許文献1においては、各走査線の空間位置を既知とし、また、同一の受信時刻を同一深さと近似する事で、空間位置を与えている。
したがって、走査線方向が屈折により変化すると、各格子点の空間位置が方位方向にも、深さ方向にも変化してしまい、算出される局所音速に誤差が生じるという問題がある。
例えば、測定対象が肝臓の場合、それよりも浅い腹壁における屈折により各走査線方向が変化する。
また、特許文献2に記載の方法においても、超音波画像におけるフォーカスは音速のみでなく屈折にも依存するため、屈折を求めずに音速分布を求めることはできない。
特許文献3に記載の方法は、以下の課題がある。
(1)所望の角度で送波・受波できるための専用の装置構成が必要である。
(2)多大な送波・受波回数、及び仮想音速分布を探索するための多大な処理時間が必要になる。着目領域のみの音速を求めるためでも、他の全領域を伝播した結果の時間の計測値と計算値とを比較する原理のため、多大な送・受回数及び処理時間が必要なことに変わりない。
(3)伝播経路を算出するために、音速場を線型近似した三角形に展開するモデルにて与えているが、実際の音速場はより複雑であり、特許文献3に記載のような近似モデルにて解を得ることは困難である。また、もし複雑な音速場を再現できるモデルを与えたとしても、特許文献3に記載の発明のように限られた送波・受波データから解を得ることは困難である。
本発明はこのような事情に鑑みてなされたもので、順次走査される超音波走査線のうち屈折のない走査線を抽出することができ、また、その抽出した屈折のない走査線のみを使用して着目領域における音速を精度よく算出することができる超音波診断装置及び方法を提供することを目的とする。
尚、横波音速計測や横方向のスペックルトラッキング、画像歪み補正など、方位方向の変位情報を要する何れの機能に対しても、屈折のない走査線のみを用いる事により精度改良できるため、これも目的とする。
前記目的を達成するために本発明の一の態様に係る超音波診断装置は、超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して超音波検出信号を出力する複数の素子を含む超音波探触子と、前記超音波検出信号に基づいて複数の反射点からの受信信号を取得する受信信号取得手段と、前記取得した受信信号に基づいて走査線が屈折しているか否かを判定する判定手段と、を備えている。
本発明の一の態様によれば、複数本の走査線上の各反射の受信信号に基づいて前記複数本の走査線のうちの屈折のない走査線を抽出するようにしている。
本発明の他の態様に係る超音波診断装置において、前記受信信号は、複数の素子が受信した信号である。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、複数の異なる深さの反射点からの受信信号を取得し、前記判定手段は、該受信信号に基づいて該複数の反射点における各局所音速を算出する局所音速算出手段を有し、前記算出した局所音速の深さ方向の変化に基づいて走査線が屈折しているか否かを判定する。
走査線上の深さ方向の各領域における局所音速が、深さにかかわらず変動しない場合、その走査線は屈折がないと認められる。そこで、本発明の更に他の態様では、複数の異なる深さの反射点からの受信信号に基づいて複数の反射点が存在する各領域における局所音速を算出し、これらの算出した局所音速の深さ方向の変化に基づいて屈折のない走査線を判定するようにしている。
本発明の更に他の態様に係る超音波診断装置において、前記判定手段は、前記複数の走査線毎に前記算出した局所音速の深さ方向に対する傾きを算出し、前記算出した傾きによりゼロ近傍の所定の閾値以内となる走査線を、屈折のない走査線として判定する。
本発明の更に他の態様に係る超音波診断装置において、前記判定手段により判定された走査線上の格子点毎に算出された局所音速の平均値を算出する算出手段を備え、前記算出された平均値を前記屈折のない走査線上の反射点を含む着目領域の局所音速とする。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、複数の異なる深さの反射点からの受信信号を取得し、前記判定手段は、該受信信号に基づいて該複数の反射点における各環境音速を算出する環境音速算出手段を有し、前記算出した環境音速の深さ方向の変化に基づいて走査線が屈折しているか否かを判定する。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、複数の異なる走査線位置からの受信信号を取得し、前記判定手段は、該受信信号に基づいて該複数の走査線位置における各環境音速を算出する環境音速算出手段を有し、前記算出した環境音速の走査方向の変化に基づいて走査線が屈折しているか否かを判定する。
連続している複数の走査線に屈折がない場合、これらの屈折のない走査線上の同一の受信時刻に対応する反射点から超音波探触子までの領域の平均音速である環境音速は、一定になる。そこで、本発明の更に他の態様では、複数の異なる走査線位置からの受信信号を取得し、これらの受信信号に基づいて各走査線位置の環境音速をそれぞれ算出する。そして、算出した環境音速の走査線の走査方向の変動に基づいて屈折のない走査線群を判定するようにしている。
本発明の更に他の態様に係る超音波診断装置において、前記判定手段は、前記算出した環境音速の前記走査線の走査方向の変動が、予め設定された閾値以内となる連続する走査線群を、屈折のない走査線群として判定する。
本発明の更に他の態様に係る超音波診断装置において、前記判定手段は、前記複数の走査線のうちの連続する所定数の走査線群毎に、前記環境音速算出手段により算出した環境音速の標準偏差、分散又は最大値と最小値の差分を算出し、その算出結果に基づいて屈折のない走査線群を判定する。連続する所定数の走査線群毎に算出した環境音速の標準偏差等がゼロに近い場合、その走査線群は屈折のない走査線群として判定される。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、複数の走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、前記局所音速算出手段は、前記下格子点での反射の受信信号に基づいて、該下格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの第1の伝播時間を算出する手段と、スネルの法則により前記上格子点から下格子点に入射する超音波の入射角と、前記着目領域の仮定音速と前記下格子点と前記超音波探触子との間の領域の環境音速とに基づいて前記下格子点から出射する超音波の出射角を算出する手段と、前記下格子点から前記算出した出射角で出射する超音波が入射する前記超音波探触子の素子の位置と該素子に入射するまでの第2の伝播時間とを算出する手段と、前記超音波探触子の素子の位置における超音波の受信時刻を、前記第1の伝播時間と第2の伝播時間とを加算して算出する手段と、前記上格子点での反射の前記超音波探触子の素子の位置における受信時刻と前記算出した受信時刻との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、を含む。
本発明の更に他の態様は、着目領域に上格子点を設定するとともに、前記上格子点と前記超音波探触子との間に下格子点を設定し、着目領域の音速(仮定音速)を仮定することにより、上格子点から下格子点までの第1の伝播時間を算出する。一方、前記仮定音速と、格子点から超音波探触子までの領域の平均音速である環境音速とから、上格子点から所定の入射角で下格子点に入射する超音波の出射角(屈折角)をスネルの法則により算出する。尚、スネルの法則は、2つの媒質中における音波のそれぞれの伝播速度と2つの媒質の境界面での入射角・出射角とは一定の関係にあることを表した法則であり、屈折の法則とも呼ばれている。このようにして下格子点での屈折角を求めることができるため、下格子点から超音波が入射する前記超音波探触子の素子の位置と該素子に入射するまでの第2の伝播時間とを算出することができる。そして、前記第1の伝播時間と第2の伝播時間とを加算することにより求めた前記超音波探触子の素子の位置における受信時刻と、前記上格子点での反射の前記超音波探触子の素子の位置における実際の受信時刻との誤差が最小となるときの仮定音速を、前記着目領域における局所音速として判定するようにしている。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、複数の走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、前記局所音速算出手段は、前記上格子点及び下格子点での反射の受信信号に基づいて、各格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、前記上格子点を反射点としたときの第1の受信波を、該上格子点に対応して算出した環境音速に基づいて算出する第1の算出手段と、前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの伝播時間を算出する手段と、前記下格子点からの第2の受信波を、該下格子点に対応して算出した環境音速及び前記算出した伝播時間に基づいて算出する第2の算出手段と、前記第1の算出手段により算出された第1の受信波と前記第2の算出手段により算出された第2の受信波との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、を含む。
本発明の更に他の態様は、着目領域に上格子点を設定するとともに、前記上格子点と前記超音波探触子との間に下格子点を設定し、着目領域の音速(仮定音速)を仮定し、また、上格子点及び下格子点での反射の受信信号に基づいて各格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する。そして、前記上格子点を反射点としたときの第1の受信波を該上格子点に対応した算出した環境音速に基づいて演算する。一方、前記着目領域における仮定音速を仮定し、前記上格子点から各下格子点までの伝播時間を算出し、この伝播時間と下格子点に対応して算出した環境音速とに基づいて、下格子点からの第2の受信波を算出する。そして、前記算出された第1の受信波と第2の受信波との誤差が最小となるときの仮定音速を、前記着目領域における局所音速として判定するようにしている。これは、ホイヘンスの原理により、上格子点からの受信波と、下格子点からの受信波とが一致することを利用している。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する上格子点と下格子点の反射の受信信号を取得し、前記局所音速算出手段は、前記取得した受信信号に基づいて前記着目領域における局所音速を算出する。
本発明の更に他の態様によれば、屈折のない走査線上の格子点のみの反射の受信信号を取得することにより、前記局所音速の算出時に誤差のない局所音速を算出することができる。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、前記下格子点での反射の受信信号に基づいて、該下格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの第1の伝播時間を算出する手段と、スネルの法則により前記上格子点から下格子点に入射する超音波の入射角と、前記着目領域の仮定音速と前記下格子点と前記超音波探触子との間の領域の環境音速とに基づいて前記下格子点から出射する超音波の出射角を算出する手段と、前記下格子点から前記算出した出射角で出射する超音波が入射する前記超音波探触子の素子の位置と該素子に入射するまでの第2の伝播時間とを算出する手段と、前記超音波探触子の素子の位置における超音波の受信時刻を、前記第1の伝播時間と第2の伝播時間とを加算して算出する手段と、前記上格子点での反射の前記超音波探触子の素子の位置における受信時刻と前記算出した受信時刻との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、を更に備えている。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、前記上格子点及び下格子点での反射の受信信号に基づいて、各格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、前記上格子点を反射点としたときの第1の受信波を、該上格子点に対応して算出した環境音速に基づいて算出する第1の算出手段と、前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの伝播時間を算出する手段と、前記下格子点からの第2の受信波を、該下格子点に対応して算出した環境音速及び前記算出した伝播時間に基づいて算出する第2の算出手段と、前記第1の算出手段により算出された第1の受信波と前記第2の算出手段により算出された第2の受信波との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、を更に備えている。
本発明の更に他の態様に係る超音波診断装置において、前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、前記判定手段により判定された屈折のない走査線上の上格子点と下格子点との間の超音波の伝播時間を、前記受信信号取得手段により取得した受信信号に基づいて算出する伝播時間算出手段と、前記判定手段により判定された屈折のない走査線上の上格子点及び下格子点の位置と、前記伝播時間算出手段により算出した前記上格子点と下格子点との間の超音波の伝播時間とに基づいて前記上格子点と下格子点との間の着目領域における局所音速を算出する局所音速算出手段と、を更に備えている。
本発明の更に他の態様によれば、屈折のない走査線上の格子点のみの反射の受信信号を取得し、これらの格子点間の位置及び伝播時間を求めることにより前記局所音速を精度よく算出できるようにしている。
本発明の更に他の態様に係る超音波診断装置において、前記伝播時間算出手段は、着目する走査線上の上格子点から前記超音波探触子の各素子までの第1の伝播時間を算出する第1の伝播時間算出手段と、前記上格子点と前記超音波探触子との間の領域に設定された下格子点から前記超音波探触子の各素子までの第2の伝播時間を算出する第2の伝播時間算出手段と、前記超音波探触子の各素子上で、前記算出された第1の伝播時間と第2の伝播時間との時間差が最大となる素子上の時間差を、前記上格子点から下格子点までの超音波の伝播時間として算出する手段と、からなる。
本発明の更に他の態様によれば、前記上格子点から前記超音波探触子の各素子までの第1の伝播時間と前記下格子点から前記超音波探触子の各素子までの第2の伝播時間とに基づいて、前記上格子点から下格子点までの超音波の伝播時間を簡単に算出することができる。
本発明の更に他の態様に係る超音波診断装置において、前記超音波探触子から送受信される走査線のステア角を調整するステア角調整手段を備え、前記判定手段は、前記ステア角調整手段によりステア角が調整される毎に前記取得した受信信号に基づいて走査線が屈折しているか否かを判定することを特徴としている。走査線が異なる媒質の境界面に対してほぼ直交するように入射する場合には、その走査線の屈折は小さくなる。従って、走査線のステア角を調整することにより、走査線の屈折が小さくなる走査線を送信することができる。
本発明の更に他の態様に係る超音波診断装置において、前記判定手段により判定された屈折のない走査線を表示する表示手段を更に備えている。
本発明の更に他の態様に係る超音波診断方法は、複数の素子を含む超音波探触子から超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して超音波検出信号を取得する工程と、前記取得した超音波検出信号に基づいて複数の反射点からの受信信号を取得する受信信号取得工程と、前記取得した受信信号に基づいて走査線が屈折しているか否かを判定する判定工程と、を含む。
本発明によれば、通常の装置構成により取得可能な超音波探触子の各素子の受信信号に基づいて屈折のない走査線を簡単に判定することができ、また、屈折のない走査線のみを使用することにより、着目領域における音速(局所音速)を精度よく算出することができるという効果、また横波音速計測や横方向のスペックルトラッキングなどの機能において方位方向の変位検出精度を改良できるという効果がある。
以下、添付図面に従って本発明に係る超音波診断装置及び方法の好ましい実施の形態について説明する。
[装置構成]
図1は本発明に係る超音波診断装置の実施形態を示すブロック図である。
図1は本発明に係る超音波診断装置の実施形態を示すブロック図である。
図1に示す超音波診断装置10は、超音波探触子300から被検体OBJに超音波ビームを送信して、被検体OBJによって反射された超音波ビーム(超音波エコー)を受信し、超音波エコーの検出信号から超音波画像を作成・表示する装置である。
CPU(Central Processing Unit:中央処理装置)100は、操作入力部200からの操作入力に応じて超音波診断装置10の各ブロックの制御を行う。
操作入力部200は、オペレータからの操作入力を受け付ける入力デバイスであり、操作卓202とポインティングデバイス204とを含んでいる。操作卓202は、文字情報(例えば、患者情報)の入力を受け付けるキーボードと、振幅画像(Bモード画像)を単独で表示するモードと局所音速の判定結果を表示するモードとの間で表示モードを切り替える表示モード切り替えボタンと、ライブモードとフリーズモードとの切り替えを指示するためのフリーズボタンと、シネメモリ再生を指示するためのシネメモリ再生ボタンと、超音波画像の解析・計測を指示するための解析・計測ボタンとを含んでいる。ポインティングデバイス204は、表示部104の画面上における領域の指定の入力を受け付けるデバイスであり、例えば、トラックボール又はマウスである。尚、ポインティングデバイス204としては、タッチパネルを用いることも可能である。
格納部102は、CPU100により超音波診断装置10の各ブロックの制御を制御するための制御プログラム、パラメータ及び本発明に係る屈折のない走査線を抽出し、局所音速を算出するためのプログラムを格納する記憶装置であり、例えば、ハードディスク又は半導体メモリである。
表示部104は、例えば、CRT(Cathode Ray Tube)ディスプレイ又は液晶ディスプレイであり、超音波画像(動画及び静止画)の表示、本発明に係る走査線方向、局所音速マップ、及び各種の設定画面を表示する。
超音波探触子300は、被検体OBJに当接させて用いるプローブであり、1次元又は2次元の超音波トランスデューサアレイを構成する複数の素子302を備えている。複数の素子302は、送信回路402から印加される駆動信号に基づいて超音波ビームを被検体OBJに送信するとともに、被検体OBJから反射される超音波エコーを受信して検出信号を出力する。
超音波探触子300の各素子302は、圧電性を有する材料(圧電体)の両端に電極が形成されて構成された振動子を含んでいる。上記振動子を構成する圧電体としては、例えば、PZT(チタン酸ジルコン酸鉛:Pb (lead) zirconate titanate)のような圧電セラミック、PVDF(ポリフッ化ビニリデン:polyvinylidene difluoride)のような高分子圧電素子を用いることができる。上記振動子の電極に電気信号を送って電圧を印加すると圧電体が伸縮し、この圧電体の伸縮により各振動子において超音波が発生する。例えば、振動子の電極にパルス状の電気信号を送るとパルス状の超音波が発生し、振動子の電極に連続波の電気信号を送ると連続波の超音波が発生する。そして、各振動子において発生した超音波が合成されて超音波ビームが形成される。また、各振動子により超音波が受信されると、各振動子の圧電体が伸縮して電気信号を発生する。各振動子において発生した電気信号は、超音波の検出信号として受信回路404に出力される。
尚、超音波探触子300の素子302としては、超音波変換方式の異なる複数種類の素子を用いることも可能である。例えば、超音波を送信する素子として上記圧電体により構成される振動子を用いて、超音波を受信する素子として光検出方式の超音波トランスデューサを用いるようにしてもよい。ここで、光検出方式の超音波トランスデューサとは、超音波信号を光信号に変換して検出するものであり、例えば、ファブリーペロー共振器又はファイバブラッググレーティングである。
次に、ライブモード時における超音波診断処理について説明する。ライブモードは、被検体OBJに超音波探触子300を当接させて超音波の送受信を行うことによって得られた超音波画像(動画)の表示、解析・計測を行うモードである。
超音波探触子300が被検体OBJに当接されて、操作入力部200からの指示入力により超音波診断が開始されると、CPU100は、送受信部400に制御信号を出力して、超音波ビームの被検体OBJへの送信、及び被検体OBJからの超音波エコーの受信を開始させる。CPU100は、素子302毎に超音波ビームの送信方向と超音波エコーの受信方向とを設定する。
また、CPU100は、超音波ビームの送信方向に応じて送信遅延パターンを選択するとともに、超音波エコーの受信方向に応じて受信遅延パターンを選択する。ここで、送信遅延パターンとは、複数の素子302から送信される超音波によって所望の方向に超音波ビームを形成するために駆動信号に与えられる遅延時間のパターンデータであり、受信遅延パターンとは、複数の素子302によって受信される超音波によって所望の方向からの超音波エコーを抽出するために検出信号に与えられる遅延時間のパターンデータである。上記送信遅延パターン及び受信遅延パターンは予め格納部102に格納されている。CPU100は、格納部102に格納されているものの中から送信遅延パターン及び受信遅延パターンを選択し、選択した送信遅延パターン及び受信遅延パターンに従って、送受信部400に制御信号を出力して超音波の送受信制御を行う。
送信回路402は、CPU100からの制御信号に応じて駆動信号を生成して、該駆動信号を素子302に印加する。また、送信回路402は、図2に示すように素子302毎に遅延回路τ1~τNを有し、CPU100によって選択された送信遅延パターンに基づいて、各素子302に印加する駆動信号を遅延させる。ここで、送信回路402は、複数の素子302から送信される超音波が超音波ビームを形成するように、各素子302に駆動信号を印加するタイミングを調整(遅延)したり、図2に示すように超音波ビームの方向(ステア角α)を調整するように、各素子302に駆動信号を印加するタイミングを調整(遅延)する。尚、複数の素子302から一度に送信される超音波が被検体OBJの撮像領域全体に届くように、駆動信号を印加するタイミングを調節するようにしてもよい。
受信回路404は、超音波探触子300の各素子302から出力される超音波検出信号を受信して増幅する。上記のように、各素子302と被検体OBJ内の超音波反射源との間の距離がそれぞれ異なるため、各素子302に反射波が到達する時間が異なる。受信回路404は遅延回路を備えており、CPU100によって選択された音速(以下、「仮定音速」という)又は音速の分布に基いて設定される受信遅延パターンに従って、反射波の到達時刻の差(遅延時間)に相当する分、各検出信号を遅延させる。
次に、受信回路404は、遅延時間を与えた検出信号を整合加算することにより受信フォーカス処理を行う。超音波反射源XROIと異なる位置に別の超音波反射源がある場合には、別の超音波反射源からの超音波検出信号は到達時刻が異なるので、上記加算回路で加算することにより、別の超音波反射源からの超音波検出信号の位相が打ち消し合う。これにより、超音波反射源XROIからの受信信号が最も大きくなり、フォーカスが合う。上記受信フォーカス処理によって、超音波エコーの焦点が絞り込まれた音線信号(以下、「RF信号」という)が形成される。
A/D変換器406は、受信回路404から出力されるアナログのRF信号をデジタルRF信号(以下、「RFデータ」という)に変換する。ここで、RFデータは、受信波(搬送波)の位相情報を含んでいる。A/D変換器406から出力されるRFデータは、信号処理部502とシネメモリ602にそれぞれ入力される。
シネメモリ602は、A/D変換器406から入力されるRFデータを順次格納する。また、シネメモリ602は、CPU100から入力されるフレームレートに関する情報(例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータ)を上記RFデータに関連付けて格納する。
信号処理部502は、上記RFデータに対して、STC(Sensitivity Time gain Control)によって、超音波の反射位置の深度に応じて距離による減衰の補正をした後、包絡線検波処理を施し、Bモード画像データ(超音波エコーの振幅を点の明るさ(輝度)により表した画像データ)を生成する。
信号処理部502によって生成されたBモード画像データは、通常のテレビジョン信号の走査方式と異なる走査方式によって得られたものである。このため、DSC(Digital Scan Converter)504は、上記Bモード画像データを通常の画像データ(例えば、テレビジョン信号の走査方式(NTSC方式)の画像データ)に変換(ラスター変換)する。画像処理部506は、DSC504から入力される画像データに、各種の必要な画像処理(例えば、階調処理)を施す。
画像メモリ508は、画像処理部506から入力される画像データを格納する。D/A変換器510は、画像メモリ508から読み出された画像データをアナログの画像信号に変換して表示部104に出力する。これにより、超音波探触子300によって撮影された超音波画像(動画)が表示部104に表示される。
尚、本実施形態では、受信回路404において受信フォーカス処理が施された検出信号をRF信号としたが、受信フォーカス処理が施されていない検出信号をRF信号としてもよい。この場合、複数の素子302から出力される複数の超音波検出信号が、受信回路404において増幅され、増幅された検出信号、即ち、RF信号が、A/D変換器406においてA/D変換されることによってRFデータが生成される。そして、上記RFデータは、信号処理部502に供給されるとともに、シネメモリ602に格納される。受信フォーカス処理は、信号処理部502においてデジタル的に行われる。
次に、シネメモリ再生モードについて説明する。シネメモリ再生モードは、シネメモリ602に格納されているRFデータに基づいて超音波診断画像の表示、解析・計測を行うモードである。
操作卓202のシネメモリ再生ボタンが押下されると、CPU100は、超音波診断装置10の動作モードをシネメモリ再生モードに切り替える。シネメモリ再生モード時には、CPU100は、オペレータからの操作入力により指定されたRFデータの再生をシネメモリ再生部604に指令する。シネメモリ再生部604は、CPU100からの指令に従って、シネメモリ602からRFデータを読み出して、画像信号生成部500の信号処理部502に送信する。シネメモリ602から送信されたRFデータは、信号処理部502、DSC504及び画像処理部506において所定の処理(ライブモード時と同様の処理)が施されて画像データに変換された後、画像メモリ508及びD/A変換器510を経て表示部104に出力される。これにより、シネメモリ602に格納されたRFデータに基づく超音波画像(動画又は静止画)が表示部104に表示される。
ライブモード又はシネメモリ再生モード時において、超音波画像(動画)が表示されているときに操作卓202のフリーズボタンが押下されると、フリーズボタン押下時に表示されている超音波画像が表示部104に静止画表示される。これにより、オペレータは、着目領域(ROI:Region of Interest)の静止画を表示させて観察することができる。
操作卓202の計測ボタンが押下されると、オペレータからの操作入力により指定された解析・計測が行われる。データ解析計測部106は、各動作モード時に計測ボタンが押下された場合に、A/D変換器406又はシネメモリ602から、画像処理が施される前のRFデータを取得し、当該RFデータを用いてオペレータ指定の解析・計測(例えば、組織部の歪み解析(硬さ診断)、血流の計測、組織部の動き計測、又はIMT(内膜中膜複合体厚:Intima-Media Thickness)値計測)を行う。データ解析計測部106による解析・計測結果は、画像信号生成部500のDSC504に出力される。DSC504は、データ解析計測部106による解析・計測結果を超音波画像の画像データに挿入して表示部104に出力する。これにより、超音波画像と解析・計測結果とが表示部104に表示される。
また、表示モード切り替えボタンが押下されると、Bモード画像を単独で表示するモード、Bモード画像に局所音速の判定結果を重畳して表示するモード(例えば、局所音速に応じて色分け又は輝度を変化させる表示、又は局所音速が等しい点を線で結ぶ表示)、Bモード画像と局所音速値の判定結果の画像を並べて表示するモードの間で表示モードが切り替わる。これにより、オペレータは、局所音速の判定結果を観察することで、例えば、病変を発見することができる。尚、局所音速の判定結果に基づいて、送信フォーカス処理及び受信フォーカス処理の少なくとも一方を施すことにより得られたBモード画像を表示部104に表示してもよい。
また、超音波診断装置10は、上記局所音速を精度よく算出するための前提として、超音波探触子300の各素子302での受信信号に基づいて、超音波ビーム(以下、「走査線」という)のうちの屈折のない走査線を抽出する。尚、屈折のない走査線の抽出方法については、後述する。また、表示部104には、屈折のない走査線を表示することができる。
[局所音速測定の第1の実施形態]
図3は、被検体の着目領域における局所音速を算出する処理手順の第1の実施形態を示すフローチャートである。
図3は、被検体の着目領域における局所音速を算出する処理手順の第1の実施形態を示すフローチャートである。
図3に示すように、まず被検体の着目領域を設定する(ステップS1)。この着目領域は、表示部104に表示される超音波画像の静止画上で、オペレータがポインティングデバイスにより設定してもよいし、制御プログラムが自動的に所定位置、所定サイズにて設定してもよいし、超音波画像を二値化処理するとともに、白の部分(又は黒の部分)が連続した画素に同じ番号を割り振るラベリング処理を行い、ラベリングした番号順に自動的に設定してもよい。
続いて、前記設定した着目領域内に格子点(上格子点及び下格子点を含む)を設定し、各格子点における環境音速を求める(ステップS2)。各格子点は、走査線位置と受信時刻によって、その位置が定義される。即ち、図4に示すように、超音波探触子300から被検体OBJの着目領域に出射される、所定の間隔の走査線1,2,…,n上の深さの異なる反射点を格子点として設定する。ここで、下格子点A1,A2,A3,…,Anは、各走査線1,2,…,n上の受信時刻が同一の反射点であり、同様に上格子点B1,B2,B3,…,Bn、上格子点C1,C2,C3,…,Cn、…もそれぞれ各走査線1,2,…,n上の受信時刻が同一の反射点である。
尚、図4上では、下格子点A、上格子点B,Cは、同じ深さの格子点として図示されているが、実際には各格子点と超音波探触子300との間の領域の音速は均一でないため、空間上では異なる深さの反射点となり、また、リニア走査される各走査線1,2,…,nも走査線の伝播領域の音速の違いにより屈折するため、すべての走査線は必ずしも平行にはならない。
また、各格子点の範囲及び個数は予め決めておく。ここで、局所音速演算に使用する格子点の範囲が広いと局所音速値の誤差が大きくなり、狭いと仮想受信波との誤差が大きくなるため、格子点の範囲はこれらの兼ね合いで決める。各格子点のx方向の間隔は、分解能と処理時間の兼ね合いで決定される。格子点のx方向の間隔は、一例で1mmから1cmである。また、下格子点と上格子点のy方向の間隔が狭いと誤差計算における誤差が大きくなり、広いと局所音速の誤差が大きくなる。格子点のy方向の間隔は、超音波画像の画像分解能の設定に基づいて決定され、一例で1cmである。
次に、上記のように設定した各格子点における環境音速を以下のようにして算出する。
<環境音速の算出>
図12(a)に示したように、ある反射点(格子点)X1ROIから超音波探触子300Aまでの距離をLとすると、格子点X1ROIで超音波が反射されてから格子点X1ROIの直下の素子302A0で受信されるまでの経過時間Tは、T=L/Vである。素子302A0からX方向(素子302Aの配列方向)に距離X離れた位置にある素子302Aiで受信されるまで経過時間をT+ΔTとすると、素子302A0と302Aiとの間の遅延時間ΔTは、下記の式(1)により表される。
図12(a)に示したように、ある反射点(格子点)X1ROIから超音波探触子300Aまでの距離をLとすると、格子点X1ROIで超音波が反射されてから格子点X1ROIの直下の素子302A0で受信されるまでの経過時間Tは、T=L/Vである。素子302A0からX方向(素子302Aの配列方向)に距離X離れた位置にある素子302Aiで受信されるまで経過時間をT+ΔTとすると、素子302A0と302Aiとの間の遅延時間ΔTは、下記の式(1)により表される。
従って、超音波が送波されて格子点X1ROIで時間T後に反射された後、各素子により受信されるまでの経過時間[2T,2T+ΔT]を測定することにより、格子点X1ROIまでの距離Lと速度Vを一意に求めることができる。
ここで、ある格子点の環境音速とは、その格子点から超音波探触子までの領域の平均音速であり、画像のコントラスト、シャープネスが最も高くなる音速である。従って、環境音速の判定方法としては、例えば、画像のコントラスト、スキャン方向の空間周波数、分散などから判定する方法(例えば、特開平8-317926号公報)を適用することができる。また、環境音速を十分な精度で求められるように、着目領域内にて細密深度間隔にて送信焦点を形成するようにフォーカスをかけることが好ましい。
この様にして求めた環境音速から、その格子点からの反射の各素子の受信時刻を得る事ができる。つまり、その格子点からの反射の各素子受信信号に対して、ある音速を仮定して遅延を決定し、その遅延を使用して生成した画像のコントラスト、シャープネスが最も高くなる場合、その遅延が各素子受信時刻に最も近づく事を意味しているため、その音速(環境音速)、つまり遅延を以って各素子受信時刻と見做す事ができる。環境音速の代わりに、位相収差解析などの手法により、各素子の受信時刻を求めて、以降に使用してもよい。
<局所音速の算出>
次に、各下格子点A1,A2,A3,…,Anを通る走査線1,2,…,nが屈折しておらず、各走査線の走査方向の位置が既知であるという前提の下で、局所音速を求める(ステップS3)。
次に、各下格子点A1,A2,A3,…,Anを通る走査線1,2,…,nが屈折しておらず、各走査線の走査方向の位置が既知であるという前提の下で、局所音速を求める(ステップS3)。
<局所音速の第1の算出方法>
まず、着目する走査線上の上格子点から、周囲走査線上の各下格子点への伝播時間を求め、その後、各走査線が屈折していないという前提で、着目走査線上の上格子点に対して周囲走査線上の各下格子点の相対的な空間座標を与える事で、上格子点から各下格子点への伝播時間から局所音速を求める。
まず、着目する走査線上の上格子点から、周囲走査線上の各下格子点への伝播時間を求め、その後、各走査線が屈折していないという前提で、着目走査線上の上格子点に対して周囲走査線上の各下格子点の相対的な空間座標を与える事で、上格子点から各下格子点への伝播時間から局所音速を求める。
図5に示すように、着目する走査線上の格子点(上格子点)Bから周囲走査線上の格子点(下格子点)A1,A2,A3への伝播時間ΔT1,ΔT2,ΔT3は、次式により算出することができる。
[数2]
ΔTn=max(TBi-TAni)
上記[数2]式において、
ΔTn:格子点Bから格子Anまでの伝播時間
TBi:格子点Bから超音波探触子の素子iまでの伝播時間(格子点B反射の素子iでの受信時刻)
TAni:格子点Anから超音波探触子の素子iまでの伝播時間(格子点An反射の素子iでの受信時刻)
である。
ΔTn=max(TBi-TAni)
上記[数2]式において、
ΔTn:格子点Bから格子Anまでの伝播時間
TBi:格子点Bから超音波探触子の素子iまでの伝播時間(格子点B反射の素子iでの受信時刻)
TAni:格子点Anから超音波探触子の素子iまでの伝播時間(格子点An反射の素子iでの受信時刻)
である。
但し、TBi及びTAniは、伝播経路の片道の伝播時間を示しており、超音波探触子の各素子の受信時刻又は環境音速から往路分の伝播時間(着目走査線上の素子における受信時刻または最短受信時刻の半分)を引くことにより求められる。
図5には、格子点Bの反射の超音波探触子の各素子iでの受信時刻を示す曲線と、格子点A1,A2,A3の反射の超音波探触子の各素子iでの受信時刻を示す曲線とが示されている。ホイヘンスの原理によれば、格子点Bからの受信波(受信時刻を示す曲線)と、格子点A1,A2,A3…からの受信波を、格子点Bからそれぞれの格子点までの伝播時間だけ遅延させて仮想的に合成した仮想合成受信波(各受信時刻を示す曲線の包絡線)とは一致する。上記[数2]式により算出されるΔTnは、格子点Bからの受信波と、格子点A1,A2,A3…からの仮想合成受信波とが一致するために必要な格子点Bから格子点Anまでの伝播時間を示す。又は、[数2]式の様に各伝播時間ΔTnを独立に求める代わりに、下格子点A1,A2,A3・・・の位置を仮想的に既知として、特許文献1に開示されている方法によって上格子点Bと下格子点A1,A2,A3・・・との間の局所音速を仮想的に求めても良い。仮想的に決めた下格子点A1,A2,A3・・・の位置および求めた局所音速に基き、各伝播時間ΔTnを求める事ができる。
ここで、着目走査線および周囲走査線が屈折していないという前提の下、上格子点Bおよび各下格子点A1,A2,A3・・・の走査方向の位置が与えられている。そして着目走査線において上格子点から下格子点までの伝播時間が与えられているため、上格子点と下格子点との間の局所音速を仮定する事によって上格子点と下格子点の深さ方向の距離が与えられ、従って上格子点Bから各下格子点A1,A2,A3・・・までの伝播時間を求める事ができる。この各伝播時間と上記で求めた各伝播時間ΔTnを比較し、誤差が最小となるときの仮定した局所音速を真の音速(局所音速)として判定する。
走査線毎に上記の処理を繰り返すことにより、走査線毎に局所音速を算出することができる。また、上格子点Bの替わりに、深さ方向の異なる上格子点C、D、…からの受信信号を使用し(図4参照)、上記と同様にして深さ方向の異なる格子点における局所音速を算出する。
<局所音速の第2の算出方法>
特許文献1に開示されている手法において、各走査線が屈折していないという前提で下格子点に空間座標を与えた上で、着目走査線上の各上格子点について演算を繰り返して局所音速を求める。
特許文献1に開示されている手法において、各走査線が屈折していないという前提で下格子点に空間座標を与えた上で、着目走査線上の各上格子点について演算を繰り返して局所音速を求める。
図6は、特許文献1に開示された屈折モデル計算により局所音速を算出する方法を模式的に示した図である。
以下の説明では、超音波探触子300の各素子302が配置された素子面S2に平行な方向をX方向とし、X方向に垂直な方向(被検体OBJの深さ方向)をY方向とする。
図6に示すように、被検体OBJJ内の領域A内の着目領域ROIを代表する上格子点をBROIとし、下格子点をA1,A2,…,An,…とする。これらの格子点の空間座標は、各走査線が屈折していないという前提で与えられている。
下格子点A1,A2,…,An,…を連結した境界面S1と被検体OBJ内の上格子点BROIとの間の領域を領域Aとし、境界面S1と超音波探触子300の素子面S2との間の領域を領域Bとする。領域Aと領域Bの中における音速はそれぞれ一定と仮定する。
下格子点A1,A2,…から超音波探触子300の素子面S2に至る領域の音速(環境音速)が略同じ場合、又は格子点A1,A2,…からの受信波が互いに同じ場合、又は該受信波が近似的に同じと見なせる場合、又は該受信波がゆるやかに変化する場合に、図6に示すように、上格子点BROIと下格子点A1,A2,…との間の領域Aにおける音速(局所音速)と領域Bにおける環境音速に基づいて、スネルの法則に従って領域AとBとの境界面で屈折する音線を追跡することにより各素子302における受信時刻を求める。
具体的には、着目領域の仮定音速をVAとし、上格子点BROIから或る下格子点X’に入射する音線の入射角をΘとすると、下格子点X’を通る音線の出射角(屈折角)Θ’は、スネルの法則により、次式で表すことができる。
[数3]
sinΘ/sinΘ’=VA/VB
各格子点の空間座標は既知であるため、入射角Θも既知であり、仮定音速VAを仮定することにより、下格子点X’を通る音線の屈折角Θ’を、上記[数3]式により求めることができる。
sinΘ/sinΘ’=VA/VB
各格子点の空間座標は既知であるため、入射角Θも既知であり、仮定音速VAを仮定することにより、下格子点X’を通る音線の屈折角Θ’を、上記[数3]式により求めることができる。
これにより、下格子点X’を通る音線が入射する、超音波探触子300上の素子302の位置X”と、下格子点X’から素子302の位置X”までの音線の伝播時間を算出することができる。
また、上格子点BROIから下格子点A1,A2,…までの伝播時間を算出する。この伝播時間は、各格子点間の距離を求めることができるため、仮定音速VAを仮定することにより算出することができる。
上記のように音線を追跡することにより、上格子点BROIから各下格子点A1,A2,…を通過し、超音波探触子300のどの素子の位置に、どの伝播時間(受信時刻)で受信されるかをそれぞれ算出することができる。
一方、格子点BROIでの反射の超音波探触子300の各素子302の位置における実際の受信時刻は、ステップS2で測定してあるため、前記算出した受信時刻と測定した受信時刻との誤差が最小となるときの仮定音速を、着目領域における真の音速(局所音速)として判定する。
走査線毎に上記の処理を繰り返すことにより、走査線毎に局所音速を算出することができる。また、上格子点Bの替わりに、深さ方向の異なる上格子点C、D、…からの受信信号を使用し(図4参照)、上記と同様にして深さ方向の異なる格子点における局所音速を算出する。
<局所音速の第3の算出方法>
図7は、特許文献1に開示されたホイヘンスの原理を利用して局所音速を算出する方法を模式的に示した図である。
図7は、特許文献1に開示されたホイヘンスの原理を利用して局所音速を算出する方法を模式的に示した図である。
図7(b)に示すように下格子点A1,A2,…からの受信波(それぞれWA1,WA2,…)の(伝播時間T及び遅延時間ΔT)を既知として、格子点BROIと格子点A1,A2,…の位置関係から格子点BROIにおける局所音速を求める。具体的には、ホイヘンスの原理により、上格子点BROIからの受信波WXと下格子点A1,A2,…からの受信波を仮想的に合成した受信波WSUMとが一致することを利用する。
ここで、上格子点BROI、下格子点A1,A2,…の空間座標は、各走査線が屈折していないという前提により与えられている。
図7に示すように、上格子点BROIにおける環境音速Vに基づいて上格子点BROIを反射点としたときの受信波WXの波形を算出する。また、上格子点BROIから各下格子点A1,A2,…までの伝播時間をそれぞれ算出する。これらの伝播時間は、各走査線が屈折していないという前提により、下格子点A1,A2,…の空間座標が与えられているため、仮定音速VAを仮定することにより算出することができる。
各格子点A1,A2,…における環境音速に基づいて下格子点A1,A2,…を反射点としたときの受信波WA1,WA2,…を算出する。そして、これらの受信波WA1,WA2,…を、各格子点A1,A2,…毎に算出した上格子点BROIからの伝播時間だけ遅延させて合成することにより、仮想的な合成受信波WSUMを算出する。
次に、上記受信波WXと合成受信波WSUMの誤差を算出する。受信波WXと合成受信波WSUMの誤差は、互いの相互相関をとる方法、受信波WXに合成受信波WSUMから得られる遅延を掛けて位相整合加算する方法、又は逆に合成受信波WSUMに受信波WXから得られる遅延を掛けて位相整合加算する方法により算出される。ここで、受信波WXから遅延を得るには、格子点BROIを反射点とし、音速Vで伝播した超音波が各素子に到着する時刻を遅延とすればよい。また、合成受信波WSUMから遅延を得るには、隣り合う素子間での合成受信波の位相差から等位相線を抽出し、その等位相線を遅延とするか、又は単に各素子の合成受信波の最大(ピーク)位置の位相差を遅延としてもよい。また、各素子からの合成受信波の相互相関ピーク位置を遅延としてもよい。位相整合加算時の誤差は、整合加算後の波形のpeak to peakとする方法、又は包絡線検波した後の振幅の最大値とする方法により求められる。
上記受信波WXと合成受信波WSUMの誤差は、仮定音速VAによって変化する。そして、誤差が最小(位相整合加算時には最大)となるときの仮定音速を、着目領域における真の音速(局所音速)として判定する。
走査線毎に上記の処理を繰り返すことにより、走査線毎に局所音速を算出することができる。また、上格子点Bの替わりに、深さ方向の異なる上格子点C、D、…からの受信信号を使用し(図4参照)、上記と同様にして深さ方向の異なる格子点における局所音速を算出する。
尚、第1~3の算出方法において算出した上格子点B、C、D、・・・における局所音速は、それぞれの上格子点と下格子点A1、A2、…との間の領域の局所音速に相当し、つまり下格子点を共通として各領域に重なりがあるが、各領域が重ならない様に、それぞれの上格子点に対して下格子点を別に設定しても良い。この時、それぞれの上格子点に対して基準となる深さにそれぞれの上格子点共通の格子点を別途設定し、この格子点における環境音速が一致するように、それぞれの上格子点における局所音速を算出しても良い。
具体的には、以下の手順により算出する。図8に示す様に、上格子点をB、下格子点を点A1、A2、…とした時に、仮の下格子点A’1、A’2、…を別の基準深さに設定する。次に、仮の下格子点A’1、A’2、…における環境音速を仮定する。すると上格子点B、および下格子点A1、A2、…、それぞれに対し第1の算出方法または第2の算出方法によって、下格子点A’1、A’2、…までの領域の局所音速を算出することができる。ここで、上格子点Bに対する局所音速と下格子点A1、A2、…に対する局所音速が一致する場合の局所音速および下格子点A’1、A’2、…における環境音速を真値と見做す事ができる。この様にして、それぞれの上格子点B、C、D、・・・における局所音速を算出することができる。
<着目領域の局所音速算出>
次に、各走査線について、深さの異なる格子点(領域)における局所音速の深さ方向の変化に基づいて屈折のない走査線を抽出し、着目領域の局所音速を求める(図3のステップS4)。
次に、各走査線について、深さの異なる格子点(領域)における局所音速の深さ方向の変化に基づいて屈折のない走査線を抽出し、着目領域の局所音速を求める(図3のステップS4)。
図9(a)~(c)は、生体ファントムの各深さにおいて、着目走査線および周囲走査線が屈折していない前提で局所音速を測定した結果を示すグラフであり、それぞれ測定点より浅い領域において異なる音速媒質を、その境界面の形状を変えて設置している。尚、これらのグラフにおいて、振幅が大きく振れているのはノイズによる影響である。
生体ファントムと異なる音速媒質の境界面を探触子面と平行にし各走査線を屈折させない場合、図9(a)に示すように、局所音速の深さ方向に対する傾きはゼロ近傍となる。
一方、境界面形状を変えて各走査線を深さと共に閉じてゆく様に屈折させた場合、または開いてゆく様に屈折させた場合、それぞれ図9(b)及び(c)に示すように、局所音速の深さ方向に対する傾きは単調増加又は単調減少する。
そこで、走査線毎に局所音速の深さ方向に対する傾きを求め、その傾きがゼロ近傍の走査線を抽出することにより、屈折のない走査線を抽出することができる。尚、屈折のない走査線が連続している走査線群を抽出することが好ましい。
そして、上記のように屈折がないと判定した走査線上の各局所音速の平均値を算出し、この算出した平均値を着目領域における局所音速とする。
また、屈折がないと判定した走査線上の上格子点及び下格子点からの受信信号のみを使用し、改めて特許文献1に開示された局所音速の算出方法を適用して局所音速を算出するようにしてもよい。このとき、屈折のない各走査線の環境音速を平均化し、その深さプロファイルに基づいて局所音速を求めてもよいし、受信信号又はフォーカス指標を平均化して、それに基づき環境音速を求め、更に局所音速を求めてもよい。
[局所音速測定の第2の実施形態]
図10は、被検体の着目領域における局所音速を算出する処理手順の第2の実施形態を示すフローチャートである。尚、図3に示した第1の実施形態と共通する部分には、同一のステップ番号を付し、その詳細な説明は省略する。
図10は、被検体の着目領域における局所音速を算出する処理手順の第2の実施形態を示すフローチャートである。尚、図3に示した第1の実施形態と共通する部分には、同一のステップ番号を付し、その詳細な説明は省略する。
図10に示すように、第2の実施形態は、ステップS3’、S4’の処理が、第1の実施形態のステップS3、S4の処理と相違する。
ステップS3’では、ステップS2で求めた走査線毎の環境音速に基づいて所定の走査線幅における走査方向の環境音速の変動を算出する。このとき、走査線の深さ方向に求めた複数の環境音速を平均し、その平均した環境音速の変動を算出するようにしてもよい。
図11は走査線の走査線方向に対する環境音速の変化の一例を示すグラフである。
ステップS3’の処理は、予め設定した走査線数分の連続する環境音速の標準偏差、分散、又は最大値と最小値の差分を、連続する環境音速の変動情報として算出する。このとき、着目領域の格子点における走査線方向の環境音速の変動のみでなく、浅い領域の走査線方向の環境音速の変動を含めるようにしてもよい。
ステップS4’では、ステップS3’で求めた環境音速の変動が、予め設定した閾値以下となる走査線群を抽出し、この走査線群を屈折のない走査線群として判定する。このようにして屈折のない走査線が抽出されると、屈折のない走査線上の上格子点及び下格子点からの受信信号のみを使用し、特許文献1に開示された屈折モデル計算、ホインヘンスの原理等を使用して局所音速を算出する。
<その他の実施形態>
超音波探触子300から送出される複数の走査線と着目領域との位置関係によっては、屈折のない走査線群が存在しないことが考えられる。この場合には、超音波探触子300から着目領域に向けて出射する走査線のステア角αを、図2に示すように各素子への遅延時間を調整し、着目領域への走査線の入射角を変更する。
超音波探触子300から送出される複数の走査線と着目領域との位置関係によっては、屈折のない走査線群が存在しないことが考えられる。この場合には、超音波探触子300から着目領域に向けて出射する走査線のステア角αを、図2に示すように各素子への遅延時間を調整し、着目領域への走査線の入射角を変更する。
このステア角αの調整は、送信回路402から超音波探触子300の素子iに印加する駆動信号を、次式に示す遅延時間Δτiだけ遅らせることにより行うことができる。
[数4]
Δτi=(i-1)p・sinα/V
但し、p:素子のピッチ、V:音速(例えば、皮下脂肪等における既知の音速)
また、受信時には、ステア角αの方向の各深さにおいて焦点を形成する様に受信フォーカスを実施し、各格子点の環境音速を求める。ここで、各格子点の環境音速を精度良く求めるために、送信時にもステア角αの方向の各深さに送信焦点を形成するようにフォーカスをかける事が好ましい。
Δτi=(i-1)p・sinα/V
但し、p:素子のピッチ、V:音速(例えば、皮下脂肪等における既知の音速)
また、受信時には、ステア角αの方向の各深さにおいて焦点を形成する様に受信フォーカスを実施し、各格子点の環境音速を求める。ここで、各格子点の環境音速を精度良く求めるために、送信時にもステア角αの方向の各深さに送信焦点を形成するようにフォーカスをかける事が好ましい。
この様にして、着目領域への走査線の入射角を変更する毎に、第1の実施形態又は第2の実施形態の方法を適用して、屈折のない走査線群又は屈折の小さい走査線群を求めることができる。
また、着目領域を小さくとることにより音速一様でない対象へも適用可能である。各領域において独立に局所音速を求めてもよいが、超音波探触子に近い領域(浅い領域)の結果を活用してもよい。例えば、浅い領域における走査線の屈折なしの判定結果を含めて屈折なしの走査線群を判定してもよい。
また、抽出した屈折のない走査線方向を示す各走査線を表示部104に表示させてもよい。
また、<環境音速の算出>で述べたように超音波走査線上のある反射点(格子点)からの受信信号とその格子点における環境音速とは一定の関係にあるため、本発明における「受信信号」は環境音速を含む概念である。また、屈折のない走査線とは、走査線が全く屈折しない場合に限らず、屈折の小さい走査線も含み、要求される局所音速の精度にもよるが、例えば、走査線方向の変化が1度程度以下の走査線をいう。
10…超音波診断装置、100…中央処理装置(CPU)、102…格納部、104…表示部、106…データ解析計測部、200…操作入力部、202…操作卓、204…ポインティングデバイス、300…超音波探触子、302…超音波トランスデューサ(素子)、400…送受信部、402…送信回路、404…受信回路、500…画像信号生成部、502…信号処理部、506…画像処理部、508…画像メモリ、510…D/A変換器、600…再生部
Claims (19)
- 超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して超
音波検出信号を出力する複数の素子を含む超音波探触子と、
前記超音波検出信号に基づいて複数の反射点からの受信信号を取得する受信信号取得手段と、
前記取得した受信信号に基づいて走査線が屈折しているか否かを判定する判定手段と、
を備えた超音波診断装置。 - 前記受信信号は、複数の素子が受信した信号である請求項1に記載の超音波診断装置。
- 前記受信信号取得手段は、複数の異なる深さの反射点からの受信信号を取得し、
前記判定手段は、該受信信号に基づいて該複数の反射点における各局所音速を算出する局所音速算出手段を有し、前記算出した局所音速の深さ方向の変化に基づいて走査線が屈折しているか否かを判定する請求項1又は2に記載の超音波診断装置。 - 前記判定手段は、複数の走査線毎に前記算出した局所音速の深さ方向に対する傾きを算出し、前記算出した傾きによりゼロ近傍の所定の閾値以内となる走査線を、屈折のない走査線として判定する請求項3に記載の超音波診断装置。
- 前記判定手段により判定された走査線上の反射点毎に算出された局所音速の平均値を算出する算出手段を備え、前記算出された平均値を前記屈折のない走査線上の反射点を含む着目領域の局所音速とする請求項3又は4に記載の超音波診断装置。
- 前記受信信号取得手段は、複数の異なる深さの反射点からの受信信号を取得し、
前記判定手段は、該受信信号に基づいて該複数の反射点における各環境音速を算出する環境音速算出手段を有し、前記算出した環境音速の深さ方向の変化に基づいて走査線が屈折しているか否かを判定する請求項1又は2に記載の超音波診断装置。 - 前記受信信号取得手段は、複数の異なる走査線位置からの受信信号を取得し、
前記判定手段は、該受信信号に基づいて該複数の走査線位置における各環境音速を算出する環境音速算出手段を有し、前記算出した環境音速の走査方向の変化に基づいて走査線が屈折しているか否かを判定する請求項1又は2に記載の超音波診断装置。 - 前記判定手段は、前記算出した環境音速の前記走査線の走査方向の変動が、予め設定された閾値以内となる連続する走査線群を、屈折のない走査線群として判定する請求項7に記載の超音波診断装置。
- 前記判定手段は、前記複数の走査線のうちの連続する所定数の走査線群毎に、前記環境音速算出手段により算出した環境音速の標準偏差、分散又は最大値と最小値の差分を算出し、その算出結果に基づいて屈折のない走査線群を判定する請求項7に記載の超音波診断装置。
- 前記受信信号取得手段は、複数の走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、
前記局所音速算出手段は、
前記下格子点での反射の受信信号に基づいて、該下格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、
前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの第1の伝播時間を算出する手段と、
スネルの法則により前記上格子点から下格子点に入射する超音波の入射角と、前記着目領域の仮定音速と前記下格子点と前記超音波探触子との間の領域の環境音速とに基づいて前記下格子点から出射する超音波の出射角を算出する手段と、
前記下格子点から前記算出した出射角で出射する超音波が入射する前記超音波探触子の素子の位置と該素子に入射するまでの第2の伝播時間とを算出する手段と、
前記超音波探触子の素子の位置における超音波の受信時刻を、前記第1の伝播時間と第2の伝播時間とを加算して算出する手段と、
前記上格子点での反射の前記超音波探触子の素子の位置における受信時刻と前記算出した受信時刻との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、を含む請求項3から5のいずれか1項に記載の超音波診断装置。 - 前記受信信号取得手段は、複数の走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、
前記局所音速算出手段は、
前記上格子点及び下格子点での反射の受信信号に基づいて、各格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、
前記上格子点を反射点としたときの第1の受信波を、該上格子点に対応して算出した環境音速に基づいて算出する第1の算出手段と、
前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの伝播時間を算出する手段と、
前記下格子点からの第2の受信波を、該下格子点に対応して算出した環境音速及び前記算出した伝播時間に基づいて算出する第2の算出手段と、
前記第1の算出手段により算出された第1の受信波と前記第2の算出手段により算出された第2の受信波との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、を含む請求項3から5のいずれか1項に記載の超音波診断装置。 - 前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する上格子点と下格子点の反射の受信信号を取得し、
前記局所音速算出手段は、前記取得した受信信号に基づいて前記着目領域における局所音速を算出する請求項10又は11に記載の超音波診断装置。 - 前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、
前記下格子点での反射の受信信号に基づいて、該下格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、
前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの第1の伝播時間を算出する手段と、
スネルの法則により前記上格子点から下格子点に入射する超音波の入射角と、前記着目領域の仮定音速と前記下格子点と前記超音波探触子との間の領域の環境音速とに基づいて前記下格子点から出射する超音波の出射角を算出する手段と、
前記下格子点から前記算出した出射角で出射する超音波が入射する前記超音波探触子の素子の位置と該素子に入射するまでの第2の伝播時間とを算出する手段と、
前記超音波探触子の素子の位置における超音波の受信時刻を、前記第1の伝播時間と第2の伝播時間とを加算して算出する手段と、
前記上格子点での反射の前記超音波探触子の素子の位置における受信時刻と前記算出した受信時刻との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、
を更に備えた請求項6から9のいずれか1項に記載の超音波診断装置。 - 前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、
前記上格子点及び下格子点での反射の受信信号に基づいて、各格子点から前記超音波探触子までの領域の平均音速である環境音速を算出する環境音速算出手段と、
前記上格子点を反射点としたときの第1の受信波を、該上格子点に対応して算出した環境音速に基づいて算出する第1の算出手段と、
前記着目領域における仮定音速を仮定し、前記上格子点から下格子点までの伝播時間を算出する手段と、
前記下格子点からの第2の受信波を、該下格子点に対応して算出した環境音速及び前記算出した伝播時間に基づいて算出する第2の算出手段と、
前記第1の算出手段により算出された第1の受信波と前記第2の算出手段により算出された第2の受信波との誤差が最小となる前記仮定音速を、前記着目領域における局所音速として判定する局所音速判定手段と、
を更に備えた請求項6から9のいずれか1項に記載の超音波診断装置。 - 前記受信信号取得手段は、前記判定手段により屈折のない走査線が判定されると、該走査線上の反射点に対応する格子点であって、所望の着目領域に設定される上格子点と、前記上格子点と前記超音波探触子との間に設定される下格子点の反射の受信信号を取得し、
前記判定手段により判定された屈折のない走査線上の上格子点と下格子点との間の超音波の伝播時間を、前記受信信号取得手段により取得した受信信号に基づいて算出する伝播時間算出手段と、
前記判定手段により判定された屈折のない走査線上の上格子点及び下格子点の位置と、前記伝播時間算出手段により算出した前記上格子点と下格子点との間の超音波の伝播時間とに基づいて前記上格子点と下格子点との間の着目領域における局所音速を算出する局所音速算出手段と、
を更に備えた請求項3から9のいずれか1項に記載の超音波診断装置。 - 前記伝播時間算出手段は、
着目する走査線上の上格子点から前記超音波探触子の各素子までの第1の伝播時間を算出する第1の伝播時間算出手段と、
前記上格子点と前記超音波探触子との間の領域に設定された下格子点から前記超音波探触子の各素子までの第2の伝播時間を算出する第2の伝播時間算出手段と、
前記超音波探触子の各素子上で、前記算出された第1の伝播時間と第2の伝播時間との時間差が最大となる素子上の時間差を、前記上格子点から下格子点までの超音波の伝播時間として算出する手段と、
からなる請求項15に記載の超音波診断装置。 - 前記超音波探触子から送受信される走査線のステア角を調整するステア角調整手段を備え、
前記判定手段は、前記ステア角調整手段によりステア角が調整される毎に前記取得した受信信号に基づいて走査線が屈折しているか否かを判定する請求項1から16のいずれか1項に記載の超音波診断装置。 - 前記判定手段により判定された屈折のない走査線を表示する表示手段を更に備えた請求項1から17のいずれか1項に記載の超音波診断装置。
- 複数の素子を含む超音波探触子から超音波を被検体に送信するとともに、該被検体によって反射される超音波を受信して超音波検出信号を取得する工程と、
前記取得した超音波検出信号に基づいて複数の反射点からの受信信号を取得する受信信号取得工程と、
前記取得した受信信号に基づいて走査線が屈折しているか否かを判定する判定工程と、
を含む超音波診断方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201280055913.2A CN103930040A (zh) | 2011-11-14 | 2012-11-13 | 超声波诊断装置及方法 |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011248776A JP2013102960A (ja) | 2011-11-14 | 2011-11-14 | 超音波診断装置及び方法 |
JP2011-248776 | 2011-11-14 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2013073514A1 true WO2013073514A1 (ja) | 2013-05-23 |
Family
ID=48429571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2012/079336 WO2013073514A1 (ja) | 2011-11-14 | 2012-11-13 | 超音波診断装置及び方法 |
Country Status (3)
Country | Link |
---|---|
JP (1) | JP2013102960A (ja) |
CN (1) | CN103930040A (ja) |
WO (1) | WO2013073514A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106999152A (zh) * | 2014-11-26 | 2017-08-01 | Ge医疗系统环球技术有限公司 | 超声诊断设备和程序 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6138596B2 (ja) | 2013-06-05 | 2017-05-31 | 富士フイルム株式会社 | 超音波診断装置、超音波診断装置の作動方法およびプログラム |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0595946A (ja) | 1991-10-09 | 1993-04-20 | Yokogawa Medical Syst Ltd | 超音波の音速測定法及び音速測定手段を備えた超音波診断装置 |
JPH08317926A (ja) | 1995-05-26 | 1996-12-03 | Hitachi Medical Corp | 超音波断層装置 |
JPH11235341A (ja) * | 1998-02-23 | 1999-08-31 | Toshiba Corp | 超音波診断装置 |
JP2009056140A (ja) | 2007-08-31 | 2009-03-19 | Fujifilm Corp | 超音波診断装置 |
JP2009061086A (ja) * | 2007-09-06 | 2009-03-26 | Fujifilm Corp | 超音波診断装置、並びに、画像処理方法及びプログラム |
JP2010099452A (ja) | 2008-09-25 | 2010-05-06 | Fujifilm Corp | 超音波診断装置及び超音波診断方法 |
JP2011172611A (ja) * | 2010-02-23 | 2011-09-08 | Canon Inc | 超音波画像装置および遅延制御方法 |
-
2011
- 2011-11-14 JP JP2011248776A patent/JP2013102960A/ja not_active Abandoned
-
2012
- 2012-11-13 CN CN201280055913.2A patent/CN103930040A/zh not_active Withdrawn
- 2012-11-13 WO PCT/JP2012/079336 patent/WO2013073514A1/ja active Application Filing
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0595946A (ja) | 1991-10-09 | 1993-04-20 | Yokogawa Medical Syst Ltd | 超音波の音速測定法及び音速測定手段を備えた超音波診断装置 |
JPH08317926A (ja) | 1995-05-26 | 1996-12-03 | Hitachi Medical Corp | 超音波断層装置 |
JPH11235341A (ja) * | 1998-02-23 | 1999-08-31 | Toshiba Corp | 超音波診断装置 |
JP2009056140A (ja) | 2007-08-31 | 2009-03-19 | Fujifilm Corp | 超音波診断装置 |
JP2009061086A (ja) * | 2007-09-06 | 2009-03-26 | Fujifilm Corp | 超音波診断装置、並びに、画像処理方法及びプログラム |
JP2010099452A (ja) | 2008-09-25 | 2010-05-06 | Fujifilm Corp | 超音波診断装置及び超音波診断方法 |
JP2011172611A (ja) * | 2010-02-23 | 2011-09-08 | Canon Inc | 超音波画像装置および遅延制御方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106999152A (zh) * | 2014-11-26 | 2017-08-01 | Ge医疗系统环球技术有限公司 | 超声诊断设备和程序 |
Also Published As
Publication number | Publication date |
---|---|
CN103930040A (zh) | 2014-07-16 |
JP2013102960A (ja) | 2013-05-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5404141B2 (ja) | 超音波装置及びその制御方法 | |
JP5389722B2 (ja) | 超音波診断装置及びその作動方法 | |
US10918355B2 (en) | Ultrasound diagnostic device and ultrasound diagnostic method | |
JP2010099452A (ja) | 超音波診断装置及び超音波診断方法 | |
WO2012043569A1 (ja) | 超音波画像生成装置、超音波画像生成方法およびプログラム | |
JP5623160B2 (ja) | 超音波診断装置及びその作動方法 | |
US9223011B2 (en) | Ultrasonic signal processing device and ultrasonic signal processing method | |
EP2306186A2 (en) | Ultrasound diagnostic apparatus and signal processing method thereof | |
JP2010234013A (ja) | 超音波診断装置及び超音波診断方法 | |
JP5623157B2 (ja) | 超音波診断装置及びその作動方法 | |
US10788459B2 (en) | Ultrasound diagnostic apparatus, ultrasound image generation method, and recording medium | |
JP2012071037A (ja) | 超音波画像診断装置、超音波画像生成方法およびプログラム | |
JP5247330B2 (ja) | 超音波信号処理装置及び超音波信号処理方法 | |
WO2013073514A1 (ja) | 超音波診断装置及び方法 | |
JP2013244159A (ja) | 超音波診断装置及び音速推定方法 | |
JP2013102959A (ja) | 超音波診断装置及び方法 | |
JP2014064856A (ja) | 超音波検査装置、超音波検査装置の信号処理方法およびプログラム | |
JP5346987B2 (ja) | 超音波診断装置 | |
WO2013176255A1 (ja) | 超音波診断装置およびデータ処理方法 | |
JP2013244162A (ja) | 超音波診断装置 | |
JP2013255599A (ja) | 超音波診断装置及び方法 | |
JP2013255598A (ja) | 超音波診断装置及び方法 | |
JP5797235B2 (ja) | 超音波診断装置及び超音波診断装置の作動方法 | |
JP5681755B2 (ja) | 超音波診断装置及びその作動方法 | |
JP2013244161A (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12849751 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012849751 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |