JP2022125505A - 計測装置、計測方法、およびプログラム - Google Patents

Abstract

【課題】物品の保管状態を効率的に管理する。【解決手段】本開示の一態様に係る計測装置は、少なくとも１つの対象物品を収容する空間に設置された音波送信装置と、空間に設置された音波受信装置とに接続可能な計測装置であって、音波送信装置に音波を送信させる手段と、音波受信装置による音波の受信結果を参照して空間を仮想的に分割した複数の区画に亘る物理量の分布を算出する手段と、対象物品の位置に対応する区画である保管区画の物理量が対象物品に関連付けられる基準を満たすか否かを判定する手段と、保管区画の物理量が基準を満たさないと判定された場合に、所定の動作を行う手段とを具備する。【選択図】図７

Description

本発明は、計測装置、計測方法、およびプログラムに関する。

物品の品質を維持するために、当該物品の保管に関して何らかの基準が設けられることがある。例えば、食品、飲料、試薬、薬品、ワクチンなどの物品に対して保管温度を定めることがある。

物品が保管される空間（例えば倉庫、コンテナ）内の温度を計測するために、空間内に温度計を設置することができる。しかしながら、温度計は、当該温度計の設置された位置付近の温度を計測するに過ぎない。空間内の温度分布が一様でない場合に、温度計の設置された位置の温度と、物品の設置された位置の温度とは大きく異なるおそれがある。故に、温度計による計測温度を参照するだけでは、物品の保管温度に関する基準が満たされているか否かの判定を誤ることがある。

空間内の全物品に温度計を取り付けることで、各物品の保管温度を正確に計測することが可能となる。特許文献１には、恒温槽に保管されたボトル毎に当該ボトルの温度を検出する温度計を取り付けている。

特開2018-177463号公報

特許文献１に記載の技術を用いる場合、温度管理の対象となる物品が増えるほど、必要な温度計の取り付け数が多くなり、温度計の調達コストおよび取り付けコストが増加する。

本開示の目的は、物品の保管状態を効率的に管理することである。

本開示の一態様に係る計測装置は、少なくとも１つの対象物品を収容する空間に設置された音波送信装置と、空間に設置された音波受信装置とに接続可能な計測装置であって、音波送信装置に音波を送信させる手段と、音波受信装置による音波の受信結果を参照して空間を仮想的に分割した複数の区画に亘る物理量の分布を算出する手段と、対象物品の位置に対応する区画である保管区画の物理量が対象物品に関連付けられる基準を満たすか否かを判定する手段と、保管区画の物理量が基準を満たさないと判定された場合に、所定の動作を行う手段とを具備する。

第１実施形態の保管状態管理システムの構成を示すブロック図である。 第１実施形態の保管状態管理システムの詳細な構成を示すブロック図である。 第１実施形態の音波送信装置の構成を示す概略図である。 第１実施形態の音波受信装置の構成を示す概略図である。 対象空間のメッシュ構造を示す図である。 対象空間における物品の配置例を示す図である。 第１実施形態の概要の説明図である。 第１実施形態の空間データテーブルのデータ構造を示す図である。 第１実施形態のセンサデータテーブルのデータ構造を示す図である。 第１実施形態の経路データテーブルのデータ構造を示す図である。 第１実施形態のメッシュデータテーブルのデータ構造を示す図である。 第１実施形態のフィルタの説明図である。 第１実施形態の物品データテーブルのデータ構造を示す図である。 第１実施形態の保管状態管理処理のフローチャートである。 第１実施形態のセンサ配置の一例を示す図である。 図１４の温度計測処理の詳細なフローチャートである。 図１６の経路温度の計算の詳細なフローチャートである。 図１６の処理において表示される画面例を示す図である。 第２実施形態の概要の説明図である。 第２実施形態の保管状態管理処理のフローチャートである。 変形例１の概要の説明図である。 変形例１の温度計測処理のフローチャートである。 変形例３のセンサ配置の一例を示す図である。 変形例１の音波送信装置及び音波受信装置の構成を示す概略図である。 変形例４の保管状態管理処理のフローチャートである。

以下、本発明の一実施形態について、図面に基づいて詳細に説明する。なお、実施形態を説明するための図面において、同一の構成要素には原則として同一の符号を付し、その繰り返しの説明は省略する。

（１）第１実施形態
第１実施形態について説明する。

（１－１）保管状態管理システムの構成
第１実施形態の保管状態管理システムの構成について説明する。図１は、第１実施形態の保管状態管理システムの構成を示すブロック図である。図２は、第１実施形態の保管状態管理システムの詳細な構成を示すブロック図である。

図１及び図２に示すように、保管状態管理システム１は、計測装置１０と、音波送信装置２０と、音波受信装置３０と、空調装置４０と、温度計５０と、を備える。
計測装置１０は、音波送信装置２０、音波受信装置３０、空調装置４０、及び、温度計５０に接続されている。
計測装置１０は、複数の同種の装置に接続可能である。例えば、計測装置１０は、複数の音波送信装置２０に接続されてもよいし、複数の音波受信装置３０に接続されてもよいし、複数の空調装置４０に接続されてもよい。
計測装置１０、音波送信装置２０、音波受信装置３０、空調装置４０、及び、温度計５０は温度計測の対象となる空間（以下「対象空間」という）ＳＰに配置されている。

対象空間ＳＰは、保管状態を管理されるべき対象物品が保管される空間であり、例えば、以下のいずれかであってよい。
・トラック、船、航空機、または鉄道によって輸送されるコンテナの内部空間
・倉庫の内部空間
保管状態を管理されるべき対象物品は、例えば、薬品、飲料、食品、ワクチン、及び試薬などである。これらの対象物品は、品質を保つために、保管状態（例えば保管時の温度）を管理することが求められる。
なお、本実施形態における対象空間ＳＰ及び対象物品は上記の例に限定されない。

計測装置１０は、以下の機能を備える。
・音波送信装置２０を制御する（例えば、音波送信装置２０に音波を送信させる）機能
・音波受信装置３０から受信波形データを取得する機能
・対象空間ＳＰに亘る温度（「物理量」の一例）の分布を算出する機能
・対象空間ＳＰの特定の区画の物理量が当該区画に設置された物品に関連付けられる基準を満たすか否かを判定する機能
・特定の区画の温度が基準を満たさないと判定した場合に、所定の動作を行う機能
・空調装置４０を制御する機能
・温度計５０から対象空間ＳＰの温度の測定結果に関する基準温度情報を取得する機能
本実施形態における物理量は、空気の状態に関する物理量であり、例えば温度、風速、風量、風向き、又は湿度である。ただし、計測装置１０が算出する物理量は、これらの例に限定されない。

音波送信装置２０は、計測装置１０の制御に従い、指向性を有する音波（例えば、超音波ビーム）を送信するように構成される。また、音波送信装置２０は、超音波の送信方向を変更するように構成される。

音波受信装置３０は、音波送信装置２０から送信された超音波ビームを受信し、且つ、受信した超音波ビームに応じた受信波形データを生成するように構成される。音波受信装置３０は、例えば、無指向性マイクロフォン又は指向性マイクロフォンである。

空調装置４０は、計測装置１０の制御に従い、対象空間ＳＰの温度を調整するように構成される。

温度計５０は、対象空間ＳＰの温度（以下「基準温度」という）を測定するように構成される。温度計５０は、接触式の温度計であってもよいし、非接触式の温度計（例えば、赤外線放射温度計）であってもよい。

（１－１－１）計測装置の構成
第１実施形態の計測装置１０の構成について説明する。

図２に示すように、計測装置１０は、記憶装置１１と、プロセッサ１２と、入出力インタフェース１３と、通信インタフェース１４と、を備える。

記憶装置１１は、プログラム及びデータを記憶するように構成される。記憶装置１１は、例えば、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、及び、ストレージ（例えば、フラッシュメモリ又はハードディスク）の組合せである。

プログラムは、例えば、以下のプログラムを含む。
・ＯＳ（Operating System）のプログラム
・情報処理（例えば、対象空間ＳＰに亘る物理量の分布を算出する情報処理、対象空間ＳＰに収容された物品の保管状態を管理する情報処理）を実行するアプリケーションのプログラム

データは、例えば、以下のデータを含む。
・情報処理において参照されるデータ及びデータベース
・情報処理を実行することによって得られるデータ（つまり、情報処理の実行結果）
・空間の温度に対する音波の速度に関する音波速度特性に関するデータ

プロセッサ１２は、記憶装置１１に記憶されたプログラムを起動することによって、計測装置１０の機能を実現するように構成される。プロセッサ１２は、コンピュータの一例である。記憶装置１１により記憶されるプログラム及びデータは、ネットワークを介して提供されてもよいし、コンピュータにより読み取り可能な記録媒体に記録して提供されてもよい。なお、計測装置１０の機能の少なくとも一部が、１又は複数の専用のハードウェアにより実現されていてもよい。

入出力インタフェース１３は、計測装置１０に接続される入力デバイスから信号（例えば、ユーザの指示、受信波形データ）を取得し、かつ、計測装置１０に接続される出力デバイスに信号（例えば、制御信号、画像信号）を出力するように構成される。
入力デバイスは、例えば、キーボード、ポインティングデバイス、タッチパネル、又は、それらの組合せである。また、入力デバイスは、温度計５０を含む。
出力デバイスは、例えば、ディスプレイである。また、出力デバイスは、音波送信装置２０、音波受信装置３０、および空調装置４０を含む。

通信インタフェース１４は、外部装置（例えば、サーバ）との間の通信を制御するように構成される。

（１－１－２）音波送信装置の構成
第１実施形態の音波送信装置２０の構成を説明する。図３は、第１実施形態の音波送信装置の構成を示す概略図である。

図３Ａに示すように、音波送信装置２０は、複数の超音波振動子（「振動素子」の一例）２１と、制御回路２２と、を備える。

図３Ｂに示すように、制御回路２２は、計測装置１０の制御に従って、複数の超音波振動子２１を振動させる。複数の超音波振動子２１が振動すると、送信面（ＸＹ平面）に対して直交する送信方向（Ｚ軸方向）に向かって、超音波ビームが送信される。

