"三次元画像のボリュームデータに対する画像処理を実行する画像処理装置であって、前記ボリュームデータから予め与えられた条件を満たす三次元の初期領域を抽出する抽出処理実行部と、前記初期領域に修正処理を施して三次元の修正領域を抽出する領域修正処理実行部と、前記ボリュームデータから前記三次元画像の複数の断面図を作成して当該複数の断面図の少なくとも一部を出力する可視化処理実行部とを有しており、前記可視化処理実行部は、前記初期領域を表示する場合には、当該初期領域に含まれるボクセルを有する一以上の断面図を当該可視化処理実行部が作成した複数の断面図から選択して、選択された一以上の断面図を当該初期領域に含まれるボクセルとその他の領域とを区別可能に出力し、前記修正領域を表示する場合には、当該修正領域に含まれるボクセルを有する一以上の断面図を当該可視化処理実行部が作成した複数の断面図から選択して、選択された一以上の断面図を当該修正領域に含まれるボクセルとその他の領域とを区別可能に出力することを特徴とする画像処理装置。請求項1記載の画像処理装置であって、前記可視化処理実行部は、前記初期領域を表示した後に前記修正領域を表示する場合に、当該期領域に含まれるボクセルを有する断面図の数と当該修正領域に含まれるボクセルを有する断面図の数が異なる場合には、出力する断面図の数を変更することを特徴とする画像処理装置。請求項1記載の画像処理装置であって、前記可視化処理実行部は、前記初期領域を表示した後に前記修正領域を表示する場合に、前記初期領域と前記修正領域の三次元領域におけるサイズが異なる場合には、出力する断面図の数を変更することを特徴とする画像処理装置。請求項1乃至3記載の画像処理装置であって、前記領域修正処理実行部は、前記条件として与えられるパラメータを変更することにより前記修正領域を抽出することを特徴とする画像処理装置。請求項1乃至3記載の画像処理装置であって、前記領域修正処理実行部は、予め定められたルールに沿って領域が変化するように前記修正領域を抽出することを特徴とする画像処理装置。請求項1乃至5記載の画像処理装置であって、前記領域修正処理実行部は、前記初期領域の輪郭を補正することによって前記修正領域を抽出することを特徴とする画像処理装置。画像処理装置が三次元画像のボリュームデータに対する画像処理を実行する方法であって、前記画像処理装置の抽出処理実行部が、前記ボリュームデータから予め与えられた条件を満たす三次元の初期領域を抽出するステップと、前記画像処理装置の領域修正処理実行部が、前記初期領域に修正処理を施して三次元の修正領域を抽出するステップと、前記画像処理装置の可視化処理実行部が、前記ボリュームデータから前記三次元画像の複数の断面図を作成するステップと、前記可視化処理実行部が、前記初期領域に含まれるボクセルを有する一以上の断面図を当該可視化処理実行部が作成した複数の断面図から選択して、選択された一以上の断面図を当該初期領域に含まれるボクセルとその他の領域とを区別可能に出力するステップと、前記可視化処理実行部が、前記修正領域に含まれるボクセルを有する一以上の断面図を当該可視化処理実行部が作成した複数の断面図から選択して、選択された一以上の断面図を当該修正領域に含まれるボクセルとその他の領域とを区別可能に出力するステップとを有することを特徴とする方法。請求項7記載の方法であって、前記可視化処理実行部は、前記初期領域を表示した後に前記修正領域を表示する場合に、当該期領域に含まれるボクセルを有する断面図の数と当該修正領域に含まれるボクセルを有する断面図の数が異なる場合には、出力する断面図の数を変更することを特徴とする方法。請求項7記載の方法であって、前記可視化処理実行部は、前記初期領域を表示した後に前記修正領域を表示する場合に、前記初期領域と前記修正領域の三次元領域におけるサイズが異なる場合には、出力する断面図の数を変更することを特徴とする方法。請求項7乃至9記載の方法であって、前記領域修正処理実行部は、前記条件として与えられるパラメータを変更することにより前記修正領域を抽出することを特徴とする方法。請求項7乃至9記載の方法であって、前記領域修正処理実行部は、予め定められたルールに沿って領域が変化するように前記修正領域を抽出することを特徴とする方法。請求項7乃至11記載の方法であって、前記領域修正処理実行部は、前記初期領域の輪郭を補正することによって前記修正領域を抽出することを特徴とする方法。自動で抽出された領域にたいして手動で修正を加える場合にそれを簡便に行う方法としてこれまでに提案されている技術には、輪郭線上に修正したい点が複数ある場合であっても、一箇所ずつしか修正できないことや、自由曲線として修正するため、修正後の曲線輪郭線が元の画像情報から大きく異なってしまう可能性があること、また変化の経過を保存する場合には、修正点の量によってはデータ量が大きくなってしまうということや、抽出目的の領域以外にガイドとなる領域の設定をしなければならないことなどの問題点がある。そこで、三次元領域に対する手動修正を行う際に、できるだけ簡便に修正を行うことが可能な画像処理方法、画像処理装置、及びプログラムが要求されている。三次元画像のボリュームデータに対する画像処理を実行する画像処理装置であって、ボリュームデータから予め与えられた条件を満たす三次元の初期領域を抽出する抽出処理実行部と、初期領域に修正処理を施して三次元の修正領域を抽出する領域修正処理実行部と、ボリュームデータから三次元画像の複数の断面図を作成して当該複数の断面図の少なくとも一部を出力する可視化処理実行部とを有する。可視化処理実行部は、初期領域を表示する場合には、当該初期領域に含まれるボクセルを有する一以上の断面図を当該可視化処理実行部が作成した複数の断面図から選択して、選択された一以上の断面図を当該初期領域に含まれるボクセルとその他の領域とを区別可能に出力し、修正領域を表示する場合には、当該修正領域に含まれるボクセルを有する一以上の断面図を当該可視化処理実行部が作成した複数の断面図から選択して、選択された一以上の断面図を当該修正領域に含まれるボクセルとその他の領域とを区別可能に出力する。3次元ボリュームデータに対して自動で抽出された3次元関心領域に対して、簡便にかつ短時間に修正を行うことができる。以下に、画像処理装置、画像処理方法およびプログラムについて、図を参照しながら詳細に説明する。図1は、画像処理装置の構成例を示す図である。図1に示すように、システムは、三次元画像以下ボリュームデータとも呼ぶやそれに関わるデータを記憶する外部記憶装置10、X線CT装置やMRI装置等で撮像され再構成されたボリュームデータに対して画像処理を行う画像処理装置11、画像処理された画像を表示する表示装置12、各処理の開始等の指示や、画面や抽出結果に対する操作を入力する入力装置13を備える。画像処理装置11は、ボリュームデータやそれに関わる情報を保存する内部メモリ20、画像処理アルゴリズム実行部21及び入力操作取得部22を備える。画像処理アルゴリズム実行部21は画像処理装置11に備えるCPUCentralProcessingUnitによるプログラム処理や専用回路により実現される。さらに、内部メモリ20、及び画像処理装置11の機能を実現するためのプログラムを格納する記憶部図示略は、RAM、ROM、HDD、フラッシュメモリ等の記憶媒体により構成される。画像処理アルゴリズム実行部21は、ボリュームデータに対して可視化処理を実行し表示装置に表示する可視化処理実行部31、抽出された領域の情報を格納する抽出領域格納部33、自動抽出された領域に対する修正処理を行う領域修正処理実行部32、ボリュームデータに対して自動抽出処理を行う自動抽出処理実行部34を備える。図2を用い、自動抽出処理結果を修正する際の全体の流れを説明する。まず、入力装置13からの入力、もしくはシステムからの指示で、三次元自動抽出領域修正処理を開始するステップ40。外部記憶装置10から内部メモリ20に転送されるボリュームデータが画像処理アルゴリズム実行部21に入力され、自動抽出処理実行部34は、入力ボリュームデータに対して自動抽出処理を実行するステップ41。自動抽出処理実行部34で行われる自動抽出処理とは例えば、領域拡張法や、LevelSet法などの処理を指す。ここでいう領域拡張法とは、自動または手動で指定した1ボクセル以上のシード点から、隣接ボクセルが拡張条件に適合するかどうかを判断しながら領域を拡張させてボリュームデータから所望の三次元領域を抽出する手法であり、LevelSetとは、二次元空間内に存在するオブジェクトの境界面を、三次元空間内で定義される曲面を時間発展させることにより検出することで、ボリュームデータから所望の三次元領域を抽出する手法である。抽出領域格納部33は、ステップ41で自動抽出された三次元領域抽出領域とも呼ぶをまずは初期領域として格納するステップ42。可視化処理実行部31は、抽出領域格納部33に格納された初期領域に対して可視化処理を実行するステップ43。可視化処理実行部31では例えば、Axial、Sagittal、Coronal等の規定された断面や任意断面での断面図を表示するMultiPlanarReformatMPR、任意曲面での断面図を表示するCurevedMultiPlanarReformatCMPR等の三次元可視化手法を入力ボリュームデータに対して用いる。この結果入力ボリュームデータによって現される三次元画像三次元領域の複数の断面図が作成される。この際、入力ボリュームデータのうち、自動抽出処理実行部34で抽出された三次元領域と、そうでない領域とは別の色付けを行うなどし、抽出された三次元領域を視認可能とする。抽出された三次元領域を正確に確認するためには、断面以下スライスともいう上のどのボクセルが初期領域として抽出されているかを確認することが重要である。そこでステップ43では、二次元平面上で三次元形状を正確に把握するため、初期領域として設定されるボクセル全てを二次元平面上に表示する方法をとる。具体的には例えば、入力ボリュームデータから複数のスライスを作成した際に抽出領域に含まれるボクセルを有するスライスを可視化処理実行部31が記憶しておき、この記憶しておいたスライスのデータ表示装置12に出力して表示させる方法がある。この際、表示されるスライス内で抽出領域として設定されたボクセルには色付けを行っておくと良い。また、CMPRで作成される任意に設定した平面に平行な複数の断面のうち、抽出領域として設定されるボクセルが含まれる断面を可視化処理実行部31が記憶しておき、記憶しておいた断面を表示装置12に表示させる方法でもよい。このように表示された状態で、抽出されている領域を確認しステップ44、初期領域が関心領域を過不足なく抽出している場合は、この時点で抽出処理の終了となるステップ47。過不足がある場合には、表示されている領域例えば初期領域の修正処理を行うステップ45が、ここで用いられる修正処理とは例えば、先に述べた自動抽出処理のパラメータを変更する方法や、モルフォロジカル処理や輪郭補正法による輪郭の修正等が挙げられる。ただしここでは、ユーザが自由に輪郭を修正するのではなく、初期領域または入力ボリュームデータの輝度情報にある程度忠実な修正を加えたものを新たな抽出領域とする。この新たな抽出領域は修正領域として抽出領域格納部33に格納される。ステップ45で抽出された修正領域について、可視化処理実行部31はステップ43の場合と同様に可視化処理を実行し、表示するステップ46。この後ステップ44~46を繰り返し、ステップ44にて関心領域を過不足なく抽出できたと判断された時点で、抽出処理の終了となるステップ47。三次元領域を複数の二次元断面像の集合として表示する場合、抽出領域格納部33に格納される三次元領域が複数あってかつ現在表示されている三次元領域とは別の三次元領域に表示を切り替える必要が生じることがある。たとえば図2の例だと、当初初期領域について可視化表示をしていたところ、修正領域に可視化処理の対象を切り替えて表示する場合などである。この場合、元のボリュームデータから作成される断面の角度とそのスライス間隔が同じであっても、表示を指示された三次元領域が、現在表示されている三次元領域とサイズが異なる場合、表示すべき二次元断面の数が異なることもある。そのような場合であっても、可視化処理実行部31は表示する二次元断面の数を自動的に変更し、表示すべき三次元領域全体をカバーするように二次元断面の集合を選択して表示する。これは先述の通り、可視化処理実行部31が可視化処理を実行するたびに、抽出領域を含むスライスを記憶しておき、記憶したスライスのデータを表示するようにしているからである。三次元領域を二次元断面の集合として表示する場合の詳細について、図3を用いて説明する.図3は、可視化処理実行部31で実行される三次元領域の可視化処理により表示装置12に示される画面表示の一例を示すイメージ図である。図3においてハッチングされた部分が抽出領域である。図3左は初期領域の表示を表す。初期領域は入力ボリュームデータから作成された複数スライスのうち、3つのスライスに跨っている。したがって3つのスライスのデータが表示される。ステップ44で行われる確認において、修正の必要あり,と判断された場合の修正領域1の表示を図3中に示す。ここでは修正領域1のスライス方向の大きさがが、初期領域よりも大きい場合を例示している。修正領域1は初期領域に対してスライス二枚分大きいサイズであったためすなわち修正領域1は初期領域よりサイズが大きく5つのスライスに跨っていたため、修正領域1全てが表示されるよう、自動的に表示スライス数が変更されている様子を示している。修正領域1にも修正が必要だと判断され、修正領域1を修正した次の修正領域2の表示を図3右に示す。修正領域2の大きさは修正領域1よりも小さく、修正領域2は4つのスライスに跨っているめ、表示スライス数を減じて表示が行われている。ここで,ステップ43および46における表示に際し,上記に述べたような二次元スライス表示に加えて、三次元形状を二次元画面上である程度把握することを目的とし、三次元可視化画像を同時に表示することもできる。三次元可視化画像を作成する方法としては例えば、ボクセル値から透明度を設定し、各視線上にあるボクセルに対して光の吸収と拡散を想定して光を加算して立体的に表示するVolumeRenderingVRや、各視線上にあるボクセルの、最大のボクセル値を投影するMaximumIntensityProjectionMIP等が挙げられる。この方法を用いれば、初期領域から修正されていく過程で、正確な形状と大まかな三次元形状を1つの画面でとらえることが可能である。次に、ステップ45における領域の修正処理の具体的な例について述べる。初期領域は入力ボリュームデータの画素やその部位において一般的で尤もらしいと判断される形状に抽出された画像であるが、入力ボリュームデータの画像特徴やシステムの利用場面によっては異なる抽出形状が望ましい場合もある。ここでは例の一つとして、均質なCT値を持つ臓器内に発生した腫瘍を抽出する場合について述べる。CT画像を用いた画像診断では、臓器の健全な実質部分と腫瘍に代表される病変部分を画像上で見分けるため、造影剤が利用されることも多い。造影剤とは、血管から血液内に注入する特定の薬剤で、血液に混入した造影剤が腫瘍に到達する時間と、腫瘍から臓器実質に造影剤が流れ出るまでの時間が異なるため、造影剤が注入されてから撮影するまでの時間に応じて、腫瘍と臓器実質の輝度が異なることが知られている。例えば臓器実質よりも輝度が低くなるような画像を対象とした腫瘍領域抽出処理としては、輝度が低い部分を腫瘍として抽出することが一般的であるが、腫瘍の性質によっては腫瘍内に更に低輝度な部分を生じることがある。原因としては腫瘍内の造影剤の広がりにむらがある場合や、腫瘍内に血栓ができている場合、腫瘍内に壊死領域がある場合等が考えられる。このような場合、輝度が低い部分を抽出することを基本アルゴリズムとする抽出処理を行うと、最も低輝度な領域のみが初期領域として抽出されることが多い。しかしこの場合、最も低輝度な領域のみを抽出するのではなく、画素の値としては最も低輝度な領域と高輝度な領域臓器実質との中間の値となる領域も腫瘍として抽出することが望ましい。本手法によれば、初期領域の確認時図2のステップ44に過不足ありと入力装置13から入力し、抽出領域格納部内で次に尤もらしい領域を領域修正実行部33が抽出して図2のステップ45、可視化処理実行部31が提示することで図2のステップ46、所望の抽出領域を得ることができる。この修正処理の一例を、図4を用いて説明する。図4は、自動抽出処理実行部34で抽出される初期領域が、領域修正処理実行部32で抽出アルゴリズムのパラメータを変更することで修正された例を示したイメージ図である。ここでは図4上の輝度を持つような画像を入力画像と仮定する。説明のため、領域1、2、3はそれぞれ領域内で単一の輝度を持つと仮定し、自動抽出アルゴリズムは閾値による二値化処理で、予め設定された輝度よりも小さい輝度を持つ領域を抽出領域とするというアルゴリズムであったとする。ここで、領域1は輝度1、領域2は輝度2、領域3は輝度3なる輝度を持つボクセルで構成され、輝度1は輝度2よりも大きく、輝度2は輝度3よりも大きいとする。予め設定された初期閾値が領域3よりも大きく、領域2よりも小さい場合、自動抽出処理の結果としては領域3のみ、即ち図4下左に示すように初期領域Aが抽出される。例えば領域2が臓器全体の領域、領域3が病変部分の領域であり、処理の目的としては病変部分の抽出であったとすれば、初期領域Aは目的の領域を過不足なく抽出できていることとなるので、処理はここで終了する。しかし、領域1が臓器の一部を写ったもので、領域2および領域3が病変部分であった場合、初期領域Aは病変部分の一部しか抽出できていないこととなるため、この場合は領域修正処理が必要となる。そこで、初期領域Aの次に抽出結果として尤もらしいものを抽出するため、抽出アルゴリズムのパラメータを変更する。即ち、ここでの抽出アルゴリズムは閾値による二値化処理であるので、初期閾値を輝度2よりも大きく、輝度1よりも小さい修正閾値Aへと変更することで、図4下右に示す修正領域Bを得ることができ、目的の領域を過不足なく抽出することが可能となる。このような輝度の変化は三段階以上にわたることもあるが、修正領域確認時に過不足ありと入力していくことで三段階以上の輝度の変化にも対応可能である。また、修正領域の作成および領域格納部への格納は、初期領域の表示前に全て行っていてもよいし、表示処理と並列して処理を行うことや、過不足ありと入力されてから処理を行うことも可能である。初期領域の次に尤もらしい領域の抽出とは例えば,初期領域の抽出時に利用された抽出アルゴリズムのパラメータを変更し,より外側に領域境界を設定するような方法などがあり,これは初期領域を元に修正領域を作成するため、複数の領域抽出処理を行うよりも処理は短時間で済む。また勿論これは抽出対象領域がその周辺領域よりも高輝度となる場合であっても適用可能である。修正閾値の設定方法としては、図5に示すようなパラメータ表を、領域修正処理実行部が予め持っていて、それを参照して閾値を設定する方法でもよいし、図6に示すパラメータ表の例のように、パラメータを式の形で持っておき、修正閾値とする方法もとれる。ここでf1、f2、f3、f’1、f’2、f’3は、予め設定する任意の関数であり、f1=f’1、f2=f’2、f3=f’3、でもよいし、個々に設定されていてもよい。それ以外にもパラメータそのものを規定するのではなく、アルゴリズムの実行結果を表として保持しておくこともできる。例えば、図7は修正処理項目表の例であるが、このように設定した修正項目表を参照し、各修正項目を全て満たす場合、もしくはいずれかを満たす場合、一定数以上満たす場合、等で収束条件を決めておく。次に,輝度に忠実な形での抽出領域よりも、それを覆う形での抽出領域が求められる場合について述べる。即ち,関心領域の範囲としては正しい抽出が行えているが、その利用場面によって、球に近い形や、より滑らかな抽出形状が望まれる場合などである。この場合は抽出アルゴリズムやそのパラメータを変更するのではなく,初期領域から大きく逸脱しない形で形状を変化させて表示させていく。つまり,領域輪郭に対して多角形近似や離散的正弦変換を利用した近似を行ったり,モルフォロジカルフィルタを利用して滑らかな輪郭にしたりする方法を用いる。1つの手法でパラメータを変更させる場合や,複数の手法を適用することで複数の輪郭が生成され,尤度は適用場面によって異なってくるとも考えられるが,一般的な場合の一つの方法としては,初期領域に最も近い形から球に代表される三次元幾何学図形に近い形までを、できるだけ段階的に変化させる順序を尤もらしさが低下する軸とすることが挙げられる。尤もらしさの低下を予めルールとして定義しておき、このルールに沿って抽出領域の形状が変化するよう領域修正を行う。ここでは、尤もらしさの低下は領域形状が初期領域から楕円球へと変わる順であると定義した場合の例を具体的に述べる。まず、自動抽出手法を用いて初期領域を抽出する。ここでの自動抽出手法は先述の領域拡張法やLevelSet等のアルゴリズムを用いてもよいし,もしくは画像の特徴によっては閾値を用いた二値化などの手法でもよい。こうして得られた初期領域に対する修正領域として、まずは領域内部の空白部分を塗りつぶすアルゴリズムを考える。塗りつぶしのためのアルゴリズムの例として、モルフォロジカルフィルタ等を用いる場合について図8を用いて説明する。図8は自動抽出処理実行部34で抽出される初期領域が領域修正処理実行部32で修正された領域の一例を修正の過程も含めて示したイメージ図である。まず、図8左に示す初期領域B内の空白部分を埋め尽くすまでDilation処理を行う。Dilatiob処理とは、領域輪郭を一画素分ずつ膨張させる処理である。この処理により初期領域B内部の空白部分が塗りつぶされた状態を修正途中領域として図8中に表した。この後に、Dilation処理を行った回数と同じ階数だけErosion処理を行ったものを修正領域イとして図8右に示す。Erosiob処理とはDilation処理を逆の操作を行う処理で、領域の輪郭を一画素分ずつ収縮させる。こうして、初期領域B内部の空白部分を塗りつぶした修正領域イが作成され、ここではこれが初期領域Bの次に尤もらしい領域として表示される。修正領域イよりも更に楕円球に近づける方法の一例を、図9を用いて説明する。図9は領域修正処理実行部32で修正された複数の修正領域の例を示すイメージ図である。図9左に示した修正領域イの輪郭の、凹となる部分を直線で置き換えたものが図9中の修正領域ロであり,初期領域の長径および短径を長軸および短軸とした楕円が図9右の修正領域ハである。こうして修正領域は、修正領域イ→修正領域ロ→修正領域ハの順に変化することとなる。なお、図4、図8および図9とそれに伴う説明は、抽出領域が二次元領域とした場合であるが、抽出領域が三次元領域であった場合でも同様のアルゴリズムによって修正領域を生成することが可能である。"