(57)【要約】      （修正有）【課題】異常検出システムにおいて、学習モデルを作成する入力データの質を向上させることにより、異常検出の精度を向上させることを可能とする。【解決手段】監視対象に関するデータを用いて特徴量を計算して得られた特徴量データに対して状態空間モデルへの写像手法を用いて特徴量データに対応する特徴空間を定義する基底を生成し、基底にもとづき特徴量データを写像して、特徴空間の新たな次元数を決定し、写像された特徴量データの次元数を削減して新たな次元数とし、変換後特徴量データを生成する写像部と、混合分布手法にもとづき変換後特徴量データを複数のクラスタに分類して、学習モデル候補を作成する学習モデル作成部と、学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除することにより、学習モデル候補を改変して学習モデルを作成するモデル改変部を備える。【選択図】図１
2020144626.tif 000002
【特許請求の範囲】
【請求項１】
監視対象の異常を検出する異常検出システムであって、
監視対象に関するデータを用いて特徴量を計算して得られた特徴量データに対して状態空間モデルへの写像手法を用いて前記特徴量データに対応する特徴空間を定義する基底を生成し、前記基底にもとづき前記特徴量データを写像して、前記特徴空間の新たな次元数を決定し、写像された前記特徴量データの次元数を削減して前記新たな次元数とし、変換後特徴量データを生成する写像部と、
混合分布手法にもとづき前記変換後特徴量データを複数のクラスタに分類して、学習モデル候補を作成する学習モデル作成部と、
前記学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除することにより、前記学習モデル候補を改変して学習モデルを作成するモデル改変部と、を備える
ことを特徴とする異常検出システム。
【請求項２】
前記写像部が、
監視対象に関する他のデータを用いて特徴量を計算して得られた他の特徴量データに対して状態空間モデルへの写像手法を用いて前記基底にもとづき前記他の特徴量データを写像し、写像された前記他の特徴量データの次元数を削減して前記新たな次元数として、他の変換後特徴量データを生成し、
前記モデル改変部が、
前記学習モデル候補における各クラスタについて前記他の変換後特徴量データが生起する生起確率を計算し、前記生起確率が設定した閾値よりも低い場合に正常と判定し、
前記他の特徴量データに対応する前記他のデータが本来は異常である場合、当該クラスタを前記学習モデル候補から削除する
ことを特徴とする請求項１記載の異常検出システム。
【請求項３】
前記学習モデル作成部が、
設定された複数の混合数にもとづいて、前記変換後特徴量データを前記混合数分のクラスタに分類して、前記混合数ごとに学習モデル候補を作成し、
前記モデル改変部が、
前記混合数ごとに、前記学習モデル候補における各クラスタについて前記他の変換後特徴量データが生起する生起確率を計算し、前記生起確率が設定した閾値よりも高い場合に異常と判定し、
真陽性の発生率が最も高い前記学習モデル候補を前記学習モデルに決定する
ことを特徴とする請求項２記載の異常検出システム。
【請求項４】
前記写像部が、設定された閾値にもとづき前記特徴空間の新たな次元数を決定することを特徴とする請求項１〜３のいずれかに記載の異常検出システム。
【請求項５】
前記写像部が、
監視対象に関する新たなデータを用いて特徴量を計算して得られた新たな特徴量データに対して状態空間モデルへの写像手法を用いて前記基底にもとづき前記新たな特徴量データを写像し、写像された前記新たな特徴量データの次元数を削減して前記新たな次元数として、新たな変換後特徴量データを生成し、
当該異常検出システムにおいて、
前記学習モデルにおける各クラスタについて前記新たな変換後特徴量データが生起する生起確率を計算し、前記生起確率が設定した閾値よりも高い場合に異常と判定する異常判定部をさらに備えた
ことを特徴とする請求項１〜４のいずれかに記載の異常検出システム。
【請求項６】
前記監視対象に関するデータが、監視対象について測定されたデータであることを特徴とする請求項１〜５のいずれかに記載の異常検出システム。
【請求項７】
前記状態空間モデルへの写像手法が、主成分分析、独立成分分析、カーネル主成分分析、ＬＤＡ（Latent Dirichlet Allocation）、近傍探索手法の少なくともいずれかであることを特徴とする請求項１〜６のいずれかに記載の異常検出システム。
【請求項８】
前記混合分布手法が、正規分布、t分布、ポアソン分布、ベルヌ−イ分布、ディリクレ分布の少なくともいずれかによる分布モデルであることを特徴とする請求項１〜７のいずれかに記載の異常検出システム。
【請求項９】
前記監視対象について測定されたデータが、画像、寸法、色、欠陥、排出数、透過度、重量、強度、及び表面粗さの少なくともいずれかであることを特徴とする請求項６〜８のいずれかに記載の異常検出システム。
【請求項１０】
前記監視対象について測定されたデータが、音、温度、湿度、振動、変位、速度、加速度、回転角、角速度、角加速度、電圧、電流、光量、及び圧力の少なくともいずれかであることを特徴とする請求項６〜９のいずれかに記載の異常検出システム。
【請求項１１】
監視対象の異常を検出する異常検出プログラムであって、
コンピュータを、
監視対象に関するデータを用いて特徴量を計算して得られた特徴量データに対して状態空間モデルへの写像手法を用いて前記特徴量データに対応する特徴空間を定義する基底を生成し、前記基底にもとづき前記特徴量データを写像して、前記特徴空間の新たな次元数を決定し、写像された前記特徴量データの次元数を削減して前記新たな次元数とし、変換後特徴量データを生成する写像部、
混合分布手法にもとづき前記変換後特徴量データを複数のクラスタに分類して、学習モデル候補を作成する学習モデル作成部、及び、
前記学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除することにより、前記学習モデル候補を改変して学習モデルを作成するモデル改変部として機能させる
ことを特徴とする異常検出プログラム。
【請求項１２】
前記写像部に、
監視対象に関する他のデータを用いて特徴量を計算して得られた他の特徴量データに対して状態空間モデルへの写像手法を用いて前記基底にもとづき前記他の特徴量データを写像させ、写像された前記他の特徴量データの次元数を削減させて前記新たな次元数とさせ、他の変換後特徴量データを生成させ、
前記モデル改変部に、
前記学習モデル候補における各クラスタについて前記他の変換後特徴量データが生起する生起確率を計算させ、前記生起確率が設定した閾値よりも低い場合に正常と判定させ、
前記他の特徴量データに対応する前記他のデータが本来は異常である場合、当該クラスタを前記学習モデル候補から削除させる
ことを特徴とする請求項１１記載の異常検出プログラム。
【請求項１３】
前記学習モデル作成部に、
設定された複数の混合数にもとづいて、前記変換後特徴量データを前記混合数のクラスタに分類させて、前記混合数ごとに学習モデル候補を作成させ、
前記モデル改変部に、
前記混合数ごとに、前記学習モデル候補における各クラスタについて前記他の変換後特徴量データが生起する生起確率を計算させ、前記生起確率が設定した閾値よりも高い場合に異常と判定させ、
真陽性の発生率が最も高い前記学習モデル候補を前記学習モデルに決定させる
ことを特徴とする請求項１２記載の異常検出プログラム。
【請求項１４】
前記写像部に、
監視対象に関する新たなデータを用いて特徴量を計算して得られた新たな特徴量データに対して状態空間モデルへの写像手法を用いて前記基底にもとづき前記新たな特徴量データを写像させ、写像された前記新たな特徴量データの次元数を削減させて前記新たな次元数とさせ、新たな変換後特徴量データを生成させ、
前記コンピュータを、
前記学習モデルにおける各クラスタについて前記新たな変換後特徴量データが生起する生起確率を計算し、前記生起確率が設定した閾値よりも高い場合に異常と判定する異常判定部としてさらに機能させる
ことを特徴とする請求項１１〜１３のいずれかに記載の異常検出プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、監視対象の異常を検出する異常検出技術に関し、特に異常検出の精度を向上させ得る異常検出システムに関する。
【背景技術】
【０００２】
監視対象の異常を検出する異常検出技術においては、監視対象が異常である状態を的確に識別することが難しいという問題がある。
例えば、様々な監視対象に対して汎用的に用いることのできる異常検出システムの構築を考えた場合、監視対象のどのような状態が異常であるのかを定義することは容易ではない。
【０００３】
そこで、監視対象の異常状態をあらかじめ定義することなく、監視対象の正常な状態から得られる各種の測定データから抽出される特徴量をモデル化することで、正常な状態との比較から異常を検出する手法を用いることが考えられる。
この手法による異常検出は、何が異常であるかを定義しないため、想定外の異常状態にも対応することが可能となる。
【０００４】
一方、このような異常検出では、入力するデータの質が検出精度に与える影響は、非常に大きい。
すなわち、監視対象を測定して得られる入力データとしては、例えば画像や波形データなど様々なものがあるが、画像や波形データからは多くの特徴量を抽出することが可能である。また、検査対象が複雑になればなるほど、異常状態の判定に必要となる特徴量の選択が難しくなる。また、判定に用いる特徴量の数を単に多くするのみでは、反って判定結果の質を劣化させる特徴量が含まれてしまうという問題もある。さらに、入力データには、正常又は異常な状態を示すデータの双方において、ノイズが含まれる場合がある。
このため、このような異常検出システムにおいては、入力データの質を向上させることにより、高精度な異常検出を可能にすることが望ましい。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２０１６−２０６７８４号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
このような異常検出システムにおいて、写像手法などを用いることにより、入力データの特徴量を加工して、特徴量のうち有用なものだけを使用して異常検出を行うことが考えられる。しかしながら、このような写像手法を用いる場合でも、入力データにノイズなどが残っていると、それにより検出精度が低下するが、写像手法では、入力データのノイズを適切に除去することは、困難であった。
【０００７】
そこで、本発明者らは鋭意研究し、入力データの特徴量に対して写像手法とクラスタリング手法を適用し、一定のクラスタを除去することによって、特徴量の外れ値を削除し、異常検出の精度を向上させることに成功して、本発明を完成させた。
ここで、写像手法とクラスタリング手法を用いる技術として、特許文献１に記載のデータ分析装置を挙げることができる。しかしながら、この装置では、本発明のような入力データの特徴量の外れ値を削除することはできなかった。
【０００８】
本発明は、上記事情に鑑みてなされたものであり、異常検出システムにおいて、学習モデルを作成する入力データの質を向上させることにより、異常検出の精度を向上させることが可能な異常検出システム、及び異常検出プログラムの提供を目的とする。
【課題を解決するための手段】
【０００９】
上記目的を達成するため、本発明の異常検出システムは、監視対象の異常を検出する異常検出システムであって、監視対象に関するデータを用いて特徴量を計算して得られた特徴量データに対して状態空間モデルへの写像手法を用いて前記特徴量データに対応する特徴空間を定義する基底を生成し、前記基底にもとづき前記特徴量データを写像して、前記特徴空間の新たな次元数を決定し、写像された前記特徴量データの次元数を削減して前記新たな次元数とし、変換後特徴量データを生成する写像部と、混合分布手法にもとづき前記変換後特徴量データを複数のクラスタに分類して、学習モデル候補を作成する学習モデル作成部と、前記学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除することにより、前記学習モデル候補を改変して学習モデルを作成するモデル改変部を備える構成としてある。
【００１０】
また、本発明の異常検出プログラムは、監視対象の異常を検出する異常検出プログラムであって、コンピュータを、監視対象に関するデータを用いて特徴量を計算して得られた特徴量データに対して状態空間モデルへの写像手法を用いて前記特徴量データに対応する特徴空間を定義する基底を生成し、前記基底にもとづき前記特徴量データを写像して、前記特徴空間の新たな次元数を決定し、写像された前記特徴量データの次元数を削減して前記新たな次元数とし、変換後特徴量データを生成する写像部、混合分布手法にもとづき前記変換後特徴量データを複数のクラスタに分類して、学習モデル候補を作成する学習モデル作成部、及び、前記学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除することにより、前記学習モデル候補を改変して学習モデルを作成するモデル改変部として機能させる構成としてある。
【発明の効果】
【００１１】
本発明によれば、異常検出システムにおいて、学習モデルを作成する入力データの質を向上させることにより、異常検出の精度を向上させることが可能な異常検出システム、及び異常検出プログラムの提供が可能となる。
【図面の簡単な説明】
【００１２】
【図１】本発明の第一実施形態の異常検出システム（異常検出装置）の構成を示すブロック図である。
【図２】本発明の第一実施形態、第二実施形態、及び第三実施形態で用いられる主成分分析による写像手法を示す説明図である。
【図３】本発明の第一実施形態、第二実施形態、及び第三実施形態で用いられる主成分分析にもとづいて、特徴量データを低次元空間に圧縮する流れを示す説明図である。
【図４】本発明の第一実施形態及び第三実施形態で用いられる混合ガウス分布モデルによるモデル改変処理についての説明図である。
【図５】本発明の第一実施形態の異常検出システムによる処理手順を示すフローチャートである。
【図６】本発明の第二実施形態の異常検出システムの構成を示すブロック図である。
【図７】本発明の第二実施形態の異常検出システムによる処理手順を示すフローチャートである。
【図８】本発明の第三実施形態の異常検出システムの構成を示すブロック図である。
【図９】本発明の第三実施形態の異常検出システムによる処理手順を示すフローチャートである。
【図１０】実施例に示す試験で用いた画像データを示す図である。
【図１１】実施例に示す試験で用いた画像データから抽出した特徴量の種類を示す図である。
【図１２】実施例に示す試験における主成分分析の結果を表すグラフを示す図である。
【図１３】実施例及び参考例によって得られた学習モデルの評価結果を示す図である。
【発明を実施するための形態】
【００１３】
以下、本発明の異常検出システム、及び異常検出プログラムの実施形態について詳細に説明する。ただし、本発明は、以下の実施形態及び後述する実施例の具体的な内容に限定されるものではない。
【００１４】
［第一実施形態］
まず、本発明の第一実施形態の異常検出システム、及び異常検出プログラムについて、図面を参照して説明する。
本実施形態の異常検出システムは、監視対象の異常を検出するものである。この監視対象には検査対象も含まれ、監視対象としては、例えば製品や装置、建物などにおける環境等を挙げることができる。
【００１５】
具体的には、本実施形態の異常検出システムに相当する異常検出装置１は、図１に示すように、測定データ入力部１０、測定データ記憶部１１、特徴量抽出部１２、特徴量データ記憶部１３、写像部１４、写像データ記憶部１５、学習モデル作成部１６、学習モデル記憶部１７、及びモデル改変部１８を備えている。
【００１６】
本実施形態の異常検出装置１におけるこれらの各構成は、図１に示すように一個の情報処理装置に全て備えることができる。また、これらの各構成は、複数の情報処理装置からなる異常検出システムの各装置に分散して備えてもよい。これは、後述する各実施形態についても同様である。
【００１７】
測定データ入力部１０は、監視対象について測定されたデータ（以下、測定データと称する場合がある）などを入力して、測定データ記憶部１１に記憶させる。
測定データ記憶部１１は、測定データ入力部１０により入力された測定データを記憶する。
測定データとしては、製品などから取得される検査データと、装置などから取得されるプロセスデータを挙げることができる。
【００１８】
検査データとしては、画像、寸法（サイズ、円形度など）、色、欠陥（数、サイズなど）、排出数、透過度、重量、強度、及び表面粗さ等を挙げることができる。なお、検査データは、検査する項目だけ存在するため、これ以外にも多くのものが想定され得る。
プロセスデータとしては、音、温度、湿度、振動、変位、速度、加速度、回転角、角速度、角加速度、電圧、電流、光量、圧力等を挙げることができる。
測定指標として、このような検査データとプロセスデータを組合わせて用いることで、異常検知の性能の向上を図ることができる。
【００１９】
特徴量抽出部１２は、測定データ記憶部１１から測定データを読み込んで、この測定データにもとづき特徴量を計算し、特徴量データを生成して特徴量データ記憶部１３に記憶させる。
このとき、測定データを２つの集合に分けて、特徴量抽出部１２により、特徴量データとして、学習モデル候補を作成するために用いる学習用データセットと、学習モデル候補を改変して学習モデルを完成させるために用いる検証用データセットの２種類を生成させて、特徴量データ記憶部１３に記憶させることが好ましい。
【００２０】
なお、測定データ入力部１０、測定データ記憶部１１、及び特徴量抽出部１２は、本実施形態の異常検出システム及び異常検出プログラムの外部に構成することができ、本実施形態の異常検出システム及び異常検出プログラムから、これらの構成を省略してもよい。
【００２１】
特徴量データ記憶部１３は、特徴量抽出部１２により生成された特徴量データを記憶する。
特徴量データの構造は、例えば測定データ（例えば画像）ＩＤ、複数の特徴量ＩＤ等からなるものとすることができる。
【００２２】
特徴量としては、例えば測定データが画像である場合、監視対象の形状に関するものとして、大きさ，アスペクト比，円形度，周囲長，外接四角面積比，凸数等を挙げることができる。また、監視対象の輝度に関するものとして、色空間（Ｒ空間，Ｇ空間，Ｂ空間，ＧＲＡＹ空間），平均値，標準偏差，最大値，最小値，最大値−最小値，中央値，第1四分位数，第2四分位数，1σ確率，3σ確率，四分位範囲，平均値−中央値，ヒストグラム確率密度分布（歪度，尖度，コントラスト，2次モーメント，エントロピー），差分確率密度分布（コントラスト，2次モーメント，エントロピー）等を挙げることができる。
【００２３】
また、特徴量として、時系列データなどの波形データにおける測定値、波形データのある時間的幅の平均値、最大値、最小値、中央値、標準偏差、最大値-最小値、１σ確率、３σ確率、第１四分位数、第３四分位数、波形データのヒストグラムの最頻値、尖度、歪度、波形データの周波数特性におけるピーク位置、ピーク値、重心位置、田口法による変化量、存在量等を挙げることもできる。
【００２４】
写像部１４は、特徴量データ記憶部１３から特徴量データを読み込んで、この特徴量データに対して状態空間モデルへの写像手法を用いて、特徴量データに対応する特徴空間を定義する基底を生成し、この基底にもとづき特徴量データを写像する。このとき、特徴量データとしては、学習用データセットが好適に用いられる。
また、写像部１４は、設定された閾値にもとづいて、特徴空間の新たな次元数（圧縮後の次元数）を決定する。
【００２５】
さらに、写像部１４は、写像された特徴量データの次元数を削減して、新たな次元数とする次元削減処理を行う。
このとき、写像部１４は、写像された特徴量データの主成分が持つ固有値を上位から足し合わせて、その合計値が閾値未満で、かつ最大の値になるように新たな次元数を決定し、次元削減処理を行うことができる。
そして、写像部１４は、得られた変換後特徴量データを写像データとして、写像データ記憶部１５に記憶させる。
【００２６】
写像データ記憶部１５は、写像部１４により生成された変換後特徴量データを写像データとして記憶する。
変換後特徴量データの構造は、例えば測定データ（例えば画像）ＩＤ、複数の特徴量ＩＤ等からなるものとすることができる。
【００２７】
具体的には、写像部１４は、状態空間モデルへの写像手法として、主成分分析（ＰＣＡ；Principal Component Analysis）、独立成分分析、カーネル主成分分析、ＬＤＡ（Latent Dirichlet Allocation）、近傍探索手法の少なくともいずれかを用いることができ、主成分分析を好適に用いることができる。
【００２８】
主成分分析は、主部分空間（principal subspace）と呼ばれる低次元の線形空間の上へのデータ点の直交射影である。また、主成分分析は、もともとのデータ点と射影した点の間の２乗距離の平均値で定義される射影のコスト関数の期待値を最小化するような線形射影である。
主成分分析によれば、元のデータが持つ情報をできるだけ失わずに、低次元空間に情報を縮約することができる。すなわち、異常検出に必要のない無駄な特徴量を使うことなく、最低限の情報量だけで元のデータを表現することが可能である。
【００２９】
また、このとき、写像部１４により、特徴量データを変換後特徴量データに変換する、以下の変換式が作成される。
【数１】
2020144626.tif 000003
【００３０】
例えば、図２に示すように、主成分分析により、第１主成分として８０％、第２主成分として２０％が得られた場合において、設定された閾値が８０％である場合、写像部１４によって、特徴空間の新たな次元数（圧縮後の次元数）として１次元が決定されて変換後特徴量データが生成され、特徴量データを１次元の変換後特徴量データに変換する変換式が作成される。
また、図３に示すように、特徴量データＸに対応する特徴空間がＭ次元である場合には、写像処理によりＭ次元の特徴量データＸ’が生成され、次元削減処理によりＭ’次元（Ｍ＞Ｍ’）の特徴量データＸ’が生成される。
【００３１】
さらに、写像部１４は、特徴量データ記憶部１３から検証用データセットの特徴量データを読み込んで、この特徴量データに対して状態空間モデルへの写像手法を用いて、先に生成した基底にもとづき特徴量データを写像する。
このとき、写像部１４は、上記の変換式により、特徴量データを写像して、次元数を削減し、得られた変換後特徴量データを写像データとして、写像データ記憶部１５に記憶させることができる。
【００３２】
学習モデル作成部１６は、写像データ記憶部１５から変換後特徴量データ（写像データ）を読み込んで、混合分布手法にもとづき変換後特徴量データを複数のクラスタ（塊）に分類し、学習モデル候補を作成して、これを学習モデル記憶部１７に記憶させる。
学習モデル候補は、一又は二以上のクラスタからなり、クラスタの構造は、例えばクラスタＩＤ、各クラスの平均値、共分散行列、混合比等からなるものとすることができる。
【００３３】
本実施形態の異常検出システムでは、測定データを最終的に「正常」と「異常」の２種類に分類する。また、本実施形態の異常検出システムは、教師なし学習を採用しており、教師なし学習では「正常」データのみを使用して学習モデルを構築し、外れ値を「異常」と判定する。このため、異常を定義することなく、様々な種類の異常を検知することが可能となっている。
ところで、「正常」にも様々な種類が存在する可能性（多態性）があるところ、１つの均一な集団においては外れ値とならないデータが、１つの均一な集団を３つの集団とした場合には、データの分布を細かく表現することが可能になる結果、外れ値になることがある。
このため、本実施形態の異常検出システムでは、混合分布手法を用いて変換後特徴量データを複数のクラスタに分類することにより、異常検出の精度を向上させることが可能になっている。
【００３４】
混合分布手法としては、正規分布、t分布、ポアソン分布、ベルヌ−イ分布、ディリクレ分布の少なくともいずれかによる分布モデルを用いることができ、混合ガウス分布モデル（Gaussian Mixture Model, GMM）を好適に用いることができる。
混合ガウス分布モデルは、複数の正規分布を組み合わせた重み付け線形和で表されるモデルである。
【００３５】
ここで、教師データなしで、データ間の類似度あるいは非類似度を手掛かりに、いくつかのクラスタにグループ分けをすることをクラスタリングという。正規分布は単峰性の確率分布しか表現できないが、実際に観測される測定データを特徴空間に布置していくと、様々な種類の測定データが存在するため、複数のクラスタが構成されると考えられる。この場合、単峰性の確率分布を複数組み合わせた重み付け線形和で全体の確率分布をモデル化すると整合性の高いモデルが構築される。このような確率モデルを混合分布モデルといい、特に複数の正規分布を混ぜ合わせて構成された混合分布モデルを、混合ガウス分布モデルという。
【００３６】
混合ガウス分布モデルは、xを特徴量ベクトル、jを正規分布の番号、dを次元数としたとき、j番目の正規分布でのxの確率分布は、正規分布パラメータ（平均j、共分散行列j）により、以下の式で表される。
【数２】
2020144626.tif 000004
【００３７】
また、学習モデル作成部１７は、混合分布手法にもとづき変換後特徴量データを複数のクラスタに分類するにあたり、正規分布の各パラメータを少しずつ調整することで、特徴が似ている測定データに分けた状態を近似して表現する。これは、ＥＭアルゴリズムという手法によって実現される。ＥＭアルゴリズムは、不完全データ、すなわち検出する測定データを学習しておらず、その発生確率が分からないような場合に用いられる。このような直接求められない値がある場合に、繰り返し演算を適用することによって値を推定する最尤推定の計算アルゴリズムの一つである。
【００３８】
モデル改変部１８は、学習モデル記憶部１７における学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除すること（以下、クラスタを削除することを「クラスタカット」と称する場合がある）により、学習モデル候補を改変して学習モデルを作成し、学習モデル記憶部１７における学習モデル候補を学習モデルに更新する。
すなわち、本実施形態の異常検出システムによれば、例えば、図４に示すように、混合比が異常値（閾値未満）のクラスタが存在している場合、学習モデル候補におけるクラスタから、このような異常値のクラスタを削除することで、誤検出を抑制して、異常検出の精度を向上させることが可能になっている。
【００３９】
モデル改変部１８により作成される学習モデルは、後述するように、クラスタの混合数ごとに複数であってもよい。その後、最適なクラスタの混合数にもとづいて、最も優れた学習モデルを決定し、学習モデルを完成させることができる。
また、本実施形態の異常検出システムにおいて、このようにして完成した学習モデルを記憶させるために、学習モデル記憶部１７とは別個の記憶部を備えてもよい。
【００４０】
また、モデル改変部１８により、検証用データセットの変換後特徴量データを用いて異常判定処理を行い、その結果にもとづいて、クラスタカットを行うことも好ましい。
すなわち、モデル改変部１８は、学習モデル記憶部１７の学習モデル候補における各クラスタについて検証用データセットの変換後特徴量データが外れ値として生起する確率（生起確率）を計算し、生起確率が設定した閾値よりも低い場合に正常と判定し、この変換後特徴量データに対応する検証用データが本来は異常である場合（誤判定の場合）、当該クラスタを学習モデル候補から削除することが好ましい。
【００４１】
なお、モデル改変部１８によって、検証用データセットの特徴量データを入力して、写像部１４により生成された基底を用いて特徴量データを写像し、写像された特徴量データの次元数を削減して新たな次元数とし、変換後特徴量データを生成し、学習モデル記憶部１７の学習モデル候補における各クラスタについて変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも低い場合に正常と判定し、この変換後特徴量データに対応する検証用データが本来は異常である場合、当該クラスタを学習モデル候補から削除する構成とすることも可能である。
【００４２】
ここで、モデル改変部１８による生起確率の計算手法は、同一の目的の処理を行い得る範囲において特に限定されないが、例えばマハラノビスタグチ（Maharanobis Taguchi, MT）法のマハラノビス距離を用いて好適に行うことができる。
具体的には、以下の式によりマハラノビス距離を算出することができる。
【数３】
2020144626.tif 000005
【００４３】
そして、得られたマハラノビス距離を以下の累積分布関数（ＣＤＦ）に適用することで、生起確率を計算することができる。なお、後述する試験においては、閾値（threshold）を0.99（カイ二乗分布の片側有意水準1％）としている。
【数４】
2020144626.tif 000006
【００４４】
また、マハラノビスタグチ法にかえて、ＭＴＡ（Maharanobis Taguchi Afjoint）法、ＲＴ（Recognition Taguchi）法、Taguchi法、及びTaguchi Schmidt法のいずれかを用いることもできる。
【００４５】
また、モデル改変部１８は、特徴量データと各クラスタのマハラノビス距離の計算にあたり、各クラスタの重心を用いることができる。
また、「本来は異常である場合」は、変換後特徴量データに対応する検証用データが、検証用データセットにおける異常データであることを意味する。
【００４６】
さらに、学習モデル作成部１７により、クラスタの混合数が相違する複数の学習モデル候補を作成させると共に、モデル改変部１８により、複数の学習モデル候補のうちの真陽性の発生率が最も高いものを学習モデルに決定させることも好ましい。
すなわち、学習モデル作成部１７は、設定された複数の混合数にもとづいて、変換後特徴量データ（写像データ）を混合数分のクラスタに分類して、混合数ごとに学習モデル候補を作成することが好ましい。
【００４７】
また、写像部１４により、特徴量データ（検証用データセット）を入力して、写像部１４により先に生成された基底を用いて特徴量データを写像し、写像された特徴量データの次元数を削減して新たな次元数とし、変換後特徴量データを生成することが好ましい。
そして、モデル改変部１８により、混合数ごとに、学習モデル候補における各クラスタについて変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも高い場合に外れ値となることで異常と判定し、真陽性の発生率が最も高い学習モデル候補を学習モデルに決定することも好ましい。
【００４８】
なお、モデル改変部１８によって、特徴量データ（検証用データセット）を入力して、写像部１４により生成された基底を用いて特徴量データを写像し、写像された特徴量データの次元数を削減して新たな次元数とし、変換後特徴量データを生成し、混合数ごとに、学習モデル候補における各クラスタについて変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも高い場合に外れ値となることで異常と判定し、真陽性の発生率が最も高い学習モデル候補を学習モデルに決定する構成とすることも可能である。
【００４９】
真陽性の発生率は、例えば、検証用データセットにおいて事前に正常データと異常データを準備しておき、全ての異常データにおける異常データと判定されたものの割合として算出することができる。また、偽陽性の発生率は、全ての正常データにおける異常データと誤判定されたものの割合として算出することができる。
【００５０】
なお、本実施形態の異常検出システムは、学習モデルを作成するまでの構成を備えており、測定データの異常検出を行うものではないが、測定データの異常検出を行うことができる学習モデルを作成しているため、「異常検出システム」と称している。
【００５１】
次に、本実施形態の異常検出システムによる処理手順について、図５を参照して説明する。
まず、測定データ入力部１０が、測定データを入力して、測定データ記憶部１１に記憶させる（ステップ１０）。
次いで、特徴量抽出部１２が、測定データ記憶部１１における測定データを読み込んで、特徴量を抽出し（ステップ１１）、特徴量データを生成して、特徴量データ記憶部１３に記憶させる。このとき、特徴量抽出部１２は、特徴量データを学習用データセットと検証用データセットの２種類に分けて、特徴量データ記憶部１３に記憶させることができる（ステップ１２）。
【００５２】
次に、写像部１４が、特徴量データ記憶部１３における特徴量データを読み込んで、写像部による処理を実行する（ステップ１３）。
このとき、写像部１４は、学習用データセットの写像処理と次元削減を行う。
具体的には、写像部１４は、主成分分析などにより、特徴量データに対応する特徴空間を定義する基底を生成して、この基底にもとづき特徴量データを写像する。また、設定された閾値にもとづき特徴空間の新たな次元数（圧縮後の次元数）を決定する。そして、写像された特徴量データの次元数を削減して新たな次元数とし、得られた変換後特徴量データを写像データとして、写像データ記憶部１５に記憶させる。
【００５３】
次に、学習モデル作成部１６が、学習モデル候補を作成して、学習モデル記憶部１７に記憶させる（ステップ１４）。
このとき、学習モデル作成部１６は、写像データ記憶部１５から変換後特徴量データ（写像データ）を読み込んで、混合ガウス分布モデルなどにより、変換後特徴量データを複数のクラスタに分類し、学習モデル候補を作成して、これを学習モデル記憶部１７に記憶させる。
【００５４】
次に、モデル改変部１８が、学習モデル候補に対してクラスタカットを行うことにより、学習モデルを作成し、最終的な学習モデルを決定する（ステップ１５）。
このとき、モデル改変部１８が、学習モデル候補に対して、設定された閾値よりも混合比が低い値をとるクラスタ又は混合比が最も小さい値をとるクラスタを削除することにより、学習モデル候補を改変して学習モデルを作成することが好ましい。
【００５５】
また、モデル改変部１８が、学習モデル候補における各クラスタについて変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも低い場合に正常と判定し、特徴量データに対応する測定データが本来は異常である場合、当該クラスタを学習モデル候補から削除することにより、学習モデル候補を改変して学習モデルを作成することも好ましい。
【００５６】
さらに、モデル改変部１８が、混合数ごとに、学習モデル候補における各クラスタについて変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも高い場合に外れ値となることで異常と判定し、真陽性の発生率が最も高い学習モデル候補を最終的な学習モデルとして決定することも好ましい。
【００５７】
このような本実施形態の異常検出システム及び異常検出プログラムによれば、学習モデル候補から混合比が一定以下の低い値をとるクラスタなどを削除して学習モデルを作成することができる。このため、実質的に入力データの質を向上させることができ、異常検出の精度を向上させることが可能である。
【００５８】
［第二実施形態］
次に、本発明の第二実施形態に係る異常検出システム、及び異常検出プログラムについて、図６及び図７を参照して説明する。
本実施形態は、第一実施形態により完成した学習モデルを用いて、測定データの異常判定を行う装置に関するものである。
【００５９】
本実施形態の異常検出システムに相当する異常検出装置１ａは、図６に示すように、測定データ入力部１０ａ、測定データ記憶部１１ａ、特徴量抽出部１２ａ、特徴量データ記憶部１３ａ、写像部１４ａ、写像データ記憶部１５ａ、学習モデル記憶部１７ａ、異常判定部１９ａ、及び判定結果記憶部２０ａを備えている。
【００６０】
本実施形態の異常検出装置１ａにおける測定データ入力部１０ａ、測定データ記憶部１１ａ、特徴量データ記憶部１３ａ、写像データ記憶部１５ａ、学習モデル記憶部１７ａは、第一実施形態における同一名称の構成と同様のものとすることができる。また、本実施形態の異常検出装置１ａにおけるその他の構成についても、以下に説明する点を除いて、第一実施形態における同一名称の構成と同様のものとすることができる。
【００６１】
特徴量抽出部１２ａは、測定データ記憶部１１ａから測定データを読み込んで、この測定データにもとづき特徴量を計算し、特徴量データ（測定データセット）を生成して特徴量データ記憶部１３ａに記憶させる。
【００６２】
写像部１４ａは、特徴量データ記憶部１３ａから特徴量データを読み込んで、この特徴量データに対して状態空間モデルへの写像手法を用いて、設定された基底にもとづき特徴量データを写像し、写像された特徴量データの次元数を削減して新たな次元数として、変換後特徴量データを生成する。
このとき、写像部１４ａは、第一実施形態において作成された特徴量データを変換後特徴量データに変換する変換式を用いて、特徴量データを変換後特徴量データに変換し、変換後特徴量データを写像データとして、写像データ記憶部１５ａに記憶させることができる。なお、本実施形態の異常検出システム１ａに変換式記憶部を設けて、上記の変換式を予め記憶させておくことができる。
【００６３】
異常判定部１９ａは、学習モデル記憶部１７ａの学習モデル候補における各クラスタについて、写像データ記憶部１５ａにおける変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも高い場合に外れ値となることで、対応する測定データを異常と判定して、この判定結果を判定結果記憶部２０ａに記憶させる。
異常判定部１９ａによる生起確率の計算は、第一実施形態のモデル改変部１８と同一の手法で行うことができ、例えばマハラノビスタグチ（Maharanobis Taguchi, MT）法のマハラノビス距離を用いて好適に行うことができる。
【００６４】
次に、本実施形態の異常検出システムの処理手順について、図７を参照して説明する。
まず、測定データ入力部１０ａが、測定データを入力して、測定データ記憶部１１ａに記憶させる（ステップ３０）。
次いで、特徴量抽出部１２ａが、測定データ記憶部１１ａにおける測定データを読み込んで、特徴量を抽出し（ステップ３１）、特徴量データを（測定データセット）生成して、特徴量データ記憶部１３ａに記憶させる（ステップ３２）。
【００６５】
次に、写像部１４ａが、特徴量データ記憶部１３ａにおける特徴量データを読み込んで、写像部による処理を実行する（ステップ３３）。
このとき、写像部１４ａは、測定データセットの写像処理と次元削減を行う。
具体的には、写像部１４ａは、上記の変換式を用いて、特徴量データを変換後特徴量データに変換し、変換後特徴量データを写像データとして写像データ記憶部１５ａに記憶させる。
【００６６】
次に、異常判定部１９ａが、測定データの異常判定処理を行う（ステップ３４）。
このとき、異常判定部１９ａは、学習モデル記憶部１７ａの学習モデルにおける各クラスタについて、写像データ記憶部１５ａにおける変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも高い場合に外れ値となることで、対応する測定データを異常と判定して、この判定結果を判定結果記憶部２０ａに記憶させる。
【００６７】
このような本実施形態の異常検出システム及び異常検出プログラムによれば、例えば第一実施形態の異常検出システムにより学習モデルと特徴量データを変換後特徴量データに変換する変換式とを作成し、これらを本実施形態の異常検出システムに記憶させることにより、本実施形態の異常検出システムによって測定データの異常判定を行うことが可能である。
【００６８】
［第三実施形態］
次に、本発明の第三実施形態に係る異常検出システム、及び異常検出プログラムについて、図８及び図９を参照して説明する。
本実施形態は、第一実施形態と第二実施形態を組み合わせて、第一実施形態の異常検出システムに異常判定を行うための構成を追加したものである。
【００６９】
本実施形態の異常検出システムに相当する異常検出装置１ｂは、図８に示すように、測定データ入力部１０ｂ、測定データ記憶部１１ｂ、特徴量抽出部１２ｂ、特徴量データ記憶部１３ｂ、写像部１４ｂ、写像データ記憶部１５ｂ、学習モデル作成部１６ｂ、学習モデル記憶部１７ｂ、モデル改変部１８ｂ、異常判定部１９ｂ、及び判定結果記憶部２０ｂを備えている。
【００７０】
本実施形態の異常検出装置１ｂにおける測定データ入力部１０ｂ、測定データ記憶部１１ｂ、特徴量抽出部１２ｂ、特徴量データ記憶部１３ｂ、写像部１４ｂ、写像データ記憶部１５ｂ、学習モデル作成部１６ｂ、学習モデル記憶部１７ｂ、モデル改変部１８ｂは、第一実施形態における同一名称の構成と同様のものとすることができる。また、本実施形態の異常検出装置１ｂにおける異常判定部１９ｂ、及び判定結果記憶部２０ｂは、第二実施形態における同一名称の構成と同様のものとすることができる。
【００７１】
すなわち、学習モデルの作成段階において、特徴量抽出部１２ｂは、測定データ記憶部１１ｂから測定データを読み込んで、この測定データにもとづき特徴量を計算し、特徴量データを生成して特徴量データ記憶部１３ｂに記憶させる。
このとき、測定データを２つの集合に分けて、特徴量抽出部１２ｂにより、特徴量データとして、学習モデル候補を作成するために用いる学習用データセットと、学習モデル候補を改変して学習モデルを完成させるために用いる検証用データセットの２種類を生成させて、特徴量データ記憶部１３ｂに記憶させる。
【００７２】
また、写像部１４ｂは、特徴量データ記憶部１３ｂから特徴量データ（学習用データセット）を読み込んで、この特徴量データに対して状態空間モデルへの写像手法を用いて、特徴量データに対応する特徴空間を定義する基底を生成し、この基底にもとづき特徴量データを写像する。また、写像部１４ｂは、設定された閾値にもとづいて、特徴空間の新たな次元数（圧縮後の次元数）を決定する。さらに、写像部１４ｂは、写像された特徴量データの次元数を削減して、新たな次元数とする次元削減処理を行う。そして、写像部１４ｂは、得られた変換後特徴量データを写像データとして、写像データ記憶部１５ｂに記憶させる。また、このとき、写像部１４ｂにより、特徴量データを変換後特徴量データに変換する変換式が作成される。
【００７３】
また、学習モデルの完成後、異常検出段階において、特徴量抽出部１２ｂは、測定データ記憶部１１ｂから測定データを読み込んで、この測定データにもとづき特徴量を計算し、特徴量データ（測定データセット）を生成して特徴量データ記憶部１３ｂに記憶させる。
【００７４】
写像部１４ｂは、特徴量データ記憶部１３ｂから特徴量データを読み込んで、この特徴量データに対して状態空間モデルへの写像手法を用いて、設定された基底にもとづき特徴量データを写像し、写像された特徴量データの次元数を削減して新たな次元数として、変換後特徴量データを生成する。
このとき、写像部１４ｂは、上記の変換式を用いて、特徴量データを変換後特徴量データに変換し、変換後特徴量データを写像データとして、写像データ記憶部１５ｂに記憶させることができる。
【００７５】
次に、本実施形態の異常検出システムによる処理手順について、図９を参照して説明する。
まず、測定データ入力部１０ｂが、測定データを入力して、測定データ記憶部１１ｂに記憶させる（ステップ５０）。
次いで、特徴量抽出部１２ｂが、測定データ記憶部１１ｂにおける測定データを読み込んで、特徴量を抽出し（ステップ５１）、特徴量データを生成して、特徴量データ記憶部１３ｂに記憶させる。このとき、特徴量抽出部１２ｂは、特徴量データを学習用データセットと検証用データセットの２種類に分けて、特徴量データ記憶部１３ｂに記憶させる（ステップ５２）。
【００７６】
次に、写像部１４ｂが、特徴量データ記憶部１３ｂにおける特徴量データを読み込んで、写像部による処理を実行する（ステップ５３）。このとき、写像部１４ｂは、学習用データセットの写像処理と次元削減を行う。
次いで、学習モデル作成部１６ｂが、学習モデル候補を作成して、学習モデル記憶部１７ｂに記憶させる（ステップ５４）。
【００７７】
次に、モデル改変部１８ｂが、学習モデル候補に対してクラスタカットを行うことにより、学習モデルを作成し、最終的な学習モデルを決定する（ステップ５５）。
そして、異常検出段階において、測定データ入力部１０ｂが、測定データを入力して、測定データ記憶部１１ｂに記憶させる（ステップ５６）。
次いで、特徴量抽出部１２ｂが、測定データ記憶部１１ｂにおける測定データを読み込んで、特徴量を抽出し（ステップ５７）、特徴量データを（測定データセット）生成して、特徴量データ記憶部１３ｂに記憶させる（ステップ５８）。
【００７８】
次に、写像部１４ｂが、特徴量データ記憶部１３ｂにおける特徴量データを読み込んで、写像部による処理を実行する（ステップ５９）。このとき、写像部１４ｂは、測定データセットの写像処理と次元削減を行う。
最後に、異常判定部１９ｂが、測定データの異常判定処理を行う（ステップ６０）。
このとき、異常判定部１９ｂは、学習モデル記憶部１７ｂの学習モデルにおける各クラスタについて、写像データ記憶部１５ｂにおける変換後特徴量データが生起する生起確率を計算し、生起確率が設定した閾値よりも高い場合に外れ値となることで、対応する測定データを異常と判定して、この判定結果を判定結果記憶部２０ｂに記憶させる。
【００７９】
このような本実施形態の異常検出システム及び異常検出プログラムによれば、学習モデル候補から混合比が一定以下の低い値をとるクラスタなどを削除して学習モデルを作成することができる。また、この異常検出システムをそのまま異常検出に用いることができる。
これにより、本実施形態によれば、実質的に入力データの質を向上させ、異常検出の精度を向上させることが可能となっている。
【００８０】
上記実施形態の異常検出システムは、本発明の異常検出プログラムに制御されたコンピュータを用いて実現することができる。コンピュータのＣＰＵは、異常検出プログラムにもとづいてコンピュータの各構成要素に指令を送り、異常検出システムの動作に必要となる所定の処理、例えば、測定データの入力処理、特徴量の抽出処理、データセットの作成処理、写像部による処理、学習モデルの生成処理、モデル改変処理、異常判定処理等を行わせる。このように、本発明の異常検出システムにおける各処理、動作は、プログラムとコンピュータとが協働した具体的手段により実現できるものである。
【００８１】
プログラムは予めＲＯＭ、ＲＡＭ等の記録媒体に格納され、コンピュータに実装された記録媒体から当該コンピュータにプログラムを読み込ませて実行されるが、例えば通信回線を介してコンピュータに読み込ませることもできる。
また、プログラムを格納する記録媒体は、例えば半導体メモリ、磁気ディスク、光ディスク、その他任意のコンピュータで読取り可能な任意の記録手段により構成できる。
【実施例】
【００８２】
以下、本発明の実施形態に係る異常検出システム及び異常検出プログラムを用いて学習モデルを作成し、その効果を評価した試験について、図１０〜図１３を用いて説明する。
実施例として、本発明の第一実施形態に係る異常検出システムを使用した。また、参考例として、本発明の第一実施形態に係る画像処理装置の一部を変更したものを使用した。具体的には、参考例ではモデル改変部１８を省略して、クラスタカットを行わないように変更した。
【００８３】
本試験において、画像データを測定データとして用いるため、画像データを準備した。
具体的には、図１０に示すように、形状と色の異なる９種類の画像を含む画像クラスを用意した。各画像クラスにおける画像のサイズは、５０×５０ピクセルであり、画像中央の２０×２０〜３４×３４ピクセル内に検出対象が描かれている。
【００８４】
画像クラス１〜３の検出対象の形状は円形であり、画像クラス４〜６の検出対象の形状は三角形であり、画像クラス７〜９の検出対象の形状は四角形である。
また、画像クラス１，４，７の検出対象の色は赤系であり、ＲＧＢ空間においてＲの値がＧ，Ｂより大きいことを示している。また、画像クラス２，５，８の検出対象の色は緑系であり、ＲＧＢ空間においてＧの値がＲ，Ｂより大きいことを示している。画像クラス３，６，９の検出対象の色は青系であり、ＲＧＢ空間においてＢの値がＲ，Ｇより大きいことを示している。なお、代表画像は、それぞれの画像クラスにおける代表的な画像の一例を示している。
【００８５】
また、画像クラス１，３〜９に含まれる画像は、全て正常のものであり、画像クラス２に含まれる画像は、全て異常のものであるとしている。
具体的には、学習用データセットを作成するための画像クラス１，３〜９に含まれる画像は各５００枚であり、ステータスは全て正常としている。また、学習用データセットを作成するための画像クラス２に含まれる画像は各２８０枚であり、ステータスは全て異常としている。
【００８６】
さらに、検証用データセットを作成するための画像クラス１，３〜９に含まれる画像は各１００枚であり、ステータスは全て正常としている。また、検証用データセットを作成するための画像クラス２に含まれる画像は各８００枚であり、ステータスは全て異常としている。
【００８７】
上記において正常、異常の違いは、画像クラスに対してステータスを恣意的に決定している。
すなわち、学習用データセットにおいてステータスが異常の画像２８０枚は、全４２８０枚の画像中の７％程度であり、本実施形態の異常検出システムによって、これらが異常と判定される学習モデルが作成されることを想定している。
【００８８】
各画像の背景は、黒の濃淡によるグレー色（0〜25輝度（ランダム）に1ピクセルごとに決定）を設定し、ガウシアンフィルタを加えてぼやかすことにより描かれている。
また、全ての画像クラスの画像における検出対象に対し、ノイズがランダムに加えられている。具体的には、グレー色（0〜25輝度（ランダム））が1ピクセルごとに20％の確率で加えられている。このように、正常データと異常データの双方にノイズを加えることによって、学習モデルによる異常検出の難易度を向上させている。なお、これらのノイズは、特徴量抽出部１２による特徴量の抽出処理、及び写像部１４による主成分分析における次元削減処理によって、除去又は低減され得ると推定される。
【００８９】
特徴量抽出部１２による特徴量の抽出において対象とする特徴量は、図１１に示すように、形状に関する特徴量６種、輝度の特徴量を算出する色空間が４種、輝度の特徴量が２０種であり、全部で８６種類（＝6＋4×20）とした。
【００９０】
そして、写像部１４により学習用データセットを用いて主成分分析を実施して、基底を生成した。また、生成した基底を用いて学習用データセットを写像し、写像した学習用データセットの主成分が持つ固有値を上位から足し合わせて、その合計値が９５％未満で最大の値になるように、圧縮後の次元数を決定して、次元の削減を行った。
具体的には、図１２に示す特徴量８６種類に対応する８６次元が、１１次元に圧縮された。また、検証用データセットも同様に写像した後、同じ次元数に圧縮した。
【００９１】
次に、学習モデル作成部１６により、圧縮された学習用データセット（変換後特徴量データ）を用いて、混合ガウス分布モデルによる学習モデル候補を生成した。このとき、クラスタの混合数を１〜２０の２０種類とした。
【００９２】
そして、実施例では、モデル改変部１８により、圧縮された検証用データセット（変換後特徴量データ）を用いて、混合数ごとに、学習モデル候補による異常判定処理を行った。また、モデル改変部１８により、混合数ごとに、混合比率が一番小さいクラスと、偽陰性（異常を正常と誤判定）が発生したクラスを削除することにより、クラスタカットを実施した。なお、クラスタカットは、閾値を片側有意水準0％として計算することにより行い、これによってそのクラスに存在したデータを全て外れ値とした。
一方、参考例では、このようなクラスタカットは、行わなかった。
【００９３】
以上によって得られた実施例と参考例のそれぞれの学習モデル候補について、混合数ごとに、真陽性率（異常を異常と判定した割合）と偽陽性率（正常を異常と誤判定した割合）を算出した。その結果を図１３に示す。
【００９４】
異常検出システムにおいては、異常を異常と的確に判定することが重要であるところ、これらの図に示されるように、実施例では、ほとんどの混合数について、真陽性率が１．００であり、ほぼ全ての異常データを異常と判定できていることが分かる。これに対して、参考例は、全ての混合数について、真陽性率の点で実施例に大きく劣っている。
【００９５】
一方、実施例ではクラスタカットを行うことにより、偽陽性率が低下しており、真陽性率と偽陽性率はトレードオフの関係を示している。この傾向は混合数が低い場合に顕著に表れているが、混合数が９を超えるあたりから偽陽性率の低下が抑制されており，本試験では混合数１８のときに最も良い性能が得られている。
したがって、本実施形態の異常検出システムにおけるモデル改変部１８は、この場合、混合数１８の学習モデル候補を最終的な学習モデルとして好適に決定することができる。
【００９６】
本発明は、以上の実施形態に限定されるものではなく、本発明の範囲内において、種々の変更実施が可能であることは言うまでもない。
例えば、異常検出システムにおける各構成を複数の情報処理装置に分散して持たせたり、特徴量として実施形態及び実施例に示したもの以外のものを含めて使用したりするなど適宜変更することが可能である。
【産業上の利用可能性】
【００９７】
本発明は、監視対象の装置などから取得されるプロセスデータや、製品などから取得される検査データ等を用いて監視対象の異常の有無の検査などを行う場合に好適に利用することが可能である。
【符号の説明】
【００９８】
１、１ａ、１ｂ  異常検出装置
１０、１０ａ、１０ｂ  測定データ入力部
１１、１１ａ、１１ｂ  測定データ記憶部
１２、１２ａ、１２ｂ  特徴量抽出部
１３、１３ａ、１３ｂ  特徴量データ記憶部
１４、１４ａ、１４ｂ  写像部
１５、１５ａ、１５ｂ  写像データ記憶部
１６、１６ｂ  学習モデル作成部
１７、１７ａ、１７ｂ  学習モデル記憶部
１８、１８ｂ  モデル改変部
１９ａ、１９ｂ  異常判定部
２０ａ、２０ｂ  判定結果記憶部
【図１】
2020144626.tif 000007
【図２】
2020144626.tif 000008
【図３】
2020144626.tif 000009
【図４】
2020144626.tif 000010
【図５】
2020144626.tif 000011
【図６】
2020144626.tif 000012
【図７】
2020144626.tif 000013
【図８】
2020144626.tif 000014
【図９】
2020144626.tif 000015
【図１０】
2020144626.tif 000016
【図１１】
2020144626.tif 000017
【図１２】
2020144626.tif 000018
【図１３】
2020144626.tif 000019