（１－１－３）音波受信装置の構成
第１実施形態の音波受信装置３０の構成を説明する。図４は、第１実施形態の音波受信装置の構成を示す概略図である。

図４に示すように、音波受信装置３０は、超音波振動子３１と、制御回路３２と、を備える。

超音波振動子３１は、音波送信装置２０から送信された超音波ビームを受信すると振動する。

制御回路３２は、超音波振動子３１の振動に応じた受信波形データを生成するように構成される。

（１－２）実施形態の概要
第１実施形態の概要について説明する。図５は、対象空間のメッシュ構造を示す図である。図６は、対象空間における物品の配置例を示す図である。図７は、第１実施形態の概要の説明図である。

図５に示すように、対象空間ＳＰには、音波送信装置２０ａ～２０ｂと、音波受信装置３０ａ～３０ｂと、が配置されている。計測装置１０は、音波送信装置２０及び音波受信装置３０と接続可能である。
対象空間ＳＰは、計測装置１０によって個別に物理量を計測可能な複数のメッシュ（「区画」の一例）Ｍｉ（ｉは引数）に仮想的に分割される。各メッシュＭｉは、２次元領域又は３次元領域である。
対象空間ＳＰにおいて、空間的に連続する複数のメッシュを連結することで、音波送信装置２０から音波受信装置３０に至る経路を構成することができる。例えば、経路Ｐ２００は、音波送信装置２０ａから音波受信装置３０ａに至る経路である。経路Ｐ２０１は、音波送信装置２０ｂから音波受信装置３０ｂに至る経路である。ここで、メッシュＭｔは、経路Ｐ２００に含まれると同時に、経路Ｐ２０１にも含まれる。
計測装置１０は、音波を送信させるように、音波送信装置２０ａ～２０ｂを制御する。
計測装置１０は、音波受信装置３０ａ～３０ｂから、受信された音波の波形に関する受信波形データを取得する。
計測装置１０は、受信波形データに基づいて、メッシュＭｔの温度を計算する。メッシュＭｔの温度は、経路Ｐ２００における音波の伝搬距離及び伝搬時間と、経路Ｐ２０１における音波の伝搬距離及び伝搬時間とに基づいて計算される。メッシュの温度の計算方法の詳細については後述する。
また、計測装置１０は、音波が伝搬する経路を変更して同様の計測を行うことで、他のメッシュの温度を計算する。すなわち、計測装置１０は、対象空間ＳＰの複数の区画に亘る温度分布（例えば、各メッシュの温度）を計算する。

このように、計測装置１０は、対象空間ＳＰの区画単位で温度を算出できる。つまり、図６に示すように、メッシュＴＭ１に物品ＴＩ１（「対象物品」の一例）が設置され、かつメッシュＴＭ２に物品ＴＩ２（「対象物品」の一例）が設置されている場合に、計測装置１０は、メッシュＴＭ１の温度と、メッシュＴＭ２の温度とを個別に算出することができる。メッシュＴＭ１の温度は物品ＴＩ１の保管温度に相当し、メッシュＴＭ２の温度は物品ＴＩ２の保管温度に相当する。

故に、計測装置１０は、メッシュＴＭ１の温度を参照し、物品ＴＩ１に関連付けられる基準（例えば、保管温度が閾値Ｔｈ１を超過しない）が満たされているか否かを判定することができる。計測装置１０は、メッシュＴＭ２の温度を参照し、物品ＴＩ２に関連付けられる基準（例えば、保管温度が閾値Ｔｈ２を超過しない）が満たされているか否かを判定することができる。図７の例では、時刻ｔにメッシュＴＭ１の温度が閾値Ｔｈ１を超過した。故に、計測装置１０は、メッシュＴＭ１の温度が物品ＴＩ１に関連付けられる基準を満たさないと判定する。

計測装置１０は、メッシュＴＭ１の温度が物品ＴＩ１に関連付けられる基準を満たさないと判定した場合に、所定の動作を行う。一例として、計測装置１０は、物品ＴＩ１の基準に関するアラートを報知する。

（１－３）データテーブル
第１実施形態のデータテーブルについて説明する。

（１－３－１）空間データテーブル
第１実施形態の空間データテーブルについて説明する。図８は、第１実施形態の空間データテーブルのデータ構造を示す図である。空間データテーブルは、計測装置１０が使用する空間形状情報の一例である。ただし、空間形状情報は計測対象の空間の形状（２次元形状又は３次元形状）を表す情報であればよく、図８に示す例に限定されない。空間形状情報は、例えば事前に対象空間をLiDARスキャナー等のセンサにより測定することで取得され、記憶装置１１に記憶される。ただし空間形状情報の取得方法はこれに限定されず、例えば超音波送信装置２０から送信され壁に反射して超音波受信装置３０に到達する超音波の伝搬時間を用いて空間形状を測定してもよい。

図８の空間データテーブルには、対象空間に関する空間情報が格納される。
空間データテーブルは、「座標」フィールドと、「反射特性」フィールドと、を含む。各フィールドは、互いに関連付けられている。

「座標」フィールドには、対象空間に存在する反射部材の座標（以下「反射部材座標」という）が格納される。反射部材座標は、対象空間の任意の基準点を原点とする座標系（以下「空間座標系」という）で表される。

「反射特性」フィールドには、反射部材の反射特性に関する反射特性情報が格納される。「反射特性」フィールドは、「反射種別」フィールドと、「反射率」フィールドと、「法線角」フィールドと、を含む。

「反射種別」フィールドには、反射種別に関する情報が格納される。反射種別は、以下の何れかである。
・拡散反射
・鏡面反射

「反射率」フィールドには、反射部材の反射率の値が格納される。

「法線角」フィールドには、反射部材の反射面の法線角度の値が格納される。
なお、図８の例では、空間形状情報が各３次元位置の属性情報として反射特性の情報を有するものとしているが、空間形状情報は各３次元位置に関する他の属性情報（例えば窓や換気口の通気特性など）を有していてもよい。

（１－３－２）センサデータテーブル
第１実施形態のセンサデータテーブルについて説明する。図９は、第１実施形態のセンサデータテーブルのデータ構造を示す図である。

図９に示すように、センサデータテーブルには、音波送信装置２０及び音波受信装置３０に関する情報（以下「センサ情報」という）が格納される。
センサデータテーブルは、「センサＩＤ」フィールドと、「座標」フィールドと、「センサタイプ」フィールドと、を含む。
各フィールドは、互いに関連付けられている。

「センサＩＤ」フィールドには、音波送信装置２０又は音波受信装置３０を識別するセンサ識別情報が格納される。

「座標」フィールドには、音波送信装置２０又は音波受信装置３０の位置を示す座標（以下「センサ座標」という）が格納される。センサ座標は、空間座標系で表される。

「センサタイプ」フィールドには、音波送信装置２０であることを示すタグ「送信」、又は、音波受信装置３０であることを示すタグ「受信」が格納される。

（１―３－３）経路データテーブル
第１実施形態の経路データテーブルについて説明する。図１０は、第１実施形態の経路データテーブルのデータ構造を示す図である。

図１０に示すように、経路データテーブルには、経路に関する経路情報が格納される。
経路データテーブルは、「経路ＩＤ」フィールドと、「送信センサ」フィールドと、「受信センサ」フィールドと、を含む。

「経路ＩＤ」フィールドには、メッシュを通る音波の伝搬経路を識別する経路識別情報が格納される。

「送信センサ」フィールドには、経路を構成する音波送信装置２０のセンサ識別情報が格納される。

「受信センサ」フィールドには、経路を構成する音波受信装置３０のセンサ識別情報が格納される。

（１－３－４）メッシュデータテーブル
第１実施形態のメッシュデータテーブルについて説明する。図１１は、第１実施形態のメッシュデータテーブルのデータ構造を示す図である。図１２は、第１実施形態のフィルタの説明図である。

図１１に示すように、メッシュデータテーブルには、メッシュに関するメッシュ情報が格納される。
メッシュデータテーブルは、「メッシュＩＤ」フィールドと、「座標」フィールドと、「経路ＩＤ」フィールドと、「フィルタ」フィールドと、を含む。

「メッシュＩＤ」フィールドには、メッシュを識別するメッシュ識別情報が格納される。

「座標」フィールドには、メッシュの位置を示すメッシュ座標が格納される。メッシュ座標は、空間座標系で表される。
メッシュ座標、または他のパラメータは、BIM（Building Information Modeling）、または他のCAD（Computer-Aided Design）データを参照して定義されてもよい。

「経路ＩＤ」フィールドには、経路の経路識別情報が格納される。

「フィルタ」フィールドには、音波受信装置３０によって受信された受信波形データによって再現される超音波ビームの波形から特定波形を抽出するためのフィルタに関するフィルタ情報が格納される。フィルタ情報は、「経路ＩＤ」フィールドに格納された経路識別情報に関連付けられる。「フィルタ」フィールドは、「時間フィルタ」フィールドと、「振幅フィルタ」フィールドと、を含む。

「時間フィルタ」フィールドには、時間軸に沿って特定波形を抽出するための時間フィルタに関する情報が格納される。時間フィルタは、例えば、以下の少なくとも１つである（図１２）。
・下限時間閾値ＴＨｔｂ
・上限時間閾値ＴＨｔｔ
・下限時間閾値ＴＨｔｂと上限時間閾値ＴＨｔｔとによって規定される時間ウインドウＷｔ

「振幅フィルタ」フィールドには、振幅軸に沿って特定波形を抽出するための振幅フィルタに関する情報が格納される。振幅フィルタは、例えば、以下の少なくとも１つである（図１２）。
・下限振幅閾値ＴＨａｂ
・上限振幅閾値ＴＨａｔ
・下限振幅閾値ＴＨａｂと上限振幅閾値ＴＨａｔとによって規定される振幅ウインドウＷａ

（１－３－５）物品データテーブル
第１実施形態の物品データテーブルについて説明する。図１３は、第１実施形態の物品データテーブルのデータ構造を示す図である。

