JP2010110503A - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- JP2010110503A JP2010110503A JP2008286528A JP2008286528A JP2010110503A JP 2010110503 A JP2010110503 A JP 2010110503A JP 2008286528 A JP2008286528 A JP 2008286528A JP 2008286528 A JP2008286528 A JP 2008286528A JP 2010110503 A JP2010110503 A JP 2010110503A
- Authority
- JP
- Japan
- Prior art keywords
- signal
- doppler
- frequency
- wave
- diagnostic apparatus
- 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.)
- Pending
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】連続波を利用して目標位置からドプラ情報を抽出する技術において抽出の精度を高める。
【解決手段】FM変調器20から出力されるFM連続波が遅延回路25において遅延処理され、ミキサ32には遅延処理されたFM連続波が直接供給され、一方、ミキサ34には遅延処理されたFM連続波がπ/2シフト回路26を経由して供給される。遅延回路25において、目的とする深さからの受信信号の位相と参照波の位相を一致させることにより、目的とする深さからの復調信号を選択的に抽出することができる。さらに、半周期遅延回路35,37の各々において、FM連続波の周期の半周期に対応した遅延量だけ、復調信号に対して遅延処理が実行され、ミキサ36,38の各々において、復調信号と遅延された復調信号との乗算に相当する処理が実行され、復調信号に含まれるドプラ信号の周期的な変動成分が除去される。
【選択図】図1
【解決手段】FM変調器20から出力されるFM連続波が遅延回路25において遅延処理され、ミキサ32には遅延処理されたFM連続波が直接供給され、一方、ミキサ34には遅延処理されたFM連続波がπ/2シフト回路26を経由して供給される。遅延回路25において、目的とする深さからの受信信号の位相と参照波の位相を一致させることにより、目的とする深さからの復調信号を選択的に抽出することができる。さらに、半周期遅延回路35,37の各々において、FM連続波の周期の半周期に対応した遅延量だけ、復調信号に対して遅延処理が実行され、ミキサ36,38の各々において、復調信号と遅延された復調信号との乗算に相当する処理が実行され、復調信号に含まれるドプラ信号の周期的な変動成分が除去される。
【選択図】図1
Description
本発明は、超音波診断装置に関し、特に、連続波を利用する超音波診断装置に関する。
超音波診断装置の連続波を利用した技術として、連続波ドプラが知られている。連続波ドプラでは、例えば、数MHzの正弦波である送信波が生体内へ連続的に放射され、生体内からの反射波が連続的に受波される。反射波には、生体内における運動体(例えば血流など)によるドプラシフト情報が含まれる。そこで、そのドプラシフト情報を抽出して周波数解析することにより、運動体の速度情報を反映させたドプラ波形などを形成することができる。
連続波を利用した連続波ドプラは、パルス波を利用したパルスドプラに比べて一般に高速の速度計測の面で優れている。こうした事情などから、本願発明者は、連続波ドプラに関する研究を重ねてきた。その成果の一つとして、特許文献1において、周波数変調処理を施した連続波ドプラ(FMCWドプラ)に関する技術を提案している。
一方、連続波ドプラでは、連続波を利用していることにより位置計測が困難である。例えば、従来の一般的な連続波ドプラの装置(FMCWドプラを利用しない装置)では、位置計測を行うことができなかった。これに対し、本願発明者は、特許文献2において、FMCWドプラにより生体内組織の速度に加えて生体内組織の位置を計測することができる技術を提案している。
特許文献1や特許文献2に記載されたFMCWドプラの技術は、それまでにない超音波診断の可能性を秘めた画期的な技術である。本願発明者は、この画期的な技術の改良についてさらに研究を重ねてきた。特に、連続波を利用して目標位置からドプラ情報を抽出する技術に注目して研究を重ねてきた。
本発明は、このような背景において成されたものであり、その目的は、連続波を利用して目標位置からドプラ情報を抽出する技術において抽出の精度を高めることにある。
上記目的を達成するために、本発明の好適な態様である超音波診断装置は、周波数を周期的に変化させた連続波の送信信号を出力する送信信号処理部と、前記送信信号に基づいて超音波を生体に送波して生体からの反射波を受波することにより得られる受信信号に対して、前記送信信号に実質的に等しい波形の参照信号を用いて復調処理を施す受信信号処理部と、復調処理により得られた信号からドプラ信号を抽出するドプラ信号抽出部と、を有し、前記受信信号処理部は、生体内の目標位置から得られる受信信号との間の相関関係が調整された参照信号を用いて前記復調処理を施すことにより、前記目標位置に対応した復調信号を選択的に抽出し、前記ドプラ信号抽出部は、復調信号と当該復調信号に対して周期的に反転関係にある反転復調信号とに基づいて、前記連続波の周期性に伴うドプラ周波数の変動成分を低減し、前記目標位置に対応したドプラ信号を抽出する、ことを特徴とする。
上記態様では、参照信号を用いて受信信号に対して復調処理を施しているため、参照信号との相関が比較的大きい信号成分を含んだ復調信号を得ることができる。そして、その復調処理にあたり、目標位置から得られる受信信号と参照信号との間の相関関係が調整され、例えば、参照信号との相関が比較的大きい信号成分として目標位置からの受信信号を抽出することができる。さらに、連続波の周期性に伴うドプラ周波数の変動成分を低減してドプラ信号を抽出するため、目標位置に対応したドプラ信号に関する抽出の精度を高めることができる。
なお、上記態様において、参照信号と送信信号は、完全に等しい波形であることが望ましい。但し、参照信号と送信信号は、実質的に等しい波形とみなせる程度の対応関係であってもよい。
望ましい態様において、前記送信信号処理部は、周期Tmで周波数を変化させた連続波の送信信号を出力し、前記ドプラ信号抽出部は、復調信号に対して前記周期Tmの半周期に対応した遅延量だけ遅延処理を施すことにより反転復調信号を得ることを特徴とする。
望ましい態様において、前記ドプラ信号抽出部は、復調信号と反転復調信号との乗算結果に基づいて、ドプラ周波数の前記変動成分を低減する、ことを特徴とする。
望ましい態様において、前記ドプラ信号抽出部は、復調信号と反転復調信号を乗算することにより、復調信号に含まれるドプラ周波数と反転復調信号に含まれるドプラ周波数とを加算して前記変動成分を相殺する、ことを特徴とする。
望ましい態様において、前記ドプラ信号抽出部は、復調信号と反転復調信号を乗算することにより得られる乗算信号から、前記加算されたドプラ周波数に対応した周波数成分を前記ドプラ信号として抽出する、ことを特徴とする。
望ましい態様において、前記超音波診断装置は、前記目標位置からの受信信号の位相と前記参照信号の位相が等しくなるように、目標位置の深さに応じた遅延量だけ前記参照信号を遅延処理することにより、目標位置から得られる受信信号と参照信号との間の相関を強める、ことを特徴とする。
本発明により、連続波を利用して目標位置からドプラ情報を抽出する技術において抽出の精度が高められる。例えば、本発明の好適な態様によれば、連続波の周期性に伴うドプラ周波数の変動成分が低減されてドプラ信号が抽出される。
以下、本発明の好適な実施形態を説明する。
図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1はその全体構成を示す機能ブロック図である。送信用振動子10は生体内へ送信波を連続的に送波し、また、受信用振動子12は生体内からの反射波を連続的に受波する。このように、送信および受信がそれぞれ異なる振動子で行われて、いわゆる連続波ドプラ法による送受信が実行される。なお、送信用振動子10は複数の振動素子を備えており、これら複数の振動素子が制御されて超音波の送信ビームが形成される。また、受信用振動子12も複数の振動素子を備えており、これら複数の振動素子により得られた信号が処理されて受信ビームが形成される。
送信ビームフォーマ(送信BF)14は、送信用振動子10が備える複数の振動素子に対して送信信号を出力する。送信ビームフォーマ14には、例えば正弦波によるFM変調処理が施されたFM連続波(FMCW波)が入力される。送信ビームフォーマ14は、FM連続波に対して、各振動素子に応じた遅延処理を施して各振動素子に対応した送信信号を形成する。なお、送信ビームフォーマ14において形成された各振動素子に対応した送信信号に対して、必要に応じて電力増幅処理が施されてもよい。こうして、FM連続波による送信ビームが形成される。
FM変調器20は、送信ビームフォーマ14にFM連続波を出力する。FM変調器20は、変調波生成部24から供給される変調信号を用いて、RF波発振器22から供給されるRF波(搬送波信号)に対して周波数変調を施すことにより、FM連続波を発生する。このFM連続波の波形等については後に詳述する。
受信ビームフォーマ(受信BF)16は、受信用振動子12が備える複数の振動素子から得られる複数の受波信号を整相加算処理して受信ビームを形成する。つまり、受信ビームフォーマ16は、各振動素子から得られる受波信号に対してその振動素子に応じた遅延処理を施し、複数の振動素子から得られる複数の受波信号を加算処理することにより受信ビームを形成する。なお、各振動素子から得られる受波信号に対して低雑音増幅等の処理を施してから、受信ビームフォーマ16に複数の受波信号が供給されてもよい。こうして受信ビームに沿った受信RF信号が得られる。
受信ミキサ30は受信RF信号に対して直交検波を施して複素ベースバンド信号を生成する回路であり、2つのミキサ32,34で構成される。各ミキサは受信RF信号を所定の参照信号と混合する回路である。
受信ミキサ30の各ミキサに供給される参照信号は、FM変調器20から出力されるFM連続波に基づいて生成される。つまり、FM変調器20から出力されるFM連続波が遅延回路25において遅延処理され、ミキサ32には遅延処理されたFM連続波が直接供給され、一方、ミキサ34には遅延処理されたFM連続波がπ/2シフト回路26を経由して供給される。
π/2シフト回路26は遅延処理されたFM連続波の位相をπ/2だけずらす回路である。この結果、2つのミキサ32,34の一方から同相信号成分(I信号成分)が出力され、他方から直交信号成分(Q信号成分)が出力される。なお、必要に応じて、受信ミキサ30の後段にLPF(ローパスフィルタ)を設け、同相信号成分および直交信号成分の各々の高周波数成分をカットし、検波後の必要な帯域のみの復調信号を抽出してもよい。
後に詳述するが、各ミキサで実行される受信RF信号と参照信号との混合処理の結果である受信ミキサ出力信号(復調信号)には、変調波生成部24から供給される変調波の変調波周波数fmに関する複数の第n次波成分(nは0以上の自然数)が含まれている。つまり、第0次波成分である直流成分、第1次波成分である基本波成分、さらに、nが2以上の複数の高調波成分が含まれている。つまり、これら複数の第n次波成分を含んだ復調信号が、2つのミキサ32,34の各々から出力される。
さらに、本実施形態においては、ミキサ32から出力される復調信号に対して、半周期遅延回路35とミキサ36により、FM連続波の周期性に伴うドプラ周波数の変動成分(ドプラ折り返し電力)の低減処理が実行され、同様に、ミキサ34から出力される復調信号に対して、半周期遅延回路37とミキサ38により、FM連続波の周期性に伴うドプラ周波数の変動成分(ドプラ折り返し電力)の低減処理が実行される。これらの低減処理についても後に詳述する。
FFT処理部(高速フーリエ変換処理部)42は、復調信号(同相信号成分および直交信号成分)の各々に対してFFT演算を実行する。その結果、FFT処理部42において復調信号が周波数スペクトラムに変換される。なお、FFT処理部42から出力される周波数スペクトラムは、回路の設定条件などにより周波数分解能δfの周波数スペクトラムデータとして出力される。
ドプラ情報解析部44は、周波数スペクトラムに変換された復調信号からドプラ信号を抽出する。後に詳述するが、本実施形態では、変調波生成部24における変調処理と遅延回路25における遅延処理により目標位置が設定され、ドプラ情報解析部44において目標位置からのドプラ信号が選択的に抽出される。ドプラ情報解析部44は、生体内の各深さ(各位置)ごとにドプラ信号を抽出して、例えば、超音波ビーム(音線)上の各深さごとに生体内組織の速度を算出し、リアルタイムで出力する。なお、超音波ビームを走査させて二次元的あるいは三次元的に生体内組織の各位置の速度を算出してもよい。
表示処理部46は、生体組織の深さ(位置)ごとの速度に基づいて、例えばドプラ波形や、深さ速度の情報を含むグラフなどを形成し、形成したドプラ波形やグラフなどを表示部48にリアルタイムで表示させる。なお、図1に示す超音波診断装置内の各部は、システム制御部50によって制御される。つまり、システム制御部50は、送信制御や受信制御や表示制御などを行う。
以上、概説したように、本実施形態では、連続波(CW)を変調波でFM変調した超音波(FMCW波)を送受波して受信信号が得られて、目標位置からのドプラ情報が選択的に抽出される。そこで、目標位置からのドプラ情報が選択的に抽出される原理について詳述する。なお、図1に示した部分(構成)については、以下の説明においても図1の符号を利用する。
周波数f0のRF波(搬送波)に対して、周波数fmの正弦波によりFM変調を施したFMCW送信波は次式のように表現できる。
数1式において、Δfは周波数変動幅の0−P値(ゼロピーク値:最大周波数偏移)であり、最大周波数偏移Δfと変調周波数fmの比であるβはFMの変調指数である。
また、ドプラシフトを伴う場合のFMCW受信波は、生体における往復の減衰をαとすると次式で表現できる。
なお、数2式においてfmに対するドプラシフトは、f0のシフト分fdに比較して小さいので無視している。
図2は、本実施形態における送信波と受信波を説明するための図である。図2(a)は、FMCW送信波(送信信号)の波形(数1式参照)を示しており、横軸が時間軸であり縦軸が振幅である。また、図2(b)は、FMCW送信波(送信信号)とFMCW受信波(受信信号)の各々についての瞬時周波数変化を示している。図2(b)の横軸は時間軸であり縦軸は周波数(瞬時周波数)である。なお、図2(a)と図2(b)の時間軸は互いに揃えられている。
図2(b)に示されるように、送信信号(破線)は、周期Tm=1/fmで周波数を変化させた連続波となっている。また、受信信号(実線)は、送信信号から、位相角でφmだけ遅れている。なお、図2(b)においては、数2式で示した受信信号の減衰やドプラシフトを省略している。
数2式で表される受信波形は、超音波振動子を介して受信される信号波形(受信RF信号)である。FMCWドプラでは、受信RF信号に対する復調処理において、FMCW送信波を参照信号として受信波と乗算を行う。図1を利用して説明したように、FM変調器20から出力されるFM連続波が遅延回路25において遅延処理され、参照信号として、ミキサ32には遅延処理されたFM連続波が直接供給され、一方、ミキサ34には遅延処理されたFM連続波がπ/2シフト回路26を経由して供給される。したがって、ミキサ32へ供給される参照信号vrI(t)と、ミキサ34へ供給される参照信号vrQ(t)は、次式のように表現できる。
数3式において、φmrは、遅延回路25における遅延処理により任意に設定できる参照信号の位相を示しており、φ0rは、任意に設定した参照信号の位相に対応して決まる搬送波の位相変化量を示している。
受信ミキサ30では、復調処理として直交検波が行われる。つまり、ミキサ32において、受信RF信号vR(t)と参照信号vrI(t)の乗算に相当する処理が実行され、また、ミキサ34において、受信RF信号vR(t)と参照信号vrQ(t)の乗算に相当する処理が実行される。
ミキサ32における受信RF信号vR(t)と参照信号vrI(t)の乗算vDI(t)は次式のように表現される。なお、次式の計算途中において、周波数2f0の成分が消去されている。これは、例えばローパスフィルタによって除去することができる。
ここで、ベッセル関数に関する次の公式を利用する。
数5式の公式を用いると、数4式はさらに次式のように計算される。
一方、ミキサ34における受信RF信号vR(t)と参照信号vrQ(t)の乗算vDQ(t)は次式のように表現される。なお、次式の計算途中において、周波数2f0の成分が消去されている。これは、例えばローパスフィルタによって除去することができる。
ここで、数6式のvDI(t)と数7式のvDQ(t)とに基づいて、複素ベースバンド信号を定義する。まず、vDI(t)とvDQ(t)に含まれている直流(DC)成分、変調周波数fmの偶数次高調波成分を次式のように表現する。
次に、vDI(t)とvDQ(t)に含まれている変調周波数fmの成分、変調周波数fmの奇数次高調波成分を次式のように表現する。
数8式と数9式から、直交検波後のベースバンド信号において、ドプラシフトfdを含んだドプラ信号は、DC成分と変調周波数fmの成分と変調周波数fmの高調波成分とからなる複数の成分の各々についての両側帯波として出現することがわかる。通信工学ではこの種の信号形式を両側帯波搬送波除去変調(Double-Sideband Suppressed-Carrier, DSB-SC)と呼んでいる。
ここで、受信信号と参照信号の位相を互いに揃えた場合、つまり、遅延回路25における遅延処理によりφmrを調整してφmと一致させた場合(φmr=φm)を考える。φmrとφmを一致させた場合には、数4式におけるkが0となる。この結果を数5式のベッセル関数に適用すると、次式のように、0次のベッセル関数の値のみが1となり、それ以外のベッセル関数の値は0となる。
数10式に示す結果を数8式と数9式に適用すると次式のとおりとなる。
数11式は、参照波(参照信号)の位相φmrを送受信間の位相差φmに設定すると、圧縮変換により、DC成分(直流信号成分)に対応したドプラ信号のみが抽出できることを示している。その結果として得られる複素ドプラ周波数fdの値と極性は、血流などの流体の速度とその極性を表わしている。また、ドプラ信号の振幅は、搬送波および参照波の位相に依存しないこともわかる。
この事実は以下のように解釈することもできる。図1における遅延回路25は、参照波(参照信号)における変調波の位相φmrを、送受信間の変調波の位相差φmに設定する役目を持っている。しかし、この遅延回路25は変調波の位相ばかりでなく、搬送波の位相も同時に変化させる。この値がφ0rである。搬送波の位相は、参照波における変調波の位相φmrに応じて変化するので、送受信間の変調波の位相差に応じた特定の値に定めることができない。しかし、数11式に示されているように、φ0rは、φ0と同様に、どんな値になろうとも、直交検波する限りは、ドプラ信号の振幅、周波数およびその極性に影響を与えない。
そのため、例えば、変調信号(変調波)のみを目標位置の深さに応じた遅延量だけ遅延処理して遅延変調信号を形成し、その遅延変調信号を用いて搬送波信号を変調処理することにより参照信号を形成し、その参照信号とπ/2だけ位相をずらした参照信号とを用いて復調処理を施すようにしてもよい。
そして、本実施形態においては、以下に説明するように、PWドプラ(パルスドプラ)と同様に特定位置のドプラ情報を比較的良好なSNRで得ることができる。数6式から数9式において、ドプラ信号の振幅を支配するJ0(kβ)の因数であるkβについて考察する。数4式におけるkの定義から、kβは次式のように表現できる。
数12式は、kβが深さdに依存して正弦波状に変化することを意味している。
図3は、kβが深さdに依存して正弦波状に変化する様子を示す図である。第1次ベッセル関数の性質により、kβが0のときにJ0(kβ)が最大値となる。図3において実線で示されるkβの波形は、体表からの深さdが正の範囲において0となる深さが3箇所ある。これら3箇所の深さから得られるドプラ信号の振幅が最大となることを意味している。
数12式などから、目的とする深さからの受信信号の位相φmと、参照波の位相φmrとを一致させるとkβを0とすることができ、kβが0となる深さにおいてJ0(kβ)が最大となりドプラ信号の振幅が最大となる。つまり、遅延回路25において、目的とする深さからの受信信号の位相φmと参照波の位相φmrを一致させることにより、目的とする深さからのドプラ信号の振幅が最大となるようにして、そのドプラ信号を選択的に抽出することができる。
以上のように、ドプラ信号が選択的に抽出される目標位置は、遅延回路25における遅延処理に基づいて決定される。図1のシステム制御部50は、目標位置の深さに応じて遅延回路25における遅延時間を制御し、ドプラ情報解析部44において得られるドプラ信号と、そのドプラ信号が得られた位置(深さ)とを対応付ける。その対応関係は、表示処理部46などに伝えられる。
さらに、本実施形態においては、FM連続波の周期性に伴うドプラ周波数の変動成分(ドプラ折り返し電力)の低減処理も実現されている。そこで、その低減処理について以下に詳述する。なお、図1に示した部分(構成)については、以下の説明においても図1の符号を利用する。
まず、ドプラ法の基本原理において、移動体(例えば血流)に関するドプラ周波数(ドプラシフト周波数)は、計測に利用される超音波の周波数と移動体の速度に比例する。本実施形態においては、FM連続波を利用しており、FM連続波は、図2を利用して説明したように、周波数(瞬時周波数)が周期的に変化している。そのため、移動体の速度が一定の場合においても、FM連続波を利用してその移動体のドプラ周波数を計測すると、FM連続波の周期性に伴ってドプラ周波数が周期的に変動する。
図4は、FM連続波の周期性がドプラ周波数へ与える影響を説明する図である。図4には、ドプラシフトの影響を受けていないFM連続波60と、ドプラシフトの影響を受けたFM連続波62が図示されている。なお、図4の横軸は時間軸であり、図4の縦軸にはFM連続波60,62の瞬時周波数が示されている。
本実施形態における超音波の送信信号は、その瞬時周波数がFM連続波60のように周期的に正弦波状に変化する。そのため、移動体の速度が一定の場合においても、ドプラシフトが周期的に変化し、その結果としてFM連続波62のような波形が得られる。つまり、FM連続波60の瞬時周波数が低い(小さい)時刻においては、比較的小さいドプラ周波数fdLとなり、FM連続波60の瞬時周波数が高い(大きい)時刻においては、比較的大きいドプラ周波数fdHとなる。
このように、FM連続波60を利用して得られるドプラ周波数の変動は、FM連続波60の周期性に対応した周期的なものとなる。特に、移動体の速度が大きい場合には、ドプラ周波数fdLとドプラ周波数fdHの差も大きくなり、ドプラ周波数の周期性が比較的顕著になる。一方、移動体の速度が小さい場合にはドプラ周波数fdLとドプラ周波数fdHの差が小さくなり、ドプラ周波数の周期性が比較的目立たなくなる。
図5は、ドプラ周波数の周期性を説明するための図である。図5には、FM連続波を利用して得られるドプラ周波数(ドプラシフト周波数)70の時間変化の様子が示されている。なお、図5の横軸は時間軸であり、図5の縦軸にはドプラ周波数70の瞬時周波数が示されている。
FM連続波を利用して速度が一定の移動体を計測すると、図5に示すようにドプラ周波数70の瞬時周波数が時間的に周期的に変化する。その周期は、FM連続波の周波数変化の周期Tmと同じである。ドプラ周波数70の瞬時周波数が時間的に変化するため、例えば、その時間的な変化を無視して任意の時刻の瞬時周波数を利用して移動体の速度などを計算してしまうと、時刻ごとに結果がばらついてしまい、また、本来観測対象としていないドプラ周波数の高調波成分が出現することとなる。
そこで、本実施形態においては、ドプラ周波数70の周期的な変動成分を低減させることにより、望ましくは、その変動成分を完全に除去することにより、血流などの移動体の速度などを高精度に計測できるようにしている。
図5に示すように、ドプラ周波数70の周期を反転させた信号72とドプラ周波数70とを加算することにより、互いの変動成分が相殺されて時間的に変化しない2fdの値を得ることができる。しかしながら、ドプラ周波数70は、復調信号(ベースバンド信号)に含まれるドプラ信号の瞬時周波数であるため、ドプラ信号とその反転信号とを単純に加え合わせたとしても、図5に示すドプラ周波数70と信号72の加算は実現できない。
本実施形態においては、ドプラ信号を含む復調信号と、その復調信号を周期Tmの半周期だけ遅延させた遅延復調信号(反転復調信号)とを乗算することにより、図5に示すドプラ周波数70と信号72の加算を実現して、ドプラ周波数70の周期的な変動成分を低減している。そこで、ドプラ周波数70の変動成分が低減される原理について詳述する。なお、図1に示した部分(構成)については、以下の説明においても図1の符号を利用する。
本実施形態における超音波の送信信号は、周波数f0のRF波(搬送波)に対して、周波数fmの正弦波によりFM変調を施したFMCW送信波であり、その信号は前述の数1式のとおりである。その送信信号(FMCW送信波)の瞬時角周波数は、数1式の位相項を時間微分して次式のように表現される。
ここで、ドプラシフトを音速(超音波の速度)cと移動体の速度vの比だけ、瞬時周波数が変化する量として定義する。この場合、相対速度vに対するドプラ周波数変化は往復で速度2vとして次式で表現される。
数14式で表現されるドプラ周波数変化を瞬時位相に変換すると次式となる。
数15式で表現される瞬時位相は、移動体からの受信波の瞬時位相に対して、初項で表現される搬送波f0によるドプラシフトに加え、第2項で表現される変調波によるドプラシフトが追加されることを意味している。なお、第3項は積分定数であり、ドプラ周波数の位相を意味する。一般に、血流などの速度計測では、ドプラ周波数の位相情報までは必要としない。また、時間的に変化しない位相成分であるため、速度計測において物理的に大きな意味を含んでいないと考える。
受信波は、送受信時間差(目標位置までの往復の伝播時間)τだけ送信波よりも遅れて到着するため、送受信時間差τを考慮すると、受信波は次式のように表現される。
受信ミキサ30では、送信波に実質的に同じ波形の参照波(参照信号)と受信波との乗算(次式)に相当する処理が実行される。
数17式から2f0の周波数成分をローパスフィルタで除去すると、受信ミキサ30の出力(例えばミキサ32の出力)は次式のように表現することができる。
数18式の結果について、さらに計算を進めると、次式のようになる。
数19式は、ドプラ信号が、新たに定義された変調度β´(数18式参照)と変調周波数fmにより周波数変調されていることを意味している。
図6は、周波数変調されたドプラ信号を説明するための図であり、図6には、数19式に対応したドプラ信号の周波数スペクトラムが示されている。なお、図6の横軸は周波数であり縦軸は電圧である。
図6や数19式に示されるように、変調信号の影響を受けたドプラ信号には、変調周波数fmのゼロ次成分J0(β´)に加え、1次成分J1(β´),2次成分J2(β´),3次成分J3(β´),・・・の折り返し成分が含まれている。
移動体の速度vがゼロに近い場合には、変調度β´(数18式参照)がゼロに近い値となり、ベッセル関数の性質から、ゼロ次成分J0(β´)が1に近い値となり、折り返し成分である1次成分J1(β´),2次成分J2(β´),3次成分J3(β´),・・・は全て0(ゼロ)に近い値となる。つまり、移動体の速度vがゼロに近い場合には、ドプラ信号に含まれる折り返し成分が無視できる程度に小さくなることがわかる。
しかし、移動体の速度vが比較的大きい場合には、ドプラ信号に含まれる折り返し成分を単純には無視できない。そこで、本実施形態においては、次のようにしてドプラ信号に含まれる折り返し成分を低減し、望ましくは折り返し成分を除去している。
図1の受信ミキサ30の出力(例えばミキサ32の出力)に含まれるドプラ信号、つまり、復調信号に含まれるベースバンド領域におけるドプラ信号は、前述の数18式のように表現される。受信ミキサ30の後段には、半周期遅延回路35,37が設けられており、半周期遅延回路35,37の各々において、周期Tmの半周期に対応した遅延量だけ、復調信号に対して遅延処理が実行される。
そこで、数18式のvd(t)をvd1(t)とし、数18式のvd(t)に対して、周期Tmの半周期に対応した遅延量だけ遅延処理した信号をvd2(t)とすると、次式が得られる。
例えば、ミキサ32から数20式のvd1(t)が出力され、半周期遅延回路35から数20式のvd2(t)が出力される。そして、ミキサ36において、vd1(t)とvd2(t)の乗算に相当する処理が実行される。vd1(t)とvd2(t)の乗算は次式のように表現される。
数21式の結果である最終3行においては、第1項のドプラ信号が、第2項および第3項に含まれる変調に伴う成分J0(2β´),J1(2β´),J2(2β´),J3(2β´),・・・から独立している。また、第1項のドプラ信号内において、ドプラ周波数fdが2倍に逓倍されている。ドプラ周波数fdが2倍に逓倍される結果は、図5におけるドプラ周波数70と信号72の加算に対応している。図5に示したように、ドプラ周波数70と信号72の加算により、ドプラ周波数70の変動成分である折り返し成分が除去される。
図7は、折り返し成分を除去したドプラ信号を説明するための図である。図7には、数21式のvd3(t)に関する周波数スペクトラムが示されている。なお、図7の横軸は周波数であり縦軸は電圧である。
図7の周波数スペクトラム内において、2倍に逓倍されたドプラ周波数2fdに対応したドプラ信号80が、J0(2β´),J1(2β´),J2(2β´),J3(2β´),・・・から独立している。例えば、ミキサ36から得られる信号がFFT処理部42において周波数スペクトラムに変換されることにより、図7の周波数スペクトラムが得られる。そして、ドプラ情報解析部44において、図7の周波数スペクトラムからドプラ周波数2fdに対応したドプラ信号80が抽出される。これにより、周波数変調に伴うドプラ周波数の変動成分である折り返し成分が除去されて、時間的に変化しない2fdの値を得ることが可能になる。
なお、半周期遅延回路35とミキサ36における処理と同様に、半周期遅延回路37とミキサ38においても、復調信号と遅延された復調信号の乗算に相当する処理が実行され、周波数変調に伴うドプラ周波数の変動成分である折り返し成分が除去される。
このように、本実施形態においては、FM連続波の周期性に伴うドプラ周波数の変動成分(ドプラ折り返し電力)の低減処理が実現され、ドプラ折り返し電力が低減されて、望ましくはドプラ折り返し電力が完全に除去されて、ドプラ情報として、時間的に変化しない2fdの値を得ることが可能になる。
以上、本発明の好適な実施形態を説明したが、上述した本発明の好適な実施形態等は、あらゆる点で単なる例示にすぎず、本発明の範囲を限定するものではない。本発明は、その本質を逸脱しない範囲で各種の変形形態を包含する。
例えば、上述した実施形態においては、周波数を周期的に変化させる連続波の変調送信信号を形成するにあたって、搬送波信号(RF波発振器22から供給されるRF波)に対して周波数変調処理を施している。この周波数変調処理に換えて、周波数変調処理と同じ角度変調の方式として当業者において明らかな位相変調処理(PM処理)を利用してもよい。つまり、搬送波信号(RF波発振器22から供給されるRF波)に対して位相変調処理を施すことにより、FM変調器20から出力されるFM連続波と同じ波形あるいは同等な波形を形成してもよい。なお、周波数を周期的に変化させる連続波のデータをメモリなどに記憶しておき、このメモリから読み出されるデータに基づいて、当該連続波を生成してもよい。
20 FM変調器、22 RF波発振器、24 変調波生成部、25 遅延回路、35,37 半周期遅延回路、36,38 ミキサ、42 FFT処理部、44 ドプラ情報解析部。
Claims (6)
- 周波数を周期的に変化させた連続波の送信信号を出力する送信信号処理部と、
前記送信信号に基づいて超音波を生体に送波して生体からの反射波を受波することにより得られる受信信号に対して、前記送信信号に実質的に等しい波形の参照信号を用いて復調処理を施す受信信号処理部と、
復調処理により得られた信号からドプラ信号を抽出するドプラ信号抽出部と、
を有し、
前記受信信号処理部は、生体内の目標位置から得られる受信信号との間の相関関係が調整された参照信号を用いて前記復調処理を施すことにより、前記目標位置に対応した復調信号を選択的に抽出し、
前記ドプラ信号抽出部は、復調信号と当該復調信号に対して周期的に反転関係にある反転復調信号とに基づいて、前記連続波の周期性に伴うドプラ周波数の変動成分を低減し、前記目標位置に対応したドプラ信号を抽出する、
ことを特徴とする超音波診断装置。 - 請求項1に記載の超音波診断装置において、
前記送信信号処理部は、周期Tmで周波数を変化させた連続波の送信信号を出力し、
前記ドプラ信号抽出部は、復調信号に対して前記周期Tmの半周期に対応した遅延量だけ遅延処理を施すことにより反転復調信号を得る、
ことを特徴とする超音波診断装置。 - 請求項2に記載の超音波診断装置において、
前記ドプラ信号抽出部は、復調信号と反転復調信号との乗算結果に基づいて、ドプラ周波数の前記変動成分を低減する、
ことを特徴とする超音波診断装置。 - 請求項3に記載の超音波診断装置において、
前記ドプラ信号抽出部は、復調信号と反転復調信号を乗算することにより、復調信号に含まれるドプラ周波数と反転復調信号に含まれるドプラ周波数とを加算して前記変動成分を相殺する、
ことを特徴とする超音波診断装置。 - 請求項4に記載の超音波診断装置において、
前記ドプラ信号抽出部は、復調信号と反転復調信号を乗算することにより得られる乗算信号から、前記加算されたドプラ周波数に対応した周波数成分を前記ドプラ信号として抽出する、
ことを特徴とする超音波診断装置。 - 請求項1から5のいずれか1項に記載の超音波診断装置において、
前記目標位置からの受信信号の位相と前記参照信号の位相が等しくなるように、目標位置の深さに応じた遅延量だけ前記参照信号を遅延処理することにより、目標位置から得られる受信信号と参照信号との間の相関を強める、
ことを特徴とする超音波診断装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008286528A JP2010110503A (ja) | 2008-11-07 | 2008-11-07 | 超音波診断装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2008286528A JP2010110503A (ja) | 2008-11-07 | 2008-11-07 | 超音波診断装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2010110503A true JP2010110503A (ja) | 2010-05-20 |
Family
ID=42299543
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2008286528A Pending JP2010110503A (ja) | 2008-11-07 | 2008-11-07 | 超音波診断装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2010110503A (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012070874A (ja) * | 2010-09-28 | 2012-04-12 | Kanazawa Inst Of Technology | 超音波診断装置 |
WO2012099341A2 (ko) * | 2011-01-21 | 2012-07-26 | 알피니언메디칼시스템 주식회사 | 의료진단장치 |
KR101310930B1 (ko) | 2011-01-21 | 2013-09-25 | 알피니언메디칼시스템 주식회사 | 의료진단장치 |
-
2008
- 2008-11-07 JP JP2008286528A patent/JP2010110503A/ja active Pending
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012070874A (ja) * | 2010-09-28 | 2012-04-12 | Kanazawa Inst Of Technology | 超音波診断装置 |
WO2012099341A2 (ko) * | 2011-01-21 | 2012-07-26 | 알피니언메디칼시스템 주식회사 | 의료진단장치 |
WO2012099341A3 (ko) * | 2011-01-21 | 2012-10-11 | 알피니언메디칼시스템 주식회사 | 의료진단장치 |
KR101310930B1 (ko) | 2011-01-21 | 2013-09-25 | 알피니언메디칼시스템 주식회사 | 의료진단장치 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8323200B2 (en) | Ultrasound diagnostic apparatus | |
US8100832B2 (en) | Ultrasound diagnostic apparatus | |
JP5459963B2 (ja) | 超音波診断装置 | |
JP5654198B2 (ja) | 超音波診断装置 | |
JP2006288974A (ja) | 超音波診断装置 | |
JP2010110503A (ja) | 超音波診断装置 | |
JP2011229557A (ja) | 超音波診断装置 | |
JP2010012160A (ja) | 超音波診断装置 | |
JP2011234846A (ja) | 超音波診断装置 | |
JP2011036448A (ja) | 超音波診断装置 | |
JP5235110B2 (ja) | 超音波診断装置 | |
JP2010279531A (ja) | 超音波診断装置 | |
JP2011240006A (ja) | 超音波診断装置 | |
JP5291952B2 (ja) | 超音波診断装置 | |
JP2009261749A (ja) | 超音波診断装置 | |
JP2010125246A (ja) | 超音波診断装置 | |
JP2011098162A (ja) | 超音波診断装置 | |
JP2010188115A (ja) | 超音波診断装置 | |
JP2009254526A (ja) | 超音波診断装置 | |
JP2011104241A (ja) | 超音波診断装置 | |
JP5297082B2 (ja) | 超音波診断装置 | |
JP2013017731A (ja) | 超音波診断装置 | |
JP2010022737A (ja) | 超音波診断装置 | |
JP2011036506A (ja) | 超音波診断装置 | |
JP2010178793A (ja) | 超音波診断装置 |