(57)【要約】【課題】        肺ガン等の異常陰影の判定を行う際の判定精度を向上させる。【解決手段】    医用画像から異常陰影の候補領域を検出して、候補領域の少なくとも一部を含む小領域を注目領域Mとして設定し、さらに、注目領域Mの近傍に存在する小領域を近傍領域R1,R2,・・・,R8として設定する。注目領域Mおよび近傍領域R1,R2,・・・,R8が正常な人工画像を生成し、注目領域Mの人工画像と原画像との差（最小自乗距離）と、近傍領域R1,R2,・・・,R8の人工画像と原画像との差（最小自乗距離）とを比較して、両差の違いの有無によって、注目領域Mに含まれる候補領域が異常陰影であるか否かの判定を行う。【選択図】          図２０
2006325937.tif 000002
【特許請求の範囲】
【請求項１】
被写体の所定の部位を撮影して得られた解剖学的特徴を有する構造物を含む医用画像において、該医用画像を細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類する分類手段と、
前記分類された小領域のうち少なくとも１つの小領域について、前記分類手段の分類結果に応じて前記構造物が正常な小領域の人工画像を生成する人工画像生成手段と、
前記医用画像から異常陰影の候補領域を検出する候補領域検出手段と、
該候補領域検出手段により検出された前記候補領域の少なくとも一部を含む小領域を注目領域として設定し、さらに、該注目領域の近傍に存在する小領域を近傍領域として設定する領域設定手段と、
前記分類手段により前記注目領域および前記近傍領域を分類し、分類された該各領域について前記人工画像生成手段により該両領域が正常な人工画像を生成し、該注目領域の人工画像と該注目領域の原画像との差と、前記近傍領域の人工画像と該近傍領域の原画像との差とを比較して、両差の違いの有無によって、前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行う判定手段とを備えることを特徴とする画像判定装置。
【請求項２】
前記判定手段が、前記注目領域の人工画像と該注目領域の原画像間の差分に応じた差分特徴量と、前記近傍領域の人口画像と該近傍領域の原画像間の差分に応じた差分特徴量を算出する差分特徴量算出手段を有するものであり、
該差分特徴量算出手段により算出した前記注目領域の差分特徴量と前記近傍領域の差分特徴量を比較して、両特徴量の違いの有無によって前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行うことを特徴とする請求項１記載の画像判定装置。
【請求項３】
前記差分特徴量は、前記人工画像と前記原画像との最小自乗距離または正規化相互相関値であることを特徴とする請求項２記載の画像判定装置。
【請求項４】
前記人工画像生成手段が、
前記分類された各小領域の構造物が正常であるときに該小領域毎に現れる固有の特徴を有する複数の固有画像を記憶する固有画像記憶手段を有し、
前記人工画像を生成する小領域と前記複数の固有画像それぞれとの類似度に応じた重み付けをして、該複数の固有画像を加算することにより該小領域の人工画像を生成するものであることを特徴とする請求項１から３いずれか記載の画像判定装置。
【請求項５】
前記固有画像記憶手段に記憶される固有画像が、被写体を撮影して得られた医用画像の前記構造物が正常な前記各小領域の画像に基づいて主成分分析して得た画像であることを特徴とする請求項４記載の画像判定装置。
【請求項６】
前記医用画像が、前記被写体の胸部を撮影して得た胸部画像であり、
前記分類手段が、
前記胸部画像より肺野領域を認識する肺野認識手段と、
前記胸部画像より肋骨を認識する肋骨認識手段とを有し、
肺野領域内の位置と肋骨の重なりに応じて前記多数の小領域を分類するものであることを特徴とする請求項１から５いずれか記載の画像判定装置。
【請求項７】
被写体の所定の部位を撮影して得られた解剖学的特徴を有する構造物を含む医用画像において、該医用画像を細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類する分類手段と、
前記分類された小領域のうち少なくとも１つの小領域について、前記分類手段の分類結果に応じて前記構造物が正常な小領域の人工画像を生成する人工画像生成手段とを備えた画像判定装置において、
前記医用画像から異常陰影の候補領域を検出する候補領域検出ステップと、
該候補領域検出ステップにより検出された前記候補領域の少なくとも一部を含む小領域を注目領域として設定し、さらに、該注目領域の近傍に存在する小領域を近傍領域として設定する領域設定ステップと、
前記分類手段により前記注目領域および前記近傍領域を分類し、分類された該各領域について前記人工画像生成手段により該両領域が正常な人工画像を生成し、該注目領域の人工画像と該注目領域の原画像との差と、前記近傍領域の人工画像と該近傍領域の原画像との差とを比較して、両差の違いの有無によって、前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行う判定ステップとを備えることを特徴とする画像判定方法。
【請求項８】
コンピュータを、
被写体の所定の部位を撮影して得られた解剖学的特徴を有する構造物を含む医用画像において、該医用画像を細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類する分類手段と、
前記分類された小領域のうち少なくとも１つの小領域について、前記分類手段の分類結果に応じて前記構造物が正常な小領域の人工画像を生成する人工画像生成手段と、
前記医用画像から異常陰影の候補領域を検出する候補領域検出手段と、
該候補領域検出手段により検出された前記候補領域の少なくとも一部を含む小領域を注目領域として設定し、さらに、該注目領域の近傍に存在する小領域を近傍領域として設定する領域設定手段と、
前記分類手段により前記注目領域および前記近傍領域を分類し、分類された該各領域について前記人工画像生成手段により該両領域が正常な人工画像を生成し、該注目領域の人工画像と該注目領域の原画像との差と、前記近傍領域の人工画像と該近傍領域の原画像との差とを比較して、両差の違いの有無によって、前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行う判定手段として機能させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、医用画像の読影による異常陰影候補の検出を支援する画像判定装置、画像判定装置方法、およびそのプログラムに関するものである。
【背景技術】
【０００２】
従来より、医療分野においては、ある患者についての撮影時期が異なる２つ以上の医用画像を比較読影して、画像間の差異を調べ、その差異に基づいて異常陰影を検出したり、疾患の進行状況や治癒状況を把握して治療方針を検討したりすることが行われている。
【０００３】
特に、異常陰影が初期の肺癌のような小型の陰影は見落としやすいため、このような異常陰影を検出する際には、撮影時期が異なる２つの医用画像について、画像間で画素を対応させた減算処理等を行って差異を求め、その差異を表す差異画像（経時サブトラクション画像）において、画素値が所定値以上の領域あるいは特徴的な形状、例えば略円形形状のものを異常陰影候補として検出する手法が提案されている（特許文献１等）。
【０００４】
しかし、比較対象となる撮影時期が異なる２つの時系列的な画像は、撮影時の患者の姿勢変化等により被写体の位置ずれが生じるため、差分画像を得る際に通常は比較対象となる画像の位置合わせが行なわれる。位置合わせを行う手法として、比較対象となる同一被写体についての２つの画像のうち少なくとも一方の画像全体に対して、回転、平行移動および拡大・縮小等の大局的位置合わせ（アフィン変換等の線形変換）や、局所領域ごとのテンプレートマッチングにより得られた対応位置関係に基づく非線形歪変換（例えば２次元多項式によるカーブフィッティングを用いた非線形歪み変換等）処理（ワーピング）による局所位置合わせ、またはこれら大局的位置合わせと局所位置合わせとを組み合わせて行う手法や、さらに位置ずれ程度の高い局所領域およびその近傍で再度位置合わせを行う手法が提案されている（特許文献２、特許文献３等）。
【０００５】
上記の手法は、あくまで２次元に投影された後の画像を用いて位置合わせを行うため、時系列的な画像間で患者の姿勢変化などにより大きく３次元的な変動（前傾、後傾、側旋）が生じた場合には、位置合わせしきれず、差異画像上にアーチファクト（偽像）が発生してしまい、異常陰影の検出が正確に行われなかった。また、過去に撮影された医用画像が存在しない患者の場合は、比較対象となる画像がなく、肺がんのような見落としやすい異常陰影の検出をする際に効果を発揮する差異画像による手法を上手く活用することができなかった。
【０００６】
そこで、被写体の正常な構造を表す正常構造画像が蓄積記憶されているデータベース中の各画像と、その各画像から抽出した解剖学的特徴の形状を教師データとして使用し、被写体の正常構造の統計的なモデル（形状変化モデル、テクスチャ変化モデル、形状とテクスチャの相関モデル）を作成し、このモデルを用いることにより、入力された画像に対応する正常構造画像を人工的に生成し、生成された正常構造画像と入力された医用画像とに基づいて異常陰影の候補領域を検出するものが提案されている（特許文献４等）。
【特許文献１】特開２００２‐１５８９２３号公報
【特許文献２】特開平７‐３７０７４号公報
【特許文献３】特開２００２‐３２７３５号公報
【特許文献４】特開２００４‐４１６９４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、上記の正常構造画像を人工的に生成する手法は、画像間で患者の姿勢変化などにより大きな３次元的変動が生じた場合や、過去に撮影された医用画像が存在しない患者の場合でも、見落としやすい小型の異常陰影の検出に効果を発揮するものであるが、通常、胸部CAD（computer aided diagnose）では検出された肺癌候補を部分的（ROI：region of interest）に切り出して、ROI毎に肺癌であるか否かの判別を行うことが行われていることからすると、胸部全体の人工画像を生成するのは非効率的である。
【０００８】
また、肋骨のような構造部が存在する部分と肋骨のような構造物が存在しない部分では現れる特徴が異なるため、胸部全体の人工画像を生成した場合には、部分的に元の画像とはずれた特徴を有したものが現れる場合がある。
【０００９】
そこで、胸部全体ではなく異常陰影が検出された注目する領域に対応した人工画像を精度よく生成して、肺がんのような見落としやすい異常陰影の検出を向上させることが望まれる。
【００１０】
本発明は、上記事情に鑑みてなされたものであり、比較読影による手法において従来に比べてより効率的でかつ精度よく異常陰影の判定をすることが可能となる画像判定装置、画像判定方法およびそのプログラムを提供することを目的とするものである。
【課題を解決するための手段】
【００１１】
本発明の画像判別装置は、被写体の所定の部位を撮影して得られた解剖学的特徴を有する構造物を含む医用画像において、該医用画像を細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類する分類手段と、
前記分類された小領域のうち少なくとも１つの小領域について、前記分類手段の分類結果に応じて前記構造物が正常な小領域の人工画像を生成する人工画像生成手段と、
前記医用画像から異常陰影の候補領域を検出する候補領域検出手段と、
該候補領域検出手段により検出された前記候補領域の少なくとも一部を含む小領域を注目領域として設定し、さらに、該注目領域の近傍に存在する小領域を近傍領域として設定する領域設定手段と、
前記分類手段により前記注目領域および前記近傍領域を分類し、分類された該各領域について前記人工画像生成手段により該両領域が正常な人工画像を生成し、該注目領域の人工画像と該注目領域の原画像との差と、前記近傍領域の人工画像と該近傍領域の原画像との差とを比較して、両差の違いの有無によって、前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行う判定手段とを備えることを特徴とするものである。
【００１２】
また、本発明の画像判別方法は、被写体の所定の部位を撮影して得られた解剖学的特徴を有する構造物を含む医用画像において、該医用画像を細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類する分類手段と、
前記分類された小領域のうち少なくとも１つの小領域について、前記分類手段の分類結果に応じて前記構造物が正常な小領域の人工画像を生成する人工画像生成手段とを備えた画像判定装置において、
前記医用画像から異常陰影の候補領域を検出する候補領域検出ステップと、
該候補領域検出ステップにより検出された前記候補領域の少なくとも一部を含む小領域を注目領域として設定し、さらに、該注目領域の近傍に存在する小領域を近傍領域として設定する領域設定ステップと、
前記分類手段により前記注目領域および前記近傍領域を分類し、分類された該各領域について前記人工画像生成手段により該両領域が正常な人工画像を生成し、該注目領域の人工画像と該注目領域の原画像との差と、前記近傍領域の人工画像と該近傍領域の原画像との差とを比較して、両差の違いの有無によって、前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行う判定ステップとを備えることを特徴とするものである。
【００１３】
また、本願発明のプログラムは、コンピュータを、被写体の所定の部位を撮影して得られた解剖学的特徴を有する構造物を含む医用画像において、該医用画像を細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類する分類手段と、
前記分類された小領域のうち少なくとも１つの小領域について、前記分類手段の分類結果に応じて前記構造物が正常な小領域の人工画像を生成する人工画像生成手段と、
前記医用画像から異常陰影の候補領域を検出する候補領域検出手段と、
該候補領域検出手段により検出された前記候補領域の少なくとも一部を含む小領域を注目領域として設定し、さらに、該注目領域の近傍に存在する小領域を近傍領域として設定する領域設定手段と、
前記分類手段により前記注目領域および前記近傍領域を分類し、分類された該各領域について前記人工画像生成手段により該両領域が正常な人工画像を生成し、該注目領域の人工画像と該注目領域の原画像との差と、前記近傍領域の人工画像と該近傍領域の原画像との差とを比較して、両差の違いの有無によって、前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行う判定手段として機能させることを特徴とするものである。
【００１４】
「被写体の所定の部位」とは、胸部や頭部や上腕部など人間や動物の身体などの中で、それがどの辺に位置するかという観点から見た時の部分をいう。
【００１５】
「構造物」とは、被写体の所定の部位を構成する組織をいう。
【００１６】
「構造物から得られる解剖学的特徴が類似する」とは、各組織が存在する場所や各組織の構成などから得られる特徴が類似することをいう。例えば、胸部は肋骨、軟部、心臓などの組織から構成されるが、肋骨が１本のみ撮影されている箇所と肋骨が２本重なって撮影されている箇所では画素値の濃度が異なり解剖学的特徴は類似しないが、肋骨が１本のみ撮影されている箇所同士は画素値の濃度が近い値となり解剖学的特徴が類似した画像となる。また、胸部の肺野領域のうち辺縁部であるか縦隔部であるかなどの位置の相違によっても、現れる濃度やテクスチャが異なり類似しない画像となるが、同じ辺縁部内であれば肋骨が１本重なっているか２本重なっているかにより類似した画像があらわれる。
【００１７】
また、「構造物が正常」とは、各部位を構成する組織が正常であることを意味する。
【００１８】
「候補領域の少なくとも一部を含む小領域を注目領域として設定する」とは、注目領域として設定する小領域内に、候補領域検出手段によって検出された候補領域が存在するように設定を行うことをいい、候補領域全体が小領域内に含まれず候補領域の一部が含まれるものであってもよい。
【００１９】
「注目領域の近傍に存在する小領域を近傍領域として設定する」は、近傍領域を注目領域に隣接するように小領域を設定するものであっても、注目領域と重なるように小領域を設定するものであってもよい。
【００２０】
「原画像」とは、撮影して得られた医用画像をいう。
【００２１】
また、前記判定手段は、前記注目領域の人工画像と該注目領域の原画像間の差分に応じた差分特徴量と、前記近傍領域の人口画像と該近傍領域の原画像間の差分に応じた差分特徴量を算出する差分特徴量算出手段を有するものであって、
該差分特徴量算出手段により算出した前記注目領域の差分特徴量と前記近傍領域の差分特徴量を比較して、両特徴量の違いの有無によって前記注目領域に含まれる前記候補領域が異常陰影であるか否かの判定を行うものを用いることができる。
【００２２】
また、前記差分特徴量は、前記人工画像と前記原画像との最小自乗距離または正規化相互相関値を用いるようにしてもよい。
【００２３】
また、前記人工画像生成手段は、前記分類された各小領域の構造物が正常であるときに該小領域毎に現れる固有の特徴を有する複数の固有画像を記憶する固有画像記憶手段を有すものであって、
前記人工画像を生成する小領域と前記複数の固有画像それぞれとの類似度に応じた重み付けをして、該複数の固有画像を加算することにより該小領域の人工画像を生成するものを用いることができる。
【００２４】
「固有画像」としては、例えば、複数の異なる正常構造の画像をサンプルとして、統計的な手法を用いて、構造物が正常であるときに各小領域毎に現れる特徴を有する画像を取得して固有画像とすることができる。
【００２５】
さらに、前記固有画像記憶手段に記憶される固有画像が、被写体を撮影して得られた医用画像の前記構造物が正常な前記各小領域の画像に基づいて主成分分析して得た画像が望ましい。
【００２６】
前記医用画像は、前記被写体の胸部を撮影して得た胸部画像である場合には、
前記分類手段が、
前記胸部画像より肺野領域を認識する肺野認識手段と、
前記胸部画像より肋骨を認識する肋骨認識手段とを有し、
肺野領域内の位置と肋骨の重なりに応じて前記多数の小領域を分類するものが望ましい。
【発明の効果】
【００２７】
本発明によれば、人体等の被写体の所定の部位を撮影した医用画像から、細分化した多数の小領域を前記構造物から得られる解剖学的特徴が互いに類似する小領域毎に分類して、分類された小領域に対応する人工画像を生成することにより、各領域の特徴に応じた人工画像を生成することができるので、異常陰影の候補領域を検出した際、候補領域とその近傍領域における人工画像と原画像とを比較することにより候補領域がアーティファクトなどによるものでなく異常陰影である否かの判定を正確に行うことが可能になる。さらに、胸部全体ではなく指定された小領域のみの人工画像を生成するため効率的である。
【００２８】
最小自乗距離または正規化相互相関値などの人工画像と原画像間の差分に応じた差分特徴量を用いて判別することにより、人工画像と原画像間の差分を定量的に判別を行うことが可能になる。
【００２９】
分類された各小領域の構造物が正常であるときに該小領域毎に現れる固有の特徴を有する複数の固有画像を予め記憶して、各小領域と前記複数の固有画像それぞれとの類似度に応じた重み付けで加算して人工画像を生成するようにすれば、特定の小領域に対応する正常構造の画像を人工的を正確に生成することが可能になる。
【００３０】
固有画像が、正常な前記各領域の画像に基づいて主成分分析して得た画像であれば、少ない画像を組み合わせて、様々な画像に対応した人工画像を生成することが可能である。
【００３１】
医用画像が胸部を撮影して得た胸部画像の場合には、肺野領域の認識と、肋骨の認識から肺野領域内の位置と肋骨の重なりに応じて小領域を分類するようにすれば、肺ガン等の異常陰影の判定を行う際の判定精度を向上させることが可能である。
【発明を実施するための最良の形態】
【００３２】
本発明の画像判定装置について、図に基づいて説明する。図１に示すように、本発明の画像判定装置１は、医用画像Ｐに撮影された部位を細分化した多数の小領域を解剖学的特徴が互いに類似する小領域に分類する分類手段１０と、医用画像より異常陰影の候補領域を検出する候補領域検出手段３０と、検出された異常陰影の候補領域を含む小領域の構造物が正常な人工画像を生成する人工画像生成手段２０と、検出された候補領域を含む小領域を注目領域として設定し、さらに、注目領域の近傍に存在する小領域を近傍領域として設定する領域設定手段４０と、注目領域および近傍領域の人工画像を生成し、注目領域の人工画像と原画像との差と、近傍領域の人工画像と原画像との差とを比較して、候補領域が異常陰影であるか否かの判定を行う判定手段５０を備える。
【００３３】
人工画像生成手段２０は、解剖学的特徴に応じて分類された各小領域の構造物が正常な構造であるときに現れる固有の特徴を有する固有画像を、分類された小領域毎に複数記憶する固有画像記憶手段２１を備え、医用画像Ｐに撮影された部位を細分化した多数の小領域のうち指定された小領域の人工画像を、その小領域が分類された小領域に対応する固有画像を重み付け加算して人工画像を生成する。
【００３４】
例えば、医用画像Pが胸部を撮影した胸部画像の場合、肺野領域の位置や肋骨の重なり（具体的には、肋骨の重なりが１本であるか２本であるか）によって解剖学的に異なった特徴が現れる。また、肺ガンのような異常陰影は、肋骨上にあるか、肋骨以外のところにあるかによって、その見え方が異なってくる。そこで、本実施の形態では、以下、胸部画像に肺ガンのような異常陰影が現れた小領域の人工画像を生成する場合について具体的に説明する。
【００３５】
分類手段１０は、肺野を認識する肺野認識手段１１と、肋骨を認識する肋骨認識手段１２とを備え、心胸郭の位置や肋骨の重なり具合などの解剖学的な特徴に応じて、解剖学的特長が類似する小領域に分類する。
【００３６】
候補領域検出手段３０は肺ガンのような異常陰影の候補領域を検出し、人工画像生成手段２０は、候補領域検出手段３０より異常陰影の候補領域が検出された場所の正常構造の人工画像を生成する。
【００３７】
さらに、判定手段５０は、人口画像と原画像間の差分に応じた差分特徴量を算出する差分特徴量算出手段５１を備え、人工画像生成手段２０により注目領域と近傍領域が正常な構造物からなる人工画像を生成し、注目領域の人工画像と原画像との差分特徴量と近傍領域の人工画像と原画像との差分特徴量から、注目領域に含まれる候補領域が異常陰影であるか否かの判定を行う。
【００３８】
以下、本実施の形態の画像判定装置１の処理の流れを、図２のフローチャートに基づいて説明を行う。
【００３９】
まず、画像判定装置１に被写体の胸部画像Ｐが入力されると（ステップ１）、分類手段１０は、肺野認識手段１１により胸部画像Ｐから心胸郭を検出して肺野領域の認識を行い（ステップ２）、さらに、肋骨認識手段１２により肺野領域より肋骨認識を行う（ステップ３）。次に、候補領域検出手段３０でガン化した部分の結節や腫瘤等の異常陰影を強調するような強調処理を行なって肺ガンのような異常陰影の候補領域を検出する（ステップ４）。そこで、領域設定手段４０は、候補領域を含むような小領域を注目領域として設定し、さらに、その注目領域に隣接する８近傍の小領域を隣接領域として設定する（ステップ５）。分類手段１０は、注目領域と隣接領域の肺野領域内の位置と、その各領域の位置での肋骨の重なり具合に応じて、注目領域と隣接領域を所定の画像タイプに分類する（ステップ６）。この分類された小領域に対応した人工画像を生成するために、人工画像生成手段２０は、固有画像記憶手段２１に記憶された複数の固有画像と候補領域との類似度を算出し、この類似度に応じて各固有画像を加算して各領域の人工画像を生成する（ステップ７）。さらに、判定手段５０では、注目領域の人工画像と医用画像の原画像との差と、近傍領域の人工画像と医用画像の原画像との差とを比較して、両差の違いの有無によって、注目領域に含まれる候補領域が異常陰影であるか否かの判定を行う（ステップ８）。
【００４０】
次に、上述の各ステップを以下詳細に説明する。
【００４１】
（１）肺野認識（ステップ２）
胸部画像に撮影された肺野領域を心胸郭を検出することにより認識する。まず、胸部画像からガボール関数等のエッジ検出マスクを用いて心胸郭のおおよその輪郭を抽出し、抽出した心胸郭の略中心となる位置を求め、その中心位置を基準に極座標変換する。極座標平面上で、基準となる平均的な心胸郭の輪郭と略相似形のテンプレートを用いてテンプレートマッチング処理を行なうことにより心胸郭の輪郭を自動検出する。検出した心胸郭の輪郭に囲まれた領域を肺野領域として認識し、得られた肺野領域を基準にして、図３に示すような肺尖部（２、７の部分）、辺縁部（３、８の部分）、縦隔部（４、９の部分）、横隔膜下部（５、１０の部分）の各領域を抽出する（例えば、本出願人が出願の特開２００３−６６６１号公報参照）。
【００４２】
あるいは、米国特許６５４９６４６号公報において提案されている方法を用いて、肺尖部（２、７の部分）、辺縁部（３、８の部分）、縦隔部（４、９の部分）、横隔膜下部（５、１０の部分）の各領域に分割して、辺縁部を抽出するようにしてもよい。
【００４３】
（２）肋骨認識（ステップ３）
肋骨認識は、例えば、胸部が撮影されたサンプル画像を教師データとして、異常陰影を含まない正常構造の構造物の統計学的なモデルを予め作成しておき、そのモデルから入力された胸部画像Pに対応する肋骨形状を人工的に生成して行う。
【００４４】
まず、多くの胸部画像から肋骨が明瞭に撮影されているサンプル画像としてＮ枚選択し、胸部画像を表示してマウス等のポインティングデバイスを用いて各画像の肋骨上の点を前肋骨と後肋骨とを合わせてランドマークをn点(例えば、400点)指定して教師データとして用いてモデルを予め作成する。
【００４５】
まず、Ｎ枚のサンプル画像の各々についてランドマークをｎ点指定した肋骨の形状X＝（ｘ1 ，ｙ1 ，・・・，ｘi ，ｙi ，・・・，ｘn ，ｙn ）と、Ｎ枚の肋骨の形状の平均から胸部肋骨の平均形状Xave＝（ｘave1 ，ｙave1 ，・・・，ｘavei ，ｙavei ，・・・，ｘaven ，ｙaven ）を求める（図４（ａ）参照。○は前肋骨を示し、△は後肋骨を示す）。Ｎ枚のサンプル画像の肋骨の形状Xとその平均形状Xaveとの差分ベクトルΔXｊ＝Xｊ−Xave（ｊ＝１、・・・、Ｎ）を求め、Ｎ個の差分ベクトルΔXｊ（ｊ＝１、・・・、Ｎ）について主成分分析を行ない、第１主成分から第ｍ主成分までの固有ベクトル（以下、主成分形状という）Pｓ（ｓ＝１、・・・、ｍ）を求める。
【００４６】
主成分分析を行うと図４(b）、(c)に示すように、第１主成分形状P１は肋骨を図４(b)の矢印の方向に拡げる成分として現れ、第２主成分形状P2は肋骨を図４(c)の矢印の方向に拡げる成分として現れる。肋骨の任意形状のモデルは平均形状Xaveと各主成分形状Ps（ｓ＝１、・・・、ｍ）の線形和により次式のように近似することができ、形状係数bｓを変えることにより、平均形状からワーピングして様々な肋骨形状を生成することができる。
【数１】
2006325937.tif 000003
【００４７】
そこで、入力された胸部画像Ｐと一致する肋骨形状を人工的に生成するために、上記形状係数bｓを求める。胸部画像Ｐから撮影されている肋骨の上の点を胸部画像Ｐからいくつか取得し、その肋骨の上の点の座標値を上式に代入して形状係数bｓの解を求める（具体的には、例えば、形状係数bｓと同じ数の点を代入することにより連立方程式の解として得られる）。肋骨の形状が明瞭に撮影されていない胸部画像Pであっても、求めた形状係数bｓを上記の式に代入することにより肋骨の全体の形状を生成することができる。具体的には、胸部画像では後肋骨などはエッジ検出により抽出することが容易であるため、後肋骨のエッジ上の点を抽出して形状係数ｂsを求めることが可能である（本出願人が出願の特開２００４‐４１６９４号公報等参照）。
【００４８】
あるいは、別の手法として、胸部画像から肋骨のエッジを抽出し、抽出された肋骨上の点をBスプラインなどで補間して肋骨形状を抽出する手法を用いてもよい。
【００４９】
具体的には、図５に示すように、胸部画像から検出された肋骨のエッジの曲線上の複数の点Ｐ１、Ｐ2、Ｐ3、・・・を抽出し、これらの点を補間するBスプライン曲線Ｐ(ｔ)を求める（図５及び図６参照）。ｎ次のBスプライン曲線Ｐ(ｔ)は、制御点Ｑｉ（ｉ＝1,2，・・・，n）とパラメータｔで定義され下式のように表される。
【数２】
2006325937.tif 000004
【００５０】
また、ｎ＝３の場合の3次のBスプライン曲線P(t)は、次のように表わすこともできる。
【数３】
2006325937.tif 000005
【００５１】
そこで、式（２）においてｔ＝０とすると、次式が導かれる。
【数４】
2006325937.tif 000006
【００５２】
制御点は図６に示すように与えられ、第２の制御点Q２と第３の制御点Q３は、肋骨のエッジを表す曲線の始点と終点の接線ｔ１、ｔ２上に存在する。そこで、この関係と肋骨のエッジの曲線上の点Pｉの位置と上記(３)式の関係を満足するように制御点Qｉ（i＝1,2,3,・・・）を求め（詳細は、”MEDICAL IMAGING TECHNOLOGY Vol.20 No.6 November 2002，Page694-Page701”参照）、抽出したエッジ上の点をBスプライン曲線で補間して肋骨の形状を得るようにしてもよい。
【００５３】
（３）肺ガンの候補領域の検出（ステップ４）
候補領域検出手段３０は、まず、胸部画像P上に表れる結節や腫瘤等の異常陰影候補を強調する強調処理を施し強調処理画像１１０を作成する。胸部画像P上に撮影された癌化した部分の結節や腫瘤等の異常陰影は、おおむね丸味をおびた輪郭を持ち、かつ、画像上では周囲に比べて画素値が大きく勾配ベクトルが集中する領域として観測される。このような異常陰影は、半球状で同じ濃度の画素が同心円状に広がる形状の円形凸領域として観測され、円形凸領域は、画素値（濃度値）の分布が周縁部から中心部に向かうにしたがって画素値が高く（濃度値が低く）なるという画素値の勾配が認められる。その勾配線は異常陰影の中心方向に向かって集中するもので、画素値の勾配を勾配ベクトルとして算出し、その勾配ベクトルの集中度から異常陰影を検出することができることができる。
【００５４】
勾配ベクトルの集中度は、具体的に、次のようにして求める。
【００５５】
まず、勾配ベクトルは計算対象となる画像を構成する全ての画素について、下記式（４）に示す計算式に基づいて勾配ベクトルの向きφを求める。
【数５】
2006325937.tif 000007
【００５６】
ここでｆ1１ 〜ｆ５５は、図７に示すように、その画素ｊを中心とした縦５画素×横５画素のマスクの外周上の画素に対応した画素値である。
【００５７】
そこで、対象となる画像を構成する全ての画素Ｐについて、勾配ベクトルの集中度Ｃを式（５）にしたがって算出する。
【数６】
2006325937.tif 000008
【００５８】
ここで、Ｎは注目画素を中心に半径ｌの円内に存在する画素の数、θj は、画素iとその円内の各画素ｊとを結ぶ直線と、その各画素ｊにおける上記式（４）で算出された勾配ベクトルとのなす角である（図８参照）。
【００５９】
上記式（２）で表される集中度Ｃが大きな値となるのは、各画素ｊの勾配ベクトルの向きが集中する画素で大きな値をとる。結節や腫瘤の陰影近傍の各画素ｊの勾配ベクトルは、その陰影のコントラストの大小に拘らず、その腫瘤陰影の略中心部に向くため、上記集中度Ｃが大きな値を採る画素は、異常陰影の中心部の画素である。
【００６０】
また、このような集中度を評価するものの中には、異常陰影の大きさや形状に左右されない検出力を達成するために、フィルタの大きさと形状に工夫をしたものがあり、その一例として、適応リングフィルタ（Adaptive  ring  filter）が挙げられる。
【００６１】
異常陰影は中心から周辺に向かって画素値が次第に大きくなるものだけではなく、肺がんに現われる結節等のように中央部分における画素値の変化が単調ではなくベクトル場が乱れて集中度が小さくなるケースもある。適応リングフィルタは、画素値が単調に変化するものだけでなく、中央部分における画素値が単調ではなくベクトル場が乱れて集中度が小さくなるケースにも適応できるものである。
【００６２】
適応リングフィルタのマスク領域は図９に示すリング状領域で、内側半径は適応的に決定される。内円と外円の半径をそれぞれｒとＲとすると、その関係は式（６）のように表される。
【００６３】
Ｒ＝ｒ＋ｄ                                  （６）
ただし、ｄは定数である。
【００６４】
また、適応リングフィルタの出力は、
【数７】
2006325937.tif 000009
【００６５】
となり、適応リングフィルタの出力は、円形凸領域の中心付近で極大値をとる。例えば、原画像上に現れる図１０（Ａ）に示すような円形凸領域は、白い線上で図１０（Ｂ）に示すような画素値を持ち、適応リングフィルタを施すと図１０（Ｃ）に示すように原画像より高い画素値が中心部分に現れる。（詳細は、例えば、魏  軍、荻原  義裕、小畑  秀文、“がん陰影候補抽出のための勾配ベクトル集中フィルタ”、電子情報通信学会論文誌（D-II）  Vol.J83-D-II  No.1，pp.118-125，Jan．2000．参照）。
【００６６】
そこで、この適応リングフィルタ処理を胸部画像Pに施して、結節や腫瘤が強調された強調処理画像１１０を作成する。図１１では、ｌ＝20mm、ｄ＝4mmとした適応リングフィルタを用いて結節を強調した例を示す（以下、結節の検出について説明する）。図１１（Ａ）の原画像上に現れている結節（白い矢印）の部分が適応リングフィルタ処理を施すと図１１（Ｂ）に示すように強調される。
【００６７】
しかし、心胸郭の辺縁部では、肋骨等が重なり合って撮影されるため濃度勾配の集中度が乱れ、原画像の胸部画像Pを用いて強調処理を行っても円形凸領域がうまく強調されない。そこで、辺縁部においては、背景画像の影響を除去して強調処理を行う方がよい。
【００６８】
具体的には、前述の胸郭の抽出方法を用いて得られた辺縁部（３、８の部分）に対して、原画像から背景画像を差し引いた差分画像１２０を作成し、差分画像１２０に対して強調処理を行うことにより、背景画像の影響を除去して結節を強調することができる。具体的には、例えば原画像をガウスフィルタでぼかした平滑化画像を原画像から差し引くことによって背景画像成分を除去することができる。
【００６９】
図１２は、適応リングフィルタを用いて強調処理画像１１０を作成したときの様子を示す。図１２（Ａ）は、原画像を用いて強調処理を施して強調処理画像１１０を作成したときの様子を示すが、背景画像の影響を受け結節がうまく強調されていない。一方、図１２（Ｂ）は、辺縁部の原画像１００から、ガウスフィルタでぼかした平滑化画像を差し引いた差分画像１２０に、適応リングフィルタを用いて強調処理を施して強調処理画像１１０を作成したときの様子を示すが、背景画像の影響を受けず結節が強調されている。
【００７０】
次に、強調処理画像１１０に閾値を低い値から高い値に徐々に変えて２値化処理を施して、複数の２値化画像１３０を作成する。２値化処理は、与えられた閾値以上の画素値を持つ画素を第１画素値（例えば、画素値を８ビットで表したときに２５５（白））に置き換え、閾値以下の画素の画素値を第２画素値（例えば、画素値を８ビットで表したときに０（黒））に置き換えて２値化画像１３０を作成する。２値化処理を行うと、２値化画像１３０上に構造物や結節や腫瘤等の高い画素値を持つ領域が第１画素値に置き換えられ、それ以外の領域が第２の画素値に置き換えられて第１画素値を持つ画素が連なった領域が島状の孤立領域として現れる。与えられた閾値が低いときには、２値化画像１３０上に現れる孤立領域には背景画像に現れる白い雲状の部分等を含んで抽出されるが、閾値が高くなるに従って、背景画像を含まない結節や腫瘤や構造物等の部分のみが孤立領域として抽出される。特に、適応リングフィルタを用いて強調された結節や腫瘤は、他の構造物と比較しても画素値が高く、高い閾値で２値化した２値化画像１３０上にも孤立領域として現われる。このことから、高い閾値で２値化した２値化画像１３０上に現れる孤立領域が異常陰影の可能性が高いものといえる。
【００７１】
図１３に、閾値を変えて２値化処理を行った一例を示す。図１３（Ａ）は、胸部画像Pの原画像に対して適応リングフィルタ処理を施して円形凸領域を強調した強調処理画像１１０である。この強調処理画像１１０は８ビットで量子化したもので、０〜２５５の階調を持つ。この強調処理画像１１０に対して１００を閾値として２値化処理を施すと、図１３（Ｂ）に示すような２値化画像１３０となり、白い孤立領域（第１の画素値に置き換えられた領域）が現れる。さらに、４刻みで閾値を変えていき３９段階の２値化画像１３０を作成し、各２値化画像１３０から孤立領域を抽出する。図１３（Ｃ）、（Ｄ）は、閾値、１７６、２５２で２値化した２値化画像１３０の例である。
【００７２】
結節は上述の孤立領域が円に近い形状で現れ面積が小さいという特徴がある。一方、背景画像が含まれて抽出された孤立領域や構造物が抽出された孤立領域は、円形度が低く、面積が大きいものが多い。そこで、２値化処理して得た複数の２値化画像１３０の各々から抽出した孤立領域のうち所定の値以下の円形度と所定値以下の大きさを持つものを結節の候補として抽出する。
【００７３】
円形度ｃは、例えば、抽出された孤立領域の面積Ｓとその周長Ｌから以下のようにして求める。
【００７４】
面積Ｓを持つ正円の半径Ｒは、
【数８】
2006325937.tif 000010
【００７５】
となり、面積Ｓを持つ正円の円周の長さＬ１は、
【数９】
2006325937.tif 000011
【００７６】
となる。そこで、円形度ｃを、正円のときは１になり正円から外れるほど０に近づくように求める。そこで、孤立領域の周長の長さＬと孤立領域の面積Ｓから、
【数１０】
2006325937.tif 000012
【００７７】
として求める。
【００７８】
あるいは、上記式（１０）で求めた円形度を求めるのではなく、以下のようにして円形度を求めてもよい。図２１に示めすように、孤立領域の重心を中心として半径Ｒの近似円Ｃｉｒを描いたときに、孤立領域の面積Ｓと、孤立領域と近似円の重なり部分の面積Ａから下式（１１）のようにして、円形度ｃを求めるようにしてもよい。
【数１１】
2006325937.tif 000013
【００７９】
具体的には、胸部画像で結節の陰影を判定する場合は、各２値化画像１３０上に現れた孤立領域のうち、式（１１）における円形度が０．７以上で、２．２６ｍｍ＜半径＜４．９４ｍｍとなる孤立領域を結節の候補領域として抽出するのが適当である。
【００８０】
また、適応リングフィルタで強調された結節は、図１０に示すように、原画像上に現れる結節より中心部分に高い画素値が現れるため、閾値を変えても同じ位置に孤立領域として現れる場合が多い。そこで、閾値を段階的に変えて２値化画像１３０を作成した場合、結節の候補が同じ位置で拾い上げられる回数が多いものを結節の候補領域として検出する。
【００８１】
以上、結節の検出について説明したが、適応リングフィルタの大きさや２値化画像による判定方法を変えることにより、結節に限らず腫瘤についても検出を行うことが可能である。
【００８２】
（４）領域設定（ステップ５）
そこで、領域設定手段４０は、図１９にしめすように、前述の候補領域検出手段３０で検出した異常陰影の候補領域を含むように小領域を設定し、この小領域を注目領域とし、その注目領域の近傍の小領域を近傍領域として設定する。異常陰影のうち結節などは通常30mm以下であるので、小領域を一辺が40mm程度の矩形とし、検出した候補領域を含むように小領域を設定して注目領域とする。その候補領域を含む注目領域M（以下、中心ROIという）に隣接する８近傍の隣接領域R1,R2,R3,R4,R5,R6,R7,R8（以下、隣接ROIという）を近傍領域として設定する。
【００８３】
（５）分類（ステップ６）
次に、解剖学的特徴に基づいて肺野領域を分類する方法について説明する。まず、分類手段１０は、ステップ２の肺野認識により検出した心胸郭の輪郭を基準に図３に示す肺尖部（２、７の部分）、辺縁部（３、８の部分）、縦隔部（４、９の部分）、横隔膜下部（５、１０の部分）のいずれの位置にあるかを区分けする。
【００８４】
まず、図１４に示すように、左肺（あるいは、右肺）の心胸郭の外側の位置のＸ座標を０とし内側の位置のＸ座標を１となるようにＸ方向の座標値を正規化し、心胸郭の上部の位置のＹ座標を０とし下部の位置のＹ座標を１となるようにＹ方向の座標値を正規化する。
【００８５】
例えば、候補領域検出手段３０により検出された肺ガンの候補領域の中心位置の座標を、胸部画像Ｐの左上を原点として表わしたときに、座標値が（ｘ、ｙ）のとき、肺ガンの候補領域の中心位置を通るY方向の直線が心胸郭と交わった２点のY座標がyup（最上部の点）とydown（最下部の点）の場合には、正規化したときのY座標値newyは、
newy＝(y −yup)/(ydown −yup)
となる。同様に、肺ガンの候補領域の中心位置を通るX方向の直線が心胸郭と交わった２点のY座標がxup（最上部の点）とxdown（最下部の点）の場合には、正規化したときのY座標値newxは、
newx＝(xout −x)/( xout −xin)
となる。ただし、心胸郭の縦隔部（４、９の部分）内では、常にｘ＝1.0とする。
【００８６】
この正規化された座標値に基づいて、図１５に示すように１１個に区分けする。  
(1)  0≦newx<0.1，0≦newy<0.25
(2)  0.1≦newx<0.8，0≦newy<0.25
(3)  0.8≦newx<1.0，0≦newy<0.25
(4)  0≦newx<0.25，0.25≦newy<0.75
(5)  0.25≦newx<0.6，0.25≦newy<0.75
(6)  0.6≦newx<1.0，0.25≦newy<0.75
(7)  0≦newx<0.25，0.75≦newy<1.0
(8)  0.25≦newx<0.6，0.75≦newy<1.0
(9)  0.6≦newx<1.0，0.75≦newy<1.0
(10)  newx=1.0
(11)  newy=1.0
【００８７】
次に、肋骨認識により検出した肋骨の位置から、肋骨があるところかないところかや、肋骨の重なりの数を求めて、これを表す肋骨特徴量よって、(1)から(9)の領域に属する肺ガンの候補領域を含む小領域をさらに３つに分類する。小領域には、一部に肋骨が１本重なっている画像を含み、一部には肋骨がない部分を含むような場合があるので、肋骨の特徴量は小領域に含まれる肋骨の重なりを平均したものを用いて各画像タイプに分類する。  
(a）  肋骨特徴量=0        ；  肋骨なし
(b)   0<肋骨特徴量≦1.0   ；  肋骨１本重なり
(c)  1.0<肋骨特徴量       ；  肋骨２本重なり
以上より、肺野領域は２９（＝９×３＋２）通り画像タイプの小領域に分類される。
【００８８】
（６）人工画像生成（ステップ７）
次に、人工画像生成手段２０で分類された各小領域の人工画像を生成する方法について説明する。まず、分類された各小領域の構造物が正常な構造であるときに現れる固有の特徴を有する複数の固有画像を、主成分分析により生成する方法について説明する。固有画像は、上記の分類手段１０に分類された各小領域の画像タイプ毎に生成する。例えば、上記の(4)に属する小領域で肋骨の特徴量が(b)（肋骨１本重なり）の場合（以下、カテゴリ4-１と記載する）を例に説明する。図１６に示すように、異常陰影が撮影されていない正常な胸部画像からカテゴリ4-１に属する80画素×80画素の矩形の小領域を抽出してサンプル画像としてＭ枚（同図(a)）選択して教師データとして用いる。図１６では、肋骨が１本重なっているところを主に抽出しているがその周辺の領域を含んで主成分分析を行っている。
【００８９】
まず、Ｍ枚のサンプル画像（80画素×80画素）に含まれる6400画素の画素値を成分とする画像ベクトルＱ＝（ｑ1,１，・・・,ｑ1,８０，ｑ２,１，・・・,ｑ２,８０，・・・,ｑｉ,ｊ,・・・,ｑ８０,１，・・・,ｑ８０,８０）ついて（図１７参照）主成分分析を行い、固有ベクトル（固有画像）を求める。図１６に、サンプル画像を主成分分析して、上位９つの主成分ベクトルを求めて固有画像（80画素×80画素）Q1,Q2,・・・,Q9を取得した例を示す。得られた固有画像Q1,Q2,・・・,Q9はそれぞれ一次独立となり、固有画像Q1,Q2,・・・,Q9の線形和により人工画像を生成することができる。
【００９０】
例えば、候補領域検出手段３０により検出された肺ガンの候補領域がカテゴリ4-bに属する場合には、肺ガンの候補領域を含んだ80画素×80画素の小領域の画像Rを胸部画像Pより切り取り、画像Rと各固有画像Q1,Q2,・・・,Q9との相関値a1, a2,・・・, a9を取得する。これより、小領域の画像Rに対応する人工画像R-modelは、
R-model＝a1×Q1＋ a2×Q2＋・・・＋a9×Q9
より生成される。
【００９１】
ここでは、サンプル画像を主成分分析して人工画像を生成する場合について説明したが、肋骨認識で説明した手法と同様に、まず、サンプル画像の平均画像を求め、平均画像と各サンプル画像との差分画像を主成分分析して、平均画像をワープして人工画像を生成するようにしてもよい。
【００９２】
図１８(a)（カテゴリ4-1）に、カテゴリ4-1にガンが存在する画像に対応した正常構造の人工画像を生成した一例を示す。カテゴリ4-1にガンが存在する原画像から人工画像を差し引くとガンのところが黒く現れているようすがわかる。
【００９３】
一方、図１８(b)（カテゴリ4-1）には、ガンが存在しないカテゴリ4-1の画像に対応した人工画像を生成した一例を示す。カテゴリ4-1にガンが存在しないため、原画像から人工画像を差し引くと画像全体が一様な濃度を示す。
【００９４】
（８）判定方法（ステップ８）
しかし、上述の手法で原画像になるべく近くなるような人工画像を作成した場合であっても、肋骨のエッジなどは人工画像上に再現するのは難しいため、図１８に示すように差分画像を作成して結節などのガンが現れている部位を強調しようとすると、差分画像にエッジ部分がアーティファクトとして現れ、それを結節などのガンとして判定する場合がある。また、撮影した被写体に依存して背景画像上に特徴のあるテクスチャが部分的に現れる場合などにも、そのテクスチャを結節などのガンとして判定する可能性がある。
【００９５】
しかし、結節のような異常陰影は広い範囲にわたって現れることはなく、小さな陰影が結節のある場所にのみ現れるが、アーティファクトや背景画像の特徴のあるテクスチャは広い範囲わたって同じような陰影が現れる。つまり、候補領域の存在する小領域のみでなくその周囲にも同じような陰影が現れるかによって、候補領域が異常陰影であるか否かを判定することができる。
【００９６】
そこで、判定手段５０は、領域設定手段４０によって設定された中心ROIと隣接ROIに現れる陰影を比較することにより、中心ROI内の候補領域が、結節などの異常陰影であるか否かを判定する。
【００９７】
まず、分類手段１０により中心ROIと隣接ROI各々を、対応する画像タイプに分類し、各ROIが分類された画像タイプに応じて、前述の人工画像生成手段２０を用いて、中心ROIと隣接ROIに対応する人工画像を生成する。
【００９８】
中心ROI内の候補陰影が結節である場合には、周囲の８近傍の隣接ROIには同じような結節が表れる可能性は少ない。一方、候補陰影がアーティファクトの場合には、周囲の８近傍のROIにも同じようなアーティファクトが現れている可能性が高い。つまり、候補陰影が結節である場合には、中心ROIの原画像とその人口画像との差分をとると、その差分画像上に図１８(a)に示すような結節の陰影が現れるが、周囲の８近傍の隣接ROIの原画像とその人口画像間で差分を取った差分画像は平坦な画像になる。一方、候補陰影がアーティファクトの場合には中心ROIにも８近傍の隣接ROIにも同じような陰影が連続して現れる。
【００９９】
従って、中心ROIの候補領域が異常陰影であれば、中心ROIの原画像とその人口画像間の差と、周囲の隣接ROIの原画像とその人口画像間の差を比較すると、中心ROIの差が隣接ROIに比べるとかなり大きくなる。
【０１００】
そこで、判定手段５０は、差分特徴量算出手段５１を用いて、生成した各ROIの人工画像と原画像との差分特徴量を算出して、中心ROIの候補領域が異常陰影であるか否かの判定を行う。差分特徴量は差分画像上の画像がアーティファクトなどによるものであるか否かを判定する指標となるものであればよく、具体的には、最小自乗距離や正規化相互相関などを用いることが出来る。
【０１０１】
例えば、中心ROIの原画像と人工画像の最小自乗距離と、隣接ROIの原画像と人工画像の最小自乗距離を比較すると、中心ROIに結節のようなガンによる異常陰影が存在する場合には、図２０に示すように、中心ROIの最小自乗距離が隣接ROIの最小自乗距離に比べて相対的に大きくなる（原画像と人工画像の正規化相互相関は小さくなる）。一方、中心ROIに結節のようなガンによる異常陰影が存在しない場合には、中心ROIの最小自乗距離は他の最小自乗距離と同程度になる。
【０１０２】
具体的な判定方法として、
・ ８近傍の隣接ROIの差分特徴量の平均と、中心ROIの差分特徴量との差の違い
・ ８近傍の隣接ROIの差分特徴量の分散と、中心ROIの差分特徴量との比較
・ ８近傍の隣接ROIの差分特徴量と中心ROIの差分特徴量との差の総和
などを用いることができる。
【０１０３】
あるいは、最小自乗距離や正規化相互相関を特徴量の１つとして用いて、他の特徴量と組み合わせてSVM(support vector machine)やマハラノビス距離やニューラルネットワークなどの判別器（具体例として、本出願人が出願の特開平9-167238号公報、特開2002-74325号公報などを参照）により、異常陰影であるか否かの判別を行うようにしてもよい。
【０１０４】
上述では、注目領域の近傍にある近傍領域を中心ROIと重ならないように設定した隣接ROIとして説明したが、近傍領域が注目領域と部分的に重なるように設定してもよい。
【０１０５】
以上、胸部画像について説明したが、乳房画像においても胸筋領域、乳腺領域、脂肪領域などに区分けして（例えば、本出願人出願の特開2003-334183公報参照）、人工画像を生成して、異常陰影の判定を行ってもよい。
【０１０６】
また、上述では、30mm以下の結節を検出する場合について説明したが、小領域の大きさを適宜変えることによって、腫瘤のような30mm以上の陰影を検出するようにしてもよい。
【０１０７】
また、上述の処理をコンピュータ上で実行するプログラムをＣＤ−ＲＯＭやネットワークを介してパソコン等のコンピュータにインストールしたものを画像判定装置とすることができる。
【０１０８】
以上詳細に説明したように、解剖学的特徴を有する構造物を含む医用画像を細分化し、各小領域を分類して分類された小領域毎にベースとなる固有画像を生成して、固有画像から人工画像を生成するため、その小領域にあった特長が現れた画像を生成することができ人工画像の精度を上げることが可能である。また、部分的な画像を生成するため、全体の画像を生成するより効率的に画像を生成することが可能になる。
【０１０９】
さらに、異常陰影の候補領域を検出した際、検出した候補領域を含む小領域とその近傍の近傍領域における人口画像と原画像の差から検出した候補領域が異常陰影であるか否かの判定を行うようにすれば、アーティファクトなどによる誤検出を少なくすることができる。
【図面の簡単な説明】
【０１１０】
【図１】本発明の画像判定装置の概略構成を示す図
【図２】本発明の画像判定装置の処理の流れを説明するためのフローチャート
【図３】肺野領域の分割結果の一例
【図４】主成分分析により肋骨認識をする手法を説明するための図
【図５】肋骨の形状をＢスプライン曲線を用いて抽出する手法を説明するための図
【図６】Ｂスプライン曲線の制御点を説明するための図
【図７】勾配を説明するための図
【図８】集中度の算出を説明する方法を説明するための図
【図９】適応リングフィルタの算出を算出する方法を説明するための図
【図１０】適応リングフィルタで出力した画素値を説明するための図
【図１１】適応リングフィルタの出力画像の一例
【図１２】差分画像を用いた適応リングフィルタの出力結果を説明するための図
【図１３】閾値を変えたときの２値化画像の変化の様子を表す図
【図１４】肺野領域の正規化する方法を説明するための図
【図１５】正規化された位置に基づく分類を説明するための図
【図１６】主成分分析した一例を表す図
【図１７】サンプル画像に含まれる画素の成分を説明するための図
【図１８】人工画像を生成した例を示す図
【図１９】注目領域と近傍の隣接領域の一例
【図２０】人工画像と原画像の最小自乗距離の一例
【図２１】円形度の求め方を説明するための図
【符号の説明】
【０１１１】
１  画像判定装置
10  分類手段
11  肺野認識手段
12  肋骨認識手段
20  人工画像生成手段
21  固有画像記憶手段
30  候補領域検出手段
40  領域設定手段
50  判定手段
51  差分特徴量算出手段
Ｐ  医用画像
【図１】
2006325937.tif 000014
【図２】
2006325937.tif 000015
【図３】
2006325937.tif 000016
【図４】
2006325937.tif 000017
【図６】
2006325937.tif 000018
【図７】
2006325937.tif 000019
【図８】
2006325937.tif 000020
【図９】
2006325937.tif 000021
【図１４】
2006325937.tif 000022
【図１５】
2006325937.tif 000023
【図１７】
2006325937.tif 000024
【図２０】
2006325937.tif 000025
【図２１】
2006325937.tif 000026
【図５】
2006325937.tif 000027
【図１０】
2006325937.tif 000028
【図１１】
2006325937.tif 000029
【図１２】
2006325937.tif 000030
【図１３】
2006325937.tif 000031
【図１６】
2006325937.tif 000032
【図１８】
2006325937.tif 000033
【図１９】
2006325937.tif 000034