WO2005060832A1 - 超音波診断装置および超音波イメージング方法 - Google Patents

超音波診断装置および超音波イメージング方法 Download PDF

Info

Publication number
WO2005060832A1
WO2005060832A1 PCT/JP2004/018260 JP2004018260W WO2005060832A1 WO 2005060832 A1 WO2005060832 A1 WO 2005060832A1 JP 2004018260 W JP2004018260 W JP 2004018260W WO 2005060832 A1 WO2005060832 A1 WO 2005060832A1
Authority
WO
WIPO (PCT)
Prior art keywords
ultrasonic
correction
subject
depth
dynamic range
Prior art date
Application number
PCT/JP2004/018260
Other languages
English (en)
French (fr)
Inventor
Mitsuhiro Oshiki
Akira Sasaki
Shinji Kishimoto
Original Assignee
Hitachi Medical Corporation
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 Hitachi Medical Corporation filed Critical Hitachi Medical Corporation
Priority to JP2005516451A priority Critical patent/JP4664209B2/ja
Publication of WO2005060832A1 publication Critical patent/WO2005060832A1/ja

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves

Definitions

  • the present invention relates to an apparatus and a method for capturing an ultrasonic image as a diagnostic image of a subject.
  • An ultrasonic diagnostic apparatus transmits and receives an ultrasonic beam to and from a subject using an ultrasonic probe, reconstructs an ultrasonic image based on a received signal output from the ultrasonic probe, and reconstructs an ultrasonic image.
  • the configured ultrasonic image is displayed on the display screen.
  • the ultrasonic beam is formed by a plurality of transducers arranged in an ultrasonic probe, having a focal point by transmitted or received ultrasonic waves.
  • the reflected echo generated by the subject force is attenuated in the process of propagating inside the subject.
  • the image quality of the ultrasound image affected by this attenuation is degraded, for example, a signal as diagnostic information is buried in noise. Therefore, it has been practiced to set the correction data exponentially as a function of the depth of the subject and to correct the received signal based on the set correction data, thereby improving the image quality of the ultrasonic image (for example, see Patent Document 1).
  • Patent Document 1 Japanese Patent Application Laid-Open No. Hei 4 40945
  • An object of the present invention is to realize an ultrasonic diagnostic apparatus and an ultrasonic imaging method that are more suitable for improving the image quality of an ultrasonic image.
  • An ultrasonic diagnostic apparatus of the present invention provides an ultrasonic probe that transmits and receives ultrasonic waves to and from a subject, a transmission unit that supplies a drive signal to the ultrasonic probe, and an ultrasonic probe.
  • Receiving means for processing a received signal output from the receiver, and a supersonic based on the received signal output from the receiving means.
  • Display means for displaying a wave image, and means for correcting a reception signal output from the reception means in accordance with the beam shape of the ultrasonic beam transmitted and received by the ultrasonic probe and the depth of the subject It is characterized by having.
  • the ultrasonic beam of at least one of the ultrasonic transmitting beam and the ultrasonic receiving beam may have various factors such as the position of the focal point and the number of transducers forming the diameter of the ultrasonic probe.
  • the beam shape is determined by the element.
  • Such an ultrasonic beam has a shape in which the aperture force is narrowed down to the focal point, and expands after the focal point.
  • the ultrasonic beam has a different beam cross-sectional area in the depth direction of the subject, the beam intensity varies in the depth direction of the subject. Due to this, the received signal output from the ultrasonic probe is affected by both the shape of the ultrasonic beam and the depth of the subject.
  • the reception signal after the correction becomes the same as the shape of the ultrasonic beam.
  • the signal is less affected by the depth of the subject.
  • the subject power can be reduced.
  • a faithfully displayed ultrasonic image can be displayed.
  • the correcting means corrects the dynamic range of the received signal from which the power of the receiving means is also output, based on the dynamic range calculated according to the beam shape of the ultrasonic beam and the depth of the subject.
  • Dynamic range correction means can be used.
  • the signal intensity of the reception signal output from the reception means is based on the beam intensity calculated corresponding to the beam shape of the ultrasonic beam and the depth of the subject.
  • the receiving means may include means for converting a received signal output from the ultrasonic probe into a digital signal based on a sampling clock.
  • the correction means may include one or more sensors set in the depth direction of the subject based on the sampling clock. For each sampling point, the converting means can also correct the received signal output.
  • the correction means may include a calculation means for calculating correction data.
  • the calculating means calculates the frequency of the drive signal, the wave number of the drive signal, the coordinates of the focal point of the ultrasonic beam, the number of transducers forming the aperture of the ultrasonic probe, and the focus data for forming the ultrasonic beam.
  • Correction data can be calculated based on parameters including at least one.
  • the correction means may include a calculation means for calculating correction data.
  • the calculating means calculates one correction data corresponding to the beam shape of one ultrasonic beam and the depth of the subject, and calculates one correction data while a reception signal corresponding to one ultrasonic beam is received.
  • Other correction data can be calculated according to the beam shape of another ultrasonic beam different from the ultrasonic beam and the depth of the subject.
  • the dynamic range calculated in accordance with the beam shape of the ultrasonic beam can be displayed as a graph as a function of the depth of the subject.
  • the beam intensity calculated corresponding to the beam shape of the ultrasonic beam can be displayed as a graph as a function of the depth of the subject.
  • An ultrasonic imaging method includes a step of calculating correction data corresponding to a beam shape of an ultrasonic beam formed by an ultrasonic probe and a depth of a subject; Transmitting a driving signal to the ultrasonic probe force to transmit an ultrasonic wave to the subject; and a receiving step of processing a reception signal output from the ultrasonic probe in response to the transmitting step.
  • Receiving step power The method includes a step of correcting an output received signal based on correction data, and a step of displaying an ultrasonic image based on the received signal output from the correction step.
  • a dynamic range corresponding to the beam shape of the ultrasonic beam and the depth of the subject can be calculated.
  • the dynamic range of the reception signal output from the reception step can be corrected based on the dynamic range calculated in the correction data calculation step.
  • a beam intensity corresponding to the beam shape of the ultrasonic beam and the depth of the subject can be calculated.
  • the signal strength of the received signal from which the reception step power is also output is calculated. Can be corrected.
  • correction data calculation step correction data corresponding to at least one of the ultrasonic transmission beam and the ultrasonic reception beam can be calculated.
  • the received signal output from the ultrasonic probe can be converted into a digital signal based on a sampling clock.
  • the reception signal output at the reception step can be corrected for each of one or more sampling points set in the depth direction of the subject based on the sampling clock.
  • the frequency of the drive signal, the wave number of the drive signal, the coordinates of the focal point of the ultrasonic beam, the number of transducers forming the aperture of the ultrasonic probe, and the ultrasonic beam are calculated.
  • the correction data can be calculated based on a parameter including at least one of the focus data to be formed.
  • a step of calculating one correction data corresponding to a beam shape of one ultrasonic beam and a depth of the subject, and a step of receiving a signal corresponding to one ultrasonic beam are performed. While receiving the signal, the method may include calculating another correction data corresponding to the beam shape of another ultrasonic beam different from the one ultrasonic beam and the depth of the subject. .
  • the dynamic range calculated according to the beam shape of the ultrasonic beam can be displayed as a graph as a function of the depth of the subject.
  • the beam intensity calculated corresponding to the beam shape of the ultrasonic beam can be displayed as a graph as a function of the depth of the subject.
  • FIG. 1 is a block diagram showing a configuration of an ultrasonic diagnostic apparatus according to a first embodiment to which the present invention has been applied.
  • FIG. 2 is a flowchart showing an entire process of the ultrasonic diagnostic apparatus of FIG. 1.
  • FIG. 3 is a conceptual diagram for explaining an operation of a dynamic range correction unit in FIG. 1.
  • FIG. 4 is a display example of an ultrasonic image captured by the ultrasonic diagnostic apparatus of FIG. 1 and a comparative example captured by the apparatus of the reference embodiment.
  • FIG. 5 is a diagram illustrating that the beam shape of an ultrasonic beam differs in the depth direction.
  • FIG. 6 is a block diagram showing a configuration of an ultrasonic diagnostic apparatus according to a second embodiment to which the present invention is applied.
  • FIG. 7 is a flowchart showing an entire process of the ultrasonic diagnostic apparatus of FIG. 6.
  • FIG. 8 is a conceptual diagram for explaining the operation of the signal strength correction unit in FIG. 6.
  • FIG. 9 is a block diagram showing a configuration of an ultrasonic diagnostic apparatus according to a third embodiment to which the present invention has been applied.
  • FIG. 10 is a flowchart showing an entire process of the ultrasonic diagnostic apparatus of FIG. 9.
  • FIG. 11 is a conceptual diagram for explaining the operation of the correction unit in FIG. 9.
  • FIG. 1 is a block diagram showing a configuration of an ultrasonic diagnostic apparatus 1 of the present embodiment.
  • the ultrasonic diagnostic apparatus 1 includes an ultrasonic probe 10 for transmitting and receiving ultrasonic waves to and from a subject, and a transmission unit for supplying a drive signal to the ultrasonic probe 10. And a receiving means for processing a received signal output from the ultrasonic probe 10.
  • the receiving means here includes a receiving circuit 16 that performs processing such as amplification on a received signal output from the ultrasonic probe 10, and converts the received signal output from the receiving circuit 16 into a digital signal based on a sampling clock.
  • An analog-to-digital conversion unit 18 (hereinafter referred to as an AZD conversion unit 18) and a phasing addition unit 20 as software phasing addition means for performing a phasing addition process on a received signal output from the AZD conversion unit 18.
  • a signal processing unit 22 that performs processing such as detection on the received signal output from the phasing addition unit 20.
  • the ultrasonic diagnostic apparatus 1 based on correction data calculated in correlation with the beam shape of the ultrasonic beam transmitted / received by the ultrasonic probe 10 and the depth of the subject, And a dynamic range correction unit 24 as means for correcting the received signal output from the signal processing unit 22.
  • the ultrasonic beam in the present embodiment means an ultrasonic wave transmitting beam and an ultrasonic wave receiving beam, but either one may be used.
  • the ultrasonic probe 10 is formed by arranging a plurality of transducers. Each transducer is supplied from the transmission circuit 12
  • the drive signal is converted into an ultrasonic wave, and the reflected echo generated by the subject force is converted into a received signal.
  • the drive signal here is a pulse signal for transmitting ultrasonic waves.
  • the transmission circuit 12 is provided in the transmission / reception unit 26 together with the reception circuit 16.
  • the transmission / reception unit 26 outputs a drive signal output from the transmission circuit 12 to the ultrasonic probe 10, and a transmission / reception separation circuit that outputs a reception signal output from the ultrasonic probe 10 to the reception circuit 16.
  • the ultrasonic probe 10 has a high-voltage switching switch for selecting a group of transducers that form a transmission aperture or a reception aperture of the ultrasonic probe 10.
  • the receiving circuit 16 has a receiving amplifier that amplifies a received signal output from the ultrasonic probe 10.
  • the receiving amplifier has a preamplifier that amplifies the received signal with a predetermined signal amplification factor, and a gain control amplifier that varies the signal amplification factor in accordance with the attenuation of the reflected echo caused by the depth of the ultrasonic reflection portion.
  • a focus position determination unit 28 is connected to the transmission / reception unit 26.
  • the focus position determination unit 28 outputs the set coordinates of the focal point of the ultrasonic transmission beam to the transmission / reception unit 26.
  • a keyboard or mouse for inputting the coordinates of the focal point may be provided in the focus position determining unit 28.
  • the AZD converter 18 has a clock generator that generates a sampling clock.
  • the phasing addition section 20 includes a delay amount correction section 34 for digitally phasing the received signal output from the AZD conversion section 18 based on the focus data, and a memory 36 storing focus data provided to the delay amount correction section 34.
  • an adder 38 that adds the received signals for each channel output from the delay amount corrector 34 and outputs the result to the signal processor 22.
  • the signal processing unit 22 has a function of performing processing such as detection, filtering, and logarithmic compression on the reception signal output from the phasing addition unit 20, and an ultrasonic image (based on the processed reception signal). For example, it has a function of forming an image signal that is a basis of a B-mode image, a Doppler image, and an M-mode image).
  • the image signal output from the signal processing unit 22 is sent to the dynamic range correction unit 24. Is output.
  • the image signal output from the signal processing unit 22 is also appropriately referred to as a received signal.
  • the dynamic range correction unit 24 has a dynamic range calculation unit 40 and a memory 42.
  • the dynamic range calculation unit 40 calculates a dynamic range (hereinafter referred to as a correction dynamic range) as correction data correlated with the beam shape of the ultrasonic beam transmitted and received by the ultrasonic probe 10 and the depth of the subject. It has a function to do. For example, the frequency of the drive signal, the wave number of the drive signal, the coordinates of the focal point of the ultrasonic transmission beam or the ultrasonic reception beam, the size of the aperture of the ultrasonic probe 10, the ultrasonic transmission beam or the reception beam
  • the correction dynamic range is calculated based on a parameter including at least one of the focus data for forming the correction data.
  • the size of the aperture is obtained based on the number of transducers forming the transmission aperture or the reception aperture, and the size of each transducer.
  • the dynamic range calculation unit 40 arranges the dynamic range for correction in a dynamic range table.
  • the dynamic range table is a database in which correction dynamic ranges are arranged in correspondence with the depth of the subject (for example, the depth of a sampling point or a focus stage).
  • the dynamic range table has an ultrasound beam profile as a function of the depth of focus.
  • Such a dynamic range table is stored in the memory 42.
  • the dynamic range calculator 40 corrects the dynamic range of the received signal output from the signal processor 22 based on the correction dynamic range. For example, based on the sampling clock of the AZD conversion unit 18, the reception signal output from the signal processing unit 22 is corrected for each sampling point set in the depth direction of the subject.
  • the received signal instead of correcting the received signal for each sampling point, the received signal may be corrected for each focus stage having a plurality of sampling points.
  • the dynamic range correction unit 24 calculates the first correction dynamic range based on the beam shape of the first ultrasonic beam and the depth of the subject, and corresponds to the first ultrasonic beam.
  • the second dynamic range for correction can be calculated based on the beam shape of the second ultrasonic beam different from the first ultrasonic beam and the depth of the subject. That is, the dynamic range calculation unit 40 can calculate the dynamic range for correction in advance by time division processing. Or, the shape of the ultrasonic beam set in advance On the other hand, the corresponding dynamic range for correction can be stored in the dynamic range table 42 in advance.
  • a digital scan converter 30 (hereinafter, referred to as DSC30) that converts the corrected reception signal output from the dynamic range correction unit 24 into a display signal and performs luminance matching processing on the reception signal.
  • a display unit 32 for displaying an ultrasonic image based on the received signal output from the DSC 30 on a display screen.
  • the ultrasonic probe 10 is brought into contact with the subject.
  • the drive signal is supplied from the transmission circuit 12 to the ultrasonic probe 10
  • the ultrasonic probe 10 irradiates the subject with an ultrasonic transmission beam.
  • the reflected echo generated from the subject is converted into a received signal by being received by the ultrasonic probe 10.
  • the received signal is amplified by the receiving circuit 16.
  • the amplified received signal is converted into a digital signal by the AZD converter 18.
  • the digitally received signal is formed by the phasing adder 20 as an ultrasonic wave receiving beam.
  • the ultrasonic reception beam output from the phasing addition section 20 is output as an image signal by being subjected to predetermined processing by the signal processing section 22.
  • the dynamic range of the image signal output from the signal processing unit 22 is corrected by the dynamic range correction unit 24.
  • the corrected image signal is displayed as an ultrasonic image on the display unit 32 via the DSC 30.
  • FIG. 2 is a flowchart showing the overall processing of the ultrasonic diagnostic apparatus 1 of FIG.
  • FIG. 3 is a conceptual diagram for explaining the operation of the dynamic range correction unit in FIG. Although FIG. 3 shows processing in one focus stage for convenience, the same applies to processing in other focus stages.
  • the dynamic range for correction D1 to D6 is calculated by the dynamic range calculation unit 40.
  • the dynamic range for correction D1 to D6 is correction data correlated with the beam shape of the ultrasonic beam formed by the ultrasonic probe 10 and the depth of each of the reflection parts 46-1 to 46-7.
  • the dynamic range calculation unit 40 also calculates the correction dynamic ranges D1 to D6 for the factor power related to the beam shape of the ultrasonic beam.
  • the factors are the frequency of the drive signal supplied from the transmission circuit 12, the wave number of the drive signal, the coordinates of the focal point of the ultrasonic beam, the number of transducers forming the aperture of the ultrasonic probe 10, the phasing method, and the like. Is the parameter.
  • the phasing method is transmission focus data for delaying the drive signal, and reception focus data for phasing the reception signal (for example, focus data stored in the memory 36).
  • a general formula for setting the correction dynamic range D1 to D6 is represented by the following formula (1).
  • the dynamic range is set as shown in the table.
  • the correction dynamic range Dl-D6 is a unique dynamic range of the ultrasonic beam, and depends on the beam intensity at each of the reflection portions 46-1-46-7.
  • Such correction dynamic ranges D1 to D6 are stored in the memory 42 as a dynamic range table after being associated with the sampling point or the depth of the focus step.
  • FIG. 3C shows a dynamic range for correction D1 to D6 corresponding to the received signals S1 to S6.
  • the dynamic range for correction D1 to D6 may be obtained from an actual measurement and a simulation.
  • a drive signal is supplied to the ultrasonic probe 10 to transmit ultrasonic waves from the ultrasonic probe 10 to the subject.
  • a drive signal is supplied from the wave transmitting circuit 12 to the ultrasonic probe 10 that has come into contact with the body surface of the subject.
  • a plurality of drive signals are generated corresponding to the number of transducers included in the ultrasonic probe 10, and each of the generated drive signals is delayed by using transmission focus data based on the focal position. I have.
  • the focus position is output from the focus position determination unit 28 to the transmission / reception unit 26 at a predetermined timing.
  • a predetermined transducer group is selected by a high-voltage switching switch, thereby forming a transmission aperture.
  • a drive signal is supplied to a group of transducers having a transmission aperture, ultrasonic waves are transmitted from the supplied transducers.
  • Each transmitted ultrasonic wave forms an ultrasonic transmission beam.
  • the formed ultrasonic transmission beam has a focal point at a predetermined distance in the transmission direction.
  • the vertical axis in FIG. 3A indicates the depth of the subject, and the horizontal axis indicates the scanning address of the beam.
  • a plurality of reflection portions 46-1-1 46-7 are scattered on the subject.
  • Each reflection part 46-1—46—7 has a different size and acoustic impedance.
  • the reflecting portions 46-1 to 46-7 have different depths except for the reflecting portion 46-6.
  • the reflection part 46-6 has the same depth as the reflection part 46-1.
  • the beam address of the ultrasonic probe 10 is initialized to the address A shown in FIG. 3A. Then, an ultrasonic transmission beam is emitted from the ultrasonic probe 10. In other words, the ultrasonic transmission beam is sequentially emitted to the reflecting portions 46-1 to 46-7 while changing the address in the scanning direction with the address A force.
  • FIG. 3A shows a two-dimensional distribution in the depth direction and the scanning direction, the same applies to a three-dimensional distribution.
  • the focal position of the ultrasonic transmission beam can be changed as appropriate.
  • the reception signal output from the ultrasonic probe 10 is processed in accordance with the ultrasonic wave transmitting step (S101).
  • the reflected echo generated from the reflecting portion 46-1—46-7 is received by the ultrasonic probe 10 and converted into a received signal.
  • the converted reception signal is sequentially processed by the reception circuit 16, the AZD conversion unit 18, and the phasing addition unit 20, whereby an ultrasonic wave reception beam is formed.
  • the formed ultrasonic receiving beam is focused at a predetermined depth of the subject.
  • the ultrasonic receiving beam is output to the dynamic range correction unit 24 as received signals S1 to S6 as shown in FIG. 3B.
  • the vertical axis in FIG. 3B indicates the depth of the subject, and the horizontal axis indicates the signal intensity.
  • the signal strength of the received signals S1 to S6 corresponds to, for example, luminance information when an image signal is formed.
  • FIG. 3B shows the signal intensities of the received signals S1 to S6 as projected for each depth.
  • the reflection part 46-1 and the reflection part 46-6 have the same depth, since the signal strength of the reception signal corresponding to the reflection part 46-1 is larger, the reception signal corresponding to the reflection part 46-1 is larger. As the received signal S 1.
  • the reception signal SI-S6 output from the ultrasonic receiving step (S102) is corrected based on the dynamic range for correction.
  • the received signal S1 is input to the dynamic range correction unit 24.
  • the dynamic range calculator 40 reads the dynamic range D1 for correction from the memory 42.
  • the dynamic range D1 for correction is arranged in the dynamic range table, and the dynamic range table power is also extracted corresponding to the depth of the reflection part 46-1.
  • the signal strength of the received signal S1 is corrected by the dynamic range calculation unit 40 based on the correction dynamic range D1.
  • FIG. 3D shows the corrected received signals DN1 to DN6 processed in this way.
  • the received signals DN1 to DN6 are signals in which the inherent dynamic range of the ultrasonic beam is considered.
  • Equation 2 The correction calculation in this step is performed based on, for example, Equation 2.
  • Equation 2 k indicates the number of received signals to be corrected. In this embodiment, the value of k is 1 to 6. The force may be changed as appropriate.
  • a indicates a coefficient for correcting the signal intensity attenuated when the reflected echo propagates in the subject.
  • ⁇ Ultrasonic image display process (SI 04)> An ultrasonic image based on the received signals DN1 to DN6 output from the correction step (S103) is displayed on the display screen.
  • the received signals DN1 to DN6 are converted into display signals by the DSC 30, and are subjected to luminance adjustment processing.
  • the received signals DN1 to DN6 output from the DSC 30 are displayed on the display screen of the display unit 32 as ultrasonic images.
  • the dynamic range for correction D1 to D6 is displayed as a graph alongside the ultrasonic image as a function of the depth of the reflection area 46-1 to 46-7.
  • the dynamic range D1-D6 for correction is visually grasped by the graph display.
  • the ultrasonic wave transmitting beam or the ultrasonic wave receiving beam has a shape which is narrowed down to the focal point focal point, and spreads out after the focal point.
  • a difference occurs in sound pressure or reflection intensity depending on the depth of each of the reflection portions 46-1-1-46-7.
  • the ultrasonic beam has a different beam cross-sectional area in the depth direction of the subject, the beam intensity differs in the depth direction of the subject. Due to this, the received signal output from the signal processing unit 22 is influenced by both the shape of the ultrasonic beam and the depth of the reflection part 46-1 to 46-7.
  • the reception signals S1-S6 output from the signal processing section 22 are corrected in correlation with the shape of the ultrasonic beam and the depth of the reflection portions 46-1-146-7.
  • the corrected reception signals DN1 to DN6 are signals in which the influence of the shape of the ultrasonic beam and the depth of the reflection portion 46-1 to 46-7 is reduced.
  • FIG. 4A is a display example of an ultrasonic image captured by a phantom by the ultrasonic diagnostic apparatus of the present embodiment.
  • FIG. 4B is a display example of an ultrasonic image captured by a phantom using the ultrasonic diagnostic apparatus according to the reference embodiment.
  • the dynamic range at the time of displaying an ultrasonic image is set to be constant (for example, 100 dB) regardless of the depth of the force reflection portion.
  • a reflected echo generated by a shallow reflecting portion force has a large signal strength
  • a reflected echo generated by a deep reflecting portion force has a small signal intensity due to attenuation in the subject. Therefore, if the image is displayed according to the dynamic range of the received signal corresponding to the shallow reflecting portion, the contrast of the image corresponding to the deep reflecting portion may be reduced.
  • the information corresponding to the shallow reflection part may be lost due to luminance saturation.
  • the ultrasonic beam having the widest dynamic range at each depth of the reflection portion is used as a reference.
  • the dynamic range in which all the received signals can be expressed is set, so that the display is wider than the dynamic range inherent in the received signals.
  • unnecessary components are displayed on the ultrasonic image as noise.
  • the dynamic range of the received signal is set in correlation with the beam shape of the ultrasonic beam. That is, an ultrasonic image faithful to the dynamic range of the ultrasonic beam at each depth of the reflection portions 46-1-147-7 is displayed. Therefore, the ultrasonic image shown in FIG. 4A is compared with FIG. 4B in that there is no loss of a signal as diagnostic information, and unnecessary components that become noise are reduced. Such an effect becomes more remarkable in a portion where the inherent dynamic range of the ultrasonic beam is narrowed. For example, it is remarkable in a shallow part when a beam is formed using a variable aperture.
  • the corrected received signals DN1 to DN6 are weighted by a function according to the dynamic range D1 to D6 inherent to the ultrasonic beam. Based on the received signals DN1 to DN6 after such correction, the dynamic range for displaying an ultrasonic image is set for each depth. Therefore, the ultrasonic image displayed on the display unit 32 has a reduced noise component and more faithfully represents the difference in acoustic impedance of the reflection portions 46-1 to 46-7.
  • the present embodiment has been made by paying attention to the fact that the ultrasonic beam shape changes in the depth direction of the subject. That is, the ultrasonic beam shape changes in the depth direction of the subject due to the focal position of the ultrasonic transmission beam, the size of the transmission aperture, the size of the reception aperture, the phasing method of the ultrasonic reception beam, etc. . Therefore, the intensity of the ultrasonic beam (the magnitude of the sound field) differs depending on the depth direction. If the shape of the ultrasonic beam differs in the depth direction, the dynamic range of the ultrasonic beam changes accordingly.
  • the image quality of the ultrasonic image can be improved by appropriately setting the dynamic range of the image signal in the depth direction of the subject in correlation with the ultrasonic beam.
  • the dynamic range of the received signal may be relatively wide, for example, exceeding 70 dB.
  • the dynamic range when displayed as an ultrasonic image is set for each depth of the reflection site in correlation with the ultrasonic beam shape.
  • the shape of the ultrasonic beam changes at any time as the focal position of the ultrasonic beam or the phasing method changes.
  • the dynamic range table is automatically updated by the dynamic range calculation unit 40 following the change in the ultrasonic beam shape. Therefore, the usability and operability of the device are improved.
  • the present invention has been described with reference to the first embodiment, the present invention is not limited to this.
  • the correction processing shown in FIG. 2 may be executed in a unit of a beam address or in a unit of a frame.
  • the processing is performed on a frame basis, as described with reference to FIG. 3A, when there are a plurality of reflection portions (for example, the reflection portions 46-1, 46-6) at the same depth, the received signal strength is high. Should be used as the basis.
  • the present invention is not limited to capturing a B-mode image (tomographic image), and can be applied to capturing a Doppler image or an M-mode image.
  • the received signal input to the dynamic range correction unit 24 is an output signal of the signal processing unit 22. Instead, the received signal is output from the signal processing unit 22 or the output signal of the phasing addition unit 20. Or the output signal of the AZD converter 18. In short, the received signal output from the receiving means may be output to the dynamic range correction unit 24.
  • the dynamic range of the transmission / reception system circuit of the ultrasonic diagnostic apparatus 1 be wider than the inherent dynamic range of the ultrasonic beam transmitted and received by the ultrasonic probe 10.
  • FIG. 5A is a conceptual diagram of a plurality of sampling points T1 and TP (P: a natural number of 2 or more) set in the depth direction of the subject.
  • FIG. 5B is a diagram showing a beam shape of a reception beam formed at each of the focus stages A to D.
  • Figure 5C shows the beam shape of the ultrasonic transmission beam.
  • a plurality of sampling points T1 and TP are set in the depth direction of the subject. Each interval between the sampling points T1 and TP is set equal to the sampling clock interval of the AZD converter 18.
  • the plurality of sampling points T1 to TP are divided into focus stages A to D, for example, every four. For example, sampling points T1 to T4 belong to focus stage A. Note that the number of sampling points belonging to one focus stage is appropriately determined within a range that does not hinder diagnosis.
  • Common focus data is set for each of such focus stages A to D.
  • the set focus data is stored in the memory 36.
  • the reception focus data is a delay amount or a minute delay amount for phasing each reception signal output from each transducer of the ultrasonic probe 10, and is switched for each of the focus stages A to D.
  • the focus data stored in the memory 36 is read according to a control command. Based on the read focus data, the reception signal output from the AZD conversion unit 18 is delayed and phased by the phasing addition unit 20. As a result, ultrasonic receiving beams corresponding to the respective sampling points T1 and TP are respectively formed.
  • so-called dynamic focus is performed.
  • the dynamic focus is a method of focusing a beam in a relatively wide range in the depth direction by focusing on a part at a shallow depth, a part depth, and a part.
  • focus data here, independent focus data is set for each of the focus stages A to D within a range that does not affect the diagnosis.
  • the total number of focus data can be reduced as compared with the case where each sampling point T1 is set for each TP, and the storage area of the memory 36 can be used effectively.
  • the beam shape (beam profile) of the ultrasonic receiving beam differs in the depth direction of the subject due to the difference in the focus data for each of the focus stages A to D.
  • Fig. 5B shows the received beam shape corresponding to sampling point T1 belonging to focus stage A, the received beam shape corresponding to sampling point T5 belonging to focus stage B, and the sampling point belonging to focus stage C in order from the top. Show the received beam shape corresponding to T (P-8) and the received beam shape corresponding to sampling point ⁇ ( ⁇ -4) belonging to focus stage D! / ⁇ .
  • the ultrasonic wave reception beam has a beam shape in the depth direction of the subject. Different.
  • the depth differs slightly between a plurality of sampling points (for example, sampling points T1 and T4) belonging to the same focus stage (for example, focus stage A).
  • the beam shape of the wave beam is also slightly different.
  • the beam shape of the ultrasonic wave receiving beam also differs in the depth direction due to the variable aperture.
  • the variable aperture is a method of automatically reducing the aperture width to a shallower portion. With the variable aperture, it is possible to suppress the spread of the reflected echo generated at the sampling points (for example, sampling points T1 and T2) close to the ultrasonic probe 10. In the case of performing such a variable aperture, the size of the aperture for receiving the reflected echo, that is, the number of transducer elements forming the received aperture changes. Therefore, the beam shape of the ultrasonic receiving beam differs in the depth direction.
  • the beam shape of the ultrasonic transmission beam also differs in the depth direction.
  • an ultrasonic transmission beam is formed by focusing on one point.
  • the ultrasonic transmission beam is narrowed from the ultrasonic probe 10 to the focal point 41, and has a shape that spreads after the focal point.
  • an ultrasonic transmission beam having a focus set for each of the focus stages A to D may be independently formed and transmitted plural times.
  • the beam shape of the ultrasonic transmitting beam or the ultrasonic receiving beam differs in the depth direction of the subject. Therefore, due to the difference in ultrasonic beam intensity in the depth direction, the inherent dynamic range of the ultrasonic beam also differs in the depth direction. For example, at a depth where the ultrasonic beam intensity is large, the inherent dynamic range of the ultrasonic beam is widened. Ultrasonic beam intensity is small! / At the depth, the inherent dynamic range becomes narrow.
  • the present embodiment employs a dynamic range calculated based on the beam shape of the ultrasonic beam and the depth of the subject. Then, the dynamic range of the received signal output from the signal processing unit 22 is corrected.
  • a second embodiment of the ultrasonic diagnostic apparatus to which the present invention is applied will be described with reference to FIGS.
  • This embodiment correlates the beam shape of the ultrasonic beam with the depth of the subject.
  • the second embodiment differs from the first embodiment in that the received signal output from the signal processing unit 22 is corrected based on the calculated beam intensity. Therefore, portions that correspond to those of the first embodiment are denoted by the same reference numerals, and differences will be mainly described.
  • FIG. 6 is a block diagram showing a configuration of the ultrasonic diagnostic apparatus 2 of the present embodiment.
  • the ultrasonic diagnostic apparatus 2 includes a signal strength correction unit 50 instead of the dynamic range correction unit 24 of the first embodiment (for example, FIG. 1).
  • the signal intensity correction unit 50 receives and outputs the signal output from the signal processing unit 22 based on the ultrasonic beam calculated in correlation with the beam shape of the ultrasonic beam transmitted and received by the ultrasonic probe 10 and the depth of the subject. Correct the signal strength of the signal.
  • the signal strength correction unit 50 has a signal strength calculation unit 52 and a memory 54.
  • the signal intensity calculator 52 calculates a beam intensity as correction data correlating to the beam shape of the ultrasonic beam and the depth of the subject (hereinafter, appropriately referred to as a correction beam intensity), and calculates the correction beam intensity.
  • the signal strength of the received signal output from the signal processing unit 22 is corrected based on the signal strength.
  • the memory 54 stores a beam intensity table in which correction beam intensities are arranged as a function of the depth of the subject (for example, the depth of a sampling point or a focus step).
  • FIG. 7 is a flowchart showing the overall processing of the ultrasonic diagnostic apparatus 2 in FIG.
  • FIG. 8 is a conceptual diagram for explaining the operation of the signal strength correction unit 50 of FIG.
  • the present embodiment includes a correction beam intensity calculation step S200 instead of the correction dynamic range calculation step S100 of the first embodiment, and replaces the received signal correction step S103. And a receiving signal correcting step S203. Therefore, the description will focus on the correction beam intensity calculation step S200 and the received signal correction step S203.
  • the correction beam intensities A1 to A4 are calculated by the signal intensity calculator 52.
  • the correction beam intensities A1 to A4 are correction data correlated with the beam shape of the ultrasonic beam formed by the ultrasonic probe 10 and the depth of each reflection portion 46-1 to 46-7.
  • the beam intensity A1 corresponds to the beam intensity at the focal position of the ultrasonic beam.
  • the beam intensity A2 corresponds to the reflection portions 46-4 and 46-5.
  • the beam intensity A3 corresponds to the reflection portions 46-2 and 46-7.
  • Beam intensity A4 corresponds to the reflection area 46-1, 46-3, 46-6 To do.
  • the signal intensity calculation unit 52 calculates the factor correction beam intensity A1-A4 related to the beam shape of the ultrasonic beam.
  • the factors here include the frequency of the drive signal supplied from the transmission circuit 12, the number of waves of the drive signal, the coordinates of the focal point of the ultrasonic beam, the number of transducers forming the aperture of the ultrasonic probe 10, and the integer. It is a parameter such as a phase method.
  • the phasing method is determined by transmission focus data for delaying the drive signal and reception focus data (for example, focus data stored in the memory 36) for phasing the reception signal.
  • the correction beam intensities A1 to A4 are the signal intensities inherent to the ultrasonic beam.
  • Such correction beam intensities A1 to A4 are stored in the memory 54 as a beam intensity table after being associated with the sampling point or the depth of the focus stage.
  • the correction beam intensity is calculated based on the beam shape of the first ultrasonic beam and the depth of the subject, and the reception corresponding to the first ultrasonic beam is performed.
  • the correction beam intensity may be calculated in advance based on the beam shape of the second ultrasonic beam different from the first ultrasonic beam and the depth of the subject.
  • the beam intensity A1—A4 for force correction such as actual measurement or simulation may be obtained.
  • the received signals S1 to S6 output from the ultrasonic receiving step (S102) are corrected based on the correction beam intensities A1 to A4.
  • the signal strength calculation unit 52 reads them from the correction beam strengths A1 to A4 memory 54. Then, among the correction beam intensities A1 to A4, the correction beam intensity A1 is set as the reference intensity.
  • the reference beam intensity for correction may be changed as appropriate.
  • the relative ratio (for example, A1ZA4) between the correction beam intensity A1 as the reference intensity and the correction beam intensity A4 is obtained.
  • the obtained relative ratio corresponds to the amplitude ratio between the beam amplitude at the focal position of the ultrasonic beam and the beam amplitude at the depth of the reflecting portion 46-1.
  • Such a relative ratio is corrected Set as a coefficient.
  • the signal strength of received signal S1 is corrected by multiplying the signal strength of received signal S1 by a correction coefficient (A1ZA4). Similar processing is performed on the received signals S2 and S6. In short, the received signals S1 to S6 are corrected as received signals B1 to B6 in this step.
  • the corrected received signals B1 to B6 are signals that take into account the inherent beam intensity of the ultrasonic beam.
  • an ultrasonic image is a reflection part (for example, a reflection part) scattered in the depth direction of the subject.
  • the beam intensity of the ultrasonic beam differs in the depth direction in correlation with the beam shape.
  • the displayed ultrasonic image includes an error due to the intensity distribution of the transmitted beam in the depth direction.
  • the corrected received signals B1-B6 are signals whose signal intensities have been corrected based on the intensity distribution of the ultrasonic beam in the depth direction. That is, the received signals B1 to B6 after the correction become equivalent to the received signals when the beam intensity distribution in the depth direction of the ultrasonic beam is uniform.
  • a third embodiment of the ultrasonic diagnostic apparatus to which the present invention is applied will be described with reference to FIGS.
  • This embodiment is different from the first embodiment in that the received signal output from the signal processing unit 22 is corrected in the second embodiment as a first-stage correction, and then is corrected in the first embodiment as a second-stage correction.
  • This is different from the first and second embodiments. Therefore, the same reference numerals are given to the portions that correspond to the first and second embodiments, and the differences will be mainly described. I will tell.
  • FIG. 9 is a block diagram showing a configuration of the ultrasonic diagnostic apparatus 3 of the present embodiment.
  • the ultrasonic diagnostic apparatus 3 includes a dynamic range calculator 40 of the first embodiment (for example, FIG. 1) and a signal strength calculator 52 of the second embodiment (for example, FIG. 2). And a correction unit 56 having a memory 58.
  • the memory 58 stores a dynamic range table for correction and a beam intensity table for correction.
  • the correction unit 56 corrects correction data (for example, correction beam intensity A1 to A4, correction beam intensity A1 to A4) in correlation with the beam shape of the ultrasonic beam transmitted and received by the ultrasonic probe 10 and the depth of the subject.
  • correction data for example, correction beam intensity A1 to A4, correction beam intensity A1 to A4
  • the dynamic range D1 to D6) is calculated, the received signal output from the signal processing unit 22 is corrected based on the calculated correction data, and the corrected signal is output to the DSC 30.
  • FIG. 10 is a flowchart showing the overall processing of the ultrasonic diagnostic apparatus 3 in FIG.
  • FIG. 11 is a conceptual diagram for explaining the operation of the correction unit 56 in FIG.
  • the correction data calculation process according to the present embodiment includes a correction dynamic range calculation process S100 according to the first embodiment (for example, FIG. 2) and a correction dynamic range calculation process S100 according to the second embodiment (for example, FIG. 7).
  • a beam intensity calculation process S 200 is provided.
  • the correction dynamic ranges D1 to D6 are calculated and stored in the memory 58 as a dynamic range table
  • the correction beam intensities A1 to A4 are calculated and stored in the memory 58 as a signal intensity table. .
  • the step of correcting the received signal includes the step of correcting the received signal S203 of the second embodiment and the step of correcting the received signal S103 of the first embodiment.
  • the received signals S1 to S6 output from the signal processing unit 22 are subjected to the first-stage correction using the corrected beam intensities A1 to A4, thereby obtaining the received signals B1 to B6. become.
  • the received signals B1 to B6 are subjected to the second-stage correction using the corrected dynamic ranges D1 to D6, and thus become the received signals DN1 to DN6 after the second-stage correction.
  • the received signals S1-S6 output from the signal processing unit 22 are subjected to first-stage and second-stage corrections that are correlated with the shape of the ultrasonic beam and the depth of the subject.
  • the received signals DN1 to DN6 are signals in which the influence of the ultrasonic beam shape and the depth of the subject is reduced by one layer.
  • the dynamic range correction unit 24 and the signal intensity correction unit 50 serve as means for correcting the reception signal output from the signal processing unit 22. And at least one of the correction units 56!

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

 超音波診断装置1は、被検体との間で超音波を送受する超音波探触子10と、超音波探触子に駆動信号を供給する送信回路12と、超音波探触子10から出力される受信信号を処理する受信回路16と、A/D変換部18と、信号処理部22と、信号処理部22から出力される受信信号に基づいた超音波像が表示される表示部32を備え、超音波探触子10で送受される超音波ビームのビーム形状と被検体の深度とに対応して算出される補正データに基づき、信号処理部22から出力される受信信号を補正するダイナミックレンジ補正部24、信号強度補正部50、補正部56の少なくとも1つを有する。

