WO2005087109A1 - 超音波撮像装置 - Google Patents

超音波撮像装置 Download PDF

Info

Publication number
WO2005087109A1
WO2005087109A1 PCT/JP2005/000047 JP2005000047W WO2005087109A1 WO 2005087109 A1 WO2005087109 A1 WO 2005087109A1 JP 2005000047 W JP2005000047 W JP 2005000047W WO 2005087109 A1 WO2005087109 A1 WO 2005087109A1
Authority
WO
WIPO (PCT)
Prior art keywords
transmission
signal
pulse
ultrasonic
reception
Prior art date
Application number
PCT/JP2005/000047
Other languages
English (en)
French (fr)
Inventor
Takashi Azuma
Shin-Ichiro Umemura
Hiroshi Kuribara
Tatsuya Hayashi
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 EP05703319.3A priority Critical patent/EP1726261B1/en
Priority to US10/586,665 priority patent/US8002703B2/en
Priority to JP2006510881A priority patent/JP4415011B2/ja
Publication of WO2005087109A1 publication Critical patent/WO2005087109A1/ja

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/02Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
    • G01S15/06Systems determining the position data of a target
    • G01S15/08Systems for measuring distance only
    • G01S15/10Systems for measuring distance only using transmission of interrupted, pulse-modulated waves
    • G01S15/102Systems for measuring distance only using transmission of interrupted, pulse-modulated waves using transmission of pulses having some particular characteristics
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/481Diagnostic techniques involving the use of contrast agent, e.g. microbubbles introduced into the bloodstream
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8959Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8959Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes
    • G01S15/8963Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using coded signals for correlation purposes using pulse inversion
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target
    • G01S7/52039Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target exploiting the non-linear response of a contrast enhancer, e.g. a contrast agent
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52019Details of transmitters
    • G01S7/5202Details of transmitters for pulse systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52036Details of receivers using analysis of echo signal for target characterisation
    • G01S7/52038Details of receivers using analysis of echo signal for target characterisation involving non-linear properties of the propagation medium or of the reflective target

