JP5437328B2 - 観測値信頼度評価装置、観測値信頼度評価方法及び観測値信頼度評価プログラム - Google Patents

Description

本発明は、正確さや誤差分布が様々なセンサデータから算出した仮想センサデータを用いて閾値判定を行う際に、信頼度という指標を用いる技術に関する。

一般的なユビキタスアプリケーションにおけるイベント検出処理には、センサデータに含まれる観測値の閾値判定処理がよく用いられている。

従来、閾値判定に用いるセンサデバイスは、予め決まっていることが前提であった。観測値の正確さや誤差範囲の見積りは容易に得られ、これらを活用して閾値を最適化したり、観測値が誤差範囲内に存在する確率を計算したりすることで、安定した閾値判定を実現していた（特許文献１参照）。

しかしながら、閾値判定に用いるセンサデバイスやその領域・位置が予め決まっていない環境においては、観測値の正確さや誤差範囲が刻々と変化し、安定した閾値判定が困難であるという問題を有する。

本発明が対象としている環境に、複数の移動ユーザにセンサデータを収集させる参加型センシングがある。参加型センシング環境において、ある領域を代表する仮想センサデバイスを定義し、閾値判定処理を行うことを考える（非特許文献１参照）。ある領域を代表する仮想センサデータは、各時刻におけるその領域に存在する複数のセンサデバイスについて、センサデータの中央値や平均値などの統計値を利用することができる。

従来、観測値とその誤差範囲を見積もる方法として、１回の観測に対して複数回のサンプリングを行い、それらの統計値を用いることが有効とされてきた。また、近隣の複数の観測値に相関があれば、これらの統計値を用いることも同様に有効である（非特許文献２参照）。

しかし、測定するセンサデバイスが予め決まっていない参加型センシング環境で得られる仮想センサデータは、正確さやその誤差範囲が時刻によってダイナミックに変化するという特徴がある。また、ある領域を代表する仮想センサデータは、その領域自体がユーザの要求により可変となるという特徴がある。

特開平８−２０２９７０号公報

武山政直、他１名、「参加型都市センシングによる価値共創モデルの可能性」、会誌「情報処理」、Ｖｏｌ．５１、Ｎｏ．９、情報処理学会、２０１０年 吉澤康和著、「新しい誤差論」、共立出版株式会社、１９８９年 独立行政法人製品評価技術基盤機構認定センター著、「不確かさの入門ガイド」、２００４年

前記述べたような特徴をもつ参加型センシング環境において、閾値判定を安定して動作させるにあたって、以下に述べる問題があった。

仮想センサデータは不確かさをもつデータであるため、閾値判定を行う場合、仮想センサデータが閾値を超えたかどうかだけでなく、仮想センサデータの誤差範囲と、その区間内に仮想センサデータが含まれる確率を考慮することが望ましい（非特許文献３参照）。

しかしながら、仮想センサデータを算出するためのセンサデータ集合が予め決まっていないことにより、センサデバイスの違いや観測位置や観測時刻のばらつきによって、算出された確率の確からしさが変動するという問題がある。また、ユーザ要求によって仮想センサデータを算出する領域が伸縮することで、仮想センサデータの確からしさが変動するという問題がある。

そこで、前記課題を解決するために、本発明は、仮想センサデータや閾値判定の確からしさを評価することにより、仮想センサデータや閾値判定の確からしさが変動するときであっても、より現実的な閾値判定を行うことを目的とする。

上記目的を達成するために、仮想センサデータを構成しながら、センサデバイスの違いや観測位置や観測時刻のばらつきを定量化する信頼度という指標を導入することで、仮想センサデータの信頼度を考慮した確率的な閾値判定を行うこととした。

具体的には、本発明は、複数の観測装置における観測値について、代表値を算出するとともに、前記複数の観測装置における観測位置及び観測時刻の少なくとも一方を含む観測条件について、代表値を算出し、前記観測条件についての代表値からの前記観測条件のばらつき、又は所定値からの前記観測条件についての代表値のかたよりを算出する統計処理部と、前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値の信頼度を算出する代表値信頼度算出部と、を備えることを特徴とする観測値信頼度評価装置である。

