対象血管を含む二次元断面を超音波で繰り返し走査し、複数の超音波データを取得するデータ取得手段と、前記複数の超音波データを用いて複数の超音波画像を生成する画像生成手段と、前記各超音波画像上の少なくとも二か所において、前記対象血管の走行方向に対する中心軸と前記二次元断面とのずれを判定するための特徴量を計算する計算手段と、前記少なくとも二か所における前記特徴量に基づいて、前記複数の超音波画像の中から、最適画像を判定する判定手段と、前記最適画像を用いて、前記対象血管の中膜厚を計測する計測手段と、を具備することを特徴とする超音波診断装置。前記データ取得手段は、周期的に変化する生体信号と同期しながら前記走査することで、前記複数の超音波データを前記生体信号の時相と対応付けて取得し、前記画像生成手段は、前記生体信号の時相と対応付けられた前記複数の超音波画像を生成し、前記計算手段は、前記複数の超音波画像の中から所定時相に対応する複数の超音波画像を選択し、当該選択された各超音波画像上の少なくとも二か所における前記特徴量を計算し、前記判定手段は、前記少なくとも二か所における前記特徴量に基づいて、前記選択された複数の超音波画像の中から前記最適画像を判定すること、を具備することを特徴とする請求項1記載の超音波診断装置。前記生体信号は、ECG信号であることを特徴とする請求項2記載の超音波診断装置。前記計算手段は、前記各超音波画像の画素値又は画素値変化率に基づいて被計測領域を抽出すると共に、当該被計測領域の短手方向幅を前記特徴量として計算することを特徴とする請求項1乃至3のうちいずれか一項記載の超音波診断装置。前記計算手段は、前記各超音波画像がカラーフローマッピング画像である場合には、血流信号を含む領域を被計測領域として抽出すると共に、当該被計測領域の短手方向幅を前記特徴量として計算することを特徴とする請求項1乃至3のうちいずれか一項記載の超音波診断装置。前記判定手段は、前記対象血管の長手方向に対する前記特徴量の変化を示すグラフを前記選択された超音波画像毎に生成し、前記生成した複数の前記グラフに基づいて、前記最適画像を判定すること、を特徴とする請求項1乃至5のうちいずれか一項記載の超音波診断装置。前記判定手段は、前記選択された各超音波画像上の前記対象血管の断面を楕円フィッティングして得られる楕円長軸長と楕円短軸長とを用いて、前記選択された複数の超音波画像の中から前記最適画像を判定することを特徴とする請求項1乃至5のうちいずれか一項記載の超音波診断装置。対象血管を含む二次元断面を超音波で繰り返し走査することで得られた複数の超音波データを記憶する記憶手段と、前記複数の超音波データを用いて複数の超音波画像を生成する画像生成手段と、前記各超音波画像上の少なくとも二か所において、前記対象血管の走行方向に対する中心軸と前記二次元断面とのずれを判定するための特徴量を計算する計算手段と、前記少なくとも二か所における前記特徴量に基づいて、前記複数の超音波画像の中から、最適画像を判定する判定手段と、前記最適画像を用いて、前記対象血管の中膜厚を計測する計測手段と、を具備することを特徴とする超音波画像処理装置。対象血管を含む二次元断面を繰り返し走査し、複数の画像データを取得するデータ取得手段と、前記複数の画像データを用いて複数の画像を生成する画像生成手段と、前記各超画像上の少なくとも二か所において、前記対象血管の走行方向に対する中心軸と前記二次元断面とのずれを判定するための特徴量を計算する計算手段と、前記少なくとも二か所における前記特徴量に基づいて、前記複数の画像の中から、最適画像を判定する判定手段と、前記最適画像を用いて、前記対象血管の中膜厚を計測する計測手段と、を具備することを特徴とする医用画像診断装置。コンピュータに、対象血管を含む二次元断面を超音波で繰り返し走査することで得られた複数の超音波データを用いて、複数の超音波画像を生成させる画像生成機能と、前記各超音波画像上の少なくとも二か所において、前記対象血管の走行方向に対する中心軸と前記二次元断面とのずれを判定するための特徴量を計算させる計算機能と、前記少なくとも二か所における前記特徴量に基づいて、前記複数の超音波画像の中から、最適画像を判定させる判定機能と、前記最適画像を用いて、前記対象血管の中膜厚を計測させる計測機能と、を実現させることを特徴とする超音波画像処理プログラム。超音波診断装置に用いるIMT計測では、超音波走査面が対象血管の中心軸に実質的に沿うように超音波走査面が対象血管の中心軸を実質的に含むように、的確に超音波プローブを位置決めする必要がある。しかしながら、従来の超音波診断装置では、操作者は、超音波走査面が対象血管の中心軸を実質的に含むか否かの判断を自身の感覚だけで行わなければならない。このため、例えば、超音波走査面が対象血管の中心軸に実質的に沿っていない場合超音波走査面が対象血管の中心軸を実質的に含まないには、計測結果を過大評価することがあった。また、部位によっては、IMT計測の対象とする血管の直径例えば約1cm程度と超音波プローブの超音波送受信面の縦幅とは同程度となる場合がある。係る場合、超音波走査面を正確に設定することは、経験の乏しい操作者にとって相当困難な作業であると言える。さらに、超音波プローブ位置の微調整より、超音波走査面が対象血管の中心軸を実質的に含むようにできた場合であっても、当該位置によって取得された複数の画像からIMT計測に最適な画像最適画像を選択することは、操作者の経験に依存するところが大きい。なお、例えば血管弾性計測のため、反射波強度に基づいて超音波プローブの位置を照準する手法例えば、反射波強度最大となる位置が血管中心軸に沿った位置等も提案されている例えば、特許文献1参照しかしながら、当該手法では、複数の超音波振動子を機械的に揺動させる必要があり、装置の構成が複雑になる。上記問題に鑑みてなされたものであり、操作者の経験に関わらず、客観性の高いIMT計測を比較的簡単な構成によって実行可能な超音波診断装置、超音波画像処置装置、医用画像診断装置及び超音波画像処理プログラムを提供することを目的としている。一実施形態に係る超音波診断装置は、対象血管を含む二次元断面を超音波で繰り返し走査し、複数の超音波データを取得するデータ取得手段と、前記複数の超音波データを用いて複数の超音波画像を生成する画像生成手段と、前記各超音波画像上の少なくとも二か所において、前記対象血管の走行方向に対する中心軸と前記二次元断面とのずれを判定するための特徴量を計算手段と、前記少なくとも二か所における前記特徴量に基づいて、前記複数の超音波画像の中から、最適画像を判定する判定手段と、前記最適画像を用いて、前記対象血管の中膜厚を計測する計測手段と、を具備するものである。以下、実施形態を図面に従って説明する。なお、以下の説明において、略同一の機能及び構成を有する構成要素については、同一符号を付し、重複説明は必要な場合にのみ行う。図1は、本実施形態に係る超音波診断装置1のブロック構成図を示している。同図に示すように、本超音波診断装置1は、超音波プローブ12、入力装置13、モニター14、超音波送信ユニット21、超音波受信ユニット22、Bモード処理ユニット23、血流検出ユニット24、画像生成処理ユニット25、計測処理ユニット27、表示処理ユニット28、制御プロセッサCPU29、記憶ユニット30、インタフェースユニット31を具備している。以下、個々の構成要素の機能について説明する。超音波プローブ12は、被検体に対して超音波を送信し、当該送信した超音波に基づく被検体からの反射波を受信するデバイス探触子であり、その先端に複数に配列された圧電振動子、整合層バッキング材等を有している。圧電振動子は、超音波プローブ12は、超音波送信ユニット21からの駆動信号に基づきスキャン領域内の所望の方向に超音波を送信し、当該被検体からの反射波を電気信号に変換する。整合層は、当該圧電振動子に設けられ、超音波エネルギーを効率良く伝播させるための中間層である。バッキング材は、当該圧電振動子から後方への超音波の伝播を防止する。当該超音波プローブ12から被検体Pに超音波が送信されると、当該送信超音波は、体内組織の音響インピーダンスの不連続面で次々と反射され、エコー信号として超音波プローブ12に受信される。このエコー信号の振幅は、反射することになった不連続面における音響インピーダンスの差に依存する。また、送信された超音波パルスが、移動している血流で反射された場合のエコーは、ドプラ効果により移動体の超音波送信方向の速度成分に依存して、周波数偏移を受ける。なお、本実施形態に係る超音波プローブ12は、一次元アレイプローブ複数の超音波振動子が所定の方向に沿って一次元的に配列されたプローブであるとする。しかしながら、当該例に拘泥されず、超音波プローブ12は、例えば二次元アレイプローブ複数の超音波振動子が二次元マトリックス状に配列されたプローブであってもよい。入力装置13は、装置本体11に接続され、オペレータからの各種指示、条件、関心領域ROIの設定指示、種々の画質条件設定指示等を装置本体11にとりこむための各種スイッチ、ボタン、トラックボール、マウス、キーボード等を有している。例えば、操作者が入力装置13の終了ボタンやFREEZEボタンを操作すると、超音波の送受信は終了し、当該超音波診断装置は一時停止状態となる。モニター14は、画像処理ユニット28からのビデオ信号に基づいて、生体内の形態学的情報や、血流情報を画像として表示する。超音波送信ユニット21は、図示しないトリガ発生回路、遅延回路およびパルサ回路等を有している。トリガ発生回路では、所定のレート周波数frHz周期;1/fr秒で、送信超音波を形成するためのトリガパルスが繰り返し発生される。また、遅延回路では、チャンネル毎に超音波をビーム状に集束し且つ送信指向性を決定するのに必要な遅延時間が、各トリガパルスに与えられる。パルサ回路は、このトリガパルスに基づくタイミングで、プローブ12に駆動パルスを印加する。なお、超音波送信ユニット21は、制御プロセッサ28の指示に従って所定のスキャンシーケンスを実行するために、送信周波数、送信駆動電圧等を瞬時に変更可能な機能を有している。特に送信駆動電圧の変更については、瞬間にその値を切り替え可能なリニアアンプ型の発信回路、又は複数の電源ユニットを電気的に切り替える機構によって実現される。超音波受信ユニット22は、図示していないアンプ回路、A/D変換器、加算器等を有している。アンプ回路では、プローブ12を介して取り込まれたエコー信号をチャンネル毎に増幅する。A/D変換器では、増幅されたエコー信号に対し受信指向性を決定し、受信ダイナミックフォーカスを行うのに必要な遅延時間を与え、その後加算器において加算処理を行う。この加算により、エコー信号の受信指向性に応じた方向からの反射成分が強調され、受信指向性と送信指向性とにより超音波送受信の総合的なビームが形成される。Bモード処理ユニット23は、受信ユニット22からエコー信号を受け取り、対数増幅、包絡線検波処理などを施し、信号強度が輝度の明るさで表現されるデータを生成する。血流検出ユニット24は、受信ユニット22から受け取ったエコー信号から血流信号を検出し、血流データを生成する。血流信号の検出は、通常CFMColorFlowMappingで行われる。この場合、血流信号を解析し、血流データとして平均速度、分散、パワー等の血流情報を多点について求める。画像生成処理ユニット25は、Bモード処理ユニット23、血流検出ユニット24から出力されるデータを用いて、画像データを生成する。また、画像生成処理ユニット25は、Bモード処理ユニット23、血流検出ユニット24から出力されるデータを用いてボリュームデータを生成し、ボリュームレンダリング、多断面変換表示MPR、最大値投影表示MIP:maximumintensityprojection等の所定の画像処理を行う。なお、ノイズ低減や画像の繋がりを良くすることを目的として、画像生成処理ユニット25の後に二次元的なフィルタを挿入し、空間的なスムージングを行うようにしてもよい。計測処理ユニット27は、後述するIMT計測支援機能において、所定の計測処理等を実行する。表示処理ユニット28は、画像処理ユニット28において生成処理された各種画像データに対し、ダイナミックレンジ、輝度ブライトネス、コントラスト、γカーブ補正、RGB変換等の各種を実行する。制御プロセッサ29は、情報処理装置計算機としての機能を持ち、本超音波診断装置本体の動作を制御する。制御プロセッサ29は、記憶ユニット31から後述するIMT計測支援機能を実現するための専用プログラムを読み出して自身が有するメモリ上に展開し、各種処理に関する演算制御等を実行する。記憶ユニット30は、後述するIMT計測支援機能を実現するための専用プログラムや、診断情報患者ID、医師の所見等、診断プロトコル、送受信条件、ボディマーク生成プログラムその他のデータ群が保管されている。また、必要に応じて、RAWデータメモリ中の画像の保管などにも使用される。記憶ユニット30のデータは、インタフェースユニット31を経由して外部周辺装置へ転送することも可能となっている。インタフェースユニット31は、入力装置13、ネットワーク、新たな外部記憶装置図示せずに関するインタフェースである。当該装置によって得られた超音波画像等のデータや解析結果等は、インタフェースユニット31よって、ネットワークを介して他の装置に転送可能である。IMT計測支援機能次に、本超音波診断装置1が有するIMT計測支援機能について説明する。この機能は、IMT計測を行う場合に、超音波走査面が対象血管の中心軸を実質的に含む画像を、対象血管に関する客観的な定量値特徴量に基づいて判定するものである。なお、本実施形態においては、説明を具体的にするために、超音波診断装置においてIMT計測支援機能を適用する場合を例とする。しかしながら、当該例に拘泥されず、例えば、X線コンピュータ断層撮像装置、磁気共鳴イメージング装置、X線診断装置、核医学診断装置等の医用画像診断装置や、医用ワークステーション医用画像参照装置においても本IMT計測支援機能は実現可能である。また、本実施形態では、説明を具体的にするため、本IMT計測支援処理を、画像生成処理ユニット25において処理を受けた二次元画像データBモードデータ、CFMデータ等を用いて実行する場合を例として説明する。しかしながら、処理対象となるデータは、Bモード処理ユニット23、血流検出ユニット23から出力された各RAWデータを用いて実行するようにしてもよい。図2は、本IMT計測支援機能に従う処理IMT計測支援処理の流れを示したフローチャートである。以下、各ステップにおける処理の内容について説明する。[患者情報等を入力:ステップS1]入力装置13を介して患者情報の入力、送受信条件被走査領域の大きさを決めるための画角、焦点位置、送信電圧等、対象血管を含む二次元断面を超音波走査するための撮像モードBモード、CFMモード等、IMT計測支援処理を実行するためのモードIMT計測支援処理モード等の選択が実行されるステップS1。入力、選択された各種情報条件モード等は、自動的に記憶ユニット31に記憶される。[二次元画像データの取得表示:ステップS2]超音波プローブ12が被検体表面の所望の位置に当接され、対象血管を含む二次元断面を被走査領域として、Bモード或いはCFMモードによる超音波走査が実行される。Bモードによる超音波走査によって取得されたエコー信号は、逐次超音波受信ユニット22を経由してBモード処理ユニット23、血流検出ユニット24に送られる。Bモード処理ユニット23は、対数増幅処理、包絡線検波処理等を実行し複数のBモードデータを生成する。また、血流検出ユニット24は、受信ユニット22から受け取ったエコー信号から血流信号を検出し、血流データを生成する。画像生成処理ユニット25は、Bモードデータ或いは血流データを用いて、Bモード画像或いは血流画像を生成する。生成された各画像は、表示処理ユニット28において所定の表示処理を受け、モニター14において所定の形態でリアルタイム表示される。操作者は、リアルタイムで表示される二次元画像を観察しながら図3に示す様に超音波プローブ12の位置を調整し、超音波走査面が対象血管の中心軸を実質的に含む様にするステップS2。[対象血管画像の保存:ステップS3]次に、ECG信号をトリガとして、対象血管に関する複数フレーム分の超音波画像対象血管画像が、画像生成処理ユニット25が有するシネメモリに保存される。このとき、ECG信号を基準とする画像毎の時間情報も各画像と対応付けて保存されるステップS3。[同一時相画像の選択:ステップS4]計測処理ユニット27は、ECG信号を基準とする時間情報に基づいて、シネメモリに保存された複数フレーム分の対象血管画像から、同一時相例えば、R波から所定期間後の時相の複数の対象血管画像を選択する図6参照。[特徴量の計算:ステップS5]計測処理ユニット27は、選択された同一時相の各対象血管画像における被計測領域を抽出する。この被計測領域の抽出は、次の様な手法によって実現することができる。例えば被計測領域を管腔領域とする場合、Bモード画像においては、被計測領域の輝度値或いは画素値は周辺の組織領域の輝度値より低く、被計測領域とそれ以外の領域との境界での輝度値変化率画素値変化率は他の位置での輝度値変化率よりも高い。従って、対象血管画像がBモード画像である場合には、輝度値又は輝度値変化率の適切な値による閾値処理により、管腔領域を抽出することができる。また、対象血管画像がCFM画像である場合には、血流信号を含む領域を被計測領域として抽出することができる。なお、被計測領域の抽出対象となる血管には、図4に示す様に幅が実質的に均一なもの、図5に示す様にプラークを有し幅が不均一なもの等が存在する。後者の場合、プラークの隣接領域の情報を用いた補間処理によりプラーク領域を除外する、被計測領域を管腔領域だけではなく血管内中膜も含めた領域として定義する等の手法により、プラーク形状の影響を低減するようにしてもよい。また、これらの例に拘泥されず、対象血管に関する客観的な定量値に応じて、任意の被計測領域を採用することが可能である。また、計測処理ユニット27は、抽出された対象血管画像毎の各被計測領域の長手方向に関する少なくとも二か所以上の位置において、対象血管の走行方向に対する中心軸と二次元断面とのずれを判定するための特徴量を計算するステップS5。なお、本実施形態では、図6に示す様に、上記特徴量として対象血管の短手方向幅を採用し、特徴量を被計測領域の長手方向に関する三か所の位置において計算するものとする。最適化画像の判定:ステップS6]計測処理ユニット27は、計算で得られた複数の特徴量を用いて、対象血管の長手方向に対する特徴量の変化を表すグラフ特徴量曲線を、選択された同一時相の対象血管画像毎に生成する。図7に、同一時相の対象血管画像毎の特徴量曲線の一例を示した。計測処理ユニット27は、複数の特徴量曲線のうち最大となる特徴量最大特徴量を持つ特徴量曲線を選択し、当該特徴量曲線に対応する対象血管画像をIMT計測に最適な画像最適画像として判定するステップS6。ここで、最大特徴量を持つ特徴量曲線に対応する画像を最適画像とするのは、IMT計測に最適な画像であれば、同一時相の対象血管画像の全てを通じて最大特徴量ここでは対象血管の短手方向幅の最大値、すなわち血管直径を有するからである。なお、最大特徴量を有する特徴量曲線が複数ある場合には、例えば次の手法によって最適画像を判定する。図8に示す様に、実際の血管の立体形状は円筒の形で近似することができることができる。仮に、超音波走査面が血管中心軸Aを完全に含む場合には、対象血管画像に含まれる対象血管の断面或いは被計測領域は実質的に矩形状となる。このため、三か所で計測された各特徴量は、実質的に同じ値となる。一方、断面が図8のS1~S3に示す様に血管中心軸Aを横切る場合には、その横切る位置においてのみ最大特徴量となる。従って、最大特徴量を有する複数の特徴量曲線のうち、長手方向に沿った変化率が最も少ないものに対応する対象血管画像を最適画像とすることができる。[IMT計測結果の表示:ステップS7]計測処理ユニット27は、判定された最適画像を用いてIMT計測を実行する。得られた結果は、モニター14に所定の形態で表示される。以上述べた本超音波診断装置によれば、対象血管画像毎の各被計測領域の長手方向に関する少なくとも二か所以上の位置において、対象血管の走行方向に対する中心軸と二次元断面とのずれを判定するための特徴量を計算し、当該特徴量に基づいてIMT計測に最適な画像を判定する。従って、操作者の経験を問わず、客観的な指標に基づいて最適な画像を常に用いてIMT計測を行うことができる。第2の実施形態次に、第2の実施形態について説明する。本実施形態に係る超音波診断装置は、楕円フィッティングを用いて最適画像を判定するものである。図9に示す様に、中心軸Aを包含しない平面によって切断される円筒の断面は、理論的には、図10に示すような楕円形を呈する。一方、切断面が完全に円筒の中心軸を通る場合には、切断面としての楕円の長軸長が無限大と考えることができる。これらの観点から、平面による血管断面と中心軸の位置関係を楕円フィッティングの長軸長で評価することができる。すなわち、ステップS5において選択した同一時相の血管対象画像の全てにつき、被計測領域の長手方向の長さと短手方向の幅すなわち、血管直径を計算する。その後、図10に示す様に、各血管対象画像の血管断面を楕円フィッティングし、得られる楕円短軸長と血管直径から得られる血管半径とを比較する。その結果、短軸長が血管半径より短い場合には、血管断面は血管中心軸を通らないと判定することができる。一方、楕円短軸長と血管半径とが実質的に等しい場合、或いは長軸長が事前設定された閾値より短い場合等には、血管断面は中心軸とある角度で交差すると判定することができる。また、楕円短軸長が血管半径と等しく、長軸長が事前設定された閾値より長い場合、血管断面は中心軸Aを実質的に含むものとみなすことができる。この様に、血管断面の楕円フィッティングによって得られる楕円短軸長及び長軸長を基準としても、最適画像を判定することが可能である。なお、上記各実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。例えば、本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスクフロッピー登録商標ディスク、ハードディスクなど、光ディスクCD-ROM、DVDなど、半導体メモリなどの記録媒体に格納して頒布することも可能である。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。