WO2014050797A1 - 超音波検査装置、超音波検査方法、プログラム及び記録媒体 - Google Patents
超音波検査装置、超音波検査方法、プログラム及び記録媒体 Download PDFInfo
- Publication number
- WO2014050797A1 WO2014050797A1 PCT/JP2013/075666 JP2013075666W WO2014050797A1 WO 2014050797 A1 WO2014050797 A1 WO 2014050797A1 JP 2013075666 W JP2013075666 W JP 2013075666W WO 2014050797 A1 WO2014050797 A1 WO 2014050797A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- sound speed
- element data
- unit
- ultrasonic
- sound
- Prior art date
Links
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
- A61B8/145—Echo-tomography characterised by scanning multiple planes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4444—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to the probe
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- 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/5269—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8997—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using synthetic aperture techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4483—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
- A61B8/4488—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
Definitions
- the present invention relates to an ultrasonic inspection apparatus that performs imaging of an inspection target such as an organ in a living body by transmitting and receiving an ultrasonic beam, and generates an ultrasonic image used for inspection and diagnosis of the inspection target.
- the present invention relates to an ultrasonic inspection method, a program, and a recording medium.
- an ultrasonic inspection apparatus such as an ultrasonic diagnostic imaging apparatus using an ultrasonic image
- this type of ultrasonic inspection apparatus has an ultrasonic probe (ultrasonic probe) including a plurality of elements (ultrasonic transducers), and an apparatus main body connected to the ultrasonic probe.
- the ultrasonic probe transmits the ultrasonic beam from the multiple elements of the ultrasonic probe toward the inspection object (subject), receives the ultrasonic echo from the subject, and receives the ultrasonic echo.
- An ultrasonic image is generated by electrically processing the ultrasonic echo signal thus processed in the apparatus main body.
- an ultrasonic wave is focused on a region to be inspected of a subject, for example, an organ in a living body or a lesion in the organ from a plurality of elements of the probe. Transmits a beam and receives ultrasonic echoes from the surface or interface of a reflector in the examination target area, for example, an organ or a lesion, via multiple elements, but is reflected by the same reflector.
- a reflector located at the focal position of the ultrasonic beam transmitted from the transmitting element, and reflected by the same reflector with respect to the ultrasonic echo signal received by the transmitting element. Therefore, the ultrasonic echo signal received by another element different from the transmitting element is delayed.
- the ultrasonic echo signals received by a plurality of elements are A / D (analog / digital) converted into element data, and then the element data is subjected to reception focus processing, that is, delay correction. Then, the phase is matched and phased and added to generate a sound ray signal, and an ultrasonic image is generated based on the sound ray signal thus obtained.
- a virtual point sound source is formed by focusing transmission ultrasonic waves radiated from a plurality of vibration elements constituting a transmission vibration element group on a transmission focusing point.
- Received ultrasonic waves reflected from a plurality of continuous observation points by transmitted ultrasonic waves radiated from a sound source are received by a plurality of vibration elements constituting a reception vibration element group, and the received signals for the obtained channels are received.
- An ultrasonic diagnostic apparatus that performs reception phasing and addition such that an observation point becomes a reception focusing point is disclosed. Further, in this ultrasonic diagnostic apparatus, the same reception phasing addition is applied to a reception signal obtained by using each of the reception vibration element group and the transmission vibration element group sequentially shifted in the arrangement direction of the vibration elements.
- the transmission phasing addition for correcting the transmission delay due to the difference in propagation distance from each transmission focusing point to the observation point is performed on the reception signal after the reception phasing addition.
- Patent Document 1 discloses that an ultrasonic diagnostic apparatus can generate and display image data having excellent spatial resolution, contrast resolution, and S / N.
- Patent Document 1 see the embodiment shown in FIG. 5
- three transmitting vibration element groups are focused at respective transmission convergence points (hereinafter referred to as focal points) Ft1 to Ft3.
- the focal points are regarded as virtual point sound sources (hereinafter referred to as virtual point sound sources)
- the wave front from each virtual point sound source is reflected at the same observation point (reflection point) Px position. It is shown that the reflected wave obtained by the above is received by the receiving vibration element group.
- the signal quality can be improved as compared with the conventional technique by adding together signals transmitted from a single point at a plurality of different focal points.
- An object of the present invention is to solve the above-described problems of the prior art and to superimpose element data acquired by transmission from two or more different elements when transmitting ultrasonic beams from a plurality of elements of an ultrasonic probe.
- line processing By performing line processing, the influence of the spread of the transmitted ultrasonic beam can be reduced, the S / N ratio can be increased, and the resolution can be increased.
- the rate While maintaining the rate, it is possible to obtain a sharp ultrasound image with optimal spatial resolution at a high resolution that does not depend on the width of the transmitted beam, and by correcting the sound speed after including calculation of multiline processing, Ultrasonic inspection apparatus, ultrasonic inspection method, program, and recording medium capable of obtaining high-quality ultrasonic image with reduced influence of sound velocity unevenness in living tissue It is to provide.
- the present inventor has achieved image quality by performing multi-line processing when transmitting / receiving ultrasonic beams from a plurality of elements of an ultrasonic probe.
- An ultrasonic inspection apparatus that can obtain a good ultrasonic image is proposed.
- the technique disclosed in the previous application describes the sound speed correction for correcting the unevenness of the sound speed distribution. Only the sound speed from the observation point Px) to the probe in FIG. 6 is corrected, and the sound speed from the transmitting element forming each virtual point sound source to the reflection point (Px point) is not corrected. The present inventors have found that the correction of sound speed is insufficient and have arrived at the present invention.
- the ultrasonic inspection apparatus of the present invention is an ultrasonic inspection apparatus that inspects an inspection object using an ultrasonic beam, and transmits an ultrasonic beam and is reflected by the inspection object. And outputs an analog element signal corresponding to the received ultrasonic echo, and a probe in which a plurality of elements are arranged and a plurality of elements are used as the probe to form a predetermined transmission focus
- the transmission unit that transmits the ultrasonic beam a plurality of times, and the reception of the analog element signals output from the plurality of elements in response to the transmission of the individual ultrasonic beams, and performing a predetermined process
- An A / D converter that converts the analog element signal processed by the receiving unit into first element data that is a digital element signal, and a sound speed determining unit that determines a sound speed value in the inspection object , The sound velocity value in the inspection object
- an element data processing unit that generates second element data corresponding to any one of the two or more first element data from the two or more first element data, and
- the sound speed determination unit includes an element data processing unit, and obtains an optimum sound speed by optimizing a sound speed value used for delay time correction calculation of two or more first element data in the element data processing unit.
- a sound speed determination part determines the sound speed value in each area
- an element data holding unit that holds a part of the first element data including two or more first element data or all the first element data.
- the element data holding unit transmits an ultrasonic beam for each of at least two overlapping target regions in the inspection target by the transmission unit, and at least two in the inspection target by the ultrasonic beam by the reception unit.
- Two or more first element data including reception data for each reception time in each element received for each element data generated by receiving an ultrasonic echo signal generated for each of two or more overlapping target areas is held.
- the element data processing unit has two or more second elements based on the geometric arrangement of the elements when transmitting the ultrasonic beam for each element data and the geometric arrangement of each element of the reception data for each element data. It is preferable to generate the second element data from the element data of one.
- a transmission part performs transmission of an ultrasonic beam in multiple times, changing the element used as a center.
- the element data processing unit calculates a delay time calculation unit that calculates each delay time of the two or more first element data using the sound velocity value in the inspection object, and the two or more first element data. It is preferable to have a superimposition processing unit that generates the second element data by superimposing based on the calculated delay times and the received element positions of the probe. And an image generation unit that generates an ultrasonic image based on the second element data generated by the element data processing unit. The element data processing unit It is preferable to create one line of second element data on the sonic image.
- the sound speed determination unit uses the element data processing unit, the sound speed changing unit that changes the sound speed value of the inspection object used for delay time calculation in the element data processing unit, and the sound speed value changed by the sound speed changing unit.
- An element data processing unit that generates second element data, and a second element data generated by the element data processing unit are phased and added using the sound velocity value changed by the sound velocity changing unit to generate a sound ray signal.
- An index calculation unit, a determination unit that determines the image quality index of the ultrasonic image calculated by the image quality index calculation unit, and a sound speed setting unit that sets an optimal sound speed value according to the determination result of the image quality index by the determination unit It is preferable to have.
- the sound speed determination unit uses the sound speed value changed by the element data processing unit every time the sound speed change unit changes the sound speed value from the initial sound speed value to the final sound speed value at predetermined sound speed intervals.
- the element data is created, and the phasing / adding unit generates the sound ray signal by phasing and adding the second element data using the same sound velocity value, and the image generation unit generates the ultrasonic image based on the sound ray signal.
- the image quality index calculation unit determines an optimum image quality index from among a plurality of image quality indexes calculated for all the plurality of sound speed values, and the sound speed setting unit determines the optimum image quality index determined by the determination unit. It is preferable to set the corresponding sound speed value as the optimal sound speed value.
- a second sound speed determination unit that calculates an optimal sound speed value using the first element data is provided before the sound speed determination unit, and the sound speed determination unit is a sound speed calculated by the second sound speed determination unit. The value is preferably used as the initial value.
- the ultrasonic inspection method of the present invention transmits an ultrasonic beam, receives an ultrasonic echo reflected by the inspection object, and responds to the received ultrasonic echo.
- the ultrasonic beam is transmitted so as to form a plurality of times, and in response to the transmission of the individual ultrasonic beams, a plurality of elements output analog element signals.
- 1st element data that is a digital element signal by performing D / D conversion, and corresponds to one of the first element data from two or more first element data using the sound velocity value in the inspection object Create second element data And determining as an optimal sound speed value to optimize the sound speed value to be used when.
- the second element data is created using the changed sound speed value,
- the second element data is phased and added using the same sound velocity value to generate a sound ray signal, and an ultrasonic image is generated based on the generated sound ray signal.
- a plurality of image quality indices corresponding to each of the plurality of changed sound speed values are calculated, and a plurality of image quality indices calculated for all of the plurality of sound speed values are calculated. It is preferable to determine an optimal image quality index from among them, and set a sound speed value corresponding to the determined optimal image quality index as an optimal sound speed value.
- the program of the present invention transmits an ultrasonic beam and receives an ultrasonic echo reflected by the inspection object to inspect the inspection object.
- the probe in which the elements are arranged is used a plurality of times to transmit the ultrasonic beam so as to form a predetermined transmission focal point, and the transmission of individual ultrasonic beams is performed.
- a plurality of elements output analog element signals
- the sound speed value used when creating the second element data corresponding to one of the first element data from the two or more pieces of first element data obtained is optimized to obtain the optimum sound speed value.
- the recording medium of the present invention is a computer-readable recording medium that records the above-described program.
- multi-line processing when transmitting an ultrasonic beam from a plurality of elements of an ultrasonic probe, multi-line processing is performed to superimpose element data acquired by transmission from two or more different elements, thereby transmitting ultrasonic waves.
- the width of the transmitted beam without reducing the influence of the beam spread, increasing the signal-to-noise ratio and increasing the resolution, without using a dedicated ultrasonic probe for wide beam transmission, and at the same frame rate
- High-resolution and sharp ultrasound images with optimal spatial resolution can be obtained, and by correcting the sound speed including calculation of multiline processing, the effect of uneven sound speed in living tissue can be reduced.
- a high-quality ultrasonic image can be obtained.
- FIG. 5 is an explanatory diagram showing element data obtained respectively.
- (A) And (c) is explanatory drawing in the case of transmitting an actual ultrasonic beam from the element directly above the reflection point of the subject and the element not directly above, respectively, (b) and (d) It is explanatory drawing which shows the element data obtained, respectively.
- (A) And (b) is explanatory drawing explaining the distance of the transmission path
- (A), (b) and (c) and (d), (e) and (f) are respectively the element data obtained by a plurality of elements in the case of a true signal and the case of a ghost, their delay time and It is explanatory drawing which shows the superimposition state of element data, (g) And (h) is explanatory drawing which shows the superimposition state of the element data corresponding to a some element, respectively, and its result. It is explanatory drawing explaining transmission / reception of the ultrasonic wave in the ultrasonic diagnostic apparatus of a prior art. It is a block diagram which shows notionally an example of a structure of Example 2 of the ultrasonic inspection apparatus which concerns on this invention. It is a block diagram which shows the principal part containing the sound speed correction
- FIG. 1 is a block diagram conceptually showing an embodiment of the configuration of the ultrasonic inspection apparatus of the present invention.
- the ultrasonic inspection apparatus 10 includes an ultrasonic probe 12, a transmission unit 14 and a reception unit 16 connected to the ultrasonic probe 12, an A / D conversion unit 18, and an element data holding unit 20.
- a sound speed correction unit 24 including an element data processing unit 22, an image generation unit 26, a display control unit 27, a display unit 28, a control unit 30, an operation unit 32, and a storage unit 34.
- the ultrasonic probe 12 has a probe 36 used in a normal ultrasonic inspection apparatus.
- the probe 36 has a plurality of elements arranged in a one-dimensional or two-dimensional array, that is, ultrasonic transducers. These ultrasonic transducers transmit an ultrasonic beam to a subject in accordance with a drive signal supplied from the transmission unit 14 when an ultrasonic image of an object to be examined (hereinafter referred to as a subject) is captured. Receives ultrasonic echoes from the specimen and outputs a received signal.
- each of a predetermined number of ultrasonic transducers constituting a set of the plurality of ultrasonic transducers of the probe 36 generates each component of one ultrasonic beam, and a set of a predetermined number of ultrasonic transducers.
- the ultrasonic transducer generates one ultrasonic beam that is transmitted to the subject.
- Each ultrasonic transducer is, for example, a piezoelectric ceramic represented by PZT (lead zirconate titanate), a polymer piezoelectric element represented by PVDF (polyvinylidene fluoride), or PMN-PT (magnesium niobate / lead titanate). It is constituted by an element in which electrodes are formed at both ends of a piezoelectric body made of a piezoelectric single crystal or the like typified by a solid solution, that is, a vibrator. That is, the probe 36 can be referred to as a transducer array in which a plurality of transducers are arranged in a one-dimensional or two-dimensional array as a plurality of ultrasonic elements.
- PZT lead zirconate titanate
- PVDF polyvinylidene fluoride
- PMN-PT magnesium niobate / lead titanate
- each transducer When a pulsed or continuous wave voltage is applied to the electrodes of such a vibrator, the piezoelectric material expands and contracts, and pulse or continuous wave ultrasonic waves are generated from the respective vibrators, and the synthesis of these ultrasonic waves. As a result, an ultrasonic beam is formed.
- each transducer generates an electric signal by expanding and contracting by receiving propagating ultrasonic waves, and these electric signals are output as ultrasonic reception signals.
- the transmission unit 14 includes, for example, a plurality of pulsars, and according to the sound speed or the distribution of sound speeds set based on the transmission delay pattern selected according to the control signal from the control unit 30, one of the probes 36.
- a plurality of sets formed by adjusting the delay amount of each drive signal so that ultrasonic beam components transmitted from a predetermined number of ultrasonic transducers (hereinafter referred to as ultrasonic elements) form one ultrasonic beam.
- ultrasonic elements formed by adjusting the delay amount of each drive signal so that ultrasonic beam components transmitted from a predetermined number of ultrasonic transducers (hereinafter referred to as ultrasonic elements) form one ultrasonic beam.
- the receiving unit 16 receives, from the subject, an ultrasonic echo generated by the interaction between the ultrasonic beam and the subject by each ultrasonic element of the probe 36 in accordance with a control signal from the control unit 30. Then, the reception signal, that is, the analog element signal for each ultrasonic element is amplified and output, and the amplified analog element signal is supplied to the A / D converter 18.
- the A / D converter 18 is connected to the receiver 16 and converts the analog element signal supplied from the receiver 16 into digital element data.
- the A / D conversion unit 18 supplies the A / D converted digital element data to the element data holding unit 20.
- the element data holding unit 20 sequentially stores the digital element data output from the A / D conversion unit 18. Further, the element data holding unit 20 receives information on the frame rate input from the control unit 30 (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) from the above-described digital element data (hereinafter referred to as the “element data”). , Simply referred to as element data) and stored.
- the element data holding unit 20 inspects at least two or more target regions that overlap in a target region on two or more dimensional position coordinates in accordance with a control signal from the control unit 30, the two or more target regions Each of the two or more element data generated for each two or more target regions from the ultrasonic echo received by the receiving unit 16 is stored and held, and the element data stored and held in the element data holding unit 20 is Two or more element data including reception data for each reception time in each element received for each element data. That is, the element data holding unit 20 holds a part of element data including two or more element data or all the first element data.
- the element data processing unit 22 is provided in the sound speed correction unit 24, and based on control by the control unit 30, the calculation coordinates of a predetermined calculation region in the examination target region of the subject (sampling points in the predetermined region). Or a collection of observation points) is subjected to multi-line processing using the sound speed of the calculated coordinates (hereinafter also referred to as environmental sound speed), that is, for each of two or more target regions stored and held by the element data holding unit 20
- Two or more element data (hereinafter referred to as “unprocessed element data”) generated at the same time are superimposed on the reception time in accordance with the time to generate element data after overlapping processing (hereinafter referred to as “processed element data”). is there.
- the sound speed of the calculated coordinates or the environmental sound speed is the transmission path from the ultrasonic element transmitted from the ultrasonic element of the probe 36 to the predetermined calculated coordinates in the examination target region of the subject,
- the average sound speed in the reception path until the ultrasonic echo reflected by the calculated coordinates returns to the ultrasonic element of the probe 36, that is, the average sound speed between the ultrasonic element of the probe 36 and the calculated coordinates.
- the sound speed correction unit 24 is provided between the element data holding unit 20 and the phasing addition unit 38 of the image generation unit 26, and includes an element data processing unit 22 that performs multiline processing therein. Based on the control, it is used by the phasing addition unit 38 based on the processed data obtained by performing the multiline processing by the element data processing unit 22 for each calculation coordinate of the predetermined calculation region in the examination target region of the subject.
- the optimum sound speed of transmission / reception is obtained, the ambient sound speed in the examination target area of the subject, and hence the optimum environmental sound speed for the ultrasonic image of the examination target area of the subject, and the optimum environment by the element data processing unit 22 are obtained. Processed data that has been subjected to multiline processing using the speed of sound is obtained. Details of the element data processing unit 22 and the sound speed correction unit 24 will be described later.
- the image generation unit 26 generates a sound ray signal (reception data) from the processed element data generated by the element data processing unit 22 and supplied from the sound speed correction unit 24 under the control of the control unit 30. An ultrasonic image is generated from the signal.
- the image generation unit 26 includes a phasing addition unit 38, a detection processing unit 40, a DSC 42, an image creation unit 44, and an image memory 46.
- the phasing / adding unit 38 processes the element data generated by the element data processing unit 22 in accordance with the optimum environmental sound speed obtained by the sound speed correction unit 24 or its distribution according to the reception direction set by the control unit 30.
- the reception focus processing is performed by adding the respective delays to and adding. By this reception focus processing, reception data (sound ray signal) in which the focus of the ultrasonic echo is narrowed is generated.
- reception delay pattern is selected from a plurality of reception delay patterns stored in advance, and the selected reception delay pattern is selected. It is also possible to use the ambient sound speed set based on
- the phasing addition unit 38 supplies the received data to the detection processing unit 40.
- the detection processing unit 40 corrects attenuation according to the distance according to the depth of the reflection position of the ultrasonic wave on the reception data generated by the phasing addition unit 38, and then performs envelope detection processing to perform detection.
- B-mode image data that is tomographic image information related to the tissue in the specimen is generated.
- a DSC (digital scan converter) 42 converts (raster conversion) the B-mode image data generated by the detection processing unit 40 into image data in accordance with a normal television signal scanning method.
- the image creation unit 44 performs various necessary image processing such as gradation processing on the B-mode image data input from the DSC 42 to create B-mode image data for use in inspection and display, and then creates the created inspection. Or display B-mode image data is output to the display control unit 27 for display or stored in the image memory 46.
- the image memory 46 temporarily stores the inspection B-mode image data created by the image creation unit 44.
- the inspection B-mode image data stored in the image memory 46 is read by the display control unit 27 for display on the display unit 28 as necessary.
- the display control unit 27 causes the display unit 28 to display an ultrasonic image based on the inspection B-mode image signal that has been subjected to image processing by the image creation unit 44.
- the display unit 28 includes a display device such as an LCD, for example, and displays an ultrasonic image under the control of the display control unit 27.
- the control unit 30 controls each unit of the ultrasonic inspection apparatus 10 based on a command input from the operation unit 32 by the operator.
- the control unit 30 performs various types of information, particularly information necessary for calculating the optimum environmental sound speed used in the sound speed correction unit 24, the optimum environment in the phasing addition unit 38 of the element data processing unit 22 and the image generation unit 26.
- Information necessary for calculating the delay time used together with the speed of sound and information necessary for element data processing in the element data processing unit 22 are read from the storage unit 34 and / or input by the operator via the operation unit 32.
- the above-described various information input from the operation unit 32 and / or read out from the storage unit 34 is output to the transmission unit 14, the reception unit 16, and the element data holding unit 20 as necessary.
- the operation unit 32 is for an operator to perform an input operation, and can be formed from a keyboard, a mouse, a trackball, a touch panel, and the like.
- the operation unit 32 provides various information as required by the operator, in particular, a plurality of ultrasonic elements of the probe 36 of the probe 12 used for the above-described delay time calculation, an examination target region of the subject, and calculation thereof.
- An input device is provided for inputting and operating information related to the position, transmission aperture and reception aperture of the probe 36, and information related to multi-line processing of element data such as the number of overlapping element data and the overlay processing method. .
- the storage unit 34 receives various types of information input from the operation unit 32, in particular, information related to the probe 12, calculation area, sound velocity search range, focal position, transmission aperture, reception aperture, and the like, the number of superimposing element data, and superposition.
- Information related to multi-line processing of element data such as processing method, etc., and transmission unit 14, reception unit 16, element data holding unit 20, element data processing unit 22, sound speed correction unit 24, image generation unit 26, display control unit 27, etc.
- the information required for the processing and operation of each unit controlled by the control unit 30 and the operation program and processing program for executing the processing and operation of each unit are stored.
- Hard disk, flexible disk, MO, MT Recording media such as RAM, CD-ROM, and DVD-ROM can be used.
- the element data processing unit 22, the sound speed correction unit 24, the phasing addition unit 38 of the image generation unit 26, the detection processing unit 40, the DSC 42, the image creation unit 44, and the display control unit 27 are a CPU and a CPU.
- the program is composed of operation programs for performing processing, but may be configured by a digital circuit.
- FIG. 2 is a block diagram showing a main part of the ultrasonic inspection apparatus including details of an example of the sound velocity correction unit of the ultrasonic inspection apparatus of the present invention shown in FIG. 1 along the processing flow.
- the sound speed correction unit 24 changes the entire predetermined sound speed range in which the set sound speed V is estimated as the sound speed in the inspection target region with a predetermined sound speed value change width (sound speed interval), for example, for a living body such as a human body.
- the image in the inspection target area is changed from 1400 m / s to 1650 m / s or 1700 m / s by a predetermined step (step) of 1 m / s to 50 m / s, for example, 10 m / s.
- the optimum sound speed is comprehensively obtained based on the image quality index (image quality evaluation index) such as the brightness value and the sharpness of the sound and the degree of convergence of the sound speed.
- the sound speed correcting method disclosed in Japanese Patent Application Laid-Open No. 2011-92686 can be used as the sound speed correcting method by the sound speed correcting unit 24.
- the sound speed correction method disclosed in this publication it is necessary to establish a transmission focus for each region of interest.
- processed element data that can achieve transmission focus at any depth can be obtained by superimposing processing of the superimposition processing unit of the element data processing unit 22 in the sound velocity correction unit 24. Therefore, it is not necessary to set the transmission focus for each calculation area or attention area.
- the transmission focus may be set for each calculation area or area of interest as in the sound speed correction method disclosed in this publication.
- the sound speed correction unit 24 includes a calculation area setting unit 70, an initial calculation coordinate setting unit 72, an optimum sound speed determination unit 23, a final calculation coordinate determination unit 88, and a calculation coordinate update unit 90.
- the optimum sound speed determination unit 23 includes a sound speed search range setting unit 74, an initial sound speed value setting unit 76, a phasing addition unit 38a, a detection processing unit 40a, an image creation unit 44a, and an image quality index calculation.
- the sound speed correction unit 24 changes the set sound speed V at the calculated coordinates X1 to Xend in the inspection target region, and from the element data holding unit 20 based on the set speed V at each set speed V.
- the phasing addition unit 38a uses the processed element data obtained by multi-line processing the supplied unprocessed element data by the element data processing unit 22 to perform reception focus processing to perform phasing addition to generate a sound ray signal.
- An ultrasonic image is formed in the image creation unit 44a from the sound ray signal, and image quality indices such as a luminance value and a sharpness of the ultrasonic image in the inspection target area at each set sound velocity V are calculated.
- the image quality index such as the brightness value and sharpness of the ultrasonic image at each set sound speed V is compared in the optimum sound speed determination unit 84.
- the image quality index such as the obtained brightness value and sharpness is obtained.
- the change of the sound speed value and the setting of the optimum sound speed may be performed comprehensively as shown in FIG.
- the area and calculation coordinates for setting the sound speed value are not particularly limited, and the entire subject or the entire inspection target area is set to the same sound speed value.
- the sound speed value may be set for each minute area by dividing the subject or the inspection target area for each minute area, or the sound speed value may be set for each pixel.
- the element data processing unit 22 of the sound speed correction unit 24 shown in FIG. 2 will be described in detail based on FIG.
- the element data processing unit 22 is for performing multiline processing, and includes a delay time calculation unit 48 and an overlay processing unit 50.
- the sound speed value Vst or the sound speed value V updated by the sound speed value update unit 82 (hereinafter, represented by the sound speed value X) is set in advance. Further, it is assumed that unprocessed element data related to the calculation coordinate Xi necessary for the multiline processing in the element data processing unit 22 is read from the element data holding unit 20 and supplied.
- the delay time calculation unit 48 receives a plurality of ultrasonic elements and ultrasonic beams of the probe 36 of the probe 12 input from the operation unit 32 or input from the operation unit 32 and stored in the storage unit 34. Information relating to the focal position, the transmission aperture and the reception aperture of the probe 36 is acquired in advance, and the sound speed of the calculated coordinate Xi in the examination target region of the subject is determined by the initial sound speed value setting unit 76 or the sound speed value update unit. An ultrasonic element (transmission element) of a transmission aperture that forms and transmits an ultrasonic beam using a sound velocity value V (Vst) set in advance by 82, and an ultrasonic echo by an ultrasonic beam from a subject are transmitted. The delay time of the element data received by each ultrasonic element of the reception aperture is calculated based on the geometric arrangement with the ultrasonic element of the reception aperture.
- the superposition processing unit 50 is a multi-element data unit such as the number of element data to be superposed and the superposition processing method input from the operation unit 32 or input from the operation unit 32 and stored in the storage unit 34. Based on the information regarding the line processing, two or more unprocessed element data generated for each of the two or more target areas stored and held by the element data holding unit 20 are read, and the delay times calculated by the delay time calculation unit 48, respectively. Based on the above, two or more unprocessed element data are superimposed on the reception time, that is, by matching the time and the absolute position of the elements of the received probe. Generate.
- an ultrasonic beam (hereinafter simply referred to as a transmission beam) is transmitted from a transmission ultrasonic element (hereinafter simply referred to as a transmission element) of the probe 36 of the ultrasonic probe 12 to the subject, and a mutual relationship with the subject is obtained.
- a transmission ultrasonic element hereinafter simply referred to as a transmission element
- element data it is obtained by the transmission beam from the transmitting element and the receiving element.
- ultrasonic data is received by receiving ultrasonic echoes using seven ultrasonic elements (hereinafter also simply referred to as elements) 52a to 52g and 52b to 52h as receiving elements.
- elements ultrasonic elements
- the transmission beam 56 transmitted to the inspection target area including the reflection point 54 is ideally narrowed to be equal to or smaller than the element interval, as shown in FIG.
- the element data is transmitted.
- the focal point 58 of the beam 56 is on a straight line connecting the element 54d and the reflection point 54. Since the transmission beam 56 is transmitted to the reflection point 54, an ultrasonic echo reflected from the reflection point 54 is generated. The ultrasonic echoes from the reflection point 54 are received by the receiving elements 52a to 52g through the receiving path 60 spreading at a predetermined angle, and the element data 62 as shown in FIG. 4B is obtained by the receiving elements 52a to 52g. Will be.
- the center of the transmitting element is shifted by one element with respect to the reflection point 54 in the direction of the element (right direction in the figure), and just above the reflection point 54.
- the transmission direction of the transmission beam 56 that is, the transmission element 52e and the focal point 58 are received. Therefore, the transmission beam 56 is not transmitted to the reflection point 54. For this reason, the ultrasonic echo reflected from the reflection point 54 is not generated, and the receiving elements 52b to 52h do not receive the ultrasonic echo, so that element data cannot be obtained as shown in FIG. become.
- the actual transmission beam 64 is wider than the element spacing.
- the transmission beam 64 is transmitted as in FIG. 4A. Is wide, the focal point 58 is on a straight line connecting the element 54d and the reflection point 54, and the transmission beam 64 is reflected by the reflection point 54 to generate an ultrasonic echo.
- the ultrasonic echo from the reflection point 54 is received by the receiving elements 52a to 52g through the receiving path 60 spreading to a predetermined angle, and is received by the receiving elements 52a to 52g.
- True element data 66 as shown in FIG. 5B is obtained.
- the center of the transmitting element is shifted by one element with respect to the reflection point 54 in the element direction (right direction in the figure).
- the transmission beam 64 is transmitted using the element 52e adjacent to the element 52d immediately above the reflection point 54 as a transmission element and the ultrasonic echoes are received by the reception elements 52b to 52h, the transmission beam 64 is wide. Even if the reflection point 54 does not exist in the transmission direction, that is, on the straight line connecting the transmission element 52e and the focal point 58, the transmission beam 64 is transmitted to the reflection point 54.
- an ultrasonic echo that does not exist originally that is, a so-called ghost reflection signal is generated from the reflection point 54, and the ghost reflection signal from the reflection point 54 passes through the reception path 60 that spreads to a predetermined angle to the reception elements 52b to 52h.
- the ghost element data 68 as shown in FIG. 5 (d) is obtained by the reception elements 52b to 52h.
- Such ghost element data 68 causes a decrease in the accuracy of the ultrasonic image generated from the element data.
- the transmission beam 64 shown in FIG. 5C is transmitted from the transmission element 52e via the focal point 58 to the reflection point 54, and the ghost reflection signal is received from the reflection point 54 to each of the reception elements 52b to 52h.
- the sum (propagation distance) with the path is such that the transmission beam 64 shown in FIG. 5A reaches the reflection point 54 from the transmission element 52d via the focal point 58 and the true reflected ultrasonic echo is the reflection point 54. 5 is longer than the sum (propagation distance) with the receiving path from each of the receiving elements 52a to 52g to the ghost element data 68 as shown in FIG. 5D, the true element data 68 as shown in FIG. There will be a delay with respect to the element data 66.
- the delay time calculation unit 48 of the element data processing unit 22 of the present invention the time difference between the ghost element data and the true element data, that is, the delay time is the transmission element, the focal point of the ultrasonic beam, the reflection point of the subject, and Calculated from the geometry of the receiving elements. Therefore, the calculation of the delay time requires information such as the shape of the ultrasonic probe 12 (element spacing, linear, convex, etc.), focal position, transmission aperture, reception aperture, and the like, and the sound velocity value of the examination region of the subject.
- the delay time calculation unit 48 acquires the information input by the operation unit 32 or stored in the storage unit 34, and the sound speed value V set in advance by the sound speed value update unit 82, etc. Perform time calculations.
- the delay time is calculated from, for example, the transmission element, the focal point of the ultrasonic beam, the reflection point of the subject, and the geometric arrangement of the reception element using the sound velocity value V, and reaches the reflection point from the transmission element through the focal point.
- the lengths of the transmission path and the reception path of the transmission beam in the case of the true reflected ultrasonic echo and the reflected signal of the ghost are obtained. be able to.
- the transmitting element 52d and the receiving element 52d coincide with each other, and the focal point 58 and the reflection point 54 are directly below.
- the position of the element 52d directly above the reflection point 54 is the coordinate (x0, 0) on the xy two-dimensional coordinate, the element interval is Le, the position of the focal point 58 is the coordinate (x0, df), the reflection point
- the position of 54 is the coordinate (x0, z)
- the position of the transmitting element 52d is also the same coordinate (x0, 0) as the element 52d immediately above the reflecting point 54, and the reflecting point 54 passes through the focal point 58 from the transmitting element 52d.
- the position of the transmission element 52e is one element lateral to the reflection point 54 as compared to the case of FIG. Rightward)
- the focal point 58 is arranged directly below the transmitting element 52e, but the reflection point 54 is arranged immediately below the receiving element 52d, and the position of the receiving element 52d immediately above the reflection point 54 is As in the case of FIG.
- the value obtained by dividing the ultrasonic propagation distance Lub which is the sum of the distance Ltb of the transmission path 61 and the distance Lrb of the reception path 60 obtained by the geometric arrangement shown in FIG.
- the delay time is calculated by calculating the x coordinate of the reflection point 54 and the transmission element 52 (52e) from the propagation time of the true ultrasonic echo when the reflection point 54 and the x coordinate of the transmission element 52 (52d) match.
- the delay time is obtained from the difference in propagation time of the reflected signal of the ghost when it is shifted by one element interval.
- the transmission path 61 is a model that passes through the focal point 58, but the present invention is not limited to this, for example, it passes through the focal point 58. Alternatively, a route that directly reaches the reflection point 54 may be used.
- the geometric model of Fig.6 (a) and FIG.6 (b) is a case of a linear probe, not only this but another probe can perform the same geometric calculation from the shape of a probe.
- a geometric model can be set from the radius of the probe and the angle between the elements, and the calculation can be performed in the same way.
- a geometric model (not shown) that considers information such as the transmission angle is used, and the element data of the true ultrasonic echo and its surroundings are determined from the positional relationship between the transmission element and the reflection point. The delay time of the ghost element data can be calculated.
- the delay time is not limited to the method of calculating the delay time using the geometric model, and the delay time is obtained for each measurement condition from the measurement result obtained by measuring the high-intensity reflection point according to the measurement condition of the apparatus in advance. By storing in the apparatus, the delay time of the same measurement condition may be read out.
- FIG. 6C shows element data 66 of the true signal at the center and element data 68 of the ghost around the true signal
- FIG. 6D shows element data 66 obtained from the geometric calculation described above.
- An example of the delay time of the ghost element data 68 is shown. It is shown that the element data 68 of the ghost signal is symmetrically delayed with respect to the element data 66 of the true signal. In this way, the delay time calculated by the delay time calculation unit 48 of the element data processing unit 22 can also be used for delay correction in the phasing addition unit 38.
- information on the number of superposition element data and the superposition processing method at the time of superposition is necessary, but these may be input in advance by the operation unit 32. Alternatively, it may be stored in the storage unit 34.
- the range of superimposing element data when superimposing in the superimposing processing unit 50 is desirably element data within a range that remains within the target region of the sound velocity correction calculation set by the calculation region setting unit 70. However, it is not limited by the area set by the calculation area setting unit 70. For example, overlapping may be performed in such a range that element data to be overlapped overlaps between adjacent regions.
- FIGS. 7A to 7H show a specific example of overlay processing performed by the overlay processor 50 when the number of element data is 5 and the number of overlay element data is 3.
- FIG. FIG. 7A shows five element data arranged side by side, and shows a state in which an ultrasonic beam is transmitted and a reflected signal is received for each element data.
- the horizontal axis of each element data represents a receiving element, and the respective element data are displayed with the center element at the time of transmission of the ultrasonic beam as the center.
- the vertical axis represents the reception time.
- the middle element data there is a reflection point directly below the element at the center of the element data (element at the center of the receiving element), that is, the element at the center at the time of transmission (transmitting element).
- a reflection signal from the reflection point is received. That is, the reflected signal is a true signal, and the element data in the middle represents the true signal.
- the transmission element of the middle element data A reflected signal, i.e., a ghost, which is generated when an ultrasonic beam hits a reflection point existing directly below, is reflected. Since the propagation time of the ultrasonic wave to the reflection point becomes longer as the ghost is away from the true signal, the reception time is delayed as compared with the true signal.
- the position of the receiving element where the reflected signal from the reflecting point is first received is the element immediately above the reflecting point, but the horizontal axis of the element data is centered on the central element when transmitting the ultrasonic beam.
- the central element is shifted by one element for each element data
- the absolute position of the element is shifted by one element in each element data. That is, in the middle element data, the receiving element from which the reflected signal from the reflection point is received first is the middle element, but the element data on both sides is shifted by one element from the middle element data.
- the element data is shifted one element to the left, and the left element data is shifted one element to the right. Further, the element data at both ends are shifted by two elements from the middle element data, the leftmost element data is shifted by two elements to the left, and the leftmost element data is shifted by two elements to the right.
- the ghost signal is not only delayed in reception time with respect to the true signal, but also deviated from the direction of the receiving element.
- FIG. 7B shows an example of the delay time of the reception time for the middle element data among the five element data shown in FIG.
- the overlay processing unit 50 uses the delay time shown in FIG. 7B to set the element data of the target element as the element data of the target element.
- the delay time is corrected by three element data, and each element data is horizontally shifted by one element on both sides in the illustrated example in accordance with the difference in element position from the element of interest (difference in the position of the central element).
- the unprocessed element data corresponding to the three element data is overlapped by shifting in the direction, that is, in phase, and obtained as one overlapped element data of the element data of the element of interest.
- FIG. 7C shows element data that has been subjected to the overlay processing of the element data of the element of interest obtained in this way. Since the element data of the target element shown in FIG. 7A is the element data of the true signal, the delay time correction and the lateral shift are performed on the unprocessed element data of the adjacent element data on both sides of the element data of the target element. As shown in FIG. 7C, the unprocessed element data of the adjacent element data and the unprocessed element data of the target element overlap each other at the high luminance position because the phases match. Therefore, when these element data are added, for example, the element data value shows a large value (high luminance value), and for example, even if an average value is obtained by averaging, an emphasized value (high luminance value) is shown.
- FIG. 7D shows the same element data as FIG. 7A, but shows an example in which the element data on the left side of the middle element data, that is, the ghost is the element data of the element of interest.
- FIG. 7 (e) is the same as FIG. 7 (b) and shows an example of the delay time of the reception time for the middle element data of the five element data shown in FIG. 7 (a). That is, since FIG. 7A and FIG. 7D are the same element data, the delay time of the reception time for the element data in the middle of the five element data shown in FIG. In the overlay processing unit 50, the delay time shown in FIG. 7E (that is, the same as FIG.
- FIG. 7F shows element data that has been subjected to the overlay processing of the element data of the element of interest obtained in this way. Since the element data of the target element shown in FIG. 7D is ghost element data, the delay time correction and the lateral shift are performed on the unprocessed element data of the adjacent element data on both sides of the element data of the target element. Even if phase matching is performed, as shown in FIG. 7F, the unprocessed element data of the adjacent element data and the unprocessed element data of the target element do not overlap each other because the phases do not match. For this reason, even if these three element data are added, for example, since the phases are not matched, signals that are inverted in phase cancel each other out, so the added value does not increase. When the average value is obtained, a small value is shown.
- FIG. 7H shows a result obtained by performing, for example, addition processing or averaging processing as superimposing processing on these as shown in g).
- the element data of the true signal has a high luminance value.
- the element data of the ghost adds or averages the element data whose phases are not matched with each other, so they cancel each other out. Therefore, the superposed processed element data of the ghost is smaller than the superposed processed element data having a high luminance value, which is the element data of the true signal, and is smaller than the element data of the true signal.
- the influence of the ghost element data can be reduced, or the influence can be reduced to the extent that the influence can be ignored. Therefore, by performing phasing addition and detection processing on the processed element data, generating reception data and generating an ultrasonic image, similarly, the influence of ghost is eliminated, that is, at all points on the sound ray. Since an ultrasonic image can be generated with element data equal to the focal point, it is possible to generate a high-quality ultrasonic image with high brightness and excellent sharpness.
- the generation of processed element data is referred to as multiline processing.
- the central element is an element at the center in the azimuth direction when the numerical aperture for transmission (the number of elements for transmitting ultrasonic waves) is an odd number.
- the numerical aperture is an even number
- one of the central elements in the azimuth direction is set as the central element, or the central element is assumed to be in the middle of the azimuth direction. That is, when the numerical aperture is an even number, the calculation may be performed assuming that the focal point is on the middle line of the aperture.
- the superposition processing method in the superposition processing unit 50 is not simply the addition, but may be an average value or a median value, or may be added after being multiplied (weighted) by multiplying coefficients. Good. Note that taking an average value or median value is thought to correspond to applying an averaging filter or median filter at the element data level, but is performed by normal image processing instead of the averaging filter or median filter. An inverse filter or the like may also be applied. Alternatively, the element data to be superimposed are compared, and if they are similar, the maximum value is taken, if not, the average value is taken, and if there is a distribution bias, the intermediate value is taken. The overlay process may be changed based on the feature amount of each element data to be superimposed.
- the number of element data to be superimposed on the element data of the element of interest is not limited to two in the illustrated example, and may be one or three or more. That is, the number of element data to be superimposed on the element data of the element of interest may be appropriately set according to the required processing speed (frame rate or the like), image quality, and the like. Further, it is desirable that the number of element data to be superimposed is matched to the extent of the beam width of the ultrasonic beam. Therefore, when the beam width changes depending on the depth, the number of overlapping element data may be changed depending on the depth. Further, since the beam width depends on the transmission numerical aperture, the number of overlapping element data may be changed according to the transmission numerical aperture.
- the number of overlapping element data may be changed based on a feature quantity such as the luminance value of the image, or the optimum number of overlapping element data is selected from images created by changing the number of overlapping element data. May be.
- the signal phase matches in the element data of the true signal, but the signal phase does not match in the ghost.
- superposition processing such as addition, signals of various phases are mutually connected. Cancel each other and the signal becomes weaker.
- the true signal remains as effective element data having a valid value, for example, high luminance
- the ghost signal can be obtained as element data having a reduced value, for example, low luminance.
- the element data processing unit used in the present invention is basically configured as described above.
- the calculation region setting unit 70 is a region to be subjected to sound velocity correction calculation in the inspection target region in accordance with an input from the operation unit 32 (see FIG. 1) by the operator.
- the region of interest ROI is set, all the points to be calculated on the xy coordinates of the set region are extracted and set as the calculation coordinates X1 to Xend.
- each area represented by the calculated coordinates X1 to Xend is not particularly limited as long as the environmental sound speed is considered to be the same, and may be the entire subject or the entire inspection target area.
- the initial calculation coordinate setting unit 72 sets the calculation coordinate Xi to the initial calculation coordinate X1 among the calculation coordinates X1 to Xend set by the calculation region setting unit 70.
- the data of the calculated coordinate X1 set by the initial calculated coordinate setting unit 72 is sent to the optimum sound speed determining unit 23.
- the calculation region setting unit 70 is not limited to the configuration for setting the calculation region in accordance with the input from the operation unit 32, and the ultrasonic wave generated by the image generation unit 44 of the image generation unit 26 (see FIG. 1).
- a region of interest ROI may be set by analyzing an image (B-mode image signal) and used as a calculation region.
- the calculation region setting unit 70 sets a position (pixel) where the difference in luminance value between adjacent pixels is equal to or greater than a predetermined value in the ultrasonic image generated by the image creation unit 44 as a structural region (tissue or lesion).
- the region surrounded by the edge or the region including many edges is determined as the structural region, that is, the target part P, extracted, includes the target part, and
- the region having the shape and the size may be set as the region of interest ROI and set as the calculation region.
- the structure area may be regarded as one kind of area, or two or more kinds according to the continuity of the edge, the level of the luminance value, and the like.
- Each region may be determined as a target region and extracted, and a region of interest ROI may be set for each target region to be a calculation region.
- the optimum sound speed determination unit 23 calculates the optimum environmental sound speed value (hereinafter referred to as the optimum sound speed value) at the initial calculation coordinate X1 set by the initial calculation coordinate setting unit 72 or the calculation coordinate Xi updated by the calculation coordinate update unit 90. This is a part to be calculated and set.
- the optimum sound speed value is obtained by changing the set sound speed (set environment sound speed) V and performing the reception focus processing based on each set sound speed V to form an ultrasound image, / Or the value of the set sound velocity V at which the sharpness is the highest.
- the optimum sound velocity value is determined based on the contrast of the image, the spatial frequency in the scan direction, the variance, etc. It can be performed.
- the optimum sound speed determination unit 23 uses the unprocessed element data related to the calculation coordinates Xi read from the element data holding unit 20 as an input value, and outputs an optimum environmental sound speed value (hereinafter also referred to as an optimum sound speed value) as an output value. It is for exhaustively seeking and setting as. As described above, as the optimum sound speed determination unit 23, a sound speed search range setting unit 74, an initial sound speed value setting, and the like can be exemplified. Unit 76, element data processing unit 22, phasing addition unit 38 a, detection processing unit 40 a, image creation unit 44 a, image quality index calculation unit 78, final sound speed value determination unit 80, and optimum sound speed determination unit 84. And a sound speed value updating unit 82 and a sound speed setting unit 86.
- the sound speed search range setting unit 74 is for setting the search range (initial sound speed value and final sound speed value) of the environmental sound speed and the sound speed value change width (sound speed interval) in the calculated coordinates Xi.
- the sound speed search range setting unit 74 sets, for example, the initial sound speed value Vst at the start of sound speed search, the final sound speed value Vend at the end of sound speed search, and the step width (sound speed change width: sound speed interval) ⁇ V for changing the set sound speed V. .
- the sound speed search range setting unit 74 can set the initial sound speed value Vst to 1400 m / s, the final sound speed value Vend to 1700 m / s, and the sound speed change width ⁇ V to 10 m / s. .
- the initial sound speed value Vst, the final sound speed value Vend, and the sound speed change width ⁇ V set by the sound speed search range setting unit 74 are not limited to the above-described values, and depend on the subject, the region to be examined, the calculated coordinate Xi, and the like. As long as it is determined.
- the initial sound velocity value Vst, the final sound velocity value Vend, and the sound velocity change width ⁇ V are input to the sound velocity search range setting section 74 by the operator according to the subject, the inspection target area, the calculated coordinate Xi, and the like. 32 may be input, or may be read out from those stored in the storage unit 34 in advance.
- the initial sound speed value setting unit 76 sets the set sound speed V to the initial sound speed value Vst set by the sound speed search range setting unit 74. In the case of the above-described living body, the initial sound speed value setting unit 76 sets the set sound speed V as the initial sound speed value Vst to 1400 m / s. The data of the initial sound speed value Vst set by the initial sound speed value setting unit 76 is sent to the element data processing unit 22 at the subsequent stage.
- the phasing / adding unit 38a has the same configuration and function as the phasing / adding unit 38, and with respect to the processed element data related to the calculated coordinate Xi obtained by the element data processing unit 22 using the set sound velocity V. Then, a reception focus process for adding a delay based on the optimum sound speed set by the sound speed correction unit 24 is performed to generate a sound ray signal.
- the detection processing unit 40a has the same configuration and function as the detection processing unit 40a, and performs processing such as envelope detection processing on the sound ray signal generated by the phasing addition unit 38a.
- B-mode image data that is tomographic image information related to the tissue in the specimen is generated.
- the image creation unit 44a has the same configuration and function as the image creation unit 44, and performs image processing on the B-mode image data generated by the detection processing unit 40a to create an ultrasound image.
- the image quality index calculation unit 78 is an image quality index at the set sound speed V from the ultrasonic image created by the image creation unit 44a, that is, an image quality index such as the brightness value and sharpness of the image, for example, calculated coordinates Xi at each set sound speed V. For example, the sharpness of the image in the region of interest ROI is calculated.
- the calculated image quality index is stored in the image quality index calculation unit 78 in association with the set sound speed V, for example, in a memory or the like in the image quality index calculation unit 78.
- the final sound speed value determination unit 80 compares the set sound speed V with the final sound speed value Vend to determine whether or not they are equal, and when the set sound speed V is not equal to the final sound speed value Vend (V ⁇ Vend).
- the optimum sound speed determination unit 84 reads the image quality index obtained for each sound speed value V (Vst to Vend) held in the image quality index calculation unit 78 and reads the image quality index at each set sound speed V, for example, an ultrasonic image.
- the set sound speed V is determined when the image quality is the best, for example, when the ultrasonic image has the highest sharpness (sharpness value).
- the sound speed setting unit 86 sets the set sound speed V determined to have the best image quality of the ultrasonic image generated by the optimum sound speed determination unit 84 as the optimum sound speed value. Note that the optimum sound speed determination unit 23 shown in FIG. 2 repeatedly processes the set sound speed V. However, the present invention is not limited to this.
- the optimum sound speed determination unit 23 sets the optimum sound speed value corresponding to the calculated coordinate Xi in the sound speed setting unit 86, and then proceeds to the final calculated coordinate determination unit 88.
- the optimum sound speed determination unit 23 is basically configured as described above, and calculates the optimum sound speed value related to the calculation coordinates Xi.
- the calculated coordinate update unit 90 updates the current calculated coordinate Xi to a new next calculated coordinate Xi + 1.
- the data of the new next calculated coordinate Xi + 1 updated by the calculated coordinate update unit 90 is sent to the sound speed search range setting unit 74 of the optimum sound speed determination unit 23 in the previous stage.
- the element data processing unit 22 that performs multiline processing is incorporated in the loop for searching for the optimum sound speed in the sound speed correction performed by the sound speed correction unit 24 of the ultrasonic inspection apparatus 10 of the present embodiment.
- the receiving element of the probe 36 of the ultrasonic probe 12 from the reflection point is used. It is possible to obtain an optimum environmental sound speed value including not only the sound speed up to, but also the sound speed from the transmitting element of the probe 36 that forms each virtual point sound source to the reflection point.
- the ultrasonic inspection apparatus of the present invention is basically configured as described above. Hereinafter, with reference to FIG. 1 and FIG. 2, the operation and action of the ultrasonic inspection apparatus of the present invention and the super A sound wave inspection method will be described.
- a plurality of probes 36 of the probe 36 are driven according to the drive signal supplied from the transmission unit 14.
- An ultrasonic beam is transmitted from the element, the ultrasonic echoes from the subject are received by a plurality of elements of the probe 36, and an analog element signal is output as a reception signal.
- the receiving unit 16 amplifies the analog element signal and supplies it to the A / D conversion unit 18, and the A / D conversion unit 18 converts the analog element signal into digital element data and supplies it to the element data holding unit 20. And stored as unprocessed element data.
- the sound velocity correction unit 24 shown in FIG. 2 after the calculation coordinates X1 to Xend are set by the calculation region setting unit 70 in accordance with the input from the operation unit 32 (see FIG. 1) by the operator, the initial values are set.
- the calculated coordinate Xi is set to the initial calculated coordinate X1 by the calculated coordinate setting unit 72, and the data of the set calculated coordinate X1 is sent to the optimum sound speed determining unit 23.
- the sound speed search range setting unit 74 is input by the operator via the operation unit 32 according to the subject, the inspection target region, the calculated coordinate Xi, or the like, or in advance.
- the initial sound speed value Vst, the final sound speed value Vend, and the sound speed change width ⁇ V read from the information on the sound speed search range stored in the storage unit 34 are set.
- the initial sound speed value Vst and the sound speed change width ⁇ V are set.
- the value Vst is set to 1400 m / s
- the sound speed value Vend is set to 1700 m / s
- the sound speed change width ⁇ V is set to 10 m / s.
- the initial sound speed value setting unit 76 sets the set sound speed V to the initial sound speed value Vst (for example, 1400 m / s), and proceeds to the element data processing unit 22.
- the element data processing unit 22 causes the delay time calculation unit 48 (see FIG. 3) to delay the unprocessed element data of the peripheral ghost signal with respect to the unprocessed element data of the true signal (for example, FIG. (b) and FIG. 7 (e) are the same), the geometrical arrangement of the transmitting element, the focal point, the reflection point, and the receiving element, and the set sound speed V set in advance by the initial sound speed value setting unit 76. (E.g., using the geometric model of FIG. 6).
- the element data processing unit 22 reads unprocessed element data related to the calculation coordinate X1 from the element data holding unit 20 from the element data holding unit 20, sets the element data to be processed as target element data, and the overlay processing unit 50. (See FIG. 3), using the delay time calculated by the delay time calculation unit 48, multi-line processing is performed by superimposing the target element data and the surrounding unprocessed element data in phase, and processing the processed element Ask for data. As a result, enhanced processed element data is obtained for unprocessed element data including a true signal, and attenuated processed element data is determined for ghost unprocessed element data. The element data processing unit 22 supplies the processed element data thus obtained to the phasing addition unit 38a.
- the phasing addition unit 38a performs reception focus processing using the processed element data corresponding to the calculated coordinate X1 supplied from the element data processing unit 22 based on the set sound speed V that is the initial sound speed value Vst. Then, phasing addition is performed to generate a sound ray signal.
- the detection processing unit 40a performs detection processing on the generated sound ray signal to generate a B-mode image signal. Then, the image creation unit 44a creates an ultrasonic image from the generated B-mode image signal.
- the image quality index calculation unit 78 determines the image quality index such as the luminance value and sharpness of the ultrasonic image at the set sound speed V from the created ultrasonic image, for example, the sharpness of the Xi (X1) image at each set sound speed V. Degree etc. are calculated.
- the calculated image quality index is stored in the image quality index calculation unit 78 in association with the set sound speed value V.
- the final sound speed value determination unit 80 compares the set sound speed V with the final sound speed value Vend, and determines whether or not both are equal.
- the process proceeds to the optimum sound speed determination unit 84, and the set sound speed V is not equal to the final sound speed value Vend ( If V ⁇ Vend), the procedure proceeds to the sound velocity value update unit 82.
- the process proceeds to the sound speed value update unit 82 and is set based on the new set sound speed V and the new set sound speed V by the sound speed value update unit 82.
- Multi-line processing by the element data processing unit 22, phasing addition by the phasing addition unit 38a, detection processing by the detection processing unit 40a, image creation by the image creation unit 44a, calculation and holding of the image quality index by the image quality index calculation unit 78, and The determination by the final sound velocity determination unit 80 is repeated until V Vend.
- the set sound speed V is changed from 1400 m / s to 1700 m / s in increments of 10 m / s while changing the set sound speed by ⁇ V until the final sound speed value Vend at the end of the sound speed range is reached. Repeat while changing.
- the image quality index is obtained for all the set sound speeds V within the sound speed search range.
- the optimal sound speed determination unit 84 reads the image quality index obtained for each sound speed value held in the image quality index calculation unit 78, and the optimal sound speed determination unit 84 reads the image quality index at each set sound speed V. For example, the ultrasonic image sharpness (sharpness value) is compared, and when the image quality is the best, for example, the set sound speed V when the ultrasonic image sharpness (sharpness value) is the highest is determined, and the sound speed setting unit 86 adopts and sets the set sound speed V determined by the optimal sound speed determination unit 84 as the optimal sound speed value.
- the optimum sound speed determination unit 23 calculates the optimum sound speed value corresponding to the calculated coordinates Xi (X1).
- the process proceeds to the final calculated coordinate determining unit 88.
- the final calculated coordinate determination unit 88 compares the calculated coordinates Xi with the final calculated coordinates Xend to determine whether or not they are equal. As a result of the determination, if the calculated coordinate Xi is not equal to the final calculated coordinate Xend (X ⁇ Xend), the process proceeds to the calculated coordinate update unit 90. Thereafter, the calculated coordinate update unit 90 updates the current calculated coordinate Xi to a new next calculated coordinate Xi + 1, and proceeds to the optimum acoustic velocity determination unit 23 in the previous stage.
- the optimum sound speed determination unit 23 performs an exhaustive search for the set sound speed V described above, calculates the optimum sound speed value corresponding to the calculated coordinates Xi + 1, and moves to the final calculated coordinate determination unit 88 again.
- the final calculation coordinate Xend is determined.
- the process proceeds to the calculation coordinate update unit 90 as described above, the setting of the new calculation coordinate Xi by the calculation coordinate update unit 90, and the optimum sound speed determination unit 23
- the optimum sound speed value calculated at each of the calculated coordinates X1 to Xend and the processed element data processed by the element data processing unit 22 using the optimum sound speed value are adjusted together with the coordinate data of each of the calculated coordinates X1 to Xend. It is supplied to the phase adder 38.
- the phasing addition unit 38 of the image generation unit 26 performs reception focus processing on the element data to generate reception data (sound ray signal) and supplies the reception data (sound ray signal) to the detection processing unit 40.
- the detection processing unit 40 processes the sound ray signal and generates a B-mode image signal.
- the DSC 42 performs raster conversion on the B-mode image signal, and the image creation unit 44 performs image processing to generate an ultrasonic image.
- the generated ultrasonic image is stored in the image memory 46, and the ultrasonic image is displayed on the display unit 28 by the display control unit 27.
- the ultrasonic inspection apparatus of the present invention is an ultrasonic image created using processed element data that has been subjected to multiline processing by the element data processing unit, the ghost signal is attenuated, and the true signal is enhanced.
- the image quality index is used to determine the optimal environmental sound speed, and the ultrasonic image is created using the processed element data that has been multiline processed using the determined optimal environmental sound speed.
- High-quality ultrasound images with reduced effects of unevenness can be obtained, and it depends on the width of the transmitted beam without using an ultrasonic probe dedicated to wide beam transmission and at the same frame rate as before. It is possible to obtain a sharp ultrasonic image with high resolution, high S / N ratio, and optimum spatial resolution.
- FIG. 9 is a block diagram conceptually showing the configuration of the ultrasonic inspection apparatus according to the second embodiment of the present invention.
- FIG. 10 is an example of a main part of the ultrasonic inspection apparatus shown in FIG. FIG.
- the ultrasonic inspection apparatus 100 shown in FIGS. 9 and 10 is exactly the same as the ultrasonic inspection apparatus 10 shown in FIGS. 1 and 2 except that it includes a second sound speed correction unit 102. Since they have a configuration, the same components are denoted by the same reference numerals, and detailed description thereof is omitted.
- the ultrasonic inspection apparatus 100 includes an ultrasonic probe 12, a transmission unit 14 and a reception unit 16 connected to the ultrasonic probe 12, an A / D conversion unit 18, and an element data holding unit 20.
- a storage unit 34 a storage unit 34.
- the second sound speed correction unit 102 performs the sound speed correction using the unprocessed element data (first element data) to provisionally.
- the sound speed value is obtained, and the sound speed value is used as an initial value when the sound speed correction unit 24 performs sound speed correction.
- the second sound speed correction unit 102 is connected between the element data holding unit 20 and the sound speed correction unit 24 and uses unprocessed element data (first element data) stored in the element data holding unit 20. Then, the sound speed value calculated by correcting the sound speed is supplied to the sound speed correcting unit 24.
- the sound speed correction performed by the second sound speed correction unit 102 is not particularly limited, but for example, a method described in Japanese Patent Application Laid-Open No. 2011-92686 can be used. In this method, reception focus processing is performed on unprocessed element data for each of a plurality of set sound speeds, a focus index is calculated for each set sound speed, and an optimum sound speed value is determined based on the focus index. That is, the set sound speed at which the unprocessed element data can be most appropriately received and focused is set to the optimum sound speed value.
- the sound speed correction by the second sound speed correction unit 102 may be roughly obtained as an initial value of the sound speed search at the time of sound speed correction by the sound speed correction unit 24. Therefore, a transmission focus is not necessarily set for each calculation region or region of interest. do not have to. Moreover, it is not limited to what is performed for every calculation area
- the sound speed search range setting unit 74 of the optimum sound speed determination unit 23 sets the sound speed value with the sound speed value supplied from the second sound speed correction unit 102 as the initial sound speed value Vst.
- the sound velocity value update unit 82 sequentially sets the initial sound velocity value Vst and the preceding and following sound velocity values (Vst ⁇ ⁇ V), the image quality index calculation unit 78 calculates the image quality index, and based on the results of the three sound velocity values. Then, the next set sound velocity V is set. For example, if the set sound speed V is Vst + ⁇ V, the sound speed value updating unit 82 sets Vst + 2 ⁇ ⁇ V as the next set sound speed V, and the set sound speed V is Vst.
- Vst is adopted as the optimum sound speed value. That is, the optimum sound speed determination unit 23 searches for the optimum sound speed using the sound speed value supplied from the second sound speed correction unit 102 as an initial sound speed value.
- the second sound speed correction unit 102 performs sound speed correction using the unprocessed element data, and is used as an initial value of sound speed correction in the sound speed correction unit 24, so that the time required for measurement and calculation of the optimum sound speed value is reduced. It can be shortened.
- each component shown in FIGS. 1 and 2 may be configured as hardware, or may be configured as software executed by a computer or the like.
- each component including an element data processing unit included in the above-described sound speed correction unit, and an image generation unit A program for ultrasonic inspection that operates a computer to function also includes each step of the ultrasonic inspection method, for example, each step of sound velocity correction including each step of the multiline processing described above, and each of image generation.
- a program for ultrasonic inspection that causes a computer to execute steps as a procedure is also one embodiment of the present invention.
- a computer-readable storage medium storing such a program is also one embodiment of the present invention.
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Acoustics & Sound (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Gynecology & Obstetrics (AREA)
- Physiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
マルチライン処理の計算を含めた上で音速補正を行うことにより、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができる超音波検査装置、超音波検査方法、プログラム及び記録媒体を提供する。複数の素子を備える探触子と、複数の素子から超音波ビームを検査対象物に送信させる送信部と、検査対象物からの超音波エコー信号を受信する受信部と、検査対象物内の音速値を決定する音速決定部と、その音速値を用いて受信部によって得られる2つ以上の第1の素子データから第1の素子データのいずれかに対応する第2の素子データを作成する素子データ処理部とを有し、音速決定部は、素子データ処理部において第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として求めることにより、上記課題を解決する。
Description
本発明は、超音波ビームを送受信することにより生体内の臓器等の検査対象物の撮像を行って、検査対象物の検査や診断のために用いられる超音波画像を生成する超音波検査装置、超音波検査方法、プログラム及び記録媒体に関する。
従来から、医療分野において、超音波画像を利用した超音波画像診断装置等の超音波検査装置が実用化されている。一般に、この種の超音波検査装置は、複数の素子(超音波トランスデューサ)を内蔵した超音波探触子(超音波プローブ)と、この超音波探触子に接続された装置本体とを有しており、超音波探触子の複数の素子から検査対象物(被検体)に向けて超音波ビームを送信し、被検体からの超音波エコーを超音波探触子で受信して、その受信した超音波エコー信号を装置本体で電気的に処理することにより超音波画像が生成される。
超音波検査装置においては、超音波画像を生成するとき、被検体の検査対象領域、例えば、生体内の臓器やその臓器内の病巣等に探触子の複数の素子から焦点を合わせて超音波ビームを送信し、検査対象領域の反射体、例えば臓器や病巣等の表面や界面からの超音波エコーを複数の素子を介して受信しているが、同一の反射体で反射された超音波エコーを複数の素子で受信するので、送信素子から送信された超音波ビームの焦点位置に位置する反射体で反射され、送信素子で受信された超音波エコー信号に対して、同一の反射体で反射され、送信素子と異なるその他の素子で受信された超音波エコー信号は遅延することになる。このため、超音波検査装置では、複数の素子で受信した超音波エコー信号をA/D(アナログ/デジタル)変換して素子データにとした後、素子データを受信フォーカス処理して、即ち遅延補正して位相を合わせ整相加算して音線信号を生成し、こうして得られた音線信号に基づいて超音波画像を生成している。
このような超音波検査技術において、超音波画像の画質を向上させるために、複数の異なる焦点に収束する複数の超音波ビームの送信で得られた信号を足し合わせることで、従来よりも信号の質を改善することが行われている。
例えば、特許文献1には、送信用振動素子群を構成する複数個の振動素子から放射される送信超音波を送信集束点に集束させることにより仮想的な点音源を形成し、そして、この点音源から放射される送信超音波によって連続した複数の観測点から反射する受信超音波を、受信用振動素子群を構成する複数個の振動素子によって受信し、得られたチャンネル分の受信信号に対し観測点が受信集束点となるような受信整相加算を行う超音波診断装置が開示されている。更に、この超音波診断装置では、受信用振動素子群と振動素子の配列方向に順次シフトさせた送信用振動素子群の各々を用いて得られた受信信号に対しても同様の受信整相加算を行い、これら受信整相加算後の受信信号に対し各々の送信集束点から観測点までの伝搬距離の差異に起因した送信遅延を補正する送信整相加算を行う。
例えば、特許文献1には、送信用振動素子群を構成する複数個の振動素子から放射される送信超音波を送信集束点に集束させることにより仮想的な点音源を形成し、そして、この点音源から放射される送信超音波によって連続した複数の観測点から反射する受信超音波を、受信用振動素子群を構成する複数個の振動素子によって受信し、得られたチャンネル分の受信信号に対し観測点が受信集束点となるような受信整相加算を行う超音波診断装置が開示されている。更に、この超音波診断装置では、受信用振動素子群と振動素子の配列方向に順次シフトさせた送信用振動素子群の各々を用いて得られた受信信号に対しても同様の受信整相加算を行い、これら受信整相加算後の受信信号に対し各々の送信集束点から観測点までの伝搬距離の差異に起因した送信遅延を補正する送信整相加算を行う。
特許文献1に開示の超音波診断装置では、複数の振動素子から得られた受信信号に対して受信整相加算と送信整相加算を行なうことにより、被検体の深さ方向に対してほぼ一様な細いビーム幅を有した送信ビーム及び受信ビームを高精度かつ高感度で形成することができる。このため、特許文献1は、超音波診断装置において、空間分解能、コントラスト分解能及びS/Nに優れた画像データの生成と表示が可能となることを開示している。
特に、特許文献1(その図5に示す実施例参照)には、図8に示すように、3つの送信用振動素子群がそれぞれの送信収束点(以下、焦点とする)Ft1~Ft3で焦点を結び、各焦点を仮想的な点音源(以下、仮想点音源とする)と見なした場合に、各仮想点音源からの波面が同一の観測点(反射点)Pxの位置で反射することによって得られた反射波を受信用振動素子群で受信することが示されている。特許文献1では、これにより、1つの点を複数の異なる焦点で送信した信号を足し合わせることで、従来よりも信号の質を改善させることができるとしている。
特に、特許文献1(その図5に示す実施例参照)には、図8に示すように、3つの送信用振動素子群がそれぞれの送信収束点(以下、焦点とする)Ft1~Ft3で焦点を結び、各焦点を仮想的な点音源(以下、仮想点音源とする)と見なした場合に、各仮想点音源からの波面が同一の観測点(反射点)Pxの位置で反射することによって得られた反射波を受信用振動素子群で受信することが示されている。特許文献1では、これにより、1つの点を複数の異なる焦点で送信した信号を足し合わせることで、従来よりも信号の質を改善させることができるとしている。
しかしながら、特許文献1に開示の技術では、従来技術よりも高画質な画像が得られるが、1ラインのデータを作るのに複数の送信ビームを、送信位置を変えて発生させる必要があり、従来技術よりも送信回数が増えるためフレームレートが低下し、リアルタイム性が悪くなるという問題があった。
また、特許文献1に開示の技術では、生体組織内の音速分布にムラがあると、各々の仮想点音源からの波面にずれが生じるため、精度が悪くなってしまうという問題があった。
また、特許文献1に開示の技術では、生体組織内の音速分布にムラがあると、各々の仮想点音源からの波面にずれが生じるため、精度が悪くなってしまうという問題があった。
本発明の目的は、上記従来技術の問題点を解消し、超音波プローブの複数の素子からの超音波ビームを送信する際に、異なる2以上の素子からの送信で取得した素子データを重ねるマルチライン処理を行うことで、送信超音波ビームの広がりによる影響を低減し、SN比を上げ、解像度を上げることができ、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高い解像度で、最適な空間分解能を持つシャープな超音波画像を得ることができると共に、マルチライン処理の計算を含めた上で音速補正を行うことにより、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができる超音波検査装置、超音波検査方法、プログラム及び記録媒体を提供することにある。
本発明者は、本出願人の先の出願に係る特願2012-158009号明細書において、超音波プローブの複数の素子からの超音波ビームを送受信する際にマルチライン処理を行うことで画質の良い超音波画像を得ることができる超音波検査装置を提案している。上記目的を達成するために、更に、鋭意研究を重ねた結果、先の出願に開示の技術には、音速分布のムラを補正する音速補正について記載されているが、反射点(上記特許文献1の図6の観測点Px)からプローブまでの音速しか補正されておらず、各々の仮想的な点音源を形成する送信素子から反射点(Pxの点)までの音速は補正されていないため、音速の補正が不十分であったことを知見し、本発明に至ったものである。
すなわち、本発明の超音波検査装置は、超音波ビームを用いて検査対象物を検査する超音波検査装置であって、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子と、探触子に、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行わせる送信部と、個々の超音波ビームの送信に対応して、複数の素子が出力したアナログ素子信号を受け、所定の処理を施す受信部と、受信部が処理したアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとするAD変換部と、検査対象物内の音速値を決定する音速決定部と、検査対象物内の音速値を用いて2つ以上の第1の素子データから、2つ以上の第1の素子データのいずれかに対応する第2の素子データを生成する素子データ処理部と、を有し、音速決定部は、素子データ処理部において2つ以上の第1の素子データから第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として求めることを特徴とする。
ここで、音速決定部は、素子データ処理部を有し、素子データ処理部における2つ以上の第1の素子データの遅延時間補正計算に用いる音速値を最適化して最適な音速として求めることが好ましい。
また、音速決定部は、検査対象物内を複数の領域に分割したときの各領域における音速値を決定することが好ましい。
また、更に、2つ以上の第1の素子データを含む一部の第1の素子データ、あるいは全ての第1の素子データを保持する素子データ保持部を有することが好ましい。
また、音速決定部は、検査対象物内を複数の領域に分割したときの各領域における音速値を決定することが好ましい。
また、更に、2つ以上の第1の素子データを含む一部の第1の素子データ、あるいは全ての第1の素子データを保持する素子データ保持部を有することが好ましい。
また、素子データ保持部は、送信部により、検査対象物内における少なくとも2つ以上の重なり合う対象領域毎において、超音波ビームを送信し、受信部により、超音波ビームによって検査対象物内における少なくとも2つ以上の重なり合う対象領域毎に発生した超音波エコー信号を受信して生成された、素子データ毎に受信した各素子における受信時間毎の受信データを含む2以上の第1の素子データを保持し、素子データ処理部は、素子データ毎の超音波ビームを送信したときの素子の幾何学的な配置及び素子データ毎の受信データの各素子の幾何学的な配置に基づいて、2以上の第1の素子データから第2の素子データを生成することが好ましい。
また、送信部は、超音波ビームの送信を、中心となる素子を変更して、複数回、行うことが好ましい。
また、送信部は、超音波ビームの送信を、中心となる素子を変更して、複数回、行うことが好ましい。
また、素子データ処理部は、検査対象物内の音速値を用いて2つ以上の第1の素子データの各遅延時間を算出する遅延時間算出部と、2つ以上の第1の素子データを算出された各遅延時間及び受信された探触子の素子の位置に基づいて重ね合わせ、第2の素子データを生成する重ね合わせ処理部とを有することが好ましい。
また、更に、素子データ処理部で生成された第2の素子データに基づいて超音波画像を生成する画像生成部と、を有し、素子データ処理部は、画像生成部で生成されるべき超音波画像上の1ラインの第2の素子データを作成することが好ましい。
また、更に、素子データ処理部で生成された第2の素子データに基づいて超音波画像を生成する画像生成部と、を有し、素子データ処理部は、画像生成部で生成されるべき超音波画像上の1ラインの第2の素子データを作成することが好ましい。
また、音速決定部は、素子データ処理部と、素子データ処理部における遅延時間算出に用いられる検査対象物の音速値を変更する音速変更部と、音速変更部で変更された音速値を用いて第2の素子データを作成する素子データ処理部と、素子データ処理部により作成された第2の素子データを、音速変更部で変更された音速値を用いて整相加算して音線信号を生成する整相加算部と、整相加算部で生成された音線信号に基づいて超音波画像を生成する画像生成部と、画像生成部で生成された超音波画像の画質指標を算出する画質指標算出部と、画質指標算出部で算出された超音波画像の画質指標を判定する判定部と、判定部による画質指標の判定結果に応じて最適な音速値を設定する音速設定部と、を有することが好ましい。
また、音速決定部は、音速変更部が、音速値を初期音速値から最終音速値まで所定音速間隔で変更する度毎に、素子データ処理部が変更された当該音速値を用いて第2の素子データを作成し、整相加算部が同じ当該音速値を用いて第2の素子データを整相加算して音線信号を生成し、画像生成部が音線信号に基づいて超音波画像を生成し、画質指標算出部が超音波画像の画質指標を算出することを繰り返して、音速変更部によって変更された全ての複数の音速値にそれぞれ対応する複数の画質指標を算出しておき、判定部が、画質指標算出部で全ての複数の音速値毎に算出された複数の画質指標の中から最適な画質指標を判定し、音速設定部が、判定部で判定された最適な画質指標に対応する音速値を最適な音速値として設定することが好ましい。
さらに、音速決定部よりも前段に、第1の素子データを用いて最適な音速値を算出する第2の音速決定部を有し、音速決定部は、第2の音速決定部で算出した音速値を初期値として用いることが好ましい。
さらに、音速決定部よりも前段に、第1の素子データを用いて最適な音速値を算出する第2の音速決定部を有し、音速決定部は、第2の音速決定部で算出した音速値を初期値として用いることが好ましい。
また、上記目的を達成するために、本発明の超音波検査方法は、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子を用いて検査対象物を検査する超音波検査方法であって、探触子によって、複数の素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の素子がアナログ素子信号を出力し、このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとし、検査対象物内の音速値を用いて、2つ以上の第1の素子データから、第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として決定することを特徴とする。
ここで、最適な音速を決定するステップは、音速値を初期音速値から最終音速値まで所定音速間隔で変更する毎に、変更された当該音速値を用いて第2の素子データを作成し、同じ当該音速値を用いて第2の素子データを整相加算して音線信号を生成し、生成された音線信号に基づいて超音波画像を生成し、画質指標算出部が超音波画像の画質指標を算出することを繰り返して、変更された全ての複数の音速値毎にそれぞれ対応する複数の画質指標を算出しておき、全ての複数の音速値毎に算出された複数の画質指標の中から最適な画質指標を判定し、判定された最適な画質指標に対応する音速値を最適な音速値として設定することが好ましい。
また、上記目的を達成するために、本発明のプログラムは、超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して検査対象物を検査するために、複数の素子が配列された探触子によって、複数の素子を用い、所定の送信焦点を形成するように、超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の素子がアナログ素子信号を出力し、このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データを複数、取得するステップと、検査対象物内の音速値を用いて、得られた2つ以上の第1の素子データから第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる音速値を最適化して最適な音速値として決定するステップとを手順としてコンピュータに実行させることを特徴とする。
また、上記目的を達成するために、本発明の記録媒体は、上記記載のプログラムを記録したコンピュータ読み取り可能な記録媒体である。
また、上記目的を達成するために、本発明の記録媒体は、上記記載のプログラムを記録したコンピュータ読み取り可能な記録媒体である。
本発明によれば、超音波プローブの複数の素子からの超音波ビームを送信する際に、異なる2以上の素子からの送信で取得した素子データを重ねるマルチライン処理を行うことで、送信超音波ビームの広がりによる影響を低減し、SN比を上げ、解像度を上げることができ、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高い解像度で、最適な空間分解能を持つシャープな超音波画像を得ることができると共に、マルチライン処理の計算を含めた上で音速補正を行うことにより、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができる。
本発明に係る超音波検査装置、超音波検査方法、プログラム及び記録媒体を添付の図面に示す好適実施形態に基づいて以下に詳細に説明する。
図1は、本発明の超音波検査装置の構成の一実施例を概念的に示すブロック図である。
同図に示すように、超音波検査装置10は、超音波プローブ12と、超音波プローブ12に接続される送信部14及び受信部16と、A/D変換部18と、素子データ保持部20と、素子データ処理部22を備える音速補正部24と、画像生成部26と、表示制御部27と、表示部28と、制御部30と、操作部32と、格納部34とを有する。
同図に示すように、超音波検査装置10は、超音波プローブ12と、超音波プローブ12に接続される送信部14及び受信部16と、A/D変換部18と、素子データ保持部20と、素子データ処理部22を備える音速補正部24と、画像生成部26と、表示制御部27と、表示部28と、制御部30と、操作部32と、格納部34とを有する。
超音波プローブ12は、通常の超音波検査装置に用いられる探触子36を有する。
探触子36は、1次元又は2次元アレイ状に配列された複数の素子、即ち超音波トランスデューサを有している。これらの超音波トランスデューサは、検査対象物(以下、被検体という)の超音波画像の撮像の際に、それぞれ送信部14から供給される駆動信号に従って超音波ビームを被検体に送信すると共に、被検体からの超音波エコーを受信して受信信号を出力する。本実施形態では、探触子36の複数の超音波トランスデューサの内の一組を成す所定数の超音波トランスデューサの各々は、1つの超音波ビームの各成分を発生し、一組の所定数の超音波トランスデューサは、被検体に送信する1つの超音波ビームを発生する。
探触子36は、1次元又は2次元アレイ状に配列された複数の素子、即ち超音波トランスデューサを有している。これらの超音波トランスデューサは、検査対象物(以下、被検体という)の超音波画像の撮像の際に、それぞれ送信部14から供給される駆動信号に従って超音波ビームを被検体に送信すると共に、被検体からの超音波エコーを受信して受信信号を出力する。本実施形態では、探触子36の複数の超音波トランスデューサの内の一組を成す所定数の超音波トランスデューサの各々は、1つの超音波ビームの各成分を発生し、一組の所定数の超音波トランスデューサは、被検体に送信する1つの超音波ビームを発生する。
各超音波トランスデューサは、例えば、PZT(チタン酸ジルコン酸鉛)に代表される圧電セラミックや、PVDF(ポリフッ化ビニリデン)に代表される高分子圧電素子、PMN-PT(マグネシウムニオブ酸・チタン酸鉛固溶体)に代表される圧電単結晶等からなる圧電体の両端に電極を形成した素子、即ち振動子によって構成される。即ち、探触子36は、複数の超音波素子として複数の振動子が1次元又は2次元アレイ状に配列された振動子アレイということができる。
このような振動子の電極に、パルス状又は連続波状の電圧を印加すると、圧電体が伸縮し、それぞれの振動子からパルス状又は連続波状の超音波が発生して、それらの超音波の合成により超音波ビームが形成される。また、それぞれの振動子は、伝搬する超音波を受信することにより伸縮して電気信号を発生し、それらの電気信号は、超音波の受信信号として出力される。
送信部14は、例えば、複数のパルサを含んでおり、制御部30からの制御信号に応じて選択された送信遅延パターンに基づいて設定される音速または音速の分布に従い、探触子36の一組の所定数の超音波トランスデューサ(以下、超音波素子という)から送信される超音波ビーム成分が1つの超音波ビームを形成するようにそれぞれの駆動信号の遅延量を調節して組を成す複数の超音波素子に供給する。
受信部16は、制御部30からの制御信号に応じて、探触子36の各超音波素子によって超音波ビームと被検体との間の相互作用によって発生された超音波エコーを被検体から受信して、受信信号、即ち超音波素子毎のアナログ素子信号を増幅して出力し、増幅されたアナログ素子信号をA/D変換部18に供給する。
A/D変換部18は、受信部16に接続され、受信部16から供給されたアナログ素子信号を、デジタル素子データに変換する。A/D変換部18は、A/D変換されたデジタル素子データを素子データ保持部20に供給する。
A/D変換部18は、受信部16に接続され、受信部16から供給されたアナログ素子信号を、デジタル素子データに変換する。A/D変換部18は、A/D変換されたデジタル素子データを素子データ保持部20に供給する。
素子データ保持部20は、A/D変換部18から出力されるデジタル素子データを順次格納する。また、素子データ保持部20は、制御部30から入力されるフレームレートに関する情報(例えば、超音波の反射位置の深度、走査線の密度、視野幅を示すパラメータ)を上記のデジタル素子データ(以下、単に素子データという)に関連付けて格納する。
素子データ保持部20は、制御部30からの制御信号に応じて、2次元以上の位置座標上の対象領域において重なり合う少なくとも2つ以上の対象領域を検査する際に、この2つ以上の対象領域毎に受信部16で受信した超音波エコーから2つ以上の対象領域毎に生成された2以上の素子データを記憶保持するものであり、素子データ保持部20に記憶保持される素子データは、素子データ毎に受信した各素子における受信時間毎の受信データを含む、2つ以上の素子データである。即ち、素子データ保持部20は、2つ以上の素子データを含む一部の素子データ、あるいは全ての第1の素子データを保持する。
素子データ保持部20は、制御部30からの制御信号に応じて、2次元以上の位置座標上の対象領域において重なり合う少なくとも2つ以上の対象領域を検査する際に、この2つ以上の対象領域毎に受信部16で受信した超音波エコーから2つ以上の対象領域毎に生成された2以上の素子データを記憶保持するものであり、素子データ保持部20に記憶保持される素子データは、素子データ毎に受信した各素子における受信時間毎の受信データを含む、2つ以上の素子データである。即ち、素子データ保持部20は、2つ以上の素子データを含む一部の素子データ、あるいは全ての第1の素子データを保持する。
また、素子データ処理部22は、音速補正部24の内に設けられ、制御部30による制御に基づいて、被検体の検査対象領域内の所定の計算領域の計算座標(所定領域内のサンプリング点や観測点の集まり)について、その計算座標の音速(以下、環境音速ともいう)を用いてマルチライン処理を行って、即ち、素子データ保持部20によって記憶保持された2つ以上の対象領域毎に生成された2以上の素子データ(以下、未処理素子データという)を受信時間上で時間を合わせて重ね合わせて重ね合わせ処理後素子データ(以下、処理済素子データという)を生成するものである。なお、本発明においては、計算座標の音速又は環境音速とは、探触子36の超音波素子から送信された超音波ビームが被検体の検査対象領域内の所定の計算座標までの送信経路及び計算座標で反射した超音波エコーが探触子36の超音波素子に戻るまでの受信経路における平均音速、即ち、探触子36の超音波素子と計算座標との間の平均音速を言う。
音速補正部24は、素子データ保持部20と画像生成部26の整相加算部38との間に設けられ、内部にマルチライン処理を行う素子データ処理部22を備えるもので、制御部30による制御に基づいて、被検体の検査対象領域内の所定の計算領域の計算座標毎に、素子データ処理部22によるマルチライン処理を行った処理済データに基づいて、整相加算部38で用いられる送受信の最適な音速を求め、被検体の検査対象領域内の環境音速、したがって、被検体の検査対象領域の超音波画像に最も適した環境音速を求めると共に、素子データ処理部22によって最適な環境音速を用いるマルチライン処理が行われた処理済データを得るものである。
なお、素子データ処理部22及び音速補正部24の詳細については、後述する。
なお、素子データ処理部22及び音速補正部24の詳細については、後述する。
画像生成部26は、制御部30による制御下で、素子データ処理部22で生成され、音速補正部24から供給された処理済素子データから音線信号(受信データ)を生成し、この音線信号から超音波画像を生成するものである。
画像生成部26は、整相加算部38、検波処理部40、DSC42、画像作成部44、および、画像メモリ46を有する。
画像生成部26は、整相加算部38、検波処理部40、DSC42、画像作成部44、および、画像メモリ46を有する。
整相加算部38は、制御部30において設定された受信方向に応じて、音速補正部24で求められた最適な環境音速またはその分布に従い、素子データ処理部22で生成された処理済素子データにそれぞれの遅延を与えて加算することにより、受信フォーカス処理を行う。この受信フォーカス処理により、超音波エコーの焦点が絞り込まれた受信データ(音線信号)が生成される。なお、音速補正部24において最適な環境音速またはその分布が求まっていない場合には、予め記憶されている複数の受信遅延パターンの中から1つの受信遅延パターンを選択し、選択された受信遅延パターンに基づいて設定される環境音速またはその分布を用いても良い。
整相加算部38は、受信データを検波処理部40に供給する。
整相加算部38は、受信データを検波処理部40に供給する。
検波処理部40は、整相加算部38で生成された受信データに対し、超音波の反射位置の深度に応じて距離による減衰の補正を施した後、包絡線検波処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像データを生成する。
DSC(digital scan converter)42は、検波処理部40で生成されたBモード画像データを通常のテレビジョン信号の走査方式に従う画像データに変換(ラスター変換)する。
DSC(digital scan converter)42は、検波処理部40で生成されたBモード画像データを通常のテレビジョン信号の走査方式に従う画像データに変換(ラスター変換)する。
画像作成部44は、DSC42から入力されるBモード画像データに階調処理等の各種の必要な画像処理を施して検査や表示に供するためのBモード画像データを作成した後、作成された検査用又は表示用Bモード画像データを表示のために表示制御部27に出力する、或いは画像メモリ46に格納する。
画像メモリ46は、画像作成部44で作成された検査用Bモード画像データを一旦格納する。画像メモリ46に格納された検査用Bモード画像データは、必要に応じて、表示部28で表示するために表示制御部27に読み出される。
画像メモリ46は、画像作成部44で作成された検査用Bモード画像データを一旦格納する。画像メモリ46に格納された検査用Bモード画像データは、必要に応じて、表示部28で表示するために表示制御部27に読み出される。
表示制御部27は、画像作成部44によって画像処理が施された検査用Bモード画像信号に基づいて、表示部28に超音波画像を表示させる。
表示部28は、例えば、LCD等のディスプレイ装置を含んでおり、表示制御部27の制御の下で、超音波画像を表示する。
表示部28は、例えば、LCD等のディスプレイ装置を含んでおり、表示制御部27の制御の下で、超音波画像を表示する。
制御部30は、操作者により操作部32から入力された指令に基づいて超音波検査装置10の各部の制御を行う。
ここで、制御部30は、種々の情報、特に、音速補正部24で用いられる最適環境音速の算出に必要な情報、素子データ処理部22及び画像生成部26の整相加算部38で最適環境音速と共に用いられる遅延時間算出に必要な情報及び素子データ処理部22における素子データ処理に必要な情報が、操作者によって操作部32を介して入力された際に、かつ/又は格納部34から読み出された際に、操作部32から入力された、かつ/又は格納部34から読み出された上述の種々の情報を、必要に応じて、送信部14、受信部16、素子データ保持部20、素子データ処理部22、音速補正部24、画像生成部26及び表示制御部27等の各部に供給する。
ここで、制御部30は、種々の情報、特に、音速補正部24で用いられる最適環境音速の算出に必要な情報、素子データ処理部22及び画像生成部26の整相加算部38で最適環境音速と共に用いられる遅延時間算出に必要な情報及び素子データ処理部22における素子データ処理に必要な情報が、操作者によって操作部32を介して入力された際に、かつ/又は格納部34から読み出された際に、操作部32から入力された、かつ/又は格納部34から読み出された上述の種々の情報を、必要に応じて、送信部14、受信部16、素子データ保持部20、素子データ処理部22、音速補正部24、画像生成部26及び表示制御部27等の各部に供給する。
操作部32は、操作者が入力操作を行うためのものであり、キーボード、マウス、トラックボール、タッチパネル等から形成することができる。
また、操作部32は、操作者が必要に応じて各種の情報、特に上述の遅延時間算出に用いられるプローブ12の探触子36の複数の超音波素子、被検体の検査対象領域、その計算領域(計算座標X1~Xend)、音速補正(最適環境音速決定)のための音速探索範囲(初期音速値Vst、最終音速値Vend、音速値変更幅(音速間隔)ΔV)、超音波ビームの焦点位置、探触子36の送信開口及び受信開口等に関する情報、並びに重ね合わせ素子データ数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報等を入力操作するための入力装置を備えている。
また、操作部32は、操作者が必要に応じて各種の情報、特に上述の遅延時間算出に用いられるプローブ12の探触子36の複数の超音波素子、被検体の検査対象領域、その計算領域(計算座標X1~Xend)、音速補正(最適環境音速決定)のための音速探索範囲(初期音速値Vst、最終音速値Vend、音速値変更幅(音速間隔)ΔV)、超音波ビームの焦点位置、探触子36の送信開口及び受信開口等に関する情報、並びに重ね合わせ素子データ数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報等を入力操作するための入力装置を備えている。
格納部34は、操作部32から入力された各種の情報、特に、上述のプローブ12、計算領域、音速探索範囲、焦点位置、送信開口及び受信開口等に関する情報、重ね合わせ素子データ数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報等、及び送信部14、受信部16、素子データ保持部20、素子データ処理部22、音速補正部24、画像生成部26及び表示制御部27等の制御部30で制御される各部の処理や動作に必要な情報、並びに、各部の処理や動作を実行させるための動作プログラムや処理プログラム等を格納するもので、ハードディスク、フレキシブルディスク、MO、MT、RAM、CD-ROM、DVD-ROM等の記録媒体を用いることができる。
なお、素子データ処理部22、音速補正部24、画像生成部26の整相加算部38、検波処理部40、DSC42並びに画像作成部44、及び表示制御部27は、CPUと、CPUに各種の処理を行わせるための動作プログラムから構成されるが、それらをデジタル回路で構成してもよい。
なお、素子データ処理部22、音速補正部24、画像生成部26の整相加算部38、検波処理部40、DSC42並びに画像作成部44、及び表示制御部27は、CPUと、CPUに各種の処理を行わせるための動作プログラムから構成されるが、それらをデジタル回路で構成してもよい。
次に、図1に示す素子データ処理部22を備える音速補正部24について説明する。
図2は、図1に示す本発明の超音波検査装置の音速補正部の一例の詳細含む超音波検査装置の要部を処理フローに沿って示すブロック図である。
音速補正部24は、設定音速Vを検査対象領域内の音速として推定される所定音速範囲全域を所定音速値変更幅(音速間隔)で変化させて、例えば、人体等の生体を対象とする場合、具体的には、1400m/sから1650m/s又は1700m/sまで、1m/s~50m/sの所定刻み(ステップ)、例えば、10m/s刻みで変化させて、検査対象領域内の画像の輝度値、鮮鋭度等の画質指標(画質評価指標)や音速の収束度合い等に基づいて網羅的に最適音速を求めるものである。
図2は、図1に示す本発明の超音波検査装置の音速補正部の一例の詳細含む超音波検査装置の要部を処理フローに沿って示すブロック図である。
音速補正部24は、設定音速Vを検査対象領域内の音速として推定される所定音速範囲全域を所定音速値変更幅(音速間隔)で変化させて、例えば、人体等の生体を対象とする場合、具体的には、1400m/sから1650m/s又は1700m/sまで、1m/s~50m/sの所定刻み(ステップ)、例えば、10m/s刻みで変化させて、検査対象領域内の画像の輝度値、鮮鋭度等の画質指標(画質評価指標)や音速の収束度合い等に基づいて網羅的に最適音速を求めるものである。
このような音速補正部24による音速補正方法は、例えば、特開2011-92686号公報に開示の音速補正方法を用いることができる。なお、この公報に開示の音速補正方法では、着目領域毎に送信フォーカスを結ぶ必要がある。しかしながら、本発明では、音速補正部24内の素子データ処理部22の重ね合わせ処理部の重ね合わせ処理によって擬似的に如何なる深さにおいても送信フォーカスを結ぶことができている処理済素子データが得られているので、計算領域や着目領域毎に送信フォーカスを設定する必要はない。ただし、本発明においても、この公報に開示の音速補正方法の通りに、計算領域や着目領域毎に送信フォーカスを設定するようにしても良いのは勿論である。
図2に示す実施形態においては、音速補正部24は、計算領域設定部70と、初期計算座標設定部72と、最適音速決定部23と、最終計算座標判定部88と、計算座標更新部90とを有し、最適音速決定部23は、音速探索範囲設定部74と、初期音速値設定部76と、整相加算部38aと、検波処理部40aと、画像作成部44aと、画質指標算出部78と、最終音速値判定部80と、音速値更新部82と、最適音速判定部84と、音速設定部86と、を有する。
音速補正部24では、図2に示すように、検査対象領域内の計算座標X1~Xendにおける設定音速Vを変化させ、それぞれの設定速度Vにおいて設定速度Vに基づいて、素子データ保持部20から供給された未処理素子データを素子データ処理部22でマルチライン処理した処理済素子データを用いて整相加算部38aにおいて受信フォーカス処理をして整相加算を行い音線信号を生成し、この音線信号から画像作成部44aにおいて超音波画像を形成し、それぞれの設定音速Vでの検査対象領域内の超音波画像の輝度値、鮮鋭度等の画質指標を算出する。
音速補正部24では、各設定音速Vでの超音波画像の輝度値、鮮鋭度等の画質指標を最適音速判定部84において比較して、例えば、得られた輝度値、鮮鋭度等の画質指標の値が最も高い設定音速Vを判定し、判定された設定音速Vを音速設定部86において最適な音速値(最適音速値)として設定する。
音速補正部24では、各設定音速Vでの超音波画像の輝度値、鮮鋭度等の画質指標を最適音速判定部84において比較して、例えば、得られた輝度値、鮮鋭度等の画質指標の値が最も高い設定音速Vを判定し、判定された設定音速Vを音速設定部86において最適な音速値(最適音速値)として設定する。
なお、音速補正部24においては、音速値の変更及び最適音速の設定は、図2に示すように網羅的に行っても良いし、探索的に行っても良い。
また、音速補正部24における音速値の変更及び最適音速の設定の際に、音速値を設定する領域や計算座標は、特に制限的ではなく、被写体全体や検査対象領域全体を同一の音速値として設定してもよいが、被写体や検査対象領域を微小領域ごとに区切って、微小領域ごとに音速値を設定してもよく、画素単位で1点ずつ音速値を設定しても良い。
また、音速補正部24における音速値の変更及び最適音速の設定の際に、音速値を設定する領域や計算座標は、特に制限的ではなく、被写体全体や検査対象領域全体を同一の音速値として設定してもよいが、被写体や検査対象領域を微小領域ごとに区切って、微小領域ごとに音速値を設定してもよく、画素単位で1点ずつ音速値を設定しても良い。
まず、図2に示す音速補正部24の素子データ処理部22について、図3に基づいて詳細に説明する。
同図に示すように、素子データ処理部22は、マルチライン処理を行うためのもので、遅延時間算出部48と、重ね合わせ処理部50とを有する。
なお、ここでは、素子データ処理部22でマルチライン処理を行う被検体の検査対象領域内の計算座標は、予め、初期計算座標設定部72によって初期計算座標X1に、又は計算座標更新部90によって更新された計算座標Xi(i=2~end:以下、計算座標Xi(i=1~end)で代表する)に設定されているものとし、その音速値は、初期音速値設定部76によって初期音速値Vstに、又は音速値更新部82によって更新された音速値V(以下、音速値Xで代表する)に予め設定されているものとする。また、素子データ処理部22におけるマルチライン処理に必要である、計算座標Xiに関わる未処理素子データは、素子データ保持部20から読み出されて、供給されているものとする。
同図に示すように、素子データ処理部22は、マルチライン処理を行うためのもので、遅延時間算出部48と、重ね合わせ処理部50とを有する。
なお、ここでは、素子データ処理部22でマルチライン処理を行う被検体の検査対象領域内の計算座標は、予め、初期計算座標設定部72によって初期計算座標X1に、又は計算座標更新部90によって更新された計算座標Xi(i=2~end:以下、計算座標Xi(i=1~end)で代表する)に設定されているものとし、その音速値は、初期音速値設定部76によって初期音速値Vstに、又は音速値更新部82によって更新された音速値V(以下、音速値Xで代表する)に予め設定されているものとする。また、素子データ処理部22におけるマルチライン処理に必要である、計算座標Xiに関わる未処理素子データは、素子データ保持部20から読み出されて、供給されているものとする。
遅延時間算出部48は、操作部32から入力された、もしくは、操作部32から入力されて格納部34に格納されているプローブ12の探触子36の複数の超音波素子、超音波ビームの焦点位置、探触子36の送信開口及び受信開口等に関する情報を事前に取得しておき、被検体の検査対象領域内の計算座標Xiの音速は、初期音速値設定部76又は音速値更新部82によって予め設定された音速値V(Vst)を用いて、超音波ビームを形成し送信する、送信開口の超音波素子(送信素子)と、被検体からの、超音波ビームによる超音波エコーを受信する、受信開口の超音波素子との幾何学的配置に基づいて、受信開口の各超音波素子で受信される素子データの遅延時間を算出する。
重ね合わせ処理部50は、操作部32から入力された、もしくは、操作部32から入力されて格納部34に格納されている、重ね合わせる素子データの数及び重ね合わせ処理方法等の素子データのマルチライン処理に関する情報に基づいて、素子データ保持部20によって記憶保持された2以上の対象領域毎に生成された2以上の未処理素子データを読み出し、遅延時間算出部48でそれぞれ算出された遅延時間に基づいて、2以上の未処理素子データを受信時間上で、即ち時間を合わせて、かつ、受信された探触子の素子の絶対的な位置を合わせて、重ね合わせて処理済素子データを生成する。
次に、素子データ処理部22で行う素子データのマルチライン処理について詳細に説明する。
まず、超音波プローブ12の探触子36の送信用超音波素子(以下、単に送信素子という)から超音波ビーム(以下、送信ビームという)を被検体に送信し、被検体との間の相互作用によって発生された超音波エコーを探触子36の受信用超音波素子(以下、単に受信素子という)で受信して素子データを得る場合において、送信素子からの送信ビームと受信素子で得られる素子データとの関係について説明する。
まず、超音波プローブ12の探触子36の送信用超音波素子(以下、単に送信素子という)から超音波ビーム(以下、送信ビームという)を被検体に送信し、被検体との間の相互作用によって発生された超音波エコーを探触子36の受信用超音波素子(以下、単に受信素子という)で受信して素子データを得る場合において、送信素子からの送信ビームと受信素子で得られる素子データとの関係について説明する。
図4(a)及び(c)に示すように、それぞれ7つの超音波素子(以下、単に、素子ともいう)52a~52g及び52b~52hを受信素子として超音波エコーを受信して素子データを取得する時、反射点54を含む検査対象領域に送信する送信ビーム56が理想的に素子間隔以下に絞れている理想的な場合には、図4(a)のように、検査対象領域内の反射点54の真上にある、素子52a~52gの中心の素子52dを送信素子として送信ビーム56を送信し、受信素子52a~52gで超音波エコーを受信して素子データを取得する場合、送信ビーム56の焦点58は、素子54dと反射点54とを結ぶ一直線上にあり、送信ビーム56は、反射点54まで送信されるので、反射点54から反射される超音波エコーが生成される。反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子52a~52gに受信され、受信素子52a~52gによって、図4(b)に示すような素子データ62が得られることになる。
これに対し、図4(c)に示すように、送信素子の中心を、反射点54に対して1素子分、素子の方向(図中右方向)にシフトさせ、反射点54の真上にある素子52dに隣接する素子52eを送信素子として送信ビーム56を送信し、受信素子52b~52hで超音波エコーを受信する場合には、送信ビーム56の送信方向、即ち、送信素子52eと焦点58とを結ぶ直線上に反射点54が存在しないため、送信ビーム56は、反射点54に送信されない。このため、反射点54から反射される超音波エコーは生成されず、受信素子52b~52hでは、超音波エコーを受信しないので、図4(d)に示すように、素子データが得られないことになる。
しかしながら、図5(a)及び(c)に示すように、実際の送信ビーム64は、素子間隔より幅が広い。
ここで、図5(a)のように、反射点54の真上にある素子52dを送信素子として送信ビーム64を送信した場合には、図4(a)の場合と同様に、送信ビーム64が幅広であっても、その焦点58は、素子54dと反射点54とを結ぶ一直線上にあり、送信ビーム64は、反射点54で反射され、超音波エコーが生成される。その結果、図4(a)の場合と同様に、反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子52a~52gに受信され、受信素子52a~52gによって、図5(b)に示すような真の素子データ66が得られることになる。
ここで、図5(a)のように、反射点54の真上にある素子52dを送信素子として送信ビーム64を送信した場合には、図4(a)の場合と同様に、送信ビーム64が幅広であっても、その焦点58は、素子54dと反射点54とを結ぶ一直線上にあり、送信ビーム64は、反射点54で反射され、超音波エコーが生成される。その結果、図4(a)の場合と同様に、反射点54からの超音波エコーは、所定角度に拡がる受信経路60を通って受信素子52a~52gに受信され、受信素子52a~52gによって、図5(b)に示すような真の素子データ66が得られることになる。
一方、図5(c)に示すように、図4(c)の場合と同様に、送信素子の中心を、反射点54に対して1素子分、素子の方向(図中右方向)にシフトさせ、反射点54の真上にある素子52dに隣接する素子52eを送信素子として送信ビーム64を送信し、受信素子52b~52hで超音波エコーを受信する場合、送信ビーム64は幅広であるため、その送信方向、即ち、送信素子52eと焦点58とを結ぶ直線上に反射点54が存在していなくても、送信ビーム64は、反射点54に送信されることになる。このため、反射点54から本来存在しない超音波エコー、所謂ゴーストの反射信号が発生し、反射点54からのゴーストの反射信号は、所定角度に拡がる受信経路60を通って受信素子52b~52hに受信され、受信素子52b~52hによって、図5(d)に示すようなゴーストの素子データ68が得られることになる。
このようなゴーストの素子データ68は、素子データから生成される超音波画像の精度を低下させる原因となる。
このようなゴーストの素子データ68は、素子データから生成される超音波画像の精度を低下させる原因となる。
ここで、図5(c)に示す送信ビーム64が送信素子52eから焦点58を経由して反射点54に至る送信経路とゴーストの反射信号が反射点54から各受信素子52b~52hに至る受信経路との和(伝播距離)は、それぞれ図5(a)に示す送信ビーム64が送信素子52dから焦点58を経由して反射点54に至る送信経路と真の反射超音波エコーが反射点54から各受信素子52a~52gに至る受信経路との和(伝播距離)より長くなるため、図5(d)に示すようなゴーストの素子データ68は、図5(b)に示すような真の素子データ66に対して遅延することになる。
本発明の素子データ処理部22の遅延時間算出部48においては、真の素子データに対するゴーストの素子データの時間差、即ち遅延時間は、送信素子、超音波ビームの焦点、被検体の反射点、及び受信素子の幾何学的配置から算出される。したがって、遅延時間の計算には、超音波プローブ12の形状(素子間隔、リニア、コンベックスなど)、焦点位置、送信開口、受信開口などの情報と、被検体の検査対象領域の音速値とが必要であり、遅延時間算出部48では、操作部32によって入力された、若しくは格納部34に格納されたこれらの情報と、音速値更新部82等で予め設定された音速値Vを取得して遅延時間の計算を行う。遅延時間は、音速値Vを用いて、例えば送信素子、超音波ビームの焦点、被検体の反射点、及び受信素子の幾何学的配置から算出される、送信素子から焦点を経て反射点に至る送信ビームの送信経路及び反射点から受信素子に至る真の反射超音波エコー又はゴーストの反射信号の受信経路の合計長さ(伝播距離)とその音速によって算出される伝播時間の差から算出することができる。
本発明では、例えば、図6(a)及び図6(b)に示すようにして、真の反射超音波エコーとゴーストの反射信号の場合の送信ビームの送信経路及び受信経路の長さを求めることができる。
真の反射超音波エコーの場合、図6(a)に示すように、送信素子52dと受信素子52d(受信素子52a~52gの中心)とが一致し、その真下に、焦点58及び反射点54が配置されており、反射点54の真上の素子52dの位置をxy2次元座標上の座標(x0、0)とし、素子間隔をLe、焦点58の位置を座標(x0、df)、反射点54の位置を座標(x0、z)とする時、送信素子52dの位置も反射点54の真上の素子52dと同じく座標(x0、0)となり、送信素子52dから焦点58を経て反射点54に至る送信ビームの送信経路61の長さ(送信経路距離)Ltaは、及び、反射点54から受信素子52dに至る真の反射超音波エコーの受信経路60の長さ(受信経路距離)Lraは、Lta=Lra=zによって算出することができる。
したがって、真の反射超音波エコーの場合の超音波の伝播距離Luaは、Lua=Lta+Lra=2zとなる。
真の反射超音波エコーの場合、図6(a)に示すように、送信素子52dと受信素子52d(受信素子52a~52gの中心)とが一致し、その真下に、焦点58及び反射点54が配置されており、反射点54の真上の素子52dの位置をxy2次元座標上の座標(x0、0)とし、素子間隔をLe、焦点58の位置を座標(x0、df)、反射点54の位置を座標(x0、z)とする時、送信素子52dの位置も反射点54の真上の素子52dと同じく座標(x0、0)となり、送信素子52dから焦点58を経て反射点54に至る送信ビームの送信経路61の長さ(送信経路距離)Ltaは、及び、反射点54から受信素子52dに至る真の反射超音波エコーの受信経路60の長さ(受信経路距離)Lraは、Lta=Lra=zによって算出することができる。
したがって、真の反射超音波エコーの場合の超音波の伝播距離Luaは、Lua=Lta+Lra=2zとなる。
ゴーストの反射信号の場合、図6(b)に示すように、図6(a)の場合と比べて、送信素子52eの位置が反射点54に対して1素子分横(x方向:図中右方向)にずれて、焦点58は送信素子52eの真下に配置されるが、反射点54は受信素子52dの真下に配置されており、反射点54の真上の受信素子52dの位置を、図6(a)の場合と同じくxy2次元座標上の座標(x0、0)とし、素子間隔をLe、反射点54の位置を座標(x0、z)とすると、送信素子52eの位置は座標(x0+Le、0)、焦点58の位置は座標(x0+Le、df)となるので、送信素子52eから焦点58を経て反射点54に至る送信ビームの送信経路61の長さ(送信経路距離)Ltbは、Ltb=df+√{(z-df)2+Le2}によって算出することができ、反射点54から受信素子52dに至るゴーストの反射信号の受信経路60の長さ(受信経路距離)Lrbは、Lrb=zによって算出することができる。
したがって、ゴーストの反射信号の場合の超音波の伝播距離Lubは、Lub=Ltb+Lrb=df+√{(z-df)2+Le2}+zとなる。
したがって、ゴーストの反射信号の場合の超音波の伝播距離Lubは、Lub=Ltb+Lrb=df+√{(z-df)2+Le2}+zとなる。
こうして、図6(a)に示す幾何学配置で求めた送信経路61の距離Ltaと受信経路60の距離Lraを合計した超音波の伝播距離Luaを音速で割った値が真の反射超音波エコーの伝播時間となり、図6(b)に示す幾何学配置で求めた送信経路61の距離Ltbと受信経路60の距離Lrbを合計した超音波の伝播距離Lubを音速で割った値がゴーストの反射信号の伝播時間となる。遅延時間の算出は、反射点54と送信素子52(52d)のx座標が一致している時の真の超音波エコーの伝播時間から、反射点54と送信素子52(52e)のx座標を1素子間隔ずつずらしたときのゴーストの反射信号の伝播時間の差から遅延時間を求める。
なお、図6(a)及び図6(b)の幾何学モデルでは、送信経路61が焦点58を経由したモデルになっているが、本発明はこれに限定されず、例えば、焦点58を経由せずに直接反射点54に至る経路であっても良い。
なお、図6(a)及び図6(b)の幾何学モデルでは、送信経路61が焦点58を経由したモデルになっているが、本発明はこれに限定されず、例えば、焦点58を経由せずに直接反射点54に至る経路であっても良い。
また、図6(a)及び図6(b)の幾何学モデルはリニアプローブの場合であるが、これに限らず他のプローブにおいても、プローブの形状から同様の幾何学計算を行うことができる。例えば、コンベックスプローブの場合、プローブの半径と素子間隔の角度から幾何学モデルを設定して同じように計算することができる。
また、ステア送信の場合には、送信角度などの情報を考慮した幾何学モデル(図示せず)を用い、送信素子と反射点との位置関係から真の超音波エコーの素子データ及びその周辺のゴーストの素子データの遅延時間を算出することができる。
さらに、幾何学モデルよって遅延時間を算出する方法に限らず、あらかじめ装置の計測条件に合わせて高輝度反射点を計測した計測結果から、計測条件毎に遅延時間を求めておき、その遅延時間を装置内に記憶しておくことで、同じ計測条件の遅延時間を読み出すようにしておいてもよい。
また、ステア送信の場合には、送信角度などの情報を考慮した幾何学モデル(図示せず)を用い、送信素子と反射点との位置関係から真の超音波エコーの素子データ及びその周辺のゴーストの素子データの遅延時間を算出することができる。
さらに、幾何学モデルよって遅延時間を算出する方法に限らず、あらかじめ装置の計測条件に合わせて高輝度反射点を計測した計測結果から、計測条件毎に遅延時間を求めておき、その遅延時間を装置内に記憶しておくことで、同じ計測条件の遅延時間を読み出すようにしておいてもよい。
図6(c)に、中央に真の信号の素子データ66及び真の信号の周辺のゴーストの素子データ68を示し、図6(d)に、上述の幾何学計算から得られた素子データ66に対するゴーストの素子データ68の遅延時間の一例を示す。真の信号の素子データ66を中心にゴーストの信号の素子データ68は、対称的に時間が遅れることが示されている。
なお、こうして、素子データ処理部22の遅延時間算出部48において算出された遅延時間を整相加算部38における遅延補正に用いることもできる。
なお、こうして、素子データ処理部22の遅延時間算出部48において算出された遅延時間を整相加算部38における遅延補正に用いることもできる。
次に、本発明の素子データ処理部22の重ね合わせ処理部50においては、こうして遅延時間算出部48において算出された遅延時間を用いて真の超音波エコーの素子データ及びその周辺のゴーストの素子データの重ね合わせ処理を行う。
重ね合わせ処理部50における重ね合わせ処理では、重ね合わせる時の重ね合わせ素子データ数と重ね合わせ処理方法の情報が必要になるが、これらは、予め、操作部32によって入力しておいても良いし、格納部34に格納しておいても良い。
なお、重ね合わせ処理部50において重ね合わせる時の重ね合わせ素子データの範囲は、計算領域設定部70で設定される音速補正計算の対象とする領域内にとどまる範囲にある素子データであることが望ましいが、計算領域設定部70で設定される領域によって制限されるものではない。例えば、重ね合わせる素子データが隣接する領域同士でオーバーラップするような範囲で重ね合わせてもよい。
重ね合わせ処理部50における重ね合わせ処理では、重ね合わせる時の重ね合わせ素子データ数と重ね合わせ処理方法の情報が必要になるが、これらは、予め、操作部32によって入力しておいても良いし、格納部34に格納しておいても良い。
なお、重ね合わせ処理部50において重ね合わせる時の重ね合わせ素子データの範囲は、計算領域設定部70で設定される音速補正計算の対象とする領域内にとどまる範囲にある素子データであることが望ましいが、計算領域設定部70で設定される領域によって制限されるものではない。例えば、重ね合わせる素子データが隣接する領域同士でオーバーラップするような範囲で重ね合わせてもよい。
図7(a)~(h)に、重ね合わせ処理部50で行われる、素子データ数が5つ、重ね合わせ素子データ数が3つの場合の重ね合わせ処理の一具体例を示す。
図7(a)は、5つの素子データを横に並べて表示しており、素子データ毎に、超音波ビームを送信し、反射信号を受信した様子を表している。各素子データの横軸は、受信素子を表しており、それぞれの素子データにおいて超音波ビームの送信時における中心の素子を中心にして表示している。縦軸は、受信時間を表す。
5つの素子データのうち、真中の素子データでは、素子データの中心の素子(受信素子の中心の素子)、即ち、送信時における中心の素子(送信素子)の真下に反射点が存在しており、反射点からの反射信号が受信されている。つまり、この反射信号は真の信号であり、真中の素子データは、真の信号を表す。
図7(a)は、5つの素子データを横に並べて表示しており、素子データ毎に、超音波ビームを送信し、反射信号を受信した様子を表している。各素子データの横軸は、受信素子を表しており、それぞれの素子データにおいて超音波ビームの送信時における中心の素子を中心にして表示している。縦軸は、受信時間を表す。
5つの素子データのうち、真中の素子データでは、素子データの中心の素子(受信素子の中心の素子)、即ち、送信時における中心の素子(送信素子)の真下に反射点が存在しており、反射点からの反射信号が受信されている。つまり、この反射信号は真の信号であり、真中の素子データは、真の信号を表す。
真中の素子データ以外の両側2つの素子データについては、送信時における中心の素子の真下には反射点は存在していないが、送信した超音波ビームの広がりによって、真中の素子データの送信素子の真下に存在する反射点に超音波ビームが当たることで生じた反射信号、即ちゴーストが写り込んでいる。ゴーストは、真の信号から離れるほど反射点までの超音波の伝播時間が長くなるため、真の信号よりも受信時間が遅くなる。また、反射点からの反射信号が初めに受信される受信素子の位置は、反射点の真上の素子であるが、素子データの横軸は超音波ビームの送信時における中心の素子を中心にしているため、素子データ毎にこの中心の素子を1素子ずつずらして送信していることから、各素子データにおいて素子の絶対位置は1素子ずつずれている。つまり、真中の素子データでは、反射点からの反射信号がはじめに受信される受信素子は真中の素子であるが、両隣の素子データにおいては、真中の素子データよりも1素子ずれており、右側の素子データでは左に1素子ずれ、左側の素子データでは右に1素子ずれている。更に、両端の素子データでは、真中の素子データよりも2素子ずれており、右端の素子データでは左に2素子ずれ、左端の素子データでは右に2素子ずれている。このように、ゴーストの信号は、真の信号に対して、受信時間が遅れるだけでなく、受信素子の方向に対してもずれを生じている。
図7(b)に、図7(a)に示す5つの素子データのうちの真中の素子データに対する受信時間の遅延時間の一例を示す。
重ね合わせ処理部50では、図7(b)に示す遅延時間を用いて、真中の素子データを注目素子の素子データとした場合に、注目素子の素子データを中心に、重ね合わせ素子データ数分、図示例では3素子データ分だけ遅延時間補正を行うと共に、各素子データを注目素子との素子位置の差(中心素子の位置の差)に応じて、図示例では両側に1素子分だけ横方向にシフトさせて、即ち位相を合わせて3素子データ分の未処理素子データを重ね合わせ、注目素子の素子データの1つの重ね合わせ処理済素子データとして求める。
重ね合わせ処理部50では、図7(b)に示す遅延時間を用いて、真中の素子データを注目素子の素子データとした場合に、注目素子の素子データを中心に、重ね合わせ素子データ数分、図示例では3素子データ分だけ遅延時間補正を行うと共に、各素子データを注目素子との素子位置の差(中心素子の位置の差)に応じて、図示例では両側に1素子分だけ横方向にシフトさせて、即ち位相を合わせて3素子データ分の未処理素子データを重ね合わせ、注目素子の素子データの1つの重ね合わせ処理済素子データとして求める。
こうして得られた注目素子の素子データの重ね合わせ処理済素子データを図7(c)に示す。
図7(a)に示す注目素子の素子データは、真の信号の素子データであることから、注目素子の素子データの両側の隣接素子データの未処理素子データに遅延時間補正及び横方向のシフトを行って位相合わせを行うと、図7(c)に示すように、隣接素子データの未処理素子データと、注目素子の未処理素子データとは、位相が合うので高輝度位置で重なり合う。したがって、これらの素子データを、例えば加算すると素子データ値は大きな値(高輝度値)を示し、例えば、平均して平均値を求めても強調された値(高輝度値)を示す。
図7(a)に示す注目素子の素子データは、真の信号の素子データであることから、注目素子の素子データの両側の隣接素子データの未処理素子データに遅延時間補正及び横方向のシフトを行って位相合わせを行うと、図7(c)に示すように、隣接素子データの未処理素子データと、注目素子の未処理素子データとは、位相が合うので高輝度位置で重なり合う。したがって、これらの素子データを、例えば加算すると素子データ値は大きな値(高輝度値)を示し、例えば、平均して平均値を求めても強調された値(高輝度値)を示す。
これに対し、図7(d)は、図7(a)と同じ素子データであるが、真中の素子データの左隣の素子データ、つまりゴーストを注目素子の素子データとした場合の一例を示す。
図7(e)は、図7(b)と同じものであり、図7(a)に示す5つの素子データのうちの真中の素子データに対する受信時間の遅延時間の一例を示す。即ち、図7(a)と図7(d)は同じ素子データであるので、図7(d)に示す5つの素子データの真中の素子データに対する受信時間の遅延時間とも同じである。
重ね合わせ処理部50では、図7(e)(即ち、図7(b)と同じ)に示す遅延時間を用いて、注目素子の素子データを中心に、重ね合わせ素子データ数分、図示例では3素子データ分だけ遅延時間補正を行うと共に、各素子データを注目素子との素子位置の差(中心素子の位置の差)に応じて、図示例では両側に1素子分だけ横方向にシフトさせて、3素子データ分の未処理素子データを重ね合わせ、注目素子の1つの重ね合わせ処理済素子データとして求める。
図7(e)は、図7(b)と同じものであり、図7(a)に示す5つの素子データのうちの真中の素子データに対する受信時間の遅延時間の一例を示す。即ち、図7(a)と図7(d)は同じ素子データであるので、図7(d)に示す5つの素子データの真中の素子データに対する受信時間の遅延時間とも同じである。
重ね合わせ処理部50では、図7(e)(即ち、図7(b)と同じ)に示す遅延時間を用いて、注目素子の素子データを中心に、重ね合わせ素子データ数分、図示例では3素子データ分だけ遅延時間補正を行うと共に、各素子データを注目素子との素子位置の差(中心素子の位置の差)に応じて、図示例では両側に1素子分だけ横方向にシフトさせて、3素子データ分の未処理素子データを重ね合わせ、注目素子の1つの重ね合わせ処理済素子データとして求める。
こうして得られた注目素子の素子データの重ね合わせ処理済素子データを図7(f)に示す。
図7(d)に示す注目素子の素子データは、ゴーストの素子データであることから、注目素子の素子データの両側の隣接素子データの未処理素子データに遅延時間補正及び横方向のシフトを行って位相合わせを行っても、図7(f)に示すように、隣接素子データの各未処理素子データと注目素子の未処理素子データとは、それぞれ位相が合わないので重なり合わない。このため、これらの3つの素子データを、例えば加算しても、位相が合っていないために、位相が反転している信号などは信号が打ち消しあうため、加算値は大きくならず、例えば、平均して平均値を求めると小さな値を示すことになる。
図7(d)に示す注目素子の素子データは、ゴーストの素子データであることから、注目素子の素子データの両側の隣接素子データの未処理素子データに遅延時間補正及び横方向のシフトを行って位相合わせを行っても、図7(f)に示すように、隣接素子データの各未処理素子データと注目素子の未処理素子データとは、それぞれ位相が合わないので重なり合わない。このため、これらの3つの素子データを、例えば加算しても、位相が合っていないために、位相が反転している信号などは信号が打ち消しあうため、加算値は大きくならず、例えば、平均して平均値を求めると小さな値を示すことになる。
他の素子データに関しても、注目素子の素子データとして同様の遅延時間補正及び横方向のシフトを行った結果、図示例の5素子データそれぞれについての隣接する3つの素子データの重なり状態を図7(g)に示し、これらに対して、重ね合わせ処理として、例えば、加算処理、若しくは平均処理した結果を図7(h)に示す。
図7(h)に示すように、図7(a)に示す送信素子と反射点との座標が一致している時の注目素子の素子データでは、真の信号の素子データが高輝度値を持つ重ね合わせ処理済素子データとして求められ、その両側の各2素子データの全4素子データでは、ゴーストの素子データは互いに位相が合わない素子データを加算し、又は平均するので、互いに打ち消し合うことになるため、ゴーストの重ね合わせ処理済素子データは、その値が真の信号の素子データである高輝度値を持つ重ね合わせ処理済素子データに対して小さくなり、真の信号の素子データに対してゴーストの素子データの影響を低減させることができ、又は、その影響を無視できる程、小さくすることができる。
そのため、処理済素子データに整相加算や検波処理を行って、受信データを生成して、超音波画像を生成することにより、同様に、ゴーストの影響を無くし、すなわち音線上の全ての点で焦点を結んだのに等しい素子データで超音波画像生成できるので、高輝度で、鮮鋭性に優れた、高画質な超音波画像を生成することができる。
なお、本発明においては、この処理済素子データの生成を、マルチライン処理と言う。
図7(h)に示すように、図7(a)に示す送信素子と反射点との座標が一致している時の注目素子の素子データでは、真の信号の素子データが高輝度値を持つ重ね合わせ処理済素子データとして求められ、その両側の各2素子データの全4素子データでは、ゴーストの素子データは互いに位相が合わない素子データを加算し、又は平均するので、互いに打ち消し合うことになるため、ゴーストの重ね合わせ処理済素子データは、その値が真の信号の素子データである高輝度値を持つ重ね合わせ処理済素子データに対して小さくなり、真の信号の素子データに対してゴーストの素子データの影響を低減させることができ、又は、その影響を無視できる程、小さくすることができる。
そのため、処理済素子データに整相加算や検波処理を行って、受信データを生成して、超音波画像を生成することにより、同様に、ゴーストの影響を無くし、すなわち音線上の全ての点で焦点を結んだのに等しい素子データで超音波画像生成できるので、高輝度で、鮮鋭性に優れた、高画質な超音波画像を生成することができる。
なお、本発明においては、この処理済素子データの生成を、マルチライン処理と言う。
本発明において、中心素子とは、送信の開口数(超音波の送信を行う素子数)が奇数の場合には、方位方向の中央の素子である。
他方、開口数が偶数の場合には、方位方向の中央の素子のいずれかを中心素子とし、あるいは、方位方向の真中に素子が有ると仮定して、中心素子とする。すなわち、開口数が偶数の場合には、開口の真中のライン上に焦点が有るものとして計算を行ってもよい。
他方、開口数が偶数の場合には、方位方向の中央の素子のいずれかを中心素子とし、あるいは、方位方向の真中に素子が有ると仮定して、中心素子とする。すなわち、開口数が偶数の場合には、開口の真中のライン上に焦点が有るものとして計算を行ってもよい。
なお、重ね合わせ処理部50における重ね合わせ処理方法としては、単に、加算するだけでなく、平均値や中央値をとってもよいし、係数を掛け合わせた上(重み付して)で加算してもよい。なお、平均値や中央値を取ることは、素子データレベルでの平均化フィルタやメディアンフィルタを掛けることに相当すると考えられるが、平均化フィルタやメディアンフィルタの代わりに、通常の画像処理で行われる逆フィルタなども適用してもよい。或いは、重ね合わせる各素子データ同士を比較し、類似している場合には最大値、類似していない場合には平均値、分布の偏りがある場合には中間値をとるなど、これに限らず、重ね合わせる各素子データの特徴量に基づいて重ね合わせ処理を変えてもよい。
また、注目素子の素子データに重ね合わせる素子データの数は、図示例の2つに限定はされず、1つでもよく、あるいは、3つ以上でもよい。すなわち、注目素子の素子データに重ね合わせる素子データの数は、要求される処理速度(フレームレート等)や画質などに応じて、適宜、設定すればよい。
また、重ね合わせる素子データ数は、超音波ビームのビーム幅の広がり程度に合わせた方が望ましい。従って、深さによってビーム幅が変わる場合には、重ね合わせ素子データ数も深さによって変更してもよい。また、ビーム幅は送信開口数に依存することから、送信開口数に応じて重ね合わせ素子データ数を変更してもよい。或いは、画像の輝度値などの特徴量に基づいて重ね合わせ素子データ数を変更してもよいし、重ね合わせ素子データ数を複数パターン変えて作成した画像から最適な重ね合わせ素子データ数を選択してもよい。
重ね合わせた結果、上述したように、真の信号の素子データでは信号の位相が合うが、ゴーストでは信号の位相が合わないため、加算などの重ね合わせ処理の結果、様々な位相の信号がお互いに打ち消し合い、信号が弱くなる。結果的に、真の信号は、有効な値を持つ、例えば高輝度の素子データとして残り、ゴーストの信号は、減弱した値を持つ、例えば低輝度の素子データとして得ることができる。
本発明に用いられる素子データ処理部は、基本的に以上のように構成される。
また、重ね合わせる素子データ数は、超音波ビームのビーム幅の広がり程度に合わせた方が望ましい。従って、深さによってビーム幅が変わる場合には、重ね合わせ素子データ数も深さによって変更してもよい。また、ビーム幅は送信開口数に依存することから、送信開口数に応じて重ね合わせ素子データ数を変更してもよい。或いは、画像の輝度値などの特徴量に基づいて重ね合わせ素子データ数を変更してもよいし、重ね合わせ素子データ数を複数パターン変えて作成した画像から最適な重ね合わせ素子データ数を選択してもよい。
重ね合わせた結果、上述したように、真の信号の素子データでは信号の位相が合うが、ゴーストでは信号の位相が合わないため、加算などの重ね合わせ処理の結果、様々な位相の信号がお互いに打ち消し合い、信号が弱くなる。結果的に、真の信号は、有効な値を持つ、例えば高輝度の素子データとして残り、ゴーストの信号は、減弱した値を持つ、例えば低輝度の素子データとして得ることができる。
本発明に用いられる素子データ処理部は、基本的に以上のように構成される。
再び、図2に戻って、音速値の変更及び最適音速の設定を網羅的に行う音速補正部24の各部についての説明を続ける。
図2に示す音速補正部24において、計算領域設定部70は、操作者による操作部32(図1参照)からの入力等に応じて、検査対象領域内の音速補正計算の対象とする領域、例えば、着目領域ROIを設定し、設定された領域のxy座標上の計算対象とする全ての点を抽出し、計算座標X1~Xendとして設定する。ここで、計算座標X1~Xendのそれぞれで表される各領域は、環境音速が同一であると考えられる領域であれば、特に制限的ではなく、被写体全体や検査対象領域全体であってもよいが、被写体や検査対象領域を区切った微小領域であってもよく、画素単位で1点ずつであっても良い。
初期計算座標設定部72は、計算座標Xiを、計算領域設定部70で設定された計算座標X1~Xendの中の初期計算座標X1に設定する。初期計算座標設定部72で設定された計算座標X1のデータは、最適音速決定部23に送られる。
図2に示す音速補正部24において、計算領域設定部70は、操作者による操作部32(図1参照)からの入力等に応じて、検査対象領域内の音速補正計算の対象とする領域、例えば、着目領域ROIを設定し、設定された領域のxy座標上の計算対象とする全ての点を抽出し、計算座標X1~Xendとして設定する。ここで、計算座標X1~Xendのそれぞれで表される各領域は、環境音速が同一であると考えられる領域であれば、特に制限的ではなく、被写体全体や検査対象領域全体であってもよいが、被写体や検査対象領域を区切った微小領域であってもよく、画素単位で1点ずつであっても良い。
初期計算座標設定部72は、計算座標Xiを、計算領域設定部70で設定された計算座標X1~Xendの中の初期計算座標X1に設定する。初期計算座標設定部72で設定された計算座標X1のデータは、最適音速決定部23に送られる。
なお、計算領域設定部70は、操作部32からの入力に応じて計算領域を設定する構成に限定はされず、画像生成部26(図1参照)の画像作成部44で生成された超音波画像(Bモード画像信号)を解析して、例えば着目領域ROIを設定し、計算領域とする構成としてもよい。
例えば、計算領域設定部70は、画像作成部44で生成された超音波画像において、隣接する画素同士の輝度値の差が所定の値以上となる位置(画素)を構造領域(組織や病変)の辺縁部(エッジ)と判定して、エッジに囲まれた領域、あるいは、エッジを多く含む領域を構造領域、即ち、対象部位Pと判定して抽出し、対象部位を含み、かつ、所定の形状、大きさとなる領域を、例えば着目領域ROIに設定し、計算領域に設定してもよい。
なお、超音波画像を解析して計算領域を設定する場合には、構造領域を1種類の領域とみなしても良く、あるいは、エッジの連続性や、輝度値のレベル等に応じて2種以上の領域とみなして、各領域を対象部位と判定して抽出し、対象部位毎に着目領域ROIを設定し、計算領域としても良い。
例えば、計算領域設定部70は、画像作成部44で生成された超音波画像において、隣接する画素同士の輝度値の差が所定の値以上となる位置(画素)を構造領域(組織や病変)の辺縁部(エッジ)と判定して、エッジに囲まれた領域、あるいは、エッジを多く含む領域を構造領域、即ち、対象部位Pと判定して抽出し、対象部位を含み、かつ、所定の形状、大きさとなる領域を、例えば着目領域ROIに設定し、計算領域に設定してもよい。
なお、超音波画像を解析して計算領域を設定する場合には、構造領域を1種類の領域とみなしても良く、あるいは、エッジの連続性や、輝度値のレベル等に応じて2種以上の領域とみなして、各領域を対象部位と判定して抽出し、対象部位毎に着目領域ROIを設定し、計算領域としても良い。
最適音速決定部23は、初期計算座標設定部72で設定された初期計算座標X1、又は計算座標更新部90で更新された計算座標Xiにおける最適な環境音速値(以下、最適音速値という)を算出して設定する部位である。
ここで、最適音速値とは、設定音速(設定環境音速)Vを種々変更し、それぞれの設定音速Vに基づいて、受信フォーカス処理をして超音波画像を形成した場合に、画像のコントラストおよび/またはシャープネスが最も高くなる設定音速Vの値であり、例えば、特開平8-317926号公報に記載のように、画像のコントラスト、スキャン方向の空間周波数、分散等に基づいて最適音速値の判定を行うことができる。
ここで、最適音速値とは、設定音速(設定環境音速)Vを種々変更し、それぞれの設定音速Vに基づいて、受信フォーカス処理をして超音波画像を形成した場合に、画像のコントラストおよび/またはシャープネスが最も高くなる設定音速Vの値であり、例えば、特開平8-317926号公報に記載のように、画像のコントラスト、スキャン方向の空間周波数、分散等に基づいて最適音速値の判定を行うことができる。
最適音速決定部23は、素子データ保持部20から読み出された、計算座標Xiに関わる未処理素子データを入力値として用いて最適な環境音速値(以下、最適音速値ともいう)を出力値として網羅的に求めて設定するためのものである。
最適音速決定部23としては、上述したように、所定の音速探索範囲を網羅的に探索して最適音速値を設定するものを挙げることができ、音速探索範囲設定部74と、初期音速値設定部76と、素子データ処理部22と、整相加算部38aと、検波処理部40aと、画像作成部44aと、画質指標算出部78と、最終音速値判定部80と、最適音速判定部84と、音速値更新部82と、音速設定部86と、を有する。
最適音速決定部23としては、上述したように、所定の音速探索範囲を網羅的に探索して最適音速値を設定するものを挙げることができ、音速探索範囲設定部74と、初期音速値設定部76と、素子データ処理部22と、整相加算部38aと、検波処理部40aと、画像作成部44aと、画質指標算出部78と、最終音速値判定部80と、最適音速判定部84と、音速値更新部82と、音速設定部86と、を有する。
音速探索範囲設定部74は、計算座標Xiにおける環境音速の探索範囲(初期音速値及び最終音速値)及び音速値変更幅(音速間隔)を設定するためのものである。音速探索範囲設定部74は、例えば、音速の探索開始の初期音速値Vst、音速の探索終了の最終音速値Vend、設定音速Vを変更するステップ幅(音速変更幅:音速間隔)ΔVを設定する。一例として、生体等の場合、音速探索範囲設定部74は、初期音速値Vstを1400m/sに、最終音速値Vendを1700m/sに、音速変更幅ΔVを10m/sに設定することができる。
なお、音速探索範囲設定部74が設定する初期音速値Vst、最終音速値Vend及び音速変更幅ΔVは、上述の値に限定されず、被検体やその検査対象領域やその計算座標Xi等に応じて決定されるものであればよい。また、これらの初期音速値Vst、最終音速値Vend及び音速変更幅ΔVは、音速探索範囲設定部74に、被検体やその検査対象領域やその計算座標Xi等に応じて、操作者によって操作部32を介して入力されるものであっても良いし、予め、格納部34に格納されているものの中から読み出されるものであっても良い。
初期音速値設定部76は、設定音速Vを音速探索範囲設定部74が設定した初期音速値Vstに設定するものである。上述の生体の場合、初期音速値設定部76は、設定音速Vを初期音速値Vstとして、1400m/sに設定する。初期音速値設定部76で設定された初期音速値Vstのデータは、後段の素子データ処理部22に送られる。
初期音速値設定部76は、設定音速Vを音速探索範囲設定部74が設定した初期音速値Vstに設定するものである。上述の生体の場合、初期音速値設定部76は、設定音速Vを初期音速値Vstとして、1400m/sに設定する。初期音速値設定部76で設定された初期音速値Vstのデータは、後段の素子データ処理部22に送られる。
整相加算部38aは、整相加算部38と全く同様の構成及び機能を有し、設定音速Vを用いて素子データ処理部22で求められた計算座標Xiに関わる処理済素子データに対して、音速補正部24で設定された最適音速に基づいて遅延を与えて加算する受信フォーカス処理を行って音線信号を生成する。
検波処理部40aは、検波処理部40aと全く同様の構成及び機能を有し、整相加算部38aで生成された音線信号に対して、包絡線検波処理等の処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像データを生成する。
画像作成部44aは、画像作成部44と全く同様の構成及び機能を有し、検波処理部40aで生成されたBモード画像データに画像処理を施して超音波画像を作成する。
検波処理部40aは、検波処理部40aと全く同様の構成及び機能を有し、整相加算部38aで生成された音線信号に対して、包絡線検波処理等の処理を施すことにより、被検体内の組織に関する断層画像情報であるBモード画像データを生成する。
画像作成部44aは、画像作成部44と全く同様の構成及び機能を有し、検波処理部40aで生成されたBモード画像データに画像処理を施して超音波画像を作成する。
画質指標算出部78は、画像作成部44aで作成された超音波画像から設定音速Vにおける画質指標、即ち画像の輝度値や鮮鋭度等の画質指標、例えば、各設定音速Vでの計算座標Xi(例えば、着目領域ROI)の画像の鮮鋭度(シャープネス)等を算出する。算出された画質指標は、設定音速Vに対応付けられて画質指標算出部78内に保持される、例えば、画質指標算出部78内のメモリ等に保持される。
最終音速値判定部80は、設定音速Vを最終音速値Vendと比較して、両者が等しいか否かの判定を行い、設定音速Vが最終音速値Vendに等しくない(V≠Vend)場合には、音速値更新部82に移行させ、設定音速Vが最終音速値Vendに等しい(V=Vend)場合には、最適音速判定部84に移行させる。
音速値更新部82は、現在の設定音速Vに音速変更幅ΔVを加算して新しい設定音速V(V=V+ΔV)を設定する。音速値更新部82で更新設定された設定音速V(V=V+ΔV)のデータは、前段の素子データ処理部22に送られる。
最終音速値判定部80は、設定音速Vを最終音速値Vendと比較して、両者が等しいか否かの判定を行い、設定音速Vが最終音速値Vendに等しくない(V≠Vend)場合には、音速値更新部82に移行させ、設定音速Vが最終音速値Vendに等しい(V=Vend)場合には、最適音速判定部84に移行させる。
音速値更新部82は、現在の設定音速Vに音速変更幅ΔVを加算して新しい設定音速V(V=V+ΔV)を設定する。音速値更新部82で更新設定された設定音速V(V=V+ΔV)のデータは、前段の素子データ処理部22に送られる。
最適音速判定部84は、画質指標算出部78に保持されている各音速値V(Vst~Vend)毎に求められた画質指標を読み出し、各設定音速Vでの画質指標、例えば、超音波画像の鮮鋭度(シャープネス値)を比較し、最も画質がよい時、例えば、超音波画像の鮮鋭度(シャープネス値)が最も高い時の設定音速Vを判定するものである。
音速設定部86は、最適音速判定部84において生成された超音波画像の画質が最もよいと判定された設定音速Vを最適音速値として設定するものである。
なお、図2に示す最適音速決定部23では、設定音速Vの変更を繰り返して処理しているが、本発明はこれに限定されず、GPUなどのマルチスレッドを用いて、全探索音速値を並列計算し、一度に全ての音速値における画質指標を求めることも可能である。
最適音速決定部23は、音速設定部86において計算座標Xiに対応する最適音速値を設定した後、最終計算座標判定部88に移行する。
最適音速決定部23は、基本的に以上のように構成され、計算座標Xiに係る最適音速値を算出する。
音速設定部86は、最適音速判定部84において生成された超音波画像の画質が最もよいと判定された設定音速Vを最適音速値として設定するものである。
なお、図2に示す最適音速決定部23では、設定音速Vの変更を繰り返して処理しているが、本発明はこれに限定されず、GPUなどのマルチスレッドを用いて、全探索音速値を並列計算し、一度に全ての音速値における画質指標を求めることも可能である。
最適音速決定部23は、音速設定部86において計算座標Xiに対応する最適音速値を設定した後、最終計算座標判定部88に移行する。
最適音速決定部23は、基本的に以上のように構成され、計算座標Xiに係る最適音速値を算出する。
最終計算座標判定部88は、計算座標Xiを最終計算座標Xendと比較して、両者が等しいか否かの判定を行い、計算座標Xiが最終計算座標Xendに等しくない(X≠Xend)場合には、計算座標更新部90に移行させ、計算座標Xiが最終計算座標Xendに等しい(X=Xend)場合には、画像生成部26の整相加算部38に移行させ、計算座標X1~Xendの各々において算出された最適音速値及び最適音速値を用いて素子データ処理部22で処理された処理済素子データを計算座標X1~Xendの各々の座標データと共に整相加算部38に供給する。
計算座標更新部90は、現在の計算座標Xiを新しい次の計算座標Xi+1に更新する。計算座標更新部90で更新された新しい次の計算座標Xi+1のデータは、前段の最適音速決定部23の音速探索範囲設定部74に送られる。
計算座標更新部90は、現在の計算座標Xiを新しい次の計算座標Xi+1に更新する。計算座標更新部90で更新された新しい次の計算座標Xi+1のデータは、前段の最適音速決定部23の音速探索範囲設定部74に送られる。
音速補正部24は、基本的に以上のように構成され、全ての計算座標Xi(i=1~end)各々に係る最適音速値及び処理済素子データを算出する。
なお、上述したように、本実施形態の超音波検査装置10の音速補正部24による音速補正における最適音速を探索するループの中にマルチライン処理を行う素子データ処理部22を組み込み、マルチライン処理時の遅延時間補正を計算する際の音速値に、音速補正部24で探索する音速値を用いることで、反射点(サンプリング点、観測点)から超音波プローブ12の探触子36の受信素子までの音速だけでなく、各々の仮想的な点音源を形成する探触子36の送信素子から反射点までの音速も含めて最適な環境音速値を求めることができる。
なお、上述したように、本実施形態の超音波検査装置10の音速補正部24による音速補正における最適音速を探索するループの中にマルチライン処理を行う素子データ処理部22を組み込み、マルチライン処理時の遅延時間補正を計算する際の音速値に、音速補正部24で探索する音速値を用いることで、反射点(サンプリング点、観測点)から超音波プローブ12の探触子36の受信素子までの音速だけでなく、各々の仮想的な点音源を形成する探触子36の送信素子から反射点までの音速も含めて最適な環境音速値を求めることができる。
本発明の超音波検査装置は、基本的に以上のように構成されるが、以下に、図1及び図2を参照して、本発明の超音波検査装置の動作、作用及び本発明の超音波検査方法について説明する。
操作者が、図1に示す超音波プローブ12を被検体の表面(検査対象領域)に当接させて、測定を開始すると、送信部14から供給される駆動信号に従って探触子36の複数の素子から超音波ビームが送信され、被検体からの超音波エコーを、探触子36の複数の素子が受信し、受信信号としてアナログ素子信号を出力する。
受信部16は、アナログ素子信号を増幅してA/D変換部18に供給し、A/D変換部18は、アナログ素子信号をデジタル素子データに変換して素子データ保持部20に供給して、未処理素子データとして記憶保持させる。
操作者が、図1に示す超音波プローブ12を被検体の表面(検査対象領域)に当接させて、測定を開始すると、送信部14から供給される駆動信号に従って探触子36の複数の素子から超音波ビームが送信され、被検体からの超音波エコーを、探触子36の複数の素子が受信し、受信信号としてアナログ素子信号を出力する。
受信部16は、アナログ素子信号を増幅してA/D変換部18に供給し、A/D変換部18は、アナログ素子信号をデジタル素子データに変換して素子データ保持部20に供給して、未処理素子データとして記憶保持させる。
次に、図2に示す音速補正部24において、操作者による操作部32(図1参照)からの入力等に応じて、計算領域設定部70によって計算座標X1~Xendが設定された後、初期計算座標設定部72によって計算座標Xiが、初期計算座標X1に設定され、設定された計算座標X1のデータが、最適音速決定部23に送られる。
最適音速決定部23では、まず、音速探索範囲設定部74が、被検体やその検査対象領域やその計算座標Xi等に応じて、操作者によって操作部32を介して入力された、又は予め、格納部34に格納されている音速探索範囲の情報の中から読み出された初期音速値Vst、最終音速値Vend及び音速変更幅ΔVを設定し、例えば、上述したように、生体では、初期音速値Vstを1400m/sに、音速値Vendを1700m/sに、音速変更幅ΔVを10m/sに設定する。
続いて、初期音速値設定部76が、設定音速Vを初期音速値Vst(例えば、1400m/s)に設定して、素子データ処理部22に移行する。
続いて、初期音速値設定部76が、設定音速Vを初期音速値Vst(例えば、1400m/s)に設定して、素子データ処理部22に移行する。
次に、素子データ処理部22は、遅延時間算出部48(図3参照)において、真の信号の未処理素子データに対する周辺のゴーストの信号の未処理素子データの遅延時間(例えば、図7(b)、図7(e)、どちらも同じものである)を、送信素子、焦点、反射点、及び受信素子の幾何学的配置、及び予め初期音速値設定部76で設定された設定音速V等から算出する(例えば、図6の幾何学モデルを用いて算出する)。
次に、素子データ処理部22は、素子データ保持部20から計算座標X1に関わる未処理素子データを素子データ保持部20から読み出し、処理を行う素子データを注目素子データとし、重ね合わせ処理部50(図3参照)において、遅延時間算出部48で算出された遅延時間を用いて、注目素子データとその周辺の未処理素子データとを位相合わせて重ね合わせてマルチライン処理を行い、処理済素子データを求める。これにより、真の信号を含む未処理素子データであれば、強調された処理済素子データが求まると共に、ゴーストの未処理素子データであれば減弱した処理済素子データが求まる。
素子データ処理部22は、こうして求めた処理済素子データを整相加算部38aに供給する。
素子データ処理部22は、こうして求めた処理済素子データを整相加算部38aに供給する。
次に、整相加算部38aは、初期音速値値Vstである設定音速Vに基づいて、素子データ処理部22から供給された計算座標X1に対応する処理済素子データを用いて受信フォーカス処理をして整相加算を行い音線信号を生成する。次いで、検波処理部40aは、生成された音線信号に対して検波処理を行ってBモード画像信号を生成する。そして、画像作成部44aは、生成されたBモード画像信号から超音波画像を作成する。
その後、画質指標算出部78は、作成された超音波画像から設定音速Vにおける超音波画像の輝度値や鮮鋭度等の画質指標、例えば、各設定音速VでのXi(X1)の画像の鮮鋭度等を算出する。算出された画質指標は、設定音速値Vに対応付けられて画質指標算出部78内に保持される。
次に、最終音速値判定部80は、設定音速Vを最終音速値Vendと比較して、両者が等しいか否かの判定を行う。最終音速値判定部80の判定において、設定音速Vが最終音速値Vendに等しい(V=Vend)場合には、最適音速判定部84に移行し、設定音速Vが最終音速値Vendに等しくない(V≠Vend)場合には、音速値更新部82に移行する。
次に、最終音速値判定部80は、設定音速Vを最終音速値Vendと比較して、両者が等しいか否かの判定を行う。最終音速値判定部80の判定において、設定音速Vが最終音速値Vendに等しい(V=Vend)場合には、最適音速判定部84に移行し、設定音速Vが最終音速値Vendに等しくない(V≠Vend)場合には、音速値更新部82に移行する。
音速値更新部82に移行した場合には、音速値更新部82が、現在の設定音速Vに音速変更幅ΔVを加算して新しい設定音速V(V=V+ΔV)を設定する。
この後、素子データ処理部22に戻り、新しい設定音速に基づいて素子データ処理部22によるマルチライン処理、整相加算部38aによる整相加算、検波処理部40aによる検波処理、画像作成部44aによる画像作成、画質指標算出部78による画質指標の算出及び保持、並びに最終音速値判定部80による判定が行われる。
最終音速値判定部80による判定結果がV≠Vend場合には、上述したように、音速値更新部82に移行し、音速値更新部82による新しい設定音速Vの設定、新しい設定音速Vに基づく素子データ処理部22によるマルチライン処理、整相加算部38aによる整相加算、検波処理部40aによる検波処理、画像作成部44aによる画像作成、画質指標算出部78による画質指標の算出及び保持、並びに最終音速値判定部80による判定が、V=Vendになるまで繰り返される。即ち、これらの処理を、音速範囲の終了の最終音速値Vendになるまで設定音速の値をΔVずつ変更しながら、例えば、設定音速Vを、1400m/sから1700m/sまで、10m/s刻みで変化させながら繰り返し行う。
この後、素子データ処理部22に戻り、新しい設定音速に基づいて素子データ処理部22によるマルチライン処理、整相加算部38aによる整相加算、検波処理部40aによる検波処理、画像作成部44aによる画像作成、画質指標算出部78による画質指標の算出及び保持、並びに最終音速値判定部80による判定が行われる。
最終音速値判定部80による判定結果がV≠Vend場合には、上述したように、音速値更新部82に移行し、音速値更新部82による新しい設定音速Vの設定、新しい設定音速Vに基づく素子データ処理部22によるマルチライン処理、整相加算部38aによる整相加算、検波処理部40aによる検波処理、画像作成部44aによる画像作成、画質指標算出部78による画質指標の算出及び保持、並びに最終音速値判定部80による判定が、V=Vendになるまで繰り返される。即ち、これらの処理を、音速範囲の終了の最終音速値Vendになるまで設定音速の値をΔVずつ変更しながら、例えば、設定音速Vを、1400m/sから1700m/sまで、10m/s刻みで変化させながら繰り返し行う。
一方、判定部80による判定結果がV=Vendである場合には、音速の探索範囲内の全ての設定音速Vに関して画質指標が求められたことになるので、上述したように、最適音速判定部84に移行し、最適音速判定部84は、画質指標算出部78に保持されている各音速値毎に求められた画質指標を読み出し、最適音速判定部84は、各設定音速Vでの画質指標、例えば、超音波画像の鮮鋭度(シャープネス値)を比較し、最も画質がよい時、例えば、超音波画像の鮮鋭度(シャープネス値)が最も高い時の設定音速Vを判定し、音速設定部86は、最適音速判定部84にて判定された設定音速Vを最適音速値として採用し、設定する。
こうして、最適音速決定部23は、計算座標Xi(X1)に対応する最適音速値を算出する。
こうして、最適音速決定部23は、計算座標Xi(X1)に対応する最適音速値を算出する。
次に、最適音速決定部23による計算座標Xi(X1)に対応する最適音速値の算出が終了すると、最終計算座標判定部88に移行する。
最終計算座標判定部88では、計算座標Xiを最終計算座標Xendと比較して、両者が等しいか否かの判定を行う。その判定の結果、計算座標Xiが最終計算座標Xendに等しくない(X≠Xend)場合には、計算座標更新部90に移行する。
この後、計算座標更新部90では、現在の計算座標Xiを新しい次の計算座標Xi+1に更新し、前段の最適音速決定部23に移行する。
この後、最適音速決定部23では、上述した設定音速Vの網羅的な探索を行い、計算座標Xi+1に対応する最適音速値を算出し、再び、最終計算座標判定部88に移行して、最終計算座標Xendの判定が行われる。
最終計算座標判定部88による判定結果がXi≠Xend場合には、上述したように、計算座標更新部90に移行し、計算座標更新部90による新しい計算座標Xiの設定、最適音速決定部23における設定音速Vの網羅的な探索による計算座標Xiに対応する最適音速値の算出、最終計算座標判定部88の判定が、Xi=Xendになるまで繰り返される。
最終計算座標判定部88では、計算座標Xiを最終計算座標Xendと比較して、両者が等しいか否かの判定を行う。その判定の結果、計算座標Xiが最終計算座標Xendに等しくない(X≠Xend)場合には、計算座標更新部90に移行する。
この後、計算座標更新部90では、現在の計算座標Xiを新しい次の計算座標Xi+1に更新し、前段の最適音速決定部23に移行する。
この後、最適音速決定部23では、上述した設定音速Vの網羅的な探索を行い、計算座標Xi+1に対応する最適音速値を算出し、再び、最終計算座標判定部88に移行して、最終計算座標Xendの判定が行われる。
最終計算座標判定部88による判定結果がXi≠Xend場合には、上述したように、計算座標更新部90に移行し、計算座標更新部90による新しい計算座標Xiの設定、最適音速決定部23における設定音速Vの網羅的な探索による計算座標Xiに対応する最適音速値の算出、最終計算座標判定部88の判定が、Xi=Xendになるまで繰り返される。
一方、最終計算座標判定部88における判定の結果、計算座標Xiが最終計算座標Xendに等しい(X=Xend)場合には、音速補正部24による最適音速値の設定が終了し、全ての計算座標Xi(i=1~end)各々に対応する最適音速値及び処理済素子データが算出され、画像生成部26の整相加算部38に移行する。こうして、計算座標X1~Xendの各々において算出された最適音速値及び最適音速値を用いて素子データ処理部22で処理された処理済素子データは、計算座標X1~Xendの各々の座標データと共に整相加算部38に供給される。
次に、画像生成部26の整相加算部38では、素子データに受信フォーカス処理を施して受信データ(音線信号)を生成し、検波処理部40に供給する。検波処理部40では、音線信号を処理してBモード画像信号を生成する。Bモード画像信号を、DSC42がラスター変換し、画像作成部44が画像処理を施し、超音波画像が生成される。生成された超音波画像は、画像メモリ46に格納されると共に、表示制御部27により超音波画像が表示部28に表示される。
このように、本発明の超音波検査装置は、素子データ処理部でマルチライン処理され、ゴーストの信号が弱減され、真の信号が強調された処理済素子データを用いて作成した超音波画像を評価する画質指標を用いて最適な環境音速を求め、求められた最適な環境音速を用いてマルチライン処理された処理済素子データを用いて超音波画像を作成するので、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができ、また、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高解像度で、SN比が高く、最適な空間分解能を持つシャープな超音波画像を得ることができる。
このように、本発明の超音波検査装置は、素子データ処理部でマルチライン処理され、ゴーストの信号が弱減され、真の信号が強調された処理済素子データを用いて作成した超音波画像を評価する画質指標を用いて最適な環境音速を求め、求められた最適な環境音速を用いてマルチライン処理された処理済素子データを用いて超音波画像を作成するので、生体組織内の音速ムラの影響を軽減させた高画質な超音波画像を得ることができ、また、幅広ビーム送信専用超音波プローブを使うことなく、かつ、従来と変わらないフレームレートのまま、送信ビームの幅に依存しない高解像度で、SN比が高く、最適な空間分解能を持つシャープな超音波画像を得ることができる。
次に、図9及び図10に基づいて、本発明の実施例2の超音波検査装置について説明する。
図9は、本発明の実施例2に係る超音波検査装置の構成を概念的に示すブロック図であり、図10は、図9に示す超音波検査装置の要部の一例を処理フローに沿って概念的に示すブロック図である。
なお、図9及び図10に示す超音波検査装置100は、図1及び図2に示す超音波検査装置10と、第2の音速補正部102を備えている点で異なる以外は、全く同様の構成を有するものであるので、同一の構成要素には、同一の参照符号を付し、その詳細な説明は省略する。
図9は、本発明の実施例2に係る超音波検査装置の構成を概念的に示すブロック図であり、図10は、図9に示す超音波検査装置の要部の一例を処理フローに沿って概念的に示すブロック図である。
なお、図9及び図10に示す超音波検査装置100は、図1及び図2に示す超音波検査装置10と、第2の音速補正部102を備えている点で異なる以外は、全く同様の構成を有するものであるので、同一の構成要素には、同一の参照符号を付し、その詳細な説明は省略する。
図9に示すように、超音波検査装置100は、超音波プローブ12と、超音波プローブ12に接続される送信部14及び受信部16と、A/D変換部18と、素子データ保持部20と、第2の音速補正部102と、素子データ処理部22を備える音速補正部24と、画像生成部26と、表示制御部27と、表示部28と、制御部30と、操作部32と、格納部34とを有する。
超音波検査装置100では、音速補正部24による音速補正を行う前に、第2の音速補正部102にて、未処理素子データ(第1の素子データ)を用いて音速補正を行って暫定的な音速値を求め、この音速値を、音速補正部24で音速補正を行う際の初期値として利用して最適音速値を求めるものである。
第2の音速補正部102は、素子データ保持部20と音速補正部24との間に接続されており、素子データ保持部20に保存された未処理素子データ(第1の素子データ)を用いて、音速補正を行い算出した音速値を音速補正部24に供給する。
第2の音速補正部102が行う音速補正には、特に限定はないが、例えば、特開2011-92686号に記載の方法が利用できる。この方法は、複数の設定音速毎に未処理素子データに対して受信フォーカス処理を行って設定音速毎にフォーカス指標を算出し、このフォーカス指標に基づいて最適な音速値を決定するものである。すなわち、未処理素子データを最も適切に受信フォーカス処理できる設定音速を最適な音速値とする。
第2の音速補正部102が行う音速補正には、特に限定はないが、例えば、特開2011-92686号に記載の方法が利用できる。この方法は、複数の設定音速毎に未処理素子データに対して受信フォーカス処理を行って設定音速毎にフォーカス指標を算出し、このフォーカス指標に基づいて最適な音速値を決定するものである。すなわち、未処理素子データを最も適切に受信フォーカス処理できる設定音速を最適な音速値とする。
なお、第2の音速補正部102での音速補正は、音速補正部24における音速補正の際の音速探索の初期値として大まかに求めればよいので、必ずしも計算領域や着目領域毎に送信フォーカスを設定する必要はない。また、計算領域や着目領域毎に行うものに限定はされず、画面全体を1つの着目領域として音速値を求めてもよい。
図10に示すように、第2の音速補正部102は、求めた音速値を音速補正部24の最適音速決定部23の音速探索範囲設定部74に供給する。
図10に示すように、第2の音速補正部102は、求めた音速値を音速補正部24の最適音速決定部23の音速探索範囲設定部74に供給する。
最適音速決定部23の音速探索範囲設定部74は、第2の音速補正部102から供給された音速値を初期音速値Vstとして音速値を設定する。
音速値更新部82は、初期音速値Vst及びこの前後の音速値(Vst±ΔV)を順次設定して、画質指標算出部78で画質指標を算出し、3つの音速値での結果に基づいて、次の設定音速Vを設定する。例えば、設定音速VがVst+ΔVの場合が最も良い結果であれば、音速値更新部82は、次の設定音速VとしてVst+2×ΔVを設定し、設定音速VがVstの場合が最も良い結果であれば、Vstを最適音速値として採用する。すなわち、最適音速決定部23は、第2の音速補正部102から供給された音速値を初期音速値として探索的に最適音速を求める。
第2の音速補正部102が未処理素子データを用いて音速補正を行い、音速補正部24での音速補正の初期値として用いる構成とすることにより、最適音速値の測定・演算にかかる時間を短縮することができる。
音速値更新部82は、初期音速値Vst及びこの前後の音速値(Vst±ΔV)を順次設定して、画質指標算出部78で画質指標を算出し、3つの音速値での結果に基づいて、次の設定音速Vを設定する。例えば、設定音速VがVst+ΔVの場合が最も良い結果であれば、音速値更新部82は、次の設定音速VとしてVst+2×ΔVを設定し、設定音速VがVstの場合が最も良い結果であれば、Vstを最適音速値として採用する。すなわち、最適音速決定部23は、第2の音速補正部102から供給された音速値を初期音速値として探索的に最適音速を求める。
第2の音速補正部102が未処理素子データを用いて音速補正を行い、音速補正部24での音速補正の初期値として用いる構成とすることにより、最適音速値の測定・演算にかかる時間を短縮することができる。
以上、本発明について詳細に説明したが、本発明は上記実施例に限定されず、本発明の主旨を逸脱しない範囲において、種々の改良や変更をしてもよいのはもちろんである。
例えば、図1及び図2に示す各構成要素をハードウエアとして構成いても良いし、コンピュータ等で実行されるソフトウエアとして構成しても良い。
更に、本発明の各実施形態について説明した超音波検査装置の各機能に対応して、例えば、上述した音速補正部に含まれる素子データ処理部を始めとする各構成要素、及び画像生成部として機能するようにコンピュータを動作させる超音波検査ためのプログラムも、また、上記超音波検査方法の各ステップ、例えば、上述したマルチライン処理の各ステップを含む音速補正の各ステップ、及び画像生成の各ステップを手順としてコンピュータに実行させる超音波検査のためのプログラムも、本発明の実施の形態の1つである。また、そのようなプログラムを記録したコンピュータ読み取り可能な記憶媒体も、本発明の実施の形態の1つである。
例えば、図1及び図2に示す各構成要素をハードウエアとして構成いても良いし、コンピュータ等で実行されるソフトウエアとして構成しても良い。
更に、本発明の各実施形態について説明した超音波検査装置の各機能に対応して、例えば、上述した音速補正部に含まれる素子データ処理部を始めとする各構成要素、及び画像生成部として機能するようにコンピュータを動作させる超音波検査ためのプログラムも、また、上記超音波検査方法の各ステップ、例えば、上述したマルチライン処理の各ステップを含む音速補正の各ステップ、及び画像生成の各ステップを手順としてコンピュータに実行させる超音波検査のためのプログラムも、本発明の実施の形態の1つである。また、そのようなプログラムを記録したコンピュータ読み取り可能な記憶媒体も、本発明の実施の形態の1つである。
10、100 超音波検査装置
12 超音波プローブ
14 送信部
16 受信部
18 A/D変換部
20 素子データ保持部
22 素子データ処理部
23 最適音速決定部
24 音速補正部
26 画像生成部
27 表示制御部
28 表示部
30 制御部
32 操作部
34 格納部
36 探触子
38、38a 整相加算部
40、40a 検波処理部
42 DSC
44、44a 画像作成部
46 画像メモリ
48 遅延時間算出部
50 重ね合わせ処理部
70 計算領域設定部
72 初期計算座標設定部
74 音速探索範囲設定部
76 初期音速値設定部
78 画質指標算出部
80 最終音速値判定部
82 音速値更新部
84 最適音速判定部
86 音速設定部
88 最終計算座標判定部
90 計算座標更新部
102 第2の音速補正部
12 超音波プローブ
14 送信部
16 受信部
18 A/D変換部
20 素子データ保持部
22 素子データ処理部
23 最適音速決定部
24 音速補正部
26 画像生成部
27 表示制御部
28 表示部
30 制御部
32 操作部
34 格納部
36 探触子
38、38a 整相加算部
40、40a 検波処理部
42 DSC
44、44a 画像作成部
46 画像メモリ
48 遅延時間算出部
50 重ね合わせ処理部
70 計算領域設定部
72 初期計算座標設定部
74 音速探索範囲設定部
76 初期音速値設定部
78 画質指標算出部
80 最終音速値判定部
82 音速値更新部
84 最適音速判定部
86 音速設定部
88 最終計算座標判定部
90 計算座標更新部
102 第2の音速補正部
Claims (15)
- 超音波ビームを用いて検査対象物を検査する超音波検査装置であって、
前記超音波ビームを送信し、かつ、前記検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子と、
前記探触子に、複数の前記素子を用い、所定の送信焦点を形成するように前記超音波ビームを送信させることを、複数回、行わせる送信部と、
個々の前記超音波ビームの送信に対応して、複数の前記素子が出力したアナログ素子信号を受け、所定の処理を施す受信部と、
前記受信部が処理したアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとするAD変換部と、
前記検査対象物内の音速値を決定する音速決定部と、
前記検査対象物内の音速値を用いて2つ以上の前記第1の素子データから、前記2つ以上の第1の素子データのいずれかに対応する第2の素子データを生成する素子データ処理部と、を有し、
前記音速決定部は、前記素子データ処理部において前記2つ以上の第1の素子データから前記第2の素子データを作成する時に用いる前記音速値を最適化して最適な音速値として求めることを特徴とする超音波検査装置。 - 前記音速決定部は、前記素子データ処理部を有し、前記素子データ処理部における前記2つ以上の第1の素子データの遅延時間補正計算に用いる前記音速値を最適化して最適な音速として求める請求項1に記載の超音波検査装置。
- 前記音速決定部は、前記検査対象物内を複数の領域に分割したときの各領域における音速値を決定する請求項1又は2に記載の超音波検査装置。
- 更に、前記2つ以上の第1の素子データを含む一部の第1の素子データ、あるいは全ての第1の素子データを保持する素子データ保持部を有する請求項1~3のいずれか1項に記載の超音波検査装置。
- 前記素子データ保持部は、前記送信部により、前記検査対象物内における少なくとも2つ以上の重なり合う対象領域毎において、前記超音波ビームを送信し、前記受信部により、前記超音波ビームによって前記検査対象物内における少なくとも2つ以上の重なり合う対象領域毎に発生した前記超音波エコー信号を受信して生成された、素子データ毎に受信した各素子における受信時間毎の受信データを含む前記2以上の第1の素子データを保持し、
前記素子データ処理部は、前記素子データ毎の前記超音波ビームを送信したときの素子の幾何学的な配置及び前記素子データ毎の前記受信データの各素子の幾何学的な配置に基づいて、前記2以上の第1の素子データから前記第2の素子データを生成する請求項4に記載の超音波検査装置。 - 前記送信部は、前記超音波ビームの送信を、中心となる素子を変更して、複数回、行う請求項5に記載の超音波診断装置。
- 前記素子データ処理部は、前記検査対象物内の音速値を用いて前記2つ以上の第1の素子データの各遅延時間を算出する遅延時間算出部と、前記2つ以上の第1の素子データを算出された各遅延時間及び受信された前記探触子の素子の位置に基づいて重ね合わせ、前記第2の素子データを生成する重ね合わせ処理部とを有する請求項1~6のいずれか1項に記載の超音波検査装置。
- 更に、前記素子データ処理部で生成された前記第2の素子データに基づいて超音波画像を生成する画像生成部と、を有し、
前記素子データ処理部は、前記画像生成部で生成されるべき超音波画像上の1ラインの前記第2の素子データを作成する請求項1~7のいずれか1項に記載の超音波検査装置。 - 前記音速決定部は、
前記素子データ処理部と、
前記素子データ処理部における前記遅延時間算出に用いられる前記検査対象物の音速値を変更する音速変更部と、
前記音速変更部で変更された音速値を用いて前記第2の素子データを作成する前記素子データ処理部と、
前記素子データ処理部により作成された前記第2の素子データを、前記音速変更部で変更された音速値を用いて整相加算して音線信号を生成する整相加算部と、
前記整相加算部で生成された前記音線信号に基づいて超音波画像を生成する画像生成部と、
前記画像生成部で生成された前記超音波画像の画質指標を算出する画質指標算出部と、
前記画質指標算出部で算出された前記超音波画像の画質指標を判定する判定部と、
前記判定部による前記画質指標の判定結果に応じて前記最適な音速値を設定する音速設定部と、を有する請求項1~8のいずれか1項に記載の超音波検査装置。 - 前記音速決定部は、
前記音速変更部が、前記音速値を初期音速値から最終音速値まで所定音速間隔で変更する度毎に、前記素子データ処理部が変更された当該音速値を用いて前記第2の素子データを作成し、前記整相加算部が同じ当該音速値を用いて前記第2の素子データを整相加算して前記音線信号を生成し、前記画像生成部が前記音線信号に基づいて前記超音波画像を生成し、前記画質指標算出部が前記超音波画像の前記画質指標を算出することを繰り返して、前記音速変更部によって変更された全ての複数の前記音速値にそれぞれ対応する複数の前記画質指標を算出しておき、
前記判定部が、前記画質指標算出部で全ての複数の前記音速値毎に算出された複数の前記画質指標の中から最適な画質指標を判定し、前記音速設定部が、前記判定部で判定された前記最適な画質指標に対応する音速値を前記最適な音速値として設定する請求項9に記載の超音波検査装置。 - さらに、前記音速決定部よりも前段に、前記第1の素子データを用いて最適な音速値を算出する第2の音速決定部を有し、
前記音速決定部は、前記第2の音速決定部で算出した音速値を初期値として用いる請求項1~10のいずれか1項に記載の超音波検査装置。 - 超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して、受信した超音波エコーに応じたアナログ素子信号を出力する、複数の素子が配列された探触子を用いて検査対象物を検査する超音波検査方法であって、
前記探触子によって、複数の前記素子を用い、所定の送信焦点を形成するように超音波ビームを送信させることを、複数回、行い、かつ、個々の前記超音波ビームの送信に対応して、複数の前記素子がアナログ素子信号を出力し、
このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データとし、
前記検査対象物内の音速値を用いて、2つ以上の前記第1の素子データから、前記第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる前記音速値を最適化して最適な音速値として決定することを特徴とする超音波検査方法。 - 前記最適な音速を決定するステップは、
前記音速値を初期音速値から最終音速値まで所定音速間隔で変更する毎に、変更された当該音速値を用いて前記第2の素子データを作成し、同じ当該音速値を用いて前記第2の素子データを整相加算して前記音線信号を生成し、生成された前記音線信号に基づいて前記超音波画像を生成し、前記画質指標算出部が前記超音波画像の前記画質指標を算出することを繰り返して、変更された全ての複数の前記音速値毎にそれぞれ対応する複数の前記画質指標を算出しておき、
全ての複数の前記音速値毎に算出された複数の前記画質指標の中から最適な画質指標を判定し、判定された前記最適な画質指標に対応する音速値を前記最適な音速値として設定する請求項12に記載の超音波検査方法。 - 超音波ビームを送信し、かつ、検査対象物によって反射された超音波エコーを受信して検査対象物を検査するために、複数の素子が配列された探触子によって、複数の前記素子を用い、所定の送信焦点を形成するように、超音波ビームを送信させることを、複数回、行い、かつ、個々の超音波ビームの送信に対応して、複数の前記素子がアナログ素子信号を出力し、このアナログ素子信号をA/D変換して、デジタル素子信号である第1の素子データを複数、取得するステップと、
前記検査対象物内の音速値を用いて、得られた2つ以上の前記第1の素子データから前記第1の素子データのいずれかに対応する第2の素子データを作成する時に用いる前記音速値を最適化して最適な音速値として決定するステップとを手順としてコンピュータに実行させることを特徴とするプログラム。 - 請求項14に記載のプログラムを記録したコンピュータ読み取り可能な記録媒体。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201380050604.0A CN104684487A (zh) | 2012-09-28 | 2013-09-24 | 超声波检查装置、超声波检查方法、程序及记录介质 |
US14/671,236 US10687786B2 (en) | 2012-09-28 | 2015-03-27 | Ultrasound inspection apparatus, ultrasound inspection method and recording medium |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-215721 | 2012-09-28 | ||
JP2012215721 | 2012-09-28 | ||
JP2013-125337 | 2013-06-14 | ||
JP2013125337A JP5946427B2 (ja) | 2012-09-28 | 2013-06-14 | 超音波検査装置、超音波検査方法、プログラム及び記録媒体 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/671,236 Continuation US10687786B2 (en) | 2012-09-28 | 2015-03-27 | Ultrasound inspection apparatus, ultrasound inspection method and recording medium |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014050797A1 true WO2014050797A1 (ja) | 2014-04-03 |
Family
ID=50388196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2013/075666 WO2014050797A1 (ja) | 2012-09-28 | 2013-09-24 | 超音波検査装置、超音波検査方法、プログラム及び記録媒体 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10687786B2 (ja) |
JP (1) | JP5946427B2 (ja) |
CN (1) | CN104684487A (ja) |
WO (1) | WO2014050797A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104188684A (zh) * | 2014-09-15 | 2014-12-10 | 声泰特(成都)科技有限公司 | 一种自适应医疗超声成像声速优化和信号修正方法及系统 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104303184B (zh) * | 2012-03-21 | 2018-05-15 | 皇家飞利浦有限公司 | 整合医疗成像和活检数据的临床工作站以及使用其的方法 |
JP6006249B2 (ja) * | 2014-03-24 | 2016-10-12 | 富士フイルム株式会社 | 音響波処理装置、音響波処理装置の信号処理方法およびプログラム |
JP5815784B2 (ja) | 2014-04-08 | 2015-11-17 | ファナック株式会社 | 同期加工における同期誤差を低減するサーボ制御装置 |
KR20150118732A (ko) * | 2014-04-15 | 2015-10-23 | 삼성전자주식회사 | 초음파 장치 및 그 제어 방법 |
WO2019044075A1 (ja) * | 2017-08-31 | 2019-03-07 | 富士フイルム株式会社 | 画像生成装置および作動方法 |
CN110680398A (zh) * | 2019-10-11 | 2020-01-14 | 深圳开立生物医疗科技股份有限公司 | 一种超声回波成像方法、系统及超声设备 |
CN115867202A (zh) * | 2020-07-09 | 2023-03-28 | 深圳迈瑞生物医疗电子股份有限公司 | 超声声速确定方法和超声成像系统 |
CN116529629A (zh) * | 2020-12-01 | 2023-08-01 | 皇家飞利浦有限公司 | 使用发射速度系统、设备和方法的回顾性发射聚焦 |
US20230061869A1 (en) * | 2021-08-26 | 2023-03-02 | GE Precision Healthcare LLC | System and methods for beamforming sound speed selection |
US11867807B2 (en) * | 2021-09-01 | 2024-01-09 | GE Precision Healthcare LLC | System and methods for beamforming sound speed selection |
CN114521917A (zh) * | 2022-02-15 | 2022-05-24 | 清华大学 | 超声重新波束形成方法、装置和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004522515A (ja) * | 2001-01-26 | 2004-07-29 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 空間的複合のための適応画像処理 |
JP2009090102A (ja) * | 2007-09-18 | 2009-04-30 | Fujifilm Corp | 超音波診断方法及び装置 |
JP2010082190A (ja) * | 2008-09-30 | 2010-04-15 | Fujifilm Corp | 超音波信号処理装置及び方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4893284A (en) * | 1988-05-27 | 1990-01-09 | General Electric Company | Calibration of phased array ultrasound probe |
JPH0595946A (ja) * | 1991-10-09 | 1993-04-20 | Yokogawa Medical Syst Ltd | 超音波の音速測定法及び音速測定手段を備えた超音波診断装置 |
EP1974672B9 (en) * | 2007-03-28 | 2014-04-16 | Kabushiki Kaisha Toshiba | Ultrasonic imaging apparatus and ultrasonic velocity optimization method |
JP2009090104A (ja) | 2007-09-18 | 2009-04-30 | Fujifilm Corp | 超音波診断方法及び装置 |
JP2009089940A (ja) * | 2007-10-10 | 2009-04-30 | Toshiba Corp | 超音波診断装置 |
JP5355924B2 (ja) | 2008-03-31 | 2013-11-27 | 株式会社東芝 | 超音波診断装置 |
JP2010099452A (ja) * | 2008-09-25 | 2010-05-06 | Fujifilm Corp | 超音波診断装置及び超音波診断方法 |
JP5566773B2 (ja) * | 2009-06-30 | 2014-08-06 | 株式会社東芝 | 超音波診断装置および音速設定方法 |
JP5486257B2 (ja) * | 2009-09-28 | 2014-05-07 | 富士フイルム株式会社 | 超音波診断装置及び弾性指標算出方法 |
JP2012192133A (ja) * | 2011-03-18 | 2012-10-11 | Fujifilm Corp | 超音波診断装置および超音波画像生成方法 |
JP5292440B2 (ja) * | 2011-06-03 | 2013-09-18 | 富士フイルム株式会社 | 超音波診断装置 |
-
2013
- 2013-06-14 JP JP2013125337A patent/JP5946427B2/ja active Active
- 2013-09-24 CN CN201380050604.0A patent/CN104684487A/zh active Pending
- 2013-09-24 WO PCT/JP2013/075666 patent/WO2014050797A1/ja active Application Filing
-
2015
- 2015-03-27 US US14/671,236 patent/US10687786B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004522515A (ja) * | 2001-01-26 | 2004-07-29 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 空間的複合のための適応画像処理 |
JP2009090102A (ja) * | 2007-09-18 | 2009-04-30 | Fujifilm Corp | 超音波診断方法及び装置 |
JP2010082190A (ja) * | 2008-09-30 | 2010-04-15 | Fujifilm Corp | 超音波信号処理装置及び方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104188684A (zh) * | 2014-09-15 | 2014-12-10 | 声泰特(成都)科技有限公司 | 一种自适应医疗超声成像声速优化和信号修正方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
US10687786B2 (en) | 2020-06-23 |
US20150196274A1 (en) | 2015-07-16 |
JP5946427B2 (ja) | 2016-07-06 |
CN104684487A (zh) | 2015-06-03 |
JP2014079565A (ja) | 2014-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5946427B2 (ja) | 超音波検査装置、超音波検査方法、プログラム及び記録媒体 | |
JP5905856B2 (ja) | 超音波検査装置 | |
US11439368B2 (en) | Acoustic wave processing device, signal processing method for acoustic wave processing device, and program | |
JP6000197B2 (ja) | 超音波診断装置、超音波画像生成方法およびプログラム | |
JP5905808B2 (ja) | 超音波検査装置、超音波画像データ生成方法およびプログラム | |
JP6074299B2 (ja) | 超音波診断装置、超音波診断装置の信号処理方法およびプログラム | |
US10231711B2 (en) | Acoustic wave processing device, signal processing method for acoustic wave processing device, and program | |
JP6000196B2 (ja) | 超音波診断装置、音速決定方法およびプログラム | |
US10980515B2 (en) | Acoustic wave processing apparatus, signal processing method, and program for acoustic wave processing apparatus | |
JP5964774B2 (ja) | 超音波診断装置、超音波診断装置の信号処理方法およびプログラム | |
JP5873412B2 (ja) | 超音波診断装置、音速決定方法およびプログラム | |
WO2014192466A1 (ja) | 超音波診断装置、超音波診断装置の音線信号生成方法、及び超音波診断装置の音線信号生成プログラム | |
JP6047041B2 (ja) | 超音波診断装置、超音波診断装置の信号処理方法およびプログラム | |
WO2014050847A1 (ja) | 超音波診断装置、超音波画像データ生成方法およびプログラム | |
WO2014050756A1 (ja) | 超音波検査装置、超音波画像データ生成方法およびプログラム | |
WO2014050897A1 (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: 13842757 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
122 | Ep: pct application non-entry in european phase |
Ref document number: 13842757 Country of ref document: EP Kind code of ref document: A1 |