(57)【要約】      （修正有）【課題】物質のサンプリングに適した場所におけるサンプリングを支援することができるサンプリング支援装置、システムおよび方法を提供する。【解決手段】サンプリング支援装置１は、外部環境の空間３においてサンプリングの対象とする物質に関する情報を出力する。サンプリング支援装置は、物質センサ１０と、位置センサと、制御部と、情報出力部を備える。物質センサは、物質の濃度に応じたセンサ値を測定する。位置センサは、空間中でセンサ値が測定された位置を測定する。制御部は、測定されたセンサ値及び位置に基づいて、空間中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出する。情報出力部は、算出された推定分布に基づいて、物質センサが空間中で次にセンサ値を測定する箇所に誘導する誘導情報を出力する。【選択図】図１
2020041996.tif 000002
【特許請求の範囲】
【請求項１】
外部環境の空間においてサンプリングの対象とする物質に関する情報を出力するサンプリング支援装置であって、
前記物質の濃度に応じたセンサ値を測定する物質センサと、
前記空間中で前記センサ値が測定された位置を測定する位置センサと、
測定されたセンサ値及び位置に基づいて、前記空間中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出する制御部と、
算出された推定分布に基づいて、前記物質センサが前記空間中で次にセンサ値を測定する箇所に誘導する誘導情報を出力する情報出力部と
を備えるサンプリング支援装置。
【請求項２】
前記制御部は、算出した推定分布においてセンサ値が極大となることが推定される箇所に応じて、前記誘導情報を生成する
請求項１に記載のサンプリング支援装置。
【請求項３】
前記制御部は、前記物質センサおよび前記位置センサから逐次、測定されたセンサ値及び位置を取得して、前記誘導情報を更新する
請求項１又は２に記載のサンプリング支援装置。
【請求項４】
前記情報出力部は、前記空間中で前記物質をサンプリングする箇所に誘導するサンプリング誘導情報を、更に出力する
請求項１〜３のいずれか１項に記載のサンプリング支援装置。
【請求項５】
前記制御部は、
算出した推定分布に基づき、前記空間中の位置毎に前記センサ値を測定する価値を示す価値関数を演算し、
演算した価値関数に基づき、前記誘導情報を生成する
請求項１〜４のいずれか１項に記載のサンプリング支援装置。
【請求項６】
前記誘導情報は、前記空間における前記価値関数の分布を提示する
請求項５に記載のサンプリング支援装置。
【請求項７】
前記誘導情報は、前記価値関数が最大値となる位置を提示する
請求項５に記載のサンプリング支援装置。
【請求項８】
前記制御部は、測定されたセンサ値及び位置に基づくガウス過程回帰により、前記推定分布を算出する
請求項１〜７のいずれか１項に記載のサンプリング支援装置。
【請求項９】
前記センサ値が測定された位置からの距離に応じて推定されるセンサ値の変動幅を規定する推定情報を記憶する記憶部をさらに備え、
前記制御部は、前記推定情報を参照して、前記推定分布を算出する
請求項１〜７のいずれか１項に記載のサンプリング支援装置。
【請求項１０】
前記情報出力部は、前記空間を示す背景に重畳させて、前記誘導情報を表示する
請求項１〜９のいずれか１項に記載のサンプリング支援装置。
【請求項１１】
前記情報出力部は、表示部、スピーカ、振動デバイス、及びプロジェクタのうちの少なくとも１つを含む
請求項１〜１０のいずれか１項に記載のサンプリング支援装置。
【請求項１２】
前記物質センサは、モル質量に基づく密度センサ、分光光度計、吸光光度計、屈折度計、ガスセンサ、ｐＨ計、および電気伝導度計の少なくとも１つを含む
請求項１〜１１のいずれか１項に記載のサンプリング支援装置。
【請求項１３】
請求項１〜１２のいずれか１項に記載のサンプリング支援装置と、
前記物質をサンプリングするサンプラーと
を備えるシステム。
【請求項１４】
外部環境の空間から対象とする物質をサンプリングする方法であって、
物質センサにより、前記物質の濃度に応じたセンサ値を測定するステップと、
位置センサにより、前記空間中で前記センサ値が測定された位置を測定するステップと、
制御部により、測定されたセンサ値及び位置に基づいて、前記空間中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出するステップと、
情報出力部により、算出された推定分布に基づいて、前記物質センサが前記空間中で次にセンサ値を測定する箇所に誘導する誘導情報を出力するステップと、
前記誘導情報に応じて、前記センサ値の測定及び前記位置の測定を繰り返し、前記空間中で前記物質をサンプリングする箇所を探索するステップと、
サンプラーにより、探索した箇所において前記物質をサンプリングするステップと
を含む方法。
【請求項１５】
分析装置により、前記サンプラーにサンプリングされた物質を分析するステップをさらに含む
請求項１４に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、外部環境の空間から対象とする物質をサンプリングする支援を行うサンプリング支援装置、サンプリング支援装置を備えたシステム、及びサンプリングを行う方法に関する。
【背景技術】
【０００２】
種々の環境における空間中で汚染物質等の影響を調査する技術が知られている（例えば特許文献１）。
【０００３】
特許文献１は、調査者が歩行しながら対象地の放射線分布を計測する放射線分布計測システムを開示している。放射線分布計測システムは、調査地点の放射線の物理量を放射線計測値として取得する放射線計測手段と、調査地点の位置情報を取得する位置情報取得手段と、当該調査地点の放射線計測値を表示する表示手段とを備えている。調査者は、放射線計測手段、位置情報取得手段、及び表示手段を携行しながら、複数の調査地点の放射線計測値を取得することで、計測対象地の放射線分布図を得る。引用文献１の放射線分布計測システムによると、住宅や学校において、欠測箇所を残すことなく網羅的に放射線量を調査するという目的が達成されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１５−１７５８０３号公報
【非特許文献】
【０００５】
【非特許文献１】Kota Shiba, Genki Yoshikawa, "Aero-Thermo-Dynamic Mass Analysis", Scientific Reports volume 6, Article number: 28849 (2016)
【非特許文献２】J. Engel, J. Stueckler, D. Cremers, "Large-Scale Direct SLAM with Stereo Cameras", In International Conference on Intelligent Robots and Systems (IROS), 2015
【非特許文献３】Michael A. Osborne, et al., "Gaussian Processes for Global Optimization", 3rd International Conference on Learning and Intelligent Optimization (LION3), 2009
【発明の概要】
【発明が解決しようとする課題】
【０００６】
環境中においては物質のサンプリングを行い、サンプリング結果を別途、分析装置を用いて詳細に分析するという調査のニーズが存在する。本発明者は、このような調査ではサンプリングに適した場所の探索が重要であることに着目し、サンプリングを支援する技術について鋭意検討を重ねた。
【０００７】
例えば、サンプリングされた物質の濃度が小さすぎると、前記分析装置の検出限界を下回るような事態が懸念される。この例では、サンプリングに適した場所は、環境空間における対象物質の空間分布において、対象物質の濃度が極大（または最大）である地点、またはその近傍である。
【０００８】
以上に鑑みて、本発明は、物質のサンプリングに適した場所におけるサンプリングを支援することができるサンプリング支援装置、システムおよび方法を提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明の一態様に係るサンプリング支援装置は、外部環境の空間においてサンプリングの対象とする物質に関する情報を出力する装置である。サンプリング支援装置は、物質センサと、位置センサと、制御部と、情報出力部とを備える。物質センサは、物質の濃度に応じたセンサ値を測定する。位置センサは、空間中でセンサ値が測定された位置を測定する。制御部は、測定されたセンサ値及び位置に基づいて、空間中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出する。情報出力部は、算出された推定分布に基づいて、物質センサが空間中で次にセンサ値を測定する箇所に誘導する誘導情報を出力する。
【００１０】
本発明の一態様に係るシステムは、上記のサンプリング支援装置と、物質をサンプリングするサンプラーとを備える。
【００１１】
本発明の一態様に係る方法は、外部環境の空間から対象とする物質をサンプリングする方法である。本方法は、物質センサにより、物質の濃度に応じたセンサ値を測定するステップと、位置センサにより、空間中でセンサ値が測定された位置を測定するステップとを含む。本方法は、制御部により、測定されたセンサ値及び位置に基づいて、空間中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出するステップを含む。本方法は、情報出力部により、算出された推定分布に基づいて、物質センサが空間中で次にセンサ値を測定する箇所に誘導する誘導情報を出力するステップを含む。本方法は、誘導情報に応じて、センサ値の測定及び位置の測定を繰り返し、空間中で物質をサンプリングする箇所を探索するステップと、サンプラーにより、探索した箇所において物質をサンプリングするステップとを含む。
【発明の効果】
【００１２】
本発明に係るサンプリング支援装置、システムおよび方法によると、物質のサンプリングに適した場所におけるサンプリングを支援することができる。
【図面の簡単な説明】
【００１３】
【図１】実施形態１に係るサンプリング支援装置の概要を説明するための図
【図２】実施形態１に係るサンプリング支援装置の構成を例示するブロック図
【図３】実施形態１に係る分析方法の手順を示すフローチャート
【図４】サンプリング支援装置の動作の概要を説明するための図
【図５】実施形態１に係るサンプリング支援装置の処理を例示するフローチャート
【図６】実施形態１における推定分布と価値関数を例示する図
【図７】ガウス過程回帰による推定分布を説明するための図
【図８】サンプリング支援装置の表示部における誘導情報の表示例を示す図
【図９】サンプリング支援装置の表示部におけるサンプリング誘導情報の表示例を示す図
【図１０】実施形態１における推定分布の更新を説明するための図
【図１１】実施形態２に係るサンプリング支援装置の処理を例示するフローチャート
【図１２】実施形態２における推定分布の経験則情報を説明するための図
【図１３】実施形態２に係るサンプリング支援装置の動作を説明するための図
【図１４】サンプリング支援装置による表示の変形例１を示す図
【図１５】サンプリング支援装置による表示の変形例２を示す図
【図１６】サンプリング支援装置による表示の変形例３を示す図
【図１７】サンプリング支援装置のシステムの変形例を示す図
【図１８】サンプリング支援装置のシステムの変形例の動作を例示するフローチャート
【発明を実施するための形態】
【００１４】
以下、添付の図面を参照して本発明に係るサンプリング支援装置の実施の形態を説明する。なお、以下の各実施形態において、同様の構成要素については同一の符号を付している。
【００１５】
（実施形態１）
実施形態１に係るサンプリング支援装置の構成および動作について、以下説明する。
【００１６】
１．構成
１−１．概要
実施形態１に係るサンプリング支援装置の概要について、図１を用いて説明する。
【００１７】
図１は、本実施形態に係るサンプリング支援装置１の概要を説明するための図である。図１（ａ）は、サンプリング支援装置１の使用状態を例示する。図１（ｂ）は、サンプリング支援装置１を用いたサンプリング結果の分析装置４を例示する。
【００１８】
本実施形態に係るサンプリング支援装置１は、図１（ａ）に例示するように、ユーザ２が所望の環境空間３において分析対象とする対象物質のサンプリングを行う際に利用される。対象物質は、例えばＶＯＣ（揮発性有機化合物）等の空気中に分布し得る汚染物質又は化学物質である。環境空間３は、対象物質が存在し得ると想定される環境の空間であり、例えば住宅の内部である。住宅等においては、壁材といった建材がＶＯＣ等の汚染物質及び異臭の発生元となり得る。
【００１９】
図１（ａ）の例において、ユーザ２は、サンプリング支援装置１と共にシステム２０を構成するサンプラー２１を用いて、環境空間３における空気を捕集することにより、対象物質をサンプリングする。本実施形態では、環境空間３のサンプリング結果について、図１（ｂ）に示すようにシステム２０とは別途、提供される分析装置４において、対象物質の詳細な分析が行われる。分析装置４は、例えばＧＣ／ＭＳ（ガスマトグラフ質量分析計）等である。
【００２０】
住宅等の環境空間３においては、換気による気流等から、対象物質が偏って存在することが考えられる。図１（ａ）に例示する環境空間３では、対象物質のホットスポット３０が生じている。ホットスポット３０は、環境空間３における対象物質の空間分布において、対象物質の濃度が極大である地点である。
【００２１】
ここで、対象物質のサンプリングがホットスポット３０とはかけ離れた場所で行われると、サンプリング結果において対象物質の濃度が薄くなり過ぎ、図１（ｂ）の分析時に、分析装置４による対象物質の検出限界を下回るような事態が懸念される。これに対して、環境空間３中の各々の地点のＶＯＣを測定することも考えられるが、ＧＣ／ＭＳ等によるＶＯＣ測定を空間中に網羅的に実施することは、現実的なコストでは困難である。
【００２２】
そこで、本実施形態では、例えばホットスポット３０のように環境空間３において対象物質のサンプリングに適した箇所を検出して、ユーザ２を誘導することにより、サンプリングを支援するサンプリング支援装置１を提供する。
【００２３】
１−２．サンプリング支援装置の構成
本実施形態に係るサンプリング支援装置１の構成について、図２を用いて説明する。図２は、サンプリング支援装置１の構成を例示するブロック図である。
【００２４】
サンプリング支援装置１は、例えば図２に示すように、物質センサ１０と、カメラ１１と、制御部１２と、記憶部１３と、表示部１４と、操作部１５と、スピーカ１６と、振動デバイス１７とを備える。サンプリング支援装置１は、例えばスマートフォンなどのモバイル端末を用いて構成できる。
【００２５】
物質センサ１０は、外部の空間における対象物質の濃度に相関を有する各種のセンサ値を測定可能な種々のセンサデバイスで構成できる。例えば、物質センサ１０は、非特許文献１のＡＭＡといった混合気体の平均モル質量、すなわち密度に相当するセンサ値を測定する密度センサであってもよい。この場合、物質センサ１０は、気体密度に応じて、対象物質以外の物質にも感度を有し得る。このように、物質センサ１０は、対象物質ではない物質を含めた感度に基づき、センサ値を測定してもよい。
【００２６】
物質センサ１０は、例えば２秒以下などのスループットにおいてセンサ値を測定する。これにより、ユーザ２が物質センサ１０を位置決めしながら各位置のセンサ値を測定し易い。図１の例において、物質センサ１０は、サンプリング支援装置１の筐体から突出するように棒材等に取り付けられている。
【００２７】
物質センサ１０は、モル質量に基づく密度センサに限らず、例えば分光光度計、吸光光度計、屈折度計およびガスセンサなどであってもよい。また、サンプリングする空間を満たす流体が液体の場合、物質センサ１０は、ｐＨ計又は電気伝導度計であってもよい。
【００２８】
カメラ１１は、例えばＲＧＢ−Ｄカメラ又はステレオカメラであり、自己位置推定（ＳＬＡＭ）に適用可能な撮像画像を生成する。カメラ１１は、公知のＳＬＡＭ技術（例えば非特許文献２）を実現するソフトウェアとの協働により、本実施形態におけるサンプリング支援装置１の位置センサとして機能できる。カメラ１１は、例えばサンプリング支援装置１に組み込まれることにより、物質センサ１０までの距離が固定される。
【００２９】
制御部１２は、例えばソフトウェアと協働して所定の機能を実現するＣＰＵ又はＭＰＵを含み、サンプリング支援装置１の全体動作を制御する。制御部１２は、記憶部１３に格納されたデータ及びプログラムを読み出して種々の演算処理を行い、各種の機能を実現する。例えば、制御部１２は、カメラ１１の撮像画像等に基づき、自己位置推定を行う。制御部１２が実行するプログラムは、通信ネットワーク等から提供されてもよいし、可搬性を有する記録媒体に格納されていてもよい。
【００３０】
なお、制御部１２は、所定の機能を実現するように設計された専用の電子回路又は再構成可能な電子回路などのハードウェア回路であってもよい。制御部１２は、ＣＰＵ、ＭＰＵ、ＧＰＵ、マイコン、ＤＳＰ、ＦＰＧＡ及びＡＳＩＣ等の種々の半導体集積回路で構成されてもよい。
【００３１】
記憶部１３は、サンプリング支援装置１の機能を実現するために必要なプログラム及びデータを記憶する記憶媒体である。記憶部１３は、例えばフラッシュメモリ、ＳＳＤ又はＨＤＤで構成され、所定の機能を実現するためのパラメータ、データ及び制御プログラム等を格納する。記憶部１３は、例えばＤＲＡＭ又はＳＲＡＭ等のＲＡＭであってもよく、制御部１２の作業エリアとして機能してもよい。例えば、記憶部１３は、物質センサ１０の測定結果、およびサンプリング支援装置１に関する位置を示す位置情報などを記録する。記憶部１３は、環境空間３を示すマップ情報などを記憶してもよい。
【００３２】
表示部１４は、例えば液晶ディスプレイ又は有機ＥＬディスプレイで構成される。表示部１４は、各種の情報を表示する情報出力部の一例である。例えば、本実施形態の表示部１４は、制御部１２の制御により、カメラ１１が撮像した撮像画像に重畳させて、ＭＲ（複合現実）窓型の情報表示を行う。
【００３３】
操作部１５は、ユーザ２が操作を行うユーザインタフェースである。例えば、操作部１５が表示部１４と共にタッチパネルを構成する。操作部１５はタッチパネルに限らず、例えば、キーボード、タッチパッド、ボタン及びスイッチ等であってもよい。
【００３４】
スピーカ１６は、音声により情報出力を行う出力デバイスである。振動デバイス１７は、振動により情報出力を行う出力デバイスである。スピーカ１６及び振動デバイス１７は、それぞれサンプリング支援装置１における情報出力部の一例である。サンプリング支援装置１の情報出力部は、光刺激による情報出力を行う各種プロジェクタ等であってもよい。
【００３５】
１−３．サンプラーについて
以上のようなサンプリング支援装置１は、図１の例において、サンプラー２１と共にシステム２０を構成している。サンプラー２１は、対象物質をサンプリングするために気体等の流体を捕集する装置である。サンプラー２１は、流体を吸入する吸入ダクト２２と、吸入した流体が捕集される捕集バッグ２３とを備える。
【００３６】
本システム２０において、サンプリング支援装置１とサンプラー２１とは、有線又は無線通信により通信接続されてもよい。例えば、サンプリング支援装置１の制御部１２は、サンプラー２１の吸入ダクト２２を制御してもよい。サンプラー２１の吸入ダクト２２における吸入口は、例えば物質センサ１０の近傍に配置されてもよい。なお、サンプリング支援装置１とサンプラー２１とは特にシステム２０を構成しなくてもよく、別個に提供可能である。
【００３７】
２．動作
以上のように構成されるサンプリング支援装置１の動作、及びサンプリング支援装置１を用いた分析方法について、以下説明する。
【００３８】
２−１．分析方法について
サンプリング支援装置１を用いてユーザ２が行う分析方法について、図１，３，４を用いて説明する。
【００３９】
図３は、本実施形態に係る分析方法の手順を示すフローチャートである。図４は、サンプリング支援装置１の動作の概要を説明するための図である。
【００４０】
まず、ユーザ２は、サンプリング支援装置１を用いて、環境空間３（図１（ａ））において対象物質のサンプリングを行う場所を探索する（Ｓ１）。図１（ａ）の例において、ユーザ２は、環境空間３中を移動しながら、移動した位置毎に物質センサ１０によるセンサ値の測定を繰り返す。図４に、ステップＳ１の測定結果を例示する。
【００４１】
図４では、環境空間３におけるセンサ値の分布の真値を示す等高線と、複数の測定データＤ１とを例示している。測定データＤ１は、物質センサ１０のセンサ値の測定結果と、センサ値が測定された位置とを関連付けたデータである。図４は、説明の便宜上、二次元で環境空間３における空間位置ｘ，ｙを示している。
【００４２】
ステップＳ１において、サンプリング支援装置１は、得られた測定データＤ１に基づいて、環境空間３において未測定の範囲を含めたセンサ値の分布を推定する処理を行う。サンプリング支援装置１は、例えば図４に示すように、より大きいセンサ値が見込まれる位置ｐ１を推定結果として算出し、算出した位置ｐ１にユーザ２を誘導するための情報提示を行う。ステップＳ１におけるサンプリング支援装置１の動作の詳細については後述する。
【００４３】
次に、ユーザ２は、サンプリング支援装置１から提示された情報を参照して、環境空間３中の特定の箇所において対象物質のサンプリングを行う（Ｓ２）。対象物質のサンプリングは、例えばサンプラー２１の吸入ダクト２２により環境空間３中の特定の箇所の空気を吸入して、捕集バッグ２３に捕集することによって行われる。
【００４４】
例えばユーザ２が複数の捕集バッグ２３を所持している場合、ユーザ２は環境空間３における複数の箇所からサンプリングを行える（Ｓ３）。この場合、ユーザ２は、所望の回数分、ステップＳ１，Ｓ２を繰り返す（Ｓ３でＮＯ）。
【００４５】
ユーザ２は、対象物質のサンプリングを完了すると（Ｓ３でＹＥＳ）、分析装置４（図１（ｂ））を用いて、サンプリング結果に対する対象物質の分析を行う（Ｓ４）。ユーザ２は、例えば分析装置４（図１（ｂ））が設置された場所までサンプリング結果の捕集バッグ２３を持ち帰り、ステップＳ４を実施する。
【００４６】
ステップＳ４では、分析装置４によって、例えば捕集された気体中における対象物質の濃度を精度良く計測したり、複数種類の対象物質の中で捕集された気体中に含まれた対象物質の種類を特定したりすることができる。対象物質の分析が為されることにより、本フローチャートによる手順は終了する。
【００４７】
以上の分析方法によると、サンプリング支援装置１によって提示される情報に基づき（Ｓ１）、環境空間３中で物質のサンプリングに適した場所でサンプリングを行える（Ｓ２）。このため、ステップＳ４における対象物質の分析の精度を良くすることができる。
【００４８】
さらに、サンプリング支援装置１によると、図４に示すように、環境空間３中で得られたセンサ値の測定データＤ１が疎な状態で、センサ値がより高いと見込まれる位置ｐ１がユーザ２に提示される。これにより、環境空間３におけるセンサ値の測定およびサンプリングといった調査を効率良く行うことができる。以下、サンプリング支援装置１の動作の詳細を説明する。
【００４９】
２−２．サンプリング支援装置の処理
図３のステップＳ１におけるサンプリング支援装置１の処理について、図５〜１０を用いて説明する。本実施形態のサンプリング支援装置１は、ベイズ最適化を適用した推定の処理を実行する。以下では、サンプリング支援装置１が、物質センサ１０による測定の誘導と、サンプリングの誘導とをそれぞれ行う処理例を説明する。
【００５０】
図５は、本実施形態に係るサンプリング支援装置１の処理を例示するフローチャートである。図５のフローチャートは、例えば図３のステップＳ１においてユーザ２がサンプリング支援装置１を操作することにより、開始する。本フローチャートによる各処理は、サンプリング支援装置１の制御部１２によって実行される。
【００５１】
図５のフローチャートにおいて、サンプリング支援装置１の制御部１２は、物質センサ１０から測定結果のセンサ値を取得する（Ｓ１１）。物質センサ１０は、例えば所定の時間間隔で逐次、現在の測定結果のセンサ値を制御部１２に出力する。ユーザ２の移動に応じて、環境空間３中の種々の位置におけるセンサ値を得ることができる。
【００５２】
また、制御部１２は、例えば位置センサとして機能するカメラ１１による撮像画像に基づいて、物質センサ１０がセンサ値を測定した測定位置を取得する（Ｓ１２）。例えば、制御部１２は、撮像画像に基づいて環境空間３における自己位置推定を行い、撮像画像の中の物質センサ１０の位置と、予め設定されたカメラ１１から物質センサ１０までの距離とに基づいて、物質センサ１０の測定位置を算出する。なお、ステップＳ１１，Ｓ１２の処理の順序は特に限定されず、例えば並列に実行されてもよい。
【００５３】
制御部１２は、取得したセンサ値と測定位置とを関連付けて１点分の測定データＤ１として記憶部１３に記録し、Ｎ点分の測定データＤ１が新たに取得されたか否かを判断する（Ｓ１３）。Ｎは、後段の演算効率の観点から予め設定される１以上の整数である。制御部１２は、Ｎ点分の測定データＤ１が取得されるまで（Ｓ１３でＮＯ）、ステップＳ１１〜Ｓ１３の処理を繰り返す。
【００５４】
Ｎ点分の測定データＤ１が新たに取得されると（Ｓ１３でＹＥＳ）、制御部１２は、現在までに取得された測定データＤ１に基づいて、環境空間３におけるセンサ値の空間分布を推定するための演算処理を実行する（Ｓ１４）。２以上のＮによると、ステップＳ１４がバッチ処理として実行される。
【００５５】
本実施形態において、制御部１２は、ステップＳ１４の演算処理においてベイズ最適化に従ってガウス過程の回帰を演算することにより、推定分布を算出する。推定分布は、環境空間３においてセンサ値が分布すると推定される空間分布を示す。ステップＳ１４における推定分布の一例を図６（ａ）に示す。
【００５６】
図６（ａ）は、取得された測定データＤ１に基づいて、実線で示した真値の空間分布を推定する推定分布を例示している。図６（ａ）では、説明の簡単化のために１次元の空間位置ｘを例示しているが、空間位置は適宜、２次元又は３次元に拡張可能である。図６（ａ）の横軸は空間位置ｘを示し、縦軸はセンサ値を示す。図６の例では、ホットスポット３０（図１（ａ））に対応する真値の極大点の一つが、ｘ＝２近傍に存在している。
【００５７】
図６（ａ）の例の推定分布は、現在までの測定データＤ１に基づく予測曲線Ｃ１と、予測曲線Ｃ１の上側領域Ｒ１ａおよび下側領域Ｒ１ｂとを含んでいる。予測曲線Ｃ１は、測定データＤ１に基づく推定結果における空間位置ｘ毎のセンサ値の期待値を示す。上側領域Ｒ１ａは、センサ値が予測曲線Ｃ１を上回る可能性が推定される領域であり、空間位置ｘ毎にセンサ値の幅Ｗ１ａを有する。下側領域Ｒ１ｂは、センサ値が予測曲線Ｃ１を下回る可能性が推定される領域であり、上側領域Ｒ１ａと同様に幅Ｗ１ｂを有する。
【００５８】
以上のような推定分布によると、現在までの測定データＤ１に基づき未測定の範囲までセンサ値を推定する予測曲線Ｃ１と共に、上側領域Ｒ１ａの幅Ｗ１ａと下側領域Ｒ１ｂの幅Ｗ１ｂの分の変動幅Ｗ１が得られる。変動幅Ｗ１は、各空間位置ｘにおけるセンサ値の真値が、予測曲線Ｃ１の期待値からずれ得ることが推定される範囲を示す。
【００５９】
ガウス過程回帰による推定分布は、空間位置ｘ毎のガウス分布により規定される。図７に、推定分布における特定の空間位置ｘについてのガウス分布を例示する。図７の横軸はセンサ値を示し、縦軸は確率密度を示す。
【００６０】
推定分布は、空間位置ｘ毎のガウス分布の平均値μ（ｘ）及び分散σ（ｘ）によって表される。具体的に、空間位置ｘ毎の予測曲線Ｃ１の値Ｃ１（ｘ）は、ガウス分布の平均値μ（ｘ）で表される。また、上側領域Ｒ１ａの幅Ｗ１ａは、例えばガウス分布における５０〜９０パーセンタイルに設定される。下側領域Ｒ１ｂの幅Ｗ１ｂは、例えばガウス分布における１０〜５０パーセンタイルに設定される。なお、各領域Ｒ１ａ，Ｒ１ｂの幅Ｗ１ａ，Ｗ１ｂは、別々に設定可能である。例えば、下側領域Ｒ１ｂの幅Ｗ１ｂが省略されてもよい。
【００６１】
ステップＳ１４におけるガウス過程回帰の演算は、カーネル法など公知の手法を適用して実行することができる（例えば非特許文献３）。この際、事前分布としては、例えばμ＝０．５，σ＝１等の初期値を適宜、設定してもよいし、経験に基づき適宜、設定変更が為されてもよい。また、カーネルとしては例えばガウスカーネル採用できる。
【００６２】
図５に戻り、制御部１２は、算出した推定分布に基づいて、環境空間３においてサンプラー２１でサンプリングすべき候補の箇所、即ちサンプリング候補が確定したか否かを判断する（Ｓ１５）。ステップＳ１５の判断は、例えば推定分布（図６（ａ））において、上側領域Ｒ１ａの上端が極大（又は最大）となる点ｍ１の変動幅Ｗ１が、所定のしきい値以下であるか否かに基づき行われる。当該しきい値は、例えば、推定分布における点ｍ１が、実際のセンサ値の極大点となる可能性を充分に見込める程度に小さい変動幅Ｗ１を考慮して設定される。
【００６３】
例えば、制御部１２は、サンプリング候補が確定していないと判断した場合に（Ｓ１５でＮＯ）、物質センサ１０を用いて次にセンサ値を測定するべき位置を決定する（Ｓ１６）。制御部１２は、算出した推定分布に基づき価値関数を演算することにより、ステップＳ１６の処理を実行する。本実施形態の価値関数は、ベイズ最適化におけるアクイジョン関数である。図６（ａ）の推定分布に基づく価値関数の一例を図６（ｂ）に示す。
【００６４】
価値関数は、環境空間３において物質センサ１０が新たにセンサ値を測定する価値が高い空間位置ｘほど、大きい値を有するように設定された関数である。図６（ｂ）の例の価値関数は、空間位置ｘ毎に予測曲線Ｃ１に係数０．２を乗じて、上側領域Ｒ１ａの幅Ｗ１ａと加算することによって得られる。このような価値関数によると、図６（ａ），（ｂ）に示すように、推定分布の極大点ｍ１の近傍に、価値関数が最大値となる空間位置ｐ１が生じることとなる。制御部１２は、例えば空間位置ｐ１を、次に測定すべき位置として決定する（Ｓ１６）。
【００６５】
制御部１２は、決定した位置においてユーザ２が物質センサ１０で次にセンサ値を測定することを誘導する情報である誘導情報を生成し、生成した誘導情報を表示するように表示部１４を制御する（Ｓ１７）。ステップＳ１７における表示部１４の表示例を図８に示す。
【００６６】
図８の表示例において、表示部１４は、背景画像５０に重畳して、測定結果アイコン５１と、測定目標アイコン５２と、誘導アイコン５３とを表示している。背景画像５０は、各種誘導情報の背景として環境空間３を示す画像である。本例の背景画像５０は、カメラ１１による撮像画像であり、物質センサ１０が映っている。
【００６７】
測定結果アイコン５１は、各点の測定データＤ１に対応して、センサ値の大きさおよび測定位置を示す。本例において、測定結果アイコン５１は、キューブ状を有し、背景画像５０において測定位置に対応する位置に配置される。例えば、測定結果アイコン５１の色またはサイズが、測定データＤ１のセンサ値の大きさに応じて設定される。
【００６８】
測定目標アイコン５２は、背景画像５０においてステップＳ１６で決定された位置に配置され、物質センサ１０によって次に測定するべき目標の位置を示すことにより、ユーザ２を誘導する誘導情報の一例である。誘導アイコン５３は、背景画像５０中で物質センサ１０から測定目標アイコン５２へ向かう方向を示すことにより、ユーザ２による物質センサ１０の測定を促す誘導情報の一例である。ステップＳ１７においては、測定目標アイコン５２と誘導アイコン５３とのうちの一方が表示されてもよい。また、測定結果アイコン５１の表示が省略されてもよい。
【００６９】
図５に戻り、制御部１２は、以上のように物質センサ１０で次に測定するべき誘導情報を表示部１４に表示させると（Ｓ１７）、ステップＳ１１に戻る。制御部１２は、物質センサ１０による新たな測定結果に基づいて、ステップＳ１１〜Ｓ１５の処理を再度行う。
【００７０】
一方、制御部１２は、サンプリング候補が確定されたと判断すると（Ｓ１５でＹＥＳ）、確定されたサンプリング候補に関するサンプリング誘導情報を生成して、表示部１４に表示させる（Ｓ１８）。サンプリング誘導情報は、サンプラー２１によるサンプリング候補の場所（サンプリングに適した場所の一例）にユーザ２を誘導するための情報である。図９に、ステップＳ１８の表示例を示す。
【００７１】
図９の表示例において、表示部１４は、図８の測定目標アイコン５２及び誘導アイコン５３の代わりに、それぞれサンプリング誘導情報の一例であるサンプリング候補アイコン５４及び誘導メッセージ５５を表示している。サンプリング候補アイコン５４は、背景画像５０において、確定されたサンプリング候補の位置を示す。誘導メッセージ５５は、ユーザ２にサンプリングを促す内容のテキスト情報を含む。
【００７２】
制御部１２は、サンプリング誘導情報を表示部１４に表示させると（Ｓ１８）、図５のフローチャートによる処理を終了する。この際、ユーザ２は、図３のステップＳ１からステップＳ２に進み、サンプリング誘導情報に従ってサンプリングを行うことができる。
【００７３】
以上の処理によると、ユーザ２が環境空間３においてサンプリングの場所を探索する際に（Ｓ１）、サンプリング支援装置１は逐次、得られる測定データＤ１に基づき推定分布の算出（Ｓ１４）を繰り返し、誘導情報を更新する（Ｓ１７，Ｓ１８）。図１０（ａ）〜（ｃ）を用いて、誘導情報及び推定分布の更新について説明する。
【００７４】
図１０（ａ）は、図６（ａ）から更新後の推定分布を例示する。図１０（ｂ）は、図１０（ａ）の推定分布に基づく価値関数を例示する。図１０（ｃ）は、図１０（ｂ）の価値関数に応じて、図１０（ａ）から更新後の推定分布を例示する。
【００７５】
図１０（ａ）の推定分布は、図６（ａ）の推定分布の算出（Ｓ１４）後に再度、ステップＳ１４が実行された結果を示す。この際、ユーザ２は、図６（ｂ）の価値関数に応じて決定された位置ｐ１（Ｓ１６）についての誘導情報に従って、物質センサ１０によるセンサ値の測定を行う（図８参照）。このことから、図１０（ａ）では、位置ｐ１近傍などにおいて図６（ａ）からデータ点Ｄ１が増えている。
【００７６】
図１０（ａ）によると、ガウス過程回帰（Ｓ１４）が再度、演算されることにより、推定分布の極大点ｍ１近傍などで変動幅Ｗ１が小さくなるように、ベイズ更新が為されている。図１０（ａ）の例では、未だサンプリング候補が確定していない場合（Ｓ１５でＮＯ）を想定しており、制御部１２は、図１０（ｂ）の価値関数を演算する（Ｓ１６）。
【００７７】
このとき、例えば図１０（ｂ）の価値関数が最大となる空間位置に測定目標アイコン５２（図８）が位置するように、誘導情報が更新され（Ｓ１７）、当該位置の近傍におけるセンサ値の測定が誘導される。すると、その後のステップＳ１４において、推定分布の極大点ｍ１近傍の変動幅Ｗ１がさらに小さくなる。図１０（ｃ）では、極大点ｍ１の位置ｐ１０がサンプリング候補として確定された状態を例示する。
【００７８】
図１０（ｃ）の推定分布では、極大点ｍ１の位置ｐ１０から離れると、同位置ｐ１０よりも大きい変動幅Ｗ１が生じている。本実施形態のサンプリング支援装置１によると、極大点ｍ１から離れた位置の変動幅Ｗ１の大きさに拘わらず、サンプリング候補の位置ｐ１０を確定することができる。
【００７９】
以上の処理におけるステップＳ１７では、図８に示すように、測定目標アイコン５２および誘導アイコン５３といった誘導情報が、現実の環境空間３を示す背景画像５０に合わせて表示される。また、ステップＳ１８では、図９に示すように、サンプリング候補アイコン５４および誘導メッセージ５５といったサンプリング誘導情報が、背景画像５０に合わせて表示される。これにより、ユーザ２は、現実の環境空間３と各種誘導情報とを同時に目視でき、サンプリング等の誘導の妥当性を判断し易くすることができる。
【００８０】
以上の処理において、サンプリング候補に関するステップＳ１５，Ｓ１８の処理は省略されてもよい。この場合であっても、ステップＳ１７の誘導情報が逐次、更新されることにより、例えば誘導情報が誘導する位置が収束した場合などに、ユーザ２は、当該位置においてサンプリングを行うことを自己判断できる。このように、サンプリング支援装置１は誘導情報の提示により、ユーザ２が行うサンプリングを支援することできる。
【００８１】
また、以上のステップＳ１５の判断は、上記の例に限らず、例えば、下側領域Ｒ１ｂの下端が所定のしきい値以上であるか否かに基づき行われてもよい。当該しきい値は、例えば、推定分布における点ｍ１が、真値の極大点でなかったとしても、分析装置４において充分に検出可能と想定される値に設定される。これにより、分析装置４において検出可能なサンプリング候補にユーザ２を誘導することができる。また、ステップＳ１５の判断に、価値関数が用いられてもよい。価値関数の演算時は、ステップＳ１５でＮＯの場合でなくてもよい。
【００８２】
また、ステップＳ１６において、決定される位置は１つに限らず、複数の位置が、次にセンサ値を測定すべき位置として算出されてもよい。この場合、ステップＳ１７において、例えば複数の位置に測定目標アイコン５２が配置されてもよい。これにより、ユーザ２は、例えば最寄りの測定目標アイコン５２など、センサ値の測定を行う位置を適宜選択することができる。
【００８３】
３．まとめ
以上のように、本実施形態に係るサンプリング支援装置１は、外部環境の空間である環境空間３においてサンプリングの対象とする物質である対象物質に関する情報を出力する。サンプリング支援装置１は、物質センサ１０と、位置センサの一例であるカメラ１１と、制御部１２と、情報出力部の一例である表示部１４とを備える。物質センサ１０は、物質の濃度に応じたセンサ値を測定する。位置センサは、環境空間３中でセンサ値が測定された位置を測定する。制御部１２は、測定されたセンサ値及び位置に基づいて、環境空間３中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出する。情報出力部は、算出された推定分布に基づいて、物質センサ１０が環境空間３中で次にセンサ値を測定する箇所に誘導する誘導情報を出力する。
【００８４】
以上のサンプリング支援装置１によると、センサ値が極大となることが推定される箇所を含んだ推定分布に基づいて、物質センサ１０が環境空間３中で次にセンサ値を測定する箇所に誘導する誘導情報が出力される。これにより、サンプリング支援装置１のユーザ２は、環境空間３において物質のサンプリングに適した場所（例えば、ホットスポットまたはその近傍）を探索する際に、誘導情報によって誘導される箇所を参考にして、物質センサ１０による測定を行える。よって、誘導情報を参照するユーザ２が、環境空間３の探索をし易くなり、物質のサンプリングに適した場所におけるサンプリングを支援することができる。
【００８５】
本実施形態において、制御部１２は、算出した推定分布においてセンサ値が極大となることが推定される箇所に応じて、誘導情報を生成する。これにより、センサ値が極大となるホットスポット３０近傍にユーザ２を誘導して、ホットスポット３０でのサンプリングを実現し易くすることができる。
【００８６】
また、本実施形態において、制御部１２は、物質センサ１０および位置センサから逐次、測定されたセンサ値及び位置を取得して（Ｓ１１，Ｓ１２）、例えば次にセンサ値を測定する箇所をサンプリングに適した場所に近づけるように、誘導情報を更新する（Ｓ１７，Ｓ１８）。これにより、ユーザ２の誘導先が、逐次、更新される毎にサンプリングに適した場所に近づき、ユーザ２がサンプリングに適した場所でサンプリングを行い易くすることができる。
【００８７】
また、本実施形態において、誘導情報の一例である測定目標アイコン５２は、物質センサ１０が環境空間３中で次にセンサ値を測定するべき位置を提示する（Ｓ１６，Ｓ１７）。これにより、測定目標アイコン５２の位置において物質センサ１０による測定をユーザ２に促して、サンプリングを行い易くすることができる。
【００８８】
また、本実施形態において、情報出力部は、環境空間３中で物質をサンプリングする箇所に誘導するサンプリング誘導情報を、更に出力する（Ｓ１８）。サンプリング誘導情報は、例えばサンプリング候補アイコン５４が示すサンプリング候補の位置により、物質のサンプリングに適した場所をユーザ２に提示する。サンプリング誘導情報により、物質センサ１０による測定に加えて、サンプラー２１によるサンプリングをユーザ２に促すことができる。
【００８９】
また、本実施形態において、制御部１２は、算出した推定分布に基づき、環境空間３中の位置毎にセンサ値を測定する価値を示す価値関数を演算し（Ｓ１６）、演算した価値関数に基づき、誘導情報を生成する（Ｓ１７）。制御部１２は、測定されたセンサ値及び位置に基づくガウス過程回帰により、推定分布を算出する。これにより、ベイズ最適化に従って、推定分布を精度良く算出することができる。
【００９０】
また、本実施形態において、情報出力部の一例の表示部１４は、環境空間３を示す背景に重畳させて、各種の誘導情報を表示する。これにより、ユーザ２が誘導情報の妥当性を推察し易くすることができる。
【００９１】
また、本実施形態において、物質センサ１０は、モル質量に基づく密度センサ、分光光度計、吸光光度計、屈折度計、ガスセンサ、ｐＨ計、および電気伝導度計の少なくとも１つを含む。本実施形態のサンプリング支援装置１によると、対象物質以外の物質にも反応し得る物質センサ１０を用いても、ユーザ２によるサンプリングを行い易くすることができる。
【００９２】
また、本実施形態において、サンプリング支援装置１と、対象物質をサンプリングするサンプラー２１とを備えるシステム２０が提供される。本システム２０では、サンプリング支援装置１の支援により、物質のサンプリングに適した場所でサンプリングを行い易くすることができる。
【００９３】
また、本実施形態において、環境空間３から対象とする物質をサンプリングする方法が提供される。本方法は、物質センサ１０により、物質の濃度に応じたセンサ値を測定するステップＳ１と、位置センサにより、空間中でセンサ値が測定された位置を測定するステップＳ２とを含む。本方法は、制御部１２により、測定されたセンサ値及び位置に基づいて、空間中でセンサ値が極大となることが推定される箇所を含んだ推定分布を算出するステップＳ１４を含む。本方法は、情報出力部により、算出された推定分布に基づいて、物質センサ１０が空間中で次にセンサ値を測定する箇所に誘導する誘導情報を出力するステップＳ１７を含む。本方法は、誘導情報に応じて、センサ値の測定及び位置の測定を繰り返し、空間中で物質をサンプリングする箇所を探索するステップＳ１と、サンプラー２１により、探索した箇所において物質をサンプリングするステップＳ２とを含む。
【００９４】
本方法によると、誘導情報を参照するユーザ２等が、環境空間３において物質のサンプリングに適した場所でサンプリングを行い易くすることができる。
【００９５】
本実施形態の分析方法は、上記各ステップに加えて、分析装置４により、サンプラー２１にサンプリングされた物質を分析するステップＳ４をさらに含む。本分析方法によると、物質のサンプリングに適した場所のサンプリング結果により、分析の精度を良くすることができる。
【００９６】
（実施形態２）
以下、図１１〜１３を参照して、実施形態２を説明する。実施形態１では、ベイズ最適化に基づく演算処理により、センサ値の推定分布が算出された。本実施形態では、経験則に基づく簡易な演算処理によって推定分布を算出するサンプリング支援装置１について説明する。
【００９７】
図１１は、実施形態２に係るサンプリング支援装置１の処理を例示するフローチャートである。実施形態１では、図５のステップＳ１４において、推定分布の算出にガウス過程回帰が採用された。これに代えて、本実施形態におけるサンプリング支援装置１の制御部１２は、経験則情報を用いて推定分布を算出する（Ｓ１４Ａ）。経験則情報は、センサ値が測定位置から離れるとどの程度変動し得るかの経験則を示す情報であり、例えば予め記憶部１３に格納される。経験則情報の一例を図１２に示す。
【００９８】
図１２（ａ）は、推定分布における上側領域Ｒ２ａについての経験則情報Ｄ２ａを例示する。図１２（ｂ）は、下側領域Ｒ２ｂについての経験則情報Ｄ２ｂを例示する。図１２（ａ），（ｂ）において、横軸は現在位置からの距離Δｘ（≧０）を示し、縦軸は各領域Ｒ２ａ，Ｒ２ｂの幅Ｗ２ａ，Ｗ２ｂの変動量を示す。経験則情報Ｄ２ａ，Ｄ２ｂは、本実施形態における推定情報の一例である。
【００９９】
図１２（ａ）の経験則情報Ｄ２ａは、センサ値の測定位置から距離Δｘが離れるほど、センサ値が増大し得ると経験的に考えられる傾向を示している。図１２（ｂ）の経験則情報Ｄ２ｂは、センサ値の測定位置から距離Δｘが離れるほど、センサ値が減少し得る傾向を示している。各経験則情報Ｄ２ａ，Ｄ２ｂは、例えばルックアップテーブルとして記憶部１３に格納できる。
【０１００】
図１１のステップＳ１４Ａにおいて、制御部１２は、以上のような経験則情報Ｄ２ａ，Ｄ２ｂにより、２点の測定データＤ１間に三角近似を適用して、推定分布を算出する。三角近似は、２つの測定位置の間の上側領域Ｒ２ａ（又は下側領域Ｒ２ｂ）を、例えば以下のように、略三角形状に近似する近似法である。三角近似に基づく本実施形態の推定分布の一例を図１３（ａ）に示す。
【０１０１】
図１３（ａ）の例において、予測曲線Ｃ２は、測定データＤ１間を線形補間することによって算出される。例えば、制御部１２は、各点の測定テータＤ１からの距離に応じて、図１２（ａ）の経験則情報Ｄ２ａを参照して、予測曲線Ｃ２の上側領域Ｒ２ａにおける幅Ｗ２ａを設定する。同様に、制御部１２は、図１２（ｂ）の経験則情報Ｄ２ｂを参照して、予測曲線Ｃ２の下側領域Ｒ２ｂにおける幅Ｗ２ｂを設定する。
【０１０２】
以上のような三角近似によって、図１３（ａ）に示すように推定分布における変動幅Ｗ２が算出される。図１３（ａ）の推定分布に基づく価値関数を図１３（ｂ）に例示する。
【０１０３】
図１３（ｂ）の価値関数は、図１３（ａ）の推定分布の予測曲線Ｃ１および幅Ｗ２ａに基づいて、実施形態１の例（図６（ｂ））と同様に算出される。このような価値関数によっても、図１３（ｂ）に示すように、最大値となる空間位置ｐ２が、センサ値の真値の極大点（ｘ＝２）の近傍に生じ得る。制御部１２は、実施形態１と同様に、センサ値を次に測定すべき位置として当該位置ｐ２を決定し（図１１のＳ１６）、同位置ｐ２に誘導する誘導情報を出力する（Ｓ１７）。
【０１０４】
図１３（ｃ）は、図１３（ａ）から更新後の推定分布を例示する。本実施形態の推定分布によっても、図１３（ｃ）に例示するように、上側領域Ｒ２ａの上端が極大となる点ｍ２の空間位置２０は、真値の極大点（ｘ＝２）近傍に収束し得る。よって、当該位置２０に誘導するサンプリング誘導情報の出力により（Ｓ１８）、ユーザ２のサンプリングを支援することができる。
【０１０５】
以上のように、本実施形態において、サンプリング支援装置１の記憶部１３は、経験則情報Ｄ２ａ，Ｄ２ｂを記憶する。経験則情報Ｄ２ａ，Ｄ２ｂは、センサ値が測定された位置からの距離Δｘに応じて推定されるセンサ値の変動幅Ｗ２を規定する推定情報の一例である。制御部１２は、経験則情報Ｄ２ａ，Ｄ２ｂを参照して、推定分布を算出する（Ｓ１４Ａ）。これにより、例えば測定されたセンサ値及び位置の間に三角近似を適用するような、簡易な演算処理によって推定分布を算出し、サンプリング支援装置１の処理負荷を低減することができる。サンプリング支援装置１による推定分布の算出には、以上のように種々の手法を採用することができる。
【０１０６】
（他の実施形態）
上記の実施形態１では、ＭＲ窓型の表示部１４による表示例を説明した（図８，９）。本実施形態のサンプリング支援装置１において、表示部１４は、特にＭＲに限らず、例えばＡＲ（拡張現実）あるいはＶＲ（仮想現実）による表示を行ってもよい。このような変形例について、図１４〜１６を用いて説明する。
【０１０７】
図１４は、サンプリング支援装置１による表示の変形例１を示す。本変形例において、サンプリング支援装置１の表示部１４は、背景画像５０に重畳する誘導情報として、分布画像５６を表示している。分布画像５６は、価値関数が示す価値の空間分布を示す画像である。例えば、サンプリング支援装置１の制御部１２は、図５のステップＳ１６，１７の代わりに、価値関数の演算結果に基づき分布画像５６を生成する。分布画像５６は、各種色彩の濃度分布により構成されてもよいし、点群で構成されてもよい。
【０１０８】
例えば、物質センサ１０が湿った地面の湿度と対象物質の両方に反応する場合、湿度は、地面の高さに相関して滑らかに変動しており、水平方向にはあまり変動しないことが予想される。そのような水平方向には滑らかな濃度変動の中で、急峻な濃度変動があれば、ユーザ２は、何かしら湿度要因以外の物質による分布があることを容易に推察できる。
【０１０９】
以上のように、本変形例は、環境空間３における価値関数の分布を提示する。これにより、例えば対象物質と他の物質が混合したセンサ値であっても、現実の環境空間３と照らし合わせることにより、実際の対象物質の空間分布をユーザ２に推察させることができる。
【０１１０】
図１５は、サンプリング支援装置１による表示の変形例２を示す。サンプリング支援装置１は、環境空間３を俯瞰するように、各種誘導情報を表示してもよい。図１５の例において、サンプリング支援装置１の表示部１４は、環境空間３を俯瞰する背景画像５０Ａに重畳させて、分布画像５６とユーザアイコン５７とを表示している。ユーザアイコン５７は、ユーザ２の現在位置を示す。
【０１１１】
本実施形態のサンプリング支援装置１は、図１５の俯瞰型の表示とＭＲ窓型の表示とが、ユーザ２の操作によって切り替え可能に構成されてもよい。図１５のような背景画像５０Ａは、例えばあらかじめサンプリング支援装置１の記憶部１３に記憶される。或いは、サンプリング支援装置１は、カメラ１１の撮像画像に基づき背景画像５０Ａを画像合成してもよいし、環境空間３の天井などに設置されたカメラから撮像画像を取得してもよい。
【０１１２】
図１６は、サンプリング支援装置１による表示の変形例３を示す。図１６の例において、サンプリング支援装置１の表示部１４は、ＨＭＤ（ヘッドマウントディスプレイ）で構成され、環境空間３を背景として各種誘導情報をユーザ２に提示している。表示部１４としてのＨＭＤは、シースルー型であってもよいし、カメラ合成を用いたビデオシースルー型であってもよい。また、表示部１４は、可搬型プロジェクタであってもよい。
【０１１３】
また、上記の各実施形態では、情報出力部として表示部１４を例示したが、サンプリング支援装置１の情報出力部は、表示部、スピーカ、振動デバイス、及びプロジェクタのうちの少なくとも１つを含んでもよい。サンプリング支援装置１は、例えば音、信号あるいは光刺激などの種々の形式において、各種誘導情報を出力してもよい。各種誘導情報の出力は、サンプリング支援装置１の現在位置から誘導する向きを指し示すように行われてもよい。
【０１１４】
上記の各実施形態では、サンプリング支援装置１のユーザ２が物質センサ１０を動かしながらサンプリングを行う例を説明したが、物質センサ１０の移動あるいはサンプリングは、機械的に行われてもよい。本変形例について、図１７，１８を用いて説明する。
【０１１５】
図１７は、サンプリング支援装置１のシステム２０の変形例を示す図である。本変形例において、サンプリング支援装置１およびサンプラー２１は移動体６に搭載されている。移動体６は、例えば各種ロボット、車両およびドローンなどであってもよい。移動体６は、ユーザ２の操作によって、或いは自動的に移動可能である。本変形例において、サンプリング支援装置１は、例えば移動体６を制御する制御装置６０に通信接続している。また、制御装置６０は、サンプラー２１を制御可能に構成される。
【０１１６】
サンプリング支援装置１は、制御装置６０に組み込んで構成されてもよい。また、サンプリング支援装置１は、例えば物質センサ１０を移動体６に配置して、他の構成要素は移動体６とは別体で提供されてもよい。ユーザ２は、例えば移動体６を操作することにより、環境空間３中で物質センサ１０を間接的に動かして、サンプリングの場所を探索することができる。
【０１１７】
図１８は、本変形例の動作を例示するフローチャートである。図１８のフローチャートによる各処理は、移動体６の制御装置６０によって実行される。具体的に、制御装置６０は、例えば図３と同様に、サンプリング支援装置１を用いてサンプリングの場所を探索し（Ｓ１）、サンプラー２１によって対象物質のサンプリングを行う（Ｓ２）。これによっても、サンプリング支援装置１による推定により、対象物質のサンプリングを行い易くすることができる。
【０１１８】
上記の例では、移動体６の制御装置６０が、サンプラー２１を制御してサンプリングを行った（Ｓ２）。サンプリング支援装置１の制御部１２が、ステップＳ２においてサンプラー２１を制御してもよい。環境空間３の状況によっては、対象物質の濃度分布が時系列で変動しており、可視化してからサンプリングしては遅い場合が考えられる。そこでより好ましくはセンサ値が急激に上昇した場合には自動でサンプラー２１のポンプ及び弁を操作し捕集する処理を行ってもよい。
【０１１９】
上記の各実施形態では、環境空間３の例として住宅内部を例示した。本実施形態のサンプリング支援装置１は、上記の例に限らず種々の外部環境の空間に適用可能であり、例えば、各種建物の内部および外部、各種大気中の環境、海洋、河川において適用可能である。また、本実施形態のサンプリング支援装置１は、火山ガスの調査、各種ガス又は汚染物質の漏洩の検査に適用可能である。また、対象物質は、種々の毒性物質および放射性物質であってもよい。対象物質を含む流体は、気体であってもよいし液体であってもよい。
【０１２０】
本実施形態のサンプリング支援装置１によると、副次的に検出限界に近い物質を測定したい場合にも利用できる。例えば大気や海洋河川における有毒物質の漏洩事故の場合、溶媒は検出できるが毒性物質濃度が検出限界以下であり同定が困難な場合がある。このとき溶媒濃度が高い場所をサンプリングすれば、毒性物質の濃度も比例して高まり、同定が可能である可能性が高まる。
【０１２１】
また、上記の各実施形態において、サンプリング支援装置１の位置センサとしてカメラ１１を例示した。本実施形態のサンプリング支援装置１において、位置センサはカメラ１１に限らず、種々の方式で各種信号を受信することにより、位置測定を行ってもよい。例えば、サンプリングを行う規模および求められる精度によっては、位置センサはＧＰＳ受信機を含んでもよいし、環境空間３に配置された各種マーカからの信号の受信機を含んでもよい。
【符号の説明】
【０１２２】
１    サンプリング支援装置
１０    物質センサ
１１    カメラ
１２    制御部
１３    記憶部
１４    表示部
１５    操作部
１６    スピーカ
１７    振動デバイス
２０    システム
２１    サンプラー
３    環境空間
４    分析装置
【図１】
2020041996.tif 000003
【図２】
2020041996.tif 000004
【図３】
2020041996.tif 000005
【図４】
2020041996.tif 000006
【図５】
2020041996.tif 000007
【図６】
2020041996.tif 000008
【図７】
2020041996.tif 000009
【図８】
2020041996.tif 000010
【図９】
2020041996.tif 000011
【図１０】
2020041996.tif 000012
【図１１】
2020041996.tif 000013
【図１２】
2020041996.tif 000014
【図１３】
2020041996.tif 000015
【図１４】
2020041996.tif 000016
【図１５】
2020041996.tif 000017
【図１６】
2020041996.tif 000018
【図１７】
2020041996.tif 000019
【図１８】
2020041996.tif 000020