物品データテーブルには、対象空間ＳＰに収容される物品に関する情報（「物品情報」）が格納される。

図１３に示すように、物品データテーブルは、物品ＩＤフィールドと、名称フィールドと、種別フィールドと、荷主フィールドと、位置フィールドと、基準フィールドと、状態フィールドと含む。

物品ＩＤフィールドには、物品ＩＤが格納される。物品ＩＤは、物品を識別する情報である。

名称フィールドには、名称情報が格納される。名称情報は、物品の名称に関する情報である。

種別フィールドには、種別情報が格納される。種別情報は、物品の種別に関する情報である。

荷主フィールドには、荷主情報が格納される。荷主情報は、物品の荷主に関する情報である。

位置フィールドには、位置情報が格納される。位置情報は、物品の位置に関する情報である。位置情報は、物品が設置されたメッシュ（「保管区画」と称する）を特定可能な情報である。例えば、位置情報は、物品が設置されたメッシュのメッシュＩＤ、または物品が設置された位置の座標情報である。

基準フィールドには、基準情報が格納される。基準情報は、物品の保管環境に課される基準に関する情報である。保管基準は、例えば、物品の保管温度に関する基準（閾値、上限温度、下限温度、平均温度、またはこれらの組み合わせ）を含む。

状態フィールドには、状態情報が格納される。状態情報は、物品の保管状態に関する情報である。一例として、状態情報の値は、物品情報の物品データテーブルへの登録時に「正常」をセットされるが、物品の保管基準を満たさないことが検知された場合に「異常」に書き換えられる。

物品情報は、例えば、対象空間ＳＰに物品を収容する際にユーザ操作に応じて入力され、記憶装置１１に記憶される。なお、物品情報を構成する情報要素の一部（例えば、物品ＩＤ、名称情報、種別情報、荷主情報、位置情報、及び基準情報の少なくとも何れか）は、物品に取り付けられた無線タグ（例えば、ＲＦＩＤ（Radio Frequency Identifier）タグ、又はＢＬＥ（Bluetooth（登録商標）Low Energy）タグ）に格納された情報をリーダーにより読み取ることで取得されてもよい。

（１－４）保管状態管理処理
第１実施形態の保管状態管理処理について説明する。図１４は、第１実施形態の保管状態管理処理のフローチャートである。図１５は、第１実施形態のセンサ配置の一例を示す図である。図１６は、図１４の温度計測処理の詳細なフローチャートである。図１７は、図１６の経路温度の計算の詳細なフローチャートである。図１８は、図１６の処理において表示される画面例を示す図である。

図１４の処理は、開始条件が成立したことに応じて開始する。開始条件は、例えば以下の少なくとも１つを含むことができる。
・計測装置１０に対してユーザにより保管状態管理処理の開始指示が入力された。
・所定の日時が到来した。
・物品が対象空間に収容された。

図１４に示すように、計測装置１０は、温度計測処理（Ｓ１１）を実行する。
具体的には、プロセッサ１２は、対象空間を仮想的に分割した複数の区画（例えば、メッシュ）毎に温度を算出する（つまり、温度分布を算出する）。

図１５に示すように、対象空間ＳＰには、複数の音波送信装置２０ａ～２０ｅと、複数の音波受信装置３０ａ～３０ｅと、が配置される。
複数の音波送信装置２０ａ～２０ｅは、それぞれ、複数の音波受信装置３０ａ～３０ｅに対向している。例えば、音波送信装置２０ａは、音波受信装置３０ａに対向している。これは、音波送信装置２０ａ及び音波受信装置３０ａが、センサペアを形成することを意味している。
図１５の例では、５つのセンサペアが形成される。
計測装置１０は、音波の伝搬距離と伝搬時間とに基づいて、複数の経路を構成するメッシュの温度（以下、「メッシュ温度」と称する）の計測が可能である。
図１５の例では、計測装置１０は、メッシュＭ１～Ｍ４のメッシュ温度の計測が可能である。

プロセッサ１２は、対象空間ＳＰにおいて温度計測の対象となるメッシュ（「計測対象メッシュ」）を、対象空間ＳＰに存在する物品の位置に基づいて決定してもよい。プロセッサ１２は、以下の少なくとも１つの情報を利用して、対象空間ＳＰに存在する物品の位置を特定可能である。
・計測装置１０に接続されたイメージセンサ（不図示）の出力
・物品に取り付けられた無線タグと屋内測位を行う無線通信装置（ロケータ）との間で送受信されるＲＦ信号の到来角または発信角
・物品に取り付けられた無線タグから読み取った位置情報
・物品データテーブル（図１３）に格納された位置情報
ただし、温度計測の対象となる領域の決定方法はこれに限らない。また、プロセッサ１２は、対象空間１２に含まれるすべてのメッシュを温度計測の対象となる領域として決定してもよい。

ステップＳ１１の後、計測装置１０は、保管基準の判定（Ｓ１２）を実行する。
具体的には、プロセッサ１２は、ステップＳ１１において計測したメッシュ温度に基づいて、各物品に関連付けられる基準が満たされているか否かを判定する。一例として、プロセッサ１２は、物品データテーブル（図１３）に格納された基準情報を参照し、物品毎に、物品に関連付けられる基準を特定する。プロセッサ１２は、物品データテーブルに格納された位置情報を参照し、物品毎に、物品の位置に対応するメッシュである保管区画を特定する。プロセッサ１２は、物品毎に、物品の保管区画のメッシュ温度が物品に関連付けられる基準を満たすか否かを判定する。例えば、プロセッサ１２は、物品の保管区画のメッシュ温度が、当該物品に関連付けられる閾値以下であるか否かを判定する。

ステップＳ１２において全ての物品について保管区画のメッシュ温度が基準を満たすと判定した場合に、計測装置１０は、温度計測処理（Ｓ１１）を再実行する。

ステップＳ１２において１つ以上の物品について保管区画のメッシュ温度が基準を満たさないと判定した場合に、計測装置１０は、所定の動作（Ｓ１３）を実行する。
プロセッサ１２は、基準を満たさないと判定した物品毎に動作内容を個別に決定してもよい。また、プロセッサ１２は、基準を満たさないと判定した物品の組み合わせに応じて動作内容を決定してもよい。プロセッサ１２は、例えば以下の少なくとも１つの動作を行う。
・物品の基準に関するアラートの報知
・物品データテーブル（図１３）に格納された状態情報の書き換え
・空調装置４０の制御

所定の動作（Ｓ１３）の第１の例では、プロセッサ１２は、物品の基準に関するアラートを報知する。例えば、プロセッサ１２は、外部装置（例えば、サーバ、または関係者（例えば、物品の管理者、荷主、発荷主、着荷主）の端末）に物品が基準を満たさない環境で保管されている事象を報知するための情報を送信する。一例として、プロセッサ１２は、以下の少なくとも１つを示す情報を送信できる。
・物品情報の少なくとも１つの要素
・基準を逸脱した温度の計測された日時（以下、「計測日時」と称する）
・計測日時における保管区画または他のメッシュの少なくとも１つのメッシュ温度
・保管区画または他のメッシュの少なくとも１つのメッシュ温度の時系列データ（例えば、時系列データをプロットした時間対温度グラフ）
・空調装置４０の動作ログ
・空調装置４０の推奨設定（例えば、どの空調装置４０の設定温度をいくつにすべきか）
・対象空間のマップ画像（例えば、物品の保管区画に対応する位置に当該物品に関連付けられる基準を示す情報が表示される画像）

所定の動作（Ｓ１３）の第２の例では、プロセッサ１２は、計測装置１０に接続されたデバイスを動作させてアラートを行う。一例として、プロセッサ１２は、以下の少なくとも１つを行うことができる。
・ランプを点灯する。
・ディスプレイに所定の画面を表示する。
・スピーカから所定の音声を出力する。

所定の動作（Ｓ１３）の第３の例では、プロセッサ１２は、物品データテーブル（図１３）に格納された状態情報の値を「異常」に書き換える。

所定の動作（Ｓ１３）の第４の例では、プロセッサ１２は、基準を満たさない環境で保管されている物品（以下、「異常物品」と称する）の保管状態が基準を満たすように、空調装置４０に対する制御信号、または制御信号の時系列パターン（つまり、複数時点に亘って空調装置４０に適用される制御信号のセット）を生成し、空調装置４０へ送信する。空調装置４０に対する制御信号は、以下の少なくとも１つを指示する信号であってよい。
・空調装置４０の起動
・空調装置４０の停止
・空調装置４０の運転モードの変更
・空調装置４０の設定温度の変更
・空調装置４０の設定風量の変更
・空調装置４０の設定風向の変更
・空調装置４０の設定湿度の変更
・空調装置４０に内蔵されるモータの回転数の変更
・空調装置４０の加熱冷却を行う交換機の内部温度の変更

一例として、プロセッサ１２は、対象空間ＳＰに収容される物品のうち未だ基準を満たす環境で保管されていると判定されている物品（以下、「正常物品」と称する）の保管区画のメッシュ温度を参照して、空調装置４０に対する制御内容を決定する。プロセッサ１２は、以下の少なくとも１つの動作を実行するように空調動作を制御してもよい。
・正常物品の保管区画の温度が基準を逸脱しないように、且つ異常物品の保管区画の温度が基準を満たすように、空調装置４０を起動または停止する。
・正常物品の保管区画の温度が基準を逸脱しないように、且つ異常物品の保管区画の温度が基準を満たすように、空調装置４０の運転モードを暖房モード、冷房モードまたは他のモードのいずれかに変更する。
・正常物品の保管区画の温度が基準を逸脱しないように、且つ異常物品の保管区画の温度が基準を満たすように、空調装置４０の設定温度を増加または減少させる。
・正常物品の保管区画の温度が基準を逸脱しないように、且つ異常物品の保管区画の温度が基準を満たすように、空調装置４０の設定風量を増加または減少させる。
・正常物品の保管区画の温度が基準を逸脱しないように、且つ異常物品の保管区画の温度が基準を満たすように、空調装置４０の設定風向を空調装置４０から当該保管区画への方向に近づく方向または当該保管区画から遠ざかる方向に変更する。

