JP3691874B2 - Ultrasonic diagnostic equipment - Google Patents

Ultrasonic diagnostic equipment Download PDF

Info

Publication number
JP3691874B2
JP3691874B2 JP20219595A JP20219595A JP3691874B2 JP 3691874 B2 JP3691874 B2 JP 3691874B2 JP 20219595 A JP20219595 A JP 20219595A JP 20219595 A JP20219595 A JP 20219595A JP 3691874 B2 JP3691874 B2 JP 3691874B2
Authority
JP
Japan
Prior art keywords
doppler
data
signal
range gate
time
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 - Lifetime
Application number
JP20219595A
Other languages
Japanese (ja)
Other versions
JPH0947451A (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.)
Toshiba Corp
Original Assignee
Toshiba Corp
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 Toshiba Corp filed Critical Toshiba Corp
Priority to JP20219595A priority Critical patent/JP3691874B2/en
Publication of JPH0947451A publication Critical patent/JPH0947451A/en
Application granted granted Critical
Publication of JP3691874B2 publication Critical patent/JP3691874B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Description

【0001】
【産業上の利用分野】
本発明は、受信信号を周波数解析することによりドプラスペクトラムを得る超音波診断装置に関する。
【0002】
【従来の技術】
超音波は医学的な見地から様々に応用されているが、その主流は超音波パルス反射法を用いて生体の軟部組織の断層像を得る超音波診断装置である。この超音波診断装置は無侵襲検査法で、組織の断層像を表示するものであり、X線診断装置、X線CT装置、MRIおよび核医学診断装置などの他の診断装置に比べて、リアルタイム表示が可能、装置が小型で安価、X線などの被曝がなく安全性が高い、ドプラ効果を利用してドプラスペクトラムや血流イメージングが可能である等の独占的な特徴を有している。ここでは、ドプラスペクトラムに着目する。ドプラスペクトラムの表示としては、縦軸にドプラ周波数(速度)、横軸に時間、輝度を各ドプラ周波数成分の強さ(パワー)でMモードのようにスクロール表示が使用される。
【0003】
このような超音波診断装置のドプラスペクトラム機能において、感度向上を目的として、バースト波が実用化されている。バースト波とは複数の超音波パルスの系列をいう。バースト波の時間幅をバースト長という。バースト波の場合、パルス波と違って、レンジゲートマーカ内のドプラ信号の積分処理に工夫が必要である。図7(a)に示す従前のケースでは、レンジゲートマーカの深さ及び長さに相当する期間t1 〜t2 に受信されるドプラ信号だけを積分処理に取り込む。なお、t1 は、バースト波の送信時刻t0 から、レンジゲートマーカの深さに相当する時間(往復を考慮してレンジゲートマーカの深さの2倍の距離を超音波が伝搬するのに要する時間)を経過した時刻であり、t2 はt1 から、レンジゲートマーカの長さに相当する時間(レンジゲートマーカの長さの2倍の距離を超音波が伝搬するのに要する時間)を経過した時刻である。この従前のケースでは、バースト波のバースト長τを考慮していない、つまり同一反射源からの反射波がバースト長τの期間、継続的に受信されることを考慮していないので、積分信号の感度分布を見て理解されるように、レンジゲートマーカよりτに相当する浅い領域からのドプラ信号が積分処理に取り込まれてしまう。
【0004】
この問題点を解決するように工夫したケースを図7(b)に示す。このケースでは、バースト波のバースト長τは、レンジゲートマーカの長さに相当する時間の半分に意図的に調整される。そして、レンジゲートマーカの深さ及び長さに相当する期間t1 〜t2 の後半の期間t3 〜t2 に受信されるドプラ信号だけを積分処理に取り込む。結果的に、積分期間とバースト長τとは一致する。このケースでは、積分信号の感度分布を見て理解されるように、レンジゲートマーカ外のドプラ信号は積分処理に取り込まれない。
【0005】
しかし、図7(b)のケースには次のような問題点がある。距離分解能の著しい低下を回避するために、バースト波のバースト長τは制限される。したがって、レンジゲートマーカが非常に長く設定されたとき、図7(c)に示すように、積分期間が最長のバースト長τmax より長くなってしまうケースが起こり得る。この場合、レンジゲートマーカ内にも関わらず、そこからのドプラ信号が積分処理に全く取り込まれない、いわゆる不感領域が存在してしまう。
【0006】
また、ドプラスペクトラム機能には、次のような問題もある。パルス波ドプラ(PWドプラ)では、レンジゲートマーカの長さを変化させると、積分期間が変化するが、この積分期間の変化は、実質的にゲインが変化することを意味する。また、連続波ドプラ(CWドプラ)では、ドプラ信号を高域通過型フィルタ(HPF )にかけるが、レート周波数PRF を変化させると、この高域通過型フィルタの入力帯域に対する出力帯域の比が変化するが、この入出力帯域比の変化は、実質的にゲインが変化することを意味する。さらに、PWドプラ、CWドプラの両者に言えることとして、周波数解析としての高速フーリエ変換(FFT )のデータ数の変化は、実質的にゲインが変化することを意味する。したがって、ゲインの変化に伴って、ノイズレベルも変化してしまうという問題があった。
【0007】
また、PWドプラでは、実質的にレート周波数を低くして、低速検出能を向上させるために、交互スキャンと呼ばれる機能が実用化されている。通常、ある1本のラスタに対して超音波送受信動作をレート周波数で所定回数繰り返しながら、順次次のラスタに移っていくというスキャンが行われる。これに対し、交互スキャンでは、ある1本のラスタに対して超音波送受信動作を1回行い、同じラスタに対する次回の超音波送受信動作までに間に、他の何本かのラスタに対して1回ずつ超音波送受信動作を行うものであり、1本のラスタに関して見ると、超音波送受信動作の間隔が長くなり、実質的にレート周波数が低くなるというものである。なお、同じラスタに対する前後2回の超音波送受信動作の間隙に行う超音波送受信動作の回数を交互段数として定義する。交互段数の増加は、周波数解析のデータ数が一定であれば、データ取り込み時間の延長を意味する。データ取り込み時間の延長は、ドプラスペクトラムの表示上で周波数分解能と時間分解能との比が変化することにつながる。これは好ましくないので、現在の多くの装置には、交互段数が増加すると、それに応じて周波数解析のデータ数を減少させて、周波数分解能と時間分解能との比を一定に保つ機能が装備されている。しかし、周波数解析のデータ数の減少は、S/Nを劣化(低下)させてしまう。
【0008】
【発明が解決しようとする課題】
本発明の目的は、パルス波ドプラでは積分期間の変化、連続波ドプラでは高域通過型フィルタの入力帯域に対する出力帯域の比の変化、周波数解析としての高速フーリエ変換(FFT )のデータ数の変化等により実質的にゲインが変化しても、ノイズレベルは変化しない超音波診断装置を提供することである
【0009】
【課題を解決するための手段】
本発明による超音波診断装置は、超音波パルスを被検体内に送信し、反射波を受信することにより受信信号を得る送受信手段と、レンジゲートマーカの深さ及び長さに応じた特定の期間に前記送受信手段から出力される前記受信信号を積分することにより積分信号を得る手段と、前記レンジゲートマーカ内のドプラスペクトラムを得るために前記積分信号を周波数解析する手段と、前記送受信手段と前記積分手段との間に設けられ、前記積分手段による積分期間をΔT RG として1/ΔT RG 1/2 に応じた重み係数を前記受信信号に乗算する手段とを具備する。
本発明によれば、積分期間が変化すると、それに伴って実質的にゲインが変化してノイズレベルも変化してしまうが、積分期間に応じた重み係数を受信信号に乗算するので、積分期間の変化に伴う実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を連続波として被検体内に送信し、反射波を受信することにより受信信号を得る手段と、前記受信信号の特定の周波数帯域を減衰するフィルタ手段と、前記フィルタ手段の出力帯域に対する前記フィルタ手段の入力帯域の比に応じた重み係数を前記フィルタ手段の出力信号に乗算する手段と、ドプラスペクトラムを得るために前記フィルタ手段の出力信号を周波数解析する手段とを具備する。
本発明によれば、例えばレート周波数が変化すると、それに伴ってフィルタ手段の出力帯域に対する入力帯域の比が変化して実質的にゲインが変化してノイズレベルも変化してしまうが、この比に応じた重み係数をフィルタ手段の出力信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を被検体内に送信し、反射波を受信することにより受信信号を得る手段と、ドプラスペクトラムを得るために前記受信信号をデータ数で周波数解析する手段と、前記データ数をNとして重み係数1/N 1/2 を前記受信信号に乗算する手段とを具備する。
本発明によれば、周波数解析のデータ数が変化すると実質的にゲインが変化してノイズレベルも変化してしまうが、このデータ数に応じた重み係数を受信信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
【0019】
【実施の形態】
以下に、本発明の実施の形態を図面を参照して説明する。
(第1の実施の形態)
図1は第1の実施の形態による超音波診断装置のブロック図である。本実施の形態による超音波診断装置は、超音波プローブ1、送受信系(T/R)2、スペクトラムドプラユニット3、ディジタルスキャンコンバータ(DSC)4、モニタ5、コンソール6とから構成される。
【0020】
超音波プローブ1は、圧電セラミック等の機械運動/電気信号の可逆的変換素子としての複数の圧電素子が一次元に配列された圧電素子アレイをその先端に装備してなる。
【0021】
送受信系2は、送信系と受信系とを有する。送信系は、パルス発生器、送信遅延回路、パルサとからなる。パルス発生器は例えば5KHzのレート周波数fr (周期;1/fr 秒)でレートパルスを繰り返し発生する。このレートパルスはチャンネル数に分配され、送信遅延回路に送られる。送信遅延回路は、超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間を各レートパルスに与える。パルサは、送信遅延回路からレートパルスを受けた時刻で超音波プローブ1の圧電素子にチャンネル毎に送信電圧を印加する。これにより超音波ビームが被検体に送信される。
【0022】
被検体内の音響インピーダンスの不連続面(反射源)で反射した反射波は超音波プローブ1で受信され、電気信号に変換される。超音波プローブ1から素子毎に出力される電気信号は、チャンネル毎に纏められて受信系に取り込まれる。受信系は、プリアンプ、受信遅延回路、加算器、位相検波回路から構成される。電気信号は、プリアンプで増幅され、受信遅延回路により受信指向性を決定するのに必要な遅延時間を与えられ、そして加算器で加算される。この加算により受信指向性に応じた方向からの反射成分が強調された受信信号が得られる。この受信信号は位相検波回路で位相検波される。これによりドプラ効果による偏移周波数のドプラ信号が得られる。ドプラ信号はスペクトラムドプラユニット3に供給される。
【0023】
スペクトラムドプラユニット3は、特定の期間に供給されるドプラ信号のみをレンジゲートマーカユニット8で取り出し、そして比較的動きの速い高周波の血流成分と比較的動きの遅い血管壁や弁等の低周波のクラッタ成分とを高域通過型フィルタ9で分離し、続いて高域通過型フィルタ9を通過したドプラ信号を高速フーリエ変換処理回路(FFT )10でリアルタイムで周波数解析することにより、各周波数成分の強さ(パワー)を表すドプラスペクトラムデータを生成する。レンジゲートマーカユニット8でドプラ信号を取り出す特定の期間は、コンソール6を介してBモード画像上に設定されるレンジゲートマーカマーカの深さ及び長さに基づいてコントローラ7で制御される。
【0024】
高速フーリエ変換処理回路10で得られたドプラスペクトラムデータは、ディジタルスキャンコンバータ4、図示しないディジタルアナログコンバータを介してモニタ5にドプラスペクトラムとして表示される。ドプラスペクトラムの表示としては、縦軸にドプラ周波数(速度)、横軸に時間、輝度を各ドプラ周波数成分の強さ(パワー)として、Mモードのようにスクロール表示が使用される。
【0025】
レンジゲートマーカユニット8は、ドプラ信号を乗算器(増幅器)11で振幅変調し、振幅変調されたドプラ信号を特定の期間だけ積分器12で積分する。乗算器11による振幅変調データは、関数発生ROM13から供給される。積分器12が積分を実行する特定の期間は、RG幅カウンタ14及びRGスタート位置カウンタ15により調整される。
【0026】
RGスタート位置カウンタ15は、コントローラ7からレンジゲートマーカの深さに対応するRG位置データの供給を受ける。RGスタート位置カウンタ15は、スキャンの開始からクロックパルスの計数(アップカウント)を開始し、計数値がRG位置データが表す値に達した時、スタートパルスを出力する。
【0027】
RG幅カウンタ14は、コントローラ7からレンジゲートマーカの長さに対応するRG幅データの供給を受ける。RG幅カウンタ14は、RGスタート位置カウンタ15からスタートパルスを受けた時刻でクロックパルスの計数(アップカウント)を開始し、計数値がRG幅データが表す値に達する時まで計数動作を継続する。したがって、RG幅カウンタ14からは、レンジゲートマーカの深さ及び長さに対応する特定の期間だけ、計数信号が出力される。
【0028】
積分器12は、RG幅カウンタ14から計数信号が出力されている期間に限り積分を実行する。これにより、特定の期間内の振幅変調されたドプラ信号だけが積分信号として取り出される。
【0029】
このようにコントローラ7からのRG位置データ及びRG幅データに、レンジゲートユニット8で取り出されるドプラ信号の期間(積分期間)が依存される。換言するとコントローラ7はRG位置データ及びRG幅データを調整することにより積分期間の開始時刻及びその時間幅を自由に変えることができることが理解されよう。
【0030】
関数発生ROM13は、RG幅カウンタ14から計数信号が出力されている特定の期間だけ係数を乗算器11に出力する。ドプラ信号は乗算器11で関数発生ROM13からの係数に応じて振幅変調される。
【0031】
関数発生ROM13には、予め複数種類の重み係数が記憶されている。関数発生ROM13には、RG幅カウンタ14からの計数値信号と、コントローラ7からの関数選択信号とがアドレス信号として供給される。アドレス信号に応じたアドレス部の係数が関数発生ROM13から読み出される。コントローラ7は複数種類の関数選択信号の中から、コンソール6を介して選択された感度特性に対応する1つの関数選択信号を選択的に関数発生ROM13に供給する。ここでいう関数とは関数発生ROM13から読み出される係数の時間変化を意味する。第1の関数選択信号のもとでは、関数発生ROM13からは一定の係数(=1.0 )が読み出され、つまり振幅変調処理の無いが如く振る舞われる。第2の関数選択信号のもとでは、関数発生ROM13からは、例えばハミング窓関数を満足するように経時的に変化する係数が読み出される。
【0032】
次に本実施の形態の動作を説明する。ここで、以下で使用される用語を定義する。レンジゲートマーカの深さに相当する時間とは、超音波送受信の往復を考慮して、レンジゲートマーカの深さの2倍の距離を超音波が伝搬するのに要する時間として定義する。また、レンジゲートマーカの長さに相当する時間とは、レンジゲートマーカの長さの2倍の距離を超音波が伝搬するのに要する時間として定義する。
【0033】
コントローラ7は、第1のモードと第2のモードとを有する。第1のモードはレンジゲートマーカの長さに相当する時間の半分と、バースト波のバースト長τとが一致する場合に適用され、第2のモードはレンジゲートマーカの長さに相当する時間の半分より、バースト波のバースト長τ(最長のバースト長τmax )が短い場合に適用される。なお、図7(b)を参照して説明した通り、バースト波のバースト長τは、レンジゲートマーカの長さに相当する時間の半分に積極的に調整されるのではあるが、距離分解能の著しい低下を抑える目的であえてバースト長τに制限を与えている由により、設定されるレンジゲートマーカによっては、レンジゲートマーカの長さに相当する時間の半分より、バースト波のバースト長τ(最長のバースト長τmax )が短い場合が発生する。
【0034】
第1のモードでは、図7(b)を参照して説明した通り、バースト波のバースト長τはレンジゲートマーカの長さに相当する時間の半分に意図的に調整される。時刻t3 〜t2 の期間に供給されるドプラ信号だけがレンジゲートユニット8で積分処理に取り込まれる。時刻t3 は、時刻t1 からバースト長τだけ経過した時刻であり、時刻t1 はバースト波の送信時刻t0 からレンジゲートマーカの深さに相当する時間△Td を経過した時刻である。時刻t2 は時刻t1 からレンジゲートマーカの長さに相当する時間△TRGを経過した時刻、つまりレンジゲートマーカの最深部に相当する時刻である。積分期間は、レンジゲートマーカの長さに相当する時間△TRGの半分(△TRG/2)に設定される。第1のモードでは、積分信号の感度分布を見て理解されるように、レンジゲートマーカ外のドプラ信号は好ましく全く積分処理に取り込まれない。
【0035】
図2(a),(b)は第2のモードでの積分開始時刻及び積分期間の時間幅を示すタイムチャートであり、同図(a)は第1の関数選択信号に対応し、また同図(b)は第2の関数選択信号に対応する。
【0036】
第2のモードでは、図2(a),(b)に示すように、時刻t1 ´〜t2 ´の積分期間に供給されるドプラ信号だけがレンジゲートユニット8で積分処理に取り込まれる。時刻t1 ´は、時刻t1 からバースト長τの半分(τ/2)だけ経過した時刻であり、時刻t1 はバースト波の送信時刻t0 からレンジゲートマーカの深さに相当する時間△Td を経過した時刻である。時刻t2 ´は時刻t1 ´からレンジゲートマーカの長さに相当する時間△TRGを経過した時刻である。積分期間は、レンジゲートマーカの長さに相当する時間△TRGに設定される。第2のモードでは、積分信号の感度分布を見て理解されるように、レンジゲートマーカ内に不感領域は存在しない。
【0037】
このように第2のモードでは、レンジゲートマーカの深さに相当する時刻t1 からτ/2だけ遅延した時刻t1 ´で積分を開始し、レンジゲートマーカの長さに相当する時間△TRG、当該積分を継続することにより、レンジゲートマーカより浅い側及び深い側それぞれにτ/2に相当する範囲からのドプラ信号が積分処理に取り込まれてしまうが、図7(a)のケースと比較すると、そのレンジゲートマーカ外のドプラ信号の取り込み量(面積)は絶対的に減少する。また、図7(c)と比較して理解されるように、積分期間がバースト長τより長くなるケースに直面しても、レンジゲートマーカ内の不感領域を完全に排除できる。
【0038】
さらに図2(b)の積分信号の感度分布から理解されるように、関数発生ROM13から例えばハミング窓関数を満足するように経時的に変化する係数が読み出され、ドプラ信号に乗算される場合、感度分布が中央付近で高くなり、レンジゲートマーカより浅い側及び深い側それぞれにτ/2に相当する範囲の感度分布が抑制される。
(第2の実施の形態)
図3は第2の実施の形態による超音波診断装置のブロック図である。本実施の形態による超音波診断装置は、超音波プローブ21、送受信系(T/R)22、スペクトラムドプラユニット23、ディジタルスキャンコンバータ(DSC)24、モニタ25、コンソール26とから構成される。
【0039】
超音波プローブ21は、第1の実施の形態の超音波プローブ1に同じである。送受信系22は、第1の実施の形態の送受信系2が、超音波をパルス波(PW)又は連続波(CW)として選択的に発生させることが可能に改良され、アナログのドプラ信号をドプラデータにディジタル化して出力するように最後段にADCを追加される。広域通過型フィルタ30は、第1の実施の形態の広域通過型フィルタ9に同じである。高速フーリエ変換処理回路31は、第1の実施の形態の高速フーリエ変換処理回路10に同じである。ディジタルスキャンコンバータ24は、第1の実施の形態のディジタルスキャンコンバータ4に同じである。モニタ25は、第1の実施の形態のモニタ5に同じである。
【0040】
PW時には、送受信系22からのドプラ信号はRG積分回路28に供給される。CW時には、送受信系22からのドプラ信号はCW帯域フィルタ29に供給される。図4(a)に示すように、RG積分回路28は、データ数Nのドプラデータを、加算器331 と、加算器331 への帰還経路に挿入された遅延回路332 とで積分する。勿論、遅延回路332 の遅延時間は、ADCのサンプリング周期△TS に設定される。データ数Nは、レンジゲートマーカの長さに応じた積分期間△TRGに依存し、△TRG/△TS によりコントローラ27で決定される。加算器331 の前段には、ドプラデータに、積分期間△TRGに応じた係数を乗算するための乗算器32が設けられる。係数は、データ数Nに基づいて、1/N1/2 でコントローラ27により与えられる。
【0041】
ドプラデータをデータ数Nで加算することにより、信号レベルはN倍になり、一方、ノイズレベルはN1/2 倍になる。したがって、レンジゲートマーカの長さ、つまり積分期間、すなわち加算データ数の変動に依存して、ノイズレベルが変動してしまう。本実施の形態では、ドプラデータに係数1/N1/2 を乗算することにより、レンジゲートマーカの長さの変動、つまり積分期間の変動、すなわち加算データ数の変動によらず、ノイズレベルを一定に保つことができる。
【0042】
図4(b)に示すように、CW帯域フィルタ29は、ドプラデータの周波数帯域を入力帯域finから出力帯域fout に制限するために、 FIR型の低域通過フィルタ(FIR LPF )34と、 IIR型の低域通過フィルタ(IIR LPF )35とが縦列に接続される。 IIR型の低域通過フィルタ35の出力側には、濾過されたドプラデータに係数を乗算するための乗算器36が設けられる。この係数は、入力帯域finと出力帯域fout との比に基づいて、(fin/fout )1/2 でコントローラ27により与えられる。
【0043】
入力帯域finが出力帯域fout に制限されることは、1/(fin/fout )1/2 <1のゲインがかかることに等価的であり、これは入力帯域finと出力帯域fout との比に応じてノイズレベルが変動してしまうことを意味する。本実施の形態では、 FIR型の低域通過フィルタ34及び IIR型の低域通過フィルタ35の出力データに係数(fin/fout )1/2 を乗算することにより、入力帯域finと出力帯域fout との比の変動に関わらず、ノイズレベルを一定に保つことができる。
【0044】
図4(c)に示すように、高速フーリエ変換処理回路31は、高速フーリエ変換処理部38でレンジゲートマーカの長さ、つまり積分期間に応じて、ドプラデータをデータ数Nで高速フーリエ変換処理に供する。高速フーリエ変換処理部38の前段には、ドプラデータに、係数を乗算するための乗算器37が設けられる。係数は、データ数Nに基づいて、1/N1/2 でコントローラ27により与えられる。
【0045】
ドプラデータをデータ数Nで高速フーリエ変換処理に供することにより、ノイズレベルに実質的にN1/2 のゲインがかかる。したがって、データ数の変動に依存して、ノイズレベルが変動してしまう。本実施の形態では、ドプラデータに係数1/N1/2 を乗算することにより、データ数の変動によらず、ノイズレベルを一定に保つことができる。
(第3の実施の形態)
図5は第3の実施の形態による超音波診断装置のブロック図である。なお、超音波プローブ41は、第1の実施の形態の超音波プローブ1に同じである。高速フーリエ変換処理回路50は、第1の実施の形態の高速フーリエ変換処理回路10に同じである。ディジタルスキャンコンバータ44は、第1の実施の形態のディジタルスキャンコンバータ4に同じである。モニタ45は、第1の実施の形態のモニタ5に同じである。
【0046】
本実施の形態による超音波診断装置は、超音波プローブ41、送受信系42、スペクトラムドプラユニット43、ディジタルスキャンコンバータ(DSC)44、モニタ45、システムコントローラ46、ドプラ用コントローラ47、コンソール48とから構成される。送受信系42は、超音波プローブ41を駆動して交互スキャンを実行することが可能に構成されている。交互スキャンとは、通常スキャンが、ある1本のラスタに対して超音波送受信動作をレート周波数で所定回数繰り返しながら、順次次のラスタに移っていくのに対して、ある1本のラスタに対して超音波送受信動作を1回行い、同じラスタに対する次回の超音波送受信動作までに間に、他の何本かのラスタに対して1回ずつ超音波送受信動作を行うものであり、1本のラスタに関して見ると、超音波送受信動作の間隔が長くなり、実質的にレート周波数が低くなって低速検出能が向上するというものである。なお、同じラスタに対する前後2回の超音波送受信動作の間に行う他のラスタの超音波送受信動作の回数を交互段数として定義する。また、送受信系42は、システムコントローラ46の制御にしたがって交互段数を変化させることが可能に構成されている。さらに、送受信系42は、ドプラ用コントローラ47の制御にしたがって送信電圧を変化させることが可能に構成されている。
【0047】
スペクトラムドプラユニット43は、ドプラゲイン補正用の乗算器49と、高速フーリエ変換処理回路50とから構成される。なお、送受信系42の出力がディジタルでなくアナログであれば、乗算器49は増幅器に置換される。
【0048】
システムコントローラ46には、コンソール48が接続される。コンソール48を介して、交互段数を変化させることが可能である。
図6に示すように、ドプラ用コントローラ47は、データ数調整用コントローラ51と、ドプラ感度補正用コントローラ52とから構成されている。データ数調整用コントローラ51は、交互段数が変化してレート周波数fPRF が実質的に変化しても、ドプラスペクトラムの表示上で周波数分解能と時間分解能との比を一定にするように、高速フーリエ変換処理回路50の周波数解析のデータ数を調整するためのものである。この調整方法は周知であるので、簡単に説明する。交互段数をK、レート周波数をfPRF 、スクロールスピードをS、1ピクセルに相当する周波数幅に対する時間幅の比を表すピクセル比をPn 、係数をαとする。まず、パラメータAA を(1)式にしたがって求める。
【0049】
AA =α・((fPRF ・S)/(Pn ・K)) …(1)
ここで、データ数をN´とする。N´=Ww ・Nで与えられる。Nは、基準データ数(通常、最小データ数とされる)であり、Ww は基準データ数に対するデータ数の比(N´/N)である。選択可能なデータ数が、16、32、64、128、256であるとすると、例えば、N=16であり、この場合、Ww =1、2、4、8又は16である。
【0050】
そして、(2)式にしたがってWw が求められる。
Ww =AA /N2 …(2)
求められたWw が、1、2、4、8、16のいずれにも一致しない場合、求められたWw に最も近いものが選択される。求められた、又は選択されたWw から、データ数N´(=Ww ・N)が決定される。
【0051】
なお、原理的には、交互段数がn倍になると、周波数解析のデータ数は、1/n1/2 倍に変更される。高速フーリエ変換処理回路50ではデータ数調整用コントローラ51から、決定されたデータ数Ww ・Nの情報を受け取り、ドプラ信号をデータ数Ww ・Nで周波数解析に供する。データ数調整用コントローラ51で決定されたデータ数Ww ・Nの情報は、ドプラ感度補正用コントローラ52にも供給される。ドプラ感度補正用コントローラ52には、システムコントローラ46から交互段数Kの情報も供給される。
【0052】
ドプラ感度補正用コントローラ52は、周波数解析のデータ数の変動に伴うS/Nの変動を抑制するためのものである。例えば、低速検出能を半値まで下げるように交互段数を変更すると、周波数解析のデータ数は、1/21/2 倍に減少する。このデータ数の減少は、S/Nを0.8dB程度低下させる。このS/Nの低下を抑制するために、ドプラゲイン補正用の乗算器49のドプラゲインを、変更前の交互段数のときより、0.8dB程度下げて、ノイズレベルを0.8dB程度下げる。これにより、データ数によらずノイズレベルを一定に保つことができるという効果が得られる。ただし、当然のこととして、ノイズレベルの低下と共に、信号レベルも0.8dB程度低下してしまう。この信号レベルの低下を補償するために、送信電圧を交互段数の変更前より1.09倍に増大する。このようなドプラゲイン調整及び送信電圧調整により、交互段数が増加して、データ数が減少しても、S/Nが劣化しないという効果が得られる。
本発明は上述した実施の形態に限定されることなく種々変形して実施可能である。
【0053】
【発明の効果】
本発明による超音波診断装置は、超音波パルスを被検体内に送信し、反射波を受信することにより受信信号を得る送受信手段と、レンジゲートマーカの深さ及び長さに応じた特定の期間に前記送受信手段から出力される前記受信信号を積分することにより積分信号を得る手段と、前記レンジゲートマーカ内のドプラスペクトラムを得るために前記積分信号を周波数解析する手段と、前記送受信手段と前記積分手段との間に設けられ、前記積分手段による積分期間をΔT RG として1/ΔT RG 1/2 に応じた重み係数を前記受信信号に乗算する手段とを具備する。
本発明によれば、積分期間が変化すると、それに伴って実質的にゲインが変化してノイズレベルも変化してしまうが、積分期間に応じた重み係数を受信信号に乗算するので、積分期間の変化に伴う実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を連続波として被検体内に送信し、反射波を受信することにより受信信号を得る手段と、前記受信信号の特定の周波数帯域を減衰するフィルタ手段と、前記フィルタ手段の出力帯域に対する前記フィルタ手段の入力帯域の比に応じた重み係数を前記フィルタ手段の出力信号に乗算する手段と、ドプラスペクトラムを得るために前記フィルタ手段の出力信号を周波数解析する手段とを具備する。
本発明によれば、例えばレート周波数が変化すると、それに伴ってフィルタ手段の出力帯域に対する入力帯域の比が変化して実質的にゲインが変化してノイズレベルも変化してしまうが、この比に応じた重み係数をフィルタ手段の出力信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
本発明による超音波診断装置は、超音波を被検体内に送信し、反射波を受信することにより受信信号を得る手段と、ドプラスペクトラムを得るために前記受信信号をデータ数で周波数解析する手段と、前記データ数をNとして重み係数1/N 1/2 を前記受信信号に乗算する手段とを具備する。
本発明によれば、周波数解析のデータ数が変化すると実質的にゲインが変化してノイズレベルも変化してしまうが、このデータ数に応じた重み係数を受信信号に乗算することにより、この実質的なゲインの変化を打ち消して、ノイズレベルを一定に保つことができる。
【図面の簡単な説明】
【図1】第1の実施の形態に係る超音波診断装置のブロック図。
【図2】第1の実施の形態による積分期間を示すタイムチャート。
【図3】第2の実施の形態に係る超音波診断装置のブロック図。
【図4】図3のRG積分回路、CW帯域フィルタ、FFTのブロック図。
【図5】第3の実施の形態に係る超音波診断装置のブロック図。
【図6】図5のドプラ用コントローラのブロック図。
【図7】第1の実施の形態に対応する従来の問題点の説明図。
【符号の説明】
1…超音波プローブ、
2…送受信系、
3…スペクトラムドプラユニット、
4…ディジタルスキャンコンバータ、
5…モニタ、
6…コンソール、
7…コントローラ、
8…レンジゲートユニット、
9…高域通過型フィルタ、
10…高速フーリエ変換処理回路、
11…乗算器、
12…積分器、
13…関数発生ROM、
14…RG幅カウンタ、
15…RGスタート位置カウンタ。
[0001]
[Industrial application fields]
The present invention relates to an ultrasonic diagnostic apparatus that obtains a Doppler spectrum by frequency analysis of a received signal.
[0002]
[Prior art]
Ultrasound has been applied in various ways from a medical standpoint, but the mainstream is an ultrasonic diagnostic apparatus that obtains a tomographic image of a soft tissue of a living body using an ultrasonic pulse reflection method. This ultrasonic diagnostic apparatus is a non-invasive examination method and displays a tomographic image of a tissue. Compared to other diagnostic apparatuses such as an X-ray diagnostic apparatus, an X-ray CT apparatus, an MRI, and a nuclear medicine diagnostic apparatus, the ultrasonic diagnostic apparatus is real-time. It has exclusive features such as display capability, small size and low cost, high safety without exposure to X-rays, and Doppler spectrum and blood flow imaging using the Doppler effect. Here, we focus on the Doppler spectrum. As a display of the Doppler spectrum, scroll display is used as in the M mode with the Doppler frequency (speed) on the vertical axis, the time and the luminance on the horizontal axis, and the strength (power) of each Doppler frequency component.
[0003]
In the Doppler spectrum function of such an ultrasonic diagnostic apparatus, burst waves have been put into practical use for the purpose of improving sensitivity. A burst wave refers to a series of a plurality of ultrasonic pulses. The time width of the burst wave is called the burst length. In the case of a burst wave, it is necessary to devise an integration process for Doppler signals in a range gate marker, unlike a pulse wave. In the previous case shown in FIG. 7A, only the Doppler signal received during the period t1 to t2 corresponding to the depth and length of the range gate marker is taken into the integration process. T1 is the time corresponding to the depth of the range gate marker from the transmission time t0 of the burst wave (the time required for the ultrasonic wave to propagate a distance twice the depth of the range gate marker in consideration of the round trip) ), And t2 is the time when the time corresponding to the length of the range gate marker (the time required for the ultrasonic wave to propagate the distance twice the length of the range gate marker) has elapsed from t1. It is. In this previous case, the burst length τ of the burst wave is not considered, that is, the reflected wave from the same reflection source is not considered to be continuously received during the burst length τ. As understood from the sensitivity distribution, a Doppler signal from a shallow region corresponding to τ from the range gate marker is taken into the integration process.
[0004]
A case devised to solve this problem is shown in FIG. In this case, the burst length τ of the burst wave is intentionally adjusted to half of the time corresponding to the length of the range gate marker. Then, only the Doppler signal received during the period t3 to t2 in the latter half of the period t1 to t2 corresponding to the depth and length of the range gate marker is taken into the integration process. As a result, the integration period and the burst length τ match. In this case, as understood from the sensitivity distribution of the integration signal, the Doppler signal outside the range gate marker is not taken into the integration process.
[0005]
However, the case of FIG. 7B has the following problems. In order to avoid a significant decrease in distance resolution, the burst length τ of the burst wave is limited. Therefore, when the range gate marker is set to be very long, there may occur a case where the integration period becomes longer than the longest burst length τmax as shown in FIG. In this case, there is a so-called insensitive area where the Doppler signal from the range gate marker is not taken into the integration process at all.
[0006]
The Doppler spectrum function also has the following problems. In the pulse wave Doppler (PW Doppler), when the length of the range gate marker is changed, the integration period changes. This change in the integration period means that the gain changes substantially. In continuous wave Doppler (CW Doppler), the Doppler signal is applied to a high-pass filter (HPF). When the rate frequency PRF is changed, the ratio of the output band to the input band of this high-pass filter changes. However, this change in the input / output bandwidth ratio means that the gain changes substantially. Furthermore, as can be said for both PW Doppler and CW Doppler, a change in the number of data of Fast Fourier Transform (FFT) as a frequency analysis means that the gain changes substantially. Therefore, there has been a problem that the noise level also changes as the gain changes.
[0007]
In PW Doppler, a function called alternate scan has been put to practical use in order to substantially reduce the rate frequency and improve the low-speed detection capability. Usually, a scan is performed in which the ultrasonic transmission / reception operation is repeated a predetermined number of times at a rate frequency with respect to a certain raster, and the next raster is sequentially moved. On the other hand, in the alternate scan, an ultrasonic transmission / reception operation is performed once for a certain raster, and 1 is performed for several other rasters before the next ultrasonic transmission / reception operation for the same raster. The ultrasonic transmission / reception operation is performed one by one. When one raster is viewed, the interval of the ultrasonic transmission / reception operation is increased, and the rate frequency is substantially reduced. The number of ultrasonic transmission / reception operations performed in the gap between the two ultrasonic transmission / reception operations before and after the same raster is defined as the number of alternating stages. An increase in the number of alternating stages means an extension of the data acquisition time if the number of frequency analysis data is constant. The extension of the data acquisition time leads to a change in the ratio between the frequency resolution and the time resolution on the Doppler spectrum display. Since this is not desirable, many current devices are equipped with a function to keep the ratio of frequency resolution to time resolution constant by decreasing the number of frequency analysis data as the number of alternating stages increases. Yes. However, a reduction in the number of data for frequency analysis degrades (decreases) the S / N.
[0008]
[Problems to be solved by the invention]
The purpose of the present invention, the change of the integration period is pulse wave Doppler, the number of data of the change in the ratio of the output band continuous wave Doppler for the input band of the high pass filter, fast Fourier transform of the frequency analysis (FFT) It is an object to provide an ultrasonic diagnostic apparatus in which the noise level does not change even if the gain changes substantially due to changes in the above .
[0009]
[Means for Solving the Problems]
An ultrasonic diagnostic apparatus according to the present invention includes a transmission / reception unit that transmits an ultrasonic pulse into a subject and receives a reflected wave to obtain a reception signal, and a specific period according to the depth and length of the range gate marker. Means for obtaining an integrated signal by integrating the received signal output from the transmitting / receiving means, means for analyzing the frequency of the integrated signal in order to obtain a Doppler spectrum in the range gate marker, the transmitting / receiving means, And a means for multiplying the received signal by a weighting factor corresponding to 1 / ΔT RG 1/2, where ΔT RG is an integration period by the integrating means.
According to the present invention, when the integration period changes, the gain changes substantially and the noise level changes accordingly, but the received signal is multiplied by the weighting coefficient corresponding to the integration period. By canceling the substantial gain change accompanying the change, the noise level can be kept constant.
The ultrasonic diagnostic apparatus according to the present invention includes means for transmitting an ultrasonic wave as a continuous wave into a subject and receiving a reflected wave to obtain a received signal; and filter means for attenuating a specific frequency band of the received signal; Means for multiplying the output signal of the filter means by a weighting factor corresponding to the ratio of the input band of the filter means to the output band of the filter means; and frequency analysis of the output signal of the filter means to obtain a Doppler spectrum. Means.
According to the present invention, for example, when the rate frequency changes, the ratio of the input band to the output band of the filter means changes accordingly, the gain changes substantially, and the noise level also changes. By multiplying the output signal of the filter means by a corresponding weighting factor, this substantial gain change can be canceled and the noise level can be kept constant.
An ultrasonic diagnostic apparatus according to the present invention includes a means for obtaining a received signal by transmitting an ultrasonic wave into a subject and receiving a reflected wave, and a means for frequency-analyzing the received signal with the number of data to obtain a Doppler spectrum. And means for multiplying the received signal by a weighting factor 1 / N 1/2 , where N is the number of data.
According to the present invention, when the number of data for frequency analysis changes, the gain substantially changes and the noise level also changes. By multiplying the received signal by a weighting factor corresponding to this number of data, this substantial The gain level can be canceled and the noise level can be kept constant.
[0019]
Embodiment
Embodiments of the present invention will be described below with reference to the drawings.
(First embodiment)
FIG. 1 is a block diagram of an ultrasonic diagnostic apparatus according to the first embodiment. The ultrasonic diagnostic apparatus according to this embodiment includes an ultrasonic probe 1, a transmission / reception system (T / R) 2, a spectrum Doppler unit 3, a digital scan converter (DSC) 4, a monitor 5, and a console 6.
[0020]
The ultrasonic probe 1 is provided with a piezoelectric element array at its tip, in which a plurality of piezoelectric elements as a reversible conversion element for mechanical motion / electrical signals such as piezoelectric ceramics are arranged one-dimensionally.
[0021]
The transmission / reception system 2 has a transmission system and a reception system. The transmission system includes a pulse generator, a transmission delay circuit, and a pulser. The pulse generator repeatedly generates rate pulses at a rate frequency fr (period: 1 / fr second) of 5 KHz, for example. This rate pulse is distributed to the number of channels and sent to the transmission delay circuit. The transmission delay circuit provides each rate pulse with a delay time necessary for focusing the ultrasonic wave into a beam and determining transmission directivity. The pulser applies a transmission voltage for each channel to the piezoelectric element of the ultrasonic probe 1 at the time when the rate pulse is received from the transmission delay circuit. Thereby, an ultrasonic beam is transmitted to the subject.
[0022]
The reflected wave reflected by the discontinuity surface (reflection source) of the acoustic impedance in the subject is received by the ultrasonic probe 1 and converted into an electric signal. The electrical signals output from the ultrasonic probe 1 for each element are collected for each channel and taken into the receiving system. The reception system includes a preamplifier, a reception delay circuit, an adder, and a phase detection circuit. The electric signal is amplified by a preamplifier, given a delay time necessary for determining reception directivity by a reception delay circuit, and added by an adder. By this addition, a reception signal in which the reflection component from the direction corresponding to the reception directivity is emphasized is obtained. This received signal is phase-detected by a phase detection circuit. Thereby, a Doppler signal having a deviation frequency due to the Doppler effect is obtained. The Doppler signal is supplied to the spectrum Doppler unit 3.
[0023]
The spectrum Doppler unit 3 takes out only the Doppler signal supplied during a specific period by the range gate marker unit 8, and uses a relatively fast moving high-frequency blood flow component and a relatively slow moving low-frequency blood vessel wall or valve. Each frequency component is separated by a high-pass filter 9 and then the Doppler signal that has passed through the high-pass filter 9 is analyzed in real time by a fast Fourier transform processing circuit (FFT) 10. Doppler spectrum data representing the strength of power is generated. The specific period during which the Doppler signal is extracted by the range gate marker unit 8 is controlled by the controller 7 based on the depth and length of the range gate marker marker set on the B-mode image via the console 6.
[0024]
The Doppler spectrum data obtained by the fast Fourier transform processing circuit 10 is displayed as a Doppler spectrum on the monitor 5 via the digital scan converter 4 and a digital analog converter (not shown). As a display of the Doppler spectrum, scroll display is used as in the M mode, with the Doppler frequency (speed) on the vertical axis, the time on the horizontal axis, and the luminance as the strength (power) of each Doppler frequency component.
[0025]
The range gate marker unit 8 amplitude-modulates the Doppler signal by a multiplier (amplifier) 11 and integrates the amplitude-modulated Doppler signal by an integrator 12 for a specific period. Amplitude modulation data by the multiplier 11 is supplied from the function generation ROM 13. The specific period during which the integrator 12 performs integration is adjusted by the RG width counter 14 and the RG start position counter 15.
[0026]
The RG start position counter 15 receives supply of RG position data corresponding to the depth of the range gate marker from the controller 7. The RG start position counter 15 starts counting clock pulses (up counting) from the start of scanning, and outputs a start pulse when the count value reaches the value represented by the RG position data.
[0027]
The RG width counter 14 receives supply of RG width data corresponding to the length of the range gate marker from the controller 7. The RG width counter 14 starts counting clock pulses (up-counting) at the time when the start pulse is received from the RG start position counter 15, and continues counting until the count value reaches the value represented by the RG width data. Therefore, the count signal is output from the RG width counter 14 only for a specific period corresponding to the depth and length of the range gate marker.
[0028]
The integrator 12 executes integration only during the period when the count signal is output from the RG width counter 14. Thereby, only the amplitude-modulated Doppler signal within a specific period is extracted as an integration signal.
[0029]
Thus, the period (integration period) of the Doppler signal extracted by the range gate unit 8 depends on the RG position data and the RG width data from the controller 7. In other words, it will be understood that the controller 7 can freely change the start time of the integration period and its time width by adjusting the RG position data and the RG width data.
[0030]
The function generation ROM 13 outputs the coefficient to the multiplier 11 only for a specific period during which the count signal is output from the RG width counter 14. The Doppler signal is amplitude-modulated by the multiplier 11 according to the coefficient from the function generation ROM 13.
[0031]
The function generation ROM 13 stores a plurality of types of weighting factors in advance. The function generation ROM 13 is supplied with a count value signal from the RG width counter 14 and a function selection signal from the controller 7 as an address signal. The coefficient of the address part corresponding to the address signal is read from the function generation ROM 13. The controller 7 selectively supplies one function selection signal corresponding to the sensitivity characteristic selected via the console 6 from the plurality of types of function selection signals to the function generation ROM 13. The function here means the time change of the coefficient read from the function generation ROM 13. Under the first function selection signal, a constant coefficient (= 1.0) is read from the function generation ROM 13, that is, it behaves as if there is no amplitude modulation processing. Under the second function selection signal, a coefficient that changes with time so as to satisfy, for example, a Hamming window function is read from the function generation ROM 13.
[0032]
Next, the operation of the present embodiment will be described. Here, the terms used below are defined. The time corresponding to the depth of the range gate marker is defined as the time required for the ultrasonic wave to propagate a distance twice the depth of the range gate marker in consideration of the round trip of ultrasonic transmission / reception. The time corresponding to the length of the range gate marker is defined as the time required for the ultrasonic wave to propagate a distance twice the length of the range gate marker.
[0033]
The controller 7 has a first mode and a second mode. The first mode is applied when the half of the time corresponding to the length of the range gate marker and the burst length τ of the burst wave match, and the second mode is the time corresponding to the length of the range gate marker. This is applied when the burst length τ (longest burst length τmax) of the burst wave is shorter than half. As described with reference to FIG. 7B, the burst length τ of the burst wave is positively adjusted to half the time corresponding to the length of the range gate marker. Depending on the range gate marker to be set, the burst length τ (maximum length of the burst wave) is longer than half of the time corresponding to the length of the range gate marker because the burst length τ is limited for the purpose of suppressing a significant drop. Occurs when the burst length τmax) is short.
[0034]
In the first mode, as described with reference to FIG. 7B, the burst length τ of the burst wave is intentionally adjusted to half the time corresponding to the length of the range gate marker. Only the Doppler signal supplied during the period of time t3 to t2 is taken into the integration process by the range gate unit 8. The time t3 is the time when the burst length τ has elapsed from the time t1, and the time t1 is the time when the time ΔTd corresponding to the depth of the range gate marker has elapsed from the transmission time t0 of the burst wave. The time t2 is the time when the time ΔTRG corresponding to the length of the range gate marker has elapsed from the time t1, that is, the time corresponding to the deepest part of the range gate marker. The integration period is set to half the time ΔTRG (ΔTRG / 2) corresponding to the length of the range gate marker. In the first mode, as understood by looking at the sensitivity distribution of the integration signal, the Doppler signal outside the range gate marker is preferably not taken into the integration process at all.
[0035]
FIGS. 2A and 2B are time charts showing the integration start time and the time width of the integration period in the second mode. FIG. 2A corresponds to the first function selection signal. FIG. 5B corresponds to the second function selection signal.
[0036]
In the second mode, as shown in FIGS. 2A and 2B, only the Doppler signal supplied during the integration period from time t1 'to time t2' is taken into the integration process by the range gate unit 8. The time t1 'is the time when half the burst length τ (τ / 2) has elapsed from the time t1, and the time t1 has passed the time ΔTd corresponding to the depth of the range gate marker from the burst wave transmission time t0. It's time. Time t2 'is the time when time ΔTRG corresponding to the length of the range gate marker has elapsed from time t1'. The integration period is set to a time ΔTRG corresponding to the length of the range gate marker. In the second mode, there is no dead area in the range gate marker, as can be seen by looking at the sensitivity distribution of the integrated signal.
[0037]
As described above, in the second mode, integration is started at time t1 ′ delayed by τ / 2 from time t1 corresponding to the depth of the range gate marker, and time ΔTRG corresponding to the length of the range gate marker, By continuing the integration, Doppler signals from a range corresponding to τ / 2 are taken into the integration process on the shallower side and the deeper side than the range gate marker, respectively, but compared with the case of FIG. The amount (area) of Doppler signals taken outside the range gate marker is absolutely reduced. Further, as understood from comparison with FIG. 7C, even in the case of facing a case where the integration period becomes longer than the burst length τ, the dead area in the range gate marker can be completely eliminated.
[0038]
Further, as understood from the sensitivity distribution of the integrated signal in FIG. 2B, a coefficient that changes over time so as to satisfy, for example, a Hamming window function is read from the function generation ROM 13 and multiplied by the Doppler signal. The sensitivity distribution becomes higher near the center, and the sensitivity distribution in the range corresponding to τ / 2 is suppressed on each of the shallower side and the deeper side than the range gate marker.
(Second Embodiment)
FIG. 3 is a block diagram of an ultrasonic diagnostic apparatus according to the second embodiment. The ultrasonic diagnostic apparatus according to the present embodiment includes an ultrasonic probe 21, a transmission / reception system (T / R) 22, a spectrum Doppler unit 23, a digital scan converter (DSC) 24, a monitor 25, and a console 26.
[0039]
The ultrasonic probe 21 is the same as the ultrasonic probe 1 of the first embodiment. The transmission / reception system 22 is improved so that the transmission / reception system 2 of the first embodiment can selectively generate an ultrasonic wave as a pulse wave (PW) or a continuous wave (CW), and an analog Doppler signal is generated. An ADC is added to the last stage so as to digitize and output the data. The wide-pass filter 30 is the same as the wide-pass filter 9 of the first embodiment. The fast Fourier transform processing circuit 31 is the same as the fast Fourier transform processing circuit 10 of the first embodiment. The digital scan converter 24 is the same as the digital scan converter 4 of the first embodiment. The monitor 25 is the same as the monitor 5 of the first embodiment.
[0040]
During PW, the Doppler signal from the transmission / reception system 22 is supplied to the RG integration circuit 28. During CW, the Doppler signal from the transmission / reception system 22 is supplied to the CW band filter 29. As shown in FIG. 4A, the RG integration circuit 28 integrates the Doppler data of the number of data N by the adder 331 and the delay circuit 332 inserted in the feedback path to the adder 331. Of course, the delay time of the delay circuit 332 is set to the ADC sampling period ΔTS. The number N of data depends on the integration period ΔTRG corresponding to the length of the range gate marker and is determined by the controller 27 based on ΔTRG / ΔTS. A multiplier 32 for multiplying the Doppler data by a coefficient corresponding to the integration period ΔTRG is provided before the adder 331. The coefficient is given by the controller 27 at 1 / N 1/2 based on the number of data N.
[0041]
By adding Doppler data by the number N of data, the signal level is increased N times, while the noise level is increased N 1/2 times. Therefore, the noise level varies depending on the length of the range gate marker, that is, the integration period, that is, the variation in the number of added data. In the present embodiment, by multiplying the Doppler data by a coefficient 1 / N 1/2 , the noise level is reduced regardless of the variation in the length of the range gate marker, that is, the variation in the integration period, that is, the variation in the number of added data. Can be kept constant.
[0042]
As shown in FIG. 4B, the CW band filter 29 includes an FIR type low-pass filter (FIR LPF) 34, IIR to limit the frequency band of Doppler data from the input band fin to the output band fout. A series of low pass filters (IIR LPF) 35 are connected in series. A multiplier 36 for multiplying the filtered Doppler data by a coefficient is provided on the output side of the IIR type low-pass filter 35. This coefficient is given by the controller 27 at (fin / fout) 1/2 based on the ratio of the input band fin and the output band fout.
[0043]
Limiting the input band fin to the output band fout is equivalent to applying a gain of 1 / (fin / fout) 1/2 <1, which is a ratio between the input band fin and the output band fout. This means that the noise level fluctuates accordingly. In this embodiment, by multiplying the output data of the FIR type low pass filter 34 and the IIR type low pass filter 35 by a coefficient (fin / fout) 1/2 , the input band fin and the output band fout and The noise level can be kept constant regardless of the fluctuation of the ratio.
[0044]
As shown in FIG. 4C, the fast Fourier transform processing circuit 31 performs fast Fourier transform processing on the Doppler data with the number of data N in the fast Fourier transform processing unit 38 according to the length of the range gate marker, that is, the integration period. To serve. A multiplier 37 for multiplying the Doppler data by a coefficient is provided in the preceding stage of the fast Fourier transform processing unit 38. The coefficient is given by the controller 27 at 1 / N 1/2 based on the number of data N.
[0045]
By applying the Doppler data to the fast Fourier transform process with the number of data N, a gain of N1 / 2 is applied to the noise level substantially. Therefore, the noise level varies depending on the variation in the number of data. In the present embodiment, by multiplying the Doppler data by the coefficient 1 / N 1/2 , the noise level can be kept constant regardless of the variation in the number of data.
(Third embodiment)
FIG. 5 is a block diagram of an ultrasonic diagnostic apparatus according to the third embodiment. The ultrasonic probe 41 is the same as the ultrasonic probe 1 of the first embodiment. The fast Fourier transform processing circuit 50 is the same as the fast Fourier transform processing circuit 10 of the first embodiment. The digital scan converter 44 is the same as the digital scan converter 4 of the first embodiment. The monitor 45 is the same as the monitor 5 of the first embodiment.
[0046]
The ultrasonic diagnostic apparatus according to the present embodiment includes an ultrasonic probe 41, a transmission / reception system 42, a spectrum Doppler unit 43, a digital scan converter (DSC) 44, a monitor 45, a system controller 46, a Doppler controller 47, and a console 48. Is done. The transmission / reception system 42 is configured to be able to drive the ultrasonic probe 41 and execute alternate scanning. The alternate scan is a normal scan in which an ultrasonic transmission / reception operation is repeated a predetermined number of times at a rate frequency for a certain raster, and then moves to the next raster sequentially. The ultrasonic transmission / reception operation is performed once, and the ultrasonic transmission / reception operation is performed once for several other rasters before the next ultrasonic transmission / reception operation for the same raster. Looking at the raster, the interval between ultrasonic transmission / reception operations is increased, the rate frequency is substantially lowered, and the low-speed detection capability is improved. The number of ultrasonic transmission / reception operations of other rasters performed between two ultrasonic transmission / reception operations before and after the same raster is defined as the number of alternating stages. Further, the transmission / reception system 42 is configured to be able to change the number of alternating stages according to the control of the system controller 46. Further, the transmission / reception system 42 is configured to be able to change the transmission voltage in accordance with the control of the Doppler controller 47.
[0047]
The spectrum Doppler unit 43 includes a multiplier 49 for Doppler gain correction and a fast Fourier transform processing circuit 50. If the output of the transmission / reception system 42 is not digital but analog, the multiplier 49 is replaced with an amplifier.
[0048]
A console 48 is connected to the system controller 46. The number of alternating stages can be changed via the console 48.
As shown in FIG. 6, the Doppler controller 47 includes a data number adjustment controller 51 and a Doppler sensitivity correction controller 52. The controller 51 for adjusting the number of data performs fast Fourier transform so that the ratio of the frequency resolution to the time resolution is constant on the display of the Doppler spectrum even if the rate frequency fPRF is substantially changed by changing the number of alternate stages. This is for adjusting the number of frequency analysis data of the processing circuit 50. This adjustment method is well known and will be described briefly. The number of alternating stages is K, the rate frequency is fPRF, the scroll speed is S, the pixel ratio representing the ratio of the time width to the frequency width corresponding to one pixel is Pn, and the coefficient is α. First, the parameter AA is obtained according to the equation (1).
[0049]
AA = α. ((FPRF.S) / (Pn.K)) (1)
Here, the number of data is N ′. N ′ = Ww · N is given. N is the reference data number (usually the minimum data number), and Ww is the ratio of the data number to the reference data number (N ′ / N). If the number of selectable data is 16, 32, 64, 128, 256, for example, N = 16, and in this case, Ww = 1, 2, 4, 8, or 16.
[0050]
Then, Ww is obtained according to equation (2).
Ww = AA / N 2 (2)
If the obtained Ww does not match any of 1, 2, 4, 8, and 16, the one closest to the obtained Ww is selected. The number of data N ′ (= Ww · N) is determined from the obtained or selected Ww.
[0051]
In principle, when the number of alternating stages is n times, the number of data for frequency analysis is changed to 1 / n 1/2 times. The fast Fourier transform processing circuit 50 receives information on the determined data number Ww · N from the data number adjustment controller 51, and uses the Doppler signal for frequency analysis with the data number Ww · N. The information on the data number Ww · N determined by the data number adjustment controller 51 is also supplied to the Doppler sensitivity correction controller 52. The Doppler sensitivity correction controller 52 is also supplied with information on the number of alternating stages K from the system controller 46.
[0052]
The Doppler sensitivity correction controller 52 is for suppressing S / N fluctuations due to fluctuations in the number of data of frequency analysis. For example, when the number of alternating stages is changed so as to reduce the low speed detection capability to half value, the number of data of frequency analysis is reduced to 1/2 1/2 times. This reduction in the number of data reduces the S / N by about 0.8 dB. In order to suppress this decrease in S / N, the Doppler gain of the multiplier 49 for Doppler gain correction is lowered by about 0.8 dB compared to the number of alternating stages before the change, and the noise level is lowered by about 0.8 dB. Thereby, the effect that a noise level can be kept constant irrespective of the number of data is acquired. However, as a matter of course, the signal level decreases by about 0.8 dB as the noise level decreases. In order to compensate for this decrease in signal level, the transmission voltage is increased 1.09 times before the change in the number of alternating stages. By such Doppler gain adjustment and transmission voltage adjustment, even if the number of alternating stages increases and the number of data decreases, the effect that S / N does not deteriorate can be obtained.
The present invention is not limited to the above-described embodiment, and can be implemented with various modifications.
[0053]
【The invention's effect】
An ultrasonic diagnostic apparatus according to the present invention includes a transmission / reception unit that transmits an ultrasonic pulse into a subject and receives a reflected wave to obtain a reception signal, and a specific period according to the depth and length of the range gate marker. Means for obtaining an integrated signal by integrating the received signal output from the transmitting / receiving means, means for analyzing the frequency of the integrated signal in order to obtain a Doppler spectrum in the range gate marker, the transmitting / receiving means, And a means for multiplying the received signal by a weighting factor corresponding to 1 / ΔT RG 1/2, where ΔT RG is an integration period by the integrating means.
According to the present invention, when the integration period changes, the gain changes substantially and the noise level changes accordingly, but the received signal is multiplied by the weighting coefficient corresponding to the integration period. By canceling the substantial gain change accompanying the change, the noise level can be kept constant.
The ultrasonic diagnostic apparatus according to the present invention includes means for transmitting an ultrasonic wave as a continuous wave into a subject and receiving a reflected wave to obtain a received signal; and filter means for attenuating a specific frequency band of the received signal; Means for multiplying the output signal of the filter means by a weighting factor corresponding to the ratio of the input band of the filter means to the output band of the filter means; and frequency analysis of the output signal of the filter means to obtain a Doppler spectrum. Means.
According to the present invention, for example, when the rate frequency changes, the ratio of the input band to the output band of the filter means changes accordingly, the gain changes substantially, and the noise level also changes. By multiplying the output signal of the filter means by a corresponding weighting factor, this substantial gain change can be canceled and the noise level can be kept constant.
An ultrasonic diagnostic apparatus according to the present invention includes a means for obtaining a received signal by transmitting an ultrasonic wave into a subject and receiving a reflected wave, and a means for frequency-analyzing the received signal with the number of data to obtain a Doppler spectrum. And means for multiplying the received signal by a weighting factor 1 / N 1/2 , where N is the number of data.
According to the present invention, when the number of data for frequency analysis changes, the gain substantially changes and the noise level also changes. By multiplying the received signal by a weighting factor corresponding to this number of data, this substantial The gain level can be canceled and the noise level can be kept constant.
[Brief description of the drawings]
FIG. 1 is a block diagram of an ultrasonic diagnostic apparatus according to a first embodiment.
FIG. 2 is a time chart showing an integration period according to the first embodiment.
FIG. 3 is a block diagram of an ultrasonic diagnostic apparatus according to a second embodiment.
4 is a block diagram of the RG integration circuit, CW band filter, and FFT in FIG. 3;
FIG. 5 is a block diagram of an ultrasonic diagnostic apparatus according to a third embodiment.
6 is a block diagram of the Doppler controller of FIG. 5. FIG.
FIG. 7 is an explanatory diagram of a conventional problem corresponding to the first embodiment.
[Explanation of symbols]
1 ... ultrasonic probe,
2 ... transmission and reception system,
3 ... Spectrum Doppler unit,
4 ... Digital scan converter,
5 ... Monitor,
6 ... Console,
7 ... Controller,
8 ... Range gate unit,
9 ... High-pass filter,
10: Fast Fourier transform processing circuit,
11 ... multiplier,
12 ... integrator,
13 ... Function generation ROM,
14 ... RG width counter,
15: RG start position counter.