Definitions

  • the present invention relates to an ultrasonic technology for imaging the inside of a living body by transmitting and receiving ultrasonic waves to and from a living body, and more particularly to an ultrasonic imaging technology for imaging using a microbubble based contrast agent.
  • An ultrasonic imaging apparatus that transmits and receives pulsed ultrasonic waves to a living body to image the inside of the living body is widely used for medical diagnosis.
  • contrast agents have previously been used for imaging of the vascular system and the like.
  • the purpose is to administer a contrast agent into the blood to obtain an image with enhanced contrast of the structure and distribution of the vasculature, and diagnose with high accuracy the diseases reflected in the vasculature such as malignant tumors and infarcts. It is.
  • micro bubbles having a size of micron 'order have become stable in some way
  • contrast agents with formulated formulations, they are beginning to be widely used.
  • microbubbles with a diameter of about 1 micron resonate with ultrasonic waves of several MHz used for ultrasonic diagnosis and vibrate with large amplitude, and as a result, ultrasonic waves in this frequency range are well scattered and contrasted. It takes advantage of the ability to
  • a characteristic of microbubble ultrasound contrast agents is their strength and non-linearity. This is because the micro-bubbles have the property that the volume increase under negative pressure is much larger than the volume decrease under positive pressure of the same amplitude. For this reason, the echo signal scattered by the microbubbles contains a large amount of the second harmonic component having a frequency twice that of the transmission signal. V ⁇ . Newhouse et al. First reported in 1992 the method of obtaining an enhanced blood flow dobbler signal for soft tissue based on this second harmonic component (see, for example, Non-Patent Document 1). .
  • W. Wilkening has also proposed a method of performing transmission / reception N times by rotating the phase of the transmission sound pressure pulse waveform by 360 ° ZN (see, for example, Non-Patent Document 2).
  • N 3
  • the following process is performed.
  • ultrasound imaging is performed on a phantom region divided in advance into a spatial region in which contrast agent echo is dominant and a spatial region in which tissue reflection echo is dominant, and signals are acquired.
  • the filter coefficients at the time of the above-mentioned addition are determined by the method of least squares so that the difference between the tissue reflection echo and the contrast agent signal is maximized.
  • Non-Patent Document 1 1992 IEEE Ultrasonics Symposium Proceedings, pp. 1175- 1177
  • Non Patent Literature 2 2001 IEEE Ultrasonics Symposium Proceedings, pp. 1733-1737
  • Non Patent Literature 3 2003 IEEE Ultrasonics Symposium Proceedings, pp. 429-432
  • Non Patent Literature 4 1996 IEEE Transactions on Ultrasonics, Ferroelectrics, and
  • Patent Document 1 US Patent No. 6,095, 980
  • Non-Patent Document 1 proposes a method for obtaining a blood flow dobbler signal emphasized for soft tissue based on the second harmonic component. If only the band pass filter is used to extract the echo signal power of the second harmonic component, the problem remains that the pulse of the second harmonic echo signal obtained as the output signal becomes long.
  • Patent Document 1 for the purpose of solving this problem, two echoes obtained by transmitting and receiving twice using a transmission sound pressure pulse waveform whose polarity is reversed are used.
  • Patent Document 2 For the purpose of solving this problem, two echoes obtained by transmitting and receiving twice using a transmission sound pressure pulse waveform whose polarity is reversed are used.
  • an echo signal having a harmonic component returns to the soft tissue force even in the absence of the microbubble contrast agent.
  • This harmonic component The method of forming an echo image based on the above is called Tissue Harmonic Imaging method, and it is preferred that the acoustic S / N ratio is higher than that of the echo image by the fundamental wave component, and has recently become commonly used. The This is because even if the pulse inversion method is used, it is possible to obtain an echo signal in which harmonic components generated by scattering by the microbubble contrast agent and harmonic components generated by the propagation of the transmission pulse are mixed. This means that it is difficult to distinguish the contrast agent that was originally intended from soft tissue.
  • Non-linear scattering by a microbubble based contrast agent tends to be generally observed even at low sound pressure as compared to non-linear propagation in soft tissue. Therefore, it is widely used to form an echo image mainly by the non-linear component from the micro-bubble-based agent while suppressing the generation of the Tissue Harmonic component by the pulse inversion method in which the transmission sound pressure is suppressed low. Since the signal amplitude is not sufficient, a high echo image having an S / N ratio sufficient to make a definitive diagnosis expected for contrast diagnosis can not be obtained at present.
  • Non-patent Document 2 As a method of solving this problem to some extent, in the prior art (Non-patent Document 2), as described above, a method of transmitting / receiving N times with a pulse waveform phase-rotated by 360 ° ZN is proposed. .
  • this method is a phantom-optimized filter, it is not necessarily appropriate for different signals from one test subject to another.
  • the above-mentioned least squares filter is used. It is difficult to optimize for all depths and focus conditions.
  • the present invention clearly distinguishes an echo component generated by being scattered by a microbubble-based contrast agent from a Tissue Harmonic component generated by nonlinear propagation of a transmission pulse and displays an image. Therefore, it is an object of the present invention to provide an ultrasound imaging technique for achieving a high contrast angiogram with a high S / N ratio sufficient to make a definitive diagnosis based on a contrast echo image.
  • the transmission sound pressure pulse waveform is shared by common, and the phase of its carrier wave is rotated by 120 ° at three times for transmission and reception.
  • the fundamental wave component and the second harmonic component of the echo signal scattered by the linear scatterer which is negligible, are canceled simultaneously.
  • the phase of the fundamental wave component is natural for the phase of the fundamental wave component to rotate by 120 ° at a time, and the phase of the second harmonic component is also opposite to the fundamental wave component. In order to rotate 120 degrees in the direction.
  • FIG. 1 is a diagram for explaining the principle of the pulse inversion method and the principle of vibration of an in-line four-cylinder four-stroke engine.
  • FIG. 1 (a) shows the phase relationship of the fundamental wave, and this configuration cancels out the fundamental wave components generated by each pair of pistons
  • FIG. 1 (b) shows the second harmonic wave.
  • the phase relationship is shown, and the second harmonic component grows twice, and as a result, vibration with a frequency twice the crankshaft rotational speed becomes a problem.
  • the solid line indicates the phase of vibration by the first set of pistons
  • the dotted line indicates the phase of vibration by the second set of pistons.
  • a four-stroke in-line six-cylinder engine is usually configured such that three pairs of two in-phase pistons symmetrically arranged have a crank angle of 120 °.
  • FIG. 2 is a diagram for explaining the principle of the 3-pulse method of the present invention and the principle of vibration of an in-line six-cylinder four-stroke engine.
  • FIG. 2 (a) shows the phase relationship of the fundamental wave, and with this configuration, the fundamental wave components generated by the respective pairs of pistons are generated so as to form 120 ° in terms of crank angle and strike each other. Erase each other.
  • the alternate long and short dash line indicates the phase of vibration due to the third set of pistons.
  • the configuration of the in-line four-cylinder engine in which the second harmonic vibration is emphasized in principle corresponds to the pulse inversion method, and in the in-series six-cylinder engine, not only the fundamental wave but also the second harmonic vibration is canceled in principle.
  • the configuration corresponds to the method of the present invention.
  • the phase of the echo signal generated by scattering by the microbubble contrast agent is influenced by the envelope amplitude due to its strong and non-linear resonance characteristics, and is constant with respect to the phase of the transmission signal carrier. It does not matter. Therefore, even if three echo signals obtained by rotating the phase of the transmission pulse carrier by 120 ° and performing three times of transmission and reception are added, echo signals scattered and generated by contrast bubbles are generated. There is a component that can not be canceled. Therefore, since the remaining echo signal component reflects only the presence of the microbubble-based contrast agent, it can be used to perform ultrasound imaging in which the contrast agent is clearly distinguished from soft tissue.
  • the ultrasonic imaging apparatus is an ultrasonic imaging apparatus that forms an image of the inside of a living body by transmitting and receiving ultrasonic pulses to the living body to which a microbubble for contrast has been introduced.
  • N is an integer of 3 or more
  • ultrasonic echoes from soft tissue and weave of a living body are performed by performing N times of transmission and reception using transmission pulse waves having different waveforms under the same transmission and reception focus conditions. It is configured to suppress the transmission / reception sensitivity to the component from the fundamental wave of the signal to the (N-1) -th harmonic and obtain the transmission / reception sensitivity to the ultrasonic echo signal of the microbubble force for contrast.
  • the ultrasonic imaging apparatus has a transmission amplifier for transmitting a transmission pulse wave into the living body, and the input cycle of the signal given to the transmission amplifier is N for the maximum frequency of the frequency components of the transmission pulse. Is an integral multiple of
  • a DZA converter is provided as a means to provide a signal to the transmission amplifier, and the signal output cycle of the DZA converter is This is an integral multiple of N for the largest frequency among frequency components of the loop.
  • the transmission pulse wave is a waveform obtained by adding the fundamental wave and the second harmonic to the fundamental wave.
  • the ultrasonic imaging apparatus of the present invention is an ultrasonic imaging apparatus that transmits and receives ultrasonic pulses to and from a living body into which a microbubble for contrast has been introduced, and forms an image of the inside of the living body.
  • the fundamental wave power of the ultrasonic echo signal from the soft tissue of the living body up to the second harmonic
  • An imaging sequence that suppresses transmission / reception sensitivity to components and selectively obtains transmission / reception sensitivity to an ultrasound echo signal from an imaging microbubble, and transmission of different waveforms under the same transmission / reception focus conditions as described above
  • the transmission and reception sensitivity to the fundamental wave component of the ultrasonic echo signal from the soft tissue yarn of the living body is suppressed, and the transmission and reception to the second or higher order nonlinear signal is performed.
  • the transmission amplitude is configured to be different in the sequence in which transmission / reception is performed three times and the sequence in which transmission / reception is performed twice.
  • the transmission amplitude in the sequence in which transmission / reception is performed three times is made larger than the transmission amplitude in the sequence in which transmission / reception is performed twice.
  • the ultrasonic imaging apparatus of the present invention transmission and reception of ultrasonic pulses to the living body, force to form a contrast-enhanced image by the imaging microbubbles inside the living body, same transmission and reception focus conditions.
  • the signals are added to form a biomedical non-linear image, and the two obtained images are superimposed and displayed.
  • echo components scattered and generated by the microbubble contrast agent can be clearly distinguished from those of the Tissue Harmonic component generated by the non-linear propagation of the transmission pulse, and a high contrast contrast echo image of the S / N ratio can be obtained.
  • a signal component derived only from an imaging microbubble is extracted from the echo signal without including a signal component derived from non-linear propagation in soft tissue, etc.
  • a background representing the position and morphology of soft tissue as a signal representing the spatial distribution of bubbles It can be superimposed on the top in a distinguishable color tone and displayed in an easy-to-understand manner by a two-dimensional or three-dimensional image.
  • the medical diagnostic utility of the device embodying the present invention is so great that the meaning of the present invention in the industry supporting medical diagnostics is also great.
  • FIG. 3 is a typical example of a block diagram of an ultrasonic diagnostic apparatus configured to practice the present invention.
  • Each element constituting the ultrasonic probe 1 is connected to the transmission beam former 3 and the reception beam former 10 via the switching switch group 2.
  • the transmission beam former 3 has directivity when transmitted through each element using the waveform selected and read from the transmission waveform memory 5 by the transmission waveform selection unit 4 under the control of the transmission sequence control unit 6.
  • a signal is generated that results in an ultrasonic pulse.
  • This signal is converted into an ultrasonic pulse by each element of the ultrasonic probe 1 and transmitted to the living body.
  • the ultrasound echo signal reflected or scattered in the living body and returned to the ultrasound probe 1 is received by each element and converted into an electrical signal.
  • the reception signals having directivity are generated by adding delay time to each reception signal to be generated.
  • the time-series signals obtained by the delay addition are also written to one bank in the reception waveform memory 12 selected by the reception wave memory selection unit 11 under the control of the transmission / reception sequence control unit 6 and added together.
  • the output signal of the adder 13 passes through a band pass filter A14 for removing noise components, is converted into an envelope signal in an envelope signal detector A16, and is input to a scan converter 18.
  • a part of the time-series signal written to the reception waveform memory 12 is read as it is.
  • the signal passes through a band pass filter B15 that removes noise components that are not added to each other, is converted into an envelope signal by an envelope signal detector B17, and is input to a scan converter 18.
  • the scan converter 18 appropriately generates signals and controls them so as to superimpose a plurality of input signals and display them in a two-dimensional display unit 19.
  • FIGS. 15 and 16 are diagrams for explaining the relationship between the envelope of transmission pulses and the carrier wave of each pulse in the 3-pulse method of the present invention.
  • envelope signals (indicated by dotted lines) are superimposed on FIGS. 15 (b), 15 (c) and 15 (d), respectively.
  • FIG. 16 (a), 16 (b) and 16 (c) envelope signals (indicated by dotted lines) are superimposed on FIGS. 15 (b), 15 (c) and 15 (d), respectively.
  • FIG. 4 is a diagram showing an example of echo signals by point reflectors in a soft tissue of the living body having non-linear propagation characteristics, obtained by the 3-pulse method of the present invention.
  • the transmission waveform memory 5 has a common envelope signal (Envelope) as shown in FIG. 15 (a), and is shown in FIGS. 15 (b), 15 (c) and 15 (d) respectively.
  • envelope envelope
  • the signals written to each bank in the reception waveform memory 12 in this case are shown in Fig. 4 (a), Fig. 4 (b) and Fig. 4 (c).
  • the received echo signal generated as a result of the transmission ultrasonic pulse propagating through the soft tissue is reflected by one point reflector, and is obtained by numerical simulation and shown.
  • the carrier frequency is 2 MHz.
  • FIG. 4 (d) is an output signal obtained when the signals shown in FIGS. 4 (a), 4 (b) and 4 (c) are input to the adder 13. Since the transmission ultrasonic pulse propagates nonlinearly in living soft tissue, the signals shown in FIGS. 4 (a), 4 (b) and 4 (c) include not only the fundamental wave component but also the second harmonic component. However, as predicted from the principle of the present invention, as shown in FIG. 4 (d), adder output results, the second harmonic components that cancel each other out only by the fundamental wave components cancel each other, and the signal The amplitude is almost zero. To compare with this, Fig. 5 shows the results when the pulse inversion method is used.
  • FIG. 5 is a diagram showing an example of an echo signal by a point reflector in living soft tissue with nonlinear propagation characteristics obtained by the pulse inversion method.
  • Fig. 5 (a) and Fig. 5 (b) show the signals written to each bank in the received waveform memory 12 when this sequence is executed twice while changing the waveform.
  • Fig. 5 (c) shows the output signal of the adder 13 obtained at that time, but the fundamental wave components are mutually canceled, but the signals of the result of emphasizing the second harmonic components are obtained. ing.
  • This signal is called a Tissue Harmonic signal, and it has the advantage of obtaining a high acoustic signal-to-noise ratio when depicting living soft tissue tissue, but it distinguishes only the distribution and dynamics of the contrast agent from the soft tissue. If you want to show it, it will be the biggest factor that disturbs it.
  • FIG. 6 is a diagram showing an example of the scattering coefficient signal by the imaging microbubble obtained by the three pulse method of the present invention.
  • FIG. 8 is a diagram showing an example of a scattered echo signal by a microbubble for contrast obtained by the pulse inversion method.
  • FIG. 7 is a view showing a waveform after passing through a band pass filter of a scattering coefficient signal by a microbubble for contrast obtained by the three pulse method of the present invention.
  • FIG. 9 is a diagram showing a waveform after passing through a band pass filter of a scattered echo signal by a microbubble for contrast obtained by the pulse inversion method.
  • Fig. 6 and Fig. 8 show the result of numerical calculation simulation of the reception echo signal generated by scattering by the micro bubble of radius 1.5 ⁇ m
  • Fig. 7 and Fig. 9 show the results at that time.
  • Bandpass filter A14 input / output signal.
  • the scale of the vertical axis, which is proportional to the sound pressure, is common throughout FIGS. 6, 7, 8 and 9.
  • Figure 7 (b), Figure 7 (c), Figure 7 (d), Figure 7 (e) and Figure 7 (f) respectively show the center frequency of the passband as the fundamental frequency (2 MHz), 1.5 harmonics Frequency (3 MHz), Second harmonic frequency (4 MHz), 2.5 This is the output signal of the bandpass filter A14 when it is set to the double harmonic frequency (5 MHz) and the third harmonic frequency (6 MHz).
  • FIGS. 8 (a) and 8 (b) the signal powers written to each bank in the reception waveform memory 12 are shown in FIGS. 8 (a) and 8 (b).
  • Figures 9 (b), 9 (c), 9 (d), 9 (e) and 9 (f) show the band pass filters when the pass band center frequency is set as in the case of Figure 8. It is an output signal of A14.
  • the 1.5 th harmonic power As a matter of course, it is possible to obtain a microbubble-derived signal containing a large amount of the second harmonic component, according to the invention of the present invention designed to cancel out the component caused by nonlinear propagation etc. in the second harmonic component of the echo signal. As shown in FIG. 7 (d), FIG. 7 (e) and FIG. 7 (f) according to the transmission / reception sequence, a signal of sufficient amplitude derived from microbubbles containing a large number of second to third harmonic components. It should be noted that is obtained.
  • the microbubble force is a resonator with large non-linearity, more generally speaking, it has a response characteristic that the delay time depends on the amplitude. is there. That is, if there is non-linearity between the input and output sound pressure and there is no amplitude dependency in the delay response time, as shown in FIG. 4 (d), the second harmonic component in the output signal is It will be defeated. On the other hand, in the case of a simple linear resonator, the second harmonic component itself is not generated, so it is out of the question.
  • Non-Patent Document 4 is considered to be useful also in the case of ultrasonic imaging using a microbubble for contrast.
  • FIG. 10 shows an example of echo signals by point reflectors in soft tissue of a living body obtained by intentionally superimposing a second harmonic on a transmission pulse wave in the 3-pulse method of the present invention.
  • Figure It is.
  • FIG. 11 is a view showing an example of a scattered echo signal by the imaging microbubbles obtained when the second harmonic is intentionally superimposed on the transmission pulse wave in the three pulse method of the present invention. .
  • FIG. 12 shows that after passing a band pass filter of the scattered echo signal by the imaging microbubbles obtained when the second harmonic is intentionally superimposed on the transmission pulse wave in the three pulse method of the present invention.
  • FIGS. 10, 11, and 12 show examples of echo signals obtained when the second harmonic component is intentionally superimposed on the transmission pulse waveform.
  • the transmission waveform selection unit 5 has a common envelope signal, three fundamental pulse waves as carrier waves, and two supersonic pulse waveforms with different phases of the second harmonic by 120 ° in the transmission waveform memory 5.
  • the signal to be written to each bank in the reception waveform memory 12 is determined similarly to the case of FIG. 10 (a), 10 (b) and 10 (c), and the output signal of the adder 13 obtained at that time is shown in FIG. 10 (d).
  • reception echo signals generated by scattering by the imaging microbubbles are determined as in the case of FIG. 6, and the signals written to each bank in the reception waveform memory 12 are shown in FIG. 11 (a), FIG.
  • . 12 (b), 12 (c), 12 (d), 12 (e) and 12 (f) show the band pass filter when the pass band center frequency is set as in the case of FIG. It is an output signal of A14.
  • the received echo signal generated by being scattered by the point reflector is added to the adder 1 Similar to the second harmonic components generated by the nonlinear propagation in the case of FIG. 4, the second harmonic components of the input signal of No. 3 cancel each other, and the amplitude of the output signal of the adder is almost zero.
  • the received echo signal generated by scattering by the imaging microbubbles as in the case of FIG. 6, the sufficient amplitude including many second to third harmonic components which can not be canceled out by the addition by the adder 13. The output signal of is obtained.
  • FIG. 13 shows an example of an adder output signal waveform of an echo signal by a point reflector in living soft tissue obtained when there is an error in the phase of a transmission pulse in the 3-pulse method of the present invention.
  • the output signal of the adder 13, ie, the input signal of the band pass filter A14, and the center frequency of the passband are the fundamental wave and the second harmonic.
  • the filter output signal in the case of fitting to the wave was determined as in the case of FIG. 4 and is shown in FIG. 13 (b) and FIG. 13 (c).
  • FIG. 13 (a) shows a signal in the case of FIG. 5, that is, when the pulse inversion method can be performed without a phase error.
  • the phase of the third transmit pulse wave is shifted by 10 ° corresponding to a phase shift of 20 ° of the second transmit pulse wave so that the sum signal of the three transmit pulse waves becomes zero. In the case of adjustment.
  • FIG. 13 (c) is the result of performing such correction of the third pulse, and the fundamental wave component in the output signal of the adder 13 can be canceled despite the phase error of the second pulse.
  • N reception echo signals Because (Equation 4) holds, and generally, the sum signal of N pulse waveforms used for transmission is substantially zero. The fundamental wave component in the signal is cancelled.
  • FIG. 14 shows the phase of the amplitude signal of the output signal of the adder of the echo signal by the point reflector in the soft tissue of the living body obtained when there is an error in the phase of the transmission pulse in the 3-pulse method of the present invention. It is a figure explaining error dependence.
  • FIG. 14 shows the phase error (degree) (horizontal axis) in which the peak-to-peak value (output signal relative amplitude value PP value) (vertical axis) of the amplitude of the output signal of the adder 13 is given to the second pulse. Shown as a function of axis).
  • the signal amplitude normalized by the value when the pulse inversion method can be performed without a phase error is corrected when the third pulse is not corrected ((a) in FIG. 14) and when the correction is performed (figure). It showed about (b) in 14).
  • the phase error reaches 20 °, even if the third pulse is corrected, the signal amplitude force derived from nonlinear propagation in the soft tissue and not derived from the imaging microbubbles is more than half that of the conventional pulse inversion method.
  • the effects of the present invention will not be fully exhibited. From this result, it can be said that it is desirable to set the phase error of the transmission pulse to about 10 degrees or less in order to sufficiently obtain the effects of the present invention.
  • the echo signal contains no signal component derived from non-linear propagation in soft tissue and the like, and the signal derived only from the imaging microbubble is It can be extracted.
  • Such a signal can be obtained as an output signal of the adder 13, and a signal with an improved S / N ratio can be obtained as an output signal of the band pass filter A14, and an envelope signal can be obtained as the envelope signal.
  • the output signal of the line signal detector A16 is obtained and input to the scan converter 18 as a signal representing the spatial distribution of the imaging microbubbles.
  • envelope signal detector B 17 an envelope signal is detected by envelope signal detector B 17 from a signal in which the signal written in one bank in received waveform memory 12 is passed through band pass filter B 15 and the S / N ratio is improved. This is input to the scan converter 18 as a signal representing the position and form of soft tissue.
  • the scan converter 18 uses the output signal of the envelope signal detector B 17 as a background, and the output signal of the envelope signal detector A 16 is superimposed by a different color tone so as to be convenient to distinguish it from the display portion 19 Is displayed on. In this way, the distribution of the imaging microbubbles in the patient's body to be examined can be easily understood by the two-dimensional or three-dimensional image. It can be displayed.
  • the case where the signal written in one bank in the received waveform memory 12 is used as it is as a signal representing the position and form of the soft tissue is generally described. It is possible to use summing signals with appropriate weighting to the signals written to the three banks in memory 12.
  • N integers of 3 or more having a common envelope signal are used.
  • the transmit pulse wave the phase of the carrier is rotated by 360 ° ZN at a time to transmit and receive N times, and the N echo signals obtained are once written to N banks in the reception waveform memory 12 and read out.
  • the present invention can be implemented by inputting the above signal to the adder 13.
  • FIG. 17 is a diagram showing an example of input / output characteristics of the transmission amplifier.
  • FIG. 18 is a diagram for explaining the relationship between the output waveform of the DZA converter and the sampling point.
  • FIG. 17 is a diagram showing the input / output characteristics of a general transmission amplifier. The horizontal axis is the input voltage, and the vertical axis is the output voltage. Although it is desirable that the output is proportional to the input, it is actually nonlinear as shown in FIG. The nonlinear characteristics of this amplifier have a major impact on 3-pulse transmission.
  • FIG. 18 shows the output waveforms of the DZA converter of three pulses of 0 °, 120 °, and 120 ° as solid lines, broken lines, and dotted lines, respectively.
  • the circled points indicate the sampling points of the DZA converter.
  • the vertical axis represents the output value (relative value) of the DZA converter, and the vertical axis represents a dimensionless number indicating the timing at which the DZA converter is provided with an output for control.
  • Figure 18 (a) shows the center when sampled at four times the center frequency
  • Figure 18 (b) shows the center.
  • Fig. 18 (c) shows the case where sampling is performed at a frequency eight times the center frequency.
  • the output value of DZA takes a different value for each pulse. In these cases, the output value does not become as set because of the nonlinearity of the transmission amplifier shown in Fig.17.
  • the signal is output from the DZA converter at a sampling frequency of 6 times the center frequency (in the case of N pulses, an integer multiple of N).
  • the output value of DZA does not differ greatly for each pulse, so there is no problem even if the transmission amplifier has nonlinearity.
  • FIG. 19 is a view for explaining the transmission amplitude dependency of the contrast ratio between the contrast agent signal and the biononlinear signal and the sensitivity of the contrast agent signal.
  • Fig. 19 (a) is a diagram in which the discrimination ratio of the contrast agent signal and the biological nonlinear signal in the pulse inversion method and the three pulse method is plotted against the amplitude of the transmitted signal
  • Fig. 19 (b) is the contrast agent signal Is plotted against the amplitude of the transmitted signal.
  • the horizontal axis indicates transmission pressure in terms of sound pressure amplitude X O.
  • the vertical axis indicates (a) discrimination ratio (Sensitivity, dB), (b) Sensitivity (dB),
  • the solid line indicates the 3-pulse method, and the dotted line indicates the pulse inversion method.
  • the discrimination ratio the 3-pulse method is always 20 dB better and 40 dB better than the conventional pulse inversion method.
  • the sensitivity of the 3-pulse method decreases as the transmission amplitude decreases. This is because a large transmission amplitude, that is, when the nonlinearity is large, the phase relationship between the signal from the contrast agent and the transmission signal does not collapse. However, when imaging with a large transmission amplitude to obtain sensitivity, the contrast agent is destroyed during imaging, and eventually it becomes difficult to maintain a state of high sensitivity.
  • the present invention avoids this problem in the following two ways.
  • One is by using a waveform in which the second harmonic is superimposed as shown in FIG.
  • the contrast agent is more easily broken when pulled than when pressed, the sound pressure is large by using a waveform in which the second harmonic is superimposed so that the negative pressure is smaller than the positive pressure.
  • Relatively broken contrast agent It is possible to devise to be
  • a method by sequence control is also effective. This alternates transmission and reception of pulse reverse transfer wave and 3-pulse method. Normally, the sound pressure is weak, the pulse inversion method is used for monitoring, and sometimes the high sound pressure 3-pulse method is used.
  • FIG. 20 and FIG. 21 are control flow diagrams of the switching imaging sequence of the pulse inversion method and the three pulse method.
  • Pulse inversion (reversal) imaging processing 20 performs the end of N imaging operations of 20 times or performs switching determination processing 21 by external input and performs 3-pulse imaging processing 22.
  • 3-pulse imaging processing 22 M imaging ends of end or external It performs decision processing 23 of switching by input.
  • FIG. 20 shows a method of switching between pulse inversion (inversion) and 3-pulse method for each frame.
  • switching is performed for each raster.
  • a process 104 is performed, and it is determined whether or not one frame is completed.
  • a raster movement is performed after one frame is completed. According to this method, even when there is movement of the subject, the switching time is sufficiently short, so that it is possible to regard the targets of the signals imaged by the pulse inversion method and the 3-pulse method as equivalent. Therefore, it becomes possible to superimpose the contrast agent signal by the 3-pulse method with another color code on the black-and-white image imaged by the pulse inversion method.
  • an ultrasonic imaging apparatus having a high SZN ratio sufficient to perform a definitive diagnosis based on a contrast-enhanced echo image.
  • FIG. 1 illustrates the principle of pulse inversion method and the principle of vibration of an in-line four-cylinder four-stroke engine.
  • FIG. 2 is a view for explaining the principle of the 3-pulse method of the present invention and the principle of vibration of an in-line six-cylinder four-stroke engine. 3) It is a block diagram showing the configuration of the ultrasonic diagnostic apparatus of the embodiment of the present invention.
  • FIG. 5 is a view showing an example of an echo signal by a point reflector in a living soft tissue having nonlinear propagation characteristics, obtained by a pulse inversion method.
  • FIG. 6 is a view showing an example of a scattered echo signal by an imaging microbubble obtained by the three pulse method of the present invention.
  • FIG. 7 is a view showing a waveform after passing through a band pass filter of a scattered echo signal by an imaging microbubble obtained by the three pulse method of the present invention.
  • FIG. 8 is a view showing an example of a scattered echo signal by a microbubble for contrast obtained by a pulse inversion method.
  • FIG. 9 is a view showing a waveform after passing through a band pass filter of a scattered echo signal by a microbubble for contrast obtained by a pulse inversion method.
  • FIG. 10 is a view showing an example of an echo signal by a point reflector in the soft tissue of a living body obtained when the second harmonic is intentionally superimposed on the transmission pulse wave in the 3-pulse method of the present invention. Ru.
  • ⁇ 11 In the three pulse method of the present invention, it is a figure showing an example of the scattered echo signal by the microbubble for contrast obtained when the second harmonic is intentionally superimposed on the transmission pulse wave.
  • FIG. 12 In the 3-pulse method of the present invention, the waveform after passing through the band pass filter of the scattered echo signal by the contrast microbubble obtained when the second harmonic is intentionally superimposed on the transmission pulse wave.
  • FIG. 13 An example of an adder output signal waveform of an echo signal by a point reflector in living soft tissue obtained when there is an error in the phase of transmission pulse in the 3-pulse method of the present invention
  • FIG. 14 In the 3-pulse method of the present invention, the phase of the amplitude signal of the output signal of the adder of the echo signal by the point reflector in the soft tissue of the living body obtained when there is an error in the phase of the transmission pulse. It is a figure explaining error dependence.
  • FIG. 15 is a view for explaining the relationship between the envelope of transmission pulses and the carrier wave of each pulse in the 3-pulse method of the present invention.
  • FIG. 16 is a view for explaining the relationship between the envelope of a transmission pulse and the carrier wave of each pulse in the 3-pulse method of the present invention.
  • FIG. 17 is a diagram showing an example of input / output characteristics of a transmission amplifier.
  • FIG. 18 is a diagram for explaining the relationship between an output waveform of a DZA converter and sampling points.
  • FIG. 19 is a diagram for explaining the discrimination ratio between the contrast agent signal and the biological non-linear signal, and the transmission amplitude dependency of the sensitivity of the contrast agent signal.
  • FIG. 20 A control flow diagram of a switched imaging sequence of the pulse inversion method and the three pulse method.
  • FIG. 21 is a control flow diagram of a switching imaging sequence of pulse inversion method and 3-pulse method. Explanation of sign

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Nonlinear Science (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biomedical Technology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Biophysics (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

 マイクロバブル系造影剤により散乱されて生ずるエコー成分を、送信パルスの非線形伝播により生ずるTissue Harmonic成分と明確に峻別して映像化する超音波撮像装置を提供するため、生体に超音波パルスの送受波を行ない、生体内部の造影用マイクロバブルによる造影画像を形成する超音波撮像装置であって、同一の送受波フォーカス条件下で、包絡線信号を共通とする送信パルス波を用いその搬送波の位相を(a)=0°、(b)=120°、(c)=−120°、(d)=180°の4回の送受波を行ない、(a)、(b)、(c)の3個の時系列受信エコー信号を加算し造影画像を形成し、(a)、(d)の2個の時系列受信エコー信号を加算し生体伝播非線形画像を形成し、2種類の画像を重畳表示する。

Description

明 細 書
超音波撮像装置
技術分野
[0001] 本発明は、超音波を生体に対して送受信することにより、生体内部を撮像する超音 波技術に係り、特に、マイクロバブル系造影剤を利用して撮像する超音波撮像技術 に関する。
背景技術
[0002] 生体に対してパルス状超音波を送受信し、生体内部を撮像する超音波撮像装置 は、医療診断に広く用いられている。
[0003] 画像診断モダリティの中でも X線や MRIの分野では、血管系の撮像等に造影剤が 以前力も用いられてきた。その目的は、造影剤を血中に投与して、血管系の構造や 分布のコントラストを高めた像を得て、悪性腫瘍や梗塞のような血管系に反映される 疾患を高い確度で診断することにある。
[0004] これに対し、超音波診断では、これまで造影剤が広く用いられることはな力つたが、 ここ数年、ミクロン 'オーダーのサイズをもつ微小気泡(マイクロバブル)を何らかの方 法で安定ィ匕した製剤による造影剤が出現したことにより、広く用いられはじめている。 この原理は、径 1ミクロン程度の微小気泡が、超音波診断に用いられる数 MHzの超 音波に共振して大振幅で振動し、その結果として、この周波数域の超音波をよく散乱 し、造影能が生ずることを利用するものである。
[0005] マイクロバブル系超音波造影剤の特徴は、その強 、非線形性にある。これは、マイ クロバブルには、負圧を受けた時の体積増加が、同振幅の正圧を受けた時の体積減 少よりはるかに大きくなる性質があることによる。このため、マイクロバブルにより散乱 されたエコー信号には、送信信号の 2倍の周波数をもつ第 2高調波成分が多く含ま れる。 V丄 .Newhouseらは、この第 2高調波成分をもとに軟部組織に対して強調された 血流ドッブラ信号を得る方式を、 1992年に初めて報告した (例えば、非特許文献 1を 参照)。
[0006] また、 P.N.Burnsらは、正負を反転させた送信音圧ノ ルス波形を用いて 2回の送受 信を行ない、得られた 2つのエコー信号をカ卩算するパルス反転(Pulse Inversion)法を 提案している(例えば、特許文献 1)。この加算により、動きの無視できる軟部組織か らのエコー信号の基本波成分は 180° 位相回転された信号が加算されるため打ち 消されてしまうが、第 2高調波成分は、 360° 位相回転されたものが加算されるため、 加算により 2倍に成長する。必要な送信回数が 2倍に増えるものの、原理的に帯域通 過フィルタなしに軟部組織力もの基本波成分を除くことができるので、距離分解能に 優れた第 2高調波エコー信号を得ることができる。また、血流中のマイクロバブル系 造影剤のように、 2回の送受信の間に生ずる変化が無視できな ヽ散乱体にっ 、ては 、それからの基本波エコー信号も完全には打ち消されないが、これは、軟部組織に 対して造影剤を強調したエコー像を得る今の目的に、むしろ合致している。
[0007] また、 W.Wilkeningは、送信音圧パルス波形の位相を 360° ZNずつ回転して N回 の送受信を行なう方法を提案している(例えば、非特許文献 2を参照)。例えば、 N = 3の場合は、 0° 、 120° 、 240° の 3回の送受信を行って得られたエコー信号を加 算する方法である。この方法を用いると、(N— 1)次高調波成分まで除去出来るとして いる。また、単純に加算するのでなぐ加算時にフィルタを通すことで、スペクトル特性 の異なる信号を峻別可能であるとして 、る。生体組織力ゝらの反射エコーと造影剤から の信号を峻別するには、以下のプロセスを経る。まず、予め造影剤エコーが支配的な 空間的な領域と、組織反射エコーが支配的な空間的領域が分かれたファントムに対 して超音波撮像を行ない、信号を取得する。次に、組織反射エコーと造影剤信号の 差が最大となるように、前述の加算時のフィルタ係数を最小自乗法によって決定する 。以下このフィルタを生体から信号に対しても適応することで、造影剤成分と生体反 射エコー成分の峻別が可能になるとしている。
[0008] 梅村は、 0° 、 120° 、 240° の 3回の送受信を行って得られたエコー信号を加算 し、造影剤信号と生体非線形信号を弁別する方法を報告して!/ヽる (非特許文献 3を 参照)。但し、非特許文献 2に記載の 3つのパルスを受信後にフィルタを使うか否かに 関しては記載がない。
[0009] 超音波送信波形に第 2高調波成分を意図的に重畳することにより、生体中又は液 体中におけるマイクロバブルの振動 ·成長 ·圧壊を強調したり、その反対に抑圧したり することができることが知られて!/ヽる (非特許文献 4を参照)。
[0010] 非特許文献 1: 1992 IEEE Ultrasonics Symposium Proceedings, pp.1175- 1177
非特許文献 2 : 2001 IEEE Ultrasonics Symposium Proceedings, pp.1733- 1737 非特許文献 3 : 2003 IEEE Ultrasonics Symposium Proceedings, pp.429- 432 非特許文献 4 : 1996 IEEE Transactions on Ultrasonics, Ferroelectrics, and
Frequency Control, vol.43, no.6, pp.1054— 1062
特許文献 1 :米国特許第 6, 095, 980号
発明の開示
発明が解決しょうとする課題
[0011] 上述したように、従来技術 (非特許文献 1)では、第 2高調波成分をもとに軟部組織 に対して強調された血流ドッブラ信号を得る方式を提案しているが、この第 2高調波 成分をエコー信号力も抽出するのに、帯域通過フィルタのみを用いたのでは、出力 信号として得られる第 2高調波エコー信号のパルスが長くなつてしまうという課題を残 している。
[0012] これは、エコー信号に含まれる基本波成分の振幅が第 2高調波成分よりも少なくと も 1桁以上大き ヽので、鋭!ヽ帯域遮断特性又は狭 、帯域通過特性をもつフィルタを 用いねばならないためである。この課題は、血流を 2次元画像表示しょうとする時、そ の距離方向の分解能を劣化させてしまうので、特に重大であった。
[0013] 従来技術 (特許文献 1)では、この課題を解決することを目的に、正負を反転させた 送信音圧パルス波形を用いて 2回の送受信を行な 、、得られた 2つのエコー信号を 加算するパルス反転法を提案して ヽる。
[0014] 生体軟部組織等多くの物質中の音速は、高い圧力の下では、低い圧力の下よりも 大きい。この非線形性に起因し、超音波パルスが軟部組織中を伝播する時、音圧の 高い部分が低い部分よりも早く伝播し、結果として、もともと正弦波状であった音圧波 形力 伝播するうちに、急に立ち上がりなだらかに下る N波状に変化すること、即ち、 第 2高調波等の高調波成分を持つようになることが知られて 、る。
[0015] これが、軟部組織中で散乱されると、マイクロバブル系造影剤が存在しなくても、高 調波成分を有するエコー信号が軟部組織力 戻ってくることになる。この高調波成分 をもとにエコー像を形成する方法は、 Tissue Harmonic Imaging法と呼ばれ、基本波 成分によるエコー像よりも音響 S/N比が高いことが好まれて、最近では、一般に用い られるようになった。し力し、このことは、パルス反転法を用いても、マイクロバブル系 造影剤により散乱されて生ずる高調波成分と送信パルスの伝播により生ずる高調波 成分とが混合したエコー信号が得られることになり、本来目的としていた造影剤を軟 部組織と峻別するイメージングが困難であることを意味する。
[0016] マイクロバブル系造影剤による非線形散乱は、軟部組織中の非線形伝播と比較す ると、一般に低い音圧でも観測される傾向にある。そこで、送信音圧を低く抑えたパ ルス反転法により、 Tissue Harmonic成分の生成を抑えて、主にマイクロバブル系造 影剤からの非線形成分によりエコー像を形成することが広く行われているが、信号振 幅が充分でな 、ため、造影診断に期待される確定診断を行なうに足る S/N比の高 ヽ エコー像が得にく 、のが現状である。
[0017] この課題をある程度解決する方法として、従来技術 (非特許文献 2)では、上述した ように、 360° ZNずつ位相回転したパルス波形で N回の送受信を行なう方法を提 案している。しかし、この方法は、ファントムで最適化されたフィルタであるから、被検 体毎に異なる信号に対して必ずしも適切でない。また、反射又は散乱源の深さによつ ても信号のスペクトルは変化すること、及び、非線形現象は音圧にも依存してスぺタト ルが変化することから、前述の最小自乗フィルタを全ての深さ、フォーカス条件に対 応して最適化するのは困難である。
[0018] 本発明は、このような現状に鑑み、マイクロバブル系造影剤により散乱されて生ずる エコー成分を、送信パルスが非線形伝播することにより生ずる Tissue Harmonic成分 と明確に峻別して映像ィ匕することにより、造影エコー像を基に確定診断を行うに足る S/N比の高 ヽ造影エコー像を実現する超音波撮像技術を提供することを目的とする 課題を解決するための手段
[0019] この課題を根本的に解決する手段として、本発明では、送信音圧パルス波形の包 絡線信号を共通とし、その搬送波の位相を 120° ずつ回転させて、 3回の送受信を 行ない、得られる 3つのエコー信号を加算する方法を提案する。この加算により、変 化の無視できる線形散乱体により散乱されたエコー信号の基本波成分と第 2高調波 成分とが、同時に打ち消される。これは、このような散乱体による 3つのエコー信号の 位相に着目する時、基本波成分の位相が 120° ずつ回転するのは当然として、第 2 高調波成分の位相も基本波成分とは逆方向に 120° ずつ回転するためである。
[0020] この原理は、レシプロ 4ストローク直列エンジンの振動問題を考えると理解しやすい 。クランクシャフトが一定の角速度で回転している時、レシプロエンジンを構成する各 ピストンが振動する線速度は、この角速度の基本波成分だけでなぐ無視できない振 幅の高調波成分を含む。 4ストローク直列 4気筒エンジンは、通常、対称配置された 2 つの同位相ピストン力もなる 2つの組が 180° のクランク角をなすように構成される。
[0021] 図 1は、パルス反転法の原理と直列 4気筒 4ストロークエンジンの振動の原理とを説 明する図である。
[0022] 図 1 (a)は、基本波の位相関係を示し、この構成により、それぞれの組のピストンが 発生する基本波成分は打ち消されるが、図 1 (b)は、第 2高調波の位相関係を示し、 第 2高調波成分は 2倍に成長することになり、結果として、クランクシャフト回転数の 2 倍の周波数をもつ振動が問題となる。図 1中、実線は、第 1のピストンの組による振動 の位相を、点線は、第 2のピストンの組による振動の位相を示す。クランクシャフトの 2 倍の角速度で回転するバランサーを装備する 4気筒エンジンがあるのは、この振動を 打ち消すためである。
[0023] 一方、 4ストローク直列 6気筒エンジンは、通常、対称配置された 2つの同位相ピスト ンカもなる 3つの組が 120° のクランク角をなすように構成される。
[0024] 図 2は、本発明の 3パルス法の原理と直列 6気筒 4ストロークエンジンの振動の原理 とを説明する図である。
[0025] 図 2 (a)は、基本波の位相関係を示し、この構成により、それぞれの組のピストンが 発生する基本波成分が、クランク角で数えて 120° をなすよう発生して互いに打ち消 しあう。図 2 (b)は、第 2高調波の位相関係を示し、第 2高調波成分も、図 2に示すよう に、クランク角で数えて 120° X 2 = 240° 、即ち、逆力も数えて 120° をなす位相 をもって発生するため、互いに打ち消される。これが、直列 6気筒エンジンが振動の 少ない所以である。図 2中、一点鎖線は、第 3のピストンの組による振動の位相を示 す。この第 2高調波振動が原理的に強調される直列 4気筒エンジンの構成が、パルス 反転法に相当し、基本波だけでなく第 2高調波の振動も原理的に打ち消される直列 6気筒エンジンの構成が、本発明の方法に相当する。
[0026] 一方、マイクロバブル系造影剤により散乱されて生ずるエコー信号の位相は、その 強 、非線形共振特性のために包絡線振幅の影響を受け、送波信号搬送波の位相 に対して、一定の関係にない。このため、送信パルス搬送波の位相を 120° ずつ回 転させて 3回の送受信を行って得られる 3つのエコー信号を加算しても、造影用マイ クロバブルにより散乱されて生じたエコー信号の場合には、打ち消されない成分が残 る。従って、この残ったエコー信号成分は、マイクロバブル系造影剤の存在のみを反 映するので、これを用いて、造影剤を軟部組織と明確に峻別した超音波イメージング が可能となる。
[0027] 非特許文献 2の公知例では、加算時にフィルタ処理が必要として!/、るが、本発明で は、それが不要になる理由は、上述の送波信号搬送波と造影剤エコーの位相関係 が一定の関係が崩れるような音圧の高い領域を用いることに加え、設計通りの送波 波形を送信可能なハードウェア構成を備えたことにある。具体的には DZAコンパ一 タに与える送波波形のサンプリング間隔を中心周波数の 3の倍数とする。このことによ り、包絡線共通で 120° ずつ位相を変えた 3つの超音波探触子から出る送波信号の 加算値を 0にすることが可能となる。
[0028] 本発明の超音波撮像装置は、造影用マイクロバブルが導入された生体に対して超 音波パルスの送受波を行なヽ、生体内部の画像を形成する超音波撮像装置であつ て、 Nを 3以上の整数とする時、同一の送受波フォーカス条件の下で波形の異なる送 信パルス波を用いて N回の送受波を行なうことにより、生体の軟部糸且織からの超音波 エコー信号の基本波から第 (N - 1)次高調波までの成分に対する送受波感度を抑圧 し、造影用マイクロバブル力 の超音波エコー信号に対する送受波感度を得るように 構成される。超音波撮像装置は、送信パルス波を生体内に送波するための送信用ァ ンプを備え、送信アンプに与える信号の入力周期が送波パルスの周波数成分のうち の最大の周波数に対して Nの整数倍である。また、送信用アンプに信号を与える手 段として DZAコンバータを備えており、 DZAコンバータの信号出力周期が送波パ ルスの周波数成分のうちの最大の周波数に対して Nの整数倍である。さらに、送信パ ルス波は、基本波とこの基本波に対する第二高調波を加算した波形である。
[0029] 本発明の超音波撮像装置は、造影用マイクロバブルが導入された生体に対して超 音波パルスの送受波を行な ヽ、生体内部の画像を形成する超音波撮像装置であつ て、同一の送受波フォーカス条件の下で波形の異なる送信ノ ルス波を用いて 3回の 送受波を行なうことにより、生体の軟部組織からの超音波エコー信号の基本波力 第 2次高調波までの成分に対する送受波感度を抑圧し、造影用マイクロバブルからの 超音波エコー信号に対して選択的に送受波感度を得る撮像シーケンスと、前述と同 一の送受波フォーカス条件の下で波形の異なる送信パルス波を用いて 2回の送受波 を行なうことにより、生体の軟部糸且織からの超音波エコー信号の基本波成分に対する 送受波感度を抑圧し、 2次以上の非線形信号に対して送受波感度を得る撮像シーケ ンスの二つのシーケンスを備え、二つのシーケンスを適宜切り替えて撮像を可能とし ている。 3回の送受波を行なうシーケンスと 2回の送受波を行なうシーケンスで、送波 振幅を異なるように構成する。また、 3回の送受波を行なうシーケンスでの送信振幅を 、 2回の送受波を行なうシーケンスでの送波振幅より大きくする。
[0030] 本発明の超音波撮像装置では、生体に対して超音波パルスの送受波を行な 、、生 体内部の造影用マイクロバブルによる造影画像を形成する力 同一の送受波フォー カス条件下で、包絡線信号を共通とする送信パルス波を用いその搬送波の位相を (a ) =0° 、(b) = 120° 、(c) =— 120° 、(d) = 180° の 4回の送受波を行ない、(a) 、(b)、(c)の 3個の時系列受信エコー信号を加算し造影画像を形成し、(a)、(d)の 2 個の時系列受信エコー信号を加算し生体伝播非線形画像を形成し、得られた 2種類 の画像を重畳して表示する。この結果、マイクロバブル系造影剤により散乱されて生 ずるエコー成分を、送信パルスの非線形伝播により生ずる Tissue Harmonic成分と明 確に峻別して映像ィ匕でき、 S/N比の高 ヽ造影エコー像を取得ができる。
発明の効果
[0031] 本発明によれば、エコー信号の中から、軟部組織中の非線形伝播等に由来する信 号成分を含まず、造影用マイクロバブルにのみ由来する信号を抽出し、これを造影 用マイクロバブルの空間分布を表す信号として、軟部組織の位置と形態を表す背景 上に、識別可能な色調で重畳して、 2次元又は 3次元画像によりわかりやすく表示で きる。これにより、造影エコー像を基に確定診断を行なうに足る SZN比の高い診断 用画像を提供することが可能となる。このように、本発明を実施した装置の医用診断 上の有用性はきわめて大きぐ従って、医用診断を支える工業における本発明の意 義も、また大きい。
発明を実施するための最良の形態
[0032] 以下、本発明の実施例について図を参照して説明する。なお、図 4、図 5、図 6、図
7、図 8、図 9、図 10、図 11、図 12、図 13、図 15、図 16【こお!/、て、縦軸 ίま音圧【こ it 例する信号振幅 (相対値で示す)を示し、横軸は時間( μ s)を示す。
[0033] 図 3は、本発明を実施すべく構成された超音波診断装置のブロック図の典型例で ある。
[0034] 超音波探触子 1を構成する各素子は、切り替えスィッチ群 2を介して、送波ビームフ ォーマ 3と受波ビームフォーマ 10に接続されている。送波ビームフォーマ 3では、送 受信シークェンス制御部 6による制御に従って、送信波形メモリ 5から送信波形選択 部 4により選択されて読み出された波形を用い、各素子を通じて送信された時に指向 性を持つ超音波パルスとなるような信号が生成される。この信号が、超音波探触子 1 の各素子により超音波ノ ルスに変換されて生体に送信される。生体中で反射あるい は散乱されて超音波探触子 1に戻ってきた超音波エコー信号は、各素子に受信され て、電気信号に変換される。
[0035] 受波ビームフォーマ 10では、送受信シークェンス制御部 6による制御に従って、指 向性を持つ受信感度を生成すベぐ各受波信号に遅延時間を与えて互いに加算す る。遅延加算により得られた時系列信号は、やはり送受信シークェンス制御部 6によ る制御に従って、受波メモリ選択部 11により選択された受信波形メモリ 12中の 1つの バンクへー且書き込まれ、互いに加算すべき時系列信号がそろったのちに読み出さ れて、加算器 13により互いに加算される。加算器 13の出力信号は、雑音成分を除去 する帯域通過フィルタ A14を通過したのち包絡線信号検出器 A16において包絡線 信号に変換され、スキャンコンバータ 18に入力される。
[0036] 一方、受信波形メモリ 12に書き込まれた時系列信号の一部は、読み出されたまま、 互いに加算されることなぐ雑音成分を除去する帯域通過フィルタ B15を通過し、包 絡線信号検出器 B 17にて包絡線信号に変換され、スキャンコンバータ 18に入力され る。スキャンコンバータ 18では、入力された複数の信号を適宜重畳して表示部 19〖こ て 2次元な 、し 3次元表示すべく、信号の生成 ·制御を行なう。
[0037] 図 15、図 16は、本発明の 3パルス法において、送波パルスの包絡線と各パルスの 搬送波の関係を説明する図である。図 16 (a)、図 16 (b)、図 16 (c)は、それぞれ、図 15 (b)、図 15 (c)、図 15 (d)に包絡線信号 (点線で示す)を重畳した図である。
[0038] 図 4は、本発明の 3パルス法により得られる、非線形伝播特性をもつ生体軟部組織 中の点反射体によるエコー信号の例を示す図である。
[0039] 送信波形メモリ 5に、図 15 (a)に示すように、共通の包絡線信号(Envelope)をもち、 図 15 (b)、図 15 (c)、図 15 (d)にそれぞれ示す、搬送波の位相が略 120° ずつ異な る 3つの超音波ノ ルス波形を書き込んでおき、送信波形選択部 4により、そのうち 1つ を選択して送受信するというシークェンスを、波形を代えて 3回実行した場合に、受 信波形メモリ 12中の各バンクに書き込まれる信号を、図 4 (a)、図 4 (b)、図 4 (c)に示 した。ここでは、簡単のために、送信超音波パルスが生体軟部組織中を伝播したの ち 1つの点反射体により反射されて生ずる受信エコー信号を、数値計算シミュレーシ ヨンにより求めて示した。搬送波周波数は、 2MHzとした。
[0040] 図 4 (d)は、図 4 (a)、図 4 (b)、図 4 (c)に示す信号が加算器 13に入力された時に 得られる出力信号である。送信超音波パルスが生体軟部組織中を非線形伝播する ため、図 4 (a)、図 4 (b)、図 4 (c)に示す信号には基本波成分だけでなく第 2高調波 成分も含まれているが、加算器出力結果図 4 (d)においては、本発明の原理から予 測されるように、基本波成分同士が打ち消し合うだけでなぐ第 2高調波成分同士も 打ち消し合い、信号振幅がほぼゼロとなっている。これと対比するため、図 5にパルス 反転法を用いた場合の結果を示した。
[0041] 図 5は、パルス反転法により得られる、非線形伝播特性をもつ生体軟部組織中の点 反射体によるエコー信号の例を示す図である。
[0042] この場合には、送信波形メモリ 5に、共通の包絡線信号をもち、搬送波の位相が 18 0° 異なる 2種類の超音波パルス波形を書き込んでおき、送信波形選択部 4により、 そのうち 1つを選択して送受信する。図 5 (a)、図 5 (b)には、このシークェンスを、波 形を代えて 2回実行した場合に、受信波形メモリ 12中の各バンクに書き込まれる信号 を示した。図 5 (c)は、その時得られる加算器 13の出力信号であるが、基本波成分同 士は打ち消し合っているが、第 2高調波成分同士はむしろ強調し合った結果の信号 が得られている。この信号は、 Tissue Harmonic信号と呼ばれるもので、生体軟部組 織を描出する時、高い音響 S/N比が得られるという利点はあるものの、造影剤の分布 や動態だけを軟部組織と峻別して描出したい場合には、それを妨害する最大要因と なってしまう。
[0043] 次に、図 4、図 5の場合と同じ送受信シークェンスにおいて、造影マイクロバブルに よる受信エコー信号を示す。
[0044] 図 6は、本発明の 3パルス法により得られる、造影用マイクロバブルによる散乱ェコ 一信号の例を示す図である。
[0045] 図 8は、パルス反転法により得られる、造影用マイクロバブルによる散乱エコー信号 の例を示す図である。
[0046] 図 7は、本発明の 3パルス法により得られる、造影用マイクロバブルによる散乱ェコ 一信号の帯域通過フィルタ通過後の波形を示す図である。
[0047] 図 9は、パルス反転法により得られる、造影用マイクロバブルによる散乱エコー信号 の帯域通過フィルタ通過後の波形を示す図である。
[0048] 図 6、図 8は、一例として、半径 1.5 μ mのマイクロバブルにより散乱されて生ずる受 信エコー信号を数値計算シミュレーションにより求めた結果であり、図 7、図 9は、その 時の帯域通過フィルタ A14の入出力信号である。図 6、図 7、図 8、図 9を通じ、音圧 に比例する縦軸の尺度は共通である。
[0049] それぞれの対応関係について述べれば、図 4の場合と同じく本発明を実施した時 の送受信シークェンスにおいて、受信波形メモリ 12中の各バンクに書き込まれる信 号が、図 6 (a)、図 6 (b)、図 6 (c)であり、その時得られる加算器 13の出力信号、即ち 、帯域通過フィルタ A14の入力信号が、図 6 (d)、及び、図 7中の(a)である。図 7 (b) 、図 7 (c)、図 7 (d)、図 7 (e)、図 7 (f)は、それぞれ、通過帯域中心周波数を、基本波 周波数(2MHz)、 1.5倍高調波周波数(3MHz)、第 2高調波周波数 (4MHz)、 2.5 倍高調波周波数(5MHz)、第 3高調波周波数 (6MHz)に設定した時の帯域通過フ ィルタ A14の出力信号である。
[0050] 一方、図 5の場合と同じくパルス反転法を実施した時の送受信シークェンスにおい て、受信波形メモリ 12中の各バンクに書き込まれる信号力 図 8 (a)、図 8 (b)であり、 その時得られる加算器 13の出力信号、即ち、帯域通過フィルタ A14の入力信号が、 図 8 (c)、及び、図 9 (a)である。図 9 (b)、図 9 (c)、図 9 (d)、図 9 (e)、図 9 (f)は、通過 帯域中心周波数を図 8の場合と同様に設定した時の帯域通過フィルタ A14の出力 信号である。
[0051] 元来、エコー信号中の第 2高調波成分を強調すべく考案されたパルス反転法により 、図 9 (c)、図 9 (d)に示されているように 1.5倍高調波力 第 2高調波成分を多く含む マイクロバブル由来の信号が得られるのは当然として、ェコ一信号中の第 2高調波成 分のうち非線形伝播等により生ずる成分を打ち消すべく考案された本発明の送受信 シークェンスにより、図 7 (d)、図 7 (e)、図 7 (f)に示されているように、第 2から第 3高 調波成分を多く含むマイクロバブル由来の充分な振幅の信号が得られていることは、 注目すべきである。
[0052] この特異かつ有用な現象の起源は、マイクロバブル力 大きな非線形性をもつ共振 体であること、より一般ィヒして言えば、遅延時間が振幅に依存する応答特性をもつこ とにある。即ち、入出力音圧の間に、非線形性があっても遅延応答時間に振幅依存 性がなければ、図 4 (d)に示されているように、出力信号中の第 2高調波成分は打ち 消されてしまう。一方、単なる線形共振体では、第 2高調波成分そのものが発生しな いので論外である。
[0053] 本発明による送受信シークェンスでは、以上の原理を背景に、送信パルス波形に 第 2高調波成分を意図的に重畳しても、それを打ち消しながら、造影用マイクロパブ ル由来の充分な振幅の信号が得られるという特長もある。従来技術 (非特許文献 4) に記載の第 2高調波重畳送波は、造影用マイクロバブルを用いた超音波イメージン グの場合にも、有用と考えられる。
[0054] 図 10は、本発明の 3パルス法において、送信パルス波に第 2高調波を意図的に重 畳した場合に得られる、生体軟部組織中の点反射体によるエコー信号の例を示す図 である。
[0055] 図 11は、本発明の 3パルス法において、送信パルス波に第 2高調波を意図的に重 畳した場合に得られる、造影用マイクロバブルによる散乱エコー信号の例を示す図 である。
[0056] 図 12は、本発明の 3パルス法において、送信パルス波に第 2高調波を意図的に重 畳した場合に得られる、造影用マイクロバブルによる散乱エコー信号の帯域通過フィ ルタ通過後の波形を示す図である。
[0057] 図 10、図 11、図 12に、送信パルス波形に第 2高調波成分を意図的に重畳した場 合に得られるエコー信号の例を示した。送信波形メモリ 5に、共通の包絡線信号をも ち、搬送波である基本波、及び、第 2高調波の位相が 120° ずつ異なる 3つの超音 波パルス波形を書き込んでおき、送信波形選択部 4により、そのうち 1つを選択して 送受信するというシークェンスを、波形を代えて 3回実行した場合に、受信波形メモリ 12中の各バンクに書き込まれる信号を、図 4の場合と同様に求めて、図 10 (a)、図 1 0 (b)、図 10 (c)に示し、その時得られる加算器 13の出力信号を図 10 (d)に示した。
[0058] また、造影用マイクロバブルにより散乱されて生ずる受信エコー信号を図 6の場合と 同様に求めて、受信波形メモリ 12中の各バンクに書き込まれる信号を図 11 (a)、図 1 1 (b)、図 11 (c)に示し、その時得られる加算器 13の出力信号、即ち、帯域通過フィ ルタ A14の入力信号を、図 11 (d)、及び、図 12 (a)に示した。図 12 (b)、図 12 (c)、 図 12 (d)、図 12 (e)、図 12 (f)は、通過帯域中心周波数を図 8の場合と同様に設定 した時の帯域通過フィルタ A14の出力信号である。
[0059] 図 10 (d)より明らかなように、送信パルス波形に第 2高調波成分を意図的に重畳し た場合でも、点反射体により散乱されて生ずる受信エコー信号については、加算器 1 3の入力信号の第 2高調波成分同士が、図 4の場合の非線形伝播により生ずる第 2 高調波成分と同様に、打ち消し合い、加算器の出力信号の振幅がほぼゼロとなって いる。一方、造影用マイクロバブルにより散乱されて生ずる受信エコー信号について は、図 6の場合と同様に、加算器 13による加算によって打ち消されることなぐ第 2か ら第 3高調波成分を多く含む充分な振幅の出力信号が得られている。
[0060] さらに、送信パルス波の位相について、本発明の効果を得るために必要な誤差範 囲について検討した。
[0061] 図 13は、本発明の 3パルス法において、送波パルスの位相に誤差があった場合に 得られる、生体軟部組織中の点反射体によるエコー信号の、加算器出力信号波形 の例を示す図である。
[0062] 例として、第 2送信パルス波の位相が 20° ずれた場合について、加算器 13の出力 信号、即ち、帯域通過フィルタ A14の入力信号と、通過帯域中心周波数を基本波と 第 2高調波に合わせた場合のフィルタ出力信号を、図 4の場合と同様に求めて、図 1 3(b),図 13(c)に示した。図 13 (a)には、比較のために、図 5の場合、即ち、パルス 反転法を位相誤差なしで実施できた場合の信号を示した。図 13(c)は、第 2送信パ ルス波の位相ずれ 20° に対応させて、第 3送信パルス波の位相を 10° ずらして、 3 つの送信パルス波の和信号がゼロになるように調整した場合である。
[0063] 送信ノ ルス波の包絡線信号を時間 tの関数として A (t)と書くと、本発明の実施例に おける第 1、第 2、第 3の各送波パルス信号 PI (t)、 P2 (t)、 P3 (t)は、位相誤差のな い時、(数 1)、(数 2)、(数 3)と書くことができる。この時、(数 4)の関係が成り立って いる。
[0064] Pl(t)=A(t)sincot …(数 1)
P2 (t) = A (t) sin ( ω t + 2 π /3) …(数 2)
Ρ3 (t) = A (t) sin ( ω t-2 π /3) …(数 3)
Pl(t)+P2(t)+P3(t)=0 …(数 4)
第 2パルスに位相誤差 φを生じた時、第 2パルスは、(数 5)と書くことができる。
[0065] この時、第 3パルスを (数 6)、(数 7)で書けるよう調整すれば、位相誤差 φにかかわ らず、(数 4)を成り立たせることができる。図 13(c)は、そのような第 3パルスの補正を 行った結果であり、第 2パルスの位相誤差にもかかわらず、加算器 13の出力信号中 の基本波成分を打ち消すことができている。(数 4)が成立していること、一般的には、 送信に用いる N個のパルス波形の和信号が実質的にゼロであることにより、 N個の受 信エコー信号を加算して得られる信号中の基本波成分が打ち消される。
[0066]
Figure imgf000015_0001
φ) …(数 5)
P3 (t) = A3 (t) sin (ωί~2π/3+ /2) · · . (数 6) A3 (t) = 2A (t) cos ( π /3 + /2) …(数 7)
図 14は、本発明の 3パルス法において、送波パルスの位相に誤差があった場合に 得られる、生体軟部組織中の点反射体によるエコー信号の、加算器の出力信号の振 幅の位相誤差依存性を説明する図である。
[0067] 図 14は、加算器 13の出力信号の振幅の peak-to-peak値(出力信号相対振幅値 PP 値)(縦軸)を、第 2パルスに与えた位相誤差 (度)(横軸)の関数として示している。パ ルス反転法を位相誤差なしで実施できた場合の値で規格化した信号振幅を、第 3パ ルスの補正なしの場合(図 14中の(a) )と、補正を行った場合(図 14中の(b) )につい て示した。位相誤差が 20° に達すると、第 3パルスの補正を行っても、軟部組織中の 非線形伝播に由来し、造影用マイクロバブルに由来しない信号振幅力 従来のパル ス反転法の半分以上となり、本発明の効果が充分に発揮されないことになつてしまう 。この結果から、本発明の効果を充分に得るには、送信パルスの位相誤差を 10° 程 度以下とすることが望ま 、と 、える。
[0068] 以上説明したように、本発明を実施することにより、エコー信号の中から、軟部組織 中の非線形伝播等に由来する信号成分を含まず、造影用マイクロバブルにのみ由 来する信号を抽出することができる。
[0069] そのような信号を加算器 13の出力信号として得ることができ、さらに S/N比を向上さ れた信号を帯域通過フィルタ A14の出力信号として得て、その包絡線信号として包 絡線信号検出器 A16の出力信号を得て、これが造影用マイクロバブルの空間分布 を表す信号としてスキャンコンバータ 18に入力される。
[0070] 一方、受信波形メモリ 12中の 1つのバンクに書き込まれた信号を帯域通過フィルタ B15に通して S/N比を向上させた信号から、包絡線信号検出器 B17により包絡線信 号を得て、これが軟部組織の位置と形態を表す信号としてスキャンコンバータ 18に 入力される。
[0071] スキャンコンバータ 18では、包絡線信号検出器 B17の出力信号を背景とし、これと 識別するのに便利なように異なる色調によって包絡線信号検出器 A16の出力信号 が重畳されて表示部 19に表示される。このようにして、検査対象である患者の体の中 における造影用マイクロバブルの分布を、 2次元又は 3次元画像によりわかりやすく 表示できる。
[0072] 上述した実施例では、軟部組織の位置と形態を表す信号として、受信波形メモリ 12 中の 1つのバンクに書き込まれた信号をそのまま用いる場合について説明した力 一 般的には、受信波形メモリ 12中の 3つのバンクに書き込まれた信号に適切な重みを つけた加算信号を用いることができる。
[0073] また、以上の一連の説明では、 3つの送信ノ ルス波を用いる実施例について詳しく 述べたが、原理説明において示したように、共通の包絡線信号をもつ 3以上の整数 N個の送信パルス波を用い、その搬送波の位相を 360° ZNずつ回転させて、 N回 の送受信を行ない、得られる N個のエコー信号を受信波形メモリ 12中の N個のバン クに一旦書き込み、読み出した信号を加算器 13に入力することによって、本発明を 実施することちでさる。
[0074] 次に、送波波形の与え方について、図 17、 18を用いて説明する。
[0075] 図 17は、送波アンプの入出力特性の一例を示す図である。
[0076] 図 18は、 DZAコンバータの出力波形とサンプリング点の関係を説明する図である
[0077] 送波波形の出力方法として、公知の方法としては、スィッチで正負の信号を与える 方式と、 DZAコンバータと送波アンプを組み合わせる方法がある。前者は、 0° と 18 0° の信号にしか対応しないので、 3以上の位相の異なる複数波形を送波するのに 適さない。一方、 DZAコンバータを使った場合も、そのまま複数波形送波が可能に なるというわけではない。図 17は一般的な送波アンプの入出力特性を示した図であ る。横軸は入力電圧、縦軸は出力電圧である。本来入力に対して出力は比例してい ることが望ましいが、実際には図 17に示すように非線形性がある。このアンプの非線 形特性は 3パルス送波に大きな影響を与える。
[0078] 図 18は、 0° 、 120° 、 一 120° の 3つのパルスの DZAコンバータの出力波形を、 それぞれ、実線、破線、点線で示している。丸印の箇所が DZAコンバータのサンプ リング点を示す。図 18において、縦軸は DZAコンバータの出力値 (相対値)を示し、 縦軸は、 DZAコンバータに制御用の出力を与えるタイミングを示す無次元数を示す 。図 18 (a)は中心周波数の 4倍の周波数でサンプリングした場合、図 18 (b)は中心 周波数の 6倍の周波数でサンプリングした場合、図 18 (c)は中心周波数の 8倍の周 波数でサンプリングした場合である。 4倍、 8倍の場合、ノ ルス毎に DZAの出力値が 異なる値をとる。これらの場合、図 17に示す送波アンプの非線形性のため、設定した とおりの出力値とならない。
[0079] このため、全く非線形性の無い媒質からの受信信号すら、加算後に打ち消し残りが 生じてしまうので、目的を実現することが出来なくなる。これが非特許文献 2の 5ベー ジ目二番目の段落に記載してある、 3パルスの非対称性の問題である。本発明では 図 18 (b)に示すように中心周波数の 6倍 (Nパルスの場合は Nの整数倍)のサンプリ ング周波数で DZAコンバータから出力する。その結果 DZAの出力値がパルス毎 に大きく異ならないため、送波アンプに非線形性があっても問題とならない。
[0080] 次に、造影剤の破壊に対する考慮した撮像法に関して説明する。
[0081] 図 19は、造影剤信号と生体非線形信号の弁別比と、造影剤信号の感度の、送波 振幅依存性を説明する図である。図 19 (a)は、パルス反転法と 3パルス法の造影剤 信号と生体非線形信号の弁別比を送波信号の振幅に対してプロットした図であり、図 19 (b)は、造影剤信号の感度を送波信号の振幅に対してプロットした図である。図 1 9において、横軸は、音圧振幅 X O. IMPaで送波振幅を示し、縦軸は、(a)弁別比( Sensitivity, dB)、 (b)感度(Sensitivity, dB)を示し、実線は 3パルス法、点線はパルス 反転 (インバージョン)法を示す。弁別比に関しては、 3パルス法の方が従来のパルス 反転法に比べ、常に 20dB力も 40dB程度優れて 、る。
[0082] 一方、感度に関しては、送波振幅が小さくなると、 3パルス法の感度が低下してくる 。これは、大きな送波振幅、即ち、非線形性が大きい時にし力、前述の造影剤からの 信号と送波信号の位相関係が崩れないからである。しかし、感度を得るために大きな 送波振幅で撮像すると、撮像している最中に造影剤が破壊され、結局感度の良い状 態を持続することが困難となる。
[0083] よって、本発明では、以下の二通りの方法でこの問題を回避する。一つは図 10に 示したような、第二高調波を重畳した波形を用いることによる。例えば、造影剤は押さ れた時より、引っ張られた場合に破壊し易いので、正圧に比べ負圧が小さくなるよう に、第二高調波重畳した波形を用いることで、音圧が大きくでも比較的造影剤が壊 れに《なるように工夫することが出来る。
[0084] また、別の方法としては、シーケンス制御による方法も有効である。これはパルス反 転送波と 3パルス法を交互に送受信を行ない。普段は弱い音圧で、パルス反転法で モニタリングを行な 、、時々高音圧の 3パルス法を用いる方法である。
[0085] 図 20、図 21は、パルス反転法と 3パルス法の切り替え撮像シーケンスの制御フロー 図である。
[0086] 図 20に示す制御処理フローのように、予め設定した M回と N回で、切り替える以外 に、オペレータからの入力に応じても切り替えを可能としている。パルスインバージョ ン (反転)撮像処理 20の N回の撮像終了又は外部入力により切り替えの判定処理 21 を行ない、 3パルス撮像処理 22を行なう、 3パルス撮像処理 22の M回の撮像終了又 は外部入力により切り替えの判定処理 23を行なう。
[0087] 図 20はフレーム毎にパルス反転 (インバージョン)と 3パルス法を切り替える方法で あるが、図 21に示す方法では、ラスター毎に切り替える。図 21に示すように、 (a) 18 0° パルス送受波処理 101、 (b) 0° パルス送受波処理 102、 (c) 120° パルス送受 波処理 103、(d)— 120° パルス送受波処理 104を行ない、 1フレーム終了か否かの 判定処理 105を行ない、 1フレーム終了後にラスター移動を行なう。この方法によって 、被検体の動きがある場合でも、切り替え時間が充分短いため、パルス反転法と 3パ ルス法で撮像した信号の対象が同等のものと見なすことが出来るようになる。そのた め、パルス反転法で撮像した白黒像の上に、 3パルス法による造影剤信号を別のカラ 一コードで重畳することが可能となる。 産業上の利用可能性
[0088] 本発明によれば、造影エコー像を基に確定診断を行なうに足る SZN比の高い超 音波撮像装置を提供できる。
図面の簡単な説明
[0089] [図 1]パルス反転法の原理と直列 4気筒 4ストロークエンジンの振動の原理とを説明す る図である。
[図 2]本発明の 3パルス法の原理と直列 6気筒 4ストロークエンジンの振動の原理とを 説明する図である。 圆 3]本発明の実施例の超音波診断装置の構成を示すブロック図である。
圆 4]本発明の 3パルス法により得られる、非線形伝播特性をもつ生体軟部組織中の 点反射体によるエコー信号の例を示す図である。
[図 5]パルス反転法により得られる、非線形伝播特性をもつ生体軟部組織中の点反 射体によるエコー信号の例を示す図である。
[図 6]本発明の 3パルス法により得られる、造影用マイクロバブルによる散乱エコー信 号の例を示す図である。
[図 7]本発明の 3パルス法により得られる、造影用マイクロバブルによる散乱エコー信 号の帯域通過フィルタ通過後の波形を示す図である。
[図 8]パルス反転法により得られる、造影用マイクロバブルによる散乱エコー信号の例 を示す図である。
[図 9]パルス反転法により得られる、造影用マイクロバブルによる散乱エコー信号の帯 域通過フィルタ通過後の波形を示す図である。
[図 10]本発明の 3パルス法において、送信パルス波に第 2高調波を意図的に重畳し た場合に得られる、生体軟部組織中の点反射体によるエコー信号の例を示す図であ る。
圆 11]本発明の 3パルス法において、送信パルス波に第 2高調波を意図的に重畳し た場合に得られる、造影用マイクロバブルによる散乱エコー信号の例を示す図である
[図 12]本発明の 3パルス法において、送信パルス波に第 2高調波を意図的に重畳し た場合に得られる、造影用マイクロバブルによる散乱エコー信号の帯域通過フィルタ 通過後の波形を示す図である。
[図 13]本発明の 3パルス法において、送波パルスの位相に誤差があった場合に得ら れる、生体軟部組織中の点反射体によるエコー信号の、加算器出力信号波形の例 を示す図である。
[図 14]本発明の 3パルス法において、送波パルスの位相に誤差があった場合に得ら れる、生体軟部組織中の点反射体によるエコー信号の、加算器の出力信号の振幅 の位相誤差依存性を説明する図である。 [図 15]本発明の 3パルス法において、送波パルスの包絡線と各パルスの搬送波の関 係を説明する図である。
[図 16]本発明の 3パルス法において、送波パルスの包絡線と各パルスの搬送波の関 係を説明する図である。
[図 17]送波アンプの入出力特性の一例を示す図である。
[図 18]DZAコンバータの出力波形とサンプリング点の関係を説明する図である。
[図 19]造影剤信号と生体非線形信号の弁別比と、造影剤信号の感度の送波振幅依 存性を説明する図である。
[図 20]ノ ルス反転法と 3パルス法の切り替え撮像シーケンスの制御フロー図である。
[図 21]パルス反転法と 3パルス法の切り替え撮像シーケンスの制御フロー図である。 符号の説明
1…超音波探触子、 2…切り替えスィッチ群、 3…送波ビームフォーマ、 4…送信波 形選択部、 5…送信波形メモリ、 6…送受信シークェンス制御部、 10…受波ビームフ ォーマ、 11···受波メモリ選択部、 12···受信波形メモリ、 13···加算器、 14···帯域通過 フィルタ A、 15···帯域通過フィルタ B、 16···包絡線信号検出器 A、 17···包絡線信号 検出器 B、 18…スキャンコンバータ、 19···表示器、 20···パルスインバージョン撮像処 ¾、 21···半 lj 、 22"'3ノ ノレス 、 23···半 lj 、 101···180ο ノ ノレス 送受波処理、 102···0° パルス送受波処理、 103···120° パルス送受波処理、 104 …一 120° パルス送受波処理、 105···判定処理。