所定の動作（Ｓ１３）の第５の例は、第１～第４の例のうち２つ以上の組み合わせである。

ステップＳ１３の後、計測装置１０は、温度計測処理（Ｓ１１）を再実行する。

図１４の処理は、終了条件が成立したことに応じて終了する。終了条件は、例えば以下の少なくとも１つを含むことができる。
・計測装置１０に対してユーザにより保管状態管理処理の終了指示が入力された。
・所定の日時が到来した。
・全ての物品が対象空間から搬出された。

以下、温度計測処理（Ｓ１１）の詳細を説明する。
図１６に示すように、計測装置１０は、計測対象メッシュの決定（Ｓ１１０）を実行する。
具体的には、プロセッサ１２は、図１５に示すように、対象空間ＳＰを構成する複数のメッシュの中から計測対象メッシュＭｔ（ｔ＝１～４）のメッシュ識別情報を決定する。一例として、プロセッサ１２は、対象空間ＳＰに収容される物品の保管区画に基づいて計測対象メッシュＭｔを決定する。

ステップＳ１１０の後、計測装置１０は、所定の経路温度計算モデルに従って、経路温度の計算（Ｓ１１１）を実行する。
図１７を参照して、ステップＳ１１１の詳細を説明する。

計測装置１０は、対象経路の決定（Ｓ１１１０）を実行する。
具体的には、プロセッサ１２は、メッシュデータテーブル（図１１）を参照して、ステップＳ１１０で決定したメッシュ識別情報に関連付けられた「経路ＩＤ」フィールドの情報（つまり、計測対象メッシュＭｔを通る経路（以下「対象経路」という）Ｐｉ（ｉは、経路の引数）の経路識別情報）を特定する。

ステップＳ１１１０の後、計測装置１０は、超音波ビームの出力（Ｓ１１１１）を実行する。
具体的には、プロセッサ１２は、経路データテーブル（図１０）を参照して、ステップＳ１１１０で特定した経路識別情報に関連付けられた「送信センサ」フィールドの情報（つまり、制御対象となる音波送信装置（以下、「対象音波送信装置」という）２０）と、「受信センサ」フィールドの情報（つまり、制御対象となる音波受信装置（以下「対象音波受信装置」という）３０）と、を特定する。
プロセッサ１２は、対象音波送信装置２０に超音波制御信号を送信する。

対象音波送信装置２０は、計測装置１０から送信された超音波制御信号に応じて超音波ビームを送信する。
具体的には、複数の超音波振動子２１は、超音波制御信号に応じて同時に振動する。
これにより、対象音波送信装置２０から対象音波受信装置３０に向かって、送信方向（Ｚ軸方向）に進行する超音波ビームが送信される。

ステップＳ１１１１の後、計測装置１０は、受信波形データの取得（Ｓ１１１２）を実行する。
具体的には、対象音波受信装置３０の超音波振動子３１は、ステップＳ１１１１で対象音波送信装置２０から送信された超音波ビームを受信することにより振動する。
制御回路３２は、超音波振動子３１の振動に応じた受信波形データ（図１２）を生成する。
制御回路３２は、生成した受信波形データを計測装置１０に送信する。

計測装置１０のプロセッサ１２は、音波受信装置３０から送信された受信波形データを取得する。プロセッサ１２は、取得した受信波形データに対して、増幅・帯域制限処理などの信号処理をしてもよい

ステップＳ１１１２の後、計測装置１０は、フィルタリング（Ｓ１１１３）を実行する。
具体的には、プロセッサ１２は、メッシュデータテーブル（図１１）を参照して、ステップＳ１１０で決定した対象経路Ｐｉの経路識別情報に関連付けられた「フィルタ」フィールドを特定する。
例えば、計測対象メッシュのメッシュ識別情報が「Ｍ００１」である場合、以下のフィルタ情報が特定される。
・経路識別情報「Ｐ００１」：時間閾値ＴＨｔｂ１以上時間閾値ＴＨｔｔ１未満、及び、振幅閾値ＴＨａｂ１以上振幅閾値ＴＨａｔ１未満
・経路識別情報「Ｐ００２」：時間ウインドウＷｔ２内及び振幅ウインドウＷａ２内

プロセッサ１２は、特定したフィルタ情報に基づいて、受信波形データに含まれる成分のうち、対象経路Ｐｉに沿って進行した超音波ビームの成分を抽出する。

ステップＳ１１１３の後、計測装置１０は、経路温度の計算（Ｓ１１１４）を実行する。
具体的には、プロセッサ１２は、センサデータテーブル（図９）の「座標」フィールドを参照して、センサペア毎に、センサペアを構成する音波送信装置２０の座標及び音波受信装置３０の座標を特定する。
プロセッサ１２は、特定した音波送信装置２０の座標及び音波受信装置３０の座標の組合せに基づいて、当該音波送信装置２０と当該音波受信装置３０との間の距離（以下「センサ間距離」という）Ｄｓを計算する。
プロセッサ１２は、ステップＳ１１１３で抽出された成分のピーク値に対応する時間（以下「伝搬時間」という）ｔを特定する。伝搬時間ｔは、音波送信装置２０が超音波ビームを送信してから、対象経路Ｐｉに沿って進行した超音波ビームが音波受信装置３０に到達するまでの所要時間（つまり、対象経路の始点から終点までを超音波ビームが伝搬する時間）を意味する。
プロセッサ１２は、超音波の音速Ｃ、センサ間距離Ｄｓと、伝搬時間ｔと、基準温度Ｔ０と、を用いて、対象経路Ｐｉの経路温度ＴＥＭＰｐａｔｈｉを計算する。

全ての対象経路ＰｉについてステップＳ１１１４が終了していない場合（Ｓ１１１５－ＮＯ）、計測装置１０は、ステップＳ１１１０を実行する。

全ての対象経路ＰｉについてステップＳ１１１４が終了すると（Ｓ１１１５－ＹＥＳ）、計測装置１０は、図１６のメッシュ温度の計算（Ｓ１１２）を実行する。
具体的には、プロセッサ１２は、ステップＳ１１１４（図１７）において計算された複数の対象経路Ｐｉの経路温度ＴＥＭＰｐａｔｈｉを用いて、計測対象メッシュＭｔのメッシュ温度ＴＥＭＰｍｅｓｈｔを計算する（式１）。
ＴＥＭＰｍｅｓｈｔ ＝ ＡＶＥ（ＴＥＭＰｐａｔｈｉ）…（式１）
・ＡＶＥ（ｘ）：ｘの平均値を求める関数

全ての計測対象メッシュＭｔについてステップＳ１１２が終了していない場合（Ｓ１１３－ＮＯ）、計測装置１０は、ステップＳ１１０を実行する。

全ての計測対象メッシュＭｔについてステップＳ１１２が終了すると（Ｓ１１３－ＹＥＳ）、計測装置１０は、計測結果の提示（Ｓ１１４）を実行する。
具体的には、プロセッサ１２は、画面Ｐ１０（図１８）をディスプレイに表示する。

画面Ｐ１０は、表示オブジェクトＡ１０を含む。
表示オブジェクトＡ１０には、画像ＩＭＧ１０が表示される。
画像ＩＭＧ１０は、対象空間ＳＰを構成する複数のメッシュのそれぞれについて、ステップＳ１１２で計算されたメッシュ温度ＴＥＭＰｍｅｓｈｔを示している。

（１－５）小括
第１実施形態の計測装置１０は、物品を収容する空間を仮想的に分割した複数の区画に亘る温度の分布を音波の伝搬特性に基づいて算出し、物品の保管区画の温度が当該物品に関連付けられる基準を満たすか否かを判定する。ここで、各区画の温度を算出するために必要とされる超音波送信装置および超音波受信装置の数は、空間に収容される物品の数に殆ど依存しない。また、超音波送信装置および超音波受信装置は、物品に取り付けられているのでなく対象空間に設置されるので、対象空間に収容される物品が入れ替わっても繰り返し使用することができる。故に、計測装置１０によれば、物品が所定の基準で保管されているか否かを効率的（高コストパフォーマンス）に監視することができる。

計測装置は、保管区画の温度が閾値を超える場合に、保管区画の温度が基準を満たさないと判定してもよい。これにより、物品が閾値以下の温度で保管されているか否かを効率的（高コストパフォーマンス）に監視することができる。

計測装置１０は、物品の保管区画の温度が基準を満たさないと判定した場合に、所定の動作を行う。所定の動作は、物品の基準に関するアラート（一例として、推奨する空調設定）の報知であってもよい。これにより、関係者（例えば、物品の管理者、荷主、発荷主、着荷主）に、物品を欠陥品として間引くように促したり、空調装置４０の設定を変更するように促したりすることができる。或いは、所定の動作は、空調装置４０の制御であってもよい。これにより、正常物品の保管区画の温度が基準を逸脱しないように対象空間内の温度分布を管理することができる。

（２）第２実施形態
第２実施形態について説明する。但し、上述の実施形態と同様の説明は省略する。

（２－１）実施形態の概要
第２実施形態の概要について説明する。図１９は、第２実施形態の概要の説明図である。

第２実施形態において、計測装置１０は第１実施形態の場合と同様に、対象空間ＳＰのメッシュ（区画）単位で温度を計測する。第２実施形態の計測装置１０は、計測したメッシュ温度を参照して将来に亘る物品の保管区画の予想温度を算出する。一例として、計測装置１０は、時刻ｔまでに計測したメッシュＴＭ１の温度を参照して補間処理を行うことで、時刻ｔ以降の予想温度を算出する。

計測装置１０は、メッシュＴＭ１の現在温度のみならず将来の予想温度を参照し、時刻ｔ＋δにおける予想温度が、物品ＴＩ１に関連付けられる基準（例えば、保管温度が閾値Ｔｈ１を超過しない）を満たすか否かを判定することができる。

図１９の例では、時刻ｔ＋δにおけるメッシュＴＭ１の予想温度が閾値Ｔｈ１を超過する。故に、計測装置１０は、メッシュＴＭ１の予想温度が物品ＴＩ１に関連付けられる基準を満たさないと判定する。