Description

明 細 書
超音波診断装置および超音波イメージング方法
技術分野
[0001] 本発明は、被検体の診断画像として超音波像を撮像する装置およびその方法に関 する。
[0002] 超音波診断装置は、被検体との間で超音波探触子により超音波ビームを送受し、 超音波探触子から出力される受信信号に基づき超音波像を再構成し、再構成された 超音波象を表示画面に表示する。ここで超音波ビームとは、超音波探触子に配列さ れた複数の振動子力 送波又は受波される超音波により焦点を有して形成されるも のである。
[0003] ところで、被検体力 発生した反射エコーは、被検体内を伝播する過程で減衰する 。この減衰の影響を受けた超音波像は、診断情報としての信号がノイズに埋もれるな ど、画質が低下する。そこで、被検体の深度の関数として指数関数的に補正データ を設定し、設定された補正データに基づいて受信信号を補正することにより、超音波 像の画質を向上させることが行われている(例えば、特許文献 1参照)。
[0004] 特許文献 1:特開平 4 40945号公報
[0005] しかしながら、超音波診断装置により表示される超音波像は、被検体の深度以外の 要因によっても画質が影響を受ける。したがって、被検体の深度の関数として補正す る特許文献 1のような方法では、例えば、診断情報としての信号がノイズに埋もれた 超音波像や、画像コントラストが劣化した超音波像が表示されるという問題が依然とし て生じる場合があり、超音波像の画質向上に改善すべき余地がある。
[0006] 本発明の目的は、超音波像の画質を改善するのにより好適な超音波診断装置およ び超音波イメージング方法を実現することにある。
発明の開示
[0007] 本発明の超音波診断装置は、被検体との間で超音波を送受する超音波探触子と、 超音波探触子に駆動信号を供給する送信手段と、超音波探触子から出力される受 信信号を処理する受信手段と、受信手段から出力される受信信号に基づいた超音 波像が表示される表示手段を備え、超音波探触子で送受される超音波ビームのビー ム形状と被検体の深度とに対応して、受信手段から出力される受信信号を補正する 手段を有してなることを特徴とする。
[0008] すなわち、超音波送波ビーム又は超音波受波ビームの少なくとも一方の超音波ビ ームは、焦点の位置や、超音波探触子の口径を形成する振動子の数など種々の因 子によりビーム形状が決まってくる。このような超音波ビームは、口径力も焦点に至る まで絞られ、焦点を超えた後は広がる形状を有する。又、伝搬距離に応じて減衰が 生じるから、各反射部位の深度によって音圧又は反射強度に差が生じる。例えば、 超音波ビームは、被検体の深度方向にビーム断面積が異なるため、ビーム強度が被 検体の深度方向に異なる。これに起因して、超音波探触子から出力される受信信号 は、超音波ビームの形状と被検体の深度の双方による影響を受けたものになる。
[0009] したがって、受信手段から出力される受信信号に対し、超音波ビームの形状と被検 体の深度とに対応した補正を施すことにより、補正後の受信信号は、超音波ビームの 形状と被検体の深度による影響が低減された信号になる。このような受信信号に基 づく超音波像を表示することにより、超音波像の画質を改善できる。
[0010] 要するに、超音波ビームの形状が被検体の深度方向に異なることに着目し、超音 波ビーム形状と被検体の深度とに対応させて受信信号を補正することにより、被検体 力 り忠実に現された超音波像を表示できる。
[0011] この場合において、補正手段としては、超音波ビームのビーム形状と被検体の深度 とに対応して算出されたダイナミックレンジに基づき、受信手段力も出力される受信 信号のダイナミックレンジを補正するダイナミックレンジ補正手段を用いることができる 。また、ダイナミックレンジ補正手段に代えて、又はそれと共に、超音波ビームのビー ム形状と被検体の深度とに対応して算出されたビーム強度に基づき、受信手段から 出力される受信信号の信号強度を補正する信号強度補正手段を用いることができる
[0012] また、受信手段については、超音波探触子から出力される受信信号をサンプリング クロックに基づきディジタル信号に変換する手段を有することができる。この場合の補 正手段は、サンプリングクロックに基づき被検体の深度方向に設定された一以上のサ ンプリング点ごとに、変換手段力も出力される受信信号を補正できる。
[0013] また、補正手段については、補正データを算出する計算手段を有することができる 。計算手段は、駆動信号の周波数、駆動信号の波数、超音波ビームの焦点の座標、 超音波探触子の口径を形成する振動子の数、超音波ビームを形成するためのフォ 一カスデータの少なくとも 1つを含むパラメータに基づき補正データを算出できる。
[0014] さらに、補正手段については、補正データを算出する計算手段を有することができ る。計算手段は、一の超音波ビームのビーム形状と被検体の深度とに対応して一の 補正データを算出し、一の超音波ビームに対応する受信信号が受信処理される間に 、一の超音波ビームとは異なる他の超音波ビームのビーム形状と被検体の深度とに 対応して他の補正データを算出できる。
[0015] また、表示手段については、超音波ビームのビーム形状に対応して算出されたダイ ナミックレンジが被検体の深度の関数としてグラフ表示できる。また、超音波ビームの ビーム形状に対応して算出されたビーム強度が被検体の深度の関数としてグラフ表 示できる。
[0016] 本発明の超音波イメージング方法は、超音波探触子により形成される超音波ビー ムのビーム形状と被検体の深度とに対応した補正データを算出する工程と、超音波 探触子に駆動信号を供給して超音波探触子力 被検体に超音波を送波する工程と 、送波工程に対応して超音波探触子から出力される受信信号を処理する受信工程と 、受信工程力 出力される受信信号を補正データに基づき補正する工程と、補正ェ 程から出力される受信信号に基づいた超音波像を表示する工程を有することを特徴 とする。
[0017] この場合において、補正データ算出工程では、超音波ビームのビーム形状と被検 体の深度とに対応したダイナミックレンジを算出できる。この場合の補正工程では、補 正データ算出工程で算出されたダイナミックレンジに基づき、受信工程から出力され る受信信号のダイナミックレンジを補正できる。
[0018] また、補正データ算出工程では、超音波ビームのビーム形状と被検体の深度とに 対応したビーム強度を算出できる。この場合の補正工程では、補正データ算出工程 で算出されたビーム強度に基づき、受信工程力も出力される受信信号の信号強度を 補正できる。
[0019] また、補正データ算出工程では、超音波送波ビームまたは超音波受波ビームの少 なくとも一方のビーム形状に対応した補正データを算出できる。
[0020] また、受信工程では、超音波探触子力 出力される受信信号をサンプリングクロック に基づきディジタル信号に変換できる。この場合の補正工程では、サンプリングクロッ クに基づき被検体の深度方向に設定された一以上のサンプリング点ごとに、受信ェ 程力 出力される受信信号を補正できる。
[0021] また、補正データ算出工程では、駆動信号の周波数、駆動信号の波数、超音波ビ ームの焦点の座標、超音波探触子の口径を形成する振動子の数、超音波ビームを 形成するためのフォーカスデータの少なくとも 1つを含むパラメータに基づき補正デ ータを算出できる。
[0022] さらに、補正データ算出工程では、一の超音波ビームのビーム形状と被検体の深 度とに対応して一の補正データを算出する工程と、一の超音波ビームに対応する受 信信号が受信処理される間に、一の超音波ビームとは異なる他の超音波ビームのビ ーム形状と被検体の深度とに対応して他の補正データを算出する工程を含めること ができる。
[0023] また、表示工程では、超音波ビームのビーム形状に対応して算出されたダイナミツ クレンジを被検体の深度の関数としてグラフ表示できる。また、超音波ビームのビーム 形状に対応して算出されたビーム強度を被検体の深度の関数としてグラフ表示でき る。
図面の簡単な説明
[0024] [図 1]本発明を適用した第 1の実施形態の超音波診断装置の構成を示すブロック図 である。
[図 2]図 1の超音波診断装置の全体処理を示すフローチャートである。
[図 3]図 1のダイナミックレンジ補正部の動作を説明するための概念図である。
[図 4]図 1の超音波診断装置により撮像した超音波像の表示例と参考の形態の装置 により撮像した比較例である。
[図 5]超音波ビームのビーム形状が深度方向に異なることを説明する図である。 [図 6]本発明を適用した第 2の実施形態の超音波診断装置の構成を示すブロック図 である。
[図 7]図 6の超音波診断装置の全体処理を示すフローチャートである。
[図 8]図 6の信号強度補正部の動作を説明するための概念図である。
[図 9]本発明を適用した第 3の実施形態の超音波診断装置の構成を示すブロック図 である。
[図 10]図 9の超音波診断装置の全体処理を示すフローチャートである。
[図 11]図 9の補正部の動作を説明するための概念図である。
発明を実施するための最良の形態
[0025] (第 1の実施形態)
本発明を適用した超音波診断装置の第 1の実施形態について図面を参照して説 明する。図 1は、本実施形態の超音波診断装置 1の構成を示すブロック図である。
[0026] 図 1に示すように、超音波診断装置 1は、被検体との間で超音波を送受する超音波 探触子 10と、超音波探触子 10に駆動信号を供給する送信手段としての送信回路 1 2と、超音波探触子 10から出力される受信信号を処理する受信手段など力 構成さ れる。ここでの受信手段は、超音波探触子 10から出力される受信信号に対し増幅な どの処理をする受信回路 16と、受信回路 16から出力される受信信号をサンプリング クロックに基づきディジタル信号に変換するアナログディジタル変換部 18 (以下、 AZ D変換部 18という)と、 AZD変換部 18から出力される受信信号に対し整相加算処 理をするソフトウェア整相加算手段としての整相加算部 20と、整相加算部 20から出 力される受信信号に対し検波などの処理をする信号処理部 22を有する。
[0027] そして、本発明に係る超音波診断装置 1は、超音波探触子 10で送受される超音波 ビームのビーム形状と被検体の深度に相関して算出される補正データに基づき、信 号処理部 22から出力される受信信号を補正する手段としてのダイナミックレンジ補正 部 24を有する。本実施形態での超音波ビームとは、超音波送波ビーム及び超音波 受波ビームを意味するが、いずれか一方でもよい。
[0028] さらに詳細に、超音波診断装置 1の構成について説明する。超音波探触子 10は、 複数の振動子を配列して形成されている。各振動子は、送信回路 12から供給される 駆動信号を超音波に変換すると共に、被検体力 発生する反射エコーを受信信号に 変換する。ここでの駆動信号は、超音波送波用のパルス信号である。
[0029] 送信回路 12は、受信回路 16と伴に送受信部 26に設けられている。送受信部 26は 、送信回路 12から出力される駆動信号を超音波探触子 10に出力すると共に、超音 波探触子 10から出力される受信信号を受信回路 16に出力する送受分離回路と、超 音波探触子 10の送波口径又は受波口径を形成する振動子群を選択する高圧切替 えスィッチを有する。
[0030] 受信回路 16は、超音波探触子 10から出力される受信信号を増幅する受信アンプ を有する。受信アンプは、受信信号に対し所定の信号増幅度で増幅するプリアンプ と、超音波の反射部位の深度に起因する反射エコーの減衰に応じて信号増幅度が 可変するゲインコントロールアンプを有する。このような受信回路 16の処理により、受 信信号は、診断に足りるに十分な大きさに増幅されると共に、被検体内での減衰が 補正される。
[0031] 送受信部 26にフォーカス位置決定部 28が接続している。フォーカス位置決定部 2 8は、送受信部 26に超音波送波ビームの焦点の設定座標を出力する。なお、フォー カス位置決定部 28に焦点の座標を入力するキーボードやマウスなどを設けてもよい
[0032] AZD変換部 18は、サンプリングクロックを発生するクロック発生器を有する。ここで のサンプリングクロックのクロック間隔は、等間隔でもよいし、不等間隔でもよい。整相 加算部 20は、 AZD変換部 18から出力される受信信号をフォーカスデータに基づき ディジタル整相をする遅延量補正部 34と、遅延量補正部 34に付与されるフォーカス データを格納したメモリ 36と、遅延量補正部 34から出力されるチャネル毎の受信信 号を加算して信号処理部 22に出力する加算部 38を有する。このような整相加算部 2 0により、超音波受波ビームが形成される。
[0033] 信号処理部 22は、整相加算部 20から出力される受信信号に対し、検波、フィルタリ ング、対数圧縮などの処理をする機能と、処理された受信信号に基づいて超音波像 (例えば、 Bモード像、ドプラ像、 Mモード像)の基になる画像信号を構成する機能を 有する。信号処理部 22から出力される画像信号は、ダイナミックレンジ補正部 24に 出力される。便宜上、信号処理部 22から出力される画像信号も受信信号と適宜称す る。
[0034] ダイナミックレンジ補正部 24は、ダイナミックレンジ計算部 40と、メモリ 42を有する。
まず、ダイナミックレンジ計算部 40は、超音波探触子 10で送受される超音波ビーム のビーム形状と被検体の深度に相関した補正データとしてダイナミックレンジ (以下、 補正用ダイナミックレンジと称する)を算出する機能を有する。例えば、駆動信号の周 波数、駆動信号の波数、超音波送波ビーム又は超音波受波ビームの焦点の座標、 超音波探触子 10の口径の大きさ、超音波送波ビーム又は受波ビームを形成するた めのフォーカスデータの少なくとも 1つを含むパラメータに基づき、補正用ダイナミック レンジを計算する。ここでの口径の大きさとは、送波口径又は受波口径を形成する振 動子の数と、各振動子の大きさに基づき求められる。そしてダイナミックレンジ計算部 40は、補正用ダイナミックレンジをダイナミックレンジテーブルに配列する。ダイナミツ クレンジテーブルは、被検体の深度 (例えば、サンプリング点又はフォーカス段の深 度)に対応させて補正用ダイナミックレンジが配列されたデータベースである。換言す れば、ダイナミックレンジテーブルは、フォーカス深度の関数とした超音波ビームプロ ファイルを有する。このようなダイナミックレンジテーブルは、メモリ 42に格納される。
[0035] また、ダイナミックレンジ計算部 40は、補正用ダイナミックレンジに基づき、信号処 理部 22から出力される受信信号のダイナミックレンジを補正する。例えば、 AZD変 換部 18のサンプリングクロックに基づ 、て被検体の深度方向に設定されたサンプリン グ点ごとに、信号処理部 22から出力される受信信号を補正する。ここで、サンプリン グ点ごとに受信信号を補正することに代えて、複数のサンプリング点を有するフォー カス段ごとに受信信号を補正してもよい。
[0036] このようなダイナミックレンジ補正部 24は、第 1の超音波ビームのビーム形状と被検 体の深度に基づいて第 1の補正用ダイナミックレンジを算出し、第 1の超音波ビーム に対応する受信信号が受信処理される間に、第 1の超音波ビームとは異なる第 2の 超音波ビームのビーム形状と被検体の深度に基づいて第 2の補正用ダイナミックレン ジを算出できる。すなわち、ダイナミックレンジ計算部 40は、時分割処理で補正用ダ イナミックレンジを事前に算出できる。又は、予め設定してある超音波ビームの形状 に対し、それに対応する補正用ダイナミックレンジを事前にダイナミックレンジテープ ル 42へ格納しておくこともできる。
[0037] また、ダイナミックレンジ補正部 24から出力される補正後の受信信号を表示用の信 号に変換すると共に、受信信号に輝度合わせ処理をするディジタルスキャンコンパ一 タ 30 (以下、 DSC30という)と、 DSC30から出力された受信信号に基づいた超音波 像を表示画面に表示する表示部 32が設けられて 、る。
[0038] このように構成される超音波診断装置 1の基本動作を説明する。まず、被検体に超 音波探触子 10を接触させる。超音波探触子 10に駆動信号が送信回路 12から供給 されると、超音波探触子 10から超音波送波ビームが被検体に照射される。被検体か ら発生する反射エコーは、超音波探触子 10により受波されることによって受信信号 に変換される。受信信号は、受信回路 16により増幅される。増幅された受信信号は、 AZD変換部 18によりディジタル信号に変換される。ディジタルィ匕された受信信号は 、整相加算部 20により超音波受波ビームとして形成される。整相加算部 20から出力 される超音波受波ビームは、信号処理部 22により所定の処理が施されることによって 画像信号として出力される。信号処理部 22から出力された画像信号は、ダイナミック レンジ補正部 24によりダイナミックレンジが補正される。補正された画像信号は、 DS C30を介して表示部 32に超音波像として表示される。
[0039] ここで、超音波診断装置 1の詳細動作についてダイナミックレンジ補正部 24を中心 に図 2及び図 3を参照して説明する。図 2は、図 1の超音波診断装置 1の全体処理を 示すフローチャートである。図 3は、図 1のダイナミックレンジ補正部の動作を説明す るための概念図である。なお、図 3は、便宜上、 1つのフォーカス段内における処理を 示すが、他のフォーカス段での処理も同様である。
[0040] く補正用ダイナミックレンジの算出工程 (S100) >
補正用ダイナミックレンジ D1— D6がダイナミックレンジ計算部 40により算出される 。ここで補正用ダイナミックレンジ D1— D6は、超音波探触子 10により形成される超 音波ビームのビーム形状と各反射部位 46— 1一 46— 7の深度に相関した補正データ である。なお、図 3Aに示すように、反射部位 46— 1と反射部位 46— 6は、同じ深度で あるため、双方が補正用ダイナミックレンジ D1に対応する。 [0041] より具体的には、ダイナミックレンジ計算部 40は、超音波ビームのビーム形状に関 連する因子力も補正用ダイナミックレンジ D1— D6を算出する。因子とは、送信回路 12から供給される駆動信号の周波数、その駆動信号の波数、超音波ビームの焦点 の座標、超音波探触子 10の口径を形成する振動子の数、整相方式などのパラメータ である。整相方式は、駆動信号を遅延するための送信用フォーカスデータ、受信信 号を整相するための受信用フォーカスデータ (例えば、メモリ 36に記憶されているフ オーカスデータ)である。
[0042] すなわち、補正用ダイナミックレンジ D1— D6を設定するための一般式は、数 1式の ように表される。これにより、例えば、送受波フォーカスが表 1に示す場合は、ダイナミ ックレンジは同表のように設定される。
[0043] (数 1式) 補正用ダイナミックレンジ DR=f (口径、超音波送受フォーカス深度、超 音波周波数、整相方式)
[0044] [表 1]
Figure imgf000011_0001
要するに、補正用ダイナミックレンジ Dl— D6は、超音波ビームが有する固有のダ イナミックレンジであり、各反射部位 46—1— 46— 7におけるビーム強度に依存したも のである。このような補正用ダイナミックレンジ D1— D6は、サンプリング点又はフォー カス段の深度に対応付けられた後、ダイナミックレンジテーブルとしてメモリ 42に格納 される。図 3Cは、受信信号 S1— S6に対応した補正用ダイナミックレンジ D1— D6を 示している。なお、補正用ダイナミックレンジ D1— D6については、実測ゃシミュレ一 シヨンなどから求めてもよ 、。
[0045] <超音波の送波工程(S101) > 超音波探触子 10に駆動信号を供給して超音波探触子 10から被検体に超音波を 送波する。例えば、被検体の体表に接触させた超音波探触子 10に対し、送波回路 1 2から駆動信号が供給される。なお、駆動信号は、超音波探触子 10に有する振動子 の数に対応して複数生成され、生成された各駆動信号は、焦点位置に基づいた送 信用フォーカスデータを用いてそれぞれ遅延されている。また、フォーカス位置は、フ オーカス位置決定部 28から送受信部 26に所定タイミングで出力される。
[0046] 一方、超音波探触子 10に有する複数の振動子群のうち、所定の振動子群が高圧 切替えスィッチにより選択されることにより、送波口径が形成される。送波口径の振動 子群に駆動信号が供給されると、供給された各振動子から超音波が送波される。送 波された各超音波により超音波送波ビームが形成される。形成された超音波送波ビ ームは、送波方向の所定距離に焦点を有する。
[0047] より具体的に図 3Aを参照して説明する。図 3Aの縦軸は被検体の深度を示し、横 軸はビームの走査アドレスを示している。図 3Aに示すように、複数の反射部位 46— 1 一 46— 7が被検体に散乱している。各反射部位 46-1— 46— 7は、大きさや音響イン ピーダンスが異なる。また、反射部位 46— 1— 46— 7は、反射部位 46— 6を除き、それ ぞれ深度が異なる。反射部位 46- 6は、反射部位 46- 1と同じ深度である。
[0048] 本工程が開始されると、超音波探触子 10のビームアドレスが図 3Aに示すアドレス Aに初期化される。そして、超音波探触子 10から超音波送波ビームが射出される。 すなわち、超音波送波ビームは、アドレス A力 走査方向にアドレスを変えながら順 番に反射部位 46— 1一 46— 7に射出される。なお、図 3Aでは深度方向と走査方向の 二次元分布を示すが、三次元分布の場合も同様である。また、超音波送波ビームの 焦点位置につ!、ては適宜変更できる。
[0049] <超音波の受信工程(S102) >
超音波送波工程 (S101)に対応して超音波探触子 10から出力される受信信号を 処理する。例えば、反射部位 46-1— 46- 7から発生した反射エコーは、超音波探触 子 10により受波されて受信信号に変換される。変換された受信信号は、受信回路 16 、 AZD変換部 18、整相加算部 20により順に処理されることにより、超音波受波ビー ムが形成される。形成された超音波受波ビームは、被検体の所定深度に焦点を有す る。そして超音波受波ビームは、信号処理部 22により処理された後、図 3Bに示すよ うに、受信信号 S1— S6としてダイナミックレンジ補正部 24に出力される。なお、図 3B の縦軸は被検体の深度を示し、横軸は信号強度を示す。受信信号 S1— S6の信号 強度は、画像信号を形成する際の例えば輝度情報に相当する。
[0050] また、説明の便宜上、図 3Bでは、受信信号 S1— S6の信号強度を深度ごとに射影 されたものとして表している。また、反射部位 46— 1と反射部位 46— 6は同じ深度であ るが、反射部位 46-1に対応する受信信号の信号強度の方が大きいため、反射部位 46-1に対応する受信信号を受信信号 S 1として 、る。
[0051] <受信信号の補正工程 (S103) >
超音波の受信工程 (S102)から出力される受信信号 SI— S6は、補正用ダイナミツ クレンジに基づいて補正される。例えば、受信信号 S1は、ダイナミックレンジ補正部 2 4に入力される。受信信号 S1の入力に応じ、ダイナミックレンジ計算部 40により補正 用ダイナミックレンジ D1がメモリ 42から読み出される。補正用ダイナミックレンジ D1は 、ダイナミックレンジテーブルに配列されたものであり、反射部位 46— 1の深度に対応 してダイナミックレンジテーブル力も抽出される。そして、受信信号 S1は、補正用ダイ ナミックレンジ D1に基づき、ダイナミックレンジ計算部 40により信号強度が補正される 。反射部位 46—2— 46— 7についても同様に処理される。このように処理された補正 後の受信信号 DN1— DN6を図 3Dに示す。要するに、本工程により、受信信号 DN 1一 DN6は、超音波ビームが有する固有のダイナミックレンジが考慮された信号にな る。
[0052] なお、本工程の補正計算は、例えば数 2式に基づき行われる。数 2式の kは、補正 対象の受信信号の数を示す。本形態の kの値は 1一 6である力 適宜変更すればよ い。 aは、反射エコーが被検体内を伝播する際に減衰した信号強度を補正する係数 を示す。
[0053] (数 2式)
D Nk= k O lo l 0 (¾) /Dk]
<超音波像の表示工程 (SI 04) > 補正工程 (S103)から出力される受信信号 DN1— DN6に基づいた超音波像が表 示画面に表示される。例えば、受信信号 DN1— DN6は、 DSC30により表示用の信 号に変換されると共に、輝度合わせ処理が施される。 DSC30から出力される受信信 号 DN1— DN6は、超音波像として表示部 32の表示画面に表示される。また、補正 用ダイナミックレンジ D1— D6は、反射部位 46—1— 46— 7の深度の関数として超音 波像と並べてグラフ表示される。グラフ表示により、補正用ダイナミックレンジ D1— D 6が視覚的に把握される。
[0054] 上述したように、超音波送波ビーム又は超音波受波ビームは、口径力 焦点に至る まで絞られ、焦点を超えた後は広がる形状を有する。又、伝搬距離に応じた減衰が 生じるから、各反射部位 46— 1一 46— 7の深度によって音圧又は反射強度に差が生 じる。例えば、超音波ビームは、被検体の深度方向にビーム断面積が異なるため、ビ ーム強度が被検体の深度方向に異なる。これに起因して、信号処理部 22から出力さ れる受信信号は、超音波ビームの形状と反射部位 46— 1一 46— 7の深度の双方によ る影響を受けたものになる。
[0055] したがって、信号処理部 22から出力される受信信号 S1— S6に対し、超音波ビーム の形状と反射部位 46— 1一 46— 7の深度に相関した補正を施す。これにより、補正後 の受信信号 DN1— DN6は、超音波ビームの形状と反射部位 46— 1一 46— 7の深度 による影響が低減された信号になる。このような受信信号 DN1— DN6に基づく超音 波像を表示することにより、超音波像の画質を改善できる。
[0056] 図 4Aは、本実施形態の超音波診断装置によりファントム撮像した超音波像の表示 例である。図 4Bは、参考の形態における超音波診断装置によりファントム撮像した超 音波像の表示例である。
[0057] 参考の形態の超音波診断装置では、超音波像を表示する際のダイナミックレンジ 力 反射部位の深度にかかわらず一定 (例えば、 lOOdB)に設定される。一般に、深 度の浅い反射部位力 発生した反射エコーは信号強度が大きくなる一方、深度の深 い反射部位力 発生した反射エコーは、被検体内での減衰に起因して信号強度が 小さくなる。したがって、浅い反射部位に対応する受信信号のダイナミックレンジにあ わせて表示すると、深い反射部位に対応する画像のコントラストが低下する場合があ る。また、深い反射部位に対応する受信信号のダイナミックレンジにあわせて表示す ると、浅!ヽ反射部位に対応する画像が輝度飽和して情報が欠落する場合がある。
[0058] 図 4Bの場合、反射部位の各深度で最もダイナミックレンジが広い超音波ビームを 基準としている。換言すれば、全ての受信信号を表現できるダイナミックレンジが設 定されているため、受信信号が本来有しているダイナミックレンジよりも広く表示され ている。その結果、図 4Bに示すように、不要な成分がノイズとして超音波像に表示さ れること〖こなる。
[0059] この点、本実施形態では、超音波ビームのビーム形状に相関させて受信信号のダ イナミックレンジが設定される。すなわち、反射部位 46— 1一 47— 7の各深度における 超音波ビームが有するダイナミックレンジに忠実な超音波像が表示される。したがつ て、図 4Aに示す超音波像は、図 4Bと比較すると、診断情報としての信号の欠落がな いことにカ卩え、ノイズとなる不要な成分が低減されている。このような効果は、特に、超 音波ビームが有する固有のダイナミックレンジが狭くなる部位ほど顕著になる。例え ば、可変口径を採用してビーム形成を行った場合の浅部などで顕著である。
[0060] すなわち、本実施形態によれば、補正後の受信信号 DN1— DN6は、超音波ビー ムが有する固有のダイナミックレンジ D1— D6に応じた関数で重み付けされたものと なる。このような補正後の受信信号 DN1— DN6に基づき、超音波像を表示する際の ダイナミックレンジが深度ごとに設定される。したがって、表示部 32に表示される超音 波像は、ノイズ成分が低減されると共に、反射部位 46-1— 46— 7の音響インピーダ ンスの違いをより忠実に現したものになる。
[0061] 要するに、本実施形態は、超音波ビーム形状が被検体の深度方向に変化すること に着目してなされたものである。すなわち、超音波ビーム形状は、超音波送波ビーム の焦点位置、送波口径の大きさ、受波口径の大きさ、超音波受波ビームの整相方式 などにより被検体の深度方向に変化する。したがって、超音波ビーム強度 (音場の大 きさ)が深度方向によって異なる。超音波ビーム形状が深度方向に異なると、それに 伴って超音波ビームが有するダイナミックレンジが変わる。この点、本実施形態によ れば、被検体の深度方向の画像信号のダイナミックレンジを超音波ビームに相関さ せて適切に設定することにより、超音波像の画質を向上させることができる。 [0062] また、受信信号のダイナミックレンジが例えば 70dBを超えるなど比較的広 ヽ幅を有 する場合がある。この点、本実施形態によれば、超音波像として表示される際のダイ ナミックレンジは、超音波ビーム形状に相関させて反射部位の深度ごとに設定される 。これにより、受信信号のダイナミックレンジが比較的広い幅を有するときでも、微小 なノイズ信号が画像に顕在化したり、画素間の輝度差が小さくなつたり、反射部位に 関する診断情報が欠落したりすることを抑制できる。
[0063] さらに、超音波ビームの焦点位置や整相方式などの変更に伴って超音波ビーム形 状が随時変わるのが一般的である。この点、本実施形態によれば、超音波ビーム形 状の変化に追従して、ダイナミックレンジ計算部 40によりダイナミックレンジテーブル が自動更新される。したがって、装置の使い勝手や操作性が向上する。
[0064] 第 1の実施形態により本発明を説明したが、これに限られるものではない。例えば、 図 2に示す補正処理については、ビームアドレス単位で実行してもよいし、フレーム 単位で実行してもよい。フレーム単位で実行する場合、図 3Aを参照して説明したよう に、同一の深度に複数の反射部位 (例えば、反射部位 46—1、 46— 6)があるときは、 受信信号の強度が大きい方を基準にすればよい。また、 Bモード像 (断層像)の撮像 に限らず、ドプラ像や Mモード像を撮像するときにも本発明を適用できる。
[0065] また、ダイナミックレンジ補正部 24に入力される受信信号は、信号処理部 22の出力 信号であるが、これに代えて、信号処理部 22の途中や、整相加算部 20の出力信号 としてもよいし、 AZD変換部 18の出力信号としてもよい。要するに、受信手段から出 力される受信信号をダイナミックレンジ補正部 24に出力すればよい。
[0066] なお、超音波診断装置 1の送受信系回路のダイナミックレンジは、超音波探触子 1 0により送受される超音波ビームが有する固有のダイナミックレンジよりも広いことが望 ましい。
[0067] <超音波ビーム形状に関する補足 >
超音波ビームのビーム形状が被検体の深度方向に異なる一例について補足説明 する。図 5Aは、被検体の深度方向に設定された複数のサンプリング点 T1一 TP (P : 2以上の自然数)の概念図である。図 5Bは、各フォーカス段 A— Dで形成された受波 ビームのビーム形状を示す図である。図 5Cは、超音波送波ビームのビーム形状を示 す図である。
[0068] 図 5Aに示すように、被検体の深度方向に複数のサンプリング点 T1一 TPが設定さ れている。サンプリング点 T1一 TPの各間隔は、 AZD変換部 18のサンプリングクロッ ク間隔に対応して等しく設定されている。複数のサンプリング点 T1一 TPは、例えば 4 つごとにフォーカス段 A— Dに分けられる。例えば、サンプリング点 T1一 T4は、フォ 一カス段 Aに属する。なお、 1つのフォーカス段に属するサンプリング点の数につい ては、診断に支障を来たさない範囲内で適宜決められる。
[0069] このようなフォーカス段 A— Dごとに共通のフォーカスデータが設定される。設定さ れたフォーカスデータは、メモリ 36に格納される。受信用フォーカスデータとは、超音 波探触子 10の各振動子から出力される各受信信号を整相する遅延量や微小遅延 量であり、フォーカス段 A— Dごとに切替えられる。メモリ 36に格納されたフォーカス データは、制御指令に応じて読み出される。読み出されたフォーカスデータに基づき 、AZD変換部 18から出力される受信信号が整相加算部 20により遅延され整相され る。これによつて、各サンプリング点 T1一 TPに対応した超音波受波ビームがそれぞ れ形成される。要するに、いわゆるダイナミックフォーカスが行われる。ダイナミックフ オーカスとは、深度が浅 、部位力 深 、部位まで焦点を 、くつかもつことで深度方向 に比較的広い範囲でビームを集束する手法である。
[0070] ここでのフォーカスデータは、診断に影響の無い範囲でフォーカス段 A— Dごとに 独立のフォーカスデータが設定されている。これにより、サンプリング点 T1一 TPごと に設定するときよりも、フォーカスデータの総数を減らすことができ、メモリ 36の記憶 領域を有効に利用できる。
[0071] このようにフォーカス段 A— Dごとにフォーカスデータが異なることに起因して、超音 波受波ビームのビーム形状 (ビームプロファイル)は、被検体の深度方向に異なる。 例えば、図 5Bは、上段から順に、フォーカス段 Aに属するサンプリング点 T1に対応し た受波ビーム形状、フォーカス段 Bに属するサンプリング点 T5に対応した受波ビーム 形状、フォーカス段 Cに属するサンプリング点 T(P— 8)に対応し受波ビーム形状、フ オーカス段 Dに属するサンプリング点 Τ (Ρ— 4)に対応した受波ビーム形状を示して!/ヽ る。図 5Βから分かるように、超音波受波ビームは、ビーム形状が被検体の深度方向 に異なる。
[0072] なお、同一のフォーカス段 (例えば、フォーカス段 A)に属する複数のサンプリング 点(例えば、サンプリング点 T1一 T4)間においても、深度が若干ことなるため、各サ ンプリング点における超音波受波ビームのビーム形状も若干異なる。
[0073] また、超音波受波ビームのビーム形状は、可変口径が行われることでも深度方向に 異なる。可変口径とは、口径幅を自動的に浅い部分になるほど小さくする方法である 。可変口径により、超音波探触子 10に近いサンプリング点(例えば、サンプリング点 T 1、 T2)力 発生する反射エコーの広がりを抑制できる。このような可変口径を行う場 合、反射エコーを受波する口径の大きさ、つまり受波口径を形成する振動子の素子 数が変わる。したがって、超音波受波ビームのビーム形状が深度方向に異なる。
[0074] さらに、超音波送波ビームについても、ビーム形状が深度方向に異なる。例えば、 図 5Cに示すように、超音波送波ビームは、焦点を一点に絞ることにより形成される。 換言すれば、超音波送波ビームは、超音波探触子 10から焦点 41に至るまで絞られ 、焦点を超えた後は広がる形状となる。なお、フレームレートが要求されないときは、 フォーカス段 A— Dごとに焦点を設定した超音波送波ビームを独立に形成し、複数 回送波してもよい。
[0075] このように超音波送波ビームまたは超音波受波ビームのビーム形状は、被検体の 深度方向に異なったものになる。したがって、超音波ビーム強度が深度方向に異な ることに起因して、超音波ビームが有する固有のダイナミックレンジも深度方向に異な る。例えば、超音波ビーム強度が大きい深度では超音波ビームが有する固有ダイナ ミックレンジが広くなる。超音波ビーム強度が小さ!/、深度では固有ダイナミックレンジ が狭くなる。このような超音波ビームの形状と被検体の深度による影響を低減するた めに、本実施形態は、超音波ビームのビーム形状と被検体の深度に相関して算出さ れたダイナミックレンジに基づき、信号処理部 22から出力される受信信号のダイナミ ックレンジを補正する。
[0076] (第 2の実施形態)
本発明を適用した超音波診断装置の第 2の実施形態について図 6ないし図 8を参 照して説明する。本実施形態は、超音波ビームのビーム形状と被検体の深度に相関 して算出されるビーム強度に基づき、信号処理部 22から出力される受信信号を補正 する点で第 1の実施形態と異なる。したがって、第 1の実施形態と相互に対応する箇 所については同一符号を付し、相違点を中心に説明する。
[0077] 図 6は、本実施形態の超音波診断装置 2の構成を示すブロック図である。図 2に示 すように、超音波診断装置 2は、第 1の実施形態 (例えば図 1)のダイナミックレンジ補 正部 24に代えて、信号強度補正部 50を備えている。信号強度補正部 50は、超音波 探触子 10により送受される超音波ビームのビーム形状と被検体の深度に相関して算 出される超音波ビームに基づき、信号処理部 22から出力される受信信号の信号強 度を補正する。例えば、信号強度補正部 50は、信号強度計算部 52と、メモリ 54を有 する。信号強度計算部 52は、超音波ビームのビーム形状と被検体の深度に相関し た補正データとしてのビーム強度 (以下、補正用ビーム強度と適宜称する)を算出す ると共に、補正用ビーム強度に基づき信号処理部 22から出力される受信信号の信 号強度を補正する。メモリ 54は、被検体の深度 (例えば、サンプリング点又はフォー カス段の深度)の関数として補正用ビーム強度が配列されたビーム強度テーブルを 格納する。
[0078] 図 7は、図 6の超音波診断装置 2の全体処理を示すフローチャートである。図 8は、 図 6の信号強度補正部 50の動作を説明するための概念図である。図 7に示すように 、本実施形態は、第 1の実施形態の補正ダイナミックレンジの算出工程 S100に代え て、補正用ビーム強度の算出工程 S200を備えると共に、受信信号の補正工程 S10 3に代えて、受信信号の補正工程 S203を備えている。したがって、補正用ビーム強 度の算出工程 S200と受信信号の補正工程 S203を中心に説明する。
[0079] <補正用ビーム強度の算出工程 (S200) >
補正用ビーム強度 A1— A4が信号強度計算部 52により算出される。補正用ビーム 強度 A1— A4とは、超音波探触子 10により形成される超音波ビームのビーム形状と 各反射部位 46-1— 46— 7の深度に相関した補正データである。例えば、図 8Cに示 すように、ビーム強度 A1は、超音波ビームの焦点位置のビーム強度に対応する。ビ ーム強度 A2は、反射部位 46-4、 46- 5に対応する。ビーム強度 A3は、反射部位 4 6—2、 46— 7に対応する。ビーム強度 A4は、反射部位 46— 1、 46— 3、 46— 6に対応 する。
[0080] 例えば、信号強度計算部 52は、超音波ビームのビーム形状に関連する因子力 補 正用ビーム強度 A1— A4を算出する。ここでの因子とは、送信回路 12から供給され る駆動信号の周波数、その駆動信号の波数、超音波ビームの焦点の座標、超音波 探触子 10の口径を形成する振動子の数、整相方式などのパラメータである。整相方 式は、駆動信号を遅延するための送信用フォーカスデータ、受信信号を整相するた めの受信用フォーカスデータ(例えば、メモリ 36に記憶されているフォーカスデータ) により決まってくる。要するに、補正用ビーム強度 A1— A4は、超音波ビームが固有 に有する信号強度である。このような補正用ビーム強度 A1— A4は、サンプリング点 又はフォーカス段の深度に対応付けられた後、ビーム強度テーブルとしてメモリ 54に 格納される。
[0081] 本工程の計算処理については、例えば、第 1の超音波ビームのビーム形状と被検 体の深度に基づいて補正用ビーム強度を算出し、第 1の超音波ビームに対応する受 信信号が受信処理される間に、第 1の超音波ビームとは異なる第 2の超音波ビーム のビーム形状と被検体の深度に基づいて補正用ビーム強度を事前に算出すればよ い。また、実測やシミュレーションなど力 補正用ビーム強度 A1— A4を求めてもよい
[0082] <受信信号の補正工程 (S203) >
超音波の受信工程 (S102)から出力される受信信号 S1— S6は、補正用ビーム強 度 A1— A4に基づいて補正される。例えば、受信信号 S1— S6が信号強度補正部 5 0に入力すると、信号強度計算部 52により補正用ビーム強度 A1— A4カ モリ 54か ら読み出される。そして、補正用ビーム強度 A1— A4のうち、補正用ビーム強度 A1 が基準強度として設定される。なお、基準とする補正用ビーム強度については適宜 変更してよい。
[0083] 次に、例えば、受信信号 S1の信号強度を補正する場合、基準強度としての補正用 ビーム強度 A1と、補正用ビーム強度 A4との相対比(例えば、 A1ZA4)が求められ る。求められた相対比は、超音波ビームの焦点位置におけるビーム振幅と、反射部 位 46— 1の深度におけるビーム振幅との振幅比に対応する。このような相対比が補正 係数として設定される。
[0084] そして、図 8Dに示すように、受信信号 S1の信号強度に補正係数 (A1ZA4)を乗 算することにより、受信信号 S1の信号強度が補正される。同様な処理が受信信号 S2 一 S6に対して行われる。要するに、受信信号 S1— S6は、本工程により、受信信号 B 1一 B6として補正される。補正後の受信信号 B1— B6は、超音波ビームが有する固 有のビーム強度が考慮された信号となる。
[0085] 一般に、超音波像は、被検体の深度方向に散乱する反射部位 (例えば、反射部位
46-1-46-7)の音響インピーダンスの違いを画像化したものである。したがって、 超音波ビームが深度方向に一様なビーム強度 (音場)であるとすれば、各反射部位 の音響インピーダンスの違いがそのまま超音波像に表示されることになる。しかし、第
1の実施形態でも説明したように、超音波ビームのビーム強度はビーム形状に相関し て深度方向に異なる。その結果、表示される超音波像は、送波ビームの深度方向の 強度分布による誤差を含んだものとなる。
[0086] この点、本実施形態によれば、補正後の受信信号 B1— B6は、超音波ビームの深 度方向の強度分布に基づいて信号強度が補正された信号になる。すなわち、補正 後の受信信号 B1— B6は、超音波ビームの深度方向のビーム強度分布が一様であ るとした場合の受信信号と等価になる。
[0087] したがって、補正後の受信信号 B1— B6に基づき超音波像を表示することにより、 反射部位 46— 1一 46— 7の音響インピーダンスの違いをより忠実に超音波像に現すこ とができる。例えば、本来明るく表示されるはずの反射部位が若干暗く表示されたり、 本来暗く表示されるはずの反射部位が若干明るく表示されたりすることがないなど、 反射部位が本来有する輝度情報をより忠実に表示できる。
[0088] (第 3の実施形態)
本発明を適用した超音波診断装置の第 3の実施形態について図 9ないし図 11を参 照して説明する。本実施形態は、信号処理部 22から出力される受信信号に対し、第 一段階補正として第 2の実施形態の補正をした後、第二段階補正として第 1の実施 形態の補正をする点で、第 1及び第 2の実施形態と異なる。したがって、第 1及び第 2 の実施形態と相互に対応する箇所については同一符号を付し、相違点を中心に説 明する。
[0089] 図 9は、本実施形態の超音波診断装置 3の構成を示すブロック図である。図 9に示 すように、超音波診断装置 3は、第 1の実施形態 (例えば図 1)のダイナミックレンジ計 算部 40と、第 2の実施形態 (例えば図 2)の信号強度計算部 52と、メモリ 58とを備え た補正部 56を有する。ここでのメモリ 58は、補正用ダイナミックレンジテーブルと補正 用ビーム強度テーブルを格納して 、る。
[0090] このような補正部 56は、超音波探触子 10で送受される超音波ビームのビーム形状 と被検体の深度に相関して補正データ (例えば、補正用ビーム強度 A1— A4、補正 用ダイナミックレンジ D1— D6)を算出し、算出した補正データに基づき、信号処理部 22から出力される受信信号を補正して DSC30に出力する。
[0091] 図 10は、図 9の超音波診断装置 3の全体処理を示すフローチャートである。図 11 は、図 9の補正部 56の動作を説明するための概念図である。図 10に示すように、本 実施形態の補正データの算出工程は、第 1の実施形態 (例えば図 2)の補正ダイナミ ックレンジの算出工程 S100と、第 2の実施形態 (例えば図 7)の補正ビーム強度の算 出工程 S 200を備えている。例えば、補正データの算出工程は、補正ダイナミックレ ンジ D1— D6を算出してダイナミックレンジテーブルとしてメモリ 58に格納すると共に 、補正ビーム強度 A1— A4を算出して信号強度テーブルとしてメモリ 58に格納する。
[0092] また、受信信号を補正する工程は、第 2の実施形態の受信信号の補正工程 S203 と、第 1の実施形態の受信信号の補正工程 S103を備えている。例えば、図 11に示 すように、信号処理部 22から出力された受信信号 S1— S6は、補正ビーム強度 A1 一 A4を用いた第一段階の補正が行われることにより、受信信号 B1— B6になる。そ の後、受信信号 B1— B6は、補正ダイナミックレンジ D1— D6を用いた第二段階の補 正が行われることにより、第二段階補正後の受信信号 DN1— DN6になる。
[0093] すなわち、信号処理部 22から出力される受信信号 S1— S6に対し、超音波ビーム の形状と被検体の深度に相関した第一段階及び第二段階の補正を施すことにより、 補正後の受信信号 DN1— DN6は、超音波ビームの形状と被検体の深度による影響 がー層低減された信号になる。このような受信信号 DN1— DN6に基づく超音波像を 表示することにより、超音波像の画質をより一層改善できる。 [0094] 以上、第 1ないし第 3の実施形態により本発明を説明したが、要するに、信号処理 部 22から出力される受信信号を補正する手段として、ダイナミックレンジ補正部 24、 信号強度補正部 50、補正部 56の少なくとも 1つを備えればよ!、。
[0095] 本発明は、その精神または主要な特徴力も逸脱することなぐ他のいろいろな形で 実施できる。そのため、前述の実施形態はあらゆる点で単なる例示にすぎず、限定 的に解釈されるものではない。本発明の範囲は、均等範囲に属する変形や変更を含 むものとする。

Claims

請求の範囲
[1] 被検体との間で超音波を送受する超音波探触子と、前記超音波探触子に駆動信 号を供給する送信手段と、前記超音波探触子から出力される受信信号を処理する受 信手段と、前記受信手段から出力される前記受信信号に基づいた超音波像が表示 される表示手段を備えた超音波診断装置において、
前記超音波探触子で送受される超音波ビームのビーム形状と前記被検体の深度と に対応して、前記受信手段から出力される前記受信信号を補正する手段を有してな ることを特徴とする超音波診断装置。
[2] 請求項 1に記載の超音波診断装置において、前記補正手段は、前記超音波ビー ムのビーム形状と前記被検体の深度に対応して算出されたダイナミックレンジに基づ き、前記受信手段から出力される前記受信信号のダイナミックレンジを補正するダイ ナミックレンジ補正手段を備えたことを特徴とする超音波診断装置。
[3] 請求項 1に記載の超音波診断装置において、前記補正手段は、前記超音波ビー ムのビーム形状と前記被検体の深度に対応して算出されたビーム強度に基づき、前 記受信手段から出力される前記受信信号の信号強度を補正する信号強度補正手段 を備えたことを特徴とする超音波診断装置。
[4] 請求項 1に記載の超音波診断装置にお!、て、前記補正手段は、超音波送波ビー ムまたは超音波受波ビームの少なくとも一方のビーム形状に対応した補正データを 算出し、該補正データに基づいて前記補正を行うことを特徴とする超音波診断装置
[5] 請求項 1に記載の超音波診断装置において、前記受信手段は、前記超音波探触 子から出力される受信信号をサンプリングクロックに基づきディジタル信号に変換す る手段を有し、
前記補正手段は、前記サンプリングクロックに基づき前記被検体の深度方向に設 定された一以上のサンプリング点ごとに、前記受信手段から出力される前記受信信 号を補正することを特徴とする超音波診断装置。
[6] 請求項 1に記載の超音波診断装置にお!、て、前記補正手段は、前記補正データを 算出する計算手段を有し、 前記計算手段は、前記駆動信号の周波数、前記駆動信号の波数、前記超音波ビ ームの焦点の座標、前記超音波探触子の口径を形成する振動子の数、前記超音波 ビームを形成するためのフォーカスデータの少なくとも 1つを含むパラメータに基づき 前記補正データを算出することを特徴とする超音波診断装置。
[7] 請求項 1に記載の超音波診断装置において、前記補正手段は、前記補正データを 算出する計算手段を有し、
前記計算手段は、一の前記超音波ビームのビーム形状と前記被検体の深度とに対 応して一の前記補正データを算出し、前記一の超音波ビームに対応する受信信号 が受信処理される間に、前記一の超音波ビームとは異なる他の前記超音波ビームの ビーム形状と前記被検体の深度とに対応して他の前記補正データを算出することを 特徴とする超音波診断装置。
[8] 請求項 1に記載の超音波診断装置にお!、て、前記表示手段は、前記超音波ビー ムのビーム形状に対応して算出されたダイナミックレンジが前記被検体の深度の関 数としてグラフ表示されることを特徴とする超音波診断装置。
[9] 請求項 1に記載の超音波診断装置にお!、て、前記表示手段は、前記超音波ビー ムのビーム形状に対応して算出されたビーム強度が前記被検体の深度の関数として グラフ表示されることを特徴とする超音波診断装置。
[10] 超音波探触子により形成される超音波ビームのビーム形状と被検体の深度とに対 応した補正データを算出する工程と、
前記超音波探触子に駆動信号を供給して前記超音波探触子から前記被検体に超 音波を送波する工程と、
前記送波工程に対応して前記超音波探触子から出力される受信信号を処理する 受信工程と、
前記受信工程から出力される受信信号を前記補正データに基づき補正する工程と 前記補正工程から出力される前記受信信号に基づいた超音波像を表示する工程 を有することを特徴とする超音波イメージング方法。
[11] 請求項 10に記載の超音波イメージング方法において、前記補正データ算出工程 は、前記超音波ビームのビーム形状と前記被検体の深度とに対応したダイナミックレ ンジを算出し、前記補正工程は、前記補正データ算出工程で算出されたダイナミック レンジに基づき、前記受信工程から出力される前記受信信号のダイナミックレンジを 補正することを特徴とする超音波イメージング方法。
[12] 請求項 10に記載の超音波イメージング方法において、前記補正データ算出工程 は、前記超音波ビームのビーム形状と前記被検体の深度とに対応したビーム強度を 算出し、前記補正工程は、前記補正データ算出工程で算出されたビーム強度に基 づき、前記受信工程から出力される前記受信信号の信号強度を補正することを特徴 とする超音波イメージング方法。
[13] 請求項 10に記載の超音波イメージング方法において、前記補正データ算出工程 は、超音波送波ビームまたは超音波受波ビームの少なくとも一方のビーム形状に対 応した補正データを算出することを特徴とする超音波イメージング方法。
[14] 請求項 10に記載の超音波イメージング方法にぉ 、て、前記受信工程は、前記超 音波探触子から出力される受信信号をサンプリングクロックに基づきディジタル信号 に変換する工程を有し、
前記補正工程は、前記サンプリングクロックに基づき前記被検体の深度方向に設 定された一以上のサンプリング点ごとに、前記受信工程から出力される前記受信信 号を補正することを特徴とする超音波イメージング方法。
[15] 請求項 10に記載の超音波イメージング方法において、前記補正データ算出工程 は、前記駆動信号の周波数、前記駆動信号の波数、前記超音波ビームの焦点の座 標、前記超音波探触子の口径を形成する振動子の数、前記超音波ビームを形成す るためのフォーカスデータの少なくとも 1つを含むパラメータに基づき前記補正データ を算出することを特徴とする超音波イメージング方法。
[16] 請求項 10に記載の超音波イメージング方法において、前記補正データ算出工程 は、一の前記超音波ビームのビーム形状と前記被検体の深度とに対応して一の前記 補正データを算出する工程と、前記一の超音波ビームに対応する受信信号が受信 処理される間に、前記一の超音波ビームとは異なる他の前記超音波ビームのビーム 形状と前記被検体の深度とに対応して他の前記補正データを算出する工程を含む ことを特徴とする超音波イメージング方法。
[17] 請求項 10に記載の超音波イメージング方法において、前記表示工程は、前記超 音波ビームのビーム形状に対応して算出されたダイナミックレンジを前記被検体の深 度の関数としてグラフ表示する工程を含むことを特徴とする超音波イメージング方法
[18] 請求項 10に記載の超音波イメージング方法において、前記表示工程は、前記超 音波ビームのビーム形状に対応して算出されたビーム強度を前記被検体の深度の 関数としてグラフ表示する工程を含むことを特徴とする超音波イメージング方法。
PCT/JP2004/018260 2003-12-19 2004-12-08 超音波診断装置および超音波イメージング方法 WO2005060832A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005516451A JP4664209B2 (ja) 2003-12-19 2004-12-08 超音波診断装置およびその撮像を実行する超音波イメージングプログラム

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2003-422297 2003-12-19
JP2003422297 2003-12-19

Publications (1)

Publication Number Publication Date
WO2005060832A1 true WO2005060832A1 (ja) 2005-07-07

Family

ID=34708731

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2004/018260 WO2005060832A1 (ja) 2003-12-19 2004-12-08 超音波診断装置および超音波イメージング方法

Country Status (2)

Country Link
JP (1) JP4664209B2 (ja)
WO (1) WO2005060832A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009022656A (ja) * 2007-07-23 2009-02-05 Aloka Co Ltd 超音波診断装置

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102188148B1 (ko) * 2014-01-17 2020-12-07 삼성메디슨 주식회사 광음향 영상 장치 및 광음향 영상 디스플레이 방법

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS612207U (ja) * 1984-06-08 1986-01-08 株式会社日立メディコ 超音波断層装置
JPH05245137A (ja) * 1992-03-06 1993-09-24 Yokogawa Medical Syst Ltd 超音波診断装置
JPH1128210A (ja) * 1997-05-07 1999-02-02 General Electric Co <Ge> 超音波散乱物質の三次元イメージング・システムおよび方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02177948A (ja) * 1988-12-28 1990-07-11 Olympus Optical Co Ltd 超音波診断装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS612207U (ja) * 1984-06-08 1986-01-08 株式会社日立メディコ 超音波断層装置
JPH05245137A (ja) * 1992-03-06 1993-09-24 Yokogawa Medical Syst Ltd 超音波診断装置
JPH1128210A (ja) * 1997-05-07 1999-02-02 General Electric Co <Ge> 超音波散乱物質の三次元イメージング・システムおよび方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009022656A (ja) * 2007-07-23 2009-02-05 Aloka Co Ltd 超音波診断装置

