(57)【特許請求の範囲】
【請求項１】
超音波パルスを被検体内に送信し、反射波を受信することにより受信信号を得る送受信手段と、
レンジゲートマーカの深さ及び長さに応じた特定の期間に前記送受信手段から出力される前記受信信号を積分することにより積分信号を得る手段と、
前記レンジゲートマーカ内のドプラスペクトラムを得るために前記積分信号を周波数解析する手段と、
前記送受信手段と前記積分手段との間に設けられ、前記積分手段による積分期間をΔＴRGとして１／ΔＴRG１／２に応じた重み係数を前記受信信号に乗算する手段とを具備したことを特徴とする超音波診断装置。
【請求項２】
超音波を連続波として被検体内に送信し、反射波を受信することにより受信信号を得る手段と、
前記受信信号の特定の周波数帯域を減衰するフィルタ手段と、
前記フィルタ手段の出力帯域に対する前記フィルタ手段の入力帯域の比に応じた重み係数を前記フィルタ手段の出力信号に乗算する手段と、
ドプラスペクトラムを得るために前記フィルタ手段の出力信号を周波数解析する手段とを具備したことを特徴とする超音波診断装置。
【請求項３】
超音波を被検体内に送信し、反射波を受信することにより受信信号を得る手段と、
ドプラスペクトラムを得るために前記受信信号をデータ数で周波数解析する手段と、
前記データ数をＮとして重み係数１／Ｎ１／２を前記受信信号に乗算する手段とを具備したことを特徴とする超音波診断装置。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、受信信号を周波数解析することによりドプラスペクトラムを得る超音波診断装置に関する。
【０００２】
【従来の技術】
超音波は医学的な見地から様々に応用されているが、その主流は超音波パルス反射法を用いて生体の軟部組織の断層像を得る超音波診断装置である。この超音波診断装置は無侵襲検査法で、組織の断層像を表示するものであり、Ｘ線診断装置、Ｘ線ＣＴ装置、ＭＲＩおよび核医学診断装置などの他の診断装置に比べて、リアルタイム表示が可能、装置が小型で安価、Ｘ線などの被曝がなく安全性が高い、ドプラ効果を利用してドプラスペクトラムや血流イメージングが可能である等の独占的な特徴を有している。ここでは、ドプラスペクトラムに着目する。ドプラスペクトラムの表示としては、縦軸にドプラ周波数（速度）、横軸に時間、輝度を各ドプラ周波数成分の強さ（パワー）でＭモードのようにスクロール表示が使用される。
【０００３】
このような超音波診断装置のドプラスペクトラム機能において、感度向上を目的として、バースト波が実用化されている。バースト波とは複数の超音波パルスの系列をいう。バースト波の時間幅をバースト長という。バースト波の場合、パルス波と違って、レンジゲートマーカ内のドプラ信号の積分処理に工夫が必要である。図７（ａ）に示す従前のケースでは、レンジゲートマーカの深さ及び長さに相当する期間ｔ1 〜ｔ2 に受信されるドプラ信号だけを積分処理に取り込む。なお、ｔ1 は、バースト波の送信時刻ｔ0 から、レンジゲートマーカの深さに相当する時間（往復を考慮してレンジゲートマーカの深さの２倍の距離を超音波が伝搬するのに要する時間）を経過した時刻であり、ｔ2 はｔ1 から、レンジゲートマーカの長さに相当する時間（レンジゲートマーカの長さの２倍の距離を超音波が伝搬するのに要する時間）を経過した時刻である。この従前のケースでは、バースト波のバースト長τを考慮していない、つまり同一反射源からの反射波がバースト長τの期間、継続的に受信されることを考慮していないので、積分信号の感度分布を見て理解されるように、レンジゲートマーカよりτに相当する浅い領域からのドプラ信号が積分処理に取り込まれてしまう。
【０００４】
この問題点を解決するように工夫したケースを図７（ｂ）に示す。このケースでは、バースト波のバースト長τは、レンジゲートマーカの長さに相当する時間の半分に意図的に調整される。そして、レンジゲートマーカの深さ及び長さに相当する期間ｔ1 〜ｔ2 の後半の期間ｔ3 〜ｔ2 に受信されるドプラ信号だけを積分処理に取り込む。結果的に、積分期間とバースト長τとは一致する。このケースでは、積分信号の感度分布を見て理解されるように、レンジゲートマーカ外のドプラ信号は積分処理に取り込まれない。
【０００５】
しかし、図７（ｂ）のケースには次のような問題点がある。距離分解能の著しい低下を回避するために、バースト波のバースト長τは制限される。したがって、レンジゲートマーカが非常に長く設定されたとき、図７（ｃ）に示すように、積分期間が最長のバースト長τmax より長くなってしまうケースが起こり得る。この場合、レンジゲートマーカ内にも関わらず、そこからのドプラ信号が積分処理に全く取り込まれない、いわゆる不感領域が存在してしまう。
【０００６】
また、ドプラスペクトラム機能には、次のような問題もある。パルス波ドプラ（PWドプラ）では、レンジゲートマーカの長さを変化させると、積分期間が変化するが、この積分期間の変化は、実質的にゲインが変化することを意味する。また、連続波ドプラ（CWドプラ）では、ドプラ信号を高域通過型フィルタ（HPF ）にかけるが、レート周波数PRF を変化させると、この高域通過型フィルタの入力帯域に対する出力帯域の比が変化するが、この入出力帯域比の変化は、実質的にゲインが変化することを意味する。さらに、PWドプラ、CWドプラの両者に言えることとして、周波数解析としての高速フーリエ変換（FFT ）のデータ数の変化は、実質的にゲインが変化することを意味する。したがって、ゲインの変化に伴って、ノイズレベルも変化してしまうという問題があった。
【０００７】
また、PWドプラでは、実質的にレート周波数を低くして、低速検出能を向上させるために、交互スキャンと呼ばれる機能が実用化されている。通常、ある１本のラスタに対して超音波送受信動作をレート周波数で所定回数繰り返しながら、順次次のラスタに移っていくというスキャンが行われる。これに対し、交互スキャンでは、ある１本のラスタに対して超音波送受信動作を１回行い、同じラスタに対する次回の超音波送受信動作までに間に、他の何本かのラスタに対して１回ずつ超音波送受信動作を行うものであり、１本のラスタに関して見ると、超音波送受信動作の間隔が長くなり、実質的にレート周波数が低くなるというものである。なお、同じラスタに対する前後２回の超音波送受信動作の間隙に行う超音波送受信動作の回数を交互段数として定義する。交互段数の増加は、周波数解析のデータ数が一定であれば、データ取り込み時間の延長を意味する。データ取り込み時間の延長は、ドプラスペクトラムの表示上で周波数分解能と時間分解能との比が変化することにつながる。これは好ましくないので、現在の多くの装置には、交互段数が増加すると、それに応じて周波数解析のデータ数を減少させて、周波数分解能と時間分解能との比を一定に保つ機能が装備されている。しかし、周波数解析のデータ数の減少は、Ｓ／Ｎを劣化（低下）させてしまう。
【０００８】
【発明が解決しようとする課題】
本発明の目的は、パルス波ドプラでは積分期間の変化、連続波ドプラでは高域通過型フィルタの入力帯域に対する出力帯域の比の変化、周波数解析としての高速フーリエ変換（FFT ）のデータ数の変化等により実質的にゲインが変化しても、ノイズレベルは変化しない超音波診断装置を提供することである。
【０００９】
【課題を解決するための手段】
本発明による超音波診断装置は、超音波パルスを被検体内に送信し、反射波を受信することにより受信信号を得る送受信手段と、レンジゲートマーカの深さ及び長さに応じた特定の期間に前記送受信手段から出力される前記受信信号を積分することにより積分信号を得る手段と、前記レンジゲートマーカ内のドプラスペクトラムを得るために前記積分信号を周波数解析する手段と、前記送受信手段と前記積分手段との間に設けられ、前記積分手段による積分期間をΔＴRGとして１／ΔＴRG１／２に応じた重み係数を前記受信信号に乗算する手段とを具備する。
本発明によれば、積分期間が変化すると、それに伴って実質的にゲインが変化してノイズレベルも変化してしまうが、積分期間に応じた重み係数を受信信号に乗算するので、積分期間の変化に伴う実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を連続波として被検体内に送信し、反射波を受信することにより受信信号を得る手段と、前記受信信号の特定の周波数帯域を減衰するフィルタ手段と、前記フィルタ手段の出力帯域に対する前記フィルタ手段の入力帯域の比に応じた重み係数を前記フィルタ手段の出力信号に乗算する手段と、ドプラスペクトラムを得るために前記フィルタ手段の出力信号を周波数解析する手段とを具備する。
本発明によれば、例えばレート周波数が変化すると、それに伴ってフィルタ手段の出力帯域に対する入力帯域の比が変化して実質的にゲインが変化してノイズレベルも変化してしまうが、この比に応じた重み係数をフィルタ手段の出力信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を被検体内に送信し、反射波を受信することにより受信信号を得る手段と、ドプラスペクトラムを得るために前記受信信号をデータ数で周波数解析する手段と、前記データ数をＮとして重み係数１／Ｎ１／２を前記受信信号に乗算する手段とを具備する。
本発明によれば、周波数解析のデータ数が変化すると実質的にゲインが変化してノイズレベルも変化してしまうが、このデータ数に応じた重み係数を受信信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
【００１９】
【実施の形態】
以下に、本発明の実施の形態を図面を参照して説明する。
（第１の実施の形態）
図１は第１の実施の形態による超音波診断装置のブロック図である。本実施の形態による超音波診断装置は、超音波プローブ１、送受信系（Ｔ／Ｒ）２、スペクトラムドプラユニット３、ディジタルスキャンコンバータ（ＤＳＣ）４、モニタ５、コンソール６とから構成される。
【００２０】
超音波プローブ１は、圧電セラミック等の機械運動／電気信号の可逆的変換素子としての複数の圧電素子が一次元に配列された圧電素子アレイをその先端に装備してなる。
【００２１】
送受信系２は、送信系と受信系とを有する。送信系は、パルス発生器、送信遅延回路、パルサとからなる。パルス発生器は例えば５ＫＨｚのレート周波数ｆr （周期；１／ｆr 秒）でレートパルスを繰り返し発生する。このレートパルスはチャンネル数に分配され、送信遅延回路に送られる。送信遅延回路は、超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間を各レートパルスに与える。パルサは、送信遅延回路からレートパルスを受けた時刻で超音波プローブ１の圧電素子にチャンネル毎に送信電圧を印加する。これにより超音波ビームが被検体に送信される。
【００２２】
被検体内の音響インピーダンスの不連続面（反射源）で反射した反射波は超音波プローブ１で受信され、電気信号に変換される。超音波プローブ１から素子毎に出力される電気信号は、チャンネル毎に纏められて受信系に取り込まれる。受信系は、プリアンプ、受信遅延回路、加算器、位相検波回路から構成される。電気信号は、プリアンプで増幅され、受信遅延回路により受信指向性を決定するのに必要な遅延時間を与えられ、そして加算器で加算される。この加算により受信指向性に応じた方向からの反射成分が強調された受信信号が得られる。この受信信号は位相検波回路で位相検波される。これによりドプラ効果による偏移周波数のドプラ信号が得られる。ドプラ信号はスペクトラムドプラユニット３に供給される。
【００２３】
スペクトラムドプラユニット３は、特定の期間に供給されるドプラ信号のみをレンジゲートマーカユニット８で取り出し、そして比較的動きの速い高周波の血流成分と比較的動きの遅い血管壁や弁等の低周波のクラッタ成分とを高域通過型フィルタ９で分離し、続いて高域通過型フィルタ９を通過したドプラ信号を高速フーリエ変換処理回路（FFT ）１０でリアルタイムで周波数解析することにより、各周波数成分の強さ（パワー）を表すドプラスペクトラムデータを生成する。レンジゲートマーカユニット８でドプラ信号を取り出す特定の期間は、コンソール６を介してＢモード画像上に設定されるレンジゲートマーカマーカの深さ及び長さに基づいてコントローラ７で制御される。
【００２４】
高速フーリエ変換処理回路１０で得られたドプラスペクトラムデータは、ディジタルスキャンコンバータ４、図示しないディジタルアナログコンバータを介してモニタ５にドプラスペクトラムとして表示される。ドプラスペクトラムの表示としては、縦軸にドプラ周波数（速度）、横軸に時間、輝度を各ドプラ周波数成分の強さ（パワー）として、Ｍモードのようにスクロール表示が使用される。
【００２５】
レンジゲートマーカユニット８は、ドプラ信号を乗算器（増幅器）１１で振幅変調し、振幅変調されたドプラ信号を特定の期間だけ積分器１２で積分する。乗算器１１による振幅変調データは、関数発生ＲＯＭ１３から供給される。積分器１２が積分を実行する特定の期間は、ＲＧ幅カウンタ１４及びＲＧスタート位置カウンタ１５により調整される。
【００２６】
ＲＧスタート位置カウンタ１５は、コントローラ７からレンジゲートマーカの深さに対応するRG位置データの供給を受ける。ＲＧスタート位置カウンタ１５は、スキャンの開始からクロックパルスの計数（アップカウント）を開始し、計数値がRG位置データが表す値に達した時、スタートパルスを出力する。
【００２７】
ＲＧ幅カウンタ１４は、コントローラ７からレンジゲートマーカの長さに対応するRG幅データの供給を受ける。ＲＧ幅カウンタ１４は、ＲＧスタート位置カウンタ１５からスタートパルスを受けた時刻でクロックパルスの計数（アップカウント）を開始し、計数値がRG幅データが表す値に達する時まで計数動作を継続する。したがって、ＲＧ幅カウンタ１４からは、レンジゲートマーカの深さ及び長さに対応する特定の期間だけ、計数信号が出力される。
【００２８】
積分器１２は、ＲＧ幅カウンタ１４から計数信号が出力されている期間に限り積分を実行する。これにより、特定の期間内の振幅変調されたドプラ信号だけが積分信号として取り出される。
【００２９】
このようにコントローラ７からのRG位置データ及びRG幅データに、レンジゲートユニット８で取り出されるドプラ信号の期間（積分期間）が依存される。換言するとコントローラ７はRG位置データ及びRG幅データを調整することにより積分期間の開始時刻及びその時間幅を自由に変えることができることが理解されよう。
【００３０】
関数発生ＲＯＭ１３は、ＲＧ幅カウンタ１４から計数信号が出力されている特定の期間だけ係数を乗算器１１に出力する。ドプラ信号は乗算器１１で関数発生ＲＯＭ１３からの係数に応じて振幅変調される。
【００３１】
関数発生ＲＯＭ１３には、予め複数種類の重み係数が記憶されている。関数発生ＲＯＭ１３には、ＲＧ幅カウンタ１４からの計数値信号と、コントローラ７からの関数選択信号とがアドレス信号として供給される。アドレス信号に応じたアドレス部の係数が関数発生ＲＯＭ１３から読み出される。コントローラ７は複数種類の関数選択信号の中から、コンソール６を介して選択された感度特性に対応する１つの関数選択信号を選択的に関数発生ＲＯＭ１３に供給する。ここでいう関数とは関数発生ＲＯＭ１３から読み出される係数の時間変化を意味する。第１の関数選択信号のもとでは、関数発生ＲＯＭ１３からは一定の係数（＝1.0 ）が読み出され、つまり振幅変調処理の無いが如く振る舞われる。第２の関数選択信号のもとでは、関数発生ＲＯＭ１３からは、例えばハミング窓関数を満足するように経時的に変化する係数が読み出される。
【００３２】
次に本実施の形態の動作を説明する。ここで、以下で使用される用語を定義する。レンジゲートマーカの深さに相当する時間とは、超音波送受信の往復を考慮して、レンジゲートマーカの深さの２倍の距離を超音波が伝搬するのに要する時間として定義する。また、レンジゲートマーカの長さに相当する時間とは、レンジゲートマーカの長さの２倍の距離を超音波が伝搬するのに要する時間として定義する。
【００３３】
コントローラ７は、第１のモードと第２のモードとを有する。第１のモードはレンジゲートマーカの長さに相当する時間の半分と、バースト波のバースト長τとが一致する場合に適用され、第２のモードはレンジゲートマーカの長さに相当する時間の半分より、バースト波のバースト長τ（最長のバースト長τmax ）が短い場合に適用される。なお、図７（ｂ）を参照して説明した通り、バースト波のバースト長τは、レンジゲートマーカの長さに相当する時間の半分に積極的に調整されるのではあるが、距離分解能の著しい低下を抑える目的であえてバースト長τに制限を与えている由により、設定されるレンジゲートマーカによっては、レンジゲートマーカの長さに相当する時間の半分より、バースト波のバースト長τ（最長のバースト長τmax ）が短い場合が発生する。
【００３４】
第１のモードでは、図７（ｂ）を参照して説明した通り、バースト波のバースト長τはレンジゲートマーカの長さに相当する時間の半分に意図的に調整される。時刻ｔ3 〜ｔ2 の期間に供給されるドプラ信号だけがレンジゲートユニット８で積分処理に取り込まれる。時刻ｔ3 は、時刻ｔ1 からバースト長τだけ経過した時刻であり、時刻ｔ1 はバースト波の送信時刻ｔ0 からレンジゲートマーカの深さに相当する時間△Ｔd を経過した時刻である。時刻ｔ2 は時刻ｔ1 からレンジゲートマーカの長さに相当する時間△ＴRGを経過した時刻、つまりレンジゲートマーカの最深部に相当する時刻である。積分期間は、レンジゲートマーカの長さに相当する時間△ＴRGの半分（△ＴRG／２）に設定される。第１のモードでは、積分信号の感度分布を見て理解されるように、レンジゲートマーカ外のドプラ信号は好ましく全く積分処理に取り込まれない。
【００３５】
図２（ａ），（ｂ）は第２のモードでの積分開始時刻及び積分期間の時間幅を示すタイムチャートであり、同図（ａ）は第１の関数選択信号に対応し、また同図（ｂ）は第２の関数選択信号に対応する。
【００３６】
第２のモードでは、図２（ａ），（ｂ）に示すように、時刻ｔ1 ´〜ｔ2 ´の積分期間に供給されるドプラ信号だけがレンジゲートユニット８で積分処理に取り込まれる。時刻ｔ1 ´は、時刻ｔ1 からバースト長τの半分（τ／２）だけ経過した時刻であり、時刻ｔ1 はバースト波の送信時刻ｔ0 からレンジゲートマーカの深さに相当する時間△Ｔd を経過した時刻である。時刻ｔ2 ´は時刻ｔ1 ´からレンジゲートマーカの長さに相当する時間△ＴRGを経過した時刻である。積分期間は、レンジゲートマーカの長さに相当する時間△ＴRGに設定される。第２のモードでは、積分信号の感度分布を見て理解されるように、レンジゲートマーカ内に不感領域は存在しない。
【００３７】
このように第２のモードでは、レンジゲートマーカの深さに相当する時刻ｔ1 からτ／２だけ遅延した時刻ｔ1 ´で積分を開始し、レンジゲートマーカの長さに相当する時間△ＴRG、当該積分を継続することにより、レンジゲートマーカより浅い側及び深い側それぞれにτ／２に相当する範囲からのドプラ信号が積分処理に取り込まれてしまうが、図７（ａ）のケースと比較すると、そのレンジゲートマーカ外のドプラ信号の取り込み量（面積）は絶対的に減少する。また、図７（ｃ）と比較して理解されるように、積分期間がバースト長τより長くなるケースに直面しても、レンジゲートマーカ内の不感領域を完全に排除できる。
【００３８】
さらに図２（ｂ）の積分信号の感度分布から理解されるように、関数発生ＲＯＭ１３から例えばハミング窓関数を満足するように経時的に変化する係数が読み出され、ドプラ信号に乗算される場合、感度分布が中央付近で高くなり、レンジゲートマーカより浅い側及び深い側それぞれにτ／２に相当する範囲の感度分布が抑制される。
（第２の実施の形態）
図３は第２の実施の形態による超音波診断装置のブロック図である。本実施の形態による超音波診断装置は、超音波プローブ２１、送受信系（Ｔ／Ｒ）２２、スペクトラムドプラユニット２３、ディジタルスキャンコンバータ（ＤＳＣ）２４、モニタ２５、コンソール２６とから構成される。
【００３９】
超音波プローブ２１は、第１の実施の形態の超音波プローブ１に同じである。送受信系２２は、第１の実施の形態の送受信系２が、超音波をパルス波（PW）又は連続波（CW）として選択的に発生させることが可能に改良され、アナログのドプラ信号をドプラデータにディジタル化して出力するように最後段にＡＤＣを追加される。広域通過型フィルタ３０は、第１の実施の形態の広域通過型フィルタ９に同じである。高速フーリエ変換処理回路３１は、第１の実施の形態の高速フーリエ変換処理回路１０に同じである。ディジタルスキャンコンバータ２４は、第１の実施の形態のディジタルスキャンコンバータ４に同じである。モニタ２５は、第１の実施の形態のモニタ５に同じである。
【００４０】
PW時には、送受信系２２からのドプラ信号はＲＧ積分回路２８に供給される。CW時には、送受信系２２からのドプラ信号はCW帯域フィルタ２９に供給される。図４（ａ）に示すように、ＲＧ積分回路２８は、データ数Ｎのドプラデータを、加算器３３1 と、加算器３３1 への帰還経路に挿入された遅延回路３３2 とで積分する。勿論、遅延回路３３2 の遅延時間は、ＡＤＣのサンプリング周期△ＴS に設定される。データ数Ｎは、レンジゲートマーカの長さに応じた積分期間△ＴRGに依存し、△ＴRG／△ＴS によりコントローラ２７で決定される。加算器３３1 の前段には、ドプラデータに、積分期間△ＴRGに応じた係数を乗算するための乗算器３２が設けられる。係数は、データ数Ｎに基づいて、１／Ｎ1/2 でコントローラ２７により与えられる。
【００４１】
ドプラデータをデータ数Ｎで加算することにより、信号レベルはＮ倍になり、一方、ノイズレベルはＮ1/2 倍になる。したがって、レンジゲートマーカの長さ、つまり積分期間、すなわち加算データ数の変動に依存して、ノイズレベルが変動してしまう。本実施の形態では、ドプラデータに係数１／Ｎ1/2 を乗算することにより、レンジゲートマーカの長さの変動、つまり積分期間の変動、すなわち加算データ数の変動によらず、ノイズレベルを一定に保つことができる。
【００４２】
図４（ｂ）に示すように、CW帯域フィルタ２９は、ドプラデータの周波数帯域を入力帯域ｆinから出力帯域ｆout に制限するために、 FIR型の低域通過フィルタ（FIR LPF ）３４と、 IIR型の低域通過フィルタ（IIR LPF ）３５とが縦列に接続される。 IIR型の低域通過フィルタ３５の出力側には、濾過されたドプラデータに係数を乗算するための乗算器３６が設けられる。この係数は、入力帯域ｆinと出力帯域ｆout との比に基づいて、（ｆin／ｆout ）1/2 でコントローラ２７により与えられる。
【００４３】
入力帯域ｆinが出力帯域ｆout に制限されることは、１／（ｆin／ｆout ）1/2 ＜１のゲインがかかることに等価的であり、これは入力帯域ｆinと出力帯域ｆout との比に応じてノイズレベルが変動してしまうことを意味する。本実施の形態では、 FIR型の低域通過フィルタ３４及び IIR型の低域通過フィルタ３５の出力データに係数（ｆin／ｆout ）1/2 を乗算することにより、入力帯域ｆinと出力帯域ｆout との比の変動に関わらず、ノイズレベルを一定に保つことができる。
【００４４】
図４（ｃ）に示すように、高速フーリエ変換処理回路３１は、高速フーリエ変換処理部３８でレンジゲートマーカの長さ、つまり積分期間に応じて、ドプラデータをデータ数Ｎで高速フーリエ変換処理に供する。高速フーリエ変換処理部３８の前段には、ドプラデータに、係数を乗算するための乗算器３７が設けられる。係数は、データ数Ｎに基づいて、１／Ｎ1/2 でコントローラ２７により与えられる。
【００４５】
ドプラデータをデータ数Ｎで高速フーリエ変換処理に供することにより、ノイズレベルに実質的にＮ1/2 のゲインがかかる。したがって、データ数の変動に依存して、ノイズレベルが変動してしまう。本実施の形態では、ドプラデータに係数１／Ｎ1/2 を乗算することにより、データ数の変動によらず、ノイズレベルを一定に保つことができる。
（第３の実施の形態）
図５は第３の実施の形態による超音波診断装置のブロック図である。なお、超音波プローブ４１は、第１の実施の形態の超音波プローブ１に同じである。高速フーリエ変換処理回路５０は、第１の実施の形態の高速フーリエ変換処理回路１０に同じである。ディジタルスキャンコンバータ４４は、第１の実施の形態のディジタルスキャンコンバータ４に同じである。モニタ４５は、第１の実施の形態のモニタ５に同じである。
【００４６】
本実施の形態による超音波診断装置は、超音波プローブ４１、送受信系４２、スペクトラムドプラユニット４３、ディジタルスキャンコンバータ（ＤＳＣ）４４、モニタ４５、システムコントローラ４６、ドプラ用コントローラ４７、コンソール４８とから構成される。送受信系４２は、超音波プローブ４１を駆動して交互スキャンを実行することが可能に構成されている。交互スキャンとは、通常スキャンが、ある１本のラスタに対して超音波送受信動作をレート周波数で所定回数繰り返しながら、順次次のラスタに移っていくのに対して、ある１本のラスタに対して超音波送受信動作を１回行い、同じラスタに対する次回の超音波送受信動作までに間に、他の何本かのラスタに対して１回ずつ超音波送受信動作を行うものであり、１本のラスタに関して見ると、超音波送受信動作の間隔が長くなり、実質的にレート周波数が低くなって低速検出能が向上するというものである。なお、同じラスタに対する前後２回の超音波送受信動作の間に行う他のラスタの超音波送受信動作の回数を交互段数として定義する。また、送受信系４２は、システムコントローラ４６の制御にしたがって交互段数を変化させることが可能に構成されている。さらに、送受信系４２は、ドプラ用コントローラ４７の制御にしたがって送信電圧を変化させることが可能に構成されている。
【００４７】
スペクトラムドプラユニット４３は、ドプラゲイン補正用の乗算器４９と、高速フーリエ変換処理回路５０とから構成される。なお、送受信系４２の出力がディジタルでなくアナログであれば、乗算器４９は増幅器に置換される。
【００４８】
システムコントローラ４６には、コンソール４８が接続される。コンソール４８を介して、交互段数を変化させることが可能である。
図６に示すように、ドプラ用コントローラ４７は、データ数調整用コントローラ５１と、ドプラ感度補正用コントローラ５２とから構成されている。データ数調整用コントローラ５１は、交互段数が変化してレート周波数ｆPRF が実質的に変化しても、ドプラスペクトラムの表示上で周波数分解能と時間分解能との比を一定にするように、高速フーリエ変換処理回路５０の周波数解析のデータ数を調整するためのものである。この調整方法は周知であるので、簡単に説明する。交互段数をＫ、レート周波数をｆPRF 、スクロールスピードをＳ、１ピクセルに相当する周波数幅に対する時間幅の比を表すピクセル比をＰn 、係数をαとする。まず、パラメータＡA を（１）式にしたがって求める。
【００４９】
ＡA ＝α・（（ｆPRF ・Ｓ）／（Ｐn ・Ｋ））  …（１）
ここで、データ数をＮ´とする。Ｎ´＝Ｗw ・Ｎで与えられる。Ｎは、基準データ数（通常、最小データ数とされる）であり、Ｗw は基準データ数に対するデータ数の比（Ｎ´／Ｎ）である。選択可能なデータ数が、１６、３２、６４、１２８、２５６であるとすると、例えば、Ｎ＝１６であり、この場合、Ｗw ＝１、２、４、８又は１６である。
【００５０】
そして、（２）式にしたがってＷw が求められる。
Ｗw ＝ＡA ／Ｎ2                     …（２）
求められたＷw が、１、２、４、８、１６のいずれにも一致しない場合、求められたＷw に最も近いものが選択される。求められた、又は選択されたＷw から、データ数Ｎ´（＝Ｗw ・Ｎ）が決定される。
【００５１】
なお、原理的には、交互段数がｎ倍になると、周波数解析のデータ数は、１／ｎ1/2 倍に変更される。高速フーリエ変換処理回路５０ではデータ数調整用コントローラ５１から、決定されたデータ数Ｗw ・Ｎの情報を受け取り、ドプラ信号をデータ数Ｗw ・Ｎで周波数解析に供する。データ数調整用コントローラ５１で決定されたデータ数Ｗw ・Ｎの情報は、ドプラ感度補正用コントローラ５２にも供給される。ドプラ感度補正用コントローラ５２には、システムコントローラ４６から交互段数Ｋの情報も供給される。
【００５２】
ドプラ感度補正用コントローラ５２は、周波数解析のデータ数の変動に伴うＳ／Ｎの変動を抑制するためのものである。例えば、低速検出能を半値まで下げるように交互段数を変更すると、周波数解析のデータ数は、１／２1/2 倍に減少する。このデータ数の減少は、Ｓ／Ｎを０．８ｄＢ程度低下させる。このＳ／Ｎの低下を抑制するために、ドプラゲイン補正用の乗算器４９のドプラゲインを、変更前の交互段数のときより、０．８ｄＢ程度下げて、ノイズレベルを０．８ｄＢ程度下げる。これにより、データ数によらずノイズレベルを一定に保つことができるという効果が得られる。ただし、当然のこととして、ノイズレベルの低下と共に、信号レベルも０．８ｄＢ程度低下してしまう。この信号レベルの低下を補償するために、送信電圧を交互段数の変更前より１．０９倍に増大する。このようなドプラゲイン調整及び送信電圧調整により、交互段数が増加して、データ数が減少しても、Ｓ／Ｎが劣化しないという効果が得られる。
本発明は上述した実施の形態に限定されることなく種々変形して実施可能である。
【００５３】
【発明の効果】
本発明による超音波診断装置は、超音波パルスを被検体内に送信し、反射波を受信することにより受信信号を得る送受信手段と、レンジゲートマーカの深さ及び長さに応じた特定の期間に前記送受信手段から出力される前記受信信号を積分することにより積分信号を得る手段と、前記レンジゲートマーカ内のドプラスペクトラムを得るために前記積分信号を周波数解析する手段と、前記送受信手段と前記積分手段との間に設けられ、前記積分手段による積分期間をΔＴRGとして１／ΔＴRG１／２に応じた重み係数を前記受信信号に乗算する手段とを具備する。
本発明によれば、積分期間が変化すると、それに伴って実質的にゲインが変化してノイズレベルも変化してしまうが、積分期間に応じた重み係数を受信信号に乗算するので、積分期間の変化に伴う実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を連続波として被検体内に送信し、反射波を受信することにより受信信号を得る手段と、前記受信信号の特定の周波数帯域を減衰するフィルタ手段と、前記フィルタ手段の出力帯域に対する前記フィルタ手段の入力帯域の比に応じた重み係数を前記フィルタ手段の出力信号に乗算する手段と、ドプラスペクトラムを得るために前記フィルタ手段の出力信号を周波数解析する手段とを具備する。
本発明によれば、例えばレート周波数が変化すると、それに伴ってフィルタ手段の出力帯域に対する入力帯域の比が変化して実質的にゲインが変化してノイズレベルも変化してしまうが、この比に応じた重み係数をフィルタ手段の出力信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を被検体内に送信し、反射波を受信することにより受信信号を得る手段と、ドプラスペクトラムを得るために前記受信信号をデータ数で周波数解析する手段と、前記データ数をＮとして重み係数１／Ｎ１／２を前記受信信号に乗算する手段とを具備する。
本発明によれば、周波数解析のデータ数が変化すると実質的にゲインが変化してノイズレベルも変化してしまうが、このデータ数に応じた重み係数を受信信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
【図面の簡単な説明】
【図１】第１の実施の形態に係る超音波診断装置のブロック図。
【図２】第１の実施の形態による積分期間を示すタイムチャート。
【図３】第２の実施の形態に係る超音波診断装置のブロック図。
【図４】図３のＲＧ積分回路、ＣＷ帯域フィルタ、ＦＦＴのブロック図。
【図５】第３の実施の形態に係る超音波診断装置のブロック図。
【図６】図５のドプラ用コントローラのブロック図。
【図７】第１の実施の形態に対応する従来の問題点の説明図。
【符号の説明】
１…超音波プローブ、
２…送受信系、
３…スペクトラムドプラユニット、
４…ディジタルスキャンコンバータ、
５…モニタ、
６…コンソール、
７…コントローラ、
８…レンジゲートユニット、
９…高域通過型フィルタ、
１０…高速フーリエ変換処理回路、
１１…乗算器、
１２…積分器、
１３…関数発生ＲＯＭ、
１４…ＲＧ幅カウンタ、
１５…ＲＧスタート位置カウンタ。
【図１】
0003691874.tif 000002
【図２】
0003691874.tif 000003
【図３】
0003691874.tif 000004
【図４】
0003691874.tif 000005
【図５】
0003691874.tif 000006
【図６】
0003691874.tif 000007
【図７】
0003691874.tif 000008