計測装置１０は、メッシュＴＭ１の予想温度が物品ＴＩ１に関連付けられる基準を満たさないと判定した場合に、所定の動作を行う。一例として、計測装置１０は、物品ＴＩ１の基準に関するアラートを報知する。

（２－２）保管状態管理処理
第２実施形態の保管状態管理処理について説明する。図２０は、第２実施形態の保管状態管理処理のフローチャートである。

図２０の処理は、開始条件が成立したことに応じて開始する。開始条件は、例えば以下の少なくとも１つを含むことができる。
・計測装置１０に対してユーザにより保管状態管理処理の開始指示が入力された。
・所定の日時が到来した。
・物品が対象空間に収容された。

図２０に示すように、計測装置１０は図１４と同様に、温度計測処理（Ｓ１１）を実行する。

ステップＳ１１の後、計測装置１０は、温度の予測（Ｓ２２）を実行する。
具体的には、プロセッサ１２は、ステップＳ１１において計測したメッシュ温度に基づいて、各物品の保管区画の将来に亘る予想温度を算出する。一例として、プロセッサ１２は、物品データテーブル（図１３）に格納された位置情報を参照し、物品毎に、保管区画を特定する。プロセッサ１２は、保管区画毎に、現在および過去の複数時点に亘るメッシュ温度を取得する。プロセッサ１２は、複数のメッシュ温度を参照して補間処理を行うことで、保管区画の将来に亘る予想温度を算出する。補間は、線形補間、３次スプライン補間、多項式補間、または他の任意の種類の補間であってもよい。

ステップＳ２２の後、計測装置１０は、保管基準の判定（Ｓ２３）を実行する。
具体的には、プロセッサ１２は、ステップＳ１１において計測したメッシュ温度、およびＳ２２において算出した予想温度に基づいて、各物品に関連付けられる基準が満たされているか否かを判定する。一例として、プロセッサ１２は、物品データテーブル（図１３）に格納された基準情報を参照し、物品毎に、物品に関連付けられる基準を特定する。プロセッサ１２は、物品データテーブルに格納された位置情報を参照し、物品毎に、保管区画を特定する。プロセッサ１２は、物品毎に、物品の保管区画の現在温度（ステップＳ１１において計測したメッシュ温度）および予想温度が物品に関連付けられる基準を満たすか否かを判定する。例えば、プロセッサ１２は、物品の保管区画の現在温度および予想温度が、当該物品に関連付けられる閾値以下であるか否かを判定する。

ステップＳ２３において全ての物品について保管区画の現在温度および予想温度が基準を満たすと判定した場合に、計測装置１０は、温度計測処理（Ｓ１１）を再実行する。

ステップＳ２３において１つ以上の物品について保管区画の現在温度または予想温度が基準を満たさないと判定した場合に、計測装置１０は、所定の動作（Ｓ２４）を実行する。
プロセッサ１２は、基準を満たさないと判定した物品毎に動作内容を個別に決定してもよい。また、プロセッサ１２は、基準を満たさないと判定した物品の組み合わせに応じて動作内容を決定してもよい。プロセッサ１２は、例えば以下の少なくとも１つの動作を行う。
・物品の基準に関するアラートの報知
・物品データテーブル（図１３）に格納された状態情報の書き換え
・空調装置４０の制御

プロセッサ１２は、保管区画の現在温度が基準を満たさないと判定した場合に、第１～第５の例に従う動作を行ってもよい。所定の動作（Ｓ２４）の第１～第５の例は、所定の動作（Ｓ１３）の第１～第５の例と同一である。

プロセッサ１２は、保管区画の現在温度が基準を満たさないと判定した場合に、以下に説明する第６～第９の例に従う動作を行ってもよい。

所定の動作（Ｓ２４）の第６の例では、プロセッサ１２は、物品の基準に関するアラートを報知する。例えば、プロセッサ１２は、外部装置（例えば、サーバ、または関係者（例えば、物品の管理者、荷主、発荷主、着荷主）の端末）に物品が将来的に基準を満たさない可能性がある環境で保管されている事象を報知するための情報を送信する。一例として、プロセッサ１２は、以下の少なくとも１つを示す情報を送信できる。
・物品情報の少なくとも１つの要素
・基準を逸脱すると判定された予想温度に対応する日時（以下、「予想日時」と称する）
・予想日時の到来までの猶予時間
・予想日時における保管区画または他のメッシュの少なくとも１つの予想温度
・保管区画または他のメッシュの少なくとも１つの計測温度または予測温度の少なくとも１つの時系列データ（例えば、時系列データをプロットした時間対温度グラフ）
・空調装置４０の動作ログ
・空調装置４０の推奨設定（例えば、どの空調装置４０の設定温度をいくつにすべきか）
・対象空間のマップ画像（例えば、物品の保管区画に対応する位置に当該物品に関連付けられる基準を示す情報が表示される画像）

所定の動作（Ｓ２４）の第７の例では、プロセッサ１２は、計測装置１０に接続されたデバイスを動作させてアラートを行う。一例として、プロセッサ１２は、以下の少なくとも１つを行うことができる。
・ランプを点灯する。
・ディスプレイに所定の画面を表示する。
・スピーカから所定の音声を出力する。
所定の動作（Ｓ２４）の第７の例において、ランプの点灯パターン、ディスプレイの表示画面、またはスピーカの出力音声の少なくとも１つが、上記第２の例と異なっていてもよい。

所定の動作（Ｓ２４）の第８の例では、プロセッサ１２は、空調装置４０に対する制御信号、または制御信号の時系列パターンを生成し、空調装置４０へ送信する。空調装置４０に対する制御信号は、以下の少なくとも１つを指示する信号であってよい。
・空調装置４０の起動
・空調装置４０の停止
・空調装置４０の運転モードの変更
・空調装置４０の設定温度の変更
・空調装置４０の設定風量の変更
・空調装置４０の設定風向の変更
・空調装置４０の設定湿度の変更
・空調装置４０に内蔵されるモータの回転数の変更
・空調装置４０の加熱冷却を行う交換機の内部温度の変更

一例として、プロセッサ１２は、基準を逸脱する予想温度が算出された物品（以下、「要注意物品」と称する）の保管区画の現在温度または予想温度を参照して、空調装置４０に対する制御内容を決定する。プロセッサ１２は、以下の少なくとも１つの動作を実行するように空調動作を制御してもよい。
・要注意物品の保管区画の予想温度が基準を逸脱しないように、空調装置４０を起動または停止する。
・要注意物品の保管区画の予想温度が基準を逸脱しないように、空調装置４０の運転モードを暖房モード、冷房モードまたは他のモードのいずれかに変更する。
・要注意物品の保管区画の予想温度が基準を逸脱しないように、空調装置４０の設定温度を増加または減少させる。
・要注意物品の保管区画の予想温度が基準を逸脱しないように、空調装置４０の設定風量を増加または減少させる。
・要注意物品の保管区画の予想温度が基準を逸脱しないように、空調装置４０の設定風向を空調装置４０から当該保管区画への方向に近づく方向または当該保管区画から遠ざかる方向に変更する。

所定の動作（Ｓ２４）の第９の例は、第６～第８の例のうち２つ以上の組み合わせである。

ステップＳ２４の後、計測装置１０は、温度計測処理（Ｓ１１）を再実行する。

図２０の処理は、終了条件が成立したことに応じて終了する。終了条件は、例えば以下の少なくとも１つを含むことができる。
・計測装置１０に対してユーザにより保管状態管理処理の終了指示が入力された。
・所定の日時が到来した。
・全ての物品が対象空間から搬出された。
（２－３）小括
第２実施形態の計測装置１０は、複数時点に亘る保管区画の温度を参照して将来に亘る当該保管区画の予想温度を推定する。計測装置１０は、予想温度が閾値を超える場合に、保管区画の温度が基準を満たさないと判定する。故に、計測装置１０によれば、物品が将来的に基準を満たさない可能性がある環境で保管されているか否かを効率的（高コストパフォーマンス）に監視することができる。

計測装置１０は、複数時点に亘る物品の保管区画の温度に対して補間処理を行うことにより保管区画の予測温度を算出してもよい。これにより、保管区画の予測温度を数値計算により求めることができる。

計測装置１０は、物品の保管区画の予測温度の時間変化に関する情報（例えば、予測温度の時系列データ）を出力してもよい。これにより、関係者（例えば、物品の管理者、荷主、発荷主、着荷主）に、物品の保管環境が現実に基準を逸脱する前に、予防手段を講じるよう支援することが可能となる。

（３）変形例
本実施形態の変形例について説明する。

（３－１）変形例１
変形例１について説明する。変形例１は、時系列フィルタを用いた温度計測アルゴリズムの例である。

（３－１－１）変形例１の概要
変形例１の概要について説明する。図２１は、変形例１の概要の説明図である。
図２１に示すように、変形例１のプロセッサ１２は、経路温度計算モデルＭｐｔ（ｔ）と、時系列フィルタＦＩＬと、を実行するように構成される。

経路温度計算モデルＭｐｔ（ｔ）は、時刻ｔの受信波形データＲＷ（ｔ｜ｘ，ｙ，ｚ）に応じて、時刻ｔの経路温度ＰＤ（ｔ｜ｘ，ｙ，ｚ）を出力するように構成される。

時系列フィルタＦＩＬは、経路温度計算モデルＭｐｔ（ｔ）の出力（経路温度ＰＤ（ｔ｜ｘ，ｙ，ｚ））、温度計５０によって計測された時刻ｔの基準温度Ｔｒｅｆ（ｔ）、及び、時刻ｔ－１の温度分布Ｄ（ｔ－１）の組合せに応じて、時刻ｔの温度分布Ｄｔ（ｔ）を出力するように構成される。
時系列フィルタＦＩＬは、例えば、以下の少なくとも１つを含む。
・カルマンフィルタ
・拡張カルマンフィルタ
・無香カルマンフィルタ
・パーティクルフィルタ

（３－１－２）温度計測処理
変形例１の温度計測処理について説明する。図２２は、変形例１の温度計測処理のフローチャートである。

