"被検体の所定部位を撮像することで取得された複数の時相に対応する複数のボリュームデータを記憶する記憶ユニットと、前記複数のボリュームデータ間の位置合わせを行うことで、前記各ボリュームデータに含まれる各ボクセルにつき、空間的な移動ベクトルを算出する算出ユニットと、前記各ボクセルについての前記移動ベクトルを用いて、前記診断部位の局所的な動きを示す画像を生成する画像生成ユニットと、前記診断部位の局所的な動きを示す画像を表示する表示ユニットと、を具備する画像処理装置。前記算出ユニットは、前記複数のボリュームデータのうちから選択された所定のボリュームデータを基準とした位置合わせを行うことで、前記移動ベクトルを算出する請求項1記載の画像処理装置。前記算出ユニットは、時間的に隣り合うボリューム間において位置合わせを行うことで、前記移動ベクトルを算出する請求項1記載の画像処理装置。前記算出ユニットは、前記移動ベクトルを規格化する請求項1乃至3のうちいずれか一項記載の画像処理装置。前記画像生成ユニットは、前記移動ベクトルの各成分の大きさに応じて、成分毎に異なる色彩及び輝度を割り当てることで、前記診断部位の局所的な動きを示す画像を生成する請求項1乃至4のうちいずれか一項記載の画像処理装置。記前記画像生成ユニットは、前記移動ベクトルの各成分をRG,Bに割り当てることで、前記診断部位の局所的な動きを示す画像を生成する請求項5記載の画像処理装置。前記画像生成ユニットは、前記各移動ベクトルを空間位置毎に矢印で図示することで、前記診断部位の局所的な動きを示す画像を生成する請求項1乃至4のうちいずれか一項記載の画像処理装置。前記画像生成ユニットは、第1の時相における前記所定部位の輪郭と、前記第1の時相から所定期間経過後の第2の時相における前記所定部位の輪郭と、を映像化することで、前記診断部位の局所的な動きを示す画像を生成する請求項1乃至4のうちいずれか一項記載の画像処理装置。前記画像生成ユニットは、前記診断部位の局所的な動きを示す画像をMPR画像として生成する請求項1乃至8のうちいずれか一項記載の画像処理装置。被検体の所定部位に関するボリュームデータを複数の時相に亘って撮像する撮像ユニットと、前記複数のボリュームデータ間の位置合わせを行うことで、前記各ボリュームデータに含まれる各ボクセルにつき、時相間における空間的な移動ベクトルを算出する算出ユニットと、前記各ボクセルについての前記移動ベクトルを用いて、前記診断部位の局所的な動きを示す画像を生成する画像生成ユニットと、前記診断部位の局所的な動きを示す画像を表示する表示ユニットと、を具備するX線コンピュータ断層撮像装置。前記診断部位は胸膜、気管支、細気管支の少なくともいずれかである請求項1乃至10のうちいずれか一項記載の画像処理装置。コンピュータに、画像処理装置を用いて、被検体の所定部位を撮像することで取得された複数の時相に対応する複数のボリュームデータ間の位置合わせを実行させる位置合わせ機能と、前記各ボリュームデータに含まれる各ボクセルにつき、時相間における空間的な移動ベクトルを算出させる算出機能と、前記各ボクセルについての前記移動ベクトルを用いて、前記診断部位の局所的な動きを示す画像を生成させる生成機能と、前記診断部位の局所的な動きを示す画像を表示させる表示機能と、を実現させる画像処理プログラム。しかしながら、前述したカラーマップは、全ての時間情報を基に作られたデータであり、表示時にはその時間毎の情報が失われている。また、MPRや3D画像の動画では、被検体の各部位それぞれの動きがどのように動いているのか、把握しづらいという問題があった。上記実情に鑑みてなされたものであり、その目的は、各部位それぞれの動きがどのように動いているのか、容易に把握することのできる医用画像を用いた臓器動態の画像処理装置、X線コンピュータ断層撮像装置及び画像処理プログラムを提供することである。本発明は、上記目的を達成するため、次のような手段を講じている。請求項1に記載の発明は、被検体の所定部位を撮像することで取得された複数の時相に対応する複数のボリュームデータを記憶する記憶ユニットと、前記複数のボリュームデータ間の位置合わせを行うことで、前記各ボリュームデータに含まれる各ボクセルにつき、空間的な移動ベクトルを算出する算出ユニットと、前記各ボクセルについての前記移動ベクトルを用いて、前記診断部位の局所的な動きを示す画像を生成する画像生成ユニットと、前記診断部位の局所的な動きを示す画像を表示する表示ユニットと、を具備する画像処理装置である。各部位それぞれの動きがどのように動いているのか、容易に把握することのできる医用画像を用いた臓器動態の画像処理装置、X線コンピュータ断層撮像装置及び画像処理プログラムを提供することができる。以下、図面を参照して実施形態を説明する。なお、以下の説明においては、X線コンピュータ断層撮影装置を例とする実施形態について説明する。しかしながら、当該例に拘泥されず、他の医用画像診断装置例えば、磁気共鳴イメージング装置、超音波診断装置、X線診断装置、核医学診断装置、医用画像診断装置によって取得される医用画像を用いた画像処理装置に関する実施形態も、実現可能である。図1は、一実施形態に係るX線コンピュータ断層撮影装置の構成を示したブロック図である。尚、X線コンピュータ断層撮影装置1には、X線管とX線検出器とが一体として被検体の周囲を回転する回転/回転方式と、リング上に多数のX線検出器が配置され、X線管のみが被検体の周囲を回転する固定/回転方式とがあり、何れの方式にも適用可能である。回転/回転方式に関して、一対のX線管とX線検出器とが回転フレームに搭載された一管球型と、X線管とX線検出器との対が回転フレームに複数搭載されたいわゆる多管球型とがあるが、何れの型でも適用可能である。X線検出器に関して、被検体を透過したX線をシンチレータ等の蛍光体で光に変換してからフォトダイオード等の光電変換素子で電荷に変換する間接変換形と、X線による半導体内の電子正孔対の発生及びその電極への移動すなわち光導電現象を利用した直接変換形とがあるが、いずれの形を採用してもよい。図1に於いて、架台10は、X線管11を有する。X線管11は、高電圧発生装置からスリップリング機構18を経由して管電圧の印加及びフィラメント電流の供給を受け、コーンビーム形のX線を発生する。X線管11は、回転軸Z軸を中心に回転自在に支持される回転フレーム12に、X線検出器13と共に搭載される。呼吸センサ17は、被検体の呼吸の動きを検出するために設けられている。尚、呼吸センサ17は、検査対象に応じて、被検体の心臓の動き心拍の時相を検出する心電計または拍動センサに置換されることができる。X線検出器13は、X線管11から被検体を透過したX線を検出する。X線検出器13は、コーンビームに対応するマルチスライス型または2次元アレイ型である。すなわち、X線検出器13は、回転軸にそって並列される複数のX線検出素子列を有する。各X線検出素子列は、回転軸RAに直交する方向に沿って一列に配列された複数のX線検出素子を有する。X線検出器13の出力は、チャンネル毎にデータ収集回路DAS14によって増幅され、デジタル信号に変換され、そして例えば非接触型のデータ伝送装置15を介して前処理装置16に送られ、そこで感度補正等の補正処理を受け、再構成処理の直前段階にあるいわゆる投影データとして投影データ/画像データ記憶部23に、それを収集した時期に対応する呼吸時相コードと共に記憶される。スキャンコントローラ20は、データ収集スキャンのために、回転駆動部、高電圧発生装置19、データ収集回路14、投影データ/画像データ記憶部23等を制御する。再構成部24は、ダイナミックスキャニングにより繰り返し収集された投影データに基づいて、呼吸時相の異なる複数の2次元または3次元の画像のデータを再構成する。呼吸時相の異なる複数の2次元又は3次元の画像のデータは、投影データ/画像データ記憶部23に、その再構成処理に用いた投影データのセットの例えば中心時相に対応する呼吸時相コードと共に記憶される。3次元の画像再構成処理法としては、典型的には、フェルドカンプ法Feldkampmethodである。フェルドカンプ法は、周知の通り、ファンビームコンボリューションバックプロジェクション法を基にした近似的再構成法であり、コンボリューション処理は、コーン角が比較的小さいことを前提として、データをファン投影データと見なして行われる。しかし、バックプロジェクション処理は、実際のレイrayに沿って行われる。つまり、投影データにZ座標に依存した重みをかけ、その重みをかけた投影データに、ファンビーム再構成と同じ再構成関数をコンボリューションし、そのデータをコーン角を持つ斜めの実際のレイに沿って逆投影する、という手順で画像が再構成される。前述したように、本実施形態に係るX線CT装置1は画像処理装置を有する。画像処理装置は、上記投影データ/画像データ記憶部23と共に、特定位相決定部25、表示処理部26、基準点決定部27、関心点決定部28、ベクトル処理部29、肺機能指標計算部30、画像処理部31を有して構成される。特定位相決定部25は、投影データと共に記憶された流量時間曲線の例えば極大点と極小点の特定により最大吸気時相と最大呼気時相とを決定する。基準点決定部27は、前記最大吸気時相の画像と最大呼気時相の画像それぞれに対して、解剖学上同一部位に基準点を設定する。また、位置合わせの基準となるデータを設定する機能を有している。関心点決定部28は、例えば、肺野上例えば、胸膜、気管支、細部気管支等の肺組織に複数の関心点を設定する。呼吸に伴う基準点に対する関心点の移動距離等から肺機能指標が得られる。複数の関心点は、最大吸気時相の画像と、最大呼気時相の画像それぞれに対して設定される。また、複数の関心点は、肺の壁輪郭上、結節部上、腫瘍上に設定される。関心点決定部28は、最大吸気時相の画像と最大呼気時相の画像それぞれから、閾値処理、例えば領域成長処理regiongrowingにより肺領域を抽出する。抽出した肺領域の壁輪郭上に、基準点から一定角度毎に関心点を設定する。ベクトル処理部29は、最大吸気時相の画像上の複数の関心点各々についてベクトルを計算する。同様に、ベクトル処理部29は、最大呼気時相の画像上の複数の関心点各々についてもベクトルを計算する。更に、ベクトル処理部29は、最大吸気時相の画像上の複数の関心点についての複数のベクトルと、最大呼気時相の画像上の複数の関心点についての複数のベクトルとの間で、同じ角度同士でベクトル差分を計算する。つまり、呼吸運動による各関心点の移動距離が、比較的固定的な基準点を基準として定量的に求められる。また、ベクトル処理部29は、非線形位置合わせ時の変形量から、各ボクセルに於ける3次元の動きベクトルx→,y→,z→を算出する機能を有する。肺機能指標計算部30は、計算された呼吸運動による複数の関心点各々の移動距離から、各時相の肺体積の定量値、肺体積の変化率、その変化体積の定量値等の肺機能指標を計算する。表示処理部26は、計算された肺機能指標を画像と共に数値として、または画像の対応位置に指標値に応じた色相や輝度に対応付けて表示させるために必要な処理を行い、これを表示する。画像処理部31は、医用画像や部位モデル等の各種画像を処理する。この画像処理部31は、図示されないが、ソフトウェアまたはハードウェア回路、或いはそれらの両方により構成されているもので、画像及びモデルの位置合わせ機能を有している。更に、画像処理部31は、前記ベクトル処理部29で算出された動きベクトルの成分それぞれを正規化し、R,G,Bカラーに割り付ける機能と、基準となるデータについて非線形位置合わせを行う機能を有している。次に、図2のフローチャートを参照して、造影された血管について処理する場合の動作について説明する。図1に示されるX線CT装置1を用いて得られた造影された画像が読み込まれる。尚、ここで読み込まれるデータは、少なくとも2位相以上のデータとする。続くステップS2にて、前記ステップS1で読み込まれたデータから、位置合わせの基準となるボリュームが選択される。これは、操作者による手動の選択でもよいし、各位相間の動きの差分等から、動きが最も少ないボリュームを自動で検出する方法等が用いられてもよい図3参照。なお、以下の説明においては、この基準ボリュームに対応する時相を基準時相と呼ぶことにする。次に、ステップS3に於いて、前記ステップS2で設定された基準ボリュームを基に、各位相について位置合わせ処理が行われる。位相間の位置合わせには、周知の線形位置合わせや、非線形位置合わせが用いられる。したがって、これらの詳細な説明は省略する。次に、ステップS4にて、位置合わせする際の変形量から、各ボクセルに於ける動きベクトル成分x→y→,z→がベクトル処理部29にて算出される。例えば、図4Aに示されるようなボクセル351から、所定時間経過後の図4Bに示されるボクセル352との間で、あるボクセル35a1の移動量に注目する。すると、ボクセル35a1から35a2までの移動量Vn→については、図4Cに示されるような、動きベクトル成分Vn→=x→ny→n,z→nが算出される。その後、ステップS5にて、前記ステップS4で得られた各ボクセルの動きベクトルの成分について正規化が行われ、画像処理部31にて、3次元上のxyz軸に対してR,G,Bカラーの割り付けすなわち、方向に応じて異なる色彩の割り当てると共に、各成分の大きさに対応した輝度値を割り当てる処理が行われる。割り付けられた画像は、新規ボリュームとして保存される。R,G,Bへの割り付けには、絶対値|x→||y→|,|z→|を算出してもよいし、符号を持った値をそのまま割り付けてもよい。符号を持っている場合は、x→y→,z→ベクトルの最小値から最大値までの中間値が、R,G,Bでも中間値となる。次いで、ステップS6に於いて、前述したステップS3~S5の処理動作が全ての元ボリュームについて行われ、最終ボリュームでない場合はステップS7に移行して、次のボリュームへシフトする。そして、位置合わせが行われるべく前記ステップS3に移行する。このようにして、動きベクトルの成分をRGB化したボリュームが作成される。最後に、ステップS8にて、前述したステップS3~S6で得られたボリュームを用いてMPR画像が生成され、表示される。これは、一位相で表示しても構わないし、動画で表示しても構わない。図5A、図5Bは、本実施形態に於けるX線CT装置1による画像処理を行った画像の例を示した図である。尚、図5A、図5Bに於いては、便宜上、R成分を○印、G成分を□印、B成分を△印で表しており、それぞれの成分が密になっている部分はその色の成分が強い濃い状態を表し、当該成分が疎になっている部分はその色の成分が弱い薄い状態を表しているものとする。例えば、図5Aに示される肺の基準画像411に対して入力画像412を合わせると、結果画像413は、全体としてRの成分が多いことがわかる。つまり、Rの成分はx軸方向に相当するので、この結果画像413から、肺は、上側が主にx軸方向に動いていると判断される。また、図5Bに示される肺の基準画像431に対して入力画像432を合わせると、結果画像433は、全体としてRの成分が多いが、右下側がややBの成分が多く表示されていることがわかる。つまり、全体としてはx軸方向R成分に動いているものの、右下側はz軸方向B成分に動いていると判断される。このように、3次元上の画像の動きをR,G,Bカラーに割り当てて、動きベクトル成分をその方向に分けてR,G,B成分で表示するようにしたので、被検体の各部位それぞれの動きが容易に把握することができる。また、例えば癌部位は動きにくいため、前述したR,G,Bカラー割り付けした場合、比較的黒く表示され、他の部位との識別ができると考えられる。第2の実施形態次に、第2の実施形態について説明する。本実施形態では、診断部位の局所的な動きを、各位置における矢印によってベクトル場として表示ベクトル場表示するものである。図6は、第2の実施形態に係る動き情報の映像化処理の流れを示したフローチャートである。ステップS1~ステップS4までの処理は、図2に示した第1の実施形態における処理を同様である。以下、ステップS5a、5bの各処理について説明する。画像処理部31は、算出された時相毎の各ボクセルの動きベクトルを用いて、予め選択された時相選択時相に対応する診断部位の局所的な動きを示す画像を生成する。具体的には、画像処理部31は、選択時相のボリュームデータに対して所定のMPR断面を設定し、算出された選択時相のボクセル毎の動きベクトルを用いて、当該MPR断面上の各位置における移動方向及び移動量を矢印で示したベクトル場画像を生成する。表示処理部26は、生成された選択時相に対応するベクトル場画像を所定の形態で表示するステップS5a。図7、図8は、それぞれステップS5aにおいて生成され表示されるベクトル場画像の例を示した図である。各図に示すように、MPR画像の各位置につき、基準時相から選択時相までの期間において移動が発生した場合には、移動方向が矢印の向きによって、移動量が矢印の長さによって明示されている。また、動きが発生していない場合には、移動方向及び移動量を示す矢印が割り当てられていない。従って、ユーザは、表示されたベクトル場画像を観察することで、被検体の診断部位の各位置の動きを容易に把握することができる。続いて、他の時相のベクトル場画像の生成指示他の時相の選択指示が入力された場合には、新たに選択された時相につき、ステップS3~S5aの処理が繰り返し実行される。また、他の時相の選択指示が入力されない場合には、動き情報の映像化処理を終了するステップS6a。第3の実施形態次に、第3の実施形態について説明する。本実施形態では、診断部位の局所的な動きを、ワイヤーフレーム、サーフェスモデル等によって表示するものである。図9は、第3の実施形態に係る動き情報の映像化処理の流れを示したフローチャートである。ステップS1~ステップS4までの処理は、図2に示した第1の実施形態における処理を同様である。以下、ステップS5b、6bの各処理について説明する。画像処理部31は、基準ボリュームを用いて、診断部位の輪郭をワイヤーフレームで示すワイヤーフレーム画像或いは、診断部位の表面を示すサーフェスモデル画像を生成する。同様に、画像処理部31は、予め選択された時相選択時相に対応するボリュームを用いて、診断部位の輪郭をワイヤーフレームで示すワイヤーフレーム画像或いは、診断部位の表面を示すサーフェスモデル画像を生成する。また、画像処理部31は、生成された各ワイヤーフレーム画像と、基準時相から選択時相までの期間における各ボクセルについての移動ベクトルとを用いて、診断部位の局所的な動きを示す画像を生成する。すなわち、画像処理部31は、基準時相に対応するワイヤーフレームと選択時相に対応するワイヤーフレームとを位置合わせすると共に、選択時相に対応するワイヤーフレームに動きを示す情報を付加することで、診断部位の局所的な動きを示す画像を生成する。表示処理部26は、生成された選択時相に対応するワイヤーフレーム画像を所定の形態で表示するステップS5b。なお、選択時相に対応するワイヤーフレームに付加される動きを示す情報は、どのようなものであってもよい。例えば、ワイヤーフレーム上の各位置について、第1の実施形態において説明したR,G,Bの割り当てを行ってもよい。また、ワイヤーフレーム上の各位置について、第2の実施形態において説明したベクトル場表示を行うようにしてもよい。さらに、各位置毎の動き表示ではなく、所定の範囲について平均的な動きの大きさ、動きの方向を矢印等によって表示するようにしてもよい。図10は、ステップS5bにおいて生成され表示されるワイヤーフレーム画像の例を示した図である。同図に示すように、ワイヤーフレーム画像を観察することで、輪郭上の各位置が基準時相から選択時相までの期間においてどのように移動したのかを容易且つ迅速に視認することができる。続いて、他の時相のワイヤーフレーム画像の生成指示他の時相の選択指示が入力された場合には、新たに選択された時相につき、ステップS3~S5bの処理が繰り返し実行される。また、他の時相の選択指示が入力されない場合には、動き情報の映像化処理を終了するステップS6b。以上述べた各実施形態において、各部位に対して1つの結果画像を表示するようにしている。しかしながら、これに限られるものではなく、例えば、オーバーレイ表示されるものであってもよい。また、各実施形態では静止画を例として説明したが、これに限られるものではなく、動画でも適用が可能である。また、上記実施形態においては、ステップS2において基準ボリュームを設定し、当該基準ボリュームに対する位置合わせを行うことで、各ボクセルについて基準ボリュームに対応する時相からの動きベクトルを算出した。しかしながら、当該例に拘泥されず、例えば時間的に隣り合うボリューム間において位置合わせを行うことで、各ボクセルについての時相間の移動ベクトルを算出するようにしてもよい。なお、本発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合わせにより、種々の発明を形成できる。例えば、本実施形態に係る各機能は、当該処理を実行するプログラムをワークステーション等のコンピュータにインストールし、これらをメモリ上で展開することによっても実現することができる。このとき、コンピュータに当該手法を実行させることのできるプログラムは、磁気ディスクフロッピー登録商標ディスク、ハードディスクなど、光ディスクCD-ROM、DVDなど、半導体メモリなどの記録媒体に格納して頒布することも可能である。また、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態にわたる構成要素を適宜組み合わせてもよい。"