"患者の所定の周期的運動を行う部位を表す3次元動画像および前記部位を表す超音波動画像を取得する画像取得手段と、前記取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から前記周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出する特徴部抽出手段と、前記3次元動画像および前記超音波動画像の前記周期的運動のフェーズを取得するフェーズ取得手段であって、前記3次元動画像および前記超音波動画像の少なくとも一方については前記抽出された特徴部の形態に基づいて前記フェーズを取得するフェーズ取得手段と、前記抽出した特徴部と前記取得したフェーズに基づいて、前記3次元動画像および前記超音波動画像を構成する各フレーム画像の前記特徴部の位置および前記フェーズを対応付ける対応付け手段と、前記対応付けられた特徴部の位置およびフェーズに基づいて、前記3次元動画像および前記超音波動画像の前記対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成する画像生成手段と、前記生成された各比較画像を対比可能に表示装置に表示する表示制御手段とを備えたことを特徴とする画像処理装置。前記特徴部抽出手段は、前記特徴部を自動抽出するものであることを特徴とする請求項1記載の画像処理装置。前記周期的運動を行う部位は心臓であり、前記特徴部は、心室、心房、心筋、心弁、心尖のいずれかであることを特徴とする請求項1または2記載の画像処理装置。前記フェーズ取得手段が前記心弁の開閉状態から前記フェーズを取得するものであることを特徴とする請求項3記載の画像処理装置。前記フェーズ取得手段が前記心弁のうち僧帽弁および/または大動脈弁の開閉状態から前記心臓の拡張期および/または収縮期の終了時を特定することにより前記フェーズを取得するものであることを特徴とする請求項4記載の画像処理装置。前記フェーズ取得手段は、前記3次元動画像および前記超音波動画像のうち少なくとも一方について、該動画像の付帯情報からフェーズを取得するものであることを特徴とする請求項1から5いずれか1項記載の画像処理装置。前記画像生成手段は、前記3次元動画像および前記超音波動画像の付帯情報から画素間隔を取得し、取得した前記画素間隔に基づいて前記両動画像の画素間隔を一致させて前記比較画像を生成するものであることを特徴とする請求項1から6のいずれか1項記載の画像処理装置。前記画像生成手段は、血流のドップラーシフトに基づくカラードップラー法により表された前記比較画像を生成するものであることを特徴とする請求項1から7のいずれか1項記載の画像処理装置。前記超音波動画像が、前記特徴部を含む断面を表す動画像であることを特徴とする請求項1から8のいずれか1項記載の画像処理装置。患者の所定の周期的運動を行う部位を表す3次元動画像および前記部位を表す超音波動画像を取得し、前記取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から前記周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出し、前記3次元動画像および前記超音波動画像の前記周期的運動のフェーズを取得するとともに、前記3次元動画像および前記超音波動画像の少なくとも一方については前記抽出された特徴部の形態に基づいて前記フェーズを取得し、前記抽出した特徴部と前記取得したフェーズに基づいて、前記3次元動画像および前記超音波動画像を構成する各フレーム画像の前記特徴部の位置および前記フェーズを対応付けし、前記対応付けられた特徴部の位置およびフェーズに基づいて、前記3次元動画像および前記超音波動画像の前記対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成し、前記生成された各比較画像を対比可能に表示装置に表示することを特徴とする画像処理方法。コンピュータを、患者の所定の周期的運動を行う部位を表す3次元動画像および前記部位を表す超音波動画像を取得する画像取得手段と、前記取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から前記周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出する特徴部抽出手段と、前記3次元動画像および前記超音波動画像の前記周期的運動のフェーズを取得するフェーズ取得手段であって、前記3次元動画像および前記超音波動画像の少なくとも一方については前記抽出された特徴部の形態に基づいて前記フェーズを取得するフェーズ取得手段と、前記抽出した特徴部と前記取得したフェーズに基づいて、前記3次元動画像および前記超音波動画像を構成する各フレーム画像の前記特徴部の位置および前記フェーズを対応付ける対応付け手段と、前記対応付けられた特徴部の位置およびフェーズに基づいて、前記3次元動画像および前記超音波動画像の前記対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成する画像生成手段と、前記生成された各比較画像を対比可能に表示装置に表示する表示制御手段として機能させることを特徴とする画像処理プログラム。しかし、特許文献1に記載された技術によれば、静止画についてCT装置により撮影された3次元画像と超音波画像を空間的に位置合わせすることができても、CT装置により撮影された3次元の動画像と超音波動画像を心臓の拍動のフェーズを対応付けることができないため、両動画像間で対応するフェーズおよび空間的位置を一致させて比較可能に表示することはできない。また、特許文献2に記載した技術によれば、拍動のフェーズを一致させるために両動画像の心電図データを取得する必要があるため、いずれかの画像または両動画像の心電図データが利用できない場合には、両動画像の対応するフェーズを対応付けることが難しい。そこで、本発明では、CT装置またはMR装置により撮影された3次元動画像と超音波動画像を、周期的運動のフェーズおよび位置を一致させた比較画像を生成して表示することを容易にする画像処理装置、方法およびプログラムを提供することを目的とする。本発明による画像処理装置は、患者の所定の周期的運動を行う部位を表す3次元動画像および前記部位を表す超音波動画像を取得する画像取得手段と、前記取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から前記周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出する特徴部抽出手段と、前記3次元動画像および前記超音波動画像の前記周期的運動のフェーズを取得するフェーズ取得手段であって、前記3次元動画像および前記超音波動画像の少なくとも一方については前記抽出された特徴部の形態に基づいて前記フェーズを取得するフェーズ取得手段と、前記抽出した特徴部と前記取得したフェーズに基づいて、前記3次元動画像および前記超音波動画像を構成する各フレーム画像の前記特徴部の位置および前記フェーズを対応付ける対応付け手段と、前記対応付けられた特徴部の位置およびフェーズに基づいて、前記3次元動画像および前記超音波動画像の前記対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成する画像生成手段と、前記生成された各比較画像を対比可能に表示装置に表示する表示制御手段とを備えたことを特徴とするものである。本発明による画像処理方法は、患者の所定の周期的運動を行う部位を表す3次元動画像および前記部位を表す超音波動画像を取得し、前記取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から前記周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出し、前記3次元動画像および前記超音波動画像の前記周期的運動のフェーズを取得するとともに、前記3次元動画像および前記超音波動画像の少なくとも一方については前記抽出された特徴部の形態に基づいて前記フェーズを取得し、前記抽出した特徴部と前記取得したフェーズに基づいて、前記3次元動画像および前記超音波動画像を構成する各フレーム画像の前記特徴部の位置および前記フェーズを対応付けし、前記対応付けられた特徴部の位置およびフェーズに基づいて、前記3次元動画像および前記超音波動画像の前記対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成し、前記生成された各比較画像を対比可能に表示装置に表示することを特徴とするものである。本発明による画像処理プログラムは、コンピュータを、患者の所定の周期的運動を行う部位を表す3次元動画像および前記部位を表す超音波動画像を取得する画像取得手段と、前記取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から前記周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出する特徴部抽出手段と、前記3次元動画像および前記超音波動画像の前記周期的運動のフェーズを取得するフェーズ取得手段であって、前記3次元動画像および前記超音波動画像の少なくとも一方については前記抽出された特徴部の形態に基づいて前記フェーズを取得するフェーズ取得手段と、前記抽出した特徴部と前記取得したフェーズに基づいて、前記3次元動画像および前記超音波動画像を構成する各フレーム画像の前記特徴部の位置および前記フェーズを対応付ける対応付け手段と、前記対応付けられた特徴部の位置およびフェーズに基づいて、前記3次元動画像および前記超音波動画像の前記対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成する画像生成手段と、前記生成された各比較画像を対比可能に表示装置に表示する表示制御手段として機能させることを特徴とするものである。ここで、本発明における部位は、所定の周期的運動を行う部位であれば何でもよく、代表的には心臓があげられる。また、所定の周期的運動とは、部位に含まれる各部分がそれぞれ所定の方向に所定の範囲で動作する再現性のある運動であれば何でもよく、心臓の拍動でも、肺の呼吸でも、関節の屈伸運動でもよい。また、本発明における前記周期的運動を行う部位が心臓である場合、前記特徴部は、心室、心房、心筋、心弁、心尖のいずれかであってもよい。この場合、例えば、前記フェーズ取得手段が前記心弁の開閉状態から前記フェーズを取得するものであることが好ましく、フェーズ取得手段が前記心弁のうち僧帽弁および/または大動脈弁の開閉状態から前記心臓の拡張期および/または収縮期の終了時を特定することにより前記フェーズを取得するものであることがさらに好ましい。また、本発明における3次元動画像は、被写体の部位の形態を表す3次元動画像であればよく、例えば、CT画像やMRI画像により撮影された3次元動画像があげられる。また、本発明における前記超音波動画像は、3次元動画像であってもよく、あるいは前記特徴部を含む断面を表す動画像であってもよい。また、本発明における周期的運動のフェーズとは、周期的運動を複数の段階で表すもので、周期的運動の一周期を何段階に分けたものであってもよい。また、本発明における対比可能に表示装置に表示するとは、比較画像を対比させて表示させるものであればどのように表示してもよく、例えば、左右に並列表示してもよく、上下に並列表示してもよく、一方の比較画像中に他方の比較画像を子画像としてピクチャーインピクチャー方式で表示してもよい。また、1つのディスプレイに両比較画像を表示してもよく、複数のディスプレイに両比較画像を表示してもよい。また、両比較画像を同じ大きさで表示してもよく、異なる大きさで表示してもよい。また、本発明における前記特徴部抽出手段は、前記特徴部を自動抽出するものであることが好ましい。また、本発明における前記フェーズ取得手段は、前記3次元動画像および前記超音波動画像のうち少なくとも一方について、該動画像の付帯情報からフェーズを取得するものであってもよい。また、本発明における前記画像生成手段は、前記3次元動画像および前記超音波動画像の付帯情報から画素間隔を取得し、取得した前記画素間隔に基づいて前記両動画像の画素間隔を一致させて前記比較画像を生成するものであってもよい。また、本発明における前記画像生成手段は、血流のドップラーシフトに基づいてカラードップラー法により表された前記超音波動画像を比較画像として生成するものであることが好ましい。本発明の画像処理装置、方法およびプログラムによれば、患者の所定の周期的運動を行う部位を表す3次元動画像およびその部位を表す超音波動画像を取得し、取得した3次元動画像および超音波動画像を構成する複数のフレーム画像から周期的運動に応じて形態が変化する所定の特徴部をそれぞれ抽出し、3次元動画像および超音波動画像のそれぞれについて周期的運動のフェーズを取得するとともに両動画像の少なくとも一方については抽出された特徴部の形態に基づいてフェーズを取得し、抽出した特徴部と取得したフェーズに基づいて、3次元動画像および超音波動画像を構成する各フレーム画像の特徴部の位置およびフェーズを対応付けし、対応付けられた特徴部の位置およびフェーズに基づいて、3次元動画像および超音波動画像の対応付けられた特徴部の位置およびフェーズを一致させた比較画像をそれぞれ生成し、生成された比較画像を表示装置に表示するものであるため、心電図データが利用できない場合であっても特徴部の形態に基づいて容易に3次元動画像及び超音波動画像の比較画像を生成することができる。また、3次元動画像に基づいて生成された比較画像により高空間解像度で対象部位を把握するとともに必要に応じて超音波動画像に基づいて生成された比較画像より超音波動画像からのみ得られる情報を容易に把握することができる。このため、医師らの診断効率と診断精度の向上に資する。以下、本発明の画像処理装置、画像処理プログラムおよび画像処理方法の実施形態について、図面を参照しながら詳細に説明する。図1に、本発明の一実施形態における画像処理装置6を含む病院システム1の概略構成を示す。病院システム1は、ローカルエリアネットワークLAN2を介して互いに接続された検査室システム3、データサーバ4、および診断用ワークステーションWS6により構成される。検査室システム3は、被検体を撮影する各種のモダリティ32と、各モダリティから出力された画像の確認や調整を行う検査室ワークステーションWS31により構成される。モダリティ32としては、心臓の形態情報を表す形態画像を取得可能なモダリティであるCTComputedTomography装置、MRIMagneticResonanceImaging装置が備えられ、さらに超音波診断装置などが備えられているものとする。これらのモダリティ32のうちCTComputedTomography装置、MRIMagneticResonanceImaging装置は、いずれもDICOMDigitalImagingandCommunicationinMedicine規格に準拠した装置であり、取得したボリュームデータに付帯情報を付し、DICOMファイルとして出力する。モダリティ32から出力されたファイルは、検査室WS31により、データサーバ4に転送される。データサーバ4は、高性能プロセッサと大容量メモリを備えた比較的処理能力の高いコンピュータに、データベースマネージメントサーバDBMS:DatabaseManagementServerの機能を提供するソフトウェアプログラムを実装したものである。プログラムはストレージに記憶され、起動時にメモリにロードされ、プロセッサにより実行される。データサーバ4は、検査室WS31から転送されたファイルを、大容量ストレージ5に記憶させる。また、データサーバ4は、診断用WS6からの検索要求に応じて、大容量ストレージ5に記憶されている複数のファイルの中から、検索条件に適ったファイルを選出し、診断用WS6に送信する。診断用WS6は、標準的なプロセッサ、メモリおよびストレージを備えた汎用のワークステーションに、診断支援するための画像処理プログラムを実装したものである。画像処理プログラムは、DVD等の記録媒体から、またはネットワーク上のサーバコンピュータからのダウンロードにより、診断用WS6にインストールされる。また、診断用WS6には、ディスプレイ7と、マウス、キーボード等の入力装置8が接続されている。診断用WS6に実装される画像処理プログラムは、各種機能を実現するプログラムモジュール群により構成され、その中には、画像処理機能を実現するためのプログラムモジュール群が含まれている。これらのプログラムはストレージに記憶され、起動時にメモリにロードされ、プロセッサにより実行される。これにより、診断用WS6は、患者の所定の周期的運動を行う部位心臓を表す3次元動画像V1およびその部位を表す超音波動画像V2を取得する画像取得手段61と、取得した3次元動画像V1および超音波動画像V2を構成する複数のフレーム画像から周期的運動に応じて形態が変化する所定の特徴部僧帽弁MVをそれぞれ抽出する特徴部抽出手段62と、3次元動画像V1および超音波動画像V2の周期的運動のフェーズを取得するフェーズ取得手段であって、3次元動画像V1および超音波動画像V2の少なくとも一方については抽出された特徴部の形態に基づいてフェーズを取得するフェーズ取得手段63と、抽出した特徴部と取得したフェーズに基づいて、3次元動画像V1および超音波動画像V2を構成する各フレーム画像の特徴部の位置およびフェーズを対応付ける対応付け手段64と、対応付けられた特徴部の位置およびフェーズに基づいて、3次元動画像V1および超音波動画像V2の対応付けられた特徴部の位置およびフェーズを一致させた比較画像I1、I2をそれぞれ生成する画像生成手段65と、生成された比較画像I1、I2を表示装置7に表示する表示制御手段66として動作するようになる。図2は、本実施形態の画像処理の流れを表すフローチャートである。図3に、比較画像I1、I2の表示例を示す。図2および図3を用いて、本実施形態におけるWS6画像処理装置の各機能処理の流れについて、以下詳細に説明する。本実施形態では、心臓の検査を行う場合を例に説明する。本実施形態の処理に先立って、心臓の検査では、CT装置等を用いて1周期の拍動を含む被写体の胸部の動画撮影が行われ、撮影された3次元動画像V1ボリュームデータが付帯情報を付されDICOMファイルとしてデータサーバ4に転送され、大容量ストレージ5に記憶される。ボリュームデータは、3次元空間内の濃度や密度の分布を表す多数のボクセルデータの集合で構成され、各ボクセルデータ内にX線の吸収量等が画素値として表われる。また、超音波プローブを口から食道まで挿入して超音波撮影を行う経食道心エコー法Transesophagealechocardiography:TEEにより同じ被写体の胸部の動画撮影が行われ、撮影された3次元超音波動画像V2がデータサーバ4に転送され、大容量ストレージ5に記憶される。まず、初期画面において心臓の画像処理機能が選択され、所定の入力画面において患者の識別番号や検査番号などが入力されると、画像取得手段61は、入力された情報をデータサーバ4へと送信し、大容量ストレージ5に保存されているファイルの検索および転送を要求する。上記要求を受けたデータサーバ4は、大容量ストレージ5内のファイルを検索し、画像取得手段61に対し、要求されたファイルを転送する。画像取得手段61は、データサーバ4から転送されたファイルに含まれる3次元動画像V1および3次元超音波動画像V2を取得し、メモリに保存するS01。続いて特徴部抽出手段62は、3次元動画像V1および3次元超音波動画像V2からそれぞれ所定の特徴部として左心室LVと左心房LDの間に位置する心臓弁である僧帽弁MVを抽出するS02。ここでは、3次元動画像V1および3次元超音波動画像V2の特徴部抽出処理およびフェーズ取得処理に、RazvanIoanIonasecetalPatientSpecificModelingandQuantificationoftheAorticandMitralValvesFrom4DCardiacCTandTEEMedicalImagingIEEEvol292010pp16361651に示される方法を適用する。特徴部抽出手段62は、上記非特許文献に記載された方法により両動画像V1、V2について、少なくとも1回の拍動周期にわたって時系列に心臓の僧帽弁MVのセグメンテーションを行い、両動画像V1、V2を構成する各フレーム画像において僧帽弁MVの輪郭上の各サンプル点の位置を特定する情報を抽出する。そして、フェーズ取得手段63は、両動画像V1、V2について、僧帽弁MVの輪郭上の各サンプル点の位置に基づいて各フレーム画像の心臓の拍動のフェーズを取得するS03。ここで、心臓の拍動は、その1周期が収縮期と拡張期から構成される。そして、収縮期が終了時に、大動脈弁AVが開いた状態から閉じた状態になるとともに僧帽弁MVが閉じた状態から開き始める。また、拡張期の終了時に、僧帽弁MVが閉じた状態になるとともに大動脈弁が閉じた状態から開き始める。フェーズ取得手段63は、この性質を利用して、拡張期の終了時と収縮期の終了時を特定することにより、心臓の拍動のフェーズを特定する。本実施形態では、上記非特許文献に記載された方法により所定の時間間隔ごとに両動画像V1、V2の僧帽弁MVの形態を所定のパラメータによって取得し、この僧帽弁MVの形態に基づいて僧帽弁MVの開閉状態を特定し、僧帽弁MVの開閉状態に対応する心臓の拍動のフェーズを取得する。ここで、本実施形態では、両動画像V1、V2においてそれぞれ、僧帽弁MVの開閉状態と、僧帽弁MVの形態を表す所定のパラメータが対応付けられて記憶されている。フェーズ取得手段63は、各動画像V1、V2ごとに、僧帽弁MVが開いた状態から閉じた状態に変化したフレームを心臓の拍動の拡張期の終了時に対応するフレームとして特定する。また、フェーズ取得手段63は、各動画像V1、V2ごとに、僧帽弁MVが閉じた状態から開いた状態開き始めた状態に変化したフレームを心臓の拍動の収縮期の終了時に対応するフレームとして特定する。なお、例えば、僧帽弁MVの形態を表す所定のパラメータとして、僧帽弁MVの輪郭上の特定の複数のサンプル点間の距離などが考えられる。そして、対応付け手段64は、両動画像V1、V2について、両動画像V1、V2の収縮期および拡張期の終了時のフェーズがそれぞれ一致するように、必要に応じて時間軸方向の補間処理を行って、両動画像を構成する各フレームを時間的に互いに対応付ける。ここでは、同じフェーズを表すフレーム画像同士を互いに対応付けるとともに、対応付けられたフレーム画像同士で同一の特徴の空間的位置を対応付ける。S04。対応付け手段64は、上記処理において3次元動画像V1と3次元超音波動画像V2の1周期を表すフレーム数が異なる場合は、1周期を表すフレーム数が少ない方の動画像を基準として対応付けを行う。例えば、必要に応じて、1周期を表すフレーム数が多い動画像のフレームを適宜間引く処理を行う。また、対応付けるフレーム間でフェーズが少しずれている場合に、周知の方法により両動画像の対応するフレーム画像のフェーズを一致させるように補間処理を行ってもよい。例えば、一方の動画像を構成するフレーム画像のフェーズを取得し、他方の動画像のフレーム画像からこのフェーズの前後のフレーム画像を用いて、取得したフェーズに応じた形態の変形を周知の方法により補間したフレーム画像を生成し、フェーズを一致させて一方の動画像のフレーム画像と他方の動画像から生成したフレーム画像を対応付けることが考えられる。画像生成手段65は、上述の処理により3次元動画像V1から抽出した一連のフレーム画像についてボリュームレンダリング画像を生成するとともに、3次元超音波動画像V2から抽出した一連のフレーム画像について、対応付け手段64により対応付けられた特徴部の位置、方向、サイズが一致するように3次元超音波動画像V2の座標系を3次元動画像V1の座標系に変換した画像を生成し、周知の方法により生成した両動画像V1、V2の比較画像I1、I2を生成してストレージ5に記憶するS05。詳細には、画像生成手段65は、対応付け手段64が対応付けた位置に基づいて、同一の特徴の空間的位置を一致させるように一方の座標系を他方の座標系に変換し、変換した座標系において、同一の特徴の空間的位置、方向、サイズを一致させるように適宜補正を行って、空間的位置合わせを行う。なお、対応付け手段64は、空間的位置の位置合わせの際に、3次元動画像V1および3次元超音波動画像V2の画素範囲PixelSpacing情報をそれぞれDICOMのヘッダ情報から取得し、画素範囲情報に基づいて両動画像V1、V2を適宜拡大または縮小して、3次元動画像V1および3次元超音波動画像V2から抽出した一連のフレーム画像の画素範囲を一致させる。また、本実施形態における画像生成手段65は、図3に示すように、3次元動画像V1に基づく画素値を所定の透明度でボリュームレンダリング法により表す比較画像I1を生成するとともに、図3の矢印Cに示すように3次元超音波動画像V2に基づいて画素値および血流の方向を周知のドップラーカラー法により表す比較画像I2を生成する。なお、画像生成手段65は、3次元動画像V1および3次元超音波動画像V2から抽出した一連のフレーム画像の比較画像I1、I2の生成方法として、同一の特徴の空間的位置、方向、サイズを一致させるように両動画像V1、V2から抽出した一連のフレーム画像の比較表示を行えるものであれば、各種の周知の生成法を適用可能である。表示制御手段66は、画像生成手段65に生成された比較画像I1、I2を取得し、図3に示すように、比較画像I1、I2をディスプレイ7に表示させるS06。以上、本実施形態によれば、3次元動画像V1および超音波動画像V2の特徴部の形態に基づいて、所定の動作を行う部位の周期的運動のフェーズを取得し、両動画像のフェーズを一致させるとともに、両動画像の特徴部の位置に基づいて空間的な位置合わせを行うため、両動画像のいずれかが心電図データなどを利用できない場合であっても、両動画像を好適に対応付けることができる。また、対応付けた両動画像の比較画像I1、I2を対比可能に表示することにより、ユーザらは、3次元動画像V1に基づいた比較画像I1により高空間解像度で対象部位を把握するとともに、必要に応じて比較画像I2より超音波動画像からのみ得られる情報を容易に把握することができる。このため、効率よく正確に画像診断を行うことができる。また、本実施形態では、両動画像の比較画像I1、I2を同じ大きさで並列表示しているため、さらに両動画像を対比して観察しやすく、両動画像に表された観察対象をより容易かつ正確に把握することができる。本実施形態では、超音波動画像V2に基づいて血流のドップラーシフトに基づいてカラードップラー法により表された比較画像I2を生成表示するため、超音波動画像からのみ得られる血流情報を比較画像I2から容易に直感的に把握することができ、好ましい。また、本実施形態では、3次元動画像V1および3次元超音波動画像V2から自動認識処理により特徴部の形態を抽出するために、ユーザのマニュアル操作による特徴部の抽出の労力を要さず、効率よく容易に特徴部の形態が抽出できる。また、周期的運動を行う部位が心臓であり、所定の特徴部が心弁であるため、心臓の周期的運動を心弁の開閉状態により精度良く特定でき、フェーズを好適に取得することができる。また、フェーズ取得手段が心弁のうち僧帽弁の開閉状態から心臓の拡張期および収縮期の終了時を特定することによりフェーズを取得するため、心臓の拍動に応じた特徴部の形態の変化の性質を好適に利用して、心臓の周期的運動をより精度良く特定できる。また、本実施形態においては、対応付け手段64が、両動画像V1、V2について、複数の心臓の拍動のフェーズ収縮期の終了時および拡張期の終了時を一致させているため、より正確に両動画像V1、V2の対応付けができる。また、本実施形態では、3次元動画像V1および超音波動画像V2の両方について自動認識処理によりフェーズを取得するため、容易に精度良くフェーズを取得して対応付けることができる。また、3次元動画像および3次元超音波動画像の一方について、特徴部の形態により周期的運動のフェーズを特定し、他方についてDICOMなどのヘッダ情報によりフェーズを取得してもよい。この場合には、一方についてのみ自動認識処理を適用することにより計算負荷の増加を抑制できるとともに効率よく両動画像のフェーズを取得することができる。また、画像生成手段65は、3次元動画像V1および超音波動画像V2の付帯情報から画素間隔を取得し、取得した画素間隔に基づいて両動画像の画素間隔を一致させて比較画像I1、I2を生成するものであるため、両動画像について画素間隔を容易に取得できるため、容易に精度良く両動画像のサイズを一致させて比較画像I1、I2を生成することができる。また、本実施形態では、CT装置またはMRI装置により撮影された3次元動画像V1と3次元の超音波動画像V2とに基づいて、上記画像処理を行うため観察対象をより詳細に把握できる。また、本実施形態による特徴部抽出処理には、YefengZhengetalFourChamberHeartModelingandAutomaticSegmentationfor3DCardiacCTVolumesUsingMarginalSpaceLearningandSteerableFeaturesMedicalImagingIEEEvol272008pp16681681に示される方法を適用してもよい。なお、特徴部抽出手段62は、構造物の特徴部を2つの3次元動画像V1、V2から抽出できる方法であれば、周知の種々の方法を適用することができる。例えば、ユーザが手動で心弁部などの特徴部の位置や形状を各3次元動画像V1、V2ごとにマウス等により入力して、かかる入力を取得して、特徴部の位置や形状を抽出してもよい。また、フェーズ取得手段63は、収縮期が終了する際に、大動脈弁AVが開いた状態から閉じた状態になるとともに僧帽弁MVが閉じた状態から開き始め、拡張期の終了時には、僧帽弁MVが閉じた状態になるとともに大動脈弁が閉じた状態から開き始めるという性質を利用するいかなる方法で心臓の拍動のフェーズを判断してもよい。例えば、僧帽弁MVが閉じた状態から開き始めてから収縮期の終了時、再度閉じた状態から開き始める収縮期の終了時までを心臓の拍動の1周期として検出して、両画像V1、V2を収縮期の終了時のフェーズが一致するように対応付けてもよく、僧帽弁MVが開いた状態から閉じた状態になってから拡張期の終了時、再度僧帽弁MVが開いた状態から閉じた状態になる拡張期の終了時までを心臓の拍動の1周期として検出して、両画像V1、V2を拡張期の終了時のフェーズが一致するように対応付けてもよい。また、例えば、僧帽弁MVに替えて大動脈弁AVの開閉状態により、拍動のフェーズを判断してもよく、僧帽弁MVの開閉状態および大動脈弁AVの開閉状態の両方の情報を重み付けして用いることにより拍動のフェーズを判断してもよい。また、所定の特徴部として、図3に示すような左心室LV、左心房LA、右心室RV、右心房RA、心弁MV、AV、PV、TV、心尖AC、などの種々の特徴のいずれかまたはこれらの任意の組合せにより、本実施形態と同様に心臓の拍動のフェーズに応じた形態の周期的変化に基づいて、心臓の周期的運動を特定した場合には、複数の情報に基づくためより精度良くフェーズを取得することができる。また、フェーズ取得手段63は、両動画像V1、V2のうち、複数周期の周期的運動を含むものに関しては、複数の周期のうち、両動画像V1、V2の対応付けに用いる周期を任意に指定可能である。本実施形態によるフェーズ取得手段63は、マウスやキーボードによりユーザの入力を受け付けて、周知の任意の方法により複数の周期のうち1つの周期を特定するものとする。例えば、3次元動画像V1または3次元超音波動画像V2に含まれる複数の周期に対応する周期選択ボタンを表示し、ユーザの周期選択を受け付けてもよく、3次元動画像V1または3次元超音波動画像V2に含まれる複数の周期の開始時間をキーボード等が入力されるように促し、入力を受け付けて周期選択を行ってもよい。図4は、上記実施形態の変形例による比較画像I1、I2の表示例である。上記実施形態では、3次元超音波動画像V2を例に説明したが、当業者にとっては明らかであるように、図4に示すような断面のように、心室LV、RVや心房LA、RAや心弁MV、AV、PV、TVや心尖ACなどの部位に含まれる認識可能な特徴を表す断面を表すものであれば、超音波診断装置により撮影された2次元動画像でも同様に本発明を適用可能である。特徴部を含む所定の断面において、CT装置またはMR装置により撮影された高空間解像度の画像を観察できると同時に血流の情報など超音波診断装置により撮影された2次元動画像にのみ得られる情報も同時に把握できるため、容易かつ正確に画像診断を行うことができる。なお、本実施形態に限られず、所定の部位は、所定の周期的な運動を行う部位であれば何でもよく、例えば、膝の屈伸運動のような所定の周期的運動を行う関節のような部位であってもよい。膝などの屈伸運動に本発明を適用する場合には、膝関節を構成する1つ以上の部位をセグメント化し、例えば、大腿骨と頸骨の所定点間の距離など、セグメント化された各部位から膝の屈伸状態を表すパラメータを取得し、屈伸状態を表すパラメータに基づいて曲げた状態から伸ばした状態までのフェーズを取得することが考えられる。また、表示制御手段66は、比較画像を対比させて表示させるものであればどのように表示してもよく、例えば、左右に並列表示してもよく、上下に並列表示してもよく、一方の比較画像中に他方の比較画像を子画像としてピクチャーインピクチャー方式で表示してもよい。また、1つのディスプレイに両比較画像を表示してもよく、複数のディスプレイに両比較画像を表示してもよい。また、両比較画像を同じ大きさで表示してもよく、異なる大きさで表示してもよい。上記各実施形態では、両比較画像I1,I2を同じ大きさで並列表示したため、医師らは比較画像I1、I2を対比しやすく、効率よく正確に画像診断を行うことができる。なお、CT装置またはMRI装置により撮影された3次元動画像V1と超音波動画像V2との位置合わせは、超音波動画像V2の座標系にCT装置またはMRI装置により撮影された3次元動画像V1の座標系を変換するものであってもよい。また、対応付け手段64は、CT装置またはMRI装置により撮影された3次元動画像V1と超音波動画像V2とのフェーズの対応付けを、1回の周期的運動の一部についてのみ行ってもよく、1回の周期的運動について行ってもよく、複数の周期的運動について行ってもよい。また、各実施形態において、本発明の画像処理プログラムが一台の診断用WSに実装されて画像処理装置として機能する場合について説明したが、画像処理プログラムを複数台のコンピュータに分散させてインストールして、複数台のコンピュータが画像処理装置として機能するようにしてもよい。"