また、本発明は、複数の観測装置における観測値について、代表値を算出するとともに、前記複数の観測装置における観測位置及び観測時刻の少なくとも一方を含む観測条件について、代表値を算出し、前記観測条件についての代表値からの前記観測条件のばらつき、又は所定値からの前記観測条件についての代表値のかたよりを算出する統計処理手順と、前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値の信頼度を算出する代表値信頼度算出手順と、を順に備えることを特徴とする観測値信頼度評価方法である。

この構成によれば、複数の観測装置における観測条件にばらつきやかたよりがあるときでも、ばらつきやかたよりに基づいて複数の観測装置における観測値についての代表値の信頼度を算出するため、仮想センサデータの信頼度を考慮することができる。

また、本発明は、前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出する閾値判定信頼度算出部、をさらに備えることを特徴とする観測値信頼度評価装置である。

また、本発明は、前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出する閾値判定信頼度算出手順、を前記代表値信頼度算出手順の後にさらに備えることを特徴とする観測値信頼度評価方法である。

この構成によれば、複数の観測装置における観測条件にばらつきやかたよりがあるときでも、ばらつきやかたよりに基づいて複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出するため、仮想センサデータの信頼度を考慮した確率的な閾値判定を行うことができる。

また、本発明は、前記閾値判定信頼度算出部は、前記複数の観測装置における観測値についての代表値からの観測値のばらつきに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出することを特徴とする観測値信頼度評価装置である。

また、本発明は、前記閾値判定信頼度算出手順では、前記複数の観測装置における観測値についての代表値からの観測値のばらつきに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出することを特徴とする観測値信頼度評価方法である。

この構成によれば、仮想センサデータが閾値を超えたかどうかだけでなく、仮想センサデータの誤差範囲内に仮想センサデータが含まれる確率を考慮することができる。

また、本発明は、前記複数の観測装置における観測条件に応じて、前記複数の観測装置が観測したデータをグループに分類するデータグループ化部、をさらに備え、前記統計処理部及び前記代表値信頼度算出部は、前記データグループ化部が各グループに分類した各データについて処理を実行することを特徴とする観測値信頼度評価装置である。

この構成によれば、複数の観測装置における観測条件ごとに細密に、複数の観測装置が観測したデータについて仮想センサデータを算出することができる。

また、本発明は、上記の観測値信頼度評価方法をコンピュータに実行させるための観測値信頼度評価プログラムである。

この構成によれば、観測値信頼度評価プログラムを、記録媒体に記録することも、ネットワークを通して提供することも可能である。

本発明は、仮想センサデータや閾値判定の確からしさを評価することにより、仮想センサデータや閾値判定の確からしさが変動するときであっても、より現実的な閾値判定を行うことができる。

仮想センサデータ構成装置の構成を示す図である。 センサデータの内容を例示する図である。 グループ化の条件を例示する図である。 センサデータグループ化機能部の処理を示す図である。 信頼度付き仮想センサデータ構成機能部の処理を示す図である。 センサデータグループの内容を例示する図である。 メタデータの統計情報の内容を例示する図である。 メタデータの統計情報付きのセンサデータグループの内容を例示する図である。 観測値の統計情報の内容を例示する図である。 観測値の統計情報付きのセンサデータグループの内容を例示する図である。 仮想センサデータの内容を例示する図である。 信頼度付きの仮想センサデータの内容を例示する図である。 信頼度付き閾値判定機能部の処理を示す図である。 信頼度付き閾値判定機能部の処理を例示する図である。 信頼度付きの仮想センサデータの内容を例示する図である。

添付の図面を参照して本発明の実施形態を説明する。以下に説明する実施形態は本発明の実施の例であり、本発明は以下の実施形態に制限されるものではない。

