JP3453415B2 - 超音波エコーグラフィによる動脈の弾力性測定装置及び方法 - Google Patents
超音波エコーグラフィによる動脈の弾力性測定装置及び方法Info
- Publication number
- JP3453415B2 JP3453415B2 JP32005293A JP32005293A JP3453415B2 JP 3453415 B2 JP3453415 B2 JP 3453415B2 JP 32005293 A JP32005293 A JP 32005293A JP 32005293 A JP32005293 A JP 32005293A JP 3453415 B2 JP3453415 B2 JP 3453415B2
- Authority
- JP
- Japan
- Prior art keywords
- artery
- measuring
- stimulation
- instantaneous
- pressure
- 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
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/02007—Evaluating blood vessel condition, e.g. elasticity, compliance
-
- 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/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Vascular Medicine (AREA)
- Hematology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Description
【0001】
【産業上の利用分野】本発明は、血圧Pの影響下で動脈
の生理学的パラメタを測定する装置であって、プロファ
イロメータ<profilometer>モード(Mモード)で使用さ
れる超音波エコーグラフ<ultrasonic echograph>内で、
且つチャネルを受信モードで形成するデバイスを含む送
信/受信手段を具えて、瞬間的な血液の流速Q(t)と、動
脈の半径の瞬間的変動Δr(t)と、1番目の刺激ラインA
に沿っての平均半径r0とを測定し計算する1番目のアセ
ンブリーを有して成る動脈の生理学的パラメタを測定す
る装置及び方法に関する。
の生理学的パラメタを測定する装置であって、プロファ
イロメータ<profilometer>モード(Mモード)で使用さ
れる超音波エコーグラフ<ultrasonic echograph>内で、
且つチャネルを受信モードで形成するデバイスを含む送
信/受信手段を具えて、瞬間的な血液の流速Q(t)と、動
脈の半径の瞬間的変動Δr(t)と、1番目の刺激ラインA
に沿っての平均半径r0とを測定し計算する1番目のアセ
ンブリーを有して成る動脈の生理学的パラメタを測定す
る装置及び方法に関する。
【0002】
【従来の技術】循環系は長い間にわたり非常に良く研究
され記述されて来ている。データは動物又は人間につい
て生体内及び生体外の測定により収集されているが、そ
の時間の大部分は侵害的か外傷性のやり方でなされてい
た。ところが新しい超音波のツールが、非侵害的な生体
内測定を可能とするようになっている、例えばパルス化
ドップラー技術がそれである。
され記述されて来ている。データは動物又は人間につい
て生体内及び生体外の測定により収集されているが、そ
の時間の大部分は侵害的か外傷性のやり方でなされてい
た。ところが新しい超音波のツールが、非侵害的な生体
内測定を可能とするようになっている、例えばパルス化
ドップラー技術がそれである。
【0003】CVI(Colour Velocity Imaging - 色速度画
像化) 技術で用いられる時間相関方法は、動脈の機械的
特性に適合する測定の定義及び測定ツールを可能とす
る。この分野の先行技術の背景は主として循環系を記述
する科学業績及び以前にその研究に使われた手段、例え
ば:Londonで刊行された Edward Arnold編集W.W.Nichol
及びM.F.O'Rourke著の "McDonald's Blood flow in Art
eries"により形成されている。
像化) 技術で用いられる時間相関方法は、動脈の機械的
特性に適合する測定の定義及び測定ツールを可能とす
る。この分野の先行技術の背景は主として循環系を記述
する科学業績及び以前にその研究に使われた手段、例え
ば:Londonで刊行された Edward Arnold編集W.W.Nichol
及びM.F.O'Rourke著の "McDonald's Blood flow in Art
eries"により形成されている。
【0004】この種の装置は特に本出願人によるヨーロ
ッパ特許出願第0 458 384 A1号から既知であり、茲でも
それを参考文献として挙げて置く。このような装置は血
管内の血流のエコーグラフィ試験に用いて好適であり、
特にこれらの血流の診断測定や特性的な生理学的パラメ
タの表示に用いて好適である。
ッパ特許出願第0 458 384 A1号から既知であり、茲でも
それを参考文献として挙げて置く。このような装置は血
管内の血流のエコーグラフィ試験に用いて好適であり、
特にこれらの血流の診断測定や特性的な生理学的パラメ
タの表示に用いて好適である。
【0005】上掲のヨーロッパ特許出願による装置の実
施例は、血流のスピードV(t,z)を時間t及び走査幅zの
関数として超音波エコーグラフィを用いて測定する第1
ユニットを含み、また、スピードV(t,z)の標本を記憶す
るメモリも含み、このスピードV(t,z)の測定は使用され
る超音波とは独立である。この実施例は更に、スピード
V(t,z)の上記標本から瞬間流速Q(t)を計算する第1回路
と上記血流を正反対方向に躍動させる2つの血管壁のシ
フトの半径方向スピードV1(t,z) 及びV2(t,z)を測定す
る第2ユニットとの組合せを含み、固定エコーE2(t,z)
が未だそこから除去されていない局所エネルギーを計算
する回路を有する第2計算回路を含み、これに続いて上
記壁のそれぞれの厚み d1 =z4 −z3 及びd2 =z6 −z5 を計算するためz3,z4,z5,z6 の値を定めるしきい値
E0'のしきい値検出器で形成される回路を含み、各時間
値t0に対し上記壁のそれぞれの平均スピード
施例は、血流のスピードV(t,z)を時間t及び走査幅zの
関数として超音波エコーグラフィを用いて測定する第1
ユニットを含み、また、スピードV(t,z)の標本を記憶す
るメモリも含み、このスピードV(t,z)の測定は使用され
る超音波とは独立である。この実施例は更に、スピード
V(t,z)の上記標本から瞬間流速Q(t)を計算する第1回路
と上記血流を正反対方向に躍動させる2つの血管壁のシ
フトの半径方向スピードV1(t,z) 及びV2(t,z)を測定す
る第2ユニットとの組合せを含み、固定エコーE2(t,z)
が未だそこから除去されていない局所エネルギーを計算
する回路を有する第2計算回路を含み、これに続いて上
記壁のそれぞれの厚み d1 =z4 −z3 及びd2 =z6 −z5 を計算するためz3,z4,z5,z6 の値を定めるしきい値
E0'のしきい値検出器で形成される回路を含み、各時間
値t0に対し上記壁のそれぞれの平均スピード
【数7】
を計算するため
Σd1V1(t,z) 又はΣd2V2(t,z)
を与える加算器とM1及びM2をセグメント〔z3,z4 〕又
は〔z5,z6 〕上の測定標本の数としM1又はM2による除
算器とで形成される2つの回路を含み、各壁のシフトを
計算するため
は〔z5,z6 〕上の測定標本の数としM1又はM2による除
算器とで形成される2つの回路を含み、各壁のシフトを
計算するため
【数8】
を与える加算器と、壁の対称シフトを計算するため減算
器及び2による除算器で形成されて時間標本の形で上記
血管の瞬間的半径変動 Δr(t)=(D2(t)-D1(t))/2 を与える回路とで形成される回路を含み、また曲線Q(t)
及びΔr(t)を時間tの関数として表示する手段を含むこ
とを特徴とする。
器及び2による除算器で形成されて時間標本の形で上記
血管の瞬間的半径変動 Δr(t)=(D2(t)-D1(t))/2 を与える回路とで形成される回路を含み、また曲線Q(t)
及びΔr(t)を時間tの関数として表示する手段を含むこ
とを特徴とする。
【0006】流速Q(t)及び動脈の血管の半径変動Δr(t)
は同時に既知となるから、種々の計算が実行でき、幾つ
かの極めて有益な表現が実現できる。実際に関数Δr(t)
は、血管の膨張が圧力の直接の関数であるという意味に
おいて血管中に行き渡る圧力P(t)の像と考えることがで
き、これら2つの関数の増加及び減少フェーズは心臓の
サイクルの経過と同一である。従って一次近似としては
P(t)とΔr(t)とは比例因数を除いて等しいとすることが
許される。
は同時に既知となるから、種々の計算が実行でき、幾つ
かの極めて有益な表現が実現できる。実際に関数Δr(t)
は、血管の膨張が圧力の直接の関数であるという意味に
おいて血管中に行き渡る圧力P(t)の像と考えることがで
き、これら2つの関数の増加及び減少フェーズは心臓の
サイクルの経過と同一である。従って一次近似としては
P(t)とΔr(t)とは比例因数を除いて等しいとすることが
許される。
【0007】流速及び圧力変動の表現の特に興味のある
組合せモードが時間の関数としてパラメタ化され<param
eterised> それ自身へ逆結合された<retrocoupled>心臓
サイクルの曲線の計算及び表示に在り、それは縦軸には
関数Δr(t)の標本値を、横軸には瞬間的流速関数Q(t)の
標本値をプロットして得られた点により形成されて、積
Q×Δr は表示平面の一部であり、そのとき一次近似に
おいては1乗に対して均一<homogeneous> である。
組合せモードが時間の関数としてパラメタ化され<param
eterised> それ自身へ逆結合された<retrocoupled>心臓
サイクルの曲線の計算及び表示に在り、それは縦軸には
関数Δr(t)の標本値を、横軸には瞬間的流速関数Q(t)の
標本値をプロットして得られた点により形成されて、積
Q×Δr は表示平面の一部であり、そのとき一次近似に
おいては1乗に対して均一<homogeneous> である。
【0008】サイクルの面<surface of the cycle>とサ
イクルにより張られた<subtended>面との測定に基づい
て、分析される血管の軸上の厚みのスライスdxに対し動
脈の効率の計算がこうして顕著に推論され得る。
イクルにより張られた<subtended>面との測定に基づい
て、分析される血管の軸上の厚みのスライスdxに対し動
脈の効率の計算がこうして顕著に推論され得る。
【0009】力比により表される効率計算、及び1サイ
クル中又はサイクルの一部中ではエネルギー比により表
される効率計算は、1サイクル中に分析される血管の一
部で失われるエネルギーを直接反映し、従って血管壁の
粗さや局所収縮を直接反映する:これは狭窄症<stenose
s>の検出に対する貴重なツールである。更にまた、上述
のようなループ中の表現は、動脈の各タイプに対し「署
名」に似た極めて特徴的な形状を示して、そのような形
状又は形状の変則は経験に基づいて放射線医学者により
そのような分析される心臓や動脈又はその病理を示すも
のと翻訳できるであろうことは想像に難くない。
クル中又はサイクルの一部中ではエネルギー比により表
される効率計算は、1サイクル中に分析される血管の一
部で失われるエネルギーを直接反映し、従って血管壁の
粗さや局所収縮を直接反映する:これは狭窄症<stenose
s>の検出に対する貴重なツールである。更にまた、上述
のようなループ中の表現は、動脈の各タイプに対し「署
名」に似た極めて特徴的な形状を示して、そのような形
状又は形状の変則は経験に基づいて放射線医学者により
そのような分析される心臓や動脈又はその病理を示すも
のと翻訳できるであろうことは想像に難くない。
【0010】上述の表現モード及びその他の所与のモー
ドの使用には、関数Q(t)及びΔr(t)の高度に正確な決定
が必要である。
ドの使用には、関数Q(t)及びΔr(t)の高度に正確な決定
が必要である。
【0011】そのように高度な精確さが達成され、超音
波周波数とは独立の血流スピードの測定を上流にむかっ
て顕著に与えることのできる超音波エコーグラフィのう
ちには、本出願人によるヨーロッパ特許出願第 0 225 6
67号に既に記載の時間相関原理に従って動作するものが
あり、それはその血流スピード測定用のユニットは2つ
の引き続くエコーに基づき相関関数をあたえる相関回路
を含み、また、該相関関数値に基づきスピードV(t,z)の
推定値を与える多重化/内挿回路をも含む。
波周波数とは独立の血流スピードの測定を上流にむかっ
て顕著に与えることのできる超音波エコーグラフィのう
ちには、本出願人によるヨーロッパ特許出願第 0 225 6
67号に既に記載の時間相関原理に従って動作するものが
あり、それはその血流スピード測定用のユニットは2つ
の引き続くエコーに基づき相関関数をあたえる相関回路
を含み、また、該相関関数値に基づきスピードV(t,z)の
推定値を与える多重化/内挿回路をも含む。
【0012】上述のような技術の状態の場合には動脈中
に生じる流体力学的現象の極めて特徴的なパラメタを時
間の関数として決定することが望ましいであろう。もっ
と詳しく云えば、動脈の任意の点における弾力性γ及び
圧力Pを決定すること、それは先ず始めはよしんば1心
臓サイクル期間に亙るこれらのパラメタの平均値γ0及
びP0 であっても、もし可能ならその時間変動、すなわ
ち関数γ(t) 及び関数P(t) を決定することが提案され
る。動脈の弾力性は微分関係 γ=dr/dP で定義することができ、順応性<compliance>とも呼ばれ
るものである。それについての説明は、EPFL(Ecole Pol
ytechnique Federale de Lausanne), Presse etinforma
tion 所載の Georges Nicod著 "Measure de la complia
nce arterielle"という文献に、動脈中の圧力変動に関
連して見出される。この文献には圧力に関連して順応性
を判定する方法も記載されている。動脈の圧力は、指の
先端の圧力の差を測定する光血量計<photoplethysmogra
ph> により連続的に測定される。更に、分析点における
動脈の直径は超音波トランスデューサによりパルス化エ
コー(すなわちエコー・フォロワー)モードで判定され
る。心臓サイクルに関して圧力と直径センサーとの間の
適切な非位相化が適用されれば、得られた2つの曲線間
の位相差の修正及び比の形成が可能になる。この方法に
より、圧力に関しては単位圧力当たりの直径の増加をパ
ーセントで表す順応性曲線を求めることができる。従っ
て動脈の任意の点の行動を特徴付けて、それによりその
点における動脈の弾力性を推定することができる。しか
しこの判定は近似的なものに留まる、という訳は血圧の
判定は局所的になされるのではなく、常に同じ場所、例
えば指の先でなされるからである。
に生じる流体力学的現象の極めて特徴的なパラメタを時
間の関数として決定することが望ましいであろう。もっ
と詳しく云えば、動脈の任意の点における弾力性γ及び
圧力Pを決定すること、それは先ず始めはよしんば1心
臓サイクル期間に亙るこれらのパラメタの平均値γ0及
びP0 であっても、もし可能ならその時間変動、すなわ
ち関数γ(t) 及び関数P(t) を決定することが提案され
る。動脈の弾力性は微分関係 γ=dr/dP で定義することができ、順応性<compliance>とも呼ばれ
るものである。それについての説明は、EPFL(Ecole Pol
ytechnique Federale de Lausanne), Presse etinforma
tion 所載の Georges Nicod著 "Measure de la complia
nce arterielle"という文献に、動脈中の圧力変動に関
連して見出される。この文献には圧力に関連して順応性
を判定する方法も記載されている。動脈の圧力は、指の
先端の圧力の差を測定する光血量計<photoplethysmogra
ph> により連続的に測定される。更に、分析点における
動脈の直径は超音波トランスデューサによりパルス化エ
コー(すなわちエコー・フォロワー)モードで判定され
る。心臓サイクルに関して圧力と直径センサーとの間の
適切な非位相化が適用されれば、得られた2つの曲線間
の位相差の修正及び比の形成が可能になる。この方法に
より、圧力に関しては単位圧力当たりの直径の増加をパ
ーセントで表す順応性曲線を求めることができる。従っ
て動脈の任意の点の行動を特徴付けて、それによりその
点における動脈の弾力性を推定することができる。しか
しこの判定は近似的なものに留まる、という訳は血圧の
判定は局所的になされるのではなく、常に同じ場所、例
えば指の先でなされるからである。
【0013】更にもう1つの文献:スイスで刊行された
Journal of Hypertension誌1988年第6号 (付録4)179
-181ページ所載のMOOSER他による"Diameter of periphe
ralarteries in cardic cycle" には、動脈の瞬間的直
径の生体内の高精度測定を非侵害的やり方でいかにして
実行するか、が明らかにされている。
Journal of Hypertension誌1988年第6号 (付録4)179
-181ページ所載のMOOSER他による"Diameter of periphe
ralarteries in cardic cycle" には、動脈の瞬間的直
径の生体内の高精度測定を非侵害的やり方でいかにして
実行するか、が明らかにされている。
【0014】
【発明が解決しようとする課題】本発明は、分析される
動脈の予め定められた領域で完全に実行されるエコーグ
ラフィックな測定に基づき、動脈の弾力性(順応性)を
決定することを目的とする。
動脈の予め定められた領域で完全に実行されるエコーグ
ラフィックな測定に基づき、動脈の弾力性(順応性)を
決定することを目的とする。
【0015】本発明のもう1つの目的は、動脈の所与の
領域で動脈の順応性を更に正確に決定することである。
領域で動脈の順応性を更に正確に決定することである。
【0016】
【課題を解決するための手段】これらの目的を達成し、
先行技術の不都合を軽減若しくは消去するためには、上
に規定した生理学的パラメタを測定する装置が更に、上
記動脈の弾力性γ(t)及びその内部の圧力P(t)を決定す
るために、2番目の刺激ラインBに沿ってのパラメタQ'
(t) 及びΔr'(t) を超音波エコーグラフィによって測定
し計算する2番目のアセンブリーを有し、該2番目の刺
激ラインBは1番目の刺激ラインAから距離eだけ離れ
た位置に所在し、その距離は数十ミリメートルと数ミリ
メートルとの間であり、この2つの刺激ラインA及びB
は動脈の経線平面の1つに沿って動脈を横切って通過
し、距離eは動脈の軸の方向に取られたものであり、ま
た、血流に関する信号及び動脈の壁に関する信号のディ
ジタル標本を記憶する1番目及び2番目の記憶手段が、
該ディジタル信号標本について演算する計算手段と共
に、上記1番目及び2番目のアセンブリーの出力に接続
されていることを特徴とする。
先行技術の不都合を軽減若しくは消去するためには、上
に規定した生理学的パラメタを測定する装置が更に、上
記動脈の弾力性γ(t)及びその内部の圧力P(t)を決定す
るために、2番目の刺激ラインBに沿ってのパラメタQ'
(t) 及びΔr'(t) を超音波エコーグラフィによって測定
し計算する2番目のアセンブリーを有し、該2番目の刺
激ラインBは1番目の刺激ラインAから距離eだけ離れ
た位置に所在し、その距離は数十ミリメートルと数ミリ
メートルとの間であり、この2つの刺激ラインA及びB
は動脈の経線平面の1つに沿って動脈を横切って通過
し、距離eは動脈の軸の方向に取られたものであり、ま
た、血流に関する信号及び動脈の壁に関する信号のディ
ジタル標本を記憶する1番目及び2番目の記憶手段が、
該ディジタル信号標本について演算する計算手段と共
に、上記1番目及び2番目のアセンブリーの出力に接続
されていることを特徴とする。
【0017】Rを血液の流体力学的抵抗とし、また、ρ
が血液の体積とするとき L=L0=ρ/(Πr0 2) とし、ΔP/Δx =RQ0 を動脈に沿った血圧勾配とし、更
にまた可測な ∂r/∂x を距離eに亙る動脈の内径の瞬間的変動とするとき、次
の基礎方程式:
が血液の体積とするとき L=L0=ρ/(Πr0 2) とし、ΔP/Δx =RQ0 を動脈に沿った血圧勾配とし、更
にまた可測な ∂r/∂x を距離eに亙る動脈の内径の瞬間的変動とするとき、次
の基礎方程式:
【数9】
から、順応性γの値を計算することが顕著に可能であ
り、それは最小自乗法を適用するのが好適である。
り、それは最小自乗法を適用するのが好適である。
【0018】最初のステップで1心臓サイクルに亙って
の平均値R0, R0Q0及びγ0 が計算され、その後で継起的
反復により、R(t)及びγ(t) が切片ごとに、従ってこの
サイクル中の2つの関数の数個の標本により定数関数の
形で計算される。
の平均値R0, R0Q0及びγ0 が計算され、その後で継起的
反復により、R(t)及びγ(t) が切片ごとに、従ってこの
サイクル中の2つの関数の数個の標本により定数関数の
形で計算される。
【0019】更に、圧力P(t)が偏微分方程式:
【数10】
を積分することにより計算できる。
【0020】この関数P(t)は加算定数を除いて求められ
るが、次に述べるようにP(t)の正しい値は心臓サイクル
中の平均値P0を定めることにより正確に分かるようにす
ることができる。
るが、次に述べるようにP(t)の正しい値は心臓サイクル
中の平均値P0を定めることにより正確に分かるようにす
ることができる。
【0021】本発明による1番目の好適実施例では該装
置が、時間相関原理に従って動作し且つ2次元モードで
超音波画像を生成することを意図する超音波エコーグラ
フの送信/受信手段を有し、その検出は、距離eだけ離
れた2つの隣接する刺激ラインA及びBに限定され且つ
該2つの刺激ラインの相同<homologous>点を一対として
行われ、またそれと同時に、上記2つの隣接する刺激ラ
インに対し受信される信号を再位相化する<rephasing>
ための、上記計算手段に組み込まれた手段を有すること
を特徴とする。
置が、時間相関原理に従って動作し且つ2次元モードで
超音波画像を生成することを意図する超音波エコーグラ
フの送信/受信手段を有し、その検出は、距離eだけ離
れた2つの隣接する刺激ラインA及びBに限定され且つ
該2つの刺激ラインの相同<homologous>点を一対として
行われ、またそれと同時に、上記2つの隣接する刺激ラ
インに対し受信される信号を再位相化する<rephasing>
ための、上記計算手段に組み込まれた手段を有すること
を特徴とする。
【0022】被検患者のスライス画像を生成するための
超音波エコーグラフィでよく知られている2次元モード
すなわち2Dモードは、探針<probe> の隣接する初等トラ
ンスデューサ<elementary transducers>の送信/受信の
ための電子走査を利用して、ライン方向の刺激によって
互いに相次ぐ画像ラインを形成する。本発明を実行する
ためには、2つの隣接ライン(刺激A及びB)で十分で
ある;これらの2つのラインは2Dモードの画像ラインの
間隔により分離して例えば e0 = 0.45 mm とすることができ、或いはこの間隔の数倍で分離して9
画像ラインのジャンプに対しては例えば e0 = 4.5 mm とすることができる。2つのライン間の空間シフトとリ
ンクする時間シフトは、時間関数∂r/∂x で顕著である
ように所与の標本の測定に対し正確に補償しなければな
らない。この補償は、刺激A及びBにより受信した信号
を再位相化するにつれて明白になる;これは刺激A、B
を分離する正確に既知の周期τから容易に導くことがで
きる。
超音波エコーグラフィでよく知られている2次元モード
すなわち2Dモードは、探針<probe> の隣接する初等トラ
ンスデューサ<elementary transducers>の送信/受信の
ための電子走査を利用して、ライン方向の刺激によって
互いに相次ぐ画像ラインを形成する。本発明を実行する
ためには、2つの隣接ライン(刺激A及びB)で十分で
ある;これらの2つのラインは2Dモードの画像ラインの
間隔により分離して例えば e0 = 0.45 mm とすることができ、或いはこの間隔の数倍で分離して9
画像ラインのジャンプに対しては例えば e0 = 4.5 mm とすることができる。2つのライン間の空間シフトとリ
ンクする時間シフトは、時間関数∂r/∂x で顕著である
ように所与の標本の測定に対し正確に補償しなければな
らない。この補償は、刺激A及びBにより受信した信号
を再位相化するにつれて明白になる;これは刺激A、B
を分離する正確に既知の周期τから容易に導くことがで
きる。
【0023】本発明による2番目の好適実施例では該装
置が、時間相関原理に従って動作し且つプロファイロメ
ータとして使用される超音波エコーグラフによって形成
され、該プロファイロメータ内の送信手段は、通常のプ
ロファイロメータのそれに較べてより強力で僅かに焦点
をぼかした超音波ビームを送出するのに適しており、ま
た該装置は、距離eだけ離れた2つの刺激ラインA及び
Bの信号を並列に且つ2つのチャネル形成デバイスの出
力で同時に供給する二重受信チャネルを形成するための
デバイスを、その中に含むことを特徴とする。
置が、時間相関原理に従って動作し且つプロファイロメ
ータとして使用される超音波エコーグラフによって形成
され、該プロファイロメータ内の送信手段は、通常のプ
ロファイロメータのそれに較べてより強力で僅かに焦点
をぼかした超音波ビームを送出するのに適しており、ま
た該装置は、距離eだけ離れた2つの刺激ラインA及び
Bの信号を並列に且つ2つのチャネル形成デバイスの出
力で同時に供給する二重受信チャネルを形成するための
デバイスを、その中に含むことを特徴とする。
【0024】この形のときは動作が二重Mモードと考え
られるので、2つの同時の刺激A、Bから受信する信号
を位相ロックすることは最早必要でない。
られるので、2つの同時の刺激A、Bから受信する信号
を位相ロックすることは最早必要でない。
【0025】小さい距離eにより顕著に物理的に分離さ
れている2つの関数r(t)に対し所要の精度を達成するた
めに、通常のドップラー・システム上では、その貧弱な
軸解像の故に、時間相関原理に基づくエコーグラフが好
適である。時間相関エコーグラフについては、例えば前
に引用した本出願人によるヨーロッパ特許出願第 0 225
667号に記載されている。
れている2つの関数r(t)に対し所要の精度を達成するた
めに、通常のドップラー・システム上では、その貧弱な
軸解像の故に、時間相関原理に基づくエコーグラフが好
適である。時間相関エコーグラフについては、例えば前
に引用した本出願人によるヨーロッパ特許出願第 0 225
667号に記載されている。
【0026】メモリ中に時間標本の形で記憶されている
刺激A及びBからのエコー信号である周期関数γ(t) 及
びP(t)は、再位相化をするならばその後で、データ処理
計算出力を用いて定めることができる。
刺激A及びBからのエコー信号である周期関数γ(t) 及
びP(t)は、再位相化をするならばその後で、データ処理
計算出力を用いて定めることができる。
【0027】血圧Pの影響下で動脈の弾力性γを測定す
る1番目の方法であって、これを上述の実施例と一致す
る装置により実行するために、経線平面の1つに沿って
動脈を横切って通過する2つの隣接する平行な刺激ライ
ンA及びBに対する瞬間的な血液の流速Q(t)と瞬間的な
半径変動Δr(t)と平均半径r0との値で、上記1番目及び
2番目の記憶手段に記憶されているところの値Q(t),Δ
r(t)及びr0を生成するのに適する超音波エコーグラフを
用いる動脈の弾力性γを測定する方法は、次の諸ステッ
プ、すなわち i - 2つの関数Q(t),Q'(t)のうちの1つに基づいて、平
均値Q0及びその関数の導関数∂Q(t)/∂t を計算するス
テップ、 ii - 刺激ラインA及びBに対する瞬間的な値Δr(t),
Δr'(t) に基づいて、瞬間値 δr(t)=Δr(t)−Δr'(t) を計算するステップ、及び iii - 弾力性パラメタγと流体力学的抵抗パラメタRと
圧力勾配ΔP/Δx とのそれぞれ平均値、すなわちγ0,R0
及びR0Q0を、1心臓サイクルに亙って、
る1番目の方法であって、これを上述の実施例と一致す
る装置により実行するために、経線平面の1つに沿って
動脈を横切って通過する2つの隣接する平行な刺激ライ
ンA及びBに対する瞬間的な血液の流速Q(t)と瞬間的な
半径変動Δr(t)と平均半径r0との値で、上記1番目及び
2番目の記憶手段に記憶されているところの値Q(t),Δ
r(t)及びr0を生成するのに適する超音波エコーグラフを
用いる動脈の弾力性γを測定する方法は、次の諸ステッ
プ、すなわち i - 2つの関数Q(t),Q'(t)のうちの1つに基づいて、平
均値Q0及びその関数の導関数∂Q(t)/∂t を計算するス
テップ、 ii - 刺激ラインA及びBに対する瞬間的な値Δr(t),
Δr'(t) に基づいて、瞬間値 δr(t)=Δr(t)−Δr'(t) を計算するステップ、及び iii - 弾力性パラメタγと流体力学的抵抗パラメタRと
圧力勾配ΔP/Δx とのそれぞれ平均値、すなわちγ0,R0
及びR0Q0を、1心臓サイクルに亙って、
【数11】
のように置いたときの関係式:
【数12】
に適用される最小自乗法に基づいて判定するステップを
有する。
有する。
【0028】上述の1番目の方法から演繹される動脈の
弾力性γを測定する2番目の方法は、関数γ(t) 及び関
数R(t)を、少なくとも1心臓サイクルに亙ってその心臓
サイクルの切片内で判定するために、最小自乗法が再び
適用され、それは心臓サイクルを幾つかに分割したその
コンプリメンタリな一部分での判定を継起的に反復し、
最初の反復では1心臓サイクルを2分割し、以下反復毎
に順次4分割、8分割、…、というふうに分割数を次々
と増加させることによるものであり、また、該心臓サイ
クルに亙って幾つかの標本の標本化された関数をγ及び
Rに対して供給する関係式
弾力性γを測定する2番目の方法は、関数γ(t) 及び関
数R(t)を、少なくとも1心臓サイクルに亙ってその心臓
サイクルの切片内で判定するために、最小自乗法が再び
適用され、それは心臓サイクルを幾つかに分割したその
コンプリメンタリな一部分での判定を継起的に反復し、
最初の反復では1心臓サイクルを2分割し、以下反復毎
に順次4分割、8分割、…、というふうに分割数を次々
と増加させることによるものであり、また、該心臓サイ
クルに亙って幾つかの標本の標本化された関数をγ及び
Rに対して供給する関係式
【数13】
によるものであることを特徴とする。
【0029】この2番目の方法に対して実行される計算
は1番目の方法に対するものより更に複雑で精巧であ
る。その理由は、1心臓サイクルに亙って平均すれば0
となる所与の時間の関数が、サイクルの一部分での平均
は最早0にならない、従って最小自乗法を含む計算を考
慮しなければならないからである。
は1番目の方法に対するものより更に複雑で精巧であ
る。その理由は、1心臓サイクルに亙って平均すれば0
となる所与の時間の関数が、サイクルの一部分での平均
は最早0にならない、従って最小自乗法を含む計算を考
慮しなければならないからである。
【0030】こうして順応性γ(t) は、測定領域の動脈
内にある圧力P(t)を知ることなく定められる、というこ
とに注意されたい。しかし、P(t)を判定することはそれ
自体でも又関数Q(t)及びγ(t) に関しても極めて望まし
い。
内にある圧力P(t)を知ることなく定められる、というこ
とに注意されたい。しかし、P(t)を判定することはそれ
自体でも又関数Q(t)及びγ(t) に関しても極めて望まし
い。
【0031】本発明のもう1つの目的は、分析される動
脈の予め定められた任意の領域で瞬間的な圧力P(t)を、
該領域の動脈の順応性γ(t) 及び半径r(t)の瞬間的関数
に基づいて決定することである。
脈の予め定められた任意の領域で瞬間的な圧力P(t)を、
該領域の動脈の順応性γ(t) 及び半径r(t)の瞬間的関数
に基づいて決定することである。
【0032】この目的のために、値r0及びΔr(t)に基づ
いて刺激A又はBのうちの1つに対してr(t)の値が定め
られるところの、動脈の内部の圧力 P(t)=P0+p(t) を測定する方法であって、前に掲げた方法と両立し得る
方法においては、圧力の変数部分p(t)は、偏導関数P1:
いて刺激A又はBのうちの1つに対してr(t)の値が定め
られるところの、動脈の内部の圧力 P(t)=P0+p(t) を測定する方法であって、前に掲げた方法と両立し得る
方法においては、圧力の変数部分p(t)は、偏導関数P1:
【数14】
を持つ関数を時間に関して積分して、更にこれに続いて
P1(t) の連続コンポネント PC : p(t)=P1(t) − PC を抑圧することにより決定されること、及び、p(t)に加
算される固定部分P0は、 r=r(t), T=T(t)=r(t)P(t) とし、またT(t)を血液による動脈の外周の物理的な意味
における張力とするときに、次の関係式:
P1(t) の連続コンポネント PC : p(t)=P1(t) − PC を抑圧することにより決定されること、及び、p(t)に加
算される固定部分P0は、 r=r(t), T=T(t)=r(t)P(t) とし、またT(t)を血液による動脈の外周の物理的な意味
における張力とするときに、次の関係式:
【数15】
に最小自乗法を適用し、それによって、記号< >が1
心臓サイクルに亙っての関数の平均値を示すものとすれ
ば、上記から
心臓サイクルに亙っての関数の平均値を示すものとすれ
ば、上記から
【数16】
が導かれ、また Pe を動脈の外の圧力(定数)とすれ
ば、 P0=T0/r0 + Pe であることにより決定されること、を特徴とする。
ば、 P0=T0/r0 + Pe であることにより決定されること、を特徴とする。
【0033】
【実施例】以下、本発明を更によく理解するために図面
により詳しく説明する。
により詳しく説明する。
【0034】図1は動脈2中の血流1の生理学的パラメ
タを測定する装置の概略図である。この装置は圧電トラ
ンスデューサ11を有し、それは例えば複数のエレメント
を含むアレイである。トランスデューサ11に接続される
送信/受信手段12は超音波走査ビームの形成を保証し、
第1測定ユニット21はトランスデューサ11に返されたエ
コーグラフィック信号を処理して、血流1のスピードV
(t,z)の推定値を、時間t及び方向
タを測定する装置の概略図である。この装置は圧電トラ
ンスデューサ11を有し、それは例えば複数のエレメント
を含むアレイである。トランスデューサ11に接続される
送信/受信手段12は超音波走査ビームの形成を保証し、
第1測定ユニット21はトランスデューサ11に返されたエ
コーグラフィック信号を処理して、血流1のスピードV
(t,z)の推定値を、時間t及び方向
【外1】
の1番目の刺激ラインAに対する走査深度
【外2】
の関数として与える。
【0035】送信/受信手段は通常少なくとも送信段1
3、及び送信段13と測定段21とを分離する分離器14を有
する。送信段13は発振器と選定された循環周波数 1/Tr
で信号発生器を制御する周波数分割器とを含み、該発生
器の電気的刺激信号は複数(M個)の初等<elementary>
トランスデューサから成るトランスデューサ11に与えら
れ、トランスデューサ11はその信号を超音波パルス信号
の周期的な列に変換する。分離器14は測定回路が送信信
号により過負荷になるのを防止し、また受信チャネルを
形成するためM個の初等トランスデューサの受信したM
個の信号をデバイス10に供給する。トランスデューサ11
の初等トランスデューサ用の電子走査制御回路(図示さ
れていない)が、トランスデューサ11の置換を伴わずに
すべての刺激ラインの選定を可能にする、それらのライ
ンは扇形の配置か或いは茲で仮定しているように平行に
配置してある。チャネル形成デバイス10は、分離器14か
ら受信した信号に課された特定の遅延配置を経て、トラ
ンスデューサ11で選択された刺激ラインに沿って最も近
い点から最遠隔の点までに逐次焦点を合わせることによ
り、与えられた数だけの点を取得することを可能とす
る。チャネル形成用の回路10は下流の信号処理回路に出
力信号S(t,z)を供給する。
3、及び送信段13と測定段21とを分離する分離器14を有
する。送信段13は発振器と選定された循環周波数 1/Tr
で信号発生器を制御する周波数分割器とを含み、該発生
器の電気的刺激信号は複数(M個)の初等<elementary>
トランスデューサから成るトランスデューサ11に与えら
れ、トランスデューサ11はその信号を超音波パルス信号
の周期的な列に変換する。分離器14は測定回路が送信信
号により過負荷になるのを防止し、また受信チャネルを
形成するためM個の初等トランスデューサの受信したM
個の信号をデバイス10に供給する。トランスデューサ11
の初等トランスデューサ用の電子走査制御回路(図示さ
れていない)が、トランスデューサ11の置換を伴わずに
すべての刺激ラインの選定を可能にする、それらのライ
ンは扇形の配置か或いは茲で仮定しているように平行に
配置してある。チャネル形成デバイス10は、分離器14か
ら受信した信号に課された特定の遅延配置を経て、トラ
ンスデューサ11で選択された刺激ラインに沿って最も近
い点から最遠隔の点までに逐次焦点を合わせることによ
り、与えられた数だけの点を取得することを可能とす
る。チャネル形成用の回路10は下流の信号処理回路に出
力信号S(t,z)を供給する。
【0036】スピードV(t,z)を測定するユニット21は、
上述の特許出願第0 458 384 A1号に詳しく説明されてい
る。該ユニット21は(図示されてはいないが)増幅器
と、分析される動脈の壁から発するエコーの主要部分を
消去する固定エコー消去器とを有する。信号は、上述の
ヨーロッパ特許出願第 0 225 667号に記述された相関方
法に従って逐次処理されるのが好適であり、該相関方法
は相関回路と好適には多重内挿回路とを具えて既知のや
り方で時間相関原理を利用する。
上述の特許出願第0 458 384 A1号に詳しく説明されてい
る。該ユニット21は(図示されてはいないが)増幅器
と、分析される動脈の壁から発するエコーの主要部分を
消去する固定エコー消去器とを有する。信号は、上述の
ヨーロッパ特許出願第 0 225 667号に記述された相関方
法に従って逐次処理されるのが好適であり、該相関方法
は相関回路と好適には多重内挿回路とを具えて既知のや
り方で時間相関原理を利用する。
【0037】通常のドップラー速度計とは反対に、こう
して定められた流速V(t,z)には用いられた超音波周波数
の分散に鋭敏であるという利点があり、従って結果が遙
かに広い用途を持つということが可能になる。スピード
V(t,z)として見出された数値は、後で処理するためにユ
ニット21でメモリに記憶される。
して定められた流速V(t,z)には用いられた超音波周波数
の分散に鋭敏であるという利点があり、従って結果が遙
かに広い用途を持つということが可能になる。スピード
V(t,z)として見出された数値は、後で処理するためにユ
ニット21でメモリに記憶される。
【0038】ユニット21では固定エコー消去器の出力信
号が血流の直径D(t)を決定するのにも利用される。瞬間
的な関数V(t,z)及びD(t)に基づいて、瞬間流速Q(t)を測
定する回路22は、ユニット21の2つの出力信号に基づい
て血流の速度に関する情報をディジタルな時間標本の形
でディジタル標本を記憶する第1手段23に供給する。
号が血流の直径D(t)を決定するのにも利用される。瞬間
的な関数V(t,z)及びD(t)に基づいて、瞬間流速Q(t)を測
定する回路22は、ユニット21の2つの出力信号に基づい
て血流の速度に関する情報をディジタルな時間標本の形
でディジタル標本を記憶する第1手段23に供給する。
【0039】信号D(t)は平均化回路24にも与えられ、該
回路は動脈の平均内径r0(t) を表すディジタル標本を回
路23に供給する。
回路は動脈の平均内径r0(t) を表すディジタル標本を回
路23に供給する。
【0040】信号S(t,z)に基づいて、更に別の生理学的
パラメタすなわち瞬間的な半径変動Δr(t)の高精度測定
を実行することも可能である、ということはヨーロッパ
特許出願第 0 225 667号に記載の通りである。Δr(t)の
測定のために、血流を直径方向に弾ませる動脈の2つの
壁のシフトの半径方向スピードV1(t,z) 及びV2(t,z)用
の第2測定ユニット25が設けられ、またその他の記憶及
び計算用の回路も設けられる。ユニット25では最初の連
鎖が、2番目の増幅器、メモリ及びそれに続いて相関連
鎖により(図には示されていないやり方で)形成され、
それはV(t,z)の測定に対し上述したのと類似である。本
実施例には固定エコー消去器はない、その理由はこの場
合実際には壁から発するエコー信号は処理されて、壁の
動きは血液の方向にほぼ直角に起き、刺激Aは動脈の経
線平面内に斜めに起きるからである。
パラメタすなわち瞬間的な半径変動Δr(t)の高精度測定
を実行することも可能である、ということはヨーロッパ
特許出願第 0 225 667号に記載の通りである。Δr(t)の
測定のために、血流を直径方向に弾ませる動脈の2つの
壁のシフトの半径方向スピードV1(t,z) 及びV2(t,z)用
の第2測定ユニット25が設けられ、またその他の記憶及
び計算用の回路も設けられる。ユニット25では最初の連
鎖が、2番目の増幅器、メモリ及びそれに続いて相関連
鎖により(図には示されていないやり方で)形成され、
それはV(t,z)の測定に対し上述したのと類似である。本
実施例には固定エコー消去器はない、その理由はこの場
合実際には壁から発するエコー信号は処理されて、壁の
動きは血液の方向にほぼ直角に起き、刺激Aは動脈の経
線平面内に斜めに起きるからである。
【0041】ユニット25では動脈のその他の特性パラメ
タも計算され、その結果として瞬間的な半径変動関数Δ
r(t)の計算となる。それは先ず最初に、直径方向に反対
側に位置する壁の厚みd1(t) 及びd2(t) に関するもの
で、その測定には上記2番目の増幅器の出力における、
局所エネルギーを計算する回路と、調整可能なしきい値
を持つ検出器により形成される瞬間的な壁の厚みd1及び
d2を計算する回路とを利用する。壁の平均スピード
タも計算され、その結果として瞬間的な半径変動関数Δ
r(t)の計算となる。それは先ず最初に、直径方向に反対
側に位置する壁の厚みd1(t) 及びd2(t) に関するもの
で、その測定には上記2番目の増幅器の出力における、
局所エネルギーを計算する回路と、調整可能なしきい値
を持つ検出器により形成される瞬間的な壁の厚みd1及び
d2を計算する回路とを利用する。壁の平均スピード
【数17】
は、1つの加算器と1つの除算器とで各々が形成される
回路を用いて定められる。最終的に、平均スピードの時
間積分用の2番目の積分器(加算器)が、各壁の変位:
回路を用いて定められる。最終的に、平均スピードの時
間積分用の2番目の積分器(加算器)が、各壁の変位:
【数18】
を高い精度で求めることを可能とする。
【0042】壁のスピードは極めて遅い、それは 0.5cm
/sのオーダーである。
/sのオーダーである。
【0043】次いで、求めようとする関数Δr(t)すなわ
ち瞬間的な半径変動は、減算器及び2で割り算する除算
器を使って、 Δr(t)=1/2(D2(t)-D1(t)) という関係式に基づいて得られ、この関係式は動脈壁の
対称シフト(膨張又は収縮)を特徴付けるものである。
ディジタル標本の形で得られたこの結果も、それ自身が
やはり第1記憶手段23に与えられる。
ち瞬間的な半径変動は、減算器及び2で割り算する除算
器を使って、 Δr(t)=1/2(D2(t)-D1(t)) という関係式に基づいて得られ、この関係式は動脈壁の
対称シフト(膨張又は収縮)を特徴付けるものである。
ディジタル標本の形で得られたこの結果も、それ自身が
やはり第1記憶手段23に与えられる。
【0044】エレメント21から25までは、刺激Aにより
動脈1内に生成された超音波エコーに基づき、Q(t),r
0(t)及びΔr(t)を測定するための第1アセンブリーを形
成する。
動脈1内に生成された超音波エコーに基づき、Q(t),r
0(t)及びΔr(t)を測定するための第1アセンブリーを形
成する。
【0045】本発明によれば、2番目の刺激ラインB
が、トランスデューサ11から出発して刺激ラインAと同
じ方向
が、トランスデューサ11から出発して刺激ラインAと同
じ方向
【外3】
且つ刺激ラインAに極く接近して形成される。ラインA
とBを動脈の軸方向に分離する距離eは典型的には 0.1
ミリメートルの数倍ないし数ミリメートルである。eの
値の選択は動脈に沿った心臓収縮の伝播波長λに支配さ
れ、λは0.5mのオーダーであり、eはλよりかなり小さ
くなければならないので 100倍から1000倍小さい値を距
離eとして選ぶのが正しいとされる。
とBを動脈の軸方向に分離する距離eは典型的には 0.1
ミリメートルの数倍ないし数ミリメートルである。eの
値の選択は動脈に沿った心臓収縮の伝播波長λに支配さ
れ、λは0.5mのオーダーであり、eはλよりかなり小さ
くなければならないので 100倍から1000倍小さい値を距
離eとして選ぶのが正しいとされる。
【0046】使用されるエコーグラフィ装置に依存し
て、以下、図3及び図4を引用して説明するように、刺
激Aと刺激Bとは同時であってもよいし、或いは、心臓
サイクルの期間 TC に比較して小さい周期τだけ僅かに
間隔を置いてもよい。それによって送信/受信手段12の
2番目の出力信号、もっと精確に云えばチャネル形成デ
バイス10からの2番目の出力信号S'(t,z) が発せられる
ことになり、該2番目の信号は、Q'(t) 及びΔr'(t) を
測定するための第2アセンブリーに与えられ、それは上
述の測定用第1アセンブリーと同様のものである(第1
アセンブリーで既に計算されているr0の測定はこの場合
にはもはや無用である)。ブロック21, 22, 25と全く同
一のブロック31, 32, 35で処理した後、血液の流速Q'
(t) 及び瞬間的半径変動Δr'(t) はディジタル標本の形
で第2記憶手段33に与えられる。
て、以下、図3及び図4を引用して説明するように、刺
激Aと刺激Bとは同時であってもよいし、或いは、心臓
サイクルの期間 TC に比較して小さい周期τだけ僅かに
間隔を置いてもよい。それによって送信/受信手段12の
2番目の出力信号、もっと精確に云えばチャネル形成デ
バイス10からの2番目の出力信号S'(t,z) が発せられる
ことになり、該2番目の信号は、Q'(t) 及びΔr'(t) を
測定するための第2アセンブリーに与えられ、それは上
述の測定用第1アセンブリーと同様のものである(第1
アセンブリーで既に計算されているr0の測定はこの場合
にはもはや無用である)。ブロック21, 22, 25と全く同
一のブロック31, 32, 35で処理した後、血液の流速Q'
(t) 及び瞬間的半径変動Δr'(t) はディジタル標本の形
で第2記憶手段33に与えられる。
【0047】データ処理を実行する計算手段(計算器)
41は、記憶手段23及び33に記憶されているディジタル標
本に基づいて、分析される動脈を特徴付ける生理学的パ
ラメタの判定を可能とする、その顕著なものは
41は、記憶手段23及び33に記憶されているディジタル標
本に基づいて、分析される動脈を特徴付ける生理学的パ
ラメタの判定を可能とする、その顕著なものは
【数19】
で定義される弾力性γ及び以下に説明する内圧P(t)であ
る。AとBとが非同時刺激である場合には、メモリ23及
び33に記憶されている標本間の再位相化が、計算器41の
実行すべき最初の演算であって、修正されるべき位相の
ずれは 2πτ/TCである(tC は心臓サイクルの周期) と
いうことに留意されたい。
る。AとBとが非同時刺激である場合には、メモリ23及
び33に記憶されている標本間の再位相化が、計算器41の
実行すべき最初の演算であって、修正されるべき位相の
ずれは 2πτ/TCである(tC は心臓サイクルの周期) と
いうことに留意されたい。
【0048】分かり易くするため、本発明による動脈の
エコーグラフィ分析を行う領域を拡大したのが図2であ
る。この領域の境界は、刺激ラインA及びBが動脈の軸
53と交わる交点51及び52を通る2つの半径方向の平面x
及び x+dxであって、平行なラインA,Bが動脈の経線
平面を定義し、それがこの図の断面として選ばれてい
る。これら2つの平面は距離 e0 /sinθ だけ離れている、但し茲で、e0 は2つのラインA,B
間の距離であり、θは動脈の方向と2つの刺激ラインの
方向
エコーグラフィ分析を行う領域を拡大したのが図2であ
る。この領域の境界は、刺激ラインA及びBが動脈の軸
53と交わる交点51及び52を通る2つの半径方向の平面x
及び x+dxであって、平行なラインA,Bが動脈の経線
平面を定義し、それがこの図の断面として選ばれてい
る。これら2つの平面は距離 e0 /sinθ だけ離れている、但し茲で、e0 は2つのラインA,B
間の距離であり、θは動脈の方向と2つの刺激ラインの
方向
【外4】
とのなす角である。e0 として選定された値とθの値と
に依存して、距離eは、0.1ミリメートルの数倍ないし
数ミリメートルの間にある。この図はまた、動脈内で測
定された点血液スピードV(t,z)及びV'(t,z) を示し、該
動脈の通過面すなわち内部断面は S=πr2 (平均断面は So =πr0 2 ) である。
に依存して、距離eは、0.1ミリメートルの数倍ないし
数ミリメートルの間にある。この図はまた、動脈内で測
定された点血液スピードV(t,z)及びV'(t,z) を示し、該
動脈の通過面すなわち内部断面は S=πr2 (平均断面は So =πr0 2 ) である。
【0049】図3及び図4は本発明による2つの特定の
実施例を示す。
実施例を示す。
【0050】図3に示すエコーグラフの構造の本質は、
CVI(色速度画像化) エコーグラフとして知られるエコー
グラフからそれが簡単に演繹できるものだから既知であ
ると考えられるかも知れない。 CVIの適用が係わるのは
128の初等トランスデューサを用いたエコーグラフィ画
像の形成であって、該エコーグラフィ画像は、例えば12
8の共平面平行ラインから成り、その平行ラインの点で
スピードが予め定められたしきい値を正又は負の方向に
超えているものは異なる色によって示されるのである
(静脈及び動脈内の血液で超音波探針に近づき又は離れ
て行く場合)。分析されるラインの数はM個の初等トラ
ンスデューサの数と等しいのが好適である。画像ライン
はその一連番号の順序に従って逐次形成される。しかし
各ラインの継起的な点の測定周波数を測定すべき(血液
の)スピード範囲に適合させるために、2,3又は4本
の継起的画像ラインで同時に進める(続いて次の2,3
又は4本のラインを形成する)のが最適アプローチであ
る。数回の刺激を必要とする2つの隣接するライン(典
型的には0.45mmの間隔で離れている)を形成するために
は、1番目のラインの点が定められ、次いでそれに対応
する2番目のラインの点が期間τの後に定められ、続い
て1番目のラインのその次の点が期間 2τの後に定めら
れ、等々、これら2つのラインのすべての点に及ぶので
ある。この目的のために、送信段63により制御される分
離器64経由でトランスデューサ61からM個のエコー信号
を受信する受信チャネル形成デバイス60は、各時点τで
の遅延を今回は最初のラインへ次回は次のラインへとい
うふうに再配置する。2つの上記ラインに対し得られる
信号は本発明を実行するのに十分である。チャネル形成
デバイス60内で既知のやり方で多重化されたこれらの信
号は、チャネル形成デバイス60の出力ではS(t,z)及びS'
(t+τ,z) と記される。出力信号は分析される一対のラ
インの対応する点では厳密に同時でなければならない。
その結果、大きさが2πτ/TC の再位相化が、信号の処
理の後で、(図1の)メモリ33に含まれるディジタル標
本とメモリ23に含まれるディジタル標本との間で実行さ
れる。
CVI(色速度画像化) エコーグラフとして知られるエコー
グラフからそれが簡単に演繹できるものだから既知であ
ると考えられるかも知れない。 CVIの適用が係わるのは
128の初等トランスデューサを用いたエコーグラフィ画
像の形成であって、該エコーグラフィ画像は、例えば12
8の共平面平行ラインから成り、その平行ラインの点で
スピードが予め定められたしきい値を正又は負の方向に
超えているものは異なる色によって示されるのである
(静脈及び動脈内の血液で超音波探針に近づき又は離れ
て行く場合)。分析されるラインの数はM個の初等トラ
ンスデューサの数と等しいのが好適である。画像ライン
はその一連番号の順序に従って逐次形成される。しかし
各ラインの継起的な点の測定周波数を測定すべき(血液
の)スピード範囲に適合させるために、2,3又は4本
の継起的画像ラインで同時に進める(続いて次の2,3
又は4本のラインを形成する)のが最適アプローチであ
る。数回の刺激を必要とする2つの隣接するライン(典
型的には0.45mmの間隔で離れている)を形成するために
は、1番目のラインの点が定められ、次いでそれに対応
する2番目のラインの点が期間τの後に定められ、続い
て1番目のラインのその次の点が期間 2τの後に定めら
れ、等々、これら2つのラインのすべての点に及ぶので
ある。この目的のために、送信段63により制御される分
離器64経由でトランスデューサ61からM個のエコー信号
を受信する受信チャネル形成デバイス60は、各時点τで
の遅延を今回は最初のラインへ次回は次のラインへとい
うふうに再配置する。2つの上記ラインに対し得られる
信号は本発明を実行するのに十分である。チャネル形成
デバイス60内で既知のやり方で多重化されたこれらの信
号は、チャネル形成デバイス60の出力ではS(t,z)及びS'
(t+τ,z) と記される。出力信号は分析される一対のラ
インの対応する点では厳密に同時でなければならない。
その結果、大きさが2πτ/TC の再位相化が、信号の処
理の後で、(図1の)メモリ33に含まれるディジタル標
本とメモリ23に含まれるディジタル標本との間で実行さ
れる。
【0051】ラインに対して選択的に CVI動作モードを
用いることにより、ライン間ピッチが2倍(0.9mm) で分
離されている2つの隣接ラインを組合せることもでき、
その時には修正すべき再位相化は 4πτ/TC に等しく、
またピッチが3倍(1.35mm)で分離されている2つの隣接
ラインを組合せることもでき、その時には修正すべき再
位相化は 6πτ/TC に等しい。粥腫<atheroma>又は狭窄
症を特徴付けるため、本発明によれば探針の変位を伴わ
ずに幾つかの連続したライン対を分析することができ、
それによって動脈の、軸距離で1ないし3cmに亙って相
互に近く位置する複数個の点で、順応性γ及び圧力Pを
判定することができる。
用いることにより、ライン間ピッチが2倍(0.9mm) で分
離されている2つの隣接ラインを組合せることもでき、
その時には修正すべき再位相化は 4πτ/TC に等しく、
またピッチが3倍(1.35mm)で分離されている2つの隣接
ラインを組合せることもでき、その時には修正すべき再
位相化は 6πτ/TC に等しい。粥腫<atheroma>又は狭窄
症を特徴付けるため、本発明によれば探針の変位を伴わ
ずに幾つかの連続したライン対を分析することができ、
それによって動脈の、軸距離で1ないし3cmに亙って相
互に近く位置する複数個の点で、順応性γ及び圧力Pを
判定することができる。
【0052】図4に示す実施例の動作は図3に示す実施
例のそれよりも単純である。それは本質的にはプロファ
イロメータに係わり、最初は単一の刺激ラインを実現す
ることを考えて、M個の初等トランスデューサを持つト
ランスデューサ71、送信段73、分離器74及び信号S(t,z)
を供給するチャネル形成デバイス70を用いる。このプロ
ファイロメータは、受信モードでは1番目の分析ライン
に隣接しそれに平行に延びる2番目の分析ラインS'(t,
z) の供給に適合する。この目的のために、先ず最初に
必要なことは、僅かに幅の広いビームを送出することで
ある、換言すれば、通常のプロファイロメータの場合よ
りも僅かに焦点をぼかし、僅かに強力でないビームを送
出することである。焦点をぼかすことは、パルス刺激中
にトランスデューサ71における初等トランスデューサの
(図示してない)電子走査の適切な制御により簡単に達
成される。上述の改訂された送信を使い、デバイス70と
同じ信号を受信しデバイス70の配置より僅かにシフトし
た遅延配置を利用する2番目のチャネル形成デバイス7
0' が、受信モードで、2番目の分析ラインの形成を可
能にし、該2番目のラインは1番目のラインに隣接しそ
れに平行に延びて、出力信号S'(t,z) によりそれを識別
できるものである。この実施例では、2つのラインの対
応する点の各1対が同じ瞬間にデバイス70及び70' の出
力に供給されるという意味で、2番目のラインは1番目
のラインと同時に得られる。チャネル形成デバイス70及
び70' 中で適用される遅延配置間の時間シフトは調整可
能とすることができ、この時間シフトから直接演繹され
る距離eの変動として明らかになる。チャネル形成デバ
イスを二重化することは、この場合には2つのラインの
対応する点の時間ロックがもはや必要ないにも拘わら
ず、図4に示すデバイスを図3のそれよりも高価にす
る。
例のそれよりも単純である。それは本質的にはプロファ
イロメータに係わり、最初は単一の刺激ラインを実現す
ることを考えて、M個の初等トランスデューサを持つト
ランスデューサ71、送信段73、分離器74及び信号S(t,z)
を供給するチャネル形成デバイス70を用いる。このプロ
ファイロメータは、受信モードでは1番目の分析ライン
に隣接しそれに平行に延びる2番目の分析ラインS'(t,
z) の供給に適合する。この目的のために、先ず最初に
必要なことは、僅かに幅の広いビームを送出することで
ある、換言すれば、通常のプロファイロメータの場合よ
りも僅かに焦点をぼかし、僅かに強力でないビームを送
出することである。焦点をぼかすことは、パルス刺激中
にトランスデューサ71における初等トランスデューサの
(図示してない)電子走査の適切な制御により簡単に達
成される。上述の改訂された送信を使い、デバイス70と
同じ信号を受信しデバイス70の配置より僅かにシフトし
た遅延配置を利用する2番目のチャネル形成デバイス7
0' が、受信モードで、2番目の分析ラインの形成を可
能にし、該2番目のラインは1番目のラインに隣接しそ
れに平行に延びて、出力信号S'(t,z) によりそれを識別
できるものである。この実施例では、2つのラインの対
応する点の各1対が同じ瞬間にデバイス70及び70' の出
力に供給されるという意味で、2番目のラインは1番目
のラインと同時に得られる。チャネル形成デバイス70及
び70' 中で適用される遅延配置間の時間シフトは調整可
能とすることができ、この時間シフトから直接演繹され
る距離eの変動として明らかになる。チャネル形成デバ
イスを二重化することは、この場合には2つのラインの
対応する点の時間ロックがもはや必要ないにも拘わら
ず、図4に示すデバイスを図3のそれよりも高価にす
る。
【0053】血圧Pの影響下で動脈の弾力性γを、(図
1の)メモリ23及び33に記憶されたディジタル標本に基
づいて測定することは、本発明の目的である。
1の)メモリ23及び33に記憶されたディジタル標本に基
づいて測定することは、本発明の目的である。
【0054】血液の流れに対し循環系に当てはめられた
動力学の基本方程式は: ρ:血液の容積 V=V(t,x):血液のスピード P=P(t,x):血圧 FR :摩擦力 とするとき
動力学の基本方程式は: ρ:血液の容積 V=V(t,x):血液のスピード P=P(t,x):血圧 FR :摩擦力 とするとき
【数20】
と書くことができる。
【0055】Sについて積分し、Sで割り算し、dxで割
り算すると
り算すると
【数21】
が求められる。
【数22】
であるから、式(5) は、
Q=Q(t,x):血液の流速
R:動脈の壁及び血液の粘度による摩擦を考慮した流体
力学的抵抗 L=ρ/S
力学的抵抗 L=ρ/S
【数23】
とするとき
【数24】
と書くことができる。
【0056】計算を続けるために関数
【外5】
が仮定され、血液の場合には定数になるという極めて良
い近似により
い近似により
【数25】
とすれば、簡単化した式(6) は
【数26】
となる。
【0057】この簡単化を利用して、計算を続けるため
には、Q(t)又はQ'(t) のうちから唯1つの関数で十分で
あることに注意されたい。
には、Q(t)又はQ'(t) のうちから唯1つの関数で十分で
あることに注意されたい。
【0058】圧力Pは動脈の壁にも作用しこれらの壁を
膨張させたり収縮させたりするが、それは動脈システム
の行動を記述する重要なパラメタである弾力性すなわち
順応性γにも関係する。動脈の壁に関しては次のように
書くことができる:
膨張させたり収縮させたりするが、それは動脈システム
の行動を記述する重要なパラメタである弾力性すなわち
順応性γにも関係する。動脈の壁に関しては次のように
書くことができる:
【数27】
従って、γ=∂r/∂P に対して、且つ
∂P/∂x =ΔP/Δx
と仮定すれば、式(7) は
【数28】
となる。
【0059】本発明は量r(x)及びr(x+dx) を、その差が
有意であるのに十分な精度で測定する物理的可能性に基
づいている。実行上は、この瞬間的な差の精度は次式: δr(t)=Δr(t)−Δr'(t) に基づいて求めることができる。また、式(8) では次の
恒等性が成り立つ:
有意であるのに十分な精度で測定する物理的可能性に基
づいている。実行上は、この瞬間的な差の精度は次式: δr(t)=Δr(t)−Δr'(t) に基づいて求めることができる。また、式(8) では次の
恒等性が成り立つ:
【数29】
量 dx/dt=e/dtは、動脈に沿った心臓収縮波の (6m/sの
オーダーの) 伝播スピードであって、動脈内の血液のス
ピードではない、従って、0.45mmというeの値に対し
て、dt=75μs という値は値δr の存在を正当化し、こ
れはそのように小さい値dxで測定され得る、ということ
に注意されたい。
オーダーの) 伝播スピードであって、動脈内の血液のス
ピードではない、従って、0.45mmというeの値に対し
て、dt=75μs という値は値δr の存在を正当化し、こ
れはそのように小さい値dxで測定され得る、ということ
に注意されたい。
【0060】式(8) では、量∂Q/∂t が先ず(既知の関
数Q(t)に基づいて)計算され、そして(等式(9) に従っ
て)−∂r/∂x が計算される。
数Q(t)に基づいて)計算され、そして(等式(9) に従っ
て)−∂r/∂x が計算される。
【0061】茲で先ず最初にパラメタγ,R,ΔP/Δx
の平均値を推定することとし、それらをそれぞれ γ0,R0,〔ΔP/Δx 〕0 と書くこととする。ΔP/Δx は関数 dP(x,t)/dx の連続コンポネントを表す、換言すれば平均流速 Q0=<Q(t)> を用いて生成される動脈に沿った平均圧力勾配であっ
て、関数Q(t)に基づいて計算することができ、それはR0
Q0である。更にまた、距離eに亙ってのrの変動は小さ
いから
の平均値を推定することとし、それらをそれぞれ γ0,R0,〔ΔP/Δx 〕0 と書くこととする。ΔP/Δx は関数 dP(x,t)/dx の連続コンポネントを表す、換言すれば平均流速 Q0=<Q(t)> を用いて生成される動脈に沿った平均圧力勾配であっ
て、関数Q(t)に基づいて計算することができ、それはR0
Q0である。更にまた、距離eに亙ってのrの変動は小さ
いから
【数30】
と仮定することができる。
【0062】この推定は最小自乗法を用いるのが好適で
あって、それは
あって、それは
【数31】
とするとき、ノルム:
【数32】
の平方を最小化するという意味である。そこから、偏導
関数を持つ次の3つの方程式、すなわち記号var をそれ
ぞれ var=γ,R,RQ0 とするときの
関数を持つ次の3つの方程式、すなわち記号var をそれ
ぞれ var=γ,R,RQ0 とするときの
【数33】
が出てくる。すべての計算がなされた時に、そこから
(心臓サイクルの期間 TCに亙ってとられた平均値)次
式が求められる:
(心臓サイクルの期間 TCに亙ってとられた平均値)次
式が求められる:
【数34】
順応性は、水銀柱1ミリメートル当たりミクロン (μ/m
mHg)で表すのが好適であり、従ってρ=10-3g/mm3 に対
して γ0 =419.14×103r0 2b μ/mmHg (但し r0,b:in mm,s として)である。γ0 に対して見
出される値は水銀40μに対し1μのオーダーである。上
記最小自乗法は更に次の関係式:
mHg)で表すのが好適であり、従ってρ=10-3g/mm3 に対
して γ0 =419.14×103r0 2b μ/mmHg (但し r0,b:in mm,s として)である。γ0 に対して見
出される値は水銀40μに対し1μのオーダーである。上
記最小自乗法は更に次の関係式:
【数35】
を与え、それに続いて
【数36】
及び
【数37】
そして更に
【数38】
を与える。健康な動脈に対しては、圧力勾配R0Q0は全動
脈に亙って、すなわち心臓から各器官に至るまでほぼ一
定である。しかし狭窄症の領域では動脈の他の点よりも
圧力勾配がかなり高い。上述のような勾配の判定によっ
て、放射線医学者は狭窄症の場所を正確に確認し、また
動脈の異なる部分に沿って数ミリメートル間隔で上述の
測定を数回繰り返すことによりその分布図を作り、圧力
勾配のプロフィールを求めることができるようになる。
またそれは、非侵害過程を用いて上記狭窄症の治療のた
めにはどんなタイプの看護や介入 (血管形成<angioplas
ty>)が適当であるかを考察するのにも有用である。従っ
て動脈瘤も高い精度で検出できる。
脈に亙って、すなわち心臓から各器官に至るまでほぼ一
定である。しかし狭窄症の領域では動脈の他の点よりも
圧力勾配がかなり高い。上述のような勾配の判定によっ
て、放射線医学者は狭窄症の場所を正確に確認し、また
動脈の異なる部分に沿って数ミリメートル間隔で上述の
測定を数回繰り返すことによりその分布図を作り、圧力
勾配のプロフィールを求めることができるようになる。
またそれは、非侵害過程を用いて上記狭窄症の治療のた
めにはどんなタイプの看護や介入 (血管形成<angioplas
ty>)が適当であるかを考察するのにも有用である。従っ
て動脈瘤も高い精度で検出できる。
【0063】ところが、圧力勾配の平均値が既に極めて
正確な徴候を示しているなら、心臓サイクルの中の関数
R(t)及びγ(t) の、時間の関数としての変動を知ること
が望ましくなろう。それは、以前に定めた平均値 L0,γ
0,R0, R0Q0に基づいて再び最小自乗法を適用することに
より達成される、但しそれは1心臓サイクルの全期間に
ついてではなく、まず2つの半サイクルについて、続い
て1サイクルの4分の1、8分の1、…、というふうに
逐次繰り返すのである。この目的のため、出発点は式
(8) とし、但し茲で
正確な徴候を示しているなら、心臓サイクルの中の関数
R(t)及びγ(t) の、時間の関数としての変動を知ること
が望ましくなろう。それは、以前に定めた平均値 L0,γ
0,R0, R0Q0に基づいて再び最小自乗法を適用することに
より達成される、但しそれは1心臓サイクルの全期間に
ついてではなく、まず2つの半サイクルについて、続い
て1サイクルの4分の1、8分の1、…、というふうに
逐次繰り返すのである。この目的のため、出発点は式
(8) とし、但し茲で
【数39】
と仮定すれば、
【数40】
となる。この方法は、セグメントγ及びRでは定数、す
なわち該心臓サイクルの全期間にわたる数個の標本の形
での関数γ(t) 及びR(t)を与える。
なわち該心臓サイクルの全期間にわたる数個の標本の形
での関数γ(t) 及びR(t)を与える。
【0064】この場合には、心臓サイクルの全期間にわ
たる平均は0であったγ0 及びR0計算用の所与の関数
が、心臓サイクルの一部分ではその平均が最早0ではな
いという事実によって、計算が複雑になる、ということ
に注意されたい。
たる平均は0であったγ0 及びR0計算用の所与の関数
が、心臓サイクルの一部分ではその平均が最早0ではな
いという事実によって、計算が複雑になる、ということ
に注意されたい。
【0065】時間の関数としての順応性はその近似関数
γ(t) により既知であり、動脈の内圧は、順応性γの定
義に基づいて時間の関数P(t)として以下の関係式より計
算できる:
γ(t) により既知であり、動脈の内圧は、順応性γの定
義に基づいて時間の関数P(t)として以下の関係式より計
算できる:
【数41】
記号P1(t) は、求めようとしていた圧力P(t)が積分の後
で付加定数すなわち連続コンポネント PC を除いて得ら
れたという事実を表し、 PC については PC =<P1(t) > (1心臓サイクルに亙って計算された平均値)である。
次にもし P(t)=P0+p(t) と仮定すれば、圧力P(t)の変数部分p(t)は差:
で付加定数すなわち連続コンポネント PC を除いて得ら
れたという事実を表し、 PC については PC =<P1(t) > (1心臓サイクルに亙って計算された平均値)である。
次にもし P(t)=P0+p(t) と仮定すれば、圧力P(t)の変数部分p(t)は差:
【数42】
p(t)=P1(t) −<P1(t) > (12)
として求められる。
【0066】固定部分P0は、語の物理的意味において
(医学的意味においてではなく)血液が動脈の周囲に及
ぼす張力T(t)を利用して求められる:
(医学的意味においてではなく)血液が動脈の周囲に及
ぼす張力T(t)を利用して求められる:
【数43】
P(t)−Pe =T(t)/r(t) (壁側間<transmural>圧力) (13)
但し茲でPe は動脈の外側の圧力で、定数であり、その
大きさは大気圧のオーダーである。
大きさは大気圧のオーダーである。
【0067】時間について微分すれば、P=P(t), T=
T(t), r=r(t)と書いて、
T(t), r=r(t)と書いて、
【数44】
が得られる。これにrを掛け算して式(2) を使えば
【数45】
が得られる。
【数46】
とするとき、ノルム:
【数47】
の平方を最小にするために茲で再び最小自乗法を用い
る。すると var=∂T/∂t 及び var=T に対する、偏導関数を持つ2つの方程式:
る。すると var=∂T/∂t 及び var=T に対する、偏導関数を持つ2つの方程式:
【数48】
が得られる。
【0068】一旦すべての計算が完了すれば、そこから
Tの平均値T0が導かれ、この段階では簡単のため γ=γ0 とすれば
Tの平均値T0が導かれ、この段階では簡単のため γ=γ0 とすれば
【数49】
となる。続いて
P0= Pe +T0/r0
そして最終的にこれから
【数50】
導かれる。
【0069】従って(x側の)半径r(t)又は(x+dx側
の) 半径r'(t) に依存して、(図2の)点51又は52にお
ける圧力P(t)が求められる。なお、関係式(15)による平
均値の計算に関する限り、関数∂r/∂t は、関数Δr(t)
又はΔr'(t) から直接微分することにより求められる、
ということに注意されたい。
の) 半径r'(t) に依存して、(図2の)点51又は52にお
ける圧力P(t)が求められる。なお、関係式(15)による平
均値の計算に関する限り、関数∂r/∂t は、関数Δr(t)
又はΔr'(t) から直接微分することにより求められる、
ということに注意されたい。
【図1】図1は、本発明による装置の実施例の概略回路
図である。
図である。
【図2】図2は、動脈の2つの隣接する超音波刺激ライ
ンの部分の径方向の断面図である。
ンの部分の径方向の断面図である。
【図3】図3は、本発明による装置の1番目の実施例の
部分的概略回路図である。
部分的概略回路図である。
【図4】図4は、本発明による装置の2番目の実施例の
部分的概略回路図である。
部分的概略回路図である。
1 血流
2 動脈
10 チャネル形成デバイス
11 トランスデューサ
12 送信/受信手段
13 送信段
14 分離器
21, 31 第1測定ユニット
22, 32 瞬間流速測定回路
23, 33 記憶手段
24, 34 平均化回路
25, 35 第2測定ユニット
41 計算手段(計算器)
─────────────────────────────────────────────────────
フロントページの続き
(72)発明者 オデュ ボンフス
フランス国 94130 ノジャン−シュル
−マルヌ リュ ド ボテ 12 ビス
(56)参考文献 特開 昭60−259250(JP,A)
(58)調査した分野(Int.Cl.7,DB名)
A61B 8/00 - 8/15
Claims (7)
- 【請求項1】 瞬時の血圧P(t)の影響下で動脈
(2)の生理学的パラメタを測定する装置であって、該
装置は、プロファイロメータ・モード(Mモード)で使
用される超音波エコーグラフ内で、且つチャネルを受信
モードで形成するデバイス(10)を含む送信/受信手
段(12)を具えて、動脈内の瞬間的な血液の流速Q(t)
と、動脈の瞬間的な半径変動Δr(t)と、1番目の刺激ラ
インAに沿っての平均半径r0とを測定し計算する1番目
のアセンブリー(22)(25)(24)を有して成る
動脈の生理学的パラメタを測定する装置において、該装
置は更に、上記動脈の弾力性γ(t) 及びその内部の圧力
P(t)を決定するために、2番目の刺激ラインBに沿って
のパラメタ、瞬間的な血液の流速Q'(t) 及び動脈の瞬間
的な半径変動Δr'(t) を超音波エコーグラフィによって
測定し計算する2番目のアセンブリー(32)(35)
を有し、該2番目の刺激ラインBは1番目の刺激ライン
Aから距離eだけ離れた位置に所在し、その距離は動脈
に沿った心臓収縮伝播波長と比べて小さく、この2つの
刺激ラインA及びBは動脈の経線平面の1つに沿って動
脈を横切って通過し、距離eは動脈の軸の方向に取られ
たものであり、また、血流に関する信号及び動脈の壁に
関する信号のディジタル標本を記憶する1番目及び2番
目の記憶手段が、該ディジタル信号標本について演算す
る計算手段(41)と共に、上記1番目及び2番目のア
センブリーの出力に接続されていることを特徴とする動
脈の生理学的パラメタを測定する装置。 - 【請求項2】 請求項1に記載の動脈の生理学的パラメ
タを測定する装置において、 該装置は、時間相関原理に従って動作し且つ2次元モー
ドで超音波画像を生成することを意図する超音波エコー
グラフの送信/受信手段を有し、 その検出は、距離eだけ離れた2つの隣接する刺激ライ
ンA及びBに限定され且つ該2つの刺激ラインの相同点
を一対として行われ、 またそれと同時に、上記2つの隣接する刺激ラインに対
し受信される信号を再位相化するための、上記計算手段
に組み込まれた手段を有することを特徴とする動脈の生
理学的パラメタを測定する装置。 - 【請求項3】 請求項1に記載の動脈の生理学的パラメ
タを測定する装置において、 該装置は、時間相関原理に従って動作し且つプロファイ
ロメータとして使用される超音波エコーグラフによって
形成され、該プロファイロメータ内の送信手段は、通常
のプロファイロメータのそれに較べてより強力で僅かに
焦点をぼかした超音波ビームを送出するのに適してお
り、 また該装置は、距離eだけ離れた2つの刺激ラインA及
びBの信号を並列に且つ2つのチャネル形成デバイスの
出力で同時に供給する二重受信チャネルを形成するため
のデバイスを、その中に含むことを特徴とする動脈の生
理学的パラメタを測定する装置。 - 【請求項4】 血圧Pの影響下で動脈の弾力性γを測定
する方法であって、これを請求項1ないし3のうちのい
ずれか1項に記載の装置により実行するために、経線平
面の1つに沿って動脈を横切って通過する2つの隣接す
る平行な刺激ラインA及びBに対する瞬間的な血液の流
速Q(t)と瞬間的な半径変動Δr(t)と平均半径r0との値
で、上記1番目及び2番目の記憶手段に記憶されている
ところの値Q(t),Δr(t)及びr0を生成するのに適する超
音波エコーグラフを用いる方法において、 該方法は次の諸ステップ、すなわち i - 2つの関数Q(t),Q'(t)のうちの1つに基づいて、平
均値Q0及びその関数の導関数∂Q(t)/∂t を計算するス
テップ、 ii - 刺激ラインA及びBに対する瞬間的な値Δr(t),
Δr'(t) に基づいて、瞬間値 δr(t)=Δr(t)−Δr'(t) を計算するステップ、及び iii - 弾力性パラメタγと流体力学的抵抗パラメタRと
圧力勾配ΔP/Δx とのそれぞれ平均値、すなわちγ0,R0
及びR0Q0を、1心臓サイクルに亙って、 【数1】 のように置いたときの関係式: 【数2】 に適用される最小自乗法に基づいて判定するステップを
有することを特徴とする動脈の弾力性γを測定する方
法。 - 【請求項5】 請求項4に記載の動脈の弾力性γを測定
する方法において、関数γ(t) 及び関数R(t)を、少なく
とも1心臓サイクルに亙ってその心臓サイクルの切片内
で判定するために、最小自乗法が再び適用され、それは
心臓サイクルを幾つかに分割したそのコンプリメンタリ
な一部分での判定を継起的に反復し、最初の反復では1
心臓サイクルを2分割し、以下反復毎に順次4分割、8
分割、…、というふうに分割数を次々と増加させること
によるものであり、また、該心臓サイクルに亙って幾つ
かの標本の標本化された関数をγ及びRに対して供給す
る関係式 【数3】 によるものであることを特徴とする動脈の弾力性γを測
定する方法。 - 【請求項6】 請求項5に記載の方法により、値r0及び
Δr(t)に基づいて刺激A又はBのうちの1つに対してr
(t)の値が定められるところの、動脈の内部の圧力 P(t)=P0+p(t) を測定する方法において、 圧力の変数部分p(t)は、偏導関数P1: 【数4】 を持つ関数を時間に関して積分して、更にこれに続いて
P1(t) の連続コンポネント PC : p(t)=P1(t) − PC を抑圧することにより決定されること、及び、 p(t)に加算される固定部分P0は、 r=r(t), T=T(t)=r(t)P(t) とし、またT(t)を血液による動脈の外周の物理的な意味
における張力とするときに、次の関係式: 【数5】 に最小自乗法を適用し、それによって、記号< >が1
心臓サイクルに亙っての関数の平均値を示すものとすれ
ば、上記から 【数6】 が導かれ、また Pe を動脈の外の圧力(定数)とすれ
ば、 P0=T0/r0 + Pe であることにより決定されること、を特徴とする動脈の
内部の圧力を測定する方法。 - 【請求項7】 請求項1ないし3のうちのいずれか1項
に記載の測定装置を有するか或いは請求項4ないし6の
うちのいずれか1項に記載の方法を実行する超音波エコ
ーグラフ。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR9215506 | 1992-12-22 | ||
FR9215506 | 1992-12-22 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH06237929A JPH06237929A (ja) | 1994-08-30 |
JP3453415B2 true JP3453415B2 (ja) | 2003-10-06 |
Family
ID=9436941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP32005293A Expired - Fee Related JP3453415B2 (ja) | 1992-12-22 | 1993-12-20 | 超音波エコーグラフィによる動脈の弾力性測定装置及び方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US5411028A (ja) |
EP (1) | EP0603967B1 (ja) |
JP (1) | JP3453415B2 (ja) |
DE (1) | DE69329416T2 (ja) |
Families Citing this family (44)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07241288A (ja) * | 1994-03-04 | 1995-09-19 | Hitachi Ltd | 超音波装置 |
EP0674185A1 (fr) * | 1994-03-25 | 1995-09-27 | Laboratoires D'electronique Philips S.A.S. | Procédé et appareil de détection et de caractérisation d'un segment d'artère par échographie ultrasonore |
US6027452A (en) * | 1996-06-26 | 2000-02-22 | Vital Insite, Inc. | Rapid non-invasive blood pressure measuring device |
EP0832604A1 (en) * | 1996-09-30 | 1998-04-01 | Koninklijke Philips Electronics N.V. | Method and device for measuring the elasticity of an artery by ultrasonic echography |
US5919144A (en) * | 1997-05-06 | 1999-07-06 | Active Signal Technologies, Inc. | Apparatus and method for measurement of intracranial pressure with lower frequencies of acoustic signal |
JP3857788B2 (ja) * | 1997-09-01 | 2006-12-13 | テルモ株式会社 | 循環器情報計測システム |
US6135957A (en) * | 1998-01-23 | 2000-10-24 | U.S. Philips Corporation | Method of and apparatus for echographic determination of the viscosity and the pressure gradient in a blood vessel |
EP0931508A1 (fr) * | 1998-01-23 | 1999-07-28 | Koninklijke Philips Electronics N.V. | Procédé et échographe pour la détermination de la viscosité et du gradient de pression dans un vaisseau sanguin |
US6066095A (en) * | 1998-05-13 | 2000-05-23 | Duke University | Ultrasound methods, systems, and computer program products for determining movement of biological tissues |
JP2003517912A (ja) * | 1999-12-21 | 2003-06-03 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 動脈の超音波複合映像シーケンスを表示するための超音波映像処理方法及び検査システム |
US6589189B2 (en) | 2000-01-07 | 2003-07-08 | Rice Creek Medical, Llc | Non-invasive method and apparatus for monitoring intracranial pressure |
US7374538B2 (en) * | 2000-04-05 | 2008-05-20 | Duke University | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
WO2002043564A2 (en) | 2000-11-28 | 2002-06-06 | Allez Physionix Limited | Systems and methods for making non-invasive physiological assessments |
US20060079773A1 (en) * | 2000-11-28 | 2006-04-13 | Allez Physionix Limited | Systems and methods for making non-invasive physiological assessments by detecting induced acoustic emissions |
US7547283B2 (en) * | 2000-11-28 | 2009-06-16 | Physiosonics, Inc. | Methods for determining intracranial pressure non-invasively |
US6817982B2 (en) * | 2002-04-19 | 2004-11-16 | Sonosite, Inc. | Method, apparatus, and product for accurately determining the intima-media thickness of a blood vessel |
US7147605B2 (en) * | 2002-07-08 | 2006-12-12 | Uab Vittamed | Method and apparatus for noninvasive determination of the absolute value of intracranial pressure |
US6835177B2 (en) * | 2002-11-06 | 2004-12-28 | Sonosite, Inc. | Ultrasonic blood vessel measurement apparatus and method |
JP3694019B2 (ja) * | 2003-04-03 | 2005-09-14 | 松下電器産業株式会社 | 超音波診断装置および超音波診断装置の制御方法 |
US7727153B2 (en) * | 2003-04-07 | 2010-06-01 | Sonosite, Inc. | Ultrasonic blood vessel measurement apparatus and method |
JP4263943B2 (ja) * | 2003-05-07 | 2009-05-13 | テルモ株式会社 | 超音波診断装置 |
EP1633234A4 (en) | 2003-06-03 | 2009-05-13 | Physiosonics Inc | SYSTEMS AND METHODS FOR DETERMINING INTRACRANIAL FA NON-INVASIVE PRESSURE AND ACOUSTIC TRANSDUCER ASSEMBLIES FOR USE THEREIN |
US7425199B2 (en) | 2003-12-30 | 2008-09-16 | General Electric Company | Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring |
US7125383B2 (en) * | 2003-12-30 | 2006-10-24 | General Electric Company | Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring |
US9066679B2 (en) | 2004-08-31 | 2015-06-30 | University Of Washington | Ultrasonic technique for assessing wall vibrations in stenosed blood vessels |
JP5529378B2 (ja) * | 2004-08-31 | 2014-06-25 | ユニヴァーシティ オブ ワシントン | 狭窄血管における壁振動を評価するための超音波技法 |
US8062224B2 (en) * | 2004-10-28 | 2011-11-22 | Uab Vittamed | Method and apparatus for non-invasive continuous monitoring of cerebrovascular autoregulation state |
CA2661128C (en) * | 2005-08-18 | 2015-02-10 | Christoffel Leendert De Korte | A method and apparatus for generating hardness and/or strain information of a tissue |
US8109880B1 (en) | 2006-12-26 | 2012-02-07 | Osvaldas Pranevicius | Noninvasive method to measure intracranial and effective cerebral outflow pressure |
US7938780B2 (en) * | 2008-01-31 | 2011-05-10 | Uab Vittamed Technologijos | Ultrasonic method and apparatus for measuring intracranial contents volume change |
RU2372023C1 (ru) * | 2008-05-06 | 2009-11-10 | Геннадий Константинович Пилецкий | Устройство для измерения внутричерепного давления у новорожденных и детей грудного возраста и опорный элемент для этого устройства |
US20100241001A1 (en) * | 2009-03-20 | 2010-09-23 | Palmeri Mark L | Ultrasound Methods, Systems and Computer Program Products for Imaging Fluids |
JP5494958B2 (ja) * | 2010-05-19 | 2014-05-21 | セイコーエプソン株式会社 | 画像形成装置 |
JP5773281B2 (ja) * | 2010-07-14 | 2015-09-02 | 株式会社 東北テクノアーチ | 血管疾患を判定するためのプログラム、媒体および装置 |
WO2013090850A1 (en) * | 2011-12-14 | 2013-06-20 | California Institute Of Technology | Noninvasive systems for blood pressure measurement in arteries |
WO2014084382A1 (ja) | 2012-11-30 | 2014-06-05 | 株式会社 東芝 | 医用画像処理装置 |
WO2015112512A1 (en) | 2014-01-21 | 2015-07-30 | California Institute Of Technology | Portable electronic hemodynamic sensor systems |
TWI562759B (en) * | 2014-05-09 | 2016-12-21 | Prolific Technology Inc | Method for measuring blood vessel elasticity |
CN107635472A (zh) | 2015-06-19 | 2018-01-26 | 神经系统分析公司 | 经颅多普勒探测器 |
US10617388B2 (en) | 2016-01-05 | 2020-04-14 | Neural Analytics, Inc. | Integrated probe structure |
US11589836B2 (en) | 2016-01-05 | 2023-02-28 | Novasignal Corp. | Systems and methods for detecting neurological conditions |
JP2019504670A (ja) | 2016-01-05 | 2019-02-21 | ニューラル アナリティクス、インコーポレイテッド | 臨床指標を決定するためのシステム及び方法 |
EP3278735A1 (en) | 2016-08-03 | 2018-02-07 | PI-Harvest Holding AG | A system and method for non-invasive measurement of pressure inside a body including intravascular blood pressure |
CN109561877A (zh) | 2016-08-03 | 2019-04-02 | 皮埃-哈韦斯特控股有限公司 | 一种用于无创测量包括血管内血压的体内压力的系统及方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS5849137A (ja) * | 1981-09-18 | 1983-03-23 | 株式会社東芝 | 超音波血流測定装置 |
US4630612A (en) * | 1983-05-25 | 1986-12-23 | Aloka Co., Ltd. | Ultrasonic diagnostic apparatus |
GB2156985B (en) * | 1984-04-02 | 1987-06-24 | Teltec Electronic Equip | Apparatus for measuring movable part-structures, eg blood vessels, within a living body |
US4803990A (en) * | 1985-12-03 | 1989-02-14 | U.S. Philips Corporation | Examining moving objects by ultrasound echograpy |
FR2662348A1 (fr) * | 1990-05-22 | 1991-11-29 | Philips Electronique Lab | Dispositif de mesure et de visualisation par echographie ultrasonore de debit d'un ecoulement sanguin et de dilatation du vaisseau associe. |
US5203337A (en) * | 1991-05-08 | 1993-04-20 | Brigham And Women's Hospital, Inc. | Coronary artery imaging system |
DE4227800C2 (de) * | 1991-08-21 | 1996-12-19 | Toshiba Kawasaki Kk | Thrombuslösende Behandlungsvorrichtung |
US5265612A (en) * | 1992-12-21 | 1993-11-30 | Medical Biophysics International | Intracavity ultrasonic device for elasticity imaging |
-
1993
- 1993-12-20 JP JP32005293A patent/JP3453415B2/ja not_active Expired - Fee Related
- 1993-12-21 US US08/173,330 patent/US5411028A/en not_active Expired - Lifetime
- 1993-12-21 EP EP93203600A patent/EP0603967B1/fr not_active Expired - Lifetime
- 1993-12-21 DE DE69329416T patent/DE69329416T2/de not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH06237929A (ja) | 1994-08-30 |
EP0603967A1 (fr) | 1994-06-29 |
EP0603967B1 (fr) | 2000-09-13 |
US5411028A (en) | 1995-05-02 |
DE69329416D1 (de) | 2000-10-19 |
DE69329416T2 (de) | 2001-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3453415B2 (ja) | 超音波エコーグラフィによる動脈の弾力性測定装置及び方法 | |
Kanai et al. | Transcutaneous measurement and spectrum analysis of heart wall vibrations | |
US4754760A (en) | Ultrasonic pulse temperature determination method and apparatus | |
Overbeck et al. | Vector Doppler: Accurate measurement of blood velocity in two dimensions | |
Kanai et al. | Noninvasive evaluation of local myocardial thickening and its color-coded imaging | |
US6464637B1 (en) | Automatic flow angle correction by ultrasonic vector | |
Eik-Nes et al. | Methodology and basic problems related to blood flow studies in the human fetus | |
US7125383B2 (en) | Method and apparatus for ultrasonic continuous, non-invasive blood pressure monitoring | |
Hasegawa et al. | Improving accuracy in estimation of artery-wall displacement by referring to center frequency of RF echo | |
US5183046A (en) | Ultrasonic plethysmograph | |
US20090012399A1 (en) | Ultrasonic diagnostic apparatus | |
US8298143B2 (en) | Ultrasonograph that determines tissue properties utilizing a reference waveform | |
JPS6226051A (ja) | 超音波血流量自動測定装置 | |
Jensen | Directional transverse oscillation vector flow estimation | |
US6293914B1 (en) | Ultrasonic system and method for measurement of fluid flow | |
Hoskins | A comparison of single-and dual-beam methods for maximum velocity estimation | |
US5383464A (en) | Ultrasonic doppler diagnostic system | |
US20030163043A1 (en) | Method and apparatus for high strain rate rejection filtering | |
KR20000038847A (ko) | 초음파영상화시스템에서의 집속방법 | |
Tortoli et al. | Accurate Doppler angle estimation for vector flow measurements | |
Kanai et al. | A new method for measuring small local vibrations in the heart using ultrasound | |
EP3493744A1 (en) | System and method for determining cardiac output | |
US20100312110A1 (en) | Ultrasonograph | |
US9918698B2 (en) | System and method for gradient-based k-space search for shear wave velocity dispersion estimation | |
US20090299179A1 (en) | Method For Detecting Cardiac Transplant Rejection |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
LAPS | Cancellation because of no payment of annual fees |