(57)【要約】【課題】知りたい移動情報として、より正確な情報を得ることができる計測装置を提供する。【解決手段】一の時相のＢモード画像ＢＩと他の時相のＢモード画像ＢＩとの間で、第一関心領域ＲＯＩ１が設定された第一部分Ｐ１の移動ベクトルを算出する第一移動ベクトル算出部と、前記一の時相のＢモード画像ＢＩと前記他の時相の超音波画像ＢＩとの間で、第二関心領域ＲＯＩ２が設定された第二部分Ｐ２の移動ベクトルを算出する第二移動ベクトル算出部と、前記第一部分Ｐ１の移動ベクトルに対して、前記第二部分Ｐ２の移動ベクトルによる補正を行なって、補正済移動ベクトルを算出する補正済移動ベクトル算出部と、を備えることを特徴とする。【選択図】図４
2014087507.tif 000002
【特許請求の範囲】
【請求項１】
被検体の超音波画像における観察対象に第一関心領域を設定する第一関心領域設定部と、
前記超音波画像において、前記第一関心領域を含まない部分に第二関心領域を設定する第二関心領域設定部と、
一の時相の超音波画像と他の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルを算出する第一移動ベクトル算出部と、
前記一の時相の超音波画像と前記他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動ベクトルを算出する第二移動ベクトル算出部と、
前記第一部分の移動ベクトルに対して、前記第二部分の移動ベクトルによる補正を行なって、補正済移動ベクトルを算出する補正済移動ベクトル算出部と、
を備えることを特徴とする計測装置。
【請求項２】
被検体の超音波画像における観察対象に第一関心領域を設定する第一関心領域設定部と、
前記超音波画像において、前記第一関心領域を含まない部分に第二関心領域を設定する第二関心領域設定部と、
一の時相の超音波画像と他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動ベクトルを算出する第二移動ベクトル算出部と、
前記第二部分の移動ベクトルの逆ベクトルによって前記他の時相の超音波画像を位置補正した位置補正済超音波画像と、前記一の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルを算出する第一移動ベクトル算出部と、
を備えることを特徴とする計測装置。
【請求項３】
前記第二関心領域は、知りたい移動情報以外の移動を検出するための部分に設定されることを特徴とする請求項１又は２に記載の計測装置。
【請求項４】
前記一の時相の超音波画像と前記他の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動を追跡する第一追跡部と、
前記一の時相の超音波画像と前記他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動を追跡する第二追跡部と、
を備え、
前記第一移動ベクトル算出部は、前記第一追跡部による前記第一部分の移動の追跡に基づいて、前記第一部分の移動ベクトルを算出し、
前記第二移動ベクトル算出部は、前記第二追跡部による前記第二部分の移動の追跡に基づいて、前記第二部分の移動ベクトルを算出する
ことを特徴とする請求項１に記載の計測装置。
【請求項５】
前記位置補正済超音波画像と、前記一の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動を追跡する第一追跡部と、
前記一の時相の超音波画像と前記他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動を追跡する第二追跡部と、
を備え、
前記第二移動ベクトル算出部は、前記第二追跡部による前記第二部分の移動の追跡に基づいて、前記第二部分の移動ベクトルを算出し、
前記前記第一移動ベクトル算出部は、前記第一追跡部による前記第一部分の移動の追跡に基づいて、前記第一部分の移動ベクトルを算出する
ことを特徴とする請求項２に記載の計測装置。
【請求項６】
前記第一関心領域の移動情報として、前記補正済移動ベクトルの大きさ及び方向に応じた表示形態を有する移動情報が表示される表示部を備えることを特徴とする請求項１又は４に記載の計測装置。
【請求項７】
前記第一関心領域の移動情報として、前記第一部分の移動ベクトルの大きさ及び方向に応じた表示形態を有する移動情報が表示される表示部を備えることを特徴とする請求項２又は５に記載の計測装置。
【請求項８】
コンピュータに、
被検体の超音波画像における観察対象に第一関心領域を設定する第一関心領域設定機能と、
前記超音波画像において、前記第一関心領域を含まない部分に第二関心領域を設定する第二関心領域設定機能と、
一の時相の超音波画像と他の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルを算出する第一移動ベクトル算出機能と、
前記一の時相の超音波画像と前記他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動ベクトルを算出する第二移動ベクトル算出機能と、
前記第一部分の移動ベクトルに対して、前記第二部分の移動ベクトルによる補正を行なって、補正済移動ベクトルを算出する補正済移動ベクトル算出機能と、
を実行させることを特徴とする計測装置の制御プログラム。
【請求項９】
コンピュータに、
被検体の超音波画像における観察対象に第一関心領域を設定する第一関心領域設定機能と、
前記超音波画像において、前記第一関心領域を含まない部分に第二関心領域を設定する第二関心領域設定機能と、
一の時相の超音波画像と他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動ベクトルを算出する第二移動ベクトル算出機能と、
前記第二部分の移動ベクトルの逆ベクトルによって前記他の時相の超音波画像を位置補正した位置補正済超音波画像と、前記一の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルを算出する第一移動ベクトル算出機能と、
を実行させることを特徴とする計測装置の制御プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、超音波画像に設定された関心領域が設定された部分の移動を追跡して移動量を算出する計測装置及びその制御プログラムに関する。
【背景技術】
【０００２】
脳梗塞や心筋梗塞などの循環器系疾患を予防するには、動脈硬化の兆候を早期に察知することが有効である。動脈硬化の検査としては、超音波診断装置を用いて血管を観察することが行われている。例えば、特許文献１，２では、動脈硬化を診断するため、超音波画像において血管壁の動きを追跡して計測を行なうことが記載されている。
【０００３】
また、血管内膜の動脈硬化による部分的な肥厚であるプラークが破綻すると、心筋梗塞などを引き起こすといわれている。プラークの破綻しやすさを決める因子としては、脂質コアーの大きさ、脂質コアーを覆う繊維性皮膜の厚さなどがある。これら脂質コアーの大きさや繊維性皮膜の厚さは、プラーク内部の移動の様子を観察することによって推定することが可能である。そこで、特許文献３には、血管のプラーク性状を知るために、拍動に伴うプラーク内部における各部分の移動を追跡することが開示されている。具体的には、この特許文献３では、関心領域を分割して分割領域を設定し、各分割領域が設定された部分の移動を追跡して分割領域の移動ベクトルを表示している。この移動ベクトルにより、観察対象であるプラーク内部の移動の様子を知ることができる。
【０００４】
プラークの性状判定について具体的に説明すると、各分割領域の移動量や移動方向が一様ではない場合、そのプラークは軟らかいものであることが推定される。一方、各分割領域の移動量や移動方向が一様である場合、そのプラークは硬いものであることが推定される。拍動によるプラーク内部の移動（変形）を検出するため、特許文献３によれば、各分割領域が設定された部分の移動を追跡し、その移動量を計測してベクトル表示するので、各分割領域の移動量や移動方向を知ることができ、プラークの性状判定を行なうことができる。
【０００５】
このように、プラークの性状判定を行なうためには、拍動によるプラーク内部の移動を知ることが有用である。しかし、血管は、粘弾性体である血液が血管内を移動することに伴う慣性力や、超音波プローブを体表面上において煽る動作などによって、並進運動する。従って、前記分割領域の追跡を行なって得られた移動距離には、拍動による移動の要素のほか、前記並進運動による移動の要素も含まれてしまう。そこで、この並進運動による移動の要素を排除するために、特許文献３では、分割領域の移動量の平均を、関心領域全体の移動量として算出し、これを差し引いた移動量で各分割領域のベクトルを表示することが記載されている（特許文献３における第９ページの段落００４９）。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００２−２３８９０３号公報
【特許文献２】特開２０１０−１１０３７３号公報
【特許文献３】特開２０１２−９０８１９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、分割領域を追跡して算出された移動距離には、拍動による移動の要素も含まれるため、分割領域の移動量の平均には、本来知りたい情報である拍動によるプラーク内部の移動の要素も含まれている。従って、分割領域の移動量の平均を差し引くと、本来知りたい拍動による移動の要素まで差し引かれてしまい、拍動による移動情報を正確に得ることができないおそれがある。そこで、観察対象について知りたい移動情報以外の移動要素をより正確に排除し、知りたい移動情報として、より正確な情報を得ることができる計測装置が望まれている。
【課題を解決するための手段】
【０００８】
上述の課題を解決するためなされた一の観点の発明は、被検体の超音波画像における観察対象に第一関心領域を設定する第一関心領域設定部と、前記超音波画像において、前記第一関心領域を含まない部分に第二関心領域を設定する第二関心領域設定部と、一の時相の超音波画像と他の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルを算出する第一移動ベクトル算出部と、前記一の時相の超音波画像と前記他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動ベクトルを算出する第二移動ベクトル算出部と、前記第一部分の移動ベクトルに対して、前記第二部分の移動ベクトルによる補正を行なって、補正済移動ベクトルを算出する補正済移動ベクトル算出部と、を備えることを特徴とする計測装置である。
【０００９】
また、他の観点の発明は、被検体の超音波画像における観察対象に第一関心領域を設定する第一関心領域設定部と、前記超音波画像において、前記第一関心領域を含まない部分に第二関心領域を設定する第二関心領域設定部と、一の時相の超音波画像と他の時相の超音波画像との間で、前記第二関心領域が設定された第二部分の移動ベクトルを算出する第二移動ベクトル算出部と、前記第二部分の移動ベクトルの逆ベクトルによって前記他の時相の超音波画像を位置補正した位置補正済超音波画像と、前記一の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルを算出する第一移動ベクトル算出部と、を備えることを特徴とする計測装置である。
【発明の効果】
【００１０】
上記一の観点の発明によれば、前記第二関心領域は、観察対象に設定される前記第一関心領域を含んでおらず、前記第二関心領域が設定された第二部分の移動ベクトルによって、第一関心領域が設定された第一部分の移動ベクトルが補正されるので、観察対象について知りたい移動情報以外の移動要素をより正確に排除し、知りたい移動情報として、より正確な情報を得ることができる。
【００１１】
上記他の観点の発明によれば、前記第二関心領域は、観察対象に設定される前記第一関心領域を含んでおらず、前記他の時相の超音波画像が、前記第二関心領域が設定された第二部分の移動ベクトルによって位置補正された前記位置補正済超音波画像は、観察対象について知りたい移動情報以外の移動要素が正確に排除された画像である。この位置補正済超音波画像と前記一の時相の超音波画像との間で、前記第一関心領域が設定された第一部分の移動ベクトルが算出されるので、観察対象について知りたい移動情報以外の移動要素をより正確に排除し、知りたい移動情報として、より正確な情報を得ることができる
【図面の簡単な説明】
【００１２】
【図１】本発明の実施形態における超音波診断装置の概略構成の一例を示すブロック図である。
【図２】第一実施形態の制御部で実行される機能のブロック図である。
【図３】第一実施形態の超音波診断装置の作用の一例を示すフローチャートである。
【図４】第一関心領域及び第二関心領域が設定されたＢモード画像を示す図である。
【図５】時相Ｔ１及び時相Ｔ２におけるＢモード画像を示す概念図である。
【図６】表示部に表示された時相Ｔ１のＢモード画像と時相Ｔ２のＢモード画像とを示す図である。
【図７】観察対象指定領域の移動を説明する図である。
【図８】補正済移動ベクトルの算出を説明する図である。
【図９】補正済移動ベクトルを示すベクトル表示を示す図である。
【図１０】第一実施形態の変形例における第二関心領域の設定を説明する図である。
【図１１】第二実施形態の制御部で実行される機能のブロック図である。
【図１２】第二実施形態の超音波診断装置の作用の一例を示すフローチャートである。
【図１３】位置補正済Ｂモード画像の作成を説明する図である。
【図１４】二つの第二関心領域が設定されたＢモード画像を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明の実施形態について詳細に説明する。
（第一実施形態）
先ず、第一実施形態について図１〜図８に基づいて説明する。図１に示す超音波診断装置１は、超音波プローブ２、送受信ビームフォーマ３、エコーデータ処理部４、表示制御部５、表示部６、操作部７、制御部８及び記憶部９を備える。前記超音波診断装置１は、本発明における計測装置の実施の形態の一例である。
【００１４】
前記超音波プローブ２は、アレイ状に配置された複数の超音波振動子（図示省略）を有して構成され、この超音波振動子によって被検体に対して超音波を送信し、そのエコー信号を受信する。
【００１５】
前記送受信ビームフォーマ３は、前記超音波プローブ２から所定の走査条件で超音波を送信するための電気信号を、前記制御部８からの制御信号に基づいて前記超音波プローブ２に供給する。また、前記送受信ビームフォーマ３は、前記超音波プローブ２で受信したエコー信号について、Ａ／Ｄ変換、整相加算処理等の信号処理を行ない、信号処理後のエコーデータを前記エコーデータ処理部４へ出力する。
【００１６】
前記エコーデータ処理部４は、前記送受信ビームフォーマ３から出力されたエコーデータに対し、超音波画像を作成するための信号処理を行なう。例えば、前記エコーデータ処理部４は、対数圧縮処理、包絡線検波処理等を含むＢモード処理を行ない、Ｂモードデータを作成する。
【００１７】
前記表示制御部５は、前記Ｂモードデータに対し、スキャンコンバータ（Ｓｃａｎ  Ｃｏｎｖｅｒｔｅｒ）による走査変換を行なってＢモード画像データを作成し、このＢモード画像データに基づくＢモード画像を前記表示部６に表示させる。また、前記表示制御部５は、後述のベクトル表示ｖｉ等その他の表示を前記表示部６に表示させる。
【００１８】
前記表示部６は、ＬＣＤ（Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）やＣＲＴ（Ｃａｔｈｏｄｅ  Ｒａｙ  Ｔｕｂｅ）などで構成される。前記操作部７は、操作者が指示や情報を入力するためのキーボード及びポインティングデバイス（図示省略）などを含んで構成されている。
【００１９】
前記制御部８は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）であり、前記記憶部９に記憶された制御プログラムを読み出し、前記超音波診断装置１の各部における機能を実行させる。例えば、前記送受信ビームフォーマ３、前記エコーデータ処理部４、前記表示制御部５の機能は、前記制御プログラムによって実行されてもよい。
【００２０】
さらに、前記制御部８は、図２に示された第一関心領域設定部８１、第二関心領域設定部８２、第一追跡部８３、第二追跡部８４、第一移動ベクトル算出部８５、第二移動ベクトル算出部８６、補正済移動ベクトル算出部８７の機能を実行させる。詳細は後述する。前記第一関心領域設定部８１は、本発明における第一関心領域設定部の実施の形態の一例であり、その機能は本発明における第一関心領域設定機能の実施の形態の一例である。前記第二関心領域設定部８２は、本発明における第二関心領域設定部の実施の形態の一例であり、その機能は本発明における第二関心領域設定機能の実施の形態の一例である。前記第一追跡部８３は、本発明における第一追跡部の実施の形態の一例である。前記第二追跡部８４は、本発明における第二追跡部の実施の形態の一例である。前記第一移動ベクトル算出部８５は、本発明における第一移動ベクトル算出部の実施の形態の一例であり、その機能は本発明における第一移動ベクトル算出機能の実施の形態の一例である。前記第二移動ベクトル算出部８６は、本発明における第二移動ベクトル算出部の実施の形態の一例であり、その機能は本発明における第二移動ベクトル算出機能の実施の形態の一例である。前記補正済移動ベクトル算出部８７は、本発明における補正済移動ベクトル算出部の実施の形態の一例であり、その機能は本発明における補正済移動ベクトル算出機能の実施の形態の一例である。
【００２１】
前記記憶部９は、例えばＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）や半導体メモリ（ｍｅｍｏｒｙ）などである。
【００２２】
さて、本例の超音波診断装置１の作用について図３のフローチャートに基づいて説明する。先ず、ステップＳ１では、被検体に対する超音波の送受信によって得られたエコー信号に基づいて前記表示部６に表示されたＢモード画像ＢＩに、操作者は、図４に示すように第一関心領域ＲＯＩ１と第二関心領域ＲＯＩ２とを設定する。前記第一関心領域ＲＯＩ１及び前記第二関心領域ＲＯＩ２は、フリーズ（ｆｒｅｅｚｅ）させたＢモード画像ＢＩにおいて設定される。ちなみに、このＢモード画像ＢＩは、血管ＢＬの長軸断面の画像である。
【００２３】
前記第一関心領域ＲＯＩ１の設定について説明する。前記第一関心領域ＲＯＩ１は、観察対象指定領域Ｒを分割した分割領域である。先ず、操作者が、前記操作部７を用いてＢモード画像ＢＩ上で観察対象であるプラークＸを含む領域を指示する入力を行なうと、前記第一関心領域設定部８１は、指示された領域に観察対象指定領域Ｒを設定する。次に、操作者が前記操作部７を用いて前記観察対象指定領域Ｒを分割する入力を行なうと、前記第一関心領域設定部８１は前記観察対象指定領域Ｒを分割して、前記第一関心領域ＲＯＩを設定する。前記観察対象指定領域Ｒの分割数は、操作者によって指定することができてもよい。この場合、操作者は縦及び横の分割数を入力することにより、分割数を指定することができてもよい。
【００２４】
前記第二関心領域ＲＯＩ２は、操作者が、前記操作部７を用いてＢモード画像ＢＩ上で所定の領域を指示する入力を行なうと、前記第二関心領域設定部８２により、前記所定の領域に設定される。例えば、前記第二関心領域ＲＯＩ２は、血管ＢＬの下方に設定される。
【００２５】
前記第一関心領域ＲＯＩ１は、プラークＸの内部の変位を観察するために設定される領域である。また、前記第二関心領域ＲＯＩ２は、拍動による移動の要素（プラークＸの内部について操作者が知りたい移動情報）以外の排除したい移動要素を検出するための領域に設定される。ここでは、血管ＢＬの下方に設定された前記第二関心領域ＲＯＩ２により、粘弾性体である血液が前記血管ＢＬ内を移動することに伴う慣性力などによる血管ＢＬの並進運動を検出することができる。
【００２６】
次に、ステップＳ２では、図５に示すように、時相Ｔ１のＢモード画像ＢＩ１と、時相Ｔ１よりも後の時相Ｔ２のＢモード画像ＢＩ２との間で、前記第一追跡部８３は、前記第一関心領域ＲＯＩ１が設定された第一部分Ｐ１の追跡を行なう。また、前記第二追跡部８４は、同様に前記Ｂモード画像ＢＩ１と前記Ｂモード画像ＢＩ２との間で、前記第二関心領域ＲＯＩ２が設定された第二部分Ｐ２の追跡を行なう。前記各第一部分Ｐ１及び前記第二部分Ｐ２の追跡は、例えばオプティカルフロー法などの公知の手法を用いて行われる。前記Ｂモード画像ＢＩ１は、本発明における一の時相の超音波画像の実施の形態の一例である。また、前記Ｂモード画像ＢＩ２は、本発明における他の時相の超音波画像の実施の形態の一例である。
【００２７】
また、ステップＳ２では、前記各第一部分Ｐ１の追跡結果に基づいて、前記第一移動ベクトル算出部８５は、前記Ｂモード画像ＢＩ１及び前記Ｂモード画像ＢＩ２の間における前記各第一部分Ｐ１の第一移動ベクトルｐ１を算出する。また、前記第二部分Ｐ２の追跡結果に基づいて、前記第二移動ベクトル算出部８６は、前記Ｂモード画像ＢＩ１及び前記Ｂモード画像ＢＩ２の間における前記第二部分Ｐ２の第二移動ベクトルｐ２を算出する。
【００２８】
例えば、前記時相Ｔ１における超音波画像Ｂ１が、図６の上方に示された状態の画像であり、前記時相Ｔ２における超音波画像Ｂ２が、図６の下方に示された状態の画像であったとする。前記超音波画像Ｂ１における血管ＢＬに対し、前記超音波画像Ｂ２における血管ＢＬは拍動によって縮径し、なおかつ図において左方向へ並進運動した状態になっている。従って、観察対象指定領域Ｒは、図７に示すように斜め左上方向へ移動している。
【００２９】
次に、ステップＳ３では、前記補正済移動ベクトル算出部８７が、前記各第一移動ベクトルｐ１に対し、前記第二移動ベクトルｐ２による補正を行なって、前記各第一関心領域ＲＯＩ１について、補正済移動ベクトルｐｃを算出する。具体的には、図８に示すように、下記（式１）のベクトル演算（減算）を行ない、補正済移動ベクトルｐｃを算出する。
（補正済移動ベクトルｐｃ）
＝（前記第一移動ベクトルｐ１）−（前記第二移動ベクトルｐ２）
・・・（式１）
【００３０】
図８では、前記第二移動ベクトルｐ２は、水平方向のベクトルで図示されているが、前記第二部分Ｐ２が、拍動によって血管の短軸方向（縦方向）にも移動する場合、前記第二移動ベクトルｐ２は、特に図示しないがこの短軸方向の要素が加わったベクトルになる。従って、拍動による短軸方向の移動の要素も削除された補正済移動ベクトルｐｃを得ることができる。
【００３１】
ステップＳ３において補正済移動ベクトルｐｃが算出されると、ステップＳ４では、前記表示制御部５は前記表示部６に、前記第一関心領域ＲＯＩ１の移動情報ｉを表示させる。この移動情報ｉとして、図９に示すように、前記各第一関心領域ＲＯＩ１に、補正済移動ベクトルｐｃを示すベクトル表示ｖｉが表示されてもよい。このベクトル表示ｖｉは、前記補正済移動ベクトルｐｃの大きさ及び方向に応じた表示形態を有する。前記ベクトル表示ｖｉは、前記Ｂモード画像ＢＩに設定された第一関心領域ＲＯＩとは別に表示されてもよい。前記ベクトル表示ｖｉは、本発明における移動情報の実施の形態の一例である。
【００３２】
また、移動情報ｉとして、特に図示しないが、前記各第一関心領域ＲＯＩ１に、前記補正済移動ベクトルｐｃに応じた色（ベクトルの方向と大きさに応じた色）を表示してもよい。
【００３３】
以上説明した本例によれば、第一移動ベクトルｐ１から第二移動ベクトルｐ２が引かれて補正済移動ベクトルｐｃが算出され、この補正済移動ベクトルｐｃに基づく移動情報ｉが表示されるので、並進運動の要素が除かれた第一関心領域ＲＯＩの移動情報を表示することができる。しかも、前記第二移動ベクトルｐ２は、前記第一関心領域ＲＯＩ１を含まない部分（前記第二関心領域ＲＯＩ２が設定された第二部分Ｐ２）の移動ベクトルであるので、観察対象であるプラーク内部について知りたい移動情報以外の移動要素が正確に排除された移動情報を得ることができる。
【００３４】
次に、第一実施形態の変形例について説明する。前記第二関心領域設定部８２は、第二関心領域ＲＯＩ２′を所定の位置に自動的に設定してもよい。例えば、図１０に示すように、前記第二関心領域ＲＯＩ２′は被検体の体表面の近傍に設定されてもよい。このように被検体の体表面の近傍に設定された前記第二関心領域ＲＯＩ２′により、前記超音波プローブ２を体表面上において煽る動作などによる生体組織の移動を検出することができる。
【００３５】
前記ステップＳ２では、前記時相Ｔ１のＢモード画像ＢＩ１及び前記時相Ｔ２のＢモード画像ＢＩ２との間で前記第二関心領域ＲＯＩ２′が設定された第二部分Ｐ２′の追跡が行われ第二移動ベクトルｐ２が算出される。
【００３６】
なお、上記実施形態において、操作者が前記操作部７を用いて前記第二関心領域ＲＯＩ２を設定する場合においても、第一変形例と同じ位置、すなわち被検体の体表面の近傍に設定されてもよい。
【００３７】
（第二実施形態）
次に、第二実施形態について説明する。第二実施形態において、前記制御部８は、図１１に示された第一関心領域設定部８１、第二関心領域設定部８２、第一追跡部８３、第二追跡部８４、第一移動ベクトル算出部８５、第二移動ベクトル算出部８６、画像移動部８８の機能を実行させる。本例の制御部８は、第一実施形態とは異なり、画像移動部８８を有する一方で、補正済移動ベクトル算出部８７を有さない。
【００３８】
次に、第二実施形態の超音波診断装置１の作用について図１２のフローチャートに基づいて説明する。先ず、ステップＳ１１では、第一実施形態のステップＳ１と同様に、前記第一関心領域ＲＯＩ１及び前記第二関心領域ＲＯＩ２を設定する。
【００３９】
次に、ステップＳ１２では、時相Ｔ１のＢモード画像ＢＩ１と、時相Ｔ１よりも後の時相Ｔ２のＢモード画像ＢＩ２との間で、前記第二追跡部８４は、前記第二関心領域が設定された第二部分Ｐ２の追跡を行なう。また、前記第二移動ベクトル算出部８６が、前記第二部分Ｐ２の追跡結果に基づいて、前記Ｂモード画像ＢＩ１及び前記Ｂモード画像ＢＩ２の間における前記第二部分Ｐ２の第二移動ベクトルｐ２を算出する。
【００４０】
次に、ステップＳ１３では、前記画像移動部８８は、図１３に示すように、前記第二移動ベクトルｐ２の逆ベクトル（−ｐ２）によって前記Ｂモード画像ＢＩ２を位置補正した位置補正済Ｂモード画像ＢＩ２′のデータを作成する。これにより、前記Ｂモード画像ＢＩ２に対し、前記第二移動ベクトルｐ２と同じ大きさで逆方向に移動した位置補正済Ｂモード画像ＢＩ２′のデータが得られる。これにより、並進運動前の位置に戻った位置補正済Ｂモード画像ＢＩ２′のデータを得ることができる。この位置補正済Ｂモード画像ＢＩ２′は、拍動によるプラーク内部の移動以外の移動要素が排除された画像である。
【００４１】
次に、ステップＳ１４では、前記第一追跡部８３は、時相Ｔ１のＢモード画像ＢＩ１と前記位置補正済Ｂモード画像ＢＩ２′との間で、前記第一関心領域ＲＯＩ１が設定された第一部分Ｐ１の追跡を行なう。また、前記第一部分Ｐ１の追跡結果に基づいて、前記Ｂモード画像ＢＩ１及び前記位置補正済Ｂモード画像ＢＩ２′の間における前記第一部分Ｐ１の第一移動ベクトルｐ１を算出する。
【００４２】
ステップＳ１４において前記第一移動ベクトルｐ１が算出されると、ステップＳ１５では、前記表示制御部５は前記表示部６に、前記第一関心領域ＲＯＩ１の移動情報ｉを表示させる。本例においても、第一実施形態と同様に、前記移動情報ｉとしてベクトル表示ｖｉ（図８参照）が表示されるが、このベクトル表示ｖｉは、前記ステップＳ１４において算出された前記第一移動ベクトルｐ１の大きさ及び方向を示す表示である。
【００４３】
また、移動情報ｉとして、特に図示しないが前記各第一関心領域ＲＯＩ１に、前記ステップＳ１４において算出された前記第一移動ベクトルｐ１に応じた色（ベクトルの方向と大きさに応じた色）を表示してもよい。
【００４４】
ちなみに、時相Ｔ２とこの時相Ｔ２の後の時相Ｔ３との間における第一移動ベクトルｐ１を算出する場合には、前記ステップＳ１２では、前記第二移動ベクトル算出部８６は、前記位置補正済Ｂモード画像ＢＩ２′と前記時相Ｔ３のＢモード画像ＢＩ３との間で、前記第二部分Ｐ２の追跡を行なって第二移動ベクトルｐ２を算出する。そして、この第二移動ベクトルｐ２によって前記Ｂモード画像ＢＩ３を位置補正した位置補正済Ｂモード画像ＢＩ３′を得て、この位置補正済Ｂモード画像ＢＩ３′と前記位置補正済Ｂモード画像ＢＩ２′との間で、前記第一部分Ｐ１の追跡を行ない、第一移動ベクトルｐ１を算出する。
【００４５】
以上説明した本例によれば、並進運動がキャンセルされた位置補正済Ｂモード画像ＢＩ２′と第一時相のＢモード画像ＢＩ２との間で、前記第一関心領域ＲＯＩ１が設定された第一部分Ｐ１の第一移動ベクトルｐ１が算出される。そして、この第一移動ベクトルｐ１に基づく移動情報ｉが表示されるので、並進運動の要素が除かれた第一関心領域ＲＯＩの移動情報を表示することができる。しかも、第一実施形態と同様に、前記第二移動ベクトルｐ２は、前記第一関心領域ＲＯＩ１を含まない部分（前記第二関心領域ＲＯＩ２が設定された第二部分Ｐ２）の移動ベクトルであるので、観察対象であるプラーク内部について知りたい移動情報以外の移動要素が正確に排除された移動情報を得ることができる。
【００４６】
なお、第二実施形態においても、第一実施形態の変形例と同様に、前記第二関心領域ＲＯＩ２′が所定の位置に自動的に設定されてもよい。
【００４７】
以上、本発明を前記各実施形態によって説明したが、本発明はその主旨を変更しない範囲で種々変更実施可能なことはもちろんである。例えば、上記第二実施形態において、前記第二関心領域ＲＯＩ２が複数設定されてもよい。この場合、複数の第二部分Ｐ２の移動ベクトルの逆ベクトルによって、時相Ｔ２のＢモード画像を位置補正した位置補正済超音波画像が作成される。
【００４８】
具体的には、図１４に示すように、Ｂモード画像ＢＩにおいて、血管ＢＬの下方及び被検体の体表面の近傍に、二つの前記第二関心領域ＲＯＩ２−１，ＲＯＩ２−２が設定されてもよい。この場合、前記第二移動ベクトル算出部８６は、時相Ｔ１のＢモード画像ＢＩ１と時相Ｔ２のＢモード画像ＢＩ２との間で、各第二関心領域ＲＯＩ２−１，ＲＯＩ２−２が設定された第二部分Ｐ２−１，Ｐ２−２の第二移動ベクトルｐ２−１，ｐ２−２をそれぞれ算出する。そして、前記第二移動ベクトルｐ２−１の逆ベクトルによって前記Ｂモード画像ＢＩ２が位置補正された位置補正済Ｂモード画像ＢＩ２′を作成し、さらにこの位置補正済Ｂモード画像ＢＩ２′が前記第二移動ベクトルｐ２−２の逆ベクトルによって位置補正された位置補正済Ｂモード画像ＢＩ２′′を作成する。
【００４９】
また、前記第二関心領域ＲＯＩ２−１，ＲＯＩ２−２が設定された場合に、時相Ｔ１のＢモード画像ＢＩ１と時相Ｔ２のＢモード画像ＢＩ２との間で、先ず、前記第二移動ベクトル算出部８６は、前記第二関心領域ＲＯＩ２−２が設定された第二部分Ｐ２−２の第二移動ベクトルｐ２−２を算出し、この第二移動ベクトルｐ２−２の逆ベクトルによって前記Ｂモード画像ＢＩ２が位置補正された位置補正済Ｂモード画像ＢＩ２′を作成してもよい。この場合、次に、前記第二移動ベクトル算出部８６は、前記位置補正済Ｂモード画像ＢＩ２′と前記Ｂモード画像ＢＩ１との間で、前記第二関心領域ＲＯＩ２−１が設定された第二部分Ｐ２−１の第二移動ベクトルｐ２−１を算出し、この第二移動ベクトルｐ２−１の逆ベクトルによって前記位置補正済Ｂモード画像ＢＩ２′が位置補正された位置補正済Ｂモード画像ＢＩ２′′を作成する。この場合、位置補正済Ｂモード画像ＢＩ２′は、他の時相の超音波画像の実施の形態の一例である。
【００５０】
位置補正済Ｂモード画像ＢＩ２′′が作成されると、この位置補正済Ｂモード画像ＢＩ２′′と前記時相Ｔ１のＢモード画像ＢＩ１との間で前記第一部分Ｐ１の追跡を行ない、前記第一移動ベクトルｐ１を算出する。
【００５１】
また、本発明に係る計測装置は、超音波診断装置以外の機器において実施されてもよい。例えば、本発明に係る計測装置は、パーソナルコンピュータ（ｐｅｒｓｏｎａｌ  ｃｏｍｐｕｔｅｒ）などの汎用のコンピュータにおいて実施されてもよい。この場合、超音波診断装置から、Ｂモードデータなどのローデータ又はＢモード画像データなどの画像データを、例えば汎用のコンピュータに取り込んで、この汎用のコンピュータにおいて、上記実施形態で説明した処理を行なう。
【符号の説明】
【００５２】
１  超音波診断装置（計測装置）
６  表示部                        
８１  第一関心領域設定部
８２  第二関心領域設定部
８３  第一追跡部
８４  第二追跡部
８５  第一移動ベクトル算出部
８６  第二移動ベクトル算出部
８７  補正済移動ベクトル算出部
【図１】
2014087507.tif 000003
【図２】
2014087507.tif 000004
【図３】
2014087507.tif 000005
【図４】
2014087507.tif 000006
【図５】
2014087507.tif 000007
【図６】
2014087507.tif 000008
【図７】
2014087507.tif 000009
【図８】
2014087507.tif 000010
【図９】
2014087507.tif 000011
【図１０】
2014087507.tif 000012
【図１１】
2014087507.tif 000013
【図１２】
2014087507.tif 000014
【図１３】
2014087507.tif 000015
【図１４】
2014087507.tif 000016