JP2004041382A - Ultrasonograph - Google Patents
Ultrasonograph Download PDFInfo
- Publication number
- JP2004041382A JP2004041382A JP2002201799A JP2002201799A JP2004041382A JP 2004041382 A JP2004041382 A JP 2004041382A JP 2002201799 A JP2002201799 A JP 2002201799A JP 2002201799 A JP2002201799 A JP 2002201799A JP 2004041382 A JP2004041382 A JP 2004041382A
- Authority
- JP
- Japan
- Prior art keywords
- blood pressure
- local
- time
- blood
- diameter
- 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.)
- Granted
Links
Images
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
【0001】
【発明の属する技術分野】
本発明は超音波診断装置に関し、特に、超音波を利用した局所血圧の測定に関する。
【0002】
【従来の技術】
血管内の特定の局所部位について血圧波形を得るために、カテーテル型の血圧計が用いられる。つまり、カテーテルを血管内に挿入し、その先端部に設けられた圧力センサによって特定部位の局所血圧を測るものである。しかし、この手法では、侵襲的かつ観血的で、患者への負担が大きい。なお、非侵襲的な血圧計測の手法として、とう骨動脈の圧力脈波を圧力センサで測るトノメトリ法が知られている。しかし、絶対血圧が得られないため、上腕部に巻かれるカフ型血圧計で測定された血圧値を用いて、圧力センサの検出値を較正する必要がある。
【0003】
ところで、超音波を利用した非侵襲的血圧測定法が知られている(文献1:菅原 基晃 ほか,血圧波形の非侵襲的計測法の開発,医用電子と生体工学,第21巻,pp429,1983.など参照)。つまり、血圧変化と血管直径変化が近似的に線形関係にあることから、超音波エコートラッキング法を用いて測定部位の血管直径変化を測定し、カフ型血圧計で最高血圧と最低血圧を測定し、最大血管直径を最高血圧で、最小血管直径を最低血圧で較正し、これにより、血管直径変化を血圧波形とみなすものである。
【0004】
【発明が解決しようとする課題】
血圧変化と血管直径変化との関係は、厳密には直線ではなく、非線形特性をもつことが知られている(文献2:Hansen F, et al., Diameter and compliance in the human common carotid artery−variations with age and sex, Ultrasound Med & Bilo 21, pp 1−9,1995.及び文献3:Hayashi K, et al., Mechanical properties of human arteries, Biorheology 17, pp 211−218,1980.など参照)。すなわち、血圧が高くなると動脈が硬くなり、血管直径変化が小さくなる。また、加齢や疾患による動脈硬化もある。それらの要因があるにもかかわらず、血管直径変化と血圧変化の関係を線形とみなすと誤差が大きくなる場合がある。したがって、線形近似で求めた血圧波形を用いて循環動態指標、例えば公知のwaveintensityを演算すると、これにも計測誤差を生じることになる。
【0005】
本発明の目的は、超音波を利用して局所血圧をより正確に測定することにある。
【0006】
【課題を解決するための手段】
本発明は、被検体に対して超音波の送受波を行う送受波手段と、前記超音波の送受波により得られた受信信号に基づいて、前記被検体内の血管における特定の局所部位について、そのサイズの時間変化を表す局所サイズ波形情報を求める血管サイズ波形情報演算手段と、前記被検体の体表面に装着され、最高血圧及び最低血圧を測定する血圧計と、前記血圧計によって測定された最高血圧及び最低血圧に基づいて前記血管サイズ波形情報を換算することにより、前記特定の局所部位について、その局所血圧の時間変化を表す局所血圧波形情報を推定する推定手段と、を含み、前記推定手段は、前記局所サイズ波形情報から前記局所血圧波形情報を推定する非線形関数に対して、前記最高血圧及び最低血圧を与えることにより、前記局所血圧波形情報を求める演算を実行することを特徴とする。
【0007】
上記構成によれば、血管サイズ波形情報から局所血圧波形情報を推定する非線形関数を用いて、局所血圧波形情報が求められる。よって、従来の線形関数を用いた場合よりも、局所血圧の推定精度をより向上することができる。
【0008】
望ましくは、前記局所サイズは血管直径である。望ましくは、前記最高血圧をPsとし、このときの血管直径をDsとし、一方、前記最低血圧をPdとし、このときの血管直径をDdとし、ある時刻における血管直径をDとした場合に、当該時刻における局所血圧Pが、
P = Pd*exp[β(D/Dd−1)]
但し、β= ln(Ps/Pd)/(Ds/Dd−1)
によって求められる。
【0009】
望ましくは、前記局所サイズは血管面積である。望ましくは、前記最高血圧をPsとし、このときの血管面積をSsとし、一方、前記最低血圧をPdとし、このときの血管面積をSdとし、ある時刻における血管面積をSとした場合に、当該時刻における局所血圧Pが、
P= Pd*exp[βs(S/Sd−1)]
但し、βs = ln(Ps/Pd)/(Ss/Sd−1)
によって求められる。
【0010】
局所部位についての血管直径は、例えば、超音波ビーム上において、血管の前壁と後壁についてエコートラッキングを行うことにより、2つのトラッキングポイント(ゲート)間の距離として求めることができる。血管断面(短軸断面)の面積は、Bモード画像から心臓壁をトレースして左室面積を求める公知のKI法やA−SMA法などを用いることができる。また、カラードプラ法の血流表示から血管断面積を求めてもよい。
【0011】
上記によって求められた局所血圧波形情報は、例えばwave intensityの演算で用いられる。もちろん、その情報を数値としてあるいはグラフとして表示するようにしてもよい。
【0012】
【発明の実施の形態】
まず、本発明に係る実施形態の原理について説明する。
【0013】
上腕部に巻き付けられるカフ型血圧計で得られた最高血圧(収縮期圧力)をPsとし、このときの血管直径をDsとする。同様に、最低血圧(拡張期圧力)をPdとし、このときの血管直径をDdとする。
【0014】
最初に、(A):局所血圧−血管直径(局所直径)の関係について検討する。いま、血管における局所部位についての局所血圧と局所直径の非線形関係が2つのパラメータa,bを有する次の指数関数で表されるとする。
【0015】
P = P(D,a,b)
= a*exp(b*D) (1)
上記(1)式から,最高血圧と最低血圧は次式で表される。
【0016】
Ps= a*exp(b*Ds) (2)
Pd= a*exp(b*Dd) (3)
上記(2)式、(3)式の関係から、a,bは次式で表される。
【0017】
b= ln(Ps/Pd)/(Ds−Dd) (4)
a= Pd/exp(b*Dd) (5)
したがって、局所血圧Pと血管直径Dとの関係は次式で表される。
【0018】
P = Pd*exp[b*(D−Dd)]
= Pd*exp[β(D/Dd−1)] (6)
ここで、βは、
β= ln(Ps/Pd)/(Ds/Dd−1) (7)
で定義されるstiffness parameterと呼ばれる血管弾性指標である(文献3参照)。つまり、βが大きいほど、血管が硬くなる。
【0019】
以上の結果は、局所血圧−血管直径の特性を指数関数で近似したとき、その指数がstiffness parameterであることを意味している。
【0020】
次に、上記と同様に、(B):局所血圧−血管面積(短軸断面積)の関係について検討する。カフ型血圧計で得られた最高血圧をPsとし、このときの血管面積をSsとする。また、最低血圧をPdとし、このときの血管面積をSdとする。
【0021】
いま、局所血圧Pと血管短軸断面積Sの非線形関係が2つのパラメータA,Bを有する以下の指数関数で表されるとする。
【0022】
P = P(S,A,B)
= A*exp(B*S) (8)
上記(8)式から、最高血圧と最低血圧は次式で表される。
【0023】
Ps= A*exp(B*Ss) (9)
Pd= A*exp(B*Sd) (10)
上記(9)式、(10)式の関係から、A,Bは次式で表される。
【0024】
B= ln(Ps/Pd)/(Ss−Sd) (11)
A= Pd/exp(B*Sd) (12)
したがって、局所血圧Pと血管面積Sとの関係は次式で表される。
【0025】
P= Pd*exp[B*(S−Sd)]
= Pd*exp[βs(S/Sd−1)] (13)
ここで、stiffness parameterは、以下の通りである。
【0026】
βs = ln(Ps/Pd)/(Ss/Sd−1) (14)
血管の断面を円形とすると、その面積Sは血管の直径Dを用いて
S = πD2/4
で表される。また、最大血管直径Dsと血管直径Dは、最小血管直径Ddとそれからの変化分ΔDs,ΔDの和として、それぞれ次式で表される。
【0027】
Ds= Dd+ΔDs (15)
D = Dd+ΔD (16)
血管直径変化ΔDs,ΔDは、最小血管直径Ddに比べて十分小さいので、つまり、ΔDs/Dd<<1, ΔD/Dd<<1なので、
S/Sd−1 = (D/Dd)2−1 ≒2(D/Dd−1) (17)
Ss/Sd−1 =( Ds/Dd)2−1≒2(Ds/Dd−1) (18)
が成り立つ。したがって、
βs≒β/2 (19)
P = Pd*exp[βs(S/Sd−1)] (20)
≒ Pd*exp[β(D/Dd−1)] (21)
と近似できる。上記(20)式、(21)式から、局所血圧−血管面積の関係を指数関数で近似した場合と、血圧−血管直径の関係を指数関数で近似した場合とではほぼ同じ結果が得られることが理解される。
【0028】
次に、図面に基づいて本発明の実施形態について説明する。
【0029】
図1には、本発明に係る超音波診断装置の好適な実施形態が示されており、図1はその全体構成を示すブロック図である。
【0030】
図1において、プローブ10は、図1に示す実施形態において、体表面12上に当接して用いられ、このプローブ10によって超音波が送受波される。具体的には、プローブ10内には複数の振動素子からなるアレイ振動子が設けられており、そのアレイ振動子によって超音波ビームBが形成される。その超音波ビームBを電子走査することにより、走査面Sが形成される。図1に示す例では、超音波ビームBが電子リニア走査されているが、電子走査方式としては、その他に電子セクタ走査などをあげることができる。
【0031】
生体内には、血管14が存在しており、プローブ10から見て血管14は前壁14A及び後壁14Bを有する。血管14における特定の局所部位について計測を行う場合には、例えばBモード画像(二次元断層画像)を観察しながら、プローブ10の位置や姿勢が調整され、これによって特定の局所部位が含まれるように走査面Sが位置決めされる。そして、さらにその走査面S上において特定のビーム方位が指定され、そのビーム方位上においてエコートラッキング技術を用いて血管14の直径が計測される。これについては後に詳述するが、その場合においてはその特定のビーム方位上において前壁14A及び後壁14Bをそれぞれトラッキングするためのトラッキングゲートが設定される。
【0032】
送信部16は送信ビームフォーマーとして機能し、複数の振動素子に対して所定の遅延関係をもって複数の送信信号を供給する。一方、受信部18は複数の振動素子から出力される複数の受信信号に対して整相加算を実行する。すなわち受信部18は受信ビームフォーマーとして機能する。
【0033】
制御部20は、上記の送信部16及び受信部18の動作制御の他、装置内における各構成の動作制御を行っている。制御部20には操作パネル22が接続されている。この操作パネル22はトラックボールやキーボードなどによって構成される。受信部18から出力される整相加算後の受信信号は、本実施形態において、超音波画像形成部24、血管変位演算部26、血流速度演算部28に出力されている。
【0034】
超音波画像形成部24は、公知のBモード画像を形成するために必要な信号処理を実行する。その信号処理には、例えば検波、対数圧縮などの公知の信号処理が含まれる。
【0035】
血管変位演算部26は、上述したエコートラッキング技術を用いて、血管における局所部位の直径を演算する。具体的には、2つのトラッキングゲート内において前壁14A及び後壁14Bの拍動を検知し、すなわち、それぞれの壁の変位から直径を演算している。これ自体は公知の技術である。本実施形態においては、この血管変位演算部26が1心拍内における最高血圧に対応した直径Dsと最低血圧に対応した直径Ddとを求める機能を有している。もちろん、各時刻における直径Dについても演算されている。最高血圧及び最低血圧に対応する直径については後述する圧力換算部32において特定するようにしてもよい。
【0036】
血流速度演算部28はドプラ法に基づいて走査面内における各点の血流の速度を演算している。
【0037】
カフ型血圧計30は例えば被検者の上腕部に巻き付けられる公知の血圧計であり、その血圧計によって最高血圧Ps及び最低血圧Pdが測定される。それらの情報は圧力換算部32に出力されている。
【0038】
画像表示処理部34は、例えばDSC(デジタルスキャンコンバータ)などを含む表示処理ユニットとして構成されており、超音波画像形成部24から出力される画像データに対して座標変換や補間処理などを実行する。また、この表示処理部34によって血管における局所部位の血圧波形が作成され、また血流速度演算部28から出力される速度情報に基づいて二次元血流画像(カラードプラ画像)を作成する。さらに、後述する圧力換算部32から出力される換算後の圧力Pの情報に基づいて局所圧力波形を作成する機能を有している。もちろん、この画像表示処理部34が例えばwave intensityなどの指標を演算する機能していてもよい。
【0039】
圧力換算部32は、最大血圧Ps及び最低血圧Pdと、それらに対応する最大直径Ds及び最小直径Ddとに基づいて、各時刻において入力される直径Dに対して上述した(6)式の非線形演算を実行し、これによって局所部位の圧力Pを推定演算している。その情報は画像表示処理部34へ出力されており、上述したように画像表示処理部34は各時刻の局所圧力Pに基づいて局所圧力波形を形成している。
【0040】
表示部36には、二次元断層画像、カラーフローマッピング画像などが表示され、また本実施形態においては上記の局所圧力波形がグラフあるいは数値として表示される。またそのような局所圧力波形の情報に基づいて演算される各種の情報についても数値あるいはグラフで表示部36に表示される。
【0041】
図2は、図1に示される装置の動作例がフローチャートとして概念的に示されている。
【0042】
まず、S101では、プローブ10によって超音波の送受波が開始され、これと共にカフ型血圧計30によって血圧測定が開始される。S102では、カフ型血圧計30により最高血圧Ps及び最低血圧Pdが取得される。一方、S103では、最高血圧Psに対応する最大直径Dsと、最低血圧Pdに対応する最小直径Ddが取得される。この場合においては、望ましくは1心拍内においてそれらの情報Ds,Ddを特定するようにしてもよいし、複数心拍分のアンサンブル平均をとることにより、それらの情報を取得するようにしてもよい。ただし、血圧の測定と直径の演算はできる限り時間的に近接しているのが望ましい。
【0043】
例えば、S102の工程においては、1心拍内における最高血圧Psと最低血圧Pdを取得し、当該心拍内において最大直径Ds及び最小直径Ddを取得するようにしてもよい。
【0044】
S104では、上述した(6)式を実行することにより、すなわち非線形関数に上述した各パラメータ値Ps,Pd,Ds,Ddを与えて、入力される直径Dを換算することにより、結果として局所血圧Pが求められる。ちなみに、この演算はリアルタイムで実行されており、各時刻において換算された局所血圧Pが求められることになる。
【0045】
S105では、そのように求められた局所血圧Pの時間変化を表す局所血圧波形が表示部36上に表示され、あるいはそのような情報が数値表示される。S106では、そのような局所血圧波形に基づいて演算される他の情報、望ましくはwave intensity波形が表示部36に表示される。
【0046】
S107において上記演算を続行する場合には、再び上述した各工程が繰り返し実行されることになる。その場合においてはS102及びS103の工程を再び実行するようにしてもよいし、それらの工程によって得られたパラメータ値をそのまま維持して上記の非線形関数の演算を連続的に実行するようにしてもよい。
【0047】
図3には、本実施形態の特性202と従来例における特性200とがグラフとして示されている。
【0048】
このグラフの横軸はある時刻における局所部位の血管直径であり、その縦軸はその局所部位における局所圧力である。従来例200においては線形関数が利用されていたが、本実施形態においては上記のような非線形関数が用いられているため、直径Dに対して従来とは異なる圧力Pが求められる。
【0049】
図4には、従来例と本実施形態との特性の対比を行った場合における圧力差分がグラフとして示されている。例えば両者の圧力差分は最大で7%程度となり、本実施形態によれば真の圧力値に従来よりも近づけることが可能となる。
【0050】
図5には本実施形態による演算で求められた圧力波形が示されており、その横軸は時間を表しており、その縦軸は圧力を表している。ここで符号204は従来例であり、符号206は本実施形態である。このように両者が異なることから、例えばwave intensityを演算する場合においても従来よりもその値をより正確に求めることが可能である。本発明者の実験によれば、従来よりも本実施形態の方がwave intensityの波形において第1ピークが低くなり、第2ピークが高くなっている。
【0051】
図6には、他の実施形態の構成が示されている。なお、図1に示す構成と同様の構成には同一符号を付しその説明を省略する。
【0052】
図6に示す実施形態においては、血管14の短軸断面が利用される。図1においては走査面Sが血管14の長軸断面をとるようにプローブ10が位置決めされたが、図6に示す実施形態では、プローブ10が血管14の短軸断面を取得できるように血管軸と直交する方向に走査面Sが設定される。そして、画像表示処理部34によって例えばBモード断層画像が形成された場合、その画像情報が面積演算部50に出力され、その面積演算部50は断層画像上において血管14内に存在する画素数をカウントすることにより、あるいは公知の各種の手法を利用することにより、血管14の断面積を演算する。この場合に、二次元血流画像を利用してもよい。
【0053】
以上により、最大血圧Psに対応する最大面積Ssと、最低血圧Pdに対応する最小面積Sdとが特定され、それらの情報が圧力換算部52に出力される。
【0054】
圧力換算部52は上記の(20)式を実行することにより、上述した(6)式と同様に、特定されたパラメータ値に従って、入力される面積Sを換算することにより局所部位の圧力Pを求めている。この場合においても、面積演算にあたっては1心拍内における最大面積と最小面積を求めるようにしてもよいし、複数心拍分にわたってアンサンブル平均を求めることにより、最大面積及び最小面積を求めるようにしてもよい。
【0055】
【発明の効果】
以上説明したように、本発明によれば局所血圧を従来よりも正確に求めることが可能となる。
【図面の簡単な説明】
【図1】本発明に係る超音波診断装置の好適な実施形態を示すブロック図である。
【図2】図1に示す装置の動作例を説明するためのフローチャートである。
【図3】直径と圧力との関係を示すグラフである。
【図4】従来例の特性と本実施形態の特性との間における差分を説明するためのグラフである。
【図5】局所圧力波形を示すグラフである。
【図6】他の実施形態に係る超音波診断装置の構成を示すブロック図である。
【符号の説明】
10 プローブ、14 血管、16 送信部、18 受信部、20 制御部、24 超音波画像形成部、26 血管変位演算部、28 血流速度演算部、30
カフ型血圧計、32 圧力換算部、34 画像表示処理部。[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to an ultrasonic diagnostic apparatus, and more particularly, to measurement of local blood pressure using ultrasonic waves.
[0002]
[Prior art]
A catheter-type sphygmomanometer is used to obtain a blood pressure waveform for a specific local site in a blood vessel. That is, a catheter is inserted into a blood vessel, and a local blood pressure at a specific site is measured by a pressure sensor provided at a distal end thereof. However, this technique is invasive and invasive, and places a heavy burden on patients. As a method of noninvasive blood pressure measurement, a tonometry method in which a pressure pulse wave of a radial artery is measured by a pressure sensor is known. However, since the absolute blood pressure cannot be obtained, it is necessary to calibrate the detection value of the pressure sensor using the blood pressure value measured by the cuff type sphygmomanometer wound around the upper arm.
[0003]
By the way, a non-invasive blood pressure measurement method using ultrasound is known (Reference 1: Motoaki Sugawara et al., Development of non-invasive measurement method of blood pressure waveform, Medical Electronics and Biotechnology, Vol. 21, pp. 429, 1983). Etc.). That is, since the blood pressure change and the blood vessel diameter change are approximately in a linear relationship, the blood vessel diameter change at the measurement site is measured using the ultrasonic echo tracking method, and the systolic blood pressure and the diastolic blood pressure are measured with the cuff type blood pressure monitor. The maximum blood vessel diameter is calibrated with the systolic blood pressure, and the minimum blood vessel diameter is calibrated with the diastolic blood pressure, whereby the change in blood vessel diameter is regarded as a blood pressure waveform.
[0004]
[Problems to be solved by the invention]
It is known that the relationship between the blood pressure change and the blood vessel diameter change is not strictly linear but has a non-linear characteristic (Hansen F, et al., Diameter and compliance in the human common artery-variations-reference). with age and sex, Ultrasound Med & Biro 21, pp 1-9, 1995. and Reference 3: Hayashi K, et al., Mechanical properties, human arteries, et al. That is, as the blood pressure increases, the artery becomes stiff, and the change in blood vessel diameter decreases. There is also arteriosclerosis due to aging and disease. Despite these factors, errors may increase if the relationship between the change in blood vessel diameter and the change in blood pressure is regarded as linear. Therefore, if a circulatory dynamic index, for example, a known wave intensity is calculated using the blood pressure waveform obtained by the linear approximation, a measurement error also occurs.
[0005]
It is an object of the present invention to more accurately measure local blood pressure using ultrasound.
[0006]
[Means for Solving the Problems]
The present invention is a wave transmitting and receiving means for transmitting and receiving ultrasonic waves to and from a subject, based on a received signal obtained by transmitting and receiving the ultrasonic waves, for a specific local site in a blood vessel in the subject, Blood vessel size waveform information calculating means for obtaining local size waveform information representing a time change of the size; a blood pressure monitor mounted on the body surface of the subject to measure systolic blood pressure and diastolic blood pressure; Estimating means for estimating local blood pressure waveform information representing a temporal change of the local blood pressure for the specific local site by converting the blood vessel size waveform information based on systolic blood pressure and diastolic blood pressure. Means for giving the systolic blood pressure and the diastolic blood pressure to a nonlinear function for estimating the local blood pressure waveform information from the local size waveform information, And executes the calculation for obtaining the information.
[0007]
According to the above configuration, the local blood pressure waveform information is obtained using the nonlinear function for estimating the local blood pressure waveform information from the blood vessel size waveform information. Therefore, the accuracy of estimating the local blood pressure can be further improved as compared with the case where the conventional linear function is used.
[0008]
Preferably, the local size is a blood vessel diameter. Preferably, the systolic blood pressure and P s, the vessel diameter at this time is D s, whereas, the minimum blood pressure and P d, the vessel diameter at this time is D d, the vessel diameter at a certain time is represented by D In this case, the local blood pressure P at the time is
P = Pd * exp [β (D / Dd- 1)]
However, β = ln (P s / P d) / (D s / D d -1)
Required by
[0009]
Preferably, the local size is a blood vessel area. Preferably, the systolic blood pressure and P s, the vessel area at this time is S s, whereas, the minimum blood pressure and P d, the vessel area at this time is S d, the vessel area at a certain time was S In this case, the local blood pressure P at the time is
P = P d * exp [β s (S / S d −1)]
Where β s = ln (P s / P d ) / (S s / S d −1)
Required by
[0010]
For example, the blood vessel diameter for a local site can be obtained as the distance between two tracking points (gates) by performing echo tracking on the front wall and the back wall of the blood vessel on the ultrasonic beam. As the area of the blood vessel cross section (short axis cross section), a known KI method, A-SMA method, or the like, which traces a heart wall from a B-mode image to obtain a left ventricular area, can be used. Further, the blood vessel cross-sectional area may be obtained from the blood flow display by the color Doppler method.
[0011]
The local blood pressure waveform information obtained as described above is used, for example, in the calculation of the wave intensity. Of course, the information may be displayed as a numerical value or as a graph.
[0012]
BEST MODE FOR CARRYING OUT THE INVENTION
First, the principle of the embodiment according to the present invention will be described.
[0013]
Systolic blood pressure obtained by the cuff-type sphygmomanometer is wound around the upper arm portion (systolic pressure) and P s, the vessel diameter at this time is D s. Similarly, the lowest blood pressure (diastolic pressure) and P d, the vessel diameter at this time is D d.
[0014]
First, the relationship between (A): local blood pressure and blood vessel diameter (local diameter) will be examined. Now, it is assumed that the nonlinear relationship between the local blood pressure and the local diameter at a local site in a blood vessel is represented by the following exponential function having two parameters a and b.
[0015]
P = P (D, a, b)
= A * exp (b * D) (1)
From the above equation (1), the systolic blood pressure and the diastolic blood pressure are expressed by the following equations.
[0016]
P s = a * exp (b * D s) (2)
Pd = a * exp (b * Dd ) (3)
From the relationship between the above expressions (2) and (3), a and b are expressed by the following expressions.
[0017]
b = ln (P s / P d) / (D s -D d) (4)
a = Pd / exp (b * Dd ) (5)
Therefore, the relationship between the local blood pressure P and the blood vessel diameter D is represented by the following equation.
[0018]
P = Pd * exp [b * (D- Dd )]
= Pd * exp [β (D / Dd- 1)] (6)
Where β is
β = ln (P s / P d) / (D s / D d -1) (7)
Is a blood vessel elasticity index called a stiffness parameter defined by (see Document 3). That is, the blood vessels become harder as β becomes larger.
[0019]
The above results indicate that when the characteristic of local blood pressure-diameter of blood vessel is approximated by an exponential function, the exponent is a stiffness parameter.
[0020]
Next, similarly to the above, the relationship of (B): local blood pressure-vascular area (short-axis cross-sectional area) will be examined. The systolic blood pressure obtained by the cuff-type sphygmomanometer and P s, the vessel area at this time is S s. Also, the minimum blood pressure and P d, the vessel area at this time is S d.
[0021]
Now, it is assumed that the nonlinear relationship between the local blood pressure P and the blood vessel short-axis cross-sectional area S is represented by the following exponential function having two parameters A and B.
[0022]
P = P (S, A, B)
= A * exp (B * S) (8)
From the above equation (8), the systolic blood pressure and the diastolic blood pressure are expressed by the following equations.
[0023]
P s = A * exp (B * S s) (9)
Pd = A * exp (B * Sd ) (10)
From the relationship between the above equations (9) and (10), A and B are expressed by the following equations.
[0024]
B = ln (P s / P d) / (S s -S d) (11)
A = P d / exp (B * S d ) (12)
Therefore, the relationship between the local blood pressure P and the blood vessel area S is represented by the following equation.
[0025]
P = Pd * exp [B * (S- Sd )]
= P d * exp [β s (S / S d -1)] (13)
Here, the stiffness parameters are as follows.
[0026]
β s = ln (P s / P d) / (S s / S d -1) (14)
When the blood vessel cross-section is circular, the area S by using the diameter D of the vessel S = [pi] D 2/4
Is represented by The maximum vessel diameter D s and vessel diameter D is the minimum blood vessel diameter D d and the change amount [Delta] D s therefrom, as the sum of [Delta] D, are respectively expressed by the following equation.
[0027]
D s = D d + ΔD s (15)
D = D d + ΔD (16)
Since the blood vessel diameter changes ΔD s and ΔD are sufficiently smaller than the minimum blood vessel diameter D d , that is, ΔD s / D d << 1, ΔD / D d << 1,
S / S d -1 = (D / D d) 2 -1 ≒ 2 (D / D d -1) (17)
S s / S d −1 = (D s / D d ) 2 −1 ≒ 2 (D s / D d −1) (18)
Holds. Therefore,
β s ≒ β / 2 (19)
P = P d * exp [β s (S / S d -1)] (20)
≒ P d * exp [β (D / D d -1)] (21)
Can be approximated. From the above equations (20) and (21), almost the same results are obtained when the relationship between local blood pressure and blood vessel area is approximated by an exponential function and when the relationship between blood pressure and blood vessel diameter is approximated by an exponential function. Is understood.
[0028]
Next, an embodiment of the present invention will be described with reference to the drawings.
[0029]
FIG. 1 shows a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention, and FIG. 1 is a block diagram showing the entire configuration.
[0030]
In FIG. 1, a
[0031]
A
[0032]
The
[0033]
The
[0034]
The ultrasonic
[0035]
The blood
[0036]
The blood
[0037]
Cuff-
[0038]
The image
[0039]
The
[0040]
The
[0041]
FIG. 2 conceptually shows an operation example of the apparatus shown in FIG. 1 as a flowchart.
[0042]
First, in S101, transmission and reception of ultrasonic waves are started by the
[0043]
For example, in S102 of the process to obtain the systolic blood pressure P s and the minimum blood pressure P d in the one heartbeat, may acquire a maximum diameter D s and the minimum diameter D d at the relevant heart.
[0044]
In S104, by executing the above-mentioned (6), that is, each parameter value P s described above to a nonlinear function, by P d, D s, giving D d, converts the diameter D to be input, As a result, the local blood pressure P is determined. Incidentally, this calculation is executed in real time, and the local blood pressure P converted at each time is obtained.
[0045]
In S105, the local blood pressure waveform indicating the time change of the local blood pressure P thus obtained is displayed on the
[0046]
When the above calculation is continued in S107, each of the above-described steps is repeatedly executed. In that case, the steps of S102 and S103 may be executed again, or the calculation of the above-mentioned nonlinear function may be executed continuously while maintaining the parameter values obtained in those steps. Good.
[0047]
FIG. 3 is a graph showing the characteristic 202 of the present embodiment and the characteristic 200 of the conventional example.
[0048]
The horizontal axis of this graph is the blood vessel diameter at a local site at a certain time, and the vertical axis is the local pressure at the local site. In the conventional example 200, a linear function is used. However, in the present embodiment, since a non-linear function as described above is used, a pressure P different from the conventional pressure is obtained for the diameter D.
[0049]
FIG. 4 is a graph showing a pressure difference when the characteristics of the conventional example and the present embodiment are compared. For example, the pressure difference between the two is about 7% at the maximum, and according to the present embodiment, it is possible to approach the true pressure value more than before.
[0050]
FIG. 5 shows a pressure waveform obtained by the calculation according to the present embodiment. The horizontal axis represents time, and the vertical axis represents pressure. Here,
[0051]
FIG. 6 shows a configuration of another embodiment. The same components as those shown in FIG. 1 are denoted by the same reference numerals, and description thereof will be omitted.
[0052]
In the embodiment shown in FIG. 6, a short-axis cross section of the
[0053]
Thus, the maximum area S s corresponding to the maximum blood pressure P s, is identified with minimum area S d corresponding to the minimum blood pressure P d, the information is output to the
[0054]
The
[0055]
【The invention's effect】
As described above, according to the present invention, the local blood pressure can be obtained more accurately than before.
[Brief description of the drawings]
FIG. 1 is a block diagram showing a preferred embodiment of an ultrasonic diagnostic apparatus according to the present invention.
FIG. 2 is a flowchart for explaining an operation example of the device shown in FIG. 1;
FIG. 3 is a graph showing the relationship between diameter and pressure.
FIG. 4 is a graph for explaining a difference between a characteristic of a conventional example and a characteristic of the present embodiment.
FIG. 5 is a graph showing a local pressure waveform.
FIG. 6 is a block diagram illustrating a configuration of an ultrasonic diagnostic apparatus according to another embodiment.
[Explanation of symbols]
DESCRIPTION OF
Cuff type sphygmomanometer, 32 pressure conversion unit, 34 image display processing unit.
Claims (5)
前記超音波の送受波により得られた受信信号に基づいて、前記被検体内の血管における特定の局所部位について、そのサイズの時間変化を表す局所サイズ波形情報を求める血管サイズ波形情報演算手段と、
前記被検体の体表面に装着され、最高血圧及び最低血圧を測定する血圧計と、前記血圧計によって測定された最高血圧及び最低血圧に基づいて前記局所サイズ波形情報を換算することにより、前記特定の局所部位について、その局所血圧の時間変化を表す局所血圧波形情報を推定する推定手段と、
を含み、
前記推定手段は、前記局所サイズ波形情報から前記局所血圧波形情報を推定する非線形関数に対して、前記最高血圧及び最低血圧を与えることにより、前記局所血圧波形情報を求める演算を実行することを特徴とする超音波診断装置。Transmitting and receiving means for transmitting and receiving ultrasonic waves to and from the subject,
Based on a received signal obtained by the transmission and reception of the ultrasonic waves, for a specific local site in a blood vessel in the subject, a blood vessel size waveform information calculating means for obtaining local size waveform information representing a temporal change in the size,
A sphygmomanometer which is attached to the body surface of the subject and measures systolic and diastolic blood pressures, and converts the local size waveform information based on the systolic and diastolic blood pressures measured by the sphygmomanometer, whereby the identification is performed. Estimating means for estimating local blood pressure waveform information representing a time change of the local blood pressure,
Including
The estimating means performs an operation for obtaining the local blood pressure waveform information by giving the systolic blood pressure and the diastolic blood pressure to a nonlinear function for estimating the local blood pressure waveform information from the local size waveform information. Ultrasound diagnostic device.
前記局所サイズは血管直径であることを特徴とする超音波診断装置。The device of claim 1,
An ultrasonic diagnostic apparatus, wherein the local size is a blood vessel diameter.
前記最高血圧をPsとし、このときの血管直径をDsとし、一方、前記最低血圧をPdとし、このときの血管直径をDdとし、ある時刻における血管直径をDとした場合に、当該時刻における局所血圧Pが、
P = Pd*exp[β(D/Dd−1)]
但し、β= ln(Ps/Pd)/(Ds/Dd−1)
によって求められることを特徴とする超音波診断装置。The device according to claim 2,
The systolic blood pressure and P s, the vessel diameter at this time as D s, whereas, the minimum blood pressure and P d, the vessel diameter at this time is a D d, the vessel diameter at a certain time when as D, The local blood pressure P at the time is
P = Pd * exp [β (D / Dd- 1)]
However, β = ln (P s / P d) / (D s / D d -1)
An ultrasonic diagnostic apparatus characterized by the following.
前記局所サイズは血管面積であることを特徴とする超音波診断装置。The device of claim 1,
An ultrasonic diagnostic apparatus, wherein the local size is a blood vessel area.
前記最高血圧をPsとし、このときの血管面積をSsとし、一方、前記最低血圧をPdとし、このときの血管面積をSdとし、ある時刻における血管面積をSとした場合に、当該時刻における局所血圧Pが、
P= Pd*exp[βs(S/Sd−1)]
但し、βs = ln(Ps/Pd)/(Ss/Sd−1)
によって求められることを特徴とする超音波診断装置。The device according to claim 4,
When the systolic blood pressure is P s and the vascular area at this time is S s , while the diastolic blood pressure is P d , the vascular area at this time is S d, and the vascular area at a certain time is S, The local blood pressure P at the time is
P = P d * exp [β s (S / S d −1)]
Where β s = ln (P s / P d ) / (S s / S d −1)
An ultrasonic diagnostic apparatus characterized by the following.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002201799A JP4091365B2 (en) | 2002-07-10 | 2002-07-10 | Ultrasonic diagnostic equipment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002201799A JP4091365B2 (en) | 2002-07-10 | 2002-07-10 | Ultrasonic diagnostic equipment |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004041382A true JP2004041382A (en) | 2004-02-12 |
JP4091365B2 JP4091365B2 (en) | 2008-05-28 |
Family
ID=31708225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002201799A Expired - Fee Related JP4091365B2 (en) | 2002-07-10 | 2002-07-10 | Ultrasonic diagnostic equipment |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4091365B2 (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010214013A (en) * | 2009-03-18 | 2010-09-30 | Fujifilm Corp | Ultrasonograph, method for measuring pressure gradient, and method for measuring blood vessel elasticity |
CN102247169A (en) * | 2010-05-19 | 2011-11-23 | 精工爱普生株式会社 | Blood pressure measuring device and blood pressure measuring method |
JP2012065747A (en) * | 2010-09-22 | 2012-04-05 | Seiko Epson Corp | Device and method for measuring fluid pulsation pressure, and device and method for measuring pulse pressure |
JP2013123529A (en) * | 2011-12-14 | 2013-06-24 | Seiko Epson Corp | Device and method of measuring blood pressure |
US20130281852A1 (en) * | 2012-04-18 | 2013-10-24 | Seiko Epson Corporation | Method for correcting blood pressure estimation parameter and blood pressure measurement apparatus |
JP2013220243A (en) * | 2012-04-18 | 2013-10-28 | Seiko Epson Corp | Central blood pressure estimation parameter calibration method, central blood pressure measurement method and central blood pressure measurement device |
JP2014104163A (en) * | 2012-11-28 | 2014-06-09 | Seiko Epson Corp | Apparatus and method for measuring blood pressure |
JP5701439B1 (en) * | 2014-08-20 | 2015-04-15 | 医療法人 駿東育愛会 望星第一クリニック | Subcutaneous diagnostic apparatus and blood pressure measurement method |
US20150112214A1 (en) * | 2013-10-23 | 2015-04-23 | Seiko Epson Corporation | Blood pressure measuring device and blood pressure measuring method |
CN104887195A (en) * | 2014-03-06 | 2015-09-09 | 日本光电工业株式会社 | Vascular endothelial function inspection apparatus |
JP2017522063A (en) * | 2014-05-19 | 2017-08-10 | クアルコム,インコーポレイテッド | Non-interfering blood pressure measurement |
JP2017217474A (en) * | 2016-06-07 | 2017-12-14 | 東芝メディカルシステムズ株式会社 | Medical image diagnostic apparatus and medical image processing system |
JP2017538539A (en) * | 2014-12-22 | 2017-12-28 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | System and method for measuring arterial parameters |
JP2022121340A (en) * | 2021-02-08 | 2022-08-19 | 阿部 真一 | Calculating blood flow rate, blood vessel diameter, blood vessel flow resistance, heart load factor, and the like according to measurement result of blood pressure manometer and displaying them to blood pressure manometer |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9826959B2 (en) | 2008-11-04 | 2017-11-28 | Fujifilm Corporation | Ultrasonic diagnostic device |
-
2002
- 2002-07-10 JP JP2002201799A patent/JP4091365B2/en not_active Expired - Fee Related
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010214013A (en) * | 2009-03-18 | 2010-09-30 | Fujifilm Corp | Ultrasonograph, method for measuring pressure gradient, and method for measuring blood vessel elasticity |
US8444565B2 (en) | 2009-03-18 | 2013-05-21 | Fujifilm Corporation | Ultrasonic diagnostic apparatus, method of measuring pressure gradient, and method of measuring blood vessel elasticity |
CN102247169A (en) * | 2010-05-19 | 2011-11-23 | 精工爱普生株式会社 | Blood pressure measuring device and blood pressure measuring method |
JP2011239972A (en) * | 2010-05-19 | 2011-12-01 | Seiko Epson Corp | Blood pressure measuring device and method |
JP2012065747A (en) * | 2010-09-22 | 2012-04-05 | Seiko Epson Corp | Device and method for measuring fluid pulsation pressure, and device and method for measuring pulse pressure |
JP2013123529A (en) * | 2011-12-14 | 2013-06-24 | Seiko Epson Corp | Device and method of measuring blood pressure |
JP2013220243A (en) * | 2012-04-18 | 2013-10-28 | Seiko Epson Corp | Central blood pressure estimation parameter calibration method, central blood pressure measurement method and central blood pressure measurement device |
JP2013220242A (en) * | 2012-04-18 | 2013-10-28 | Seiko Epson Corp | Method for calibrating blood pressure estimation parameter, and blood pressure measurement apparatus |
US20130281852A1 (en) * | 2012-04-18 | 2013-10-24 | Seiko Epson Corporation | Method for correcting blood pressure estimation parameter and blood pressure measurement apparatus |
JP2014104163A (en) * | 2012-11-28 | 2014-06-09 | Seiko Epson Corp | Apparatus and method for measuring blood pressure |
US20150112214A1 (en) * | 2013-10-23 | 2015-04-23 | Seiko Epson Corporation | Blood pressure measuring device and blood pressure measuring method |
CN104887195A (en) * | 2014-03-06 | 2015-09-09 | 日本光电工业株式会社 | Vascular endothelial function inspection apparatus |
JP2017522063A (en) * | 2014-05-19 | 2017-08-10 | クアルコム,インコーポレイテッド | Non-interfering blood pressure measurement |
JP5701439B1 (en) * | 2014-08-20 | 2015-04-15 | 医療法人 駿東育愛会 望星第一クリニック | Subcutaneous diagnostic apparatus and blood pressure measurement method |
JP2016042903A (en) * | 2014-08-20 | 2016-04-04 | 医療法人 駿東育愛会 望星第一クリニック | Subcutaneous diagnostic system, and blood pressure measuring method |
JP2017538539A (en) * | 2014-12-22 | 2017-12-28 | コーニンクレッカ フィリップス エヌ ヴェKoninklijke Philips N.V. | System and method for measuring arterial parameters |
JP2017217474A (en) * | 2016-06-07 | 2017-12-14 | 東芝メディカルシステムズ株式会社 | Medical image diagnostic apparatus and medical image processing system |
JP2022121340A (en) * | 2021-02-08 | 2022-08-19 | 阿部 真一 | Calculating blood flow rate, blood vessel diameter, blood vessel flow resistance, heart load factor, and the like according to measurement result of blood pressure manometer and displaying them to blood pressure manometer |
Also Published As
Publication number | Publication date |
---|---|
JP4091365B2 (en) | 2008-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3268396B2 (en) | Ultrasound diagnostic equipment | |
JP4667394B2 (en) | Ultrasonic diagnostic equipment | |
US10716538B2 (en) | Hemodynamic ultrasound medical monitoring device | |
JP4884476B2 (en) | Ultrasonic diagnostic equipment | |
Seo et al. | Noninvasive arterial blood pressure waveform monitoring using two-element ultrasound system | |
JP4766546B2 (en) | Ultrasonic diagnostic equipment | |
JP5158880B2 (en) | Ultrasonic diagnostic equipment | |
WO2008023618A1 (en) | Ultrasonograph | |
JP4091365B2 (en) | Ultrasonic diagnostic equipment | |
JP2003010183A (en) | Ultrasonograph | |
WO2003077765A1 (en) | Ultrasonographic system and ultrasonography | |
WO2006129545A1 (en) | Ultrasonograph | |
WO2004103185A1 (en) | Ultrasonograph | |
WO2006043528A1 (en) | Ultrasonographic device and ultrasonographic device control method | |
JP4918369B2 (en) | Ultrasonic diagnostic equipment | |
JP2009039277A (en) | Ultrasonic diagnostic apparatus | |
JP5346555B2 (en) | Ultrasound diagnostic device with arteriosclerosis risk display function | |
WO2007080870A1 (en) | Ultrasonograph | |
CN115844452B (en) | Pulse wave detection method, device and storage medium | |
JP2008161546A (en) | Ultrasonic diagnostic apparatus | |
JP2004159672A (en) | Ultrasonograph and ultrasonic measuring method | |
JP3668687B2 (en) | Pulse wave velocity measuring device and ultrasonic diagnostic device | |
JP5462474B2 (en) | Ultrasonic diagnostic equipment | |
WO2007011000A1 (en) | Ultrasonographic device | |
JP6098641B2 (en) | Ultrasonic diagnostic apparatus, control method for ultrasonic diagnostic apparatus, and controller for ultrasonic diagnostic apparatus |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050614 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20071122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071211 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080125 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080226 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080228 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140307 Year of fee payment: 6 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |