(57)【特許請求の範囲】
【請求項１】
超音波プローブにより被検体から取得された反射超音波に基づいて前記被検体の検査対象血管を特定する超音波診断装置であって、
前記反射超音波に基づいて、前記被検体の断層画像を生成する断層画像生成部と、
前記反射超音波に基づいて、前記断層画像上での前記被検体の血流領域を示す血流情報を生成する血流情報生成部と、
前記血流情報生成部が生成した前記血流情報を解析することで、前記血流領域の代表点を示す血流点を抽出し、同じ特徴を持った血流点どうしを１つの血流領域群として抽出し、前記血流領域群を単位として前記血流領域が前記検査対象血管に対応する血流領域であるか否かの判定を行う血流領域判定部とを備える
超音波診断装置。
【請求項２】
前記血流領域判定部は、
前記血流情報として、前記断層画像における前記血流領域の位置を解析することで前記判定を行う
請求項１に記載の超音波診断装置。
【請求項３】
前記血流領域判定部は、
前記血流領域を含む前記断層画像が取得された前記被検体上の位置に対応させて前記血流領域を配置した場合に、配置した前記血流領域に対応する血流点の間の距離が閾値以下である複数の血流領域をまとめて１つの前記血流領域群として抽出し、抽出した前記血流領域群の属性に基づいて、前記血流領域群に含まれる血流領域に対して前記判定を行う
請求項１又は２に記載の超音波診断装置。
【請求項４】
前記血流領域判定部は、
抽出した前記血流領域群の属性があらかじめ定められた前記検査対象血管の属性に一致する場合に、抽出した前記血流領域群に含まれる血流領域が前記検査対象血管に対応する血流領域であるとの前記判定を行う
請求項３に記載の超音波診断装置。
【請求項５】
前記血流領域判定部は、
抽出した前記血流領域群に含まれる血流領域の数と、前記血流領域群に含まれる血流領域のうち面積が最大である血流領域の面積と、前記血流領域群を含む複数の断層画像が取得された被検体上の範囲のうちの始端及び終端から取得された断層画像のそれぞれにおける血流領域の位置とのうちの少なくとも１つを、前記血流領域群の属性として、前記判定を行う
請求項４に記載の超音波診断装置。
【請求項６】
前記血流領域判定部は、
１つの前記血流領域群が含む前記血流領域の数が所定数以下の場合には、当該血流領域群を除く血流領域群に含まれる血流領域に対して、前記判定を行う
請求項５に記載の超音波診断装置。
【請求項７】
前記血流領域判定部は、
前記血流領域群のうち、面積が最大である血流領域を含む血流領域群に含まれる血流領域が、前記検査対象血管に対応する血流領域であるとの判定を行う
請求項５に記載の超音波診断装置。
【請求項８】
前記血流領域判定部は、さらに、
抽出した前記血流領域群のそれぞれを含む断層画像が取得された前記被検体の位置に対応させて、前記血流領域群のそれぞれを配置した場合に、第一血流領域群を補間することで形成される図形の端部と、第二血流領域群を補間することで形成される図形の一部との距離が所定値以下であるとき、前記第一血流領域群と前記第二血流領域群とをまとめて１つの血流領域群として新たに抽出し、新たに抽出した血流領域群に対して前記判定を行う
請求項５に記載の超音波診断装置。
【請求項９】
前記血流領域判定部は、さらに、
抽出した前記血流領域群のそれぞれを含む前記断層画像が取得された前記被検体上の位置に対応させて前記血流領域群のそれぞれを配置した場合に、前記血流領域群に対応する前記被検体上の範囲のうちの始端又は終端から取得された断層画像に含まれる第一血流領域の前記断層画像上における第一位置と、前記被検体上の範囲のうちの始端又は終端から取得された断層画像に含まれる血流領域であって前記第一血流領域と異なる血流領域である第二血流領域の前記断層画像上における第二位置との差が所定値以下であり、かつ、
前記第一血流領域が取得された前記被検体の位置と前記第二血流領域が取得された前記被検体の位置とが所定距離以下であるときに、
前記第一血流領域を含む前記血流領域群と、前記第二血流領域を含む前記血流領域群とをまとめて１つの血流領域群として新たに抽出し、新たに抽出した血流領域群に対して前記判定を行う
請求項５に記載の超音波診断装置。
【請求項１０】
前記血流領域判定部は、
前記血流情報として、前記断層画像における前記血流領域の面積を解析することで前記判定を行う
請求項１に記載の超音波診断装置。
【請求項１１】
前記血流領域判定部は、
取得された複数の断層画像のそれぞれにおける前記血流領域の面積の変動を解析し、拍動様の変動の有無を判定することで前記判定を行う
請求項１０に記載の超音波診断装置。
【請求項１２】
前記血流領域判定部は、
前記血流情報として、前記血流領域の血流の方向を解析することで前記判定を行う
請求項１に記載の超音波診断装置。
【請求項１３】
前記血流領域判定部は、
前記血流領域の前記血流の方向が、あらかじめ定められた前記検査対象血管の血流の方向と一致する場合に、前記血流領域が前記検査対象血管に対応する前記血流領域であるとの前記判定を行う
請求項１２に記載の超音波診断装置。
【請求項１４】
前記超音波診断装置は、さらに、
前記血流領域判定部における前記判定の結果を示す情報を表示する表示部と、
前記表示部が表示している前記判定の結果を示す情報に対するユーザによる修正要求に基づいて前記血流領域判定部における前記判定の結果を修正する修正部とを備える
請求項１〜１３のいずれか１項に記載の超音波診断装置。
【請求項１５】
前記超音波診断装置は、さらに、
前記血流領域判定部が判定した血流領域に基づいて、仮血管輪郭を設定する仮血管輪郭設定部と、
前記仮血管輪郭設定部が設定した前記仮血管輪郭を用いて、前記断層画像生成部が生成した前記断層画像内において、前記検査対象血管の輪郭を抽出する血管輪郭抽出部とを備える
請求項１〜１４のいずれか１項に記載の超音波診断装置。
【請求項１６】
前記超音波診断装置は、さらに、
前記超音波プローブの位置又は姿勢を示す位置姿勢情報を取得するプローブ位置姿勢取得部と、
前記プローブ位置姿勢取得部が取得した前記位置姿勢情報と、前記血流情報生成部が生成した前記血流情報とに基づいて、前記血流領域を３次元空間上で示す血流３次元情報を生成する血流３次元構成部とを備え、
前記血流領域判定部は、
前記血流３次元構成部が生成した前記血流３次元情報を解析することで、前記判定を行う
請求項１〜１５のいずれか１項に記載の超音波診断装置。
【請求項１７】
前記検査対象血管は、頸動脈である
請求項１〜１６のいずれか１項に記載の超音波診断装置。
【請求項１８】
超音波プローブにより被検体から取得された反射超音波に基づいて前記被検体の検査対象血管を特定する血管特定方法であって、
前記反射超音波に基づいて、前記被検体の断層画像を生成する断層画像生成ステップと、
前記反射超音波に基づいて、前記断層画像上での前記被検体の血流領域を示す血流情報を生成する血流情報生成ステップと、
前記血流情報生成ステップで生成した前記血流情報を解析することで、前記血流領域の代表点を示す血流点を抽出し、同じ特徴を持った血流点どうしを１つの血流領域群として抽出し、前記血流領域群を単位として前記血流領域が前記検査対象血管に対応する血流領域であるか否かの判定を行う血流領域判定ステップとを含む
血管特定方法。
【請求項１９】
請求項１８に記載の血管特定方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波診断装置および血管特定方法に関する。特に、超音波プローブにより被検体から取得された反射超音波に基づいて検査対象血管を特定する超音波診断装置および血管特定方法に関する。
【背景技術】
【０００２】
近年、動脈硬化及び血管疾患等を早期に発見するために、超音波診断装置が用いられている。具体的には、超音波診断装置を用いて、血管壁の内膜と中膜とをあわせた膜の厚さである内膜中膜複合体厚（Ｉｎｔｉｍａ  Ｍｅｄｉａ  Ｔｈｉｃｋｎｅｓｓ、以下ＩＭＴ）が計測されている。ＩＭＴの厚みやプラークの有無を判断する際には、検査者は、超音波画像内の血管外膜及び血管内膜の形状を手動で指定（スケッチ）し、このスケッチを元に診断を行う。具体的には、検査者は、Ｂモード画像に対し、外膜の辺縁形状を外膜輪郭線としてスケッチし、さらに内腔の辺縁形状を内腔輪郭線としてスケッチする。そして最終的にこのスケッチの形状をもとに、検査者は、プラークの有無などの診断を行う（非特許文献１参照）。
【先行技術文献】
【非特許文献】
【０００３】
【非特許文献１】Ａｉｎｓｗｏｒｔｈ  ＣＤ、  Ｂｌａｋｅ  ＣＣ、  Ｔａｍａｙｏ  Ａ、  Ｂｅｌｅｔｓｋｙ  Ｖ、  Ｆｅｎｓｔｅｒ  Ａ、  Ｓｐｅｎｃｅ  ＪＤ、  “３Ｄ  ｕｌｔｒａｓｏｕｎｄ  ｍｅａｓｕｒｅｍｅｎｔ  ｏｆ  ｃｈａｎｇｅ  ｉｎ  ｃａｒｏｔｉｄ  ｐｌａｑｕｅ  ｖｏｌｕｍｅ：  ａ  ｔｏｏｌ  ｆｏｒ  ｒａｐｉｄ  ｅｖａｌｕａｔｉｏｎ  ｏｆ  ｎｅｗ  ｔｈｅｒａｐｉｅｓ．”、  Ｓｔｒｏｋｅ  ２００５、  ３６（９）：１９０４−１９０９．
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、検査者が超音波画像上の血管を手動で指定する必要があるので、正確に血管が検出されない場合があるという問題がある。【０００５】
そこで、本発明は、より正確に検査対象の血管を検出する超音波診断装置を提供する。
【課題を解決するための手段】
【０００６】
本発明の一態様に係る超音波診断装置は、超音波プローブにより被検体から取得された反射超音波に基づいて前記被検体の対象血管を検出する超音波診断装置であって、前記反射超音波に基づいて、前記被検体の断層画像を生成する断層画像生成部と、前記反射超音波に基づいて、前記断層画像上での前記被検体の血流領域を示す血流情報を生成する血流情報生成部と、前記血流情報生成部が生成した前記血流情報を解析することで、前記血流領域が前記対象血管に対応する血流領域であるか否かの判定を行う血流領域判定部と備える。【０００７】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムおよび記録媒体の任意な組み合わせで実現されてもよい。
【発明の効果】
【０００８】
本発明の超音波診断装置は、より正確に検査対象の血管を特定することができる。
【図面の簡単な説明】
【０００９】
【図１】図１は、実施の形態１における超音波診断装置の概略構成を示す図である。
【図２】図２は、実施の形態１における超音波診断装置の詳細構成を示すブロック図である。
【図３】図３は、実施の形態１における超音波診断装置の血管輪郭の抽出処理のフローチャートの第一例である。
【図４】図４は、（ａ）血流画像の血流領域を示す図、及び、（ｂ）血流点を示す図である。
【図５】図５は、（ａ）複数フレームにおいて血流点を抽出した結果を示す図、（ｂ）血流点をグループ化した結果を示す図、（ｃ）抽出された対象血管の血流領域相当の血流グループを示す図、及び、（ｄ）抽出された対象血管の血流領域相当の血流点を示す図である。
【図６】図６は、実施の形態１における超音波診断装置の血管輪郭の抽出処理のフローチャートの第二例である。
【図７】図７は、実施の形態２における超音波診断装置の概略構成を示すブロック図である。
【図８】図８は、実施の形態３における超音波診断装置の概略構成を示すブロック図である。
【図９】図９は、実施の形態３おける超音波診断装置の血管輪郭の抽出処理のフローチャートである。
【図１０】図１０は、各実施の形態に係る超音波診断装置の外観図である。
【図１１】図１１は、想定される超音波診断装置の第一例を示す図である。
【図１２】図１２は、想定される超音波診断装置の第二例を示す図である。
【発明を実施するための形態】
【００１０】
（本発明の基礎となった知見）
本発明者は、「背景技術」の欄において記載した超音波診断装置に関し、以下の問題が生じることを見出した。
【００１１】
近年、動脈硬化及び血管疾患等を早期に発見するために、超音波診断装置が用いられている。具体的には、超音波診断装置を用いて、血管壁の内膜と中膜とをあわせた膜の厚さである内膜中膜複合体厚（Ｉｎｔｉｍａ  Ｍｅｄｉａ  Ｔｈｉｃｋｎｅｓｓ、以下ＩＭＴ）が計測されている。また、超音波診断装置を用いて、血管の内腔が狭窄することによって発生するプラークの有無が確認されている。これは動脈硬化が進行するほどＩＭＴが厚くなったり、プラークが形成されたりすることが明らかになってきたからである。なお、動脈硬化は全身的に進行すると考えられており、主に表在性の頸動脈が動脈硬化の進行度合いを判断する際の計測対象となっている。ここでプラークとは、血管の内壁が局所的に血管の内側（内腔）に突出した隆起性病変を意味している。プラークは、血栓、脂肪性、及び繊維性など様々な形態をとり、頸動脈の狭窄及び閉塞、並びに、脳梗塞及び脳虚血を起こす原因になる恐れがある。
【００１２】
超音波診断装置によるＩＭＴの計測、および、プラーク形状の検出は、血管の超音波画像を用いて行われる。超音波診断装置は、探触子を介して被検体表より被検体内に超音波を送波し、被検体から発生する反射波に基づき超音波画像（例えばＢモード画像）を形成する。そして、検査者が、この超音波画像を見てプラーク有無を判断する。
【００１３】
ここで、検査者がＩＭＴの厚みやプラークの有無を判断する際には、検査者は、超音波画像内の血管外膜及び血管内膜の形状を手動で指定（スケッチ）し、このスケッチを元に診断を行う（例えば、図１１の構成）。具体的には、検査者は、Ｂモード画像に対し、外膜の辺縁形状を外膜輪郭線としてスケッチし、さらに内腔の辺縁形状を内腔輪郭線としてスケッチする。そして最終的にこのスケッチの形状をもとに、検査者は、プラークの有無などの診断を行う（例えば、非特許文献１参照）。
【００１４】
しかしながら、非特許文献１の手法では、画像取得後にオフライン状態において、検査者は、超音波画像上に血管壁の外膜輪郭及び内腔輪郭の位置を手動で指定する必要がある。つまり、検査者は、従来の診断方法によってプラーク検出を行う場合には手動で輪郭を指定する必要があり、手間がかかる。その結果、検査者に依存して発生する検査結果の揺らぎが生ずる。よって、正確に血管が検出されない場合があるという問題がある。
【００１５】
そこで、本発明は、より正確に検査対象の血管を検出する超音波診断装置を提供する。
【００１６】
上記の問題を解決するために、本発明の一態様に係る超音波診断装置は、超音波プローブにより被検体から取得された反射超音波に基づいて前記被検体の対象血管を検出する超音波診断装置であって、前記反射超音波に基づいて、前記被検体の断層画像を生成する断層画像生成部と、前記反射超音波に基づいて、前記断層画像上での前記被検体の血流領域を示す血流情報を生成する血流情報生成部と、前記血流情報生成部が生成した前記血流情報を解析することで、前記血流領域が前記対象血管に対応する血流領域であるか否かの判定を行う血流領域判定部とを備える。
【００１７】
これによれば、超音波診断装置は、血流領域に関する血流情報に基づいて、超音波画像において対象血管に対応する血流領域を判定することができる。超音波画像上には、対象血管に対応する血流領域を含む複数の血流領域が存在している場合がある。そのような場合に、血流情報に基づいて、複数の血流領域のうち対象血管に対応する血流領域を特定することができる。よって、超音波診断装置は、より正確に検査対象の血管を検出することができる。
【００１８】
例えば、前記血流領域判定部は、前記血流情報として、前記断層画像における前記血流領域の位置を解析することで前記判定を行う。
【００１９】
これによれば、超音波診断装置は、超音波画像における血流領域の位置に基づいて、血流領域を特定することができる。検査者が超音波プローブを移動（スキャン）させながら順次超音波画像を取得するので、超音波画像における血流領域の位置を解析することで、複数の超音波画像において、同一の血管に対応する血流領域を検出することができる。よって、より正確に検査対象の血管を検出することができる。
【００２０】
例えば、前記血流領域判定部は、前記血流領域を含む前記断層画像が取得された前記被検体上の位置に対応させて前記血流領域を配置した場合に、配置した前記血流領域の間の距離が閾値以下である複数の血流領域をまとめて１つの血流領域群として抽出し、抽出した前記血流領域群の属性に基づいて、前記血流領域群に含まれる血流領域に対して前記判定を行う。
【００２１】
これによれば、超音波診断装置は、複数の超音波画像のそれぞれにおける血流領域の位置が閾値以下である血流領域を同一の血管に対応する血流領域（血流領域群）として抽出することができる。検査者が超音波プローブ移動させながら順次超音波画像を取得するので、比較的近い時刻に取得された複数の超音波画像において、同一の血管に対応する血流領域は、比較的近い位置に存在する。超音波画像における血流領域の位置に基づいて血流領域を特定することにより、より正確に検査対象の血管を検出することができる。
【００２２】
例えば、前記血流領域判定部は、抽出した前記血流領域群の属性があらかじめ定められた前記対象血管の属性に一致する場合に、抽出した前記血流領域群に含まれる血流領域が前記対象血管に対応する血流領域であるとの前記判定を行う。
【００２３】
これによれば、超音波診断装置は、同一の血管に対応する血流領域群の属性が、あらかじめ定められた対象血管の属性と一致する場合に、当該血流領域群を対象血管に対応する血流領域として特定することができる。
【００２４】
例えば、前記血流領域判定部は、抽出した前記血流領域群に含まれる血流領域の数と、前記血流領域群に含まれる血流領域のうち面積が最大である血流領域の面積と、前記血流領域群を含む複数の断層画像が取得された被検体上の範囲のうちの始端及び終端から取得された断層画像のそれぞれにおける血流領域の位置とのうちの少なくとも１つを、前記血流領域群の属性として、前記判定を行う。
【００２５】
これによれば、これらの属性を利用することにより、より正確にＹ字形状の血管を特定することができる。血流領域群に含まれる血流領域のうち面積が最大である血流領域は、Ｙ字形状の分岐部に対応し、始端と終端とは、それぞれ、超音波プローブによるスキャンの開始時と終了時とに取得された超音波画像に対応する。つまり、超音波診断装置は、Ｙ字形状の分岐部と端部とを特定することにより、より正確にＹ字形状の血管を特定することができる。
【００２６】
例えば、前記血流領域判定部は、１つの前記血流領域群が含む前記血流領域の数が所定数以下の場合には、当該血流領域群を除く血流領域群に含まれる血流領域に対して、前記判定を行う。
【００２７】
これによれば、超音波診断装置は、含まれる血流領域の数が少ない血流領域群を血流ノイズと特定し、血流ノイズと特定された血流領域群を、血流領域の判定対象から除外することができる。よって、血流ノイズなど対象血管以外の血流情報がある場合にも正確に対象血管を検出することができる。
【００２８】
例えば、前記血流領域判定部は、前記血流領域群のうち、面積が最大である血流領域を含む血流領域群に含まれる血流領域が、前記対象血管に対応する血流領域であるとの判定を行う。
【００２９】
これによれば、特にＹ字形状の分岐部に対応する部分を有する血管を対象血管として特定することができる。
【００３０】
例えば、前記血流領域判定部は、さらに、抽出した前記血流領域群のそれぞれを含む断層画像が取得された前記被検体の位置に対応させて、前記血流領域群のそれぞれを配置した場合に、第一血流領域群を補間することで形成される図形の端部と、第二血流領域群を補間することで形成される図形の一部との距離が所定値以下であるとき、前記第一血流領域群と前記第二血流領域群とをまとめて１つの血流領域群として新たに抽出し、新たに抽出した血流領域群に対して前記判定を行う。
【００３１】
例えば、前記血流領域判定部は、さらに、抽出した前記血流領域群のそれぞれを含む前記断層画像が取得された前記被検体上の位置に対応させて前記血流領域群のそれぞれを配置した場合に、前記血流領域群に対応する前記被検体上の範囲のうちの始端又は終端から取得された断層画像に含まれる第一血流領域の前記断層画像上における第一位置と、前記被検体上の範囲のうちの始端又は終端から取得された断層画像に含まれる血流領域であって前記第一血流領域と異なる血流領域である第二血流領域の前記断層画像上における第二位置との差が所定値以下であり、かつ、前記第一血流領域が取得された前記被検体の位置と前記第二血流領域が取得された前記被検体の位置とが所定距離以下であるときに、前記第一血流領域を含む前記血流領域群と、前記第二血流領域を含む前記血流領域群とをまとめて１つの血流領域群として新たに抽出し、新たに抽出した血流領域群に対して前記判定を行う。
【００３２】
これによれば、超音波診断装置は、異なる血流領域群として抽出された複数の血流領域群を、まとめて１つの血流領域群とし、対象血管であるか否かの判定を行うことができる。
【００３３】
例えば、前記血流領域判定部は、前記血流情報として、前記断層画像における前記血流領域の面積を解析することで前記判定を行う。
【００３４】
これによれば、超音波診断装置は、超音波画像における血流領域の面積に基づいて、血流領域を特定することができる。検査者が超音波プローブを移動（スキャン）させながら順次超音波画像を取得するので、超音波画像における血流領域の面積を解析することで、複数の超音波画像において、同一の血管に対応する血流領域を検出することができる。よって、より正確に検査対象の血管を検出することができる。
【００３５】
例えば、前記血流領域判定部は、取得された複数の断層画像のそれぞれにおける前記血流領域の面積の変動を解析し、拍動様の変動の有無を判定することで前記判定を行う。
【００３６】
これによれば、超音波診断装置は、拍動様に面積が変動する動脈を、対象血管として特定することができる。
【００３７】
例えば、前記血流領域判定部は、前記血流情報として、前記血流領域の血流の方向を解析することで前記判定を行う。
【００３８】
例えば、前記血流領域判定部は、前記血流領域の前記血流の方向が、あらかじめ定められた前記対象血管の血流の方向と一致する場合に、前記血流領域が前記対象血管に対応する前記血流領域であるとの前記判定を行う。
【００３９】
これによれば、超音波診断装置は、血流の方向に基づいて対象血管を特定することができる。対象血管の血流の方向又は向きと、対象血管のおおまかな位置とは従来知られているので、超音波診断装置は、解析された血流の方向に基づいて、対象血管を特定することができる。
【００４０】
例えば、前記超音波診断装置は、さらに、前記血流領域判定部における前記判定の結果を示す情報を表示する表示部と、前記表示部が表示している前記判定の結果を示す情報に対するユーザによる修正要求に基づいて前記血流領域判定部における前記判定結果を修正する修正部とを備える。
【００４１】
これによれば、超音波診断装置は、血流領域判定部による判定結果に誤りがある場合に、検査者（ユーザ）からの修正要求に基づいて判定結果を修正し、より正確に対象血管を検出することができる。
【００４２】
例えば、前記超音波診断装置は、さらに、前記血流領域判定部が判定した血流領域に基づいて、仮血管輪郭を設定する仮血管輪郭設定部と、前記仮血管輪郭設定部が設定した前記仮血管輪郭を用いて、前記断層画像生成部が生成した前記断層画像内において、前記対象血管の輪郭を抽出する血管輪郭抽出部とを備える。
【００４３】
これによれば、超音波診断装置は、血流領域判定部が検出した対象血管の血管壁の輪郭を超音波画像上でトレースすることができる。
【００４４】
例えば、前記超音波診断装置は、さらに、前記超音波プローブの位置又は姿勢を示す位置姿勢情報を取得するプローブ位置姿勢取得部と、前記プローブ位置姿勢取得部が取得した前記位置姿勢情報と、前記血流情報生成部が生成した前記血流情報とに基づいて、前記血流領域を３次元空間上で示す血流３次元情報を生成する血流３次元構成部とを備え、前記血流領域判定部は、前記血流３次元構成部が生成した前記血流３次元情報を解析することで、前記判定を行う。
【００４５】
これによれば、超音波診断装置は、３次元空間における血流領域の形状に基づいて、対象血管を検出することができる。
【００４６】
例えば、前記対象血管は、頸動脈である。
【００４７】
これによれば、超音波診断装置は、頸動脈を検出することができる。
【００４８】
本発明の一態様に係る血管抽出方法は、超音波プローブにより被検体から取得された反射超音波に基づいて前記被検体の対象血管を検出する血管検出方法であって、前記反射超音波に基づいて、前記被検体の断層画像を生成する断層画像生成ステップと、前記反射超音波に基づいて、前記断層画像上での前記被検体の血流領域を示す血流情報を生成する血流情報生成ステップと、前記血流情報生成ステップで生成した前記血流情報を解析することで、前記血流領域が前記対象血管に対応する血流領域であるか否かの判定を行う血流領域判定ステップとを含む。
【００４９】
これにより、上記の超音波診断装置と同様の効果を奏する。
【００５０】
なお、これらの包括的または具体的な態様は、システム、方法、集積回路、コンピュータプログラムまたはコンピュータ読み取り可能なＣＤ−ＲＯＭなどの記録媒体記録媒体で実現されてもよく、システム、方法、集積回路、コンピュータプログラムまたは記録媒体の任意な組み合わせで実現されてもよい。
【００５１】
以下、実施の形態について、図面を参照しながら説明する。
【００５２】
なお、以下で説明する実施の形態は、いずれも包括的または具体的な例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、本発明を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。
【００５３】
（実施の形態１）
従来、手動の手間を減らして検査者の負荷を低減し、また、検査時間を短縮することを目的として、血流情報が存在する位置をもとに仮血管壁輪郭を設定し、設定した仮血管壁輪郭を初期輪郭として血管壁輪郭を自動的に抽出する手法がある（図１２）。しかしながら、多様な血流情報が存在する場合に、検査対象とする血管（以下、対象血管、又は検査対象血管とも呼ぶ）（例えば頸動脈）のみを抽出することが難しいという問題がある。
【００５４】
例えば、計測範囲に頸動脈以外の血管（頸静脈、又は、椎骨動脈など）が存在する場合、その血管の位置に血流が存在するので、その血管を誤って抽出することがある。また、組織の動きなどが誤って血流情報として捉えられることもある（血流ノイズ）。このように、確実に対象血管のみを抽出することが難しいという問題がある。
【００５５】
本実施の形態に係る超音波診断装置の概略構成を説明する。図１は、本実施の形態に係る超音波診断装置１５０の概略構成を示すブロック図である。
【００５６】
図１に示される超音波診断装置１５０は、超音波プローブにより被検体から取得された反射超音波２０１に基づいて対象血管を検出する。なお、検査者が被検体上で超音波プローブを直線的に移動させながら、順次、反射超音波を取得する場合を例として説明する。ただし、超音波プローブの移動は、直線的なものに限られない。つまり、操作者が曲線的に超音波プローブを移動させる場合にも適用することができる。
【００５７】
超音波診断装置１５０は、Ｂモード画像生成部１０４と、血流画像生成部１０５と、血流領域判定部１０６と、仮血管輪郭設定部１０７と、血管輪郭抽出部１０８とを備える。
【００５８】
Ｂモード画像生成部１０４は、反射超音波２０１に基づいてＢモード画像２０２を生成する。
【００５９】
血流画像生成部１０５は、血流情報生成部に相当する。血流画像生成部１０５は、反射超音波２０１に基づいて、血流がある領域を示す血流画像２０３を生成する。血流画像２０３は、血流情報に相当する。以降、血流情報の具体例として血流画像２０３を用いて説明するが、血流情報は、血流画像２０３に限られない。血流情報は、被検体において血流がある領域を示す情報であればなんでもよい。
【００６０】
血流領域判定部１０６は、血流画像２０３に基づいて、対象血管の血流領域２０４を抽出する。
【００６１】
仮血管輪郭設定部１０７は、血流領域２０４を基に仮血管輪郭２０５を設定する。
【００６２】
血管輪郭抽出部１０８は、仮血管輪郭２０５を初期輪郭として、Ｂモード画像２０２を用いて血管輪郭（血管の外膜輪郭）２０６を抽出する。具体的には、例えば、血管輪郭抽出部１０８は、仮血管輪郭２０５を初期輪郭とした探索処理をＢモード画像２０２上で行うことにより、血管輪郭２０６を抽出する。
【００６３】
以下、本実施の形態に係る超音波診断装置１５０の構成を詳細に説明する。
【００６４】
図２は、本実施の形態に係る超音波診断装置１５０の詳細構成を示すブロック図である。
【００６５】
図２に示す超音波診断装置１５０は、制御部１０２と、送受信部１０３と、Ｂモード画像生成部１０４と、血流画像生成部１０５と、血流領域判定部１０６と、仮血管輪郭設定部１０７と、血管輪郭抽出部１０８と、血管輪郭画像生成部１０９と、データ格納部１１０とを備える。また、超音波診断装置１５０の外部には、プローブ１０１と、表示部１１１とがあり、それぞれ、超音波診断装置１５０に接続されている。なお、プローブ１０１と、表示部１１１とは、超音波診断装置１５０の内部にあってもよい。また、プローブ１０１と、表示部１１１とがなくてもよい。
【００６６】
プローブ１０１は、超音波を送受信する超音波振動子を備える超音波プローブである。プローブ１０１は、送受信部１０３の指示に従い超音波を送受信する。また、プローブ１０１は、被検体からの反射超音波２０１（超音波反射信号）をエコー信号として受信する。なお、プローブ１０１は、超音波振動子が１次元方向に配列されているプローブであってもよいし、超音波振動子がマトリックス状に配置された２次元アレイプローブであってもよい。
【００６７】
制御部１０２は、超音波診断装置１５０に含まれる各処理部の制御を行う。以降、特に明記しないが、制御部１０２が各処理部の動作を制御する。例えば、制御部１０２は、動作タイミングなどを制御しながら各処理部に処理を実行させる。
【００６８】
送受信部１０３は、プローブ１０１の超音波振動子を駆動させて超音波を発生させる。また、送受信部１０３は、プローブ１０１が被検体から受信した反射超音波２０１を受信する。
【００６９】
Ｂモード画像生成部１０４は、送受信部１０３が受信した反射超音波２０１を基に、Ｂモード画像２０２を生成する。具体的には、Ｂモード画像生成部１０４は、反射超音波２０１に対しフィルタ処理を行ったのち、包絡線検波を行う。さらに、Ｂモード画像生成部１０４は、包絡線検波により取得した信号に対して対数変換及びゲイン調整を行うことでＢモード画像２０２を生成する。
【００７０】
血流画像生成部１０５は、送受信部１０３が受信した反射超音波２０１を基に、血流画像２０３を生成する。ここで、血流画像２０３とは、血流がある領域を示す画像である。具体的には、血流画像生成部１０５は、超音波が血流に反射されることにより生ずる周波数の変化を用いて、血管内の血流の速度（血流速度）を検出する。そして、血流画像生成部１０５は、検出した血流速度を、カラーデータとして画像化することで血流画像２０３を生成する。なお、血流速度を画像化する方法としては、例えばカラードプラ法又はパワードプラ法を用いることができる。
【００７１】
血流領域判定部１０６は、血流画像生成部１０５が生成した血流画像に基づいて、対象となる血管の血流領域２０４を抽出する。
【００７２】
なお、ここでは、血流画像生成部１０５は、血流画像２０３を生成するとしたが、必ずしも画像を生成する必要はない。つまり、血流画像生成部１０５が、血流がある領域を示す情報（血流情報）を生成し、血流領域判定部１０６がこの血流情報を用いて対象となる血管の血流領域２０４を抽出してもよい。この血流領域の判定方法については、後ほど詳述する。
【００７３】
仮血管輪郭設定部１０７は、血流領域判定部１０６が抽出した血流領域２０４に基づいて仮血管輪郭２０５を設定する。そして、仮血管輪郭設定部１０７は設定した仮血管輪郭２０５を示す仮血管輪郭情報を血管輪郭抽出部１０８に送る。
【００７４】
血管輪郭抽出部１０８は、仮血管輪郭情報で示される仮血管輪郭２０５を、Ｂモード画像２０２に設定した上で、この仮血管輪郭２０５を初期輪郭として用い、Ｂモード画像２０２からより詳細な血管輪郭２０６を示す情報を抽出する。ここでは、血管輪郭抽出部１０８は、血管の外膜の輪郭（外膜輪郭）に相当する輪郭を抽出する。
【００７５】
なお、本実施の形態では、血管輪郭抽出部１０８が、外膜輪郭を抽出するようにＢモード画像からの抽出を行っているが、仮血管輪郭情報を血流画像２０３に設定した上で、血流画像２０３から内腔輪郭を抽出するようにしてもよい。また、その両方を行って外膜輪郭と内腔輪郭とを抽出するようにしてもよい。
【００７６】
血管輪郭画像生成部１０９は、血管輪郭抽出部１０８が抽出した血管輪郭を示す情報を、Ｂモード画像２０２に重畳する形で合成し、血管輪郭画像を生成する。
【００７７】
データ格納部１１０は、Ｂモード画像生成部１０４が生成するＢモード画像２０２、血流画像生成部１０５が生成する血流画像２０３、及び、血管輪郭抽出部１０８が生成する血管輪郭２０６を格納する。
【００７８】
表示部１１１は、Ｂモード画像２０２、血流画像２０３若しくは血管輪郭２０６、又は、それらのデータなどを表示する。表示部１１１は、ＬＣＤ（Liquid crystal display、液晶ディスプレイ）などの表示装置である。なお、本実施の形態は、対象とする血管の輪郭をより正確に求める輪郭抽出方法を特徴とする。よって、超音波診断装置１５０が、血管輪郭画像生成部１０９及びデータ格納部１１０を具備するかどうかは任意である。
【００７９】
図３は、本実施の形態における超音波診断装置１５０の血管輪郭２０６の抽出処理のフローチャートの第一例である。
【００８０】
まず、ステップＳ２０１において、Ｂモード画像生成部１０４がＢモード画像２０２を生成するとともに、血流画像生成部１０５が血流画像２０３を生成する。具体的には、送受信部１０３は、プローブ１０１を通じて被検体に超音波を発信し、プローブ１０１を通じて反射超音波２０１を受信する。Ｂモード画像生成部１０４及び血流画像生成部１０５は、送受信部１０３が受信したデータを処理することで、Ｂモード画像２０２及び血流画像２０３を生成し、生成したＢモード画像２０２及び血流画像２０３をデータ格納部１１０に保存する。なお、Ｂモード画像生成、および、血流画像生成は時系列的に連続で行ってもよい（以降便宜上、Ｂモード画像、および、血流画像の１生成単位を、フレームと呼ぶこととする）。
【００８１】
次に、ステップＳ２０２において、血流領域判定部１０６は、血流画像２０３を解析し、対象血管の血流領域２０４を抽出する。以下、この工程について述べる。
【００８２】
本実施の形態では、血流領域判定部１０６は、まず、予め定めた値より大きい面積を有する血流領域を抽出する。血流画像２０３には、血流画像２０３を生成する過程において誤って血流領域と判定されうる小領域（小さい面積を有する領域）が存在する場合がある。よって、このような処理を行うことで、血流画像２０３の中から効率よく頸動脈以外の血流領域を取り除くことができる。さらに、血流領域判定部１０６は、血流領域の重心の位置を示す座標を算出し、これをその血流領域を示す代表点（以降、血流点と呼ぶ）とする。ここまでで、あるフレームにおける血流点を抽出したことになる。なお、血流点は１点も抽出されないこともあり、また、複数抽出されることもある。また、ここでは血流領域の重心位置を代表点としたが、血流領域範囲の最大座標値と最小座標値との中央値で示される点を代表点としてもよい。さらに、各血流点の情報として、代表点の座標の代わりに、血流領域の面積を示す値を算出してもよい。また、それらを併用してもよい。
【００８３】
図４を用いて、血流画像から血流点を抽出する処理について説明する。図４の（ａ）は、血流画像の血流領域４０１及び４０２を表示した図である。図４の（ｂ）は、血流領域４０１及び４０２のそれぞれに対応した血流点４１１及び４１２を示す図である。ここでは、それぞれの血流領域の重心の位置を血流点としている。
【００８４】
次に、血流領域判定部１０６は、血流点のグループ化を行う。血流点のグループ化とは、同じ特徴を有する血流点同士を同じラベルをもつグループとすることを意味する。なお、ラベルとはグループを識別する指標のことであり、名前、数字、色など、グループを一意に識別できるのであれば、どのような指標でも構わない。本実施の形態では、血流領域判定部１０６は、あるフレームにおける血流点に対し、当該フレームの前後の複数フレームにおいて、フレーム内における血流位置が近い血流点に同じ数字（以降、血流グループ番号と呼ぶ）を付けることで、グループ化を実現する。なお、グループ化された血流点の集まりを血流グループとよぶことがある。
【００８５】
具体的には、血流領域判定部１０６は、あるフレームにおけるある血流点を注目血流点とすると、当該フレームより過去のフレームであって、当該フレームが取得された時刻から１秒以内に取得されたフレームに着目する。血流領域判定部１０６は、着目したフレームにおいて、注目血流点との座標値の距離が１０ｍｍ相当以内である血流点が存在するならば、この過去の血流点につけられた血流グループ番号と同じ番号を、注目血流点に対して番号付ける。また、前記範囲内に血流点が存在しないならば、血流領域判定部１０６は、これまでにつけられていない、新しい血流グループ番号を、注目血流点に対して番号付ける。血流領域判定部１０６は、この処理を、取得した全てのフレームに対して繰り返し適用することで、抽出した全ての血流点に対して血流グループ番号を番号付けることができる。このようにして、血流領域判定部１０６は、全ての血流点をグループ化することができる。なお、上記で「１秒」又は「１０ｍｍ」というように具体的な数値を示したが、あくまで例示であり、それぞれ、他の時間、他の長さであってもよい。
【００８６】
次に、血流領域判定部１０６は、血流点の解析を行い、対象血管血流領域に相当する血流点を抽出する。本実施の形態では、血流グループを単位とした解析、および、抽出を行う。ここでの解析は、例えばパターンマッチングである。例えば、総頸動脈から内頸動脈及び外頸動脈に至る頸動脈の形状は、分岐部をもつＹ字型であるという形状的特徴を有することが知られている。血流領域判定部１０６は、この形状的特徴に当てはまる血流グループの組を抽出する。血流領域判定部１０６が抽出した血流グループに属する血流点が、対象血管の血流領域２０４に相当する。
【００８７】
図５は、縦軸及び横軸をそれぞれ各フレームの縦軸及び横軸とし、奥行き方向の軸をフレームの前後の方向（以降、フレーム方向とよぶ）の軸として、模式的に立体図とした説明図である。
【００８８】
図５の(ａ)は、血流領域判定部１０６が各フレームで血流点を抽出した結果を示している。図５の(ａ)において、一番手前に配置されたフレームでは、３つの血流点５０１、５０２及び５０３が抽出されている。この時点では、いずれの血流点が対象血管に相当する血流点であるかはわからない。
【００８９】
図５の（ｂ）は、血流領域判定部１０６が図５の（ａ）のように抽出した血流点をフレーム方向にグループ化した結果の例である。図５の（ｂ）に示される各円柱は、同じ血流グループ番号をもつ血流点の集まりを示したもので、血流グループに相当する。図５の（ｂ）には、４つの血流グループ５１１、５１２、５１３及び５１４が存在している。血流領域判定部１０６は、頸動脈を対象血管とする場合には、例えばＹ字型の特徴をもつパターンをリファレンスパターンとして、パターンマッチングを実行する。このパターンマッチング処理により、図５の（ａ）に示される４つの血流グループから、リファレンスパターンと同様の特徴を有する血流グループ５１１、５１２及び５１４を選ぶことができる（図５の（ｃ））。その結果、血流領域判定部１０６は、対象血管血流領域として、血流グループ５１１、５１２及び５１４に属する血流点を抽出することができる。
【００９０】
図５の（ｄ）は、対象血流領域として、血流点５０１及び５０２が抽出されたことを示している。図５の（ａ）において血流領域判定部１０６が抽出した血流点５０３は、血流グループ５１３に属するので、対象血管ではないとして適切に除去することができている。
【００９１】
なお、ここでは、血流領域判定部１０６がパターンマッチング技術を用いる例を示したが、対象血管の特徴を学習させた機械学習ベースの抽出方式を採用してもよい。また、血流領域判定部１０６は、各グループの位置関係により形作られる形状的な特徴だけでなく、血流面積などの情報を用いてもよい。たとえば、頸動脈の分岐部では血流領域面積が非常に大きくなる傾向がある。そのため、血流領域判定部１０６は、複数のグループのそれぞれにおいて血流領域の面積が最大となる血流点を求めた後、複数のグループのうち最大血流面積を有するグループが分岐部を含んでいる領域であるとして抽出してもよい。さらに、血流領域判定部１０６は、分岐部を含む血流グループを抽出した後、その血流グループとの位置関係に基づいて、さらに複数の血流グループを抽出してもよい。これにより、基準となる血流グループが先に決まるので、より効率よく対象血流グループを抽出することができる。
【００９２】
血流領域判定部１０６は、他の例として、血流面積を用いて、拍動様の面積変化を特徴として対象血管の血流グループを抽出してもよい。一般に、動脈の方が血管径の変化が大きいので、血流面積の変化が大きくなるという特徴がある。血流領域判定部１０６は、動脈と静脈とを区別することを意図する場合には、この特徴を利用することより、より適切に対象血管を抽出することができる。
【００９３】
また、血流領域判定部１０６は、パターンマッチング等の処理を行う前に、血流グループ自体を解析してもよい。例えば、血流のノイズは、フレーム方向に間欠的に描出されることが多く、血流点数の少ない個別のグループに分かれることがある。そこで、血流グループのうち、そのグループに所属する血流点が閾値より少ない場合には、その血流グループは血流ノイズであるとして、以降の解析対象から外すようにしてもよい。これによって、不要な情報を処理する必要がなくなり、より効率よく対象血流グループを抽出することができる。
【００９４】
図３に戻り、次に、ステップＳ２０３では、仮血管輪郭設定部１０７は、ステップＳ２０２で抽出した対象血管の血流領域２０４を基に仮血管輪郭２０５を設定する。血管輪郭は、抽出した対象血流領域周辺に存在すると想定されるので、対象血管の血流領域２０４の位置に仮の血管輪郭を設定する。本実施の形態では、血管を輪切りにした像が描出されるようにプローブをスキャンして円形状の血管像が取得されることを想定し、仮血管輪郭は円形状で設定する。
【００９５】
まず、仮血管輪郭設定部１０７は、各フレームにおいて、対象血流領域に相当する血流点があるかどうかを判定する。仮血管輪郭設定部１０７は、対象血流領域に相当する血流点がある場合には、その血流点を中心とした円を仮血管輪郭とする。このとき、仮血管輪郭設定部１０７は、仮血管輪郭とした円が対象血流領域を内包するように、円の半径を決定することが望ましい。また、仮血管輪郭設定部１０７は、対象とする血管の半径の統計的な平均値をもとに円の半径を決定してもよい。
【００９６】
ステップＳ２０４では、血管輪郭抽出部１０８は、ステップＳ２０３で設定した仮血管輪郭２０５を基に、Ｂモード画像２０２から血管輪郭２０６を抽出する。本実施の形態では、血管輪郭抽出部１０８は、仮血管輪郭を初期輪郭として使用して、Ｂモード画像２０２に対する動的輪郭探索処理（スネークスなど）を行った結果得られる輪郭を、血管輪郭２０６に決定する。ここで、動的輪郭探索処理とは、エネルギー最小化処理を行うことにより初期輪郭の輪郭点を移動させることで輪郭を抽出する処理である。スネークスアルゴリズムでは、例えば下記（式１）、（式２）及び（式３）で定義されるエネルギーＥｓｎａｋｅｓを最小化するように輪郭を決定する。
【００９７】
【数１】
0006020470.tif 000002
【００９８】
【数２】
0006020470.tif 000003
【００９９】
【数３】
0006020470.tif 000004
【０１００】
ここで、Ｅｉｎｔは輪郭線の内部変形エネルギーであり、Ｅｉｍａｇｅは輪郭線と画像との適合度を表す画像エネルギーである。ｖは輪郭線のパラーメータ表現であり、ｖｓはｖの１階微分、ｖｓｓはｖの２階微分を示している。α、β、ｗ１、ｗ２は重みを示す定数である。Ｇσはガウシアンフィルタであり、∇２はラプラシアンフィルタであり、「＊」は畳み込み（Convolution）の演算子であり、Ｉは画像の輝度値である。より具体的には、スネークスアルゴリズムは、輪郭線を、当該輪郭線を離散化した輪郭点として表現し、輪郭点ごとにエネルギーＥｓｎａｋｅｓが最小化されるような点を決定する。例えば、α＝０．８、β＝０．２とすることで、もとの円形形状を保ちながら輪郭を探索することができる。
【０１０１】
ステップＳ２０５では、血管輪郭画像生成部１０９は、血管輪郭を示す情報をＢモード画像に重畳することで血管輪郭画像を生成する。
【０１０２】
ステップＳ２０６では、表示部１１１は、生成された画像を表示する。
【０１０３】
なお、前述したようにステップＳ２０５（血管輪郭画像生成ステップ）、及び、ステップＳ２０６（表示ステップ）を採用するかどうかは任意である。
【０１０４】
図６は、本実施の形態における超音波診断装置１５０の血管輪郭２０６の抽出処理のフローチャートの第二例である。超音波診断装置１５０は、図６に示されるように、ステップＳ２０１〜Ｓ２０４の処理のみを行なってもよい。
【０１０５】
以上より、本実施の形態に係る超音波診断装置１５０は、血流領域の分布から対象血管の血流領域を求め、この血流領域に基づいて対象血管の輪郭を抽出する。これにより、超音波診断装置１５０は、より安定、かつより正確に輪郭位置情報を得ることができる。その結果、超音波診断装置１５０は、正しく抽出対象の血管壁の輪郭をトレースすることができる。
【０１０６】
（実施の形態２）
本実施の形態では、プローブの位置又は姿勢を示す情報を用いることにより、さらに正確に対象血管の形状を抽出することができる超音波診断装置の例について説明する。
【０１０７】
図７は、本実施の形態に係る超音波診断装置１５１の概略構成を示すブロック図である。なお、図７において、図２と同じ構成要素については同じ符号を用い、説明を省略する。
【０１０８】
図７に示される超音波診断装置１５１は、図２に示す超音波診断装置１５０の構成に加え、プローブ位置姿勢取得部１１２と、血流３次元構成部１１３とを備える。
【０１０９】
プローブ位置姿勢取得部１１２は、プローブ１０１の位置又は姿勢を示す位置姿勢情報を取得する。
【０１１０】
血流３次元構成部１１３は、プローブ１０１の位置姿勢情報と、血流画像２０３とをもとに、血流情報を３次元空間上で示す血流３次元情報を生成する。
【０１１１】
実施の形態１における超音波診断装置１５０においては、血流点の形状的特徴を解析することで、対象血管に相当する血流領域を特定した。プローブを不規則に移動させて血管をスキャンする場合には、プローブの不規則な動きが血流点の形状的特徴に、意図に反して反映される可能性がある。特に、血管の走行方向に沿う方向の、プローブの不規則な動きによる影響が大きい。例えば、血管の走行方向に沿って総頸動脈から内頸動脈及び外頸動脈が描出されるようにプローブを移動させた後、プローブの動きを反転させ、内頸動脈及び外頸動脈から総頸動脈に向けてプローブを移動させた場合には、抽出された血流点の形状は、Ｙ字を２つつなげた形となる。つまり、抽出された血流点の形状は、頸動脈の形状を正しく示すＹ字型にはならない。このように、プローブの動きが反映された形状からは、頸動脈としての形状特徴を捉えることが困難であり、正しく対象血管を抽出することができない。
【０１１２】
そこで、本実施の形態における超音波診断装置１５１は、プローブの動きの影響を除去するために、プローブの位置姿勢情報を取得する。そして、取得した位置姿勢情報に基づいて血流３次元情報を生成する。さらに、生成した血流情報をもとに、実施の形態１と同様に血流解析を行う。
【０１１３】
かかる構成によれば、プローブの動きが不規則な場合においても、対象血管の血流に対して血流情報以外の影響を受けにくくすることができる。よって、より正確に抽出対象の血管の位置を得ることができる。
【０１１４】
（実施の形態３）
本実施の形態では、上述した実施の形態１の変形例について説明する。
【０１１５】
図８は、本実施の形態に係る超音波診断装置１５２の構成を示すブロック図である。なお、図８において、図２と同じ構成要素については同じ符号を用い、説明を省略する。
【０１１６】
図８に示される超音波診断装置１５２は、図２に示される超音波診断装置１５０の構成に加え、ユーザ指定入力部１２１と、修正部１２２とを備える。
【０１１７】
ユーザ指定入力部１２１は、ユーザによる修正要求を受け付ける。
【０１１８】
修正部１２２は、ユーザ指定入力部１２１が受け付けた修正要求をもとに、血流領域判定部１０６による血流領域判定結果を修正する。
【０１１９】
図９は、本実施の形態における超音波診断装置１５２の血管輪郭の抽出処理のフローチャートである。
【０１２０】
超音波診断装置１５２は、実施の形態１と同様にステップＳ２０１からステップＳ２０６までの処理を実行する。
【０１２１】
次に、ステップＳ２１１において、ユーザ指定入力部１２１は、ユーザによる修正要求を受け付ける。
【０１２２】
次に、ステップＳ２２２において、修正部１２２は、ユーザ指定入力部１２１が受け付けた修正要求をもとに、血流領域判定部１０６による血流領域判定結果を修正する。
【０１２３】
実施の形態１においては、血流点の形状的特徴を解析することで、対象血管に相当する血流領域を特定していたが、解析結果が誤った場合に正しく対象血管を抽出することができないという課題があった。例えば、頸静脈の血流が誤って対象血流として特定された場合には、頸静脈の血管壁を抽出してしまう。このように、解析が誤るような場合には正しく対象血管を抽出することができなくなる。
【０１２４】
このため、解析が誤った場合にも正しく対象血管を抽出するために、ユーザによる修正指示を入力する手段を設け、入力された修正指示に基づいて、血流領域判定結果を修正する。さらに、修正された血流領域判定結果をもとに、実施の形態１と同様に仮血管輪郭設定、血管輪郭抽出、および、血管輪郭画像生成を行う。
【０１２５】
ユーザによる修正指示は、例えば前記血流グループを単位とし、血流グループの選択指示あるいは選択解除指示により行う。つまり、誤った血流グループが対象血流領域として選ばれている場合には、選択解除指示を行うことで、対象血流領域修正部では該血流グループを対象血流領域から除外するように前記対象血流領域を修正する。また、対象血流領域として抽出されなかった血流グループに対して選択指示を行うことで、対象血流領域修正部では該血流グループを対象血流領域として抽出するように前記対象血流領域を修正する。
【０１２６】
なお、前記対象血流領域修正部における対象血流領域の修正は、血流グループ単位でなくてもよい。
【０１２７】
かかる構成によれば、血流画像の解析が困難で対象血流領域の抽出が誤るような場合においても、誤りを修正することが可能となる。よって、ユーザの判断を反映し、より正しく抽出対象の血管を抽出することができる。
【０１２８】
以上、各実施の形態に係る超音波診断装置について説明した。各実施の形態に係る超音波診断装置の外観は、例えば、図１０に示される。
【０１２９】
なお、本発明は、各実施の形態に限定されるものではない。例えば、各実施の形態における超音波診断装置に含まれる処理部の一部又は全部が、プローブ１０１に含まれてもよい。
【０１３０】
また、各実施の形態では、Ｂモード画像及び血流画像として、走行方向に垂直な血管の断面である所謂短軸像を用いる場合を例に説明したが、本発明は、走行方向に平行な血管の断面である所謂長軸像を用いる場合にも適用できる。この場合、仮血管輪郭設定部１０７は、四角形形状の仮血管輪郭を決定する。ここで四角形形状とは、長方形、平行四辺形及び略四角形形状を意味する。
【０１３１】
また、各実施の形態に係る超音波診断装置に含まれる各処理部は典型的には集積回路であるＬＳＩとして実現される。これらは個別に１チップ化されてもよいし、一部又は全てを含むように１チップ化されてもよい。
【０１３２】
また、集積回路化はＬＳＩに限るものではなく、専用回路又は汎用プロセッサで実現してもよい。ＬＳＩ製造後にプログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、又はＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【０１３３】
また、各実施の形態に係る、超音波診断装置の機能の一部又は全てを、ＣＰＵ等のプロセッサがプログラムを実行することにより実現してもよい。
【０１３４】
さらに、本発明は上記プログラムであってもよいし、上記プログラムが記録された非一時的なコンピュータ読み取り可能な記録媒体であってもよい。また、上記プログラムは、インターネット等の伝送媒体を介して流通させることができるのは言うまでもない。
【０１３５】
また、各実施の形態に係る超音波診断装置、及びその変形例の機能のうち少なくとも一部を組み合わせてもよい。
【０１３６】
また、上記で用いた数字は、全て本発明を具体的に説明するために例示するものであり、本発明は例示された数字に制限されない。
【０１３７】
また、ブロック図における機能ブロックの分割は一例であり、複数の機能ブロックを一つの機能ブロックとして実現したり、一つの機能ブロックを複数に分割したり、一部の機能を他の機能ブロックに移してもよい。また、類似する機能を有する複数の機能ブロックの機能を単一のハードウェア又はソフトウェアが並列又は時分割に処理してもよい。
【０１３８】
また、上記のステップが実行される順序は、本発明を具体的に説明するために例示するためのものであり、上記以外の順序であってもよい。また、上記ステップの一部が、他のステップと同時（並列）に実行されてもよい。
【０１３９】
さらに、本発明の主旨を逸脱しない限り、本実施の形態に対して当業者が思いつく範囲内の変更を施した各種変形例も本発明に含まれる。
【産業上の利用可能性】
【０１４０】
本発明にかかる超音波診断装置は、血管壁の形状を正しく獲得する手段を有し、動脈硬化の診断に有用である。
【符号の説明】
【０１４１】
１０１    プローブ
１０２    制御部
１０３    送受信部
１０４    Ｂモード画像生成部
１０５    血流画像生成部
１０６    血流領域判定部
１０７    仮血管輪郭設定部
１０８    血管輪郭抽出部
１０９    血管輪郭画像生成部
１１０    データ格納部
１１１    表示部
１１２    プローブ位置姿勢取得部
１１３    血流３次元構成部
１５０、１５１、１５２    超音波診断装置
２０１    反射超音波
２０２    Ｂモード画像
２０３    血流画像
２０４    血流領域
２０５    仮血管輪郭
２０６    血管輪郭
４０１、４０２    血流領域
４１１、４１２、５０１、５０２、５０３    血流点
５１１、５１２、５１３、５１４    血流グループ
【図１】
0006020470.tif 000005
【図２】
0006020470.tif 000006
【図３】
0006020470.tif 000007
【図４】
0006020470.tif 000008
【図５】
0006020470.tif 000009
【図６】
0006020470.tif 000010
【図７】
0006020470.tif 000011
【図８】
0006020470.tif 000012
【図９】
0006020470.tif 000013
【図１０】
0006020470.tif 000014
【図１１】
0006020470.tif 000015
【図１２】
0006020470.tif 000016