(57)【特許請求の範囲】
【請求項１】
動脈上の少なくとも２ヶ所の血管の断面の形状を同時に検出する断面形状検出手段と、
該断面形状検出手段により検出された断面形状に基づいて、前記２ヶ所の断面の断面積を算出する断面積算出手段と、
前記断面積算出手段により求められた断面積に基づいて血管内の血液の平均流速を算出する流速算出手段と、
前記断面形状検出手段により検出された２ヶ所の血管の断面形状の時間変化から前記２ヶ所の断面の間の脈波の伝搬速度を算出する脈波速度算出手段と、
前記流速算出手段により算出された血液の流速と、前記脈波速度算出手段により算出された脈波の伝搬速度とに基づいて、特定の時相における血圧値を基準とする相対血圧値を算出する相対血圧算出手段と、
あらかじめ測定されている基準となる前記特定の時相における血圧値と前記相対血圧値とに基づいて、絶対血圧値を算出する絶対血圧算出手段とを具備し、
前記流速算出手段は、前記２ヶ所の断面の断面積の時間微分と前記２ヶ所の断面の距離とから、前記２ヶ所の断面と血管壁により囲まれた部分の血液の平均流速を算出することを特徴とする循環器情報計測システム。
【請求項２】
前記断面形状検出手段は、前記動脈に向けて超音波を放射する超音波放射手段と、生体内で反射された超音波を検出する超音波検出手段と、前記動脈の管壁の断面形状を抽出するために、前記超音波検出手段からの出力信号を適切な閾値を基準として２値化する２値化手段とを備えることを特徴とする請求項１に記載の循環器情報計測システム。
【請求項３】
前記断面形状検出手段は、前記超音波検出手段からの出力信号を画像信号に変換する変換手段をさらに具備することを特徴とする請求項２に記載の循環器情報計測システム。
【請求項４】
前記２値化手段は、前記超音波検出手段からの出力信号を前記変換手段により画像信号に変換した後に、適切な画像濃度値を閾値として２値化することを特徴とする請求項３に記載の循環器情報計測システム。
【請求項５】
前記２値化手段は、超音波の入射方向と前記動脈の血管壁のなす角度に応じて前記閾値を設定する設定手段を備えることを特徴とする請求項２に記載の循環器情報計測システム。
【請求項６】
前記２値化手段は、前記断面形状検出手段により検出された血管の断面形状と前記超音波の入射方向とに基づいて、該超音波の入射方向と前記血管壁のなす角度を演算する演算手段をさらに備えることを特徴とする請求項５に記載の循環器情報計測システム。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は非侵襲で連続的に生体の循環器情報を計測することができる循環器情報計測システムに関するものである。
【０００２】
【従来の技術】
従来より、患者の循環器の状態を把握するための循環器指標として、血流量、血管の拡張・収縮速度による粘弾性などの諸指標が利用されている。
血流量の計測としては、生体の血管に対して超音波を照射し、血流により生じるドップラー信号を用いて血流速度を求め、血管の断面積との関係から、血流量を求める方法が知られている。
【０００３】
また、血管の拡張・収縮速度を求める方法としては、連続した画像データの輝度勾配を用いて速度ベクトルを求めるオプティカルフロー法が用いられている。
【０００４】
【発明が解決しようとする課題】
しかしながら、ドップラー信号を用いて血流速度を求め、血流量を求める場合、血流量を求めるために、その部位の血管の断面積を容易かつ連続的に算出する技術が望まれていた。
また、上記のオプティカルフロー法においては、複雑な逐次計算が必要となるという問題点があった。
【０００５】
また、従来、血圧を計測する方法としては、主にカフを用いて生体の一部を圧迫し、得られる振動や音から血圧を計測する方法がとられてきた。しかしながら、手術後など数日にわたって血圧を計測する場合には、生体を圧迫することがかなりの精神的かつ肉体的な負担となっている。
この問題を解決するために、超音波により血管内を伝搬する脈波の伝搬速度を計測し、血圧を計測する方法が、特開平７−２４１２８８号公報に開示されている。この方法では、
ΔＰ＝ρＣΔｖ
ここで、ΔＰ：圧力変化    ρ：血液密度
Δｖ：流速変化    Ｃ：脈波伝搬速度
なる理論式を用いて血圧を求めている。この理論式は（（社）日本油空圧学会編：新版油空圧便覧、第１編、第３章、２４−２５、オーム社、１９８８）に開示されている。しかし、この理論式は血管壁に変化のないことを仮定した式であり、管壁の影響が考慮されていないという問題点があった。
【０００６】
したがって、本発明は上述した課題に鑑みてなされたものであり、その目的は、非侵襲で血管の断面形状を計測でき、それに基づいて循環器指標を算出することができる循環器情報計測システムを提供することである。
また、本発明の他の目的は、非侵襲で血管の断面形状を計測でき、血管壁の脈動を考慮して正確な血圧を測定することができる循環器情報計測システムを提供することである。
【０００７】
【課題を解決するための手段】
上述した課題を解決し、目的を達成するために、本発明に係わる循環器情報計測システムは、動脈上の少なくとも２ヶ所の血管の断面の形状を同時に検出する断面形状検出手段と、該断面形状検出手段により検出された断面形状に基づいて、前記２ヶ所の断面の断面積を算出する断面積算出手段と、前記断面積算出手段により求められた断面積に基づいて血管内の血液の平均流速を算出する流速算出手段と、前記断面形状検出手段により検出された２ヶ所の血管の断面形状の時間変化から前記２ヶ所の断面の間の脈波の伝搬速度を算出する脈波速度算出手段と、前記流速算出手段により算出された血液の流速と、前記脈波速度算出手段により算出された脈波の伝搬速度とに基づいて、特定の時相における血圧値を基準とする相対血圧値を算出する相対血圧算出手段と、あらかじめ測定されている基準となる前記特定の時相における血圧値と前記相対血圧値とに基づいて、絶対血圧値を算出する絶対血圧算出手段とを具備し、前記流速算出手段は、前記２ヶ所の断面の断面積の時間微分と前記２ヶ所の断面の距離とから、前記２ヶ所の断面と血管壁により囲まれた部分の血液の平均流速を算出することを特徴としている。
【０００８】
また、この発明に係わる循環器情報計測システムにおいて、前記断面形状検出手段は、前記動脈に向けて超音波を放射する超音波放射手段と、生体内で反射された超音波を検出する超音波検出手段と、前記動脈の管壁の断面形状を抽出するために、前記超音波検出手段からの出力信号を適切な閾値を基準として２値化する２値化手段とを備えることを特徴としている。
【０００９】
また、この発明に係わる循環器情報計測システムにおいて、前記断面形状検出手段は、前記超音波検出手段からの出力信号を画像信号に変換する変換手段をさらに具備することを特徴としている。
また、この発明に係わる循環器情報計測システムにおいて、前記２値化手段は、前記超音波検出手段からの出力信号を前記変換手段により画像信号に変換した後に、適切な画像濃度値を閾値として２値化することを特徴としている。
【００１０】
また、この発明に係わる循環器情報計測システムにおいて、前記２値化手段は、超音波の入射方向と前記動脈の血管壁のなす角度に応じて前記閾値を設定する設定手段を備えることを特徴としている。
また、この発明に係わる循環器情報計測システムにおいて、前記２値化手段は、前記断面形状検出手段により検出された血管の断面形状と前記超音波の入射方向とに基づいて、該超音波の入射方向と前記血管壁のなす角度を演算する演算手段をさらに備えることを特徴としている。
【００１６】
【発明の実施の形態】
以下、本発明の実施形態について、添付図面を参照して詳細に説明する。
（第１の実施形態）
図１は本発明の循環器情報計測システムの第１の実施形態としての血圧測定装置の構成を示した図である。
【００１７】
図１に示すように、血管断層像収集部１は同一の血管２の近接した２点の超音波断層像３ａ，３ｂを収集するための超音波振動子４ａ，４ｂから構成される。収集された断層像４ａ，４ｂは血管断面積算出部５に伝達され、血管の断面の抽出処理がおこなわれる。血管断面積算出部５で算出された血管の断面積は、血流速算出部６と脈波伝搬速度算出部７に伝達される。血流速算出部６は、断面積と２点間の距離から、血液の流速を求める。脈波伝搬速度検出部７は、２点間の断面積の時間変化の位相差から脈波伝搬速度を求める。血流速算出部６と脈波伝搬速度算出部７によって算出された血流速と脈波伝搬速度が伝達されて、以下の換算式に代入される。
【００１８】
Ｐ−Ｐ0＝ρＣ0ｖ＋（１／８）ρｖ2
ここで、Ｐ：血圧、Ｐ0＝最低血圧（拡張期血圧）、ρ：血液の比重、Ｃ0：拡張期末期の脈波伝搬速度、ｖ：血管の断面内の平均流速である。
上記の式により、最低血圧を基準とする相対値が算出され、血圧変化算出部８に出力される。一方で基準圧測定部９で計られたある時点における血圧の絶対値が血圧算出部１０に送られ、血圧算出部１０において、血圧値が連続的に算出される。
【００１９】
図２に血管断層像収集部の構成を示す。超音波振動子群４aはa1からan、超音波振動子群４bはb1からbnのそれぞれｎ個（nは自然数）の超音波振動子で構成されている。超音波振動子群４aと超音波振動子群４ｂは距離Ｌだけ離れて平行に配置されている。まず、超音波振動子４aのa1から超音波が動脈２に向けて発射され、動脈壁などで反射した信号が同一の振動子によって検出される。この一連の動作をa1からanまで続けることにより、超音波断層像３aが得られる。動脈断層像３aが得られた後、動脈断層像３aを得たのと同様の操作をおこなうことにより動脈断層像４bを得る。
【００２０】
なお、ここでは簡単のため、a1,a2,・・・,anと１個ずつの超音波振動子が送受信をおこなって断層像を得ると説明したが、アレイプローブを用いる一般的な超音波診断装置と同様、複数個の振動子によるビームフォーミングによって断層像を得ることも可能である。ビームフォーミングとはアレイ内の位相もしくは遅延の分布を制御して目的の性状の超音波ビームを形成することである。
【００２１】
動脈断層像を収集する方法に関しては、動脈断層像が得られる限り、任意である。また、本実施形態では超音波振動子群４ａで断層像が得られた後に、超音波振動子群４ｂで断層像を得たが、断層像を得ることができる限り、断層像を同時に求めてもよいし、断層像の収集順序、振動子群の数、断層像の形状などの条件が任意であることはあきらかである。また、超音波振動子群を構成する超音波振動子についても、個数、配置の仕方、大きさ、発射される周波数などの条件に関しても、本発明の趣旨に反しない限り、任意の値に設定することが可能である。
【００２２】
得られた断層像４ａ，４ｂは、血管断面積算出部５に伝達される。図３は血管断面積算出部５の構成を示した図である。
伝達された断層像はフレームメモリ１１ａ，１１ｂに蓄えられている。この画像を２値化処理部１２で２値化することで、血管壁を抽出する。血管に超音波を照射した場合、血管壁と血管内部との境界部分では強い反射が起こる。一方、血管内は成分が一様であるため、音響インピーダンスの変化が少なく、反射波が弱い。よって、反射波の強いほど、白黒の高い階調値を割り当て、反射波の弱いものを、白黒の低い階調値を割り当てて画像化する場合、血管壁と血管内部との境界部分が白っぽく、血管内は階調が低い、つまり黒っぽい画像となる。この性質を利用して、ある階調値を閾値として、それ以上のものを白、つまり最高の階調値とし、それ以下のものを黒、つまり最低の階調値とする処理（２値化処理）を行なう。この２値化処理は、画像処理の分野で一般的に用いられている。
【００２３】
２値化処理によって血管内部が黒く描出された映像から、血管形状算出部１３で血管内部の面積、血管径を求め、血流速算出部６と、脈波伝搬速度算出部７に伝達する。なお、血管の断層像を時間経過に伴って連続的に求めることにより、血管の断面形状の時間変化を求めることができる。
図４の実線で示した図形のように、血流速算出部６では、超音波振動子群４ａ,４ｂによって血管の断面１４ａ,１４ｂの断面積（時間経過とともに変化する測定データ）Ｓa（ｔ）,Ｓb（ｔ）を求め、そのデータからＡ点、Ｂ点における断面内の平均流速ｖを、
【００２４】
【数１】
0003857788.tif 000002
【００２５】
ここで、ＬはＡ点とＢ点の間の距離
により求める。なお、記号の上の「・」は、時間微分を示す。
ここで、上式が有効であることを説明する。
動脈中での流れでは、血管の半径方向の速度勾配や圧力勾配は無視してさしつかえない。この場合には、流れとしては管軸方向の平均速度だけを考えればよい。ただし、血管の断面積（従って半径）の変化は考慮する。このようなモデルを準１次元モデルとよぶ。
【００２６】
図５に準１次元モデルの座標系を示す。ここで血管の管軸をｘ軸とする。また血管の断面積をＳ（ｘ,ｔ）とする。前述の測定点（２断面）の位置をｘ＝ｘ１,ｘ＝ｘ2とする（ｘ1＜ｘ2）。従って、Ｘ2が下流側である。また、ｘ＝０は、心臓の駆出部と考える（ただし、他の任意の点としてもよい）。
準１次元流れでは、連続の式
∂Ｓ／∂ｔ＋∂（ｖＳ）／∂ｘ＝０                    （１）
が成り立つ。（１）式から、
∂（ｖＳ）／∂ｘ＝−∂Ｓ／∂ｔ                      （２）
となる。（２）式を、ｘについて区間[ｘ,＋∞]で積分すると、
【００２７】
【数２】
0003857788.tif 000003
【００２８】
であり、かつ、
【００２９】
【数３】
0003857788.tif 000004
【００３０】
が有限の値であることから、
【００３１】
【数４】
0003857788.tif 000005
【００３２】
となる。
流量Ｑ（ｘ,ｔ）は、ｖ,Ｓの定義より、
Ｑ（ｘ,ｔ）＝ｖ（ｘ,ｔ）Ｓ（ｘ,ｔ）                  （４）
であるから、（３）式を代入して、
【００３３】
【数５】
0003857788.tif 000006
【００３４】
となる。（５）式より、ｘ＝ｘ1における流量Ｑ（ｘ1,ｔ）は、
【００３５】
【数６】
0003857788.tif 000007
【００３６】
となる。
Ｔを拍動の周期より短い時間として、測定したい時間が、区間[０,Ｔ]だとすると、今考えている流れでは、心臓が１回拍動して、次の拍動までの期間では単一波とみなせるので、
ｘ2≦ｘ＜∞  ,  ０≦ｔ≦Ｔ
において、Ｓは脈波が到達していないため時間変化せず、Ｓはｘのみの関数である。従って、（５）式の第２項は、
【００３７】
【数７】
0003857788.tif 000008
【００３８】
となる。よって、
【００３９】
【数８】
0003857788.tif 000009
【００４０】
となる。
区間[ｘ1,ｘ2]の間隔が十分小さいため、区間[ｘ1,ｘ2]において、Ｓはｘについて単調増加または単調減少であるから、
Ｓ（ｘ,ｔ）＝（ｘ＋ａ）ｆ（ｔ）                        （９）
ここで、ａは定数
と直線近似できる。従って、
【００４１】
【数９】
0003857788.tif 000010
【００４２】
となる。
（１０）式を（８）式に代入して、
【００４３】
【数１０】
0003857788.tif 000011
【００４４】
となる。（４）式より、
ｖ（ｘ1,ｔ）＝Ｑ（ｘ1,ｔ）／Ｓ（ｘ1,ｔ）
であるから、この式に（９）式、（１１）式を代入して、
【００４５】
【数１１】
0003857788.tif 000012
【００４６】
となる。
ｘ1,ｘ2に比べて、ｘ2−ｘ1は十分小さいので、
ｘ1＋ｘ2＋２ａ≒２（ｘ1＋ａ）                  （１３）
と近似できるため、（１２）式は、
【００４７】
【数１２】
0003857788.tif 000013
【００４８】
となる。
Ｌ＝ｘ2−ｘ1                                  （１５）
であるから、（１４）式は、
【００４９】
【数１３】
0003857788.tif 000014
【００５０】
となる。
（９）式より、
Ｓ（ｘ1,ｔ）＝（ｘ1＋ａ）ｆ（ｔ）                （１７）
Ｓ（ｘ2,ｔ）＝（ｘ2＋ａ）ｆ（ｔ）                （１８）
である。上式の差をとると、
Ｓ（ｘ2,ｔ）−Ｓ（ｘ1,ｔ）＝ｆ（ｔ）Ｌ                （１９）
また、
Ｓ（ｘ1,ｔ）＝Ｓa（ｔ）                          （２０）
Ｓ（ｘ2,ｔ）＝Ｓb（ｔ）                          （２１）
であるから、（２０）式、（２１）式を（１９）式に代入して、
Ｓb（ｔ）−Ｓa（ｔ）＝ｆ（ｔ）Ｌ                      （２２）
となる。（２２）式を両辺微分して、
【００５１】
【数１４】
0003857788.tif 000015
【００５２】
となる。（２２）式、（２３）式より、（１６）式は、
【００５３】
【数１５】
0003857788.tif 000016
【００５４】
となる。
（２４）式により、ｘ＝ｘ1における平均流速が求められた。また、ｘ＝ｘ1での流量Ｑ（ｘ1,ｔ）は、（４）式より、
【００５５】
【数１６】
0003857788.tif 000017
【００５６】
である。
図４に示すように実際の血管においては、点線１５ａ，１５ｂのような形状をとっている可能性があり、直線近似による誤差を減らす改良として、２断面間の距離を短くする方法、あるいは超音波振動子をマトリックスアレイとすることで正確な値を得る方法や、実際の曲線つまり、点線によく合致する曲線で側面部分を仮定し、積分をおこなうことなどの改良が容易に思いつくが、これらはすべて、本発明の範囲であることは明らかである。
【００５７】
脈波伝搬速度算出部７では、Ａ点、Ｂ点で生じる断面積の時間変化の時間遅れ（位相差）τで、 Ａ点、Ｂ点間の距離を割ることにより、
Ｃ0＝Ｌ／τ
として脈波伝搬速度Ｃ0を求める。
血圧変化算出部８では、理論式
Ｐ−Ｐ0＝ρＣ0ｖ＋（１／８）ρｖ2
に、伝達された流速、脈波伝搬速度を代入する。この式は、（菅原基晃、桜井泰久ら：“動脈中の脈波の非線形理論”、医用電子と生体工学、１１巻、３号、１８０−１８９、１９７３）に記載されている。
【００５８】
以上の処理により、基準圧測定部９で計られた時相の血圧値を元に血圧算出部１０で連続的に血圧を算出することが可能である。この基準圧の測定には一例として、カフ式自動血圧計（ある時間内の血圧を１回のみ測定する１点測定）などにより最低血圧をあらかじめ求めておいて設定する方法がある。また、基準圧は最低血圧に限らず、同一時刻での圧力値とそのときの血圧変化算出部出力が得られれば設定可能である。
【００５９】
算出された血圧値は使用者の必要に応じて、値を表示する液晶モニターなどの表示装置や、記録用の記憶媒体、プリンタなどに適宜出力することが可能となる。
なお、上記の説明では、生体内で反射された超音波を超音波振動子４ａ,４ｂが受信した信号を画像信号に変換した後に２値化するように説明したが、超音波振動子４ａ,４ｂから出力された受信信号を対数増幅し、検波した信号を２値化してもよい。
【００６０】
（第２の実施形態）
第１の実施形態では、血管断面積算出部５はある固定の閾値で２値化し、血管を２値画像として抽出したが、血管断面積の算出に関してはこの方法に限らない。
図６は超音波信号と血管の壁面１６の関係を示している。この図に示すように、血管の中心の直上にあるCm番目の超音波振動子から発せられる超音波は血管壁の接線に対して垂直に入射するため、同じ振動子Cmによって検出される反射波が大きい。しかしながら、直上から離れた場合には、振動子から発せられた超音波が血管壁の接線に対して斜めに入射するため、検出される反射波はCmに比べ弱くなってしまう。なお、複数個の振動子によるビームフォーミングをおこなう場合には、超音波ビームの中心軸が血管の中心に近い場合に反射波が大きく、血管の中心から、離れるほど反射波が小さくなる。
【００６１】
血管の中心の直上にあるCm番目の超音波振動子での反射波の強度が、図６（ｂ）で示す形状であるとする。血管壁に対応する部分がＹａ，Ｙｂであり、Ｙｋの部分がノイズであるとする。この場合、血管壁をうまく抽出するためには、閾値をノイズのピークより高く、血管壁のピークより低い、α2程度に設定するのが望ましい。
【００６２】
しかしながら、図６（ｃ）では超音波振動子が血管の中心の直上にないため、血管壁の部分の反射波の強度が弱く、上述の閾値α2に達しておらず、閾値がα1からα3の間の値でなければならないという問題が生じ、ノイズのある場合には、固定の閾値を用いた２値化ではうまく血管壁を描出できないことがある。
この問題を解決する方法として、本実施形態においては、超音波の血管壁の接線に対する入射角に応じて閾値を変化させる方法を示す。
【００６３】
図７は本実施形態における血管断面積算出部の概略構成を示す図である。血管直上検出部１８は、どの素子が血管の直上にきているかを検出する。
図８に血管直上検出部の構成を示す。図６の超音波振動子群３ｃにより得られた画像が、表示用のモニター２１に映し出される。モニターには超音波振動子群３ｃの中心あるいはその付近にある振動子より得られた画像が画面上のどの部分に表示されているかを指し示すマーカ２２が表示されており、操作者はこのマーカ２２に血管像の中心が一致するように、プローブを移動させる。中心からずれた血管像２３ａが、中心のマーカと一致する血管像２３ｂの位置に来るようにプローブの位置を調整し終えた後、操作者が検出終了確認ボタンを押すと、計測を開始する。、
換算部１９はフレームメモリ１７から呼び出されたデータが、超音波振動子群のどの振動子から得られたデータであるかを、フレームメモリ１７のアドレスより計算し、閾値変化２値化部２０に伝達する。閾値変化２値化部２０は上述の換算部から出力されるデータにより、閾値を変化させて２値化する。
【００６４】
閾値変化２値化部２０により設定される閾値は、必ずしも血管壁に対する超音波信号の入射角のみによって決められるものではない。実際には、測定条件や被験者などによっても変化する。よって、閾値変化に適当な補正項をつけることは本発明の趣旨に反するものではない。
抽出した血管壁から、血管の断面積Ｓ、血管内半径ｒ、あるいはその時間変化から血管状態などを求めることが可能である。
【００６５】
血管壁のヤング率を反映し、動脈硬化の指標とされる脈波伝搬速度を求める方法を示す。図９は脈波伝搬速度を求める際の構成図（ａ）及び得られる波形（ｂ）を説明する図である。超音波振動子群を距離Ｌだけ離して、動脈上に配置する。それぞれの超音波振動子群から得られた断層像から本実施形態の上記の手法により、断面積を求める。得られた波形は、図９（ｂ）に示すように、脈動を反映した波形となる。２つの断層像３ａ，３ｂにみられる断面積変化の時間遅れ（位相差）をτとすると、脈波伝搬速度ＣはＣ＝Ｌ／τで求めることができる。通常、脈波伝搬速度は頚動脈波と股動脈波の時間差により、求められるが、頚動脈と股動脈の間に存在するすべての血管の硬化度を反映したものである。従って、本方法によって、局所的な脈波伝搬速度を計測することが可能となり、動脈硬化の進行の進んだ部位の特定などのより詳細な診断が可能となる。
【００６６】
また、図９（ｂ）の断層像３ａの断面積の変化の周期τ1から、一分あたりに換算した脈拍数ＲＲをＲＲ＝６０／τ1という換算式で求めることができる。
さらに、２つの断層像３ａ，３ｂから血流量や断面内平均流速を求めることが可能である。これらを求める方法は第１の実施形態で既に述べた。
また、血圧計測技術を組み合わせると、さらに直接的に血管の弾性率を求めることが可能となり、動脈硬化度を評価する上で有効な計測が可能となる。動脈硬化の指標として、圧力弾性率Epが、
Ep＝ΔP／（Δｒ／ｒ ）
で定義されている。この指標を特開平４−３２９９３８号公報では、血管内に挿入するカテーテルで圧力、血管径を求めているが、本実施形態に記載された方法を用いれば、非侵襲に計測することができる。
上記の脈拍、脈波伝搬速度、圧力弾性率以外にも、血管壁の厚み、血管壁の脈動速度など、抽出した血管壁の断面形状の時間変化から様々な循環器の情報を得ることが可能であると考えられるが、これらはすべて本発明の範疇であることは明らかである。
【００６７】
（第３の実施形態）
第２の実施形態においては、血管の中心の直上を検出する直上検出を、操作者が手動でおこなう方法を示したが、これ以外の方法で血管の直上を特定することも可能である。本実施形態では、感圧素子を用いた直上検出の方法を示す。
図１０は、血管２の中心を求めるための感圧素子群の概略構成を示す図である。感圧素子群３１をカフ等の押圧部で押し込むことで、血管を圧迫し、図のように平坦に変形させることで、圧力値が感圧素子群３１に伝わる。伝わる圧力の最大を示す部分を血管の直上とする。
【００６８】
ここでは、伝わる圧力の最大値をとる部分を直上としたが、これは圧迫法、圧迫強度、血管の位置などによっても異なるので、直上の判断を示す方法を条件に応じて改良することは、任意である。
図１１は上記の直上検出手段の構成の一例を示す図である。超音波振動子群４ｅと感圧素子群３１を平行に配置している。まず、血管像を計測する前に、感圧素子群３１を適度な圧力で生体に押しつけ、圧波形をもとに感圧素子群３１のどの圧力素子が、血管の直上にあるかを判定する。この直上の判定は血管の脈動によって時間変化する各圧力素子の圧波形を比較し、もっとも大きくなる感圧素子を直上とする。直上であるとされた感圧素子に近接している超音波振動子を、直上の超音波振動子と判断する。直上の検出が終了した後は、血管の扁平をさけるために圧迫をゆるめることも可能である。
【００６９】
第２の実施形態においては、超音波の入射角により閾値を変化させるとしたが、これは必ずしも、入射角のみ（たとえば入射角の正弦との比例関係でのみ）で閾値が決定されると限定しているものではない。各振動子から得られる超音波の映像に対して、逐次、入射角を計算することが機器の構成上煩雑で困難である場合には、血管の直上の超音波振動子により受信された超音波信号からなる映像部分を入射角０度と設定し、これを基準として、その振動子から離れるに従って、直線的、曲線的など任意の波形で閾値を変化させるという方法をとることも可能である。つまり、１あるいは数箇所の入射角の情報のみを利用するという簡易的な方法を用いることも、本発明の範囲内である。
【００７０】
図１２は閾値の設定の例を示した図である。図１２（ａ）では直上の血管部分の閾値をβ１とし、直線状に閾値を変化させ、直上からｒaだけ離れたところの閾値がβ２となるように閾値を設定している。しかしながら、同じ血管においても脈動がある場合には、拡張期と収縮期においては、設定すべき閾値を変えるのが適当である。従って、図１２（ｂ）では３４ａで示すように血管の拡張している場合には実線のように、３４ｂで示すように収縮している場合には点線のようにと、脈動に応じて閾値を変化させる方法をとっている。具体的には、心電図、心音などのデータから、拡張、収縮期を知る拡張・収縮検出部を設けて、そのデータにより閾値を変化させる方法が考えられる。
【００７１】
なお、上記のように閾値を入射角に応じて直線状あるいは曲線状に変化させ、且つ、この直線状あるいは曲線状の閾値の変化曲線全体を、ボリューム等により平行に上下動できるようにしておき、操作者が画面上の血管断面の２値化画像を見ながら、血管壁がもっとも鮮明に見えるようにボリュームを調整するようにしても良い。
【００７２】
以上の実施形態に示したように、血管の状態に応じて適切な閾値を設定することにより、動脈の血管壁を抽出し、血管壁の映像から種々の循環情報の計測をおこなうことが可能となる。
なお、本発明による装置を、他の装置の一部として利用することは、本発明の趣旨に反しない限り、任意である。たとえば、超音波診断装置に本発明による装置を連動させ、抽出した血管壁を超音波像に重ねてモニタ表示したり、得られた数値をモニタの画面や用紙に出力することは本発明の範囲であることは明らかである。
【００７３】
【発明の効果】
以上説明したように、本発明によれば、血管壁の脈動による管壁の動きを考慮した理論式をもちいて、血圧を連続的かつ非侵襲に測定することができる。
また、超音波断層像から血管壁部分を抽出できる構成を用いることにより、非侵襲に血管、血流状態などの循環器指標を計測できる循環器情報計測システムを提供することができる。
【図面の簡単な説明】
【図１】本発明の循環器情報計測システムの第１の実施形態としての血圧測定装置の構成を示した図である。
【図２】血管断層像収集部の構成を示す図である。
【図３】血管断面積算出部の構成を示した図である。
【図４】血管の一部を拡大して示した模式図である。
【図５】準１次元モデルの座標系を示す図である。
【図６】超音波信号と血管の壁面の関係を示した図である。
【図７】第２の実施形態における血管断面積算出部の概略構成を示す図である。
【図８】血管直上検出部の構成を示す図である。
【図９】脈波伝搬速度を求めるための構成を示す図である。
【図１０】血管の中心を求めるための感圧素子群の概略構成を示す図である。
【図１１】血管の直上検出手段の構成の一例を示す図である。
【図１２】閾値の設定の例を示した図である。
【符号の説明】
１  血管断層像収集部
２  血管
３ａ,３ｂ  超音波断層像
４ａ,４ｂ  超音波振動子群
【図１】
0003857788.tif 000018
【図２】
0003857788.tif 000019
【図３】
0003857788.tif 000020
【図４】
0003857788.tif 000021
【図５】
0003857788.tif 000022
【図６】
0003857788.tif 000023
【図７】
0003857788.tif 000024
【図８】
0003857788.tif 000025
【図９】
0003857788.tif 000026
【図１０】
0003857788.tif 000027
【図１１】
0003857788.tif 000028
【図１２】
0003857788.tif 000029