(57)【要約】【課題】短時間で正確に、しかも被検体に大きな負担を課さずに直接的な動脈硬化の診断を行なうことが可能な動脈硬化検出装置を提供する。【解決手段】血管径の数値シミュレーションを行なう数値シミュレーション部８１２と、被検体に送信された超音波のエコーデータに基づいて血管径を計測する径計測部８２と、前記数値シミュレーション部８１２で得られた血管径と前記径計測部８２で得られた血管径との差を算出する差分算出部８３と、この差分算出部８３で算出された差に基づいて動脈硬化に関する評価を行なう評価部８４と、この評価部８４による評価結果を表わす画像が表示される表示部と、を備えることを特徴とする。【選択図】図２
2012055621.tif 000002
【特許請求の範囲】
【請求項１】
血管径の数値シミュレーションを行なう数値シミュレーション部と、
被検体に送信された超音波のエコーデータに基づいて血管径を計測する径計測部と、
前記数値シミュレーション部で得られた血管径と前記径計測部で得られた血管径との差を算出する差分算出部と、
該差分算出部で算出された差に基づいて動脈硬化に関する評価を行なう評価部と、
該評価部による評価結果を表わす画像が表示される表示部と、
を備えることを特徴とする動脈硬化検出装置。
【請求項２】
前記数値シミュレーション部は、血管壁の壁面モデルとしてフォークトモデルを適用して数値シミュレーションを行なうことを特徴とする請求項１に記載の動脈硬化検出装置。
【請求項３】
前記数値シミュレーション部による数値シミュレーションは、血管壁に血流の圧力がかかった場合の前記フォークトモデルにおける血管壁の変位を算出し、該変位を初期の血管径に加算して、血管径の算出を行なうものであることを特徴とする請求項２に記載の動脈硬化検出装置。
【請求項４】
前記フォークトモデルにおいて並列に接続されたバネとダンパーにおける弾性係数及び粘性係数は、被検体の血管壁について、年代別に実験的に得られた値であることを特徴とする請求項２又は３に記載の動脈硬化検出装置。
【請求項５】
前記数値シミュレーション部は、前記エコーデータに基づく超音波画像における血管に設定された計算格子において数値シミュレーションを行なうことを特徴とする請求項１〜４のいずれか一項に記載の動脈硬化検出装置。
【請求項６】
前記超音波画像において特定される血管に前記計算格子を設定する計算格子設定部を備えることを特徴とする請求項５に記載の動脈硬化検出装置。
【請求項７】
前記径計測部は、超音波画像において特定された血管の径を計測することを特徴とする請求項１〜６のいずれか一項に記載の動脈硬化検出装置。
【請求項８】
前記評価結果を表わす画像は、前記エコーデータに基づく超音波画像上に表示されることを特徴とする請求項１〜７のいずれか一項に記載の動脈硬化検出装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、動脈硬化を検出する動脈硬化検出装置に関する。
【背景技術】
【０００２】
従来、動脈硬化を診断する手法として血液検査による手法があるが（例えば、特許文献１参照）、この手法は間接的な手法に過ぎないため、直接的に動脈硬化を診断する手法が試みられている。この直接的な動脈硬化の診断手法としては、例えば脈波伝播速度ＰＷＶ（Ｐｕｌｓｅ  Ｗａｖｅ  Ｖｅｌｏｃｉｔｙ）を計測して診断を行なう手法や、一心拍の血圧変動に対応する血管径の変化に基づいて算出されるＳｔｉｆｆｎｅｓｓ  Ｐａｒａｍｅｔｅｒ  βに基づいて診断を行なう手法や、被検体の上腕動脈においてカフなどによって血流を一時遮断し、開放後の血管径増加率を％ＦＭＤ（Ｆｌｏｗ−Ｍｅｄｉａｔｅｄ  Ｄｉｌａｔｉｏｎ）として計測し診断を行なう手法などが知られている（例えば、特許文献２参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１０−１６９６９６号公報
【特許文献２】特開２００４−１０５５４９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかし、ＰＷＭを計測して診断する手法では、伝播速度が互いに異なる弾性血管と筋性血管とを区別せずに計測を行なっているため正確な診断を行なうことが困難であるなどの問題がある。また、Ｓｔｉｆｆｎｅｓｓ  Ｐａｒａｍｅｔｅｒ  βに基づいて診断を行なう手法では、計測に要する時間が２０〜３０分と長いことなどが問題である。また、％ＦＭＤを計測する手法においても、検査に時間がかかるという問題や、検査者間、検査施設間において再現性が低いといった問題がある。さらに、上述の三つのいずれの手法も被検体の動脈を遮断する必要があり、被検体に対する負担が重い。
【０００５】
このようなことから、短時間で正確に、しかも被検体に大きな負担を課さずに直接的な動脈硬化の診断を行なうことが可能な装置が望まれている。
【課題を解決するための手段】
【０００６】
この発明は、前記課題を解決するためになされたもので、第１の観点の発明は、  血管径の数値シミュレーションを行なう数値シミュレーション部と、被検体に送信された超音波のエコーデータに基づいて血管径を計測する径計測部と、前記数値シミュレーション部で得られた血管径と前記径計測部で得られた血管径との差を算出する差分算出部と、この差分算出部で算出された差に基づいて動脈硬化に関する評価を行なう評価部と、この評価部による評価結果を表わす画像が表示される表示部と、を備えることを特徴とする動脈硬化検出装置である。
【０００７】
第２の観点の発明は、第１の観点の発明において、前記数値シミュレーション部は、血管壁の壁面モデルとしてフォークトモデルを適用して数値シミュレーションを行なうことを特徴とする動脈硬化検出装置である。
【０００８】
第３の観点の発明は、第２の観点の発明において、前記数値シミュレーション部による数値シミュレーションは、血管壁に血流の圧力がかかった場合の前記フォークトモデルにおける血管壁の変位を算出し、該変位を初期の血管径に加算して、血管径の算出を行なうものであることを特徴とする動脈硬化検出装置である。
【０００９】
第４の観点の発明は、第２、３の観点の発明において、前記フォークトモデルにおいて並列に接続されたバネとダンパーにおける弾性係数及び粘性係数は、被検体の血管壁について、年代別に実験的に得られた値であることを特徴とする動脈硬化検出装置である。
【００１０】
第５の観点の発明は、第１〜４のいずれか一の観点の発明において、前記数値シミュレーション部は、前記エコーデータに基づく超音波画像における血管に設定された計算格子において数値シミュレーションを行なうことを特徴とする動脈硬化検出装置である。
【００１１】
第６の観点の発明は、第５の観点の発明において、前記超音波画像において特定される血管に前記計算格子を設定する計算格子設定部を備えることを特徴とする動脈硬化検出装置である。
【００１２】
第７の観点の発明は、第１〜６のいずれか一の観点の発明において、前記径計測部は、超音波画像において特定された血管の径を計測することを特徴とする動脈硬化検出装置である。
【００１３】
第８の観点の発明は、第１〜７のいずれか一の観点の発明において、前記評価結果を表わす画像は、前記エコーデータに基づく超音波画像上に表示されることを特徴とする動脈硬化検出装置である。
【発明の効果】
【００１４】
上記観点の発明によれば、前記数値シミュレーション部で得られた血管径と前記径計測部で得られた血管径との差に基づく動脈硬化に関する評価結果を表わす画像が表示されるので、医師等の診断者は、短時間で正確に、しかも被検体に大きな負担を課さずに直接的な動脈硬化の診断を行なうことができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の実施形態における超音波診断装置の概略構成の一例を示すブロック図である。
【図２】エコーデータ処理部と制御部の詳細構成を示すブロック図である。
【図３】動脈硬化に関する評価を表わす画像の作成及び表示と、血流の速度及び圧力の算出のフローチャートである。
【図４】１心周期における血圧の変動を示す図である。
【図５】計算格子の設定の説明図である。
【図６】フォークトモデルを示す図である。
【図７】数値シミュレーション部で得られた血管径と径計測部で得られた血管径との差が閾値を超えた部分にカラー画像が表示された表示部の一例を示す図である。
【図８】変形例におけるエコー処理部と制御部の詳細構成を示すブロック図である。
【図９】変形例における作用を示すフローチャートである。
【発明を実施するための形態】
【００１６】
以下、本発明の実施形態について図１〜図６に基づいて詳細に説明する。図１に示す超音波診断装置１は、超音波プローブ２、送受信部３、エコーデータ処理部４、表示制御部５、表示部６、操作部７及び制御部８を備える。前記超音波診断装置は、本発明における動脈硬化検出装置の実施の形態の一例である。
【００１７】
前記超音波プローブ２は、生体組織に対して超音波の送受信を行なう。また、前記送受信部３は、前記超音波プローブ２を所定のスキャンパラメータで駆動させてスキャン面を走査させる。そして、前記送受信部３は、前記超音波プローブ２で得られたエコー信号について、整相加算処理等の信号処理を行なう。
【００１８】
前記エコーデータ処理部４は、図２に示すようにＢモード処理部４１及びドプラ処理部４２を有する。前記Ｂモード処理部４１は、前記送受信部３から出力されたエコー信号に対し、対数圧縮処理、包絡線検波処理等のＢモード処理を行ない、Ｂモードデータを作成する。また、前記ドプラ処理部４２は、前記送受信部３から出力されたエコー信号に基づいてドプラデータを作成する。ドプラデータには、流速データ、分散データ及びパワーデータが含まれる。
【００１９】
具体的には、前記ドプラ処理部４２は、前記送受信部３から出力されたエコー信号に対し、直交検波処理を行ない、ＭＴＩフィルタ（Ｍｏｖｉｎｇ  Ｔａｒｇｅｔ  Ｉｎｄｉｃａｔｉｏｎ  Ｆｉｌｔｅｒ）でＭＴＩ処理してエコーのドプラ信号を求め、さらにＭＴＩ処理後の信号に対し、自己相関演算を行なう。そして、自己相関演算結果から平均流速と、流速の分散と、パワーとを求める。
【００２０】
前記表示制御部５は、前記Ｂモードデータに対し、スキャンコンバータ（Ｓｃａｎ  Ｃｏｎｖｅｒｔｅｒ）による走査変換を行なってＢモード画像データを作成し、このＢモード画像データに基づくＢモード画像を前記表示部６に表示させる。ただし、前記表示制御部５は、Ｂモード画像データの作成に加え、前記ドプラデータに対してスキャンコンバータによる走査変換を行なってドプラ画像データを作成し、このドプラ画像データに基づくドプラ画像を表示させるようになっていてもよい。
【００２１】
また、前記表示制御部５は、後述するように動脈硬化の指標を表わす画像を前記表示部６に表示させるようになっている。
【００２２】
前記表示部６は、ＬＣＤ（Ｌｉｑｕｉｄ  Ｃｒｙｓｔａｌ  Ｄｉｓｐｌａｙ）やＣＲＴ（Ｃａｔｈｏｄｅ  Ｒａｙ  Ｔｕｂｅ）などで構成される。前記操作部７は、操作者が指示や情報を入力するためのキーボード及びポインティングデバイス（図示省略）などを含んで構成されている。
【００２３】
前記制御部８は、特に図示しないがＣＰＵ（ＣｅｎｔＲａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を有して構成される。この制御部８は、前記ＨＤＤ９に記憶された制御プログラムを読み出し、前記超音波診断装置１の各部における機能を実行させる。
【００２４】
前記制御部８の詳細構成について説明すると、この制御部８は、計測融合シミュレーション部８１、径計測部８２、差分算出部８３及び評価部８４を有している。そして、前記計測融合シミュレーション部８１は、計算格子設定部８１１、数値シミュレーション部８１２及びフィードバック部８１３を有している。これら各部の説明は後述する。ちなみに、前記径計測部８２は本発明における径計測部の実施の形態の一例であり、前記差分算出部８３は、本発明における差分算出部の実施の形態の一例である。また、前記評価部８４は本発明における評価部の実施の形態の一例であり、前記数値シミュレーション部８１２は、本発明における数値シミュレーション部の実施の形態の一例である。
【００２５】
さて、本例の超音波診断装置１の作用について説明する。本例の超音波診断装置１においては、超音波の送受信を行なってエコーデータを取得し、前記Ｂモードデータ及び前記ドプラデータを作成する。これらＢモードデータ及びドプラデータは、図示しないＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）やＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）等のメモリ、ＨＤＤ（ハードディスクドライブ：Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）などに記憶される。そして、これらＢモードデータ及びドプラデータを用いて、動脈硬化に関する評価を表わす画像を作成して表示し、さらには血流の速度や圧力を算出する。
【００２６】
動脈硬化に関する評価を表わす画像の作成及び表示と、血流の速度及び圧力の算出について、図３に示すフローチャートに基づいて説明する。この図３に示すフローチャートでは、１心周期につき、Ｓ１〜Ｓ４の各ステップをＮループ繰り返す。図４は、１心周期における血圧の変動を表わすグラフであり、前記Ｓ１〜Ｓ４の各ステップをＮループ繰り返すことにより、１心周期における各時点の血圧の血管径、血流の速度及び圧力が算出され、１心周期で拡張及び収縮する血管における血管径、血流の速度及び圧力が、１心周期についてＮ回算出される。
【００２７】
各ステップについて具体的に説明すると、先ずステップＳ１では、前記数値シミュレーション部８１２が血管径の数値シミュレーションを行なって血管径Ｘｓを算出し、前記径計測部８３がＢモード画像における血管径Ｘｍの計測を行なう。ここでは、前記血管径Ｘｓ，Ｘｍとして、血管の直径を算出及び計測するものとする。
【００２８】
前記数値シミュレーション部８１２による数値シミュレーションについて具体的に説明する。前記数値シミュレーション部８１２による数値シミュレーションを行なうにあたり、先ず前記計算格子設定部８１１は、図５に示すようにＢモード画像ＢＧにおいて特定される血管ｂｌの上半分における所定の領域Ｒに計算格子Ｃｉ（ｉ＝１〜ｎ（ｎ：自然数））を設定する。血管ｂｌの特定はＢモードデータに基づいて行なわれる。ちなみに、Ｂモード画像ＢＧにおいて血管ｂｌは周囲よりも輝度が低いため、Ｂモードデータにおいて血管ｂｌの特定が可能である。なお、Ｂモード画像データに基づいて血管ｂｌの特定が行なわれてもよい。
【００２９】
前記数値シミュレーション部８１２は、血流の圧力によって血管壁の変位ｘ（ｉ）が生じ、初期半径ｒｓ０（ｉ）の状態から半径ｒｓ（ｉ）となった場合のこの半径ｒｓ（ｉ）を算出する。すなわち、下記の（式１）により、半径ｒｓ（ｉ）を算出する。
ｒｓ（ｉ）＝ｒｓ０（ｉ）＋ｘ（ｉ）  ・・・（式１）
半径ｒｓ（ｉ）の算出は、各計算格子Ｃｉ毎に行なう。ただし、本例では血管の上半分又は下半分における各計算格子Ｃｉについて半径ｒｓ（ｉ）の算出を行なうものとする。従って、前記計算格子Ｃｉは上半分又は下半分のみ設定されてもよい。そして、前記数値シミュレーション部８１２は、算出された半径ｒｓ（ｉ）を２倍して直径を算出し、血管径Ｘｓ（ｉ）を得る。これにより、ｎ個の計算格子Ｃｉについての血管径Ｘｓ（ｉ）が得られる。
【００３０】
前記（式１）において、初期半径ｒｓ０（ｉ）は、計測部位における血管部分の標準的な半径（図４の時刻０における半径）とする。
【００３１】
次に、前記（式１）における変位ｘ（ｉ）の算出について説明する。血管壁ｗの壁面モデルとして、図６に示すようにバネｓｐとダンパーｄａが並列に接続されたモデルであるフォークトモデルを適用した場合、壁面における力学の釣り合いは、下記（式２）となる。
【数１】
2012055621.tif 000003
（式２）において、ｆ（ｉ）は血管壁ｗにかかる力、ｋは血管壁ｗの弾性係数、ｃは血管壁ｗの粘性係数である。変位ｘ（ｉ）は、（式２）のｆ（ｉ），ｋ，ｃに具体的な数値を代入して算出される。
【００３２】
前記（式２）において、ｋ，ｃは、被検体の年代に応じて実験的に得られた値であって、年代別に得られた健常者の平均値である。また、（式２）において、ｆ（ｉ）は下記（式３）で得られる値である。
ｆ（ｉ）＝Ｐ（ｉ）×Ｓ（ｉ）  ・・・（式３）
この（式３）において、Ｐ（ｉ）は所定点（例えば後述のステップＳ３で圧力を算出する際に設定される計算格子）における血流の圧力（血圧）、Ｓ（ｉ）は計算格子Ｃｉにおける血管壁ｗの面積であり、ｆ（ｉ）として計算格子Ｃｉにおける血管壁ｗ全体にかかる圧力が算出される。前記Ｐ（ｉ）としては、前のループにおけるステップＳ３で算出された値を近似的な値として用いる。なお、最初のループにおいては予め設定された値を用いる。
【００３３】
前記径計測部８３による血管径Ｘｍ（ｉ）の計測について具体的に説明する。例えば、Ｂモード画像において、血管は周囲の組織よりも輝度が低くなっている。従って、前記径計測部８３は、Ｂモード画像の輝度に基づいて特定される血管の直径を算出する。血管径は、計算格子ごとに算出される。
【００３４】
次に、ステップＳ２では、前記差分算出部８１３が、前記数値シミュレーション部８１２で得られた血管径Ｘｓ（ｉ）と前記径計測部８３で得られた血管径Ｘｍ（ｉ）の差Ｄ（ｉ）を求める。すなわち、
Ｄ（ｉ）＝Ｘｓ（ｉ）−Ｘｍ（ｉ）  ・・・（式４）
そして、この（式４）で算出された差Ｄ（ｉ）に基づいて、前記評価部８４が動脈硬化に関する評価を行なう。
【００３５】
ここで、動脈硬化が進行していると血管内が閉塞し、血管径が健常者と比べて小さくなる。前記数値シミュレーション部８１２で得られた血管径Ｘｓ（ｉ）は、年代別の健常者の血管径であり、前記差Ｄ（ｉ）が大きくなるほど動脈硬化が進行していることになる。従って、前記評価部８４は、動脈硬化に関する評価として、前記差Ｄ（ｉ）が所定の閾値ＤＴＨを超えているか否かを判定する。この閾値ＤＴＨは、動脈硬化と診断される値に設定されるものとする。前記評価部８４は、前記差Ｄ（ｉ）が前記閾値ＤＴＨを超えていると判定した場合、前記表示制御部５へ信号を出力する。
【００３６】
ステップＳ２において、前記差Ｄ（ｉ）が前記閾値ＤＴＨを超えていると判定された場合、図７に示すように、前記表示部６に表示されたＢモード画像ＢＧには、前記差Ｄ（ｉ）が前記閾値ＤＴＨを超えている部分に、所定の色からなるカラー画像ＣＧを表示する。このカラー画像ＣＧは、血管ｂｌの外側にこの血管ｂｌに沿うようにして帯状に表示される。前記カラー画像ＣＧは、動脈硬化に関する評価結果を表わす画像の実施の形態の一例である。
【００３７】
次に、ステップＳ３では、前記計測融合シミュレーション部８１が計測融合シミュレーションを行なって血流の速度及び圧力の算出を行なう。具体的には、Ｂモード画像における血管に設定された計算格子について、前記数値シミュレーション部８１２が、血流の速度及び圧力についての数値シミュレーションを、ナビエ・ストークス方程式と圧力方程式とを用いて行ない、この数値シミュレーションによる算出値と、前記ドプラ処理部４２で得られたドプラデータとの差を前記フィードバック部８１３で算出して前記数値シミュレーション部８１２にフィードバックし、数値シミュレーションで得られる速度との差を補償するように計算結果を実際の血流場に収束させて、血流の速度及び圧力を算出する（詳細は、例えば特許４２６９６２３号公報参照）。
【００３８】
ちなみに、ステップＳ３における計算格子は、ステップＳ１において設定される計算格子とは異なり、前記特許４２６９６２３号公報に記載されているように、格子点（縦線と横線の交差点）であってもよく、格子点における血流の速度及び圧力を算出するようになっていてもよい。
【００３９】
ちなみに、特に図示しないが、ステップＳ３において算出された血流の速度や圧力を示す画像を前記表示部６に表示してもよい。
【００４０】
次に、ステップＳ４では、Ｎループ目であるか否かの判定を行なう。Ｎループ目に達していない場合（ステップＳ４においてＮＯ）、ステップＳ１の処理へ戻り、再度このステップＳ１以降の処理を行なう。一方、Ｎループ目に達した場合（ステップＳ４においてＹＥＳ）、処理を終了する。
【００４１】
以上説明した本例の超音波診断装置１によれば、前記数値シミュレーション部８１２で得られた血管径Ｘｓ（ｉ）と前記径計測部８２で得られた血管径Ｘｍ（ｉ）との差Ｄ（ｉ）に基づいて、動脈硬化に関する評価結果を表わす画像として、前記カラー画像ＣＧが表示されるので、被検体の動脈を遮断する必要がなく大きな負担を課さずに、短時間で正確に直接的な動脈硬化の診断を行なうことができる。
【００４２】
次に、上記実施形態の変形例について説明する。この変形例においては、図８に示すように前記制御部８は差分算出部８３を有しておらず、前記フィードバック部８１３からの入力に基づいて、前記評価部８４が動脈硬化に関する評価を行なう。具体的に図９のフローチャートに基づいて説明する。
【００４３】
先ず、図９のステップＳ１′において、前記計測融合シミュレーション部８１が前記所定の領域Ｒ（図５参照）内において計測融合シミュレーションを行なって、血管径Ｘｒ（ｉ）を算出する。具体的には、先ず上述のステップＳ１と同様にして前記数値シミュレーション部８１２が血管径Ｘｓ（ｉ）を算出し、前記径計測部８２により血管径Ｘｍ（ｉ）を得る。そして、前記フィードバック部８１３がこれら血管径Ｘｓ（ｉ）と血管径Ｘｍ（ｉ）との差Ｄ（ｉ）を求めてこの差Ｄ（ｉ）に応じた仮想的外力ｆｓ（ｉ）を数値シミュレーションにフィードバックすることにより、前記差Ｄ（ｉ）を補償するように数値シミュレーションの結果を実際の血管径の値に収束させる。これにより前記血管径Ｘｒ（ｉ）が得られる。この場合において、前記フィードバック部８１３は、本発明における差分算出部の実施の形態の一例である。
【００４４】
フィードバックは、仮想的外力ｆｓ（ｉ）を上記（式２）の左辺に与えて行なう。具体的には、以下の（式５）である。
【数２】
2012055621.tif 000004
また、フィードバックに用いる仮想的外力ｆｓ（ｉ）は、下記（式６）となる。
ｆｓ（ｉ）＝Ｋ（Ｘｓ（ｉ）−Ｘｍ（ｉ））  ・・・（式６）
（式６）において、Ｋは計算の加速係数であり圧力の単位を含む係数である。
【００４５】
以上の計測融合シミュレーションにより得られた血管径Ｘｒ（ｉ）は、前記表示部６に表示されるようになっていてもよい。
【００４６】
前記ステップＳ１′において前記フィードバック部８１３で算出された前記差Ｄ（ｉ）は、前記評価部８４へ出力される。そして、ステップＳ２′において、前記評価部８４は、前記差Ｄ（ｉ）に基づいて、上述のステップＳ２と同様にして動脈硬化の評価を行なう。そして、前記差Ｄ（ｉ）が前記閾値ＤＴＨを超えている部分に、所定の色からなるカラー画像ＣＧを表示する。
【００４７】
ただし、前記評価部８４は、前記仮想的外力ｆｓ（ｉ）に基づいて動脈硬化の評価を行なってもよい。前記差Ｄ（ｉ）が大きくなるほど、前記仮想的外力ｆｓ（ｉ）も大きくなるので、前記仮想的外力ｆｓ（ｉ）が所定の閾値を超えていれば動脈硬化が進行しているものとする。前記仮想的外力ｆｓ（ｉ）が前記閾値を超えている場合、閾値を超えていると判定された部分に、所定色からなるカラー画像ＣＧを表示する。
【００４８】
ステップＳ３，Ｓ４については、上述と同様であり説明を省略する。
【００４９】
以上、本発明を前記各実施形態によって説明したが、本発明はその主旨を変更しない範囲で種々変更実施可能なことはもちろんである。例えば、前記評価部８４による動脈硬化に関する評価として、前記差Ｄ（ｉ）の大きさに応じた色からなるカラー画像ＣＧを表示するようにしてもよい。
【符号の説明】
【００５０】
１  超音波診断装置（動脈硬化検出装置）
６  表示部                        
８２  径計測部
８３  差分算出部
８４  評価部
８１１  計算格子設定部
８１２  数値シミュレーション部
８１３  フィードバック部（差分算出部）
【図１】
2012055621.tif 000005
【図２】
2012055621.tif 000006
【図３】
2012055621.tif 000007
【図４】
2012055621.tif 000008
【図５】
2012055621.tif 000009
【図６】
2012055621.tif 000010
【図７】
2012055621.tif 000011
【図８】
2012055621.tif 000012
【図９】
2012055621.tif 000013