JP3691874B2 - Ultrasonic diagnostic equipment - Google Patents
Ultrasonic diagnostic equipment Download PDFInfo
- 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
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
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
[0020]
The
[0021]
The transmission /
[0022]
The reflected wave reflected by the discontinuity surface (reflection source) of the acoustic impedance in the subject is received by the
[0023]
The
[0024]
The Doppler spectrum data obtained by the fast Fourier
[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
[0026]
The RG start position counter 15 receives supply of RG position data corresponding to the depth of the range gate marker from the
[0027]
The RG width counter 14 receives supply of RG width data corresponding to the length of the range gate marker from the
[0028]
The
[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
[0030]
The
[0031]
The
[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
[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
(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
[0039]
The
[0040]
During PW, the Doppler signal from the transmission /
[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
[0042]
As shown in FIG. 4B, the
[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
[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
[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
(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
[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
[0047]
The spectrum Doppler unit 43 includes a
[0048]
A console 48 is connected to the
As shown in FIG. 6, the
[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
[0052]
The Doppler
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
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.
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)
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 |
-
1995
- 1995-08-08 JP JP20219595A patent/JP3691874B2/en not_active Expired - Lifetime
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 |