仮想センサデータ構成装置の構成を示す図を図１に示す。信頼度付き閾値判定を行う信頼度付き仮想センサデータ構成装置１００は、センサデータ読み出し機能部１０１、仮想センサデータ蓄積機能部１０２、設定受付機能部１０３、応答機能部１０４、センサデータグループ化機能部１０５、信頼度付き仮想センサデータ構成機能部１０６及び信頼度付き閾値判定機能部１１１を具備する。信頼度付き仮想センサデータ構成機能部１０６は、メタデータ統計処理機能部１０７、観測値統計処理機能部１０８、仮想センサデータ構成機能部１０９及び信頼度付与機能部１１０を含む。

信頼度付き閾値判定を行う信頼度付き仮想センサデータ構成装置１００は、センサデータ読み出し機能部１０１及び仮想センサデータ蓄積機能部１０２において、センサデータ記憶装置１１２と接続され、設定受付機能部１０３及び応答機能部１０４において、アプリケーション装置１１３に接続される。センサデータ記憶装置１１２は、リレーショナルデータベース管理システム（ＲＤＢＭＳ）などを用いて容易に作成可能なものである。

信頼度付き仮想センサデータ構成装置１００の概要を説明する。センサデータグループ化機能部１０５は、データグループ化部として、複数の観測装置における観測条件に応じて、複数の観測装置が観測したデータをグループに分類する。観測装置として、温度センサなどが挙げられ、観測条件として、観測位置や観測時刻などが挙げられる。

メタデータ統計処理機能部１０７は、統計処理部として、複数の観測装置における観測条件について、代表値を算出し、代表値からの観測条件のばらつき又は所定値からの代表値のかたよりを算出する。観測条件の代表値として、観測位置や観測時刻の平均値などが挙げられ、所定値として、グループ領域における中心値などが挙げられる。観測値統計処理機能部１０８は、統計処理部として、複数の観測装置における観測値について、代表値を算出する。観測値の代表値として、温度の平均値などが挙げられる。

信頼度付与機能部１１０は、代表値信頼度算出部として、上記ばらつき又は上記かたよりに基づいて、複数の観測装置における観測値についての代表値の信頼度を算出する。例えば、信頼度付与機能部１１０は、グループにおいて観測条件のばらつきが一様であるほど、観測値の代表値の信頼度が高いと評価してもよく、グループにおいて観測条件の代表値が観測条件の中心値に近いほど、観測値の代表値の信頼度が高いと評価してもよい。

信頼度付き閾値判定機能部１１１は、閾値判定信頼度算出部として、上記ばらつき又は上記かたよりに基づいて、複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出する。例えば、信頼度付き閾値判定機能部１１１は、グループにおいて観測条件のばらつきが一様であるほど、閾値判定の信頼度が高いと評価してもよく、グループにおいて観測条件の代表値が観測条件の中心値に近いほど、閾値判定の信頼度が高いと評価してもよい。

信頼度付き閾値判定機能部１１１は、閾値判定信頼度算出部として、複数の観測装置における観測値についての代表値からの観測値のばらつきに基づいて、複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出する。つまり、信頼度付き閾値判定機能部１１１は、代表値が閾値を超えたかどうかだけでなく、代表値の誤差範囲内に代表値が含まれる確率を考慮している。

図１に説明した各構成は、コンピュータとプログラムによっても実現でき、プログラムを記録媒体に記録することも、ネットワークを通して提供することも可能である。

センサデータの内容を例示する図を図２に示す。センサデータ２００は、属性と値の組で記述されるデータであり、その属性には観測位置２０１、２０２と観測時刻２０３と観測値２０４を含んでいる。観測位置は少なくとも１つの属性からなり、緯度と経度のように複数あってもよい。また、緯度経度に限るものではなく、例えば住所区画でもよい。

センサデータ２００の観測条件はこれに限るものでなく、観測位置及び観測時刻の両方が含まれてもよく、どちらか片方が含まれなくてもよい。また、その他のメタデータとして、例えば、センサの種別、センサの型番、ネットワークアドレスなどを含んでもよい。

グループ化の条件を例示する図を図３に示す。グループ化条件３００は、グループ化するセンサの種類３０１、グループ化する観測開始位置３０２、グループ化する観測終了位置３０３、グループ化する観測開始時刻３０４及びグループ化する観測終了時刻３０５を含み、時空間で表される領域を定義している。

本実施形態では、グループ化する位置範囲の指定に、緯度経度の矩形を表す開始位置と終了位置を用いているが、例えば、住所区画、建築物内の区画、緯度経度の中央と半径などを用いることができる。本実施形態では、グループ化する時刻範囲の指定に、開始時刻と終了時刻を用いているが、例えば、中央時刻と前後の時間、あるいは、開始時刻と終了時刻とグループ化間隔、などで指定してもよい。

例えば、
開始時刻＝｛２０１１−０５−０１ ００：００：００｝
終了時刻＝｛２０１１−０５−０１ ００：３５：００｝
グループ化間隔＝｛００：１０：００｝
と指定した場合、２０１１−０５−０１ ００：００：００から１０分間隔でグループ化する。結果として、
２０１１−０５−０１ ００：００：００−００：０９：５９、
２０１１−０５−０１ ００：１０：００−００：１９：５９、
２０１１−０５−０１ ００：２０：００−００：２９：５９、
の３グループを作成し、残り時間のグループはグループ化されてもよくされなくてもよい。また、複数のグループを指定する場合は、複数個のグループ化条件を列挙すればよい。

センサデータグループ化機能部の処理を示す図を図４に示す。センサデータグループ化機能部１０５は、設定受付機能部１０３を介して、グループ化条件を受け取り、読み出し要求を生成する（手順Ｓ４０１）。センサデータグループ化機能部１０５は、読み出し要求をセンサデータ読み出し機能部１０１に送信し、センサデータ記憶装置１１２から読み出したセンサデータを受信する（手順Ｓ４０２）。

センサデータグループ化機能部１０５は、手順Ｓ４０１で受信したセンサデータのグループ化条件に基づいてグループを作成する（手順Ｓ４０３）。手順Ｓ４０１〜Ｓ４０３により、図３に示したセンサデータグループ化条件３００に基づいて、図６に示すようなセンサデータグループ６００を作成する。

センサデータグループ化機能部１０５は、グループごとに、手順Ｓ４０５〜Ｓ４０８を繰り返す（手順Ｓ４０４）。グループが空である場合は（手順Ｓ４０５において「はい」）、手順Ｓ４０６〜Ｓ４０８を省略する。グループが空でない場合は（手順Ｓ４０５において「いいえ」）、手順Ｓ４０６〜Ｓ４０８を実行する。

信頼度付き仮想センサデータ構成機能部１０６の実行（手順Ｓ４０６）の詳細は後述する。手順Ｓ４０６が完了すると、信頼度付き仮想センサデータ構成機能部１０６は、信頼度付き閾値判定機能部１１１に仮想センサデータを送信し（手順Ｓ４０７）、仮想センサデータ蓄積機能部１０２に仮想センサデータを送信する（手順Ｓ４０８）。仮想センサデータ蓄積機能部１０２は、センサデータ記憶装置１１２に仮想センサデータを保存する。

信頼度付き仮想センサデータ構成機能部の処理を示す図を図５に示す。信頼度付き仮想センサデータ構成機能部１０６は、センサデータグループ６００を受け取ると、欠損データの補正を行う（手順Ｓ５０１）。欠損データの補正は、観測値とメタデータの属性のうち１つ以上が欠損しているデータを削除することで行う。図６の欠損値６０５は観測値のみが欠損した例である。この行は欠損データの補正によって削除される。

補正方法はこれに限らず、例えば、センサデバイスが測定しうる値の範囲外の値を含むデータを削除する、あるいは、カルマンフィルタ、線形補完、スプライン補完、クリギングなどの手段を用いて欠損データの補完を行う、などの手段を用いることができる。