Also Published As

Publication number Publication date
JP4664209B2 (ja) 2011-04-06
JPWO2005060832A1 (ja) 2007-12-13

Similar Documents

Publication Publication Date Title
EP1832894B1 (en) System and method for automatic gain compensation based on image processing
JP4585326B2 (ja) 超音波撮像装置及び超音波撮像方法
EP1426010A1 (en) Ultrasonographic apparatus, ultrasonographic data processing method, and ultrasonographic data processing program
JP2004215987A (ja) 超音波診断装置および超音波診断方法
US10231709B2 (en) Ultrasound diagnostic apparatus, signal processing method for ultrasound diagnostic apparatus, and recording medium
JP3645347B2 (ja) 超音波診断装置及び遅延時間最適化方法
CN108024798B (zh) 超声波诊断装置以及延迟数据生成方法
JP2007236740A (ja) 超音波診断装置及びその制御プログラム
JP5498551B2 (ja) 超音波診断装置及び超音波送受信条件最適化プログラム
JP5132089B2 (ja) 超音波診断装置、超音波送受信条件適正化プログラム及び超音波送受信条件適正化方法
JP2011212440A (ja) 超音波信号取得及び処理の方法及び装置
KR20160099877A (ko) 초음파 영상 화질 개선 방법 및 이를 이용한 초음파 영상 장치
JP2003010180A (ja) 超音波撮像装置
JP4664209B2 (ja) 超音波診断装置およびその撮像を実行する超音波イメージングプログラム
JP2010201110A (ja) 超音波診断装置及び超音波診断装置制御方法
JP4599408B2 (ja) 超音波診断装置
JP3784109B2 (ja) 超音波映像装置
JP2003339698A (ja) 超音波診断装置
US10925578B2 (en) Acoustic wave image generating apparatus and control method thereof
JP5190248B2 (ja) 超音波診断装置
US9291601B2 (en) Ambient sound velocity obtaining method and apparatus
US11583247B2 (en) Acoustic wave image generating apparatus and control method thereof
JP2004195091A (ja) 超音波診断装置
JP3388954B2 (ja) 超音波診断装置
JP5393568B2 (ja) 超音波撮像方法および超音波撮像装置

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BW BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE EG ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NA NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): BW GH GM KE LS MW MZ NA SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LT LU MC NL PL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

DPEN Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed from 20040101)
121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2005516451

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

WWW Wipo information: withdrawn in national office

Country of ref document: DE

122 Ep: pct application non-entry in european phase