(57)【要約】【課題】或る項目に係るデータの集合に含まれる一部のデータの発生時刻周辺において、当該或る項目又は他の項目のデータの時間変化の態様を容易に把握する。【解決手段】或る項目のデータの分布を示したグラフ上から入力装置によって指定された一部のデータに対応付けられた時刻を各時刻の間隔に基づいて複数のグループに分け、当該複数のグループがそれぞれ規定する時間範囲に含まれる所望の項目に係るデータを記憶装置から検索し、当該検索結果に基づいて当該所望の項目に係るデータの時間波形をグループごとに作成する。【選択図】図８
2015072515.tif 000002
【特許請求の範囲】
【請求項１】
複数の項目に係るデータがそれぞれ時刻に対応付けて記憶されている記憶装置と、
前記複数の項目のうち１つの項目に係るデータの分布を示した第１グラフが表示される表示装置と、
前記第１グラフに示された前記データの分布に含まれる一部のデータの指定が可能な入力装置と、
当該入力装置によって指定された一部のデータに対応付けられた時刻に基づいて前記記憶装置のデータを検索する時間範囲を複数決定し、当該各時間範囲に含まれるデータを前記記憶装置から検索した結果に基づいて当該データの時間変化を示す第２グラフを前記時間範囲ごとに前記表示装置に表示する処理装置とを備えることを特徴とするデータ表示システム。
【請求項２】
請求項１に記載のデータ表示システムにおいて、
前記処理装置は、前記第２グラフを表示するために、
当該入力装置によって指定された一部のデータに対応付けられた時刻を複数のグループに分ける第１処理と、
当該複数のグループのそれぞれに含まれる時刻で規定される時間範囲に含まれる所望の項目に係るデータを前記記憶装置から検索する第２処理と、
当該検索結果に基づいて前記所望の項目に係るデータの時間変化を示す第２グラフを前記複数の時間範囲ごとに前記表示装置に表示する第３処理とを実行することを特徴とするデータ表示システム。
【請求項３】
請求項１に記載のデータ表示システムにおいて、
前記処理装置は、さらに、前記複数の時間範囲ごとの前記第２グラフの形状を比較して、類似する形状を有するグラフをグループ化して前記表示装置に表示する処理を行うことを特徴とするデータ表示システム。
【請求項４】
請求項２に記載のデータ表示システムにおいて、
前記処理装置は、前記第１グラフを前記表示装置に表示する前に、前記複数の項目の中から当該第１グラフを表示する項目を選択するための第１選択画面を前記表示装置に表示する処理を行い、前記入力装置により当該第１選択画面を介して選択された項目に係るデータの分布を前記第１グラフとして表示する処理を行うことを特徴とするデータ表示システム。
【請求項５】
請求項２に記載のデータ表示システムにおいて、
前記処理装置は、
前記第２処理の前に、前記第２処理で前記記憶装置からデータが検索される項目を前記複数の項目の中から選択するための第２選択画面を表示する第４処理を行い、
前記第２処理として、前記４処理で選択された項目に係るデータのうち前記複数の時間範囲に含まれるものを前記記憶装置から検索する処理を行うことを特徴とするデータ表示システム。
【請求項６】
請求項１に記載のデータ表示システムにおいて、
前記入力装置によって、前記第１グラフに示された前記データの分布から当該データの一部のデータとして指定されるのは、他のデータから外れた値を有する外れ値であることを特徴とするデータ表示システム。
【請求項７】
請求項１に記載のデータ表示システムにおいて、
前記記憶装置は、前記複数の項目に係るデータとして、機械に搭載された複数のセンサの検出値を記憶しており、
当該複数のセンサに係るセンサ検出値は、それぞれ、その検出時刻に対応付けて記憶されていることを特徴とするデータ表示システム。
【請求項８】
請求項１に記載のデータ表示システムにおいて、
前記第１グラフは、前記１つの項目に係るデータの分布を前記複数の項目に含まれる他の項目に係るデータの分布に対応させて示した散布図、または、前記１つの項目に係るデータの分布を示すヒストグラムであることを特徴とするデータ表示システム。
【請求項９】
複数の項目に係るデータがそれぞれ時刻に対応付けて記憶されている記憶装置と、
前記複数の項目のうち１つの項目に係るデータの分布を示した第１グラフが表示される表示装置と、
前記第１グラフに示された前記データの分布に含まれる一部のデータの指定が可能な入力装置と、
当該入力装置によって指定された一部のデータに対応付けられた時刻を各時刻の間隔に基づいて複数のグループに分ける第１処理と、当該複数のグループのそれぞれに含まれる最初と最後の時刻で規定される時間範囲に含まれる所望の項目に係るデータを前記記憶装置から検索する第２処理と、当該検索結果に基づいて前記所望の項目に係るデータの時間変化を示す第２グラフを前記複数の時間範囲ごとに前記表示装置に表示する第３処理とを行う処理装置とを備えることを特徴とするデータ表示装置。
【請求項１０】
複数の項目に係るデータがそれぞれ時刻に対応付けて記憶されている記憶装置を利用して、前記複数の項目のうち１つの項目に係るデータの分布を示す第１グラフを作成する第１ステップと、
前記第１グラフに示された前記データの分布に含まれる一部のデータを指定する第２ステップと、
当該第２ステップで指定された前記一部のデータに対応付けられた時刻を各時刻の間隔に基づいて複数のグループに分ける第３ステップと、
当該複数のグループのそれぞれに含まれる最初と最後の時刻で規定される時間範囲に含まれる所望の項目に係るデータを前記記憶装置から検索する第４ステップと、
当該検索結果に基づいて前記所望の項目に係るデータの時間変化を示す第２グラフを前記複数の時間範囲ごとに作成する第５ステップとを備えることを特徴とするデータ表示方法。
【請求項１１】
請求項１０に記載のデータ表示方法において、
前記第５ステップで作成される前記複数の時間範囲ごとの前記第２グラフの形状を比較して、類似する形状を有するグラフをグループ化する第６ステップをさらに備えることを特徴とするデータ表示方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は或る項目に係るデータの時間変化を表示するデータ表示システムに関する。
【背景技術】
【０００２】
ガスエンジンやエレベータ、採掘・建築機械といった機械を常に動作させるためには、機械の保守作業が必須である。保守作業で有効な技術の１つに、機械の各部に取り付けられた複数のセンサのセンサデータ（検出値）を収集し、収集したセンサデータから機械の異常診断を行い、異常があった場合はその原因分析を行うものがある。
【０００３】
該技術における異常診断の方法として、機械のセンサデータの分布や出現頻度を散布図やヒストグラム等のグラフで表現し、そのグラフにおいて他の値から大きく外れた外れ値に基づいて機械の異常を調べるものがある。本方法では外れ値出現時に機械に異常が発生していると推定するが、実際に何が原因で外れ値が発生しているのかを調査して異常の原因を特定する必要がある。外れ値の発生原因の特定のためには、外れ値出現時のセンサデータの時間変化を分析する必要がある。しかしながらヒストグラムや散布図から外れ値の出現時の時間波形を得るためには、まず所望の外れ値の発生時刻を特定し、その発生時刻の時間波形を読み出す作業が必要である。この作業を外れ値の調査の度に行っていたのでは、調査の効率が著しく低下するおそれがある。
【０００４】
この種の問題の解決を試みたデータ表示装置として、例えば特開平７−２８２２７７号公報がある。該文献では所定の日付（例えば、５月１６日）における機械の不良Ａ〜Ｅの発生数と不良総数の中で各不良Ａ〜Ｅの占める割合をパレート図で画面表示し、所定期間の発生数の時系列をグラフ表示したい不良Ａ〜Ｅを分析者に選ばせる。そして、分析者の選択という単純な操作のみに基づいて、当該選択した不良の発生数の所定期間における時系列を棒グラフ（トレンド図）で表示する。当該棒グラフにおける横軸は日付を示し、当該日付としては最初にパレート図で利用した所定の日を含む所定の日数（例えば、５月１日から１６日）が設定されている。これにより、時間軸における断面的な特性を示すパレート図と、時間軸に沿う特性を示すトレンド図が有機的に結合し、個々のデータの時系列変化の容易な把握が可能になっている。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開平７−２８２２７７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
上記文献に係る発明は、所定の日に係る複数の不良の発生数を示す棒グラフ（パレート図）から所望の不良を選択し、当該選択した不良について当該所定の日を含む数日間分における不良発生数を示す棒グラフを表示している。そのため、日付ごとの不良発生数の変化の把握は容易になるものの、１日よりも短い時間（例えば、数秒から数時間程度の短期間）で所定の傾向を有する外れ値が出現した場合には、日付ごとの不良発生数を示したグラフを利用しても分からないことが多く、原因分析が困難になることがある。例えば、不良によっては、短期間で外れ値がまとまって発生することがあり、さらに当該外れ値の集合と同様の傾向を有するものが間隔を空けて複数回発生することもあるため、この種の不良の分析には上記文献の技術は不向きである。
【０００７】
また、短時間で発生する外れ値に対応するために、上記技術における不良発生数の時系列のグラフの横軸（時間軸）の単位を１日未満の値（例えば、分、時間）に設定して上記問題に対応することが考えられるが、そのままでは当該時系列のグラフから外れ値が発生した時刻を特定する作業や、当該部分の抽出・比較作業が必要となり、なおも不良分析作業が滞るおそれがある。
【０００８】
さらに、上記で例示した外れ値だけに限らず、時刻に対応付けられた或る項目のデータの集合に含まれる一部のデータ（部分集合）の時間変化の態様を容易に把握することが望まれる。
【０００９】
本発明の目的は、或る項目に係るデータの集合に含まれる一部のデータ（部分集合）が数秒から数時間といった短時間に発生する場合に、当該一部のデータの発生時刻周辺における当該或る項目または他の項目に係るデータの時間変化の態様を容易に把握できるデータ表示システムを提供することにある。
【課題を解決するための手段】
【００１０】
上記目的を達成するために、本発明は、複数の項目に係るデータがそれぞれ時刻に対応付けて記憶されている記憶装置と、前記複数の項目のうち１つの項目に係るデータの分布を示した第１グラフが表示される表示装置と、前記第１グラフに示された前記データの分布に含まれる一部のデータの指定が可能な入力装置と、当該入力装置によって指定された一部のデータに対応付けられた時刻を各時刻の間隔に基づいて複数のグループに分ける第１処理と、当該複数のグループのそれぞれに含まれる最初と最後の時刻で規定される時間範囲に含まれる所望の項目に係るデータを前記記憶装置から検索する第２処理と、当該検索結果に基づいて前記所望の項目に係るデータの時間変化を示す第２グラフを前記複数の時間範囲ごとに前記表示装置に表示する第３処理とを行う処理装置とを備える。
【発明の効果】
【００１１】
本発明によれば、或る項目に係るデータの集合の一部（部分集合）の発生時刻周辺において時間変化を知りたい項目を指定することで、当該指定項目に係るデータの当該発生時刻周辺の時間変化の態様を容易に把握することができる。これにより、例えば、数秒〜数時間といった短期間に外れ値が発生している場合でも、その発生時刻周辺の短期間の時間波形グラフを表示し、外れ値発生時の波形から外れ値の主な原因を分析できる。
【図面の簡単な説明】
【００１２】
【図１】本発明の目的の説明図。
【図２】本発明の課題の説明図。
【図３】本発明の原理の説明図。
【図４】本発明の実施の形態に係るデータ表示システムの全体構成図。
【図５】センサデータベース４１０に格納されたテーブルＴ１の構造を示す図。
【図６】センサデータベース４１０に格納された他のテーブルＴ７の構造を示す図。
【図７】本発明の実施の形態に係るデータ表示システムの機能ブロック図。
【図８】図７に示した各部位が行う処理の流れを説明したシーケンス図。
【図９】図８のＳ５０５で表示装置４４０に表示される第１選択画面を示す図。
【図１０】図８のＳ５１０でテーブルＴ１から記憶装置４５０にロードされたセンサデータが格納されたテーブルＴ２の構造を示す図。
【図１１】図８のＳ５１５で作成される散布図。
【図１２】テーブルＴ３の構造を示す図。
【図１３】テーブルＴ４の構造を示す図。
【図１４】図８のＳ５３５で行われる内部処理を説明したフローチャート。
【図１５】図８のＳ５３８で表示装置４４０に表示される第２選択画面を示す図。
【図１６】テーブルＴ５の構造を示す図。
【図１７】テーブルＴ６の構造を示す図。
【図１８】図８のＳ５４５で行われる内部処理Ｓ５４５＿ＳＵＢのフローチャート。
【図１９】図１８に示した内部処理で呼び出されるサブルーチンＳ５４５＿ＳＵＢ２のフローチャート。
【図２０】図８のＳ５５０で作成されＳ５５５で表示される画面の一例を示す図。
【図２１】図８のＳ５５５で表示される画面に含まれる時間波形の他の例を示す図。
【発明を実施するための形態】
【００１３】
まず、実施の形態を説明する前に、本発明の基本概念について説明する。図１は、エンジンを搭載した機械において、圧力センサおよび回転数センサによって所定期間内に検出されたエンジン圧力とエンジン回転数の関係を示した散布図の一例であり、当該散布図にはセンサデータの部分集合として外れ値の集合１００が存在する。本方法では外れ値の出現時に機械に異常が発生している可能性があることを前提としており、外れ値の発生原因の特定のために、図１中に破線の領域１１０で示したように、集合１００に含まれる外れ値が出現した時刻についてのセンサデータ（エンジン回転数）を当該センサデータの時間変化を示すグラフ（時間波形）に基づいて分析する。
【００１４】
特開平７−２８２２７７号公報に開示された技術では、所定の日に係る複数の不良の発生数を示す棒グラフ（パレート図）から所望の不良を選択し、当該選択した不良について当該所定の日を含む数日間分における不良発生数を示す棒グラフを表示している。しかし、例えば図２に示すように、外れ値の集合２００が破線の領域２１０のように数秒〜数分程度の短期間しか継続していない場合には、日付ごとの不良発生数を示したグラフを利用しても原因が分かりづらいまたは分からないことがある。
【００１５】
（１）そこで、本発明に係るデータ表示システムでは、複数の項目（例えば、エンジン圧力と、エンジン回転数）に係るデータがそれぞれ時刻に対応付けて記憶されている記憶装置（例えば、ハードディスクドライブ等の磁気記憶装置や、フラッシュメモリ等の半導体メモリ）と、前記複数の項目のうち１つの項目に係るデータの分布を示した第１グラフ（例えば、散布図やヒストグラム）が表示される表示装置（例えば、モニタ）と、前記第１グラフに示された前記データの分布に含まれる一部のデータの指定が可能な入力装置（例えば、マウス、キーボード、タッチパネル）と、当該入力装置によって指定された一部のデータに対応付けられた時刻を各時刻の間隔に基づいて複数のグループに分ける第１処理と、当該複数のグループのそれぞれに含まれる最初と最後の時刻で規定される時間範囲に含まれる所望の項目（例えば、エンジン回転数）に係るデータを前記記憶装置から検索する第２処理と、当該検索結果に基づいて前記所望の項目に係るデータの時間変化を示す第２グラフ（例えば、エンジン回転数の時間波形）を前記複数の時間範囲ごとに前記表示装置に表示する第３処理とを行う処理装置（例えば、ＣＰＵ）とを備えた。
【００１６】
これにより、前記入力装置によって指定された一部のデータに対応付けられた時刻が、各時刻の間隔に基づいて自動的に複数のグループに分けられる（例えば、間隔が所定の閾値未満の時刻同士が同じグループに分類される）。そして、各グループに含まれる最初と最後の時刻によって規定される時間範囲について、所望の項目に係るデータの時間変化を示すグラフ（第２グラフ）が作成・表示されるので、前記入力装置を利用して第１グラフ上で指定した一部のデータが出現した時刻に係る各種データの時間変化の態様や傾向を前記表示装置上で容易に把握できる。したがって、例えば、数秒〜数時間といった短期間に外れ値が発生している場合でも、その発生時刻周辺の短期間について所望の項目のデータの時間波形が表示されるので、当該時間波形に基づいて外れ値の主な原因を分析できる。
【００１７】
さらに、本発明の作用および効果を図３に示した例を利用して説明すると次のようになる。図３の例では、回転数センサおよび圧力センサの検出値であるエンジン回転数およびエンジン圧力がその検出時刻と対応づけて記憶装置に記憶されており、エンジン回転数のデータの分布をエンジン圧力のデータの分布に対応させた散布図（図３中の左端のグラフ参照）が表示装置に表示されている。分析者（オペレータ）が、当該散布図を読み取って外れ値の集合３００をマウス等の入力装置を介して選択すると、処理装置は、当該外れ値に対応づけられた時刻を時系列順にソートし、各時刻の間隔に基づいて当該時刻を複数のグループに分ける（図３中の右上のグラフ参照）。そして、処理装置は、各グループに含まれる最初と最後の時刻に基づいて各グループの時間範囲を定義し、その各時間範囲に含まれるエンジン回転数を記憶装置から検索する。さらに、処理装置は、当該検索結果に基づいてエンジン回転数の時間変化を前記時間範囲ごと（前記グループごと）に作成する。図３の例では、時系列順にソートされた時刻は間隔が近いもの同士が集約されて３つのグループに分けられており、各グループの時間範囲に含まれるエンジン回転数データに基づいて３つの時間波形３１０，３２０，３３０が作成されている。これにより分析者は、第１グラフ上でエンジン回転数データの外れ値の集合を選択するだけで、当該選択した外れ値の検出時刻に基づいて自動的にグループ化された時間波形を取得できるので、外れ値が発生した時刻に機械にどのような変化が発生しているかを容易に把握できる。
【００１８】
（２）また、本発明は、さらに、前記グループごとの前記第２グラフの形状を比較して、類似する形状を有するグラフをグループ化して前記表示装置に表示する処理を前記処理装置により行うことが好ましい。これにより、複数の第２グラフのうち類似した形状を有するものがグループ化されて分析者に表示されるので、複数の第２グラフのうち繰り返して出現するものを容易に特定することができ、当該表示画面は、第１グラフ上で選択したデータの発生原因の推察に役立つ。
【００１９】
この特徴に関して、図３の例では、エンジン圧力とエンジン回転数の外れ値の集合に基づいて得たエンジン回転数の３つの時間波形３１０、３２０，３３０のうち、２つの波形３２０，３３０が類似している事を示すため、当該２つの波形の背景と、残りの１つの波形の背景を異ならせて表示している。すなわち、波形３１０に係る表示波形３１５の背景は白で、波形３２０，３３０に係る表示波形３２５，３３５の背景にはハッチングが施されている。例えば、この場合、分析者は、波形３１５が１件に対し、類似する波形３２５と波形３３５が２件発生しているので、外れ値の主な原因は波形３２５と波形３３５が関係していると考えることができる。また、分析者は波形３２５，３３５とその発生時間帯から、これらは機械の起動直後に発生しており、通常稼働より低いエンジン圧力・回転数になっていることが外れ値の発生原因であると推測することができる。
【００２０】
（３）また、前記第１グラフを前記表示装置に表示する前に、前記複数の項目の中から当該第１グラフを表示する項目を選択するための画面（第１選択画面）を前記表示装置に表示する処理を行い、前記入力装置により当該画面（第１選択画面）を介して選択された項目に係るデータの分布を前記第１グラフとして表示する処理を前記処理装置によって行っても良い。これにより分析者の所望の項目に係る第１グラフを表示装置に表示することができる。
【００２１】
また、その際、第１グラフとして表示するグラフの種類の選択と、当該選択した種類のグラフ表示に必要な指標の設定などを第１選択画面上で実施可能に設定しても良い。例えば、第１グラフとして表示するグラフの種類としては散布図やヒストグラムがあり、散布図の場合にはグラフの縦軸と横軸に設定する項目を第１選択画面上で選択可能にする必要がある。
【００２２】
（４）また、前記第２処理の前に、前記第２処理で前記記憶装置からデータが検索される項目を前記複数の項目の中から選択するための画面（第２選択画面）を表示する処理（第４処理）を行い、前記第２処理として、当該第４処理で選択された項目に係るデータのうち前記複数の時間範囲に含まれるものを前記記憶装置から検索する処理を前記処理装置によって行っても良い。
【００２３】
これにより、第１グラフで表示する項目と、第２グラフで表示する項目を異ならせることができる。例えば、機械の同じ部分について異なる状態を計測している複数のセンサが存在している場合には、外れ値の発生した時刻周辺における当該複数のセンサの検出値の時間変化をそれぞれ調べることで、１つのセンサの検出値だけでは分からないことが明らかになることがある。具体例としては、エンジン温度センサの外れ値についての時間波形（第２グラフ）を見ても当該外れ値の原因が分からない場合には、次に、エンジン温度と関連の深いエンジン圧力を検出する圧力センサの検出値の時間波形を表示することで、外れ値の原因をさらに探すことができる。なお、第２グラフで表示する“項目”は、その都度選択画面（第２表示画面）を表示して前記複数の項目の中から１つを選択するように構成しても良いし、予め設定しておいても良い。
【００２４】
以下、本発明の実施の形態について図面を用いて説明する。図４は本発明の実施の形態に係るデータ表示システムの全体構成図である。この図に示すデータ表示システムは、センサデータベース４１０と、処理装置４４５と、入力装置４２５と、記憶装置４５０と、表示装置４４０を備えている。
【００２５】
センサデータベース４１０は、鉄道や建設機械等の機械に搭載した各種センサにより計測したセンサデータ（例えば、エンジン圧力や回転数）を格納したデータベースである。
【００２６】
図５はセンサデータベース４１０に格納されたテーブルＴ１の構造を示す図である。図５に示すテーブルＴ１は、複数のセンサの検出値（センサデータ）を格納しており、複数のセンサデータ８１５，８２０，８２５と、その計測時刻８１０とを対応づけて格納しており、任意の時刻範囲のセンサデータが検索可能に構成されている。
【００２７】
図６はセンサデータベース４１０に格納された他のテーブルＴ７の構造を示す図である。テーブルＴ７は、互いに関連の深いセンサの組の情報を格納するテーブルである。関連の深いセンサとは、例えばエンジンの温度センサや圧力センサといった同じ部品の状態を計測しているセンサの組や、気温とエンジン温度のように連動する値を計測するセンサの組を示す。
【００２８】
例えばＴ７の１３２０の１行目のエンジン圧力と１３３０の排気温度は同じエンジンの状態を示すセンサであり関連が深い。関連の深いセンサの波形を調べることで、１つのセンサだけではわからない外れ値の原因を探ることができる。例えば温度センサの時間波形を見ても外れ値の原因が分からない場合、次に温度センサと関連の深い圧力センサの時間波形を見ることで、外れ値の原因をさらに探すことができる。関連の深いセンサは機械の設計者からヒアリングや、設計書の情報を元に決定する。
【００２９】
処理装置４４５は、ＣＰＵ等で構成されており、本実施の形態に係る各種処理に関係する演算を実行する装置である。表示装置４４０は、液晶ディスプレイなどで構成されており、後述する発生頻度情報作成部４１５や時間波形データ作成部４２０で作成された散布図（第１グラフ）や時間波形（第２グラフ）が表示される。記憶装置４５０は、処理装置４４５が各種処理を実行するためのプログラムをはじめ各種データを一次的または継続的に記憶するための装置であり、例えば、ＲＯＭ、ＲＡＭおよびフラッシュメモリ等の半導体メモリや、ハードディスクドライブ等の磁気記憶装置によって構成されている。記憶装置４５０には、例えば、後述するテーブルＴ２〜Ｔ６が格納される。入力装置４２５は、マウスやキーボード、タッチパネル等の装置で構成され、表示装置４４０に表示される散布図（第１グラフ）や表からデータやセンサをユーザが指定できるようになっている。
【００３０】
なお、処理装置４４５、表示装置４４０、記憶装置４５０、入力装置４２５と、データベース４１０とは、同じコンピュータ（計算機）に搭載しても、異なるコンピュータに搭載しても良い。前者の場合には、全ての装置４４５，４４０，４５０，４２５，４１０または当該装置の機能を搭載したコンピュータ（データ表示装置）が本発明の実施の形態となる。一方、後者の場合には、データベース４１０が搭載されたコンピュータと、他の装置が搭載されたコンピュータとをＬＡＮやインターネットなどのネットワーク４１６を介して接続し、データの送受信が可能に構成する必要がある。さらに、後述するデータ表示に係る主な処理をネットワークに接続されたサーバー側で実行し、当該サーバーへの演算指示や当該演算結果の取得をクライアントコンピュータで実行するシステム構成としても良い。すなわち、各装置４４５，４４０，４５０，４２５，４１０をどこに設置しても、本発明はその効果を発揮することができる。
【００３１】
図７は本発明の実施の形態に係るデータ表示システムの機能ブロック図である。この図に示すように本実施の形態に係るデータ表示システムは、発生頻度情報作成部４１５と、時間波形データ作成部４２０と、データ集約部４３０と、類似波形検索部４３５として機能する。
【００３２】
発生頻度情報作成部４１５は、指定した時間範囲・種類のセンサデータをセンサデータベース４１０から読み出して散布図（第１グラフ）を作成する処理を行う部分である。時間波形データ作成部４２０は、指定した時間範囲・種類のセンサデータをセンサデータベース４１０から読み出して時間波形（第２グラフ）を作成する部分である。データ集約部４３０は、テーブルＴ３の飛び飛びの時刻情報９１３〜９１８をグループ化（集約）して、その結果の時間範囲をテーブルＴ４に格納する部分である。データ集約部４３０の詳細な処理内容については図１４を用いて後述する。類似波形検索部４３５は、テーブルＴ５からセンサデータの類似した時間波形をグルーピングし、その結果をテーブルＴ６に格納する処理を行う部分である。類似波形検索部４３５の詳細な処理内容については図１８および図１９を用いて後述する。
【００３３】
次に、図８のシーケンス図を用いて、本発明で行われる処理の流れについて説明する。図８は、図７に示した各部位が行う全体的な処理の流れを説明したシーケンス図であり、各部位で行われる内部処理と、各部位間のデータ入出力を時系列順に説明している。また、図８は、本実施の形態で行われる処理の全体的な流れをＳ５０５〜Ｓ５５５のステップで説明している。各ステップで行う処理を以下で説明する。
【００３４】
Ｓ５０５では、散布図（第１グラフ）を作成するセンサデータをユーザ（分析者）に選択させる画面（第１選択画面）が表示装置４４０に表示される。図９はＳ５０５で表示装置４４０に表示される第１選択画面を示す図である。この図に示す第１選択画面は、散布図の横軸として表示するセンサを選択するためのセンサ選択部６０５と、散布図の縦軸として表示するセンサを表示するためのセンサ選択部６１０と、散布図に表示するセンサデータの計測時間範囲のうち開始時刻を決定するための開始時間設定部６１５と、当該計測時間範囲のうち終了時刻を決定するための終了時間設定部６２０とを備えている。
【００３５】
センサ選択部６０５，６１０によるセンサの選択方法としては、センサデータベース４１０にセンサデータが格納されたセンサ名の一覧が表示されるプルダウンメニューを用意しておき、当該プルダウンメニュー内から選択するものがある。また、センサ選択部６０５，６１０にキーボードでセンサ名を直接入力することで行っても良い。また、開始時間設定部６１５と終了時間設定部６２０による時刻の設定方法としては、キーボードで時刻を直接入力するものがある。センサ名と開始・終了時刻の入力を完了したら、ユーザが散布図表示のボタン６２５を押すことで散布図を作成するＳ５１０の処理に移る。
【００３６】
Ｓ５１０では、データベース４１０のテーブルＴ１（図５）に係るセンサデータのうち、Ｓ５０５でユーザが指定したセンサと時間範囲のセンサデータが記憶装置４５０にロードされる。図１０は、Ｓ５１０でテーブルＴ１から記憶装置４５０にロードされたセンサデータが格納されたテーブルＴ２の構造を示す図である。この図に示すように、記憶装置４５０にロードされたデータは、Ｓ５０５で指定した２つのセンサ（エンジン圧力およびエンジン回転数）８３５，８４０と、これに対応するセンサ計測時刻８３０とがセットでテーブルＴ２に格納されている。  
【００３７】
Ｓ５１５では、発生頻度情報作成部４１５は、テーブルＴ２から散布図を作成する。図１１はＳ５１５で作成される散布図を示している。Ｓ５１５では、具体的には、テーブルＴ２の１行目から最終行までの２つのセンサ８３５，８４０の値を図１１の散布図の上に描画する。図１１の例ではエンジン圧力８３５の値を縦軸に、エンジン回転数８４０を横軸に描画して散布図を作成している。
【００３８】
Ｓ５２０では、Ｓ５１５で作成した散布図を表示装置４４０に表示し、これをユーザに提示する。本処理はＳ５１５で発生頻度情報作成部４１５が液晶ディスプレイなどに散布図を描画することで自動的に行われる。
【００３９】
Ｓ５２５では、入力装置４２５の操作によって画面上を移動可能なカーソル７０５をＳ５２０の散布図上に表示し、当該カーソル７０５でユーザに複数のデータ点を選択させる処理を行う。ユーザは、例えば点線７２０のような軌跡に沿ってカーソル７０５を動かすことで当該軌跡に囲まれた複数のデータ点を一度に選択できる。選択した複数のデータ点は、散布図上の外れ値など、発生原因を調べたいとユーザが考えるものである。なお、散布図上のデータ点の他の選択方法としては、マウスをドラッグした方向を対角線に有する矩形を散布図上に描画し、当該矩形の内部に存在する複数のデータ点を一度に選択するものがある。また、選択したい全てのデータ点をマウスでクリックする方法など、公知の選択方法が利用可能である。
【００４０】
Ｓ５３０では、Ｓ５２５で選択した複数のデータ点に対応付けられた時刻（センサ計測時刻）をテーブルＴ２（図１０）から検索し、その結果をテーブルＴ３にまとめ、当該テーブルＴ３を記憶装置４５０に格納する処理を行う。図１２はテーブルＴ３の構造を示す図である。この図に示すように、テーブルＴ３は時刻データのみから構成されている。テーブルＴ３を作成する際には、図１１の散布図上で選択された複数のデータ点に係るエンジン圧力とエンジン回転数の値を検索キーとして、テーブルＴ２（図１０）の時刻８３０を検索する。その検索結果の時刻一覧をテーブルＴ３の時刻９１０に時系列順に格納する。
【００４１】
Ｓ５３５では、後にＳ５３８で選択されるセンサデータの時間波形を作成する準備として、テーブルＴ３に格納された複数の時刻データを時刻の間隔の大きさに応じてグループ分けする。各グループに含まれる時刻データには、最も早い時刻を示すものと、最も遅い時刻を示すものが含まれることになるが、当該２つのデータに係る時刻がＳ５５５で表示される各時間波形の時間範囲を決めることになる。
【００４２】
本発明の特徴の１つである、形状が類似した時間波形（類似波形）のグルーピングには、或るセンサデータの時間波形の中から、分析対象とする一部の時間波形を抽出することが必須である。当該一部の時間波形を抽出するには、当該一部の時間波形の時間範囲を決める必要がある。そのために、本実施の形態では、テーブルＴ３に格納された時刻データのうち、時刻が近いもの同士をグループ化（集約）することで当該一部の時間波形の時間範囲を決定している。例えば、テーブルＴ３では、10:20:43〜10:23:05の時刻データを第１グループ９１３としてグループ化（集約）している。ここでグループ化した各時間範囲の開始時刻と終了時刻はテーブルＴ４に格納する。
【００４３】
図１３はテーブルＴ４の構造を示す図である。この図に示すように、各グループの開始時刻９２０と終了時刻９３０が組となってテーブルＴ４に格納される。テーブルＴ４の第１行９２３に格納されたデータは、テーブルＴ３の第１グループ９１３の時間範囲を示している。テーブルＴ３の時刻データから第１グループ９１３を決定したように、データ値（時刻）が近いもの同士をグループ化する場合には既知の技術を利用する。既知の技術としては、例えば、データマイニング技術のクラスタリングを用いることができる。本実施の形態ではより単純な例を利用しており、具体的には、時系列順にソートされた時刻データにおいて隣接する２つの時刻データの時刻間隔が閾値以下の場合に、当該２つの時刻データを同じグループとする方法を利用している。次にこの方法について図１４を用いて説明する。なお、時刻間隔の閾値は本システムの設計時に決定しておくものとし、本実施の形態では単に閾値と呼ぶことにする。
【００４４】
図１４は、図８中のＳ５３５で行われる内部処理を説明したフローチャートである。換言すれば、テーブルＴ３の時刻情報をグループ化し、その結果をテーブルＴ４に格納するフローとなっている。まず、Ｓ１４１０では、カウンタ変数ｎの記憶領域を記憶装置４４０内に用意する。このカウンタ変数ｎはテーブルＴ３のデータの行数を示す変数であり、ｎ＝１から開始して、テーブルＴ３の最終行に係る行数値に達するまでｎは１ずつ増加する。
【００４５】
Ｓ１４２０では、テーブルＴ３からｎ行目の時刻を読み出し、テーブルＴ４の開始時刻９２０に格納する。具体例を言えば、ｎ＝１の場合は、テーブルＴ３の第１行に係るデータ（2013-03-03 10:20:43）を、テーブルＴ４の第１行９２３の開始時刻の列９２０に格納することになる。
【００４６】
Ｓ１４３０では、テーブルＴ３のｎ行目とｎ＋１行目の時間間隔ΔＴを計算する。例えばｎ＝１の場合はＴ３の第１行目と２行目が「2013-03-03 10:20:43」と「2013-03-03 10:22:05」になるので両者の時間間隔ΔＴは２２秒となる。
【００４７】
Ｓ１４４０では、Ｓ１４３０で計算した時間間隔ΔＴがシステム設計時に決めた閾値以上かどうか判定する。時間間隔ΔＴが閾値以上ならば、テーブルＴ３のｎ行目とｎ＋１行目の時刻が２つの時間波形の境界であると判断し、Ｓ１４５０に進む。反対に閾値より小さいなら１つの時間波形に含まれる時刻であるとみなし、グループ化の処理を継続するためにＳ１４８０に進む。例えば、閾値が１０分の場合において、上記のように時間間隔ΔＴが２２秒のときには、時間間隔ΔＴは閾値より小さいためＳ１４８０に進むことになる。
【００４８】
Ｓ１４８０では、ｎ行目の時刻を、ｎ−１行目の時刻と同じグループとし、次のｎ＋１行目の時刻を参照するためにカウンタ変数ｎをｎ＋１に更新してからＳ１４３０に戻る。
【００４９】
Ｓ１４５０では、Ｓ１４４０でｎ行目の時刻が波形と波形の境界であると判断したので、当該ｎ行目の時刻を波形の終了時刻としてテーブルＴ４の終了時刻の列９３０に格納する。例えばテーブルＴ３の第１グループ９１３の中では、次の時刻まで１日以上の時間間隔がある時刻「2013-03-03 10:23:05」が終了時刻となるので、テーブルＴ４の第１行９２３の列９３０に該時刻を格納する。
【００５０】
Ｓ１４６０では、本サブルーチンの終了条件を満たすかどうかの判断処理を行う。すなわち、これまでのＳ１４１０〜Ｓ１４５０の処理でテーブルＴ３の全データを参照したかを判断する。全データを参照していれば本ルーチンを完了する。一方、参照していなければＳ１４７０に移行する。
【００５１】
Ｓ１４７０では、次の時間波形を見つけるためにカウンタ変数ｎを１だけ増加してＳ１４２０に戻る。以上でＳ５３５のサブルーチンが完了し、その後、図８のＳ５３７に移行する。
【００５２】
Ｓ５３７では、Ｓ５５５で時間波形（第２グラフ）を表示するセンサの候補を決定する処理を行う。本実施の形態では、Ｓ５０５で散布図（図１１）の作成時に選択した２つのセンサ（散布図の縦軸と横軸に設定されるセンサ）と関連のあるセンサを、センサデータベース４１０に格納されたテーブルＴ７（図６）から検索し、その検索結果に係るセンサをセンサ候補としている。具体的には、図１１の散布図の縦軸と横軸であるエンジン圧力とエンジン回転数を検索キーにして、エンジン圧力またはエンジン回転数が含まれる行をテーブルＴ７の第１列から検索し、該当する行の第２列に格納されているセンサ名を関連センサとして取得する。図６に示したテーブルＴ７からは、エンジン圧力とエンジン回転数に関連するセンサとして、排気温度、排気圧、冷却水の３つが取得される。そして、当該３つのセンサに加え、検索キーとして利用した２つのセンサ（エンジン圧力とエンジン回転数）を加えた５つのセンサを、Ｓ５５５で時間波形を表示するセンサの候補とする。
【００５３】
Ｓ５３８では、Ｓ５３７で決定したセンサ候補の中から、時間波形（第２グラフ）を表示するものを１つだけユーザに選択させる画面（第２選択画面）が表示装置４４０に表示される。図１５はＳ５３８で表示装置４４０に表示される第２選択画面を示す図である。この図に示すように第２選択画面には、Ｓ５３７で取得したセンサの名前が格納された列１２３０を有するテーブルが表示される。当該テーブルの第１列にはチェックボックス１２２０が設けられており、ユーザは時間波形の表示を希望するセンサに係るチェックボックスにチェックを入れられるようになっている。チェックボックス１２２０を介してセンサの選択が完了したら、表示ボタン１２１０がアクティブになり、当該表示ボタン１２１０がユーザにより押下されることでＳ５４０に移行する。ここでは、図１５に示すようにＳ５３８でエンジン圧力が選択されたものとして説明を続ける。
【００５４】
Ｓ５４０では、Ｓ５３８で選択したセンサデータをデータベース４１０からロードして新たなテーブルＴ５を作成する処理が行われる。図１６はテーブルＴ５の構造を示す図である。本実施の形態では、Ｓ５３８で「エンジン圧力」が選択されているので、エンジン圧力をセンサデータベース４１０のテーブルＴ１（図５）からロードする。ロードするデータはテーブルＴ４（図１３）に基づいて決定され、具体的には、テーブルＴ４に規定された時間範囲（各グループの開始時刻９２０から終了時刻９３０）に含まれるセンサデータがロードされる。
【００５５】
本実施の形態では、図５に示したように、１秒ごとに検出したエンジン圧力がテーブルＴ１に格納されているため、テーブルＴ４で規定された各時間範囲に含まれる１秒ごとのエンジン圧力データがテーブルＴ５に格納される。例えば、まず、テーブルＴ４の第１行９２３に係る時間範囲（2013-03-03 10:20:43〜2013-03-03 10:23:05）のエンジン圧力データをテーブルＴ１からロードする。ロードしたセンサデータ１０１３は時系列順にテーブルＴ５に格納される。
【００５６】
テーブルＴ５の第３列１０４０の波形ＩＤに格納されるデータ（整数）は、Ｓ５３５でのグループのＩＤと一致しており、テーブルＴ５への格納の際に各センサデータに付与される。テーブルＴ５における波形ＩＤ１０４０は１から順に１、２、３、・・・と採番する。つまり、テーブルＴ４の第１行９２３に係る時間範囲のセンサデータ１０１３の波形ＩＤを１としてテーブルＴ５に格納し、テーブルＴ４の第２行９２６に係る時間範囲のセンサデータ１０１６の波形ＩＤを２としてテーブルＴ５に格納し、第３行９２７に係る時間範囲のセンサデータ１０１８の波形ＩＤを３としてテーブルＴ５に格納する。これ以降は波形ＩＤの値を１ずつ増加しながら同様にテーブルＴ５にセンサデータを格納する。波形ＩＤは、次のＳ５４５で類似した波形のグルーピングに使用される。
【００５７】
Ｓ５４５では、テーブルＴ５で同じ波形ＩＤが付されたセンサデータ１０１３，１０１６，１０１８、・・・によって作成される時間波形（第２グラフ）を、その形状に応じてグルーピングする処理が行われる。時間波形のグルーピングの結果はテーブルＴ６に格納される。図１７はテーブルＴ６の構造を示す図である。テーブルＴ６の第１列１０６０には波形ＩＤが格納され、第２列１０８０には波形のグルーピング結果である波形グループのＩＤが格納される。図１７の例では、波形ＩＤ＝１の波形は波形グループＩＤ１０００に、波形ＩＤ２，３は波形グループＩＤ２０００に属している。
【００５８】
次にＳ５４５で行われる波形のグルーピングの具体的な手順をＳ５４５のサブルーチンとして図１８，１９を用いて以下で説明する。図１８は図８中のＳ５４５で行われる内部処理Ｓ５４５＿ＳＵＢのフローチャートであり、図１９は図１８に示した内部処理で呼び出されるサブルーチンＳ５４５＿ＳＵＢ２（類似した波形をグルーピングする処理）のフローチャートである。
【００５９】
図１８のＳ１５１０ではテーブルＴ５の波形ＩＤを指定するため、波形ＩＤを示す変数Ｘの領域を記憶装置４５０に確保する。波形ＩＤ変数Ｘの初期値は波形ＩＤの最小値と同じく１とする。
【００６０】
Ｓ１５２０では波形グループＩＤを示す変数Ｇの領域を記憶装置４５０に確保する。変数Ｇは形状の類似した波形が属する波形グループＩＤを示す。変数Ｇの初期値は本実施の形態では１０００としておく。
【００６１】
Ｓ１５３０ではテーブルＴ５（図１６）から変数Ｘの値を検索キーにして第３列１０４０の波形ＩＤを検索して、波形ＩＤ＝Ｘのセンサデータを見つける。例えば、変数Ｘ＝１なら、テーブルＴ５で波形ＩＤ＝１のセンサデータ１０１３が見つけるべきセンサデータになる。
【００６２】
Ｓ１５４０ではＳ１５３０でレコードが見つかったか否かを判定する。Ｓ１５３０でテーブルＴ５からセンサデータが１行以上見つかればＳ１５５０に進む。見つからなければ全ての波形ＩＤを検索完了したと判断して本サブルーチンを完了する。
【００６３】
Ｓ１５５０では、Ｓ１５３０で見つけた波形ＩＤ＝Ｘの時間波形の形状と類似する形状を有する時間波形（類似波形）の有無をテーブルＴ５から検索し、類似波形が存在する場合には両者をグルーピングする処理を行う。例えば、変数Ｘ＝１（すなわち、波形ＩＤ＝１）ならセンサデータ１０１３の時間波形と、他の波形ＩＤに係るセンサデータ（例えば、センサデータ１０１６，１０１８）の時間波形とを比較し、類似波形が存在する場合には、テーブルＴ６において、センサデータ１０１３の時間波形と類似する時間波形の波形ＩＤに対して同じ波形グループＩＤである１０００（変数Ｇの値）を格納する。ここでは、当該グルーピング処理に当たって図１９に示したサブルーチンを呼び出す。
【００６４】
図１９のサブルーチン（Ｓ５４５＿ＳＵＢ２）の処理内容を説明する。図１９のＳ１６１０では、呼び出し元の図１８のサブルーチンから波形グループＩＤを示す変数Ｇと波形ＩＤを示す変数Ｘを受け取る。初回の呼び出しなら変数Ｘ＝１かつ変数Ｇ＝１０００である。Ｓ１６２０以降のステップでは波形ＩＤ＝Ｘの波形と類似したエンジン圧力の時間波形をテーブルＴ５から見つけて、変数Ｇの値と紐づけてテーブルＴ６に格納していく。
【００６５】
Ｓ１６２０では、波形ＩＤ＝Ｘの時間波形と形状を比較する時間波形の波形ＩＤを示す変数Ｙの領域を記憶装置４５０に確保する。変数Ｙの初期値は変数Ｘと同じ値である。変数Ｘ＝１ならＹも１となる。図１９のサブルーチン内では、変数Ｘの値は一定のままで、比較対象の波形ＩＤ＝Ｙの値を１，２，３、と増やしながら波形ＩＤ＝Ｘと比較して類似波形か確認していく。
【００６６】
Ｓ１６３０では、波形ＩＤ＝ＹのセンサデータをテーブルＴ５から検索する。Ｙ＝１ならテーブルＴ５の１０１３が検索される。なお、Ｙの初期値はＸなので、波形ＩＤ＝Ｘの波形も検索されることになるが、グルーピング結果として波形ＩＤ＝Ｘの波形の波形グループＩＤもテーブルＴ６に格納する必要があるので、当該処理は正しい処理となる。
【００６７】
Ｓ１６４０ではＳ１６３０でレコードが見つかったか否かを判定する。波形ＩＤ＝Ｙのセンサデータが見つかればＳ１６５０に、見つからなければ全ての波形を比較完了したと判断し本サブルーチンを完了、図１８のＳ１５６０に進む。
【００６８】
Ｓ１６５０では波形ＩＤ＝Ｘの波形と波形ＩＤ＝Ｙの波形を比較して両者が類似する波形かどうかを判断する処理を行う。波形ＩＤ＝Ｘ，Ｙの波形が類似しているかの判断に関して、本実施の形態では次のように定義される類似度Ｓ（Ｘ、Ｙ）の大小に基づいて決定している。Ｓ（Ｘ、Ｙ）は、波形ＩＤ＝ＸのセンサデータであるＤｘ（ｋ）（ｋ＝1,2,3,…,Nx）と、波形ＩＤ＝ＹのセンサデータであるＤｙ（ｋ）（ｋ＝1,2,3,…,Ny）から次式Ａで計算するものとする。ただし、次式Ａでは、Ｄｘ（ｋ）とＤｙ（ｋ）に含まれるデータ点数に差異が生じる場合を考慮して、m=MIN(Nx, Ny)とする。すなわち、Ｄｘ（ｋ）とＤｙ（ｋ）のセンサデータ点数Nx,Nyのうち、データ点数の少ない方をmとする。なお、当該対応は、データ点数に差異がある場合の対応の一例に過ぎず、データ点数の少ない方の時間波形上にデータ点を加えることで双方のデータ点数を等しくする等、公知の対応が可能である。また、２つの波形の比較に際して、両者の時間範囲が大きく異なる場合には、そもそも両波形は類似しないと判定しても良いし、一方の波形の時間範囲に合わせて他方の波形をスケーリングして比較しても良い。
【００６９】
【数１】
2015072515.tif 000003
【００７０】
図１６のテーブルＴ５の例では、Ｄｘ（ｋ）、Ｄｙ（ｋ）はエンジン圧力のセンサ値である。波形ＩＤ：Ｘ＝１なら、Ｄｘ（ｋ）はセンサデータ１０１３に含まれるエンジン圧力の数値であり、Ｄｘ（１）＝0.5628、Ｄｘ（２）＝0.5727、となる。そして、比較対象が波形ＩＤ：Ｙ＝２なら、Ｄｙ（ｋ）はセンサデータ１０１６に含まれるエンジン圧力の数値であり、Ｄｙ（１）＝0.5699、Ｄｙ（２）＝0.5621、となる。Ｄｘ（ｋ）とＤｙ（ｋ）から波形ＩＤ＝Ｘ，Ｙの波形の類似度Ｓ（Ｘ，Ｙ）計算したらＳ１６６０に移行する。
【００７１】
Ｓ１６６０では、Ｓ１６５０で算出した類似度Ｓ（Ｘ、Ｙ）が本システム設計時に決定する閾値Ｓ０より大きければ、波形ＩＤ＝ＸとＹの波形が類似していると判断しＳ１６７０に進む。そうでなければＳ１６８０に進む。なお、閾値Ｓ０は事後的に変更可能としても良い。
【００７２】
Ｓ１６７０では、Ｓ１６６０で類似度が高いと判断した波形のＩＤをテーブルＴ６（図１７）に格納する。具体的には、波形ＩＤ＝Ｙの値をテーブルＴ６の第１列１０６０に波形ＩＤとして格納し、同じ行の第２列１０８０に変数Ｇの値を波形グループＩＤとして格納する。変数Ｇは、Ｓ１６１０において呼び出し元から受信した波形グループＩＤの値である。例えばテーブルＴ５の波形ＩＤ＝２の波形と波形ＩＤ＝３の波形が類似しており、両者の波形の類似度を計算した時の変数Ｇ（波形グループＩＤ）の値が２０００の場合、テーブルＴ６の２，３行目のようにデータを格納する。
【００７３】
これで波形ＩＤ＝ＸとＹの波形比較が完了したので、Ｓ１６８０では変数Ｙの値を＋１して次の波形比較の準備をする。例えばＸ＝１でＹ＝１の波形比較が完了したら、Ｙを２にして波形ＩＤ＝１と波形ＩＤ＝２の波形比較の準備をする。その後Ｓ１６３０に戻る。
【００７４】
以上で図１９のサブルーチンおよび図１８のＳ１５５０は完了し、図１８のＳ１５６０に移行する。
【００７５】
図１８のＳ１５６０では、波形比較元の波形ＩＤを示す変数Ｘを＋１して次の類似波形を探す準備をする。例えば波形ＩＤ：Ｘ＝１の類似波形を探し終わったら、波形ＩＤ：Ｘ＝２に更新して次の類似波形を探す準備を行う。
【００７６】
Ｓ１５７０では波形グループＩＤを示す変数Ｇを＋１０００して次の類似波形が属する波形グループＩＤを更新し、Ｓ１５３０に戻る。例えば波形ＩＤ＝１の類似波形を探し終わったら、波形グループＩＤ：Ｇ=１０００をＧ＝２０００に更新する。なお本実施の形態では一例として変数の値を１０００ずつ増加しているが増加量は任意の値で良い。
【００７７】
以上で図１８のサブルーチンと図８のＳ５４５の処理が完了し、Ｓ５５０に移行する。
【００７８】
Ｓ５５０では、テーブルＴ５（図１６）とテーブルＴ６（図１７）に基づいて、Ｓ５３８（図１５）で選択したセンサデータの時間波形（第２グラフ）とそのグルーピング結果を示す画面を作成する。図２０は、Ｓ５５０で作成された後にＳ５５５で表示される画面の一例を示す図である。この図の例では、テーブルＴ６に格納された波形グループＩＤが３つ（１０００，２０００，３０００）の場合を示しており、波形グループＩＤの数に合わせて画面内が３つの表示部１１００、１２００、１３００に分割されている。表示部１１００は、波形グループＩＤ＝１０００に属する波形を示す部分であり、テーブルＴ６で波形グループＩＤ＝１０００に属する波形ＩＤ＝１，２２の波形（第２グラフ）をそれぞれ波形１１１０、１１２０として表示している。波形１１１０，１１２０を表示するためのセンサデータは、テーブルＴ５から波形ＩＤ＝１，２２を検索キーとして検索できる。
【００７９】
なお、各波形ＩＤの時間波形を作成するに際して検索するセンサデータは、テーブルＴ５内で各波形ＩＤを有するセンサデータ、すなわちテーブルＴ４で定義された時間範囲に含まれるセンサデータだけでなく、当該時間範囲を前後に規定の時間だけ拡張して、当該拡張後の時間範囲に含まれるセンサデータを加えて時間波形を作成しても良い。この場合には、テーブルＴ４で定義された時間範囲と拡張分の時間範囲の判別がつくように、時間波形の表示を工夫することが好ましい。この種の表示の具体例としては、時間波形の背景を変えるもの（図２１中の波形１１１０Ａ参照）や、テーブルＴ４の時間範囲に含まれるセンサデータのみ時間波形上に点を描画するもの（図２１中の波形１１１０Ｂ参照）等がある。このように時間波形を表示する時間を拡張すると、散布図（図１１）上で選択したデータの前後でセンサデータがどのように変化しているかを把握でき、データの分析に資する場合がある。なお、上記の処理に代えて、それ以前の手順で時間範囲を拡張したセンサデータの時間波形が表示される処理をしても良い。
【００８０】
表示部１１００と同様なので詳細な説明は省略するが、表示部１２００は、波形グループＩＤ＝２０００に属する波形を示す部分であり、３つの波形１２１０，１２２０，１２３０が表示されている。表示部１３００は、波形グループＩＤ＝３０００に属する波形を示す部分であり、２つの波形１３１０，１３２０が表示されている。
【００８１】
なお、テーブルＴ６に格納された全ての波形グループＩＤに係る波形を表示できない場合には、画面をスクロール可能または他の画面に遷移可能に構成にするなどして、全てのグループＩＤの波形を表示可能にするものとする。また、波形グループＩＤの数だけを表示して、所望のＩＤをマウス等でクリックする等して選択することで、当該ＩＤに属する波形が全て表示されるように構成しても良い。すなわち、テーブルＴ６に含まれる波形グループＩＤの総数（グループの数）と、各グループに含まれる波形の形状が確認可能であれば、各波形の表示方法に特に限定は無い。
【００８２】
以上の画面作成処理が完了したらＳ５５５に移行する。Ｓ５５５ではＳ５５０で作成した図２０を表示装置４４０に表示しユーザに提供する。これにより図８に示した一連の処理が完了する。なお、他のセンサデータの散布図を表示したい場合には最初のＳ５０５に戻り、散布図上のデータ点の選択に変更は無いが、他のセンサデータの時間波形を表示したい場合には、Ｓ５３８に戻れば良い。
【００８３】
以上のように構成したデータ表示システムによれば、或る項目に係るデータの集合の一部の発生時刻周辺において時間変化を知りたい項目を指定することで、当該指定項目に係るデータの当該発生時刻周辺の時間変化の態様を容易に把握することができる。これにより、例えば、数秒〜数時間といった短期間に外れ値が発生している場合でも、その発生時刻周辺の短期間の時間波形グラフを表示し、外れ値発生時の波形から外れ値の主な原因を分析できる。
【００８４】
なお、上記で図８，１４，１８，１９等を利用して説明した各処理の順番は一例に過ぎず、上記効果が発揮できる範囲であれば順番の変更が適宜可能である。
【００８５】
また、上記では、図９、図１１、図１５、図２０の順番に画面が遷移する場合について説明したが、散布図と時間波形を表示すべきセンサが予め決まっている場合等には図９と図１５は省略可能である。また、図１１上においてカーソル７０５でデータを選択した後に、予め決められたショートカット操作をマウス、キーボード、タッチパネル等の入力装置４２５で行うことで図１５によるセンサの選択を代替して行い、図１５の表示を省略することも可能である。
【００８６】
また、上記では、図１１として散布図を表示したが、ヒストグラムやパレート図をはじめとするデータベース４１０内のセンサデータの傾向が把握できるグラフであれば、他のグラフを表示しても良い。その際には、図９の画面は、図１１で表示するグラフを規定するために必要な指標が適宜入力可能に作成することは言うまでもない。
【００８７】
また、上記では、図１１上でいわゆる外れ値を選択することで当該外れ値のトレンドを時間波形から把握して異常の原因を追究する場面を例に挙げて説明したが、本発明はこれだけに限らず、或る項目に係るデータの分布について複数の部分集合が把握された場合、各部分集合に含まれるデータが検出された時刻に係る所定のデータの変化を表示することで、当該部分集合に含まれるデータの傾向を視覚的に把握する場面で広く有用である。
【００８８】
ところで、図７に示した本システムに係る機能および当該機能を発揮するための実行処理等は、それらの一部又は全部をハードウェア（例えば各機能を実行するロジックを集積回路で設計する等）で実現しても良い。また、上記のデータ表示システムに係る構成は、処理装置（例えばＣＰＵ）によって読み出し・実行されることで当該システムの構成に係る各機能が実現されるプログラム（ソフトウェア）としてもよい。当該プログラムに係る情報は、例えば、半導体メモリ（フラッシュメモリ、ＳＳＤ等）、磁気記憶装置（ハードディスクドライブ等）及び記録媒体（磁気ディスク、光ディスク等）等に記憶することができる。
【００８９】
また、上記の実施の形態の説明では、制御線や情報線は、当該実施の形態の説明に必要であると解されるものを示したが、必ずしも製品に係る全ての制御線や情報線を示しているとは限らない。実際には殆ど全ての構成が相互に接続されていると考えて良い。
【００９０】
また、本発明は、上記の実施の形態に限定されるものではなく、その要旨を逸脱しない範囲内の様々な変形例が含まれる。例えば、本発明は、上記の実施の形態で説明した全ての構成を備えるものに限定されず、その構成の一部を削除したものも含まれる。
【符号の説明】
【００９１】
４１０…センサデータベース、４１５…発生頻度情報作成部、４２０…時間波形データ作成部、４２５…入力装置、４３０…データ集約部、４３５…類似波形検索部、４４０…表示装置、４４５…処理装置、４５０…記憶装置、Ｓ…類似度
【図１】
2015072515.tif 000004
【図２】
2015072515.tif 000005
【図３】
2015072515.tif 000006
【図４】
2015072515.tif 000007
【図５】
2015072515.tif 000008
【図６】
2015072515.tif 000009
【図７】
2015072515.tif 000010
【図８】
2015072515.tif 000011
【図９】
2015072515.tif 000012
【図１０】
2015072515.tif 000013
【図１１】
2015072515.tif 000014
【図１２】
2015072515.tif 000015
【図１３】
2015072515.tif 000016
【図１４】
2015072515.tif 000017
【図１５】
2015072515.tif 000018
【図１６】
2015072515.tif 000019
【図１７】
2015072515.tif 000020
【図１８】
2015072515.tif 000021
【図１９】
2015072515.tif 000022
【図２０】
2015072515.tif 000023
【図２１】
2015072515.tif 000024