メタデータ統計処理機能部１０７は、手順Ｓ５０１で作成したセンサデータグループ６００を受け取り、観測位置および観測時刻のそれぞれの統計値を計算する。そして、それぞれの統計値の組からなるセンサデータグループ６００のメタデータの統計情報７００と、各センサデータに対してメタデータの統計処理結果を与えたセンサデータグループ８００を、観測値統計処理機能部１０８に送信する（手順Ｓ５０２）。

メタデータの統計情報の内容を例示する図を図７に示す。メタデータの統計情報７００には、平均値７０１、標準偏差７０２、歪度７０３、尖度７０４、正規化平均値７０５及び正規化標準偏差７０６を用いている。

正規化平均値７０５及び正規化標準偏差７０６は、空間範囲及び時間範囲で定義された領域の各軸において、観測開始点を−１、観測終了点を１とする正規化を行った単位領域における、メタデータの平均値及び標準偏差である。また、観測時刻の統計値７０９は、グループ化開始時刻２０１１／５／１ ０１：５５：００（３０４）を０とした経過秒数で算出している。メタデータの統計情報７００は、これらの値に限らず、例えば、中央値、不偏分散、値の重心などを含んでもよい。

メタデータの統計情報付きのセンサデータグループの内容を例示する図を図８に示す。図６のセンサデータグループ６００に、緯度の偏差８０１、経度の偏差８０２及び時刻の偏差８０３を付与している。各センサデータに与える統計処理結果は、これらの属性に限るものではなく、例えば、重心からの距離などを用いてもよい。

観測値統計処理機能部１０８は、手順Ｓ５０２で作成したセンサデータグループ８００とメタデータの統計情報７００を受け取り、観測値の統計値を計算する。そして、手順Ｓ５０２で作成したメタデータの統計情報７００と、観測値の統計値の組からなるセンサデータグループ８００の観測値の統計情報９００と、手順Ｓ５０２でメタデータの統計処理結果を与えた各センサデータに対しさらに観測値の統計処理結果を与えたセンサデータグループ１０００を、仮想センサデータ構成機能部１０９に送信する（手順Ｓ５０３）。

観測値の統計情報の内容を例示する図を図９に示す。観測値の統計情報９００には、平均値９０１、標準偏差９０２、歪度９０３及び尖度９０４を用いている。

観測値の統計情報９００は、これらの値に限らず、例えば、中央値、不偏分散などを含んでもよい。これらの統計値の算出に、手順Ｓ５０２で作成した、メタデータの統計情報７００や、各センサデータの統計情報８０１−８０３を用いてもよい。例えば、観測値の統計情報９００は、メタデータの偏差を用いた重みを付与した荷重平均値、メタデータの重心からの距離を用いた重みを付与した荷重平均値、などを含んでもよい。

観測値の統計情報付きのセンサデータグループの内容を例示する図を図１０に示す。図８のセンサデータグループ８００に、観測値の偏差１００１を付与している。

本実施形態では、メタデータ統計処理機能部１０７と観測値統計処理機能部１０８の順序で計算している（手順Ｓ５０２、Ｓ５０３）。しかし、これらの計算順序を入れ替えてもよい。そして、観測値の統計情報９００や、各センサデータの統計情報１００１を用いて、メタデータの統計情報７００を算出してもよい。

仮想センサデータ構成機能部１０９は、手順Ｓ５０２と手順Ｓ５０３で算出した統計情報から、観測値の代表値、観測値の誤差範囲及びメタデータの代表値を選択し、これらを結合した仮想センサデータを構成する（手順Ｓ５０４）。

仮想センサデータの内容を例示する図を図１１に示す。仮想センサデータ１１００は、グループＩＤ１１０１、メタデータの代表値１１０２−１１０４、観測値の代表値１１０５及び観測値の誤差範囲１１０６で構成する。仮想センサデータ１１００は、この構成に限らず、例えば、メタデータの代表値１１０２−１１０４をすべては含まなくてもよい。

