(57)【要約】【課題】過抽出を補正する際に，領域拡張（RG）のステップを戻す方法が適切であるかどうかが結果の表示でわかるようにする。【解決手段】Region Growingの抽出ステップで抽出領域130を色別表示し，隣接ボクセルとの差が大きい場合，そこに物質の境界がある可能性が高いと考え，色別表示の際の色の輝度をその部分については変更する。過抽出部分を手動修正する際には，画面内に色別表示の凡例をカラーバー131-133などで表示し，カラーバー131-133への操作を抽出領域130の修正として反映させる。【選択図】図１３
2010057671.tif 000002
【特許請求の範囲】
【請求項１】
記憶部に記憶された，複数の組織が存在するボリュームデータから特定の組織および部分のみを処理部で抽出する画像処理装置であって，
前記処理部は，自動抽出処理部と抽出結果表示部を有し，
前記自動抽出処理部は，
指定された点から特定の拡張条件に従って領域を拡張することで，目的の領域を抽出する領域自動抽出ステップと，
前記領域自動抽出ステップにおける拡張の過程を，拡張の各段階と前記各段階で拡張された拡張領域とを対応付けて前記記憶部に記録する拡張過程記録ステップと，
複数種類の色を持つ少なくとも１種類の色リストを用意し，前記領域自動抽出処理ステップにおける前記各段階に前記色リストの異なる色を割り当て，前記各段階と前記色リストの各色を対応付けた少なくとも1種類の色凡例を前記記憶部に記録する段階別色割当ステップと，
を実行し，
前記抽出結果表示部は，
前記拡張過程記録ステップで記録された全ての前記拡張領域を，抽出領域として表示する抽出領域表示ステップと，
前記拡張過程記録ステップと前記段階別色割当ステップで記録された情報から，前記抽出領域表示ステップにおける表示に，前記各段階での前記拡張領域に前記色リストの各色を反映させて表示する抽出過程表示ステップと，
を実行する，
ことを特徴とする画像処理装置。
【請求項２】
請求項１に記載の画像処理装置であって，
前記抽出結果表示部の前記抽出領域表示ステップは，前記抽出領域をそれ以外の領域とは区別して表示する，
ことを特徴とする画像処理装置。
【請求項３】
請求項１に記載の画像処理装置であって，
前記自動抽出処理部の前記段階別色割当ステップは，
前記各段階での拡張過程において，一つ前の段階の前記拡張領域に属する隣接ボクセルとの差を，前記色リストの輝度に反映させる段階別輝度割当ステップ，
を有することを特徴とする画像処理装置。
【請求項４】
請求項１に記載の画像処理装置であって，
前記抽出結果表示部の前記抽出過程表示ステップは更に，
前記段階別色割当ステップおよび前記段階別輝度割当ステップで記録された前記色凡例を，前記抽出領域の表示と同一画面内に色リストとして表示する色凡例表示ステップと，
ユーザまたはシステムからの指示によって前記色凡例の中から一つを選択する色凡例選択ステップ，
を有することを特徴とする画像処理装置。
【請求項５】
請求項１に記載の画像処理装置であって，
前記抽出結果表示部の前記抽出領域表示ステップは更に，
ユーザまたはシステムからの指示によって前記ボリュームデータの可視化表示方法を変更する表示法変更ステップ，
を有することを特徴とする画像処理装置。
【請求項６】
請求項１に記載の画像処理装置であって，
モニタを更に有し，前記モニタは，前記処理部の出力に基づき，前記抽出領域を，前記各段階での前記拡張領域に対応する前記色リストの各色を反映して表示する，
ことを特徴とする画像処理装置。
【請求項７】
請求項４に記載の画像処理装置であって，
モニタを更に有し，前記モニタは，表示画面中に，前記処理部の出力に基づき，前記抽出領域を，前記各段階での前記拡張領域に対応する前記色リストの各色を反映して表示し，且つ前記表示画面中の別の領域に，少なくとも１個の前記色凡例を表示する，
ことを特徴とする画像処理装置。
【請求項８】
記憶部に記憶された，複数の組織が存在するボリュームデータから特定の組織および部分のみを処理部で抽出する画像処理方法であって，
前記処理部は，自動抽出処理ステップと抽出結果表示ステップを実行し，
前記自動抽出処理ステップは，
指定された点から特定の拡張条件に従って自動的に領域を拡張することで，目的の領域を抽出する領域自動抽出ステップと，
前記領域自動抽出ステップにおける拡張の過程を，拡張の各段階と前記各段階で拡張された拡張領域とを対応付けて記憶領域に記録する拡張過程記録ステップと，
複数種類の色を持つ一種類または複数種類の色リストを用意し，前記領域自動抽出処理ステップにおける前記各段階に前記色リストの異なる色をそれぞれ割り当て，前記各段階と前記色リストの各色を対応付けた1種類または複数種類の色凡例を記憶領域に記録する段階別色割当ステップとを有し，
前記抽出結果表示ステップは，
前記拡張過程記録ステップで記録された全ての前記拡張領域を，抽出領域としてそれ以外の領域とは区別して表示する抽出領域表示ステップと，
前記拡張過程記録ステップと前記段階別色割当ステップで記録された情報から，前記抽出領域表示ステップにおける表示に，前記各段階での前記拡張領域に前記色リストの各色を反映させて表示する抽出過程表示ステップとを有する，
ことを特徴とする画像処理方法。
【請求項９】
請求項８に記載の画像処理方法であって，
前記段階別色割当ステップは，前記各段階での拡張過程において，一つ前の段階の前記拡張領域に属する隣接ボクセルとの差を前記色リストの輝度に反映させる段階別輝度割当ステップを有する，
ことを特徴とする画像処理方法
【請求項１０】
請求項８に記載の画像処理方法であって，
前記抽出過程表示ステップは更に，
前記段階別色割当ステップおよび前記段階別輝度割当ステップで記録された前記色凡例を前記抽出領域の表示と同一画面内に色リストとして表示する色凡例表示ステップと，
ユーザまたはシステムからの指示によって複数種類の前記色凡例の中から一意または任意に一つを選択する色凡例選択ステップを有する，
ことを特徴とする画像処理方法。
【請求項１１】
請求項８に記載の画像処理方法であって，
前記抽出領域表示ステップは更に，
ユーザまたはシステムからの指示によって前記ボリュームデータの可視化表示方法を変更する表示法変更ステップを有する，
ことを特徴とする画像処理方法。
【請求項１２】
複数の組織が存在するボリュームデータから特定の組織および部分のみを抽出するためのコンピュータで実行されるプログラムであって，
自動抽出処理ステップと抽出結果表示ステップとをコンピュータに実行させ，
前記自動抽出処理ステップは，
指定された点から特定の拡張条件に従って自動的に領域を拡張することで，目的の領域を抽出する領域自動抽出ステップと，
前記領域自動抽出ステップにおける拡張の過程を，拡張の各段階と前記各段階で拡張された拡張領域とを対応付ける形で記憶領域に記録する拡張過程記録ステップと，
複数種類の色を持つ１種類または複数種類の色リストを用意し，前記領域自動抽出処理ステップにおける前記各段階に色リストの異なる色をそれぞれ割り当て，前記各段階と前記色リストの各色を対応付けた1種類または複数種類の色凡例を記憶領域に記録する段階別色割当ステップとを有し，
前記抽出結果表示ステップは，
前記拡張過程記録ステップで記録された全ての前記拡張領域を，抽出領域としてそれ以外の領域とは区別して表示する抽出領域表示ステップと，
前記拡張過程記録ステップと前記段階別色割当ステップで記録された情報から，前記抽出領域表示ステップにおける表示に，前記各段階での前記拡張領域に前記色リストの各色を反映させて表示する抽出過程表示ステップとを有する，
ことを特徴とする画像処理プログラム。
【請求項１３】
請求項１２に記載の画像処理プログラムであって，
前記段階別色割当ステップは，前記各段階での拡張過程において，一つ前の段階の前記拡張領域に属する隣接ボクセルとの差を前記色リストの輝度に反映させる段階別輝度割当ステップを有する，
ことを特徴とする画像処理プログラム。
【請求項１４】
請求項１２に記載の画像処理プログラムであって，
前記抽出過程表示ステップは更に，
前記段階別色割当ステップおよび前記段階別輝度割当ステップで記録された前記色凡例を前記抽出領域の表示と同一画面内に色リストとして表示する色凡例表示ステップと，
ユーザまたはシステムからの指示によって複数種類の前記色凡例の中から一意または任意に一つを選択する色凡例選択ステップとを有する，
ことを特徴とする画像処理プログラム。
【請求項１５】
請求項１２に記載の画像処理プログラムであって，
前記抽出領域表示ステップは更に，
ユーザまたはシステムからの指示によってボリュームデータの可視化表示方法を変更する表示法変更ステップを有する，
ことを特徴とする画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は，3次元情報を有するボリュームデータから必要な情報を抽出または不必要な情報を除去する画像処理技術に関するものである。
【背景技術】
【０００２】
X線CT(X-ray Computed Tomography)やMRI(Magnetic Resonance Imaging)等で撮像されたボリュームデータを用いた診断においては，3次元情報を有するボリュームデータを可視化処理することで2次元画像を生成し，その2次元画像を観察しながら診断を行う場合がある。医用画像においては一般的に，画像内にピクセル値として記録される値は，臓器や組織によってその範囲が異なることが一般的であるので，フィルタリング処理を行うことで，注目したい臓器や組織のみを表示することが可能となる。
【０００３】
しかし，例えば血管に造影剤を注入して撮像するとその造影血管領域の値の範囲は骨領域の値の範囲と重なる場合がある。このように値の範囲が重なる場合には，単純なフィルタリング処理では注目したい臓器や組織のみを表示することは難しい。
【０００４】
このような場合には従来，CT画像一枚一枚から，医師や技師が手作業で目的とする領域以外を除去することで，目的とする領域のみを表示した画像を作成していた。他方、画像一枚だけでなくボリュームデータとしても連続的なデータを作成するためには，操作が煩雑になりコストも高く，また高精度で効率的な抽出のためには高度の習熟が必要となる。これらのことから，複数枚のCT画像データを1つのボリュームデータとし，自動的に目的の領域以外を除去する技術の開発が行われてきた。
【０００５】
しかしある程度高精度な自動抽出処理であったとしても，最終的に医師や技師による確認や手作業による修正は必要不可欠なものであり，この際の医師や技師にかかる負担やコストは，依然として残っている。
【０００６】
自動抽出処理アルゴリズムの代表的なものとして，領域拡張法（Region Growing）があるが，このアルゴリズムによる処理結果には，過抽出と呼ばれる拡張の際の漏れが問題として起こりがちである。これを修正する手法として，特許文献１では逆Region Growing法を開示している。
【０００７】
その他のアルゴリズムの例として，非特許文献1で用いられている，Watershedと呼ばれる領域分割手法がある。この手法による領域分割の過程には，細分化されてしまった小領域を，パラメータを段階的に変動させることで統合していく処理がある。
【０００８】
【特許文献１】特許2874445公報
【非特許文献１】Alan P. Mangan, Partitioning 3D Surface Meshes Using Watershed Segmentation, IEEE Transaction on Visualization and Computer Graphics, Vol.5, No.4, 1999
【発明の開示】
【発明が解決しようとする課題】
【０００９】
例えばRegion Growingの場合，上述した逆Region Growing法により修正できる部分は，拡張過程の新しい段階の部分のみであり，Region Growingの結果のみからでは，拡張がどのように進んできたかは判断できない。逆Region Growingの他には，拡張しすぎた部分に対して拡張段階を数段階遡って修正する方法もとることができるが，この場合にも，結果の領域だけ見てもどのような過程で成長したのか判別できない。以前は計算機の処理が遅く，計算時間が長くかかる場合には履歴を表示しながら拡張することが可能だったが，現在は計算機の処理速度も上がり，拡張し終わった結果しか見られない場合も多くある。
【００１０】
また，Watershedの統合処理の場合，本来は別の領域として分割されるべき領域が，パラメータによっては統合されてしまう。それを回避するため，フィルタリング処理等を用いて予め画像を平滑化してから処理を行う方法なども提案されているが，統合された後の結果から修正しようとする場合，どのような過程で統合処理が進んだのかはわからない。
【００１１】
本発明の目的は，上記の点に鑑み，自動抽出処理の結果を手動で修正する際の，煩雑さや困難さを軽減することが可能な画像処理装置，処理方法，及び処理プログラムを提供することにある。
【課題を解決するための手段】
【００１２】
上記目的を達成するために，本発明は，抽出結果を表示する際に抽出ステップを色別表示し，各抽出ステップで変わっていく変化量を，色別表示の際の色の輝度，に反映させる。
【００１３】
即ち，本発明では上記課題を解決するため，記憶する複数の組織が存在するボリュームデータから，処理部で特定の組織および部分のみを抽出する画像処理装置であって，処理部は自動抽出処理ステップと抽出結果表示ステップとを実行し，自動抽出処理ステップは，指定された点から特定の拡張条件に従って領域を拡張することで，目的の領域を抽出する領域自動抽出ステップと，領域自動抽出ステップにおける拡張の過程を，拡張の各段階とその段階で拡張された領域とを対応付ける形で逐次的に記憶領域に記録する拡張過程記録ステップと複数種類の色を持つ一種類または複数種類の色リストを用意し，領域自動抽出処理ステップにおける各段階に色リストの異なる色をそれぞれ割り当て，各段階と色リストの各色を対応付けた1種類または複数種類の凡例を記憶領域に記録する段階別色割当ステップとを有し，抽出結果表示ステップは，拡張過程記録ステップで記録された全拡張領域を，抽出領域として表示する抽出領域表示ステップと拡張過程記録ステップと段階別色割当ステップで記録された情報から，抽出領域表示ステップにおける表示に，各段階での拡張領域に色リストの各色を反映させ表示する抽出過程表示ステップとを有することを特徴とする画像処理装置を提供する。
【００１４】
また，本発明は上記課題を解決するために，記憶部に記憶された，複数の組織が存在するボリュームデータから特定の組織および部分のみを処理部で抽出する画像処理方法であって，処理部は，自動抽出処理ステップと抽出結果表示ステップを実行し，自動抽出処理ステップは，指定された点から特定の拡張条件に従って自動的に領域を拡張することで，目的の領域を抽出する領域自動抽出ステップと，領域自動抽出ステップにおける拡張の過程を，拡張の各段階と各段階で拡張された拡張領域とを対応付けて記憶領域に記録する拡張過程記録ステップと，複数種類の色を持つ一種類または複数種類の色リストを用意し，領域自動抽出処理ステップにおける各段階に色リストの異なる色をそれぞれ割り当て，各段階と色リストの各色を対応付けた1種類または複数種類の色凡例を記憶領域に記録する段階別色割当ステップとを有し，抽出結果表示ステップは，拡張過程記録ステップで記録された全ての拡張領域を，抽出領域としてそれ以外の領域とは区別して表示する抽出領域表示ステップと，拡張過程記録ステップと段階別色割当ステップで記録された情報から，抽出領域表示ステップにおける表示に，各段階での拡張領域に色リストの各色を反映させ表示する抽出過程表示ステップとを有する画像処理方法を提供する。
【００１５】
更に，本発明は上記課題を解決するために，複数の組織が存在するボリュームデータから，コンピュータが特定の組織および部分のみを抽出することを実行させることを目的とするプログラムであって，自動抽出処理ステップと抽出結果表示ステップとを備え，自動抽出処理ステップは，コンピュータが，指定された点から特定の拡張条件に従って自動的に領域を拡張することで，目的の領域を抽出する領域自動抽出ステップと，領域自動抽出ステップにおける拡張の過程を，拡張の各段階とその段階で拡張された領域とを対応付ける形で逐次的に記憶領域に記録する拡張過程記録ステップと複数種類の色を持つ一種類または複数種類の色リストを用意し，領域自動抽出処理ステップにおける各段階に色リストの異なる色をそれぞれ割り当て，各段階と色リストの各色を対応付けた1種類または複数種類の凡例を記憶領域に記録する段階別色割当ステップとを有し，抽出結果表示ステップは，拡張過程記録ステップで記録された全拡張領域を，抽出領域としてそれ以外の領域とは区別して表示する抽出領域表示ステップと拡張過程記録ステップと段階別色割当ステップで記録された情報から，抽出領域表示ステップにおける表示に，各段階での拡張領域に色リストの各色を反映させ表示する抽出過程表示ステップとを有することを特徴とする画像処理プログラムを提供する。
【発明の効果】
【００１６】
本発明によれば，過抽出を手動で補正する際に拡張の段階を戻す方法が適切であるかどうかが結果の表示をみることでわかり，その結果に基づいた手動補正の支援も可能となる。このことにより，ユーザの手間や負担を軽減することができる。
【発明を実施するための最良の形態】
【００１７】
以下に，図を用いて実施の形態について説明する。
【００１８】
図1は，以下の説明においても用いる，ユーザが注目したい領域を抽出しようとするボリュームデータの概念図である。ハッチングの種類の違いは，ボクセル値の違いを表す。
【実施例】
【００１９】
図２は，画像処理装置の一実施例のシステム構成図である。画像処理装置は中央処理部（Central Processing Unit：ＣＰＵ）と記憶部（メモリ）と入出力部とからなるコンピュータ60と，入力装置11，モニタ12，外部記憶装置50とからなる。コンピュータ60内には，内部メモリ10とＣＰＵからなる画像処理部20を図示し，画像処理部20は，更に機能ブロックとして，自動抽出処理部30と抽出結果表示部40を図示した。これらの機能ブロックは，後で説明するように，ＣＰＵで実行されるプログラムで構成される。
【００２０】
内部メモリ10は，外部記憶装置50から図示されない入出力部を介して，ボリュームデータを入力として受信する。画像処理部20の自動抽出処理部30は，入力装置11を介して抽出処理の起動指示を受信した場合に，内部メモリ10の入力ボリュームデータに対して自動抽出処理を開始する。抽出結果表示部40は自動抽出処理が終了したという情報を自動抽出処理部30から受信した場合に，抽出結果表示処理を開始し，入力装置11からの指示を受け，内部メモリ10の入力ボリュームデータを用いてモニタ12への表示を行う。  
図３に本実施例の自動抽出処理部30の機能構成の一例を示すシステム構成図を，図４に自動抽出処理部30で実行される処理の流れを示す。
【００２１】
領域自動抽出部31は，内部メモリ10から受け取った例えば図１に示すようなボリュームデータを入力データとし(S100)，それに対して自動で領域拡張を行う(S101)。ここでは領域拡張手法の例として，Region Growing法を用いる場合について述べる。Region Growing手法ではまず，自動または手動で特定の領域を指定する(t=0)。ここでtは拡張段階を表す。この領域は抽出したい領域内に存在する部分であればどこでもよく，領域の大きさは1ボクセルから，抽出したい領域の大きさ以下となる。また，ここで領域として指定される各ボクセルをシードと呼ぶ。
【００２２】
特定の条件に従って拡張を行うが，拡張に使用する特定の条件とは例えば，予め手動または自動で設定されたLocal ThresholdとGlobal Threholdと呼ばれる二種類の閾値と，シードとして指定されたボクセルのボクセル値の最小値，最大値，平均値等の値をシード値として用いる。当該ボクセルとそれに隣接する六つのボクセルとのボクセル値の差を算出し，その差がLocal Threshold以下で，かつシード値と当該ボクセルの隣接ボクセルのボクセル値の差を算出し，その差がGlobal Threshold以下である場合には，隣接ボクセルへと領域を拡張する，などの条件を設定する。このとき，その段階までに一度拡張された領域については参照しない。このようにして拡張していき，特定の条件に適うボクセルが隣接ボクセル中にない場合，または別に設定された終了条件に適合した場合，拡張が終了する。別に設定された終了条件とは例えば，Limitという特定の最大拡張ボクセル数を用いるもので，拡張領域として拡張されたボクセルがLimitボクセル目となった時点で拡張が終了する，などの条件を設定できる。
【００２３】
図5のt=0〜t=Tのハッチングが施された部分は，その段階までに拡張された領域全てを表し，ここでの拡張の段階はt=0からt=Tまでの(T+1)段階となる。図5で示すように拡張する場合，各段階での拡張領域は図6のようになる。図６にはt=0〜t=Tにおける拡張領域R0〜RTを示す。
【００２４】
拡張過程記録部32は，例えば図7に示す拡張過程記録表320に，各拡張段階と，各拡張段階における拡張領域を対応付け，拡張が一段階進むごとにそれらを逐次的に記録する(S102)。拡張領域R0〜RTは例えば式(1)のように，その領域に属するボクセルの3次元座標値のリストとして表される。
【００２５】
【数１】
2010057671.tif 000003
段階別色割当部33が記録する色凡例は，例えば図8に示す色対応表330のように，各拡張段階と複数種類の色リストを記録する(S103)。色リストの各色は例えばR(赤),G(緑),B(青)の三色組み合わせとして記録され，実際に色を表示に使う際には，この三色の加法混色を用いる。まずCt は，拡張が終了するまでの全段階数であるT種類の色を持つ色リストであり，式(2)で表すことができる。
【００２６】
【数２】
2010057671.tif 000004
C'tは式(3)で表されるように，CT0〜CT4の5種類の色を利用し，全拡張段階が5段階以上になる場合には6段階目にCT0を用いるというように，循環的に限られた数の色を利用する色リストの例である。尚，式(3)においてt%5は，t/5の剰余とする。
【００２７】
【数３】
2010057671.tif 000005
またC''tは式(4)に表されるように，段階t0以前は全て同じ色を設定し，t0よりもあとの段階では，(T- t0)種類の色を利用する色リストの例である。
【００２８】
【数４】
2010057671.tif 000006
これらの色凡例は，拡張過程記録部32と同様，拡張が一段階進むごとに逐次的に色を割り当てていってもよいし，全拡張過程が終了した時点で，全段階に各々色を割り当ててもよい。
【００２９】
段階別輝度割当部34で色対応表に記録する色凡例としては，例えば式(5)に表されるC'''tのような色リストを考えることができる。
【００３０】
【数５】
2010057671.tif 000007
式(４)において，LはLocal Threshold，Dは隣接ボクセルとのボクセル値の差の平均を表す。
【００３１】
C'''tは領域自動抽出部31で用いられる特定の条件の中に，隣接ボクセルのボクセル値との差がLocal Threshold以下であれば拡張させる，という条件がある場合の例となるが， C'''tはその段階での拡張領域に属する全ボクセルの各々の一つ前の拡張段階のボクセルとのボクセル値の差の平均等の値を輝度に反映させた色リストとなっている。C'''tでは，R:G:Bの割合は拡張段階によって決定されるが，R+G+Bの合計値は，隣接するボクセルと当該ボクセルのボクセル値との差Dが大きければ小さく，Dが小さければ大きくなる。即ち，Dが大きければ大きいほどその領域の色は黒に近くなり，Dが小さければ小さいほど，その領域の色は白に近くなる。尚，拡張過程記録表320に登録されない情報を色リストに反映させる場合で，更に拡張が終了してから色を割り当てる場合には，拡張過程記録表320に記載項目を追加する必要がある。拡張過程記録表320に記載項目として、拡張時の隣接ボクセルとの差の平均を追加して例を図９にB321として示す。
【００３２】
このようにして，自動抽出処理部30で記録された拡張過程記録表320，段階別色割当部33で記録された色対応表330の情報を元に，画像処理部20の抽出結果表示部40が起動される(S105)。
【００３３】
図１０に，本実施例における抽出結果表示部40の構成の一例を示すシステム構成図，図１１に抽出結果表示部40にて実行される処理の流れを示す。
【００３４】
まず抽出領域表示部41は，入力ボリュームデータに対し，例えば抽出されなかった領域は(R,G,B)=(0,0,0)即ち黒とし，抽出された領域はそれ以外の色とするなど，抽出された領域と抽出されなかった領域とを明示的に区別できるような色を割り当てて表示を行う(S200)。
【００３５】
色凡例選択部44はユーザからの指示もしくは予め定められた初期値を用いて，色凡例を選択する(S201)。ここでの指示とは，ユーザが色リスト自体へ直接指定を行うことや，ユーザが選択した作業に対し，予め内部メモリに保存されていた対応表を参照し，それに適した色リストをシステムが自動的に選択することである。そして抽出過程表示部42は，拡張過程記録部32で記録された拡張過程記録表320やB321，および段階別色割当部33と段階別輝度割当部34で記録された色対応表330を参照し，各段階で抽出された領域に，各々異なる色を割り当て，例えば図１２のように抽出領域の表示に抽出過程の情報を付与した表示を行う(S202)。 
また，色凡例表示部43は，図１３に示すように，抽出領域の表示と同時に抽出領域130を表示する画面の同一画面内に色凡例を表示を行う(S203)。なお,
図13,14の矢印は拡張段階の進む方向を示すものである。
【００３６】
ここで，色対応表330に複数の色リストが存在する場合，色凡例選択部44は、ユーザからの指示による色凡例131-133の選択に従って，抽出領域の色表示をそれに従って変更することもできる。ここでの指示とは，S201と同様に，ユーザが色リスト自体へ直接指定を行うことや，ユーザが選択した作業に対し，予め内部メモリに保存されていた対応表を参照し，それに適した色リストをシステムが自動的に選択することである。図１３では，現在表示に使用している色リスト131を，縁を太線にすることで示している。別の色リスト132，133が選択されると，選択された色リストの縁が太線になる。
【００３７】
また，段階別輝度割当部34において，前記の例では輝度の割当はその段階での拡張領域全体に設定したが，隣接ボクセルのボクセル値との差から，ボクセル単位で輝度を設定してもよい。ボクセル単位に設定する場合は，拡張領域への色の割当は逐次的に行われ，色リストへは輝度の変化は反映されないか，またはその段階における拡張領域に属するボクセルの輝度の平均などを色リストの輝度とする。
【００３８】
次に，自動抽出処理部30にてRegion Growingにより自動抽出を行った結果を上記のように表示した結果過抽出があり(S204)，それを手動修正する場合について述べる。 
まず，色凡例が修正に適した状態になっていない場合(S205)，色凡例選択部44は色凡例を変更する(S206)。ここでの指示とは，S201の場合と同じように，ユーザが色リストを直接選択することや，またはユーザが修正作業を開始させる指示したのちに，システムが予め内部メモリに保存されている対応表を参照して，その修正に適した色リストを選択することなどで実行される。変更された色凡例に従って，抽出過程表示部42は抽出領域の表示色を変更し(S207)，また色凡例表示部43も，それに伴って色凡例の表示を変更する(S208)。
【００３９】
例として，色リストに使用できる色の種類が限られている場合を図１４を用いて説明する。例えば色リストに使用できる色が三種類である場合，拡張終了段階ではBのように全体を三色に分けた色リストを用い，ユーザが，拡張領域の表示部分をクリックするなどの操作により，手動修正の開始を指示した場合に色リストDのように新しい段階が細かく観察できるような色リストへと変更する。
【００４０】
次にユーザは各拡張段階で拡張された領域に割り当てられ表示された彩度や輝度の情報を元に，図１５左側に示すように遡りたい領域と拡張したい領域の境界付近周辺をドラッグなどの操作で指定する(S209)。表示法変更部45は，抽出領域の表示に並列して，または重ねて，例えば指定された領域の断面表示などを行う(S210)。表示法の選択は，S201やS202と同様に，ユーザが直接変更する表示法を指示したり，もしくはユーザが修正の指示を出し，それを受け，システムがその場合に適した表示法を内部メモリに予め保存された対応表を参照することでシステムが選択したりすることで実行される。この断面表示では，ボリュームデータ内部の拡張過程を示すこともできるため，過抽出部分の拡張過程に関して三次元的な付加情報を提示できる。
【００４１】
ユーザは拡張領域の表示や表示法変更部45による断面表示を参照しながら，その場合に適した修正方法を検討し，拡張段階を戻すことで過抽出が修正できると判断された場合には(S211)，図１６に示すように，切断したい曲面の一点をクリックなどの操作で指定する(S212)。
【００４２】
抽出領域表示部41および抽出過程表示部42は，指定された一ボクセルが属する段階を判断したのち，拡張領域とその拡張段階を対応付けた表を参照し，指定された段階以降の段階に対応する拡張領域を，抽出領域から除外して抽出領域を表示する(S213)。
【００４３】
このような手動修正を繰り返し行う場合，拡張過程記録表320に記載されている情報に加え，修正した結果現時点でどこまでを抽出領域として表示しているかという情報を保存する必要がある。その場合の修正後拡張過程記録表322を，図１７に示す。修正後拡張過程記録表322において，修正後表示フラグが1となっている拡張段階での拡張領域は抽出領域として表示し，0となっている場合には表示しない。図１７に示した場合であれば，拡張段階2までの拡張領域を，抽出領域として表示していることになる。前記のような手動修正を行うたびに，この修正後表示フラグを変更することで，モニタに表示されている抽出領域と，修正後拡張過程記録表322を対応付けることが可能となる。
【００４４】
また，拡張段階を戻す場合には，色凡例に対して操作を行うこともできる。切断したい断面に付された色と，同じ色の部分を色凡例上で指定し，ドラッグなどの操作を行うことで，拡張段階を戻したり，進めたりという操作を繰り返しながら過抽出を修正することもできる。
【図面の簡単な説明】
【００４５】
【図１】ボリュームデータを説明するための図である。
【図２】は実施例の画像処理装置の全体のシステム構成図である。
【図３】実施例の自動抽出処理部のシステム構成図である。
【図４】実施例の自動抽出処理部で行われる処理の流れ図である。
【図５】実施例のボリュームデータに対する拡張の過程の説明図である。
【図６】実施例の各段階での拡張領域を示す図である。
【図７】実施例の拡張段階とその段階での拡張領域の対応を表す拡張過程記録表を示す図である。
【図８】実施例の拡張段階と色リストの色対応表を示す図である。
【図９】実施例の拡張段階とその段階での拡張領域の対応を表す拡張過程記録表を示す図である。
【図１０】実施例の抽出結果表示部のシステム構成図である。
【図１１】第１の実施例の抽出結果表示部の処理の流れを示す図である。
【図１２】実施例の拡張過程と色対応表を抽出領域表示に反映させた例を示す図である。
【図１３】実施例の抽出領域表示と色リストを同時に表示する例を示す図である。
【図１４】実施例の色リストの例を示す図である。
【図１５】実施例の表示法変更の例を示す図である。
【図１６】実施例の手動修正の例を示す図である。
【図１７】実施例の正後拡張過程記録表を示す図である。
【符号の説明】
【００４６】
10…内部メモリ
11…入力装置
12…モニタ
20…画像処理部
30…自動抽出処理部
31…領域自動抽出部
32…拡張過程記録部
320…拡張過程記録表
321…拡張過程記録表B
322…修正後拡張過程記録表
33…段階別色割当部
330…色対応表
34…段階別輝度割当部
40…抽出結果表示部
41…抽出領域表示部
42…抽出過程表示部
43…色凡例表示部
44…色凡例選択部
45…表示法変更部
50…外部記憶装置
60…コンピュータ。
【図１】
2010057671.tif 000008
【図２】
2010057671.tif 000009
【図３】
2010057671.tif 000010
【図４】
2010057671.tif 000011
【図５】
2010057671.tif 000012
【図６】
2010057671.tif 000013
【図７】
2010057671.tif 000014
【図８】
2010057671.tif 000015
【図９】
2010057671.tif 000016
【図１０】
2010057671.tif 000017
【図１１】
2010057671.tif 000018
【図１２】
2010057671.tif 000019
【図１３】
2010057671.tif 000020
【図１４】
2010057671.tif 000021
【図１５】
2010057671.tif 000022
【図１６】
2010057671.tif 000023
【図１７】
2010057671.tif 000024