JPH031849A - 血液流の速度測定装置 - Google Patents
血液流の速度測定装置Info
- Publication number
- JPH031849A JPH031849A JP2118836A JP11883690A JPH031849A JP H031849 A JPH031849 A JP H031849A JP 2118836 A JP2118836 A JP 2118836A JP 11883690 A JP11883690 A JP 11883690A JP H031849 A JPH031849 A JP H031849A
- Authority
- JP
- Japan
- Prior art keywords
- velocity
- blood flow
- correlation
- measuring
- correlation functions
- 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.)
- Pending
Links
- 230000017531 blood circulation Effects 0.000 title claims description 14
- 238000005314 correlation function Methods 0.000 claims abstract description 24
- 238000005259 measurement Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 10
- 230000004044 response Effects 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 4
- 210000004204 blood vessel Anatomy 0.000 description 4
- 238000001914 filtration Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 238000002604 ultrasonography Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 238000002592 echocardiography Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 241000070918 Cima Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 210000000265 leukocyte Anatomy 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/52—Discriminating between fixed and moving objects or between objects moving at different speeds
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/50—Systems of measurement, based on relative movement of the target
- G01S15/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
- G01S15/582—Velocity or trajectory determination systems; Sense-of-movement determination systems using transmission of interrupted pulse-modulated waves and based upon the Doppler effect resulting from movement of targets
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Hematology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
め要約のデータは記録されません。
Description
【発明の詳細な説明】
（産業上の利用分野）
本発明は、Ｎ個の連続するエコグラフィ信号の系列に基
づく血液流の速度測定装置であって、Ｍ点固定エコー消
去装置を具えており、それが該固定エコー消去装置によ
り供給されるＮ−Ｍ＋ １個の個別信号の相関／加算／
補間により前記速度を測定するためのユニットに追随さ
れている血液流の速度測定装置に関するものである。
づく血液流の速度測定装置であって、Ｍ点固定エコー消
去装置を具えており、それが該固定エコー消去装置によ
り供給されるＮ−Ｍ＋ １個の個別信号の相関／加算／
補間により前記速度を測定するためのユニットに追随さ
れている血液流の速度測定装置に関するものである。
本発明は、血管内の血液流のエコグラフィ試験の一般的
な分野で、特にそのような流れの画像化のために特に有
利に用いられる。
な分野で、特にそのような流れの画像化のために特に有
利に用いられる。
（従来の技術）
前述の種類の測定装置は特開昭６２−１３３９４４号公
報によって既知であり、その装置では前記の測定ユニッ
トが、前記Ｎ−Ｍ＋ １個の個別信号が記憶される記憶
装置と、記憶されたＭ−Ｎ＋１個の信号に対するＭ−Ｎ
個の相関関数を供給する相関装置と、平均相関関数を供
給する加算装置と、及びそれの位置が流れの関連速度に
直接関連する、相関頂点（ｃｏｒ−ｒｅｌａｔｉｏｎ
ｐｅａｋ）とも呼ばれる平均相関関数の最大の決定を可
能にする補間回路とによって形成されている。実際に、
この既知の装置が回帰周期Ｔによる回帰送信の場合には
、移動する目標により反射される連続した超音波信号が
次の式 ８式％（） によってつながれると言う事実を利用していることが思
い出される。これは、信号ｎ＋１が時間偏倚でかなけれ
ば先の信号ｎの複製であることを意味する。後者は、変
換器から目標までの通路を進行し一つの励起から別の励
起まで前記変換器へ戻るために、超音波によって必要な
補充時間を意味する。言い換えれば、 τ＝ ２ＶＴ／Ｃ であり、ここでＶは目標の速度であり、Ｃは音の速度で
ある。時間偏倚τの測定が関連速度■の測定を可能にす
ることが明らかである。
報によって既知であり、その装置では前記の測定ユニッ
トが、前記Ｎ−Ｍ＋ １個の個別信号が記憶される記憶
装置と、記憶されたＭ−Ｎ＋１個の信号に対するＭ−Ｎ
個の相関関数を供給する相関装置と、平均相関関数を供
給する加算装置と、及びそれの位置が流れの関連速度に
直接関連する、相関頂点（ｃｏｒ−ｒｅｌａｔｉｏｎ
ｐｅａｋ）とも呼ばれる平均相関関数の最大の決定を可
能にする補間回路とによって形成されている。実際に、
この既知の装置が回帰周期Ｔによる回帰送信の場合には
、移動する目標により反射される連続した超音波信号が
次の式 ８式％（） によってつながれると言う事実を利用していることが思
い出される。これは、信号ｎ＋１が時間偏倚でかなけれ
ば先の信号ｎの複製であることを意味する。後者は、変
換器から目標までの通路を進行し一つの励起から別の励
起まで前記変換器へ戻るために、超音波によって必要な
補充時間を意味する。言い換えれば、 τ＝ ２ＶＴ／Ｃ であり、ここでＶは目標の速度であり、Ｃは音の速度で
ある。時間偏倚τの測定が関連速度■の測定を可能にす
ることが明らかである。
幅Ｗを有する窓内で、次式
％式％（）
により定義されるＳ、（ｔ）とＳ、＋＋（ｔ）との間の
相関関数が、 Ｃ，、、＋＋（ｔＯ，ｕ）＝ ｃａｎ（ｔｏ、ｕ−ｒ）
を照合する。
相関関数が、 Ｃ，、、＋＋（ｔＯ，ｕ）＝ ｃａｎ（ｔｏ、ｕ−ｒ）
を照合する。
時刻ｔｏはｔｏ＝２ｚ／Ｃとして走査深さ２に関連する
。
。
この間数Ｃ，，、（ｔｏ、ｕ）は自己相関関数であり、
且つそれ故にＵ＝Ｏに対して最大である。
且つそれ故にＵ＝Ｏに対して最大である。
従って、時間偏倚τの測定、及びそれ故速度■の測定は
、関数Ｃ，，、。。＋（ｔＯ，ｕ）が最大であるパラメ
ータＵを探すことにより得ることができる。この最大が
Ｕの値ｕＯに対して得られた場合には、時間偏倚τはｕ
Ｏ＝及びＶ＝ｕＯＣ／２Ｔの等しいことを用いることに
より、それから演鐸される。
、関数Ｃ，，、。。＋（ｔＯ，ｕ）が最大であるパラメ
ータＵを探すことにより得ることができる。この最大が
Ｕの値ｕＯに対して得られた場合には、時間偏倚τはｕ
Ｏ＝及びＶ＝ｕＯＣ／２Ｔの等しいことを用いることに
より、それから演鐸される。
特開昭６２−１３３９４４号公報に記載された装置は、
血液流の速度の測定とそれらの映像化に決定的な貢献を
した。しかしながら、得られる結果の精密度はそれに対
して平均値が形成される相関関数の数、言い換えれば既
知の装置に対するＮ−Ｍに依存する。測定の変動を減少
させるためち、従って、同時に処理されるべき信号の数
を増加させること、あるいは固定エコー消去装置の点の
数Ｍを減少させることが有益である。実際にはこの点の
数を３あるいは４より実質的に小さくすることはほとん
ど不可能である。血液（白血球）の反射率を４０ｄＢも
越え得る血管の壁の高い反射率のせいで、速度輪郭のな
んらかの評定の前には、固定エコー消去装置の使用が不
可欠であることは既知である。最も単純な固定エコー消
去装置は、２点装置とも呼ばれるが、零遅延線と平行な
１回帰周期に匹敵する遅延線から構成される装置である
。それぞれ＋１と−１との重み付は係数がこれらの遅延
線に取り付けられ、それらは重み付けの後に加算装置に
よって加算される。従ってこの既知のフィルターは、原
理的には、固定された組織により生じたエコーの半完全
な減少に導く２つの連続的なエコグラフィ遅延線間の着
を形成する。しかしながら、この技術は低い流速に相当
する信号をも減衰させるという重大な欠点を有する。例
えば、流速の関数としての上述のフィルターの応答は、
５ｋＨｚの回復周波数と５ＭＨｚの送信周波数に対して
、ｖ＝５ｃｍ／ｓに相当する信号が３０ｄＢだけ減衰さ
れるものであることが実演され得る。これが最も低い、
言うならば血管の壁に近い流速を測定することを困難あ
るいは不可能にする。しかしながら、例えば動脈の研究
と臨床診断とについて、これらの速度の知識は非常に重
要である。それ故に、特開昭６４−３１０８１号公報に
記載されたような少なくとも３点を具えるエコー消去装
置が好適に用いられ、その装置が固定組織により生じる
エコーの完全な除去を、しかし低速の流れから生じる信
号を過度に減少させることなく可能にする。
血液流の速度の測定とそれらの映像化に決定的な貢献を
した。しかしながら、得られる結果の精密度はそれに対
して平均値が形成される相関関数の数、言い換えれば既
知の装置に対するＮ−Ｍに依存する。測定の変動を減少
させるためち、従って、同時に処理されるべき信号の数
を増加させること、あるいは固定エコー消去装置の点の
数Ｍを減少させることが有益である。実際にはこの点の
数を３あるいは４より実質的に小さくすることはほとん
ど不可能である。血液（白血球）の反射率を４０ｄＢも
越え得る血管の壁の高い反射率のせいで、速度輪郭のな
んらかの評定の前には、固定エコー消去装置の使用が不
可欠であることは既知である。最も単純な固定エコー消
去装置は、２点装置とも呼ばれるが、零遅延線と平行な
１回帰周期に匹敵する遅延線から構成される装置である
。それぞれ＋１と−１との重み付は係数がこれらの遅延
線に取り付けられ、それらは重み付けの後に加算装置に
よって加算される。従ってこの既知のフィルターは、原
理的には、固定された組織により生じたエコーの半完全
な減少に導く２つの連続的なエコグラフィ遅延線間の着
を形成する。しかしながら、この技術は低い流速に相当
する信号をも減衰させるという重大な欠点を有する。例
えば、流速の関数としての上述のフィルターの応答は、
５ｋＨｚの回復周波数と５ＭＨｚの送信周波数に対して
、ｖ＝５ｃｍ／ｓに相当する信号が３０ｄＢだけ減衰さ
れるものであることが実演され得る。これが最も低い、
言うならば血管の壁に近い流速を測定することを困難あ
るいは不可能にする。しかしながら、例えば動脈の研究
と臨床診断とについて、これらの速度の知識は非常に重
要である。それ故に、特開昭６４−３１０８１号公報に
記載されたような少なくとも３点を具えるエコー消去装
置が好適に用いられ、その装置が固定組織により生じる
エコーの完全な除去を、しかし低速の流れから生じる信
号を過度に減少させることなく可能にする。
測定精度を改善する唯一の方法は、速度の評定を得るた
めに用いられる連続した信号の数Ｎを増加させることで
あると結論することができる。しかしながら、この解決
法はそれが測定速度を減少させ、あるいは表示の場合に
１／ＮＴで与えられる映像化速度を減少させるという直
接の影響を有することは明らかである。
めに用いられる連続した信号の数Ｎを増加させることで
あると結論することができる。しかしながら、この解決
法はそれが測定速度を減少させ、あるいは表示の場合に
１／ＮＴで与えられる映像化速度を減少させるという直
接の影響を有することは明らかである。
（発明が解決しようとする課題）
本発明によって解決されるべき技術的な問題点は、速度
に影響を与えることなくより高い測定精度を達成するこ
とを可能にし、あるいは同じことではあるが、精度を低
下させることなく測定速度を増大することを可能にする
、前述の種類の測定装置を実現することである。
に影響を与えることなくより高い測定精度を達成するこ
とを可能にし、あるいは同じことではあるが、精度を低
下させることなく測定速度を増大することを可能にする
、前述の種類の測定装置を実現することである。
（課題を解決するための手段）
本発明に従って提出された技術的な問題点に対する解決
法は、前記測定ユニットが一方ではＮＦ個の並列処理チ
ャンネルを備え、その各チャンネルはフィルターと、Ｎ
−Ｍ＋１個の濾波された信号を記憶するための記憶装置
及びＮ−Ｍ個の相関関数を供給する相関装置により形成
されており、他方では得られたＮＦ（Ｎ−Ｍ）個の相関
関数の平均値を形成する加算装置を具えており、関連す
る速度の評定は平均相関関数に作用する補間回路により
供給される。本発明が基礎としている本質的な技術的効
果は、それ故に基礎として平均相関関数が計算される個
別の相関関数の数である係数ＮＦによる乗法に存する。
法は、前記測定ユニットが一方ではＮＦ個の並列処理チ
ャンネルを備え、その各チャンネルはフィルターと、Ｎ
−Ｍ＋１個の濾波された信号を記憶するための記憶装置
及びＮ−Ｍ個の相関関数を供給する相関装置により形成
されており、他方では得られたＮＦ（Ｎ−Ｍ）個の相関
関数の平均値を形成する加算装置を具えており、関連す
る速度の評定は平均相関関数に作用する補間回路により
供給される。本発明が基礎としている本質的な技術的効
果は、それ故に基礎として平均相関関数が計算される個
別の相関関数の数である係数ＮＦによる乗法に存する。
これは濾波動作が信号間の時間関係を維持し、且つ従っ
て同じ相関結果に導くという事実によって可能にされる
。種々の並列処理チャンネルにより得られた相関関数の
独立性は、Ｎｐ個のフィルター間の相関の不在によって
保証される。このようにして、同じ測定周波数１／ＮＴ
に対して精度が係数四によって改善された速度評定が得
られる。逆に、１／ＮＯＴが本発明による装置と同じ精
度を生じる前述の既知の装置の測定速度である場合には
、ＮとＮ。どの間の関係は ＮＦ（Ｎ−Ｍ）＝ＮＦＭ によって与えられ、それで Ｎ＝（Ｎｏ−Ｍ）／ＮＰ＋Ｍ であり、ここで、Ｎ０＝１５． Ｍ＝３で且つＮＦ＝６
であり、それでＮ＝５であり、これが測定速度に関して
は係数３の利得に相当する。しかしながら、並列処理チ
ャンネルの数ＮＦが増加した場合には、限界値がＭに等
しいＮに上昇する。
て同じ相関結果に導くという事実によって可能にされる
。種々の並列処理チャンネルにより得られた相関関数の
独立性は、Ｎｐ個のフィルター間の相関の不在によって
保証される。このようにして、同じ測定周波数１／ＮＴ
に対して精度が係数四によって改善された速度評定が得
られる。逆に、１／ＮＯＴが本発明による装置と同じ精
度を生じる前述の既知の装置の測定速度である場合には
、ＮとＮ。どの間の関係は ＮＦ（Ｎ−Ｍ）＝ＮＦＭ によって与えられ、それで Ｎ＝（Ｎｏ−Ｍ）／ＮＰ＋Ｍ であり、ここで、Ｎ０＝１５． Ｍ＝３で且つＮＦ＝６
であり、それでＮ＝５であり、これが測定速度に関して
は係数３の利得に相当する。しかしながら、並列処理チ
ャンネルの数ＮＦが増加した場合には、限界値がＭに等
しいＮに上昇する。
Ｎ２個のフィルターの選択は、信号分解能の退化を敢え
てすることなく、信号変換器への異なる濾波動作を適用
することにより、Ｎ１個の異なる圧電変換器を実現する
思想に基づいている。それ故に、Ｆｅがこの装置のサン
プリング周波数である本発明に従った特別の実施例では
、各フィルターが１に等しい係数応答とランダム位相と
により、及び周波数零に対する応答零により区間 ］
０． Ｆｅ／２］上に定義される。ＮＦ個のフィルター
の各々が従ってそのランダム位相分布によって特徴付け
られる。
てすることなく、信号変換器への異なる濾波動作を適用
することにより、Ｎ１個の異なる圧電変換器を実現する
思想に基づいている。それ故に、Ｆｅがこの装置のサン
プリング周波数である本発明に従った特別の実施例では
、各フィルターが１に等しい係数応答とランダム位相と
により、及び周波数零に対する応答零により区間 ］
０． Ｆｅ／２］上に定義される。ＮＦ個のフィルター
の各々が従ってそのランダム位相分布によって特徴付け
られる。
この種類のフィルターは信号の周波数帯域を維持するこ
とも注意されるべきである。
とも注意されるべきである。
最後に、既知の方法では“１ビツト”相関と呼ばれる相
関を使用する方法も用いられ得て、その方法では各記憶
装置は相当するフィルタ・−により供給されたＮ−Ｍ＋
１個の信号の符号のみを記憶する。
関を使用する方法も用いられ得て、その方法では各記憶
装置は相当するフィルタ・−により供給されたＮ−Ｍ＋
１個の信号の符号のみを記憶する。
この場合には相関関数の頂点が二等辺三角形として形成
される。この形状の認識が、最高点とその２つの隣接点
とから出発する、線型補間による相関頂点の完全な再建
を可能にし、従ってそれの最大の位置ｕＯの正確な決定
を可能にする。
される。この形状の認識が、最高点とその２つの隣接点
とから出発する、線型補間による相関頂点の完全な再建
を可能にし、従ってそれの最大の位置ｕＯの正確な決定
を可能にする。
（実施例）
以下、添付の図面を参照しつつ例を用いて本発明の詳細
な説明する。
な説明する。
第１図は０９月の歩みで１からＮまで変化するＮ個の連
続したエコグラフィ信号Ｓ、（ｔ）の系列に基づいた血
液流の速度の測定装置を略図的に示している。これらの
エコグラフィ信号は送信／受信ユニット１００（図示せ
ず）から発生し、その送信／受信ユニットは、普通の方
法で、送信段からそれが受信する電気的励起信号を、回
帰周期Ｔを有する超音波パルスの周期的な列に変換する
圧電変換器を具えている。受信段はサンプリング周波数
Ｆｅでサンプリングされたエコグラフィ信号Ｓゎ（１）
を出力し、その信号が研究されている媒体を通って圧電
変換器へ戻される。第１図の速度測定装置は、Ｎ個のエ
コグラフィ信号Ｓ、（ｔ）系列から、血管の固定壁から
生じる高振幅成分を含まないＮ−Ｍ＋ １個の個別の信
号ａｔ（ｔ）の系列を製作する、Ｍ点固定エコー消去装
置２００を具えている。先に示したように、Ｍ＝３点を
具え且つ本発明による装置に適する固定エコー消去装置
の例は、特開昭６４−３１０８１号公報に記載されてい
る。
続したエコグラフィ信号Ｓ、（ｔ）の系列に基づいた血
液流の速度の測定装置を略図的に示している。これらの
エコグラフィ信号は送信／受信ユニット１００（図示せ
ず）から発生し、その送信／受信ユニットは、普通の方
法で、送信段からそれが受信する電気的励起信号を、回
帰周期Ｔを有する超音波パルスの周期的な列に変換する
圧電変換器を具えている。受信段はサンプリング周波数
Ｆｅでサンプリングされたエコグラフィ信号Ｓゎ（１）
を出力し、その信号が研究されている媒体を通って圧電
変換器へ戻される。第１図の速度測定装置は、Ｎ個のエ
コグラフィ信号Ｓ、（ｔ）系列から、血管の固定壁から
生じる高振幅成分を含まないＮ−Ｍ＋ １個の個別の信
号ａｔ（ｔ）の系列を製作する、Ｍ点固定エコー消去装
置２００を具えている。先に示したように、Ｍ＝３点を
具え且つ本発明による装置に適する固定エコー消去装置
の例は、特開昭６４−３１０８１号公報に記載されてい
る。
第１図から明らかなように、固定エコー除去装置２００
は、特開昭６２−１３３９４４に従った固定エコー消去
装置により供給されるＮ−Ｍ＋１個の個別の信号ａｔ（
ｔ）の相関／加算／補間によって、関連する血液流の速
度■を測定するためのユニット３００に引き継がれる。
は、特開昭６２−１３３９４４に従った固定エコー消去
装置により供給されるＮ−Ｍ＋１個の個別の信号ａｔ（
ｔ）の相関／加算／補間によって、関連する血液流の速
度■を測定するためのユニット３００に引き継がれる。
第１図から明らかなように、測定ユニット３００は最初
に３１０ｊとして示したＮ１個の並列処理チャンネルを
具えており、そこではｊ＝１．２゜・・・ＮＦである。
に３１０ｊとして示したＮ１個の並列処理チャンネルを
具えており、そこではｊ＝１．２゜・・・ＮＦである。
これらの各処理チャンネルは、周波数ｆの関数としての
それらの応答Ｈ（ｆ）が第２ａ図及び第２ｂ図に略図的
に示しであるフィルター３１１ｊを具えており、係数Ｉ
Ｈ（ｆ） ｌはｆ＝ｏに対しては０であり、周波数Ｏ
から周波数Ｆｅ／２までは１であって、ここでＦｅはサ
ンプリング周波数であり、位相φ（ｆ）は同一区間内で
この周波数のランダム関数φ（ｆ） ＝ Ｒａｎ（ｆ） である。この種類のフィルターの時間応答ｈ（ｔ）の幾
つかの例を第３図に示す。この図に示した４個のフィル
ターはそれらの間の相関を完全に除去された特性を有し
ており、相互に心細の処理装置の完全な独立性を保証し
ていて、この特性が本発明による装置の有効性を保証す
るために絶対に必要である。フィルター３１１ｊにより
実行される動作は、入力信号に対して出力端子上に ｄ；（ｔ） ＝ ａｔ（ｔ）ｘｈ＋（ｔ）に等しい信号
ｄｒ（ｔ）を生じる旋回である。続いて、チャンネル３
１０ｊは、相関装置３１３ｊを用いて、Ｃｆ＋（ｔｏ、
ｕ−ｒ ）＝Ｃｆ、 ＋＋＋（ｔｏ、ｕ）＝ｆｄｉ−＋
（ｔ＋ｕ）ｄｉ（ｔ）ｄｔ によって定義されるＮ−Ｍ個の相関関数の計算に役立つ
Ｎ−Ｍ＋ １個の濾波された信号ｄｉ（Ｄの値を記憶す
るための記憶装置３１２ｊを具えている。相関関数を決
定するために用いられる信号ｄｔ（ｔ）は信号自身の符
号のみによって制限され得る。処理される信号がそこで
１ビツトのみを専有するだけであるから、この方法が相
関計算及び記憶装置３１２ｊのサイズをも単純化する。
それらの応答Ｈ（ｆ）が第２ａ図及び第２ｂ図に略図的
に示しであるフィルター３１１ｊを具えており、係数Ｉ
Ｈ（ｆ） ｌはｆ＝ｏに対しては０であり、周波数Ｏ
から周波数Ｆｅ／２までは１であって、ここでＦｅはサ
ンプリング周波数であり、位相φ（ｆ）は同一区間内で
この周波数のランダム関数φ（ｆ） ＝ Ｒａｎ（ｆ） である。この種類のフィルターの時間応答ｈ（ｔ）の幾
つかの例を第３図に示す。この図に示した４個のフィル
ターはそれらの間の相関を完全に除去された特性を有し
ており、相互に心細の処理装置の完全な独立性を保証し
ていて、この特性が本発明による装置の有効性を保証す
るために絶対に必要である。フィルター３１１ｊにより
実行される動作は、入力信号に対して出力端子上に ｄ；（ｔ） ＝ ａｔ（ｔ）ｘｈ＋（ｔ）に等しい信号
ｄｒ（ｔ）を生じる旋回である。続いて、チャンネル３
１０ｊは、相関装置３１３ｊを用いて、Ｃｆ＋（ｔｏ、
ｕ−ｒ ）＝Ｃｆ、 ＋＋＋（ｔｏ、ｕ）＝ｆｄｉ−＋
（ｔ＋ｕ）ｄｉ（ｔ）ｄｔ によって定義されるＮ−Ｍ個の相関関数の計算に役立つ
Ｎ−Ｍ＋ １個の濾波された信号ｄｉ（Ｄの値を記憶す
るための記憶装置３１２ｊを具えている。相関関数を決
定するために用いられる信号ｄｔ（ｔ）は信号自身の符
号のみによって制限され得る。処理される信号がそこで
１ビツトのみを専有するだけであるから、この方法が相
関計算及び記憶装置３１２ｊのサイズをも単純化する。
旋回は信号間の時間関係を修正しないから、フィルター
３１１ｊによる濾波は相関結果には影響しない。言い換
えれば、ｊには無関係に １Ｑｌ＃ ｆ ｄｆ＋＋（ｔ＋ｕ）ｄｆ（ｔ）ｄｔ＝ｆｄ＋＋
＋（ｔ＋Ｕ）ｄｔ（ｔ）ｄｔｌｏ
１である。
３１１ｊによる濾波は相関結果には影響しない。言い換
えれば、ｊには無関係に １Ｑｌ＃ ｆ ｄｆ＋＋（ｔ＋ｕ）ｄｆ（ｔ）ｄｔ＝ｆｄ＋＋
＋（ｔ＋Ｕ）ｄｔ（ｔ）ｄｔｌｏ
１である。
この特性がＮＦ（Ｎ−Ｍ）個の個別の相関関数を効率的
に得ることを可能にし、その平均値は加算装置３２０に
よって形成される。相関関数のインピーダンスによって
、相関関数 がＮＦによって割り算された分散σ２を提供する。
に得ることを可能にし、その平均値は加算装置３２０に
よって形成される。相関関数のインピーダンスによって
、相関関数 がＮＦによって割り算された分散σ２を提供する。
最後に、関数Ｃ（ｔｏ、 ｕ−τ）は補間装置３３０に
より普通の方法で処理され、その補間装置が相関頂点の
位置ｕＯ＝τ＝２ＶＴ／Ｃを計算することにより、関連
する速度の評定Ｖを供給し、その高さＣＭａｘも速度の
評定値■を確認するための回路３４０に用いられるため
に測定される。
より普通の方法で処理され、その補間装置が相関頂点の
位置ｕＯ＝τ＝２ＶＴ／Ｃを計算することにより、関連
する速度の評定Ｖを供給し、その高さＣＭａｘも速度の
評定値■を確認するための回路３４０に用いられるため
に測定される。
速度測定の確認は不可欠である。何故ならばこれは流れ
の領域の外側では、固定エコー除去装置２００の出力信
号が本質的に雑音であるからである。
の領域の外側では、固定エコー除去装置２００の出力信
号が本質的に雑音であるからである。
従って、本発明による装置によって供給された結果は、
速度零の表示ではなく、この結果を零ではないかどうか
確認することが必要である。従って、２つの比較が実行
される。一方では、固定エコー消去装置の出力信号の局
所エネルギーＥか計算される。
速度零の表示ではなく、この結果を零ではないかどうか
確認することが必要である。従って、２つの比較が実行
される。一方では、固定エコー消去装置の出力信号の局
所エネルギーＥか計算される。
続いて、Ｅがしきい値Ｅ０と比較される。この比較の結
果が正の場合には、相関頂点の最大ＣＭａｘが第２のし
きい値Ｅ１と比較される。速度の評定■はこのＣＩＭａ
ｘが第２のしきい値Ｅ、より大きい場合のみ確認される
。
果が正の場合には、相関頂点の最大ＣＭａｘが第２のし
きい値Ｅ１と比較される。速度の評定■はこのＣＩＭａ
ｘが第２のしきい値Ｅ、より大きい場合のみ確認される
。
従ってこの確認された速度は、記憶と走査変換及び色符
号化のための装置を具えるユニット４００（図示せず）
により、普通の方法で表示するために処理される。
号化のための装置を具えるユニット４００（図示せず）
により、普通の方法で表示するために処理される。
第１図は本発明による測定装置のブロック図を示し、
第２ａ図及び第２ｂ図は第１図の装置に用いられたフィ
ルターの周波数応答と位相とを示し、第３図は第２図に
示したような周波数応答を有するスイルターの時間応答
の幾つかの例を示す。 １００・・・送信／受信ユニット ２００・・・Ｍ点固定エコー消去装置 ３００・・・測定ユニット ３１０ｊ・・・並列処理チャンネル ３１１ｊ・・・フィルター ３１２ｊ・・・記憶装置 ３１３ｊ・・・相関装置 ３２０・・・加算装置 ３３０・・・補間装置 ３４０・・・確認回路 ４００・・・ユニット
ルターの周波数応答と位相とを示し、第３図は第２図に
示したような周波数応答を有するスイルターの時間応答
の幾つかの例を示す。 １００・・・送信／受信ユニット ２００・・・Ｍ点固定エコー消去装置 ３００・・・測定ユニット ３１０ｊ・・・並列処理チャンネル ３１１ｊ・・・フィルター ３１２ｊ・・・記憶装置 ３１３ｊ・・・相関装置 ３２０・・・加算装置 ３３０・・・補間装置 ３４０・・・確認回路 ４００・・・ユニット
Claims (1)
- 【特許請求の範囲】 １、Ｎ個の連続するエコグラフィ信号の系列に基づく血
液流の速度測定装置であって、Ｍ点固定エコー消去装置
（２００）を具えており、それが該固定エコー消去装置
（２００）により供給されるＮ−Ｍ＋１個の個別信号の
相関／加算／補間により前記速度を測定するためのユニ
ット（３００）に追随されている血液流の速度測定装置
において、 前記測定ユニット（３００）が、一方ではＮ＿Ｆ個の並
列処理チャンネル（３１０ｊ）を具え、その各チャンネ
ルはフィルター（３１１ｊ）と、Ｎ−Ｍ＋１個の濾波さ
れた信号を記憶するための記憶装置（３１２ｊ）及びＮ
−Ｍ個の相関関数を供給する相関装置（３１３ｊ）によ
り形成されており、他方では得られたＮ＿Ｆ（Ｎ−Ｍ）
個の相関関数の平均値を形成する加算装置を具えており
、関連する速度の評定は平均相関関数に作用する補間回
路（３３０）により供給されることを特徴とする血液流
の速度測定装置。 ２、Ｆｅがこの装置のサンプリング周波数であって、各
フィルター（３１１ｊ）は１に等しい計数応答とランダ
ム位相とにより、及び周波数零に対する応答零により区
間 ］０，Ｆｅ／２］上に定義されることを特徴とする
請求項１記載の血液流の速度測定装置。 ３、各記憶装置（３１２ｊ）が、相当するフィルター（
３１１ｊ）により供給されるＮ−Ｍ＋１個の信号の符号
のみを記憶することを特徴とする請求項１又は２のうち
いずれか１項記載の血液流の速度測定装置。 ４、前記速度に対して得られた評定を確認するための回
路（３４０）を具えることを特徴とする請求項１〜３の
うちいずれか１項記載の血液流の速度測定装置。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| FR8906289A FR2646918B1 (fr) | 1989-05-12 | 1989-05-12 | Dispositif de mesure de la vitesse d'ecoulements sanguins par echographie ultrasonore a cadence de mesure amelioree |
| FR8906289 | 1989-05-12 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH031849A true JPH031849A (ja) | 1991-01-08 |
Family
ID=9381634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2118836A Pending JPH031849A (ja) | 1989-05-12 | 1990-05-10 | 血液流の速度測定装置 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US5062430A (ja) |
| EP (1) | EP0402968B1 (ja) |
| JP (1) | JPH031849A (ja) |
| CN (1) | CN1047202A (ja) |
| DE (1) | DE69013234T2 (ja) |
| FR (1) | FR2646918B1 (ja) |
| IL (1) | IL94333A (ja) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2662265A1 (fr) * | 1990-05-18 | 1991-11-22 | Philips Electronique Lab | Dispositif eliminateur d'echos fixes pour echographe ultrasonore. |
| US5291892A (en) * | 1991-11-04 | 1994-03-08 | General Electric Company | Ultrasonic flow imaging |
| DE69222401T2 (de) * | 1991-12-11 | 1998-04-02 | Koninkl Philips Electronics Nv | Ultraschallechograph zur Messung hoher Geschwindigkeiten von Blutströmungen |
| GB9321379D0 (en) * | 1993-10-15 | 1993-12-08 | Atomic Energy Authority Uk | Flow measuement |
| US5462059A (en) * | 1994-05-25 | 1995-10-31 | The Regents Of The University Of California | Method for assessing and displaying vascular architecture using ultrasound |
| US5390677A (en) * | 1994-05-31 | 1995-02-21 | The Regents Of The University Of California | Method for assessing and displaying the true three dimensional magnitude of blood velocity |
| NO944736D0 (no) * | 1994-12-07 | 1994-12-07 | Vingmed Sound As | Fremgangsmåte for bestemmelse av blodhastighet |
| US5947901A (en) | 1997-09-09 | 1999-09-07 | Redano; Richard T. | Method for hemodynamic stimulation and monitoring |
| ATE436032T1 (de) * | 2004-08-30 | 2009-07-15 | Koninkl Philips Electronics Nv | Einstellbare verfolgung von strömungsgeschwindigkeiten in doppler- geschwindigkeits-spektren |
Family Cites Families (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3882444A (en) * | 1971-06-18 | 1975-05-06 | Automation Ind Inc | Steerable parametric conversion array |
| FR2562675B1 (fr) * | 1984-04-06 | 1989-10-13 | Cgr Ultrasonic | Procede de levee d'ambiguite de la mesure par effet doppler de la vitesse d'un mobile |
| ATE116741T1 (de) * | 1985-05-09 | 1995-01-15 | Met Flow Sa | Verfahren und gerät zum bestimmen von geschwindigkeiten mittels ultraschall-doppler- echographie im inneren eines sich bewegenden fluids. |
| US4803990A (en) * | 1985-12-03 | 1989-02-14 | U.S. Philips Corporation | Examining moving objects by ultrasound echograpy |
| ES2033799T3 (es) * | 1986-04-21 | 1993-04-01 | North American Philips Corporation | Circuito de correlacion y dispositivo para cartografiar flujo de ultrasonidos que comprende dicho circuito. |
| US4809249A (en) * | 1986-04-21 | 1989-02-28 | North American Philips Corporation | Apparatus for ultrasound flow mapping |
| US4761752A (en) * | 1986-04-21 | 1988-08-02 | North American Philips Corporation | Fractional step correlator |
| FR2604081A1 (fr) * | 1986-09-19 | 1988-03-25 | Labo Electronique Physique | Dispositif d'exploration par echographie ultrasonore d'organes en mouvement et d'ecoulements sanguins |
| FR2617982B1 (fr) * | 1987-07-09 | 1989-10-27 | Labo Electronique Physique | Dispositif d'elimination d'echos fixes pour echographe ultrasonore |
| FR2629997B1 (fr) * | 1988-04-19 | 1990-08-17 | Labo Electronique Physique | Dispositif de mesure par correlation de la vitesse d'organes en mouvement et d'ecoulements sanguins |
-
1989
- 1989-05-12 FR FR8906289A patent/FR2646918B1/fr not_active Expired - Lifetime
-
1990
- 1990-05-04 US US07/518,896 patent/US5062430A/en not_active Expired - Fee Related
- 1990-05-07 EP EP90201145A patent/EP0402968B1/fr not_active Expired - Lifetime
- 1990-05-07 DE DE69013234T patent/DE69013234T2/de not_active Expired - Fee Related
- 1990-05-09 CN CN90102795.2A patent/CN1047202A/zh active Pending
- 1990-05-09 IL IL94333A patent/IL94333A/xx not_active IP Right Cessation
- 1990-05-10 JP JP2118836A patent/JPH031849A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| IL94333A0 (en) | 1991-03-10 |
| DE69013234D1 (de) | 1994-11-17 |
| FR2646918B1 (fr) | 1991-10-11 |
| CN1047202A (zh) | 1990-11-28 |
| DE69013234T2 (de) | 1995-05-11 |
| EP0402968A1 (fr) | 1990-12-19 |
| US5062430A (en) | 1991-11-05 |
| FR2646918A1 (fr) | 1990-11-16 |
| EP0402968B1 (fr) | 1994-10-12 |
| IL94333A (en) | 1993-01-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lai et al. | Interpolation methods for time-delay estimation using cross-correlation method for blood velocity measurement | |
| US4803990A (en) | Examining moving objects by ultrasound echograpy | |
| US5662115A (en) | Method for determining the velocity-time spectrum of blood flow | |
| JP2849827B2 (ja) | パルスドップラ信号のフーリエプロセッシング方法および装置 | |
| US5938611A (en) | Method and apparatus for color flow imaging using coded excitation with single codes | |
| JPS6377437A (ja) | 超音波エコーグラフィによる可動物体検査装置 | |
| JPH0246213B2 (ja) | ||
| WO1998020361A1 (en) | Finite amplitude distortion-based inhomogeneous pulse echo ultrasonic imaging | |
| CN111796259B (zh) | 一种基于lfm互相关模特性逼近的子样本时延估计方法 | |
| Lai et al. | An extended autocorrelation method for estimation of blood velocity | |
| US5579771A (en) | Method of and apparatus for the detection and characterization of a segment of an artery by ultrasonic echography | |
| CN104665877A (zh) | 颈动脉血管局部脉搏波传播速度测量方法 | |
| JPH031849A (ja) | 血液流の速度測定装置 | |
| Jensen | Stationary echo canceling in velocity estimation by time-domain cross-correlation | |
| US4809249A (en) | Apparatus for ultrasound flow mapping | |
| US5226420A (en) | Ultrasonic color flow imaging using autoregressive processing | |
| US5320105A (en) | Ultrasonic echograph for measuring high velocities of blood flows | |
| JP2763155B2 (ja) | 運動する器官や血液の運動速度を超音波エコーグラフで測定する装置 | |
| JP3182419B2 (ja) | 血流の測定及び表示装置 | |
| JP2006520642A (ja) | ４−ｄ超音波画像処理および４−ｄドップラー超音波画像処理のための感度増大 | |
| JP3061430B2 (ja) | 速度測定装置 | |
| JP3281435B2 (ja) | 超音波ドプラ診断装置 | |
| JP2901048B2 (ja) | 画像信号処理方法及び装置 | |
| EP1006369B1 (en) | Method and apparatus for adaptive filtering by counting acoustic sample zeroes in ultrasound imaging | |
| US4761752A (en) | Fractional step correlator |