観測値の代表値１１０５と誤差範囲１１０６として、観測値の平均値９０１と標準偏差９０２の２倍の値（９５％信頼区間）を用いる。観測値の代表値と誤差範囲に用いる統計値として、観測値の平均値と標準偏差の組に限らず、例えば、観測値の荷重平均値と荷重標準偏差の組などを用いることができる。また、標準偏差９０２の２倍の値（９５％信頼区間）に限らず、標準偏差９０２の１倍の値（６８％信頼区間）を用いてもよい。

メタデータの代表値１１０２−１１０４として、メタデータの平均値７０１を用いる。メタデータの代表値には、メタデータの平均値だけではなく、例えば、メタデータの中央値、あるいは、メタデータの重み付き平均値などを用いることができる。

信頼度付与機能部１１０は、手順Ｓ５０４で構成した仮想センサデータ１１００に信頼度を付与する（手順Ｓ５０５）。信頼度は、観測位置及び観測時刻のかたより具合及びばらつき具合という、２つの指標を定量化した値とする。

本実施形態においては、上記の正規化後の単位領域に観測条件が一様に散らばっているほど、信頼度が高いと仮定する。なお、上記の正規化後の単位領域の中心部に観測条件が集中しているほど、信頼度が高いと仮定してもよい。

２つの指標の定量化には、メタデータの統計情報７００及び観測値の統計情報９００を用いる。本実施形態においては、単位領域の中心点から正規化平均値７０５までの距離を用いて、かたより値を算出し、正規化標準偏差７０６を用いて、ばらつき値を算出する。

まず、単位領域の中心点と単位領域の各軸における正規化平均値７０５の距離ｒを用いて、中心点からのかたより値を算出する。本実施形態においては、距離ｒが取りうる範囲は０≦ｒ≦√３である。メタデータの統計情報７００から算出される距離ｒは０．２５４であり、これをｒが取りうる値の最大値で割るとかたより値は１４．７％である。

次に、単位領域における各軸の正規化標準偏差σ＿ｘ（ｘ＝ｉ，ｊ，ｋ）を用いて、ばらつき値を算出する。一辺が２σ＿ｘ（ｘ＝ｉ，ｊ，ｋ）の直方体を仮定し、直方体の単位領域に占める割合をばらつき値とする。本実施形態においては、体積ｖの取りうる範囲は０≦ｖ≦８である。メタデータの統計情報７００から算出される直方体の体積ｖは０．０６１９であり、これをｖが取りうる最大値で割るとばらつき値は０．７７３％である。

信頼度付きの仮想センサデータの内容を例示する図を図１２に示す。かたより値１２０１とばらつき値１２０２の２値で信頼度を表す。また、かたより値とばらつき値の算出には、正規化平均値と正規化標準偏差に限らず、例えば、センサデータグループの単位領域における中心の周りのモーメント、メタデータの歪度と尖度、センサデータグループの単位領域におけるエントロピーなどを用いてもよい。

また、本発明の範囲は上記のものに限らず、例えば、メタデータの緯度、経度、時刻及び観測値からなる４次元分布を仮定し、その信頼区間に入る確率を信頼度としてもよい。また、センサの個数の密度を算出して重み付けを行ってもよい。また、信頼度は、１値に合成してもよく、複数の値からなるベクトル表現のままでもよい。

信頼度付き仮想センサデータ構成機能部１０６は、信頼度を付与した仮想センサデータ１２００を、信頼度付き閾値判定機能部１１１に送信する（手順Ｓ５０６）。そして、信頼度を付与した仮想センサデータ１２００を、仮想センサデータ蓄積機能部１０２を介し、センサデータ記憶装置１１２に蓄積する（手順Ｓ５０７）。

信頼度付き閾値判定機能部の処理を示す図を図１３に示す。信頼度付き閾値判定機能部の処理を例示する図を図１４に示す。信頼度付きの仮想センサデータの内容を例示する図を図１５に示す。図１４には、時刻ｔ１，ｔ２，ｔ３における、仮想センサデータの値、仮想センサデータの誤差範囲、仮想センサデータの信頼度及びアプリケーションの閾値を図示している。図１５には、図１４の仮想センサデータ１４００を例示している。

