JP2849159B2 - 超音波診断装置 - Google Patents
超音波診断装置Info
- Publication number
- JP2849159B2 JP2849159B2 JP12236090A JP12236090A JP2849159B2 JP 2849159 B2 JP2849159 B2 JP 2849159B2 JP 12236090 A JP12236090 A JP 12236090A JP 12236090 A JP12236090 A JP 12236090A JP 2849159 B2 JP2849159 B2 JP 2849159B2
- Authority
- JP
- Japan
- Prior art keywords
- depth
- diagnostic apparatus
- ultrasonic diagnostic
- dimensional
- shear
- 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
Description
【発明の詳細な説明】 〔概要〕 生体組織の弾性特性を算出表示するための手段を備え
た超音波診断装置に関し、 パルスドプラ法によって得られた各深さにおける組織
の変位速度を深さ方向に微分することによって、組織の
伸び縮みの度合いを示す「ずり速度」を算出表示するこ
とを目的とし、 被検体内の同一方向に複数回超音波パルスを送信し、
該同一方向から受信された複数本を受信信号を使って、
各深さに於ける被検体内組織の変位速度を検出する速度
検出手段と、該検出した速度を深さ方向に微分する深さ
方向微分手段と、前記深さ方向微分手段によって得られ
たずり速度の実効値を算出する実効値算出手段と、前記
深さ方向微分手段によって得られたずり速度を時間と深
さの二次元関数、又は、空間的な二次元分布像として、
また、前記実効値算出手段によって得られた結果を数値
として表示するための表示手段を有する構成とする。
た超音波診断装置に関し、 パルスドプラ法によって得られた各深さにおける組織
の変位速度を深さ方向に微分することによって、組織の
伸び縮みの度合いを示す「ずり速度」を算出表示するこ
とを目的とし、 被検体内の同一方向に複数回超音波パルスを送信し、
該同一方向から受信された複数本を受信信号を使って、
各深さに於ける被検体内組織の変位速度を検出する速度
検出手段と、該検出した速度を深さ方向に微分する深さ
方向微分手段と、前記深さ方向微分手段によって得られ
たずり速度の実効値を算出する実効値算出手段と、前記
深さ方向微分手段によって得られたずり速度を時間と深
さの二次元関数、又は、空間的な二次元分布像として、
また、前記実効値算出手段によって得られた結果を数値
として表示するための表示手段を有する構成とする。
本発明は、パルスドプラ法によって得られた各深さに
おける組織の変位速度を深さ方向に微分することによっ
て、組織の伸び縮みの度合いを示す弾性的組織特性を算
出表示する超音波診断装置に関するものである。
おける組織の変位速度を深さ方向に微分することによっ
て、組織の伸び縮みの度合いを示す弾性的組織特性を算
出表示する超音波診断装置に関するものである。
従来、組織の弾性的特性を表示する超音波診断装置と
して、第７図のように生体外部から低周波振動を与え、
生体内部を伝播する振動波の振幅と位相分布を計測しよ
うとするもの（第52回日本超音波医学会論文集52−143,
P287−288、及び、第53回日本超音波医学会論文集53−8
3,P271−272）や、第８図のように解析信号の空間相関
関数を用いた不均一組織の微小変位の計測をしようとす
るもの（第54回日本超音波医学会論文集54−116,P359−
360）が提案されている。
して、第７図のように生体外部から低周波振動を与え、
生体内部を伝播する振動波の振幅と位相分布を計測しよ
うとするもの（第52回日本超音波医学会論文集52−143,
P287−288、及び、第53回日本超音波医学会論文集53−8
3,P271−272）や、第８図のように解析信号の空間相関
関数を用いた不均一組織の微小変位の計測をしようとす
るもの（第54回日本超音波医学会論文集54−116,P359−
360）が提案されている。
第７図は、低周波振動を加えた時に、生体内部を伝播
する振動波の振動振幅と位相分布の計測するためのプロ
ーブ位置と機械的振動子の位置関係を示したものであ
る。機械的振動子20は、平板を通して、媒体19に周波数
10kHzから1kHz程度を振動を加える。また、プローブ11
から送波された超音波パルスは媒体内部で反射し、同プ
ローブ11で受信される。この受信信号を使って、振動波
22の振動振幅と位相分布を求めようとするが、第７図で
ある。23は、プローブと体表との間のコンタクト剤、24
は、機械的振動子20の振動を吸収するためのスプリング
である。
する振動波の振動振幅と位相分布の計測するためのプロ
ーブ位置と機械的振動子の位置関係を示したものであ
る。機械的振動子20は、平板を通して、媒体19に周波数
10kHzから1kHz程度を振動を加える。また、プローブ11
から送波された超音波パルスは媒体内部で反射し、同プ
ローブ11で受信される。この受信信号を使って、振動波
22の振動振幅と位相分布を求めようとするが、第７図で
ある。23は、プローブと体表との間のコンタクト剤、24
は、機械的振動子20の振動を吸収するためのスプリング
である。
第８図は、解析信号の二次元位相関関数を用いた組織
変位を計測する方法である。第８図において、11は、超
音波を送受信するための超音波プローブで、送受信回路
10につながっている。10は、ビームフォーマ等を含む送
受信回路で、プローブ11の構成要素である振動子群を駆
動するための信号を生成したり、各振動子で受信された
信号から特定方向にフォーカスした受信信号１を生成す
る。受信信号１は、次段の解析信号生成手段12に送られ
る。12では、受信信号をxr（ｔ）とすると、xr（ｔ）と
直交する信号xi（ｔ）を生成し、複素信号ｚ（ｔ）＝xr
（ｔ）＋jxi（ｔ）をメモリ13に格納する。13は、全て
の走査線に対応する受信複素信号ｚ（ｔ）を複数フレー
ム分以上格納する。
変位を計測する方法である。第８図において、11は、超
音波を送受信するための超音波プローブで、送受信回路
10につながっている。10は、ビームフォーマ等を含む送
受信回路で、プローブ11の構成要素である振動子群を駆
動するための信号を生成したり、各振動子で受信された
信号から特定方向にフォーカスした受信信号１を生成す
る。受信信号１は、次段の解析信号生成手段12に送られ
る。12では、受信信号をxr（ｔ）とすると、xr（ｔ）と
直交する信号xi（ｔ）を生成し、複素信号ｚ（ｔ）＝xr
（ｔ）＋jxi（ｔ）をメモリ13に格納する。13は、全て
の走査線に対応する受信複素信号ｚ（ｔ）を複数フレー
ム分以上格納する。
ここで、あるフレームに於ける、走査線ｍ、深さｎに
於ける二次元複素信号データを Ａ（m,n）＝Ar（m,n）＋jAi（m,n） 式（１） と表現する。また、別のフレームに於ける、走査線ｍ、
深さｎに於ける二次元複素信号データを Ｂ（m,n）＝Br（m,n）＋jBi（m,n） 式（２） と表現すると、（m,n）の周りでの二次元相互相関関数
Ｃ（τ，ρ|m,n）は次式で定義できる。
於ける二次元複素信号データを Ａ（m,n）＝Ar（m,n）＋jAi（m,n） 式（１） と表現する。また、別のフレームに於ける、走査線ｍ、
深さｎに於ける二次元複素信号データを Ｂ（m,n）＝Br（m,n）＋jBi（m,n） 式（２） と表現すると、（m,n）の周りでの二次元相互相関関数
Ｃ（τ，ρ|m,n）は次式で定義できる。
ここで、上式中＜ ＞は、二次元空間的な平均操作を
示す。上式の最大値を計算し、複素空間で、（1,0）か
らの距離を計算すると変位ベクトルを計算できる。
示す。上式の最大値を計算し、複素空間で、（1,0）か
らの距離を計算すると変位ベクトルを計算できる。
14は、（３）式の二次元相互相関関数を計算するため
のもので、15は、任意の位置で算出した二次元相互関関
数を格納するためのメモリである。微小変位検出手段16
は、変位ベクトルを計算し、結果を表示手段17に送る。
各位置で計算された変位ベクトルは、第８図（ｂ）に示
すごとく、ベクトルの傾きを移動方向、ベクトルの長さ
を移動量として表示する。
のもので、15は、任意の位置で算出した二次元相互関関
数を格納するためのメモリである。微小変位検出手段16
は、変位ベクトルを計算し、結果を表示手段17に送る。
各位置で計算された変位ベクトルは、第８図（ｂ）に示
すごとく、ベクトルの傾きを移動方向、ベクトルの長さ
を移動量として表示する。
しかし、第７図の方法においては、外部から振動を与
えるための機械系が必要になり、実用的でないこと、ま
た、第８図の方法では、二次元相互相関関数を計算する
のに十分に広い領域を必要とするため高い空間分解能を
得ることが難しく、計算処理に時間がかかるという問題
があった。
えるための機械系が必要になり、実用的でないこと、ま
た、第８図の方法では、二次元相互相関関数を計算する
のに十分に広い領域を必要とするため高い空間分解能を
得ることが難しく、計算処理に時間がかかるという問題
があった。
ところで、生体内組織では、拍動等の動きによって組
織が変位している。この組織の変位速度は、場所によっ
て硬さが異なるために、一様ではなく、ばらついてい
る。変位の速度の場所によるばらつきは、組織の伸び縮
みを意味する。例えば、ある深さに注目し、この深さよ
りプローブ寄りのところで、プローブに向かう変位速度
が、その深さより離れた位置での変位速度より速けれ
ば、その深さでは、組織が伸びていることを意味する。
もし、正常組織よりも硬い異常組織が、正常組織の中に
存在した場合、硬い異常組織内部での変位速度の場所に
よるばらつきは、正常組織のそれに比べて小さいと考え
られる。
織が変位している。この組織の変位速度は、場所によっ
て硬さが異なるために、一様ではなく、ばらついてい
る。変位の速度の場所によるばらつきは、組織の伸び縮
みを意味する。例えば、ある深さに注目し、この深さよ
りプローブ寄りのところで、プローブに向かう変位速度
が、その深さより離れた位置での変位速度より速けれ
ば、その深さでは、組織が伸びていることを意味する。
もし、正常組織よりも硬い異常組織が、正常組織の中に
存在した場合、硬い異常組織内部での変位速度の場所に
よるばらつきは、正常組織のそれに比べて小さいと考え
られる。
本発明は、パルスドプラ法を用いて各場所の組織の変
位速度を検出し、場所による変位速度の変化の度合いを
算出することによって、弾性的特性を表示することを目
的としている。
位速度を検出し、場所による変位速度の変化の度合いを
算出することによって、弾性的特性を表示することを目
的としている。
上記課題を解決するための具体的手段を第２図と第３
図を使って説明する。
図を使って説明する。
従来、パルスドプラ法を使って血流の速度を検出する
方法が知られている。本発明においては、組織の変位速
度を検出するのにこの手法を用いる。第３図中２は、こ
の検出方法を説明するブロック図である。
方法が知られている。本発明においては、組織の変位速
度を検出するのにこの手法を用いる。第３図中２は、こ
の検出方法を説明するブロック図である。
第３図において、２−１〜２−６は直交検波器を構成
する。従って、ローパスフィルタ（LPF）２−５の出力
は、受信信号１の直交検波出力の実部成分である。ま
た、ローパスフィルタ（LPF）２−６の出力は、受信信
号１の直交検波出力の虚部成分である。
する。従って、ローパスフィルタ（LPF）２−５の出力
は、受信信号１の直交検波出力の実部成分である。ま
た、ローパスフィルタ（LPF）２−６の出力は、受信信
号１の直交検波出力の虚部成分である。
同一走査線方向に、同期Ｔの間隔でＮ回超音波パルス
を送波した時の受信信号の直交検波出力の実部成分を第
２図（ａ）中、R1,R2,……,RNで示す。同様に、虚部成
分を第２図（ａ）中、I1,I2,……,INで示す。これらの
信号系列は、それぞれメモリ２−９、２−10に一旦格納
される。ここで、パルスドプラ法を使うと、ある深さｚ
に於ける深さ方向の変位速度は、以下のように算出する
ことができる。
を送波した時の受信信号の直交検波出力の実部成分を第
２図（ａ）中、R1,R2,……,RNで示す。同様に、虚部成
分を第２図（ａ）中、I1,I2,……,INで示す。これらの
信号系列は、それぞれメモリ２−９、２−10に一旦格納
される。ここで、パルスドプラ法を使うと、ある深さｚ
に於ける深さ方向の変位速度は、以下のように算出する
ことができる。
第２図（ａ）中の点線矢印で示すように、深さｚにお
いて時間方向に読みだした複素信号系列Vn（ｚ）を Vn（ｚ）＝（R1（ｚ），……,RN-1（ｚ）） ＋ｊ（I1（ｚ），……,IN-1（ｚ）） 式（４） と表現する。また、周期Ｔだけずらして、深さｚにおい
て時間方向に読みだした複素信号系列Vn+1（ｚ）を Vn+1（ｚ）＝（R2（ｚ），……,RN（ｚ）） ＋ｊ（I2（ｚ），……,IN（ｚ）） 式（５） と表現する。
いて時間方向に読みだした複素信号系列Vn（ｚ）を Vn（ｚ）＝（R1（ｚ），……,RN-1（ｚ）） ＋ｊ（I1（ｚ），……,IN-1（ｚ）） 式（４） と表現する。また、周期Ｔだけずらして、深さｚにおい
て時間方向に読みだした複素信号系列Vn+1（ｚ）を Vn+1（ｚ）＝（R2（ｚ），……,RN（ｚ）） ＋ｊ（I2（ｚ），……,IN（ｚ）） 式（５） と表現する。
第３図中、複素自己相関器２−11は、上記２つの複素
信号系列の相関係数を算出するものである。即ち、相関
係数Yn（ｚ）は、 Yn（ｚ）＝Yn+1（ｚ）Vn ＊（ｚ） ＝|Yn（ｚ）|exp（ｊΔθｎ）（ｚ）） 式（６） ≡YnR（ｚ）＋jYnI（ｚ） 式（７） となる。ここで、＊は複素共役、YnR（ｚ）は相関係数
の実数成分、YnI（ｚ）は相関係数の虚数成分、Δθｎ
（ｚ）は、ドプラ効果によって生じた位相差である。
信号系列の相関係数を算出するものである。即ち、相関
係数Yn（ｚ）は、 Yn（ｚ）＝Yn+1（ｚ）Vn ＊（ｚ） ＝|Yn（ｚ）|exp（ｊΔθｎ）（ｚ）） 式（６） ≡YnR（ｚ）＋jYnI（ｚ） 式（７） となる。ここで、＊は複素共役、YnR（ｚ）は相関係数
の実数成分、YnI（ｚ）は相関係数の虚数成分、Δθｎ
（ｚ）は、ドプラ効果によって生じた位相差である。
上式より、ドプラシフト周波数fdn（ｚ）は、 fdn（ｚ）＝tan-1［YnI（ｚ）/YnR（ｚ）］/T/（２π） ＝Δθｎ（ｚ）/T/（２π） 式（８） となる。従って、組織の深さ方向の変位速度は、 vn（ｚ）＝fdn（ｚ）C/（2f0） 式（９） となる。但し、f0は、送信超音波パルスの中心周波数、
Ｃは音速である。
Ｃは音速である。
第３図中、速度計算部２−12は、（９）式を計算する
ためのものである。全ての深さについて計算した組織の
深さ方向の変位速度の概念図を第２図（ｂ）に示す。
ためのものである。全ての深さについて計算した組織の
深さ方向の変位速度の概念図を第２図（ｂ）に示す。
組織の伸び縮みの度合いを示す弾性的組織パラメータ
は、変位速度の場所による変化の度合いで表現できるか
ら、第２図（ｂ）の変位速度を深さ方向に微分すること
によって得ることができる。我々は、これを組織の伸び
縮みを表す「ずり速度」Sn（ｚ）として以下のように定
義する。
は、変位速度の場所による変化の度合いで表現できるか
ら、第２図（ｂ）の変位速度を深さ方向に微分すること
によって得ることができる。我々は、これを組織の伸び
縮みを表す「ずり速度」Sn（ｚ）として以下のように定
義する。
Sn（ｚ）＝dvn（ｚ）/dz 式（10） 第２図（ｂ）の変位速度vn（ｚ）を（10）式のように
深さ方向に微分したずり速度Sn（ｚ）を概念的に第２図
（ｃ）に示す。
深さ方向に微分したずり速度Sn（ｚ）を概念的に第２図
（ｃ）に示す。
同一走査線方向に送波する回数を増やし（即ち、Ｎ回
以上にする）、ずり速度Sn（ｚ）をn,n＋1,n＋2,…につ
いて繰り返し算出し、ある走査線についての時間と深さ
の二次元関数として、Sn（ｚ）,Sn+1（ｚ），…を表示
することができる。また、走査線をスキャンすることに
よって、各走査線に対して得られたSn（ｚ）を空間的二
次元分布像を表示することができる。更に、従来の超音
波診断装置のＢモード像、Ｍモード像と組み合わせて表
示することによって、対象部位の伸び縮みの様子を的確
に判断することができるようになる。
以上にする）、ずり速度Sn（ｚ）をn,n＋1,n＋2,…につ
いて繰り返し算出し、ある走査線についての時間と深さ
の二次元関数として、Sn（ｚ）,Sn+1（ｚ），…を表示
することができる。また、走査線をスキャンすることに
よって、各走査線に対して得られたSn（ｚ）を空間的二
次元分布像を表示することができる。更に、従来の超音
波診断装置のＢモード像、Ｍモード像と組み合わせて表
示することによって、対象部位の伸び縮みの様子を的確
に判断することができるようになる。
また、ある走査線において、深さｚ＝ｊΔz,j＝0,1,
2,……（Δｚは深さ方向のサンプル間隔）について計算
された｛Sn（ｊ）,Sn+1j），…｝について、深さの範囲
｛z1,z2］、時間の範囲［t1,t2］を指定すると、ずり速
度の実効値Ｐは次式で定義できる。
2,……（Δｚは深さ方向のサンプル間隔）について計算
された｛Sn（ｊ）,Sn+1j），…｝について、深さの範囲
｛z1,z2］、時間の範囲［t1,t2］を指定すると、ずり速
度の実効値Ｐは次式で定義できる。
avは、深さの範囲［z1,z2］、時間の範囲［t1,t2］内
の平均値で、 となる。また、簡単のためのav＝０として計算しても良
い。但し、t1＝Tn1、t2＝Tn2、z1＝j1Δｚ、z2＝j2Δｚ
である。
の平均値で、 となる。また、簡単のためのav＝０として計算しても良
い。但し、t1＝Tn1、t2＝Tn2、z1＝j1Δｚ、z2＝j2Δｚ
である。
（11）式のようなずり速度の実効値は、組織の弾性的
特性の１つのパラメータとなる。また、（11）式におい
て、特定の深さｚについてのみ、計算しても良いことは
いうまでもない。
特性の１つのパラメータとなる。また、（11）式におい
て、特定の深さｚについてのみ、計算しても良いことは
いうまでもない。
本発明によると、組織の弾性的特性であるずり速度
を、ある走査線における時間と深さの二次元的関数、あ
る関心領域における空間的な二次元分布像、或いは、実
効値として、従来の超音波診断装置のＢモード像、Ｍモ
ード像と組み合わせて表示されるようになる。これによ
って、組織の位置と硬さの度合いの情報を提供できるよ
うになる。
を、ある走査線における時間と深さの二次元的関数、あ
る関心領域における空間的な二次元分布像、或いは、実
効値として、従来の超音波診断装置のＢモード像、Ｍモ
ード像と組み合わせて表示されるようになる。これによ
って、組織の位置と硬さの度合いの情報を提供できるよ
うになる。
次に、第１図〜第６図を使って本発明を順次説明す
る。
る。
第１図に、本発明の構成を示す。11は、媒体19に対し
て、超音波を送受信するための超音波プローブで、送受
信回路10につながっている。10は、ビームフォーマ等を
含む送受信回路で、プローブ11の構成要素である振動子
群を駆動するための信号を生成したり、各振動子で受信
された信号から特定方向にフォーカスした受信信号１を
生成する。受信信号１は、速度検出手段２、Ｂモード像
生成手段８、Ｍモード像生成手段９の入力になる。10、
11、２、８及び９は、従来の超音波診断装置に備わって
いるものである。速度検出手段２は、（９）式を算出す
るためのもので、直交検波回路等から構成され、任意の
深さに於ける深さ方向の組織の変位速度を検出する。メ
モリ３は、指定した走査線における深さ方向の組織の変
位速度、或いは、ある関心領域内の複数走査線における
組織の変位速度を格納するためのものである。深さ方向
微分手段４は、（10）式を算出するためのもので、手段
２によって算出した変位速度を深さ方向に微分する。実
効値算出手段５は、（11）式を算出するためのもので、
手段４によって算出したずり速度の実効値を算出する。
表示手段６は、指定した走査線におけるずり速度の時間
と深さの二次元的関数、走査線をスキャンすることによ
って得られるある関心領域における空間的なずり速度の
二次元分布像、または、手段５によって算出したずり速
度の実効値を従来の超音波診断装置のＢモード像、Ｍモ
ード像と組み合わせて表示するためのものである。走査
線またはROI指定手段７は、ずり速度を計算するための
走査線、関心領域を指定するためのROI、ずり速度の実
効値を算出するための深さの範囲や時間区間を指定する
ための手段である。第２図は、第１図の本発明を概念的
に説明するものである。
て、超音波を送受信するための超音波プローブで、送受
信回路10につながっている。10は、ビームフォーマ等を
含む送受信回路で、プローブ11の構成要素である振動子
群を駆動するための信号を生成したり、各振動子で受信
された信号から特定方向にフォーカスした受信信号１を
生成する。受信信号１は、速度検出手段２、Ｂモード像
生成手段８、Ｍモード像生成手段９の入力になる。10、
11、２、８及び９は、従来の超音波診断装置に備わって
いるものである。速度検出手段２は、（９）式を算出す
るためのもので、直交検波回路等から構成され、任意の
深さに於ける深さ方向の組織の変位速度を検出する。メ
モリ３は、指定した走査線における深さ方向の組織の変
位速度、或いは、ある関心領域内の複数走査線における
組織の変位速度を格納するためのものである。深さ方向
微分手段４は、（10）式を算出するためのもので、手段
２によって算出した変位速度を深さ方向に微分する。実
効値算出手段５は、（11）式を算出するためのもので、
手段４によって算出したずり速度の実効値を算出する。
表示手段６は、指定した走査線におけるずり速度の時間
と深さの二次元的関数、走査線をスキャンすることによ
って得られるある関心領域における空間的なずり速度の
二次元分布像、または、手段５によって算出したずり速
度の実効値を従来の超音波診断装置のＢモード像、Ｍモ
ード像と組み合わせて表示するためのものである。走査
線またはROI指定手段７は、ずり速度を計算するための
走査線、関心領域を指定するためのROI、ずり速度の実
効値を算出するための深さの範囲や時間区間を指定する
ための手段である。第２図は、第１図の本発明を概念的
に説明するものである。
受信信号１は、速度検出手段２内において先ず直交検
波される。ここで、第２図を使って、手段２の動作を説
明する。
波される。ここで、第２図を使って、手段２の動作を説
明する。
第２図（ａ）は、同一方向にＮ回送受信した時の直交
検波された複素信号系列を示している。ここで、複素信
号系列の内、ある深さｚに於ける複素信号系列Vn（ｚ）
（（４）式）を得る。次に、送波周期Ｔずれた、同じ深
さｚに於ける信号複素信号系列Vn+1（ｚ）（（５）式）
を得る。次に、（６）式に従って、複素自己相関係数を
計算し、（９）式から、深さｚにおける組織の変位速度
Vn（ｚ）を算出する。全深さについて、vn（ｚ）を計算
したのが、図２図（ｂ）である。
検波された複素信号系列を示している。ここで、複素信
号系列の内、ある深さｚに於ける複素信号系列Vn（ｚ）
（（４）式）を得る。次に、送波周期Ｔずれた、同じ深
さｚに於ける信号複素信号系列Vn+1（ｚ）（（５）式）
を得る。次に、（６）式に従って、複素自己相関係数を
計算し、（９）式から、深さｚにおける組織の変位速度
Vn（ｚ）を算出する。全深さについて、vn（ｚ）を計算
したのが、図２図（ｂ）である。
次に、深さ方向微分手段４においは、一旦メモリ３に
蓄えられた組織の変位速度Vn（ｚ）を（10）式に従って
深さ方向に微分する。第２図（ｃ）は、この時のずり速
度を概念的に示したものである。
蓄えられた組織の変位速度Vn（ｚ）を（10）式に従って
深さ方向に微分する。第２図（ｃ）は、この時のずり速
度を概念的に示したものである。
第３図は、速度検出手段２及び深さ方向微分手段４を
より詳細に説明したものである。第３図において、速度
検出手段２は、直交検波器２−１〜２−６、直交検波信
号の実部成分をAD変換するためのADC2−７、直交検波信
号の虚部成分をAD変換するためのADC2−８、実数成分を
格納するメモリ２−９、虚数成分を格納するメモリ２−
10、複素自己相関器２−11、速度計算部２−12から構成
される。２−１は、受信信号１に参照信号（cosω0t）
２−３を乗算するためのミキサ、２−５は、前記ミキサ
２−１の出力の上側波帯成分をカットするためのローパ
スフィルタである。同様に、２−２は、受信信号１に参
照信号（−sinω0t）２−４を乗算するためのミキサ、
２−６は、前記ミキサ２−２の出力の上側波帯成分をカ
ットするためのローパスフィルタである。２−５及び２
−６の出力は、直交検波信号の実数部と虚数部に対応す
る。直交検波信号の実数部と虚数部は、一旦、メモリ２
−９、２−10に格納される。２−11は、同一深さの複素
信号系列に対して、複素自己相関係数を（７）式に従っ
て計算するものである。また、２−11において算出され
た複素自己相関係数から組織の変位速度vn（ｚ）が
（９）式に従って算出される。
より詳細に説明したものである。第３図において、速度
検出手段２は、直交検波器２−１〜２−６、直交検波信
号の実部成分をAD変換するためのADC2−７、直交検波信
号の虚部成分をAD変換するためのADC2−８、実数成分を
格納するメモリ２−９、虚数成分を格納するメモリ２−
10、複素自己相関器２−11、速度計算部２−12から構成
される。２−１は、受信信号１に参照信号（cosω0t）
２−３を乗算するためのミキサ、２−５は、前記ミキサ
２−１の出力の上側波帯成分をカットするためのローパ
スフィルタである。同様に、２−２は、受信信号１に参
照信号（−sinω0t）２−４を乗算するためのミキサ、
２−６は、前記ミキサ２−２の出力の上側波帯成分をカ
ットするためのローパスフィルタである。２−５及び２
−６の出力は、直交検波信号の実数部と虚数部に対応す
る。直交検波信号の実数部と虚数部は、一旦、メモリ２
−９、２−10に格納される。２−11は、同一深さの複素
信号系列に対して、複素自己相関係数を（７）式に従っ
て計算するものである。また、２−11において算出され
た複素自己相関係数から組織の変位速度vn（ｚ）が
（９）式に従って算出される。
ここで、血流速度を検出する場合、複素自己相関器２
−11の前段にMTIフィルタを入れるのが通常であるが、
本発明では必要としない。
−11の前段にMTIフィルタを入れるのが通常であるが、
本発明では必要としない。
算出した組織の変位速度は一旦メモリ３に格納され
る。一旦格納した変位速度は、深さ方向の時系列として
読みだされ、FIRフィルタ４−１によって、深さ方向の
微分処理が行われる。FIRフィルタの係数列の１例とし
て第３図（ｂ）のようにすると、（10）式に相当する演
算をすることができる。この係数を第３図（ｃ）のFIR
フィルタの係数列４−１−２−０〜４−１−２−ｍとす
る。第３図（ｃ）において、データシフトレジスタ４−
１−１は、入力の変位速度の信号系列vn（ｉ）（ｉ＝0,
1,2,……）を図示しないシステムクロックに従って１段
ずつシフトするもので、シフトされた信号系列｛v
n（ｉ）,vn（ｉ−１）,vn（ｉ−２），……｝にそれぞ
れ係数｛am,am-1,……｝を乗算した結果を加算器４−１
−３にて加算する。加算結果は、変位速度vn（ｉ）を深
さ方向に微分した結果に相当し、ずり速度Sn（ｊ）とな
る。深さ方向のデータ位置を示すｊは例えばｊ＝ｉ−m/
2である。
る。一旦格納した変位速度は、深さ方向の時系列として
読みだされ、FIRフィルタ４−１によって、深さ方向の
微分処理が行われる。FIRフィルタの係数列の１例とし
て第３図（ｂ）のようにすると、（10）式に相当する演
算をすることができる。この係数を第３図（ｃ）のFIR
フィルタの係数列４−１−２−０〜４−１−２−ｍとす
る。第３図（ｃ）において、データシフトレジスタ４−
１−１は、入力の変位速度の信号系列vn（ｉ）（ｉ＝0,
1,2,……）を図示しないシステムクロックに従って１段
ずつシフトするもので、シフトされた信号系列｛v
n（ｉ）,vn（ｉ−１）,vn（ｉ−２），……｝にそれぞ
れ係数｛am,am-1,……｝を乗算した結果を加算器４−１
−３にて加算する。加算結果は、変位速度vn（ｉ）を深
さ方向に微分した結果に相当し、ずり速度Sn（ｊ）とな
る。深さ方向のデータ位置を示すｊは例えばｊ＝ｉ−m/
2である。
第４図は、ある走査線について算出したずり速度の時
間と深さの二次元関数の表示例を示したものである。
間と深さの二次元関数の表示例を示したものである。
第４図（ａ）は、従来の超音波診断装置で得られるＢ
モード像である。ここで、ずり速度を算出するために、
走査線またはROI指定手段７によって指定した走査線方
向を示す線を表示手段６によってＢモード像上に重ねて
表示する。この走査線に於けるずり速度のデータセット
（Sn（ｚ）,Sn+1（ｚ），…）を深さｚ方向、時間ｎ方
向の関数として表示した例が、第４図（ｂ）、及び
（ｃ）である。（ｂ）は、ずり速度の値を振幅として、
また、（ｃ）はずり速度をグレースケール或いはカラー
画像として表示したもので、同図（ｄ）は、（ｃ）のず
り速度値に対応するグレースケール、または、カラーを
表示したものである。また、第４図（ｅ）のように、対
応するＭモード像を同時に表示してもよい。
モード像である。ここで、ずり速度を算出するために、
走査線またはROI指定手段７によって指定した走査線方
向を示す線を表示手段６によってＢモード像上に重ねて
表示する。この走査線に於けるずり速度のデータセット
（Sn（ｚ）,Sn+1（ｚ），…）を深さｚ方向、時間ｎ方
向の関数として表示した例が、第４図（ｂ）、及び
（ｃ）である。（ｂ）は、ずり速度の値を振幅として、
また、（ｃ）はずり速度をグレースケール或いはカラー
画像として表示したもので、同図（ｄ）は、（ｃ）のず
り速度値に対応するグレースケール、または、カラーを
表示したものである。また、第４図（ｅ）のように、対
応するＭモード像を同時に表示してもよい。
ここで、第４図（ｂ）内に示すように、走査線または
ROI指定手段７によって、深さの範囲［z1,z2］、時間の
範囲［t1,t2］を指定する、（11）式に従って、ずり速
度の実効値を手段５によって算出することができる。
ROI指定手段７によって、深さの範囲［z1,z2］、時間の
範囲［t1,t2］を指定する、（11）式に従って、ずり速
度の実効値を手段５によって算出することができる。
第５図は、肝臓のずり速度を第４図（ｂ）の如く、時
間と深さの二次元関数として表示した例である。第５図
（ａ）は、正常肝臓例、第５図（ｂ）は、深さＡの位置
に転移癌（転移した腫瘍）がある場合である。正常な肝
臓の例（ａ）では、各深さにおいて、ずり速度が時間的
に変化していく様子が観測できる。また、転移癌（ｂ）
の場合は、腫瘍内部の深さＡにおけるずり速度の時間的
変化が非常に小さくなる。これは、腫瘍内部が硬く、伸
び縮みが小さいことを示しているもの考えられる。更
に、腫瘍と正常組織との境界に当たる深さＢでは、速度
変化が非常に大きいことを示している。このように、本
発明によって、ずり速度を表示すると、正常組織と腫瘍
組織の違いを明瞭に判断することができる。
間と深さの二次元関数として表示した例である。第５図
（ａ）は、正常肝臓例、第５図（ｂ）は、深さＡの位置
に転移癌（転移した腫瘍）がある場合である。正常な肝
臓の例（ａ）では、各深さにおいて、ずり速度が時間的
に変化していく様子が観測できる。また、転移癌（ｂ）
の場合は、腫瘍内部の深さＡにおけるずり速度の時間的
変化が非常に小さくなる。これは、腫瘍内部が硬く、伸
び縮みが小さいことを示しているもの考えられる。更
に、腫瘍と正常組織との境界に当たる深さＢでは、速度
変化が非常に大きいことを示している。このように、本
発明によって、ずり速度を表示すると、正常組織と腫瘍
組織の違いを明瞭に判断することができる。
第６図は、従来の超音波診断装置で得られるＢモード
像上に、ずり速度の空間的二次元分布像を重ねて表示し
たものである。ずり速度を算出する範囲は、図中、点線
のROIで示した。第６図の例では、Ｂモード像に重ねて
表示したが、Ｂモード像と、ずり速度の空間的二次元像
を分けて表示しても良いことは言うまでもない。
像上に、ずり速度の空間的二次元分布像を重ねて表示し
たものである。ずり速度を算出する範囲は、図中、点線
のROIで示した。第６図の例では、Ｂモード像に重ねて
表示したが、Ｂモード像と、ずり速度の空間的二次元像
を分けて表示しても良いことは言うまでもない。
以上説明したように、本発明によれば、組織の弾性的
特性であるずり速度を、ある走査線における時間と深さ
の二次元的関数、ある関心領域における空間的な二次元
分布像、或いは、実効値として、従来の超音波診断装置
のＢモード像、Ｍモード像と組み合わせて表示されるよ
うになる。これによって、組織の位置と硬さの度合いの
情報を提供できるようになる。
特性であるずり速度を、ある走査線における時間と深さ
の二次元的関数、ある関心領域における空間的な二次元
分布像、或いは、実効値として、従来の超音波診断装置
のＢモード像、Ｍモード像と組み合わせて表示されるよ
うになる。これによって、組織の位置と硬さの度合いの
情報を提供できるようになる。
第１図は、本発明の１実施例、第２図は、本発明の１実
施例を概念的に説明する図、第３図は、第１図の構成を
より詳細に説明する図、第４図は、表示例、第５図は、
正常肝臓と転移肝癌に本手法を適用した場合の例、第６
図は、表示の別の例、第７図及び第８図は従来の実施例
である。 １……受信信号 ２……速度検出手段 ３……メモリ ４……深さ方向微分手段 ５……実効値算出手段 ６……表示手段 ７……ROI指定手段 ８……Ｂモード像生成手段 ９……Ｍモード像生成手段 11……超音波プローブ 10……送受信回路
施例を概念的に説明する図、第３図は、第１図の構成を
より詳細に説明する図、第４図は、表示例、第５図は、
正常肝臓と転移肝癌に本手法を適用した場合の例、第６
図は、表示の別の例、第７図及び第８図は従来の実施例
である。 １……受信信号 ２……速度検出手段 ３……メモリ ４……深さ方向微分手段 ５……実効値算出手段 ６……表示手段 ７……ROI指定手段 ８……Ｂモード像生成手段 ９……Ｍモード像生成手段 11……超音波プローブ 10……送受信回路
Claims (6)
- 【請求項１】被検体内の同一方向に複数回超音波パルス
を送信し、該同一方向から受信された複数本の受信信号
から、各深さに於ける被検体内組織の変位速度を検出す
る速度検出手段（２）と、 該検出した速度を深さ方向に微分し、ずり速度を算出す
る深さ方向微分手段（４）と、 前記深さ方向微分手段（４）によって得られたずり速度
を時間と深さの二次元関数、又は、空間的な二次元分布
像として、表示するための表示手段（６）を有すること
を特徴とする超音波診断装置。 - 【請求項２】前記深さ方向微分手段（４）によって得ら
れたずり速度の空間的な二次元分布像をＢモード像と同
時に同一画面上に、又は、Ｂモード像に重ねて表示する
ことを特徴とする特許請求の範囲第１項記載の超音波診
断装置。 - 【請求項３】前記深さ方向微分手段（４）によって得ら
れたずり速度を時間と深さの二次元関数として、Ｂモー
ド像、或いは、Ｍモード像と同時に表示することを特徴
とする特許請求の範囲第一項記載の超音波診断装置。 - 【請求項４】走査線またはROI指定手段（７）によっ
て、走査線を指定し、 前記深さ方向微分手段は、前記走査線又はROI指定手段
（７）によって指定された超音波走査線上で、ずり速度
を算出し、前記算出されたずり速度を時間と深さの二次
元関数として表示することを特徴とする特許請求の範囲
第一項記載の超音波診断装置。 - 【請求項５】走査線またはROI指定手段（７）によっ
て、二次元的な関心領域（ROI）を指定し、該関心領域
内の上記手段（４）によって算出されたずり速度を空間
的な二次元分布像として表示することを特徴とする特許
請求の範囲第一項記載の超音波診断装置。 - 【請求項６】前記深さ方向微分手段（４）によって得ら
れたずり速度のうち、指定された時間範囲と深さ範囲の
ずり速度において、実効値を算出する実効値算出手段
（５）を有することを特徴とする特許請求の範囲第一項
記載の超音波診断装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12236090A JP2849159B2 (ja) | 1990-05-11 | 1990-05-11 | 超音波診断装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12236090A JP2849159B2 (ja) | 1990-05-11 | 1990-05-11 | 超音波診断装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0417842A JPH0417842A (ja) | 1992-01-22 |
| JP2849159B2 true JP2849159B2 (ja) | 1999-01-20 |
Family
ID=14833968
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP12236090A Expired - Lifetime JP2849159B2 (ja) | 1990-05-11 | 1990-05-11 | 超音波診断装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2849159B2 (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100450446C (zh) * | 2002-10-18 | 2009-01-14 | 株式会社日立医药 | 超声波诊断装置 |
| CN109470883A (zh) * | 2018-10-15 | 2019-03-15 | 同济大学 | 一种gps漂流式流速测量仪 |
Families Citing this family (25)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2291969A (en) * | 1993-04-19 | 1996-02-07 | Commw Scient Ind Res Org | Tissue characterisation using intravascular echoscopy |
| JP4179596B2 (ja) * | 2002-08-30 | 2008-11-12 | 株式会社日立メディコ | 超音波診断装置 |
| WO2005122907A1 (ja) * | 2004-06-22 | 2005-12-29 | Hitachi Medical Corporation | 超音波診断装置および弾性画像表示方法 |
| WO2008051639A2 (en) | 2006-10-25 | 2008-05-02 | Maui Imaging, Inc. | Method and apparatus to produce ultrasonic images using multiple apertures |
| JP4615528B2 (ja) * | 2007-01-22 | 2011-01-19 | 株式会社日立メディコ | 超音波診断装置 |
| US9282945B2 (en) * | 2009-04-14 | 2016-03-15 | Maui Imaging, Inc. | Calibration of ultrasound probes |
| US9788813B2 (en) | 2010-10-13 | 2017-10-17 | Maui Imaging, Inc. | Multiple aperture probe internal apparatus and cable assemblies |
| JP4498451B2 (ja) * | 2008-09-05 | 2010-07-07 | 株式会社日立メディコ | 超音波診断装置 |
| JP2009050720A (ja) * | 2008-11-04 | 2009-03-12 | Hitachi Medical Corp | 超音波診断装置 |
| JP2010131426A (ja) * | 2010-02-10 | 2010-06-17 | Hitachi Medical Corp | 超音波診断装置 |
| KR102322776B1 (ko) | 2010-02-18 | 2021-11-04 | 마우이 이미징, 인코포레이티드 | 초음파 이미지를 구성하는 방법 및 이를 위한 다중-개구 초음파 이미징 시스템 |
| JP6092109B2 (ja) | 2010-10-13 | 2017-03-08 | マウイ イマギング，インコーポレーテッド | 凹面超音波トランスデューサ及び３ｄアレイ |
| KR20140098843A (ko) | 2011-12-01 | 2014-08-08 | 마우이 이미징, 인코포레이티드 | 핑-기반 및 다수 개구부 도플러 초음파를 이용한 모션 검출 |
| WO2013101988A1 (en) | 2011-12-29 | 2013-07-04 | Maui Imaging, Inc. | M-mode ultrasound imaging of arbitrary paths |
| KR102134763B1 (ko) * | 2012-02-21 | 2020-07-16 | 마우이 이미징, 인코포레이티드 | 다중의 어퍼처 초음파를 사용한 물질 강성의 결정 |
| WO2013148673A1 (en) | 2012-03-26 | 2013-10-03 | Maui Imaging, Inc. | Systems and methods for improving ultrasound image quality by applying weighting factors |
| WO2013153968A1 (ja) * | 2012-04-11 | 2013-10-17 | 株式会社日立製作所 | 超音波診断装置 |
| IN2015DN00556A (ja) | 2012-08-10 | 2015-06-26 | Maui Imaging Inc | |
| WO2014031642A1 (en) | 2012-08-21 | 2014-02-27 | Maui Imaging, Inc. | Ultrasound imaging system memory architecture |
| WO2014038812A1 (en) | 2012-09-06 | 2014-03-13 | Samsung Electronics Co., Ltd. | Method and apparatus for displaying stereoscopic information related to ultrasound sectional plane of target object |
| US9510806B2 (en) | 2013-03-13 | 2016-12-06 | Maui Imaging, Inc. | Alignment of ultrasound transducer arrays and multiple aperture probe assembly |
| JP5603971B2 (ja) * | 2013-05-24 | 2014-10-08 | 株式会社日立メディコ | 超音波装置 |
| US9883848B2 (en) | 2013-09-13 | 2018-02-06 | Maui Imaging, Inc. | Ultrasound imaging using apparent point-source transmit transducer |
| WO2016028787A1 (en) | 2014-08-18 | 2016-02-25 | Maui Imaging, Inc. | Network-based ultrasound imaging system |
| WO2017132517A1 (en) | 2016-01-27 | 2017-08-03 | Maui Imaging, Inc. | Ultrasound imaging with sparse array probes |
-
1990
- 1990-05-11 JP JP12236090A patent/JP2849159B2/ja not_active Expired - Lifetime
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100450446C (zh) * | 2002-10-18 | 2009-01-14 | 株式会社日立医药 | 超声波诊断装置 |
| CN109470883A (zh) * | 2018-10-15 | 2019-03-15 | 同济大学 | 一种gps漂流式流速测量仪 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH0417842A (ja) | 1992-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP2849159B2 (ja) | 超音波診断装置 | |
| EP2816958B1 (en) | Determining material stiffness using multiple aperture ultrasound | |
| EP2504716B1 (en) | Ultrasonic shear wave imaging with focused scanline beamforming | |
| EP0909395B1 (en) | Apparatus for determining movements and velocities of moving objects | |
| EP1273267B9 (en) | Ultrasonic diagnostic apparatus | |
| CN107049361B (zh) | 使用剪切波的声音速度成像 | |
| JP2006325704A (ja) | 超音波診断装置 | |
| Maulik | Spectral Doppler: basic principles and instrumentation | |
| JP2001061840A (ja) | 超音波診断装置 | |
| EP1003051A2 (en) | Ultrasonic imaging apparatus | |
| KR20210014081A (ko) | 환경에서 유체 흐름을 감지하기 위한 초음파 시스템 | |
| CN110893103A (zh) | 用于基于超声的剪切波成像的角度 | |
| JPS62224334A (ja) | 超音波パルスを用いた流速ベクトル分布測定方法及び装置 | |
| JP2005058533A (ja) | 超音波診断装置 | |
| JPH0368694B2 (ja) | ||
| JPH0767451B2 (ja) | 超音波組織変位計測装置 | |
| JP4581596B2 (ja) | 超音波診断装置 | |
| JPH05161647A (ja) | 超音波組織変位計測装置 | |
| JPH0741036B2 (ja) | 超音波装置 | |
| JPH02195948A (ja) | 超音波ドプラ診断装置 | |
| JPH04138146A (ja) | 相関型超音波流速測定装置 | |
| JPH024349A (ja) | 超音波診断装置 | |
| JPH043221B2 (ja) | ||
| JPH04187144A (ja) | 超音波ドプラ診断装置 | |
| JPH0332653A (ja) | 超音波診断装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081106 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081106 Year of fee payment: 10 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091106 Year of fee payment: 11 |
|
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101106 Year of fee payment: 12 |
|
| EXPY | Cancellation because of completion of term | ||
| FPAY | Renewal fee payment (prs date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101106 Year of fee payment: 12 |