Claims

請求の範囲
[1] 造影用マイクロバブルが導入された生体に対して超音波パルスの送受波を行な ヽ 、前記生体内部の画像を形成する超音波撮像装置であって、 Nを 3以上の整数とす る時、同一の送受波フォーカス条件の下で波形の異なる送信パルス波を用いて N回 の送受波を行なうことにより、前記生体の軟部組織力 の超音波エコー信号の基本 波から第 (N - 1)次高調波までの成分に対する送受波感度を抑圧し、前記造影用マ イクロバブルからの超音波エコー信号に対する送受波感度を得るように構成され、前 記送信パルス波を生体内に送波するための送信用アンプを備え、前記送信アンプに 与える信号の入力周期が送波パルスの周波数成分のうちの最大の周波数に対して Nの整数倍であることを特徴とする超音波撮像装置。
[2] 請求項 1に記載の超音波診断装置において、前記送信用アンプに信号を与える手 段として DZAコンバータを備え、前記 DZAコンバータの信号出力周期が前記送波 パルスの周波数成分のうちの最大の周波数に対して Nの整数倍であることを特徴と する超音波撮像装置。
[3] 造影用マイクロバブルが導入された生体に対して超音波パルスの送受波を行な ヽ 、前記生体内部の画像を形成する超音波撮像装置であって、 Nを 3以上の整数とす る時、同一の送受波フォーカス条件の下で波形の異なる送信パルス波を用いて N回 の送受波を行なうことにより、前記生体の軟部組織力 の超音波エコー信号の基本 波から第 (N - 1)次高調波までの成分に対する送受波感度を抑圧し、前記造影用マ イクロバブルからの超音波エコー信号に対する送受波感度を得るように構成され、前 記送信パルス波は、基本波とこの基本波に対する第二高調波を加算した波形である ことを特徴とする超音波撮像装置。
[4] 造影用マイクロバブルが導入された生体に対して超音波パルスの送受波を行な ヽ 、前記生体内部の画像を形成する超音波撮像装置であって、同一の送受波フォー カス条件の下で波形の異なる送信パルス波を用いて 3回の送受波を行なうことにより 、前記生体の軟部組織力 の超音波エコー信号の基本波から第 2次高調波までの 成分に対する送受波感度を抑圧し、前記造影用マイクロバブルからの超音波エコー 信号に対して選択的に送受波感度を得る撮像シーケンスと、前記同一の送受波フォ 一カス条件の下で波形の異なる送信パルス波を用いて 2回の送受波を行なうことによ り、前記生体の軟部組織からの超音波エコー信号の基本波成分に対する送受波感 度を抑圧し、 2次以上の非線形信号に対して送受波感度を得る撮像シーケンスの二 つのシーケンスを備え、前記二つのシーケンスを適宜切り替えて撮像可能であること を特徴とする超音波撮像装置。
[5] 請求項 4に記載の超音波撮像装置において、前記 3回の送受波を行なうシーケン スと前記 2回の送受波を行なうシーケンスで、送波振幅を異なるように構成したことを 特徴とする超音波撮像装置。
[6] 請求項 5に記載の超音波撮像装置において、前記 3回の送受波を行なうシーケン スでの送信振幅が、前記 2回の送受波を行なうシーケンスでの送波振幅より大きいこ とを特徴とする超音波撮像装置。
PCT/JP2005/000047 2004-03-12 2005-01-06 超音波撮像装置 WO2005087109A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP05703319.3A EP1726261B1 (en) 2004-03-12 2005-01-06 Ultrasonographic device
US10/586,665 US8002703B2 (en) 2004-03-12 2005-01-06 Ultrasonic imaging device
JP2006510881A JP4415011B2 (ja) 2004-03-12 2005-01-06 超音波撮像装置

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2004-069984 2004-03-12
JP2004069984 2004-03-12

