(57)【特許請求の範囲】
【請求項１】
生体に超音波センサを接触させて該生体内の血管の血管径を連続的に測定する生体の血管径連続測定装置であって、
前記超音波センサにより検出された超音波反射信号に基づいて前記血管の短軸画像を繰り返し生成する短軸画像生成部と、
該短軸画像生成部により繰り返し生成された短軸画像間における複数の特徴点領域の変位を、該特徴点領域のパターンと一致する領域を検出するパターンマッチングを用いてそれぞれ逐次算出する変位算出部と、
該変位算出部により算出された前記特徴点領域の血管径方向成分の変位から前記血管の径変位量を算出する血管径変位量算出部と、
前記血管径変移量算出部により逐次求められることで脈動を示す前記血管の径変位量を微分することで周期的に拡張末期を示すスパイク状の拡張末期心拍同期信号を発生する拡張末期タイミング信号発生部と、
前記血管径変位量算出部により算出された前記血管の径変位量のうちの前記スパイク状の拡張末期心拍同期信号の発生時の径変位量に基づいて血管径の時系列的変化を測定する血管径測定部と、
を、含むことを特徴とする血管径連続測定装置。
【請求項２】
前記超音波センサは前記血管に直交するように配置される短軸用超音波アレイ探触子と該血管に平行な方向に配置される長軸用超音波アレイ探触子と、前記長軸用超音波アレイ探触子からの超音波反射信号に基づいて長軸での血管径を逐次算出する長軸血管径算出部とを含み、
前記短軸画像生成部は、該短軸用超音波アレイ探触子により検出された超音波反射信号に基づいて前記血管の短軸画像を繰り返し生成するものであり、
前記血管径測定部は、前記長軸血管径算出部により算出された血管径の時系列的変化を血管径測定値として用いるが、該血管径の時系列的変化が得られない空白区間がある場合またはすべてが空白区間の場合には、前記血管径変位量算出部により前記短軸画像中の特徴点領域の血管径方向成分の変位から算出された前記血管の径変位量を用いて該空白区間内の血管径の時系列的変化を生成する
ことを特徴とする請求項１の血管径連続測定装置。
【請求項３】
前記長軸用超音波アレイ探触子からの超音波反射信号に基づいて該超音波反射信号の検波波形のピーク値を結ぶ包絡線信号であるＡモード波形を生成するＡモード波形生成部を含み、
前記長軸血管径算出部は、該Ａモード波形生成部から生成されたＡモード波形が示すピーク波形間の時間差から長軸での血管径を逐次算出するものである
ことを特徴とする請求項２の血管径連続測定装置。
【請求項４】
前記長軸血管径算出部は、該Ａモード波形生成部から生成されたＡモード波形が示すピーク波形間の時間差のうち、前記拡張末期タイミング信号発生部からの拡張末期心拍同期信号に同期した時間差から長軸での血管径を逐次算出するものである
ことを特徴とする請求項３の血管径連続測定装置。
【請求項５】
前記血管径変位量算出部は、オプティカルフロー解析を用いて、前記短軸画像生成部により繰り返し生成された短軸画像間における複数の特徴点領域の変位を、該特徴点領域のパターンと一致する領域を検出することによりそれぞれ逐次算出する
ことを特徴とする請求項１乃至４のいずれか１の血管径連続測定装置。
【請求項６】
前記超音波センサは前記血管に直交するように配置される短軸用超音波アレイ探触子と該血管に平行な方向に配置される長軸用超音波アレイ探触子と、前記長軸用超音波アレイ探触子からの超音波反射信号に基づいて長軸での血管径を逐次算出する長軸血管径算出部とを含み、
前記短軸画像生成部は、該短軸用超音波アレイ探触子により検出された超音波反射信号に基づいて前記血管の短軸画像を繰り返し生成するものであり、
前記血管径測定部は、前記長軸血管径算出部により算出された血管径の時系列的変化を血管径測定値として用いるが、該血管径の時系列的変化が得られない空白区間がある場合またはすべてが空白区間の場合には、前記血管径変位量算出部により前記短軸画像中の特徴点領域の血管径方向成分の変位から算出された前記血管の径変位量を用いて該空白区間内の血管径の時系列的変化を生成する
ことを特徴とする請求項１又は２の血管径連続測定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、生体に超音波センサを接触させてその生体内に位置する血管の血管径を連続的に測定する生体の血管径連続測定装置に関するものである。
【背景技術】
【０００２】
生体に超音波センサを接触させてそれから超音波放射信号を発信し、その生体内の血管の血管径、内腔径などの血管径を連続的に測定することが行われている。この血管径の測定は、血管と他の組織との伝播速度差によりそれらの境界から反射される超音波反射信号間の時間差処理、或いはその反射信号から合成される超音波画像上における距離測定などにより実行される。【０００３】
上記生体内に位置する血管径の測定を正確に行なうためには、生体たとえば皮膚上に接触させた超音波センサに対して血管の相対移動を防止する必要があるため、その生体の一部たとえば上肢を所定の姿勢で固定できる載置台あるいは上肢保持装置が用いられる。たとえば、特許文献１に示す血管内皮機能検査装置の上肢保持装置がそれである。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００７−６１１８２号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、上記のような血管径の連続測定においては、一般に、血管の長手方向に沿って位置させた超音波アレイを用いて得られた、検波後の血管からの超音波反射信号のピーク値を結ぶ包絡線信号であるＡモード波形あるいは超音波反射信号を明暗として表わしたＢモード画像が用いられ、そのＡモード波形のピーク値又はその立上りの位置間あるいはＢモード画像における血管のエッジ位置間が示す時間差と音速とに基づいて血管径が逐次計測される。血管の超音波の反射点が線状である長軸断面での測定は、血管の超音波の反射点が点状である短軸断面の測定に対して、計測可能な点が多くとれるためである。しかし、安定に測定できる状態が保持できる前記上肢保持装置のような生体の一部を固定する装置を用いたとしても、被測定者の体動が発生したりすると、超音波センサと血管との相対位置にずれが発生し、少なくとも血管径の連続測定が不可能となる区間が発生する。仮に、超音波センサの位置を血管に対して自動的に位置制御を行なったとしても、体動発生時点から正常位置への復帰までは計測が不可能となる。【０００６】
本発明は以上の事情を背景としてなされたものであり、その目的とするところは、血管に対する超音波センサのずれが発生しても血管径の連続測定が可能な生体の血管径連続測定装置を提供することにある。【０００７】
本発明者等は、上記事情を背景として種々検討を重ねた結果、超音波センサと血管との相対位置がずれて長軸画像が得られなかったとしても断面画像内に得られる血管短軸画像において、オプティカルフロー解析を用いてその画像中の特徴点をパターンマッチングにて追跡するとその特徴点の移動量を一拍毎に計測でき、その特徴点の移動量から血管径の変化量を逐次求めることができることを見いだした。本発明はこの知見に基づいて為されたものである。
【課題を解決するための手段】
【０００８】
すなわち、上記目的を達成するための請求項１に係る発明の生体の血管径連続測定装置は、（ａ）生体に超音波センサを接触させて該生体内の血管の血管径を連続的に測定する生体の血管径連続測定装置であって、（ｂ）前記超音波センサにより検出された超音波反射信号に基づいて前記血管の短軸画像を繰り返し生成する短軸画像生成部と、（ｃ）該短軸画像生成部により繰り返し生成された短軸画像間における複数の特徴点領域の変位を、該特徴点領域のパターンと一致する領域を検出するパターンマッチングを用いてそれぞれ逐次算出する変位算出部と、（ｄ）該変位算出部により算出された前記特徴点領域の血管径方向成分の変位から前記血管の径変位量を算出する血管径変位量算出部と、（ｅ）前記血管径変移量算出部により逐次求められることで脈動を示す前記血管の径変位量を微分することで周期的に拡張末期を示すスパイク状の拡張末期心拍同期信号を発生する拡張末期タイミング信号発生部と、（ｆ）前記血管径変位量算出部により算出された前記血管の径変位量のうちの前記スパイク状の拡張末期同期信号の発生時の前記血管の径変位量に基づいて血管径の時系列的変化を測定する血管径測定部と、を、含むことを特徴とする。
【発明の効果】
【０００９】
本発明の生体の血管径連続測定装置によれば、超音波反射信号に基づいて繰り返し生成された短軸画像間における複数の特徴点領域の変位がパターンマッチングにより逐次求められ、それら複数の特徴点領域の血管径方向成分の変位から前記血管の径変位量が逐次算出され、その血管の径変位量に基づいて血管径の時系列的変化が測定される。このため、体動などによる超音波センサのずれに影響され難い短軸画像から生体の血管の径の時系列的変化が得られる。また、本発明の生体の血管径連続測定装置によれば、血管径変位量算出部は、前記変位算出部により算出された前記特徴点領域の変位のうち前記スパイク状の拡張末期心拍同期信号に同期した血管径方向成分の変位から前記血管の径変位量を算出し、血管径測定部は、上記拡張末期心拍同期信号に同期して１拍毎の径変化量を決定し、その１拍毎の径変化量から前記血管径を時系列的に決定するため、拡張末期タイミング信号としてＥＣＧ信号を検出する心電誘導信号検出装置を用いることが不要となる利点がある。
【００１１】
ここで、好適には、前記超音波センサは前記血管に直交するように配置される短軸用超音波アレイ探触子と該血管に平行な方向に配置される長軸用超音波アレイ探触子と、前記長軸用超音波アレイ探触子からの超音波反射信号に基づいて長軸での血管径を逐次算出する長軸血管径算出部とを含み、前記短軸画像生成部は、該短軸用超音波アレイ探触子により検出された超音波反射信号に基づいて前記血管の短軸画像を繰り返し生成するものであり、前記血管径測定部は、該長軸血管径算出部により算出された血管径の時系列的変化を血管径測定値として用いるが、該血管径の時系列的変化が得られない空白区間がある場合またはすべてが空白区間の場合には、前記血管径変位量算出部により前記短軸画像中の特徴点領域の血管径方向成分の変位から算出された前記血管の径変位量を用いて該空白区間内の血管径の時系列的変化を生成するので、体動などによる長軸用超音波アレイ探触子のずれにより血管を見失った場合でも、時系列的に連続した血管の径変化量を測定できる。
【００１２】
また、好適には、前記長軸用超音波アレイ探触子からの超音波反射信号に基づいて該超音波反射信号の検波波形のピーク値を結ぶ包絡線信号であるＡモード波形を生成するＡモード波形生成部を含み、前記長軸血管径算出部は、該Ａモード波形生成部から生成されたＡモード波形が示すピーク波形間の時間差から長軸での血管径を逐次算出するものである。このようにすれば、長軸画像から血管径を測定する場合に比較して最小計測単位が１画素とならないため、長軸での血管径の測定精度が高められる。【００１３】
また、好適には、前記長軸血管径算出部は、該Ａモード波形生成部から生成されたＡモード波形が示すピーク波形間の時間差のうち、前記拡張末期タイミング信号発生部からの拡張末期心拍同期信号に同期した時間差から長軸での血管径を逐次算出するものである。このようにすれば、長軸での血管径の算出のために、拡張末期タイミング信号としてＥＣＧ信号を検出する心電誘導信号検出装置を用いることが不要となる利点がある。【００１４】
また、好適には、前記血管径変位量算出部は、オプティカルフロー解析を用いて、前記短軸画像生成部により繰り返し生成された短軸画像間における複数の特徴点領域の変位を、該特徴点領域のパターンと一致する領域を検出することによりそれぞれ逐次算出するものである。このようにすれば、オプティカルフロー解析が用いられるので、高い精度で速度ベクトルが得られ、得られる特徴点間の変位の精度が高められる。
【図面の簡単な説明】
【００１５】
【図１】本発明の一実施例である生体の血管径連続測定装置の全体的構成を説明する図である。
【図２】図１の生体の血管径連即測定装置の超音波プローブおよび制御装置の制御機能を説明する図である。
【図３】図１の生体の血管径連続測定装置において、位置決め装置により位置決めされる超音波プローブと血管の関係を示す図である。
【図４】超音波プローブの第１短軸用超音波アレイ探触子、長軸用超音波アレイ探触子、第２短軸用超音波アレイ探触子によりそれぞれ検出された超音波反射信号からそれぞれ生成された超音波画像が図２の画像表示装置に表示された例であって、超音波プローブにずれがない場合を示す図である。
【図５】図２の電子制御装置によって測定され、画像表示装置に表示される血管径の連続測定値を示グラフを示す図である。
【図６】図４と同様に、図２の画像表示装置に表示される超音波画像の例であって、体動により超音波プローブにずれが生じた場合を示す図である。
【図７】図２の電子制御装置によって逐次測定され、画像表示装置に表示される血管径の連続測定値を示グラフであって、図６に示す超音波プローブにずれ期間に対応する空白期間を示す図である。
【図８】血管変位を算出するための特徴点および特徴点領域を短軸画像中に設定した例を説明する図である。
【図９】図８の特徴点領域の変位から求められた血管径の変化量を示す図であって、図７の空白区間の補完するための情報を例示する図である。
【図１０】図９に示す血管径の変化量の微分に相当する変化波形を、拡張末期タイミング信号Ｍと共に示す図である。
【図１１】図２の電子制御装置の制御作動の要部を説明するフローチャートである。
【図１２】本発明の他の実施例の電子制御装置の制御機能の要部を説明する機能ブロック線図であって、図２に相当する図である。
【図１３】図１２の電子制御装置の制御作動の要部を説明するフローチャートである。
【発明を実施するための形態】
【００１６】
以下、本発明の一実施例を図面を参照して詳細に説明する。なお、以下の実施例において図は適宜簡略化或いは変形されており、各部の寸法比および形状等は必ずしも正確に描かれていない。
【実施例１】
【００１７】
図１は、本発明の一実施例である生体の血管径連続測定装置１０の全体を示す図である。図２は、図１の血管径連続測定装置１０の全体的な構成、すなわちセンサ保持器２２、それに保持されたハイブリッドプローブユニット２４、および電子制御装置３８の制御機能を主に説明する図である。この血管径連続測定装置１０は、上肢１１に皮膚上に超音波センサとして機能するハイブリッドプローブユニット２４を配置して皮膚下の動脈などの血管１２の血管径、内膜厚を連続的に測定するものであり、ベッド１４上に仰臥になった被測定者（生体）１６の側方へ突き出した上肢１１が載置され測定テーブル２０上に設置されて、血管内皮機能検査装置としても機能する。上記血管径連続測定装置１０は、センサ保持器２２に保持されたハイブリッドプローブユニット２４を用いて上肢１１の皮膚２６の上からその皮膚２６直下に位置する血管１２の横断面画像（短軸画像）或いは縦断面画像（長軸画像）を測定する。【００１８】
ハイブリッドプローブユニット２４は、血管１２に関連する生体情報すなわち血管パラメータを検出するための超音波センサとして機能するものであって、図３に示すように互いに平行な２列の第１短軸用超音波アレイ探触子Ａおよび第２短軸用超音波アレイ探触子Ｂとそれらの長手方向中央部を連結する長軸用超音波アレイ探触子Ｃとを一平面すなわち平坦な探触面上に有して成るＨ型の超音波プローブ３２と、その超音波プローブ３２を位置決めするための多軸駆動装置（位置決め装置）３４とを備えている。図３は、超音波プローブ３２と血管１２の関係を示す図である。上記第１短軸用超音波アレイ探触子Ａ、第２短軸用超音波アレイ探触子Ｂ、および長軸用超音波アレイ探触子Ｃは、たとえば圧電セラミックスから構成された多数個の超音波振動子（超音波発振子）が直線的に配列されることにより長手状にそれぞれ構成されている。【００１９】
図２に示す、血管径連続測定装置１０は、所謂マイクロコンピュータから構成された電子制御装置３８と、画像表示装置４０と、Ｈ型の超音波プローブ３２の各探触子Ａ、Ｂ、Ｃとに駆動信号を供給し、超音波プローブ３２で検出された超音波反射信号を出力する超音波受信送信回路４２と、駆動モータ制御回路４４とを備えている。上記電子制御装置３８によって超音波受信送信回路４２から駆動信号が供給されると、ハイブリッドプローブユニット２４の超音波プローブ３２の各探触子Ａ、Ｂ、Ｃから超音波が放射される。ついで、その各探触子Ａ、Ｂ、Ｃにより検知された超音波反射信号を受けてその超音波反射信号の処理が行われることによって、皮膚２６下の超音波画像が発生させられ画像表示装置４０に図４に示すように表示されるとともに、連続測定された血管径の変化が図５に示すように表示される。ＦＭＤ（血流依存性血管拡張反応）の評価に際しては、血管径としては血管１２の内膜の径(内口径)が用いられる。【００２０】
画像表示装置４０は、図４に示すように、第１短軸用超音波アレイ探触子Ａによる超音波画像を表示する第１短軸画像表示領域Ｓ１と、第２短軸用超音波アレイ探触子Ｂによる超音波画像を表示する第２短軸画像表示領域Ｓ２と、長軸用超音波アレイ探触子Ｃによる超音波画像を表示する長軸画像表示領域Ｓ３とを有している。第１短軸画像表示領域Ｓ１および第２短軸画像表示領域Ｓ２および長軸画像表示領域Ｓ３は、皮膚２６からの深さ寸法を示す共通の縦軸を備えたものである。血管１２の超音波画像が生成されるに際して、超音波プローブ３２は、血管１２に対して所定の位置となるよう電子制御装置３８によって駆動モータ制御回路４４から駆動信号を供給された多軸駆動装置３４が駆動することにより位置決めさせられる。上記所定の位置とは、図３に示すように第１短軸用超音波アレイ探触子Ａおよび第２短軸用超音波アレイ探触子Ｂが血管１２に対して直交する位置し、長軸用超音波アレイ探触子Ｃが血管１２に対して平行且つ真上となる位置である。【００２１】
センサ保持器２２(図２）は、三次元空間内の所望の位置すなわち所定の位置において生体１６の上腕３６の皮膚２６の上からその皮膚２６直下に位置する血管１２を変形させない程度に軽く接触させる状態でハイブリッドプローブユニット２４を所望の姿勢で保持する。上記ハイブリッドプローブユニット２４の超音波プローブ３２の端面と皮膚２６との間には、通常、超音波の減衰、境界面における反射や散乱を抑制して超音波画像を明瞭とするためのよく知られたゼリー（超音波ゼリー）等のカップリング剤が介在させられる。このゼリーは、たとえば寒天等の高い割合で水を含むゲル状の吸水性高分子であって、空気よりは固有インピーダンス（＝音速×密度）が十分に高く大きく超音波送受信信号の減衰を抑制するものである。また、そのゼリーに換えて、水を樹脂製袋内に閉じ込めた水袋、オリーブ油、グリセリン等が用いられ得る。【００２２】
上記センサ保持器２２は、たとえば磁気的吸着力により机、台座等に固定されるマグネット台４８と、前記ハイブリッドプローブユニット２４が固定されるユニット固定具５０と、マグネット台４８およびユニット固定具５０に一端が固定され且つ球状に形成された先端部５２を備えた連結部材５４、５６と、それら連結部材５４、５６を介して、マグネット台４８とユニット固定具５０とを相対移動可能に連結し支持する自在アーム５８とを備えている。上記自在アーム５８は、相互に回動可能に連結された２つのリンク６０、６２と、そのリンク６０、６２の一端にて前記各先端部５２に所定の抵抗が付勢されつつその先端部５２に対して回曲可能に嵌め入れられた勘合穴６４をそれぞれ有する回曲関節部６６、６８と、各リンク６０、６２の他端にてその他端を相互に相対回動可能に連結し且
つその連結箇所を貫設するねじ穴に螺合されたおねじ付き固定ノブ７０が締め付けられることで得られる締着力により相対回動不能にされる回動関節部７２とを、有する。【００２３】
血管パラメータの測定に際しては、図２において、超音波駆動制御部８０は、電子制御装置３８からの指令に従って、たとえば第１短軸用超音波アレイ探触子Ａを構成する一列に配列された多数個の超音波振動子のうち、その端の超音波振動子ａ１ら、一定数の超音波振動子群たとえば１５個の毎に所定の位相差を付与しつつ１０ＭＨｚ程度の周波数で同時駆動するビームフォーミング駆動することにより超音波振動子の配列方向において収束性の超音波ビームを血管１２に向かって順次放射させ、超音波振動子を１個ずつずらしながらその超音波ビームをスキャン（走査）させたときの放射毎の反射波を受信して電子制御装置３８へ入力させる。また、上記第１短軸用超音波アレイ探触子Ａの放射面には、その超音波振動子の配列方向に直交する方向に超音波ビームを収束させるための図示しない音響レンズが設けられている。上記のようなビームフォーミング駆動および音響レンズによって収束性とされた超音波ビームには、超音波振動子の配列方向に対して直交する方向に長手状の収束断面が形成される。この収束断面の長手方向は、平面視において超音波振動子の配列方向、およびビームの放射方向に対して、それぞれ直交する方向である。電子制御装置３８は、上記反射波に基づいて画像を合成し、皮膚２６下における血管１２の横断面画像（短軸画像）、あるいは縦断面画像（長軸画像）を生成させて、画像表示装置４０に表示させる。【００２４】
電子制御装置３８には、超音波受信送信回路４２に一定の周期で駆動指令を繰り返し出力する超音波制御部８０と、超音波受信送信回路４２からの超音波反射信号を検波処理する検波処理部８２と、長軸用超音波アレイ探触子Ｃからの超音波反射信号の検波波形のピーク値を結ぶ包絡線信号であるＡモード波形を生成するＡモード信号処理部８４と、Ａモード信号処理部８４により生成されたＡモード波形に血管壁位置に対応して表れるピーク波形間の時間差から長軸血管径を算出する長軸血管径算出部８６と、Ｂモード信号処理部として機能し、超音波反射信号に基づいて前記血管１２の短軸画像を繰り返し生成する短軸画像生成部８８と、その短軸画像生成部８８により繰り返し生成された短軸画像間における複数の特徴点Ｐを含む特徴点領域Ｐの変位を、オプティカルフロー解析を用いてその特徴点領域Ｒのパターンと一致する領域を検出するパターンマッチングを用いてそれぞれ逐次算出する短軸画像血管断面のオプティカルフローによる変位算出部９０と、その変位算出部９０により算出された前記特徴点領域Ｒの変位のうち血管径方向成分の変位から血管１２の径変位量を算出する血管径変位量算出部９２と、その血管径変位量算出部９２により算出された血管１２の径変位量と長軸血管径算出部８６により算出された長軸血管径とに基づいて、血管径補完情報である短軸画像に由来する補完用血管径を算出する血管径補完情報処理部９４と、常時は、長軸血管径算出部８６により算出された長軸血管径の時系列的変化を血管径測定値として用いるが、図６に示すように体動による超音波センサの移動によって血管１２を見失って追従不可能となり図７に示すように長軸血管径の時系列的変化が得られない空白区間が形成される場合には、血管径補完情報処理部９４により生成された補完用血管径を上記空白区間の血管径として用いて血管径の時系列的変化を生成(測定)し、記録する血管径測定部９６と、その血管径測定部９６により測定された血管径(絶対値)の時系列的変化を画像表示装置４０に図５に示すように表示させる表示制御部９８と、第１短軸用超音波アレイ探触子Ａおよび第２短軸用超音波アレイ探触子Ｂが血管１２と直交し且つ長軸用超音波アレイ探触子Ｃが血管１２の真上でそれと平行となるように、たとえば表示制御部９８で表示される図４に示す短軸画像および長軸画像に基づいて制御する駆動モータ制御部１００とを、機能的に備えている。【００２５】
変位算出部９０において短軸画像間における特徴点領域Ｒ間の変位を導くオプティカルフロー解析では、各速度ベクトルの血管断面中心方向の成分は血管壁の拍動や拡大、縮小減少による移動を示している。このうち対峙する(たとえば血管壁の近位、遠位)のベクトル成分の差が血管径の変化量となる。ある時点の血管径にこの変化量を積算していくと、拍動による変化と拡大、縮小減少による変化が加わった血管径の経時的変化として計測される。図９の実線は、脈拍に同期して脈動する血管径の変化量を示し、破線は上記拡張末期のタイミングで特定された径変化量を示している。この破線が、前記空白区間の血管径の変化量として用いられる。【００２６】
電子制御装置３８の機能に含まれる拡張末期タイミング発生部１０２は、変位算出部９０により算出される、図９の実線の血管径の変化量を微分することで拡張末期のスパイク状の拡張末期タイミング信号Ｍを拡張末期心拍同期信号として生成する。図１０は、図９の実線に示される血管径の変化量の微分波形を示している。この拡張末期タイミング信号Ｍは、ＥＣＧ信号のＲ波と同様に、明瞭に表れるので信頼性があり、容易に利用され得る。【００２７】
この拡張末期タイミング信号Ｍは、長軸血管径算出部８６において長軸血管径算出の基礎となる、Ａモード信号処理部８４により生成されたＡモード波形のピーク波形間の周期的に脈動する時間差のうちから拡張末期の時間差を取り出す拡張末期タイミング信号として用いられ、取り出された時間差から長軸血管径が算出される。また、血管径変位量算出部９２は、変位算出部９０により算出された特徴点領域Ｒの変位のうち上記拡張末期タイミング信号Ｍに同期した変位から血管１２の径変化量を算出するので、血管径測定部９６では、上記心拍同期信号に同期して１拍毎の径変化量が決定され、その１拍毎の径変化量から血管径が時系列的に決定され、記憶される。このため、拡張末期タイミング信号としてＥＣＧ信号を検出する心電誘導信号検出装置を用いることが不要となる。【００２８】
電子制御装置３８では、上述のようにして得られた血管径連続測定値から、たとえば血管内皮機能を評価するために、腕に巻回された図示しない圧迫帯を用いて阻血虚血反応性充血後のＦＭＤ（血流依存性血管拡張反応）を表す血管径の変化率（％）［＝１００×（ｄmax −ｄ）／ｄ］（但し、ｄは安静時の血管径、ｄmax は阻血解放後の最大血管径）を算出する。【００２９】
図１１は、電子制御装置３８の制御作動の要部を説明するフローチャートである。図１１のステップＳ１( 以下、ステップを省略する)では、圧迫帯による阻血前から阻血後の一定時間経過後までの区間内の超音波反射信号が読み込まれる。ついで、Ａモード信号処理部８４に対応するＳ２では、長軸用超音波アレイ探触子Ｃからの超音波反射信号の検波波形のピーク値を結ぶ包絡線信号であるＡモード波形が生成される。長軸血管径算出部８６に対応するＳ３では、Ｓ２により生成されたＡモード波形に血管壁位置に対応して表れるピーク波形間の時間差から長軸血管径が算出される。次に、短軸画像生成部８８に対応するＳ４では、
Ｂモード信号処理が実行され、短軸用超音波アレイ探触子ＡあるいはＢからの超音波反射信号に基づいて血管の短軸画像が繰り返し生成される。そして、変位算出部９０に対応するＳ５では、Ｓ４により繰り返し生成された複数の短軸画像間における複数の特徴点領域Ｒの変位が、その特徴点領域Ｒのパターンと一致する領域を検出するパターンマッチングを用いてそれぞれ逐次算出される。ついで、血管径変位量算出部９２に対応するＳ６では、Ｓ５により算出された特徴点領域Ｒの変位のうちの血管径方向成分の変位から血管１２の径変位量が算出される。また、血管径補完情報処理部９４に対応するＳ７では、そのＳ６により算出された血管１２の径変位量とＳ３により算出された長軸血管径とに基づいて、血管径補完情報である短軸画像に由来する補完用血管径が算出される。そして、血管径測定部９６に対応するＳ８では、常には、Ｓ３により算出された長軸血管径の時系列的変化を血管径測定値として用いるが、図６に示すように体動による超音波センサの移動によって血管を見失って追従不可能となり図７に示すように長軸血管径の時系列的変化が得られない空白区間が形成される場合には、Ｓ７により生成された補完用血管径を上記空白区間の血管径として用いて血管径の時系列的変化が生成(測定)され、記録される。【００３０】
上述のように、本実施例の生体の血管径連続測定装置１０によれば、超音波反射信号に基づいて繰り返し生成された短軸画像間における複数の特徴点Ｐの周囲の特徴点領域Ｒの変位がパターンマッチングにより逐次求められ、それら複数の特徴点領域Ｒの変位から血管１２の径変化量が逐次算出され、その血管１２の径変化量に基づいて血管径の時系列的変化が測定されるので、体動などによる超音波センサのずれに影響され難い短軸画像から生体の血管の径の時系列的変化が得られる。【００３１】
また、本実施例の生体の血管径連続測定装置１０によれば、変位算出部９０により算出される、図９の実線の血管径方向成分の特徴点Ｐの変化量を微分することでスパイク状に得られた、拡張末期タイミング信号Ｍを心拍同期信号として生成する拡張末期タイミング発生部１０２が設けられているので、たとえば、血管径変位量算出部９２は、その変位算出部９０により算出された特徴点領域Ｒの変位のうち上記心拍同期信号に同期した血管径方向成分の変位から血管１２の径変化量を算出し、血管径測定部９６は、拡張末期タイミング信号Ｍに同期して１拍毎の径変化量を決定し、その１拍毎の径変化量から前記血管径を時系列的に決定する。このため、拡張末期タイミング信号としてＥＣＧ信号を検出する心電誘導信号検出装置を用いることが不要となる利点がある。【００３２】
また、本実施例の生体の血管径連続測定装置１０によれば、ハイブリッドプローブユニット２４（超音波センサ）は、血管１２に対して直交するように配置される第１短軸用超音波アレイ探触子Ａおよび第２短軸用超音波アレイ探触子Ｂと血管１２の直上にそれと平行な方向に配置される長軸用超音波アレイ探触子Ｃと、長軸用超音波アレイ探触子Ｃからの超音波反射信号に基づいて長軸での血管径を逐次算出する長軸血管径算出部８６を含む。短軸画像生成部８８は、その第１短軸用超音波アレイ探触子Ａまたは第２短軸用超音波アレイ探触子Ｂにより検出された超音波反射信号に基づいて前記血管１２の短軸画像を繰り返し生成するものである。血管径測定部９６は、長軸血管径算出部８６により算出された血管径の時系列的変化を血管径測定値として用いるが、その血管径の時系列的変化が得られない空白区間がある場合には、血管径変化量算出部９２により短軸画像中の特徴点領域Ｒの変位のうちの径方向成分の変位から算出された血管１２の径変化量を用いてその空白区間内の血管径の時系列的変化を生成する。このため、体動などによる長軸用超音波アレイ探触子Ｃのずれにより血管１２を見失った場合でも、時系列的に連続した血管１２の径変化量を測定できる。【００３３】
また、本実施例の生体の血管径連続測定装置１０によれば、長軸用超音波アレイ探触子Ｃからの超音波反射信号に基づいてその超音波反射信号の検波波形のピーク値を結ぶ包絡線信号であるＡモード波形を生成するＡモード信号処理部８４を含み、長軸血管径算出部８６は、そのＡモード信号処理部８４から生成されたＡモード波形が示すピーク波形間の時間差から長軸での血管径を逐次算出するものである。このため、長軸画像から血管径を測定する場合に比較して最小計測単位が１画素とならないため、長軸での血管径の測定精度が高められる。【００３４】
また、本実施例の生体の血管径連続測定装置１０によれば、長軸血管径算出部８６は、Ａモード信号処理部８４から生成されたＡモード波形が示すピーク波形間の時間差のうち、拡張末期タイミング信号発生部１０２からの拡張末期タイミング信号Ｍ（心拍同期振動）に同期した時間差から長軸での血管径を逐次算出するものである。このため、長軸での血管径の算出のために、拡張末期タイミング信号としてＥＣＧ信号を検出する心電誘導信号検出装置を用いることが不要となる利点がある。【００３５】
また、本実施例の生体の血管径連続測定装置１０によれば、血管径変位量算出部９２は、オプティカルフロー解析を用いて、短軸画像生成部８８により繰り返し生成された短軸画像間における複数の特徴点領域Ｒの変位を、その特徴点領域のパターンと一致する領域を検出することによりそれぞれ逐次算出するものであるので、高い精度でフロー推定および速度ベクトルが得られ、得られる特徴点Ｐ間の変位の精度が高められる。
【実施例２】
【００３６】
次に、本発明の他の実施例について説明する。なお、以下の説明において、前述の実施例と重複する部分については、同一の符号を付してその説明を省略する。【００３７】
図１２は、本発明の他の実施例の電子制御装置３８’の制御機能の要部を説明する機能ブロック線図であり、図１３はその電子制御装置３８’の制御作動の要部を説明するフローチャートである。【００３８】
本実施例の電子制御装置３８’は、前述の電子制御装置３８と比較して、Ａモード信号処理部８４、Ａモード波形ピーク位置による長軸血管径算出部８６、および血管径補完情報処理部９４が除去されるとともに、血管径初期値算出部８４’が新たに加えられている点で相違し、他は同様に構成されている。以下では、相違点を中心に説明する。【００３９】
図１２において、血管径変位量算出部９２は、血管径の相対的変位を逐次算出するが、血管径の絶対値を算出することができない。このため、血管径初期値算出部８４’は、前述のＡモード波形ピーク位置による長軸血管径算出部８６のようにＡモード波形ピーク位置間の時間差に基づいて、あるいはＡモード信号から生成される長軸画像上に表れる血管壁の間隔に基づいて、たとえば圧迫帯による圧迫前の血管１２の血管径の絶対値を算出する。血管径変化測定部９６は、その血管径初期値算出部８４’により求められた血管径(絶対値)と、血管径変位量算出部９２により求められた血管径の相対的変位量とに基づいて、血管径の経時的変化を測定し、記憶するとともに、図５に示すような血管径連続測定値を画像表示装置４０に表示させる。【００４０】
図１３において、ステップＳ１１では、前述のＳ１と同様に超音波反射信号が読み込まれる。ついで、血管径初期値算出部８４’に対応するＳ１２では、前述のＡモード波形ピーク位置による長軸血管径算出部８６のようにＡモード波形ピーク位置間の時間差に基づいて、あるいはＡモード信号から生成される長軸画像上に表れる血管壁の間隔に基づいて、たとえば圧迫帯による圧迫前の血管１２の血管径の初期値(絶対値)が算出される。次に、短軸画像生成部８８に対応するＳ１３では、前述のＳ４と同様に短軸画像が繰り返し生成される。続く、変位算出部９０に対応するＳ１４では、前述のＳ５と同様に、特徴点Ｐの変位が算出され、血管径変位量算出部９２に対応するＳ１５では、前述のＳ６と同様に、血管径の変位量が算出される。ついで、血管径測定部９６に対応するＳ１６では、Ｓ１２で求められた血管１２の血管径のの初期値(絶対値)と、Ｓ１５で求められた血管径の変位量とに基づいて血管径の連続測定値が求められる。そして、表示制御部９８に対応するＳ１７では、前述のＳ９と同様に、血管径の連続測定値が画像表示装置４０に表示される。【００４１】
本実施例によれば、体動によるずれの影響の殆どない短軸画像から、パターンマッチングを用いたオグテイカルフロー解析により短軸画像内の特徴点領域Ｒ間の変位が求められ、その特徴点領域Ｒ間の変位から血管径方向成分の血管系変位量が求められ、その血管径変位量と血管径初期値(絶対値)とから、血管１２の血管径連続測定値が求められるので、前述の実施例と同様の効果が得られる。【００４２】
以上、本発明の一実施例を図面を参照して詳細に説明したが、本発明はこの実施例に限定されるものではなく、別の態様でも実施され得る。【００４３】
たとえば、前述の実施例では、繰り返し生成される短軸画像間の特徴点領域Ｒの変位検出をオプティカルフロー解析を行なうことで実現していたが、必ずしもオプティカルフロー解析を用いる必要はなく、短軸画像間の特徴点領域Ｒの変位をパターンマッチングで求めればよい。【００４４】
また、前述の実施例において、繰り返し生成される短軸画像間の特徴点領域Ｒの変位をオプティカルフロー解析を用いて求められていたが、拘束直線の交点群のクラスターリングや最小自乗法を用いてフロー推定を行なって短軸画像間の特徴点領域Ｒの変位が求められてもよい。【００４５】
また、前述の実施例では、圧迫帯により圧迫前から圧迫解放後の所定時間経過後までの超音波反射信号が読み込まれた後に演算処理が行なわれて血管１２の連続測定値が算出されるというパッチ処理であったが、超音波反射信号が入力される毎に演算処理が実行されるリアルタイム処理であってもよい。【００４６】
また、前述の実施例では、血管径が測定される血管１２は皮膚直下に位置していたが、ハイブリッドプローブユニット２４を体腔内へ挿入することにより、体腔内の臓器内に位置するものであってもよい。【００４７】
なお、上述したのはあくまでも一実施形態であり、その他一々例示はしないが、本発明は、その主旨を逸脱しない範囲で当業者の知識に基づいて種々変更、改良を加えた態様で実施することができる。
【符号の説明】
【００４８】
１０：血管径連続測定装置
１１：上肢(生体)
１２：血管
２４：ハイブリッドプローブユニット（超音波センサ）
８４：Ａモード信号処理部
８６：長軸血管径算出部
８８：短軸画像生成部
９０：変位算出部
９２：血管径変位量算出部
９６：血管径測定部
１０２：拡張末期タイミング発生部
Ａ：第１短軸用超音波アレイ探触子
Ｂ：第２短軸用超音波アレイ探触子
Ｃ：長軸用超音波アレイ探触子
【図１】
0005998197.tif 000002
【図２】
0005998197.tif 000003
【図３】
0005998197.tif 000004
【図４】
0005998197.tif 000005
【図５】
0005998197.tif 000006
【図６】
0005998197.tif 000007
【図７】
0005998197.tif 000008
【図８】
0005998197.tif 000009
【図９】
0005998197.tif 000010
【図１０】
0005998197.tif 000011
【図１１】
0005998197.tif 000012
【図１２】
0005998197.tif 000013
【図１３】
0005998197.tif 000014