JP3544722B2 - 超音波診断装置 - Google Patents
超音波診断装置 Download PDFInfo
- Publication number
- JP3544722B2 JP3544722B2 JP31310294A JP31310294A JP3544722B2 JP 3544722 B2 JP3544722 B2 JP 3544722B2 JP 31310294 A JP31310294 A JP 31310294A JP 31310294 A JP31310294 A JP 31310294A JP 3544722 B2 JP3544722 B2 JP 3544722B2
- Authority
- JP
- Japan
- Prior art keywords
- frame
- feature
- feature amount
- dimensional
- image
- 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 - Fee Related
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Description
【産業上の利用分野】
本発明は、被検体の2次元領域を超音波で走査し、得られた受信信号を検波して組織断層像を生成する超音波診断装置に係り、特に、正常心筋と差別して梗塞心筋を画像化する新規な技術に関する。
【0002】
【従来の技術】
超音波は、音響インピーダンスが異なる2つの媒質の境界で反射、屈折を繰り返しながら伝搬する。反射強度は2つの媒質の音響インピーダンスの差に依存する。超音波イメージングは反射信号の振幅を輝度変換して表示するもので、反射信号の振幅が小さい、つまり音響インピーダンスの差が小さいと2つの媒質間で輝度変化(濃度変化)が小さくなり、したがってこのような音響インピーダンスの差が小さい2つの媒質を超音波断層像上で識別する精度は、診断医師の知識や経験に依存する場合が多く、客観性に乏しい。
【0003】
このようなケースは例えば心筋梗塞の超音波診断に当て嵌まる。つまり、梗塞心筋と正常心筋との音響インピーダンスの差は小さく、両者を超音波断層像上で識別することは困難であり且つ客観性に乏しい。
【0004】
【発明が解決しようとする課題】
本発明では、正常心筋と差別して梗塞心筋を特徴的に画像化することのできる超音波診断装置を提供しようとするものである。
【0005】
【課題を解決するための手段】
本発明の第1局面は、被検体の2次元又は3次元領域を超音波で走査することにより、1フレーム分の受信信号を繰り返し得、この受信信号を検波することにより画像を生成する超音波診断装置において、第1フレームの第1局所領域内の受信信号と、前記第1フレームと走査タイミングの異なる第2フレーム内における前記第1局所領域に対して空間的に近傍する複数の第2局所領域内の受信信号各々との間で相互相関係数を計算し、前記計算した相互相関係数の中の最大値を用いて前記第1局所領域に関する前記第2局所領域に対する類似性を表す特徴量を定義し、前記第1局所領域の位置を変えながら前記特徴量の2次元分布を生成する特徴量2次元分布生成手段と、前記2次元分布を特徴量画像として表示する表示手段とを具備する。
本発明の第2局面は、被検体の2次元又は3次元領域を超音波で走査することにより、1フレーム分の受信信号を繰り返し得、この受信信号を検波することにより画像を生成する超音波診断装置において、第1フレームの第1局所領域内の受信信号と、前記第1フレームと走査タイミングの異なる第2フレーム内における前記第1局所領域に対して空間的に近傍する複数の第2局所領域内の受信信号各々との間での信号差分の絶対値の総和を計算し、前記計算した総和の中の最小値を用いて前記第1局所領域に関する前記第2局所領域に対する類似性を表す特徴量を定義し、前記第1局所領域の位置を変えながら前記特徴量の2次元分布を生成する特徴量2次元分布生成手段と、前記2次元分布を特徴量画像として表示する表示手段とを具備することを特徴とする。
【0006】
【作用】
本発明によれば、正常心筋の周期的な収縮に従って、組織が変型しないで単に位置が変化(転移)する梗塞部分を、第1フレームの第1局所領域内の受信信号と第2フレームの第2局所領域内の受信信号との間の類似性を高いとして、とらえることができる。
【0007】
【実施例】
以下、図面を参照して本発明の一実施例を説明する。なお、ここでは心臓診断に好適な浅部視野が狭く深部視野の広いセクタ電子走査法が採用されるものとして説明するが、被検体の2次元領域を走査可能であれば、リニア電子走査法やコンベックス電子走査法さらに他の走査法であってもかまわない。
音響信号と電気信号とを可逆的に変換する複数の振動子を配列したセクタ走査用の振動子アレイを先端に装備してなる超音波プローブ1は、送信時には送信系2に接続され、受信時には受信系3に接続される。送信系2は、クロック発生器、レートパルス発生器、振動子毎(チャンネル毎)に設けられた遅延回路、振動子毎(チャンネル毎)に設けられたパルサから構成される。クロック発生器から発生するクロックパルスはレートパルス発生器で例えば5KHzのレートパルスに分周される。このレートパルスは分配され各遅延回路を介してパルサに送られる。各遅延回路は、超音波をビーム状に集束し且つこの超音波ビームの方向(走査線方向または方位方向と呼ばれる)を決定するのに必要な遅延時間をレートパルスに与える。パルサはレートパルスを受けたタイミングで、対応する振動子にパルス電圧を印加する。これによりプローブ1から超音波パルスが遅延時間に応じた方位方向に発射される。被検体内の音響インピーダンスの境界で反射した反射波は、プローブ1の各振動子で電気信号に変換され、受信器3で取り込まれる。
【0008】
受信系3としては、特公平6−14934号に示されているようなディジタルビームフォーマ方式が好例として挙げられ、前置増幅器、アナログディジタル変換器、ディジタル遅延回路、加算器から構成される。各振動子からのアナログの受信信号は、前置増幅器で増幅された後、アナログディジタル変換器でディジタル信号に変換され、ディジタル遅延回路に送られ、そこで受信指向性を得るために、送信時とは逆の遅延時間をチャンネル毎に与えられ、さらに加算器で加算される。この加算ディジタル信号を以下、RF受信信号と称する。上述した送受信時の遅延時間制御を超音波パルスの繰り返し毎に順次変化させることにより、被検体の2次元領域をセクタ走査して、1フレーム分のRF受信信号群を得ることができる。さらに、このようなセクタ走査を繰り返すことにより、RF受信信号群を繰り返し得ることができる。
【0009】
受信系3からの出力は、信号処理部4で検波され、さらにダイナミックレンジを実質的に拡大する対数増幅にかけられる。これによりBモード画像(組織断層像)データが生成される。このBモード画像データはディジタルスキャンコンバータ(DSC)5を介して表示器6に送られ、アナログに変換された後、ビジュアルなBモード画像として濃淡表示される。
【0010】
また、受信系3からの出力は、特徴量画像生成部7に送られる。特徴量画像生成部7は、走査タイミングが異なり、且つ心拍時相の異なる2フレームのRF受信信号群を用いて、心筋の梗塞程度に応じて値が変化する特徴量の2次元分布(特徴量画像)を生成することができるように構成される。この特徴量画像データは表示器6に送られ、アナログに変換された後、ビジュアルな特徴量画像として濃淡またはカラー表示される。
【0011】
図2は特徴量画像生成部7のブロック図である。フレームメモリ11は、一連のセクタ走査の繰り返しにより得られた時間的に連続する複数フレーム分のRF受信信号群を記憶できる容量を有している。また、フレームメモリ11には2つの出力ポートが設けられている。一方の出力ポートには、第1フレームのRF受信信号群が装填され、他方の出力ポートには、第1フレームと心拍時相の異なる、例えば第1フレームから1フレーム分遅れた第2フレームのRF受信信号群が装填される。
【0012】
例えば複数フレーム分のRF受信信号群はフレーム単位でフレーム番号nで管理される。フレーム番号nは、走査順序にしたがって連番として付される。同様に、DSC5では画像はフレーム番号Nで管理される。フレーム番号nとフレーム番号Nとは対応させておく。表示器6に表示された特定の画像を図示しない入力器を介してオペレータが指定すると、当該特定の画像のフレーム番号Nに対応する第1フレームとしてのフレーム番号nのRF受信信号群が一方の出力ポートに装填され、第2フレームとしてのフレーム番号n+1のRF受信信号群が他方の出力ポートに装填される。このようにオペレータは、特徴量画像を作成するための心拍時相を、当該所望とする心拍時相の画像を指定することで特定することが可能となる。
【0013】
一方の出力ポートからは所定のサイズの局所領域(カーネルROIという)内のRF受信信号群が、また他方の出力ポートからは同じサイズであって、カーネルROIの近傍の局所領域(オブジェクトROIという)内のRF受信信号群がそれぞれ特徴量演算部12に取り込まれる。特徴量演算部12は、両RF受信信号群どうしの類似性を表す特徴量、例えば相互相関係数を計算する。特徴量演算部12は、カーネルROIの位置を固定したままで、オブジェクトROIの位置を所定の範囲のサーチROI内で転移しながら、特徴量の計算を繰り返す。これにより、或る位置のカーネルROIに対して、複数の特徴量が計算され、特徴量演算部12はこれら複数の特徴量の中から最大値(または最小値)を選択し、この選択した最大値を当該カーネルROIに対応する特徴量として認識する。特徴量演算部12は、カーネルROIの位置を変えながら、このような特徴量の計算を繰り返す。これにより、特徴量の2次元分布、つまり特徴量画像データが得られる。この特徴量画像データはディジタルスキャンコンバータ(DSC)13を介して表示器6に送られ、アナログに変換された後、ビジュアルな特徴量画像として濃淡またはカラー表示される。
【0014】
なお、制御系10は、1フレーム分の特徴量画像データの生成処理及びそのDSC13への格納が終了して特徴量画像の表示が可能となったとき、特徴量画像表示が可能となった旨のメッセージをDSC13を介して表示器6に表示させるようにすることが好ましい。このメッセージが表示された後、図示しない入力器を介してオペレータが特徴量画像を表示させるための指示をしたとき、制御系10は、これを受けてDSC13から特徴量画像データを表示器6に読み出させ表示させる。
【0015】
図3はDSC13のブロック図である。特徴量演算部12で求められた特徴量画像データは、輝度変換部14、座標変換部15、フレームメモリ16を介して適当な例えばオペレータに指定されたタイミングでディジタルビデオ信号として表示器6に送られる。輝度変換部14は、特徴量を輝度データ(またはRGB色相データ)に変換するものであり、例えば、入力を特徴量とし、出力を輝度としたROMで構成される。この特徴量と輝度の対応関係(いわゆる設定)は割り付けテーブルに纏められて輝度変換部14に記憶される。輝度変換部14は、図5に示すような設定の異なる複数種類のテーブルを装備し、オペレータの指定により任意のテーブルを選択可能に構成されていることが好ましい。図5(a)は輝度変換の設定を示し、図5(b)は色相変換の設定を示している。座標変換部15は、特徴量が表現されている座標系(θ,d)を、表示器6に応じた適当な座標系、例えば直交座標系(x,y)に変換する。θとは超音波ビームの振り角,dは深さを示している。なお、DSC13では特徴量の空間的な補間を実行することが好ましい。
【0016】
図4は特徴量画像の生成過程の概略的説明図である。上述したように特徴量とは、心筋の梗塞程度を数値化した梗塞心筋を正常心筋と差別化可能な計算値である。例えば、特徴量が大きいほど梗塞程度が強く、特徴量が小さいほど正常に近くなることを意味する。ところで、梗塞心筋は、正常心筋に対して収縮率が小さく、心臓の拍動に関係なく、その部分の組織が固形化(梗塞)している。したがって、梗塞部分は、正常心筋の周期的な収縮に従って、単に位置が変化(転移)するのみである。本発明はこの点に着目したものであり、つまり、第1フレームの梗塞心筋部分と、その部分が転移した第2フレームの梗塞心筋部分とは、組織変化が少ないので、その部分内のRF受信信号群の分布は両フレーム間で非常に類似(酷似)し、このRF受信信号群どうしの類似性を数値化したものが特徴量として与えられる。しかし、梗塞心筋部分はどの位置に転移するか不明であるのが問題である。この問題は後述するように、カーネルROIの位置を固定したままで、オブジェクトROIの位置をサーチROI内で転移しながら、その各位置で特徴量を計算し、複数の特徴量の中の最大値を当該カーネルROIに対応する特徴量として決定する、つまり梗塞心筋部分が転移する可能性のある広い範囲でオブジェクトROIを移動して特徴量をサーチすることで解決する。以下に特徴量画像の生成過程を具体的に説明する。
【0017】
図4(a)に示すように、走査タイミングが異なり、且つ心拍時相の異なる2フレームのRF受信信号群、好ましくはBモード像を観察しながらオペレータにより指定されたフレーム番号nの第1フレームのRF受信信号群と、第1フレームから1フレーム分遅れたフレーム番号n+1の第2フレームのRF受信信号群とから、特徴量画像が生成される。
【0018】
まず、図4(b),(c),(d)を参照して、或るカーネルROI(中心点(xm,ym,n ))の特徴量を求める手順を説明する。フレームメモリ11の一方の出力ポートからは第1フレームのカーネルROI内のRF受信信号群が、また他方の出力ポートからは、カーネルROIの近傍の第2フレームのオブジェクトROI内のRF受信信号群が、特徴量演算部12に取り込まれる。両RF受信信号群の分布の類似性(特徴量)が、例えば相互相関係数として計算される。特徴量の計算方法の詳細は後述する。このカーネルROIの位置を固定したままで、オブジェクトROIの位置をサーチROI内で転移し、この特徴量を計算する。さらに、オブジェクトROIの位置を転移し、同様に特徴量を計算する。このようにオブジェクトROIをサーチROI内で転移しながら、カーネルROIの位置が同じで、オブジェクトROIの位置が異なる複数の特徴量を求める。これら複数の特徴量の中の最大値(または最小値)が、当該カーネルROIに対応する特徴量として認識され、この最大値に当該カーネルROIの中心点の位置情報が与えられる。
【0019】
そして、カーネルROIを例えば中心点(xm+1,ym+1,n )に移動し、同様に特徴量を求める。カーネルROIの位置を変えながら、その各位置で特徴量を求める。これにより、特徴量の2次元分布、つまり特徴量画像が求められる。
【0020】
こうして求められた特徴量画像は、DSC13の輝度変換部14、座標変換部15、フレームメモリ16を順に介して、適当な例えばオペレータに指定されたタイミングでディジタルビデオ信号として表示器6に送られ、図6(a)に示すようにBモード画像と並列に、または図6(b)に示すようにBモード画像に重ね合わして表示される。
【0021】
次に特徴量の計算方法を説明する。特徴量の計算方式としては、2つの方式を提供する。特徴量演算部12を両方式実行可能に構成して、オペレータの指示で選択的に使用するようにしてもよいし、予め一方の方式を固定的に採用してもよい。第1の方式は特徴量として相互相関係数を計算する方式であり、第2の方式は特徴量として両フレームのカーネルとROIオブジェクトROI間での差分の絶対値の総和を計算する方式である。
【0022】
さらに、第1の方式として、直接法と2次元フーリエ変換法とが提供される。直接法は相互相関係数の一般式(定義式)に基づいて計算するものであり、2次元フーリエ変換法は相互相関関数が原関数のパワースペクトラムのフーリエ逆変換であるというウィナー・ヒンチンの定理を利用して計算するものである。
【0023】
図7は第1の方式の全体の流れの説明図である。c(x,y)が求める相互相関係数分布である。c(x,y)は複素信号なので、最終的に画像化するには絶対値処理が必要である。ここで、演算点X,Yループのステップ(カーネルROIの転移ピッチ)は、図8に示したカーネルROI(オブジェクトROIに同じ)のサイズの縦横(rx、ry)に対して、それぞれ1/4程度とすることは、c(x,y)の誤差が許容範囲に収まり、且つ計算時間の短縮化という観点から好ましい。この場合、計算しなかったc(x,y)は、補間する。
【0024】
図9は直接法の計算手順の説明図であり、図10は2次元フーリエ変換法の計算手順の説明図である。なお、前提条件として図8のように定義する。 *c はc の複素共役、 cc(τx,τy)は相互相関係数、csumは相互相関係数の x=τx 、 y=τy の要素値、ccとcsumは複素信号とする。
【0025】
図11は第2の方式の説明図である。或るカーネルROI(中心点(x,y) )に対してオブジェクトROIをサーチROI内で移動しながら各位置で、対応位置間での差分の総和の絶対値ε(xi,yj) 、及びこのε(xi,yj) のカーネルROI内の総和に対する比率ε´(xi,yj) を計算し、ε´(xi,yj) の最小値(または最大値)を当該カーネルROI(中心点(x,y) )における特徴量とする。このような計算をカーネルROIを転移しながら繰り返すことにより特徴量の2次元分布を生成する。
【0026】
上述したように本実施例によると、心筋の変形の程度を、心時相の異なる2フレームを用いて各フレームの局所領域内のRF受信信号の分布間での類似度(特徴量)として定量化(数値化)して、梗塞心筋を客観的に画像化できる。しかも、この類似度は、サーチ領域をサーチするので、梗塞心筋の転移を追従できる。
【0027】
本実施例は、次のように変形可能であり、上述の説明では2次元の局所領域間のRF受信信号の分布に基づいて類似性を見たものであるが、これを1次元追加して3次元の局所領域間での類似性に発展させたものである。これは、RF受信信号の位置情報に奥行き情報を加えることで、容易に実現できる。つまり、図12に示すように、図1の構成にプローブ1の位置を検出する位置検出器20を追加し、特徴量画像生成部7にこのプローブ1の位置情報を奥行き情報として与えることで、プローブ1を例えば平行移動しながら得られる多段層のRF受信信号群を1つの3次元データとして取り扱うことが可能となる。図13に、3次元画像処理を施したBモード画像と並列に表示した3次元の特徴量画像の一例を示す。この変型例によれば、不明な梗塞心筋の3次元的な転移を捕捉可能である。
本発明はその要旨を逸脱しない範囲において種々変形して実施可能である。なお上述ではBモード、カラーフローマッピングにも適用可能である。
【0028】
【発明の効果】
本発明は、被検体の2次元領域を超音波で走査することにより、1フレーム分の受信信号を繰り返し得、この受信信号を検波することにより画像を生成する超音波診断装置において、第1フレームの第1局所領域内の受信信号と、前記第1フレームと走査タイミングの異なる第2フレームの第2局所領域内の受信信号との間の類似性を表す特徴量の2次元分布を生成する特徴量2次元分布生成手段と、前記2次元分布を特徴量画像として表示する表示手段とを具備し、正常心筋の周期的な収縮に従って、組織が変型しないで単に位置が変化(転移)する梗塞部分を、第1フレームの第1局所領域内の受信信号と第2フレームの第2局所領域内の受信信号との間の類似性を高いとして、とらえることができるので、正常心筋と差別して梗塞心筋を特徴的に画像化することのできる超音波診断装置を提供できる。
【図面の簡単な説明】
【図1】本発明による超音波診断装置の一実施例の構成を示すブロック図。
【図2】図1の特徴量画像生成部のブロック図。
【図3】図2のDSCのブロック図。
【図4】特徴量画像の生成過程の概略的説明図。
【図5】輝度及び色相変換の設定変化を示す図。
【図6】特徴量画像の表示例を示す図。
【図7】特徴量を求める第1の方式の説明図。
【図8】相互相関係数計算に使われるパラメータの説明図。
【図9】相互相関係数を計算する直接法の計算手順の説明図。
【図10】相互相関係数を計算する2次元フーリエ変換法の計算手順の説明図。
【図11】特徴量を求める第2の方式の説明図。
【図12】本発明の超音波診断装置の変型例の構成を示すブロック図。
【図13】変型例による特徴量画像の表示例を示す図。
【符号の説明】
1…超音波プローブ、 2…送信系、
3…受信系、 4…信号処理部、
5…DSC、 6…表示器、
7…特徴量画像生成部。
Claims (3)
- 被検体の2次元又は3次元領域を超音波で走査することにより、1フレーム分の受信信号を繰り返し得、この受信信号を検波することにより画像を生成する超音波診断装置において、
第1フレームの第1局所領域内の受信信号と、前記第1フレームと走査タイミングの異なる第2フレーム内における前記第1局所領域に対して空間的に近傍する複数の第2局所領域内の受信信号各々との間で相互相関係数を計算し、前記計算した相互相関係数の中の最大値を用いて前記第1局所領域に関する前記第2局所領域に対する類似性を表す特徴量を定義し、前記第1局所領域の位置を変えながら前記特徴量の2次元分布を生成する特徴量2次元分布生成手段と、
前記2次元分布を特徴量画像として表示する表示手段とを具備することを特徴とする超音波診断装置。 - 被検体の2次元又は3次元領域を超音波で走査することにより、1フレーム分の受信信号を繰り返し得、この受信信号を検波することにより画像を生成する超音波診断装置において、
第1フレームの第1局所領域内の受信信号と、前記第1フレームと走査タイミングの異なる第2フレーム内における前記第1局所領域に対して空間的に近傍する複数の第2局所領域内の受信信号各々との間での信号差分の絶対値の総和を計算し、前記計算した総和の中の最小値を用いて前記第1局所領域に関する前記第2局所領域に対する類似性を表す特徴量を定義し、前記第1局所領域の位置を変えながら前記特徴量の2次元分布を生成する特徴量2次元分布生成手段と、
前記2次元分布を特徴量画像として表示する表示手段とを具備することを特徴とする超音波診断装置。 - 前記表示手段は前記2次元分布を輝度変換処理と色相変換処理の一方を介して表示することを特徴とする請求項1又は2記載の超音波診断装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31310294A JP3544722B2 (ja) | 1994-12-16 | 1994-12-16 | 超音波診断装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31310294A JP3544722B2 (ja) | 1994-12-16 | 1994-12-16 | 超音波診断装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH08164139A JPH08164139A (ja) | 1996-06-25 |
JP3544722B2 true JP3544722B2 (ja) | 2004-07-21 |
Family
ID=18037184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP31310294A Expired - Fee Related JP3544722B2 (ja) | 1994-12-16 | 1994-12-16 | 超音波診断装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3544722B2 (ja) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4659974B2 (ja) * | 2000-12-12 | 2011-03-30 | 株式会社東芝 | 超音波診断装置 |
JP4907798B2 (ja) * | 2001-08-24 | 2012-04-04 | 株式会社東芝 | 超音波診断装置 |
JP4733938B2 (ja) * | 2004-07-16 | 2011-07-27 | 株式会社東芝 | 超音波診断装置および超音波画像処理装置 |
JP4754838B2 (ja) * | 2005-02-03 | 2011-08-24 | 株式会社日立メディコ | 超音波診断装置 |
JP4890554B2 (ja) | 2006-08-21 | 2012-03-07 | 国立大学法人東北大学 | 超音波診断装置 |
WO2010024023A1 (ja) * | 2008-08-25 | 2010-03-04 | 株式会社 日立メディコ | 超音波診断装置及び超音波画像表示方法 |
JP4709937B2 (ja) * | 2010-10-01 | 2011-06-29 | 株式会社東芝 | 超音波診断装置及び画像処理装置 |
CN103906473B (zh) * | 2011-10-28 | 2016-01-06 | 日立阿洛卡医疗株式会社 | 超声波成像装置、超声波成像方法 |
US11475669B2 (en) * | 2020-07-30 | 2022-10-18 | Ncr Corporation | Image/video analysis with activity signatures |
-
1994
- 1994-12-16 JP JP31310294A patent/JP3544722B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH08164139A (ja) | 1996-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5715594B2 (ja) | フローパラメータイメージングのための方法及び装置 | |
JP4570115B2 (ja) | コヒーレンス・イメージングのための方法および装置 | |
US5785654A (en) | Ultrasound diagnostic apparatus | |
JP4942237B2 (ja) | 関心領域を画像内に配置する表示方法及びイメージングシステム | |
JP4958348B2 (ja) | 超音波撮像装置 | |
KR100749973B1 (ko) | Prf 조절 방법 및 장치, 및 초음파 촬상 장치 | |
US8740798B2 (en) | Three-dimensional ultrasonic diagnosis apparatus | |
WO2001066014A1 (fr) | Dispositif d'imagerie a ultrasons | |
US20070038103A1 (en) | Apparatus for obtaining ultrasonic image and method of obtaining ultrasonic image | |
JP3144819B2 (ja) | 超音波診断装置 | |
JP2005152648A (ja) | 運動適応型空間合成の方法及びシステム | |
JPH08173422A (ja) | 超音波診断装置 | |
JP3544722B2 (ja) | 超音波診断装置 | |
JP3763924B2 (ja) | 超音波診断装置 | |
JP2012010943A (ja) | 超音波診断装置及び超音波診断方法 | |
JPH078492A (ja) | 超音波診断装置 | |
JP7366221B2 (ja) | 解剖学的、機能的及び血行動態イメージングのためのトリプルモード超音波イメージング | |
JP2006000421A (ja) | 超音波血流イメージング装置 | |
JP2003339698A (ja) | 超音波診断装置 | |
JP2007020915A (ja) | 超音波診断装置 | |
JP3695807B2 (ja) | 2次元ドプラ超音波診断装置 | |
JP3558586B2 (ja) | 超音波診断装置 | |
JPH08117227A (ja) | 超音波診断装置 | |
JP2702983B2 (ja) | 超音波診断装置 | |
JP2009112491A (ja) | 超音波診断装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20040330 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040406 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080416 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090416 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100416 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100416 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110416 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130416 Year of fee payment: 9 |
|
LAPS | Cancellation because of no payment of annual fees |