図２２に示すように、変形例１の計測装置１０は、図１６と同様に、ステップＳ１１０～Ｓ１１３を実行する。

全ての計測対象メッシュＭｔについてステップＳ１１２が終了していない場合（Ｓ１１３－ＮＯ）、計測装置１０は、ステップＳ１１０を実行する。

全ての計測対象メッシュＭｔについてステップＳ１１２が終了すると（Ｓ１１３－ＹＥＳ）、計測装置１０は、時系列フィルタリング（Ｓ３１０）を実行する。
具体的には、プロセッサ１２は、温度計５０から時刻ｔの基準温度Ｔｒｅｆ（ｔ）を取得する。
プロセッサ１２は、ステップＳ１１１で得られた時刻ｔの経路温度Ｔｐ（ｔ｜ｘ，ｙ，ｚ）と、基準温度Ｔｒｅｆ（ｔ）と、を時刻ｔ－１の温度分布（ｔ－１）と、を時系列フィルタＦＩＬに入力することにより、時刻ｔの温度分布Ｄ（ｔ）を計算する。
温度分布Ｄ（ｔ）は、時刻ｔ＋１の温度分布Ｄ（ｔ＋１）の計算において参照される。

変形例１によれば、時系列フィルタリングを実行することにより、空間の温度の計測結果のＳ／Ｎ比を更に向上させることができる。

なお、変形例１の時系列フィルタＦＩＬは、更に、時刻ｔ－１の外部環境情報を参照して、時刻ｔの温度分布Ｄ（ｔ）を計算しても良い。時刻ｔ－１の外部環境情報は、例えば、以下の情報を含む。
・空調装置４０の熱量に関する情報
・対象空間ＳＰの周辺の外気温に関する情報
・対象空間ＳＰの周辺の気温に関する情報
・対象空間ＳＰの３次元形状に関する情報
・対象空間ＳＰの断熱性能に関する情報
・対象空間ＳＰ内に存在する人の数に関する情報
・対象空間ＳＰ内に存在する人の動きに関する情報
・空調装置４０の風に関する情報
・対象空間ＳＰ内の風に関する情報

（３－２）変形例２
変形例２について説明する。変形例２は、温度に加えて、風ベクトルの分布を計測する例である。

変形例２の対象空間ＳＰには、例えば、図５の計測装置１０（不図示）と、少なくとも２つの音波送信装置２０ａ～２０ｂと、少なくとも２つの音波受信装置３０ａ～３０ｂと、が配置されている。

計測装置１０は、所定の送信周波数Ｆｓを有する超音波ビームを送信するように、音波送信装置２０ａを制御する。
音波送信装置２０ａは、計測装置１０の制御に従い、送信周波数Ｆｓを有する超音波ビームを送信する。
音波受信装置３０ａは、超音波ビームを受信すると、受信波形データを生成する。音波受信装置３０ａが受信する超音波ビームには、音波送信装置２０ａと音波受信装置３０ａとの間の風に起因するドップラー効果が発生する。したがって、音波受信装置３０ａが受信する超音波ビームの受信周波数Ｆｒａは、送信周波数Ｆｓとは異なる。
計測装置１０は、音波受信装置３０ａから受信波形データを取得し、且つ、受信波形データを参照して受信周波数Ｆｒａを特定する。
ドップラー効果を考慮すると、送信周波数Ｆｓと受信周波数Ｆｒａとの間には、式２．１の関係が成立する。
Ｆｓ ＝ （Ｃ＋Ｖｗａ）／Ｃ×Ｆｒａ …（式２．１）
・Ｃ：超音波の音速
・Ｖｗａ：音波送信装置２０ａと音波受信装置３０ａとの間の経路上の理論上の風速

式２．１を展開すると、理論上の風速Ｖｗａは、式２．２のように表すことができる。プロセッサ１２は、式２．１を用いて、理論上の風速Ｖｗａを計算する。
Ｖｗａ ＝ Ｃ×Ｆｓ／Ｆｒａ－Ｃ …（式２．２）

プロセッサ１２は、式２．３を用いて、音波送信装置２０ｂと音波受信装置３０ｂとの間の経路上の理論上の風速Ｖｗｂを計算する。
Ｖｗｂ ＝ Ｃ×Ｆｓ／Ｆｒｂ－Ｃ …（式２．３）
・Ｆｒｂ：音波受信装置３０ｂが受信した超音波ビームの受信周波数

記憶装置１１には、温度と風速の相関関数Ｇ（ｘ）が予め格納されている。
プロセッサ１２は、ステップＳ１１２で得られたメッシュ温度ＴＥＭＰｍｅｓｈ及び相関関数を用いて、式２．４～式２．５のように、補正風速Ｖｒｗａ及びＶｒｗｂを計算する。
Ｖｒｗａ ＝ Ｇ（Ｔｍｅｓｈ）×Ｖｗａ …（式２．４）
Ｖｒｗｂ ＝ Ｇ（Ｔｍｅｓｈ）×Ｖｗｂ …（式２．５）

計測装置１０は、音波送信装置２０ａ～音波受信装置３０ａの進行経路及び音波送信装置２０ｂ～音波受信装置３０ｂの進行経路の成す角度を参照して、式２．４～式２．５から得られた補正風速Ｖｒｗａ及びＶｒｗｂを合成することにより、風ベクトルを計算する。

変形例２によれば、超音波の進行経路が交差する交差点を含むメッシュの風ベクトルが得られる。したがって、複数の交差点が形成される場合、対象空間ＳＰの風ベクトルの分布（つまり、風速及び風向きの分布）が得られる。

変形例２では、２本の進行経路が交差点を形成する例を示した。この場合、得られる風ベクトルの次元は２次元である。
なお、変形例２は、３本の進行経路が交差点を形成する場合にも適用可能である。この場合、得られる風ベクトルの次元は３次元である。

（３－３）変形例３
変形例３について説明する。変形例３は、超音波を用いて風ベクトルの分布を計測する例である。図２３は、変形例３のセンサ配置の一例を示す図である。図２４は、変形例３の音波送信装置及び音波受信装置の構成を示す概略図である。

図２３に示すように、変形例３の対象空間ＳＰには、少なくとも４つのセンサユニットＳＵａ～ＳＵｄが、センサユニットＳＵａ～ＳＵｂの超音波の進行経路と、センサユニットＳＵｃ～ＳＵｄの超音波の進行経路とが互いに交わる（好ましくは、直交する）ように配置される。
図２４に示すように、各センサユニットＳＵａ～ＳＵｄは、それぞれ、音波送信装置２０と、音波受信装置３０と、を備える（つまり、送受兼用である）。

図２３に示すように、音波送信装置２０ａから音波受信装置３０ｂに向かう経路（以下「往路」という）を進行する超音波ビームの実速度の絶対値｜Ｖａｂ｜については、式３ａの関係が成立する。
｜Ｖａｂ｜ ＝ Ｄａｂ／ｔａｂ ＝ Ｃ＋Ｖｗａｂ …（式３ａ）
・Ｄａｂ：センサユニットＳＵａ及びＳＵｂの間のセンサ間距離
・ｔａｂ：往路の伝搬時間
・Ｖｗａｂ：センサユニットＳＵａ及びＳＵｂの間の風速成分

音波送信装置２０ｂから音波受信装置３０ａに向かう経路（以下「復路」という）は、往路の逆方向である。したがって、復路を進行する超音波ビームの実速度の絶対値｜Ｖｂａ｜については、式３ｂの関係が成立する。
｜Ｖｂａ｜ ＝ Ｄａｂ／ｔｂａ ＝ Ｃ－Ｖｗａｂ …（式３ｂ）
・ｔｂａ：復路の伝搬時間

計測装置１０は、式３ｃを用いて、センサユニットＳＵａ及びＳＵｂの平均速度｜Ｖａ：ｂ｜を計算する。
｜Ｖａ：ｂ｜ ＝ （｜Ｖａｂ｜＋｜Ｖｂａ｜）／２ …（式３ｃ）

式３ｃにより、往路の風速成分Ｖａｂと復路の風速成分Ｖｂａが互いに相殺される。したがって、平均速度｜Ｖａ：ｂ｜は、風速成分Ｖｗａｂを含まない。

変形例３の計測装置１０は、式４．１を用いて、センサユニットＳＵａ～ＳＵｂ間の往路を進行する超音波ビームの速度の絶対値｜Ｖａｂ｜から温度因子を除去することにより、センサユニットＳＵａ及びＳＵｂの間の風速成分Ｖｗａｂを計算する。
Ｖｗａｂ ＝ ｜Ｖａｂ｜－｜Ｖａ：ｂ｜ …（式４．１）

計測装置１０は、式４．２を用いて、センサユニットＳＵｃ～ＳＵｄ間の往路を進行する超音波ビームの速度の絶対値｜Ｖｃｄ｜から温度因子を除去することにより、センサユニットＳＵｃ及びＳＵｄの間の風速成分Ｖｗｃｄを計算する。
Ｖｗｃｄ ＝ ｜Ｖｃｄ｜－｜Ｖｃ：ｄ｜ …（式４．２）
・｜Ｖｃ：ｄ｜：センサユニットＳＵｃ及びＳＵｄの平均速度

計測装置１０は、センサユニットＳＵａ～ＳＵｂの進行経路及びセンサユニットＳＵｃ～ＳＵｄの進行経路の成す角度を参照して、式４．１から得られた風速成分Ｖｗａｂと、式４．２から得られた風速成分Ｖｗｃｄと、を合成することにより、風ベクトルを計算する。

変形例３によれば、対象空間ＳＰに配置されたセンサユニットの超音波の進行経路が交差点を含むメッシュの風ベクトルが得られる。したがって、複数の交差点が形成される場合、対象空間ＳＰの温度分布に代えて、風ベクトル（つまり、風速及び風向き）の分布が得られる。

変形例３では、往路を進行する超音波ビームの絶対値｜Ｖａｂ｜及び｜Ｖｃｄ｜に代えて、復路を進行する超音波ビームの絶対値｜Ｖｂａ｜及び｜Ｖｄｃ｜を用いて、風ベクトルを推定しても良い。