Publications (1)

Publication Number Publication Date
WO2005087109A1 true WO2005087109A1 (ja) 2005-09-22

Family

ID=34975275

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2005/000047 WO2005087109A1 (ja) 2004-03-12 2005-01-06 超音波撮像装置

Country Status (5)

Country Link
US (1) US8002703B2 (ja)
EP (1) EP1726261B1 (ja)
JP (1) JP4415011B2 (ja)
CN (1) CN100502790C (ja)
WO (1) WO2005087109A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011235018A (ja) * 2010-05-13 2011-11-24 Hitachi Aloka Medical Ltd 超音波診断装置
GB2488390A (en) * 2010-12-13 2012-08-29 Bosch Gmbh Robert Object detection
JP2014180363A (ja) * 2013-03-19 2014-09-29 Konica Minolta Inc 超音波画像診断装置
WO2023132088A1 (ja) 2022-01-04 2023-07-13 ソニア・セラピューティクス株式会社 超音波画像処理装置

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009022462A (ja) * 2007-07-18 2009-02-05 Toshiba Corp 超音波診断装置および超音波診断装置の制御プログラム
US20100256498A1 (en) * 2007-11-16 2010-10-07 Hiroki Tanaka Ultrasonic imaging device
JP5735496B2 (ja) * 2010-06-04 2015-06-17 株式会社日立メディコ 超音波診断装置
WO2014009442A1 (en) 2012-07-11 2014-01-16 Sanofi-Aventis Deutschland Gmbh Arrangement and method for determining a stopper position
US9274215B2 (en) * 2013-03-08 2016-03-01 Chison Medical Imaging, Inc. Ultrasound fusion harmonic imaging systems and methods
CN104224232A (zh) * 2014-09-25 2014-12-24 飞依诺科技(苏州)有限公司 超声谐波成像方法及装置
GB201513024D0 (en) * 2015-07-23 2015-09-09 Univ Heriot Watt A method of, and apparatus for, determination of position in ultrasound imaging
US10804942B2 (en) * 2018-05-24 2020-10-13 Analog Devices, Inc. State-machine based body scanner imaging system
CN113117261B (zh) * 2019-12-30 2023-06-02 重庆融海超声医学工程研究中心有限公司 用于检测空化效应的装置及超声治疗设备
CN117918887A (zh) * 2024-03-20 2024-04-26 中国医学科学院北京协和医院 基于脉冲反转二次谐波成像的压力确定方法及装置

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08628A (ja) * 1994-04-21 1996-01-09 Hitachi Medical Corp 超音波断層装置
JP2001061841A (ja) * 1999-08-30 2001-03-13 Toshiba Corp 超音波診断装置及び超音波画像の生成方法
JP2003038490A (ja) * 2001-05-31 2003-02-12 Novasonics Inc 位相反転超音波画像処理システムおよび方法

