"周期的運動をする心臓を撮像した時系列の動画像を入力し、前記動画像から前記心臓の位置、輪郭、長軸角度、又は、スケールの少なくとも一つを含む構造情報を、前記心臓の特定時相のタイミングで検出する構造情報検出部と、前記特定時相における前記構造情報と、前記特定時相における1周期前、又は、1周期後の特定時相における前記構造情報との差であるずれ量を算出するずれ量算出部と、前記動画像が診断に有用であるかどうかを示すように、前記ずれ量を正規化した有用度を算出する有用度算出部と、を有することを特徴とする医用画像診断装置。前記特定時相における前記構造情報の位置、輪郭、長軸角度、又は、スケールの近傍に、前記構造情報を検出するための第1探索範囲を設定し、前記特定時相における前記1周期前、又は、前記1周期後の前記特定時相における前記第2探索範囲を、前記第1探索範囲より狭い領域に設定する探索範囲設定部をさらに有する、ことを特徴とする請求項1に記載の医用画像診断装置。前記心臓の心電図を測定する心電計、又は、前記動画像から前記特定時相を検出する時相検出部をさらに有する、ことを特徴とする請求項1に記載の医用画像診断装置。前記構造情報検出部が、前記1周期前、又は、前記1周期後の前記特定時相における前記構造情報を検出できなかった場合、前記ずれ量算出部は、前記ずれ量を無限大に設定する、ことを特徴とする請求項1に記載の医用画像診断装置。前記構造情報検出部が、前記特定時相において前記構造情報を初めて検出する場合、又は、検出できなかった場合、前記探索範囲設定部は、予め設定されている領域を前記探索範囲とする、ことを特徴とする請求項2に記載の医用画像診断装置。前記医用画像診断装置が、超音波プローブによって前記動画像を取得する画像取得部を有する超音波診断装置である、ことを特徴とする請求項1に記載の医用画像診断装置。前記特定時相検出部は、複数組以上の前記特定時相を検出し、前記構造情報検出部は、前記各組の前記特定時相間の前記構造情報をそれぞれ算出し、前記ずれ量算出部は、前記各組の前記特定時相間で検出された前記構造情報の前記ずれ量をそれぞれ算出し、前記有用度算出部は、前記各組の前記特定時相間で前記有用度をそれぞれ算出し、前記各組の前記有用度の中で最も低い前記有用度を選択する、ことを特徴とする請求項1に記載の医用画像診断装置。前記構造情報の前記位置とは、前記心臓の心尖位置、心房の中心位置、心室の中心位置、大動脈の中心位置、又は、弁輪位置であり、前記輪郭とは、前記心臓の左心室内膜境界、左心室外膜境界、左心房内膜境界、右心室内膜境界、又は、右心房内膜境界である、ことを特徴とする請求項1に記載の医用画像診断装置。前記有用度を、ユーザに通知する出力部をさらに有する、ことを特徴とする請求項1に記載の医用画像診断装置。しかし、撮影中に心臓の断面がずれてしまった場合には、得られた画像の有用性を正しく評価できないという問題点があった。本発明が解決しようとする課題は、心臓の1周期分の動画像において、有用性を正しく算出できる医用画像診断装置の提供を目的とする。実施形態の医用画像診断装置は、構造情報検出部と、ずれ量算出部と、有用度算出部とを持つ。構造情報検出部は周期的運動をする心臓を撮像した時系列の動画像を入力し、前記動画像から前記心臓の位置、輪郭、長軸角度、又は、スケールの少なくとも一つを含む構造情報を、前記心臓の特定時相のタイミングで検出する。ずれ量算出部は前記特定時相における前記構造情報と、前記特定時相における1周期前、又は、1周期後の特定時相における前記構造情報との差であるずれ量を算出する。有用度算出部は前記動画像が診断に有用であるかどうかを示すように、前記ずれ量を正規化した有用度を算出する。以下、実施形態である医用画像診断装置以下、単に診断装置という10について図面を参照して説明する。以下、第1の実施形態の診断装置10について図1、図2を参照して説明する。本実施形態の診断装置10の構成について図1を参照して説明する。図1は、診断装置10のブロック図である。図1に示すように、診断装置10は、構造情報検出部1、ずれ量算出部2、有用度算出部3を有する。周期運動をする心臓の1周期分の時系列の動画像が、診断装置10に入力する。この動画像は、X線CT、MRI、超音波診断装置等の医用画像撮影装置から入力するが、画像サーバ、CD、DVD等のメディア、又は、ネットワークストレージ等から入力してもよい。構造情報検出部1は、入力された動画像は、1周期分の動画像であるので、最初と最後の画像フレームを特定時相に設定し、最初のフレーム1周期前の特定時相と最後のフレーム現在の特定時相における心臓の構造情報を検出し、その構造情報は、ずれ量算出部2に出力する。ずれ量算出部2は、入力された最初のフレーム1周期前の特定時相の構造情報と、最後のフレーム現在の特定時相の構造情報との差からずれ量を算出し、そのずれ量を有用度算出部3に出力する。有用度算出部3は、入力されたずれ量から、2つの前記特定時相間の有用度を算出する。診断装置10の動作について図1と図2を参照して説明する。なお、図2は、診断装置10の動作を示すフローチャートである。まず、ステップS1では、構造情報検出部1が、入力された周期運動をする心臓の1周期分の動画像から、最初のフレームと最後のフレームにおける心臓の構造情報を検出する。最初のフレームと最後のフレームとは、心臓の周期運動の同一特定時相の1周期ずれた画像である。構造情報とは、心臓の位置、輪郭、姿勢、特徴点、境界等を意味し、詳しくは、心房の中心位置、心室の中心位置、大動脈の中心位置、心臓の長軸角度、心臓のスケール、心尖位置、弁輪位置、又は、前記心臓の輪郭を表す左心室内膜境界、左心室外膜境界、左心房内膜境界、右心室内膜境界、右心房内膜境界の少なくとも一つを意味する。また、弁輪位置とは僧帽弁の付け根であって、大動脈弁の付け根である。本実施形態では、左心室中心位置、長軸角度、スケールを構造情報とした場合と、左心室内膜境界を構造情報とした場合について説明する。スケールとは長軸や短軸などの軸の長さを示す。構造情報検出部1が、左心室中心位置、長軸角度、スケールに関する構造情報を検出するときは、左心室周辺の画像から切り出した特徴量を用いる。特徴量としては、画素の輝度値や輝度勾配を用いる。これら構造情報を検出する場合は、構造情報検出部1は、非特許文献1TatsuoKozakayaetalFullyAutomaticFeatureLocalizationforMedicalImagesusingaGlobalVectorConcentrationApproachCVPR2008に記載されている識別器を予め学習しておき、その識別器を用いて画像中から構造情報を検出する。また、左心室内膜境界に関する構造情報は、複数の制御点によって表わされる。構造情報検出部1が、この境界を検出するときは、非特許文献2TECootesetalActiveshapemodels、TheirtrainingandapplicationCVIU1995に記載された動的輪郭モデルを用いる。すなわち、構造情報検出部1は、予め学習された輪郭形状モデルと境界周辺の輝度情報とを利用して第1形状を変形させていき、境界を求める。次に、ステップS2では、ずれ量算出部2が、検出された最初のフレーム1周期前の特定時相の構造情報と、最後のフレーム現在の特定時相の構造情報との差から、2つの特定時相で検出した構造情報間のずれ量を算出する。例えば、構造情報が左心室中心位置Cx,y、長軸角度θ、スケールsで表される場合、ずれ量算出部2は、ずれ量Dを式1で算出する。なお、x,yは、画像に設定された2次元のx-y直交座標軸上の座標値である。また、構造情報が境界の制御点Xx,yに関する点群で表される場合は、ずれ量算出部2は、ずれ量Dを式2で算出する。なお、ずれ量Dは式2で示すような各点の二乗誤差の総和に限定されず、例えば、絶対値誤差の総和、制御点と境界の距離の総和等の2つの構造情報の差を表す指標であればよい。次に、ステップS3では、有用度算出部3が、算出されたずれ量に基づいて、2つの特定時相間の動画像の有用度を算出する。ずれ量は物理的な量であり、直感的に理解しづらいため、ユーザ例えば、医師に通知するためにずれ量Dを正規化した値有用度に変換する。有用度算出部3は、ずれ量Dが小さければ高い有用度を与え、大きければ低い有用度を与える。有用度算出部3は、ずれ量Dから有用度への変換には予め用意しておいた変換関数を用いてもよいし、ずれ量Dの値と有用度を対応付けておいたテーブルを用いてもよい。また、有用度は連続値でもよいし、閾値を設けてそれ以上ならば白丸又は使える、以下ならば×又は使えないといった2値の表現でもよい。例えば、変換関数には式3を用いる。但し、Eは有用度、Dはずれ量である。式3では、構造情報が一致していればずれ量Dが0となり、有用度は最大のE=1となる。また、ずれ量Dが大きくなるほど有用度は低くなり、未検出の場合は有用度E=0となる。このように、ずれ量Dを0<=有用度E<=1で正規化できる。本実施形態によれば、撮影中に操作者の操作ミスや被験者の動きによって、撮影している心臓の断面がずれてしまった場合に、そのずれ量を検出でき、入力された動画像が診断に有用であるかどうかをユーザに通知できる。すなわち、診断装置10は、複数枚からなる1周期分の動画像の診断の有用性を算出できる。以下、第2の実施形態の診断装置10について図3~図5を参照して説明する。本実施形態の診断装置10の構成について図3を参照して説明する。図3は、本実施形態に係わる診断装置10を示すブロック図である。診断装置10は、構造情報検出部1、ずれ量算出部2、有用度算出部3、探索範囲設定部5、特定時相検出部4、出力部6を有する。特定時相検出部4は、画像情報、又は、生体情報ECG、心音、脈拍等から特定時相のタイミング時間を検出し、その特定時相のタイミングを構造情報検出部1と探索範囲設定部5へ出力する。構造情報検出部1は、第1の実施形態と同様に入力された動画像から、特定時相検出部4で検出された特定時相における心臓の構造情報を、後から説明する探索範囲設定部5で設定された探索範囲内から検出し、その構造情報を、ずれ量算出部2、探索範囲設定部5に出力する。ずれ量算出部2は、現在の特定時相の前記構造情報と1周期前の特定時相の構造情報の差からずれ量を算出し、そのずれ量を有用度算出部3に出力する。有用度算出部3は、入力されたずれ量から有用度を算出し、その有用度を出力部6に出力する。探索範囲設定部5は、特定時相検出部4で検出された特定時相において、構造情報検出部1で検出された心臓の位置等の近傍を、1周期前の特定時相の探索範囲として設定し、その探索範囲を構造情報検出部1に出力する。ここで位置の近傍とは、その位置を中心、又は、基準点として、予め定められた範囲を意味し、長軸角度の近傍とは、その角度を中心として、予め定められた角度の範囲を意味し、スケールの近傍とは、そのスケールを中心として、予め定められた範囲を意味し、輪郭例えば、左心室内膜境界の近傍とは、境界の各制御点を中心、又は、基準点として予め定められた範囲を意味する。出力部6は、液晶表示装置やCRT等のディスプレイであって、有用度、心電図波形、動画像等を表示する。次に、診断装置10の動作について図3と図4を参照して説明する。図4は、診断装置10の動作を示すフローチャートである。なお、本実施形態では、構造情報としては、左心室内膜境界の場合について説明する。左心室内膜境界は、複数の制御点によって表わされる。まず、ステップS11では、特定時相検出部4が、入力された周期運動を繰り返す心臓の特定時相のタイミングを検出し、構造情報検出部1へ出力する。特定時相とは、左心室の拡張末期、又は、左心室の収縮末期のタイミングを意味し、心電計より推定可能な時相であって、心電計より推定可能な時相から固定時間進めた時相である。本実施形態では、特定時相を左心室の拡張末期の時相とした場合について説明する。また、特定時相は、心電計に代えて、心臓の動画像から推定してもよい。次に、ステップS12では、探索範囲設定部5が、心臓の構造情報を検出するときの探索範囲を設定し、構造情報検出部1へ出力する。すなわち、探索範囲設定部5は、現在の特定時相拡張末期の構造情報に関する第1探索範囲を決定し、1周期前の特定時相1周期前の拡張末期の構造情報に関する探索範囲を決定する。なお、探索範囲設定部5は、初めて探索を行う場合は、予め設定されている領域を第1探索範囲として設定する。また、探索範囲設定部5は、1周期前の拡張末期の特定時相において左心室内膜境界が検出されている場合は、その境界の各制御点の近傍で、かつ、第1探索範囲より狭い領域を第2探索範囲として設定する。例えば、第2探索範囲は式4のように設定する。さらに、探索範囲設定部5は、1周期前の拡張末期の特定時相において左心室内膜境界が検出されていない場合は、予め設定されている領域を第2探索範囲として設定する。このように、現在の特定時相の1周期前の第2探索範囲を第1探索範囲より狭くすることにより、探索の処理量削減と検出ミス削減に寄与できる。次に、ステップS13では、構造情報検出部1が、第1の実施形態で説明したように、設定された探索範囲内から心臓の構造情報を検出し、ずれ量算出部2へ出力する。次に、ステップS14では、ずれ量算出部2が、検出された拡張末期の特定時相の左心室内膜境界と、1周期前の拡張末期の特定時相における左心室内膜境界の差から、2つの特定時相間の構造情報のずれ量Dを算出し、有用度算出部3へ出力する。ずれ量Dの計算方法は第1の実施形態のステップS2と同様であるが、設定された探索範囲内から左心室内膜境界が検出されなかった場合は、ずれ量算出部2は、未検出としてずれ量Dを無限大に設定する。次に、ステップS15では、有用度算出部3が、ずれ量Dに基づいて、2つの前記特定時相間の有用度を算出し、出力部6へ出力する。有用度の計算方法は第1の実施形態のステップS3と同様であるので説明は省略する。最後に、ステップS16では、出力部6が、算出された有用度をユーザに通知する。出力部6は、図5に示す画像を表示する。図5aに示す心臓の動きの周期を表す心電波形と、それに算出された有用度を重ねて表示したものであって、有用度が高い場合は白丸、低い場合は×で表している。図5bは、入力された動画像を表示したものである。また、出力部6は、有用度に応じてそれぞれの周期の心電波形を色分けして表示してもよい。本実施形態によれば、撮影中に操作者の操作ミスや被験者の動きによって撮影している心臓の断面がずれてしまった場合に、そのずれ量を検出でき、入力された動画像が診断に有用であるかどうかをユーザに通知できる。また、現在の特定時相の構造情報を用いて、1周期前の特定時相における探索範囲を狭くすることにより、探索の処理時間削減と検出ミスの削減ができる。以下、第3の実施形態の診断装置10について図6~図10を参照して説明する。本実施形態の診断装置10の構成について図6を参照して説明する。図6は、診断装置10を示すブロック図である。図6に示すように、本実施形態は、第2の実施形態の診断装置10の特定時相検出部4とずれ量算出部2との接続を追加した構成となる。この構成にすることで、ずれ量算出部2で現在注目しているフレームの特定時相を判断できる。そのため、1周期間に複数組の特定時相を設定した場合でも、それぞれの1周期前の構造情報のずれ量を算出できる。他の構成については第2の実施形態と同様であるため説明は省略する。次に、診断装置10の動作について図6~図10を参照して説明する。図7は、診断装置10の動作を示すフローチャートである。まず、ステップS21では、特定時相検出部4が、入力された周期運動を繰り返す心臓の特定時相のタイミングを1周期間から2つ組以上検出し、構造情報検出部1、ずれ量算出部2へ出力する。本実施形態では、構造情報を左心室内膜境界として、特定時相は図8に示すように1周期間からA、Bの2組の特定時相を選択した場合について説明する。次に、ステップS22では、探索範囲設定部5が、心臓の構造情報を検出するときの探索範囲を設定し、構造情報検出部1へ出力する。次に、ステップS23では、構造情報検出部1が、設定された探索範囲内から心臓の構造情報を検出し、ずれ量算出部2へ出力する。ステップS22、ステップS23は、第2の実施形態におけるステップS12、ステップS13の動作と同様であるため、説明は省略する。次に、ステップ24では、ずれ量算出部2が、現在の特定時相における左心室内膜境界と、1周期前の特定時相現在注目しているフレームが図8の特定時相A2ならば特定時相A1、特定時相B2ならば特定時相B1における左心室内膜境界とから、2組の特定時相間の構造情報のそれぞれのずれ量を算出し、有用度算出部3へ出力する。それぞれのずれ量の算出方法は第1の実施形態のステップS2と同様であるため、説明は省略する。次に、ステップS25では、有用度算出部3が、それぞれのずれ量に基づいて、2組の特定時相間の有用度をそれぞれ算出する。ステップS25は、第2の実施形態におけるステップS15の動作と同様であるため、説明を省略する。次に、ステップS’25では、有用度算出部3が、2組の特定時相間に対する有用度を統合し、出力部6へ出力する。有用度の統合は、例えば図9に示すように行う。具体的には、有用度算出部3は、A2-A3間が白丸であった場合、A2-A3間と一部重複するB1-B2間、B2-B3間のどちらか一方が×であれば統合された有用度を×と算出し図9a参照、双方とも白丸であれば統合された有用度を白丸と算出する図9b参照。また、有用度が連続値である場合は、有用度算出部3は有用度の低い値を採用する。すなわち、有用度算出部3は、各組の特定時相間で有用度をそれぞれ算出し、各組の有用度の中で最も低い前記有用度を選択する。このように、1周期間に2組の特定時相を設定し、それらの有用度を統合することによって、1周期中で断面が一旦ずれた後に同一位置に復帰した場合であっても断面のずれを検出できる。最後に、ステップS26では、出力部6は、算出された有用度をユーザに通知する。出力部6は、図10に示すような画像を表示する。図10aは、心臓の動きの周期を表す心電波形と、それに算出された特定時相A間、特定時相B間の有用度と、統合された有用度をそれぞれ重ねて表示したものである。図10bは、入力された動画像を表示したものである。また、出力部6は、有用度に応じてそれぞれの周期の心電波形を色分けして表示してもよい。本実施形態によれば、撮影中の心臓の断面のずれをより精度よく検出でき、よりロバストな動画像の評価を実現できる。なお、上記実施形態では、2組の特定時相間で説明したが、これに限らず3組以上の特定時相間で有用度を求めてもよい。上記各実施形態では、ずれ量Dは、動画像における任意の特定時相の構造情報と、前記任意の特定時相の1周期前の特定時相の構造情報との差で算出した。しかし、前記動画像を録画している場合は、ずれ量Dは、前記動画像における任意の特定時相の構造情報と、前記任意の特定時相の1周期後の特定時相の構造情報との差で算出してもよい。上記各実施形態の診断装置10は、例えば、汎用のコンピュータを基本ハードウェアとして用いることでも実現することが可能である。すなわち、構造情報検出部1、ずれ量算出部2、有用度算出部3、特定時相検出部4、及び、探索範囲設定部5は、上記のコンピュータに搭載されたプロセッサにプログラムを実行させることにより実現することができる。このとき、診断装置10は、上記のプログラムをコンピュータに予めインストールすることで実現してもよいし、CD-ROM等の記憶媒体に記憶して、あるいはネットワークを介して上記のプログラムを配布して、このプログラムをコンピュータに適宜インストールすることで実現してもよい。また、B及びCは、上記のコンピュータに内蔵あるいは外付けされたメモリ、ハードディスクもしくはCD-R、CD-RW、DVD-RAM、DVD-R等の記憶媒体等を適宜利用して実現することができる。上記では本発明の一実施形態を説明したが、この実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の主旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。"