(57)【特許請求の範囲】
【請求項１】
被検体に対して超音波を送信し、得られたエコー信号を周波数解析して、周波数スペクトラムを生成し、この周波数スペクトラムを時間波形として縦軸を周波数軸、横軸を時間軸として表示する超音波診断装置において、
前記周波数スペクトラムの周波数分解能と時間分解能とを両者の比を一定に保ったままで変化させる比一定モードと、時間分解能を固定したままで周波数分解能を変化させる時間分解能優先モードと、周波数分解能を固定したままで時間分解能を変化させる周波数分解能優先モードとを選択することができるように構成され、
前記時間分解能ΔＸと、前記周波数分解能ΔＹはそれぞれ、
ΔＸ＝（ＤX ×Ｎ×Ｋ×Ｍ）／（ｆPRF ×Ｓ）
ΔＹ＝ＤY ／（Ｎ×Ｗw ）
ＤX ；ピリオドグラムの表示範囲の横軸（時間軸）方向の画素数ＤY ；ピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数Ｎ；周波数解析のデータ数Ｋ；交互段数Ｍ；間引き率ｆPRF ；レート周波数Ｓ；スクロール速度Ｗw ；窓重みにより定義され、
前記時間分解能優先モードでは、周波数分解能に関わるピリオドグラムの表示範囲の周波数軸方向の画素数ＤY 、周波数解析のデータ数Ｎ、窓重みＷw のうち、少なくとも１つが変わって、周波数分解能が変動しても、データ数Ｎ以外のピリオドグラムの表示範囲の時間軸方向の画素数ＤX 、交互段数Ｋ、間引き率Ｍ、レート周波数ｆPRF 、スクロール速度Ｓのうち少なくとも１つを変えて、時間分解能を所定値に維持することを特徴とする超音波診断装置。
【請求項２】
被検体に対して超音波を送信し、得られたエコー信号を周波数解析して、周波数スペクトラムを生成し、この周波数スペクトラムを時間波形として縦軸を周波数軸、横軸を時間軸として表示する超音波診断装置において、
前記周波数スペクトラムの周波数分解能と時間分解能とを両者の比を一定に保ったままで変化させる比一定モードと、時間分解能を固定したままで周波数分解能を変化させる時間分解能優先モードと、周波数分解能を固定したままで時間分解能を変化させる周波数分解能優先モードとを選択することができるように構成され、
前記時間分解能ΔＸと、前記周波数分解能ΔＹはそれぞれ、
ΔＸ＝（ＤX ×Ｎ×Ｋ×Ｍ）／（ｆPRF ×Ｓ）
ΔＹ＝ＤY ／（Ｎ×Ｗw ）
ＤX ；ピリオドグラムの表示範囲の横軸（時間軸）方向の画素数ＤY ；ピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数Ｎ；周波数解析のデータ数Ｋ；交互段数Ｍ；間引き率ｆPRF ；レート周波数Ｓ；スクロール速度Ｗw ；窓重みにより定義され、
前記周波数分解能優先モードでは、時間分解能に関わるピリオドグラムの表示範囲の時間軸方向の画素数ＤX 、交互段数Ｋ、間引き率Ｍ、レート周波数ｆPRF 、スクロール速度Ｓのうち少なくとも１つが変わって、時間分解能が変動しても、ピリオドグラムの表示範囲の周波数軸方向の画素数ＤY 、データ数Ｎ及び窓重みＷw は変えないで、また両分解能に関わるデータ数Ｎが変わった場合には、ピリオドグラムの表示範囲の周波数軸方向の画素数ＤY と窓重みＷw との少なくとも一方を変えて周波数分解能を一定値に維持することを特徴とする超音波診断装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、ドプラ効果により偏移した周波数のスペクトラムを時間波形（ピリオドグラム）として表示する超音波診断装置に関する。
【０００２】
【従来の技術】
上述の周波数スペクトラムの描画性、つまりピリオドグラムの表示の細かさを定義する周波数分解能と時間分解能により評価され得る。これら周波数分解能と時間分解能は、レート周波数ｆPRF やＦＦＴのデータ数Ｎ、さらにはスクロール速度Ｓ等のパラメータにより決まるもので、従来では、例えばレート周波数ｆPRF が変えられたときには、時間分解能が変動するので、ピリオドグラムとしては非常に見難いという問題があった。
【０００３】
この問題を解決するために、あるパラメータが変わった或いは変えられたときには、ＦＦＴのデータ数Ｎや窓重みＷw を調整して、時間分解能ΔＸに対する周波数分解能ΔＹの比（縦横比；ΔＹ／ΔＸ）を一定に維持して、ピリオドグラムの表示の安定化を図るような技術が開発され、現在では主流派になっている。
【０００４】
しかし、この縦横比一定化技術では、例えば低流速検出能を向上して腹部等の比較的遅い血流を精度よく観測するために、レート周波数ｆPRF を下げると、時間分解能が大きくなってしまうが、それに伴ってＦＦＴのデータ数Ｎや窓重みＷw を変えて周波数分解能を大きくするように機能するため、ピリオドグラムを構成する最小単位（縦方向の画素数が周波数分解能、横方向の画素数が時間分解能で定義される）が非常に大きく粗くなって、微妙な速度変化や時間的な細かな応答等が捉えにくくなってしまうという問題があった。
【０００５】
【発明が解決しようとする課題】
本発明の目的は、観察者の要望にできる限り沿って周波数分解能や時間分解能を調整できる超音波診断装置を提供することにある。
【０００６】
【課題を解決するための手段】
本発明の第１局面は、被検体に対して超音波を送信し、得られたエコー信号を周波数解析して、周波数スペクトラムを生成し、この周波数スペクトラムを時間波形として縦軸を周波数軸、横軸を時間軸として表示する超音波診断装置において、前記周波数スペクトラムの周波数分解能と時間分解能とを両者の比を一定に保ったままで変化させる比一定モードと、時間分解能を固定したままで周波数分解能を変化させる時間分解能優先モードと、周波数分解能を固定したままで時間分解能を変化させる周波数分解能優先モードとを選択することができるように構成され、前記時間分解能ΔＸと、前記周波数分解能ΔＹはそれぞれ、ΔＸ＝（ＤX ×Ｎ×Ｋ×Ｍ）／（ｆPRF ×Ｓ）、ΔＹ＝ＤY ／（Ｎ×Ｗw ）、ＤX ；ピリオドグラムの表示範囲の横軸（時間軸）方向の画素数ＤY ；ピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数Ｎ；周波数解析のデータ数Ｋ；交互段数Ｍ；間引き率ｆPRF ；レート周波数Ｓ；スクロール速度Ｗw ；窓重みにより定義され、前記時間分解能優先モードでは、周波数分解能に関わるピリオドグラムの表示範囲の周波数軸方向の画素数ＤY 、周波数解析のデータ数Ｎ、窓重みＷw のうち、少なくとも１つが変わって、周波数分解能が変動しても、データ数Ｎ以外のピリオドグラムの表示範囲の時間軸方向の画素数ＤX 、交互段数Ｋ、間引き率Ｍ、レート周波数ｆPRF 、スクロール速度Ｓのうち少なくとも１つを変えて、時間分解能を所定値に維持することを特徴とする。
本発明の第２局面は、被検体に対して超音波を送信し、得られたエコー信号を周波数解析して、周波数スペクトラムを生成し、この周波数スペクトラムを時間波形として縦軸を周波数軸、横軸を時間軸として表示する超音波診断装置において、前記周波数スペクトラムの周波数分解能と時間分解能とを両者の比を一定に保ったままで変化させる比一定モードと、時間分解能を固定したままで周波数分解能を変化させる時間分解能優先モードと、周波数分解能を固定したままで時間分解能を変化させる周波数分解能優先モードとを選択することができるように構成され、前記時間分解能ΔＸと、前記周波数分解能ΔＹはそれぞれ、ΔＸ＝（ＤX ×Ｎ×Ｋ×Ｍ）／（ｆPRF ×Ｓ）、ΔＹ＝ＤY ／（Ｎ×Ｗw ）、ＤX ；ピリオドグラムの表示範囲の横軸（時間軸）方向の画素数ＤY ；ピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数Ｎ；周波数解析のデータ数Ｋ；交互段数Ｍ；間引き率ｆPRF ；レート周波数Ｓ；スクロール速度Ｗw ；窓重みにより定義され、前記周波数分解能優先モードでは、時間分解能に関わるピリオドグラムの表示範囲の時間軸方向の画素数ＤX 、交互段数Ｋ、間引き率Ｍ、レート周波数ｆPRF 、スクロール速度Ｓのうち少なくとも１つが変わって、時間分解能が変動しても、ピリオドグラムの表示範囲の周波数軸方向の画素数ＤY 、データ数Ｎ及び窓重みＷw は変えないで、また両分解能に関わるデータ数Ｎが変わった場合には、ピリオドグラムの表示範囲の周波数軸方向の画素数ＤY と窓重みＷw との少なくとも一方を変えて周波数分解能を一定値に維持することを特徴とする。
【０００７】
【発明の実施の形態】
以下、図面を参照して、本発明による超音波診断装置を好ましい実施形態により説明する。図１に本実施形態に係る超音波診断装置の構成を示す。超音波プローブ１の先端部分には、電気信号と音響信号とを相互変換するための複数の圧電素子が配列されている。この超音波プローブ１には、送受信回路（Ｔ＆Ｒ）２が接続されている。この送受信回路２は、送信部分と受信部分と直交検波部分とからなり、送信部分では、まず、クロック発生器からのクロックを分周して、４ｋＨｚ乃至５０ｋＨｚ程度のレートパルスを生成する。このレートパルスの周波数は、レート周波数とかパルス繰り返し周波数とか呼ばれている。ここでは、レート周波数と称するものとし、記号“ｆPRF ”で表すものとする。このレートパルスを適当に遅延させた後、これをトリガとしてパルサからプローブ１の圧電素子に電圧パルスが印加されるようになっている。この印加によりプローブ２の圧電素子は機械的に振動し、超音波パルスを発生する。
【０００８】
この超音波は被検体内部を伝播し、その途中にある音響インピーダンスの不連続面で反射し、プローブ１に返ってきて、圧電素子を機械的に振動する。これにより、圧電素子からは微弱な電気信号が発生する。この電気信号は、受信部分に取り込まれ、プリアンプで増幅され、アナログディジタル変換器でディジタル信号に変換されて、いわゆるディジタルビームフォーマで遅延相当の処理を経て加算される。この受信部分からの出力信号には、血球等の移動体のドプラ効果により偏移した周波数成分が含まれており、これを直交検波部分で検波する。
【０００９】
この送受信回路２の出力信号、つまりドプラ直交検波信号には、主に血球等の速い移動体での反射により周波数偏移を受けた高周波成分と、主に心臓壁等の遅い移動体での反射により周波数偏移を受けた低周波成分、一般的にはクラッタ成分とが含まれているので、高域通過型のクラッタフィルタ３で低周波成分を減衰し、高周波成分を抽出し、このクラッタフィルタ３で抽出された高周波成分をオーディオ回路６でオーディオ処理して、スピーカ７からドプラ音として観察者が聞けるようになっている。
【００１０】
また、このクラッタフィルタ３で抽出された高周波成分を高速フーリエ変換プロセッサ（ＦＦＴ）４で周波数分析し、その周波数スペクトラムを表示装置５で時間波形、つまりピリオドグラムとして表示するようになっている。
【００１１】
ところで制御装置（ＣＰＵ）８は、周波数分解能や時間分解能を、３種類のモードで選択的に調整することができるようになっている。この３種類のモードとは、図２（ａ）に示すような時間分解能に対する周波数分解能の比を一定に維持したままでこれら分解能を調整することのできる縦横比一定モードと、図２（ｂ）に示すような時間分解能を所定値に保ったままで、周波数分解能を調整できる時間分解能優先モードと、図２（ｃ）に示すような周波数分解能を所定値に保ったままで、時間分解能を調整できる周波数分解能優先モードとがある。
【００１２】
これら各モードについて順番に説明する。まず、時間分解能ΔＸ、周波数分解能ΔＹの計算方法について説明する。時間分解能ΔＸ、周波数分解能ΔＹは、次の式（１），（２）により計算される。
【００１３】
ΔＸ＝（ＤX ×Ｎ×Ｋ×Ｍ）／（ｆPRF ×Ｓ）  …（１）
ΔＹ＝ＤY ／（Ｎ×Ｗw ）                    …（２）
ＤX ；ピリオドグラムの表示範囲の横軸（時間軸）方向の画素数
ＤY ；ピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数
Ｎ；ＦＦＴのデータ数
Ｋ；交互段数
Ｍ；間引き率
ｆPRF ；レート周波数
Ｓ；スクロール速度
Ｗw ；窓重み
なお、窓重みＷw とは、次の通りである。周知の通り、ＦＦＴのデータ数Ｎは、２のべき乗という制約があるが、多くの場合、目的の周波数分解能を実現するための計算上のデータ数は２のべき乗にはならない。そこで、ＦＦＴのデータ数Ｎには、計算上のデータ数以上でそれに最も近い２のべき乗を選び、その選んだＦＦＴのデータ数Ｎと、計算上のデータ数との不一致を、ＦＦＴの窓をかけることで埋めるようにしているが、この窓の大きさ（窓幅）を決めるパラメータが、窓重みＷw である。例えば、計算上のデータ数が、３２０であるとき、ＦＦＴのデータ数Ｎとしては、５１２が選ばれ、この場合、窓重みＷw は、５／８に設定される
図３には、縦横比一定モードを実現するための制御方法の概念を示しており、レート周波数ｆPRF 、スクロール速度Ｓ、交互段数Ｋ、間引き計数Ｍ、ＤY が与えられると、観察者の所望の縦横比Ｐr に維持するためのＦＦＴのデータ数Ｎと窓重みＷw を計算するようになっている。つまり、例えばレート周波数ｆPRF が変わったとき（時間分解能が変わったとき）、それに応じてＦＦＴのデータ数ＡA を調整して周波数分解能を変えて、縦横比を一定にするようにることになる。このデータ数ＡA は、次の式（３）により、求められる。
【００１４】
ＡA ＝（ｆPRF ×ＤY ×Ｓ）／（Ｐr ×Ｋ×Ｍ×ＤX ）    …（３）
こうして計算されたＦＦＴのデータ数ＡA は、上述したように、ＦＦＴで扱える２のべき乗には必ずしもならない。そこで、上述したように、計算上のデータ数ＡA 以上でそれに最も近い２のべき乗を選び、その選んだＦＦＴのデータ数Ｎと、両者間の不一致を埋めるための窓重みＷw とを選択する。この処理は、実際には、ＲＯＭテーブルにより実現されていることが多い。
【００１５】
このような処理で、レート周波数等のパラメータが変わっても、時間分解能ΔＸに対する周波数分解能ΔＹの比を一定に保つことができる。
次に、時間分解能優先モードでは、周波数分解能に関わるパラメータ、つまりピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数ＤY 、ＦＦＴのデータ数Ｎ、窓重みＷw のうち、少なくとも１つが変わって、周波数分解能が変動しても、時間分解能を、観察者が所望する一定値に固定するような処理が行われ、つまり、両分解能に関わるデータ数Ｎが変わったとき、上記式（１）において、データ数Ｎ以外のパラメータ、つまりピリオドグラムの表示範囲の横軸（時間軸）方向の画素数ＤX 、交互段数Ｋ、間引き率Ｍ、レート周波数ｆPRF 、スクロール速度Ｓのうち少なくとも１つ、例えばスクロール速度Ｓを変えて、時間分解能を一定値に維持するものである。
【００１６】
また、周波数分解能優先モードでは、時間分解能に関わるパラメータ、つまり両分解能に関わるデータ数Ｎ以外のピリオドグラムの表示範囲の横軸（時間軸）方向の画素数ＤX 、交互段数Ｋ、間引き率Ｍ、レート周波数ｆPRF 、スクロール速度Ｓのうち少なくとも１つが変わって、時間分解能が変動しても、周波数分解能を、観察者が所望する一定値に固定するように、ピリオドグラムの表示範囲の縦軸（周波数軸）方向の画素数ＤY 、データ数Ｎ及び窓重みＷw は変えないで固定したままにして、また両分解能に関わるデータ数Ｎが変わった場合には、ピリオドグラムの表示範囲の周波数軸方向の画素数ＤY と窓重みＷw との少なくとも一方を変えて周波数分解能を一定値に維持するものである。
【００１７】
図４には、真のピリオドグラムに対して、実際に表示されるピリオドグラムがモードによって見え方がどのように変化するかについて示している。
このような３種類のモードが用意されているので、例えば当初、縦横比一定モードを使っていて、時間分解能が拡大又は縮小過多で観測し難いときには、時間分解能優先モードを選択し、また周波数分解能が拡大又は縮小過多で観測し難いときには、周波数分解能優先モードを選択することができるので、微妙な速度変化を観測し難くなったり、また逆に単位領域が細か過ぎて時間的な応答を観測し難くなったといった問題を解決することができる。
本発明は、上述してきたような実施形態に限定されることなく、種々変形して実施可能であることは言うまでもない。
【００１８】
【発明の効果】
本発明によれば、３種類のモードを使い分けて、例えば縦横比一定モードでは、ピリオドグラムが見え難い、例えば微妙な速度変化を観測し難くなった場合には時間分解能を固定するモードや周波数分解能を固定するモードを選択して、ピリオドグラムを観察し易くすることができる。
【図面の簡単な説明】
【図１】本発明の好ましい実施形態に係る超音波診断装置の構成を示すブロック図。
【図２】本実施形態による３種類の表示モードに関する説明図。
【図３】縦横比一定モードのローカル処理の概念図。
【図４】図２の各モードによるピリオドグラムの見え方を真のピリオドグラムと比較して示す図。
【符号の説明】
１…超音波プローブ、
２…送受信回路（Ｔ＆Ｒ）、
３…クラッタフィルタ、
４…高速フーリエ変換プロセッサ（ＦＦＴ）、
５…表示装置、
６…オーディオ回路、
７…スピーカ、
８…制御装置、
９…コンソール。
【図１】
0004191812.tif 000002
【図２】
0004191812.tif 000003
【図３】
0004191812.tif 000004
【図４】
0004191812.tif 000005