Claims (3)

超音波パルスを被検体内に送信し、反射波を受信することにより受信信号を得る送受信手段と、
レンジゲートマーカの深さ及び長さに応じた特定の期間に前記送受信手段から出力される前記受信信号を積分することにより積分信号を得る手段と、
前記レンジゲートマーカ内のドプラスペクトラムを得るために前記積分信号を周波数解析する手段と、
前記送受信手段と前記積分手段との間に設けられ、前記積分手段による積分期間をΔT RG として1/ΔT RG 1/2 に応じた重み係数を前記受信信号に乗算する手段とを具備したことを特徴とする超音波診断装置。
Transmitting / receiving means for transmitting an ultrasonic pulse into the subject and receiving a reflected wave to obtain a received signal;
Means for obtaining an integrated signal by integrating the received signal output from the transmitting / receiving means in a specific period according to the depth and length of the range gate marker;
Means for frequency analysis of the integrated signal to obtain a Doppler spectrum in the range gate marker;
Provided between the transmitting / receiving means and the integrating means, and means for multiplying the received signal by a weighting coefficient corresponding to 1 / ΔT RG 1/2, where ΔT RG is the integration period of the integrating means. A characteristic ultrasonic diagnostic apparatus.
超音波を連続波として被検体内に送信し、反射波を受信することにより受信信号を得る手段と、
前記受信信号の特定の周波数帯域を減衰するフィルタ手段と、
前記フィルタ手段の出力帯域に対する前記フィルタ手段の入力帯域の比に応じた重み係数を前記フィルタ手段の出力信号に乗算する手段と、
ドプラスペクトラムを得るために前記フィルタ手段の出力信号を周波数解析する手段とを具備したことを特徴とする超音波診断装置。
Means for transmitting an ultrasonic wave as a continuous wave into a subject and receiving a reflected wave to obtain a received signal;
Filter means for attenuating a specific frequency band of the received signal;
Means for multiplying the output signal of the filter means by a weighting factor according to the ratio of the input band of the filter means to the output band of the filter means;
An ultrasonic diagnostic apparatus comprising: means for frequency analysis of an output signal of the filter means in order to obtain a Doppler spectrum.
超音波を被検体内に送信し、反射波を受信することにより受信信号を得る手段と、
ドプラスペクトラムを得るために前記受信信号をデータ数で周波数解析する手段と、
前記データ数をNとして重み係数1/N 1/2 を前記受信信号に乗算する手段とを具備したことを特徴とする超音波診断装置。
Means for obtaining a received signal by transmitting an ultrasonic wave into a subject and receiving a reflected wave;
Means for frequency analysis of the received signal by the number of data to obtain a Doppler spectrum;
An ultrasonic diagnostic apparatus comprising: means for multiplying the received signal by a weighting factor 1 / N 1/2 , where N is the number of data.
JP20219595A 1995-08-08 1995-08-08 Ultrasonic diagnostic equipment Expired - Lifetime JP3691874B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP20219595A JP3691874B2 (en) 1995-08-08 1995-08-08 Ultrasonic diagnostic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP20219595A JP3691874B2 (en) 1995-08-08 1995-08-08 Ultrasonic diagnostic equipment

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2004203963A Division JP3663206B2 (en) 2004-07-09 2004-07-09 Ultrasonic diagnostic equipment

Publications (2)

Publication Number Publication Date
JPH0947451A JPH0947451A (en) 1997-02-18
JP3691874B2 true JP3691874B2 (en) 2005-09-07

Family

ID=16453549

Family Applications (1)

Application Number Title Priority Date Filing Date
JP20219595A Expired - Lifetime JP3691874B2 (en) 1995-08-08 1995-08-08 Ultrasonic diagnostic equipment

Country Status (1)

Country Link
JP (1) JP3691874B2 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4711498B2 (en) * 2000-10-12 2011-06-29 株式会社日立メディコ Ultrasound diagnostic imaging equipment
US7645236B2 (en) * 2005-06-28 2010-01-12 Siemens Medical Solutions Usa, Inc. Ultrasound imaging system having motion adaptive gain
JP6113892B2 (en) * 2016-05-30 2017-04-12 東芝メディカルシステムズ株式会社 Ultrasonic diagnostic apparatus, image processing apparatus, and image processing program

Also Published As

Publication number Publication date
JPH0947451A (en) 1997-02-18

Similar Documents

Publication Publication Date Title
US6937883B2 (en) System and method for generating gating signals for a magnetic resonance imaging system
JP2005536309A (en) System and method for improved harmonic imaging
JP6218400B2 (en) Ultrasonic diagnostic apparatus and control program for ultrasonic diagnostic apparatus
CN110301936B (en) Frequency scanning for acoustic radiation force pulses
JPH0975348A (en) Ultrasonic diagnostic equipment
WO2014013839A1 (en) Ultrasonic diagnostic device and image processing device
US11039777B2 (en) Ultrasonic diagnostic apparatus and control method
US6599248B1 (en) Method and apparatus for ultrasound diagnostic imaging
JP4334032B2 (en) Ultrasonic diagnostic equipment
JPH0542138A (en) Ultrasonic diagnostic device and its calibrating method
US6039692A (en) Method and apparatus for processing ultrasound signals
JP3691874B2 (en) Ultrasonic diagnostic equipment
JP5247958B2 (en) Ultrasonic diagnostic apparatus and ultrasonic echo signal processing method
JP3943653B2 (en) Ultrasonic diagnostic equipment
JP3663206B2 (en) Ultrasonic diagnostic equipment
JP4427139B2 (en) Ultrasonic diagnostic apparatus and ultrasonic transmission method
JP2011136224A (en) Ultrasonograph and ultrasonic-echo-signal processing method
US11129598B2 (en) Calibration for ARFI imaging
CN114727807A (en) Reduction of reverberation artifacts in ultrasound images and associated devices, systems, and methods
JPH03261466A (en) Ultrasonic diagnostic device
KR101809358B1 (en) Ultrasonic doppler imaging apparatus using new plane wave synthesis method and controlling method thereof
JP4077913B2 (en) Ultrasonic diagnostic equipment
US11896428B2 (en) Adaptive selection of ultrasound frequency
Maeda et al. Sensitivity time control for chirp transmission
US20240138681A1 (en) Acoustic imaging and measurements using windowed nonlinear frequency modulation chirp

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20040518

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20040525

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040709

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20050322

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050518

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: 20050614

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050617

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

Free format text: PAYMENT UNTIL: 20090624

Year of fee payment: 4

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

Free format text: PAYMENT UNTIL: 20100624

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20100624

Year of fee payment: 5

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

Free format text: PAYMENT UNTIL: 20110624

Year of fee payment: 6

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

Free format text: PAYMENT UNTIL: 20120624

Year of fee payment: 7

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

Free format text: PAYMENT UNTIL: 20130624

Year of fee payment: 8

EXPY Cancellation because of completion of term