(57)【要約】【課題】本発明の目的は、従来法のリージョングローイングで問題とされた、恣意的なパラメータや拡張の停止条件をなくし、自動的なセグメンテーションを実現できるようにすることを課題とする。【解決手段】複数フレームよりなる生体フルカラー画像（断面画像）について、関心組織領域の抽出対象フレームに対して１つ前の連続するフレームで抽出された関心組織領域と同じ位置に仮領域を与え、拡張リージョングローイング法により、その仮領域の境界近傍画素位置を中心として抽出対象フレームよりも前に抽出した複数フレームから注目画素を中心とした３次元の局所領域を与え、その局所領域について関心組織領域内外の濃度中央値を求めてそれぞれ注目画素の濃度値との関係に基づいて抽出対象フレームの関心組織領域を抽出する。【選択図】      図５
2004234579.tif 000002
【特許請求の範囲】
【請求項１】
生体断面画像の連続フレームから関心組織領域を抽出する方法であって、
前記連続フレームの一部に対して前記関心組織領域における初期の判定基準を与え、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記判定基準とに基づいて前記関心組織領域に帰属させる領域と帰属させない領域とを判断し、前記フレームに続くフレーム以降の判断では、前記関心組織領域に前記帰属させる領域を含めた領域の画素値に基づいて前記判定基準を設定することを特徴とする関心組織領域抽出方法。
【請求項２】
生体断面画像の連続フレームから関心組織領域を抽出する方法であって、
前記連続フレームの一部に対して前記関心組織領域内における初期の判定基準と前記関心組織領域外における初期の判定基準の２種類を与え、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記２種類の判定基準とに基づいて前記関心組織領域内に帰属させる領域と前記関心組織領域外に帰属させる領域とを判断し、前記フレームに続くフレーム以降の判断では、前記関心組織領域内に前記帰属させる領域を含めた領域の画素値と前記関心組織領域外に前記帰属させる領域を含めた領域の画素値に基づいて前記２種類の判定基準を設定することを特徴とする関心組織領域抽出方法。
【請求項３】
前記連続フレームの一部は１枚のフレーム又は複数枚のフレームであり、前記フレームに対して手操作により設定された前記関心組織領域から前記判定基準を求めることを特徴とする請求項１又は２に記載の関心組織領域抽出方法。
【請求項４】
前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレームの１つ前のフレームで判断された関心組織領域と同じ位置に仮領域を設け、当該仮領域の境界に沿って画素単位に帰属の判断を行うことを特徴とする請求項１，２又は３に記載の関心組織領域抽出方法。
【請求項５】
前記仮領域において、前記関心組織領域に帰属させる画素が前記仮領域外にある場合には、当該画素を前記判断対象のフレームの関心組織領域の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して拡張し、前記関心組織領域に帰属させない画素が前記仮領域内にある場合には、当該画素を前記判断対象のフレームの関心組織領域外の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して縮小することを特徴とする請求項４に記載の関心組織領域抽出方法。
【請求項６】
前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレーム以前の連続する複数枚のフレームから画素単位に対応させて前記判断対象となるフレームの各画素を含む所定サイズの局所領域を取り出しながら、前記各局所領域に含まれる前記関心組織領域の一部に対応する前記設定された判定基準を使用することを特徴とする請求項１〜５のいずれか１つに記載の関心組織領域抽出方法。
【請求項７】
前記関心組織領域の画素のみ表示対象として記憶することを特徴とする請求項１〜６のいずれか１つに記載の関心組織領域抽出方法。
【請求項８】
生体断面画像の連続フレームから関心組織領域を抽出するための処理をコンピュータに実行させるプログラムであって、
前記コンピュータに、
前記連続フレームの一部に対応する前記関心組織領域における初期の判定基準を算出する第１ステップと、
前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記算出された初期の判定基準とに基づいて前記関心組織領域に帰属させる領域と帰属させない領域とを判断する第２ステップと、
前記フレームに続くフレーム以降の判断では、前記関心組織領域に前記帰属させる領域を含めた領域の画素値に基づいて前記判定基準を設定する第３ステップと、
を実行させることを特徴とする関心組織領域抽出プログラム。
【請求項９】
生体断面画像の連続フレームから関心組織領域を抽出するための処理をコンピュータに実行させるプログラムであって、
前記連続フレームの一部に対して前記関心組織領域内における初期の判定基準と前記関心組織領域外における初期の判定基準の２種類を算出する第１ステップと、
前記連続フレームの一部に続くフレームについて当該フレームの各画素の濃度値と前記算出された２種類の判定基準とに基づいて前記関心組織領域内に帰属させる領域と前記関心組織領域外に帰属させる領域とを判断する第２ステップと、
前記フレームに続くフレーム以降の判断では、前記関心組織領域内に前記帰属させる領域を含めた領域の画素値と前記関心組織領域外に前記帰属させる領域を含めた領域の画素値に基づいて前記２種類の判定基準を設定する第３ステップと、
を実行させることを特徴とする関心組織領域抽出プログラム。
【請求項１０】
前記連続フレームの一部は１枚のフレーム又は複数枚のフレームであり、前記フレームに対して手操作により設定された前記関心組織領域から前記初期の判定基準を算出させる前記第１ステップを有したことを特徴とする請求項８又は９に記載の関心組織領域抽出プログラム。
【請求項１１】
前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレームの１つ前のフレームで判断された関心組織領域と同じ位置に仮領域を設け、当該仮領域の境界に沿って画素単位に帰属の判断を行わせる前記第２ステップを有したことを特徴とする請求項８，９又は１０に記載の関心組織領域抽出プログラム。
【請求項１２】
前記仮領域において、前記関心組織領域に帰属させる画素が前記仮領域外にある場合には、当該画素を前記判断対象のフレームの関心組織領域の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して拡張させ、前記関心組織領域に帰属させない画素が前記仮領域内にある場合には、当該画素を前記判断対象のフレームの関心組織領域外の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して縮小させる前記第２ステップを有したことを特徴とする請求項１１に記載の関心組織領域抽出プログラム。
【請求項１３】
前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレーム以前の連続する複数枚のフレームから画素単位に対応させて前記判断対象となるフレームの各画素を含む所定サイズの局所領域を取り出させながら、前記各局所領域に含まれる前記関心組織領域の一部に対応する前記設定された判定基準を使用させる前記第２ステップを有したことを特徴とする請求項８〜１２のいずれか１つに記載の関心組織領域抽出プログラム。
【請求項１４】
前記コンピュータに、さらに、前記関心組織領域の画素のみ表示対象として記憶する第４ステップを実行させることを特徴とする請求項８〜１３のいずれか１つに記載の関心組織領域抽出プログラム。
【請求項１５】
生体断面画像の連続フレームから関心組織領域を抽出する画像処理装置であって、
前記連続フレームの画素値を記憶する画素値記憶手段と、
前記画素値記憶手段を参照して前記連続フレームの一部に対応する前記関心組織領域における初期の判定基準を算出する算出手段と、
前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記算出された初期の判定基準とに基づいて前記関心組織領域に帰属させる領域と帰属させない領域とを判断する判断手段と、
前記フレームに続くフレーム以降の判断では、前記関心組織領域に前記帰属させる領域を含めた領域の画素値に基づいて前記判定基準を設定する設定手段と、
を備えることを特徴とする画像処理装置。
【請求項１６】
生体断面画像の連続フレームから関心組織領域を抽出する画像処理装置であって、
前記連続フレームの画素値を記憶する画素値記憶手段と、
前記画素値記憶手段を参照して前記連続フレームの一部に対して前記関心組織領域内における初期の判定基準と前記関心組織領域外における初期の判定基準の２種類を算出する算出手段と、
前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記算出された２種類の判定基準とに基づいて前記関心組織領域内に帰属させる領域と前記関心組織領域外に帰属させる領域とを判断する判断手段と、
前記フレームに続くフレーム以降の判断では、前記関心組織領域内に前記帰属させる領域を含めた領域の画素値と前記関心組織領域外に前記帰属させる領域を含めた領域の画素値に基づいて前記２種類の判定基準を設定する設定手段と、
を備えることを特徴とする画像処理装置。
【請求項１７】
さらに、前記関心組織領域の各画素に対応させて表示対象を示すフラグ情報を記憶する表示対象記憶手段を有することを特徴とする請求項１５又は１６に記載の画像処理装置。
【請求項１８】
さらに、前記画素値記憶手段を参照し、前記表示対象記憶手段に記憶されたフラグ情報に基づいて前記関心組織領域を表示する表示手段を有することを特徴とする請求項１７に記載の画像処理装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、例えば、関心組織領域抽出方法、関心組織領域抽出処理プログラム及び画像処理装置に関する。
【０００２】
【従来の技術】
今日、期待される生体内組織のディジタルモデルとは、生体内の各組織の構造、及び力学的特性の詳細な情報を保持する３次元モデルを意味する。もし医療に応用できれば、診断や治療、手術の事前シミュレーションが可能となり、その効果は計り知れない。
【０００３】
そういったモデル構築のためには生体内の各組織の詳細な構造データと、その力学的特性の収集が必要であり、構築過程は以下の３つに大別される。
▲１▼生体内情報を画像データとして収集
▲２▼認識（各組織のセグメンテーション）
▲３▼各組織の力学的特性を付加
【０００４】
これらの過程を経て構築される力学的モデルに、有限要素法（ＦＥＭ）などによる力学的シミュレーションが実行されるが、そのシミュレーション結果は特に画像データの精度に大きく左右される。そのため、画像データ収集のための生体内撮影装置の撮影精度が大きな問題となる。現在、主な撮影装置としてＸ線ＣＴ（Ｃｏｍｐｕｔｅｄ  Ｔｏｍｏｇｒａｐｈｙ）、ＭＲＩ（Ｍａｇｎｅｔｉｃ  Ｒｅｓｏｎａｎｃｅ  Ｉｍａｇｉｎｇ）が用いられているが、その撮影特徴上、データ収集が可能な組織が限られる。つまり、詳細部位や軟組織のデータ収集は不可能である。また、撮影された組織の解像度も十分とは言えず、現在のところ意図するモデル構築には至らない。
【０００５】
生体内の詳細なモデルを構築するためには、構成要素の組成が似ている軟組織の情報も入手し、認識することが必要不可欠であるが、現在の主に用いられている撮影装置からでは要求を満たすデータが入手できていないのが現状であった。
【０００６】
ところが近年、生体内組織の詳細データの入手に、実際の生体組織が持つカラー情報を利用しようとする動きが出てきた。米国Ｎａｔｉｏｎａｌ  Ｌｉｂｒａｒｙ  ｏｆ  Ｍｅｄｉｃｉｎｅとコロラド大学によるＶｉｓｉｂｌｅ  Ｈｕｍａｎ  Ｐｒｏｊｅｃｔ  （以下、ＶＨＰ）では人体全身を約０．３３ｍｍ間隔でスライスし、そのスライス面を撮影することで生体内フルカラー連続断面画像の入手を行った。
【０００７】
Ｘ線ＣＴ、ＭＲＩが数ｍｍの解像度であることや、入手データがカラー情報を持つことを考慮すると、その有用性は言うまでもない。しかし、このプロジェクトでは人体を一旦特殊冷却保存した後に、まず頭部から足先の方向に０．３ｍｍ間隔で削り取り、その度毎に分解能２０００×２０００画素のカラーディジタルカメラで撮影してデータを得ているため、その手間を考慮するとデータを量産出来るものではない。
【０００８】
一方、本件出願人である理化学研究所において、生体内カラー情報の収集が可能な装置として三次元内部構造顕微鏡が開発された。この装置は数十μｍの極小間隔で生体試料を全自動でスライスすることが出来る装置であり、スライス面の直上からＣＣＤカメラで直接スライス面を撮影することで、軸のずれない高精彩のフルカラー連続断面画像の入手を可能としている。僅か１時間程度で１万枚程度の連続画像が入手可能である。
【０００９】
さて、上述の装置で入手した生体内フルカラーデータを利用し、生体内ディジタルモデルを構築したいが、ここで新しく問題が生じた。従来は白黒画像であったものが、フルカラーになったことで画像の特徴が大きく変わったこと、さらに、生体内のデータを詳細に入手したため情報量が莫大に増加し、前述の過程▲２▼において従来のセグメンテーションでは良い認識結果を示すことが出来なくなったのである。
【００１０】
もともと連続断面画像のような３次元画像からの領域抽出（セグメンテーション）法は、それだけで一つの研究分野が出来るほど奥が深い。従来多く用いられているＸ線ＣＴ、ＭＲＩ画像からのセグメンテーション法でも未だ確立された手法はなく、多くの議論がなされている。
【００１１】
一般的にセグメンテーションでは、「画像内の、一つの対象物（関心組織）に対応した部分（関心組織領域）では、特徴（濃度値、テクスチャなど）がほぼ均一であり、異なる対象物間（組織間）の境界部分では特徴が急変している」という仮定が利用されている。この仮定に基づき、セグメンテーション手法は２つに大別できる。すなわち、（１）関心組織の輪郭に相当するエッジを線として抽出する手法（エッジ抽出法）および、（２）画像を関心組織に対応した部分領域に分割する手法（領域分割法）である。それぞれの代表的なものとして、（１）ではスネークやレベルセットなど、（２）では領域成長法（リージョングローイング）などが挙げられる。
【００１２】
【非特許文献１】
鳥脇純一郎著「３次元ディジタル画像処理」昭晃堂
【非特許文献２】
周藤安造著「医学における三次元画像処理基礎から応用まで」コロナ社
【００１３】
【発明が解決しようとする課題】
上述の（１）のエッジ抽出法の場合、局所処理を用いて特徴の急変性を強調し境界を検出するため、画像内に雑音（ノイズ）が含まれる場合や、関心組織形状が複雑化した場合は、偽のエッジやエッジの途切れが生じる。そのため、生体フルカラー画像の様に保持する情報が多く、エッジ特定のための必要特徴とそれ以外を判別することが難しいことや、生体組織形状の複雑さを考慮すると、最適なセグメンテーション手法には挙げられない。
【００１４】
一方、前述の（２）の領域分割法では、関心組織領域内での特徴の均一性を評価して関心組織に対応すると考えられる領域を直接抽出するため、雑音などの局所的な変動に左右されにくい。さらに、抽出結果は領域としてまとまっているために画像の大局的な構造を理解しやすく、新たに構造理解の必要が減るケースが多いなどの利点がある。しかし、ごく狭い領域（線状領域）を抽出する場合や、組織間の特徴変動がなだらかな場合に、分離すべき領域が統合される（抽出のあふれ）がしばしば生じる欠点もあり、多くの場合、何らかの処理を恣意的に加えなければならないケースが多い。
【００１５】
また、前述の（１）と（２）双方の対照的、相補的な特徴を利用し、それらを組み合わせることで性能を良くしたセグメンテーション手法として領域・エッジ併用法が提案されている。  この併用方法はリージョングローイングによる領域分割の判定条件にエッジ情報を使用するため、エッジ抽出法だけのときのエッジの途切れをそれほど問題としない。さらに、リージョングローイングだけのときよりも、エッジ情報が加わったために正しい判別が行える利点がある。
【００１６】
しかし、リージョングローイングは領域判別の制御条件としてあらかじめパラメータを与える必要がある。さらに、生成される領域の拡張量に停止条件を設定しないと処理を終了することが出来ない。このパラメータや条件を自動的に画像の特徴に合わせて調整することは難しく実現されていないため、複雑な画像からの正しいセグメンテーション結果は期待出来ない。
【００１７】
また他に、各判別段階で対話的な修正処理を組み込むことで、画像特徴に合わないパラメータの役割を補間する手法が提案されている。これは、領域抽出過程のモニタリングを行い、領域拡張の停止条件適当でないために生じた抽出のあふれを、解剖学的知識を持った人が指定し修正することで正しいセグメンテーションを実現する手法である。
【００１８】
前述のように、データの情報量が莫大に増加した画像において、人の手作業による修正処理を施す必要のあるこの方法では一般性は得られない。
【００１９】
以上のように、リージョングローイングはエッジ抽出法と比較し利点も多いが、その効果が発揮されるのは、理想的なデータに対し、特定のパラメータや条件を与えられた場合のみであり、通常の処理の場合は人の介在なくしては良い結果が得られない。今日の研究で扱う生体内部のフルカラー画像は、形状や濃度変化も複雑でデータの総容量も多く、従来法をそのまま適用しても良いセグメンテーション結果は得られない。
【００２０】
そこで、本件発明者は鋭意研究を通じて、生体内部のフルカラー画像に適したセグメンテーション手法を提案するに至った。
【００２１】
本発明の目的は、従来法のリージョングローイングで問題とされた、恣意的なパラメータや拡張の停止条件をなくし、自動的なセグメンテーションを実現することが可能な関心組織領域抽出方法、関心組織領域抽出処理プログラム及び画像処理装置を提供することにある。
【００２２】
【課題を解決するための手段】
上述した課題を解決し、上記目的を達成するため、請求項１の発明に係る関心組織領域抽出方法は、生体断面画像の連続フレームから関心組織領域を抽出する方法であって、前記連続フレームの一部に対して前記関心組織領域における初期の判定基準を与え、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記判定基準とに基づいて前記関心組織領域に帰属させる領域と帰属させない領域とを判断し、前記フレームに続くフレーム以降の判断では、前記関心組織領域に前記帰属させる領域を含めた領域の画素値に基づいて前記判定基準を設定することを特徴とするものである。
【００２３】
また、請求項２の発明に係る関心組織領域抽出方法は、生体断面画像の連続フレームから関心組織領域を抽出する方法であって、前記連続フレームの一部に対して前記関心組織領域内における初期の判定基準と前記関心組織領域外における初期の判定基準の２種類を与え、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記２種類の判定基準とに基づいて前記関心組織領域内に帰属させる領域と前記関心組織領域外に帰属させる領域とを判断し、前記フレームに続くフレーム以降の判断では、前記関心組織領域内に前記帰属させる領域を含めた領域の画素値と前記関心組織領域外に前記帰属させる領域を含めた領域の画素値に基づいて前記２種類の判定基準を設定することを特徴とするものである。
【００２４】
上記請求項１又は２の発明において、前記連続フレームの一部は１枚のフレーム又は複数枚のフレームであり、前記フレームに対して手操作により設定された前記関心組織領域から前記判定基準を求めるようにしてもよい。
【００２５】
上記請求項１又は２の発明において、前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレームの１つ前のフレームで判断された関心組織領域と同じ位置に仮領域を設け、当該仮領域の境界に沿って画素単位に帰属の判断を行うようにしてもよい。
【００２６】
上記請求項１又は２の発明において、前記仮領域において、前記関心組織領域に帰属させる画素が前記仮領域外にある場合には、当該画素を前記判断対象のフレームの関心組織領域の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して拡張し、前記関心組織領域に帰属させない画素が前記仮領域内にある場合には、当該画素を前記判断対象のフレームの関心組織領域外の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して縮小するようにしてもよい。
【００２７】
上記請求項１又は２の発明において、前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレーム以前の連続する複数枚のフレームから画素単位に対応させて前記判断対象となるフレームの各画素を含む所定サイズの局所領域を取り出しながら、前記各局所領域に含まれる前記関心組織領域の一部に対応する前記設定された判定基準を使用するようにしてもよい。
【００２８】
上記請求項１又は２の発明において、前記関心組織領域の画素のみ表示対象として記憶するようにしてもよい。
【００２９】
また、請求項８の発明に係る関心組織領域抽出プログラムは、生体断面画像の連続フレームから関心組織領域を抽出するための処理をコンピュータに実行させるプログラムであって、前記コンピュータに、前記連続フレームの一部に対応する前記関心組織領域における初期の判定基準を算出する第１ステップと、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記算出された初期の判定基準とに基づいて前記関心組織領域に帰属させる領域と帰属させない領域とを判断する第２ステップと、前記フレームに続くフレーム以降の判断では、前記関心組織領域に前記帰属させる領域を含めた領域の画素値に基づいて前記判定基準を設定する第３ステップと、を実行させることを特徴とするものである。
【００３０】
また、請求項８の発明に係る関心組織領域抽出プログラムは、生体断面画像の連続フレームから関心組織領域を抽出するための処理をコンピュータに実行させるプログラムであって、前記連続フレームの一部に対して前記関心組織領域内における初期の判定基準と前記関心組織領域外における初期の判定基準の２種類を算出する第１ステップと、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記算出された２種類の判定基準とに基づいて前記関心組織領域内に帰属させる領域と前記関心組織領域外に帰属させる領域とを判断する第２ステップと、前記フレームに続くフレーム以降の判断では、前記関心組織領域内に前記帰属させる領域を含めた領域の画素値と前記関心組織領域外に前記帰属させる領域を含めた領域の画素値に基づいて前記２種類の判定基準を設定する第３ステップと、を実行させることを特徴とするものである。
【００３１】
上記請求項８又は９の発明において、前記連続フレームの一部は１枚のフレーム又は複数枚のフレームであり、前記フレームに対して手操作により設定された前記関心組織領域から前記初期の判定基準を算出させる前記第１ステップを有したものでもよい。
【００３２】
上記請求項８又は９の発明において、前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレームの１つ前のフレームで判断された関心組織領域と同じ位置に仮領域を設け、当該仮領域の境界に沿って画素単位に帰属の判断を行わせる前記第２ステップを有したものでもよい。
【００３３】
上記請求項８又は９の発明において、前記仮領域において、前記関心組織領域に帰属させる画素が前記仮領域外にある場合には、当該画素を前記判断対象のフレームの関心組織領域の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して拡張させ、前記関心組織領域に帰属させない画素が前記仮領域内にある場合には、当該画素を前記判断対象のフレームの関心組織領域外の画素として当該関心組織領域を前記１つ前のフレームの関心組織領域に対して縮小させる前記第２ステップを有したものでもよい。
【００３４】
上記請求項８又は９の発明において、前記連続フレームの一部に続くフレーム及び当該フレームに続くフレーム以降の帰属の判断では、判断対象となるフレーム以前の連続する複数枚のフレームから画素単位に対応させて前記判断対象となるフレームの各画素を含む所定サイズの局所領域を取り出させながら、前記各局所領域に含まれる前記関心組織領域の一部に対応する前記設定された判定基準を使用させる前記第２ステップを有したものでもよい。
【００３５】
上記請求項８又は９の発明において、前記コンピュータに、さらに、前記関心組織領域の画素のみ表示対象として記憶する第４ステップを実行させるものでもよい。
【００３６】
また、請求項１５の発明に係る画像処理装置は、生体断面画像の連続フレームから関心組織領域を抽出する画像処理装置であって、前記連続フレームの画素値を記憶する画素値記憶手段と、前記画素値記憶手段を参照して前記連続フレームの一部に対応する前記関心組織領域における初期の判定基準を算出する算出手段と、前記連続フレームの一部に続くフレームについて当該フレームの各画素の濃度値と前記算出された初期の判定基準とに基づいて前記関心組織領域に帰属させる領域と帰属させない領域とを判断する判断手段と、前記フレームに続くフレーム以降の判断では、前記関心組織領域に前記帰属させる領域を含めた領域の画素値に基づいて前記判定基準を設定する設定手段と、を備えることを特徴とするものである。
【００３７】
また、請求項１６の発明に係る画像処理装置は、生体断面画像の連続フレームから関心組織領域を抽出する画像処理装置であって、前記連続フレームの画素値を記憶する画素値記憶手段と、前記画素値記憶手段を参照して前記連続フレームの一部に対して前記関心組織領域内における初期の判定基準と前記関心組織領域外における初期の判定基準の２種類を算出する算出手段と、前記連続フレームの一部に続くフレームについて当該フレームの各画素値と前記算出された２種類の判定基準とに基づいて前記関心組織領域内に帰属させる領域と前記関心組織領域外に帰属させる領域とを判断する判断手段と、前記フレームに続くフレーム以降の判断では、前記関心組織領域内に前記帰属させる領域を含めた領域の画素値と前記関心組織領域外に前記帰属させる領域を含めた領域の画素値に基づいて前記２種類の判定基準を設定する設定手段と、を備えることを特徴とするものである。
【００３８】
上記請求項１６又は１７の発明において、さらに、前記関心組織領域の各画素に対応させて表示対象を示すフラグ情報を記憶する表示対象記憶手段を有するようにしてもよい。
【００３９】
上記請求項１６又は１７の発明において、さらに、前記画素値記憶手段を参照し、前記表示対象記憶手段に記憶されたフラグ情報に基づいて前記関心組織領域を表示する表示手段を有するようにしてもよい。
【００４０】
【発明の実施の形態】
以下に添付図面を参照して、本発明にかかる一実施の形態を詳細に説明する。
【００４１】
まず、本発明の原理について説明する。本発明に係る関心組織領域抽出方法、関心組織領域抽出プログラム及び画像処理装置において共通する原理は、生体断面画像の連続フレームの一部に対して関心組織領域内における初期の判定基準（濃度中央値又は濃度平均値）と関心組織領域外における初期の判定基準（濃度中央値又は濃度平均値）の２種類を与えることろからはじめる。
【００４２】
連続フレームの一部に連続するフレームについて、そのフレームの各画素値（画素の濃度値）は関心組織領域内における初期の判定基準、関心組織領域外における初期の判定基準にそれぞれ比較される。画素値が関心組織領域内における初期の判定基準に近い場合には、関心組織領域内に帰属させる領域として判断される。
【００４３】
この場合、上記の関心組織領域とその領域内に帰属させると判断した領域とを含めた新たな関心組織領域（３次元構造）の画素値に基づいて判定基準を設定する。
【００４４】
一方、画素値が関心組織領域外における初期の判定基準に近い場合には、関心組織領域外に帰属させる領域として判断される。この場合、上記の関心組織領域外とその領域外に帰属させると判断した領域とを含めた新たな関心組織領域外（３次元構造）の画素値に基づいて判定基準を設定する。
【００４５】
さらに連続する次のフレームについては、連続する前処理済みフレームによって更新された判定基準を用いて関心組織領域内か外かを判断する。
【００４６】
つづいて、フルカラー連続断面画像からの自動領域抽出についての原理を説明する。図１はＨＳＶカラー空間を説明する図、図２は本発明に適用される極座標を説明する図、図３及び図４は本発明に適用される領域抽出方法を説明する図、そして、図５は本発明に適用される３次元空間における関心組織領域のリージョングローイング方法を説明する図である。
【００４７】
フルカラー画像は各ピクセルにＲＧＢ表色系による３つのカラーベクトル値をもっている。各２次元画像において、関心組織領域を抽出するために、組織の境界を判別する必要がある。そのためには、組織間のわずかな色の差を認識しなくてはならない。つまり、近隣画素間での微小な色変化に注目をする。
【００４８】
関心組織領域を抽出する際には、画素値に閾値を適用する方法を用いることもできるが、生体内では組織境界付近に明らかな色変化があることが少なく、また場所ごとに複雑な色分布をしているため、一意的な閾値では関心組織領域を抽出することは困難である。
【００４９】
そこで本発明では、以下の２点を抽出手法に取り入れる。
（１）ＨＳＶ（Ｈ：色相、Ｓ：色彩、Ｖ：明度）表色法で表現した画素値を利用した画素間濃度差の算出
（２）画像全体に同一閾値を適用するのではなく、局所的領域に注目した領域内外判別（拡張リージョングローイング法）
【００５０】
ＨＳＶ表色法を用いた理由は以下の２点となる。
（１）ＲＧＢ表色法よりも生体組織境界を判別しやすい点
（２）同一組織内では特にＨ値、Ｓ値に大きな色変化がおこらず値がまとまっている点
【００５１】
さて、ＨＳＶ表色法で表された画素値を用いてセグメンテーションを行うが、特にＨ値、Ｓ値を用いてセグメンテーションを行う。しかし、ＨＳＶ表色法で表されるＨ，Ｓ，Ｖの各ベクトル値は互いに等価ではなく、図１のような特殊な空間を生成しているため、画素間濃度差を計算する際、通常のユークリッド距離では正しい値が得られない。
【００５２】
そこで、ある２点の画素間濃度差を算出する際は、ＨＳＶ空間の特徴を損なわない計算方法として、Ｈ，Ｓによる極座標上での距離を計算した。例として画素Ｐ１，Ｐ２の画素値をそれぞれ
【数１】
2004234579.tif 000003
としたとき、画素  Ｐ１，Ｐ２  間の濃度差δは以下の式（１），（２），（３）で計算する。
【数２】
2004234579.tif 000004
【００５３】
本発明は、以上の画素間の濃度差を求める手法を３次元の局所領域に適用して関心組織領域に帰属するか否かを判断するものである。そこで、図３、図４及び図５を用いて簡単に説明する。
【００５４】
連続する生体断面画像（生体フルカラー画像情報の連続フレーム）から初期画像（フレームＮ−１）を選択する際は、ある程度視認できる程度のサイズを有した関心組織領域が写っていることが好ましい（図３参照）。この場合には、手作業による領域抽出となるが、関心組織領域を含むことが精度上好ましい。なお、図３の例では、フレームＮ−１，Ｎ，Ｎ＋１・・・の方向がスライス方向であり、関心組織領域の中心に向けてスライス方向が進むものとする。
【００５５】
フレームＮ−１の画像に手作業で抽出対象と思われる任意の領域が指定されると、その領域は関心組織領域ＲＯＩ０として抽出される（図４（Ａ））。そして、連続する生体断面画像では、連続する次のフレームＮの画像に対して上フレームＮ−１の画像で抽出された領域とまったく同じ位置の領域をまずは仮領域Ｓとして選択される（図４（Ｂ）参照）。
【００５６】
その後、後述する本発明の拡張リージョングローイング法により、正しい領域すなわち実際の関心組織領域ＲＯＩ１をが仮の状態から再抽出される（図４（Ｃ）参照）。このようにして関心組織領域ＲＯＩ１の抽出が完了する（図４（Ｄ）参照）。
【００５７】
拡張リージョングローイング法とは、フレームＮ−１の画像から抽出された関心組織領域ＲＯＩ０と全く同じ抽出対象の関心組織領域ＲＯＩ１をＨ値、Ｓ値に基づく濃度値を用いてフレームＮの画像から抽出し、さらに連続する次のフレームＮ＋１、Ｎ＋２、Ｎ＋３・・・については少なくとも過去フレームの関心組織領域内外の濃度値から抽出する手法である。
【００５８】
すなわち、各フレームにおいては、１つ前のフレームで抽出された関心組織領域と同じ位置の領域をまずは仮領域として与え、その仮領域の境界画素に一つ一つ注目し、断面（連続するフレーム）間で領域内と領域外のどちらかに変化したかを判定する処理が施される。
【００５９】
図５において、まず境界画素Ｐｘを中心とした７ｘ７ｘ４ピクセルの局所領域に注目する。この場合、Ｚ方向はＮ，Ｎ−１，Ｎ−２，Ｎ−３で示すが、領域抽出の処理方向は図中下方向となる。
【００６０】
その局所領域中、この時点までに領域内と判定されている画素のＨ値，Ｓ値の濃度中央値（メディアン値）Ｍｉｎを計算する。同時にこの時点までに領域外と判定されている画素のＨ値、Ｓ値の中央値Ｍｏｕｔも計算する。
【数３】
2004234579.tif 000005
【００６１】
画素Ｐｘの画素値が、この局所領域の領域内外どちらの特徴に近いかを検討し領域内，外を決定する。実際には、領域内外のＨ値、Ｓ値による濃度代表値Ｍと画素Ｐｘとの濃度差ｄをそれぞれ計算し、差が小さかった方と画素Ｐｘの特徴が類似しているとみなす。これにより、Ｐｘを領域内、領域外のどちらかに取り込む。これを計算式で表したものが以下の式（５）〜（６）である。
【数４】
2004234579.tif 000006
【００６２】
以上の処理をすべての境界画素に対して行う。この処理は領域内外への取り込みが終了するまで自動的に繰り返す。このときの取り込みが終了するとは、画素の領域内への取り込み、または除外を繰り返すことで正しい領域境界の設定が完了したことを意味する。
【００６３】
このようにして、フレームＮの画像に対する抽出処理が終了した後、さらに下に連続するフレームＮ＋１，Ｎ＋２，・・・の画像においても同様の処理を行い、連続自動抽出を実行する。
【００６４】
以上の原理に基づいて以下に構成、動作を説明する。まず、構成について説明する。図６は本発明の一実施の形態による画像処理装置の一構成を示すブロック図である。図６の画像処理装置は本発明に係る画像処理装置を構成しており、フロッピー（登録商標）ディスク、ＣＤ、ＤＶＤ等の記録媒体や通信回線を通じて図６の画像処理装置に供給されるプログラムは本発明に係る関心組織領域抽出プログラムを含むものとする。
【００６５】
本実施の形態による画像処理装置は、例えば図６に示したように、内部バス１１に、通信インタフェース１２、ＣＰＵ１３、ＲＯＭ１４、ＲＡＭ１５、ディスプレイ１６、キーボード／マウス１７、ドライブ１８、ハードディスク１９を接続させ、アドレス信号、制御信号、データ等を伝送させ、本実施の形態による領域抽出を実現する構成を備えている。
【００６６】
図６において、通信インタフェース１２は、例えばインターネット等の通信網に接続する機能を有しており、本発明に係る関心組織領域抽出プログラムをダウンロードすることも可能である。ＣＰＵ１３は、ＲＯＭ１４に格納されたＯＳにより装置全体の制御を行うとともにハードディスク１９に格納された各種のアプリケーションプログラムに基づいて処理を実行する機能を司る。
【００６７】
ＲＯＭ１４は、ＯＳ等のように装置全体の制御を行うためのプログラムを格納しており、これらをＣＰＵ１３に供給する機能を有している。ＲＡＭ１５は、ＣＰＵ１３による各種プログラムの実行時にワークエリアとして利用されるメモリ機能を有している。このＲＡＭ１５は、領域抽出の際に例えばスタックＡ，Ｂを設定するものとする。
【００６８】
ディスプレイ１６は、ＣＰＵ１３の各種の処理に伴うメニュー、ステータス、表示遷移等を表示する機能を有している。キーボード／マスス１７は、文字、数字、記号等のデータを入力するためのキー、画面上のポイント位置を操作するためのカーソルキー、ファンクションキー等を備えている。なお、ポイント位置を操作するために、マウスが備わっている。
【００６９】
ドライブ１８は、各種のプログラム、データを記録したＣＤ、ＤＶＤ等の記録媒体からインストール作業を実行するための駆動ユニットである。なお、フロッピー（登録商標）ディスクのドライブ機能を付加してもよい。
【００７０】
ハードディスク１９は、プログラム１９Ａ、メモリ１９Ｂ、生体フルカラー画像情報１９Ｃ、フラグ情報１９Ｄ等を記憶する外部記憶装置である。プログラム１９Ａは、前述した通信インタフェース１２、ドライブ１８等からインストールされたプログラムを実行形式で記憶したものに相当する。メモリ１９Ｂは、各種プログラムの実行結果等のファイルを保存する記憶部である。
【００７１】
生体フルカラー画像情報１９Ｃは、通信インタフェース１２、ドライブ１８等を介して読み込んだデータファイルである。この生体フルカラー画像情報１９Ｃは、連続する複数フレームにより構成される生体の断面画像である。フラグ情報１９Ｄは、生体フルカラー画像情報１９Ｃの各フレームに座標で対応させてＲＧＢをひとつのフラグで表すものである。
【００７２】
つづいて、本実施の形態を機能的に説明する。図７は本実施の形態による機能を示すブロック図であり、図８は本実施の形態による表示遷移の一例を示す図である。本発明の機能は、初期関心組織領域設定部１０１、仮領域設定部１０２、濃度中央値算出部１０３、境界近傍画素値取得部１０４、領域内外判定部１０５、関心組織領域抽出部１０６、表示画像形成部１０７等により構成される。
【００７３】
以上の構成によれば、まず初期関心組織領域設定部１０１により生体フルカラー画像情報１９ＣからフレームＮ−１の画像が取り出され、手作業による任意の領域指定により初期の関心組織領域が設定される。濃度中央値算出部１０３では、この関心組織領域設定部１０１で設定された領域内の各画素の濃度ちから濃度中央値が算出される。
【００７４】
仮領域設定部１０２により、フレームＮ−１に連続する次のフレームＮに対して初期関心組織領域設定部１０１で設定された関心組織領域と同じ領域の位置に仮領域が設定される。そして、境界近傍画素値取得部１０４により仮領域内外の境界近傍画素の画素値が取得される。
【００７５】
領域内外判定部１０５では、境界近傍画素と画素値と濃度中央値算出部１０３で算出された濃度中央値とに基づいて前述したように拡張リージョングローリング法による領域内外の判定が行われる。なお、７×７×４の局所領域での判定が開始されるのは、過去に４フレーム分の領域抽出が完了してからとなる。よって、抽出開始のレームをＮとした場合には、フレームＮ＋３から７×７×４の局所領域が採用される。
【００７６】
そして、関心組織領域抽出部１０６により、領域内外判定に基づいて関心組織領域の拡張又は伸縮による抽出が行われ、その抽出結果がフラグ情報１９Ｄに反映される。
【００７７】
連続する次のフレームＮ＋１においては、フレームＮ−１及びＮによる３次元の局所領域が濃度中央値算出部１０３で濃度中央ちを算出する際の対象となる。このとき、仮領域設定部１０２では、フラグ情報１９ＤよりフレームＮの関心組織領域を取り込み、その関心組織領域にしたがって仮領域が設定される。以降、同様に、拡張リージョングローリング法により領域抽出が実施される。
【００７８】
そして、抽出済の関心組織領域を表示させる場合には、表示画像形成部１０７によりフラグ情報１９Ｄ及び生体フルカラー画像情報１９Ｃに基づいて表示データが生成され、図８（Ａ）〜（Ｆ）の如く関心組織領域ＲＯＩの表示遷移を観察することができる。
【００７９】
つぎに、動作について簡単に説明する。図９は本実施の形態による動作を説明するフローチャートであり、図１０は本実施の形態によるフラグ情報の管理について説明する図である。なお、この図９による処理は、プログラム１９Ａのうち本発明に係る関心組織領域抽出プログラムに相当するものである。ただし、表示については、図７の機能ブロックに従うものである。
【００８０】
まず、１枚目の断面画像が入力され（ステップＳ１）、ＨＳＶ変換処理が実行される（ステップＳ２）。そして、抽出対象の領域が指定される（ステップＳ３）。これにより手作業部分での初期の関心組織領域が設定される。
【００８１】
この段階では次の断面画像が存在することから（ステップＳ４）、連続する次の断面画像が入力される（ステップＳ５）。このとき、入力された断面画像に対応するフラグ情報がリセットされる（ステップＳ６）。
【００８２】
入力断面画像に対しては、１つ前の断面画像に指定した関心組織領域と同じ領域の位置に仮領域として抽出開始領域が指定される（ステップＳ７）。この抽出開始領域境界の近傍画素を１つ選択され（ステップＳ８）、スタックＡ（図６参照）に格納される（ステップＳ９）。
【００８３】
この段階では、スタックＡは空ではないので（ステップＳ１０）、スタックＡから注目画素Ｐを１つ取り出し（ステップＳ１１）、局所領域による前述した拡張リージョングローリング処理が実行される（ステップＳ１２）。
【００８４】
そして、注目画素Ｐが関心組織領域内の場合には（ステップＳ１３）、注目画Ｐに対してフラグがセットされ（ステップＳ１４）、その注目画素Ｐ近傍で抽出開始領域に対して外側の画素が１つ選択される（ステップＳ１５）。この後、選択画素はスタックＡに格納され（ステップＳ９）、同様の処理が繰り返される。
【００８５】
フラグ情報は、図１０に示したように、核フレームに対して画素毎に座標で管理される。例えば、フレームＮｏ．１００の座標（Ｘ１，Ｙ１，Ｚ１）が関心組織領域内と判断された場合には、フラグ“１”にセットされる。フレームＮｏ．１０１の座標（Ｘ２，Ｙ２，Ｚ２）はフラグ“１”のため、関心組織領域内と判断されたことになる。
【００８６】
ステップＳ１３において注目画素Ｐが領域外の場合には、注目画素Ｐの近傍で領域内の画素が１つ選択され（ステップＳ１６）、その選択画素はスタックＢに格納される（ステップＳ１７）。そして、スタックＡの空状態が確認される（ステップＳ１８）。このスタックＡに関する処理は領域内を処理対象としている。
【００８７】
スタックＡが空でないときには（ステップＳ１８のＮＯルート）、処理はステップＳ１１に戻り、同様の処理が繰り返される。一方、スタックＡが空の場合には（ステップＳ１８のＹＥＳルート）、さらにスタックＢの空の状態が確認される（ステップＳ１９）。
【００８８】
スタックＢが空でないときには（ステップＳ１９のＮＯルート）、処理はステップＳ２０に進み、一方、スタックＢが空の場合には（ステップＳ１９のＹＥＳルート）、入力断面画像１枚分の処理が完了したので、処理はステップＳ４に戻り、次の断面画像の有無が判断される。
【００８９】
ステップＳ２０においては、スタックＢから注目画素Ｐとして画素が１つ取り出され、その注目画素Ｐに対して前述した拡張リージョングローリング処理が施される（ステップＳ２１）。そして、その注目画素Ｐが領域外であれば（ステップＳ２２）、さらに内側の判断を仰ぐため、注目画素Ｐの近傍で領域内の画素が１つ選択される（ステップＳ２３）。その選択画素はスタックＢに格納され（ステップＳ１７）、以下は同様の処理が繰り返される。このスタックＢに関する処理は領域外を処理対象としているので、ステップＳ２２で注目画素が領域内と判断された場合には処理はステップＳ１９に戻る。
【００９０】
以上説明したように、本実施の形態によれば、拡張リージョングローリングにより、注目画素周囲の３次元局所空間を使用し、注目画素周辺の３次元局所空間での関心組織領域内外の濃度中央値を判定基準とした領域抽出を行うようにしたので、抽出対象のフレームよりも前の連続する数フレーム分の判定結果を加味することができ、その結果、拡張リージョングローリングによる自動的なセグメンテーションを実現することが可能である。このようにすれば、例えば同一臓器の色の変化にも対応して正確に領域抽出を行い、また、同一臓器内に異物が混在しても影響されずに正確に領域抽出を行うことができる。また、従来法のリージョングローイングで問題とされた恣意的なパラメータや拡張の停止条件をなくし、判定基準の固定化を避けることができる。
【００９１】
さて、上述した実施の形態では、判定基準のデータとして濃度中央値を採用したが、本発明はこれに限定されるものではなく、濃度中央値のほうが好ましいが、濃度平均値からも同様の効果を得ることが可能である。
【００９２】
また、上述した実施の形態では、３次元局所領域を７×７×４のサイズとして説明したが、本発明はこれに限定されるものではなく、例えば３×３×２のように上記サイズよりも小さくてもよく、あるいは、大きくてもよい。また、大中小等の複数のサイズをもつことで、それぞれに重み付けをすることも可能である。
【００９３】
さらに、上述した実施の形態では、カラー画像を例に挙げているが、白黒画像への適用も可能である。
【００９４】
また、上述した実施の形態では、初期画像にあたるフレームＮ−１に対して手作業で関心組織領域を抽出していたが、本発明はこれに限定されるものではなく、画像処理技術により自動的に抽出するようにしてもよい。
【００９５】
つぎに、実験例を２つ挙げる。
（実験１：水晶体部の抽出）
図１１及び図１２は実験２の結果を示す図である。人眼全体から水晶体領域を自動抽出した。合計枚数１２０枚を利用した。水晶体付近は、周辺組織が透明であるため、断面画像では次断面での領域がすけて見える現象がおこっていた。そのため、手作業ではどこまでが正しい領域であるのかを判別するのが難しい。
【００９６】
図１２には、連続する複数フレームについて、本発明を適用して抽出された領域内外（ｉｎｓｉｄｅ  ａｒｅａ／ｏｕｔｓｉｄｅ  ａｒｅａ）の画素数（Ｒｅｓｕｌｔ）と実際に手作業で抽出した領域内外（ｉｎｓｉｄｅ  ａｒｅａ／ｏｕｔｓｉｄｅ  ａｒｅａ）の画素数（Ａｎｓｗｅｒ）との比較から、エラー画素数（Ｅ）と回答率（Ｒａｔｅ％）とを示している。フレームＮｏ．は１，５，１０，３０，５０，１００，１３０，１４０，１５０をサンプルとして挙げている。
【００９７】
例えば、フレームＮｏ．１０では、領域内は１６６３画素であるのに対して実際は１７２６画素となり、９６．３５％の正解率を得ることができた。このときのエラー画素数は、６３となった。一方、領域外は７５０１５画素であるのに対して実際は７５０７４画素となり、９９．９２％の正解率を得ることができた。このときのエラー画素数は、５２である。
【００９８】
このように、枚数が進んでも高い正解率が得られ、フレームＮｏ．１３０まで進んでもほぼ９５％（９４．８１％）の正解率を得ることができた。これは、５％のエラーを示すものであり、手作業との間で遜色のない領域抽出を達成できたことを示すものである。
【００９９】
ところが、画素値のＳ値は実際の水晶体領域と、透けて見えている領域の間に大きな値の変化があったため、本手法ではそれらを分離することに成功した。
【０１００】
（実験２：眼球の包埋剤からの抽出）
図１３は実験１の結果を示す図である。人眼の連続断面画像８４０枚を使用し、背景（包埋剤）から人眼だけを抽出した。包埋剤は青一色だが、同じ青でも場所ごとに色に差があり、一つの閾値で青だけを除外する手法では難しい。また、８４０枚もの多数枚から手作業で抽出するのは困難であった。本発明の手法により、図１３に示したように、自動的に連続抽出が可能となった。
【０１０１】
以上の実験例から示されるように、Ｈ値、Ｓ値を利用したセグメンテーションの有用性が示された。
【０１０２】
【発明の効果】
以上説明したように本発明によれば、従来法のリージョングローイングでは関心組織領域の判定基準が固定されていたのに対して、関心組織領域の判定基準を３次元の空間で行うことができ、これによって、自動的なセグメンテーションを実現することが可能な関心組織領域抽出方法、関心組織領域抽出処理プログラム及び画像処理装置を提供できるという効果を奏する。
【図面の簡単な説明】
【図１】ＨＳＶカラー空間を説明する図である。
【図２】本発明に適用される極座標を説明する図である。
【図３】本発明に適用される領域抽出方法を説明する図である。
【図４】本発明に適用される領域抽出方法を説明する図である。
【図５】本発明に適用される３次元空間における関心組織領域のリージョングローイング方法を説明する図である。
【図６】本発明の一実施の形態による画像処理装置の一構成を示すブロック図である。
【図７】本実施の形態による機能を示すブロック図である。
【図８】本実施の形態による表示遷移の一例を示す図である。
【図９】本実施の形態による動作を説明するフローチャートである。
【図１０】本実施の形態によるフラグ情報の管理について説明する図である。
【図１１】実験１の結果を示す図である。
【図１２】実験１の結果を示す図である。
【図１３】実験２の結果を示す図である。
【符号の説明】
１２  通信インタフェース
１３  ＣＰＵ
１４  ＲＯＭ
１５  ＲＡＭ
１６  ディスプレイ
１７  キーボード／マウス
１８  ドライブ
１９  ハードディスク
１９Ａ  プログラム
１９Ｂ  メモリ
１９Ｃ  生体フルカラー画像情報
１９Ｄ  フラグ情報
１０１  初期関心組織領域設定部
１０２  仮領域設定部
１０３  濃度中央値算出部
１０４  境界近傍画素値取得部
１０５  領域内外判定部
１０６  関心組織領域抽出部
１０７  表示画像形成部
【図１】
2004234579.tif 000007
【図２】
2004234579.tif 000008
【図３】
2004234579.tif 000009
【図４】
2004234579.tif 000010
【図５】
2004234579.tif 000011
【図６】
2004234579.tif 000012
【図７】
2004234579.tif 000013
【図８】
2004234579.tif 000014
【図９】
2004234579.tif 000015
【図１０】
2004234579.tif 000016
【図１１】
2004234579.tif 000017
【図１２】
2004234579.tif 000018
【図１３】
2004234579.tif 000019