アプリケーション装置１１３は、仮想センサデータが閾値以上であるかどうかの判定結果を要求する。時刻ｔ１においては、誤差範囲の下限からの１／４点で、アプリケーションの閾値と交差している。時刻ｔ２においては、誤差範囲の下限からの３／４点で、アプリケーションの閾値と交差している。時刻ｔ３においては、交差が起こっていない。

信頼度付き閾値判定処理機能部１１１は、第一段の判定として、誤差範囲と閾値から閾値判定の計算を行う（手順Ｓ１３０１）。グループ内の観測値の誤差分布は、正規分布に従うとみなす。測定値ｘの誤差範囲の幅は、９５％信頼区間を用いる。つまり、平均値μ、標準偏差σのとき、（μ−２σ≦ｘ≦μ＋２σ）の区間を用いる。

閾値ａを上回るという判定を行う場合、測定値ｘが閾値以上誤差範囲上限以下（ａ≦ｘ≦μ＋２σ）に存在する確率を計算する。このとき、（μ＋２σ≦ａ）の場合は０％、（ａ≦μ−２σ）の場合は１００％とみなす。

閾値ａを下回るという判定を行う場合、測定値ｘが閾値以下誤差範囲下限以上（μ−２σ≦ｘ≦ａ）に存在する確率を計算する。このとき、（μ＋２σ≦ａ）の場合は１００％、（ａ≦μ−２σ）の場合は０％とみなす。

本実施形態では、閾値を上回るという判定を行っている。時刻ｔ１においては、閾値以上かつ誤差範囲上限以下に存在する確率は８１．８５％である。時刻ｔ２においては、閾値以上かつ誤差範囲上限以下に存在する確率は１３．５９％である。時刻ｔ３においては、閾値以上かつ誤差範囲上限以下に存在する確率は０％である。

信頼度付き閾値判定処理機能部１１１は、第二段の判定として、信頼度（かたより値及びばらつき値）を用いて、第一段閾値判定結果の正しさを評価する（手順Ｓ１３０２）。つまり、代表値の信頼度が低い場合、その誤差範囲の信頼度も同様に低いとみなされるため、第一段閾値判定結果が正しく無い可能性がある。

本実施形態においては、上記の正規化後の単位領域に観測条件が一様に散らばっている場合に、観測値の代表値の信頼度が高いとみなすこととする。つまり、観測条件が単位領域に一様に散らばっている場合は、観測値の代表値の信頼度は１００％とし、センサの配置が局所的である場合は、観測値の代表値の信頼度は低いとみなす。

本実施形態においては、仮想センサデータに付与されている信頼度を用いて、ばらつき値×（１−かたより値）を算出する。（１−かたより値）を乗算することにより、観測条件が単位領域に一様に散らばっており、観測条件の代表値が単位領域の中心値に近いほど、観測値の代表値の信頼度が高いと、より確実にみなすことができる。

ばらつき値×（１−かたより値）を評価する。時刻ｔ１における第一段閾値判定結果は、０．０２％の信頼度である。時刻ｔ２における第一段閾値判定結果は、５１．２％の信頼度である。時刻ｔ３における第一段閾値判定結果は、０．５９％の信頼度である。

信頼度付き閾値判定処理機能部１１１は、第一段の判定結果及び第二段の判定結果を用いて、閾値判定を統合する（手順Ｓ１３０３）。

本実施形態では、第一段判定結果の確率が０％より大きいものは「閾値を超えている可能性がある」とみなし、その確率を以て第一段判定結果の不確かさを表現する。また、第二段判定結果を以て第一段判定結果の正しさを表現する。

「第一段判定結果の確率が０％であり、その信頼度が低い」場合は、閾値を超えている可能性があるとみなす。「閾値を超えている確率が０％であり、その信頼度が高い」場合のみ、閾値を超えていないという判定となる。

