JP5043623B2 - Ultrasonic diagnostic equipment - Google Patents

Ultrasonic diagnostic equipment Download PDF

Info

Publication number
JP5043623B2
JP5043623B2 JP2007318814A JP2007318814A JP5043623B2 JP 5043623 B2 JP5043623 B2 JP 5043623B2 JP 2007318814 A JP2007318814 A JP 2007318814A JP 2007318814 A JP2007318814 A JP 2007318814A JP 5043623 B2 JP5043623 B2 JP 5043623B2
Authority
JP
Japan
Prior art keywords
bubble
transmission
ultrasonic
signal
frequency
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
JP2007318814A
Other languages
Japanese (ja)
Other versions
JP2009136626A (en
Inventor
和徳 射谷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Aloka Medical Ltd
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 Aloka Medical Ltd filed Critical Hitachi Aloka Medical Ltd
Priority to JP2007318814A priority Critical patent/JP5043623B2/en
Publication of JP2009136626A publication Critical patent/JP2009136626A/en
Application granted granted Critical
Publication of JP5043623B2 publication Critical patent/JP5043623B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Description

本発明は、超音波診断装置に関し、特に、バブルを利用した画像形成技術に関する。   The present invention relates to an ultrasonic diagnostic apparatus, and more particularly to an image forming technique using bubbles.

マイクロバブル(またはナノバブル)は、液体などに注入された微細な気泡を意味している。このマイクロバブルは、様々な優れた特性を備えているため、多くの分野で利用されている。例えば、医療分野への応用として、マイクロバブルが超音波の好適な反射体となることから、超音波画像を形成する際の造影剤として利用されている。   Microbubbles (or nanobubbles) mean fine bubbles injected into a liquid or the like. Since these microbubbles have various excellent characteristics, they are used in many fields. For example, as an application to the medical field, microbubbles are suitable as ultrasonic wave reflectors, and are therefore used as a contrast agent when forming an ultrasonic image.

造影剤を利用した超音波画像の形成においては、バブルから得られる高調波成分に注目した技術が知られており、特に、フェーズインバージョン法やパワーモジュレーション法などが有名である。もちろん、高調波成分に注目した画像形成技術に関する文献も、例えば、特許文献1,2などのように多数存在する。ちなみに、特許文献1,2に記載された技術は、バブルから得られる第3次高調波を積極的に画像化するものではない。   In the formation of an ultrasonic image using a contrast agent, a technique that pays attention to harmonic components obtained from bubbles is known, and in particular, a phase inversion method, a power modulation method, and the like are well known. Of course, there are a large number of documents relating to image forming techniques that focus on harmonic components, such as Patent Documents 1 and 2, for example. Incidentally, the techniques described in Patent Documents 1 and 2 do not actively image third-order harmonics obtained from bubbles.

比較的太い血管などに投与されたマイクロバブル(造影剤)を画像化する場合には、マイクロバブルの周囲が散乱などの少ない血液のため、従来から知られている技術により、比較的高いコントラストを得ることができる。これに対し、組織や腫瘍などの毛細血管やリンパ節などに取り込まれたマイクロバブルを観測する場合には、マイクロバブルの周囲が組織や腫瘍などの実組織であるため、周囲が血液である場合に比べて、高いコントラストを得ることが難しい。   When imaging microbubbles (contrast media) administered to relatively thick blood vessels, the surroundings of the microbubbles are blood with little scattering, so a relatively high contrast can be achieved using a conventionally known technique. Obtainable. On the other hand, when observing microbubbles taken into capillaries or lymph nodes such as tissues and tumors, the surroundings are blood because the surroundings of the microbubbles are real tissues such as tissues and tumors. Compared to, it is difficult to obtain high contrast.

また、近年になり、低中音圧の超音波による画像化に適した造影剤(例えば「ソナゾイド」(登録商標))が注目されており、低中音圧の超音波による高感度の画像化技術も望まれている。   In recent years, contrast agents suitable for imaging with low and medium sound pressure ultrasound (for example, “Sonazoid” (registered trademark)) have been attracting attention, and high sensitivity imaging with low and medium sound pressure ultrasound has been attracting attention. Technology is also desired.

特開2005−152177号公報JP 2005-152177 A 特開2001−258882号公報JP 2001-258882 A

このような状況のもと、本願の発明者は、バブル(マイクロバブルなど)を利用した超音波による画像化技術について研究を重ねてきた。特に、超音波を受けたバブルの挙動に注目して研究を重ねてきた。   Under such circumstances, the inventors of the present application have conducted research on imaging technology using ultrasonic waves using bubbles (such as microbubbles). In particular, research has been conducted focusing on the behavior of bubbles subjected to ultrasonic waves.

本発明は、その研究の過程において成されたものであり、その目的は、バブルの挙動に基づいた適切な画像処理技術を提供することにある。   The present invention has been made in the course of its research, and an object thereof is to provide an appropriate image processing technique based on the behavior of bubbles.

上記目的を達成するために、本発明の好適な態様である超音波診断装置は、バブルを含む診断領域に対して超音波を送受するプローブと、前記バブルの共振周波数を基準として設定される送信周波数かつ前記バブルの膨張倍率を基準として設定される送信音圧で超音波を送波するようにプローブを制御する送信制御部と、プローブにおける超音波の受波により前記診断領域から得られる受信信号を処理することにより受信信号に含まれる高調波成分を抽出する受信処理部と、抽出された高調波成分に基づいて超音波画像の画像データを形成する画像形成部と、を有することを特徴とする。   In order to achieve the above object, an ultrasonic diagnostic apparatus according to a preferred aspect of the present invention includes a probe that transmits / receives ultrasonic waves to / from a diagnostic region including a bubble, and a transmission that is set based on the resonance frequency of the bubble. A transmission control unit for controlling the probe to transmit an ultrasonic wave at a frequency and a transmission sound pressure set on the basis of the expansion rate of the bubble, and a received signal obtained from the diagnostic region by receiving the ultrasonic wave at the probe A reception processing unit that extracts a harmonic component contained in a received signal by processing the image, and an image forming unit that forms image data of an ultrasonic image based on the extracted harmonic component. To do.

上記態様によれば、バブルの共振周波数を基準として送信周波数が設定され、バブルの膨張倍率を基準として送信音圧が設定されるため、例えば、バブルから得られる高調波成分を高めると共に実組織から得られる高調波成分を抑制することができ、その結果として、例えば、非常に高いCTR(Contrast to Tissue Ratio)を得ることができる。   According to the above aspect, the transmission frequency is set on the basis of the resonance frequency of the bubble, and the transmission sound pressure is set on the basis of the expansion rate of the bubble. As a result, for example, a very high CTR (Contrast to Tissue Ratio) can be obtained.

望ましい態様において、前記送信制御部は、前記バブルの膨張倍率と前記バブルから得られる高調波成分の信号量との対応関係に基づいて、信号量が比較的小さく且つ膨張倍率の変化に対する信号量の変化が比較的大きい対応関係を示す膨張倍率が小さな範囲と、信号量が比較的大きく且つ膨張倍率の変化に対する信号量の変化が比較的小さい対応関係を示す膨張倍率が大きな範囲と、からなる二つの範囲のうち、膨張倍率が大きな範囲に収まるように前記送信音圧を設定することを特徴とする。   In a desirable mode, the transmission control unit has a relatively small signal amount based on the correspondence relationship between the expansion rate of the bubble and the signal amount of the harmonic component obtained from the bubble, and the signal amount with respect to the change in the expansion rate. A range in which the expansion ratio indicating a relatively large change relationship is small, and a range in which the signal amount is relatively large and a change amount of the signal amount relative to the change in expansion magnification is relatively small. Of the two ranges, the transmission sound pressure is set so that the expansion magnification is within a large range.

望ましい態様において、前記送信制御部は、前記バブルの共振周波数以下となるように前記送信周波数を設定し、前記バブルの膨張倍率が1.5以上となるように前記送信音圧を設定することを特徴とする。   In a preferred aspect, the transmission control unit sets the transmission frequency so as to be equal to or lower than a resonance frequency of the bubble, and sets the transmission sound pressure so that an expansion ratio of the bubble is 1.5 or more. Features.

望ましい態様において、前記送信制御部は、共に前記送信周波数かつ前記送信音圧に設定されて互いに位相が反転関係にある超音波の組を送波するようにプローブを制御し、前記受信処理部は、前記超音波の組のうちの一方に対応した第1受信信号と他方に対応した第2受信信号の差分を算出することにより奇数次の高調波成分を抽出し、さらに、周波数フィルタを用いて奇数次の高調波成分から第3次の高調波成分を抽出し、前記画像形成部は、抽出された第3次の高調波成分に基づいて超音波画像の画像データを形成することを特徴とする。   In a desirable aspect, the transmission control unit controls the probe so as to transmit a set of ultrasonic waves that are set to the transmission frequency and the transmission sound pressure and have phases inverted to each other, and the reception processing unit The odd harmonic component is extracted by calculating the difference between the first received signal corresponding to one of the ultrasonic sets and the second received signal corresponding to the other, and using a frequency filter A third-order harmonic component is extracted from odd-order harmonic components, and the image forming unit forms image data of an ultrasonic image based on the extracted third-order harmonic component. To do.

望ましい態様において、前記受信処理部は、前記第3次の高調波成分に加えて、前記第1受信信号と第2受信信号を加算することにより第2次の高調波成分を抽出し、前記画像形成部は、抽出された第2次の高調波成分と第3次の高調波成分に基づいて超音波画像の画像データを形成することを特徴とする。   In a preferred aspect, the reception processing unit extracts the second harmonic component by adding the first received signal and the second received signal in addition to the third harmonic component, and the image The forming unit forms image data of an ultrasonic image based on the extracted second harmonic component and third harmonic component.

望ましい態様において、前記プローブは、前記送信周波数から前記送信周波数の3倍の周波数までを含む周波数帯域を備えることを特徴とする。   In a preferred aspect, the probe includes a frequency band including the transmission frequency to a frequency three times the transmission frequency.

本発明により、バブルの挙動に基づいた適切な画像処理技術が提供される。例えば、本発明の好適な態様によれば、バブルの共振周波数を基準として送信周波数が設定され、バブルの膨張倍率を基準として送信音圧が設定されるため、バブルから得られる高調波成分を高めると共に実組織から得られる高調波成分を抑制することができ、その結果として、例えば非常に高いCTRを得ることができる。   According to the present invention, an appropriate image processing technique based on the behavior of bubbles is provided. For example, according to a preferred aspect of the present invention, the transmission frequency is set based on the resonance frequency of the bubble, and the transmission sound pressure is set based on the expansion rate of the bubble, so that the harmonic component obtained from the bubble is increased. In addition, harmonic components obtained from the real tissue can be suppressed, and as a result, for example, a very high CTR can be obtained.

以下、本発明の好適な実施形態を図面に基づいて説明する。図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1は、その全体構成を示す機能ブロック図である。図1に示す超音波診断装置は、造影用のバブル(マイクロバブルやナノバブルなどの微小気泡)を含んだ造影剤を利用して画像を形成するのに適している。   DESCRIPTION OF EXEMPLARY EMBODIMENTS Hereinafter, preferred embodiments of the invention will be described with reference to the drawings. FIG. 1 shows a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention, and FIG. 1 is a functional block diagram showing the overall configuration thereof. The ultrasonic diagnostic apparatus shown in FIG. 1 is suitable for forming an image using a contrast agent including a contrast bubble (microbubbles such as microbubbles and nanobubbles).

本実施形態においては、例えば、低中音圧の超音波による画像化に適した造影剤(「ソナゾイド」(登録商標)など)を利用することが望ましいものの、本発明における造影剤は特定のものに限定されない。造影剤は、例えば生体内の血管や腫瘍などの診断部位に投与される。そして、造影剤が投与されてから、例えば生体内におけるバブルの集積や取り込みが起こるまでの一定時間が経過した後に、図1に示す超音波診断装置によって診断が行われる。   In the present embodiment, for example, it is desirable to use a contrast agent (such as “Sonazoid” (registered trademark)) suitable for imaging with ultrasound of low and medium sound pressures, but the contrast agent in the present invention is a specific one. It is not limited to. The contrast agent is administered to a diagnostic site such as a blood vessel or a tumor in a living body. A diagnosis is performed by the ultrasonic diagnostic apparatus shown in FIG. 1 after a certain period of time has elapsed from the administration of the contrast agent to, for example, the accumulation and uptake of bubbles in the living body.

信号発生部10は、図示しない制御部などによって制御され、送信パルスを形成するための駆動信号を生成して送信回路12へ出力する。本実施形態では、生体内に投与されるバブルの共振周波数を基準として設定される送信周波数で、さらに、当該バブルの膨張倍率を基準として設定される送信音圧で超音波を送波するように送信制御される。信号発生部10は、例えば、中心周波数が1.5MHz程度、送信音圧が200〜300kPa程度に対応した波形の送信パルスを出力する。   The signal generation unit 10 is controlled by a control unit (not shown) or the like, generates a drive signal for forming a transmission pulse, and outputs the drive signal to the transmission circuit 12. In the present embodiment, ultrasonic waves are transmitted at a transmission frequency set with reference to the resonance frequency of the bubble administered into the living body, and with a transmission sound pressure set with reference to the expansion ratio of the bubble. Transmission is controlled. For example, the signal generator 10 outputs a transmission pulse having a waveform corresponding to a center frequency of about 1.5 MHz and a transmission sound pressure of about 200 to 300 kPa.

送信回路12は、信号発生部10から出力された送信パルスに基づいて、プローブ14が備える図示しない複数の振動素子を制御して送信ビームを形成し、形成した送信ビームを電子的に走査することにより、走査領域の全域に亘って複数の送信ビームを形成する。   The transmission circuit 12 controls a plurality of vibration elements (not shown) included in the probe 14 based on the transmission pulse output from the signal generator 10 to form a transmission beam, and electronically scans the formed transmission beam. Thus, a plurality of transmission beams are formed over the entire scanning region.

本実施形態においては、複数の送信ビームの各ビーム方向ごとに2回の送信が行われる。つまり、一つのビーム方向に対して送信パルスに基づいて1回目の送信が行われて受信信号を取得した後に、続けてその同じビーム方向に対して、1回目の送信と同じ波形で位相だけを反転させた送信パルスに基づいて2回目の送信が行われて受信信号が取得される。そして、一つのビーム方向に対して2回の送信が行われた後に、ビーム方向を変えて次のビーム方向に対して2回の送信が行われる。こうして、走査領域の全域に亘って各ビーム方向ごとに2回の送信が実行される。   In the present embodiment, transmission is performed twice for each beam direction of a plurality of transmission beams. In other words, after the first transmission is performed for one beam direction based on the transmission pulse and the received signal is acquired, the phase is continuously changed with the same waveform as the first transmission for the same beam direction. A second transmission is performed based on the inverted transmission pulse to obtain a received signal. Then, after two transmissions are performed for one beam direction, the transmission is performed twice for the next beam direction by changing the beam direction. In this way, two transmissions are executed for each beam direction over the entire scanning region.

プローブ14は、造影剤が投与された生体内の診断領域に対して超音波を送受する。プローブ14は、超音波を送受する複数の振動素子を備えており、複数の振動素子が送信回路12によって送信制御されて送信ビームが走査される。また、複数の振動素子が生体から反射された超音波を受波し、これにより得られた信号が受信回路16へ出力される。なお、送信と受信を異なる振動子で行うようにしてもよい。   The probe 14 transmits and receives ultrasonic waves to the in-vivo diagnostic region to which the contrast agent is administered. The probe 14 includes a plurality of vibration elements that transmit and receive ultrasonic waves, and the transmission circuit 12 controls transmission of the plurality of vibration elements and scans the transmission beam. Further, the plurality of vibration elements receive the ultrasonic waves reflected from the living body, and a signal obtained thereby is output to the receiving circuit 16. Note that transmission and reception may be performed by different vibrators.

プローブ14は、送信周波数(例えば中心周波数が1.5MHz程度)から、この送信周波数の3倍の周波数までを含む周波数帯域を備えている。そのため、送信における送信周波数を基本波の周波数とすると、受信において基本波と第2次高調波と第3次高調波の成分を含んだ受信信号が得られる。   The probe 14 has a frequency band including a transmission frequency (for example, the center frequency is about 1.5 MHz) to a frequency three times the transmission frequency. Therefore, when the transmission frequency in transmission is the fundamental frequency, a reception signal including fundamental, second harmonic, and third harmonic components is obtained in reception.

受信回路16は、プローブ14の複数の振動素子から得られる信号を整相加算処理することにより、走査領域内の複数の送信ビームの各々に対応した受信信号を形成する。本実施形態においては、複数の送信ビームの各ビーム方向ごとに2回の送信が行われるため、受信回路16は、各ビーム方向ごとに、1回目の送信に対応した1回目の受信信号と2回目の送信に対応した2回目の受信信号を形成する。そして、各ビーム方向ごとに、1回目の受信信号(受信信号1)が受信信号1メモリ18に記憶され、2回目の受信信号(受信信号2)が受信信号2メモリ20に記憶される。   The reception circuit 16 forms a reception signal corresponding to each of the plurality of transmission beams in the scanning region by performing phasing addition processing on the signals obtained from the plurality of vibration elements of the probe 14. In this embodiment, since transmission is performed twice for each beam direction of a plurality of transmission beams, the reception circuit 16 receives the first received signal corresponding to the first transmission and 2 for each beam direction. A second received signal corresponding to the second transmission is formed. For each beam direction, the first received signal (received signal 1) is stored in the received signal 1 memory 18, and the second received signal (received signal 2) is stored in the received signal 2 memory 20.

加算処理部22は、各ビーム方向ごとに、受信信号1メモリ18と受信信号2メモリ20から受信信号1と受信信号2を読み出し、そして、受信信号1と受信信号2を加算処理する。受信信号1と受信信号2は、互いに位相が反転関係にある送信パルスから得られる信号であるため、フェーズインバージョン法(またはパルスインバージョン法)の原理により、受信信号1と受信信号2の加算処理の結果として、奇数次の信号がキャンセルされて偶数次の信号のみが残る。つまり、受信信号に含まれる基本波と第2次高調波と第3次高調波のうち、奇数次の基本波と第3次高調波がキャンセルされて偶数次の第2次高調波が抽出される。   The addition processing unit 22 reads the reception signal 1 and the reception signal 2 from the reception signal 1 memory 18 and the reception signal 2 memory 20 for each beam direction, and adds the reception signal 1 and the reception signal 2 to each other. Since reception signal 1 and reception signal 2 are signals obtained from transmission pulses whose phases are inverted with respect to each other, addition of reception signal 1 and reception signal 2 is performed according to the principle of phase inversion method (or pulse inversion method). As a result of the processing, the odd order signal is canceled and only the even order signal remains. That is, out of the fundamental wave, the second harmonic, and the third harmonic included in the received signal, the odd-order fundamental wave and the third-order harmonic are canceled, and the even-order second-order harmonic is extracted. The

バンドパスフィルタ(BPF)26は、第2次高調波に対応した通過帯域を備えたフィルタであり、バンドパスフィルタ26によって、ノイズ等が除去されて第2次高調波のみがより確実に抽出される。なお、原理的には加算処理部22において第2次高調波のみが確実に抽出されるため、バンドパスフィルタ26を省略することも可能である。   The band pass filter (BPF) 26 is a filter having a pass band corresponding to the second harmonic, and noise or the like is removed by the band pass filter 26 so that only the second harmonic is extracted more reliably. The Note that, in principle, only the second harmonic is reliably extracted in the addition processing unit 22, and therefore the bandpass filter 26 can be omitted.

差分処理部24は、各ビーム方向ごとに、受信信号1メモリ18と受信信号2メモリ20から受信信号1と受信信号2を読み出し、そして、受信信号1と受信信号2の差分を算出する。フェーズインバージョン法(またはパルスインバージョン法)の原理により、受信信号1と受信信号2の差分処理の結果として、偶数次の信号がキャンセルされて奇数次の信号のみが残る。つまり、受信信号に含まれる基本波と第2次高調波と第3次高調波のうち、偶数次の第2次高調波がキャンセルされて奇数次の基本波と第3次高調波が抽出される。   The difference processing unit 24 reads the reception signal 1 and the reception signal 2 from the reception signal 1 memory 18 and the reception signal 2 memory 20 for each beam direction, and calculates the difference between the reception signal 1 and the reception signal 2. Due to the principle of the phase inversion method (or pulse inversion method), as a result of differential processing between the received signal 1 and the received signal 2, the even-order signal is canceled and only the odd-order signal remains. That is, out of the fundamental wave, the second harmonic, and the third harmonic included in the received signal, the even-order second harmonic is canceled and the odd-order fundamental wave and the third harmonic are extracted. The

ハイパスフィルタ(HPF)28は、第3次高調波以上の帯域を備えたフィルタであり、ハイパスフィルタ28によって、基本波やノイズ等が除去されて第3次高調波のみが抽出される。   The high-pass filter (HPF) 28 is a filter having a band equal to or higher than the third harmonic, and the high-pass filter 28 removes the fundamental wave, noise, and the like and extracts only the third harmonic.

画像形成部30は、ハイパスフィルタ28から得られる第3次高調波のみを用いて、または、第3次高調波に加えてバンドパスフィルタ26から得られる第2次高調波を用いて、超音波画像の画像データを形成する。画像形成部30は、複数のビーム方向の各々から得られる高調波成分(第3次高調波のみ、または、第2次高調波と第3次高調波)に基づいて、例えば、2次元超音波画像の画像データを形成する。もちろん、超音波ビームを3次的に走査させて3次元超音波画像の画像データを形成してもよい。画像形成部30において形成された画像データに対応する超音波画像は表示部32に表示される。   The image forming unit 30 uses only the third harmonic obtained from the high-pass filter 28 or uses the second harmonic obtained from the bandpass filter 26 in addition to the third harmonic. Image data of an image is formed. The image forming unit 30 uses, for example, two-dimensional ultrasonic waves based on harmonic components (third harmonic only, or second harmonic and third harmonic) obtained from each of a plurality of beam directions. Image data of an image is formed. Of course, image data of a three-dimensional ultrasonic image may be formed by tertiary scanning with an ultrasonic beam. An ultrasonic image corresponding to the image data formed in the image forming unit 30 is displayed on the display unit 32.

本実施形態では、生体内に投与されるバブルの共振周波数を基準として送信周波数が設定され、バブルの膨張倍率を基準として送信音圧が設定されるため、例えば、バブルから得られる高調波成分を高めると共に実組織から得られる高調波成分を抑制することができ、その結果として、例えば、非常に高いCTR(Contrast to Tissue Ratio)の超音波画像を得ることができる。   In this embodiment, the transmission frequency is set with reference to the resonance frequency of the bubble administered into the living body, and the transmission sound pressure is set with reference to the expansion ratio of the bubble. As a result, it is possible to suppress the harmonic component obtained from the real tissue, and as a result, for example, an ultrasound image with a very high CTR (Contrast to Tissue Ratio) can be obtained.

そこで、以下において、送信条件等とバブルの挙動との関係について説明する。図2から図6は、超音波造影剤(バブル)の挙動について、バブル膜の振動を数値解析で得るシミュレータを作成して解析した結果を示している。   Therefore, in the following, the relationship between transmission conditions and the like and bubble behavior will be described. FIGS. 2 to 6 show the results of analyzing the behavior of the ultrasonic contrast agent (bubble) by creating a simulator that obtains the vibration of the bubble film by numerical analysis.

図2は、バブルの膜振動のシミュレーション結果を示す図である。図2(A)は、送信周波数(中心周波数)1.5MHz、波数6波、音圧300kPaの超音波振動を半径1.25μmのバブルに与えた場合の膜振動であり、図2(B)は、送信周波数(中心周波数)4.5MHz、波数18波、音圧300kPaの超音波振動を同バブルに与えた場合の膜振動を示している。図2(A)(B)は共に、横軸を時間軸として縦軸にバブル半径を示している。なお、バブルに与える超音波エネルギー(音圧×パルス長)を図2(A)と(B)で等しい条件として解析した。   FIG. 2 is a diagram showing a simulation result of bubble membrane vibration. FIG. 2A shows membrane vibration when ultrasonic vibration having a transmission frequency (center frequency) of 1.5 MHz, a wave number of 6 waves, and a sound pressure of 300 kPa is applied to a bubble having a radius of 1.25 μm. Shows membrane vibration when ultrasonic vibration having a transmission frequency (center frequency) of 4.5 MHz, a wave number of 18 waves, and a sound pressure of 300 kPa is applied to the bubble. 2A and 2B both show the bubble radius on the vertical axis with the horizontal axis as the time axis. The ultrasonic energy (sound pressure × pulse length) given to the bubble was analyzed under the same conditions in FIGS. 2 (A) and 2 (B).

図2(A)に示すように低周波1.5MHzでは3倍以上にバブルは膨張するが、図2(B)に示すように高周波数の4.5MHzでは1.4倍以下の膨張に抑えられている。バブルの条件として、市販造影剤ソナゾイド(登録商標)を想定し、バブルの直径2.5μm、バブルの膜厚3nm、それ以外のバブルのパラメータは、文献“Experimental and Theoretical Evaluation of Microbubble Behavior:IEEE Transaction on Ultrasonics,Ferroelectrics,and Frequency Control,Vol 47,No6,Nov 2000”で示された条件を用いた。この条件でのバブルの共振周波数は、ほぼ2.5MHzである。   As shown in FIG. 2 (A), the bubble expands more than 3 times at a low frequency of 1.5 MHz, but is suppressed to less than 1.4 times expansion at a high frequency of 4.5 MHz as shown in FIG. 2 (B). It has been. As a bubble condition, a commercially available contrast medium Sonazoid (registered trademark) is assumed. on Ultrasonics, Ferroelectrics, and Frequency Control, Vol 47, No6, Nov 2000 ”. The resonant frequency of the bubble under this condition is approximately 2.5 MHz.

図3は、図2のバブルの膜振動によって発生する超音波信号を周波数解析した結果を示す図である。図3(A)の解析結果は図2(A)の膜振動に対応しており、図3(B)の解析結果は図2(B)の膜振動に対応している。   FIG. 3 is a diagram showing the result of frequency analysis of the ultrasonic signal generated by the membrane vibration of the bubble of FIG. The analysis result of FIG. 3 (A) corresponds to the membrane vibration of FIG. 2 (A), and the analysis result of FIG. 3 (B) corresponds to the membrane vibration of FIG. 2 (B).

基本波の信号量は、図3(A)の1.5MHzと図3(B)の4.5MHzでほぼ同じである。これに対し、第3高調波量は、1.5MHzでは基本波の−2dB程度で殆ど変わらないが、4.5MHzでは基本波の−20dB以上と大きく減少する。   The signal amount of the fundamental wave is substantially the same at 1.5 MHz in FIG. 3A and 4.5 MHz in FIG. On the other hand, the third harmonic amount hardly changes at about −2 dB of the fundamental wave at 1.5 MHz, but greatly decreases to −20 dB or more of the fundamental wave at 4.5 MHz.

図4は、送信周波数と信号量との対応関係を示す図であり、図4には、音圧300kPaで異なる周波数で一定信号量(=音圧×時間)の振動を与えた場合の、バブル振動最大半径と、その時発生する1次(基本波=送信周波数)、2次、3次の高調波の解析結果のグラフが示されている。   FIG. 4 is a diagram showing a correspondence relationship between the transmission frequency and the signal amount. FIG. 4 shows a bubble in the case of applying a constant signal amount (= sound pressure × time) with different frequencies at a sound pressure of 300 kPa. A graph of the analysis result of the maximum vibration radius and the first-order (fundamental wave = transmission frequency), second-order and third-order harmonics generated at that time is shown.

バブルは、低い周波数で大きく振動し、その結果発生する超音波信号量は、その基本周波数(送信周波数)ではそれほど差は無いが、第2、第3高調波では、低い周波数で送信したほうが、強い信号が得られることとなる。   The bubble vibrates greatly at a low frequency, and the amount of ultrasonic signals generated as a result is not so different at the fundamental frequency (transmission frequency), but it is better to transmit at the lower frequency in the second and third harmonics, A strong signal will be obtained.

図5は、バブルの膨張倍率とバブルから得られる高調波成分の信号量との対応関係を示す図であり、半径1.25μmのバブルに周波数1.5、2.5、3.5、4.5MHz、音圧100、200、300、400kPaの条件で加える超音波エネルギーを一定にしてバブルに振動を与えた場合の、バブルの膨張倍率を横軸に、3次高調波信号量を縦軸にしてプロットしたものである。   FIG. 5 is a diagram illustrating a correspondence relationship between the expansion ratio of the bubble and the signal amount of the harmonic component obtained from the bubble. The frequency of 1.5, 2.5, 3.5, 4 When the bubble is vibrated with a constant ultrasonic energy applied under the conditions of 0.5 MHz, sound pressure 100, 200, 300, and 400 kPa, the expansion factor of the bubble is plotted on the horizontal axis and the third harmonic signal amount is plotted on the vertical axis. And plotted.

図5に示すように、バブルの膨張倍率が0.5から1.5の小さな範囲では、3次高調波の信号量が比較的小さく且つ膨張倍率の変化に対する信号量の変化が比較的大きい。また、バブルの膨張倍率が1.5以上の大きな範囲では、3次高調波の信号量が比較的大きく且つ膨張倍率の変化に対する信号量の変化が比較的小さい。このように、バブルを1.5倍以上に膨張させることで、強い3次高調波信号が得られることが推察される。また、バブルを2.5倍以上に膨張させても、高調波信号量の増加はほぼ横ばいとなる。   As shown in FIG. 5, when the expansion ratio of the bubble is a small range of 0.5 to 1.5, the signal amount of the third harmonic is relatively small and the change of the signal amount with respect to the change of the expansion ratio is relatively large. Further, in the large range where the expansion ratio of the bubble is 1.5 or more, the signal amount of the third harmonic is relatively large and the change of the signal amount with respect to the change of the expansion ratio is relatively small. Thus, it is guessed that a strong third harmonic signal can be obtained by expanding the bubble by 1.5 times or more. Moreover, even if the bubble is expanded 2.5 times or more, the increase in the harmonic signal amount is almost flat.

以上のことから、バブル(マイクロバブル)を1.5倍から2.0倍程度に膨張させることで、バブルから非常に強い高調波信号を発生させることができる。また、バブルの共振周波数よりも低い周波数の超音波を与えることで、より低音圧の超音波信号でバブルを大きく膨張させることが可能であることが分かる。   From the above, it is possible to generate a very strong harmonic signal from a bubble by expanding the bubble (microbubble) from about 1.5 times to about 2.0 times. It can also be seen that by applying an ultrasonic wave having a frequency lower than the resonance frequency of the bubble, the bubble can be greatly expanded by an ultrasonic signal having a lower sound pressure.

そこで、バブルを高感度に検出する条件を考える。超音波は、ティシュハーモニックで知られるように伝搬中に高調波成分を発生する。その高調波信号量は、ほぼ音圧の二乗に比例する。従って、共振周波数以下の低い周波数の超音波で送信することで、小さい音圧でバブルを1.5倍以上大きく膨張させることが可能となり、バブルからの高調波信号を増加し、組織からの高調波信号を低減することが可能となり、非常に高いCTR(Contrast to Tissue Ratio)を得ることができる。更に、バブルを1.5倍以上に膨張させることで、バブルからは第3次以上の高調波成分も強くなる。組織伝搬中の高調波信号は、高次になるほど大きく減少するため、バブルを1.5倍以上に膨張させ3次以上の高調波成分を利用して画像化することで、より高いCTRを得ることが可能となる。   Therefore, a condition for detecting bubbles with high sensitivity is considered. Ultrasound generates harmonic components during propagation, as is known in tissue harmonics. The amount of harmonic signal is approximately proportional to the square of the sound pressure. Therefore, by transmitting with ultrasonic waves having a low frequency lower than the resonance frequency, it becomes possible to expand the bubble by 1.5 times or more with a small sound pressure, increasing the harmonic signal from the bubble and increasing the harmonic from the tissue. Wave signals can be reduced, and a very high CTR (Contrast to Tissue Ratio) can be obtained. Furthermore, by expanding the bubble by 1.5 times or more, third-order or higher harmonic components are also strengthened from the bubble. Since harmonic signals during tissue propagation greatly decrease as the order increases, bubbles are expanded 1.5 times or more and imaged using harmonic components of the third or higher order to obtain a higher CTR. It becomes possible.

図6は、腹部用コンベックスプローブを用いて、組織から発生する高調波とバブルからの高調波をシミュレーションした結果を示す図である。組織中の非線形伝搬は、KZK方程式を用いて数値解析した。シミュレーションでは、送信周波数1.6MHz、深さ9cmで組織中の音圧が200kPaの条件で、組織伝搬中とバブルからの高次高調波信号を求めた。   FIG. 6 is a diagram showing a result of simulating harmonics generated from a tissue and harmonics from bubbles using an abdominal convex probe. Non-linear propagation in the tissue was numerically analyzed using the KZK equation. In the simulation, high-order harmonic signals from tissue propagation and bubbles were obtained under the conditions of a transmission frequency of 1.6 MHz, a depth of 9 cm, and a sound pressure in the tissue of 200 kPa.

図6(A)は、この時のバブル膜の振動結果である。バブルの半径は、1.25μmから2.2μmに約1.75倍に膨張している。図6(B)は、組織およびバブルの高調波信号を比較したものである。組織からの信号は、伝搬中に発生した信号が、バブルと同じ位置で組織中の反射体から反射し信号として解析した。図6(B)では、組織及びバブルからの基本波信号を0dBとして、それぞれの高調波信号量を相対的に比較した。バブルからの高調波信号の信号量は、第2次、第3次とも基本波信号とほぼ同じである。一方、組織からの高調波信号の信号量は、基本波信号に比べ、第2次高調波では、−30dB、第3次高調波では、−50dB以下に減少する。   FIG. 6A shows the vibration result of the bubble film at this time. The bubble radius expands about 1.75 times from 1.25 μm to 2.2 μm. FIG. 6B compares the harmonic signals of tissue and bubbles. The signal from the tissue was analyzed by reflecting the signal generated during propagation from the reflector in the tissue at the same position as the bubble. In FIG. 6B, the fundamental wave signals from the tissue and the bubble are set to 0 dB, and the respective harmonic signal amounts are relatively compared. The signal amount of the harmonic signal from the bubble is almost the same as that of the fundamental wave signal in both the second and third orders. On the other hand, the signal amount of the harmonic signal from the tissue decreases to −30 dB for the second harmonic and to −50 dB or less for the third harmonic compared to the fundamental signal.

以上のことから、バブル共振周波数以下の低周波で励振することで、低音圧送信でもバブルを1.5倍以上に膨張することが可能となり、強い高調波信号が発生するともに、組織中の高調波信号を抑制することができることが分かる。その結果、非常に高い検出感度とCTRが達成される。   From the above, by exciting at a low frequency below the bubble resonance frequency, it becomes possible to expand the bubble to 1.5 times or more even with low sound pressure transmission, generating a strong harmonic signal and generating harmonics in the tissue. It can be seen that the wave signal can be suppressed. As a result, very high detection sensitivity and CTR are achieved.

超音波診断装置では、周波数2.5MHz〜12MHzのプローブが標準的に利用されている。コントラストエコーで高い検出感度を得るためには、バブル共振周波数より低い周波数の送信が可能な、低い周波数のプローブを利用することが適している。例えば、中心周波数3MHz、比帯域80パーセントのプローブであれば、1.4MHzで送信し、その第3次高調波4.2MHzを受信することができる。送信周波1.4MHzは、プローブの帯域のかなり低下した領域での利用となるが、高い音圧を発生する必要がないため、十分に送信可能である。バブルの高調波信号は、第3次以降の高次でも強く発生するが、プローブの帯域から受信感度は得られないため第3次までの利用となる。もちろん、第3次以上の高調波に対応したプローブで第3次以上の高調波を利用してもよい。   In the ultrasonic diagnostic apparatus, a probe having a frequency of 2.5 MHz to 12 MHz is typically used. In order to obtain high detection sensitivity with contrast echo, it is suitable to use a low-frequency probe capable of transmitting a frequency lower than the bubble resonance frequency. For example, if the probe has a center frequency of 3 MHz and a relative bandwidth of 80 percent, it can transmit at 1.4 MHz and receive the third harmonic 4.2 MHz. Although the transmission frequency of 1.4 MHz is used in a region where the bandwidth of the probe is considerably reduced, it is not necessary to generate a high sound pressure, and can be sufficiently transmitted. The harmonic signal of the bubble is strongly generated even at higher orders after the third order. However, since reception sensitivity cannot be obtained from the band of the probe, it is used up to the third order. Of course, third and higher harmonics may be used with a probe corresponding to the third and higher harmonics.

図7は、図1の超音波診断装置による第3次の高調波信号を利用した画像形成処理を説明するためのフローチャートである。まず、生体に対して造影剤を注入後(S701)、生体内でのバブルの集積や取り込みが起こるまでの適当な時間が経過してから(S702)、観測が開始される。   FIG. 7 is a flowchart for explaining an image forming process using a third harmonic signal by the ultrasonic diagnostic apparatus of FIG. First, after injecting a contrast medium into a living body (S701), an appropriate time elapses until bubble accumulation or uptake occurs in the living body (S702), and observation is started.

観測が開始されると、第1番目の走査線上(送信ビーム方向)が選択され(S703)、位相が0°に設定された1回目の低音圧パルス信号が送信され(S704)、そのエコー信号である受信信号1が受信される(S705)。さらに、第1番目の走査線上において、1回目の送信信号の位相を180°変えた2回目のパルス信号が送信され(S706)、そのエコー信号である受信信号2が受信される(S707)。なお、S704とS706の送信においては、例えば、中心周波数が1.5MHz程度、送信音圧が200〜300kPa程度の超音波パルスが送波される。   When the observation is started, the first scanning line (transmission beam direction) is selected (S703), the first low sound pressure pulse signal whose phase is set to 0 ° is transmitted (S704), and the echo signal is transmitted. The received signal 1 is received (S705). Further, the second pulse signal obtained by changing the phase of the first transmission signal by 180 ° is transmitted on the first scanning line (S706), and the reception signal 2 that is the echo signal is received (S707). In the transmission in S704 and S706, for example, an ultrasonic pulse having a center frequency of about 1.5 MHz and a transmission sound pressure of about 200 to 300 kPa is transmitted.

第3次の高調波信号のみを利用する場合には、2つの受信信号の差分を取ることで(S708)、奇数時の信号が同相で加算されて偶数時の信号がキャンセルされる。なお、第5次以上の奇数次の高調波は、プローブの帯域外のため受信信号に含まれず、受信信号に含まれる第1(基本周波数)と第3次の信号成分がS708の処理で取り出される。次に、基本周波数は、組織信号にも多く含まれるので、ハイパスフィルタを通過させることで(S709)、第3次高調波成分のみを取り出す。そして、走査線を一つ移動させて(S710)、S704からS709までの処理を繰り返すことにより、複数の走査線から得られる第3次高調波によりコントラスト画像を構築することができる。   When only the third-order harmonic signal is used, by taking the difference between the two received signals (S708), the odd-numbered signals are added in phase and the even-numbered signals are canceled. The fifth and higher odd harmonics are not included in the received signal because they are out of the band of the probe, and the first (fundamental frequency) and third signal components included in the received signal are extracted in the process of S708. It is. Next, since many fundamental frequencies are also included in the tissue signal, only the third-order harmonic component is extracted by passing the high-pass filter (S709). Then, by moving one scanning line (S710) and repeating the processing from S704 to S709, a contrast image can be constructed from the third harmonics obtained from the plurality of scanning lines.

第3次高調波信号のみを用いて画像化することで、組織からの第3次が非常に小さいため、CTRの高い画像が得られる。なお、第3次高調波信号のみを用いる装置を構成する場合には、図1における加算処理部22とバンドパスフィルタ26を省略してもよい。   By imaging using only the third harmonic signal, since the third order from the tissue is very small, an image with a high CTR can be obtained. In the case of configuring an apparatus that uses only the third harmonic signal, the addition processing unit 22 and the bandpass filter 26 in FIG. 1 may be omitted.

なお、受信信号には、第2次の高調波も多く含まれている。そこで、第2次と第3次を合わせて利用することで、信号強度を増加することが出来る。そのため、エコーの弱い深部からの情報も画像化することが可能となる。また、信号量の帯域も広がるため、解像度の高い画像を得ることが可能となる。   The received signal contains a lot of second-order harmonics. Therefore, the signal intensity can be increased by using the second order and the third order together. Therefore, it is possible to image information from the deep part where the echo is weak. In addition, since the signal amount band is widened, an image with high resolution can be obtained.

図8は、図1の超音波診断装置による第2次と第3次の高調波信号を利用した画像形成処理を説明するためのフローチャートである。なお、図8におけるS801からS809までの処理は、図7におけるS701からS709までの処理と同じであるため、説明を省略する。   FIG. 8 is a flowchart for explaining an image forming process using the second and third harmonic signals by the ultrasonic diagnostic apparatus of FIG. Note that the processing from S801 to S809 in FIG. 8 is the same as the processing from S701 to S709 in FIG.

図8においては、第3次の高調波に加えて、2つの受信信号を加算することで(S810)、受信信号に含まれる偶数時の信号が同相で加算され、奇数時の信号がキャンセルされる。つまり、基本波と第3次高調波が取り除かれて2次の高調波成分が取り出される。そして、第2次と第3次の高調波信号を加算することで(S811)、広帯域で強い強度の信号を得ることが可能になる。なお、図8においても、走査線を一つ移動させて(S812)、S804からS811までの処理を繰り返すことで、複数の走査線から得られる第2次と第3次高調波によりコントラスト画像を構築することができる。   In FIG. 8, by adding two received signals in addition to the third harmonic (S810), even signals included in the received signals are added in phase, and odd signals are canceled. The That is, the fundamental wave and the third harmonic are removed, and the second harmonic component is extracted. Then, by adding the second-order and third-order harmonic signals (S811), it becomes possible to obtain a strong signal in a wide band. Also in FIG. 8, by moving one scanning line (S812) and repeating the processing from S804 to S811, a contrast image is obtained by the second and third harmonics obtained from a plurality of scanning lines. Can be built.

多数の反射体が波長に比べて密に、かつランダムに存在する場合、エコーと干渉パターンが現れる。このパターンはスペックルと呼ばれ、アーテイファクトとなる。バブルが多数存在する場合も、同様にスペックルが現れる。これらの干渉パターンは、周波数によって異なるパターンを呈することが知られており、それを低減する方法として、周波数コンパウンド法が知られている。同様に、秋山らはBモード画像に現れるスペックルノイズを低減する目的で、高次のハーモニックイメージを加算する方法を提案している(山本、小沢、秋山「スペックル軽減のための高次ハーモニックイメージの加算方法について」電子情報通信学会技術研究報告Vol.102,No.411(20021018)pp.19-24参照)。   When many reflectors exist densely and randomly compared to the wavelength, echoes and interference patterns appear. This pattern is called speckle and is an artifact. Speckle appears in the same way when there are many bubbles. These interference patterns are known to exhibit different patterns depending on the frequency, and a frequency compound method is known as a method for reducing the interference pattern. Similarly, Akiyama et al. Have proposed a method of adding higher-order harmonic images for the purpose of reducing speckle noise appearing in B-mode images (Yamamoto, Ozawa, Akiyama “Higher-order harmonics for speckle reduction” About the method of adding images ”, IEICE Technical Report Vol.102, No.411 (20021018) pp.19-24).

図9は、インコヒーレントな重ね合わせによりスペックルを低減させる処理を説明するためのフローチャートである。つまり、図9は、バブルからの第2次と第3次の高調波によるインコヒーレントな重ね合わせによる、コントラストエコーのスペックル低減を目的とした処理のフローを示している。なお、図9におけるS901からS910までの処理は、図8におけるS801からS810までの処理と同じであるため、説明を省略する。   FIG. 9 is a flowchart for explaining processing for reducing speckles by incoherent superposition. That is, FIG. 9 shows a flow of processing for the purpose of reducing speckles of contrast echo by incoherent superimposition by the second and third harmonics from the bubble. Note that the processing from S901 to S910 in FIG. 9 is the same as the processing from S801 to S810 in FIG.

図9においては、S909の処理によって抽出される第3次高調波成分に対して振幅演算処理が施され(S911)、その結果として、第3次高調波成分のエンベローブが得られる。同様に、S910の処理によって抽出される第2次高調波成分に対して振幅演算処理が施され(S912)、第2次高調波成分のエンベローブが得られる。そして、第3次高調波成分のエンベローブと第2次高調波成分のエンベローブが加算処理されることにより(S913)、周波数コンパウンドが実現される。   In FIG. 9, the amplitude calculation process is performed on the third harmonic component extracted by the process of S909 (S911), and as a result, the envelope of the third harmonic component is obtained. Similarly, an amplitude calculation process is performed on the second harmonic component extracted by the process of S910 (S912), and an envelope of the second harmonic component is obtained. The envelope of the third harmonic component and the envelope of the second harmonic component are added (S913), thereby realizing frequency compound.

なお、図9においても、走査線を一つ移動させて(S914)、S904からS913までの処理が繰り返されて超音波画像が形成される。   Also in FIG. 9, one scanning line is moved (S914), and the processing from S904 to S913 is repeated to form an ultrasonic image.

図10は、コヒーレントな重ね合わせによりスペックルを低減させる処理を説明するためのフローチャートである。つまり、図10は、バブルからの第2次と第3次の高調波によるコヒーレントな重ね合わせによる、コントラストエコーのスペックル低減を目的とした処理のフローを示している。なお、図10におけるS1001からS1010までの処理は、図9におけるS901からS910までの処理と同じであるため説明を省略する。   FIG. 10 is a flowchart for explaining processing for reducing speckles by coherent superposition. That is, FIG. 10 shows a flow of processing for the purpose of reducing speckles of contrast echoes by coherent superposition by the second and third harmonics from the bubble. Note that the processing from S1001 to S1010 in FIG. 10 is the same as the processing from S901 to S910 in FIG.

図10においては、S1009の処理によって抽出される第3次高調波成分に対して複素解析信号処理が施され(S1011)、その結果として、第3次高調波成分の実数部と虚数部が得られる。同様に、S1010の処理によって抽出される第2次高調波成分に対して複素解析信号処理が施され(S1012)、第2次高調波成分の実数部と虚数部が得られる。   In FIG. 10, complex analysis signal processing is performed on the third harmonic component extracted by the processing of S1009 (S1011), and as a result, a real part and an imaginary part of the third harmonic component are obtained. It is done. Similarly, complex analysis signal processing is performed on the second harmonic component extracted by the processing of S1010 (S1012), and a real part and an imaginary part of the second harmonic component are obtained.

そして、第3次高調波成分の実数部と第2次高調波成分の実数部とが加算され(S1013)、また、第3次高調波成分の虚数部と第2次高調波成分の虚数部とが加算される(S1014)。さらに、加算後の実数部と加算後の虚数部とに基づいて、振幅演算処理が施される(S1015)ことにより、第2次と第3次の高調波によるコヒーレントな重ね合わせが実現される。なお、図10においても、走査線を一つ移動させて(S1016)、S1004からS1015までの処理が繰り返されて超音波画像が形成される。   Then, the real part of the third harmonic component and the real part of the second harmonic component are added (S1013), and the imaginary part of the third harmonic component and the imaginary part of the second harmonic component are added. Are added (S1014). Further, the amplitude calculation process is performed based on the real part after the addition and the imaginary part after the addition (S1015), thereby realizing coherent superposition by the second and third harmonics. . Also in FIG. 10, one scanning line is moved (S1016), and the processing from S1004 to S1015 is repeated to form an ultrasound image.

以上、本発明の好適な実施形態を説明したが、上述した実施形態は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。   As mentioned above, although preferred embodiment of this invention was described, embodiment mentioned above is only a mere illustration in all the points, and does not limit the scope of the present invention. The present invention includes various modifications without departing from the essence thereof.

本発明に係る超音波診断装置の全体構成を示す機能ブロック図である。1 is a functional block diagram showing an overall configuration of an ultrasonic diagnostic apparatus according to the present invention. バブルの膜振動のシミュレーション結果を示す図である。It is a figure which shows the simulation result of the membrane vibration of a bubble. バブルの膜振動によって発生する超音波信号を示す図である。It is a figure which shows the ultrasonic signal generated by the membrane vibration of a bubble. 送信周波数と信号量との対応関係を示す図である。It is a figure which shows the correspondence of a transmission frequency and signal amount. バブルの膨張倍率とバブルから得られる高調波成分の信号量との対応関係を示す図である。It is a figure which shows the correspondence of the expansion ratio of a bubble and the signal amount of the harmonic component obtained from a bubble. 組織から発生する高調波とバブルからの高調波をシミュレーションした結果を示す図である。It is a figure which shows the result of having simulated the harmonic generated from a structure | tissue, and the harmonic from a bubble. 第3次の高調波信号を利用した画像形成処理を説明するためのフローチャートである。It is a flowchart for demonstrating the image formation process using a 3rd harmonic signal. 第2次と第3次の高調波信号を利用した画像形成処理を説明するためのフローチャートである。It is a flowchart for demonstrating the image formation process using the 2nd and 3rd harmonic signal. インコヒーレントな重ね合わせによりスペックルを低減させる処理を説明するためのフローチャートである。It is a flowchart for demonstrating the process which reduces a speckle by incoherent superimposition. コヒーレントな重ね合わせによりスペックルを低減させる処理を説明するためのフローチャートである。It is a flowchart for demonstrating the process which reduces a speckle by coherent superimposition.

符号の説明Explanation of symbols

10 信号発生部、14 プローブ、18 受信信号1メモリ、20 受信信号2メモリ、22 加算処理部、24 差分処理部、30 画像形成部。   DESCRIPTION OF SYMBOLS 10 Signal generation part, 14 Probe, 18 Received signal 1 memory, 20 Received signal 2 memory, 22 Addition processing part, 24 Difference processing part, 30 Image formation part

Claims (6)

バブルを含む診断領域に対して超音波を送受するプローブと、
前記バブルの共振周波数を基準として設定される送信周波数かつ前記バブルの膨張倍率を基準として設定される送信音圧で超音波を送波するようにプローブを制御する送信制御部と、
プローブにおける超音波の受波により前記診断領域から得られる受信信号を処理することにより受信信号に含まれる高調波成分を抽出する受信処理部と、
抽出された高調波成分に基づいて超音波画像の画像データを形成する画像形成部と、
を有する、
ことを特徴とする超音波診断装置。
A probe for transmitting and receiving ultrasound to and from a diagnostic region including a bubble;
A transmission control unit for controlling the probe so as to transmit an ultrasonic wave with a transmission frequency set with reference to the resonance frequency of the bubble and a transmission sound pressure set with reference to the expansion rate of the bubble;
A reception processing unit that extracts a harmonic component included in the reception signal by processing the reception signal obtained from the diagnostic region by reception of ultrasonic waves in the probe;
An image forming unit that forms image data of an ultrasonic image based on the extracted harmonic component;
Having
An ultrasonic diagnostic apparatus.
請求項1に記載の超音波診断装置において、
前記送信制御部は、前記バブルの膨張倍率と前記バブルから得られる高調波成分の信号量との対応関係に基づいて、信号量が比較的小さく且つ膨張倍率の変化に対する信号量の変化が比較的大きい対応関係を示す膨張倍率が小さな範囲と、信号量が比較的大きく且つ膨張倍率の変化に対する信号量の変化が比較的小さい対応関係を示す膨張倍率が大きな範囲と、からなる二つの範囲のうち、膨張倍率が大きな範囲に収まるように前記送信音圧を設定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 1,
The transmission control unit has a relatively small signal amount based on the correspondence relationship between the expansion ratio of the bubble and the signal amount of the harmonic component obtained from the bubble, and the change in the signal amount with respect to the change in the expansion ratio is relatively small. Of the two ranges consisting of a range in which the expansion ratio indicating a large correspondence is small and a range in which the signal amount is relatively large and a change in the signal amount relative to a change in expansion magnification is relatively small , The transmission sound pressure is set so that the expansion ratio is within a large range,
An ultrasonic diagnostic apparatus.
請求項2に記載の超音波診断装置において、
前記送信制御部は、前記バブルの共振周波数以下となるように前記送信周波数を設定し、前記バブルの膨張倍率が1.5以上となるように前記送信音圧を設定する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 2,
The transmission control unit sets the transmission frequency so as to be equal to or lower than the resonance frequency of the bubble, and sets the transmission sound pressure so that an expansion ratio of the bubble is 1.5 or more.
An ultrasonic diagnostic apparatus.
請求項1から3のいずれか1項に記載の超音波診断装置において、
前記送信制御部は、共に前記送信周波数かつ前記送信音圧に設定されて互いに位相が反転関係にある超音波の組を送波するようにプローブを制御し、
前記受信処理部は、前記超音波の組のうちの一方に対応した第1受信信号と他方に対応した第2受信信号の差分を算出することにより奇数次の高調波成分を抽出し、さらに、周波数フィルタを用いて奇数次の高調波成分から第3次の高調波成分を抽出し、
前記画像形成部は、抽出された第3次の高調波成分に基づいて超音波画像の画像データを形成する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to any one of claims 1 to 3,
The transmission control unit controls the probe so as to transmit a set of ultrasonic waves that are both set to the transmission frequency and the transmission sound pressure and have a phase inversion relationship with each other,
The reception processing unit extracts an odd-order harmonic component by calculating a difference between a first received signal corresponding to one of the ultrasonic sets and a second received signal corresponding to the other, and Extract third-order harmonic components from odd-order harmonic components using a frequency filter,
The image forming unit forms image data of an ultrasonic image based on the extracted third harmonic component;
An ultrasonic diagnostic apparatus.
請求項4に記載の超音波診断装置において、
前記受信処理部は、前記第3次の高調波成分に加えて、前記第1受信信号と第2受信信号を加算することにより第2次の高調波成分を抽出し、
前記画像形成部は、抽出された第2次の高調波成分と第3次の高調波成分に基づいて超音波画像の画像データを形成する、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 4,
The reception processing unit extracts a second harmonic component by adding the first reception signal and the second reception signal in addition to the third harmonic component,
The image forming unit forms image data of an ultrasonic image based on the extracted second harmonic component and third harmonic component.
An ultrasonic diagnostic apparatus.
請求項4または5に記載の超音波診断装置において、
前記プローブは、前記送信周波数から前記送信周波数の3倍の周波数までを含む周波数帯域を備える、
ことを特徴とする超音波診断装置。
The ultrasonic diagnostic apparatus according to claim 4 or 5,
The probe includes a frequency band including from the transmission frequency to a frequency three times the transmission frequency.
An ultrasonic diagnostic apparatus.
JP2007318814A 2007-12-10 2007-12-10 Ultrasonic diagnostic equipment Expired - Fee Related JP5043623B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007318814A JP5043623B2 (en) 2007-12-10 2007-12-10 Ultrasonic diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007318814A JP5043623B2 (en) 2007-12-10 2007-12-10 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JP2009136626A JP2009136626A (en) 2009-06-25
JP5043623B2 true JP5043623B2 (en) 2012-10-10

Family

ID=40867891

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007318814A Expired - Fee Related JP5043623B2 (en) 2007-12-10 2007-12-10 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP5043623B2 (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6186737B2 (en) * 2013-02-04 2017-08-30 セイコーエプソン株式会社 Ultrasonic measuring device, ultrasonic imaging device, and ultrasonic measuring method
JP6318594B2 (en) * 2013-12-11 2018-05-09 コニカミノルタ株式会社 Ultrasound diagnostic imaging equipment
CN105307574B (en) * 2014-01-30 2018-03-30 皇家飞利浦有限公司 For detecting the method and apparatus of bubble using ultrasound
CN107205722A (en) * 2015-01-29 2017-09-26 皇家飞利浦有限公司 The fundamental wave harmonic frequency ultrasound diagnosing image of broadband mixing
JP2023059145A (en) * 2021-10-14 2023-04-26 Hoya株式会社 electronic endoscope system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1199152A (en) * 1997-09-26 1999-04-13 Ge Yokogawa Medical Systems Ltd Ultrasonograph
JP4744727B2 (en) * 2001-06-06 2011-08-10 株式会社東芝 Ultrasound diagnostic imaging equipment
JP3908555B2 (en) * 2002-02-08 2007-04-25 株式会社東芝 Ultrasonic diagnostic equipment
US8055036B2 (en) * 2005-05-27 2011-11-08 Hitachi Medical Corporation Ultrasonic diagnostic apparatus and ultrasonic image display method
JP4464432B2 (en) * 2007-08-23 2010-05-19 Geヘルスケア・ジャパン株式会社 Ultrasonic imaging device

Also Published As

Publication number Publication date
JP2009136626A (en) 2009-06-25

Similar Documents

Publication Publication Date Title
JP4627556B2 (en) Ultrasonic diagnostic equipment
JP4667733B2 (en) Method and apparatus for enhancing ultrasound contrast imaging using stepped chirp waveform
CA2659645C (en) Pulse inversion sequences for nonlinear imaging
JP5043623B2 (en) Ultrasonic diagnostic equipment
Saris et al. A comparison between compounding techniques using large beam-steered plane wave imaging for blood vector velocity imaging in a carotid artery model
Lediju et al. A motion-based approach to abdominal clutter reduction
JP5247322B2 (en) Ultrasonic imaging device
Shen et al. Pulse-inversion-based fundamental imaging for contrast detection
Emmer et al. Radial modulation of single microbubbles
JP5325670B2 (en) Ultrasonic diagnostic equipment
Lin et al. Ultrasound contrast imaging: influence of scatterer motion in multi-pulse techniques
Shen et al. Third harmonic transmit phasing for tissue harmonic generation
Nasholm et al. Transmit beams adapted to reverberation noise suppression using dual-frequency SURF imaging
Vos et al. Parametric array technique for microbubble excitation
JP2010063493A (en) Ultrasonic diagnostic apparatus
Frijlink et al. A simulation study on tissue harmonic imaging with a single-element intravascular ultrasound catheter
Du et al. Feasibility of non-linear simulation for Field II using an angular spectrum approach
JP5186155B2 (en) Ultrasonic diagnostic equipment
Novell et al. Optimization of multi-pulse sequences for nonlinear contrast agent imaging using a cMUT array
de Jong et al. 2C-4 Chirp Reversal Ultrasound Contrast Imaging
Zhu et al. Detection of scatters motion induced by mechanical vibrator using 7-chip barker-coded excitation
Saris et al. Does large beam-steered angle velocity compounding using plane wave transmission improve blood vector velocity estimation in a 3D carotid artery flow field?
Boulos Ultrasound imaging of the ultrasound thrombolysis
Karlinsky Frequency Mixing as a Platform for Enhanced Contrast Harmonic Imaging with Microbubbles
Basset et al. Ultrasound imaging using CMUT–Techniques developed in the frame of the ANR BBMUT project

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100927

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20120702

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20120710

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120712

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150720

Year of fee payment: 3

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees