(57)【特許請求の範囲】
【請求項１】
画像中の所定個の特徴量間の類似度を算出する算出手段と、
前記所定個の特徴量の各々について、自身及び他の特徴量との前記類似度を要素とするベクトルを生成し、前記特徴量間の類似度に基づいて前記ベクトルに対応する点が原点から等距離に並ぶように当該ベクトルの次元数を削減する次元削減手段と、
前記特徴量間の類似度を基に、前記所定個の特徴量のそれぞれに対応する軸の配置を決定する形状決定手段と、
前記形状決定手段により決定した軸の配置に従って前記画像から抽出した特徴量を描画した図を出力する出力手段とを有することを特徴とする情報処理装置。
【請求項２】
前記算出手段は、前記画像に含まれる各特徴量が示すデータの分布を基に前記特徴量間の類似度を算出することを特徴とする請求項１記載の情報処理装置。
【請求項３】
前記形状決定手段は、前記所定個の特徴量のそれぞれに対応する軸を、前記特徴量間の類似度に応じた角度をなすように配置することを特徴とする請求項１又は２記載の情報処理装置。
【請求項４】
前記所定個の特徴量のそれぞれに対応する軸の配置をユーザからの入力に応じて変更することを特徴とする請求項１〜３の何れか１項に記載の情報処理装置。
【請求項５】
前記特徴量間の類似度が所定の値よりも高い特徴量に対応する軸を１つにまとめて配置することを特徴とする請求項１〜４の何れか１項に記載の情報処理装置。
【請求項６】
画像中の所定個の特徴量間の類似度を算出する算出工程と、
前記所定個の特徴量の各々について、自身及び他の特徴量との前記類似度を要素とするベクトルを生成し、前記特徴量間の類似度に基づいて前記ベクトルに対応する点が原点から等距離に並ぶように当該ベクトルの次元数を削減する次元削減工程と、
前記特徴量間の類似度を基に、前記所定個の特徴量のそれぞれに対応する軸の配置を決定する形状決定工程と、
決定した軸の配置に従って前記画像から抽出した特徴量を描画した図を出力する出力工程とを有することを特徴とする情報処理方法。
【請求項７】
画像中の所定個の特徴量間の類似度を算出する算出ステップと、
前記所定個の特徴量の各々について、自身及び他の特徴量との前記類似度を要素とするベクトルを生成し、前記特徴量間の類似度に基づいて前記ベクトルに対応する点が原点から等距離に並ぶように当該ベクトルの次元数を削減する次元削減ステップと、
前記特徴量間の類似度を基に、前記所定個の特徴量のそれぞれに対応する軸の配置を決定する形状決定ステップと、
決定した軸の配置に従って前記画像から抽出した特徴量を描画した図を出力する出力ステップとをコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、情報処理装置、情報処理方法、及びプログラムに関する。
【背景技術】
【０００２】
人がデータ解析を行う際、数字の羅列をそのまま眺めるのではなくデータを人が把握しやすい形態で表示するための視覚表現として多様なグラフ（graph）又はチャート（chart）による表現方法を使うのが一般的である。対象となるデータや把握したい情報の性質によって適切に採用された可視化方法でデータを把握しやすくすることにより単なる数字の羅列では気付かないデータの傾向や情報を読み取ることが可能となる。ビッグデータの活用によって新しい問題解決方法や付加価値を提供することが昨今の技術領域のメインストリームになりつつある。例えば、多変量のデータを効果的に可視化して人の目による分析をサポートすることが重要となっている。
【０００３】
一般的に、データを複数の項目で分析する場合、円グラフ（Pie Chart）や棒グラフ（Bar Chart、Bar Graph）等が用いられる。また、レーダチャート（Radar Chart）や平行座標プロット（Parallel Coordinate Plot）等が用いられる。これらのうち、棒グラフやレーダチャート、平行座標プロットは、多変量データの関係をそのまま知ることができる方法である。これら３種のグラフは、さまざまな評価に使用できるデータ可視化方法であり、一般的には多変量を比較するのに使用される。また、グラフ内の複数の項目同士を比較することもでき、同じ項目で評価値を出した他のグラフを比較することもできる。
【０００４】
棒グラフや平行座標プロット、レーダチャートは、多変量解析に利用されるデータ可視化法である一方、変量の数の増加とともに軸数が増えてグラフ形状が複雑化していくため、人が一度に全体を把握することが困難となるということが知られている。認知神経科学で良く知られた事実として『マジック・ナンバー  ７±２  説』（非特許文献１）というものがある。これは、人が情報を認識する際、７±２のチャンクに納まる情報であればうまく認識できるが、それ以上になると難しくなるという説である。この説によると、例えばグラフ形状の直感的な比較に優れているレーダチャートであっても、複数のグラフ形状比較をする際、隣り合う軸同士の意味的なまとまり（チャンク）のない一般的なレーダチャートの場合は各軸が１チャンクとなる。そのため、軸数が７±２を超えたレーダチャートは人の形状認識能力を超えている可能性が高いことになる。
【０００５】
下記特許文献１には、多変量データの可視化の際にデータを表現する特徴量の数の増加とともに横方向に伸びていく平行座標プロットにおいて、部分的な平行座標プロットを取り出して配置し、人の認識しやすいグラフを作成する技術が提案されている。また、下記非特許文献２には、レーダチャートの軸に割りつける変数の順番が重要であると説き、これらをＰＣＡによって近くなったものをレーダチャート上での軸で近くになるように設定することが説明されている。
【先行技術文献】
【特許文献】
【０００６】
【特許文献１】特開２０１３−１６１２２６号公報
【非特許文献】
【０００７】
【非特許文献１】“The Magical Number Seven, Plus or Minus Two: Some Limits on Our Capacity for Processing Information", George A. Miller, The Psychological Review , 1956
【非特許文献２】『多変量解析事例集第１集』（吉澤・芳賀編，日科技連，1992）
【発明の概要】
【発明が解決しようとする課題】
【０００８】
ここで、レーダチャート、平行座標プロット、棒グラフや棒グラフを折れ線表示したものすべてにおいて言える課題は、データ可視化を形状として把握しやすく可視化している傍ら、軸の相対的な位置や間隔（幅や角度）に情報価値がないことである。これらのグラフは、多変量の項目間を連結した線によって（棒グラフの場合は棒の並びによって）形状が構成されるにも関わらず、偶然隣り合った項目を結ぶことによって形状が確定する。そのため、可視化によって得られる情報において人の認識にとっては隣り合う軸同士の関係性を重点的に認識することになる。例えば、棒グラフでは遠くの棒同士の高さの差は近くの棒同士の差に比べて直感的な大小把握が難しくなり、平行座標プロットでは隣同士の変数間に相関があるかどうかしか判らない。よって、詳細にスコアの差を比較する必要のない項目を離したり、微小な差異を詳細に比較すべき項目を近くにしたりするといった、多変量データを効果的に可視化する際の工夫が必要となってくる。
【０００９】
前記特許文献１に記載の平行座標プロットを部分的に取り出して配置する方法では、人が注目すべき部分が絞られるため、部分的な平行座標プロットに関して比較しやすくなる。しかし、グラフ全体を俯瞰してデータを把握するような場合には、多変量の部分的組合せを複数可視化して一つのグラフとするため、一瞥して情報を直感的に把握するのには不向きなグラフになったり、見落としが発生しやすくなったりする可能性がある。また、前記非特許文献２のレーダチャートの変数の順序決定方法によっても、近くにまとめるべき変数成分の見当がつくだけであり、それら変数を表す軸間の類似度等を反映したチャートを得ることができない。
【００１０】
本発明は、このような事情に鑑みてなされたものであり、構成する軸同士の関連性や類似性が考慮された比較認識を容易に行うことが可能になるチャート（グラフ）を生成できるようにすることを目的とする。
【課題を解決するための手段】
【００１１】
本発明に係る情報処理装置は、画像中の所定個の特徴量間の類似度を算出する算出手段と、前記所定個の特徴量の各々について、自身及び他の特徴量との前記類似度を要素とするベクトルを生成し、前記特徴量間の類似度に基づいて前記ベクトルに対応する点が原点から等距離に並ぶように当該ベクトルの次元数を削減する次元削減手段と、前記特徴量間の類似度を基に、前記所定個の特徴量のそれぞれに対応する軸の配置を決定する形状決定手段と、前記形状決定手段により決定した軸の配置に従って前記画像から抽出した特徴量を描画した図を出力する出力手段とを有することを特徴とする。
【発明の効果】
【００１２】
本発明によれば、一瞥しただけでは把握しにくい多変量データであっても直感的に把握しやすい表現形態にて表示することができ、データ分析がしやすくなる効果が得られる。
【図面の簡単な説明】
【００１３】
【図１】本実施形態における情報処理装置としての多変量データ可視化装置の構成例を示す図である。
【図２】本実施形態における処理動作の例を示すフローチャートである。
【図３】異常パターン毎にレーダチャートの形状が異なる例を示す図である。
【図４】異なる性質のデータの判別に不向きなレーダチャートの例を示す図である。
【図５】図４に示したレーダチャートを、異なる性質のデータを判別しやすく変換を施した例を示す図である。
【図６】第１の実施形態に係るレーダチャートの生成手順の例を示す図である。
【図７】第１の実施形態によるグラフの例を示す図である。
【図８】第１の実施形態に係るレーダチャートの生成手順の他の例を示す図である。
【図９】第１の実施形態に係るレーダチャートの生成手順の他の例を示す図である。
【図１０】第１の実施形態に係る表示形式を選択可能なＧＵＩの例を示す図である。
【図１１】本実施例における情報処理装置を実現可能なコンピュータ機能を示す図である。
【発明を実施するための形態】
【００１４】
以下、本発明の実施形態を図面に基づいて説明する。
【００１５】
（第１の実施形態）
本発明の第１の実施形態について説明する。第１の実施形態は、複数の特徴量で表現される複数のデータが有り、各々のデータがそれらの中でどのような傾向にあるかを分析したい場合に、人が見て直感的に判り易い形でのレーダチャートやそれに類するチャートを生成する方法である。
【００１６】
一例として、入力画像から複数の特徴量を抽出して自動で検査を行う外観検査でのユースケースを挙げて説明する。入力画像の正常・異常を自動で判定する外観検査では、入力画像から判定に必要な特徴量を複数抽出し、それらのスコア（評価値）を総合的に判断して正常画像か欠陥画像かを判定する。抽出される特徴量の示す各スコアと検出される異常の種類には相関があり、生産現場では判定結果である正常・異常のラベル以外にも判定スコアや元の抽出特徴軸におけるスコアにより測定の難しい欠陥の傾向を分析することができる。これによって、例えば生産ライン設計にフィードバックをかけることができる。
【００１７】
特に、人の目による官能目視検査を自動化した外観検査装置では複雑なパターンかつ多様なレベルの異常データが存在するため、異常と判定した画像だけを見ても何を根拠に異常と判定したのか分析するのが難しい。そこで、抽出特徴を一覧性の高いチャート（グラフ）によって可視化することで、画像を見るだけではわかりにくかった傾向の把握を容易にすることができる。
【００１８】
図３に一例として、各々の入力画像から抽出した特徴量のスコアの傾向の違いを直感的に知ることができる例としてレーダチャートの形状の違いが異常パターン毎に異なる様子を示した。抽出特徴の数が５個である場合を示しており、各々の抽出特徴にはＩＤが番号によって付与されている。図３（Ａ）は、入力画像としての正常画像及びそれに対応するレーダチャートの例を示しており、図３（Ｂ）は、入力画像としてのムラ欠陥画像及びそれに対応するレーダチャートの例を示している。また、図３（Ｃ）は、入力画像としてのキズ欠陥画像及びそれに対応するレーダチャートの例を示しており、図３（Ｄ）は、入力画像としての異物欠陥画像及びそれに対応するレーダチャートの例を示している。図３（Ａ）〜図３（Ｄ）に示すように、それぞれのレーダチャートの形状は異なっている。
【００１９】
このように、適切なレーダチャートが設定されると、ユーザがレーダチャートの形状を見るだけで欠陥の傾向を直感的に把握するのに役立つ。ここでは見易さのために各欠陥の事例は画像だけを見て判る入力画像を示したが、生産現場で発生する欠陥サンプルには人の目には正常サンプルとほとんど外観の変わらないようなものも存在する。そのような場合に、そのサンプルが何の異常の傾向を持つのか、といった確認や分析の際にもレーダチャートの形状を見るだけで直感的に理解することが可能となる。その他、こういった情報の可視化によって複雑な判定処理アルゴリズム自体をユーザが理解していなくとも、アルゴリズムがどのように画像を判定しているかを直感的にチェックすることができるため、検査装置の信頼度の向上や不具合の早期発見に役立つ。
【００２０】
しかし、図３に示したような抽出特徴の種類が５個程度と少なく、それぞれが各欠陥の傾向を把握するのに充分であるならば良いが、一般的に多様な異常パターンに対応することを前提とした識別器では抽出特徴の数は数十から数百まで増加することが多い。このような場合に何の工夫も無くレーダチャートを作ると、図３に示したように一瞥して異常の傾向を把握できるようなチャートとは異なるチャートが生成されることが多い。
【００２１】
例として、抽出特徴の数が３５個である場合に、単純に抽出した特徴量を抽出特徴順に時計回りに軸を配しただけのレーダチャートを描画したことによって欠陥の傾向によるグラフの形状の違いを把握しにくくなる典型例を図４に示す。図４（Ａ）〜図４（Ｄ）に示すレーダチャートは、それぞれ異なる理由で欠陥とされた検査対象の画像を入力としたときに生成されたグラフである。生産現場において１日に大量の検査対象がラインを流れていく中で人がこれら４つの形状の違いを一瞥して把握し、瞬時にそれぞれがどういった欠陥を意味しているのかを理解するのは困難である。
【００２２】
そこで、本実施形態では、グラフに対して以降で説明する処理を施すことによって直感的に形状の違いを把握しやすいグラフに変換する。図５（Ａ）〜図５（Ｄ）に、図４（Ａ）〜図４（Ｄ）を各軸が示すスコアはそのまま保存しながら、それぞれ認識しやすい形状になるように軸の順を並べ替えたレーダチャートを示す。これによると図４（Ａ）〜図４（Ｄ）においては、不明瞭であった性質がよく判るようになる。例えば、図５（Ａ）、図５（Ｂ）、図５（Ｃ）はそれぞれ異なる欠陥の種類を示す画像に対応するチャートであったこと、図５（Ｄ）は、図５（Ａ）及び図５（Ｂ）に示される２つの欠陥種が複合している可能性が高いこと等が即座に理解しやすい。
【００２３】
つまり、たとえ抽出特徴の数が増えたとしても、図５（Ａ）〜図５（Ｄ）に示すように同じ欠陥信号に反応する特徴量を表す軸をそれぞれ近くに配置するだけで、形状を見るだけで欠陥の傾向を把握できるレーダチャートを作成することができるようになる。以下では、見やすいレーダチャート生成アルゴリズムの例を詳細に説明する。
【００２４】
前述したとおり、見やすいレーダチャートを生成するためには、関連性の高い軸を近くに配置することが好ましい。決まった傾向の入力欠陥信号に対して決まった反応をする特徴量があれば、それらの軸が近くに配置されたレーダチャートでは該当の箇所がグラフの形状としていつも連動して変形する。そのため、ひとまとまりの領域に同じ意味付けをすることで数少ない部分領域を人は意識してグラフを読み取れば良いことになるためである。逆に、連動して反応する軸が散在したレーダチャートでは同じような意味を持つ特徴量の集合をまとめて認識することが人にとって困難となる。
【００２５】
ここで処理の流れに関して、図１の本実施形態における情報処理装置としての多変量データ可視化装置の構成例を示すブロック図、及び図２のフローチャートを用いて説明する。多変量データ入力部１０１にて、画像等の対象となるサンプルからｎ個の特徴量が抽出されたデータが入力される（Ｓ２０１）。そして、入力されたｎ個の特徴量を基に特徴量間距離算出部１０２にて各々の特徴量同士の類似度（距離）を計算し（Ｓ２０２）、得られた特徴量間の類似度を尺度とする特徴空間で表現された各特徴量を代表するｎ個のベクトルが獲得される。
【００２６】
このｎ個のベクトルを次元削減部１０３にて原点から等距離に点が並ぶように２次元（又は３次元）に次元削減を行う（Ｓ２０３）。このとき、データ可視化時に必要となる次元削減数やどの特徴量同士を近付けたい（又は遠ざけたい）等のパラメータ類をユーザが指定したい場合、これらの変更を結果に反映させるため、データ可視化パラメータ入力部１０４から入力することも可能である。
【００２７】
そして、求まった原点から等距離に並んだｎ個の特徴量を代表するベクトルやユーザがデータ可視化パラメータ入力部１０４により入力したパラメータを基にしてチャート形状決定部１０６にてチャートの軸の配置が決定（Ｓ２０４）される（Ｓ２０４）。この配置にならって多変量データ入力部１０１から入力されたデータのそれぞれが各軸毎に評価値をプロットされてチャートが完成し、ディスプレイ等の出力部１０５にて出力される（表示される）。
【００２８】
また、チャート形状決定部１０６にて決定したチャートの軸配置をチャート形状記憶部１０８で記憶しておく。そして、チャートの軸配置に使われなかった多変量データを多変量データ追加入力部１０７にて受付け、チャート形状記憶部１０８で保持されているチャートの軸の上にプロットすることで新たなチャートを作成して出力部１０５に出力する。
【００２９】
次に、図２に示したステップＳ２０２での具体的な特徴量間の類似度（距離）の算出例について説明する。特徴量同士の関連が高いかどうかは、例えばカルバック・ライブラー情報量で定義することができる。カルバック・ライブラー情報量によって特徴量同士の類似性を計算する場合、各特徴量が示すデータの分布をカーネル密度推定（ＫＤＥ）等のアルゴリズムによって推定し、推定された分布間の距離によって定義する。特徴量ＩＤがｉとｊの２つの特徴量ＦiとＦjのそれぞれが示すデータバラつきから推定した確率分布をＰiとＰjとすると、２つの分布間の類似度は下記式（１）で定義される。
【００３０】
【数１】
0006602013.tif 000002
【００３１】
ただし、分布間の距離は式（１）で示した以外にも、２つの確率分布が多次元正規分布であると仮定して、Ｂｈａｔｔａｃｈａｒｙｙａ距離を用いたり、相互情報量によって定義する等、他の指標で算出しても良いことは言うまでもない。また、各確率分布推定と分布間の距離算出を同時に行う密度比推定と呼ばれる方法によっても同様の結果を得ることができる。その他、分布間の距離算出に使える指標としては、特徴量毎が示すデータのスコアランキングの差の２乗和の総和を使う方法がある。この場合にもスコアが大きい程、特徴量同士が似ていないものとして扱える。
【００３２】
以上のようにして、全ｎ個の入力特徴量Ｆ＝｛Ｆi｝ni=1同士の類似度を算出する。これにより、全ｎ個の入力特徴量は同じく自身を含めてｎ個の特徴量との類似度を要素とするベクトルが獲得される。このとき得られたベクトルを式（２）で表す。
【００３３】
【数２】
0006602013.tif 000003
【００３４】
次に、ステップＳ２０３に続くステップＳ２０３で行う処理について具体的な算出法の例について説明する。ステップＳ２０３では、２次元のレーダチャートの軸としてｎ個の特徴量を表す軸は２次元へ次元削減を行うため、次元削減後に対応する軸の描画のための点の座標を式（３）で表す。式（３）ではｍ次元に次元削減する場合を示しているが、通常のレーダチャートではｍ＝２とすれば良い。また、式（３）においてＢはレーダチャート生成のための埋め込み行列であり、定義は式（４）に示したとおりである。
【００３５】
【数３】
0006602013.tif 000004
【００３６】
【数４】
0006602013.tif 000005
【００３７】
ここで、埋め込み行列Ｂを求めるためのいくつかの方法を以下に示す。求めるべき埋め込み行列ＢをＢ*と呼ぶことにし、また、近付けるべき特徴量同士の規則を記述した類似度行列Ｗによって式（５）で埋め込み行列を定義する。
【００３８】
【数５】
0006602013.tif 000006
【００３９】
ここで、類似度行列Ｗの要素であるＷi,jは、特徴量同士を近付ける規則を記述した行列であり、ｉ番目の特徴とｊ番目の特徴を近付ける場合には１に、遠ざける場合には０になるような関数として設定されれば何でも良い。つまり、類似度Ｗi,jが１の場合には式（５）で次元削減後の距離が最小化対象となり、０の場合には無視されるようにしてＢ*を求めることになる。当然のことながら、近付ける優先度を対象によって変えたい場合には、Ｗi,jとして０〜１の間の値を設定しても良い。類似度Ｗi,jの例として使える式の例を式（６）、式（７）に示す。
【００４０】
【数６】
0006602013.tif 000007
【００４１】
【数７】
0006602013.tif 000008
【００４２】
式（６）は、ｎ次元特徴空間でのＧiがＧjのｋ近傍にあるかどうかで０か１かを決定する方式であり、式（７）は、定数γによって定義される距離ベースで算出される値をセットする方式である。その他のＷi,jの例として、抽出特徴量に関する事前知識を反映させた値をセットしても良い。また、抽出特徴の性質上、ユーザの都合によりデータ可視化の際に近付けたくない特徴量が存在する場合には、ユーザによる該当特徴量ＩＤ指定の後、該特徴量同士の類似度Ｗi,jを０にセットすることでユーザ希望の結果を得ることができる。逆に、ユーザの希望で近付けたい特徴があれば同様の入力により、該当する特徴量間のＷi,jを１にセットすれば良い。
【００４３】
これらの指標に基き次元削減を行うが、式（５）にある次元削減後の原点からの距離が１になるという制約条件により、各特徴量ＩＤを代表する点が原点からの距離１の円上に並ぶ結果を得る。このときの処理前と処理後のイメージ図、及びそのとき得られるチャートの例を図６に示す。
【００４４】
図６（Ａ）は、前述した全入力特徴量をカルバック・ライブラー情報量によって類似度ベースの距離により定義された１０次元空間内に１０個の各特徴量を代表するベクトルＧiが分布している様子を表している。図６（Ｂ）は、前述した方法にて算出された次元削減用の埋め込み行列によって２次元に次元削減された各特徴量を代表するベクトルｇiが原点周りの距離１の円上に分布していることを示している。これが、レーダチャートで各軸をどの位置関係に描画するかの指標になる。
【００４５】
この円上の点ｇiが各特徴量を表す特徴量ＩＤのｉ番目の軸を設定するための位置指標となるため、原点から全１０個の点に向かって軸を設定し描画することによって、図６（Ｃ）に示したようなチャートを得ることができる。このチャートはレーダチャートのように軸を等間隔に配していないが、各軸の類似度が各軸の間の角度として表現されているため、意味的なまとまりを把握しやすくなっており、一般的なレーダチャートよりも理解しやすくなる場合がある。このように、軸の相対的な位置と角度に情報価値が無かったレーダチャートの欠点を補った、新しいチャートをＣｌｕｍｓｙ  Ｓｐｉｄｅｒ  Ｃｈａｒｔ（ＣＳＣ）と称する。
【００４６】
このチャート上に任意のデータの各軸スコアをプロットして描画したＣＳＣの例を図６（Ｄ）に示す。図６（Ｄ）は説明を簡易にするために軸の数が少ない場合のＣＳＣとして表現しているが、ＣＳＣが便利になるのはより高次元の多変量を扱うときである。高次元になっても、連動してスコアが変わる軸が近く配置されているために、局所に注目する場合は近い軸同士のスコア変動をまとめて観察すればよく、全体のグラフ形状を見ると局所のまとまりを持ったスコアの動きのバランスを俯瞰して把握することができる。
【００４７】
以上説明したように、軸同士の関連性や類似性を考慮して軸の配置（並びや間隔）を設定してチャート（グラフ）を生成することによって、全体を俯瞰した際に複数の軸同士の意味的なまとまりを直感的に把握することが可能となる。これにより、巨視的な視点で見たときのグラフ全体の形状比較がし易くかつ微視的な視点では軸の近い関連性の高い物同士の差異に注目して比較することのできる、人の認識にやさしいチャート（グラフ）を生成することが可能となる。
【００４８】
図７（Ａ）、図７（Ｂ）に高次元（１０８次元）の多変量データの比較が直感的に優れていることがわかるＣＳＣの例を示した。図７（Ａ）、図７（Ｂ）によると、一般的なレーダチャートよりも、どの軸が相関が強く連動して動くことが直感的にわかるため、グラフの注目ポイントが暗に示されており、高次元になっても形状の認識がしやすいことがわかる。図７（Ａ）と図７（Ｂ）に示したチャートを比較するときも全体的な形状の傾向は変わらないが、局所的な形状としての右上部分、左上部分の形状に違いがあるということが判り易い。軸の並びが適切な順になっていない一般的なレーダチャートであると、１０８本もの軸の上のグラフを一瞥して傾向を把握するのは困難になるのは言うまでもない。
【００４９】
また、通常のレーダチャートのように軸の間の角度が等しい表示を好むユーザのために外観だけを一般的なレーダチャートに合わせるように軸の並びはそのままにしながら軸間角度を等間隔になるようにしても良い。その場合の変換過程の概念は、図８（Ａ）〜図８（Ｄ）に示しており、図６（Ｄ）に対応するグラフを図８（Ｅ）に示している。
【００５０】
詳細な分析ではなく、より短時間で直感的にデータの傾向を把握するためには、できるだけ冗長な情報になる可能性の高い部分は取り除いてユーザに提示した方が良い場合がある。そのため、さらに簡易チャート表示モード等を選択できるようにし、特徴量同士の類似度が所定の値より高い軸同士を１つの軸にまとめて表示しても良い。その際の処理手順を図９（Ａ）〜図９（Ｄ）に示す。例えば、次元削減後の軸の間の角度が設定した閾値以下であった場合にそれらの軸をまとめて１つの軸で表し、その軸で表されるスコアは元の複数の軸のスコアの平均値にする等により算出する方法がある。または次元削減前に近傍関係にある点の平均を取る等して１つの点にまとめ、その後に次元削減をしても良い。図９（Ｅ）は、図６（Ｄ）や図８（Ｅ）に示したデータを同様に図９（Ｄ）のチャート上に表示した例であるが、図９（Ｅ）をさらに一般的なレーダチャートの形式にして見たいユーザのために全軸間の角度を等間隔になるように調整して表示しても良い。
【００５１】
なお、前述した各種チャートの表示切り替えは、例えば図１０に示すように複数の画像に対応するチャートを表示させ、チャートを見ながらユーザが好みの表示形式をボタンをクリックするだけで簡単に選択できるようにしてもよい。これによりデータによって適した分析が容易になる。
【００５２】
また、前述した説明では、すべて一般的なレーダチャートに近い表示形態を想定したため、次元削減時の次元数を２次元とした例を用いて説明したが、より複雑な分析が必要な場合等は次元削減時の次元数を３次元としても良い。削減後の次元数が３である場合、前述したチャートの３次元表現となり、よりそれぞれの特徴量の関係が直感的に表現されたグラフを作成することができる。ただし、３次元のチャートを確認するディスプレイが２次元表示の場合等は２次元のチャートよりも形状の把握が難しくなる場合がある。そのために３次元化したチャートの際は、緯度や経度によってＲＧＢ（表示色）の値を変化させながらカラーで表示するなどすれば良い。なお、本実施形態は外観検査を例として説明したが複数の特徴量を扱った場合のデータ可視化全般に利用できる方法であることは言うまでもない。
【００５３】
（第２の実施形態）
次に、本発明の第２の実施形態について説明する。第２の実施形態は、第１の実施形態で説明したチャート生成のための次元削減時において次元削減する前後において特徴量間の距離関係の保持度合いを基準にして次元削減を行い、人が見て直感的に判り易い形でのレーダチャート等を生成する方法である。以下では、第２の実施形態において、前述した第１の実施形態と異なる点について説明する。
【００５４】
次元削減する前後において特徴量間の距離関係の保持度合いを基準にして次元削減を行う場合、最小化すべきは次元削減する前後で距離関係が変わることによる誤差である。このときの誤差２乗和を例えば式（８）に示した方法で定義すると良い。
【００５５】
【数８】
0006602013.tif 000009
【００５６】
ただし、次元削減前の特徴量Ｇやｇは、第１の実施形態で示したものと同じであり、Δijは式（９）のように定義し、Ｇi、Ｇjの間の距離を表す。
【００５７】
【数９】
0006602013.tif 000010
【００５８】
また、δijは同様に次元削減後のｇi、ｇjの距離を表す。ただし、δijには２次元平面に次元削減した後に原点からの距離が等しくなるように式（１０）に示す制約を加える。
【００５９】
【数１０】
0006602013.tif 000011
【００６０】
以上のようにした後、最急降下法により解を求めることができる。これによって求まった解は、２次元平面の上で原点周辺の同心円上に配置された点として次元削減後の特徴量の代表点を得る。この点を基にして以降のチャートを獲得する手順は、第１の実施形態と同様である。なお、式（８）以外の誤差２乗和関数の例を式（１１）、式（１２）に示す。これらはどれも多次元尺度法等で利用されている関数である。
【００６１】
【数１１】
0006602013.tif 000012
【００６２】
【数１２】
0006602013.tif 000013
【００６３】
その他、式（８）、式（１１）、式（１２）に記した最小化すべき誤差２乗和関数は、次元削減前の特徴量Ｇのすべての関係性を保持するための誤差となっている。ｋ近傍までを近付けるといったパラメータを追加するとＩＳＯＭＡＰと同様な手順で次元削減を行うことができる。以上手順によって得られた次元削減後の各特徴量を表す代表点であるｇを基にしてチャートを生成する手順は、第１の実施形態と同様手順によって実現でき、適切なチャートを得ることができる。
【００６４】
（第３の実施形態）
次に、本発明の第３の実施形態について説明する。第３の実施形態は、高次元の特徴量で表現される複数のデータがあり、各々のデータがそれらの中でどのような傾向にあるかを分析したい場合に人が見て直感的に判り易い表示に平行座標プロットや棒グラフを表示する方法である。第１の実施形態又は第２の実施形態で説明した方法によって、レーダチャートの軸同士の間隔や順番を決定するアルゴリズムを用いることによって、同様に平行座標プロットや棒グラフの項目の順番と間隔を決定することができる。また、レーダチャートでは円形状になるように軸を配することができたため、平行座標プロットや棒グラフでもｎ個の項目にてグラフを生成する際、末端であるｎ個目の項目の隣に１個目の項目からのグラフをさらに描画しても良い。
【００６５】
（本発明の他の実施形態）
また、本発明は、以下の処理を実行することによっても実現される。即ち、前述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【００６６】
例えば、前述した各実施形態における情報処理装置は、図１１に示すようなコンピュータ機能１１００を有し、そのＣＰＵ１１０１により各実施形態での動作が実施される。コンピュータ機能１１００は、図１１に示すように、ＣＰＵ１１０１と、ＲＯＭ１１０２と、ＲＡＭ１１０３とを備える。また、操作部（ＣＯＮＳ）１１０９のコントローラ（ＣＯＮＳＣ）１１０５と、ＬＣＤ等の表示部としてのディスプレイ（ＤＩＳＰ）１１１０のディスプレイコントローラ（ＤＩＳＰＣ）１１０６とを備える。さらに、ハードディスク（ＨＤ）１１１１、及びフレキシブルディスク等の記憶デバイス（ＳＴＤ）１１１２のコントローラ（ＤＣＯＮＴ）１１０７と、ネットワークインタフェースカード（ＮＩＣ）１１０８とを備える。それら機能部１１０１、１１０２、１１０３、１１０５、１１０６、１１０７、１１０８は、システムバス１１０４を介して互いに通信可能に接続された構成としている。
【００６７】
ＣＰＵ１１０１は、ＲＯＭ１１０２又はＨＤ１１１１に記憶されたソフトウェア、又はＳＴＤ１１１２より供給されるソフトウェアを実行することで、システムバス１１０４に接続された各構成部を総括的に制御する。すなわち、ＣＰＵ１１０１は、前述したような動作を行うための処理プログラムを、ＲＯＭ１１０２、ＨＤ１１１１、又はＳＴＤ１１１２から読み出して実行することで、各実施形態での動作を実現するための制御を行う。ＲＡＭ１１０３は、ＣＰＵ１１０１の主メモリ又はワークエリア等として機能する。
【００６８】
ＣＯＮＳＣ１１０５は、ＣＯＮＳ１１０９からの指示入力を制御する。ＤＩＳＰＣ１１０６は、ＤＩＳＰ１１１０の表示を制御する。ＤＣＯＮＴ１１０７は、ブートプログラム、種々のアプリケーション、ユーザファイル、ネットワーク管理プログラム、及び各実施形態における動作を実現するための処理プログラム等を記憶するＨＤ１１１１及びＳＴＤ１１１２とのアクセスを制御する。ＮＩＣ１１０８はネットワーク１１１３上の他の装置と双方向にデータをやりとりする。
【００６９】
なお、前記実施形態は、何れも本発明を実施するにあたっての具体化のほんの一例を示したものに過ぎず、これらによって本発明の技術的範囲が限定的に解釈されてはならないものである。すなわち、本発明はその技術思想、又はその主要な特徴から逸脱することなく、様々な形で実施することができる。
【符号の説明】
【００７０】
１０１：多変量データ入力部  １０２：特徴量間距離算出部  １０３：次元削減部  １０４：データ可視化パラメータ入力部  １０５：出力部  １０６：チャート形状決定部  １０７：多変量データ追加入力部  １０８：チャート形状記憶部
【図１】
0006602013.tif 000014
【図２】
0006602013.tif 000015
【図３】
0006602013.tif 000016
【図４】
0006602013.tif 000017
【図５】
0006602013.tif 000018
【図６】
0006602013.tif 000019
【図７】
0006602013.tif 000020
【図８】
0006602013.tif 000021
【図９】
0006602013.tif 000022
【図１０】
0006602013.tif 000023
【図１１】
0006602013.tif 000024