時刻ｔ１においては、閾値を超えている可能性がある。その確率は第一段で判定したように８１．８５％であり、第一段の信頼度は第二段で判定したように低い。時刻ｔ２においては、閾値を超えている可能性がある。その確率は第一段で判定したように１３．５９％であり、第一段の信頼度は第二段で判定したように高い。時刻ｔ３においては、閾値を超えている可能性がある。その確率は第一段で判定したように０％であるが、第一段の信頼度は第二段で判定したように低いため、閾値を超えている可能性があると考えられる。

信頼度付き閾値判定処理機能部１１１は、手順Ｓ１３０３で得られた判定結果に信頼度を付与し、応答機能部１０４を介しアプリケーション装置１１３に送信する（手順Ｓ１３０４）。本実施形態では、信頼度付き閾値判定処理機能部１１１は、信頼度として、かたより値１２０１及びばらつき値１２０２の２値を付与する。

本発明に係る観測値信頼度評価装置、観測値信頼度評価方法及び観測値信頼度評価プログラムは、特に、閾値判定に用いるセンサデバイスやその領域・位置が予め決まっていない環境において、例えば、複数の移動ユーザにセンサデータを収集させる参加型センシング環境において、安定した閾値判定を行うときに適用することができる。

１００：仮想センサデータ構成装置
１０１：センサデータ読み出し機能部
１０２：仮想センサデータ蓄積機能部
１０３：設定受付機能部
１０４：応答機能部
１０５：センサデータグループ化機能部
１０６：信頼度付き仮想センサデータ構成機能部
１０７：メタデータ統計処理機能部
１０８：観測値統計処理機能部
１０９：仮想センサデータ構成機能部
１１０：信頼度付与機能部
１１１：信頼度付き閾値判定機能部
１１２：センサデータ記憶装置
１１３：アプリケーション装置

Claims (8)

1. 複数の観測装置における観測値について、代表値を算出するとともに、前記複数の観測装置における観測位置及び観測時刻の少なくとも一方を含む観測条件について、代表値を算出し、前記観測条件についての代表値からの前記観測条件のばらつき、又は所定値からの前記観測条件についての代表値のかたよりを算出する統計処理部と、
   前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値の信頼度を算出する代表値信頼度算出部と、
   を備えることを特徴とする観測値信頼度評価装置。
2. 前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出する閾値判定信頼度算出部、
   をさらに備えることを特徴とする、請求項１に記載の観測値信頼度評価装置。
3. 前記閾値判定信頼度算出部は、前記複数の観測装置における観測値についての代表値からの観測値のばらつきに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出することを特徴とする、請求項２に記載の観測値信頼度評価装置。
4. 前記複数の観測装置における観測条件に応じて、前記複数の観測装置が観測したデータをグループに分類するデータグループ化部、をさらに備え、
   前記統計処理部及び前記代表値信頼度算出部は、前記データグループ化部が各グループに分類した各データについて処理を実行することを特徴とする、請求項１から請求項３のいずれかに記載の観測値信頼度評価装置。
5. 複数の観測装置における観測値について、代表値を算出するとともに、前記複数の観測装置における観測位置及び観測時刻の少なくとも一方を含む観測条件について、代表値を算出し、前記観測条件についての代表値からの前記観測条件のばらつき、又は所定値からの前記観測条件についての代表値のかたよりを算出する統計処理手順と、
   前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値の信頼度を算出する代表値信頼度算出手順と、
   を順に備えることを特徴とする観測値信頼度評価方法。
6. 前記ばらつき又は前記かたよりに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出する閾値判定信頼度算出手順、
   を前記代表値信頼度算出手順の後にさらに備えることを特徴とする、請求項５に記載の観測値信頼度評価方法。
7. 前記閾値判定信頼度算出手順では、前記複数の観測装置における観測値についての代表値からの観測値のばらつきに基づいて、前記複数の観測装置における観測値についての代表値及び閾値の大小判定の信頼度を算出することを特徴とする、請求項６に記載の観測値信頼度評価方法。
8. 請求項５から請求項７のいずれかに記載の観測値信頼度評価方法をコンピュータに実行させるための観測値信頼度評価プログラム。

Images