(57)【要約】【課題】超音波測定の対象となる血管の位置を検出する新たな技術の提案。【解決手段】超音波測定装置は、超音波測定によって取得したＢモード画像の中から複数の特徴点を検出する。これを測定フレーム毎に繰り返し実行する。そして、時間的に隣接する（連続する）測定フレーム間での特徴点それぞれの変位ベクトルを算出する。次いで、特徴点毎に変位ベクトルに基づいて変位方向を算出する。この変位方向の交点を、血管３の断面中心位置４として検出する。【選択図】図２
2015159933.tif 000002
【特許請求の範囲】
【請求項１】
生体の表面に超音波を送信し、受信信号に基づいて、生体の断面を超音波で測定した測定データを取得する測定データ取得部と、
前記測定データに基づいて、生体内組織の変位を検出する変位検出部と、
前記検出された変位に基づいて、前記生体内の血管位置を検出する血管位置検出部と、
を備える超音波測定装置。
【請求項２】
前記変位検出部は、複数の測定データを時系列に比較することで前記変位を検出する、
請求項１に記載の超音波測定装置。
【請求項３】
前記変位検出部は、前記生体内組織の変位方向を検出し、
前記血管位置検出部は、前記変位方向に基づいて、前記血管位置を検出する、
請求項１又は２に記載の超音波測定装置。
【請求項４】
前記測定データに基づいて、生体内組織に係る複数の特徴点を検出する特徴点検出部、
を更に備え、
前記変位検出部は、前記複数の特徴点それぞれの変位方向を検出し、
前記血管位置検出部は、前記複数の特徴点それぞれの変位方向の交点を求めて前記血管位置を検出する、
請求項１〜３の何れかに記載の超音波測定装置。
【請求項５】
前記特徴点検出部は、血管或いは血管周辺の生体内組織に係る特徴点を検出する、
請求項４に記載の超音波測定装置。
【請求項６】
前記変位検出部は、血管の収縮及び弛緩の少なくとも一方に係る前記変位を検出する、
請求項１〜５の何れかに記載の超音波測定装置。
【請求項７】
生体の表面に超音波を送信し、受信信号に基づいて、生体の断面を超音波で測定した測定データを取得することと、
前記測定データに基づいて、生体内組織の変位を検出することと、
前記検出された変位に基づいて、前記生体内の血管位置を検出することと、
を含む超音波測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波を用いて血管位置を検出する超音波測定装置等に関する。
【背景技術】
【０００２】
超音波を用いて生体情報を測定する一例として、血管疾患の判断を含む血管機能の評価が行われている。例えば、動脈硬化の指標となる頸動脈のＩＭＴ（Intima Media Thickness：内膜中膜複合体厚）を測定することもその１つである。このような測定では、生体内組織の血管の位置や形状が測定される。
【０００３】
特許文献１には、頸動脈の断面画像であるＢモード画像に基づいて、頸動脈の位置や形状を推定してモデル化する技術が開示されている。係る技術では、心拍動による動脈の動きに着目し、フレーム毎に、モデルの評価関数を生成し、最適化を行って次フレームの頸動脈の位置や形状を推定し、モデル化することを繰り返している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００９−６６２６８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述の特許文献１に開示されている技術では、フレーム毎に、評価関数の生成、最適化、そしてモデル化を繰り返し行っているため、測定に係る演算処理が煩雑となり、演算量が増加するという問題があった。
【０００６】
本発明は、上記事情に鑑みてなされたものであり、その目的とするところは、超音波測定の対象となる血管の位置を検出する新たな技術を提案することである。
【課題を解決するための手段】
【０００７】
上記課題を解決するための第１の発明は、生体の表面に超音波を送信し、受信信号に基づいて、生体の断面を超音波で測定した測定データを取得する測定データ取得部と、前記測定データに基づいて、生体内組織の変位を検出する変位検出部と、前記検出された変位に基づいて、前記生体内の血管位置を検出する血管位置検出部と、を備える超音波測定装置である。
【０００８】
また、他の発明として、生体の表面に超音波を送信し、受信信号に基づいて、生体の断面を超音波で測定した測定データを取得することと、前記測定データに基づいて、生体内組織の変位を検出することと、前記検出された変位に基づいて、前記生体内の血管位置を検出することと、を含む超音波測定方法を構成しても良い。
【０００９】
この第１の発明等によれば、生体の断面を超音波で測定した測定データをもとに検出された生体内組織の変位に基づいて、生体内の血管位置が検出される。血管は、拍動によって周期的に収縮及び弛緩を繰り返しており、これに同期して血管の周囲の生体内組織も変位している。この生体内組織の変位を検出することで、その変位の原因となっている血管位置を検出するという新たな技術を実現できる。
【００１０】
また、第２の発明として、第１の発明の超音波測定装置であって、前記変位検出部が、複数の測定データを時系列に比較することで前記変位を検出する、超音波測定装置を構成しても良い。
【００１１】
この第２の発明によれば、複数の測定データを時系列に比較することで、生体内組織の変位を検出することができる。
【００１２】
また、第３の発明として、第１又は第２の発明の超音波測定装置であって、前記変位検出部が、前記生体内組織の変位方向を検出し、前記血管位置検出部が、前記変位方向に基づいて、前記血管位置を検出する、超音波測定装置を構成しても良い。
【００１３】
この第３の発明によれば、生体内組織の変位方向が検出され、変位方向に基づいて血管位置が検出される。血管の短軸方向断面において、血管は等方的に収縮及び弛緩を繰り返すため、血管やその周囲の生体内組織の変位方向は、血管の短軸断面の略中心位置からの放射方向となる。このため、検出した生体内組織の変位方向から、血管位置を検出することができる。
【００１４】
また、第４の発明として、第１〜第３の何れかの発明の超音波測定装置であって、前記測定データに基づいて、生体内組織に係る複数の特徴点を検出する特徴点検出部、を更に備え、前記変位検出部が、前記複数の特徴点それぞれの変位方向を検出し、前記血管位置検出部が、前記複数の特徴点それぞれの変位方向の交点を求めて前記血管位置を検出する、超音波測定装置を構成しても良い。
【００１５】
この第４の発明によれば、生体内組織に係る複数の特徴点それぞれの変位方向の交点を求めて、血管位置が検出される。血管は短軸方向断面において等方的に収縮及び弛緩を繰り返すため、複数の特徴点それぞれの変位方向の交点を、血管の断軸断面の略中心位置として検出することができる。
【００１６】
また、第５の発明として、第４の発明の超音波測定装置であって、前記特徴点検出部が、血管或いは血管周辺の生体内組織に係る特徴点を検出する、超音波測定装置を構成しても良い。
【００１７】
この第５の発明によれば、特徴点として、血管やその周辺の生体内組織に係る特徴点が検出される。これにより、血管の収縮及び弛緩に伴う特徴点の変位を検出することができ、血管位置の検出が容易となる。
【００１８】
また、第６の発明として、第１〜第５の何れかの発明の超音波測定装置であって、前記変位検出部が、血管の収縮及び弛緩の少なくとも一方に係る前記変位を検出する、超音波測定装置を構成しても良い。
【００１９】
この第６の発明によれば、生体内組織に係る変位として、血管の収縮及び弛緩の少なくとも一方に係る変位が検出される。血管は収縮及び弛緩を繰り返しているため、検出される生体内組織に係る変位方向は、略直線状となる。このため、血管の収縮及び弛緩の少なくとも一方に係る変位を検出することで、生体内組織に係る変位方向を検出することが可能となる。
【図面の簡単な説明】
【００２０】
【図１】超音波測定装置の全体構成図。
【図２】血管の収縮及び弛緩を説明する図。
【図３】変位ベクトルの算出を説明する図。
【図４】変位方向の定め方の説明図。
【図５】Ｂモード画像中の特徴点の検出を説明する図。
【図６】超音波測定装置の機能構成図。
【図７】超音波測定処理のフローチャート。
【発明を実施するための形態】
【００２１】
［全体構成］
図１は、本実施形態における超音波測定装置１０の構成例を示す図である。超音波測定装置１０は、超音波を用いて被検者の生体情報を測定する装置である。本実施形態では、生体情報として、頸動脈のＩＭＴ（Intima Media Thickness：血管の内膜中膜複合体厚）といった血管機能情報を測定する。勿論、ＩＭＴ以外にも、血管径や、血管径から血圧を測定するといった別の血管機能情報を測定することとしても良い。
【００２２】
超音波測定装置１０は、タッチパネル１２と、キーボード１４と、超音波プローブ１６と、本体装置２０とを備える。本体装置２０には、制御基板２２が搭載されており、タッチパネル１２やキーボード１４、超音波プローブ１６等の各部と信号送受信可能に接続されている。
【００２３】
制御基板２２には、ＣＰＵ（Central Processing Unit）２４や、ＡＳＩＣ（Application Specific Integrated Circuit）、各種集積回路の他、ＩＣメモリーやハードディスクによる記憶媒体２６と、外部装置とのデータ通信を実現する通信ＩＣ２８とが搭載されている。本体装置２０は、ＣＰＵ２４等が記憶媒体２６に記憶されている制御プログラムを実行することにより、超音波測定をはじめとする本実施形態に係る各種機能を実現する。
【００２４】
具体的には、本体装置２０は、超音波プローブ１６から、被検者２の生体内組織に向けて超音波ビームを送信・照射し、反射波を受信する。そして、反射波の受信信号を増幅・信号処理することにより、被検者の生体内構造に係る測定データを生成することができる。測定データには、いわゆるＡモード、Ｂモード、Ｍモード、カラードップラーの各モードの画像がこれに含まれる。超音波を用いた測定は、所定周期で繰り返し実行される。測定単位を「フレーム」と呼ぶ。
【００２５】
超音波プローブ１６は、複数の超音波振動子を配列して構成されている。本実施形態では、配列数を一列とするが、複数列として面状の配列構成としてもよい。そして、各超音波振動子からの超音波が被検者の頸動脈の短軸方向に横断する相対姿勢で超音波プローブ１６を被検者２の頸部に固定し、生体情報の測定が行われる。
【００２６】
［原理］
血管機能情報の測定にあたり、先ずは血管位置の検出が行われる。図２は、血管３の短軸断面の模式図である。血管３は、拍動にあわせて等方的に収縮及び拡張を繰り返す。その収縮及び拡張方向は、血管３の短軸断面の中心位置４からの放射方向となる。このことを利用し、本実施形態では、測定データの一つであるＢモード画像を用いて、血管３の周辺の生体内組織の変位を検出し、その変位方向から血管３の断面中心位置４を検出する。図２中の破線丸印は、図３に示した特徴点の一部であり、破線丸印の中心位置の変位方向を実線両矢印で示している。
【００２７】
また、本実施形態における座標系として、皮膚面（体表面）からの深さ方向をＺ軸、血管３の長手方向（血流方向）をＹ軸、Ｚ軸及びＹ軸に直交する血管３の短軸方向をＸ軸とする直交３次元座標を用いることとする。
【００２８】
より具体的に説明する。図３に示すように、まず、連続する２つの測定フレームで測定されたＢモード画像Ｐ１，Ｐ２それぞれにおいて画像中の特徴点（図３中の破線丸印の中心位置）を抽出し、２つのＢモード画像Ｐ１，Ｐ２間で特徴点の対応付け推定（追跡とも呼ばれる）を行う。対応付けを行うことができなかった特徴点は、特徴量が小さい特徴点であるとして削除する。これにより、時間経過に対する位置の特定、すなわち追跡できる特徴点のみを検出する。特徴点の抽出方法及び対応付けの推定方法は、適宜公知技術を用いることで実現できるので、ここでの説明は省略する。
【００２９】
次いで、追跡できる特徴点として検出された特徴点毎に、連続するフレーム間での位置変化を表す変位ベクトルＶを算出する。図４は、検出されたある１つの特徴点に着目して、第ｔフレームの位置と、その次の第ｔ＋１フレームの位置とに基づいて、変位ベクトルＶｔ＋１を求めた例を示す図である。検出された特徴点毎に、変位ベクトルＶを算出する。
【００３０】
そして、変位ベクトルＶに沿った方向を、当該特徴点の変位方向とする。このとき、血管位置の検出精度を上げるため、特徴点１つ１つについて次のようにして変位方向を決定する。すなわち、１つの特徴点について位置を追跡し、連続する複数フレーム間それぞれにおける複数の変位ベクトルＶを求める。求めた複数の変位ベクトルＶそれぞれの方向のうち、出現頻度の最多の方向を当該特徴点の変位方向として決定する。多数決方式で方向を決定するのである。これを各特徴点について行い、各特徴点の変位方向を決定する。
【００３１】
変位方向は、例えば、図５に示すように、Ｘ軸に対して変位ベクトルＶがなす角度θとして表される。図２に示したように、特徴点の位置は血管３の中心から放射方向に往復するように変化するため、角度θは「０°〜１８０°」の範囲とする。つまり、図５における変位ベクトルのように、向きが互いに逆となる変位ベクトルＶは、その変位方向が同じとみなす。具体的には、変位ベクトルＶは、位置座標と方向（方位ともいえる）と大きさとで定義される。そのため、変位ベクトルＶの方向とＸ軸正方向との反時計回りの成す角度（時計回りでの成す角度でも勿論よい）が「０°〜１８０°」であれば当該角度を変位方向の角度θとし、「１８０°〜３６０°」であれば当該角度から１８０°を減算した角度を変位方向の角度θとする。
【００３２】
複数の特徴点それぞれの変位方向を決定すると、図２に示したように、これらの複数の特徴点それぞれの位置を通り、当該特徴点の変位方向に沿った方向の直線（図２中の点線）を求める（以下この直線を「変位方向線」という）。そして、変位方向線の交点を算出して、この変位方向線の交点を血管３の断面中心位置４とする。なお、このときに求められる交点は、理論上は１点となるが、複数の変位方向線の中から２つの変位方向線の組み合わせを選択して交点を求めることを各組み合わせについて行うと、１点とはならない場合がある。そこで、変位方向線の組み合わせ別に求めた交点のうち、位置が重複する数が最多の交点を、血管３の断面中心位置４とすることができる。また、別の例として、求めた交点の位置の集合が、およそ一定の近傍位置に収まることを利用する方法も考えられる。例えば、この近傍位置の中心を断面中心位置４とすることができる。近傍位置を求める際に、求めた交点の位置の集合の８割程度が含まれる位置範囲を近傍位置とするなど、測定上の誤差を取り除く処理を導入してもよいことは勿論である。
【００３３】
［機能構成］
図６は、超音波測定装置１０の機能構成図である。図６に示すように、超音波測定装置１０は、操作入力部１１０と、表示部１２０と、音出力部１３０と、通信部１４０と、処理部２００と、記憶部３００とを備えて構成される。
【００３４】
操作入力部１１０は、ボタンスイッチやタッチパネル、各種センサー等の入力装置によって実現され、なされた操作に応じた操作信号を処理部２００に出力する。図１では、タッチパネル１２やキーボード１４がこれに該当する。
【００３５】
表示部１２０は、ＬＣＤ（Liquid Crystal Display）等の表示装置によって実現され、処理部２００からの表示信号に基づく各種表示を行う。図１では、タッチパネル１２がこれに該当する。音出力部１３０は、スピーカー等の音出力装置によって実現され、処理部２００からの音信号に基づく各種音出力を行う。
【００３６】
通信部１４０は、無線ＬＡＮ（Local Area Network）やＢｌｕｅｔｏｏｔｈ（登録商標）等の無線通信装置、モデム、有線用の通信ケーブルのジャックや制御回路等の通信装置によって実現され、所与の通信回線と接続して外部装置との通信を行う。図１では、制御基板２２に搭載された通信ＩＣ２８がこれに該当する。
【００３７】
処理部２００は、ＣＰＵ（Central Processing Unit）やＧＰＵ（Graphics Processing Unit）等のマイクロプロセッサー、ＡＳＩＣ（特定用途向け集積回路：Application Specific Integrated Circuit）、ＩＣ（Integrated Circuit）メモリー等の電子部品によって実現され、記憶部３００に記憶されたプログラムやデータ、操作入力部１１０からの操作信号等に基づいて各種演算処理を実行して、超音波測定装置１０の動作を制御する。また、処理部２００は、超音波測定制御部２１０と、血管位置検出部２３０と、血管機能測定部２５０とを有する。
【００３８】
超音波測定制御部２１０は、超音波プローブ１６における超音波の送受信を制御する。具体的には、所定周期の送信タイミングで、超音波プローブ１６から超音波を送信させる。また、超音波プローブ１６において受信された超音波の反射波の信号の増幅等を行う。
【００３９】
測定データ生成部２２０は、超音波プローブ１６による反射波の受信信号をもとに、ＡモードやＢモード、Ｍモードといった各モードの画像データを含む測定データを生成する。この測定データ生成部２２０が、測定データ取得部に相当する。
【００４０】
血管位置検出部２３０は、特徴点検出部２３２と、特徴点対応付部２３４と、変位ベクトル算出部２３６と、変位方向算出部２３８と、変位方向交点算出部２４０とを有し、測定データ生成部２２０によって生成された測定データに基づく血管位置の検出を行う。
【００４１】
特徴点検出部２３２は、測定フレーム毎のＢモード画像中の特徴点を検出する。このとき、所定時間（例えば１〜２秒間）の間のＢモード画像それぞれに現れており、特徴量が所定以上の特徴点、すなわち位置の追跡が可能な特徴点を検出する。
【００４２】
特徴点対応付部２３４は、特徴点検出部２３２によって検出された特徴点それぞれのフレーム間での対応付けを行う。各特徴点の画像上の特徴（当該特徴点の輝度値や周辺輝度値）等から対応付けを行うことができる。
【００４３】
変位ベクトル算出部２３６は、特徴点対応付部２３４によって対応付けされた各特徴点のフレーム間毎の変位ベクトルＶを算出する。
【００４４】
変位方向算出部２３８は、変位ベクトル算出部２３６によって算出された変位ベクトルＶから、特徴点検出部２３２によって検出された特徴点それぞれの変位方向を算出する。具体的には、特徴点それぞれについて、変位ベクトルＶの方向の出現頻度を求め、最多の方向を当該特徴点の変位方向として決定する。この変位ベクトル算出部２３６、及び、変位方向算出部２３８が、変位検出部に相当する。
【００４５】
変位方向交点算出部２４０は、各特徴点について、当該特徴点の位置座標を通り、変位方向算出部２３８によって算出された変位方向に沿った方向の変位方向線を求める。そして、各特徴点の変位方向線の交点を算出する。交点が１つに定まらない場合の決定方法は上述した通りである。この交点が血管の断面中心位置として検出され、Ｂモード画像中の血管中心位置座標３４０として記憶される。
【００４６】
血管機能測定部２５０は、所与の血管機能情報の測定を行う。具体的には、検出された血管位置によって特定された血管の血管径やＩＭＴ等を測定したり、血管前壁及び血管後壁をトラッキングすることで血管径変動から血圧を推定演算したり、脈拍数を算出するといった血管機能情報の測定を行う。
【００４７】
記憶部３００は、ＲＯＭやＲＡＭ、ハードディスク等の記憶装置によって実現され、処理部２００が超音波測定装置１０を統合的に制御するためのプログラムやデータ等を記憶しているとともに、処理部２００の作業領域として用いられ、処理部２００が実行した演算結果や、操作入力部１１０からの操作データ等が一時的に格納される。図１では、制御基板２２に搭載される記憶媒体２６がこれに該当する。本実施形態では、記憶部３００には、超音波測定プログラム３１２と、Ｂモードデータ３２０と、特徴点データ３３０と、血管中心位置座標３４０とが記憶される。
【００４８】
Ｂモードデータ３２０は、測定フレーム毎に生成され、フレームＩＤと、Ｂモード画像データと、Ｂモード画像中の各特徴点の位置座標のリストである特徴点座標リストとが格納される。
【００４９】
特徴点データ３３０は、特徴点検出部２３２によって検出された特徴点毎に生成され、特徴点ＩＤと、測定フレーム毎のＢモード画像中の位置座標である位置座標データと、方向別頻度データと、変位方向とを格納している。方向別頻度データは、変位方向算出部２３８によって使用されるデータであり、変位ベクトル算出部２３６によって各測定フレームにおいて算出された変位ベクトルＶの方向を、方向別に集計した頻度（度数）のデータである。各測定フレームの変位ベクトルＶの数を、方向毎に累積したデータともいえる。なお、特徴点データ３３０に、各測定フレームにおいて算出された変位ベクトルＶのデータを含めることとしてもよい。
【００５０】
［処理の流れ］
図７は、超音波測定処理の流れを説明するフローチャートである。この処理は、処理部２００が超音波測定プログラム３１２を実行することで実現される。
【００５１】
処理部２００は、先ず、超音波プローブ１６による超音波の送受信を開始させる（ステップＳ１）。そして、測定フレーム毎にループＡの処理を繰り返す。ループＡでは、測定データ生成部２２０が、超音波プローブ１６による超音波反射波の受信信号をもとに、Ｂモード画像を生成する（ステップＳ３）。次いで、特徴点検出部２３２が、生成されたＢモード画像から特徴点を検出する（ステップＳ５）。
【００５２】
続いて、検出された特徴点それぞれを対象とするループＢの処理を行う。ループＢでは、処理対象の特徴点を、直前フレームのＢモード画像中の特徴点と対応付けし（ステップＳ７）、直前フレーム間との変位ベクトルＶを算出する（ステップＳ９）。次いで、処理対象の特徴点の特徴点データ３３０に含まれる方向別頻度データのうち、ステップＳ９で算出した変位ベクトルＶの方向の頻度を更新（例えば、「１」加算）する（ステップＳ１１）。そして、方向別頻度データに集計された最多の方向を、処理対象の特徴点の変位方向として決定する（ステップＳ１３）。
【００５３】
全ての特徴点を対象としたループＢの処理を行うと、続いて、特徴点それぞれの変位方向線を求め、その変位方向線の交点を算出することで、血管の断面中心位置を算出する（ステップＳ１５）。このとき、変位方向線２つずつの組み合わせの交点を求めるが、所定数の交点の座標が一致する（すなわち所定数の変位方向線が一点で交差する（詳細には、所定の狭範囲を通過する））場合に、当該交点を、血管の断面中心位置とする。血管の断面中心位置が算出されないならば（ステップＳ１７：ＮＯ）、次の測定フレームを対象としたループＡの処理を行う。測定開始から時間があまり経過していない場合には、方向別頻度データに蓄積された値が少ないため、ステップＳ１７ではＮＯと判定される可能性が高い。そして、血管の断面中心位置が算出されたならば（ステップＳ１７：ＹＥＳ）、ループＡの処理を終了する。
【００５４】
ループＡの処理を終了すると、血管位置の検出がなされたことを報知する（ステップＳ１９）。続いて、血管機能測定部２５０が、超音波プローブ１６による超音波の送受信結果を用いた所与の血管機能情報の測定を行い、測定結果の記憶・表示を行う（ステップＳ２１）。以上の処理を行うと、超音波測定処理は終了となる。
【００５５】
［作用効果］
本実施形態によれば、超音波測定装置１０は、超音波測定によって取得した時系列の複数のＢモード画像から、複数の特徴点それぞれの変位を算出し、これらの変位方向線の交点を、血管の断面中心位置として検出することができる。よって、超音波測定の対象とする血管位置を検出する新たな手法を実現することができる。また、その手法も超音波反射波の受信信号を用いる、比較的簡単な方法であるため、演算処理量が少なく、高速に演算することが可能である。
【００５６】
［変形例］
なお、本発明の適用可能な実施形態は上述の実施形態に限定されることなく、本発明の趣旨を逸脱しない範囲で適宜変更可能なのは勿論である。
【００５７】
例えば、Ｂモード画像に基づく特徴点それぞれの変位方向の算出は、血管の収縮及び弛緩の少なくとも一方の期間において行うこととしてもよい。血管は心臓の拍動に応じて収縮及び弛緩を繰り返し、この収縮及び弛緩の一周期の間に特徴点が一往復するように変位する。このため、血管の収縮及び弛緩の少なくとも一方の期間における特徴点の変位を検出することで、速やかに血管位置（断面中心位置）を検出することができる。
【符号の説明】
【００５８】
１０  超音波測定装置、１６  超音波プローブ、３０  本体装置、１１０  入力部、１２０  表示部、１３０  音出力部、１４０  通信部、２００  処理部、２１０  超音波測定制御部、２２０  測定データ生成部、２３０  血管位置検出部、２３２  特徴点検出部、２３４  特徴点対応付部、２３６  変位ベクトル算出部、２３８  変位方向算出部、２４０  変位方向交点算出部、２５０  血管機能測定部、３００  記憶部、３１０  システムプログラム、３１２  超音波測定プログラム、３２０  Ｂモードデータ、３３０  特徴点データ、３４０  血管中心位置座標、２  被検者、３  血管（頸動脈）
【図１】
2015159933.tif 000003
【図２】
2015159933.tif 000004
【図４】
2015159933.tif 000005
【図５】
2015159933.tif 000006
【図６】
2015159933.tif 000007
【図７】
2015159933.tif 000008
【図３】
2015159933.tif 000009