WO2015008534A1 - 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム - Google Patents

超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム Download PDF

Info

Publication number
WO2015008534A1
WO2015008534A1 PCT/JP2014/063395 JP2014063395W WO2015008534A1 WO 2015008534 A1 WO2015008534 A1 WO 2015008534A1 JP 2014063395 W JP2014063395 W JP 2014063395W WO 2015008534 A1 WO2015008534 A1 WO 2015008534A1
Authority
WO
WIPO (PCT)
Prior art keywords
ultrasonic
unit
transmission
frequency
correction
Prior art date
Application number
PCT/JP2014/063395
Other languages
English (en)
French (fr)
Inventor
弘孝 江田
Original Assignee
オリンパスメディカルシステムズ株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by オリンパスメディカルシステムズ株式会社 filed Critical オリンパスメディカルシステムズ株式会社
Priority to JP2015503710A priority Critical patent/JP5788624B2/ja
Priority to CN201480028486.8A priority patent/CN105392427B/zh
Priority to EP14826109.2A priority patent/EP3023058B1/en
Publication of WO2015008534A1 publication Critical patent/WO2015008534A1/ja
Priority to US14/942,163 priority patent/US9662090B2/en

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/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/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/4477—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device using several separate ultrasound transducers or probes
    • 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/54—Control of the diagnostic device
    • 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/895—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum
    • 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/52033—Gain control of receivers
    • 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/5205—Means for monitoring or calibrating
    • 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/4438—Means for identifying the diagnostic device, e.g. barcodes
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B06—GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS IN GENERAL
    • B06B—METHODS OR APPARATUS FOR GENERATING OR TRANSMITTING MECHANICAL VIBRATIONS OF INFRASONIC, SONIC, OR ULTRASONIC FREQUENCY, e.g. FOR PERFORMING MECHANICAL WORK IN GENERAL
    • B06B1/00—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency
    • B06B1/02—Methods or apparatus for generating mechanical vibrations of infrasonic, sonic, or ultrasonic frequency making use of electrical energy
    • B06B1/0207—Driving circuits