変形例３では、一対のセンサユニットが交差点を形成する例を示した。この場合、得られる風ベクトルの次元は２次元である。
なお、変形例３は、３つのセンサユニットが交差点を形成する場合にも適用可能である。この場合、得られる風ベクトルの次元は３次元である。

（３－４）変形例４
変形例４について説明する。変形例４は、温度を含む複数種の物理量の計測値を参照して将来の温度を予測する例である。

変形例４の保管状態管理処理について説明する。図２５は、変形例４の保管状態管理処理のフローチャートである。

図２５の処理は、開始条件が成立したことに応じて開始する。開始条件は、例えば以下の少なくとも１つを含むことができる。
・計測装置１０に対してユーザにより保管状態管理処理の開始指示が入力された。
・所定の日時が到来した。
・物品が対象空間に収容された。

図２５に示すように、計測装置１０は、物理量計測処理（Ｓ３１）を実行する。
具体的には、プロセッサ１２は、プロセッサ１２は、対象空間を仮想的に分割した区画（例えば、メッシュ）毎に温度を含む複数種の物理量を計測する（つまり、複数種の物理量の分布を算出する）。一例として、プロセッサ１２は、温度、風速、および風向を区画毎に計測する。

プロセッサ１２は、図１６または図２２に示した温度計測処理、または他の技法により、各区画の温度を計測してもよい。プロセッサ１２は、変形例２または変形例３において説明した技法、または他の技法により、各区画の風速、および風向を計測してもよい。

ステップＳ３１の後、計測装置１０は、温度の予測（Ｓ３２）を実行する。
具体的には、プロセッサ１２は、ステップＳ３１において計測した各区画の温度を含む複数種の物理量に基づいて、各物品の保管区画の将来に亘る予想温度を算出する。一例として、プロセッサ１２は、空間データテーブル（図８）およびメッシュデータテーブル（図１１）を参照し、対象空間の形状情報を取得する。プロセッサ１２は、物品データテーブル（図１３）に格納された位置情報を参照し、物品毎に、保管区画を特定する。プロセッサ１２は、計測日時または位置の少なくとも１つが異なる複数の区画の複数種の物理量を取得する。プロセッサ１２は、取得した対象空間の形状情報と、複数種の物理量とを参照して保管区画の温度をシミュレートすることで、保管区画の将来に亘る予想温度を算出する。

一例として、プロセッサ１２は、ＣＦＤ（Computational Fluid Dynamics）シミュレーションを用いて予想を行うことができる。ＣＦＤシミュレーションは、温度、風速、境界条件が与えられている場合に、連続の式、ナビエ・ストークス方程式、またはエネルギー方程式を用いて、空間内の温度／風速の分布を求める手法である。ただし、プロセッサ１２は、ＣＦＤシミュレーションとは異なるシミュレーションを行ってもよい。プロセッサは、ある位置または時間の物理量に基づいて他の位置または時間の物理量を算出可能な任意のシミュレーションを採用可能である。シミュレーションでは、以下のパラメータの少なくとも１つを参照可能である。
・対象空間ＳＰの少なくとも１つの区画の物理量
・対象空間ＳＰの内部環境に関する動的なパラメータ（例えば、温度計５０によって測定された基準温度、空調装置４０の熱量、対象空間ＳＰ内に存在する人の数、位置または動きに関する情報）
・対象空間ＳＰの内部環境に関する静的なパラメータ（例えば、対象空間ＳＰの３次元形状、対象空間ＳＰと外界とを隔てる構造体（例えば、壁、床、および天井）の断熱性能（例えば、素材および厚み））
・対象空間ＳＰの外部環境に関する動的または静的なパラメータ（例えば、日照条件、外気温、気圧）

ステップＳ３２の後、計測装置１０は図２０と同様に、保管基準の判定（Ｓ２３）を実行する。

ステップＳ２３において全ての物品について保管区画の現在温度および予想温度が基準を満たすと判定した場合に、計測装置１０は、物理量計測処理（Ｓ３１）を再実行する。

ステップＳ２３において１つ以上の物品について保管区画の現在温度または予想温度が基準を満たさないと判定した場合に、計測装置１０は図２０と同様に、所定の動作（Ｓ２４）を実行する。

ステップＳ２４の後、計測装置１０は、物理量計測処理（Ｓ３１）を再実行する。

図２５の処理は、終了条件が成立したことに応じて終了する。終了条件は、例えば以下の少なくとも１つを含むことができる。
・計測装置１０に対してユーザにより保管状態管理処理の終了指示が入力された。
・所定の日時が到来した。
・全ての物品が対象空間から搬出された。

変形例４の計測装置１０は、複数の区画に亘る複数種の物理量（温度を含む）の分布を算出し、計測日時または位置の異なる複数の区画の複数種の物理量を参照して将来に亘る保管区画の予想温度を推定する。計測装置１０は、予想温度が閾値を超える場合に、保管区画の温度が基準を満たさないと判定する。これにより、保管区画の予測温度を高精度に推定することができる。

上記複数種の物理量は、風速、または風向の少なくとも１つを含んでもよい。これにより、流体による熱の移動を考慮して、保管区画の予測温度を高精度に推定することができる。

（４）その他の変形例
その他の変形例について説明する。

記憶装置１１は、ネットワークＮＷを介して、計測装置１０と接続されてもよい。

上記説明では、各種の入力デバイスまたは出力デバイスが入出力インタフェースを介してプロセッサに接続される例を示した。しかしながら、各種の入力デバイスまたは出力デバイスは通信インタフェースを介してプロセッサに接続されてもよい。

図４の例では、超音波振動子３１を備える音波受信装置３０の例を示した。しかし、音波受信装置３０は、音波送信装置２０と同様に、複数の超音波振動子３１を備えても良い。
音波送信装置２０は、アレイ状に配列された複数の超音波振動子２１を備えていてもよい。複数の超音波振動子２１は、それぞれ個別の駆動制御信号によって制御されてもよいし、グループ単位で同一の駆動制御信号によって制御されてもよいし、全体で同一の駆動制御信号によって制御されてもよい。
音波受信装置３０は、アレイ状に配列された複数の超音波振動子３１を備えていてもよい。複数の超音波振動子３１は、それぞれ個別の駆動制御信号によって制御されてもよいし、グループ単位で同一の駆動制御信号によって制御されてもよいし、全体で同一の駆動制御信号によって制御されてもよい。

図５の例では、１個の音波送信装置２０が複数の経路に沿った超音波ビームを送信し、且つ、１個の音波受信装置３０が複数の経路に沿った超音波ビームを受信する例を示した。しかし、本実施形態はこれに限られない。ｎ（ｎは２以上の整数）個の音波送信装置２０のそれぞれが１本の経路に沿った超音波ビーム（つまり、ｎ個の音波送信装置２０がｎ本の経路に沿った超音波ビーム）を送信し、且つ、ｎ個の音波受信装置３０のそれぞれが各経路に沿った超音波ビームを受信しても良い（つまり、ｎ個の音波受信装置３０がｎ本の経路に沿った超音波ビームを受信しても良い）。

上記の実施形態では、メッシュ温度ＴＥＭＰｍｅｓｈｔの計算に平均値を求める関数を用いる例を示したが、本実施形態のメッシュ温度ＴＥＭＰｍｅｓｈｔの計算方法はこれに限られるものではない。

上記の実施形態では、計測装置１０が、音波の伝搬特性に基づいて、当該音波の伝搬路上にある区画の温度を計測する例を示した。しかしながら、計測装置１０は、計測した区画の温度に基づいてシミュレートすることにより、音波の伝搬路上にない区画の温度を推定してもよい。計測装置１０は、シミュレーションにおいて、温度とともに計測した風速、風向、もしくは他の物理量、または過去に計測または推定した物理量をさらに参照してもよい。

上記の実施形態では、計測装置１０が、音波の伝搬特性（例えば、伝搬時間、振幅の変化、位相の変化、および周波数の変化、など）に基づいて、複数の区画に亘る温度または風ベクトル（つまり、風速および風量）の分布を算出する例を示した。しかしながら、計測装置１０は、対象空間ＳＰの環境に関する任意の物理量の分布を算出してもよい。一例として、計測装置１０は、以下の物理量の分布を計測することも可能である。
・空気中の化学物質（例えば、ＣＯ２）の濃度の分布
・湿度の分布
・臭気の分布
・有毒ガスの分布

上記の実施形態では、計測装置１０が、物品の保管区画の温度が当該物品に関連付けられる基準を満たすか否かを判定する例を示した。しかしながら、上記説明における「温度」は、「物理量」として適宜読み替え可能である。例えば、物品の保管区画の湿度に関する基準が定められている場合に、計測装置１０は、物品の保管区画の現在湿度または予想湿度の少なくとも１つを算出し、基準が満たされているか否かを判定してもよい。

上記の実施形態では、計測装置１０が、空調装置４０の推奨設定を報知する例を示した。図示されないコントローラが、計測装置１０によって報知された推奨設定に基づいて、空調装置４０の動作を制御してもよい。

上記の変形例４では、計測装置１０が、保管区画の温度をシミュレートする例を示した。空調装置４０の動作を制御する例、および空調装置４０の推奨設定を報知する例を示した。計測装置１０は、空調装置４０に対する制御内容、または推奨設定の少なくとも１つを、以下の少なくとも１つを参照して決定してもよい。
・対象空間ＳＰの少なくとも１つの区画の物理量
・対象空間ＳＰの内部環境に関する動的なパラメータ（例えば、温度計５０によって測定された基準温度、空調装置４０の熱量、対象空間ＳＰ内に存在する人の数、位置または動きに関する情報）
・対象空間ＳＰの内部環境に関する静的なパラメータ（例えば、対象空間ＳＰの３次元形状、対象空間ＳＰと外界とを隔てる構造体（例えば、壁、床、および天井）の断熱性能（例えば、素材および厚み））
・対象空間ＳＰの外部環境に関する動的または静的なパラメータ（例えば、日照条件、外気温、気圧）

