医療用画像処理装置及び医療用画像処理プログラム
【課題】動き解析情報及びジオメトリ情報に基づき、複数のボリュームデータからそれぞれジオメトリ情報を作成する処理の労力を軽減できる医療用画像処理装置及び医療用画像処理プログラムを提供する。
【解決手段】医療用画像処理装置は、時刻を有する複数のボリュームデータを可視化する医療用画像処理装置であって、前記複数のボリュームデータについて、ジオメトリ情報を生成する第１生成手段と、前記一のボリュームデータと他の少なくとも一つのボリュームデータとの間の動き解析情報を生成する第２生成手段と、前記他の少なくとも一つのボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する新たなジオメトリ情報を生成する第３生成手段と、前記他の少なくとも一つのボリュームデータについて、前記対応する新たなジオメトリ情報を含む画像を表示する処理を行う表示処理手段と、を備える。
【解決手段】医療用画像処理装置は、時刻を有する複数のボリュームデータを可視化する医療用画像処理装置であって、前記複数のボリュームデータについて、ジオメトリ情報を生成する第１生成手段と、前記一のボリュームデータと他の少なくとも一つのボリュームデータとの間の動き解析情報を生成する第２生成手段と、前記他の少なくとも一つのボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する新たなジオメトリ情報を生成する第３生成手段と、前記他の少なくとも一つのボリュームデータについて、前記対応する新たなジオメトリ情報を含む画像を表示する処理を行う表示処理手段と、を備える。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医療用画像処理装置及び医療用画像処理プログラムに関する。
【背景技術】
【０００２】
Ｘ線ＣＴ装置およびＭＲＩ装置等の各種画像診断機器は、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行うものである。例えば、心臓や血管等の循環器系およびその他の動きのある臓器の場合、それらを構成する組織の動きを断層像により観察して、それら臓器等の機能を診断することが行なわれている。このような画像診断機器により、高解像度の画像が短時間に得られるようになったことに伴い、取得したボリュームデータから、臓器や腫瘍等の関心領域（Ｒｅｇｉｏｎ ｏｆ Ｉｎｔｅｒｅｓｔ）を抽出し、これらを見やすく可視化し、あるいは定量化して面積や体積を量ることにより病変部の診断に役立てている。
【０００３】
従来例１として、特許文献１に開示されている医療用画像処理装置は、各フェイズの複数のボリュームデータから新たに一の画像を生成し、各フェイズ現れる組織の変化を一の画像に表示することができる。しかし、各フェイズ間を補うボリュームデータを生成しているものではない。ここで、フェイズ（Ｐｈａｓｅ）とは、同一の対象を同一手法で短時間のうちに撮像した複数のボリュームデータのうちの１つのボリュームデータを指す。
【０００４】
従来例２として、特許文献２に開示されている医療用画像処理装置は、各フェイズの抽出領域から構造情報（たとえば、血管の中心線）を求めて、各フェイズの異常個所を比較、検出している。そのため、少ない計算量で領域抽出を行うことができる。しかし、各フェイズ間を補うボリュームデータを生成し、各フェイズの領域抽出を行っていない。
【０００５】
従来例３として、非特許文献１では、いわゆる、モーフィングにより、医療画像とは異なる通常の画像で構成される動画像のフレーム間隔（例えば、１／６０秒）の中間位置で新たに画像を補間している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００８−３５８９５号公報
【特許文献２】特開２００７−６８６５８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、従来例１では、複数のボリュームデータから高画質な一の画像を生成しているが、複数のボリュームデータそれぞれから作成する個別の画像の画質は改善しない。これは特に動画となったときに顕著である。また、従来例２は複数のボリュームデータに含まれる組織の対応関係を取得するに留まる。また、従来例１及び従来例２で開示されている技術では、各フェイズ間を補うボリュームデータを生成しているものではないので、各フェイズの位置合わせの精度を高めるには別の方法を検討する必要がある。また、従来例３では、新たに画像を補間するタイミングは、動画像（たとえば、テレビ映像）のフレーム間隔（例えば、１／６０秒）の中間位置に固定されており、この位置を変更することはまでは示唆されていない。さらに、新たに画像を補間しても画質を改善しているわけではない。また、観察対象の領域を画像としてよりも、その形状や大きさなどのより客観的な指標で把握したい場合がある。この場合は観察対象のジオメトリ情報を作成するが、補間ボリュームを含む多数のボリュームデータそれぞれにジオメトリ情報を作成するには多大な労力がかかる。また、異なるボリュームデータのジオメトリ情報の間の対応関係を設定するのも多大な労力がかかる。
【０００８】
本発明の目的は、動き解析情報及びジオメトリ情報に基づき、複数のボリュームデータからそれぞれジオメトリ情報を作成するときの処理やユーザの労力を軽減できる医療用画像処理装置及び医療用画像処理プログラムを提供することである。
【課題を解決するための手段】
【０００９】
本発明は、時刻を有する複数のボリュームデータのそれぞれを可視化する医療用画像処理装置であって、前記複数のボリュームデータのうち、一のボリュームデータについて、ジオメトリ情報を生成する第１生成手段と、前記複数のボリュームデータのうち、前記一のボリュームデータと他の少なくとも一のボリュームデータとの間の動き解析情報を生成する第２生成手段と、前記複数のボリュームデータのうち、前記他の少なくとも一のボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する新たなジオメトリ情報を生成する第３生成手段と、前記他の少なくとも一のボリュームデータについて、前記対応する新たなジオメトリ情報を含む画像を表示するための処理を行う表示処理手段と、を備える医療用画像処理装置を提供する。
【００１０】
上記医療用画像処理装置では、前記動き解析情報に基づき、前記一のボリュームデータに含まれる特定の領域の変異又は変形を取得し、当該取得した変異又は変形、前記ジオメトリ情報、及び前記対応する新たなジオメトリ情報に基づき、ボクセル値の変化を予測し、前記一のボリュームデータのボクセル値を補正する。
【００１１】
上記医療用画像処理装置では、さらに、前記一のボリュームデータと前記他の少なくとも一のボリュームデータとの間を補間する補間ボリュームデータを生成する第４生成手段と、前記補間ボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する補間ジオメトリ情報を生成する第５生成手段とを備え、前記表示処理手段は、前記補間ボリュームデータについて、前記補間ジオメトリ情報を含む画像を表示するための処理を行う。
【００１２】
上記医療用画像処理装置では、前記一のボリュームデータの時刻情報及び前記補間ボリュームデータの時刻情報の間隔と、前記他の少なくとも一のボリュームデータ及び前記補間ボリュームデータの時刻情報の間隔とは非等間隔である。
【００１３】
上記医療用画像処理装置では、前記表示処理手段は、前記複数のボリュームデータについて描画する時に、前記ジオメトリ情報と前記新たなジオメトリ情報とを利用して、前記複数のボリュームデータで共通する特定の領域を画像として表示するための処理を行う。
【００１４】
また、本発明は、コンピュータを、上記医療用画像処理装置の各手段として機能させるための医療用画像処理プログラムを提供する。
【発明の効果】
【００１５】
本発明に係る医療用画像処理装置及び医療用画像処理プログラムによれば、動き解析情報及びジオメトリ情報に基づき、複数のボリュームデータからそれぞれジオメトリ情報を作成するときの処理やユーザの労力を軽減できる。また、客観性を増すことができる。また、ジオメトリ情報に含まれる情報を表示することによって診断に寄与できる。
【図面の簡単な説明】
【００１６】
【図１】実施の形態１に係る医療用画像処理装置１００の構成を示すブロック図
【図２】動き解析情報生成部１０３の動作を概念的に説明するための示す図
【図３】補間ボリュームデータ生成部１０５の動作を概念的に説明するための図
【図４】画像生成部１０７が生成する画像を模式的に示す図
【図５】時刻ｔｋ−１と時刻ｔｋと間の任意の時刻で画像生成部１０７が生成する画像を模式的に示す図
【図６】心室容量の大きさの経時変化を示すモデルと画像生成部１０７が作成した補間画像ｔａｕｘとを比較した図
【図７】画像生成部１０７が生成した２組の時系列に並んだ３次元画像セットを示す図
【図８】図７に示す２組みの３次元画像セットを同期させた画像のシーケンスを示す図
【図９】等間隔に複数の補間画像を実画像に補った場合の画像のシーケンスを示す図
【図１０】実画像に対して一時的に複数の補間画像補い、表示した例を示す図
【図１１】医療用画像処理装置１００の動作処理フロー
【図１２】実施の形態２に係る医療用画像処理装置３００の構成を示すブロック図
【図１３】補正処理前後の実画像の例（１）を示す図
【図１４】医療用画像処理装置３００補正処理フローの例（１）
【図１５】医療用画像処理装置３００補正処理フローの例（２）
【図１６】実画像の補正処理の例（２）を説明するための図
【図１７】実画像の補正処理の例（３）を説明するための図
【図１８】実画像及び補間画像の補正処理の例（４）を説明するための図
【図１９】実画像の補正処理の例（５）を説明するための図
【図２０】実施の形態３に係る医療用画像処理装置５００の構成を示すブロック図
【図２１】（ａ）ジオメトリ情報生成部５０６での処理１を模式的に示す図、（ｂ）ジオメトリ情報生成部５０６での処理２を模式的に示す図
【図２２】医療用画像処理装置５００の画像処理の例（１）を示す図
【図２３】医療用画像処理装置５００の画像処理の例（２）を示す図
【図２４】医療用画像処理装置５００の画像処理の例（３）を示す図
【図２５】医療用画像処理装置５００の画像処理の例（４）を示す図
【図２６】医療用画像処理装置５００の画像処理の例（５）を示す図
【発明を実施するための形態】
【００１７】
以下、本発明の実施形態について、図面を参照して説明する。
【００１８】
図１を参照して、実施の形態１に係る医療用画像処理装置について説明する。図１に示す医療用画像処理装置１００は、ボリュームデータ生成部１０１と、動き解析情報生成部１０３と、補間ボリュームデータ生成部１０５と、画像生成部１０７と、表示処理部１０９と、を備える。医療用画像処理装置１００には、医療用画像処理装置１００で生成された画像を表示する表示部１１１と、医療用画像処理装置１００の各部を操作するためのユーザインターフェースである操作部１１３と、が設けられている。
【００１９】
図１に示すように、医療用画像処理装置１００は、コンピュータ断層撮影装置に接続されている。コンピュータ断層撮影装置は、被検体の組織等を可視化するものである。Ｘ線源１からは同図に鎖線で示す縁部ビームを有するピラミッド状のＸ線ビーム束２が放射される。Ｘ線ビーム束２は、例えば患者３である被検体を透過しＸ線検出器４に照射される。Ｘ線源１及びＸ線検出器４は、本実施形態の場合にはリング状のガントリー５に互いに対向配置されている。リング状のガントリー５は、このガントリーの中心点を通るシステム軸線６に対して、同図に示されていない保持装置に回転可能（矢印ａ参照）に支持されている。
【００２０】
患者３は、Ｘ線が透過するテーブル７上に寝ている。このテーブル７は、図示されていない支持装置によりシステム軸線６に沿って移動可能（矢印ｂ参照）に支持されている。従って、Ｘ線源１及びＸ線検出器４は、システム軸線６に対して回転可能でありかつシステム軸線６に沿って患者３に対して相対的に移動可能である測定システムを構成するので、システム軸線６に関して種々の投影角及び種々の位置のもとで患者３にＸ線を投射する。その際に発生するＸ線検出器４の出力信号は、ボリュームデータ生成部１０１に供給され、ボリュームデータに変換される。
【００２１】
シーケンス走査の場合、患者３の層毎の走査が行なわれる。その際に、Ｘ線源１及びＸ線検出器４はシステム軸線６を中心に患者３の周りを回転し、Ｘ線源１及びＸ線検出器４を含む測定システムは患者３の２次元断層を走査するために多数の投影を撮影する。その際に取得された測定値から、走査された断層を表示する断層像が再構成される。相連続する断層の走査の間に、患者３はその都度システム軸線６に沿って移動される。この過程は全ての関心断層が捕捉されるまで繰り返される。
【００２２】
スパイラル走査中は、Ｘ線源１及びＸ線検出器４を含む測定システムはシステム軸線６を中心に回転し、テーブル７は連続的に矢印ｂの方向に移動する。すなわち、Ｘ線源１及びＸ線検出器４を含む測定システムは、患者３に対して相対的に連続的にスパイラル軌道上を、患者３の関心領域が全部捕捉されるまで移動する。本実施形態の場合、同図に示されたコンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部１０１に供給される。
【００２３】
ボリュームデータ生成部１０１は、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部１０３及び画像生成部１０７へ出力する。なお、ボリュームデータ生成部１０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータを保持することができる。
【００２４】
ここで、本実施の形態では、ボリュームデータ生成部１０１は、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
【００２５】
以下、説明のため、ボリュームデータ生成部１０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部１０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部１０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。なお、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
【００２６】
動き解析情報生成部１０３は、ボリュームデータ生成部１０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成し（図２参照）、補間ボリュームデータ生成部１０５に出力する。
【００２７】
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。
【００２８】
図２は、動き解析情報生成部１０３の動作を概念的に示す図である。図２に示すように、動き解析情報生成部１０３は、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを、動き解析情報として生成している。
【００２９】
補間ボリュームデータ生成部１０５は、動き解析情報生成部１０３から出力された動き解析情報に基づき、ボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、画像生成部１０７に出力する。
【００３０】
ここで、補間ボリュームデータとは、時刻ｔｋ−１と時刻ｔｋとの間の特定の時刻ｔａ
ｕｘにおけるボリュームデータである。つまり、補間ボリュームデータ生成部１０５が補間ボリュームデータｔａｕｘを生成することは、実際はコンピュータ断層撮影装置で撮像されていない時刻ｔａｕｘにおいて、あたかもコンピュータ断層撮影装置により撮像して、ボリュームデータ生成部１０１がボリュームデータを生成したことと、同等となる。
【００３１】
図３は、補間ボリュームデータ生成部１０５の動作を概念的に示す図である。図３に示すように、補間ボリュームデータ生成部１０５は、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成する。
【００３２】
ここで、本実施の形態では、図３に示すように補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない。言
い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、後述する。
【００３３】
画像生成部１０７は、ボリュームデータ生成部１０１から出力されるボリュームデータｔｋ−１、ｔｋ、及び時刻情報ｔｋ−１、ｔｋ、並びに補間ボリュームデータ生成部１０５から出力される補間ボリュームデータｔａｕｘ、及びその時刻情報ｔａｕｘを入力とする。画像生成部１０７は、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成する一方、補間ボリュームデータｔａｕｘから補間３次元画像を生成し、時系列に並べる。
【００３４】
表示処理部１０９は、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部１０７で生成された３次元画像のシーケンスを生成して、表示部１１１で動画像として再生するための処理（再生フレームレートの設定など）を行うことができる。
【００３５】
表示部１１１は、表示処理部１０９で適宜処理された時系列となった実画像ｔｋ−１、
ｔｋ及び補間画像ｔａｕｘを表示する。表示部１１１は、時系列となった実画像ｔｋ−１
、ｔｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部１１１は、画像生成部１０７で生成された３次元画像のシーケンスを、表示処理部１０９での処理の後、動画像として再生することができる。
【００３６】
操作部１１３は、表示部１１１に表示された画像を加工するための各種のインターフェースを備えており、例えば、補間画像ｔａｕｘを生成するための補間パラメータ（時刻情報ｔａｕｘなど）を設定する。そして、ユーザは、表示部１１１を見ながら、操作部１１３により、医療用画像処理装置１００の各部を操作することができる。そのため、ユーザは、インタラクティブに補間パラメータを変更できる。したがって、本実施の形態の医療用画像処理装置において、画像生成部１０７は、操作部１１３から指示に基づき、インタラクティブに補間画像ｔａｕｘを修正することができる。
【００３７】
ここで、補間パラメータとして、たとえば、（１）位置補間アルゴリズム（リニア／キュービック／スプライン）、（２）ボクセル値の補間のアルゴリズム（リニア／キュービック／スプライン）、（３）補間により生成される３Ｄ画像のピッチ（時間間隔）や数、並びに（４）４Ｄフィルタの方法として、平滑化、エッジ強調、最大値・最小値投影、差分、累積加算、ヒストグラムマッチング、およびそれらの組み合わせなどがある。
【００３８】
図４を参照して、画像生成部１０７が生成する３次元画像について説明する。図４は、画像生成部１０７が生成した各画像を模式的に示す図である。図４に示すように、実画像ｔｋ−１と実画像ｔｋとの間で、時刻情報ｔａｕｘを有する補間画像ｔａｕｘが、画像生成部１０７によって生成される。以下、説明のためボリュームデータｔｋ−１、ｔｋから生成された、時刻ｔｋ−１、ｔｋにおける実際の３次元画像を、それぞれ実画像ｔｋ−１と実画像ｔｋと称し、補間ボリュームデータｔａｕｘから生成された、時刻ｔａｕｘにお
ける３次元の補間画像を、補間画像ｔａｕｘと称する。
【００３９】
ここで、図４に示すように、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミング
に生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻
情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが非等間隔となっている。本実施の形態において、ボリュームデータｔｋ−１と、時刻情報ｔｋ−１の次のタイミングで生成されるボリュームデータｔｋとの間隔が、およそ０．１秒から数年単位である。これは、画像生成部１０７で生成される医療画像（実画像ｔｋ−１、実画像ｔｋ、及び補間画像ｔａｕｘ）とは異なる通常の画像で構成される動画像（たとえば、テレビなどの映像）のフレーム間隔（例えば、１／６０秒）と比べて非常に大きいことに起因する。
【００４０】
つまり、本実施の形態では、実画像ｔｋ−１、と実画像ｔｋとの間の間隔（Ｃ＋Ｄ）が、例えば年単位であって、実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の関心領域の経年変化を観察するのに、その中間位置での補間画像では十分ではないからである。そのため、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが非等間隔とし、実画像ｔｋ−１、と実画像ｔｋとの間で被検体である患者３の検査部位（関心領域）の経年変化を中間位置とは異なる、実画像ｔｋ−１、と実画像ｔｋとの間の様々な位置で補間画像を生成することを可能としている。また、撮影に要する放射線量を最小にする目的で、関心の最も高い時刻では撮影の時間間隔を密にし、それ以外では疎に撮影することが行われている。また、心電図など被検体のリアルタイム情報を取得する別途の機器で指定されるタイミングで撮影することもある。これらによって撮影間隔が非等間隔となりうる。リアルタイム情報には他には例えば、造影剤の注入のタイミングもある。
【００４１】
上述のように、本実施の形態では、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイ
ミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘ
の時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが、非等間隔であるが、これに限らない。実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の関心領域の経時変化を観察するのに適したタイミングであれば、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが、等間隔に設定されてもよい。すわなち、実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の検査部位（関心領域）の経時変化を観察するのに適したタイミングであれば、時刻情報ｔｋ−１と時刻情報ｔｋ−１との間の任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成しても良い。言い換えると、本実施の形態において、ボリュームデータｔａｕｘの時刻ｔａｕｘは、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。これにより、本実施の形態では、Ｘ線の線量を削減し、被検体である患者３への影響を極力抑えることができる。
【００４２】
（任意の時刻の補間画像）
図５は、時刻ｔｋ−１と時刻ｔｋ−１との任意の時刻で画像生成部１０７が補間画像ｔ
ａｕｘを生成する一例を模式的に示す図である。図５に示すように、画像生成部１０７は、実画像ｔｋ−１を生成した時刻０．０ｓｅｃと、実画像ｔｋを生成した時刻情報１．５ｓｅｃとの間の任意の時刻情報１．０ｓｅｃの補間画像ｔａｕｘを生成する。そのため、実画像を生成した時刻の間の任意の時刻で補間画像を生成することができ、実画像ｔｋ−１、と実画像ｔｋとの間で発生している患者３の診断範囲の経時変化を観察することができる。
【００４３】
（補間画像とモデルとの比較）
図６は、心室容量の大きさの経時変化を示すモデルと画像生成部１０７が作成した補間画像ｔａｕｘとを比較した図である。上述のように、時刻情報ｔｋ−１と時刻情報ｔｋ−１との間の時刻であれば、任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成することが可能なので、たとえば図６に示すように、心室容量の大きさの経時変化を示すモデルの最大容量を示す時刻で、画像生成部１０７が補間画像ｔａｕｘを作成し、心室容量の大きさの経時変化を示すモデルと補間画像ｔａｕｘにあらわれる心室の大きさとを比較することができる。
【００４４】
心室容量の大きさの経時変化を示すモデルは、図６に示すように線形的に変化するものではないので、２つの実画像ｔｋ−１、ｔｋ間の中間位置では、心室容量は最大とならない、したがって、本実施の形態では、モデル上、心室容量が最大となる特徴を有する時刻で、実画像ｔｋ−１、ｔｋから補間画像ｔａｕｘを生成する。図６に示すように、モデル上、心室容量が最大となる時刻はｔ＝０．７８ｓｅｃである。したがって、画像生成部１０７は、時刻情報ｔａｕｘをモデル上、心室容量が最大となる時刻はｔ＝０．７８ｓｅｃとし、補間画像ｔａｕｘを生成する。
【００４５】
上述のように、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴が現れると推定される時刻を補間画像の時刻情報とすることで、画像生成部１０７は、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴を画像として生成し、表示部１１１が可視化することができる。
【００４６】
なお、心室容量の大きさの経時変化を示すモデルから、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴が現れると推定しているが、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴を推定する手段は当該モデルだけに限らない。例えば、被検体である患者３に対して造影剤などが注入された外的要因が発生した時刻を、画像上に特徴が現れる時刻であると推定することも可能できる。そのため、補間画像ｔａｕｘを作成し、実画像ｔｋ−１、ｔｋ、及び補間画像ｔａｕｘの一連の画像の変化から、被検体（患者３）の検査部位（関心領域）の外的要因による変化を観察することができる。また、時刻情報ｔａｕｘモデルのみならず外部機器からの入力によって定めることが出来る。これは心電図など被検体のリアルタイム情報を実画像の取得と同時に行い、これを参照することによって求めることが出来る。例えば、心電図上の波形の特徴的な時刻や、造影剤の注入を行ったタイミングなどが考えられる。
【００４７】
（異なる画像セットの同期）
また、画像生成部１０７は、同一の被検体から同一の検査部位を対象として、異なるモードで動作するコンピュータ断層撮影装置又は異なるコンピュータ断層撮影装置により供給された複数のボリュームデータから、時系列の複数の３次元画像セットを生成し、複数の３次元画像セットから正確に同期した画像のシーケンスを生成することができる。図７は、画像生成部１０７が生成した２組の時系列に並んだ３次元画像セットを示す。また、図８は、図７に示す２組みの３次元画像セットを同期させた画像のシーケンスを示す図である。
【００４８】
図７に示す３次元画像セットＢでは、第１の時刻情報ｔｋ−１（＝０．２ｓｅｃ）を有する第１実画像ｔｋ−１、及び第２の時刻情報ｔｋ（＝１．７ｓｅｃ）を有する第２実画像ｔｋ、並びに、第１実画像ｔｋ−１及び第２実画像ｔｋから生成され、時刻情報ｔａｕｘ２を有する第２補間画像ｔａｕｘ２が時系列に並んでいる。また、図７に示す３次元画像セットＡでは、第３の時刻情報ｔｋ−１（＝０．０ｓｅｃ）を有する第３実画像ｔｋ−１、及び第４の時刻情報ｔｋ（＝１．５ｓｅｃ）を有する第４実画像ｔｋ、並びに第３実画像ｔｋ−１及び第４実画像ｔｋから生成され、時刻情報ｔａｕｘ１を有する第１補間画像ｔａｕｘ１が時系列に並んでいる。
【００４９】
そして、図８に示すように、画像生成部１０７は、各画像の保持する時刻情報に基づき、３次元画像セットＢの第２補間画像ｔａｕｘ２を有する時刻情報ｔａｕｘ２を、３次元画像セットＡに時間的に同期させて、新たな画像のシーケンスを生成する。
【００５０】
したがって、本実施の形態では、画像生成部１０７は、同一の被検体から同一の検査部位を対象として、異なるモードで動作するコンピュータ断層撮影装置又は異なるコンピュータ断層撮影装置により供給された複数のボリュームデータから、時系列の複数の３次元画像セットを生成することができる。そして、表示処理部１０９は、複数の３次元画像セットから正確に時間的に同期した画像のシーケンスを生成することができる。
【００５１】
また、本実施の形態では、画像生成部１０７は、図８に示す３次元画像のシーケンスを用いれば、観察対象部位（心臓）のマルチフェイズデータでの同一フェイズ間のずれの補正も可能となる。また、本実施の形態では、画像生成部１０７は、図８に示す３次元画像のシーケンスを用いれば、「コンピュータ断層撮影装置から取得した検査の画像取得時間情報」と「表示処理部１０９で生成される再生フレームレート」を元に、表示部１１１で正確な実時間で動画再生することができる。
【００５２】
さらに、本実施の形態では、画像生成部１０７は、図７に示す一方の３次元画像セットＢにだけ存在する特定の取得時刻（例えば、ｔｋ＝１．７ｓｅｃ）の実画像ｔｋに合わせて、他方の３次元画像セットＡの第３実画像ｔｋ−１及び第４実画像ｔｋから新たな補間画像を作成することにより、より正確な３次元画像間の比較を可能にする。
【００５３】
図９に、画像生成部１０７が時間的に等間隔に複数の補間画像を実画像に補った場合の画像のシーケンスを示す。また、画像生成部１０７は、図９に示すように、同一の被検体から同一の検査部位を対象として、撮影間隔の異なるコンピュータ断層撮影装置により供給された複数のボリュームデータｔ１（＝０．００ｓｅｃ）、ｔ２（＝０．１７ｓｅｃ）、ｔ３（＝０．３０ｓｅｃ）から生成された複数の実画像ｔ１、ｔ２、ｔ３のうち実画像ｔ１、ｔ３に対して、時間的に等間隔に補う、時刻情報ｔ４（＝０．１０ｓｅｃ）を有する補間画像ｔａｕｘ４、及び時刻情報ｔ５（＝０．２０ｓｅｃ）を有する補間画像ｔａｕｘ５を生成することができる。そのため、本実施の形態では、時間的に等間隔でない実画像に対して、補間画像を等時間間隔に生成するので、表示部１１１が表示する医療用画像の表示や解析の質・能力を高めることができる。
【００５４】
図１０に、実画像に対して一時的に複数の補間画像補い、表示部１１１に表示した例を示す。画像生成部１０７は、図１０に示すように、表示部１１１に３次元画像のシーケンスを表示している時のみ、一時的に、画像生成部１０７が複数の補間画像ｔｅｍｐ１から補間画像ｔｅｍｐ３を生成し、各フレームごとに、レンダリング後に複数の補間画像ｔｅｍｐ１から補間画像ｔｅｍｐ３を破棄することで、メモリの使用量を削減でき、本実施の形態の医療用画像処理装置１００への負荷を軽減することができる。
【００５５】
次に、図１１を参照して、本実施の形態に係る医療用画像処理装置１００の動作処理フローを説明する。図１１は、医療用画像処理装置１００の動作処理フローである。
【００５６】
まず、ステップＳＴ１で、ボリュームデータ生成部１０１は、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータ、及び各ボリュームデータの時刻情報ｔｎを生成する。そして、ステップＳＴ２へ遷移する。
【００５７】
次に、ステップＳＴ２で、動き解析情報生成部１０３は、ボリュームデータ生成部１０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。そして、ステップＳＴ３へ遷移する。
【００５８】
次に、ステップＳＴ３で、補間ボリュームデータ生成部１０５は、動き解析情報生成部１０３から出力された動き解析情報に基づき、ボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成する。そして、ステップＳＴ４へ遷移する。
【００５９】
次に、ステップＳＴ４で、画像生成部１０７は、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成する一方、補間ボリュームデータｔａｕｘから補間３次元画像を生成し、時系列に並べる。そして、ステップＳＴ５へ遷移する。
【００６０】
次に、ステップＳＴ５で、表示処理部１０９は、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う。そして、ステップＳＴ６へ遷移する。
【００６１】
最後に、ステップＳＴ６で、表示部１１１は、表示処理部１０９で適宜処理された時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。そして、動作処理が終了する。
【００６２】
以上、本実施の形態に係る医療用画像処理装置１００によれば、動き解析情報に基づき補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示することができる。したがって、本実施の形態に係る医療用画像処理装置１００によれば、被検体である患者３の検査部位（関心領域）の観察を精度良く行うことができる。
【００６３】
例えば、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる。また、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。
【００６４】
さらに、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。さらに、また、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、セグメンテーションにより隣接する関心領域（臓器）を設定し、動き解析情報により複数の臓器の動きが同一の部分を推定、臓器間癒着の可能性あり、として表示部１１１に表示することができる。
【００６５】
また、本実施の形態に係る医療用画像処理装置１００によれば、動き解析情報に基づき補間画像を実画像間に補間すれば、表示部１１１のノイズ軽減、ノイズのばらつき調整が可能である。
【００６６】
例えば、医療用画像処理装置１００が生成した動き解析情報によって、実画像及び補間画像の対応画素を特定することにより、本来の詳細な画像情報を失うことなくノイズを軽減できる。また、医療用画像処理装置１００により、ＣＴマルチフェイズ検査において、Ｘ線照射量の多いフェイズと少ないフェイズが混在しても、Ｘ線照射量の少ないフェイズについては、前後のフェイズの画素値を利用することで、良好な補間画像を得られる。
【００６７】
さらに、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、重要なフェイズのみＸ線照射量を多くして実画像を生成し、それ以外のフェイズについては、補間画像を生成することで、患者３のＸ線の被ばく量を低減することができる。
【００６８】
なお、本実施の形態に係る医療用画像処理装置１００では、図１に示すコンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部１０１に供給されているが、これに限らない。ボリュームデータ生成部１０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
【００６９】
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。
【００７０】
なお、本実施の形態において、表示処理部１０９が、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部１１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
【００７１】
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻ｔに基づき、生成しているが、これに限らない。
【００７２】
なお、本実施の形態の説明に用いた各機能ブロック（図１参照）は、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。
【００７３】
なお、本実施の形態では、心臓の観察に用いる例を示したが、これに限らない。例えば、腫瘍組織の経時変化の観察、血流灌流の観察などに用いることも出来る。
【００７４】
（実施の形態２）
上述した実施の形態１に係る医療用画像処理装置１００によれば、動き解析情報に基づき、補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示する。しかし、補間画像を実画像間に補間する場合、実画像の画質によっては、単純に補間画像を補間するだけでは、画像内の特定の領域の変化を正確に表示できない場合もある。
【００７５】
そこで、実施の形態２に係る医療用画像処理装置３００では、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正して、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
【００７６】
実施の形態２に係る医療用画像処理装置３００では、必ずしも補間画像を生成する必要はない。動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質を補正して、２以上の実画像内の特定の領域の変化を正確に表示することができれば、補間画像を生成する必要はない。
【００７７】
図１２を参照して、本実施の形態２に係る医療用画像処理装置３００の構成について説明する。図１２に示す医療用画像処理装置３００は、ボリュームデータ生成部３０１と、動き解析情報生成部３０３と、補間ボリュームデータ生成部３０５と、補正処理部３０６と、画像生成部３０７と、表示処理部３０９と、を備える。さらに、医療用画像処理装置３００には、医療用画像処理装置３００で生成された画像を表示する表示部３１１と、医療用画像処理装置１００の各部を操作するためのユーザインターフェースである操作部３１３と、が設けられている。
【００７８】
ボリュームデータ生成部３０１は、実施の形態１と同様、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部３０３と、補正処理部３０６と、画像生成部３０７とへ出力する。
【００７９】
また、ボリュームデータ生成部３０１は、コンピュータ断層撮影装置から、装置の撮像条件（例えば、Ｘ線量、撮影時間、造影剤の注入の時刻）を取得し、画像生成部３０７へ出力する。撮像条件により、のちに生成される実画像の画質が大きく変化する。なお、コンピュータ断層撮影装置の撮像条件以外にも、コンピュータ断層撮影装置の撮影対象の条件（患者の動きや、患者の体内に含まれる金属など）により、のちに生成される実画像の画質が大きく変化する。
【００８０】
なお、ボリュームデータ生成部３０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータ、コンピュータ断層撮影装置の撮像条件、コンピュータ断層撮影装置の撮影対象の条件を保持することができる。
【００８１】
ここで、本実施の形態では、ボリュームデータ生成部３０１は、実施の形態１と同様、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
【００８２】
以下、説明のため、ボリュームデータ生成部３０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部３０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部３０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。また、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
【００８３】
動き解析情報生成部３０３は、ボリュームデータ生成部３０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。
【００８４】
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。
【００８５】
そして、動き解析情報生成部３０３は、生成した動き解析情報を、補間ボリュームデータ生成部３０５と画像生成部３０７とに出力する。動き解析情報生成部３０３が生成する動き解析情報は、実施の形態１の動き解析情報生成部１０３と同様（図２参照）、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを示す。
【００８６】
補間ボリュームデータ生成部３０５は、操作部３１３の入力指示に基づき、必要に応じて、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、補正処理部３０６に出力する。
【００８７】
ここで、本実施の形態では、実施の形態１と同様、補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない。言い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、実施の形態１で上述した理由と同じであるため、詳細な説明は省略する。
【００８８】
本実施の形態に係る医療用画像処理装置３００の一つの特徴として、補正処理部３０６は、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像（入力画像）の画質、又は２つの実画像間を補間する補間画像の画質を補正する。
【００８９】
＜補正処理１：画質の均一化＞
補正処理の一例として、コンピュータ断層撮影装置の撮影対象の条件などから、補正処理部３０６の補正処理前の時刻ｔｋ−１、ｔｋにおける画像の画質が均一でない場合、補正処理部３０６は、動き解析情報生成部３０３からの入力された動き解析情報、ボリュームデータ生成部３０１から入力されたボリュームデータｔｋのボクセル及びボリュームデータｔｋ−１のボクセルに基づき、画像生成部３０７で生成する実画像ｔｋ−１、ｔｋ、の画素値を求める。詳細は、図１３、図１４を参照して後述する。
【００９０】
図１３に、補正処理前後の実画像ｔｋ−１、ｔｋの例（１）を示す。なお、補正処理前の実画像（入力画像）は、説明のため示しているが、画像生成部３０７が実際に生成する画像（出力画像）は、補正処理部３０６で補正処理が行われた後の実画像ｔｋ−１、ｔｋである。図１３に示すように、補正処理前の時刻ｔｋ−１、ｔｋにおける画像の画質が均一ではないが、補正処理部３０６で補正処理が行われた後は、実画像ｔｋ−１、ｔｋの画質は均一となる。そのため、本実施の形態に係る医療用画像処理装置３００は、動き解析情報及び画像の画質情報に基づき、複数のボリュームデータからそれぞれ画像を作成するときの画質を向上させることができる。さらに、本実施の形態に係る医療用画像処理装置３００は、補正処理後の各画像の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の経時間変化を、正確に表示することができる。なお、必要に応じて、補正処理部３０６は、補間ボリュームデータ生成部３０５から出力された補間ボリュームデータｔａｕｘを、実画像ｔｋ−１、ｔｋと同様に補正処理を行う。
【００９１】
次に、図１４を参照して、実画像ｔｋ−１、ｔｋの画質を均一とするための補正処理について説明する。図１４は、本実施の形態に係る医療用画像処理装置３００補正処理フローの例（１）である。
【００９２】
図１４に示すように、ステップＳＴ３１では、出力フェイズの前後に相当する、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１を取得する。そして、ステップＳＴ３２へ遷移する。
【００９３】
ステップＳＴ３２では、ボリュームデータｔｋ、ｔｋ−１それぞれについて、ボクセル位置ごとに、信頼度（０〜１）を決定する。ここで、ボリュームデータｔｋ、ｔｋ−１についての信頼度（信頼性の高さ）は、出力フェイズとの時間的距離、ボリュームデータｔｋ、ｔｋ−１を撮影した時の撮影条件（Ｘ線照射量など）、当該ボリュームデータｔｋ、ｔｋ−１のＳ／Ｎ比、近傍の領域の移動速度や画素値の特徴などにより、決定することが考えられる。そして、ステップＳＴ３３へ遷移する。
【００９４】
ステップＳＴ３３では、出力フェイズの全ボクセルについて、ボリュームデータｔｋ、ｔｋ−１の画素値を、動き解析情報及びボリュームデータｔｋ、ｔｋ−１より求める。
そして、ステップＳＴ３４へ遷移する。
ステップＳＴ３４では、出力フェイズの全ボクセルについて、前後の入力フェイズでの対応する位置を、動き解析情報により求める。そして、ステップＳＴ３５へ遷移する。
【００９５】
ステップＳＴ３５では、各対応位置の画素値を入力画像上でサンプリングし、サンプリングした値を、入力フェイズに対応する入力画素値とする。そして、ステップＳＴ３６へ遷移する。
ステップＳＴ３６では、入力フェイズに対応する入力画素値の信頼度を取得する。
そして、ステップＳＴ３７へ遷移する。
【００９６】
ステップＳＴ３７では、入力フェイズに対応する入力画素値に、取得した信頼度を乗じ、当該入力画素値について、各入力フェイズで加算したものを、出力フェイズで、補正後のボリュームデータｔｋ、ｔｋ−１の画素値とする。そして、処理が終了する。
【００９７】
ここで、ステップＳＴ３２における、ボリュームデータｔｋ、ｔｋ−１の信頼度（信頼性の高さ）について、画素値の信頼性の高さに応じ、４Ｄフィルタのパラメータを各実画像ｔｋ、ｔｋ−１毎、又はボクセル毎に変えることで、全体の適正なＳ／Ｎを保ちつつ、元画像の空間分解能を最大限活用することができる。
【００９８】
画素値の信頼性は、時刻、ボクセルの位置、画素値、収縮度、速度、加速度、近傍の画素の変位から求めた回転量、変形の度合いなどの関数である。また、近傍の領域の特徴（高画素値の領域の有無、など）が、画素値の信頼性に影響することもある。
【００９９】
４Ｄフィルタの方法として、平滑化、エッジ強調、最大値・最小値投影、差分、累積加算、ヒストグラムマッチング、およびそれらの組み合わせなどがある。
【０１００】
＜補正処理２：画像の特性の均一＞
また、補正処理部３０６は、動き解析情報生成部３０３からの入力された動き解析情報、並びにボリュームデータ生成部３０１から入力されたボリュームデータｔｋ、のボクセル及びボリュームデータｔｋ−１のボクセルに基づき、画像生成部３０７で生成する実画像ｔｋ−１、ｔｋの画素値の値的な歪み量を求める。詳細は、図１５を参照して後述する。
【０１０１】
なお、補正処理部３０６が画像生成部３０７で生成する実画像ｔｋ−１、ｔｋの画素値の値的な歪み量を求めるのは、たとえば、コンピュータ断層撮影装置の撮影対象の条件(画質情報の元となる条件)などから、生成される画像の特性（明度、ダイナミックレンジ）が均一でない場合が考えられる。
【０１０２】
ここで、画素値の値的な歪み量とは、対応時刻、変位後の位置、対応位置の入力画素値、変位場の収縮度、変位の大きさ（速度）やその加速度、近傍での最大画素値や、その最大画素値のボクセルから当該ボクセルまでの距離、などの関数である。
【０１０３】
ボリュームデータ生成部３０１から入力されたボリュームデータｔｋ、のボクセル及びボリュームデータｔｋ−１のボクセルと、動き解析情報から画素値の歪の量の関数を求める方法としては、対応する位置の画素値の変化量を画像内で統計処理する、などが考えられる。
【０１０４】
そして、補正処理部３０６で求めた画素値の値的な歪み量に基づき、実画像ｔｋ−１、ｔｋを補正処理すると、生成される画像の特性（明度、ダイナミックレンジ）が均一な画像を作成することができる。
【０１０５】
＜補正処理フロー（２）：画像の特性（明度、ダイナミックレンジ）の均一化＞
ここで、図１５を参照して、実画像ｔｋ−１、ｔｋの画質、における特性（明度、ダイナミックレンジ）を均一とするための、補正処理部３０６で補正処理について説明する。図１５は、本実施の形態に係る医療用画像処理装置３００補正処理フローの例（２）である。
【０１０６】
図１５に示すように、ステップＳＴ４１では、出力フェイズの前後に相当する、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１を取得する。そして、ステップＳＴ４２へ遷移する。
【０１０７】
ステップＳＴ４２では、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１と動き解析情報から、画素値の値的な歪の量の関数を算出する。そして、ステップＳＴ４３へ遷移する。
ステップＳＴ４３では、歪の量の関数から、ボリュームデータｔｋ、ｔｋ−１の各ボクセル位置について、画素値の補正係数を算出する。そして、ステップＳＴ４４へ遷移する。
【０１０８】
ステップＳＴ４４では、出力フェイズの全ボクセルについて、動き解析情報及びボリュームデータｔｋ、ｔｋ−１から、ボリュームデータｔｋ、ｔｋ−１の画素値を算出する。
そして、ステップＳＴ４５へ遷移する。
ステップＳＴ４５では、出力フェイズの全ボクセルについて、前後の入力フェイズでの対応する位置を、動き解析情報から算出する。そして、ステップＳＴ４６へ遷移する。
【０１０９】
ステップＳＴ４６では、各対応位置の画素値を入力画像上でサンプリングし、サンプリングした値を、入力フェイズに対応する入力画素値とする。そして、ステップＳＴ４７へ遷移する。
ステップＳＴ４７では、上記入力画素値に、対応ボクセル位置の補正係数による変換をし、これをそれぞれの入力フェイズについて加算したものを、出力フェイズの画素値とする。そして、処理が終了する。
【０１１０】
＜補正処理３：変形又は変異による画素値の変化＞
また、本実施の形態に係る医療用画像処理装置３００の一つの特徴として、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報に基づき、動き解析情報から得られる入力画像内の特定の領域に関する変形又変位による画素値の変化を見込んだ補正処理を行うこともできる。図１６を参照して、説明する。
【０１１１】
図１６は、実画像ｔｋ−１、ｔｋの補正処理の例（２）を説明するための図である。説明のため、図１６では、入力画像ｔｋ−１、ｔｋ内の特定の領域を領域Ｒ１とする。図１６に示すように、領域Ｒ１は、実画像ｔｋ−１から実画像ｔｋへ遷移する間に縮小する。領域Ｒが縮小すると、領域Ｒ１では造影剤などが濃縮され、何も補正処理を実行しなければ、実画像ｔｋ上、領域Ｒ１の画素値が大きくなってしまう。これは、造影剤の注入の時刻情報を加味すると更に効果的である。
【０１１２】
そこで、補正処理部３０６は、動き解析情報生成部３０３から得られる動き解析情報に基づき、領域Ｒ１の変形変異を取得し、実画像ｔｋに対して変形変異を見込んだ補正処理を行う。
【０１１３】
上述のように、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報などに基づき、動き解析情報から得られる入力画像内の特定の領域Ｒ１に関する変形又変位による、画素値の変化を見込んだ補正処理を行うことで、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
【０１１４】
＜補正処理４：代表値での置き換え＞
図１７は、実画像ｔｋ−１、ｔｋ、ｔｋ＋１の補正処理の例（３）を説明するための図
である。説明のため、図１７では、実画像（入力画像）ｔｋ−１、ｔｋ、ｔｋ＋１内の特
定の領域を領域Ｒ２とする。図１７に示すように、領域Ｒ２は、実画像ｔｋ−１から実画像ｔｋ＋１へ遷移する間に縮小する。領域Ｒが縮小すると、領域Ｒ２では造影剤などが濃縮され、何も補正処理を実行しなければ、実画像（入力画像）ｔｋ−１、よりも、実画像ｔｋ、ｔｋ＋１の領域Ｒ２の画素値が大きくなってしまう。
【０１１５】
そこで、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋ間、及び実画像ｔｋ、ｔｋ＋１間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、実画像ｔｋについて実画像ｔｋ−１、ｔｋ、ｔｋ＋１で対応する画素値を比較し
、もっともＳ／Ｎの高い画素の画素値を、実画像ｔｋにおける画素値とする。
【０１１６】
また、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋ間、及び実画像ｔｋ、ｔｋ＋１間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、実画像ｔｋに対して、実画像ｔｋ−１、ｔｋ、ｔｋ＋１で対応する画素値で、Ｓ／Ｎを加味した加重平均を取り、それを新しい画素値とすることも可能である。
【０１１７】
上述のように、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報などに基づき、動き解析情報から得られる入力画像内の特定の領域Ｒ２に関する変形又変位による、画素値の変化を見込んだ補正処理を行うことで、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
【０１１８】
＜補正処理５：２以上実画像の画質の均一化＞
図１８は、実画像及び補間画像の補正処理の例（４）を説明するための図である。補正処理の他の例として、図１８に示すように、上述した補正処理を組み合わせて、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ−１、ｔｋ間の補間画像tａｕｘ２について、実画像ｔｋ−１、ｔｋに関する動き解析情報、及び画質を均一化する補正処理（図１３、図１４参照）後の実画像ｔｋ−１、ｔｋの画質情報などに基づき、補間画像tａｕｘ２の補正処理を行うことができる。
【０１１９】
補間画像tａｕｘ２の補正処理を行う場合、たとえば、図１６を参照して説明した方法
を用いると、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋに関する動き解析情報に基づき、図１８に示す領域Ｒ３の変形変異を取得し、補間画像tａｕｘ２に対して変形変異を見込んだ補正処理を行う。
【０１２０】
また、補間画像tａｕｘ２の補正処理を行う場合、たとえば、図１７を参照して説明し
た方法を用いると、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、補間画像ｔａｕｘ２間、及び補間像ｔａｕｘ２、実画像ｔｋ間の動き解析情報
に基づき、領域Ｒ２の変形変異を取得し、補間実画像ｔａｕｘ２について実画像ｔｋ−１
、ｔｋ、及び補間画像ｔａｕｘ２で対応する画素値を比較し、もっともＳ／Ｎの高い画素
の画素値を、実画像ｔａｕｘ２における画素値とする。この場合、比較する補間画像ｔａ
ｕｘ２の画素値をゼロとすれば良い。
【０１２１】
＜補正処理６：２番目以降の画像の画素値の差分の保存＞
図１９は、実画像の補正処理の例（５）を説明するための図である。補正処理の他の例として、図１９に示すように、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ−１、ｔｋについて、実画像ｔｋ−１、ｔｋに関する動き解析情報、及び実画像ｔｋ−１、ｔｋの画質情報などに基づき、実画像ｔｋの補正処理を行い、実画像ｔｋ−１から数えて２番目の実画像ｔｋについては、実画像ｔｋ−１の画素値との差分のみをデータとして保存しておくで、情報量（量子化ビット数）を減らしつつ、画質を維持することができる。同様に、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ、ｔｋ＋１について、実画像ｔｋ、ｔｋ＋１に関する動き解析情報、及び実画像ｔｋ、ｔｋ＋１の画質情報などに基づき、実画像ｔｋ＋１の補正処理を行い、実画像ｔｋ＋１については、実画像ｔｋの画素値との差分のみをデータとして保存しておくで、情報量（量子化ビット数）を減らしつつ、画質を維持することができる。
【０１２２】
画像生成部３０７は、補正処理部３０６から出力される補正処理後のボリュームデータｔｋ−１、ｔｋ、とその時刻情報ｔｋ−１、ｔｋ、補正処理後の補間ボリュームデータｔａｕｘ、及びその時刻情報ｔａｕｘを入力とする。そして、画像生成部３０７は、補正処理後の、補間ボリュームデータｔａｕｘ、各ボリュームデータｔｋ−１、ｔｋから実際の
３次元画像を生成し、時系列に並べる。
【０１２３】
表示処理部３０９は、画像生成部３０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部３０７で生成された３次元画像のシーケンスを生成して、表示部３１１で動画像として再生するための処理（再生フレームレートの設定など）を行う ことができる。
【０１２４】
表示部３１１は、表示処理部３０９で適宜処理された、時系列の実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。表示部１１１は、時系列となった実画像ｔｋ−１、ｔ
ｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部１１１は、画像生成部３０７で生成された３次元画像のシーケンスを、表示処理部３０９での処理の後、動画像として再生することができる。
【０１２５】
以上、実施の形態２に係る医療用画像処理装置３００によれば、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正して、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。したがって、本実施の形態に係る医療用画像処理装置３００によれば、被検体である患者３の検査部位（関心領域）の観察を精度良く行うことができる。
【０１２６】
例えば、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる。また、医療用画像処理装置３００により、動き解析情報に基づき補間画像を実画像間に補間すれば、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。
【０１２７】
さらに、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。さらに、また、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、複数のボリュームデータに対してセグメンテーション処理を行うときに処理結果を安定させることが出来る。また、臓器間癒着の判断を自動処理で行う場合も医師等が動画を直接観察する場合も容易になる。
【０１２８】
また、本実施の形態に係る医療用画像処理装置３００によれば、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、表示部３１１のノイズ軽減、ノイズのばらつき調整が可能である。
【０１２９】
例えば、医療用画像処理装置３００が生成した動き解析情報によって、実画像又は補間画像の対応画素を特定することにより、本来の詳細な画像情報を失うことなくノイズを軽減できる。また、医療用画像処理装置３００により、ＣＴマルチフェイズ検査において、Ｘ線照射量の多いフェイズと少ないフェイズが混在しても、Ｘ線照射量の少ないフェイズについては、前後のフェイズの画素値を利用することで、良好な補間画像を得られる。
【０１３０】
さらに、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、重要なフェイズのみＸ線照射量を多くして実画像を生成し、それ以外のフェイズについては、補間画像を生成することで、患者３のＸ線の被ばく量を低減することができる。
【０１３１】
なお、本実施の形態に係る医療用画像処理装置３００では、コンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部３０１に供給されているが、これに限らない。ボリュームデータ生成部３０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
【０１３２】
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。
【０１３３】
なお、本実施の形態において、表示処理部３０９が、画像生成部３０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部３１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
【０１３４】
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻に基づき、生成しているが、これに限らない。
【０１３５】
なお、本実施の形態の説明に用いた各機能ブロック（図１２参照）は、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。
【０１３６】
なお、本実施の形態では、心臓の観察に用いる例を示したが、これに限らない。例えば、腫瘍組織の経時変化の観察、血流灌流の観察などに用いることも出来る。
【０１３７】
（実施の形態３）
上述した実施の形態１に係る医療用画像処理装置１００によれば、動き解析情報に基づき、補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示する。その一方で、観察対象の領域を画像としてよりも、その形状や大きさなどのより客観的な指標で把握したい場合がある。この場合は観察対象のジオメトリ情報を作成するが、補間ボリュームを含む多数のボリュームデータそれぞれにジオメトリ情報を作成するには多大な労力がかかる。また、異なるボリュームデータのジオメトリ情報の間の対応関係を設定するのも多大な労力がかかる。
【０１３８】
そこで、実施の形態３に係る医療用画像処理装置５００では、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成することで、個別にジオメトリ情報を生成するのに比較して、セグメンテーションに係る操作の手間、補正処理を減らすことができる。さらに、実施の形態３に係る医療用画像処理装置５００では、実画像と他の実画像との画質に差異によって、セグメンテーションの結果の揺らぎを抑止できる。そのため、ユーザは、より客観的な画像解析を行うことができる。
【０１３９】
さらには、実施の形態３に係る医療用画像処理装置５００では、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報から他の実画像のジオメトリ情報を生成しているので、他の実画像のジオメトリ情報は、内在的に実画像のボリュームデータから生成したジオメトリ情報との対応関係を有している。そのため、他の実画像のジオメトリ情報と実画像のジオメトリ情報との対応関係を改めて取得、計算する必要がない。
【０１４０】
図２０を参照して、本実施の形態３に係る医療用画像処理装置５００の構成について説明する。図２０に示す医療用画像処理装置５００は、ボリュームデータ生成部５０１と、動き解析情報生成部５０３と、補間ボリュームデータ生成部５０５と、ジオメトリ情報生成部５０６と、画像生成部５０７と、表示処理部５０９と、を備える。さらに、医療用画像処理装置５００には、医療用画像処理装置５００で生成された画像を表示する表示部５１１と、医療用画像処理装置５００の各部を操作するためのユーザインターフェースである操作部５１３と、が設けられている。
【０１４１】
ボリュームデータ生成部５０１は、実施の形態１と同様、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部５０３と、ジオメトリ情報生成部５０６と、画像生成部５０７とへ出力する。
【０１４２】
また、ボリュームデータ生成部５０１は、コンピュータ断層撮影装置から、装置の撮像条件（例えば、Ｘ線量、撮影時間、造影剤の注入の時刻）を取得し、画像生成部５０７へ出力する。なお、コンピュータ断層撮影装置の撮像条件以外にも、コンピュータ断層撮影装置の撮影対象の条件（患者の動きや、患者の体内に含まれる金属など）により、のちに生成される実画像の画質が大きく変化する。
【０１４３】
なお、ボリュームデータ生成部５０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータ、コンピュータ断層撮影装置の撮像条件、コンピュータ断層撮影装置の撮影対象の条件を保持することができる。
【０１４４】
ここで、本実施の形態では、ボリュームデータ生成部５０１は、実施の形態１と同様、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
【０１４５】
以下、説明のため、ボリュームデータ生成部５０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部５０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部５０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。また、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
【０１４６】
動き解析情報生成部５０３は、ボリュームデータ生成部５０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。
【０１４７】
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。
【０１４８】
そして、動き解析情報生成部５０３は、生成した動き解析情報を、補間ボリュームデータ生成部５０５と、ジオメトリ情報生成部５０６と、画像生成部５０７とに出力する。動き解析情報生成部５０３が生成する動き解析情報は、実施の形態１の動き解析情報生成部１０３と同様（図２参照）、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを示す。
【０１４９】
補間ボリュームデータ生成部５０５は、操作部５１３の入力指示に基づき、必要に応じて、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、画像生成部５０７とに出力する。
【０１５０】
ジオメトリ情報生成部５０６は、ボリュームデータ生成部５０１から出力された、撮影時刻の異なる複数のボリュームデータとその時刻情報に基づき、必要に応じて、１又は複数のボリュームデータのジオメトリ情報を生成し、画像生成部５０７に出力する。以下、時刻情報ｔｋにおけるボリュームデータｔｋのジオメトリ情報を、ジオメトリ情報Ｇｔｋと称する。
【０１５１】
ここで、ジオメトリ情報とは、１つの３次元実画像について、中実、表面、線、位置、又はこれらの要素の組み合わせで表現した組織の情報である。特に、ボクセルデータのマスクで表現される領域の中実や表面であることが多く、セグメンテーションなどの従来手法を用いることができる。また、トラッキングなどの従来手法を用いて取得される線であることもある。これらはマスク・パス・点として表現することができる。
【０１５２】
図２１（ａ）、図２１（ｂ）を参照して、ジオメトリ情報生成部５０６での処理について説明する。図２１（ａ）は、ジオメトリ情報生成部５０６での処理１を模式的に示す図であり、図２１（ｂ）は、ジオメトリ情報生成部５０６での処理２を模式的に示す図である。
【０１５３】
図２１（ａ）に示すように、ジオメトリ情報生成部５０６は、ボリュームデータｔｋ−１における、３次元的に表現されている特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１を生成する。次に、図２１（ｂ）に示すように、ジオメトリ情報生成部５０６は、特定の領域Ｒのジオメトリ情報Ｇｔｋ−１及び、あらかじめ動き解析情報生成部５０３から取得したボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報から、特定の領域Ｒのジオメトリ情報Ｇｔｋを生成する。ジオメトリ情報Ｇｔｋはジオメトリ情報Ｇｔｋ−１に含まれる座標情報をその座標における動き解析情報による動きを適用して、ジオメトリ情報Ｇｔｋにおける座標を取得し、この時に座標情報(例えばノード)の属するマスク・パス・点の関係(例えばリンク)を維持することによって求めることが出来る。
【０１５４】
また、ジオメトリ情報生成部５０６は、操作部５１３の入力指示に基づき、必要に応じて、生成した各ボリュームデータのジオメトリ情報を、補間ボリュームデータ生成部５０５に出力する。
【０１５５】
また、補間ボリュームデータ生成部５０５は、操作部５１３の入力指示に基づき、必要に応じて、動き解析情報生成部５０３で生成した動き解析情報、及びジオメトリ情報生成部５０６で生成したジオメトリ情報Ｇｔｋに基づき、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘを生成することができる。
【０１５６】
ここで、本実施の形態では、実施の形態１と同様、補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報とは限らない。言い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、実施の形態１で上述した理由と同じであるため、詳細な説明は省略する。
【０１５７】
画像生成部５０７は、ボリュームデータ生成部５０１から出力されるボリュームデータｔｋ−１、ｔｋ、とその時刻情報ｔｋ−１、ｔｋ、補間ボリュームデータ生成部５０５から出力される補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘ、を入力とする。そして、画像生成部５０７は、補間ボリュームデータｔａｕｘ、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像、及び、補間ジオメトリ情報Ｇｔａｕｘ、各ジオメトリ情報Ｇｔｋと重ね合わせた画像を生成し、時系列に並べる。この時に、補間ジオメトリ情報Ｇｔａｕｘ、各ジオメトリ情報Ｇｔｋに含まれる数値を可視化する。数値は例えば、組織の体積(心室の体積、肺の体積、腫瘍の体積)、表面積(腫瘍が浸潤した面の面積)、長さ(血管の長さ、心壁の厚さ)、位置(弁の座標)が挙げられる。また、補間ジオメトリ情報Ｇｔａｕｘ、各ジオメトリ情報Ｇｔｋに含まれる数値の変化を可視化する。例えば、肺活量、腫瘍の成長などがある。
【０１５８】
表示処理部５０９は、画像生成部５０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部５０７で生成された３次元画像のシーケンスを生成して、表示部５１１で動画像として再生するための処理（再生フレームレートの設定など）を行うことができる。
【０１５９】
表示部５１１は、表示処理部５０９で適宜処理された、時系列の実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。表示部５１１は、時系列となった実画像ｔｋ−１、ｔ
ｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部５１１は、画像生成部５０７で生成された３次元画像のシーケンスを、表示処理部５０９での処理の後、動画像として再生することができる。
【０１６０】
＜画像処理の例（１）＞
次に、図２２を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２２は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（１）を示す図である。
【０１６１】
図２２に示すように、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集をすると、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報に基づき、インタラクティブに、その編集内容を他の３次元画像ｔｋに伝播させることはできる。これは、あらかじめ、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報と、ジオメトリ情報Ｇｔｋ−１とにより、他の３次元画像ｔｋの特定の領域Ｒについて、ジオメトリ情報Ｇｔを生成しているからである。また、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集・保存したものに、後から、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報を適用し、その編集内容を他の３次元画像ｔｋに伝播させることができる。これによってジオメトリ情報Ｇｔｋ−１作成の処理やユーザの労力が軽減される。
【０１６２】
また、ジオメトリ情報Ｇｔに含まれるマスク・パス・点などの編集をすると、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報を修正する事が出来る。これは、ジオメトリ情報Ｇｔとジオメトリ情報Ｇｔｋ−１との対応関係が定められているからである。これによって動き解析情報を容易に修正することが出来る。
【０１６３】
＜画像処理の例（２）＞
次に、図２３を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２３は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（２）を示す図である。
【０１６４】
図２３では、まず、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集をすると、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報に基づき、インタラクティブに、その編集内容を他の３次元画像ｔｋに伝播させることができる。さらに、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集した後に、同じ動き解析情報を用いて、補間ジオメトリ情報Ｇｔａｕｘ、及び補間画像ｔａｕｘ３を生成成することができる。このとき、補間画像ｔａｕｘ３の時刻情報ｔａｕｘ３は、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない、非等間隔な時刻で生成することできる。
なお、補間ジオメトリ情報Ｇｔａｕｘとは、ジオメトリ情報と同様、１つの３次元補間画像について、中実、表面、線、位置、又はこれらの要素の組み合わせで表現した組織の情報である。
【０１６５】
そして、補間画像ｔａｕｘ３の時刻情報ｔａｕｘ３を決める際に、あらかじめ生成していたボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報と、ジオメトリ情報Ｇｔｋ−１を利用する。たとえば、特定の領域Ｒが心臓である場合、その容積が最大となる時刻を、補間画像ｔａｕｘ３の時刻情報ｔａｕｘ３に設定すれば、特定の領域Ｒの時間変化を容易に画像化することができる。
【０１６６】
ここで、補間ボリュームデータとは、時刻ｔｋ−１と時刻ｔｋとの間の特定の時刻ｔａｕｘにおけるボリュームデータである。つまり、補間ボリュームデータ生成部５０５が補間ボリュームデータｔａｕｘを生成することは、実際はコンピュータ断層撮影装置で撮像されていない時刻ｔａｕｘにおいて、あたかもコンピュータ断層撮影装置により撮像して、ボリュームデータ生成部５０１がボリュームデータを生成したことと、同等となる。
【０１６７】
＜画像処理の例（３）＞
次に、図２４を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２４は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（３）を示す図である。
【０１６８】
図２４に示すように、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報から、特定の領域Ｒについて、時刻情報ｔｋ−１から時刻情報ｔｋまでの変形変異を取得し、この変形変異を見込んで、ボリュームデータｔｋ−１について特定の領域Ｒのボクセル値の変化を予測し補正する補正処理を行う際に、ジオメトリ情報Ｇｔｋ−１を利用することができる。つまり、図２４に示すように、特定の領域Ｒは、実画像ｔｋ−１から実画像ｔｋへ遷移する間に縮小する。特定の領域Ｒが縮小すると、特定の領域Ｒでは造影剤などが濃縮され、何も補正処理を実行しなければ、実画像ｔｋ上、特定の領域Ｒの画素値が大きくなってしまう。これは、造影剤の注入の時刻情報を加味すると更に効果的である。
【０１６９】
そこで、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報から、領域Ｒの変形変異を取得し、実画像ｔｋ−１に対して、ジオメトリ情報Ｇｔｋ−１に基づき、変形変異を見込んだ補正処理を行う。
【０１７０】
なお、上述した変形変異を見込んだ補正処理は、特定の領域Ｒに含まれる組織によって、その補正の度合を適宜変更するのが良い。補正の度合とは、例えば、ＣＴ検査における、「臓器の収縮」、「近接している高信号の部分の影響」であり、ＭＲ検査における、「位置の変化」、超音波検査における、「速度の変化」などがある。
【０１７１】
＜画像処理の例（４）＞
次に、図２５を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２５は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（４）を示す図である。
【０１７２】
図２５に示すように、ボリュームレンダリング法（ＶＲ）、最大値投影法（ＭＩＰ）、ＭＰＲ（Ｍｕｌｔｉ−Ｐｌａｎａｒ Ｒｅｎｄｅｒｉｎｇ）、ＣＰＲ（Ｃｕｒｖｅｄ Ｐｌａｎａｒ Ｒｅｆｏｒｍａｔｉｏｎ）、仮想内視鏡法（ＶＥ）など、医療用３ＤＷＳの既存の種々のレンダリングでの実画像ｔｋ−１において、特定の領域Ｒの注目点を画面上で固定する。そして、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報、及び実画像ｔｋ−１における特定の領域Ｒのジオメトリ情報Ｇｔｋ−１を利用することで、特定の領域Ｒ内のある領域又は点に着目して、表示部５１１上でなるべく動かないように、視点の位置・方向・回転角度・拡大率を、実画像毎に設定することができる。
【０１７３】
＜画像処理の例（５）＞
次に、図２６を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２６は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（５）を示す図である。
【０１７４】
図２６に示すように、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報、及び実画像ｔｋ−１における特定の領域Ｒのジオメトリ情報Ｇｔｋ−１を利用することで、仮想内視鏡法（ＶＥ）で得られる３次元実画像ｔｋ−１上で、視点を特定の領域Ｒ内で固定し、その視点の位置及び方向を、動き解析情報に従い、各３次元実画像ごとに設定することができる。
【０１７５】
以上、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成することで、観察対象の形状や大きさなどのより客観的な指標で把握することができる。
【０１７６】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成することで、個別にジオメトリ情報を生成するのに比較して、
セグメンテーションに係る操作の手間、補正処理を減らすことができる。さらに、実画像と他の実画像との画質に差異によって、セグメンテーションの結果の揺らぎを抑止できる。そのため、ユーザは、より客観的な画像解析を行うことができる。
【０１７７】
さらには、本実施の形態に係る医療用画像処理装置５００では、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報から他の実画像のジオメトリ情報を生成しているので、他の実画像のジオメトリ情報は、内在的に実画像のボリュームデータから生成したジオメトリ情報との対応関係を有している。そのため、他の実画像のジオメトリ情報と実画像のジオメトリ情報との対応関係を改めて取得、計算する必要がない。
【０１７８】
本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる
【０１７９】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。
【０１８０】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。
【０１８１】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、セグメンテーションにより隣接する臓器を設定し、動き解析情報により動きが一緒の部分を推定、癒着の可能性あり、として表示部５１１に表示することができる。
【０１８２】
なお、本実施の形態に係る医療用画像処理装置５００では、コンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部５０１に供給されているが、これに限らない。ボリュームデータ生成部５０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
【０１８３】
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。
【０１８４】
なお、本実施の形態において、表示処理部５０９が、画像生成部５０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部５１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部５１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
【０１８５】
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻に基づき、生成しているが、これに限らない。
【０１８６】
なお、本実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。
【符号の説明】
【０１８７】
１００、３００、５００ 医療用画像処理装置
１０１、３０１、５０１ ボリュームデータ生成部
１０３、３０３、５０３ 動き解析情報生成部
１０５、３０５、５０５ 補間ボリュームデータ生成部
１０７、３０７、５０７ 画像生成部
１０９、３０９、５０９ 表示処理部
１１１、３１１、５１１ 表示部
１１３、３１３、５１３ 操作部
３０６ 補正処理部
５０６ ジオメトリ情報生成部
【技術分野】
【０００１】
本発明は、医療用画像処理装置及び医療用画像処理プログラムに関する。
【背景技術】
【０００２】
Ｘ線ＣＴ装置およびＭＲＩ装置等の各種画像診断機器は、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行うものである。例えば、心臓や血管等の循環器系およびその他の動きのある臓器の場合、それらを構成する組織の動きを断層像により観察して、それら臓器等の機能を診断することが行なわれている。このような画像診断機器により、高解像度の画像が短時間に得られるようになったことに伴い、取得したボリュームデータから、臓器や腫瘍等の関心領域（Ｒｅｇｉｏｎ ｏｆ Ｉｎｔｅｒｅｓｔ）を抽出し、これらを見やすく可視化し、あるいは定量化して面積や体積を量ることにより病変部の診断に役立てている。
【０００３】
従来例１として、特許文献１に開示されている医療用画像処理装置は、各フェイズの複数のボリュームデータから新たに一の画像を生成し、各フェイズ現れる組織の変化を一の画像に表示することができる。しかし、各フェイズ間を補うボリュームデータを生成しているものではない。ここで、フェイズ（Ｐｈａｓｅ）とは、同一の対象を同一手法で短時間のうちに撮像した複数のボリュームデータのうちの１つのボリュームデータを指す。
【０００４】
従来例２として、特許文献２に開示されている医療用画像処理装置は、各フェイズの抽出領域から構造情報（たとえば、血管の中心線）を求めて、各フェイズの異常個所を比較、検出している。そのため、少ない計算量で領域抽出を行うことができる。しかし、各フェイズ間を補うボリュームデータを生成し、各フェイズの領域抽出を行っていない。
【０００５】
従来例３として、非特許文献１では、いわゆる、モーフィングにより、医療画像とは異なる通常の画像で構成される動画像のフレーム間隔（例えば、１／６０秒）の中間位置で新たに画像を補間している。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２００８−３５８９５号公報
【特許文献２】特開２００７−６８６５８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００７】
しかし、従来例１では、複数のボリュームデータから高画質な一の画像を生成しているが、複数のボリュームデータそれぞれから作成する個別の画像の画質は改善しない。これは特に動画となったときに顕著である。また、従来例２は複数のボリュームデータに含まれる組織の対応関係を取得するに留まる。また、従来例１及び従来例２で開示されている技術では、各フェイズ間を補うボリュームデータを生成しているものではないので、各フェイズの位置合わせの精度を高めるには別の方法を検討する必要がある。また、従来例３では、新たに画像を補間するタイミングは、動画像（たとえば、テレビ映像）のフレーム間隔（例えば、１／６０秒）の中間位置に固定されており、この位置を変更することはまでは示唆されていない。さらに、新たに画像を補間しても画質を改善しているわけではない。また、観察対象の領域を画像としてよりも、その形状や大きさなどのより客観的な指標で把握したい場合がある。この場合は観察対象のジオメトリ情報を作成するが、補間ボリュームを含む多数のボリュームデータそれぞれにジオメトリ情報を作成するには多大な労力がかかる。また、異なるボリュームデータのジオメトリ情報の間の対応関係を設定するのも多大な労力がかかる。
【０００８】
本発明の目的は、動き解析情報及びジオメトリ情報に基づき、複数のボリュームデータからそれぞれジオメトリ情報を作成するときの処理やユーザの労力を軽減できる医療用画像処理装置及び医療用画像処理プログラムを提供することである。
【課題を解決するための手段】
【０００９】
本発明は、時刻を有する複数のボリュームデータのそれぞれを可視化する医療用画像処理装置であって、前記複数のボリュームデータのうち、一のボリュームデータについて、ジオメトリ情報を生成する第１生成手段と、前記複数のボリュームデータのうち、前記一のボリュームデータと他の少なくとも一のボリュームデータとの間の動き解析情報を生成する第２生成手段と、前記複数のボリュームデータのうち、前記他の少なくとも一のボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する新たなジオメトリ情報を生成する第３生成手段と、前記他の少なくとも一のボリュームデータについて、前記対応する新たなジオメトリ情報を含む画像を表示するための処理を行う表示処理手段と、を備える医療用画像処理装置を提供する。
【００１０】
上記医療用画像処理装置では、前記動き解析情報に基づき、前記一のボリュームデータに含まれる特定の領域の変異又は変形を取得し、当該取得した変異又は変形、前記ジオメトリ情報、及び前記対応する新たなジオメトリ情報に基づき、ボクセル値の変化を予測し、前記一のボリュームデータのボクセル値を補正する。
【００１１】
上記医療用画像処理装置では、さらに、前記一のボリュームデータと前記他の少なくとも一のボリュームデータとの間を補間する補間ボリュームデータを生成する第４生成手段と、前記補間ボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する補間ジオメトリ情報を生成する第５生成手段とを備え、前記表示処理手段は、前記補間ボリュームデータについて、前記補間ジオメトリ情報を含む画像を表示するための処理を行う。
【００１２】
上記医療用画像処理装置では、前記一のボリュームデータの時刻情報及び前記補間ボリュームデータの時刻情報の間隔と、前記他の少なくとも一のボリュームデータ及び前記補間ボリュームデータの時刻情報の間隔とは非等間隔である。
【００１３】
上記医療用画像処理装置では、前記表示処理手段は、前記複数のボリュームデータについて描画する時に、前記ジオメトリ情報と前記新たなジオメトリ情報とを利用して、前記複数のボリュームデータで共通する特定の領域を画像として表示するための処理を行う。
【００１４】
また、本発明は、コンピュータを、上記医療用画像処理装置の各手段として機能させるための医療用画像処理プログラムを提供する。
【発明の効果】
【００１５】
本発明に係る医療用画像処理装置及び医療用画像処理プログラムによれば、動き解析情報及びジオメトリ情報に基づき、複数のボリュームデータからそれぞれジオメトリ情報を作成するときの処理やユーザの労力を軽減できる。また、客観性を増すことができる。また、ジオメトリ情報に含まれる情報を表示することによって診断に寄与できる。
【図面の簡単な説明】
【００１６】
【図１】実施の形態１に係る医療用画像処理装置１００の構成を示すブロック図
【図２】動き解析情報生成部１０３の動作を概念的に説明するための示す図
【図３】補間ボリュームデータ生成部１０５の動作を概念的に説明するための図
【図４】画像生成部１０７が生成する画像を模式的に示す図
【図５】時刻ｔｋ−１と時刻ｔｋと間の任意の時刻で画像生成部１０７が生成する画像を模式的に示す図
【図６】心室容量の大きさの経時変化を示すモデルと画像生成部１０７が作成した補間画像ｔａｕｘとを比較した図
【図７】画像生成部１０７が生成した２組の時系列に並んだ３次元画像セットを示す図
【図８】図７に示す２組みの３次元画像セットを同期させた画像のシーケンスを示す図
【図９】等間隔に複数の補間画像を実画像に補った場合の画像のシーケンスを示す図
【図１０】実画像に対して一時的に複数の補間画像補い、表示した例を示す図
【図１１】医療用画像処理装置１００の動作処理フロー
【図１２】実施の形態２に係る医療用画像処理装置３００の構成を示すブロック図
【図１３】補正処理前後の実画像の例（１）を示す図
【図１４】医療用画像処理装置３００補正処理フローの例（１）
【図１５】医療用画像処理装置３００補正処理フローの例（２）
【図１６】実画像の補正処理の例（２）を説明するための図
【図１７】実画像の補正処理の例（３）を説明するための図
【図１８】実画像及び補間画像の補正処理の例（４）を説明するための図
【図１９】実画像の補正処理の例（５）を説明するための図
【図２０】実施の形態３に係る医療用画像処理装置５００の構成を示すブロック図
【図２１】（ａ）ジオメトリ情報生成部５０６での処理１を模式的に示す図、（ｂ）ジオメトリ情報生成部５０６での処理２を模式的に示す図
【図２２】医療用画像処理装置５００の画像処理の例（１）を示す図
【図２３】医療用画像処理装置５００の画像処理の例（２）を示す図
【図２４】医療用画像処理装置５００の画像処理の例（３）を示す図
【図２５】医療用画像処理装置５００の画像処理の例（４）を示す図
【図２６】医療用画像処理装置５００の画像処理の例（５）を示す図
【発明を実施するための形態】
【００１７】
以下、本発明の実施形態について、図面を参照して説明する。
【００１８】
図１を参照して、実施の形態１に係る医療用画像処理装置について説明する。図１に示す医療用画像処理装置１００は、ボリュームデータ生成部１０１と、動き解析情報生成部１０３と、補間ボリュームデータ生成部１０５と、画像生成部１０７と、表示処理部１０９と、を備える。医療用画像処理装置１００には、医療用画像処理装置１００で生成された画像を表示する表示部１１１と、医療用画像処理装置１００の各部を操作するためのユーザインターフェースである操作部１１３と、が設けられている。
【００１９】
図１に示すように、医療用画像処理装置１００は、コンピュータ断層撮影装置に接続されている。コンピュータ断層撮影装置は、被検体の組織等を可視化するものである。Ｘ線源１からは同図に鎖線で示す縁部ビームを有するピラミッド状のＸ線ビーム束２が放射される。Ｘ線ビーム束２は、例えば患者３である被検体を透過しＸ線検出器４に照射される。Ｘ線源１及びＸ線検出器４は、本実施形態の場合にはリング状のガントリー５に互いに対向配置されている。リング状のガントリー５は、このガントリーの中心点を通るシステム軸線６に対して、同図に示されていない保持装置に回転可能（矢印ａ参照）に支持されている。
【００２０】
患者３は、Ｘ線が透過するテーブル７上に寝ている。このテーブル７は、図示されていない支持装置によりシステム軸線６に沿って移動可能（矢印ｂ参照）に支持されている。従って、Ｘ線源１及びＸ線検出器４は、システム軸線６に対して回転可能でありかつシステム軸線６に沿って患者３に対して相対的に移動可能である測定システムを構成するので、システム軸線６に関して種々の投影角及び種々の位置のもとで患者３にＸ線を投射する。その際に発生するＸ線検出器４の出力信号は、ボリュームデータ生成部１０１に供給され、ボリュームデータに変換される。
【００２１】
シーケンス走査の場合、患者３の層毎の走査が行なわれる。その際に、Ｘ線源１及びＸ線検出器４はシステム軸線６を中心に患者３の周りを回転し、Ｘ線源１及びＸ線検出器４を含む測定システムは患者３の２次元断層を走査するために多数の投影を撮影する。その際に取得された測定値から、走査された断層を表示する断層像が再構成される。相連続する断層の走査の間に、患者３はその都度システム軸線６に沿って移動される。この過程は全ての関心断層が捕捉されるまで繰り返される。
【００２２】
スパイラル走査中は、Ｘ線源１及びＸ線検出器４を含む測定システムはシステム軸線６を中心に回転し、テーブル７は連続的に矢印ｂの方向に移動する。すなわち、Ｘ線源１及びＸ線検出器４を含む測定システムは、患者３に対して相対的に連続的にスパイラル軌道上を、患者３の関心領域が全部捕捉されるまで移動する。本実施形態の場合、同図に示されたコンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部１０１に供給される。
【００２３】
ボリュームデータ生成部１０１は、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部１０３及び画像生成部１０７へ出力する。なお、ボリュームデータ生成部１０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータを保持することができる。
【００２４】
ここで、本実施の形態では、ボリュームデータ生成部１０１は、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
【００２５】
以下、説明のため、ボリュームデータ生成部１０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部１０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部１０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。なお、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
【００２６】
動き解析情報生成部１０３は、ボリュームデータ生成部１０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成し（図２参照）、補間ボリュームデータ生成部１０５に出力する。
【００２７】
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。
【００２８】
図２は、動き解析情報生成部１０３の動作を概念的に示す図である。図２に示すように、動き解析情報生成部１０３は、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを、動き解析情報として生成している。
【００２９】
補間ボリュームデータ生成部１０５は、動き解析情報生成部１０３から出力された動き解析情報に基づき、ボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、画像生成部１０７に出力する。
【００３０】
ここで、補間ボリュームデータとは、時刻ｔｋ−１と時刻ｔｋとの間の特定の時刻ｔａ
ｕｘにおけるボリュームデータである。つまり、補間ボリュームデータ生成部１０５が補間ボリュームデータｔａｕｘを生成することは、実際はコンピュータ断層撮影装置で撮像されていない時刻ｔａｕｘにおいて、あたかもコンピュータ断層撮影装置により撮像して、ボリュームデータ生成部１０１がボリュームデータを生成したことと、同等となる。
【００３１】
図３は、補間ボリュームデータ生成部１０５の動作を概念的に示す図である。図３に示すように、補間ボリュームデータ生成部１０５は、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成する。
【００３２】
ここで、本実施の形態では、図３に示すように補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない。言
い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、後述する。
【００３３】
画像生成部１０７は、ボリュームデータ生成部１０１から出力されるボリュームデータｔｋ−１、ｔｋ、及び時刻情報ｔｋ−１、ｔｋ、並びに補間ボリュームデータ生成部１０５から出力される補間ボリュームデータｔａｕｘ、及びその時刻情報ｔａｕｘを入力とする。画像生成部１０７は、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成する一方、補間ボリュームデータｔａｕｘから補間３次元画像を生成し、時系列に並べる。
【００３４】
表示処理部１０９は、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部１０７で生成された３次元画像のシーケンスを生成して、表示部１１１で動画像として再生するための処理（再生フレームレートの設定など）を行うことができる。
【００３５】
表示部１１１は、表示処理部１０９で適宜処理された時系列となった実画像ｔｋ−１、
ｔｋ及び補間画像ｔａｕｘを表示する。表示部１１１は、時系列となった実画像ｔｋ−１
、ｔｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部１１１は、画像生成部１０７で生成された３次元画像のシーケンスを、表示処理部１０９での処理の後、動画像として再生することができる。
【００３６】
操作部１１３は、表示部１１１に表示された画像を加工するための各種のインターフェースを備えており、例えば、補間画像ｔａｕｘを生成するための補間パラメータ（時刻情報ｔａｕｘなど）を設定する。そして、ユーザは、表示部１１１を見ながら、操作部１１３により、医療用画像処理装置１００の各部を操作することができる。そのため、ユーザは、インタラクティブに補間パラメータを変更できる。したがって、本実施の形態の医療用画像処理装置において、画像生成部１０７は、操作部１１３から指示に基づき、インタラクティブに補間画像ｔａｕｘを修正することができる。
【００３７】
ここで、補間パラメータとして、たとえば、（１）位置補間アルゴリズム（リニア／キュービック／スプライン）、（２）ボクセル値の補間のアルゴリズム（リニア／キュービック／スプライン）、（３）補間により生成される３Ｄ画像のピッチ（時間間隔）や数、並びに（４）４Ｄフィルタの方法として、平滑化、エッジ強調、最大値・最小値投影、差分、累積加算、ヒストグラムマッチング、およびそれらの組み合わせなどがある。
【００３８】
図４を参照して、画像生成部１０７が生成する３次元画像について説明する。図４は、画像生成部１０７が生成した各画像を模式的に示す図である。図４に示すように、実画像ｔｋ−１と実画像ｔｋとの間で、時刻情報ｔａｕｘを有する補間画像ｔａｕｘが、画像生成部１０７によって生成される。以下、説明のためボリュームデータｔｋ−１、ｔｋから生成された、時刻ｔｋ−１、ｔｋにおける実際の３次元画像を、それぞれ実画像ｔｋ−１と実画像ｔｋと称し、補間ボリュームデータｔａｕｘから生成された、時刻ｔａｕｘにお
ける３次元の補間画像を、補間画像ｔａｕｘと称する。
【００３９】
ここで、図４に示すように、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミング
に生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻
情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが非等間隔となっている。本実施の形態において、ボリュームデータｔｋ−１と、時刻情報ｔｋ−１の次のタイミングで生成されるボリュームデータｔｋとの間隔が、およそ０．１秒から数年単位である。これは、画像生成部１０７で生成される医療画像（実画像ｔｋ−１、実画像ｔｋ、及び補間画像ｔａｕｘ）とは異なる通常の画像で構成される動画像（たとえば、テレビなどの映像）のフレーム間隔（例えば、１／６０秒）と比べて非常に大きいことに起因する。
【００４０】
つまり、本実施の形態では、実画像ｔｋ−１、と実画像ｔｋとの間の間隔（Ｃ＋Ｄ）が、例えば年単位であって、実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の関心領域の経年変化を観察するのに、その中間位置での補間画像では十分ではないからである。そのため、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが非等間隔とし、実画像ｔｋ−１、と実画像ｔｋとの間で被検体である患者３の検査部位（関心領域）の経年変化を中間位置とは異なる、実画像ｔｋ−１、と実画像ｔｋとの間の様々な位置で補間画像を生成することを可能としている。また、撮影に要する放射線量を最小にする目的で、関心の最も高い時刻では撮影の時間間隔を密にし、それ以外では疎に撮影することが行われている。また、心電図など被検体のリアルタイム情報を取得する別途の機器で指定されるタイミングで撮影することもある。これらによって撮影間隔が非等間隔となりうる。リアルタイム情報には他には例えば、造影剤の注入のタイミングもある。
【００４１】
上述のように、本実施の形態では、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイ
ミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘ
の時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが、非等間隔であるが、これに限らない。実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の関心領域の経時変化を観察するのに適したタイミングであれば、補間画像ｔａｕｘの時刻情報ｔａｕｘと直前のタイミングに生成された実画像ｔｋ−１の時刻情報ｔｋ−１との間隔Ｃと、補間画像ｔａｕｘの時刻情報ｔａｕｘと直後のタイミングに生成された実画像ｔｋの時刻情報ｔｋとの間隔Ｄとが、等間隔に設定されてもよい。すわなち、実画像ｔｋ−１と実画像ｔｋとの間で発生している被検体である患者３の検査部位（関心領域）の経時変化を観察するのに適したタイミングであれば、時刻情報ｔｋ−１と時刻情報ｔｋ−１との間の任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成しても良い。言い換えると、本実施の形態において、ボリュームデータｔａｕｘの時刻ｔａｕｘは、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。これにより、本実施の形態では、Ｘ線の線量を削減し、被検体である患者３への影響を極力抑えることができる。
【００４２】
（任意の時刻の補間画像）
図５は、時刻ｔｋ−１と時刻ｔｋ−１との任意の時刻で画像生成部１０７が補間画像ｔ
ａｕｘを生成する一例を模式的に示す図である。図５に示すように、画像生成部１０７は、実画像ｔｋ−１を生成した時刻０．０ｓｅｃと、実画像ｔｋを生成した時刻情報１．５ｓｅｃとの間の任意の時刻情報１．０ｓｅｃの補間画像ｔａｕｘを生成する。そのため、実画像を生成した時刻の間の任意の時刻で補間画像を生成することができ、実画像ｔｋ−１、と実画像ｔｋとの間で発生している患者３の診断範囲の経時変化を観察することができる。
【００４３】
（補間画像とモデルとの比較）
図６は、心室容量の大きさの経時変化を示すモデルと画像生成部１０７が作成した補間画像ｔａｕｘとを比較した図である。上述のように、時刻情報ｔｋ−１と時刻情報ｔｋ−１との間の時刻であれば、任意の時刻で画像生成部１０７が補間画像ｔａｕｘを生成することが可能なので、たとえば図６に示すように、心室容量の大きさの経時変化を示すモデルの最大容量を示す時刻で、画像生成部１０７が補間画像ｔａｕｘを作成し、心室容量の大きさの経時変化を示すモデルと補間画像ｔａｕｘにあらわれる心室の大きさとを比較することができる。
【００４４】
心室容量の大きさの経時変化を示すモデルは、図６に示すように線形的に変化するものではないので、２つの実画像ｔｋ−１、ｔｋ間の中間位置では、心室容量は最大とならない、したがって、本実施の形態では、モデル上、心室容量が最大となる特徴を有する時刻で、実画像ｔｋ−１、ｔｋから補間画像ｔａｕｘを生成する。図６に示すように、モデル上、心室容量が最大となる時刻はｔ＝０．７８ｓｅｃである。したがって、画像生成部１０７は、時刻情報ｔａｕｘをモデル上、心室容量が最大となる時刻はｔ＝０．７８ｓｅｃとし、補間画像ｔａｕｘを生成する。
【００４５】
上述のように、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴が現れると推定される時刻を補間画像の時刻情報とすることで、画像生成部１０７は、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴を画像として生成し、表示部１１１が可視化することができる。
【００４６】
なお、心室容量の大きさの経時変化を示すモデルから、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴が現れると推定しているが、観察対象となる被検体（患者３）の検査部位（関心領域）の特徴を推定する手段は当該モデルだけに限らない。例えば、被検体である患者３に対して造影剤などが注入された外的要因が発生した時刻を、画像上に特徴が現れる時刻であると推定することも可能できる。そのため、補間画像ｔａｕｘを作成し、実画像ｔｋ−１、ｔｋ、及び補間画像ｔａｕｘの一連の画像の変化から、被検体（患者３）の検査部位（関心領域）の外的要因による変化を観察することができる。また、時刻情報ｔａｕｘモデルのみならず外部機器からの入力によって定めることが出来る。これは心電図など被検体のリアルタイム情報を実画像の取得と同時に行い、これを参照することによって求めることが出来る。例えば、心電図上の波形の特徴的な時刻や、造影剤の注入を行ったタイミングなどが考えられる。
【００４７】
（異なる画像セットの同期）
また、画像生成部１０７は、同一の被検体から同一の検査部位を対象として、異なるモードで動作するコンピュータ断層撮影装置又は異なるコンピュータ断層撮影装置により供給された複数のボリュームデータから、時系列の複数の３次元画像セットを生成し、複数の３次元画像セットから正確に同期した画像のシーケンスを生成することができる。図７は、画像生成部１０７が生成した２組の時系列に並んだ３次元画像セットを示す。また、図８は、図７に示す２組みの３次元画像セットを同期させた画像のシーケンスを示す図である。
【００４８】
図７に示す３次元画像セットＢでは、第１の時刻情報ｔｋ−１（＝０．２ｓｅｃ）を有する第１実画像ｔｋ−１、及び第２の時刻情報ｔｋ（＝１．７ｓｅｃ）を有する第２実画像ｔｋ、並びに、第１実画像ｔｋ−１及び第２実画像ｔｋから生成され、時刻情報ｔａｕｘ２を有する第２補間画像ｔａｕｘ２が時系列に並んでいる。また、図７に示す３次元画像セットＡでは、第３の時刻情報ｔｋ−１（＝０．０ｓｅｃ）を有する第３実画像ｔｋ−１、及び第４の時刻情報ｔｋ（＝１．５ｓｅｃ）を有する第４実画像ｔｋ、並びに第３実画像ｔｋ−１及び第４実画像ｔｋから生成され、時刻情報ｔａｕｘ１を有する第１補間画像ｔａｕｘ１が時系列に並んでいる。
【００４９】
そして、図８に示すように、画像生成部１０７は、各画像の保持する時刻情報に基づき、３次元画像セットＢの第２補間画像ｔａｕｘ２を有する時刻情報ｔａｕｘ２を、３次元画像セットＡに時間的に同期させて、新たな画像のシーケンスを生成する。
【００５０】
したがって、本実施の形態では、画像生成部１０７は、同一の被検体から同一の検査部位を対象として、異なるモードで動作するコンピュータ断層撮影装置又は異なるコンピュータ断層撮影装置により供給された複数のボリュームデータから、時系列の複数の３次元画像セットを生成することができる。そして、表示処理部１０９は、複数の３次元画像セットから正確に時間的に同期した画像のシーケンスを生成することができる。
【００５１】
また、本実施の形態では、画像生成部１０７は、図８に示す３次元画像のシーケンスを用いれば、観察対象部位（心臓）のマルチフェイズデータでの同一フェイズ間のずれの補正も可能となる。また、本実施の形態では、画像生成部１０７は、図８に示す３次元画像のシーケンスを用いれば、「コンピュータ断層撮影装置から取得した検査の画像取得時間情報」と「表示処理部１０９で生成される再生フレームレート」を元に、表示部１１１で正確な実時間で動画再生することができる。
【００５２】
さらに、本実施の形態では、画像生成部１０７は、図７に示す一方の３次元画像セットＢにだけ存在する特定の取得時刻（例えば、ｔｋ＝１．７ｓｅｃ）の実画像ｔｋに合わせて、他方の３次元画像セットＡの第３実画像ｔｋ−１及び第４実画像ｔｋから新たな補間画像を作成することにより、より正確な３次元画像間の比較を可能にする。
【００５３】
図９に、画像生成部１０７が時間的に等間隔に複数の補間画像を実画像に補った場合の画像のシーケンスを示す。また、画像生成部１０７は、図９に示すように、同一の被検体から同一の検査部位を対象として、撮影間隔の異なるコンピュータ断層撮影装置により供給された複数のボリュームデータｔ１（＝０．００ｓｅｃ）、ｔ２（＝０．１７ｓｅｃ）、ｔ３（＝０．３０ｓｅｃ）から生成された複数の実画像ｔ１、ｔ２、ｔ３のうち実画像ｔ１、ｔ３に対して、時間的に等間隔に補う、時刻情報ｔ４（＝０．１０ｓｅｃ）を有する補間画像ｔａｕｘ４、及び時刻情報ｔ５（＝０．２０ｓｅｃ）を有する補間画像ｔａｕｘ５を生成することができる。そのため、本実施の形態では、時間的に等間隔でない実画像に対して、補間画像を等時間間隔に生成するので、表示部１１１が表示する医療用画像の表示や解析の質・能力を高めることができる。
【００５４】
図１０に、実画像に対して一時的に複数の補間画像補い、表示部１１１に表示した例を示す。画像生成部１０７は、図１０に示すように、表示部１１１に３次元画像のシーケンスを表示している時のみ、一時的に、画像生成部１０７が複数の補間画像ｔｅｍｐ１から補間画像ｔｅｍｐ３を生成し、各フレームごとに、レンダリング後に複数の補間画像ｔｅｍｐ１から補間画像ｔｅｍｐ３を破棄することで、メモリの使用量を削減でき、本実施の形態の医療用画像処理装置１００への負荷を軽減することができる。
【００５５】
次に、図１１を参照して、本実施の形態に係る医療用画像処理装置１００の動作処理フローを説明する。図１１は、医療用画像処理装置１００の動作処理フローである。
【００５６】
まず、ステップＳＴ１で、ボリュームデータ生成部１０１は、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータ、及び各ボリュームデータの時刻情報ｔｎを生成する。そして、ステップＳＴ２へ遷移する。
【００５７】
次に、ステップＳＴ２で、動き解析情報生成部１０３は、ボリュームデータ生成部１０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。そして、ステップＳＴ３へ遷移する。
【００５８】
次に、ステップＳＴ３で、補間ボリュームデータ生成部１０５は、動き解析情報生成部１０３から出力された動き解析情報に基づき、ボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成する。そして、ステップＳＴ４へ遷移する。
【００５９】
次に、ステップＳＴ４で、画像生成部１０７は、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像を生成する一方、補間ボリュームデータｔａｕｘから補間３次元画像を生成し、時系列に並べる。そして、ステップＳＴ５へ遷移する。
【００６０】
次に、ステップＳＴ５で、表示処理部１０９は、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う。そして、ステップＳＴ６へ遷移する。
【００６１】
最後に、ステップＳＴ６で、表示部１１１は、表示処理部１０９で適宜処理された時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。そして、動作処理が終了する。
【００６２】
以上、本実施の形態に係る医療用画像処理装置１００によれば、動き解析情報に基づき補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示することができる。したがって、本実施の形態に係る医療用画像処理装置１００によれば、被検体である患者３の検査部位（関心領域）の観察を精度良く行うことができる。
【００６３】
例えば、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる。また、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。
【００６４】
さらに、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。さらに、また、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、セグメンテーションにより隣接する関心領域（臓器）を設定し、動き解析情報により複数の臓器の動きが同一の部分を推定、臓器間癒着の可能性あり、として表示部１１１に表示することができる。
【００６５】
また、本実施の形態に係る医療用画像処理装置１００によれば、動き解析情報に基づき補間画像を実画像間に補間すれば、表示部１１１のノイズ軽減、ノイズのばらつき調整が可能である。
【００６６】
例えば、医療用画像処理装置１００が生成した動き解析情報によって、実画像及び補間画像の対応画素を特定することにより、本来の詳細な画像情報を失うことなくノイズを軽減できる。また、医療用画像処理装置１００により、ＣＴマルチフェイズ検査において、Ｘ線照射量の多いフェイズと少ないフェイズが混在しても、Ｘ線照射量の少ないフェイズについては、前後のフェイズの画素値を利用することで、良好な補間画像を得られる。
【００６７】
さらに、医療用画像処理装置１００により、動き解析情報に基づき補間画像を実画像間に補間すれば、重要なフェイズのみＸ線照射量を多くして実画像を生成し、それ以外のフェイズについては、補間画像を生成することで、患者３のＸ線の被ばく量を低減することができる。
【００６８】
なお、本実施の形態に係る医療用画像処理装置１００では、図１に示すコンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部１０１に供給されているが、これに限らない。ボリュームデータ生成部１０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
【００６９】
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。
【００７０】
なお、本実施の形態において、表示処理部１０９が、画像生成部１０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部１１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部１１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
【００７１】
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻ｔに基づき、生成しているが、これに限らない。
【００７２】
なお、本実施の形態の説明に用いた各機能ブロック（図１参照）は、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。
【００７３】
なお、本実施の形態では、心臓の観察に用いる例を示したが、これに限らない。例えば、腫瘍組織の経時変化の観察、血流灌流の観察などに用いることも出来る。
【００７４】
（実施の形態２）
上述した実施の形態１に係る医療用画像処理装置１００によれば、動き解析情報に基づき、補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示する。しかし、補間画像を実画像間に補間する場合、実画像の画質によっては、単純に補間画像を補間するだけでは、画像内の特定の領域の変化を正確に表示できない場合もある。
【００７５】
そこで、実施の形態２に係る医療用画像処理装置３００では、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正して、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
【００７６】
実施の形態２に係る医療用画像処理装置３００では、必ずしも補間画像を生成する必要はない。動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質を補正して、２以上の実画像内の特定の領域の変化を正確に表示することができれば、補間画像を生成する必要はない。
【００７７】
図１２を参照して、本実施の形態２に係る医療用画像処理装置３００の構成について説明する。図１２に示す医療用画像処理装置３００は、ボリュームデータ生成部３０１と、動き解析情報生成部３０３と、補間ボリュームデータ生成部３０５と、補正処理部３０６と、画像生成部３０７と、表示処理部３０９と、を備える。さらに、医療用画像処理装置３００には、医療用画像処理装置３００で生成された画像を表示する表示部３１１と、医療用画像処理装置１００の各部を操作するためのユーザインターフェースである操作部３１３と、が設けられている。
【００７８】
ボリュームデータ生成部３０１は、実施の形態１と同様、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部３０３と、補正処理部３０６と、画像生成部３０７とへ出力する。
【００７９】
また、ボリュームデータ生成部３０１は、コンピュータ断層撮影装置から、装置の撮像条件（例えば、Ｘ線量、撮影時間、造影剤の注入の時刻）を取得し、画像生成部３０７へ出力する。撮像条件により、のちに生成される実画像の画質が大きく変化する。なお、コンピュータ断層撮影装置の撮像条件以外にも、コンピュータ断層撮影装置の撮影対象の条件（患者の動きや、患者の体内に含まれる金属など）により、のちに生成される実画像の画質が大きく変化する。
【００８０】
なお、ボリュームデータ生成部３０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータ、コンピュータ断層撮影装置の撮像条件、コンピュータ断層撮影装置の撮影対象の条件を保持することができる。
【００８１】
ここで、本実施の形態では、ボリュームデータ生成部３０１は、実施の形態１と同様、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
【００８２】
以下、説明のため、ボリュームデータ生成部３０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部３０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部３０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。また、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
【００８３】
動き解析情報生成部３０３は、ボリュームデータ生成部３０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。
【００８４】
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。
【００８５】
そして、動き解析情報生成部３０３は、生成した動き解析情報を、補間ボリュームデータ生成部３０５と画像生成部３０７とに出力する。動き解析情報生成部３０３が生成する動き解析情報は、実施の形態１の動き解析情報生成部１０３と同様（図２参照）、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを示す。
【００８６】
補間ボリュームデータ生成部３０５は、操作部３１３の入力指示に基づき、必要に応じて、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、補正処理部３０６に出力する。
【００８７】
ここで、本実施の形態では、実施の形態１と同様、補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない。言い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、実施の形態１で上述した理由と同じであるため、詳細な説明は省略する。
【００８８】
本実施の形態に係る医療用画像処理装置３００の一つの特徴として、補正処理部３０６は、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像（入力画像）の画質、又は２つの実画像間を補間する補間画像の画質を補正する。
【００８９】
＜補正処理１：画質の均一化＞
補正処理の一例として、コンピュータ断層撮影装置の撮影対象の条件などから、補正処理部３０６の補正処理前の時刻ｔｋ−１、ｔｋにおける画像の画質が均一でない場合、補正処理部３０６は、動き解析情報生成部３０３からの入力された動き解析情報、ボリュームデータ生成部３０１から入力されたボリュームデータｔｋのボクセル及びボリュームデータｔｋ−１のボクセルに基づき、画像生成部３０７で生成する実画像ｔｋ−１、ｔｋ、の画素値を求める。詳細は、図１３、図１４を参照して後述する。
【００９０】
図１３に、補正処理前後の実画像ｔｋ−１、ｔｋの例（１）を示す。なお、補正処理前の実画像（入力画像）は、説明のため示しているが、画像生成部３０７が実際に生成する画像（出力画像）は、補正処理部３０６で補正処理が行われた後の実画像ｔｋ−１、ｔｋである。図１３に示すように、補正処理前の時刻ｔｋ−１、ｔｋにおける画像の画質が均一ではないが、補正処理部３０６で補正処理が行われた後は、実画像ｔｋ−１、ｔｋの画質は均一となる。そのため、本実施の形態に係る医療用画像処理装置３００は、動き解析情報及び画像の画質情報に基づき、複数のボリュームデータからそれぞれ画像を作成するときの画質を向上させることができる。さらに、本実施の形態に係る医療用画像処理装置３００は、補正処理後の各画像の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の経時間変化を、正確に表示することができる。なお、必要に応じて、補正処理部３０６は、補間ボリュームデータ生成部３０５から出力された補間ボリュームデータｔａｕｘを、実画像ｔｋ−１、ｔｋと同様に補正処理を行う。
【００９１】
次に、図１４を参照して、実画像ｔｋ−１、ｔｋの画質を均一とするための補正処理について説明する。図１４は、本実施の形態に係る医療用画像処理装置３００補正処理フローの例（１）である。
【００９２】
図１４に示すように、ステップＳＴ３１では、出力フェイズの前後に相当する、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１を取得する。そして、ステップＳＴ３２へ遷移する。
【００９３】
ステップＳＴ３２では、ボリュームデータｔｋ、ｔｋ−１それぞれについて、ボクセル位置ごとに、信頼度（０〜１）を決定する。ここで、ボリュームデータｔｋ、ｔｋ−１についての信頼度（信頼性の高さ）は、出力フェイズとの時間的距離、ボリュームデータｔｋ、ｔｋ−１を撮影した時の撮影条件（Ｘ線照射量など）、当該ボリュームデータｔｋ、ｔｋ−１のＳ／Ｎ比、近傍の領域の移動速度や画素値の特徴などにより、決定することが考えられる。そして、ステップＳＴ３３へ遷移する。
【００９４】
ステップＳＴ３３では、出力フェイズの全ボクセルについて、ボリュームデータｔｋ、ｔｋ−１の画素値を、動き解析情報及びボリュームデータｔｋ、ｔｋ−１より求める。
そして、ステップＳＴ３４へ遷移する。
ステップＳＴ３４では、出力フェイズの全ボクセルについて、前後の入力フェイズでの対応する位置を、動き解析情報により求める。そして、ステップＳＴ３５へ遷移する。
【００９５】
ステップＳＴ３５では、各対応位置の画素値を入力画像上でサンプリングし、サンプリングした値を、入力フェイズに対応する入力画素値とする。そして、ステップＳＴ３６へ遷移する。
ステップＳＴ３６では、入力フェイズに対応する入力画素値の信頼度を取得する。
そして、ステップＳＴ３７へ遷移する。
【００９６】
ステップＳＴ３７では、入力フェイズに対応する入力画素値に、取得した信頼度を乗じ、当該入力画素値について、各入力フェイズで加算したものを、出力フェイズで、補正後のボリュームデータｔｋ、ｔｋ−１の画素値とする。そして、処理が終了する。
【００９７】
ここで、ステップＳＴ３２における、ボリュームデータｔｋ、ｔｋ−１の信頼度（信頼性の高さ）について、画素値の信頼性の高さに応じ、４Ｄフィルタのパラメータを各実画像ｔｋ、ｔｋ−１毎、又はボクセル毎に変えることで、全体の適正なＳ／Ｎを保ちつつ、元画像の空間分解能を最大限活用することができる。
【００９８】
画素値の信頼性は、時刻、ボクセルの位置、画素値、収縮度、速度、加速度、近傍の画素の変位から求めた回転量、変形の度合いなどの関数である。また、近傍の領域の特徴（高画素値の領域の有無、など）が、画素値の信頼性に影響することもある。
【００９９】
４Ｄフィルタの方法として、平滑化、エッジ強調、最大値・最小値投影、差分、累積加算、ヒストグラムマッチング、およびそれらの組み合わせなどがある。
【０１００】
＜補正処理２：画像の特性の均一＞
また、補正処理部３０６は、動き解析情報生成部３０３からの入力された動き解析情報、並びにボリュームデータ生成部３０１から入力されたボリュームデータｔｋ、のボクセル及びボリュームデータｔｋ−１のボクセルに基づき、画像生成部３０７で生成する実画像ｔｋ−１、ｔｋの画素値の値的な歪み量を求める。詳細は、図１５を参照して後述する。
【０１０１】
なお、補正処理部３０６が画像生成部３０７で生成する実画像ｔｋ−１、ｔｋの画素値の値的な歪み量を求めるのは、たとえば、コンピュータ断層撮影装置の撮影対象の条件(画質情報の元となる条件)などから、生成される画像の特性（明度、ダイナミックレンジ）が均一でない場合が考えられる。
【０１０２】
ここで、画素値の値的な歪み量とは、対応時刻、変位後の位置、対応位置の入力画素値、変位場の収縮度、変位の大きさ（速度）やその加速度、近傍での最大画素値や、その最大画素値のボクセルから当該ボクセルまでの距離、などの関数である。
【０１０３】
ボリュームデータ生成部３０１から入力されたボリュームデータｔｋ、のボクセル及びボリュームデータｔｋ−１のボクセルと、動き解析情報から画素値の歪の量の関数を求める方法としては、対応する位置の画素値の変化量を画像内で統計処理する、などが考えられる。
【０１０４】
そして、補正処理部３０６で求めた画素値の値的な歪み量に基づき、実画像ｔｋ−１、ｔｋを補正処理すると、生成される画像の特性（明度、ダイナミックレンジ）が均一な画像を作成することができる。
【０１０５】
＜補正処理フロー（２）：画像の特性（明度、ダイナミックレンジ）の均一化＞
ここで、図１５を参照して、実画像ｔｋ−１、ｔｋの画質、における特性（明度、ダイナミックレンジ）を均一とするための、補正処理部３０６で補正処理について説明する。図１５は、本実施の形態に係る医療用画像処理装置３００補正処理フローの例（２）である。
【０１０６】
図１５に示すように、ステップＳＴ４１では、出力フェイズの前後に相当する、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１を取得する。そして、ステップＳＴ４２へ遷移する。
【０１０７】
ステップＳＴ４２では、２以上の入力フェイズについて、ボリュームデータｔｋ、ｔｋ−１と動き解析情報から、画素値の値的な歪の量の関数を算出する。そして、ステップＳＴ４３へ遷移する。
ステップＳＴ４３では、歪の量の関数から、ボリュームデータｔｋ、ｔｋ−１の各ボクセル位置について、画素値の補正係数を算出する。そして、ステップＳＴ４４へ遷移する。
【０１０８】
ステップＳＴ４４では、出力フェイズの全ボクセルについて、動き解析情報及びボリュームデータｔｋ、ｔｋ−１から、ボリュームデータｔｋ、ｔｋ−１の画素値を算出する。
そして、ステップＳＴ４５へ遷移する。
ステップＳＴ４５では、出力フェイズの全ボクセルについて、前後の入力フェイズでの対応する位置を、動き解析情報から算出する。そして、ステップＳＴ４６へ遷移する。
【０１０９】
ステップＳＴ４６では、各対応位置の画素値を入力画像上でサンプリングし、サンプリングした値を、入力フェイズに対応する入力画素値とする。そして、ステップＳＴ４７へ遷移する。
ステップＳＴ４７では、上記入力画素値に、対応ボクセル位置の補正係数による変換をし、これをそれぞれの入力フェイズについて加算したものを、出力フェイズの画素値とする。そして、処理が終了する。
【０１１０】
＜補正処理３：変形又は変異による画素値の変化＞
また、本実施の形態に係る医療用画像処理装置３００の一つの特徴として、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報に基づき、動き解析情報から得られる入力画像内の特定の領域に関する変形又変位による画素値の変化を見込んだ補正処理を行うこともできる。図１６を参照して、説明する。
【０１１１】
図１６は、実画像ｔｋ−１、ｔｋの補正処理の例（２）を説明するための図である。説明のため、図１６では、入力画像ｔｋ−１、ｔｋ内の特定の領域を領域Ｒ１とする。図１６に示すように、領域Ｒ１は、実画像ｔｋ−１から実画像ｔｋへ遷移する間に縮小する。領域Ｒが縮小すると、領域Ｒ１では造影剤などが濃縮され、何も補正処理を実行しなければ、実画像ｔｋ上、領域Ｒ１の画素値が大きくなってしまう。これは、造影剤の注入の時刻情報を加味すると更に効果的である。
【０１１２】
そこで、補正処理部３０６は、動き解析情報生成部３０３から得られる動き解析情報に基づき、領域Ｒ１の変形変異を取得し、実画像ｔｋに対して変形変異を見込んだ補正処理を行う。
【０１１３】
上述のように、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報などに基づき、動き解析情報から得られる入力画像内の特定の領域Ｒ１に関する変形又変位による、画素値の変化を見込んだ補正処理を行うことで、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
【０１１４】
＜補正処理４：代表値での置き換え＞
図１７は、実画像ｔｋ−１、ｔｋ、ｔｋ＋１の補正処理の例（３）を説明するための図
である。説明のため、図１７では、実画像（入力画像）ｔｋ−１、ｔｋ、ｔｋ＋１内の特
定の領域を領域Ｒ２とする。図１７に示すように、領域Ｒ２は、実画像ｔｋ−１から実画像ｔｋ＋１へ遷移する間に縮小する。領域Ｒが縮小すると、領域Ｒ２では造影剤などが濃縮され、何も補正処理を実行しなければ、実画像（入力画像）ｔｋ−１、よりも、実画像ｔｋ、ｔｋ＋１の領域Ｒ２の画素値が大きくなってしまう。
【０１１５】
そこで、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋ間、及び実画像ｔｋ、ｔｋ＋１間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、実画像ｔｋについて実画像ｔｋ−１、ｔｋ、ｔｋ＋１で対応する画素値を比較し
、もっともＳ／Ｎの高い画素の画素値を、実画像ｔｋにおける画素値とする。
【０１１６】
また、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋ間、及び実画像ｔｋ、ｔｋ＋１間の動き解析情報に基づき、領域Ｒ２の変形変異を取得し、実画像ｔｋに対して、実画像ｔｋ−１、ｔｋ、ｔｋ＋１で対応する画素値で、Ｓ／Ｎを加味した加重平均を取り、それを新しい画素値とすることも可能である。
【０１１７】
上述のように、補正処理部３０６は、動き解析情報及び２以上の実画像（入力画像）ｔｋ−１、ｔｋの画質情報などに基づき、動き解析情報から得られる入力画像内の特定の領域Ｒ２に関する変形又変位による、画素値の変化を見込んだ補正処理を行うことで、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。
【０１１８】
＜補正処理５：２以上実画像の画質の均一化＞
図１８は、実画像及び補間画像の補正処理の例（４）を説明するための図である。補正処理の他の例として、図１８に示すように、上述した補正処理を組み合わせて、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ−１、ｔｋ間の補間画像tａｕｘ２について、実画像ｔｋ−１、ｔｋに関する動き解析情報、及び画質を均一化する補正処理（図１３、図１４参照）後の実画像ｔｋ−１、ｔｋの画質情報などに基づき、補間画像tａｕｘ２の補正処理を行うことができる。
【０１１９】
補間画像tａｕｘ２の補正処理を行う場合、たとえば、図１６を参照して説明した方法
を用いると、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、ｔｋに関する動き解析情報に基づき、図１８に示す領域Ｒ３の変形変異を取得し、補間画像tａｕｘ２に対して変形変異を見込んだ補正処理を行う。
【０１２０】
また、補間画像tａｕｘ２の補正処理を行う場合、たとえば、図１７を参照して説明し
た方法を用いると、補正処理部３０６は、動き解析情報生成部３０３から得られる実画像ｔｋ−１、補間画像ｔａｕｘ２間、及び補間像ｔａｕｘ２、実画像ｔｋ間の動き解析情報
に基づき、領域Ｒ２の変形変異を取得し、補間実画像ｔａｕｘ２について実画像ｔｋ−１
、ｔｋ、及び補間画像ｔａｕｘ２で対応する画素値を比較し、もっともＳ／Ｎの高い画素
の画素値を、実画像ｔａｕｘ２における画素値とする。この場合、比較する補間画像ｔａ
ｕｘ２の画素値をゼロとすれば良い。
【０１２１】
＜補正処理６：２番目以降の画像の画素値の差分の保存＞
図１９は、実画像の補正処理の例（５）を説明するための図である。補正処理の他の例として、図１９に示すように、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ−１、ｔｋについて、実画像ｔｋ−１、ｔｋに関する動き解析情報、及び実画像ｔｋ−１、ｔｋの画質情報などに基づき、実画像ｔｋの補正処理を行い、実画像ｔｋ−１から数えて２番目の実画像ｔｋについては、実画像ｔｋ−１の画素値との差分のみをデータとして保存しておくで、情報量（量子化ビット数）を減らしつつ、画質を維持することができる。同様に、補正処理部３０６は、補間ボリュームデータ生成部３０５で生成した実画像ｔｋ、ｔｋ＋１について、実画像ｔｋ、ｔｋ＋１に関する動き解析情報、及び実画像ｔｋ、ｔｋ＋１の画質情報などに基づき、実画像ｔｋ＋１の補正処理を行い、実画像ｔｋ＋１については、実画像ｔｋの画素値との差分のみをデータとして保存しておくで、情報量（量子化ビット数）を減らしつつ、画質を維持することができる。
【０１２２】
画像生成部３０７は、補正処理部３０６から出力される補正処理後のボリュームデータｔｋ−１、ｔｋ、とその時刻情報ｔｋ−１、ｔｋ、補正処理後の補間ボリュームデータｔａｕｘ、及びその時刻情報ｔａｕｘを入力とする。そして、画像生成部３０７は、補正処理後の、補間ボリュームデータｔａｕｘ、各ボリュームデータｔｋ−１、ｔｋから実際の
３次元画像を生成し、時系列に並べる。
【０１２３】
表示処理部３０９は、画像生成部３０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部３０７で生成された３次元画像のシーケンスを生成して、表示部３１１で動画像として再生するための処理（再生フレームレートの設定など）を行う ことができる。
【０１２４】
表示部３１１は、表示処理部３０９で適宜処理された、時系列の実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。表示部１１１は、時系列となった実画像ｔｋ−１、ｔ
ｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部１１１は、画像生成部３０７で生成された３次元画像のシーケンスを、表示処理部３０９での処理の後、動画像として再生することができる。
【０１２５】
以上、実施の形態２に係る医療用画像処理装置３００によれば、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正して、２以上の実画像内の特定の領域（たとえば、被検体（患者３）の検査部位（関心領域））の変化を正確に表示することができる。したがって、本実施の形態に係る医療用画像処理装置３００によれば、被検体である患者３の検査部位（関心領域）の観察を精度良く行うことができる。
【０１２６】
例えば、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる。また、医療用画像処理装置３００により、動き解析情報に基づき補間画像を実画像間に補間すれば、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。
【０１２７】
さらに、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。さらに、また、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、複数のボリュームデータに対してセグメンテーション処理を行うときに処理結果を安定させることが出来る。また、臓器間癒着の判断を自動処理で行う場合も医師等が動画を直接観察する場合も容易になる。
【０１２８】
また、本実施の形態に係る医療用画像処理装置３００によれば、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、表示部３１１のノイズ軽減、ノイズのばらつき調整が可能である。
【０１２９】
例えば、医療用画像処理装置３００が生成した動き解析情報によって、実画像又は補間画像の対応画素を特定することにより、本来の詳細な画像情報を失うことなくノイズを軽減できる。また、医療用画像処理装置３００により、ＣＴマルチフェイズ検査において、Ｘ線照射量の多いフェイズと少ないフェイズが混在しても、Ｘ線照射量の少ないフェイズについては、前後のフェイズの画素値を利用することで、良好な補間画像を得られる。
【０１３０】
さらに、医療用画像処理装置３００により、動き解析情報に加えて、２以上の実画像の画質情報に基づき、２以上の実画像の画質又は２つの実画像間を補間する補間画像の画質を補正すれば、重要なフェイズのみＸ線照射量を多くして実画像を生成し、それ以外のフェイズについては、補間画像を生成することで、患者３のＸ線の被ばく量を低減することができる。
【０１３１】
なお、本実施の形態に係る医療用画像処理装置３００では、コンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部３０１に供給されているが、これに限らない。ボリュームデータ生成部３０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
【０１３２】
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。
【０１３３】
なお、本実施の形態において、表示処理部３０９が、画像生成部３０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部３１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
【０１３４】
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻に基づき、生成しているが、これに限らない。
【０１３５】
なお、本実施の形態の説明に用いた各機能ブロック（図１２参照）は、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。
【０１３６】
なお、本実施の形態では、心臓の観察に用いる例を示したが、これに限らない。例えば、腫瘍組織の経時変化の観察、血流灌流の観察などに用いることも出来る。
【０１３７】
（実施の形態３）
上述した実施の形態１に係る医療用画像処理装置１００によれば、動き解析情報に基づき、補間画像を実画像間に補間することで、画像内の特定の領域の変化を正確に表示する。その一方で、観察対象の領域を画像としてよりも、その形状や大きさなどのより客観的な指標で把握したい場合がある。この場合は観察対象のジオメトリ情報を作成するが、補間ボリュームを含む多数のボリュームデータそれぞれにジオメトリ情報を作成するには多大な労力がかかる。また、異なるボリュームデータのジオメトリ情報の間の対応関係を設定するのも多大な労力がかかる。
【０１３８】
そこで、実施の形態３に係る医療用画像処理装置５００では、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成することで、個別にジオメトリ情報を生成するのに比較して、セグメンテーションに係る操作の手間、補正処理を減らすことができる。さらに、実施の形態３に係る医療用画像処理装置５００では、実画像と他の実画像との画質に差異によって、セグメンテーションの結果の揺らぎを抑止できる。そのため、ユーザは、より客観的な画像解析を行うことができる。
【０１３９】
さらには、実施の形態３に係る医療用画像処理装置５００では、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報から他の実画像のジオメトリ情報を生成しているので、他の実画像のジオメトリ情報は、内在的に実画像のボリュームデータから生成したジオメトリ情報との対応関係を有している。そのため、他の実画像のジオメトリ情報と実画像のジオメトリ情報との対応関係を改めて取得、計算する必要がない。
【０１４０】
図２０を参照して、本実施の形態３に係る医療用画像処理装置５００の構成について説明する。図２０に示す医療用画像処理装置５００は、ボリュームデータ生成部５０１と、動き解析情報生成部５０３と、補間ボリュームデータ生成部５０５と、ジオメトリ情報生成部５０６と、画像生成部５０７と、表示処理部５０９と、を備える。さらに、医療用画像処理装置５００には、医療用画像処理装置５００で生成された画像を表示する表示部５１１と、医療用画像処理装置５００の各部を操作するためのユーザインターフェースである操作部５１３と、が設けられている。
【０１４１】
ボリュームデータ生成部５０１は、実施の形態１と同様、コンピュータ断層撮影装置により供給された断層信号から、撮影時刻の異なる複数のボリュームデータを生成し、各ボリュームデータが生成された時刻を表す時刻情報と共に、動き解析情報生成部５０３と、ジオメトリ情報生成部５０６と、画像生成部５０７とへ出力する。
【０１４２】
また、ボリュームデータ生成部５０１は、コンピュータ断層撮影装置から、装置の撮像条件（例えば、Ｘ線量、撮影時間、造影剤の注入の時刻）を取得し、画像生成部５０７へ出力する。なお、コンピュータ断層撮影装置の撮像条件以外にも、コンピュータ断層撮影装置の撮影対象の条件（患者の動きや、患者の体内に含まれる金属など）により、のちに生成される実画像の画質が大きく変化する。
【０１４３】
なお、ボリュームデータ生成部５０１は、図示しないメモリに、撮影時刻の異なる複数のボリュームデータ、コンピュータ断層撮影装置の撮像条件、コンピュータ断層撮影装置の撮影対象の条件を保持することができる。
【０１４４】
ここで、本実施の形態では、ボリュームデータ生成部５０１は、実施の形態１と同様、複数のボリュームデータの各ボリュームデータが、コンピュータ断層撮影装置により撮影された時刻ｔに基づき、ボリュームデータを生成した時刻情報ｔｎ（ｎ＝１、２、…、ｋ−１、ｋ、ｋ＋１…：自然数）を生成するものとする。
【０１４５】
以下、説明のため、ボリュームデータ生成部５０１では、時刻情報ｔｋ−１のボリュームデータと、時刻ｔｋ−１から所定時間経過後に撮影され、ボリュームデータ生成部５０１に供給された断層信号に基づき生成された時刻情報ｔｋのボリュームデータと、時刻ｔｋから所定時間経過後に撮影され、ボリュームデータ生成部５０１に供給された断層信号に基づき生成された時刻情報ｔｋ＋１のボリュームデータを保持しているものとする。また、以下、説明のため、時刻情報ｔｋのボリュームデータを、ボリュームデータｔｋと適宜省略して記載する。
【０１４６】
動き解析情報生成部５０３は、ボリュームデータ生成部５０１から出力されたボリュームデータのうち、ボリュームデータｔｋ−１及びその時刻情報ｔｋ−１、並びにボリュームデータｔｋ及びその時刻情報ｔｋから、動き解析情報を生成する。
【０１４７】
ここで、動き解析情報とは、複数のボリュームデータ上の対応する位置もしくは対応する物体の対応関係の情報、及び、前記位置及び物体が移動変化する過程の情報を指す。また、各ボリュームデータの画素が、時刻ｋ−１と時刻ｋとの間の任意の時刻での位置を示す指標となる。
【０１４８】
そして、動き解析情報生成部５０３は、生成した動き解析情報を、補間ボリュームデータ生成部５０５と、ジオメトリ情報生成部５０６と、画像生成部５０７とに出力する。動き解析情報生成部５０３が生成する動き解析情報は、実施の形態１の動き解析情報生成部１０３と同様（図２参照）、３次元画像の情報を含むボリュームデータｔｋ−１の各画素がどのように変位して、ボリュームデータｔｋでの画素となっているかを示す。
【０１４９】
補間ボリュームデータ生成部５０５は、操作部５１３の入力指示に基づき、必要に応じて、３次元画像の情報を含むボリュームデータｔｋ−１とボリュームデータｔｋとの間を補間する、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘとを生成し、画像生成部５０７とに出力する。
【０１５０】
ジオメトリ情報生成部５０６は、ボリュームデータ生成部５０１から出力された、撮影時刻の異なる複数のボリュームデータとその時刻情報に基づき、必要に応じて、１又は複数のボリュームデータのジオメトリ情報を生成し、画像生成部５０７に出力する。以下、時刻情報ｔｋにおけるボリュームデータｔｋのジオメトリ情報を、ジオメトリ情報Ｇｔｋと称する。
【０１５１】
ここで、ジオメトリ情報とは、１つの３次元実画像について、中実、表面、線、位置、又はこれらの要素の組み合わせで表現した組織の情報である。特に、ボクセルデータのマスクで表現される領域の中実や表面であることが多く、セグメンテーションなどの従来手法を用いることができる。また、トラッキングなどの従来手法を用いて取得される線であることもある。これらはマスク・パス・点として表現することができる。
【０１５２】
図２１（ａ）、図２１（ｂ）を参照して、ジオメトリ情報生成部５０６での処理について説明する。図２１（ａ）は、ジオメトリ情報生成部５０６での処理１を模式的に示す図であり、図２１（ｂ）は、ジオメトリ情報生成部５０６での処理２を模式的に示す図である。
【０１５３】
図２１（ａ）に示すように、ジオメトリ情報生成部５０６は、ボリュームデータｔｋ−１における、３次元的に表現されている特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１を生成する。次に、図２１（ｂ）に示すように、ジオメトリ情報生成部５０６は、特定の領域Ｒのジオメトリ情報Ｇｔｋ−１及び、あらかじめ動き解析情報生成部５０３から取得したボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報から、特定の領域Ｒのジオメトリ情報Ｇｔｋを生成する。ジオメトリ情報Ｇｔｋはジオメトリ情報Ｇｔｋ−１に含まれる座標情報をその座標における動き解析情報による動きを適用して、ジオメトリ情報Ｇｔｋにおける座標を取得し、この時に座標情報(例えばノード)の属するマスク・パス・点の関係(例えばリンク)を維持することによって求めることが出来る。
【０１５４】
また、ジオメトリ情報生成部５０６は、操作部５１３の入力指示に基づき、必要に応じて、生成した各ボリュームデータのジオメトリ情報を、補間ボリュームデータ生成部５０５に出力する。
【０１５５】
また、補間ボリュームデータ生成部５０５は、操作部５１３の入力指示に基づき、必要に応じて、動き解析情報生成部５０３で生成した動き解析情報、及びジオメトリ情報生成部５０６で生成したジオメトリ情報Ｇｔｋに基づき、補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘを生成することができる。
【０１５６】
ここで、本実施の形態では、実施の形態１と同様、補間ボリュームデータｔａｕｘが有する時刻情報ｔａｕｘは、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報とは限らない。言い換えると、補間ボリュームデータの時刻情報ｔａｕｘと直前のタイミングに生成されたボリュームデータｔｋ−１の時刻情報ｔｋ−１との間隔Ａと、補間ボリュームデータの時刻情報ｔａｕｘと直後のタイミングに生成されたボリュームデータの時刻情報ｔｋとの間隔Ｂとが異なる間隔、すなわち、非等間隔となっている。この理由については、実施の形態１で上述した理由と同じであるため、詳細な説明は省略する。
【０１５７】
画像生成部５０７は、ボリュームデータ生成部５０１から出力されるボリュームデータｔｋ−１、ｔｋ、とその時刻情報ｔｋ−１、ｔｋ、補間ボリュームデータ生成部５０５から出力される補間ボリュームデータｔａｕｘとその時刻情報ｔａｕｘ、を入力とする。そして、画像生成部５０７は、補間ボリュームデータｔａｕｘ、各ボリュームデータｔｋ−１、ｔｋから実際の３次元画像、及び、補間ジオメトリ情報Ｇｔａｕｘ、各ジオメトリ情報Ｇｔｋと重ね合わせた画像を生成し、時系列に並べる。この時に、補間ジオメトリ情報Ｇｔａｕｘ、各ジオメトリ情報Ｇｔｋに含まれる数値を可視化する。数値は例えば、組織の体積(心室の体積、肺の体積、腫瘍の体積)、表面積(腫瘍が浸潤した面の面積)、長さ(血管の長さ、心壁の厚さ)、位置(弁の座標)が挙げられる。また、補間ジオメトリ情報Ｇｔａｕｘ、各ジオメトリ情報Ｇｔｋに含まれる数値の変化を可視化する。例えば、肺活量、腫瘍の成長などがある。
【０１５８】
表示処理部５０９は、画像生成部５０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部３１１のディスプレイに表示するための各種処理を行う。例えば、画像生成部５０７で生成された３次元画像のシーケンスを生成して、表示部５１１で動画像として再生するための処理（再生フレームレートの設定など）を行うことができる。
【０１５９】
表示部５１１は、表示処理部５０９で適宜処理された、時系列の実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを表示する。表示部５１１は、時系列となった実画像ｔｋ−１、ｔ
ｋ及び補間画像ｔａｕｘ以外に、各画像の時刻情報も表示する。また、必要に応じて、表示部５１１は、画像生成部５０７で生成された３次元画像のシーケンスを、表示処理部５０９での処理の後、動画像として再生することができる。
【０１６０】
＜画像処理の例（１）＞
次に、図２２を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２２は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（１）を示す図である。
【０１６１】
図２２に示すように、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集をすると、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報に基づき、インタラクティブに、その編集内容を他の３次元画像ｔｋに伝播させることはできる。これは、あらかじめ、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報と、ジオメトリ情報Ｇｔｋ−１とにより、他の３次元画像ｔｋの特定の領域Ｒについて、ジオメトリ情報Ｇｔを生成しているからである。また、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集・保存したものに、後から、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報を適用し、その編集内容を他の３次元画像ｔｋに伝播させることができる。これによってジオメトリ情報Ｇｔｋ−１作成の処理やユーザの労力が軽減される。
【０１６２】
また、ジオメトリ情報Ｇｔに含まれるマスク・パス・点などの編集をすると、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報を修正する事が出来る。これは、ジオメトリ情報Ｇｔとジオメトリ情報Ｇｔｋ−１との対応関係が定められているからである。これによって動き解析情報を容易に修正することが出来る。
【０１６３】
＜画像処理の例（２）＞
次に、図２３を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２３は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（２）を示す図である。
【０１６４】
図２３では、まず、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集をすると、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報に基づき、インタラクティブに、その編集内容を他の３次元画像ｔｋに伝播させることができる。さらに、ユーザが、操作部５１３を介して、３次元実画像ｔｋ−１上で、特定の領域Ｒについて、ジオメトリ情報Ｇｔｋ−１に含まれるマスク・パス・点などの編集した後に、同じ動き解析情報を用いて、補間ジオメトリ情報Ｇｔａｕｘ、及び補間画像ｔａｕｘ３を生成成することができる。このとき、補間画像ｔａｕｘ３の時刻情報ｔａｕｘ３は、時刻情報ｔｋ−１と時刻情報ｔｋとの中間の時刻情報ではない、非等間隔な時刻で生成することできる。
なお、補間ジオメトリ情報Ｇｔａｕｘとは、ジオメトリ情報と同様、１つの３次元補間画像について、中実、表面、線、位置、又はこれらの要素の組み合わせで表現した組織の情報である。
【０１６５】
そして、補間画像ｔａｕｘ３の時刻情報ｔａｕｘ３を決める際に、あらかじめ生成していたボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報と、ジオメトリ情報Ｇｔｋ−１を利用する。たとえば、特定の領域Ｒが心臓である場合、その容積が最大となる時刻を、補間画像ｔａｕｘ３の時刻情報ｔａｕｘ３に設定すれば、特定の領域Ｒの時間変化を容易に画像化することができる。
【０１６６】
ここで、補間ボリュームデータとは、時刻ｔｋ−１と時刻ｔｋとの間の特定の時刻ｔａｕｘにおけるボリュームデータである。つまり、補間ボリュームデータ生成部５０５が補間ボリュームデータｔａｕｘを生成することは、実際はコンピュータ断層撮影装置で撮像されていない時刻ｔａｕｘにおいて、あたかもコンピュータ断層撮影装置により撮像して、ボリュームデータ生成部５０１がボリュームデータを生成したことと、同等となる。
【０１６７】
＜画像処理の例（３）＞
次に、図２４を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２４は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（３）を示す図である。
【０１６８】
図２４に示すように、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報から、特定の領域Ｒについて、時刻情報ｔｋ−１から時刻情報ｔｋまでの変形変異を取得し、この変形変異を見込んで、ボリュームデータｔｋ−１について特定の領域Ｒのボクセル値の変化を予測し補正する補正処理を行う際に、ジオメトリ情報Ｇｔｋ−１を利用することができる。つまり、図２４に示すように、特定の領域Ｒは、実画像ｔｋ−１から実画像ｔｋへ遷移する間に縮小する。特定の領域Ｒが縮小すると、特定の領域Ｒでは造影剤などが濃縮され、何も補正処理を実行しなければ、実画像ｔｋ上、特定の領域Ｒの画素値が大きくなってしまう。これは、造影剤の注入の時刻情報を加味すると更に効果的である。
【０１６９】
そこで、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報から、領域Ｒの変形変異を取得し、実画像ｔｋ−１に対して、ジオメトリ情報Ｇｔｋ−１に基づき、変形変異を見込んだ補正処理を行う。
【０１７０】
なお、上述した変形変異を見込んだ補正処理は、特定の領域Ｒに含まれる組織によって、その補正の度合を適宜変更するのが良い。補正の度合とは、例えば、ＣＴ検査における、「臓器の収縮」、「近接している高信号の部分の影響」であり、ＭＲ検査における、「位置の変化」、超音波検査における、「速度の変化」などがある。
【０１７１】
＜画像処理の例（４）＞
次に、図２５を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２５は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（４）を示す図である。
【０１７２】
図２５に示すように、ボリュームレンダリング法（ＶＲ）、最大値投影法（ＭＩＰ）、ＭＰＲ（Ｍｕｌｔｉ−Ｐｌａｎａｒ Ｒｅｎｄｅｒｉｎｇ）、ＣＰＲ（Ｃｕｒｖｅｄ Ｐｌａｎａｒ Ｒｅｆｏｒｍａｔｉｏｎ）、仮想内視鏡法（ＶＥ）など、医療用３ＤＷＳの既存の種々のレンダリングでの実画像ｔｋ−１において、特定の領域Ｒの注目点を画面上で固定する。そして、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報、及び実画像ｔｋ−１における特定の領域Ｒのジオメトリ情報Ｇｔｋ−１を利用することで、特定の領域Ｒ内のある領域又は点に着目して、表示部５１１上でなるべく動かないように、視点の位置・方向・回転角度・拡大率を、実画像毎に設定することができる。
【０１７３】
＜画像処理の例（５）＞
次に、図２６を参照して、本実施の形態３に係る医療用画像処理装置５００の画像処理の一例について説明する。図２６は、本実施の形態３に係る医療用画像処理装置５００の画像処理の例（５）を示す図である。
【０１７４】
図２６に示すように、ボリュームデータｔｋ−１とボリュームデータｔｋ間の動き解析情報、及び実画像ｔｋ−１における特定の領域Ｒのジオメトリ情報Ｇｔｋ−１を利用することで、仮想内視鏡法（ＶＥ）で得られる３次元実画像ｔｋ−１上で、視点を特定の領域Ｒ内で固定し、その視点の位置及び方向を、動き解析情報に従い、各３次元実画像ごとに設定することができる。
【０１７５】
以上、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成することで、観察対象の形状や大きさなどのより客観的な指標で把握することができる。
【０１７６】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成することで、個別にジオメトリ情報を生成するのに比較して、
セグメンテーションに係る操作の手間、補正処理を減らすことができる。さらに、実画像と他の実画像との画質に差異によって、セグメンテーションの結果の揺らぎを抑止できる。そのため、ユーザは、より客観的な画像解析を行うことができる。
【０１７７】
さらには、本実施の形態に係る医療用画像処理装置５００では、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報から他の実画像のジオメトリ情報を生成しているので、他の実画像のジオメトリ情報は、内在的に実画像のボリュームデータから生成したジオメトリ情報との対応関係を有している。そのため、他の実画像のジオメトリ情報と実画像のジオメトリ情報との対応関係を改めて取得、計算する必要がない。
【０１７８】
本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、従来のレンダリング法でシネ再生するだけでも、心機能の不全等を良好に観察できる
【０１７９】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、心臓マルチフェイズデータのような短い時間間隔のものだけでなく、病変の経年変化を観察するような長い時間間隔のデータにも応用できる。
【０１８０】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、ストレイン解析、パフュージョン解析、心機能解析などで、再現性のある全自動化された定量化により、従来よりも安価に正確で結果を誰もが得られる。
【０１８１】
また、本実施の形態に係る医療用画像処理装置５００によれば、動き解析情報に加えて、実画像のボリュームデータから生成したジオメトリ情報に基づき、他の実画像のジオメトリ情報を生成するので、セグメンテーションにより隣接する臓器を設定し、動き解析情報により動きが一緒の部分を推定、癒着の可能性あり、として表示部５１１に表示することができる。
【０１８２】
なお、本実施の形態に係る医療用画像処理装置５００では、コンピュータ断層撮影装置により、患者３の診断範囲における多数の相連続する断層信号がボリュームデータ生成部５０１に供給されているが、これに限らない。ボリュームデータ生成部５０１に供給される信号は、コンピュータ断層撮影装置以外に、たとえば、ＭＲＩ装置（Ｍａｇｎｅｔｉｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）といった、被検体のボリュームデータから検査部位の断層像を生成し、モニター画面に表示して診断を行う装置からの断層信号も含まれる。
【０１８３】
なお、本実施の形態において、ボリュームデータｔｋの時刻ｔｋは、ボリュームデータｔｋ−１から所定時間経過した後の時刻であるが、この所定時間は任意の時間であり、たとえば、被検体（患者３）の検査部位（関心領域）の変化が現れると予測される時間に依存して、適宜設定することができる。
【０１８４】
なお、本実施の形態において、表示処理部５０９が、画像生成部５０７で時系列となった実画像ｔｋ−１、ｔｋ及び補間画像ｔａｕｘを、表示部５１１のディスプレイに表示するための各種処理を行う際、以下のような処理が表示部５１１で可能なように設定されても良い。つまり、（１）動き解析、マスク・パス等の編集、補間画像の生成は、時間・場所・処理装置に関して別々に行え、（２）動き解析情報と、マスク・パス等の編集状態（ワークスペース）を、別のファイルとして保存し、後で組み合わせて補間画像を生成する、（３）動き解析情報と、マスク・パス・レンダリングの情報を含む画像ファイルを、後で組み合わせて補間（増量）画像を生成する、（４）補間画像の生成やジオメトリ情報の伝播の機能があるワークステーション上に、動き解析情報の有無を表示し、ユーザに動き解析が完了したことを即時に知らせる、（５）異なる精度、アルゴリズム、による複数の動き解析情報のマネージメント（削除・コピー・移動）ができる。
【０１８５】
なお、本実施の形態では、ボリュームデータを生成した時刻情報を、コンピュータ断層撮影装置により供給された時刻に基づき、生成しているが、これに限らない。
【０１８６】
なお、本実施の形態の説明に用いた各機能ブロックは、典型的には集積回路であるＬＳＩとして実現される。特に、実画像及び補間画像を生成する計算処理は、ＧＰＵ（Ｇｒａｐｈｉｃ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）により行うことができる。ＧＰＵは、汎用のＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）と比較して画像処理に特化した設計なされている演算処理装置で、通常のＣＰＵとは別個にコンピュータに搭載される。
【符号の説明】
【０１８７】
１００、３００、５００ 医療用画像処理装置
１０１、３０１、５０１ ボリュームデータ生成部
１０３、３０３、５０３ 動き解析情報生成部
１０５、３０５、５０５ 補間ボリュームデータ生成部
１０７、３０７、５０７ 画像生成部
１０９、３０９、５０９ 表示処理部
１１１、３１１、５１１ 表示部
１１３、３１３、５１３ 操作部
３０６ 補正処理部
５０６ ジオメトリ情報生成部
【特許請求の範囲】
【請求項１】
時刻を有する複数のボリュームデータのそれぞれを可視化する医療用画像処理装置であって、
前記複数のボリュームデータのうち、一のボリュームデータについて、ジオメトリ情報を生成する第１生成手段と、
前記複数のボリュームデータのうち、前記一のボリュームデータと他の少なくとも一のボリュームデータとの間の動き解析情報を生成する第２生成手段と、
前記複数のボリュームデータのうち、前記他の少なくとも一のボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する新たなジオメトリ情報を生成する第３生成手段と、
前記他の少なくとも一のボリュームデータについて、前記対応する新たなジオメトリ情報を含む画像を表示するための処理を行う表示処理手段と、を備える医療用画像処理装置。
【請求項２】
請求項１に記載の医療用画像処理装置であって、
前記動き解析情報に基づき、前記一のボリュームデータに含まれる特定の領域の変異又は変形を取得し、当該取得した変異又は変形、前記ジオメトリ情報、及び前記対応する新たなジオメトリ情報に基づき、ボクセル値の変化を予測し、前記一のボリュームデータのボクセル値を補正する、医療用画像処理装置。
【請求項３】
請求項１に記載の医療用画像処理装置であって、
さらに、前記一のボリュームデータと前記他の少なくとも一のボリュームデータとの間を補間する補間ボリュームデータを生成する第４生成手段と、
前記補間ボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する補間ジオメトリ情報を生成する第５生成手段と、を備え、
前記表示処理手段は、
前記補間ボリュームデータについて、前記補間ジオメトリ情報を含む画像を表示するための処理を行う、医療用画像処理装置。
【請求項４】
請求項３に記載の医療用画像処理装置であって、
前記一のボリュームデータの時刻情報及び前記補間ボリュームデータの時刻情報の間隔と、前記他のボリュームデータ及び前記補間ボリュームデータの時刻情報の間隔とは非等間隔である、医療用画像処理装置。
【請求項５】
請求項１に記載の医療用画像処理装置であって、
前記表示処理手段は、
前記複数のボリュームデータについて描画する時に、前記ジオメトリ情報と前記新たなジオメトリ情報とを利用して、前記複数のボリュームデータで共通する特定の領域を画像として表示するための処理を行う、医療用画像処理装置。
【請求項６】
コンピュータを、請求項１から請求項５のいずれか１項に記載の各手段として機能させるための医療用画像処理プログラム。
【請求項１】
時刻を有する複数のボリュームデータのそれぞれを可視化する医療用画像処理装置であって、
前記複数のボリュームデータのうち、一のボリュームデータについて、ジオメトリ情報を生成する第１生成手段と、
前記複数のボリュームデータのうち、前記一のボリュームデータと他の少なくとも一のボリュームデータとの間の動き解析情報を生成する第２生成手段と、
前記複数のボリュームデータのうち、前記他の少なくとも一のボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する新たなジオメトリ情報を生成する第３生成手段と、
前記他の少なくとも一のボリュームデータについて、前記対応する新たなジオメトリ情報を含む画像を表示するための処理を行う表示処理手段と、を備える医療用画像処理装置。
【請求項２】
請求項１に記載の医療用画像処理装置であって、
前記動き解析情報に基づき、前記一のボリュームデータに含まれる特定の領域の変異又は変形を取得し、当該取得した変異又は変形、前記ジオメトリ情報、及び前記対応する新たなジオメトリ情報に基づき、ボクセル値の変化を予測し、前記一のボリュームデータのボクセル値を補正する、医療用画像処理装置。
【請求項３】
請求項１に記載の医療用画像処理装置であって、
さらに、前記一のボリュームデータと前記他の少なくとも一のボリュームデータとの間を補間する補間ボリュームデータを生成する第４生成手段と、
前記補間ボリュームデータについて、前記動き解析情報及び前記ジオメトリ情報に基づき、対応する補間ジオメトリ情報を生成する第５生成手段と、を備え、
前記表示処理手段は、
前記補間ボリュームデータについて、前記補間ジオメトリ情報を含む画像を表示するための処理を行う、医療用画像処理装置。
【請求項４】
請求項３に記載の医療用画像処理装置であって、
前記一のボリュームデータの時刻情報及び前記補間ボリュームデータの時刻情報の間隔と、前記他のボリュームデータ及び前記補間ボリュームデータの時刻情報の間隔とは非等間隔である、医療用画像処理装置。
【請求項５】
請求項１に記載の医療用画像処理装置であって、
前記表示処理手段は、
前記複数のボリュームデータについて描画する時に、前記ジオメトリ情報と前記新たなジオメトリ情報とを利用して、前記複数のボリュームデータで共通する特定の領域を画像として表示するための処理を行う、医療用画像処理装置。
【請求項６】
コンピュータを、請求項１から請求項５のいずれか１項に記載の各手段として機能させるための医療用画像処理プログラム。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【図１９】
【図２０】
【図２１】
【図２２】
【図２３】
【図２４】
【図２５】
【図２６】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図９】
【図１０】
【図１１】
【図１２】
【図１３】
【図１４】
【図１５】
【図１６】
【図１７】
【図１８】
【図１９】
【図２０】
【図２１】
【図２２】
【図２３】
【図２４】
【図２５】
【図２６】
【公開番号】特開２０１１−１１０２８２（Ｐ２０１１−１１０２８２Ａ）
【公開日】平成２３年６月９日（２０１１．６．９）
【国際特許分類】
【出願番号】特願２００９−２７０４８１（Ｐ２００９−２７０４８１）
【出願日】平成２１年１１月２７日（２００９．１１．２７）
【出願人】（５００１０９３２０）ザイオソフト株式会社 (59)
【Ｆターム（参考）】
【公開日】平成２３年６月９日（２０１１．６．９）
【国際特許分類】
【出願日】平成２１年１１月２７日（２００９．１１．２７）
【出願人】（５００１０９３２０）ザイオソフト株式会社 (59)
【Ｆターム（参考）】
[ Back to top ]