(57)【特許請求の範囲】
【請求項１】
画像の特徴量群から当該画像の分類を判定する機械学習による識別器に入力する特徴量群を算出する特徴量群算出部に入力する変換画像データを、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造を個体として、世代毎に設けられる個体集団から、適応度に従って一部の個体群を選択し、選択した個体群に含まれる個体の生殖処理により世代交代させる進化的計算により、最適解あるいは実用解となる個体であるフィルタ構造を探索する進化型画像自動分類装置であって、以下の要素を有することを特徴とする進化型画像自動分類装置
（１）フィルタ構造を個体として、世代毎に個体集団を記憶する個体集団記憶部
（２）前記特徴量群算出部と前記識別器を有し、世代毎に、個体集団記憶部に記憶している個体集団に含まれる各個体であるフィルタ構造について、学習用画像データを変換対象として当該フィルタ構造に従って学習用の変換画像データを生成し、生成した学習用の変換画像データから特徴量群を算出し、学習用の変換画像データの特徴量群を入力して学習用画像データの正解の分類である学習用分類データを出力するように識別器を学習させ、評価用画像データを変換対象として前記フィルタ構造に従って評価用の変換画像データを生成し、生成した評価用の変換画像データから特徴量群を算出し、前記学習した識別器に、評価用の変換画像データの特徴量群を入力して分類結果を出力させ、分類結果が評価用画像データの正解の分類である評価用分類データと一致する正答率を算出し、算出した正答率を適応度として、各個体の適応度を算出する適応度計算部
（３）世代毎に、算出した適応度に従って個体集団から一部の個体群を選択する個体選択部
（４）世代毎に、選択した個体群に含まれる個体のフィルタ構造から次世代の個体となるフィルタ構造を生殖処理により生成し、生成したフィルタ構造群を次世代の個体集団として個体集団記憶部に記憶させる生殖部。
【請求項２】
フィルタ構造は、順次前段フィルタの出力画像データを入力して画像を変換する多段フィルタを示すフィルタ識別情報の列であることを特徴とする請求項１記載の進化型画像自動分類装置。
【請求項３】
進化型画像自動分類装置は、更に、実用解の条件となる適応度基準を、前記算出した適応度が越えた場合に、当該適応度のフィルタ構造を出力する学習個体出力部を有することを特徴とする請求項１記載の進化型画像自動分類装置。
【請求項４】
画像の特徴量群から当該画像の分類を判定する機械学習による識別器と、識別器に入力する特徴量群を算出する特徴量群算出部と、特徴量群算出部に入力する変換画像データを、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造を個体として、世代毎に個体集団を記憶する個体集団記憶部とを有し、世代毎に設けられる個体集団から、適応度に従って一部の個体群を選択し、選択した個体群に含まれる個体の生殖処理により世代交代させる進化的計算により、最適解あるいは実用解となる個体であるフィルタ構造を探索する進化型画像自動分類装置によるフィルタ構造生成方法であって、以下の要素を有することを特徴とするフィルタ構造生成方法
（１）世代毎に、個体集団記憶部に記憶している個体集団に含まれる各個体であるフィルタ構造について、学習用画像データを変換対象として当該フィルタ構造に従って学習用の変換画像データを生成し、生成した学習用の変換画像データから特徴量群を算出し、学習用の変換画像データの特徴量群を入力して学習用画像データの正解の分類である学習用分類データを出力するように識別器を学習させ、評価用画像データを変換対象として前記フィルタ構造に従って評価用の変換画像データを生成し、生成した評価用の変換画像データから特徴量群を算出し、前記学習した識別器に、評価用の変換画像データの特徴量群を入力して分類結果を出力させ、分類結果が評価用画像データの正解の分類である評価用分類データと一致する正答率を算出し、算出した正答率を適応度として、各個体の適応度を算出する適応度計算工程
（２）世代毎に、算出した適応度に従って個体集団から一部の個体群を選択する個体選択工程
（３）世代毎に、選択した個体群に含まれる個体のフィルタ構造から次世代の個体となるフィルタ構造を生殖処理により生成し、生成したフィルタ構造群を次世代の個体集団として個体集団記憶部に記憶させる生殖工程。
【請求項５】
フィルタ構造は、順次前段フィルタの出力画像データを入力して画像を変換する多段フィルタを示すフィルタ識別情報の列であることを特徴とする請求項４記載のフィルタ構造生成方法。
【請求項６】
フィルタ構造生成方法は、更に、実用解の条件となる適応度基準を、前記算出した適応度が越えた場合に、当該適応度のフィルタ構造を出力する学習個体出力工程を有することを特徴とする請求項４記載のフィルタ構造生成方法。
【請求項７】
画像の特徴量群から当該画像の分類を判定する機械学習による識別器と、識別器に入力する特徴量群を算出する特徴量群算出部と、特徴量群算出部に入力する変換画像データを、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造を個体として、世代毎に個体集団を記憶する個体集団記憶部とを有し、世代毎に設けられる個体集団から、適応度に従って一部の個体群を選択し、選択した個体群に含まれる個体の生殖処理により世代交代させる進化的計算により、最適解あるいは実用解となる個体であるフィルタ構造を探索する進化型画像自動分類装置となるコンピュータに、以下の手順を実行させるためのプログラム
（１）世代毎に、個体集団記憶部に記憶している個体集団に含まれる各個体であるフィルタ構造について、学習用画像データを変換対象として当該フィルタ構造に従って学習用の変換画像データを生成し、生成した学習用の変換画像データから特徴量群を算出し、学習用の変換画像データの特徴量群を入力して学習用画像データの正解の分類である学習用分類データを出力するように識別器を学習させ、評価用画像データを変換対象として前記フィルタ構造に従って評価用の変換画像データを生成し、生成した評価用の変換画像データから特徴量群を算出し、前記学習した識別器に、評価用の変換画像データの特徴量群を入力して分類結果を出力させ、分類結果が評価用画像データの正解の分類である評価用分類データと一致する正答率を算出し、算出した正答率を適応度として、各個体の適応度を算出する適応度計算手順
（２）世代毎に、算出した適応度に従って個体集団から一部の個体群を選択する個体選択手順
（３）世代毎に、選択した個体群に含まれる個体のフィルタ構造から次世代の個体となるフィルタ構造を生殖処理により生成し、生成したフィルタ構造群を次世代の個体集団として個体集団記憶部に記憶させる生殖手順。
【請求項８】
フィルタ構造は、順次前段フィルタの出力画像データを入力して画像を変換する多段フィルタを示すフィルタ識別情報の列であることを特徴とする請求項７記載のプログラム。
【請求項９】
プログラムは、更に、実用解の条件となる適応度基準を、前記算出した適応度が越えた場合に、当該適応度のフィルタ構造を出力する学習個体出力手順をコンピュータに実行させることを特徴とする請求項７記載のプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像から算出した特徴量群を用いた識別器（分類器ともいう。）による画像の分類装置に係り、特徴量群算出部の前段に設けたフィルタ列を進化的計算により最適化させることにより、特徴量群算出部のインプット画像として、当該識別器による分類に適した変換画像データを得ることにより、画像の分類の精度を向上させる技術に関する。
【背景技術】
【０００２】
近年工業製品製造ラインにおいて、機械によるオートメーション化が急速に進んでいる。検品ラインにおいても不良品の抽出を画像を用いてコンピュータで処理することにより自動化する取り組みがおこなわれている。しかし現状では疵を認識する際の条件を画像処理の熟練者が手動でプログラムを組むなど、多大なコストと時間が必要とされている。さま最近パターン認識の分野で高い汎化能力を示す学習方法が考案されているが、それらを適用してもなお人間の識別能力に及ばない。
【０００３】
そこで発明者の所属している研究室が提案している進化的画像処理とこれらのパターン認識の手法を組み合わせることにより、高い認識結果を得られる識別器の獲得を目指すのが本発明の目的である。進化的画像処理とは与えられた任意の画像処理を最適化問題として定式化し、その問題を進化的計算法によって解く考え方である。
【０００４】
パターン認識とは、認識対象がいくつかの概念に分類できるとき、観測されたパターンをそれらの概念のうちのひとつに対応させる処理である。たとえば数字の認識は入力パターンを０〜９の１０種類の数字のいずれかに対応させることが目的である。
【０００５】
パターンを認識する機械を実現するためには特徴量抽出と学習・識別のふたつの課程で行われる。特徴量抽出では認識対象から何らかの特徴量を計測・抽出し、それを数値化する。このような特徴量はまとめて特徴ベクトルｘ＝（ｘ1、ｘ2、…、ｘM）Tとして表される。もし特徴ベクトルの選び方が適切ならば、同じクラスの特徴ベクトルは互いに似ており、異なるクラスの特徴ベクトルは互いに違っていると考えられる。ここで得られた特徴ベクトルからその対象がどのクラスに属するかを判定する識別方法を学習する。そのためにはまずクラスの帰属が既知の学習用のサンプル集合から特徴ベクトルとクラスとの確率的な対応関係を知識として学習することが必要である。このような学習は教師あり学習と呼ばれている。次に学習された特徴ベクトルとクラスとの対応関係に関する確率的知識を利用して、与えられた未知の認識対象の特徴からその認識対象がどのクラスに属していたかを決定する方式が必要となる。その際間違って識別する確率をできるだけ小さくすることが望ましい。この未知のデータに対する識別能力を汎化能力とよぶ。パターン認識の手法としてサポートベクタマシン（Ｓｕｐｐｏｒｔ  Ｖｅｃｔｏｒ  Ｍａｃｈｉｎｅ：ＳＶＭ）を用いた。
【０００６】
ＳＶＭを用いた画像分類装置について説明する。図２５は、従来技術における識別器の学習実行処理に係る構成を示す図である。学習用画像データ２５０２と学習用分類データ２５０２のセットである学習用画像分類データセット２５０１をＬ個用意し、特徴量群算出部２５０４で学習用画像データ２５０２から複数の特徴量を算出する。そしてこれらの特徴量群と学習用分類データ２５０３を対として識別器２５０５に入力する処理をＬ個の学習用画像分類データセットについて行ない、学習実行処理を行なう。
【０００７】
図２６は、従来技術における識別器の分類実行処理に係る構成を示す図である。特徴量群算出部２５０４で判定用画像データ２６０１の特徴量群を算出し、その特徴量群を識別器２５０５に入力し、分類実行処理することにより、判定用画像データ２６０１の分類結果を得る。
【０００８】
このようにＳＶＭを用いた画像分類において、分類の正答率を向上させる方法として２通りのアプローチがある。ひとつはＳＶＭの入力として、目的とする分類に適した特徴量を選択すること、もうひとつは画像からより的確に特徴量を抽出することである。
【０００９】
既存の画像解析ソフトでは、１枚の画像から１００を超える特徴量を得ることができる。これらの特徴量すべてをＳＶＭに入力して分類することもできるが、実際には適切な入力特徴を選ぶことが汎化能力の向上につながる。逆に言えば特徴量を増やせば増やすほど汎化能力が向上するとは限らないのである。したがって目的とする分類に適した特徴量を用意することが重要になる。
【００１０】
また画像から特徴量を算出する際に、たとえば画像にノイズがのっていて対象とする物体を的確に抽出することができない場合などは特徴量の値も正確ではなくなり、分類に悪影響を及ぼすと考えられる。このようなことを解決するためには画像にあらかじめ画像フィルタを適用し、対象物体の特徴量を的確に抽出できるようにすることが考えられる。しかしながら現在画像処理の分野で用いられているフィルタの数は膨大であり、どのフィルタをどのような順序で適用すればよいかを人手で判断するのは非常に困難である。発明者の所属する研究室ではこのような問題をＧＡなどの進化計算法を用いて解決する進化的画像処理を提案している。進化的画像処理とはフィルタの構造をＧＡや遺伝的プログラミング（Ｇｅｎｅｔｉｃ  Ｐｒｏｇｒａｍｍｉｎｇ；  ＧＰ）に適用し、目的とする画像処理を自動構築するものである。
【特許文献１】特開２００１−１３４７６３号公報
【特許文献２】特開２００５−３０９５３５号公報
【発明の開示】
【発明が解決しようとする課題】
【００１１】
本発明では画像分類において正答率を向上させることを目標に、画像から的確な特徴量を抽出することに重きを置いて、２通りの方法を提案する。いずれも原画像をＳＶＭでの識別器に入力する前に画像フィルタを複数回適用する構造とし、その構造をＧＡによって学習するものである。
【課題を解決するための手段】
【００１２】
本発明に係る進化型画像自動分類装置は、
画像の特徴量群から当該画像の分類を判定する機械学習による識別器に入力する特徴量群を算出する特徴量群算出部に入力する変換画像データを、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造を個体として、世代毎に設けられる個体集団から、適応度に従って一部の個体群を選択し、選択した個体群に含まれる個体の生殖処理により世代交代させる進化的計算により、最適解あるいは実用解となる個体であるフィルタ構造を探索する進化型画像自動分類装置であって、以下の要素を有することを特徴とする。
（１）フィルタ構造を個体として、世代毎に個体集団を記憶する個体集団記憶部
（２）前記特徴量群算出部と前記識別器を有し、世代毎に、個体集団記憶部に記憶している個体集団に含まれる各個体であるフィルタ構造について、学習用画像データを変換対象として当該フィルタ構造に従って学習用の変換画像データを生成し、生成した学習用の変換画像データから特徴量群を算出し、学習用の変換画像データの特徴量群を入力して学習用画像データの正解の分類である学習用分類データを出力するように識別器を学習させ、評価用画像データを変換対象として前記フィルタ構造に従って評価用の変換画像データを生成し、生成した評価用の変換画像データから特徴量群を算出し、前記学習した識別器に、評価用の変換画像データの特徴量群を入力して分類結果を出力させ、分類結果が評価用画像データの正解の分類である評価用分類データと一致する正答率を算出し、算出した正答率を適応度として、各個体の適応度を算出する適応度計算部
（３）世代毎に、算出した適応度に従って個体集団から一部の個体群を選択する個体選択部
（４）世代毎に、選択した個体群に含まれる個体のフィルタ構造から次世代の個体となるフィルタ構造を生殖処理により生成し、生成したフィルタ構造群を次世代の個体集団として個体集団記憶部に記憶させる生殖部。
【００１３】
フィルタ構造は、順次前段フィルタの出力画像データを入力して画像を変換する多段フィルタを示すフィルタ識別情報の列であることを特徴とする。
【００１４】
進化型画像自動分類装置は、更に、実用解の条件となる適応度基準を、前記算出した適応度が越えた場合に、当該適応度のフィルタ構造を出力する学習個体出力部を有することを特徴とする。
【００１５】
本発明に係るフィルタ構造生成方法は、
画像の特徴量群から当該画像の分類を判定する機械学習による識別器と、識別器に入力する特徴量群を算出する特徴量群算出部と、特徴量群算出部に入力する変換画像データを、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造を個体として、世代毎に個体集団を記憶する個体集団記憶部とを有し、世代毎に設けられる個体集団から、適応度に従って一部の個体群を選択し、選択した個体群に含まれる個体の生殖処理により世代交代させる進化的計算により、最適解あるいは実用解となる個体であるフィルタ構造を探索する進化型画像自動分類装置によるフィルタ構造生成方法であって、以下の要素を有することを特徴とする。
（１）世代毎に、個体集団記憶部に記憶している個体集団に含まれる各個体であるフィルタ構造について、学習用画像データを変換対象として当該フィルタ構造に従って学習用の変換画像データを生成し、生成した学習用の変換画像データから特徴量群を算出し、学習用の変換画像データの特徴量群を入力して学習用画像データの正解の分類である学習用分類データを出力するように識別器を学習させ、評価用画像データを変換対象として前記フィルタ構造に従って評価用の変換画像データを生成し、生成した評価用の変換画像データから特徴量群を算出し、前記学習した識別器に、評価用の変換画像データの特徴量群を入力して分類結果を出力させ、分類結果が評価用画像データの正解の分類である評価用分類データと一致する正答率を算出し、算出した正答率を適応度として、各個体の適応度を算出する適応度計算工程
（２）世代毎に、算出した適応度に従って個体集団から一部の個体群を選択する個体選択工程
（３）世代毎に、選択した個体群に含まれる個体のフィルタ構造から次世代の個体となるフィルタ構造を生殖処理により生成し、生成したフィルタ構造群を次世代の個体集団として個体集団記憶部に記憶させる生殖工程。
【００１６】
フィルタ構造生成方法は、更に、実用解の条件となる適応度基準を、前記算出した適応度が越えた場合に、当該適応度のフィルタ構造を出力する学習個体出力工程を有することを特徴とする。
【００１７】
本発明に係るプログラムは、
画像の特徴量群から当該画像の分類を判定する機械学習による識別器と、識別器に入力する特徴量群を算出する特徴量群算出部と、特徴量群算出部に入力する変換画像データを、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造を個体として、世代毎に個体集団を記憶する個体集団記憶部とを有し、世代毎に設けられる個体集団から、適応度に従って一部の個体群を選択し、選択した個体群に含まれる個体の生殖処理により世代交代させる進化的計算により、最適解あるいは実用解となる個体であるフィルタ構造を探索する進化型画像自動分類装置となるコンピュータに、以下の手順を実行させることを特徴とする。
（１）世代毎に、個体集団記憶部に記憶している個体集団に含まれる各個体であるフィルタ構造について、学習用画像データを変換対象として当該フィルタ構造に従って学習用の変換画像データを生成し、生成した学習用の変換画像データから特徴量群を算出し、学習用の変換画像データの特徴量群を入力して学習用画像データの正解の分類である学習用分類データを出力するように識別器を学習させ、評価用画像データを変換対象として前記フィルタ構造に従って評価用の変換画像データを生成し、生成した評価用の変換画像データから特徴量群を算出し、前記学習した識別器に、評価用の変換画像データの特徴量群を入力して分類結果を出力させ、分類結果が評価用画像データの正解の分類である評価用分類データと一致する正答率を算出し、算出した正答率を適応度として、各個体の適応度を算出する適応度計算手順
（２）世代毎に、算出した適応度に従って個体集団から一部の個体群を選択する個体選択手順
（３）世代毎に、選択した個体群に含まれる個体のフィルタ構造から次世代の個体となるフィルタ構造を生殖処理により生成し、生成したフィルタ構造群を次世代の個体集団として個体集団記憶部に記憶させる生殖手順。
【００１８】
プログラムは、更に、実用解の条件となる適応度基準を、前記算出した適応度が越えた場合に、当該適応度のフィルタ構造を出力する学習個体出力手順をコンピュータに実行させることを特徴とする。
【発明の効果】
【００１９】
本発明によれば、画像から算出した特徴量群を用いた識別器による画像分類装置において、特徴量群算出部の前段に設けたフィルタ列を進化的計算により最適化させることにより、特徴量群算出部のインプット画像として、当該識別器による分類に適した変換画像データを得ることが出来る。これにより、画像分類の精度を向上させることができる。
【発明を実施するための最良の形態】
【００２０】
実施の形態１．
ＳＶＭでの識別器構成の前にフィルタ列を適用する構造として、まず図１及び図２に示すような構造を考えた。つまりＳＶＭの入力である全特徴量に対して、ある特定の１つのフィルタ列を対応させるものである。（以下この構造を［１：全］構造と表現する。）ここでフィルタ列をＧＡの学習対象とする。
【００２１】
識別器は、機械学習により画像の特徴量群から当該画像の分類を判定するように構成されている。そして、本装置は、この識別器による学習実行処理と分類実行処理の二段階の処理を通じてフィルタ列の適応度を算出する。図１は、識別器の学習実行処理に係る構成を示す図である。学習用画像データ１０２と学習用分類データ１０３のセットである学習用画像分類データセット１０１をＮ個用意し、学習用画像データ１０２をフィルタ列１０４で画像変換し、特徴量群算出部１０５でその学習用の変換画像データから異なる複数の特徴量を算出する。そしてこれらの特徴量群と学習用分類データ１０３を対として識別器２５０５に入力する処理をＮ個の学習用画像分類データセットについて行ない、学習実行処理を行なう。学習用分類データ１０３は、学習用画像データ１０２の正解の分類を示している。
【００２２】
このようにして識別器１０６で学習した結果を用いて、別に用意した評価用画像データを分類し、その結果の正答率を評価する。図２は、識別器の分類実行処理に係る構成を示す図である。評価用画像データ２０２と評価用分類データ２０３のセットである評価用画像分類データセット２０１をＭ個用意し、評価用画像データ２０２をフィルタ列１０４で画像変換し、特徴量群算出部１０５でその評価用の変換画像データから異なる複数の特徴量を算出する。そしてこれらの特徴量群を識別器２５０５に入力し、分類実行処理を行なう。そして分類結果が評価用分類データ２０３と一致するか否かを判定する。一致する場合には正答とし、不一致の場合には誤答とする。これらの処理をＭ個の評価用画像分類データセットについて行ない、正答率を算出する。そして、この正答率をＧＡの適応率として扱い、最適あるいは実用的なフィルタ列を導き出す。
【００２３】
評価用分類データ２０３は、評価用画像データ１０２の正解の分類を示している。また、フィルタ列は、変換対象の画像データを変換して生成するフィルタ群の構成を定義するフィルタ構造の例である。
【００２４】
このモデルは、たとえば画像が撮影された環境やカメラの性能が原因で全体的に暗くなってしまった画像やノイズがのってしまった画像に対し、それぞれガンマ補正フィルタなどで明るくしたり、ノイズ除去系のフィルタでノイズ除去する構造が学習されることを期待している。このようなフィルタ列が得られれば画像中の対象物がより鮮明になり、ＳＶＭでの分類が容易になって正答率の向上に結びつくのではないかと考えた。ここでフィルタ列の構造として様々なものが考えられる。たとえば木構造状フィルタやネットワーク状フィルタ、階層型フィルタなどである。これらのフィルタ構造は表現力が高く、高い性能を得ることが期待できる。従って、１本の直列構造のフィルタ列に限らず、他のフィルタ構造であっても構わない。
【００２５】
学習は染色体をビット列ではなく実数値をそのまま用いる実数値ＧＡを使用した。用意したフィルタ群のそれぞれに固有のフィルタ番号を割り振り、その番号を染色体に用いた。つまり、この例ではフィルタ列は、順次前段フィルタの出力画像データを入力して画像を変換する多段フィルタを示すフィルタ識別情報の列である。フィルタの長さは固定とし、画像に対して何も手を加えないＮＯＰフィルタ（ＮｏｎＯｐｅｒａｔｏｒ）を用意して擬似的に可変となるようにした。また用いるフィルタとして、画像を急激に変換するようなフィルタはその画像の特徴量が失われてしまう可能性が高く適切でない。したがってあまり急激に変化しないようなフィルタを用意した。
【００２６】
画像フィルタの例を示す。
番号    名称：内容
Ｆ０    ＮＯＰ：何もしない（Ｎｏｎ  Ｏｐｅｒａｔｏｒ）
Ｆ１    反転：階調値を反転させる
Ｆ２    最大値：注目画素の８近傍の階調値の最大値に置き換える
Ｆ３    最小値：注目画素の８近傍の階調値の最小値に置き換える
Ｆ４    中央値：注目画素の８近傍の階調値の中央値に置き換える
Ｆ５    平均値：注目画素の８近傍の階調値の平均値に置き換える
Ｆ６    レンジ：注目画素の８近傍の階調値の最大値と最小値の差に置き換える
Ｆ７    分散：注目画素の８近傍の分散に置き換える
Ｆ８    Ｐｒｅｗｉｔｔ：プレヴィットフィルタ
Ｆ９    Ｓｏｂｅｌ：ゾーベルフィルタ
Ｆ１０  縦エッジ：縦エッジを強調する
Ｆ１１  横エッジ：横エッジを強調する
Ｆ１２  Ｌａｐｌａｃｉａｎ：ラプラシアンフィルタ
Ｆ１３  鮮鋭化：エッジを強調して先鋭化する
Ｆ１４  原画像差分：原画像との階調値の差分
Ｆ１５  ＬｉｇｈｔＰｉｘｅｌ：判別分析法による二値化しきい値より暗い画素を０
Ｆ１６  ＤａｒｋＰｉｘｅｌ：判別分析法による二値化しきい値より明るい画素を２５５
Ｆ１７  二値化しきい値距離：注目画素の階調値と二値化しきい値との差に置き換える
Ｆ１８  ガンマ補正（＋）：ガンマ補正（γ＝２）
Ｆ１９  ガンマ補正（−）：ガンマ補正（γ＝０。５）
Ｆ２０  線形濃度変換：最小階調値と最大階調値の間で線形の濃度変換を行う
Ｆ２１  非線形濃度変換：階調値０と２５５の間で非線形の濃度変換を行う
特徴量の例を示す。
番号    名称：内容
Ｐ０    平均階調値：画像中に含まれる階調値の平均値
Ｐ１    標準偏差：全画素の階調値の標準偏差
Ｐ２    最大階調値：画像中に含まれる階調値の最大値
Ｐ３    最小階調値：画像中に含まれる階調値の最小値
Ｐ４    レンジ：全画素の階調値の最大値から最小値を引いたもの
Ｐ５    １σ内確率：平均階調値±標準偏差に収まっている画素数の割合
Ｐ６    ３σ外確率：平均階調値±３×標準偏差に収まっていない画素数の割合
Ｐ７    第１四分位数：階調値を昇順に並べたとき下から１／４の位置にあたる階調値
Ｐ８    中央値：階調値を昇順に並べたとき中央の位置にあたる階調値
Ｐ９    第３四分位数：階調値を昇順に並べたとき下から３／４の位置にあたる階調値
Ｐ１０  最頻値：もっとも出現頻度の高い階調値
Ｐ１１  歪度：階調値分布の非対称性をあらわす
Ｐ１２  尖度：階調値分布のとがり具合をあらわす
Ｐ１３  判別分析法しきい値：判別分析法による二値化の際のしきい値
Ｐ１４  低位クラス内分散：二値化時の階調値の低い方のグループ内の分散
Ｐ１５  高位クラス内分散：二値化時の階調値の高い方のグループ内の分散
Ｐ１６  クラス間分散：二値化時の高低のグループ間の分散
Ｐ１７  領域数：二値化後の対象領域数
Ｐ１８  平均周長：二値化後の対象領域の周長の平均値
Ｐ１９  最大周長：二値化後の対象領域の周長の最大値
Ｐ２０  平均フェレ径縦横比：二値化後の対象領域の外接矩形の縦横比の平均値
Ｐ２１  最大フェレ径縦横比：二値化後の対象領域の外接矩形の縦横比の最大値
Ｐ２２  合計領域面積：二値化後の対象領域の面積の合計値
Ｐ２３  平均領域面積：二値化後の対象領域の面積の平均値
Ｐ２４  最大領域面積：二値化後の対象領域の面積の最大値
Ｐ２５  平均円形度：二値化後の対象領域の円形度（４πＳ／Ｌ２）の平均
Ｐ２６  水平ライン差分率：二値化後の水平ラインスキャン時に上下ラインでの階調値相違画素数の割合
Ｐ２７  垂直ライン差分率：二値化後の垂直ラインスキャン時に左右ラインでの階調値相違画素数の割合
Ｐ２８  最大角度：二値化後の対象領域の傾き角度の最大値
Ｐ２９  最小角度：二値化後の対象領域の傾き角度の最小値
Ｐ３０  平均角度：二値化後の対象領域の傾き角度の平均値
続いて、本発明に係る進化型画像自動分類装置について詳述する。図３は、進化型画像自動分類装置の構成を示す図である。進化型画像自動分類装置は、セットデータ入力部３０１、学習用画像分類データセット評価用画像分類データセット記憶部３０２、初期個体入力部３０３、個体集団記憶部３０４、適応度計算部３０５、適応度記憶部３０６、個体選択部３０７、選択個体記憶部３０８、生殖部３０９、及び学習個体出力部３１０を有している。
【００２７】
図４は、進化計算アルゴリズムによる学習処理のフローを示す図である。まず、セットデータ入力部３０１による学習用画像分類データセット、評価用画像分類データセット入力処理（Ｓ４０１）を行ない、入力した学習用画像分類データセットと評価用画像分類データセットを学習用画像分類データセット評価用画像分類データセット記憶部３０２に記憶させる。
【００２８】
次に、初期個体入力部３０３による初期個体入力処理（Ｓ４０２）で、初期個体となるフィルタ列を入力し、個体集団記憶部３０４に記憶させる。フィルタ列は無作為に生成されたものでよい。
【００２９】
世代交代終了判定処理（Ｓ４０６）で、終了と判定するまで世代交代の処理を繰り返す。具体的には、適応度計算部３０５による適応度計算処理（Ｓ４０３）、個体選択部３０７による個体選択処理（Ｓ４０４）、生殖部３０９による生殖処理（Ｓ４０５）を繰り返す。当該世代の各個体（フィルタ列）の適応度を算出する適応度計算処理（Ｓ４０３）については、図５〜図１５を用いて後述する。個体選択部３０７による個体選択処理（Ｓ４０４）では、適応度に従って一部の個体群を選択する。例えば、適応度の高い個体から順に所定数の個体を選択する。生殖処理（Ｓ４０５）では、選択された所定数の個体を元に次世代の個体（フィルタ列）を生成する。この処理については、図１６を用いて後述する。
【００３０】
世代交代終了判定処理（Ｓ４０６）で終了と判定すると、最終世代の個体（フィルタ列）について適応度計算部３０５で適応度を計算し（Ｓ４０７）、学習個体出力部３１０による学習個体出力処理（Ｓ４０８）でそのフィルタ列と適応度を出力する。世代交代終了判定部（図示せず）では、最適解あるいは実用解を得られたと判断した場合に、終了とする。例えば、適応度が所定世代数を交代しても一定になった場合に最適解と判断して、終了する。あるいは、予め実用解の条件となる適応度基準を設定しておき、その適応度基準を前記算出した適応度が越えた場合に、終了する。
【００３１】
前述の適応度計算部３０５による適応度計算処理（Ｓ４０３）について詳述する。図５は、適応度計算部の構成を示す図である。適応度計算部３０５は、分類学習部５０１、評価分類部５０２、及び正答率算出部５０３を有している。
【００３２】
図６は、適応度計算処理フローを示す図である。個体（フィルタ列）毎に以下の処理を繰り返す（Ｓ６０１）。まず、分類学習部５０１による分類学習処理（Ｓ６０２）を行なう。この処理では、当該フィルタ列を識別器の前段に設けて識別器の学習を行なう。
【００３３】
図７は、分類学習部の構成を示す図である。分類学習部５０１は、フィルタ列画像変換部７０１、特徴量群算出部７０２、識別器学習実行部７０３、及び識別器学習データ記憶部７０４を有している。識別器学習実行部７０３と識別器学習データ記憶部７０４は、識別器１０６の一部を構成している。つまり、分類学習部５０１は、識別器１０６の学習実行に係る構成を利用している。
【００３４】
図８は、分類学習処理フローを示す図である。学習用画像データ毎に以下の処理を繰り返す（Ｓ８０１）。まず、フィルタ列画像変換部７０１によるフィルタ列画像変換処理（Ｓ８０２）を行なう。
【００３５】
図９は、フィルタ列画像変換部の構成を示す図である。フィルタ列画像変換部７０１は、フィルタ選択部９０１、フィルタ実行部９０２〜９０４、及び画像データバッファ９０６，９０７を有している。
【００３６】
図１０は、フィルタ列画像変換処理フローを示す図である。まず、学習用画像データを画像データバッファＡに書き込み（Ｓ１００１）、入力バッファポインタに画像データバッファＡのポインタを設定し、出力バッファポインタに画像データバッファＢのポインタを設定する（Ｓ１００２）。入力バッファポインタと出力バッファポインタは、フィルタ選択部９０１内に設けられている。そして、フィルタの列に従って、順次フィルタ番号毎に以下の処理を繰り返す（Ｓ１００３）。フィルタ番号により特定されるフィルタ実行部に入力バッファポインタと出力バッファポインタをセットし、当該フィルタ実行部を起動する（Ｓ１００４）。フィルタ実行部は、入力バッファポインタにより特定される画像データバッファから画像データを読込み、画像を変換し、変換した画像を出力バッファポインタにより特定される画像データバッファに書込む。フィルタ実行部の処理が終了すると、入力バッファポインタと出力バッファポインタを切り替える（Ｓ１００５）。この処理を、すべてのフィルタ番号について行なった時点で（Ｓ１００６）、出力バッファポインタで特定される画像データバッファから変換画像データを出力する（Ｓ１００７）。
【００３７】
次に、特徴量群算出部７０２による特徴量群算出処理（図８のＳ８０３）を行なう。図１１は、特徴量群算出部の構成を示す図である。複数の特徴量算出部１１０１〜１１０３を有している。これらの特徴量算出部は、識別器のインプットの仕様として決定された所定の特徴量に対応する。
【００３８】
図１２は、特徴量群算出処理フローを示す図である。それぞれの特徴量算出部に、変換画像データを入力して特徴量を算出する。この例では、第一特徴量算出部１１０１による第一特徴量算出処理（Ｓ１２０１）、第二特徴量算出部１１０２による第二特徴量算出処理（Ｓ１２０２）、第三特徴量算出部１１０３による第三特徴量算出処理（Ｓ１２０３）を行なっている。
【００３９】
そして、すべての学習用画像データについて特徴量群を算出した時点（図８のＳ８０４）で、識別器学習実行部７０３による識別器学習実行処理（Ｓ８０５）を行なう。特徴量群と学習用分類データの対をＮ個入力し、識別器としての学習を行なう。学習の結果は、識別器学習データ（ＳＶＭでは、分類値ということもある。）として識別器学習データ記憶部７０４に格納される。このデータは、後述する識別器の分類実行処理で用いられる。
【００４０】
次に、評価分類部５０２による評価分類処理（図６のＳ６０３）を行なう。この処理では、学習した識別器の前段に当該フィルタ列を設けて、識別器による評価用画像の分類を行なう。
【００４１】
図１３は、評価分類部の構成を示す図である。評価分類部５０２は、フィルタ列画像変換部７０１、特徴量群算出部７０２、識別器分類実行部１３０１、及び識別器学習データ記憶部７０４を有している。識別器分類実行部１３０１と識別器学習データ記憶部７０４は、識別器１０６の一部を構成している。つまり、評価分類部５０２は、識別器１０６の分類実行に係る構成を利用している。
【００４２】
図１４は、評価分類処理フローを示す図である。評価用画像データ毎に以下の処理を繰り返し（Ｓ１４０１〜Ｓ１４０４）、それぞれに分類結果を得る。分類学習処理と同様にフィルタ列画像変換部７０１によるフィルタ列画像変換処理（Ｓ１４０２）を行ない、次に特徴量群算出部７０２による特徴量群算出処理（Ｓ１４０３）を行なう。そして、識別器分類実行部１３０１による識別器分類実行処理（Ｓ１４０４）を行なう。この処理では、評価用の変換画像データの特徴量群を入力し、識別器学習データを用いて分類結果を算出する。そして、すべての評価用画像データについて処理した時点で終了する（Ｓ１４０５）。
【００４３】
次に、正答率算出部５０３による正答率算出処理（図６のＳ６０４）を行なう。この処理では、評価用画像に対する分類結果の正答率を算出する。
【００４４】
図１５は、正答率算出処理フローを示す図である。分類結果毎に以下の処理を繰り返す（Ｓ１５０１）。分類結果は、元となった評価用画像データに対応する評価用分類データと一致するか判定し（Ｓ１５０２）、一致した場合には、正答数（初期値＝０）をインクリメントしてカウントアップする（Ｓ１５０３）。そして、すべての分類結果について処理した時点で（Ｓ１５０４）、正答数を分類結果数（Ｍ）で割って、商を正答率とする（Ｓ１５０５）。
【００４５】
図６の適応度算出処理は、すべての個体について上述の処理した時点で終了する。（図６のＳ６０５）
前述の生殖部３０９による生殖処理（図４のＳ４０５）について説明する。図１６は、生殖処理フローを示す図である。交叉処理により得る個体数と、突然変異処理により得る個体数と、延命処理により得る個体数を決定する（Ｓ１６０１）。交叉処理（Ｓ１６０２）では、その個体数分の新たな個体を得る交叉を行なう。交叉は、選択個体記憶部３０８から複数の個体を選択し、その個体であるフィルタ列の交叉点を設定し、複数のフィルタ列間で交叉点以降の列や交叉点間のフィルタ番号列（染色体列）を相互に交換し、新たなフィルタ列を生成する。突然変異処理（Ｓ１６０３）では、その個体数分の新たな個体を得る突然変異を行なう。突然変異は、選択個体記憶部３０８から個体を選択し、変異させる単数あるいは部分列のフィルタ番号（染色体）を選択し、他のフィルタに置き換える処理を行なう。他のフィルタは、無作為に決定してよい。延命処理（Ｓ１６０４）では、その個体数分の個体を選択個体記憶部３０８から選択し、そのまま次世代の個体とする。このようにして得られた新たな個体群を、個体集団更新処理（Ｓ１６０５）により個体集団記憶部３０４に書き込む。つまり、現世代の個体群を次世代の個体群に更新する。
【００４６】
実施の形態２．
前述の［１：全］構造では、ＳＶＭの入力にあたる特徴量の意味を無視して１つのフィルタ列を学習する構造であった。これに対し、それぞれの特徴量に、より分類しやすくなるような値を入力できるように図１７から図１９のような構造を考えた。つまりＳＶＭの入力である特徴量それぞれに対して、１つのフィルタ列を対応させるものである。（以下この構造を［１：１］構造と表現する。）［１：全］構造と同様に、フィルタ列をＧＡの学習対象とする。
【００４７】
この構造は前述の改善に加え、ＳＶＭに入力するある特徴量に対し、その特徴量が画像の種類ごとに明確にされるようなフィルタ列が得られることを期待したものである。たとえば希望する画像の分類がエッジが多い画像とそうでない画像であったときに、エッジ検出のフィルタと領域面積の特徴量を対応させることによって分類が容易になると考えられる。ここで、全特徴量にエッジ検出フィルタを適用することはよい結果を得られるとは考えにくく、特徴量ごとにフィルタを用意することは有用であると考えた。その一方で特徴量ごとにフィルタ列を用意することから、ＧＡでの最適化において探索空間が膨大になり、最適解が得られにくい、あるいは時間がかかることが考えられる。
【００４８】
フィルタ列の学習は、それぞれのフィルタ列部分を１つにつなげて１フィルタ列とみなし、［１：全］構造と同様の手法で学習した。
【００４９】
図１７は、実施の形態２に係るフィルタ列の構成を示す図である。図に示すように、１つのフィルタ列を、複数のフィルタ列部分１７０１〜１７０３の連なりとして扱う。
【００５０】
図１８は、実施の形態２に係る識別器の学習実行処理に係る構成を示す図である。第一特徴量用フィルタ列部分１７０１により変換した第一特徴量用変換画像データを、特徴量群算出部１０５による第一特徴量算出に用い、同様に第二特徴量用フィルタ列部分１７０２により変換した第二特徴量用変換画像データを、特徴量群算出部１０５による第二特徴量算出に用い、第三特徴量用フィルタ列部分１７０３により変換した第三特徴量用変換画像データを、特徴量群算出部１０５による第三特徴量算出に用いる。
【００５１】
図１９は、実施の形態２に係る識別器の分類実行処理に係る構成を示す図である。分類実行処理でも、学習実行処理と同様にそれぞれの特徴量別変換画像データ（第一特徴量用変換画像データ、第二特徴量用変換画像データ、第三特徴量用変換画像データのこと）をそれぞれの特徴量算出に用いる。
【００５２】
図２０は、実施の形態２に係る分類学習部の構成を示す図である。フィルタ列画像変換部７０１への入力が、複数の特徴量別フィルタ列部分（第一特徴量用フィルタ列部分、第二特徴量用フィルタ列部分、第三特徴量用フィルタ列部分のこと）となる。
【００５３】
図２１は、実施の形態２に係るフィルタ列画像変換部の構成を示す図である。図２２は、実施の形態２に係るフィルタ列画像変換処理フローを示す図である。特徴量別フィルタ列部分毎に以下の処理を繰り返す（Ｓ２２０１）。Ｓ２２０４では、当該特徴量用フィルタ列部分から、順次フィルタ番号毎に以下の処理を繰り返す。特徴量別フィルタ列部分に対応するループ処理が終了する毎に（Ｓ２２０１〜Ｓ２２０７）、出力バッファポインタで特定される画像データバッファから特徴量別変換画像データを出力する（Ｓ２２０８）。そして、すべての特徴量別フィルタ列部分について処理した時点で終了する（Ｓ２２０９）。
【００５４】
図２３は、実施の形態２に係る特徴量群算出部の構成を示す図である。各特徴量算出部は、専用に変換された特徴量別変換画像データを用いてそれぞれの特徴量を算出する。
図２４は、実施の形態２に係る評価分類部の構成を示す図である。分類学習部の構成と同様に、フィルタ列画像変換部７０１への入力が、複数の特徴量別フィルタ列部分となる。
【００５５】
生殖処理では、前述の形態と同様にフィルタ列全体として交叉処理、突然変異処理、及び延命処理を行なう方法の他、特徴量別フィルタ列部分毎に交叉処理、突然変異処理、及び延命処理を行なう方法も有効である。
【００５６】
前述の例では、識別器の例としては、ＳＶＭに限らず他のクラスタリングを用いてもよい。例えば、ニューラルネットワークなどの他の機械学習による識別器を用いても良い。その場合には、特徴量群を入力し、分類データを出力するように、学習用神経回路網を誤差逆伝播法で学習させ、その結果得られた結合荷重が識別器学習データとなる。そして、学習用神経回路網と同様の構成の分類用神経回路網に、特徴量群を入力し、識別器学習データ記憶部の結合荷重に従って神経回路網演算処理を行い、処理結果として分類結果を出力する。
【００５７】
進化型画像自動分類装置は、コンピュータであり、各要素はプログラムにより処理を実行することができる。また、プログラムを記憶媒体に記憶させ、記憶媒体からコンピュータに読み取られるようにすることができる。コンピュータは、バスに、演算装置、データ記憶装置、メモリが接続されている。データ記憶装置は、例えばＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）やハードディスクである。メモリは、通常ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）である。プログラムは、通常データ記憶装置に記憶されており、メモリにロードされた状態で、順次演算装置に読み込まれ処理を行う。
【図面の簡単な説明】
【００５８】
【図１】識別器の学習実行処理に係る構成を示す図である。
【図２】識別器の分類実行処理に係る構成を示す図である。
【図３】進化型画像自動分類装置の構成を示す図である。
【図４】進化計算アルゴリズムによる学習処理のフローを示す図である。
【図５】適応度計算部の構成を示す図である。
【図６】適応度計算処理フローを示す図である。
【図７】分類学習部の構成を示す図である。
【図８】分類学習処理フローを示す図である。
【図９】フィルタ列画像変換部の構成を示す図である。
【図１０】フィルタ列画像変換処理フローを示す図である。
【図１１】特徴量群算出部の構成を示す図である。
【図１２】特徴量群算出処理フローを示す図である。
【図１３】評価分類部の構成を示す図である。
【図１４】評価分類処理フローを示す図である。
【図１５】正答率算出処理フローを示す図である。
【図１６】生殖処理フローを示す図である。
【図１７】実施の形態２に係るフィルタ列の構成を示す図である。
【図１８】実施の形態２に係る識別器の学習実行処理に係る構成を示す図である。
【図１９】実施の形態２に係る識別器の分類実行処理に係る構成を示す図である。
【図２０】実施の形態２に係る分類学習部の構成を示す図である。
【図２１】実施の形態２に係るフィルタ列画像変換部の構成を示す図である。
【図２２】実施の形態２に係るフィルタ列画像変換処理フローを示す図である。
【図２３】実施の形態２に係る特徴量群算出部の構成を示す図である。
【図２４】実施の形態２に係る評価分類部の構成を示す図である。
【図２５】従来技術における識別器の学習実行処理に係る構成を示す図である。
【図２６】従来技術における識別器の分類実行処理に係る構成を示す図である。
【符号の説明】
【００５９】
１０１  学習用画像分類データセット、１０２  学習用画像データ、１０３  学習用分類データ、１０４  フィルタ列、１０５  特徴量群算出部、１０６  識別器２０１、  評価用画像分類データセット、２０２  評価用画像データ、２０３  評価用分類データ、２０４  正答率算出部、３０１  セットデータ入力部３０２、  学習用画像分類データセット評価用画像分類データセット記憶部、３０３  初期個体入力部、３０４  個体集団記憶部、３０５  適応度計算部、３０６  適応度記憶部、３０７  個体選択部、３０８  選択個体記憶部、３０９  生殖部、３１０  学習個体出力部、５０１  分類学習部、５０２  評価分類部、５０３  正答率算出部、７０１  フィルタ列画像変換部、７０２  特徴量群算出部、７０３  識別器学習実行部、７０４  識別器学習データ記憶部、９０１  フィルタ選択部、９０２〜９０４  フィルタ実行部、９０５  画像データバッファＡ、９０６  画像データバッファＢ、１１０１  第一特徴量算出部、１１０２  第二特徴量算出部、１１０３  第三特徴量算出部、１３０１  識別器分類実行部、１７０１  第一特徴量用フィルタ列部分、１７０２  第二特徴量用フィルタ列部分、１７０３  第三特徴量用フィルタ列部分、２５０１  学習用画像分類データセット、２５０２  学習用画像データ、２５０３  学習用分類データ、２５０４  特徴量群算出部、２５０５  識別器、２６０１  判定用画像データ。
【図１】
0004660765.tif 000002
【図２】
0004660765.tif 000003
【図３】
0004660765.tif 000004
【図４】
0004660765.tif 000005
【図５】
0004660765.tif 000006
【図６】
0004660765.tif 000007
【図７】
0004660765.tif 000008
【図８】
0004660765.tif 000009
【図９】
0004660765.tif 000010
【図１０】
0004660765.tif 000011
【図１１】
0004660765.tif 000012
【図１２】
0004660765.tif 000013
【図１３】
0004660765.tif 000014
【図１４】
0004660765.tif 000015
【図１５】
0004660765.tif 000016
【図１６】
0004660765.tif 000017
【図１７】
0004660765.tif 000018
【図１８】
0004660765.tif 000019
【図１９】
0004660765.tif 000020
【図２０】
0004660765.tif 000021
【図２１】
0004660765.tif 000022
【図２２】
0004660765.tif 000023
【図２３】
0004660765.tif 000024
【図２４】
0004660765.tif 000025
【図２５】
0004660765.tif 000026
【図２６】
0004660765.tif 000027