Definitions

  • the present invention relates to an ultrasonic observation apparatus that observes a tissue of a specimen using ultrasonic waves, an operation method of the ultrasonic observation apparatus, and an operation program of the ultrasonic observation apparatus.
  • an ultrasound observation apparatus that generates an ultrasound image of a specimen based on ultrasound
  • it is desirable that the image displayed when quantitatively evaluating the ultrasound image is not affected by the type of ultrasound probe that transmits and receives ultrasound.
  • vibration that has been stored in advance when an image is generated using a received signal in order to suppress variation of each transducer.
  • a technique for correcting a received signal based on a frequency characteristic of a child is disclosed (for example, see Patent Document 1).
  • a technique for monitoring the transmission drive signal and correcting the transmission drive signal based on the ideal waveform data stored in advance in order to suppress variation of each vibrator is disclosed (for example, see Patent Document 2). reference).
  • the present invention has been made in view of the above, and an ultrasonic observation apparatus and an ultrasonic observation apparatus capable of realizing observation of an ultrasonic image that eliminates the influence of the difference in the type of ultrasonic probe. It is an object to provide an operation method and an operation program for an ultrasonic observation apparatus.
  • an ultrasonic observation apparatus converts an electrical transmission drive wave into an ultrasonic transmission echo, and transmits the transmission echo to a specimen.
  • an ultrasonic probe that receives a reception echo reflected by the specimen and converts it into an electrical reception signal, and a parameter that gives characteristics according to the type of the ultrasonic probe, the transmission drive
  • a storage unit that stores parameters necessary for wave generation and correction of the received signal, and the parameters stored in the storage unit are referred to, and the frequency spectrum of the transmission echo in a predetermined frequency band is the ultrasonic probe.
  • the transmission drive wave generation unit that generates the transmission drive wave according to the type of the ultrasound probe, and the parameter stored in the storage unit, so that the same regardless of the type of A reception signal correction unit that corrects the received signal according to the type of the ultrasound probe, an image processing unit that generates image data using the reception signal corrected by the reception signal correction unit, and It is provided with.
  • the transmission drive wave generator generates a frequency spectrum having a larger value as the attenuation amount is larger when converted into the transmission echo by the ultrasonic probe.
  • the transmission drive wave is generated, and the reception signal correction unit corrects the reception signal for each frequency based on the attenuation amount.
  • the storage unit calculates the attenuation amount for each frequency in the frequency spectrum when the transmission drive wave or the reception echo is converted by the ultrasonic probe.
  • the transmission drive wave generation unit stores the parameter in association with the type of the ultrasonic probe, and the transmission drive wave generation unit stores the frequency spectrum of a predetermined frequency band to be output as the transmission echo.
  • the transmission drive wave is generated by adding an attenuation amount for each frequency
  • the reception signal correction unit stores the attenuation stored in the storage unit with respect to the frequency spectrum of the reception signal received from the ultrasonic probe.
  • the received signal is corrected by adding a quantity for each frequency.
  • the ultrasonic observation apparatus is characterized in that, in the above invention, the ultrasonic probe can be selected from a plurality of ultrasonic probes of different types.
  • a frequency analysis unit that calculates a frequency spectrum by analyzing a frequency of an ultrasonic wave received by the ultrasonic probe, and a frequency calculated by the frequency analysis unit.
  • a feature amount extraction unit that extracts at least one feature amount from the frequency spectrum by approximating the spectrum, and the image processing unit is a feature amount image corresponding to the feature amount extracted by the feature amount extraction unit
  • a feature amount image data generation unit for generating data is provided.
  • the feature amount extraction unit contributes to attenuation generated according to an ultrasonic reception depth and frequency before or after performing the frequency spectrum approximation processing. Attenuation correction processing is performed to reduce the above.
  • the ultrasonic observation apparatus is characterized in that, in the above invention, the feature amount extraction unit approximates a frequency spectrum to be approximated by a polynomial by regression analysis.
  • the feature amount extraction unit approximates the frequency spectrum to be approximated by a linear expression, and the inclination of the linear expression, the intercept of the linear expression, and the inclination At least one of intensities determined using the intercept and a specific frequency included in a frequency region of the frequency spectrum is extracted as a feature amount.
  • the operation method of the ultrasonic observation apparatus converts an electrical transmission drive wave into an ultrasonic transmission echo, transmits the transmission echo to the specimen, and receives the reception echo reflected by the specimen.
  • a parameter that gives characteristics according to the type of ultrasonic probe that is received and converted into an electrical reception signal, and that stores parameters necessary for generation of the transmission drive wave and correction of the reception signal With reference to the parameter, referring to the transmission drive wave generation step of generating the transmission drive wave according to the type of the ultrasound probe by a transmission drive wave generation unit, the parameter stored in the storage unit, A received signal correction step for correcting the received signal according to the type of the ultrasound probe, and image data using the received signal corrected in the received signal correction step.
  • the operation program of the ultrasonic observation apparatus converts an electrical transmission drive wave into an ultrasonic transmission echo, transmits the transmission echo to a specimen, and receives a reception echo reflected by the specimen.
  • a parameter that gives characteristics according to the type of ultrasonic probe that is received and converted into an electrical reception signal, and that stores parameters necessary for generation of the transmission drive wave and correction of the reception signal With reference to the parameter, referring to the transmission drive wave generation step of generating the transmission drive wave according to the type of the ultrasound probe by a transmission drive wave generation unit, the parameter stored in the storage unit, A received signal correction step for correcting the received signal in accordance with the type of the ultrasound probe, and image reception using the received signal corrected in the received signal correction step. Characterized in that to execute an image processing step of generating a data, to the computer.
  • FIG. 1 is a block diagram showing a configuration of an ultrasonic observation apparatus according to Embodiment 1 of the present invention.
  • FIG. 2 is a diagram illustrating a change in the frequency spectrum when the signal conversion unit converts the transmission drive wave into a transmission echo.
  • FIG. 3 is a diagram illustrating a transmission drive wave generated by the transmission drive wave generation unit and a case where the signal conversion unit of the ultrasonic probe converts the transmission drive wave into a transmission echo.
  • FIG. 4 is a diagram illustrating an outline (first example) of correction processing performed by the reception signal correction unit.
  • FIG. 5 is a diagram showing an outline (second example) of correction processing performed by the reception signal correction unit.
  • FIG. 1 is a block diagram showing a configuration of an ultrasonic observation apparatus according to Embodiment 1 of the present invention.
  • FIG. 2 is a diagram illustrating a change in the frequency spectrum when the signal conversion unit converts the transmission drive wave into a transmission echo.
  • FIG. 3 is a diagram illustrating a transmission
  • FIG. 6 is a diagram illustrating a relationship between the reception depth and the amplification factor in the amplification process performed by the STC correction unit on the reception signal corrected by the reception signal correction unit.
  • FIG. 7 is a flowchart showing an outline of processing of the ultrasonic observation apparatus according to Embodiment 1 of the present invention.
  • FIG. 8 is a block diagram showing a configuration of an ultrasonic observation apparatus according to Embodiment 2 of the present invention.
  • FIG. 9 is a diagram illustrating the relationship between the reception depth and the amplification factor in the amplification process performed by the amplification correction unit.
  • FIG. 10 is a diagram illustrating an example of a frequency spectrum calculated by the frequency analysis unit.
  • FIG. 11 is a diagram illustrating a frequency spectrum (first example) before and after correction by the reception signal correction unit when an ultrasonic probe is used.
  • FIG. 12 is a diagram illustrating a frequency spectrum (second example) before and after correction by the reception signal correction unit when an ultrasonic probe is used.
  • FIG. 13 is a diagram in which the feature amounts of the frequency spectrum curves shown in FIGS. 11 and 12 are plotted in the feature amount space.
  • FIG. 14 is a diagram illustrating a straight line corresponding to the feature amount corrected by the attenuation correction unit.
  • FIG. 15 is a flowchart showing an outline of processing of the ultrasonic observation apparatus according to Embodiment 2 of the present invention.
  • FIG. 16 is a flowchart illustrating an outline of processing performed by the frequency analysis unit.
  • FIG. 17 is a diagram schematically showing a data array of one sound ray.
  • FIG. 18 is a diagram schematically illustrating an outline of the attenuation correction process performed by the attenuation correction unit of the ultrasonic observation apparatus according to another embodiment of the present invention.
  • FIG. 1 is a block diagram showing a configuration of an ultrasonic observation apparatus according to Embodiment 1 of the present invention.
  • An ultrasonic observation apparatus 1 shown in FIG. 1 is an apparatus for observing a specimen that is a diagnosis target using ultrasonic waves.
  • the ultrasonic observation apparatus 1 outputs an ultrasonic transmission echo to the outside, receives an ultrasonic reception echo reflected from the outside, and converts it into an electrical reception signal;
  • the transmitter / receiver 3 transmits / receives an electrical signal to / from the acoustic probe 2, the image processor 4 generates image data corresponding to the received signal, and an interface such as a keyboard, a mouse, and a touch panel.
  • An input unit 5 that receives input of various types of information, a display unit 6 that displays various types of information including an image generated by the image processing unit 4 and that is realized by using a display panel made of liquid crystal, organic EL, or the like;
  • a storage unit 7 that stores various kinds of information necessary for ultrasonic observation including parameters indicating characteristics according to the type of the touch element 2 and a control unit 8 that controls the operation of the ultrasonic observation apparatus 1 are provided.
  • the ultrasonic observation apparatus 1 includes a scope in which an ultrasonic probe 2 is provided at a distal end portion thereof, and a processing apparatus (processor) in which a base end of the scope is detachably connected and the above-described parts other than the ultrasonic probe 2 are provided. ).
  • the ultrasonic probe 2 converts the transmission drive wave, which is an electrical pulse signal received from the transmission / reception unit 3, into a transmission echo, which is an acoustic pulse signal, and the received echo that has been reflected and returned by an external specimen. It has a signal converter 21 for converting it into an electrical reception signal.
  • the signal conversion unit 21 functions as a filter.
  • FIG. 2 is a diagram illustrating a change in frequency spectrum when the signal conversion unit 21 converts a transmission drive wave into a transmission echo.
  • the horizontal axis f represents frequency and the vertical axis I represents intensity.
  • F ⁇ f
  • FIGS. 2B and 2C are diagrams showing the spectrums of transmission echoes generated by the signal converters 21 of the ultrasonic probes 2 of different types.
  • the spectrum C 2 shown in FIG. 2B is more attenuated than the spectrum C 3 shown in FIG.
  • the amount of attenuation ⁇ I 1 (f 1 ) corresponding to the difference between the spectrum C 1 and the spectrum C 2 at an arbitrary frequency f 1 (f L ⁇ f 1 ⁇ f H ) in the frequency band F is
  • the ultrasonic probe information is stored in the ultrasonic probe information storage unit 71 of the storage unit 7 in association with the type of the acoustic probe 2.
  • the attenuation amount ⁇ I 2 (f 1 ) corresponding to the difference between the spectrum C 1 and the spectrum C 3 at an arbitrary frequency f 1 in the frequency band F is also the kind of the ultrasound probe 2.
  • the ultrasound probe information storage unit 71 is also the kind of the ultrasound probe 2.
  • two ultrasonic probes 2 have been described as examples, but other types of ultrasonic probes 2 may be stored in the ultrasonic probe information storage unit 71 in the same manner. Is possible.
  • the ultrasonic probe 2 for outputting a spectrum C 2 shown in FIG. 2 (b) as a transmission echo ultrasonic probe 2A
  • the ultrasonic probe 2 to the spectrum C 3 is output as transmission echoes shown in FIG. 2 (c) of ultrasonic probe 2B.
  • the signal conversion unit 21 of the ultrasonic probe 2A is referred to as a signal conversion unit 21A
  • the signal conversion unit 21 of the ultrasonic probe 2B is referred to as a signal conversion unit 21B.
  • the ultrasonic probe 2 may be one that mechanically scans an ultrasonic transducer, or one that electronically scans a plurality of ultrasonic transducers. In the first embodiment, it is possible to select and use any one of a plurality of different types of ultrasound probes 2 as the ultrasound probe 2.
  • the transmitting / receiving unit 3 is electrically connected to the ultrasonic probe 2.
  • the transmission / reception unit 3 generates a transmission drive wave based on a preset waveform and transmission timing, and transmits the generated transmission drive wave to the ultrasonic probe 2.
  • a received signal correction unit 32 that corrects the received signal received from the probe 2 to eliminate the influence of the difference in characteristics depending on the type of the ultrasound probe 2, and the received signal corrected by the received signal correction unit 32
  • the STC correction unit 33 performs STC (Sensitivity Time Control) correction that corrects at a higher amplification factor as the reception depth increases.
  • the transmission drive wave generation unit 31 generates a waveform based on the characteristics of the ultrasound probe 2 as a transmission drive wave based on information stored in an ultrasound probe information storage unit 71 of the storage unit 7 to be described later.
  • FIG. 3 is a diagram illustrating a transmission drive wave generated by the transmission drive wave generation unit 31 and a case where the signal conversion unit 21 of the ultrasound probe 2 converts the transmission drive wave into a transmission echo.
  • the transmission drive wave generation unit 31 generates different transmission drive waves for each type of the ultrasound probe 2. Specifically, the transmission drive wave generation unit 31 reads out the characteristics corresponding to the type of the ultrasound probe 2 from the ultrasound probe information storage unit 71, and has a rectangular shape as shown in FIG. A transmission drive wave having a spectrum C 6 is generated.
  • the characteristics of the ultrasonic probe 2 here are the attenuation amounts ⁇ I 1 (f) and ⁇ I 2 (f) (where f ⁇ F) and the ultrasonic probe shown in FIGS. 2B and 2C, respectively. This means information unique to each ultrasonic probe 2 including the type of the probe 2.
  • the curve and the straight line are composed of a set of discrete points.
  • the transmission drive wave generator 31 generates a transmission drive wave corresponding to the type of the ultrasound probe 2 based on the spectrum of the transmission echo to be generated and the characteristics of the ultrasound probe 2.
  • the transmission drive wave generation unit 31 generates a transmission drive wave having a frequency spectrum in a predetermined frequency band of the transmission echo that is different for each type of the ultrasonic probe.
  • a spectrum C 4 of the transmission drive wave shown in FIG. 3A is for the ultrasonic probe 2A
  • a spectrum C 5 of the transmission drive wave shown in FIG. 3B is for the ultrasonic probe 2B.
  • the spectrum C 4 has a shape obtained by adding an attenuation amount ⁇ I 1 (f) (f ⁇ F) when the ultrasound probe 2A converts the spectrum C 6 of the transmission echo.
  • the spectrum C 5 has a shape obtained by adding the attenuation ⁇ I 2 (f) (f ⁇ F) when the ultrasound probe 2B converts the spectrum C 6 of the transmission echo.
  • the reception signal correction unit 32 applies correction corresponding to the type of the ultrasonic probe to the reception signal received from the ultrasonic probe 2 by the transmission / reception unit 3.
  • FIG. 4 is a diagram showing an outline of the correction process performed by the reception signal correction unit 32, and specifically shows an outline of the correction process for the reception signal sent from the ultrasonic probe 2A.
  • Spectrum C 7 shown in FIG. 4 (a) is the frequency spectrum of the received signal the signal conversion unit 21A of the ultrasonic probe 2A is generated by converting the received echo.
  • the reception signal correction unit 32 corrects the reception signal by adding attenuation ⁇ I 1 (f) (f ⁇ F) when the ultrasound probe 2A converts the spectrum C 7 .
  • a spectrum C 8 shown in FIG. 4B is obtained.
  • the spectrum C 8 is a spectrum that excludes the influence of the ultrasonic probe 2A.
  • FIG. 5 is a diagram showing an outline of the correction process performed by the reception signal correction unit 32, and specifically shows an outline of the correction process for the reception signal sent from the ultrasonic probe 2B.
  • a spectrum C 9 shown in FIG. 5A is a frequency spectrum of a reception signal generated by converting a reception echo by the signal conversion unit 21B of the ultrasonic probe 2B.
  • the reception signal correction unit 32 corrects the reception signal by adding the amount of attenuation ⁇ I 2 (f) (f ⁇ F) when the ultrasound probe 2B converts the spectrum C 9 .
  • a spectrum C 8 shown in FIG. 5B is obtained.
  • This spectrum C 8 is a spectrum that excludes the influence of the ultrasound probe 2B, and is the same as the spectrum obtained as a result of correcting the received signal converted by the ultrasound probe 2A (FIG. 4). (See (b)). In other words, the spectrum C 8 is nothing but the spectrum of the received echo.
  • FIG. 6 is a diagram illustrating a relationship between the reception depth and the amplification factor in the amplification process performed by the STC correction unit 33 on the reception signal corrected by the reception signal correction unit 32.
  • the reception depth z shown in FIG. 6 is an amount calculated based on the elapsed time from the reception start time of the ultrasonic wave.
  • the amplification factor ⁇ (dB) increases linearly from ⁇ 0 to ⁇ th (> ⁇ 0 ) as the reception depth z increases.
  • the amplification factor ⁇ takes a constant value ⁇ th when the reception depth z is equal to or greater than the threshold value z th .
  • the value of the threshold value z th is such a value that the ultrasonic signal received from the specimen is almost attenuated and the noise becomes dominant. More generally, when the reception depth z is smaller than the threshold value z th , the amplification factor ⁇ may increase monotonously as the reception depth z increases.
  • the transmission / reception unit 3 performs processing such as filtering on the echo signal amplified by the STC correction unit 33, and then performs A / D conversion to generate and output a digital RF signal in the time domain.
  • processing such as filtering on the echo signal amplified by the STC correction unit 33
  • a / D conversion to generate and output a digital RF signal in the time domain.
  • the transmission / reception unit 3 has a multi-channel circuit for beam synthesis corresponding to the plurality of ultrasonic transducers.
  • the image processing unit 4 includes a B-mode image data generation unit 41 that generates B-mode image data from the received signal.
  • the B-mode image data generation unit 41 performs signal processing using a known technique such as a bandpass filter, logarithmic conversion, gain processing, contrast processing, and the like on the digital signal, and also according to the image display range on the display unit 6.
  • B-mode image data is generated by thinning out data in accordance with the data step width determined in advance.
  • the storage unit 7 includes an ultrasound probe information storage unit 71, an amplification factor information storage unit 72, and a window function storage unit 73.
  • the ultrasound probe information storage unit 71 converts the type of the ultrasound probe 2 and the ultrasound probe 2 converts the transmission drive wave into a transmission echo or converts the reception echo into a reception signal.
  • the correlation between the amount of attenuation generated in the spectrum during conversion or the like is stored as a parameter.
  • the amplification factor information storage unit 72 stores, as amplification factor information, the relationship between the amplification factor and the reception depth that is referred to when the STC correction unit 33 performs amplification processing (for example, the relationship illustrated in FIG. 2).
  • the window function storage unit 73 stores at least one of window functions such as Hamming, Hanning, and Blackman.
  • the storage unit 7 is realized using a ROM in which an operation program of the ultrasound observation apparatus 1, a program for starting a predetermined OS, and the like are stored in advance, and a RAM in which calculation parameters and data of each process are stored.
  • the control unit 8 includes an ultrasonic probe determination unit 81 that determines the type of the connected ultrasonic probe 2.
  • the determination result of the ultrasonic probe determination unit 81 is stored in the ultrasonic probe information storage unit 71, and is referred to when the signal conversion unit 21, the transmission drive wave generation unit 31, and the reception signal correction unit 32 perform processing. Is done.
  • the scope including the ultrasonic probe 2 determines the type of the ultrasonic probe 2 as a processing device at the connection portion with the processing device at the subsequent stage. It is only necessary to provide a connection pin for this purpose. Thereby, the ultrasonic probe determination unit 81 provided on the processing apparatus side can determine the type of the ultrasonic probe 2 according to the shape of the connection pin of the connected scope.
  • the control unit 8 is realized by using a CPU having calculation and control functions.
  • the control unit 8 reads out various programs including information stored and stored in the storage unit 7 and the operation program of the ultrasonic observation apparatus 1 from the storage unit 7, thereby performing various arithmetic processes related to the operation method of the ultrasonic observation apparatus 1. To control the ultrasonic observation apparatus 1 in an integrated manner.
  • the operation program of the ultrasonic observation apparatus 1 can be recorded on a computer-readable recording medium such as a hard disk, flash memory, CD-ROM, DVD-ROM, or flexible disk and widely distributed. Recording of various programs on a recording medium or the like may be performed when the computer or the recording medium is shipped as a product, or may be performed by downloading via a communication network.
  • a computer-readable recording medium such as a hard disk, flash memory, CD-ROM, DVD-ROM, or flexible disk and widely distributed. Recording of various programs on a recording medium or the like may be performed when the computer or the recording medium is shipped as a product, or may be performed by downloading via a communication network.
  • FIG. 7 is a flowchart showing an outline of processing of the ultrasonic observation apparatus 1 having the above configuration. It is assumed that the type of the ultrasound probe 2 provided in the ultrasound observation apparatus 1 is determined in advance by the ultrasound probe determination unit 81.
  • the transmission drive wave generation unit 31 generates a transmission drive wave and outputs it to the ultrasound probe 2 (step S1).
  • the spectrum of the transmission drive wave is a spectrum C 4 shown in FIG.
  • the spectrum of the transmitted driving wave is a spectrum C 5 shown in FIG. 3 (b).
  • the signal converter 21 converts the transmission drive wave into a transmission echo, and generates and outputs the converted transmission echo (step S2).
  • the spectrum of this transmission echo is, for example, the spectrum C 6 in FIG. As described with reference to FIG. 3, the transmission echo has the same spectrum regardless of the type of the ultrasound probe 2.
  • the ultrasound probe 2 receives the reception echo returned from the transmission echo reflected by the living body (step S3).
  • the signal conversion unit 21 that has received the reception echo converts the reception echo into a reception signal and outputs it to the transmission / reception unit 3 (step S4).
  • the spectrum of this received signal is, for example, the spectrum C 7 shown in FIG. 4A in the case of the ultrasonic probe 2A. Also, when the ultrasonic probe 2B, the spectrum of the received signal is a spectrum C 9 shown in Figure 5 (a).
  • the reception signal correction unit 32 that has received the reception signal from the ultrasound probe 2 corrects the frequency spectrum (step S5).
  • the spectrum of this received signal is the same spectrum regardless of the type of the ultrasound probe 2 (see the spectrum C 8 in FIGS. 4B and 5B).
  • the STC correction unit 33 performs STC correction on the reception signal whose spectrum is corrected by the reception signal correction unit 32 (step S6).
  • the STC correction unit 33 performs STC correction based on the relationship between the amplification factor and the reception depth illustrated in FIG. 6, for example.
  • the B-mode image data generation unit 41 generates B-mode image data using the echo signal amplified by the STC correction unit 33 (step S7).
  • the display unit 6 displays a B mode image corresponding to the B mode image generated by the B mode image data generation unit 41 (step S8).
  • step S8 the ultrasound observation apparatus 1 ends a series of processes. Note that the ultrasonic observation apparatus 1 may periodically repeat the processes of steps S1 to S8.
  • the processing for eliminating the influence according to the type of the ultrasonic probe is performed on both the transmission drive wave and the reception signal, the ultrasonic wave It is possible to realize observation of an ultrasonic image that eliminates the influence of the difference in the type of probe.
  • the transmission driving wave having a frequency spectrum having a larger value as the attenuation amount when converted into a transmission echo by the ultrasonic probe is larger is generated, and the received signal is generated. Since the correction is performed for each frequency based on the amount of attenuation, it is possible to reliably eliminate the influence of the difference in characteristics depending on the type of the ultrasonic probe.
  • FIG. 8 is a block diagram showing a configuration of an ultrasonic observation apparatus according to Embodiment 2 of the present invention.
  • the ultrasonic observation apparatus 11 shown in the figure includes a predetermined calculation for an ultrasonic probe 2, a transmission / reception unit 3, an input unit 5, a display unit 6, a control unit 8, and an electrical reception signal.
  • symbol as the ultrasonic observation apparatus 1 is attached
  • the arithmetic unit 12 performs an amplification correction unit 121 for performing an amplification correction to make the amplification factor constant regardless of the reception depth with respect to the digital RF signal output from the transmission / reception unit 3, and a fast Fourier transform to the digital RF signal subjected to the amplification correction
  • Frequency analysis unit 122 that calculates a frequency spectrum by performing (FFT) and performing frequency analysis, and approximation processing based on regression analysis and ultrasonic waves propagate to the frequency spectrum of each location calculated by frequency analysis unit 122
  • a feature amount extraction unit 123 that extracts the feature amount of the specimen by performing attenuation correction processing that reduces the contribution of attenuation generated according to the reception depth and frequency of the ultrasonic wave.
  • FIG. 9 is a diagram illustrating the relationship between the reception depth and the amplification factor in the amplification process performed by the amplification correction unit 121.
  • the amplification rate ⁇ (dB) in the amplification process performed by the amplification correction unit 121 takes the maximum value ⁇ th ⁇ 0 when the reception depth z is zero, and the reception depth z is zero to the threshold value z th. Decreases linearly until reaching 0 and is zero when the reception depth z is greater than or equal to the threshold z th .
  • the amplification correction unit 121 amplifies and corrects the digital RF signal with the amplification factor determined in this way, thereby canceling the influence of the STC correction in the STC correction unit 33 and outputting a signal with a constant amplification rate ⁇ th. .
  • the relationship between the reception depth z and the amplification factor ⁇ performed by the amplification correction unit 121 is different depending on the relationship between the reception depth and the amplification factor in the STC correction unit 33.
  • the STC correction is a correction that uniformly amplifies the amplitude of the analog signal waveform over the entire frequency band. For this reason, when generating a B-mode image using the amplitude of ultrasonic waves, a sufficient effect can be obtained by performing STC correction, while in the case of calculating the frequency spectrum of ultrasonic waves. However, the influence of attenuation associated with the propagation of ultrasonic waves cannot be accurately eliminated. To solve this problem, when generating a B-mode image, output a reception signal subjected to STC correction, while generating an image based on a frequency spectrum, It is conceivable to perform a different new transmission and output a received signal not subjected to STC correction.
  • the amplification correction unit 121 corrects the amplification factor.
  • the frequency analysis unit 122 calculates a frequency spectrum at a plurality of locations (data positions) on the sound ray by performing fast Fourier transform on an FFT data group having a predetermined data amount for each sound ray (line data). The calculation result by the frequency analysis unit 122 is obtained as a complex number and stored in the storage unit 14.
  • the frequency spectrum shows different tendencies depending on the tissue properties of the specimen. This is because the frequency spectrum has a correlation with the size, density, acoustic impedance, and the like of the specimen as a scatterer that scatters ultrasonic waves.
  • the “tissue property” is, for example, any of cancer, endocrine tumor, mucinous tumor, normal tissue, vessel, and the like.
  • FIG. 10 is a diagram illustrating an example of a frequency spectrum calculated by the frequency analysis unit 122. Specifically, a frequency spectrum obtained by performing a fast Fourier transform on the FFT data group is expressed by intensity I (f, z) and phase ⁇ (f, z) with frequency f and reception depth as a function of z. The spectrum of intensity I (f, z) is shown.
  • intensity refers to any of parameters such as voltage, power, sound pressure, and acoustic energy.
  • the horizontal axis f is the frequency
  • the vertical axis I is the intensity
  • the reception depth z is constant.
  • the feature amount extraction unit 123 depends on the ultrasonic reception depth and frequency with respect to the approximation unit 124 that calculates the approximate expression of the frequency spectrum calculated by the frequency analysis unit 122 by regression analysis, and the approximate expression calculated by the approximation unit 124.
  • an attenuation correction unit 125 that extracts a characteristic amount of the frequency spectrum by performing attenuation correction processing that reduces the contribution of attenuation of ultrasonic waves.
  • the approximating unit 124 extracts a feature amount before attenuation correction (hereinafter referred to as a pre-correction feature amount) that characterizes the approximated primary equation by approximating the frequency spectrum with a linear equation (regression line) by regression analysis. Specifically, the approximating unit 124 extracts the slope a 0 and the intercept b 0 of the linear expression as the pre-correction feature quantity.
  • a straight line L 10 illustrated in FIG. 10 is a straight line corresponding to a linear expression approximated by the approximation unit 124.
  • the inclination a 0 has a correlation with the size of the ultrasonic scatterer, and it is generally considered that the larger the scatterer, the smaller the inclination.
  • the intercept b 0 has a correlation with the size of the scatterer, the difference in acoustic impedance, the density (concentration) of the scatterer, and the like. Specifically, it is considered that the intercept b 0 has a larger value as the scatterer is larger, a larger value as the acoustic impedance is larger, and a larger value as the density (concentration) of the scatterer is larger.
  • the intensity at the center frequency f M (hereinafter simply referred to as “intensity”) c 0 is an indirect parameter derived from the slope a 0 and the intercept b 0 , and gives the spectrum intensity at the center in the effective frequency band. Therefore, the intensity c 0 is considered to have a certain degree of correlation with the brightness of the B-mode image in addition to the size of the scatterer, the difference in acoustic impedance, and the density of the scatterer.
  • the approximate polynomial calculated by the feature quantity extraction unit 123 is not limited to a linear expression, and it is possible to use a quadratic or higher approximate polynomial.
  • the value of the pre-correction feature value is the same regardless of the type of the ultrasound probe 2.
  • the frequency spectrum referred to in the first embodiment is used to describe the case where two different types of ultrasonic probes 2A and 2B are used.
  • description will be made using the intercept b 0 and the intensity c 0 as the pre-correction feature quantity.
  • FIG. 11 is a diagram illustrating frequency spectra before and after correction by the reception signal correction unit 32 when the ultrasound probe 2A is used. Specifically, spectrum C 7 is before reception signal correction, and spectrum C 8 is after reception signal correction.
  • FIG. 12 is a diagram illustrating frequency spectra before and after correction by the reception signal correction unit 32 when the ultrasonic probe 2B is used. Specifically, spectrum C 9 is before reception signal correction, and spectrum C 8 is after reception signal correction.
  • FIG. 13 is a diagram in which pairs of feature amounts of the spectra C 7 to C 9 shown in FIGS. 11 and 12 are plotted in a feature amount space (b 0 , c 0 ), respectively.
  • a certain point B (b 2 , c 2 ) is generally a different point on the feature amount space.
  • the set of corrected feature values (b 0 , c 0 ) is the spectrum before correction. Regardless, the point is C (b 3 , c 3 ).
  • the pre-correction feature amounts of the frequency spectrum of the received signal are identical because the frequency spectrum depending on the type of the ultrasound probe 2 when the transmission drive wave generation unit 31 generates the transmission drive wave.
  • the reception signal correction unit 32 eliminates the influence of the frequency spectrum depending on the type of the ultrasound probe 2 on the reception signal.
  • the frequency spectrum of the received signal differs depending on the type of the ultrasound probe 2. Therefore, according to the second embodiment, it is possible to reliably eliminate the influence due to the difference in the type of the ultrasound probe 2, and it is possible to realize more accurate feature extraction.
  • ⁇ is the attenuation rate
  • z is the ultrasonic reception depth
  • f is the frequency.
  • the attenuation amount A (f, z) is proportional to the frequency f.
  • the value of the attenuation rate ⁇ may be set or changed by input from the input unit 5.
  • FIG. 14 is a diagram illustrating a straight line corresponding to the feature amount corrected by the attenuation correction unit 125.
  • the image processing unit 13 generates feature amount image data for displaying information corresponding to the feature amount extracted by the B-mode image data generation unit 41 and the feature amount extraction unit 123 according to one of a plurality of display methods.
  • An image data generation unit 131 An image data generation unit 131.
  • the feature amount image data generation unit 131 determines the information assigned to each pixel in the feature amount image data according to the data amount of the FFT data group when the frequency analysis unit 122 calculates the frequency spectrum. Specifically, for example, information corresponding to the feature amount of the frequency spectrum calculated from the FFT data group is assigned to the pixel region corresponding to the data amount of one FFT data group. In the second embodiment, the number of feature amounts used when generating feature amount image data can be arbitrarily set.
  • the storage unit 14 includes a correction information storage unit 141 in addition to the ultrasound probe information storage unit 71, the amplification factor information storage unit 72, and the window function storage unit 73.
  • the amplification factor information storage unit 72 is a relationship between the amplification factor and the reception depth that is referred to when the STC correction unit 33 performs the amplification process (for example, the relationship illustrated in FIG. 6), and when the amplification correction unit 121 performs the amplification correction process.
  • the relationship (for example, the relationship shown in FIG. 9) between the amplification factor and the reception depth referred to respectively is stored as amplification factor information.
  • the correction information storage unit 141 stores information related to attenuation correction including Expression (1).
  • the storage unit 14 is realized by using a ROM in which an operation program of the ultrasonic observation apparatus 11, a program for starting a predetermined OS, and the like are stored in advance, and a RAM in which calculation parameters and data of each process are stored.
  • FIG. 15 is a flowchart showing an outline of processing of the ultrasonic observation apparatus 11 having the above configuration.
  • the processes in steps S11 to S17 correspond to steps S1 to S7 in the flowchart shown in FIG. Hereinafter, the process after step S18 is demonstrated.
  • step S18 the amplification correction unit 121 performs amplification correction on the signal output from the transmission / reception unit 3 so that the amplification factor is constant regardless of the reception depth (step S18).
  • the amplification correction unit 121 performs amplification correction based on, for example, the relationship between the amplification factor and the reception depth shown in FIG.
  • the frequency analysis unit 122 calculates a frequency spectrum by performing frequency analysis by FFT calculation (step S19).
  • step S19 the process (step S19) performed by the frequency analysis unit 122 will be described in detail with reference to the flowchart shown in FIG.
  • the frequency analysis unit 122 sets a counter k for identifying a sound ray to be analyzed as k 0 (step S31).
  • FIG. 17 is a diagram schematically showing a data array of one sound ray.
  • a white or black rectangle means one piece of data.
  • the sound ray SR k is discretized at a time interval corresponding to a sampling frequency (for example, 50 MHz) in A / D conversion performed by the transmission / reception unit 3.
  • FIG. 17 shows the case where the first data position of the sound ray SR k is set as the initial value Z (k) 0 , the position of the initial value can be arbitrarily set.
  • the frequency analysis unit 122 acquires the FFT data group at the data position Z (k) (step S33), and applies the window function stored in the window function storage unit 73 to the acquired FFT data group (step S34). ). In this way, by applying the window function to the FFT data group, it is possible to avoid the FFT data group from becoming discontinuous at the boundary and to prevent the occurrence of artifacts.
  • the frequency analysis unit 122 determines whether or not the FFT data group at the data position Z (k) is a normal data group (step S35).
  • the FFT data group needs to have a power number of 2 data.
  • the number of data in the FFT data group is 2 n (n is a positive integer).
  • the normal FFT data group means that the data position Z (k) is the 2 n-1 th position from the front in the FFT data group.
  • the FFT data groups F 2 and F 3 are both normal.
  • step S35 If the result of determination in step S35 is that the FFT data group at the data position Z (k) is normal (step S35: Yes), the frequency analysis unit 122 proceeds to step S37 described later.
  • step S35 If the result of determination in step S35 is that the FFT data group at the data position Z (k) is not normal (step S35: No), the frequency analysis unit 122 inserts zero data as much as the deficiency into the normal FFT data group. Generate (step S36).
  • the FFT function group determined to be not normal in step S35 is subjected to a window function before adding zero data. For this reason, discontinuity of data does not occur even if zero data is inserted into the FFT data group.
  • step S36 the frequency analysis unit 122 proceeds to step S37 described later.
  • step S37 the frequency analysis unit 122 obtains a frequency spectrum composed of complex numbers by performing an FFT operation using the FFT data group (step S37). As a result, for example, a spectrum C 11 as shown in FIG. 10 is obtained.
  • the frequency analysis unit 122 changes the data position Z (k) by the step width D (step S38).
  • the step width D is stored in the storage unit 14 in advance.
  • the step width D is desirably matched with the data step width used when the B-mode image data generation unit 41 generates the B-mode image data.
  • a value larger than the data step width may be set.
  • the frequency analysis unit 122 determines whether or not the data position Z (k) is larger than the maximum value Z (k) max in the sound ray SR k (step S39). If the data position Z (k) is greater than the maximum value Z (k) max (step S39: Yes), the frequency analysis section 122 is increased by one counter k (step S40). On the other hand, when the data position Z (k) is equal to or less than the maximum value Z (k) max (step S39: No), the frequency analysis unit 122 returns to step S33.
  • the frequency analysis unit 122 performs an FFT operation on [ ⁇ (Z (k) max ⁇ Z (k) 0 ) / D ⁇ +1] FFT data groups for the sound ray SR k .
  • [X] represents the maximum integer not exceeding X.
  • step S40 the frequency analysis unit 122 determines whether the counter k is greater than the maximum value k max (step S41). When the counter k is greater than k max (step S41: Yes), the frequency analysis unit 122 ends a series of FFT processing. On the other hand, when the counter k is equal to or less than k max (step S41: No), the frequency analysis unit 122 returns to step S32.
  • the frequency analysis unit 122 performs the FFT operation a plurality of times for each of (k max ⁇ k 0 +1) sound rays.
  • the frequency analysis unit 122 performs frequency analysis processing on all the areas where the ultrasonic signal is received, a setting input of a specific region of interest is received by the input unit 5 in advance, You may make it perform a frequency analysis process only within the region of interest.
  • the approximating unit 124 extracts a pre-correction feature quantity by performing regression analysis on the frequency spectrum calculated by the frequency analyzing unit 122 as an approximating process (step S20). Specifically, the approximating unit 124 calculates a linear expression that approximates the intensity I (f, z) of the frequency spectrum in the frequency spectrum frequency band f L ⁇ f ⁇ f H by regression analysis. A characteristic gradient a 0 and an intercept b 0 (, intensity c 0 ) are extracted as pre-correction feature values. Linear L 10 shown in FIG. 10, in step S20, which is an example of a regression line obtained by performing the pre-correction feature amount extraction processing for the spectrum C 11.
  • the attenuation correction unit 125 performs an attenuation correction process on the pre-correction feature quantity extracted by the approximation unit 124 (step S21).
  • the data sampling frequency is 50 MHz
  • the data sampling time interval is 20 (nsec).
  • the sound speed is 1530 (m / sec)
  • the data position Z is 0 using the data step number n and the data step width D. 0153 nD (mm).
  • the attenuation correction unit 125 substitutes the value of the data position Z obtained in this way into the reception depth z of the above-described equations (2) to (4), so that the slope a and the intercept b ( , Intensity c) is calculated.
  • a straight line L 11 shown in FIG. 14 can be given.
  • Steps S20 and S21 described above constitute a feature quantity extraction step in which the feature quantity extraction unit 123 extracts at least one feature quantity from the frequency spectrum by approximating the frequency spectrum.
  • the feature amount image data generation unit 131 generates feature amount image data using the feature amount extracted in the feature amount extraction step (steps S20 and S21) (step S22).
  • the display unit 6 displays the B mode image generated by the B mode image data generation unit 41 and / or the feature amount image generated by the feature amount image data generation unit 131 (step S23).
  • the display unit 6 may display either one of the B mode image and the feature amount image, may display the B mode image and the feature amount image side by side, or may display the B mode image and the feature amount image. May be superimposed and displayed.
  • the mixing ratio of the B-mode image and the feature amount image may be changed by an input from the input unit 5.
  • the ultrasonic observation apparatus 11 ends a series of processes. Note that the ultrasound observation apparatus 11 may periodically repeat the processes of steps S11 to S23.
  • the processing for eliminating the influence according to the type of the ultrasonic probe is performed on both the transmission drive wave and the reception signal.
  • the processing for eliminating the influence according to the type of the ultrasonic probe is performed on both the transmission drive wave and the reception signal.
  • a transmission drive wave having a frequency spectrum having a larger value as the attenuation amount when converted into a transmission echo by the ultrasonic probe is larger is generated, and the received signal is generated. Since the correction is performed for each frequency based on the attenuation amount, the influence of the difference in characteristics depending on the type of the ultrasonic probe can be surely eliminated.
  • the feature amount is extracted by performing frequency analysis, but the feature amount is excluded from the feature amount according to the type of the ultrasound probe. Therefore, according to the second embodiment, even when a quantitative evaluation of an ultrasonic image is performed, the type of ultrasonic probe is not affected.
  • the feature quantity extraction unit 123 may calculate the approximate expression of the corrected frequency spectrum after performing the attenuation correction of the frequency spectrum.
  • FIG. 18 is a diagram schematically illustrating an outline of the attenuation correction process performed by the attenuation correction unit 125. As shown in FIG. 18, the attenuation correction unit 125 attenuates the expression (1) to the intensity I (f, z) at all frequencies f (f L ⁇ f ⁇ f H ) in the band for the spectrum C 11.
  • control unit 8 may collectively perform the amplification correction process by the amplification correction unit 121 and the attenuation correction process by the attenuation correction unit 125. This process is equivalent to performing the amplification correction process in step S18 of FIG. 15 and changing the definition of the attenuation amount of the attenuation correction process in step S21 of FIG. 15 to the following equation (6).
  • a ′ 2 ⁇ zf + ⁇ (z) (6)
  • ⁇ (z) on the right side is the difference between the amplification factors ⁇ and ⁇ 0 at the reception depth z
  • ⁇ (z) ⁇ ⁇ ( ⁇ th ⁇ 0 ) / z th ⁇ z + ⁇ th ⁇ 0 (z ⁇ z th ) (7)
  • ⁇ (z) 0 (z > z th) ⁇ (8) It is expressed.

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Radiology & Medical Imaging (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Gynecology & Obstetrics (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

　超音波観測装置は、電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子と、超音波探触子の種類に応じた特性を与えるパラメータであって送信駆動波の生成および受信信号の補正に必要なパラメータを記憶する記憶部と、記憶部が記憶する前記パラメータを参照し、送信エコーの所定の周波数帯域における周波数スペクトルが超音波探触子の種類によらず同じとなるように、超音波探触子の種類に応じた送信駆動波を生成する送信駆動波生成部と、記憶部が記憶する前記パラメータを参照し、受信信号に対して超音波探触子の種類に応じた補正を行う受信信号補正部と、受信信号補正部が補正した受信信号を用いて画像データを生成する画像処理部と、を備える。

Description

超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
　本発明は、超音波を用いて検体の組織を観測する超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラムに関する。
　超音波に基づく検体の超音波画像を生成する超音波観測装置では、超音波画像を定量評価する際に表示される画像が超音波を送受信する超音波探触子の種類に影響されないことが望ましい。従来、超音波信号に基づく検体の超音波画像を生成する超音波観測装置において、各振動子のばらつきを抑制するために、受信信号を用いて画像を生成する際に、予め記憶していた振動子の周波数特性に基づいて受信信号を補正する技術が開示されている（例えば、特許文献１を参照）。また、各振動子のばらつきを抑制するために、送信駆動信号をモニタし、予め記憶している理想波形データに基づいて送信駆動信号を補正する技術が開示されている（例えば、特許文献２を参照）。
特開２０１２－２５００８０号公報 特開２００５－２７０２４７号公報
　一般に、超音波の検体内における伝播は、非線形伝播であることが知られている。このため、同一の検体に対して異なる送信エコーを送信した場合には、受信信号でいかに補正しても超音波探触子の種類の違いによる影響を排除することができない。
　また、同一の検体に対して異なる種類の超音波探触子から同一形状の送信エコーを送信した場合、検体による非線形伝播の影響は同じであるものの、超音波探触子によって変換された受信信号は同一ではないため、超音波探触子の種類の違いによる影響を排除して受信信号を生成することは難しい。
　このように、上述した従来技術のように送信側または受信側でのみ補正を行っただけでは、超音波探触子の種類の違いによる影響を排除した超音波画像の観察を実現することは難しかった。
　本発明は、上記に鑑みてなされたものであって、超音波探触子の種類の違いによる影響を排除した超音波画像の観察を実現することができる超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラムを提供することを目的とする。
　上述した課題を解決し、目的を達成するために、本発明に係る超音波観測装置は、電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、前記検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子と、前記超音波探触子の種類に応じた特性を与えるパラメータであって前記送信駆動波の生成および前記受信信号の補正に必要なパラメータを記憶する記憶部と、前記記憶部が記憶する前記パラメータを参照し、前記送信エコーの所定の周波数帯域における周波数スペクトルが前記超音波探触子の種類によらず同じとなるように、前記超音波探触子の種類に応じた前記送信駆動波を生成する送信駆動波生成部と、前記記憶部が記憶する前記パラメータを参照し、前記受信信号に対して前記超音波探触子の種類に応じた補正を行う受信信号補正部と、前記受信信号補正部が補正した前記受信信号を用いて画像データを生成する画像処理部と、を備えたことを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記送信駆動波生成部は、前記超音波探触子によって前記送信エコーに変換される際の減衰量が大きい周波数ほど値が大きい周波数スペクトルを有する前記送信駆動波を生成し、前記受信信号補正部は、前記受信信号に対して前記減衰量をもとに周波数ごとの補正を行うことを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記記憶部は、前記超音波探触子によって前記送信駆動波または前記受信エコーが変換された場合の周波数スペクトルにおける周波数ごとの減衰量を前記超音波探触子の種類と対応づけて前記パラメータとして記憶し、前記送信駆動波生成部は、前記送信エコーとして出力すべき所定の周波数帯域の周波数スペクトルに対して、前記記憶部が記憶する前記減衰量を周波数ごとに加えることによって前記送信駆動波を生成し、前記受信信号補正部は、前記超音波探触子から受信した受信信号の周波数スペクトルに対して、前記記憶部が記憶する前記減衰量を周波数ごとに加えることによって前記受信信号の補正を行うことを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記超音波探触子は、互いに種類が異なる複数の超音波探触子から選択可能であることを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記超音波探触子が受信した超音波の周波数を解析することによって周波数スペクトルを算出する周波数解析部と、前記周波数解析部が算出した周波数スペクトルを近似することによって前記周波数スペクトルから少なくとも１つの特徴量を抽出する特徴量抽出部と、をさらに備え、前記画像処理部は、前記特徴量抽出部が抽出した特徴量に応じた特徴量画像データを生成する特徴量画像データ生成部を有することを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記特徴量抽出部は、前記周波数スペクトルの近似処理を行う前または行った後に、超音波の受信深度および周波数に応じて発生する減衰の寄与を削減する減衰補正処理を行うことを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記特徴量抽出部は、回帰分析によって近似対象の周波数スペクトルを多項式で近似することを特徴とする。
　本発明に係る超音波観測装置は、上記発明において、前記特徴量抽出部は、前記近似対象の周波数スペクトルを一次式で近似し、前記一次式の傾き、前記一次式の切片、および前記傾きと前記切片と前記周波数スペクトルの周波数域に含まれる特定の周波数とを用いて定まる強度、の少なくとも１つを特徴量として抽出することを特徴とする。
　本発明に係る超音波観測装置の作動方法は、電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、前記検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子の種類に応じた特性を与えるパラメータであって前記送信駆動波の生成および前記受信信号の補正に必要なパラメータを記憶する記憶部から前記パラメータを参照し、前記超音波探触子の種類に応じた前記送信駆動波を送信駆動波生成部によって生成する送信駆動波生成ステップと、前記記憶部が記憶する前記パラメータを参照し、前記受信信号に対して前記超音波探触子の種類に応じた補正を行う受信信号補正ステップと、前記受信信号補正ステップで補正した前記受信信号を用いて画像データを生成する画像処理ステップと、を有することを特徴とする。
　本発明に係る超音波観測装置の作動プログラムは、電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、前記検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子の種類に応じた特性を与えるパラメータであって前記送信駆動波の生成および前記受信信号の補正に必要なパラメータを記憶する記憶部から前記パラメータを参照し、前記超音波探触子の種類に応じた前記送信駆動波を送信駆動波生成部によって生成する送信駆動波生成ステップと、前記記憶部が記憶する前記パラメータを参照し、前記受信信号に対して前記超音波探触子の種類に応じた補正を行う受信信号補正ステップと、前記受信信号補正ステップで補正した前記受信信号を用いて画像データを生成する画像処理ステップと、をコンピュータに実行させることを特徴とする。
　本発明によれば、超音波探触子の種類の違いによる影響を排除した超音波画像の観察を実現することができる。
図１は、本発明の実施の形態１に係る超音波観測装置の構成を示すブロック図である。 図２は、信号変換部が送信駆動波を送信エコーに変換した場合の周波数スペクトルの変化を示す図である。 図３は、送信駆動波生成部が生成する送信駆動波と、超音波探触子の信号変換部が送信駆動波を送信エコーに変換した場合を示す図である。 図４は、受信信号補正部が行う補正処理の概要（第１例）を示す図である。 図５は、受信信号補正部が行う補正処理の概要（第２例）を示す図である。 図６は、受信信号補正部によって補正された受信信号に対してＳＴＣ補正部が行う増幅処理における受信深度と増幅率との関係を示す図である。 図７は、本発明の実施の形態１に係る超音波観測装置の処理の概要を示すフローチャートである。 図８は、本発明の実施の形態２に係る超音波観測装置の構成を示すブロック図である。 図９は、増幅補正部が行う増幅処理における受信深度と増幅率との関係を示す図である。 図１０は、周波数解析部によって算出される周波数スペクトルの例を示す図である。 図１１は、超音波探触子を用いた場合の受信信号補正部が補正する前後の周波数スペクトル（第１例）を示す図である。 図１２は、超音波探触子を用いた場合の受信信号補正部が補正する前後の周波数スペクトル（第２例）を示す図である。 図１３は、図１１および図１２に示す周波数スペクトル曲線の特徴量を特徴量空間でプロットした図である。 図１４は、減衰補正部が補正した特徴量に対応する直線を示す図である。 図１５は、本発明の実施の形態２に係る超音波観測装置の処理の概要を示すフローチャートである。 図１６は、周波数解析部が行う処理の概要を示すフローチャートである。 図１７は、１つの音線のデータ配列を模式的に示す図である。 図１８は、本発明の別な実施の形態に係る超音波観測装置の減衰補正部が行う減衰補正処理の概要を模式的に示す図である。
　以下、添付図面を参照して、本発明を実施するための形態（以下、「実施の形態」という）を説明する。
（実施の形態１）
　図１は、本発明の実施の形態１に係る超音波観測装置の構成を示すブロック図である。同図に示す超音波観測装置１は、超音波を用いて診断対象である検体を観測するための装置である。超音波観測装置１は、外部へ超音波の送信エコーを出力するとともに、外部で反射された超音波の受信エコーを受信して電気的な受信信号に変換する超音波探触子２と、超音波探触子２との間で電気信号の送受信を行う送受信部３と、受信信号に対応する画像データの生成を行う画像処理部４と、キーボード、マウス、タッチパネル等のインタフェースを用いて実現され、各種情報の入力を受け付ける入力部５と、液晶または有機ＥＬ等からなる表示パネルを用いて実現され、画像処理部４が生成した画像を含む各種情報を表示する表示部６と、超音波探触子２の種類に応じた特性を示すパラメータを含む超音波観測に必要な各種情報を記憶する記憶部７と、超音波観測装置１の動作制御を行う制御部８と、を備える。超音波観測装置１は、先端部に超音波探触子２が設けられるスコープと、スコープの基端が着脱可能に接続され、超音波探触子２以外の上記部位が設けられる処理装置（プロセッサ）とによって構成される。
　超音波探触子２は、送受信部３から受信した電気的なパルス信号である送信駆動波を音響パルス信号である送信エコーに変換するとともに、外部の検体で反射されて戻ってきた受信エコーを電気的な受信信号に変換する信号変換部２１を有する。送信駆動波から送信エコーへの変換および受信エコーから受信信号への変換を行う際、信号変換部２１は、フィルタの機能を果たす。
　図２は、信号変換部２１が送信駆動波を送信エコーに変換した場合の周波数スペクトルの変化を示す図である。図２（ａ）では、横軸ｆが周波数、縦軸Ｉが強度を表している。送信駆動波のスペクトル（周波数スペクトル曲線）Ｃ1は、所定の周波数帯域Ｆ＝｛ｆ｜ｆL≦ｆ≦ｆH｝でゼロでない一定値を取る矩形状をなす。この送信駆動波が信号変換部２１によって送信エコーに変換されると、スペクトルの形状も変化する。図２（ｂ）および図２（ｃ）は、互いに種類が異なる超音波探触子２の信号変換部２１がそれぞれ生成した送信エコーのスペクトルを示す図である。図２（ｂ）に示すスペクトルＣ2は、図２（ｃ）に示すスペクトルＣ3よりも減衰が大きい。
　スペクトルＣ2において、周波数帯域Ｆ内の任意の周波数ｆ1（ｆL≦ｆ1≦ｆH）におけるスペクトルＣ1とスペクトルＣ2との差に相当する減衰量ΔＩ1（ｆ１）は、超音波探触子２の種類と対応付けて、記憶部７が有する超音波探触子情報記憶部７１に記憶されている。同様に、スペクトルＣ3において、周波数帯域Ｆ内の任意の周波数ｆ1におけるスペクトルＣ1とスペクトルＣ3との差に相当する減衰量ΔＩ2（ｆ１）も、超音波探触子２の種類と対応付けて超音波探触子情報記憶部７１に記憶されている。なお、ここでは２つの超音波探触子２を例にとって説明したが、他の種類の超音波探触子２についても同様にして超音波探触子情報記憶部７１に記憶させておくことが可能である。
　以下では、説明の便宜上、図２（ａ）に示すスペクトルＣ1を有する送信駆動波に対し、図２（ｂ）に示すスペクトルＣ2を送信エコーとして出力する超音波探触子２を超音波探触子２Ａ、図２（ｃ）に示すスペクトルＣ3を送信エコーとして出力する超音波探触子２を超音波探触子２Ｂという。また、超音波探触子２Ａの信号変換部２１を信号変換部２１Ａ、超音波探触子２Ｂの信号変換部２１を信号変換部２１Ｂという。
　超音波探触子２は、超音波振動子をメカ的に走査させるものであってもよいし、複数の超音波振動子を電子的に走査させるものであってもよい。本実施の形態１では、超音波探触子２として、互いに異なる複数種類のいずれかの超音波探触子２を選択して使用することが可能である。
　送受信部３は、超音波探触子２と電気的に接続される。送受信部３は、予め設定された波形および送信タイミングに基づいて送信駆動波を生成し、この生成した送信駆動波を超音波探触子２へ送信する送信駆動波生成部３１と、超音波探触子２から受信した受信信号を超音波探触子２の種類に応じた特性の違いによる影響を排除する補正を行う受信信号補正部３２と、受信信号補正部３２が補正した受信信号に対して、受信深度が大きいほど高い増幅率で補正するＳＴＣ（Sensitivity　Time　Control）補正を行うＳＴＣ補正部３３と、を有する。
　送信駆動波生成部３１は、後述する記憶部７の超音波探触子情報記憶部７１が記憶する情報に基づいて、超音波探触子２の特性をふまえた波形を送信駆動波として生成する。図３は、送信駆動波生成部３１が生成する送信駆動波と、超音波探触子２の信号変換部２１が送信駆動波を送信エコーに変換した場合を示す図である。送信駆動波生成部３１は、超音波探触子２の種類毎に異なる送信駆動波を生成する。具体的には、送信駆動波生成部３１は、超音波探触子２の種類に応じた特性を超音波探触子情報記憶部７１から読み出して、図３（ｃ）に示すような矩形のスペクトルＣ6を有する送信駆動波を生成する。ここでいう超音波探触子２の特性とは、図２（ｂ）、（ｃ）にそれぞれ示した減衰量ΔＩ1（ｆ）、ΔＩ2（ｆ）（ただしｆ∈Ｆ）と超音波探触子２の種類を含む個々の超音波探触子２に固有の情報を意味する。なお、本実施の形態１および後続する実施の形態２において、曲線および直線は、離散的な点の集合からなる。
　送信駆動波生成部３１は、生成する送信エコーのスペクトルと、超音波探触子２の特性とに基づいて、超音波探触子２の種類に応じた送信駆動波を生成する。ここで、送信駆動波生成部３１は、送信エコーの所定の周波数帯域における周波数スペクトルが前記超音波探触子の種類毎に異なる送信駆動波を生成する。図３（ａ）に示す送信駆動波のスペクトルＣ4は、超音波探触子２Ａに対するものであり、図３（ｂ）に示す送信駆動波のスペクトルＣ5は、超音波探触子２Ｂに対するものである。スペクトルＣ4は、送信エコーのスペクトルＣ6に対して、超音波探触子２Ａが変換する際の減衰量ΔＩ1（ｆ）（ｆ∈Ｆ）を加えた形状をなしている。同様に、スペクトルＣ5は、送信エコーのスペクトルＣ6に対して、超音波探触子２Ｂが変換する際の減衰量ΔＩ2（ｆ）（ｆ∈Ｆ）を加えた形状をなしている。
　受信信号補正部３２は、送受信部３が超音波探触子２から受信した受信信号に対して、超音波探触子の種類に応じた補正を加える。図４は、受信信号補正部３２が行う補正処理の概要を示す図であり、具体的には超音波探触子２Ａから送られてくる受信信号に対する補正処理の概要を示す図である。図４（ａ）に示すスペクトルＣ7は、超音波探触子２Ａの信号変換部２１Ａが受信エコーを変換して生成した受信信号の周波数スペクトルである。受信信号補正部３２は、スペクトルＣ7に対して超音波探触子２Ａが変換する際の減衰量ΔＩ1（ｆ）（ｆ∈Ｆ）を加えることによって受信信号を補正する。これにより、図４（ｂ）に示すスペクトルＣ8が得られる。このスペクトルＣ8は、超音波探触子２Ａの影響を排除したスペクトルである。
　図５は、受信信号補正部３２が行う補正処理の概要を示す図であり、具体的には超音波探触子２Ｂから送られてくる受信信号に対する補正処理の概要を示す図である。図５（ａ）に示すスペクトルＣ9は、超音波探触子２Ｂの信号変換部２１Ｂが受信エコーを変換して生成した受信信号の周波数スペクトルである。受信信号補正部３２は、スペクトルＣ9に対して超音波探触子２Ｂが変換する際の減衰量ΔＩ2（ｆ）（ｆ∈Ｆ）を加えることによって受信信号を補正する。これにより、図５（ｂ）に示すスペクトルＣ8が得られる。このスペクトルＣ8は、超音波探触子２Ｂの影響を排除したスペクトルであり、超音波探触子２Ａが変換した受信信号に対して補正を行った結果得られるスペクトルと同じである（図４（ｂ）を参照）。換言すると、スペクトルＣ8は、受信エコーのスペクトルに他ならない。
　図６は、受信信号補正部３２によって補正された受信信号に対してＳＴＣ補正部３３が行う増幅処理における受信深度と増幅率との関係を示す図である。図６に示す受信深度ｚは、超音波の受信開始時点からの経過時間に基づいて算出される量である。図６に示すように、増幅率β（ｄＢ）は、受信深度ｚが閾値ｚthより小さい場合、受信深度ｚの増加に伴ってβ０からβth（＞β0）へ線型に増加する。また、増幅率βは、受信深度ｚが閾値ｚth以上である場合、一定値βthをとる。閾値ｚthの値は、検体から受信する超音波信号がほとんど減衰してしまい、ノイズが支配的になるような値である。より一般に、増幅率βは、受信深度ｚが閾値ｚthより小さい場合、受信深度ｚの増加に伴って単調増加すればよい。
　送受信部３は、ＳＴＣ補正部３３によって増幅されたエコー信号に対してフィルタリング等の処理を施した後、Ａ／Ｄ変換することによって時間ドメインのデジタルＲＦ信号を生成して出力する。なお、超音波探触子２が複数の超音波振動子を電子的に走査させるものである場合、送受信部３は、複数の超音波振動子に対応したビーム合成用の多チャンネル回路を有する。
　画像処理部４は、受信信号からＢモード画像データを生成するＢモード画像データ生成部４１を有する。
　Ｂモード画像データ生成部４１は、デジタル信号に対してバンドパスフィルタ、対数変換、ゲイン処理、コントラスト処理等の公知の技術を用いた信号処理を行うとともに、表示部６における画像の表示レンジに応じて定まるデータステップ幅に応じたデータの間引き等を行うことによってＢモード画像データを生成する。
　記憶部７は、超音波探触子情報記憶部７１と、増幅率情報記憶部７２と、窓関数記憶部７３とを有する。
　超音波探触子情報記憶部７１は、上述したように、超音波探触子２の種類と、その超音波探触子２が送信駆動波を送信エコーに変換したり受信エコーを受信信号に変換したりする際にスペクトルに生じる減衰量とを関連付けたものをパラメータとして記憶する。
　増幅率情報記憶部７２は、ＳＴＣ補正部３３が増幅処理を行う際に参照する増幅率と受信深度との関係（例えば、図２に示す関係）を増幅率情報として記憶する。
　窓関数記憶部７３は、Hamming、Hanning、Blackmanなどの窓関数のうち少なくともいずれか１つの窓関数を記憶する。
　記憶部７は、超音波観測装置１の作動プログラムや所定のＯＳを起動するプログラム等が予め記憶されたＲＯＭ、および各処理の演算パラメータやデータ等を記憶するＲＡＭ等を用いて実現される。
　制御部８は、接続されている超音波探触子２の種類を判定する超音波探触子判定部８１を有する。超音波探触子判定部８１の判定結果は超音波探触子情報記憶部７１に記憶され、信号変換部２１、送信駆動波生成部３１、および受信信号補正部３２が処理を行う際に参照される。超音波探触子判定部８１が判定を行うために、超音波探触子２を備えたスコープには、後段の処理装置との接続部分に超音波探触子２の種類を処理装置に判別させるための接続ピンを設けておけばよい。これにより、処理装置側に設けられる超音波探触子判定部８１は、接続されたスコープの接続ピンの形状に応じて超音波探触子２の種類を判定することができる。
　制御部８は、演算および制御機能を有するＣＰＵを用いて実現される。制御部８は、記憶部７が記憶、格納する情報および超音波観測装置１の作動プログラムを含む各種プログラムを記憶部７から読み出すことにより、超音波観測装置１の作動方法に関連した各種演算処理を実行することによって超音波観測装置１を統括して制御する。
　なお、超音波観測装置１の作動プログラムは、ハードディスク、フラッシュメモリ、ＣＤ－ＲＯＭ、ＤＶＤ－ＲＯＭ、フレキシブルディスク等のコンピュータ読み取り可能な記録媒体に記録して広く流通させることも可能である。各種プログラムの記録媒体等への記録は、コンピュータまたは記録媒体を製品として出荷する際に行ってもよいし、通信ネットワークを介したダウンロードにより行ってもよい。
　図７は、以上の構成を有する超音波観測装置１の処理の概要を示すフローチャートである。なお、超音波観測装置１が備える超音波探触子２の種類は、超音波探触子判定部８１が予め判定しているものとする。
　図７において、送信駆動波生成部３１は送信駆動波を生成して超音波探触子２へ出力する（ステップＳ１）。この送信駆動波のスペクトルは、例えば超音波探触子２Ａの場合には、図３（ａ）に示すスペクトルＣ4である。また、超音波探触子２Ｂの場合、送信駆動波のスペクトルは、図３（ｂ）に示すスペクトルＣ5である。
　続いて、送信駆動波を受信した超音波探触子２では、信号変換部２１が送信駆動波を送信エコーに変換し、変換した送信エコーを生成して出力する（ステップＳ２）。この送信エコーのスペクトルは、例えば図３（ｃ）のスペクトルＣ6である。図３を参照して説明したように、送信エコーは超音波探触子２の種類によらず同じスペクトルを有している。
　この後、超音波探触子２は、送信エコーが生体で反射されて戻ってきた受信エコーを受信する（ステップＳ３）。
　受信エコーを受信した信号変換部２１は、受信エコーを受信信号に変換して送受信部３へ出力する（ステップＳ４）。この受信信号のスペクトルは、例えば超音波探触子２Ａの場合には、図４（ａ）に示すスペクトルＣ7である。また、超音波探触子２Ｂの場合、受信信号のスペクトルは、図５（ａ）に示すスペクトルＣ9である。
　超音波探触子２から受信信号を受信した受信信号補正部３２は、周波数スペクトルの補正を行う（ステップＳ５）。この受信信号のスペクトルは、超音波探触子２の種類によらず同じスペクトルとなる（図４（ｂ）および図５（ｂ）のスペクトルＣ8を参照）。
　続いて、ＳＴＣ補正部３３は、受信信号補正部３２によってスペクトルが補正された受信信号に対してＳＴＣ補正を行う（ステップＳ６）。ここで、ＳＴＣ補正部３３は、例えば図６に示す増幅率と受信深度との関係に基づいてＳＴＣ補正を行う。
　この後、Ｂモード画像データ生成部４１は、ＳＴＣ補正部３３が増幅したエコー信号を用いてＢモード画像データを生成する（ステップＳ７）。表示部６は、Ｂモード画像データ生成部４１が生成したＢモード画像に対応するＢモード画像を表示する（ステップＳ８）。
　ステップＳ８の後、超音波観測装置１は、一連の処理を終了する。なお、超音波観測装置１が、ステップＳ１～Ｓ８の処理を周期的に繰り返すようにしてもよい。
　以上説明した本発明の実施の形態１によれば、送信駆動波と受信信号の両方に対して超音波探触子の種類に応じた影響を排除するための処理を行っているため、超音波探触子の種類の違いによる影響を排除した超音波画像の観察を実現することができる。
　また、本実施の形態１によれば、超音波探触子によって送信エコーに変換される際の減衰量が大きい周波数ほど値が大きい周波数スペクトルを有する前記送信駆動波を生成し、受信信号に対してその減衰量をもとに周波数ごとの補正を行うため、超音波探触子の種類に応じた特性の違いの影響を確実に排除することができる。
（実施の形態２）
　図８は、本発明の実施の形態２に係る超音波観測装置の構成を示すブロック図である。同図に示す超音波観測装置１１は、超音波探触子２と、送受信部３と、入力部５と、表示部６と、制御部８と、電気的な受信信号に対して所定の演算を施す演算部１２と、電気的なエコー信号に対応する画像データの生成を行う画像処理部１３と、記憶部１４と、を備える。なお、上述した超音波観測装置１が有する構成要素と同様の機能を有する構成要素には、超音波観測装置１と同じ符号を付してある。
　演算部１２は、送受信部３が出力したデジタルＲＦ信号に対して受信深度によらず増幅率を一定とする増幅補正を行う増幅補正部１２１と、増幅補正を行ったデジタルＲＦ信号に高速フーリエ変換（ＦＦＴ）を施して周波数解析を行うことにより周波数スペクトルを算出する周波数解析部１２２と、周波数解析部１２２が算出した各箇所の周波数スペクトルに対し、回帰分析に基づく近似処理および超音波が伝播する際に該超音波の受信深度および周波数に応じて発生する減衰の寄与を削減する減衰補正処理を行うことにより、検体の特徴量を抽出する特徴量抽出部１２３と、を有する。
　図９は、増幅補正部１２１が行う増幅処理における受信深度と増幅率との関係を示す図である。図９に示すように、増幅補正部１２１が行う増幅処理における増幅率β（ｄＢ）は、受信深度ｚがゼロのとき最大値βth－β0をとり、受信深度ｚがゼロから閾値ｚthに達するまで線型に減少し、受信深度ｚが閾値ｚth以上のときゼロである。このように定められる増幅率によって増幅補正部１２１がデジタルＲＦ信号を増幅補正することにより、ＳＴＣ補正部３３におけるＳＴＣ補正の影響を相殺し、一定の増幅率βthの信号を出力することができる。なお、増幅補正部１２１が行う受信深度ｚと増幅率βの関係は、ＳＴＣ補正部３３における受信深度と増幅率の関係に応じて異なることは勿論である。
　このような増幅補正を行う理由を説明する。ＳＴＣ補正は、アナログ信号波形の振幅を全周波数帯域にわたって均一に増幅させる補正である。このため、超音波の振幅を利用するＢモード画像を生成する際には、ＳＴＣ補正を行うことによって十分な効果を得ることができる一方で、超音波の周波数スペクトルを算出するような場合には、超音波の伝播に伴う減衰の影響を正確に排除できるわけではない。この問題を解決するには、Ｂモード画像を生成する際にＳＴＣ補正を施した受信信号を出力する一方、周波数スペクトルに基づいた画像を生成する際に、Ｂモード画像を生成するための送信とは異なる新たな送信を行い、ＳＴＣ補正を施していない受信信号を出力することが考えられる。ところがこの場合には、受信信号に基づいて生成される画像データのフレームレートが低下してしまうという問題がある。そこで、本実施の形態２では、生成される画像データのフレームレートを維持しつつ、Ｂモード画像用にＳＴＣ補正を施した信号に対して一度ＳＴＣ補正の影響を排除するために、増幅補正部１２１によって増幅率の補正を行っている。
　周波数解析部１２２は、各音線（ラインデータ）に対し、所定のデータ量からなるＦＦＴデータ群を高速フーリエ変換することによって音線上の複数の箇所（データ位置）における周波数スペクトルを算出する。周波数解析部１２２による算出結果は複素数で得られ、記憶部１４に記憶される。
　一般に、周波数スペクトルは、検体の組織性状によって異なる傾向を示す。これは、周波数スペクトルが、超音波を散乱する散乱体としての検体の大きさ、密度、音響インピーダンス等と相関を有しているためである。なお、本実施の形態２において、「組織性状」とは、例えば癌、内分泌腫瘍、粘液性腫瘍、正常組織、脈管などのいずれかである。
　図１０は、周波数解析部１２２によって算出される周波数スペクトルの例を示す図である。具体的には、ＦＦＴデータ群を高速フーリエ変換することによって得られた周波数スペクトルを、周波数ｆおよび受信深度をｚの関数として、強度Ｉ（ｆ，ｚ）と位相φ（ｆ，ｚ）によって表現したときの強度Ｉ（ｆ，ｚ）のスペクトルを示している。ここでいう「強度」とは、電圧、電力、音圧、音響エネルギー等のパラメータのいずれかを指す。図１０では、横軸ｆが周波数、縦軸Ｉが強度であり、受信深度ｚは一定である。図１０に示すスペクトルＣ11において、周波数スペクトルの下限周波数ｆLおよび上限周波数ｆHは、超音波探触子２の周波数帯域、送受信部３が送信するパルス信号の周波数帯域などをもとに決定されるパラメータであり、例えばｆL＝３ＭＨｚ、ｆH＝１０ＭＨｚである。
　特徴量抽出部１２３は、周波数解析部１２２が算出した周波数スペクトルの近似式を回帰分析によって算出する近似部１２４と、近似部１２４が算出した近似式に対し、超音波の受信深度および周波数に依存する超音波の減衰の寄与を削減する減衰補正処理を施すことによって周波数スペクトルの特徴量を抽出する減衰補正部１２５と、を有する。
　近似部１２４は、回帰分析によって周波数スペクトルを一次式（回帰直線）で近似することにより、この近似した一次式を特徴付ける減衰補正前の特徴量（以下、補正前特徴量という）を抽出する。具体的には、近似部１２４は、一次式の傾きａ0および切片ｂ0を補正前特徴量として抽出する。図１０に示す直線Ｌ10は、近似部１２４が近似した一次式に相当する直線である。なお、近似部１２４は、傾きａ0および切片ｂ0以外の補正前特徴量として、周波数帯域（ｆL＜ｆ＜ｆH）の中心周波数ｆM＝（ｆL＋ｆH）／２における強度（Mid-band　fitともいう）ｃ0＝ａ0ｆM＋ｂ0を算出してもよい。
　３つの特徴量のうち、傾きａ0は、超音波の散乱体の大きさと相関を有し、一般に散乱体が大きいほど傾きが小さな値を有すると考えられる。また、切片ｂ0は、散乱体の大きさ、音響インピーダンスの差、散乱体の密度（濃度）等と相関を有している。具体的には、切片ｂ0は、散乱体が大きいほど大きな値を有し、音響インピーダンスが大きいほど大きな値を有し、散乱体の密度（濃度）が大きいほど大きな値を有すると考えられる。中心周波数ｆMにおける強度（以下、単に「強度」という）ｃ0は、傾きａ0と切片ｂ0から導出される間接的なパラメータであり、有効な周波数帯域内の中心におけるスペクトル強度を与える。このため、強度ｃ0は、散乱体の大きさ、音響インピーダンスの差、散乱体の密度に加えて、Ｂモード画像の輝度とある程度の相関を有していると考えられる。なお、特徴量抽出部１２３が算出する近似多項式は一次式に限定されるわけではなく、二次以上の近似多項式を用いることも可能である。
　ここで、図１１～図１３を参照して、補正前特徴量の値が超音波探触子２の種類によらず同じであることを説明する。なお、図１１～図１３では、２つの異なる種類の超音波探触子２Ａ、２Ｂを用いた場合を説明するため、実施の形態１において参照した周波数スペクトルを使用する。また、ここでは、補正前特徴量として切片ｂ0および強度ｃ0を用いて説明する。
　図１１は、超音波探触子２Ａを用いた場合の受信信号補正部３２が補正する前後の周波数スペクトルを示す図である。具体的には、スペクトルＣ7が受信信号補正前であり、スペクトルＣ8が受信信号補正後のものである。図１１において、スペクトルＣ7の補正前特徴量は、ｂ0＝ｂ1（直線Ｌ1の切片）、ｃ0＝ｃ1（直線Ｌ1の強度）である。また、スペクトルＣ8の補正前特徴量は、ｂ0＝ｂ3（直線Ｌ3の切片）、ｃ0＝ｃ3（直線Ｌ3の強度）である。
　図１２は、超音波探触子２Ｂを用いた場合の受信信号補正部３２が補正する前後の周波数スペクトルを示す図である。具体的には、スペクトルＣ9が受信信号補正前であり、スペクトルＣ8が受信信号補正後のものである。図１２において、スペクトルＣ9の補正前特徴量は、ｂ0＝ｂ2（直線Ｌ2の切片）、ｃ0＝ｃ2（直線Ｌ2の強度）である。ここで、一般にｂ2≠ｂ1，ｂ3であるとともに、ｃ2≠ｃ1，ｃ3である。
　図１３は、図１１および図１２に示すスペクトルＣ7～Ｃ9の特徴量の組を特徴量空間（ｂ0，ｃ0）でそれぞれプロットした図である。図１３に示すように、受信信号補正前のスペクトルＣ7の補正前特徴量の組である点Ａ（ｂ1，ｃ1）と受信信号補正前のスペクトルＣ9の補正前特徴量の組である点Ｂ（ｂ2，ｃ2）は、一般に特徴量空間上の異なる点である。また、受信信号補正部３２が補正した後のスペクトルは、補正前のスペクトルによらず同じスペクトルＣ8となるため、補正後の特徴量の組（ｂ0，ｃ0）は、補正前のスペクトルによらず点Ｃ（ｂ3，ｃ3）となる。
　このように、受信信号の周波数スペクトルの補正前特徴量が一致するのは、送信駆動波生成部３１が送信駆動波を生成する際に超音波探触子２の機種に依存した周波数スペクトルの影響を排除するとともに、受信信号補正部３２が受信信号に対して超音波探触子２の機種に依存した周波数スペクトルの影響を排除したためである。この点につき、従来技術のように、送信駆動波または受信信号のいずれか一方のみの補正を行うだけでは、補正を行わない他方において超音波探触子２の種類の違いの影響を排除することができないため、受信信号の周波数スペクトルが超音波探触子２の種類に応じて異なるものとなってしまう。したがって、本実施の形態２によれば、超音波探触子２の種類の違いによる影響を確実に排除することが可能となり、より正確な特徴量の抽出を実現することができる。
　次に、減衰補正部１２５が行う補正について説明する。一般に、超音波の減衰量Ａ（ｆ，ｚ）は、
　　Ａ（ｆ，ｚ）＝２αｚｆ　　・・・（１）
と表される。ここで、αは減衰率であり、ｚは超音波の受信深度であり、ｆは周波数である。式（１）からも明らかなように、減衰量Ａ（ｆ，ｚ）は、周波数ｆに比例している。減衰率αの具体的な値は、観察対象が生体である場合、０．０～１．０（ｄＢ／ｃｍ／ＭＨｚ）、より好ましくは０．３～０．７（ｄＢ／ｃｍ／ＭＨｚ）であり、生体の部位に応じて定まる。例えば、観察対象が膵臓である場合には、α＝０．６（ｄＢ／ｃｍ／ＭＨｚ）と定めることがある。なお、本実施の形態において、減衰率αの値を入力部５からの入力によって設定または変更可能な構成としてもよい。
　減衰補正部１２５は、近似部１２４が抽出した補正前特徴量（傾きａ0，切片ｂ0，強度ｃ0）を、以下のように減衰補正することによって特徴量を抽出する。
　　ａ＝ａ0＋２αｚ　　・・・（２）
　　ｂ＝ｂ0　　・・・（３）
　　ｃ＝ｃ0＋２αｚｆM（＝ａｆM＋ｂ）　　・・・（４）
式（２）、（４）からも明らかなように、減衰補正部１２５は、超音波の受信深度ｚが大きいほど、補正量が大きい補正を行う。また、式（３）によれば、切片に関する補正は恒等変換である。これは、切片が周波数０（Ｈｚ）に対応する周波数成分であって減衰の影響を受けないためである。
　図１４は、減衰補正部１２５が補正した特徴量に対応する直線を示す図である。直線Ｌ11の式は、
　　Ｉ＝ａｆ＋ｂ＝（ａ0＋２αｚ）ｆ＋ｂ0　　・・・（５）
で表される。この式（５）からも明らかなように、直線Ｌ11は、直線Ｌ10と比較して、傾きが大きく、かつ切片が同じである。
　画像処理部１３は、Ｂモード画像データ生成部４１と、特徴量抽出部１２３が抽出した特徴量に対応する情報を複数の表示方法のいずれかにしたがって表示する特徴量画像データを生成する特徴量画像データ生成部１３１と、を有する。
　特徴量画像データ生成部１３１は、特徴量画像データにおいて各画素に割り当てられる情報は、周波数解析部１２２が周波数スペクトルを算出する際のＦＦＴデータ群のデータ量に応じて定められる。具体的には、例えば１つのＦＦＴデータ群のデータ量に対応する画素領域には、そのＦＦＴデータ群から算出される周波数スペクトルの特徴量に対応する情報が割り当てられる。なお、本実施の形態２において、特徴量画像データを生成する際に使用する特徴量の数は任意に設定することが可能である。
　記憶部１４は、超音波探触子情報記憶部７１、増幅率情報記憶部７２、窓関数記憶部７３に加えて、補正情報記憶部１４１を有する。
　増幅率情報記憶部７２は、ＳＴＣ補正部３３が増幅処理を行う際に参照する増幅率と受信深度との関係（例えば、図６に示す関係）および増幅補正部１２１が増幅補正処理を行う際にそれぞれ参照する増幅率と受信深度との関係（例えば、図９に示す関係）を増幅率情報として記憶する。
　補正情報記憶部１４１は、式（１）を含む減衰補正に関連した情報を記憶する。
　記憶部１４は、超音波観測装置１１の作動プログラムや所定のＯＳを起動するプログラム等が予め記憶されたＲＯＭ、および各処理の演算パラメータやデータ等を記憶するＲＡＭ等を用いて実現される。
　図１５は、以上の構成を有する超音波観測装置１１の処理の概要を示すフローチャートである。図１５において、ステップＳ１１～Ｓ１７の処理は、図７に示すフローチャートのステップＳ１～Ｓ７に順次対応している。以下、ステップＳ１８以降の処理を説明する。
　ステップＳ１８において、増幅補正部１２１は、送受信部３から出力された信号に対して受信深度によらず増幅率が一定となる増幅補正を行う（ステップＳ１８）。ここで、増幅補正部１２１は、例えば図９に示す増幅率と受信深度との関係に基づいて増幅補正を行う。
　その後、周波数解析部１２２は、ＦＦＴ演算による周波数解析を行うことによって周波数スペクトルを算出する（ステップＳ１９）。
　ここで、周波数解析部１２２が行う処理（ステップＳ１９）について、図１６に示すフローチャートを参照して詳細に説明する。まず、周波数解析部１２２は、解析対象の音線を識別するカウンタｋをｋ0とする（ステップＳ３１）。
　続いて、周波数解析部１２２は、ＦＦＴ演算用に取得する一連のデータ群（ＦＦＴデータ群）を代表するデータ位置（受信深度に相当）Ｚ(k)の初期値Ｚ(k) 0を設定する（ステップＳ３２）。図１７は、１つの音線のデータ配列を模式的に示す図である。同図に示す音線ＳＲkにおいて、白または黒の長方形は、１つのデータを意味している。音線ＳＲkは、送受信部３が行うＡ／Ｄ変換におけるサンプリング周波数（例えば５０ＭＨｚ）に対応した時間間隔で離散化されている。図１７では、音線ＳＲkの１番目のデータ位置を初期値Ｚ(k) 0として設定した場合を示しているが、初期値の位置は任意に設定することができる。
　その後、周波数解析部１２２は、データ位置Ｚ(k)のＦＦＴデータ群を取得し（ステップＳ３３）、取得したＦＦＴデータ群に対し、窓関数記憶部７３が記憶する窓関数を作用させる（ステップＳ３４）。このようにＦＦＴデータ群に対して窓関数を作用させることにより、ＦＦＴデータ群が境界で不連続になることを回避し、アーチファクトが発生するのを防止することができる。
　続いて、周波数解析部１２２は、データ位置Ｚ(k)のＦＦＴデータ群が正常なデータ群であるか否かを判定する（ステップＳ３５）。ここで、ＦＦＴデータ群は、２のべき乗のデータ数を有している必要がある。以下、ＦＦＴデータ群のデータ数を２n（ｎは正の整数）とする。ＦＦＴデータ群が正常であるとは、データ位置Ｚ(k)がＦＦＴデータ群で前から２n-1番目の位置であることを意味する。換言すると、ＦＦＴデータ群が正常であるとは、データ位置Ｚ(k)の前方に２n-1－１（＝Ｎとする）個のデータがあり、データ位置Ｚ(k)の後方に２n-1（＝Ｍとする）個のデータがあることを意味する。図１７に示す場合、ＦＦＴデータ群Ｆ2、Ｆ3はともに正常である。なお、図１７ではｎ＝４（Ｎ＝７，Ｍ＝８）の場合を例示している。
　ステップＳ３５における判定の結果、データ位置Ｚ(k)のＦＦＴデータ群が正常である場合（ステップＳ３５：Ｙｅｓ）、周波数解析部１２２は、後述するステップＳ３７へ移行する。
　ステップＳ３５における判定の結果、データ位置Ｚ(k)のＦＦＴデータ群が正常でない場合（ステップＳ３５：Ｎｏ）、周波数解析部１２２は、不足分だけゼロデータを挿入することによって正常なＦＦＴデータ群を生成する（ステップＳ３６）。ステップＳ３５において正常でないと判定されたＦＦＴデータ群は、ゼロデータを追加する前に窓関数が作用されている。このため、ＦＦＴデータ群にゼロデータを挿入してもデータの不連続は生じない。ステップＳ３６の後、周波数解析部１２２は、後述するステップＳ３７へ移行する。
　ステップＳ３７において、周波数解析部１２２は、ＦＦＴデータ群を用いてＦＦＴ演算を行うことにより、複素数からなる周波数スペクトルを得る（ステップＳ３７）。この結果、例えば図１０に示すようなスペクトルＣ11が得られる。
　続いて、周波数解析部１２２は、データ位置Ｚ(k)をステップ幅Ｄで変化させる（ステップＳ３８）。ステップ幅Ｄは、記憶部１４が予め記憶しているものとする。図１７では、Ｄ＝１５の場合を例示している。ステップ幅Ｄは、Ｂモード画像データ生成部４１がＢモード画像データを生成する際に利用するデータステップ幅と一致させることが望ましいが、周波数解析部１２２における演算量を削減したい場合には、そのデータステップ幅より大きい値を設定してもよい。
　その後、周波数解析部１２２は、データ位置Ｚ(k)が音線ＳＲkにおける最大値Ｚ(k) maxより大きいか否かを判定する（ステップＳ３９）。データ位置Ｚ(k)が最大値Ｚ(k) maxより大きい場合（ステップＳ３９：Ｙｅｓ）、周波数解析部１２２はカウンタｋを１増加させる（ステップＳ４０）。一方、データ位置Ｚ(k)が最大値Ｚ(k) max以下である場合（ステップＳ３９：Ｎｏ）、周波数解析部１２２はステップＳ３３へ戻る。このようにして、周波数解析部１２２は、音線ＳＲkに対して、［｛（Ｚ(k) max－Ｚ(k) 0）／Ｄ｝＋１］個のＦＦＴデータ群に対するＦＦＴ演算を行う。ここで、［Ｘ］は、Ｘを超えない最大の整数を表す。
　ステップＳ４０の後、周波数解析部１２２は、カウンタｋが最大値ｋmaxより大きいか否かを判定する（ステップＳ４１）。カウンタｋがｋmaxより大きい場合（ステップＳ４１：Ｙｅｓ）、周波数解析部１２２は一連のＦＦＴ処理を終了する。一方、カウンタｋがｋmax以下である場合（ステップＳ４１：Ｎｏ）、周波数解析部１２２はステップＳ３２に戻る。
　このようにして、周波数解析部１２２は、（ｋmax－ｋ0＋１）本の音線の各々について複数回のＦＦＴ演算を行う。
　なお、ここでは、周波数解析部１２２が超音波信号を受信したすべての領域に対して周波数解析処理を行うことを前提としているが、あらかじめ入力部５によって特定の関心領域の設定入力を受け付けて、その関心領域内においてのみ周波数解析処理を行うようにしてもよい。
　以上説明したステップＳ１９の周波数解析処理に続いて、近似部１２４は、近似処理として周波数解析部１２２が算出した周波数スペクトルを回帰分析することにより、補正前特徴量を抽出する（ステップＳ２０）。具体的には、近似部１２４は、周波数スペクトル周波数帯域ｆL＜ｆ＜ｆHの周波数スペクトルの強度Ｉ（ｆ，ｚ）を近似する一次式を回帰分析によって算出することにより、この一次式を特徴づける傾きａ0、切片ｂ0（、強度ｃ0）を補正前特徴量として抽出する。図１０に示す直線Ｌ10は、このステップＳ２０において、スペクトルＣ11に対して補正前特徴量抽出処理を行うことによって得られる回帰直線の一例である。
　この後、減衰補正部１２５は、近似部１２４が抽出した補正前特徴量に対して減衰補正処理を行う（ステップＳ２１）。例えば、データのサンプリング周波数が５０ＭＨｚである場合、データのサンプリングの時間間隔は２０（ｎｓｅｃ）である。ここで、音速を１５３０（ｍ／ｓｅｃ）とすると、データのサンプリング距離間隔は、１５３０（ｍ／ｓｅｃ）×２０（ｎｓｅｃ）／２＝０．０１５３（ｍｍ）となる。処理対象のＦＦＴデータ群のデータ位置までの音線の１番目のデータからのデータステップ数がｎであるとすると、そのデータ位置Ｚは、データステップ数ｎとデータステップ幅Ｄとを用いて０．０１５３ｎＤ（ｍｍ）となる。減衰補正部１２５は、このようにして求まるデータ位置Ｚの値を上述した式（２）～（４）の受信深度ｚへ代入することにより、周波数スペクトルの特徴量である傾きａ、切片ｂ（、強度ｃ）を算出する。このように算出される特徴量に相当する直線の一例として、図１４に示す直線Ｌ11を挙げることができる。
　以上説明したステップＳ２０およびＳ２１は、特徴量抽出部１２３が周波数スペクトルを近似することによってその周波数スペクトルから少なくとも１つの特徴量を抽出する特徴量抽出ステップを構成する。
　この後、特徴量画像データ生成部１３１は、特徴量抽出ステップ（ステップＳ２０，Ｓ２１）で抽出した特徴量を用いて特徴量画像データを生成する（ステップＳ２２）。
　続いて、表示部６は、Ｂモード画像データ生成部４１が生成したＢモード画像および／または特徴量画像データ生成部１３１が生成した特徴量画像を表示する（ステップＳ２３）。この際、表示部６は、Ｂモード画像および特徴量画像のいずれか一方を表示してもよいし、Ｂモード画像と特徴量画像を並べて表示してもよいし、Ｂモード画像と特徴量画像を重畳して表示してもよい。Ｂモード画像と特徴量画像を重畳して表示する際には、Ｂモード画像と特徴量画像との混合比率を入力部５からの入力によって変化させることができるようにしてもよい。
　このようにして、特徴量画像をＢモード画像とともに表示部６で表示することにより、医師等のユーザは、Ｂモード画像による情報とあわせて検体の組織性状を判別することが可能となり、より高精度の診断を行うことが可能となる。
　ステップＳ２３の後、超音波観測装置１１は、一連の処理を終了する。なお、超音波観測装置１１が、ステップＳ１１～Ｓ２３の処理を周期的に繰り返すようにしてもよい。
　以上説明した本発明の実施の形態２によれば、送信駆動波と受信信号の両方に対して超音波探触子の種類に応じた影響を排除するための処理を行っているため、実施の形態１と同様に、超音波探触子の種類の違いによる影響を排除した超音波画像の観察を実現することができる。
　また、本実施の形態２によれば、超音波探触子によって送信エコーに変換される際の減衰量が大きい周波数ほど値が大きい周波数スペクトルを有する送信駆動波を生成し、受信信号に対してその減衰量をもとに周波数ごとの補正を行うため、超音波探触子の種類に応じた特性の違いの影響を確実に排除することができる。
　また、本実施の形態２においては、周波数解析を行って特徴量の抽出を行っているが、この特徴量には超音波探触子の種類に応じた特性の違いが排除されている。したがって、本実施の形態２によれば、超音波画像の定量的な評価を行う場合であっても、超音波探触子の種類に影響されることがない。
（その他の実施の形態）
　ここまで、本発明を実施するための形態を説明してきたが、本発明は、上述した一実施の形態によってのみ限定されるべきものではない。例えば、本発明において、特徴量抽出部１２３が、周波数スペクトルの減衰補正を行ってから、補正後の周波数スペクトルの近似式を算出するようにしてもよい。図１８は、減衰補正部１２５が行う減衰補正処理の概要を模式的に示す図である。図１８に示すように、減衰補正部１２５は、スペクトルＣ11に対し、帯域内のすべての周波数ｆ（ｆL＜ｆ＜ｆH）における強度Ｉ（ｆ，ｚ）に式（１）の減衰量Ａ（ｆ，ｚ）をそれぞれ加える補正（Ｉ（ｆ，ｚ）→Ｉ（ｆ，ｚ）＋Ａ（ｆ，ｚ））を行う。これにより、超音波の伝播に伴う減衰の寄与を削減した新たなスペクトルＣ12が得られる。近似部１２４は、スペクトルＣ12に対して回帰分析を行うことにより、特徴量を抽出する。この場合に抽出される特徴量は、図１８に示す直線Ｌ11の傾きａ、切片ｂ（、強度ｃ）である。この直線は、図１４に示す直線Ｌ11と同じである。
　また、本発明において、制御部８が、増幅補正部１２１による増幅補正処理と減衰補正部１２５における減衰補正処理とを一括して行わせるようにしてもよい。この処理は、図１５のステップＳ１８における増幅補正処理を行わず、図１５のステップＳ２１における減衰補正処理の減衰量の定義を次式（６）のように変更して行うことと等価である。
　　Ａ'＝２αｚｆ＋γ（ｚ）　・・・（６）
ここで、右辺のγ（ｚ）は、受信深度ｚにおける増幅率βとβ0との差であり、
　　γ（ｚ）＝－｛（βth－β0）／ｚth｝ｚ＋βth－β0　（ｚ≦ｚth）　・・・（７）
　　γ（ｚ）＝０　（ｚ＞ｚth）　・・・（８）
と表される。
　このように、本発明は、特許請求の範囲に記載した技術的思想を逸脱しない範囲内において、様々な実施の形態を含みうるものである。
　１、１１　超音波観測装置
　２　超音波探触子
　３　送受信部
　４、１３　画像処理部
　５　入力部
　６　表示部
　７　記憶部
　８　制御部
　１２　演算部
　１４　記憶部
　２１　信号変換部
　３１　送信駆動波生成部
　３２　受信信号補正部
　３３　ＳＴＣ補正部
　４１　Ｂモード画像データ生成部
　７１　超音波探触子情報記憶部
　７２　増幅率情報記憶部
　７３　窓関数記憶部
　８１　超音波探触子判定部
　１２１　増幅補正部
　１２２　周波数解析部
　１２３　特徴量抽出部
　１２４　近似部
　１２５　減衰補正部
　１３１　特徴量画像データ生成部
　１４１　補正情報記憶部

Claims (10)

  1. 　電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、前記検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子と、
    　前記超音波探触子の種類に応じた特性を与えるパラメータであって前記送信駆動波の生成および前記受信信号の補正に必要なパラメータを記憶する記憶部と、
    　前記記憶部が記憶する前記パラメータを参照し、前記送信エコーの所定の周波数帯域における周波数スペクトルが前記超音波探触子の種類によらず同じとなるように、前記超音波探触子の種類に応じた前記送信駆動波を生成する送信駆動波生成部と、
    　前記記憶部が記憶する前記パラメータを参照し、前記受信信号に対して前記超音波探触子の種類に応じた補正を行う受信信号補正部と、
    　前記受信信号補正部が補正した前記受信信号を用いて画像データを生成する画像処理部と、
    　を備えたことを特徴とする超音波観測装置。
  2. 　前記送信駆動波生成部は、
    　前記超音波探触子によって前記送信エコーに変換される際の減衰量が大きい周波数ほど値が大きい周波数スペクトルを有する前記送信駆動波を生成し、
    　前記受信信号補正部は、
    　前記受信信号に対して前記減衰量をもとに周波数ごとの補正を行うことを特徴とする請求項１に記載の超音波観測装置。
  3. 　前記記憶部は、
    　前記超音波探触子によって前記送信駆動波または前記受信エコーが変換された場合の周波数スペクトルにおける周波数ごとの減衰量を前記超音波探触子の種類と対応づけて前記パラメータとして記憶し、
    　前記送信駆動波生成部は、
    　前記送信エコーとして出力すべき所定の周波数帯域の周波数スペクトルに対して、前記記憶部が記憶する前記減衰量を周波数ごとに加えることによって前記送信駆動波を生成し、
    　前記受信信号補正部は、
    　前記超音波探触子から受信した受信信号の周波数スペクトルに対して、前記記憶部が記憶する前記減衰量を周波数ごとに加えることによって前記受信信号の補正を行うことを特徴とする請求項１または２に記載の超音波観測装置。
  4. 　前記超音波探触子は、互いに種類が異なる複数の超音波探触子から選択可能であることを特徴とする請求項１～３のいずれか一項に記載の超音波観測装置。
  5. 　前記超音波探触子が受信した超音波の周波数を解析することによって周波数スペクトルを算出する周波数解析部と、
    　前記周波数解析部が算出した周波数スペクトルを近似することによって前記周波数スペクトルから少なくとも１つの特徴量を抽出する特徴量抽出部と、
    　をさらに備え、
    　前記画像処理部は、
    　前記特徴量抽出部が抽出した特徴量に応じた特徴量画像データを生成する特徴量画像データ生成部を有することを特徴とする請求項１～４のいずれか一項に記載の超音波観測装置。
  6. 　前記特徴量抽出部は、
    　前記周波数スペクトルの近似処理を行う前または行った後に、超音波の受信深度および周波数に応じて発生する減衰の寄与を削減する減衰補正処理を行うことを特徴とする請求項５に記載の超音波観測装置。
  7. 　前記特徴量抽出部は、
    　回帰分析によって近似対象の周波数スペクトルを多項式で近似することを特徴とする請求項５または６に記載の超音波観測装置。
  8. 　前記特徴量抽出部は、
    　前記近似対象の周波数スペクトルを一次式で近似し、前記一次式の傾き、前記一次式の切片、および前記傾きと前記切片と前記周波数スペクトルの周波数域に含まれる特定の周波数とを用いて定まる強度、の少なくとも１つを特徴量として抽出することを特徴とする請求項７に記載の超音波観測装置。
  9. 　電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、前記検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子の種類に応じた特性を与えるパラメータであって前記送信駆動波の生成および前記受信信号の補正に必要なパラメータを記憶する記憶部から前記パラメータを参照し、前記超音波探触子の種類に応じた前記送信駆動波を送信駆動波生成部によって生成する送信駆動波生成ステップと、
    　前記記憶部が記憶する前記パラメータを参照し、前記受信信号に対して前記超音波探触子の種類に応じた補正を行う受信信号補正ステップと、
    　前記受信信号補正ステップで補正した前記受信信号を用いて画像データを生成する画像処理ステップと、
    　を有することを特徴とする超音波観測装置の作動方法。
  10. 　電気的な送信駆動波を超音波の送信エコーに変換し、該送信エコーを検体に対して送信するとともに、前記検体によって反射された受信エコーを受信して電気的な受信信号に変換する超音波探触子の種類に応じた特性を与えるパラメータであって前記送信駆動波の生成および前記受信信号の補正に必要なパラメータを記憶する記憶部から前記パラメータを参照し、前記超音波探触子の種類に応じた前記送信駆動波を送信駆動波生成部によって生成する送信駆動波生成ステップと、
    　前記記憶部が記憶する前記パラメータを参照し、前記受信信号に対して前記超音波探触子の種類に応じた補正を行う受信信号補正ステップと、
    　前記受信信号補正ステップで補正した前記受信信号を用いて画像データを生成する画像処理ステップと、
    　をコンピュータに実行させることを特徴とする超音波観測装置の作動プログラム。
PCT/JP2014/063395 2013-07-18 2014-05-20 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム WO2015008534A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2015503710A JP5788624B2 (ja) 2013-07-18 2014-05-20 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
CN201480028486.8A CN105392427B (zh) 2013-07-18 2014-05-20 超声波观测装置以及超声波观测装置的动作方法
EP14826109.2A EP3023058B1 (en) 2013-07-18 2014-05-20 Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and program for operating ultrasound observation apparatus
US14/942,163 US9662090B2 (en) 2013-07-18 2015-11-16 Ultrasonic observation apparatus, operation method of ultrasonic observation apparatus, and computer readable recording medium

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2013149831 2013-07-18
JP2013-149831 2013-07-18

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/942,163 Continuation US9662090B2 (en) 2013-07-18 2015-11-16 Ultrasonic observation apparatus, operation method of ultrasonic observation apparatus, and computer readable recording medium

Publications (1)

Publication Number Publication Date
WO2015008534A1 true WO2015008534A1 (ja) 2015-01-22

Family

ID=52346009

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/063395 WO2015008534A1 (ja) 2013-07-18 2014-05-20 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム

Country Status (5)

Country Link
US (1) US9662090B2 (ja)
EP (1) EP3023058B1 (ja)
JP (1) JP5788624B2 (ja)
CN (1) CN105392427B (ja)
WO (1) WO2015008534A1 (ja)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107530057A (zh) * 2015-05-13 2018-01-02 奥林巴斯株式会社 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序
WO2018186473A1 (ja) * 2017-04-05 2018-10-11 オリンパス株式会社 超音波観測装置の作動方法、超音波観測装置および超音波観測装置の作動プログラム
EP3275376A4 (en) * 2015-03-23 2019-01-16 Olympus Corporation ULTRASONIC OBSERVATION DEVICE, ULTRASONIC OBSERVATION DEVICE OPERATING PROCEDURE AND ULTRASONIC OBSERVATION DEVICE OPERATING PROGRAM
EP3295876A4 (en) * 2015-05-13 2019-02-20 Olympus Corporation ULTRASONIC OBSERVATION DEVICE, OPERATING METHOD FOR ULTRASONIC OBSERVATION DEVICE AND OPERATING PROGRAM FOR ULTRASONIC OBSERVATION DEVICE
US11207056B2 (en) 2015-05-13 2021-12-28 Olympus Corporation Ultrasound diagnostic apparatus, method for operating ultrasound diagnostic apparatus, and computer-readable recording medium

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106999150B (zh) * 2014-11-25 2021-03-30 皇家飞利浦有限公司 用于监控医学设备的使用的系统
EP3480621B1 (en) * 2017-11-01 2022-05-11 Samsung Medison Co., Ltd. Ultrasound apparatus and control method thereof

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02279145A (ja) * 1989-04-21 1990-11-15 Olympus Optical Co Ltd 超音波診断装置
JPH0615605U (ja) * 1992-08-04 1994-03-01 横河メディカルシステム株式会社 Ｃｆｍ分散表示機能付超音波ドプラ診断装置
JP2000279408A (ja) * 1999-03-30 2000-10-10 Aloka Co Ltd 超音波診断装置における送信回路及び波形生成方法
JP2003325506A (ja) * 2002-05-16 2003-11-18 Aloka Co Ltd 超音波診断装置及び送信信号設計方法
JP2005270247A (ja) 2004-03-24 2005-10-06 Aloka Co Ltd 超音波診断装置
WO2012063929A1 (ja) * 2010-11-11 2012-05-18 オリンパスメディカルシステムズ株式会社 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP2012250080A (ja) 2012-09-24 2012-12-20 Ge Medical Systems Global Technology Co Llc 超音波診断装置、超音波探触子及び信号処理方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3809880B2 (ja) * 1997-02-27 2006-08-16 株式会社ミツトヨ 振動制御系における伝達関数の短時間内取得方法および装置
JP4342859B2 (ja) * 2002-09-30 2009-10-14 富士フイルム株式会社 超音波用探触子及びそれを用いた超音波送受信装置
US7074186B2 (en) * 2003-09-23 2006-07-11 Siemens Medical Solutions Usa, Inc. Transmit based axial whitening
US7481769B2 (en) * 2003-09-30 2009-01-27 Fujifilm Corporation Ultrasonic diagnosing apparatus
JP5159480B2 (ja) * 2008-07-11 2013-03-06 株式会社東芝 超音波診断装置および超音波診断装置の制御プログラム
CN102596048B (zh) * 2010-10-26 2014-12-10 株式会社东芝 超声波诊断装置、超声波图像处理装置、医用图像诊断装置及医用图像处理装置
JPWO2012063928A1 (ja) * 2010-11-11 2014-05-12 オリンパスメディカルシステムズ株式会社 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
WO2012063976A1 (ja) * 2010-11-11 2012-05-18 オリンパスメディカルシステムズ株式会社 超音波診断装置、超音波診断装置の作動方法および超音波診断装置の作動プログラム
JP2012161562A (ja) * 2011-02-09 2012-08-30 Fujifilm Corp 超音波診断装置および超音波画像生成方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02279145A (ja) * 1989-04-21 1990-11-15 Olympus Optical Co Ltd 超音波診断装置
JPH0615605U (ja) * 1992-08-04 1994-03-01 横河メディカルシステム株式会社 Ｃｆｍ分散表示機能付超音波ドプラ診断装置
JP2000279408A (ja) * 1999-03-30 2000-10-10 Aloka Co Ltd 超音波診断装置における送信回路及び波形生成方法
JP2003325506A (ja) * 2002-05-16 2003-11-18 Aloka Co Ltd 超音波診断装置及び送信信号設計方法
JP2005270247A (ja) 2004-03-24 2005-10-06 Aloka Co Ltd 超音波診断装置
WO2012063929A1 (ja) * 2010-11-11 2012-05-18 オリンパスメディカルシステムズ株式会社 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP2012250080A (ja) 2012-09-24 2012-12-20 Ge Medical Systems Global Technology Co Llc 超音波診断装置、超音波探触子及び信号処理方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3023058A4

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3275376A4 (en) * 2015-03-23 2019-01-16 Olympus Corporation ULTRASONIC OBSERVATION DEVICE, ULTRASONIC OBSERVATION DEVICE OPERATING PROCEDURE AND ULTRASONIC OBSERVATION DEVICE OPERATING PROGRAM
CN107530057A (zh) * 2015-05-13 2018-01-02 奥林巴斯株式会社 超声波诊断装置、超声波诊断装置的工作方法以及超声波诊断装置的工作程序
EP3295876A4 (en) * 2015-05-13 2019-02-20 Olympus Corporation ULTRASONIC OBSERVATION DEVICE, OPERATING METHOD FOR ULTRASONIC OBSERVATION DEVICE AND OPERATING PROGRAM FOR ULTRASONIC OBSERVATION DEVICE
EP3295875A4 (en) * 2015-05-13 2019-02-20 Olympus Corporation ULTRASOUND DIAGNOSIS DEVICE, METHOD FOR OPERATING AN ULTRASONIC DIAGNOSIS DEVICE AND PROGRAM FOR OPERATING AN ULTRASOUND DIAGNOSTIC DEVICE
US11207056B2 (en) 2015-05-13 2021-12-28 Olympus Corporation Ultrasound diagnostic apparatus, method for operating ultrasound diagnostic apparatus, and computer-readable recording medium
WO2018186473A1 (ja) * 2017-04-05 2018-10-11 オリンパス株式会社 超音波観測装置の作動方法、超音波観測装置および超音波観測装置の作動プログラム

Also Published As

Publication number Publication date
CN105392427A (zh) 2016-03-09
JPWO2015008534A1 (ja) 2017-03-02
CN105392427B (zh) 2018-02-13
JP5788624B2 (ja) 2015-10-07
EP3023058A4 (en) 2017-03-29
US20160074008A1 (en) 2016-03-17
US9662090B2 (en) 2017-05-30
EP3023058B1 (en) 2023-07-12
EP3023058A1 (en) 2016-05-25

Similar Documents

Publication Publication Date Title
JP5788624B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5568199B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5642910B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5659324B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5430809B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5974210B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP5881905B2 (ja) 超音波観測装置、超音波観測装置の作動方法、及び超音波観測装置の作動プログラム
JP6289772B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
US9517054B2 (en) Ultrasound observation apparatus, method for operating ultrasound observation apparatus, and computer-readable recording medium
WO2016181869A1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP2018191779A (ja) 超音波観測装置
JP5981072B1 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム
JP6138402B2 (ja) 超音波観測装置、超音波観測装置の作動方法および超音波観測装置の作動プログラム

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 201480028486.8

Country of ref document: CN

ENP Entry into the national phase

Ref document number: 2015503710

Country of ref document: JP

Kind code of ref document: A

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14826109

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2014826109

Country of ref document: EP