音波送信装置２０は、自己相関が比較的強い自己相関信号（例えば、Ｍ系列信号、Goldコードなど）を含む超音波ビームを送信しても良い。これにより、空間の温度の計測結果のＳ／Ｎ比を更に向上させることができる。

音波送信装置２０が個別に異なる自己相関信号を含む超音波ビームを送信することにより、音波受信装置３０が、超音波ビームの発信源となる音波送信装置２０を識別しても良い。
また、音波送信装置２０毎に異なる発振周波数を有する超音波ビームを送信することにより、音波受信装置３０が、超音波ビームの発信源となる音波送信装置２０を識別しても良い。

本実施形態では、音波送信装置２０及び音波受信装置３０を区別して規定したが、本実施形態の範囲は、これに限られない。本実施形態は、１つの超音波振動子が超音波を送信する機能及び超音波を受信する機能を備えても良い。

本実施形態では、ステップＳ１１１４（図１７）において経路温度ＴＥＭＰｐａｔｈｉの計算に用いる式、及び、ステップＳ１１２（図１６）においてメッシュ温度ＴＥＭＰｍｅｓｈｔの計算に用いる式の少なくとも１つは、外部環境情報（例えば、外気温、外気の湿度、及び、外気圧の少なくとも１つ）をパラメータとして含んでも良い。この場合、外部環境情報に関わらず、空間の空気特性の計測結果のＳ／Ｎ比を向上させることができる。

本実施形態では、音波送信装置２０は、指向性を有する超音波ビームを送信する例を示したが、本実施形態は、これに限られない。本実施形態は、音波送信装置２０が可聴音ビーム（つまり、超音波ビームとは異なる周波数を有する音波）を送信する場合にも適用可能である。

本実施形態において、温度分布とは、メッシュ温度ＴＥＭＰｍｅｓｈに限られない。温度分布は、以下の少なくとも１つも含む。
・経路上の複数点の温度
・経路上の平均温度

以上、本発明の実施形態について詳細に説明したが、本発明の範囲は上記の実施形態に限定されない。また、上記の実施形態は、本発明の主旨を逸脱しない範囲において、種々の改良や変更が可能である。また、上記の実施形態及び変形例は、組合せ可能である。

（５）付記
実施形態および変形例で説明した事項を、以下に付記する。

（付記１）
少なくとも１つの対象物品を収容する空間に設置された音波送信装置（２０）と、空間に設置された音波受信装置（３０）とに接続可能な計測装置（１０）であって、
音波送信装置に音波を送信させる手段（Ｓ１１１１）と、
音波受信装置による音波の受信結果を参照して空間を仮想的に分割した複数の区画に亘る物理量の分布を算出する手段（Ｓ１１）と、
対象物品の位置に対応する区画である保管区画の物理量が対象物品に関連付けられる基準を満たすか否かを判定する手段（Ｓ１２，Ｓ２３）と、
保管区画の物理量が基準を満たさないと判定された場合に、所定の動作を行う手段（Ｓ１３，Ｓ２４）と
を具備する、計測装置。

（付記２）
算出する手段は、音波受信装置による音波の受信結果を参照して複数の区画に亘る温度の分布を算出し、
判定する手段は、算出された保管区画の温度が閾値を超える場合に保管区画の物理量が基準を満たさないと判定する、
付記１に記載の計測装置。

（付記３）
将来に亘る保管区画の予想温度を推定する手段（Ｓ２２）をさらに具備し、
算出する手段は、音波受信装置による音波の受信結果を参照して複数の区画に亘る温度を含む複数種の物理量の分布を算出し、
推定する手段は、複数種の物理量の分布を参照して保管区画の予想温度を推定し、
判定する手段（Ｓ２３）は、保管区画の予想温度が閾値を超える場合に保管区画の予想温度が基準を満たさないと判定する、
付記１に記載の計測装置。

（付記４）
算出する手段によって算出される複数種の物理量は、風速または風向の少なくとも１つを含む、
付記３に記載の計測装置。

（付記５）
将来に亘る保管区画の予想温度を推定する手段（Ｓ２２）をさらに具備し、
算出する手段は、音波受信装置による音波の受信結果を参照して複数時点に亘る保管区画の温度を算出し、
推定する手段は、複数時点に亘る保管区画の温度を参照して保管区画の予想温度を推定し、
判定する手段（Ｓ２３）は、保管区画の予想温度が閾値を超える場合に保管区画の予想温度が基準を満たさないと判定する、
付記１に記載の計測装置。

（付記６）
算出する手段は、複数時点に亘る保管区画の温度に対して補間処理を行うことにより保管区画の予測温度を算出する、
付記５に記載の計測装置。

（付記７）
保管区画の予測温度の時間変化に関する情報を出力する手段をさらに具備する、
付記３乃至付記６のいずれかに記載の計測装置。

（付記８）
所定の動作は、基準に関するアラートの報知である、
付記１乃至付記７のいずれかに記載の計測装置。

（付記９）
所定の動作は、空調装置の制御である、
付記１乃至付記８のいずれかに記載の計測装置。

（付記１０）
所定の動作は、推奨する空調設定の報知である、
付記１乃至付記９のいずれかに記載の計測装置。

（付記１１）
少なくとも１つの対象物品を収容する空間に設置された音波送信装置（２０）と、空間に設置された音波受信装置（３０）とに接続可能な計測装置（１０）によって行われる計測方法であって、計測装置が、
音波送信装置に音波を送信させること（Ｓ１１１１）と、
音波受信装置による音波の受信結果を参照して空間を仮想的に分割した複数の区画に亘る物理量の分布を算出すること（Ｓ１１）と、
対象物品の位置に対応する区画である保管区画の物理量が対象物品に関連付けられる基準を満たすか否かを判定すること（Ｓ１２，Ｓ２３）と、
保管区画の物理量が基準を満たさないと判定された場合に、所定の動作を行うこと（Ｓ１３，Ｓ２４）と
を具備する、計測方法。

（付記１２）
コンピュータに、付記１～付記１０の何れかに記載の各手段を実現させるためのプログラム。

１ ：保管状態管理システム
１０ ：計測装置
１１ ：記憶装置
１２ ：プロセッサ
１３ ：入出力インタフェース
１４ ：通信インタフェース
２０ ：音波送信装置
２１ ：超音波振動子
２２ ：制御回路
３０ ：音波受信装置
３１ ：超音波振動子
３２ ：制御回路
４０ ：空調装置
５０ ：温度計

Claims (12)

1. 少なくとも１つの対象物品を収容する空間に設置された音波送信装置と、前記空間に設置された音波受信装置とに接続可能な計測装置であって、
   前記音波送信装置に音波を送信させる手段と、
   前記音波受信装置による音波の受信結果を参照して前記空間を仮想的に分割した複数の区画に亘る物理量の分布を算出する手段と、
   前記対象物品の位置に対応する区画である保管区画の物理量が前記対象物品に関連付けられる基準を満たすか否かを判定する手段と、
   前記保管区画の物理量が前記基準を満たさないと判定された場合に、所定の動作を行う手段と
   を具備する、計測装置。
2. 前記算出する手段は、前記音波受信装置による音波の受信結果を参照して前記複数の区画に亘る温度の分布を算出し、
   前記判定する手段は、算出された前記保管区画の温度が閾値を超える場合に前記保管区画の物理量が前記基準を満たさないと判定する、
   請求項１に記載の計測装置。
3. 将来に亘る前記保管区画の予想温度を推定する手段をさらに具備し、
   前記算出する手段は、前記音波受信装置による音波の受信結果を参照して前記複数の区画に亘る温度を含む複数種の物理量の分布を算出し、
   前記推定する手段は、前記複数種の物理量の分布を参照して前記保管区画の予想温度を推定し、
   前記判定する手段は、前記保管区画の予想温度が閾値を超える場合に前記保管区画の予想温度が前記基準を満たさないと判定する、
   請求項１に記載の計測装置。
4. 前記算出する手段によって算出される複数種の物理量は、風速または風向の少なくとも１つを含む、
   請求項３に記載の計測装置。
5. 将来に亘る前記保管区画の予想温度を推定する手段をさらに具備し、
   前記算出する手段は、前記音波受信装置による音波の受信結果を参照して複数時点に亘る前記保管区画の温度を算出し、
   前記推定する手段は、前記複数時点に亘る保管区画の温度を参照して前記保管区画の予想温度を推定し、
   前記判定する手段は、前記保管区画の予想温度が閾値を超える場合に前記保管区画の予想温度が前記基準を満たさないと判定する、
   請求項１に記載の計測装置。
6. 前記算出する手段は、前記複数時点に亘る保管区画の温度に対して補間処理を行うことにより前記保管区画の予測温度を算出する、
   請求項５に記載の計測装置。
7. 前記保管区画の予測温度の時間変化に関する情報を出力する手段をさらに具備する、
   請求項３乃至請求項６のいずれかに記載の計測装置。
8. 前記所定の動作は、前記基準に関するアラートの報知である、
   請求項１乃至請求項７のいずれかに記載の計測装置。
9. 前記所定の動作は、空調装置の制御である、
   請求項１乃至請求項８のいずれかに記載の計測装置。
10. 前記所定の動作は、推奨する空調設定の報知である、
    請求項１乃至請求項９のいずれかに記載の計測装置。
11. 少なくとも１つの対象物品を収容する空間に設置された音波送信装置と、前記空間に設置された音波受信装置とに接続可能な計測装置によって行われる計測方法であって、前記計測装置が、
    前記音波送信装置に音波を送信させることと、
    前記音波受信装置による音波の受信結果を参照して前記空間を仮想的に分割した複数の区画に亘る物理量の分布を算出することと、
    前記対象物品の位置に対応する区画である保管区画の物理量が前記対象物品に関連付けられる基準を満たすか否かを判定することと、
    前記保管区画の物理量が前記基準を満たさないと判定された場合に、所定の動作を行うことと
    を具備する、計測方法。
12. コンピュータに、請求項１～請求項１０の何れかに記載の各手段を実現させるためのプログラム。
    
    

Images