JP2004041382A - Ultrasonograph - Google Patents

Ultrasonograph Download PDF

Info

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
Application number
JP2002201799A
Other languages
Japanese (ja)
Other versions
JP4091365B2 (en
Inventor
Takemitsu Harada
原田 烈光
Takashi Okada
岡田 孝
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Aloka Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aloka Co Ltd filed Critical Aloka Co Ltd
Priority to JP2002201799A priority Critical patent/JP4091365B2/en
Publication of JP2004041382A publication Critical patent/JP2004041382A/en
Application granted granted Critical
Publication of JP4091365B2 publication Critical patent/JP4091365B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To accurately compute a blood pressure of a local area (local blood pressure) in a blood vessel by utilizing a transmitted/received wave of an ultrasonic wave. <P>SOLUTION: A maximum blood pressure P<SB>s</SB>and a minimum blood pressure P<SB>d</SB>are obtained by a cuff type hemodynamometer 30. Also, a maximum diameter D<SB>s</SB>and a minimum diameter D<SB>d</SB>regarding a local area are obtained corresponding to them, and those parameter values are provided to a non-linear function. By the non-linear function, a diameter D at each time to be input is reduced, and a pressure P at each time regarding the local area is computed. In the place of the diameter, the area of the blood vessel (the area on a minor axis cross section) may be utilized as well. <P>COPYRIGHT: (C)2004,JPO

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】
望ましくは、前記局所サイズは血管直径である。望ましくは、前記最高血圧をPとし、このときの血管直径をDとし、一方、前記最低血圧をPとし、このときの血管直径をDとし、ある時刻における血管直径をDとした場合に、当該時刻における局所血圧Pが、
P = P*exp[β(D/D−1)]
但し、β= ln(P/P)/(D/D−1)
によって求められる。
【0009】
望ましくは、前記局所サイズは血管面積である。望ましくは、前記最高血圧をPとし、このときの血管面積をSとし、一方、前記最低血圧をPとし、このときの血管面積をSとし、ある時刻における血管面積をSとした場合に、当該時刻における局所血圧Pが、
P=  P*exp[β(S/S−1)]
但し、β = ln(P/P)/(S/S−1)
によって求められる。
【0010】
局所部位についての血管直径は、例えば、超音波ビーム上において、血管の前壁と後壁についてエコートラッキングを行うことにより、2つのトラッキングポイント(ゲート)間の距離として求めることができる。血管断面(短軸断面)の面積は、Bモード画像から心臓壁をトレースして左室面積を求める公知のKI法やA−SMA法などを用いることができる。また、カラードプラ法の血流表示から血管断面積を求めてもよい。
【0011】
上記によって求められた局所血圧波形情報は、例えばwave intensityの演算で用いられる。もちろん、その情報を数値としてあるいはグラフとして表示するようにしてもよい。
【0012】
【発明の実施の形態】
まず、本発明に係る実施形態の原理について説明する。
【0013】
上腕部に巻き付けられるカフ型血圧計で得られた最高血圧(収縮期圧力)をPとし、このときの血管直径をDとする。同様に、最低血圧(拡張期圧力)をPとし、このときの血管直径をDとする。
【0014】
最初に、(A):局所血圧−血管直径(局所直径)の関係について検討する。いま、血管における局所部位についての局所血圧と局所直径の非線形関係が2つのパラメータa,bを有する次の指数関数で表されるとする。
【0015】
P = P(D,a,b)
= a*exp(b*D)          (1)
上記(1)式から,最高血圧と最低血圧は次式で表される。
【0016】
= a*exp(b*D)                    (2)
= a*exp(b*D)                    (3)
上記(2)式、(3)式の関係から、a,bは次式で表される。
【0017】
b= ln(P/P)/(D−D)         (4)
a= P/exp(b*D)          (5)
したがって、局所血圧Pと血管直径Dとの関係は次式で表される。
【0018】
P = P*exp[b*(D−D)]
= P*exp[β(D/D−1)]        (6)
ここで、βは、
β= ln(P/P)/(D/D−1)             (7)
で定義されるstiffness parameterと呼ばれる血管弾性指標である(文献3参照)。つまり、βが大きいほど、血管が硬くなる。
【0019】
以上の結果は、局所血圧−血管直径の特性を指数関数で近似したとき、その指数がstiffness parameterであることを意味している。
【0020】
次に、上記と同様に、(B):局所血圧−血管面積(短軸断面積)の関係について検討する。カフ型血圧計で得られた最高血圧をPとし、このときの血管面積をSとする。また、最低血圧をPとし、このときの血管面積をSとする。
【0021】
いま、局所血圧Pと血管短軸断面積Sの非線形関係が2つのパラメータA,Bを有する以下の指数関数で表されるとする。
【0022】
P = P(S,A,B)
= A*exp(B*S)                   (8)
上記(8)式から、最高血圧と最低血圧は次式で表される。
【0023】
= A*exp(B*S)                   (9)
= A*exp(B*S)                  (10)
上記(9)式、(10)式の関係から、A,Bは次式で表される。
【0024】
B= ln(P/P)/(S−S)               (11)
A= P/exp(B*S)                   (12)
したがって、局所血圧Pと血管面積Sとの関係は次式で表される。
【0025】
P= P*exp[B*(S−S)]
= P*exp[β(S/S−1)]              (13)
ここで、stiffness parameterは、以下の通りである。
【0026】
β = ln(P/P)/(S/S−1)              (14)
血管の断面を円形とすると、その面積Sは血管の直径Dを用いて
S = πD/4
で表される。また、最大血管直径Dと血管直径Dは、最小血管直径Dとそれからの変化分ΔD,ΔDの和として、それぞれ次式で表される。
【0027】
= D+ΔD                     (15)
D = D+ΔD                     (16)
血管直径変化ΔD,ΔDは、最小血管直径Dに比べて十分小さいので、つまり、ΔD/D<<1, ΔD/D<<1なので、
S/S−1 = (D/D−1 ≒2(D/D−1)        (17)
/S−1 =( D/D−1≒2(D/D−1)      (18)
が成り立つ。したがって、
β≒β/2                       (19)
P  = P*exp[β(S/S−1)]              (20)
≒ P*exp[β(D/D−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心拍内における最高血圧に対応した直径Dと最低血圧に対応した直径Dとを求める機能を有している。もちろん、各時刻における直径Dについても演算されている。最高血圧及び最低血圧に対応する直径については後述する圧力換算部32において特定するようにしてもよい。
【0036】
血流速度演算部28はドプラ法に基づいて走査面内における各点の血流の速度を演算している。
【0037】
カフ型血圧計30は例えば被検者の上腕部に巻き付けられる公知の血圧計であり、その血圧計によって最高血圧P及び最低血圧Pが測定される。それらの情報は圧力換算部32に出力されている。
【0038】
画像表示処理部34は、例えばDSC(デジタルスキャンコンバータ)などを含む表示処理ユニットとして構成されており、超音波画像形成部24から出力される画像データに対して座標変換や補間処理などを実行する。また、この表示処理部34によって血管における局所部位の血圧波形が作成され、また血流速度演算部28から出力される速度情報に基づいて二次元血流画像(カラードプラ画像)を作成する。さらに、後述する圧力換算部32から出力される換算後の圧力Pの情報に基づいて局所圧力波形を作成する機能を有している。もちろん、この画像表示処理部34が例えばwave intensityなどの指標を演算する機能していてもよい。
【0039】
圧力換算部32は、最大血圧P及び最低血圧Pと、それらに対応する最大直径D及び最小直径Dとに基づいて、各時刻において入力される直径Dに対して上述した(6)式の非線形演算を実行し、これによって局所部位の圧力Pを推定演算している。その情報は画像表示処理部34へ出力されており、上述したように画像表示処理部34は各時刻の局所圧力Pに基づいて局所圧力波形を形成している。
【0040】
表示部36には、二次元断層画像、カラーフローマッピング画像などが表示され、また本実施形態においては上記の局所圧力波形がグラフあるいは数値として表示される。またそのような局所圧力波形の情報に基づいて演算される各種の情報についても数値あるいはグラフで表示部36に表示される。
【0041】
図2は、図1に示される装置の動作例がフローチャートとして概念的に示されている。
【0042】
まず、S101では、プローブ10によって超音波の送受波が開始され、これと共にカフ型血圧計30によって血圧測定が開始される。S102では、カフ型血圧計30により最高血圧P及び最低血圧Pが取得される。一方、S103では、最高血圧Pに対応する最大直径Dと、最低血圧Pに対応する最小直径Dが取得される。この場合においては、望ましくは1心拍内においてそれらの情報D,Dを特定するようにしてもよいし、複数心拍分のアンサンブル平均をとることにより、それらの情報を取得するようにしてもよい。ただし、血圧の測定と直径の演算はできる限り時間的に近接しているのが望ましい。
【0043】
例えば、S102の工程においては、1心拍内における最高血圧Pと最低血圧Pを取得し、当該心拍内において最大直径D及び最小直径Dを取得するようにしてもよい。
【0044】
S104では、上述した(6)式を実行することにより、すなわち非線形関数に上述した各パラメータ値P,P,D,Dを与えて、入力される直径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】
以上により、最大血圧Pに対応する最大面積Sと、最低血圧Pに対応する最小面積Sとが特定され、それらの情報が圧力換算部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 probe 10 is used in contact with the body surface 12 in the embodiment shown in FIG. 1, and the probe 10 transmits and receives ultrasonic waves. Specifically, an array transducer composed of a plurality of transducers is provided in the probe 10, and the array transducer generates an ultrasonic beam B. A scanning surface S is formed by electronically scanning the ultrasonic beam B. In the example shown in FIG. 1, the ultrasonic beam B is electronically linearly scanned, but the electronic scanning method may be electronic sector scanning.
[0031]
A blood vessel 14 exists in the living body, and the blood vessel 14 has a front wall 14A and a rear wall 14B when viewed from the probe 10. When measuring a specific local part in the blood vessel 14, for example, the position and posture of the probe 10 are adjusted while observing a B-mode image (two-dimensional tomographic image) so that the specific local part is included. The scanning surface S is positioned. Then, a specific beam direction is further designated on the scanning plane S, and the diameter of the blood vessel 14 is measured on the beam direction using the echo tracking technique. As will be described in detail later, in this case, a tracking gate for tracking the front wall 14A and the rear wall 14B on the specific beam direction is set.
[0032]
The transmission unit 16 functions as a transmission beamformer, and supplies a plurality of transmission signals with a predetermined delay relationship to a plurality of vibration elements. On the other hand, the receiving unit 18 performs phasing addition on a plurality of reception signals output from the plurality of vibration elements. That is, the receiving unit 18 functions as a receiving beamformer.
[0033]
The control unit 20 controls the operation of each component in the apparatus in addition to the operation control of the transmission unit 16 and the reception unit 18 described above. An operation panel 22 is connected to the control unit 20. The operation panel 22 includes a trackball, a keyboard, and the like. The reception signal after the phasing addition output from the receiving unit 18 is output to the ultrasonic image forming unit 24, the blood vessel displacement calculating unit 26, and the blood flow velocity calculating unit 28 in the present embodiment.
[0034]
The ultrasonic image forming unit 24 performs signal processing necessary for forming a known B-mode image. The signal processing includes known signal processing such as detection and logarithmic compression.
[0035]
The blood vessel displacement calculator 26 calculates the diameter of a local part in a blood vessel using the above-described echo tracking technique. Specifically, the pulsation of the front wall 14A and the rear wall 14B is detected in the two tracking gates, that is, the diameter is calculated from the displacement of each wall. This is a known technique itself. In the present embodiment has a function of this blood vessel displacement calculation unit 26 obtains a diameter D d corresponding to the diameter D s and the lowest blood pressure corresponding to the systolic blood pressure within one heartbeat. Of course, the diameter D at each time is also calculated. The diameters corresponding to the systolic blood pressure and the diastolic blood pressure may be specified in a pressure conversion unit 32 described later.
[0036]
The blood flow velocity calculator 28 calculates the blood flow velocity at each point in the scanning plane based on the Doppler method.
[0037]
Cuff-type sphygmomanometer 30 is a known blood pressure monitor wound around the upper arm, for example, subject, systolic P s and the minimum blood pressure P d is measured by the sphygmomanometer. Those pieces of information are output to the pressure conversion unit 32.
[0038]
The image display processing unit 34 is configured as a display processing unit including, for example, a DSC (Digital Scan Converter), and performs coordinate conversion, interpolation processing, and the like on image data output from the ultrasonic image forming unit 24. . The display processing unit 34 creates a blood pressure waveform at a local site in a blood vessel, and creates a two-dimensional blood flow image (color Doppler image) based on the speed information output from the blood flow speed calculation unit 28. Further, it has a function of creating a local pressure waveform based on information of the converted pressure P output from the pressure conversion unit 32 described later. Of course, the image display processing unit 34 may have a function of calculating an index such as a wave intensity.
[0039]
The pressure conversion section 32, the maximum blood pressure P s and the diastolic blood pressure P d, based on the maximum diameter D s and the minimum diameter D d that correspond to them, the above-described with respect to the diameter D input at each time (6 ) Is performed to estimate and calculate the pressure P at the local site. The information is output to the image display processing unit 34, and the image display processing unit 34 forms a local pressure waveform based on the local pressure P at each time as described above.
[0040]
The display unit 36 displays a two-dimensional tomographic image, a color flow mapping image, and the like. In the present embodiment, the local pressure waveform is displayed as a graph or a numerical value. Various types of information calculated based on such information of the local pressure waveform are also displayed on the display unit 36 as numerical values or graphs.
[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 probe 10, and at the same time, blood pressure measurement is started by the cuff type sphygmomanometer 30. In S102, systolic P s and the minimum blood pressure P d is obtained by the cuff-type sphygmomanometer 30. On the other hand, in S103, the maximum diameter D s which corresponds to the systolic blood pressure P s, the minimum diameter D d corresponding to the minimum blood pressure P d is obtained. In this case, preferably the information D s in one heartbeat, may be identify D d, by taking the ensemble average of a plurality heartbeats, be acquired their information Good. However, it is desirable that the measurement of blood pressure and the calculation of the diameter be as close as possible in time.
[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 display unit 36, or such information is numerically displayed. In S106, other information calculated based on such a local blood pressure waveform, desirably, a wave intensity waveform is displayed on the display unit 36.
[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, reference numeral 204 is a conventional example, and reference numeral 206 is this embodiment. As described above, since the two are different from each other, it is possible to obtain the value more accurately than in the conventional case, for example, when calculating the wave intensity. According to the experiment of the inventor, in the present embodiment, the first peak is lower and the second peak is higher in the waveform of the wave intensity than in the related art.
[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 blood vessel 14 is used. In FIG. 1, the probe 10 is positioned so that the scanning plane S takes the long-axis cross section of the blood vessel 14. In the embodiment shown in FIG. 6, the blood vessel axis is set so that the probe 10 can acquire the short-axis cross section of the blood vessel 14. The scanning plane S is set in a direction orthogonal to the scanning plane S. When, for example, a B-mode tomographic image is formed by the image display processing unit 34, the image information is output to the area calculating unit 50, and the area calculating unit 50 calculates the number of pixels existing in the blood vessel 14 on the tomographic image. The cross-sectional area of the blood vessel 14 is calculated by counting or by using various known methods. In this case, a two-dimensional blood flow image may be used.
[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 pressure converting unit 52.
[0054]
The pressure conversion unit 52 executes the above equation (20) to convert the input area S in accordance with the specified parameter value to thereby calculate the pressure P of the local part, similarly to the above equation (6). I'm asking. Also in this case, in calculating the area, the maximum area and the minimum area in one heartbeat may be obtained, or the maximum area and the minimum area may be obtained by obtaining an ensemble average over a plurality of heartbeats. .
[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 SYMBOLS 10 probe, 14 blood vessel, 16 transmitting part, 18 receiving part, 20 control part, 24 ultrasonic image forming part, 26 blood vessel displacement calculating part, 28 blood flow velocity calculating part, 30
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.
請求項1記載の装置において、
前記局所サイズは血管直径であることを特徴とする超音波診断装置。
The device of claim 1,
An ultrasonic diagnostic apparatus, wherein the local size is a blood vessel diameter.
請求項2記載の装置において、
前記最高血圧をPとし、このときの血管直径をDとし、一方、前記最低血圧をPとし、このときの血管直径をDとし、ある時刻における血管直径をDとした場合に、当該時刻における局所血圧Pが、
P = P*exp[β(D/D−1)]
但し、β= ln(P/P)/(D/D−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.
請求項1記載の装置において、
前記局所サイズは血管面積であることを特徴とする超音波診断装置。
The device of claim 1,
An ultrasonic diagnostic apparatus, wherein the local size is a blood vessel area.
請求項4記載の装置において、
前記最高血圧をPとし、このときの血管面積をSとし、一方、前記最低血圧をPとし、このときの血管面積をSとし、ある時刻における血管面積をSとした場合に、当該時刻における局所血圧Pが、
P=  P*exp[β(S/S−1)]
但し、β = ln(P/P)/(S/S−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.
JP2002201799A 2002-07-10 2002-07-10 Ultrasonic diagnostic equipment Expired - Fee Related JP4091365B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9826959B2 (en) 2008-11-04 2017-11-28 Fujifilm Corporation Ultrasonic diagnostic device

Cited By (18)

* Cited by examiner, † Cited by third party
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