Family Cites Families (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5517996A (en) * 1994-04-21 1996-05-21 Hitachi Medical Corporation Ultrasonic diagnostic apparatus
US5833614A (en) * 1997-07-15 1998-11-10 Acuson Corporation Ultrasonic imaging method and apparatus for generating pulse width modulated waveforms with reduced harmonic response
US6193659B1 (en) * 1997-07-15 2001-02-27 Acuson Corporation Medical ultrasonic diagnostic imaging method and apparatus
US5944666A (en) * 1997-08-21 1999-08-31 Acuson Corporation Ultrasonic method for imaging blood flow including disruption or activation of contrast agent
US6095980A (en) 1997-10-02 2000-08-01 Sunnybrook Health Science Centre Pulse inversion doppler ultrasonic diagnostic imaging
US5902243A (en) * 1998-04-15 1999-05-11 Acuson Corporation Ultrasonic imaging method with multiple pulse cancellation
US5984869A (en) * 1998-04-20 1999-11-16 General Electric Company Method and apparatus for ultrasonic beamforming using golay-coded excitation
JP2001212144A (ja) 2000-01-31 2001-08-07 Toshiba Corp 超音波診断装置及び超音波画像化方法
US6497665B1 (en) * 2000-07-14 2002-12-24 Koninklijke Philips Electronics N.V. System and method for non-linear detection of ultrasonic contrast agents at a fundamental frequency
JP4232134B2 (ja) 2001-02-01 2009-03-04 株式会社日立メディコ 超音波造影描画装置
JP4744727B2 (ja) 2001-06-06 2011-08-10 株式会社東芝 超音波画像診断装置
JP3959257B2 (ja) 2001-11-08 2007-08-15 株式会社東芝 超音波診断装置
US7094204B2 (en) * 2002-08-23 2006-08-22 Siemens Medical Solutions Usa, Inc. Coded excitation imaging for use with bipolar, unipolar and other waveforms

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08628A (ja) * 1994-04-21 1996-01-09 Hitachi Medical Corp 超音波断層装置
JP2001061841A (ja) * 1999-08-30 2001-03-13 Toshiba Corp 超音波診断装置及び超音波画像の生成方法
JP2003038490A (ja) * 2001-05-31 2003-02-12 Novasonics Inc 位相反転超音波画像処理システムおよび方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
See also references of EP1726261A4 *
UMEMURA S.I. ET AL: "Enhancement of Sonodynamic Tissue Damage Production by Second-Harmonic Superimposition: Theoretical Analysis of Its Mechanism", IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS AND FREQUENCY CONTROL, vol. 43, no. 6, 1996, pages 1054 - 1062, XP000636974 *
UMEMURA S.I. ET AL: "Triplet Pulse sequence for superior microbubble/tissue contrast", 2003 IEEE ULTRASONICS SYMPOSIUM PROCEEDINGS, October 2003 (2003-10-01), pages 429 - 432, XP010702256 *
WILKENING W. ET AL: "Phase-Coded Pulse Sequence for Non-Linear Imaging", 2000 IEEE ULTRASONICS SYMPOSIUM, PROCEEDINGS, October 2000 (2000-10-01), pages 1559 - 1562, XP010540912 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011235018A (ja) * 2010-05-13 2011-11-24 Hitachi Aloka Medical Ltd 超音波診断装置
GB2488390A (en) * 2010-12-13 2012-08-29 Bosch Gmbh Robert Object detection
JP2014180363A (ja) * 2013-03-19 2014-09-29 Konica Minolta Inc 超音波画像診断装置
WO2023132088A1 (ja) 2022-01-04 2023-07-13 ソニア・セラピューティクス株式会社 超音波画像処理装置

Also Published As

Publication number Publication date
EP1726261B1 (en) 2014-07-16
US20080228076A1 (en) 2008-09-18
CN100502790C (zh) 2009-06-24
JP4415011B2 (ja) 2010-02-17
CN1909837A (zh) 2007-02-07
JPWO2005087109A1 (ja) 2008-01-24
EP1726261A4 (en) 2010-05-26
US8002703B2 (en) 2011-08-23
EP1726261A1 (en) 2006-11-29

Similar Documents

Publication Publication Date Title
WO2005087109A1 (ja) 超音波撮像装置
JP4365909B2 (ja) パルス反転ドップラー超音波診断画像処理方法及び装置
US6171246B1 (en) Realtime ultrasonic imaging of perfusion using ultrasonic contrast agents
JP4667733B2 (ja) ステップ状チャープ波形を使用して超音波造影イメージングを増強する方法及び装置
JP4567842B2 (ja) 超音波診断装置
JP4422421B2 (ja) 超音波撮像装置
US20040039282A1 (en) System and method for improved harmonic imaging
EP2154548A1 (en) Ultrasound diagnostic apparatus for nonlinear contrast imaging
US20080249417A1 (en) Non-Linear Ultrasonic Diagnostic Imaging Using Intermodulation Product Signals
US20080275338A1 (en) Nonlinear Ultrasonic Diagnostic Imaging Using Intermodulation Product Signals
JP2007090075A (ja) コントラスト剤を超音波によって特徴付ける方法およびコントラスト剤を超音波によって特徴付けるためにプログラムされたプロセッサによって実行可能な命令を表すデータが記憶されたコンピュータ読み出し可能記憶媒体
JP4717484B2 (ja) 超音波診断装置
KR20070038471A (ko) 적절한 mi 레벨에서 초음파 진단 콘트래스트 영상형성
JPH02307436A (ja) 超音波血流イメージング装置
JP4690537B2 (ja) 超音波診断装置
JP4157688B2 (ja) 超音波診断装置
JPH10118065A (ja) 超音波診断装置
Shen et al. Optimal transmit phasing on tissue background suppression in contrast harmonic imaging
Løvstakken Signal processing in diagnostic ultrasound: Algorithms for real-time estimation and visualization of blood flow velocity
JP2009233352A (ja) 超音波撮像装置
KR20230062418A (ko) 매체의 물리적 특성을 결정하기 위한 방법 및 시스템
Kruse et al. A new imaging strategy utilizing wideband transient response of ultrasound contrast agents
CN117918887A (zh) 基于脉冲反转二次谐波成像的压力确定方法及装置
Ma et al. Investigation on phase-coded third harmonic imaging for normal and pathological tissues in transmission mode in vitro
Goertz et al. 1J-3 High Frequency Colour Flow Imaging of Mouse Tumours using Inter-Frame Clutter Filtering

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 SM 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): 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

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: 2006510881

Country of ref document: JP

WWE Wipo information: entry into national phase

Ref document number: 10586665

Country of ref document: US

Ref document number: 2005703319

Country of ref document: EP

Ref document number: 200580002796.3

Country of ref document: CN

NENP Non-entry into the national phase

Ref country code: DE

WWW Wipo information: withdrawn in national office

Ref document number: DE

WWP Wipo information: published in national office

Ref document number: 2005703319

Country of ref document: EP