(57)【要約】超音波ドプラ血流計は、ウォールフィルタの過渡現象を抑えるために、位相検波部１０４の検波出力からクラッタ成分の振幅、初期位相、位相速度をクラッタ成分推定部１１０で推定する。推定されたデータをもとにウォールフィルタ１０５の初期設定値をフィルタ初期値生成部１１で生成し、ウォールフィルタ１０５に与え、フィルタ処理を行なう。【特許請求の範囲】
【請求項１】
電気パルス信号を発生する送信部と、前記送信部で発生した電気パルスを超音波信号に変換し体内への送受信を行なう探触子と、前記探触子で得た体内からのエコー信号を位相検波する位相検波部と、前記位相検波部の出力から周波数の低いクラッタ成分を除去するウォールフィルタと、前記ウォールフィルタを通過した信号より血流情報を演算する血流情報演算部と、前記血流情報演算部で得た血流情報を走査変換するディジタルスキャンコンバータと、前記ディジタルスキャンコンバータで走査変換された信号を表示するモニタを備えた超音波ドプラ血流計において、前記位相検波部から出力された少なくとも２つの信号をもとに前記クラッタ成分を正弦波で近似するクラッタ成分推定部と、前記クラッタ成分推定部で近似された正弦波をもとに前記ウォールフィルタの初期値データを生成するフィルタ初期値生成部を備えていることを特徴とする超音波ドプラ血流計。
【請求項２】
推定されたクラッタ成分からフィルタのカットオフ周波数を設定するカットオフ周波数設定部を備えたことを特徴とする請求項１に記載の超音波ドプラ血流計。
【請求項３】
推定されたクラッタ成分の振幅からその振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする請求項２に記載の超音波ドプラ血流計。
【請求項４】
走査する部位の近傍の領域のクラッタ成分の振幅情報、または同じ部位の以前に走査した際のクラッタ成分の振幅情報をもとに、推定されたクラッタ成分の振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする請求項３に記載の超音波診断装置。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、医用分野において、超音波のドプラ現象を利用して体内の血流を測定し、測定された血流情報を画像表示する超音波ドプラ血流計に関する。
【背景技術】
【０００２】  超音波のドプラ現象を利用し、生体中の血流分布を色に対応させた状態で白黒の２次元断層像に重ね合わせて表示を行う、超音波ドプラ血流計（カラーフロー装置）が知られている（例えば、特許文献１）。
【特許文献１】特開平１１−２９９７８６号公報
【０００３】  図６に超音波ドプラ血流計１０Ａの概略構成を示す。図示する超音波血流計１０Ａでは、送信部１が、プローブ２を介して超音波パルスを被検者の体内に照射する。照射された超音波パルスの体内からのエコーは、プローブ２により電気信号に変換され、受信部３に入力される。受信部３は、微弱な信号を増幅し、図示しないＡ／Ｄ変換器により、アナログ信号からディジタル信号に変換する。その後、ディジタル信号は、位相検波部４に入力され、ベースバンドに変換されたエコー信号となる。ベースバンドに変換されたエコー信号は、ウォールフィルタ５に入力される。ウォールフィルタ５は、一般的にはＦＩＲ型フィルタもしくはＩＩＲ型フィルタによって構成されており、通常クラッタと呼ばれる低周波信号成分である不要な体内組織からの信号を除去する。ウォールフィルタ５によりクラッタを除去されたエコー信号は、血流情報演算部（速度演算部）６に送信され、そこで血流速度が算出される。この血流速度は、包絡線検波部７からのＢモード信号とともにディジタルスキャンコンバータ（以下、本明細書では「ＤＳＣ」と表す。）８に入力される。ＤＳＣ８は、Ｂモード信号と血流信号を混合し、その混合された信号を用いてモニタ９が二次元血流像を映し出す。【０００４】  図７は、２次のＩＩＲ型フィルタにより構成されたウォールフィルタの概略構成を示す。図中、ｘ（ｍ）はウォールフィルタ５に入力される信号（入力）、ｙ（ｍ）はウォールフィルタ５から出力される信号（出力）、２０１〜２０３は加算器、２１ａ２〜２１ａ３、２１ｂ１〜２１ｂ３は掛算器、２２１〜２２２は記憶素子（レジスタ）、ａ（２）〜ａ（３）、ｂ（１）〜ｂ（３）はフィルタ係数である。【０００５】  このような超音波ドプラ血流計においては、より正確な血流情報を得るためにはウォールフィルタ５に急峻な特性を持たせることが望ましいが、そのような急峻な特性は大きな過渡応答の発生を招き、このような過渡応答を含んだデータからは正確な血流情報を得ることができない。その結果、過渡応答の影響を除去して正確に血流速度を演算するために、多くのデータを破棄しなければならず、そのため速度演算部に送られるデータ点数が減り、また平均化の処理回数が減るので、信号／ノイズ比の低下を招くという問題があった。【０００６】  また、破棄データを補うためにデータサンプル数を増やすことは、同一方向への音波の繰り返し送受信回数を増やすこととなり、フレームレートが低下するという問題があった。【０００７】  さらに、従来の超音波ドプラ血流計においては、動きの速度変化の著しい部位、例えば心臓などにおいて、動きのゆっくりしている時相においてはクラッタの周波数がフィルタのカットオフ周波数以下であり、クラッタ信号は遮断されるが、動きの速い時相においてはクラッタの周波数がフィルタのカットオフ周波数以上になり、クラッタ信号が遮断されず血流速度の計算に誤差を生じることが発生する。この事象においては、クラッタ周波数の上限以上にフィルタのカットオフ周波数を上げると、動きのゆっくりしている時相において血流速度がカットオフ周波数より低くなってしまうため、血流が良好に表示されなくなるという問題も発生する。
【発明の開示】
【発明が解決しようとする課題】
【０００８】  本発明は、これらの問題を解決するためになされたもので、急峻な特性のウォールフィルタを用いつつ、過渡応答の少ないフィルタ出力を得、良好な血流表示を行なう超音波ドプラ血流計を提供することを目的とする。【０００９】  また、本発明は、従来の問題を解決するためになされたもので、動きの速度変化の激しい部位において、クラッタの周波数が大きく変動しても良好にクラッタを除去し、血流信号を抽出できるフィルタ処理を行ない、良好な血流表示を行なう超音波ドプラ血流計を提供することを目的とするものである。
【課題を解決するための手段】
【００１０】  このような課題を解決するために、本発明に係る超音波ドプラ血流計は、被検体中に送信した複数の超音波パルスに対する複数のエコー信号から、クラッタ成分の振幅、速度、位相を単一周波数の正弦波に近似して求め、この近似クラッタ成分がウォールフィルタに定常的に入力して、前記複数のエコー信号取得開始タイミングに至ったとき、ウォールフィルタ内のレジスタに蓄積されているべきデータを求め、ウォールフィルタに設定することで、信号の不連続性を抑え、過渡応答による不要な成分を低減し、良好な血流情報を得るものである。そのため、本発明の超音波ドプラ血流計は、電気パルス信号を発生する送信部と、前記送信部で発生した電気パルスを超音波信号に変換し体内への送受信を行なう探触子と、前記探触子で得た体内からのエコー信号を位相検波する位相検波部と、前記位相検波部の出力から周波数の低いクラッタ成分を除去するウォールフィルタと、前記ウォールフィルタを通過した信号より血流情報を演算する血流情報演算部と、前記血流情報演算部で得た血流情報を走査変換するディジタルスキャンコンバータと、前記ディジタルスキャンコンバータで走査変換された信号を表示するモニタを備えた超音波ドプラ血流計において、前記位相検波部から出力された少なくとも２つの信号をもとに前記クラッタ成分を正弦波で近似するクラッタ成分推定部と、前記クラッタ成分推定部で近似された正弦波をもとに前記ウォールフィルタの初期値データを生成するフィルタ初期値生成部を備えていることを特徴とする。【００１１】  本発明の他の形態に係る超音波ドプラ血流計は、単一周波数の正弦波に近似して求めたクラッタ信号の速度から、ウォールフィルタのカットオフ周波数を、このクラッタ成分を除去できるように設定する手段を設けることで、クラッタ成分を常に除去し、血流成分を常に抽出できるようにしたもので、推定されたクラッタ成分からフィルタのカットオフ周波数を設定するカットオフ周波数設定部を備えたことを特徴とする。【００１２】  また、本発明の他の形態に係る超音波ドプラ血流計は、単一周波数の正弦波に近似して求めたクラッタ信号の振幅から、この近似したクラッタ信号がクラッタであるか、それとも血流であるかを判断し、ウォールフィルタのカットオフ周波数を設定することで、クラッタ成分を常に除去し、血流成分を常に抽出できる精度を高めたもので、推定されたクラッタ成分の振幅からその振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする。【００１３】  さらに、本発明の他の形態に係る超音波ドプラ血流計は、すでに抽出を終えている、当該部位に近接した部位、もしくは過去に得た当該部位のクラッタ信号の情報から、クラッタ信号の見極めを行ない、ウォールフィルタの設定を行なうことで、クラッタ成分を常に除去し、血流成分を常に抽出できる精度を高めたもので、走査する部位の近傍の領域のクラッタ成分の振幅情報、または同じ部位の以前に走査した際のクラッタ成分の振幅情報をもとに、推定されたクラッタ成分の振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする。
【発明の効果】
【００１４】  このような構成を採用した本発明の超音波ドプラ血流計によれば、大振幅のクラッタによる過渡応答を抑えたフィルタを実現することができ、精度の良い血流情報を得ることができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施の形態１に係る超音波ドプラ血流計の構成を示すブロック図。
【図２】図１に示す超音波ドプラ血流計のウォールフィルタの過渡応答を防ぐ構成部分を示す図。
【図３】本発明の実施の形態２に係る超音波ドプラ血流計の構成を示すブロック図。
【図４】本発明の実施の形態３に係る超音波ドプラ血流計の構成を示すブロック図。
【図５】本発明の実施の形態４に係る超音波ドプラ血流計の構成を示すブロック図。
【図６】従来の超音波ドプラ血流計の構成を示すブロック図。
【図７】図６の超音波ドプラ血流計に含まれるウォールフィルタの構成を示す図。
【符号の説明】
【００１６】１，１０１…送信部、２，１０２…探触子、３，１０３…受信部、４，１０４…位相検波部、５，１０５，１０５Ａ，１０５Ｂ…ウォールフィルタ、６、１０６…血流情報演算部、７，１０７…包絡線検波部、８，１０８…ＤＳＣ、９，１０９…モニタ、１１０…クラッタ成分推定部、１１１…フィルタ初期値生成部、１１２…カットオフ周波数設定部、１１３…フィルタ係数生成部、１１５…クラッタ判定部、１１６…クラッタ情報蓄積部。
【発明を実施するための最良の形態】
【００１７】  以下、図面を参照して本発明の実施の形態を説明する。【００１８】実施の形態１【００１９】  図１は、本発明の実施の形態１に係る超音波ドプラ血流計の構成を示すブロック図である。図示する超音波ドプラ血流計１００Ａは、概略、図６を参照して説明した従来の血流計Ａにクラッタ成分推定部１１０とフィルタ初期値生成部１１１が付加されている。【００２０】  このように構成された血流計１００Ａの動作時、送信部１０１は、プローブ１０２を介して超音波パルスを被検者の体内に照射する。照射された超音波パルスの体内からのエコーは、プローブ１０２により電気信号に変換され、受信部１０３に入力される。受信部１０３は、微弱な信号を増幅し、図示しないＡ／Ｄ変換器により、アナログ信号からディジタル信号に変換する。その後、ディジタル信号は、位相検波部１０４に入力され、ベースバンドに変換されたエコー信号となる。ベースバンドに変換されたエコー信号は、ウォールフィルタ１０５とクラッタ成分推定部１１０に入力される。クラッタ成分推定部１１０は、位相検波部１０４の出力を用いてクラッタ成分を推定する。フィルタ初期値生成部１１１は、推定されたクラッタ成分を用いてフィルタ初期値を生成し、この生成されたフィルタ初期値をウォールフィルタ１０５に提供する。ウォールフィルタ１０５は、一般的にはＦＩＲ型フィルタもしくはＩＩＲ型フィルタによって構成されており、フィルタ初期値生成部１１１から提供されたフィルタ初期値を用いて、クラッタと呼ばれる低周波信号成分である不要な体内組織からの信号を除去する。ウォールフィルタ１０５によりクラッタを除去されたエコー信号は、血流情報演算部（速度演算部）１０６に送信され、そこで血流速度が算出される。この血流速度は、包絡線検波部１０７からのＢモード信号とともにディジタルスキャンコンバータ（以下、本明細書では「ＤＳＣ」と表す。）１０８に入力される。ＤＳＣ１０８は、Ｂモード信号と血流信号を混合し、その混合された信号を用いてモニタ１０９が二次元血流像を映し出す。【００２１】  図２を参照して、血流計１００Ａにおける特徴部分の構成と動作（処理）を具体的に説明する。図示するように、血流計１００Ａは、記憶部１１２を備えており、位相検波部１０４から出力された位相検波出力ｘｉ（１）〜ｘｉ（ｎ）、ｘｑ（１）〜ｘｑ（ｎ）を記憶部１１２に記憶する。記憶部１１２は、位相検波出力ｘｉ（１）〜ｘｉ（ｎ）、ｘｑ（１）〜ｘｑ（ｎ）を対応するウォールフィルタ１０５Ａ，１０５Ｂにそれぞれ出力するとともに、第１と第２の位相検波出力ｘｉ（１），ｘｉ（２）、ｘｑ（１），ｘｑ（２）をクラッタ成分推定部１１０に出力する。【００２２】  クラッタ成分推定部１１０は、第１と第２の位相検波出力ｘｉ（１），ｘｉ（２）、ｘｑ（１），ｘｑ（２）を用い、位相検波出力に含まれるクラッタ成分を正弦波で近似して推定する。この推定プロセスにおいて、クラッタ成分推定部１１０は、数式１，２に基づいて複素数Ｅ１，Ｅ２を計算する。
【数１】
2006088094.tif 000002
【数２】
2006088094.tif 000003
【００２３】  推定する正弦波の振幅をａ、位相をθｃ（ｔ）とすると、位相θｃ（ｔ）は数式３で表される。
【数３】
2006088094.tif 000004
【００２４】  ここで、初期位相ｍと位相速度ｎは、数式４，５で表される。
【数４】
2006088094.tif 000005
【数５】
2006088094.tif 000006
【００２５】  ａｎｇｌｅ（Ｅ）は複素数Ｅの位相角、Ｅ＊は共役複素数、Ｔ１，Ｔ２はそれぞれ第１と第２の検波出力のタイミング（受信時刻に対応する）である。【００２６】  一方、振幅ａは数式６で表される。
【数６】
2006088094.tif 000007
【００２７】  クラッタ成分推定部１１０は、以上の数式に基づいて計算された振幅ａと位相θｃ（ｔ）で定義される正弦波を推定した後、この正弦波に対応する信号をフィルタ初期値生成部１１１に出力する。フィルタ初期値生成部１１１は、クラッタ成分推定部１１０から出力された信号に基づいてフィルタ初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０を作成し、これらの初期値データをウォールフィルタ６Ａ，６Ｂに出力する。【００２８】  ウォールフィルタ１０５Ａ，１０５Ｂは図２に示す構成を有する。図中、１２０１〜１２０３は加算器、１２１ａ２〜１２１ａ３、１２１ｂ１〜１２１ｂ３は掛算器、１２２１，１２２２は記憶素子（レジスタ）、１２３１，１２３２はスイッチ、ａｉ（２）〜ａｉ（３）、ｂｉ（１）〜ｂｉ（３）、ａｑ（２）〜ａｑ（３）、ｂｑ（１）〜ｂｑ（３）はフィルタ係数である。【００２９】  このように構成されたウォールフィルタ１０５Ａ，１０５Ｂに対し、フィルタ初期値生成部１１１で生成された初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０に対応する信号がフィルタ１０５Ａのスイッチ１２３１，１２３２とフィルタ１０５Ｂのスイッチ１２３１，１２３２を通じて入力される。初期値データが入力されたウォールフィルタ１０５Ａ，１０５Ｂは、それら初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０をもとに、記憶部１１２からウォールフィルタ１０５Ａ，１０５Ｂへの位相検波出力ｘｉ（１）〜ｘｉ（ｎ）、ｘｑ（１）〜ｘｑ（ｎ）を処理してクラッタ成分を除去する。【００３０】  図２に示す構成において、ウォールフィルタ１０５Ａ，１０５Ｂの出力ｙｉ（ｍ）、ｙｑ（ｍ）、また各ウォールフィルタ１０５Ａ，１０５Ｂにおける記憶素子１２２１，１２２２の入出力ｚｉ１（ｍ）、ｚｉ２（ｍ）、ｚｑ１（ｍ）、ｚｑ２（ｍ）は以下の数式７〜１２で定義される。
【数７】
2006088094.tif 000008
【数８】
2006088094.tif 000009
【数９】
2006088094.tif 000010
【数１０】
2006088094.tif 000011
【数１１】
2006088094.tif 000012
【数１２】
2006088094.tif 000013
【００３１】  そして、ウォールフィルタ１０５Ａ，１０５Ｂは、数式７〜数１２のｚ１ｉ（ｍ−１）、ｚｉ２（ｍ−１）、ｚ１ｑ（ｍ−１）、ｚｑ２（ｍ−１）に上述の初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０を入力し、これらの初期値データと記憶部１１２からウォールフィルタ１０５Ａ，１０５Ｂに入力されたデータｘｉ（１）、ｘｑ（１）をもとに、ｙｉ（１）、ｙｑ（１）を計算する。以後同様にして、（ｍ−１）番目の位相検波出力に対する計算時に計算された出力ｚ１ｉ（ｍ−１）、ｚｉ２（ｍ−１）、ｚ１ｑ（ｍ−１）、ｚｑ２（ｍ−１）を用いて、ｍ番目の位相検波出力を処理する。【００３２】  このように、ウォールフィルタ１０５Ａ，１０５Ｂでは、最初の位相検波出力ｘｉ（１）、ｘｑ（１）について、クラッタ成分を正弦波で近似して初期値データを利用して出力ｙｉ（１）、ｙｑ（１）を得る。そのため、ウォールフィルタ１０５Ａ，１０５Ｂから発生する過渡応答は低減し、血流信号を精度良く抽出できる。【００３３】  以上の説明では、第１と第２の位相検波出力ｘｉ（１），ｘｉ（２）、ｘｑ（１），ｘｑ（２）を用いて、位相検波出力に含まれるクラッタ成分を正弦波で近似したが、例えばクラッタの速度変化が小さい場合、ｎ≧ｕ＞２の関係を有する任意のｕ（整数）を選択し、第１番目と第ｕ番目の位相検波出力ｘｉ（１），ｘｉ（ｕ）、ｘｑ（１），ｘｑ（ｕ）を用い、上述と同様の方法により、近似する正弦波を推定してもよい。この場合、推定される正弦波の位相速度ｎは数式１３で定義される。
【数１３】
2006088094.tif 000014
【００３４】実施の形態２【００３５】  図３は本発明の実施の形態２に係る超音波ドプラ血流計１００Ｂのブロック図である。図示する超音波ドプラ血流計１００Ｂは、実施の形態１の超音波ドプラ血流計１００Ａ（図２参照）と比べてみると明らかなように、カットオフ周波数設定部１１３とフィルタ係数生成部１１４が付加されている点に特徴を有する。【００３６】  この超音波ドプラ血流計１００Ｂによれば、位相検波部１０４の検波出力信号をもとにクラッタ成分推定部１１０がクラッタを正弦波で近似するとともに、この正弦波の振幅と位相速度を推定する。推定された情報は、実施の形態１で説明したように、フィルタ初期値生成部１１１に入力されるとともに、カットオフ周波数設定部１１３に入力される。【００３７】  とろこで、超音波ドプラ血流計にあっては、多くの場合、心臓のようなクラッタおよび血流の位相速度の変化が激しい部位において、つねにクラッタのみを除去し、血流のみを抽出することが望まれる。したがって、カットオフ周波数設定部１１３では、クラッタの位相速度と振幅をもとに、クラッタが除去される最低のカットオフ周波数を設定し、フィルタ係数生成部１１４に対応するデータを送信する。フィルタ係数生成部１１４は、選定されたカットオフ周波数のフィルタを実現するフィルタ係数〔図２に示すフィルタ係数ａｉ（２）〜ａｉ（３）、ｂｉ（１）〜ｂｉ（３）、ａｑ（２）〜ａｑ（３）、ｂｑ（１）〜ｂｑ（３）〕を生成し、ウォールフィルタ１０５に与える。ウォールフィルタ１０５は、フィルタ係数生成部１１４から与えられたフィルタ係数を用いて出力ｙｉ（ｍ）、ｙｑ（ｍ）を求める。【００３８】  実施の形態１ではフィルタ係数は不変としているが、実施の形態２ではフィルタ係数が新たな血流情報を得る度に変化する。したがって、実施の形態２では、フィルタ係数が変化に応じて、フィルタ初期値を変化させることが好ましい。そこで、実施の形態２では、フィルタ係数生成部１１４で生成されたフィルタ係数をフィルタ初期値生成部１１１に出力する。フィルタ初期値生成部１１１は、フィルタ係数生成部１１４から受け取ったフィルタ係数と、クラッタ成分推定値１１０で推定されたクラッタ近似値の振幅、初期位相、位相速度をもとに、フィルタ初期値を生成してウォールフィルタ１０５を設定する。【００３９】  フィルタ係数生成部１１４は、予め決められた方法に基づいてカットオフ周波数からクラッタ係数を計算してもよいが、カットオフ周波数とフィルタ係数を相互に関連付けたテーブルを用意（記憶）しておき、このテーブルに基づいてフィルタ係数を選択してもよい。【００４０】  このように、実施の形態２に係る超音波ドプラ血流計１００Ｂによれば、心臓のようなクラッタおよび血流の位相速度の変化が激しい部位において、つねにクラッタのみを除去し、血流のみを抽出することができ、精度の高いデータを提供することができる。【００４１】実施の形態３【００４２】  図４は本発明の実施の形態３に係る超音波ドプラ血流計１００Ｃのブロック図である。図示する超音波ドプラ血流計１００Ｃは、実施の形態２の超音波ドプラ血流計１００Ｂ（図３参照）と比べてみると明らかなように、クラッタ判定部１１４が付加さ���ている点に特徴を有する。【００４３】  超音波診断法にあっては、体内から返ってくるエコー信号において、多くの場合はクラッタの振幅が血流の振幅の数十倍もしくは百数十倍と大きい。しかし、例えば心臓の腔内のような血液のみが存在している部位からのエコー信号では、血流の振幅がクラッタの振幅より大きくなることがある。このような場合において、実施の形態２のような方法でクラッタ信号を同定し、除去を行なうと、振幅の大きな血流信号がクラッタとみなされてしまい、除去されてしまうおそれがある。【００４４】  したがって、実施の形態３では、クラッタ判定部１１５においてクラッタと血流の判別を行なう。具体的には、クラッタ判定部１１５は、振幅の閾値を設定し、クラッタ成分推定部１１０で近似された正弦波の振幅がこの閾値を超える場合にある基準以上の振幅の場合は、推定された正弦波をクラッタとみなして以後の信号処理を行なう。しかし、推定された正弦波の振幅が閾値以下の振幅の場合は血流とみなし、上述したフィルタ初期値生成部１１１、カットオフ周波数設定部１１３、フィルタ係数生成部１１４における処理の実行を禁止する。【００４５】  このように構成された実施の形態３に係る超音波ドプラ血流計１００Ｃによれば、心腔内のような血流信号がクラッタ信号より振幅が大きい場合にも血流信号を抽出でき、精度の高いデータを提供することができる。【００４６】実施の形態４【００４７】  図５は本発明の実施の形態４に係る超音波ドプラ血流計１００Ｄのブロック図である。図示する超音波ドプラ血流計１００Ｄは、実施の形態３の超音波ドプラ血流計１００Ｃ（図４参照）と比べてみると明らかなように、クラッタ情報蓄積部１１５が付加されている点に特徴を有する。【００４８】  具体的に説明すると、上述した実施の形態３では、振幅の閾値を設けることで、エコー信号中に含まれる最も大きな振幅の信号がクラッタであるか、血流であるかを判断したが、実施の形態４では、閾値を設定するための手段としてクラッタ情報蓄積部１５が設けられている。【００４９】  超音波ドプラ血流計においては、体内を順次、走査していく過程において、近傍の部位の走査がすでに終わっており、クラッタや血流の情報が既に得られ、また、同じ部位を既に走査してその際にデータが取得されていることが通常である。そこで、実施の形態４の超音波ドプラ血流計１００Ｄでは、既に取得された情報をクラッタ情報蓄積部１１５に蓄積しておき、この蓄積された情報を用いてクラッタ判定部１１４又はクラッタ情報蓄積部１１５が閾値を設定する。【００５０】  ところで、超音波ドプラ血流計の全走査範囲において、すべての範囲でクラッタが血流より小さいということは考えにくい。クラッタが血流より大きく、かつ、そのレベルは血流信号ではありえない振幅になるのが常である。そこで、現在走査している部位の近傍においてクラッタが血流より大きい部位を探すことで、クラッタと血流の振幅の敷居値を設定することができる。また、クラッタの振幅について空間的な分布を求めることで、現在走査している場所がどの程度のクラッタがあるかを予測することもできる。同じ部位について言えば、以前その部位において、例えば心臓の弁が横切るなどし��、大振幅の信号が入力した経緯があるならば、それを基準にクラッタと血流の敷居値を設けることができる。したがって、実施の形態４によれば、蓄積された情報をもとに振幅の閾値が設定されるので、精度の良い血流情報を得ることができる。
【図１】
2006088094.tif 000015
【図２】
2006088094.tif 000016
【図３】
2006088094.tif 000017
【図４】
2006088094.tif 000018
【図５】
2006088094.tif 000019
【図６】
2006088094.tif 000020
【図７】
2006088094.tif 000021
【手続補正書】
【提出日】平成18年7月27日(2006.7.27)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】請求の範囲
【補正方法】変更
【補正の内容】
【書類名】請求の範囲
【請求項１】（補正後）  電気パルス信号を発生する送信部と、前記送信部で発生した電気パルスを超音波信号に変換し体内への送受信を行なう探触子と、前記探触子で得た体内からのエコー信号を位相検波する位相検波部と、前記位相検波部の出力から周波数の低いクラッタ成分を除去するウォールフィルタと、前記ウォールフィルタを通過した信号より血流情報を演算する血流情報演算部と、前記血流情報演算部で得た血流情報を走査変換するディジタルスキャンコンバータと、前記ディジタルスキャンコンバータで走査変換された信号を表示するモニタを備えた超音波ドプラ血流計において、
前記位相検波部から出力された少なくとも２つの信号をもとに前記クラッタ成分を正弦波で推定するクラッタ成分推定部と、
前記クラッタ成分推定部で近似された正弦波をもとに前記ウォールフィルタの初期値データを生成するフィルタ初期値生成部と、
前記推定されたクラッタ成分からフィルタのカットオフ周波数を設定するカットオフ周波数設定部と、
前記推定されたクラッタ成分の振幅情報を蓄積するクラッタ情報蓄積部と、
前記蓄積部に蓄積されている、走査する部位の近傍の領域のクラッタ成分の振幅情報、または同じ部位の以前に走査した際のクラッタ成分の振幅情報をもとに、推定されたクラッタ成分の振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする超音波診断装置。
【請求項２】（削除）
【請求項３】（削除）
【請求項４】（削除）
１３
【手続補正書】
【提出日】平成19年9月13日(2007.9.13)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
電気パルス信号を発生する送信部と、前記送信部で発生した電気パルスを超音波信号に変換し体内への送受信を行なう探触子と、前記探触子で得た体内からのエコー信号を位相検波する位相検波部と、前記位相検波部の出力から周波数の低いクラッタ成分を除去するウォールフィルタと、前記ウォールフィルタを通過した信号より血流情報を演算する血流情報演算部と、前記血流情報演算部で得た血流情報を走査変換するディジタルスキャンコンバータと、前記ディジタルスキャンコンバータで走査変換された信号を表示するモニタを備えた超音波ドプラ血流計において、
前記位相検波部から出力された少なくとも２つの信号をもとに前記クラッタ成分を正弦波で推定するクラッタ成分推定部と、
前記クラッタ成分推定部で近似された正弦波をもとに前記ウォールフィルタの初期値データを生成するフィルタ初期値生成部と、
前記推定されたクラッタ成分からフィルタのカットオフ周波数を設定するカットオフ周波数設定部と、
前記推定されたクラッタ成分の振幅情報を蓄積するクラッタ情報蓄積部と、
前記蓄積部に蓄積されている、走査する部位の近傍の領域のクラッタ成分の振幅情報、または同じ部位の以前に走査した際のクラッタ成分の振幅情報をもとに、推定されたクラッタ成分の振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする超音波診断装置。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用分野において、超音波のドプラ現象を利用して体内の血流を測定し、測定された血流情報を画像表示する超音波ドプラ血流計に関する。
【背景技術】
【０００２】
超音波のドプラ現象を利用し、生体中の血流分布を色に対応させた状態で白黒の２次元断層像に重ね合わせて表示を行う、超音波ドプラ血流計（カラーフロー装置）が知られている（例えば、特許文献１）。
【特許文献１】特開平１１−２９９７８６号公報
【０００３】
図６に超音波ドプラ血流計１０Ａの概略構成を示す。図示する超音波血流計１０Ａでは、送信部１が、プローブ２を介して超音波パルスを被検者の体内に照射する。照射された超音波パルスの体内からのエコーは、プローブ２により電気信号に変換され、受信部３に入力される。受信部３は、微弱な信号を増幅し、図示しないＡ／Ｄ変換器により、アナログ信号からディジタル信号に変換する。その後、ディジタル信号は、位相検波部４に入力され、ベースバンドに変換されたエコー信号となる。ベースバンドに変換されたエコー信号は、ウォールフィルタ５に入力される。ウォールフィルタ５は、一般的にはＦＩＲ型フィルタもしくはＩＩＲ型フィルタによって構成されており、通常クラッタと呼ばれる低周波信号成分である不要な体内組織からの信号を除去する。ウォールフィルタ５によりクラッタを除去されたエコー信号は、血流情報演算部（速度演算部）６に送信され、そこで血流速度が算出される。この血流速度は、包絡線検波部7からのＢモード信号とともにディジタルスキャンコンバータ(以下、本明細書では「ＤＳＣ」と表す。)８に入力される。ＤＳＣ８は、Ｂモード信号と血流信号を混合し、その混合された信号を用いてモニタ９が二次元血流像を映し出す。
【０００４】
図７は、２次のＩＩＲ型フィルタにより構成されたウォールフィルタの概略構成を示す。図中、ｘ（ｍ）はウォールフィルタ５に入力される信号（入力）、ｙ（ｍ）はウォールフィルタ５から出力される信号（出力）、２０１〜２０３は加算器、２１ａ２〜２１ａ３、２１ｂ１〜２１ｂ３は掛算器、２２１〜２２２は記憶素子（レジスタ）、ａ（２）〜ａ（３）、ｂ（１）〜ｂ（３）はフィルタ係数である。
【０００５】
このような超音波ドプラ血流計においては、より正確な血流情報を得るためにはウォールフィルタ５に急峻な特性を持たせることが望ましいが、そのような急峻な特性は大きな過渡応答の発生を招き、このような過渡応答を含んだデータからは正確な血流情報を得ることができない。その結果、過渡応答の影響を除去して正確に血流速度を演算するために、多くのデータを破棄しなければならず、そのため速度演算部に送られるデータ点数が減り、また平均化の処理回数が減るので、信号/ノイズ比の低下を招くという問題があった。
【０００６】
また、破棄データを補うためにデータサンプル数を増やすことは、同一方向への音波の繰り返し送受信回数を増やすこととなり、フレームレートが低下するという問題があった。
【０００７】
さらに、従来の超音波ドプラ血流計においては、動きの速度変化の著しい部位、例えば心臓などにおいて、動きのゆっくりしている時相においてはクラッタの周波数がフィルタのカットオフ周波数以下であり、クラッタ信号は遮断されるが、動きの速い時相においてはクラッタの周波数がフィルタのカットオフ周波数以上になり、クラッタ信号が遮断されず血流速度の計算に誤差を生じることが発生する。この事象においては、クラッタ周波数の上限以上にフィルタのカットオフ周波数を上げると、動きのゆっくりしている時相において血流速度がカットオフ周波数より低くなってしまうため、血流が良好に表示されなくなるという問題も発生する。
【発明の開示】
【発明が解決しようとする課題】
【０００８】
本発明は、これらの問題を解決するためになされたもので、急峻な特性のウォールフィルタを用いつつ、過渡応答の少ないフィルタ出力を得、良好な血流表示を行なう超音波ドプラ血流計を提供することを目的とする。
【０００９】
また、本発明は、従来の問題を解決するためになされたもので、動きの速度変化の激しい部位において、クラッタの周波数が大きく変動しても良好にクラッタを除去し、血流信号を抽出できるフィルタ処理を行ない、良好な血流表示を行なう超音波ドプラ血流計を提供することを目的とするものである。
【課題を解決するための手段】
【００１０】
このような課題を解決するために、本発明に係る超音波ドプラ血流計は、被検体中に送信した複数の超音波パルスに対する複数のエコー信号から、クラッタ成分の振幅、速度、位相を単一周波数の正弦波に近似して求め、この近似クラッタ成分がウォールフィルタに定常的に入力して、前記複数のエコー信号取得開始タイミングに至ったとき、ウォールフィルタ内のレジスタに蓄積されているべきデータを求め、ウォールフィルタに設定することで、信号の不連続性を抑え、過渡応答による不要な成分を低減し、良好な血流情報を得るものである。そのため、本発明の超音波ドプラ血流計は、電気パルス信号を発生する送信部と、前記送信部で発生した電気パルスを超音波信号に変換し体内への送受信を行なう探触子と、前記探触子で得た体内からのエコー信号を位相検波する位相検波部と、前記位相検波部の出力から周波数の低いクラッタ成分を除去するウォールフィルタと、前記ウォールフィルタを通過した信号より血流情報を演算する血流情報演算部と、前記血流情報演算部で得た血流情報を走査変換するディジタルスキャンコンバータと、前記ディジタルスキャンコンバータで走査変換された信号を表示するモニタを備えた超音波ドプラ血流計において、前記位相検波部から出力された少なくとも２つの信号をもとに前記クラッタ成分を正弦波で近似するクラッタ成分推定部と、前記クラッタ成分推定部で近似された正弦波をもとに前記ウォールフィルタの初期値データを生成するフィルタ初期値生成部を備えていることを特徴とする。
【００１１】
本発明の他の形態に係る超音波ドプラ血流計は、単一周波数の正弦波に近似して求めたクラッタ信号の速度から、ウォールフィルタのカットオフ周波数を、このクラッタ成分を除去できるように設定する手段を設けることで、クラッタ成分を常に除去し、血流成分を常に抽出できるようにしたもので、推定されたクラッタ成分からフィルタのカットオフ周波数を設定するカットオフ周波数設定部を備えたことを特徴とする。
【００１２】
また、本発明の他の形態に係る超音波ドプラ血流計は、単一周波数の正弦波に近似して求めたクラッタ信号の振幅から、この近似したクラッタ信号がクラッタであるか、それとも血流であるかを判断し、ウォールフィルタのカットオフ周波数を設定することで、クラッタ成分を常に除去し、血流成分を常に抽出できる精度を高めたもので、推定されたクラッタ成分の振幅からその振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする。
【００１３】
さらに、本発明の他の形態に係る超音波ドプラ血流計は、すでに抽出を終えている、当該部位に近接した部位、もしくは過去に得た当該部位のクラッタ信号の情報から、クラッタ信号の見極めを行ない、ウォールフィルタの設定を行なうことで、クラッタ成分を常に除去し、血流成分を常に抽出できる精度を高めたもので、走査する部位の近傍の領域のクラッタ成分の振幅情報、または同じ部位の以前に走査した際のクラッタ成分の振幅情報をもとに、推定されたクラッタ成分の振幅がクラッタに起因するのか血流に起因するのかを判定するクラッタ判定部を備えたことを特徴とする。
【発明の効果】
【００１４】
このような構成を採用した本発明の超音波ドプラ血流計によれば、大振幅のクラッタによる過渡応答を抑えたフィルタを実現することができ、精度の良い血流情報を得ることができる。
【発明を実施するための最良の形態】
【００１５】
以下、図面を参照して本発明の実施の形態を説明する。
【００１６】
実施の形態１
【００１７】
図１は、本発明の実施の形態１に係る超音波ドプラ血流計の構成を示すブロック図である。図示する超音波ドプラ血流計１００Ａは、概略、図６を参照して説明した従来の血流計Ａにクラッタ成分推定部１１０とフィルタ初期値生成部１１１が付加されている。
【００１８】
このように構成された血流計１００Ａの動作時、送信部１０１は、プローブ１０２を介して超音波パルスを被検者の体内に照射する。照射された超音波パルスの体内からのエコーは、プローブ１０２により電気信号に変換され、受信部１０３に入力される。受信部１０３は、微弱な信号を増幅し、図示しないＡ／Ｄ変換器により、アナログ信号からディジタル信号に変換する。その後、ディジタル信号は、位相検波部１０４に入力され、ベースバンドに変換されたエコー信号となる。ベースバンドに変換されたエコー信号は、ウォールフィルタ１０５とクラッタ成分推定部１１０に入力される。クラッタ成分推定部１１０は、位相検波部１０４の出力を用いてクラッタ成分を推定する。フィルタ初期値生成部１１１は、推定されたクラッタ成分を用いてフィルタ初期値を生成し、この生成されたフィルタ初期値をウォールフィルタ１０５に提供する。ウォールフィルタ１０５は、一般的にはＦＩＲ型フィルタもしくはＩＩＲ型フィルタによって構成されており、フィルタ初期値生成部１１１から提供されたフィルタ初期値を用いて、クラッタと呼ばれる低周波信号成分である不要な体内組織からの信号を除去する。ウォールフィルタ１０５によりクラッタを除去されたエコー信号は、血流情報演算部（速度演算部）１０６に送信され、そこで血流速度が算出される。この血流速度は、包絡線検波部１０７からのＢモード信号とともにディジタルスキャンコンバータ(以下、本明細書では「ＤＳＣ」と表す。)１０８に入力される。ＤＳＣ１０８は、Ｂモード信号と血流信号を混合し、その混合された信号を用いてモニタ１０９が二次元血流像を映し出す。
【００１９】
図２を参照して、血流計１００Ａにおける特徴部分の構成と動作（処理）を具体的に説明する。図示するように、血流計１００Ａは、記憶部１１２を備えており、位相検波部１０４から出力された位相検波出力ｘｉ（１）〜ｘｉ（ｎ）、ｘｑ（１）〜ｘｑ（ｎ）を記憶部１１２に記憶する。記憶部１１２は、位相検波出力ｘｉ（１）〜ｘｉ（ｎ）、ｘｑ（１）〜ｘｑ（ｎ）を対応するウォールフィルタ１０５Ａ，１０５Ｂにそれぞれ出力するとともに、第１と第２の位相検波出力ｘｉ（１），ｘｉ（２）、ｘｑ（１），ｘｑ（２）をクラッタ成分推定部１１０に出力する。
【００２０】
クラッタ成分推定部１１０は、第１と第２の位相検波出力ｘｉ（１），ｘｉ（２）、ｘｑ（１），ｘｑ（２）を用い、位相検波出力に含まれるクラッタ成分を正弦波で近似して推定する。この推定プロセスにおいて、クラッタ成分推定部１１０は、数式１，２に基づいて複素数Ｅ１，Ｅ２を計算する。
【数１】
2006088094.tif 000022
【数２】
2006088094.tif 000023
【００２１】
推定する正弦波の振幅をａ、位相をθｃ（ｔ）とすると、位相θｃ（ｔ）は数式３で表される。
【数３】
2006088094.tif 000024
【００２２】
ここで、初期位相ｍと位相速度ｎは、数式４，５で表される。
【数４】
2006088094.tif 000025
【数５】
2006088094.tif 000026
【００２３】
ａｎｇｌｅ（Ｅ）は複素数Ｅの位相角、Ｅ＊は共役複素数、Ｔ１，Ｔ２はそれぞれ第１と第２の検波出力のタイミング（受信時刻に対応する）である。
【００２４】
一方、振幅ａは数式６で表される。
【数６】
2006088094.tif 000027
【００２５】
クラッタ成分推定部１１０は、以上の数式に基づいて計算された振幅ａと位相θｃ（ｔ）で定義される正弦波を推定した後、この正弦波に対応する信号をフィルタ初期値生成部１１１に出力する。フィルタ初期値生成部１１１は、クラッタ成分推定部１１０から出力された信号に基づいてフィルタ初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０を作成し、これらの初期値データをウォールフィルタ６Ａ，６Ｂに出力する。
【００２６】
ウォールフィルタ１０５Ａ，１０５Ｂは図２に示す構成を有する。図中、１２０１〜１２０３は加算器、１２１ａ２〜１２１ａ３、１２１ｂ１〜１２１ｂ３は掛算器、１２２１，１２２２は記憶素子（レジスタ）、１２３１，１２３２はスイッチ、ａｉ（２）〜ａｉ（３）、ｂｉ（１）〜ｂｉ（３）、ａｑ（２）〜ａｑ（３）、ｂｑ（１）〜ｂｑ（３）はフィルタ係数である。
【００２７】
このように構成されたウォールフィルタ１０５Ａ，１０５Ｂに対し、フィルタ初期値生成部１１１で生成された初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０に対応する信号がフィルタ１０５Ａのスイッチ１２３１，１２３２とフィルタ１０５Ｂのスイッチ１２３１，１２３２を通じて入力される。初期値データが入力されたウォールフィルタ１０５Ａ，１０５Ｂは、それら初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０をもとに、記憶部１１２からウォールフィルタ１０５Ａ，１０５Ｂへの位相検波出力ｘｉ（１）〜ｘｉ（ｎ）、ｘｑ（１）〜ｘｑ（ｎ）を処理してクラッタ成分を除去する。
【００２８】
図２に示す構成において、ウォールフィルタ１０５Ａ，１０５Ｂの出力ｙｉ（ｍ）、ｙｑ（ｍ）、また各ウォールフィルタ１０５Ａ，１０５Ｂにおける記憶素子１２２１，１２２２の入出力ｚｉ１（ｍ）、ｚｉ２（ｍ）、ｚｑ１（ｍ）、ｚｑ２（ｍ）は以下の数式７〜１２で定義される。
【数７】
2006088094.tif 000028
【数８】
2006088094.tif 000029
【数９】
2006088094.tif 000030
【数１０】
2006088094.tif 000031
【数１１】
2006088094.tif 000032
【数１２】
2006088094.tif 000033
【００２９】
そして、ウォールフィルタ１０５Ａ，１０５Ｂは、数式７〜数１２のｚ１ｉ（ｍ−１）、ｚｉ２（ｍ−１）、ｚ１ｑ（ｍ−１）、ｚｑ２（ｍ−１）に上述の初期値データｚｉ１０，ｚｉ２０、ｚｑ１０，ｚｑ２０を入力し、これらの初期値データと記憶部１１２からウォールフィルタ１０５Ａ，１０５Ｂに入力されたデータｘｉ（１）、ｘｑ（１）をもとに、ｙｉ（１）、ｙｑ（１）を計算する。以後同様にして、（ｍ−１）番目の位相検波出力に対する計算時に計算された出力ｚ１ｉ（ｍ−１）、ｚｉ２（ｍ−１）、ｚ１ｑ（ｍ−１）、ｚｑ２（ｍ−１）を用いて、ｍ番目の位相検波出力を処理する。
【００３０】
このように、ウォールフィルタ１０５Ａ，１０５Ｂでは、最初の位相検波出力ｘｉ（１）、ｘｑ（１）について、クラッタ成分を正弦波で近似して初期値データを利用して出力ｙｉ（１）、ｙｑ（１）を得る。そのため、ウォールフィルタ１０５Ａ，１０５Ｂから発生する過渡応答は低減し、血流信号を精度良く抽出できる。
【００３１】
以上の説明では、第１と第２の位相検波出力ｘｉ（１），ｘｉ（２）、ｘｑ（１），ｘｑ（２）を用いて、位相検波出力に含まれるクラッタ成分を正弦波で近似したが、例えばクラッタの速度変化が小さい場合、ｎ≧ｕ＞２の関係を有する任意のｕ（整数）を選択し、第１番目と第ｕ番目の位相検波出力ｘｉ（１），ｘｉ（ｕ）、ｘｑ（１），ｘｑ（ｕ）を用い、上述と同様の方法により、近似する正弦波を推定してもよい。この場合、推定される正弦波の位相速度ｎは数式１３で定義される。
【数１３】
2006088094.tif 000034
【００３２】
実施の形態２
【００３３】
図３は本発明の実施の形態２に係る超音波ドプラ血流計１００Ｂのブロック図である。図示する超音波ドプラ血流計１００Ｂは、実施の形態１の超音波ドプラ血流計１００Ａ（図２参照）と比べてみると明らかなように、カットオフ周波数設定部１１３とフィルタ係数生成部１１４が付加されている点に特徴を有する。
【００３４】
この超音波ドプラ血流計１００Ｂによれば、位相検波部１０４の検波出力信号をもとにクラッタ成分推定部１１０がクラッタを正弦波で近似するとともに、この正弦波の振幅と位相速度を推定する。推定された情報は、実施の形態１で説明したように、フィルタ初期値生成部１１１に入力されるとともに、カットオフ周波数設定部１１３に入力される。
【００３５】
とろこで、超音波ドプラ血流計にあっては、多くの場合、心臓のようなクラッタおよび血流の位相速度の変化が激しい部位において、つねにクラッタのみを除去し、血流のみを抽出することが望まれる。したがって、カットオフ周波数設定部１１３では、クラッタの位相速度と振幅をもとに、クラッタが除去される最低のカットオフ周波数を設定し、フィルタ係数生成部１１４に対応するデータを送信する。フィルタ係数生成部１１４は、選定されたカットオフ周波数のフィルタを実現するフィルタ係数〔図２に示すフィルタ係数ａｉ（２）〜ａｉ（３）、ｂｉ（１）〜ｂｉ（３）、ａｑ（２）〜ａｑ（３）、ｂｑ（１）〜ｂｑ（３）〕を生成し、ウォールフィルタ１０５に与える。ウォールフィルタ１０５は、フィルタ係数生成部１１４から与えられたフィルタ係数を用いて出力ｙｉ（ｍ）、ｙｑ（ｍ）を求める。
【００３６】
実施の形態１ではフィルタ係数は不変としているが、実施の形態２ではフィルタ係数が新たな血流情報を得る度に変化する。したがって、実施の形態２では、フィルタ係数が変化に応じて、フィルタ初期値を変化させることが好ましい。そこで、実施の形態２では、フィルタ係数生成部１１４で生成されたフィルタ係数をフィルタ初期値生成部１１１に出力する。フィルタ初期値生成部１１１は、フィルタ係数生成部１１４から受け取ったフィルタ係数と、クラッタ成分推定値１１０で推定されたクラッタ近似値の振幅、初期位相、位相速度をもとに、フィルタ初期値を生成してウォールフィルタ１０５を設定する。
【００３７】
フィルタ係数生成部１１４は、予め決められた方法に基づいてカットオフ周波数からクラッタ係数を計算してもよいが、カットオフ周波数とフィルタ係数を相互に関連付けたテーブルを用意（記憶）しておき、このテーブルに基づいてフィルタ係数を選択してもよい。
【００３８】
このように、実施の形態２に係る超音波ドプラ血流計１００Ｂによれば、心臓のようなクラッタおよび血流の位相速度の変化が激しい部位において、つねにクラッタのみを除去し、血流のみを抽出することができ、精度の高いデータを提供することができる。
【００３９】
実施の形態３
【００４０】
図４は本発明の実施の形態３に係る超音波ドプラ血流計１００Ｃのブロック図である。図示する超音波ドプラ血流計１００Ｃは、実施の形態２の超音波ドプラ血流計１００Ｂ（図３参照）と比べてみると明らかなように、クラッタ判定部１１４が付加されている点に特徴を有する。
【００４１】
超音波診断法にあっては、体内から返ってくるエコー信号において、多くの場合はクラッタの振幅が血流の振幅の数十倍もしくは百数十倍と大きい。しかし、例えば心臓の腔内のような血液のみが存在している部位からのエコー信号では、血流の振幅がクラッタの振幅より大きくなることがある。このような場合において、実施の形態２のような方法でクラッタ信号を同定し、除去を行なうと、振幅の大きな血流信号がクラッタとみなされてしまい、除去されてしまうおそれがある。
【００４２】
したがって、実施の形態３では、クラッタ判定部１１５においてクラッタと血流の判別を行なう。具体的には、クラッタ判定部１１５は、振幅の閾値を設定し、クラッタ成分推定部１１０で近似された正弦波の振幅がこの閾値を超える場合にある基準以上の振幅の場合は、推定された正弦波をクラッタとみなして以後の信号処理を行なう。しかし、推定された正弦波の振幅が閾値以下の振幅の場合は血流とみなし、上述したフィルタ初期値生成部１１１、カットオフ周波数設定部１１３、フィルタ係数生成部１１４における処理の実行を禁止する。
【００４３】
このように構成された実施の形態３に係る超音波ドプラ血流計１００Ｃによれば、心腔内のような血流信号がクラッタ信号より振幅が大きい場合にも血流信号を抽出でき、精度の高いデータを提供することができる。
【００４４】
実施の形態４
【００４５】
図５は本発明の実施の形態４に係る超音波ドプラ血流計１００Ｄのブロック図である。図示する超音波ドプラ血流計１００Ｄは、実施の形態３の超音波ドプラ血流計１００Ｃ（図４参照）と比べてみると明らかなように、クラッタ情報蓄積部１１５が付加されている点に特徴を有する。
【００４６】
具体的に説明すると、上述した実施の形態３では、振幅の閾値を設けることで、エコー信号中に含まれる最も大きな振幅の信号がクラッタであるか、血流であるかを判断したが、実施の形態４では、閾値を設定するための手段としてクラッタ情報蓄積部１５が設けられている。
【００４７】
超音波ドプラ血流計においては、体内を順次、走査していく過程において、近傍の部位の走査がすでに終わっており、クラッタや血流の情報が既に得られ、また、同じ部位を既に走査してその際にデータが取得されていることが通常である。そこで、実施の形態４の超音波ドプラ血流計１００Ｄでは、既に取得された情報をクラッタ情報蓄積部１１５に蓄積しておき、この蓄積された情報を用いてクラッタ判定部１１４又はクラッタ情報蓄積部１１５が閾値を設定する。
【００４８】
ところで、超音波ドプラ血流計の全走査範囲において、すべての範囲でクラッタが血流より小さいということは考えにくい。クラッタが血流より大きく、かつ、そのレベルは血流信号ではありえない振幅になるのが常である。そこで、現在走査している部位の近傍においてクラッタが血流より大きい部位を探すことで、クラッタと血流の振幅の敷居値を設定することができる。また、クラッタの振幅について空間的な分布を求めることで、現在走査している場所がどの程度のクラッタがあるかを予測することもできる。同じ部位について言えば、以前その部位において、例えば心臓の弁が横切るなどして、大振幅の信号が入力した経緯があるならば、それを基準にクラッタと血流の敷居値を設けることができる。したがって、実施の形態４によれば、蓄積された情報をもとに振幅の閾値が設定されるので、精度の良い血流情報を得ることができる。
【図面の簡単な説明】
【００４９】
【図１】本発明の実施の形態１に係る超音波ドプラ血流計の構成を示すブロック図。
【図２】図１に示す超音波ドプラ血流計のウォールフィルタの過渡応答を防ぐ構成部分を示す図。
【図３】本発明の実施の形態２に係る超音波ドプラ血流計の構成を示すブロック図。
【図４】本発明の実施の形態３に係る超音波ドプラ血流計の構成を示すブロック図。
【図５】本発明の実施の形態４に係る超音波ドプラ血流計の構成を示すブロック図。
【図６】従来の超音波ドプラ血流計の構成を示すブロック図。
【図７】図６の超音波ドプラ血流計に含まれるウォールフィルタの構成を示す図。
【符号の説明】
【００５０】
１，１０１…送信部、２，１０２…探触子、３，１０３…受信部、４，１０４…位相検波部、５，１０５，１０５Ａ，１０５Ｂ…ウォールフィルタ、６、１０６…血流情報演算部、７，１０７…包絡線検波部、８，１０８…ＤＳＣ、９，１０９…モニタ、１１０…クラッタ成分推定部、１１１…フィルタ初期値生成部、１１２…カットオフ周波数設定部、１１３…フィルタ係数生成部、１１５…クラッタ判定部、１１６…クラッタ情報蓄積部。
【手続補正３】
【補正対象書類名】図面
【補正対象項目名】全図
【補正方法】変更
【補正の内容】
【図１】
2006088094.tif 000035
【図２】
2006088094.tif 000036
【図３】
2006088094.tif 000037
【図４】
2006088094.tif 000038
【図５】
2006088094.tif 000039
【図６】
2006088094.tif 000040
【図７】
2006088094.tif 000041
【国際調査報告】
2006088094.tif 000042
2006088094.tif 000043
2006088094.tif 000044
2006088094.tif 000045