特開2023-141392ガス監視方法、ガス監視装置、ガス監視システム及びガス監視プログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023141392
(43)【公開日】2023-10-05
(54)【発明の名称】ガス監視方法、ガス監視装置、ガス監視システム及びガス監視プログラム
(51)【国際特許分類】
G06T 7/00 20170101AFI20230928BHJP
【ＦＩ】
G06T7/00 640
【審査請求】未請求
【請求項の数】13
【出願形態】ＯＬ
(21)【出願番号】P 2022047701
(22)【出願日】2022-03-24
(71)【出願人】
【識別番号】000006208
【氏名又は名称】三菱重工業株式会社
(74)【代理人】
【識別番号】110000785
【氏名又は名称】ＳＳＩＰ弁理士法人
(72)【発明者】
【氏名】神原 信幸
(72)【発明者】
【氏名】荒川 宜彬
(72)【発明者】
【氏名】乾 正幸
(72)【発明者】
【氏名】渋谷 秀和
(72)【発明者】
【氏名】川添 浩平
(72)【発明者】
【氏名】森竹 崇之
【テーマコード（参考）】
5L096
【Ｆターム（参考）】
5L096AA06
5L096BA02
5L096CA04
5L096DA02
5L096FA35
5L096HA04
5L096HA09
(57)【要約】
【課題】ＵＡＶに搭載された撮像装置で取得した画像データを用いて監視対象のガスを適切に監視可能なガス監視方法、ガス監視装置、ガス監視システム及びガス監視プログラムを提供する。
【解決手段】ガス監視方法は、ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップと、前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップと、前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップと、を備える。
【選択図】図３
【解決手段】ガス監視方法は、ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップと、前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップと、前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップと、を備える。
【選択図】図３
【特許請求の範囲】
【請求項1】
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップと、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップと、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップと、
を備えるガス監視方法。
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップと、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップと、
を備えるガス監視方法。
【請求項2】
前記ＵＡＶの飛行状態を示す前記少なくとも１つの変数は、前記ＵＡＶの飛行高さＨ及び前記ＵＡＶの飛行速度ｖを含む
請求項１に記載のガス監視方法。
請求項１に記載のガス監視方法。
【請求項3】
地表における風速Ｗｇを取得するステップと、
前記ダウンウォッシュ速度ｕ及び前記風速Ｗｇに基づいて、前記風速Ｗｇを考慮したダウンウォッシュ速度である補正ダウンウォッシュ速度ｕｇを算出するステップと、を備え、
前記抽出ステップでは、前記複数の速度ベクトルから、前記補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する
請求項１又は２に記載のガス監視方法。
前記ダウンウォッシュ速度ｕ及び前記風速Ｗｇに基づいて、前記風速Ｗｇを考慮したダウンウォッシュ速度である補正ダウンウォッシュ速度ｕｇを算出するステップと、を備え、
前記抽出ステップでは、前記複数の速度ベクトルから、前記補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する
請求項１又は２に記載のガス監視方法。
【請求項4】
前記ダウンウォッシュ速度ｕのベクトルと前記地表における風速Ｗｇのベクトルとを合成することで前記補正ダウンウォッシュ速度ｕｇのベクトルを算出する
請求項３に記載のガス監視方法。
請求項３に記載のガス監視方法。
【請求項5】
地表からの高さと風速との相関関係を取得するステップと、
前記ＵＡＶの飛行高さにおける風速Ｗ１を取得するステップと、を備え、
前記地表における風速Ｗｇを取得するステップでは、前記風速Ｗ１及び前記相関関係に基づいて前記地表における風速Ｗｇを取得する
請求項３又は４に記載のガス監視方法。
前記ＵＡＶの飛行高さにおける風速Ｗ１を取得するステップと、を備え、
前記地表における風速Ｗｇを取得するステップでは、前記風速Ｗ１及び前記相関関係に基づいて前記地表における風速Ｗｇを取得する
請求項３又は４に記載のガス監視方法。
【請求項6】
前記風速Ｗ１を取得するステップでは、前記ＵＡＶのモータ電流値に基づいて前記ＵＡＶの飛行高さにおける前記風速Ｗ１を算出する
請求項５に記載のガス監視方法。
請求項５に記載のガス監視方法。
【請求項7】
前記撮像装置は、監視対象のガスが吸収する波長の赤外線を選択的に透過させるフィルタを備えた赤外線カメラを含む
請求項１乃至６の何れか一項に記載のガス監視方法。
請求項１乃至６の何れか一項に記載のガス監視方法。
【請求項8】
前記複数の画像データをオプティカルフロー推定処理して得られる前記複数の速度ベクトルから、前記複数の速度ベクトルの大きさの度数分布又は前記複数の速度ベクトルの向きの度数分布に基づいてノイズ成分を除去するステップを備え、
前記抽出ステップでは、前記ノイズ成分が除去された前記複数の速度ベクトルから前記ガス速度ベクトルを抽出する
請求項１乃至７の何れか一項に記載のガス監視方法。
前記抽出ステップでは、前記ノイズ成分が除去された前記複数の速度ベクトルから前記ガス速度ベクトルを抽出する
請求項１乃至７の何れか一項に記載のガス監視方法。
【請求項9】
前記複数の速度ベクトルの大きさの度数分布において、度数が第１閾値以上の大きさの範囲に属し、かつ、前記複数の速度ベクトルの向きの度数分布において、度数が第２閾値以下の向きの範囲に属する速度ベクトルを、前記複数の速度ベクトルから前記ノイズ成分として除去する
請求項８に記載のガス監視方法。
請求項８に記載のガス監視方法。
【請求項10】
前記ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データを取得するステップと、
前記複数の画像データについてリフレーム処理を施すステップと、
前記リフレーム処理が施された複数の画像データをオプティカルフロー推定処理して前記複数の速度ベクトルを取得するステップと、
を備える請求項１乃至９の何れか一項に記載のガス監視方法。
前記複数の画像データについてリフレーム処理を施すステップと、
前記リフレーム処理が施された複数の画像データをオプティカルフロー推定処理して前記複数の速度ベクトルを取得するステップと、
を備える請求項１乃至９の何れか一項に記載のガス監視方法。
【請求項11】
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する速度ベクトル取得部と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するダウンウォッシュ速度取得部と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出部と、
を備えるガス監視装置。
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するダウンウォッシュ速度取得部と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出部と、
を備えるガス監視装置。
【請求項12】
撮像装置を搭載したＵＡＶと、
前記撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して複数の速度ベクトルを算出するように構成されたオプティカルフロー推定処理部と、
前記複数の速度ベクトルから前記ガス速度ベクトルを抽出するように構成された請求項１１に記載のガス監視装置と、
を備えるガス監視システム。
前記撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して複数の速度ベクトルを算出するように構成されたオプティカルフロー推定処理部と、
前記複数の速度ベクトルから前記ガス速度ベクトルを抽出するように構成された請求項１１に記載のガス監視装置と、
を備えるガス監視システム。
【請求項13】
コンピュータに、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する手順と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得する手順と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する手順と、
実行させるためのガス監視プログラム。
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する手順と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得する手順と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する手順と、
実行させるためのガス監視プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本開示は、ガス監視方法、ガス監視装置、ガス監視システム及びガス監視プログラムに関する。
【背景技術】
【0002】
カメラなどの撮像装置で撮影された画像データに基づいてガスを監視したり監視したりすることがある。
【0003】
特許文献１には、オプティカルフロー推定を用いた黒煙検知システムが開示されている。この黒煙検知システムでは、フレアスタックから排出される黒煙に対して、連続して撮影された時系列の２画像を用いて、両画像間で演算したオプティカルフロー推定からガスの移動で生じた画像内局所の速度ベクトルが求められる。このようにして求められた速度ベクトルのうち、大きさ及び向きが所定範囲内であるものを、黒煙の移動を示す速度ベクトルとして抽出している。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特許第４２６６５３５号公報
【発明の概要】
【発明が解決しようとする課題】
【0005】
ところで、特許文献１の黒煙検知システムでは、固定された撮像装置（カメラ）を用いているため、撮像装置と撮影対象（煙突等）との距離が一定であり、また、撮影範囲内に速度算出の基準となるスケール（煙突等）が存在するため、オプティカルフロー推定から得られる複数の速度ベクトルの大きさ（速度の絶対値）を算出することができる。
【0006】
一方、ＵＡＶ（無人航空機；unmanned aerial vehicle）に搭載された撮像装置を用いる場合、ＵＡＶの高さや位置は変化するためＵＡＶと撮影対象の距離が一定でなく、また、速度算出の基準となるスケールが画像内に存在するとは限らないため、該複数の速度ベクトルの各々の大きさ（速度の絶対値）はオプティカルフロー推定では演算できず未知である。このため、オプティカルフロー推定で得られる複数の速度ベクトルから、監視対象のガスの移動を示すベクトルを適切に抽出することが難しい。
【0007】
上述の事情に鑑みて、本発明の少なくとも一実施形態は、ＵＡＶに搭載された撮像装置で取得した画像データを用いて監視対象のガスを適切に監視可能なガス監視方法、ガス監視装置、ガス監視システム及びガス監視プログラムを提供することを目的とする。
【課題を解決するための手段】
【0008】
本発明の少なくとも一実施形態に係るガス監視方法は、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップと、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップと、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップと、
を備える。
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップと、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップと、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップと、
を備える。
【0009】
また、本発明の少なくとも一実施形態に係るガス監視装置は、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する速度ベクトル取得部と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するダウンウォッシュ速度取得部と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出部と、
を備える。
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する速度ベクトル取得部と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するダウンウォッシュ速度取得部と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出部と、
を備える。
【0010】
また、本発明の少なくとも一実施形態に係るガス監視システムは、
撮像装置を搭載したＵＡＶと、
前記撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して複数の速度ベクトルを算出するように構成されたオプティカルフロー推定処理部と、
前記複数の速度ベクトルから前記ガス速度ベクトルを抽出するように構成された上述のガス監視装置と、
を備える。
撮像装置を搭載したＵＡＶと、
前記撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して複数の速度ベクトルを算出するように構成されたオプティカルフロー推定処理部と、
前記複数の速度ベクトルから前記ガス速度ベクトルを抽出するように構成された上述のガス監視装置と、
を備える。
【0011】
また、本発明の少なくとも一実施形態に係るガス監視プログラムは、
コンピュータに、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する手順と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得する手順と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する手順と、
実行させる。
コンピュータに、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する手順と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得する手順と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する手順と、
実行させる。
【発明の効果】
【0012】
本発明の少なくとも一実施形態によれば、ＵＡＶに搭載された撮像装置で取得した画像データを用いて監視対象のガスを適切に監視可能なガス監視方法、ガス監視装置、ガス監視システム及びガス監視プログラムが提供される。
【図面の簡単な説明】
【0013】
【図1】一実施形態に係るガス監視システムを構成するＵＡＶの概略図である。
【図2】一実施形態に係るガス監視システムの概略図である。
【図3】一実施形態に係るガス監視方法のフローチャートである。
【図4】撮像装置で撮影された画像の一例を示す図である。
【図5】撮像装置で撮影された画像にオプティカルフロー推定により算出された複数の速度ベクトルを重ねた画像の一例である。
【図6】ＵＡＶの飛行高さＨ及び飛行速度ｖとＵＡＶのダウンウォッシュ速度との相関関係の一例を示すグラフである。
【図7】ＵＡＶの飛行高さＨ及び飛行速度ｖとＵＡＶのダウンウォッシュ速度との相関関係の一例を示すチャートである。
【図8】ＵＡＶのダウンウォッシュ速度の取得手順を説明するための図である。
【図9】ＵＡＶの補正ダウンウォッシュ速度の取得手順を説明するための図である。
【図10】地表からの高さとの風速との相関関係の一例を示すグラフである。
【図11】撮像装置で撮影された画像に、抽出ステップで抽出されたガス速度ベクトルＶＧを重ねた画像の一例である。
【図12】ガス速度ベクトルＶＧの抽出の手順を説明するためのグラフである。
【図13】複数の速度ベクトルの大きさの度数分布の一例を示すヒストグラムである。
【図14】複数の速度ベクトルの向きの度数分布の一例を示すヒストグラムである。
【図15】画像データのリフレーム処理を説明するための図である。
【図16】画像データのリフレーム処理を説明するための図である。
【図17】画像データのリフレーム処理を説明するための図である。
【発明を実施するための形態】
【0014】
以下、添付図面を参照して本発明の幾つかの実施形態について説明する。ただし、実施形態として記載されている又は図面に示されている構成部品の寸法、材質、形状、その相対的配置等は、本発明の範囲をこれに限定する趣旨ではなく、単なる説明例にすぎない。
【0015】
【0016】
幾つかの実施形態に係るガス監視装置及びガス監視システムは、ＵＡＶ（無人航空機；unmanned aerial vehicle）に搭載された撮像装置（カメラ等）でフィールドを撮影して得られる画像データに基づいて、画像内のガス領域（監視対象のガスが流れている領域）を特定することが可能な装置及びシステムである。このガス監視装置／ガス監視システムにより、撮影対象となるフィールド領域における監視対象ガスの検知や監視（例えば、施設からのガス漏れ検知や、ガス流出の監視等）をすることができる。なお、監視対象のガスの種類は特に限定されない。以下においては、一例として、二酸化炭素（ＣＯ２）ガスを監視対象とする場合について説明する。
【0017】
図１及び図２に示すように、一実施形態に係るガス監視システム１００は、ＵＡＶ１０と、ＵＡＶ１０に搭載された撮像装置２０と、ＵＡＶ１０及び／又は撮像装置２０等で取得された計測データや画像データを示す信号を処理するためのガス監視装置５０と、を備える。
【0018】
図１に示すように、ＵＡＶ１０は、ＵＡＶ本体１２と、ＵＡＶ本体１２に取り付けられる複数のプロペラ１４（通常、３つ以上のプロペラ１４）と、を含む。複数のプロペラ１４の各々は、モータによって回転駆動されるように構成される。ＵＡＶ１０の飛行高さ、飛行速度、飛行方向等は、複数のプロペラ１４を駆動するための各モータの電流値を調節することによって制御可能になっている。各モータの電流値を示す信号は、無線通信によりガス監視装置５０に送られるようになっていてもよい。
【0019】
ＵＡＶ１０が飛行するとき、複数のプロペラ１４によってダウンウォッシュＤ（下方に向かう気流）が生成される。ＵＡＶ１０から下方に向かうダウンウォッシュＤは、地表Ｇに衝突して水平方向に転向され、地表におけるＵＡＶ１０の真下の位置Ｐｃから地表Ｇに沿って放射状に広がるように流れる。
【0020】
ＵＡＶ１０には、ＵＡＶ１０の飛行高さＨを計測するための高度計（不図示）、又は、ＵＡＶ１０の飛行速度ｖを計測するための速度計（不図示）が設けられていてもよい。高度計で取得されたＵＡＶ１０の飛行高さＨを示す信号、及び／又は、速度計で取得されたＵＡＶ１０の飛行速度ｖを示す信号は、無線通信によりガス監視装置５０に送られるようになっていてもよい。
【0021】
撮像装置２０は、ＵＡＶ１０が飛行している状態で、ＵＡＶ１０の下方のフィールドを時系列で連続的に撮影可能なものである。撮像装置２０によって撮像された画像を示す画像データは、無線通信によりガス監視装置５０に送られるようになっていてもよい。
【0022】
撮像装置２０は、監視対象のガスが吸収する波長（ＣＯ２ガスの場合４．３μｍ）の赤外線を選択的に透過させるフィルタを備えた赤外線カメラであってもよい。このように、特定の波長を選択的にフィルタリングしてイメージングする赤外線カメラを用いることで、撮影視野の流体（気体）については特定のガスのみが画像データにイメージングされる。
【0023】
ガス監視装置５０は、ＵＡＶ１０から送られた飛行高さＨ又は飛行速度ｖ等の計測データ又は撮像装置２０から送られた画像データを処理するように構成される。図２に示すように、一実施形態に係るガス監視装置５０は、画像取得部５２と、リフレーム処理部５４と、オプティカルフロー処理部５６と、速度ベクトル取得部５８と、ノイズ成分除去部６０と、飛行状態取得部６２と、ダウンウォッシュ速度取得部６４と、抽出部６６と、表示画像生成部６８と、を備える。
【0024】
ガス監視装置５０は、プロセッサ（ＣＰＵ又はＧＰＵ等）、記憶装置（メモリデバイス；ＲＡＭ等）、補助記憶部及びインターフェース等を備えた計算機を含む。ガス監視装置５０は、インターフェースを介して、ＵＡＶ１０又は撮像装置２０からの信号を受け取るようになっている。プロセッサは、このようにして受け取った信号を処理するように構成される。また、プロセッサは、記憶装置に展開されるプログラムを処理するように構成される。これにより、上述の各機能部（画像取得部５２～表示画像生成部６８）の機能が実現される。
【0025】
ガス監視装置５０での処理内容は、プロセッサにより実行されるプログラムとして実装される。プログラムは、補助記憶部に記憶されていてもよい。プログラム実行時には、これらのプログラムは記憶装置に展開される。プロセッサは、記憶装置からプログラムを読み出し、プログラムに含まれる命令を実行するようになっている。
【0026】
画像取得部５２は、ＵＡＶ１０に搭載された撮像装置２０で複数の時刻にそれぞれ撮像して得られる複数の画像データを取得するように構成される。
【0027】
リフレーム処理部５４は、必要に応じ、画像取得部５２によって取得された複数の画像データについてリフレーム処理を行うように構成される。
【0028】
オプティカルフロー処理部５６は、画像取得部５２によって取得され、必要に応じてリフレーム処理が施された複数の画像データについて、オプティカルフロー推定処理するように構成される。
【0029】
速度ベクトル取得部５８は、オプティカルフロー処理部５６でのオプティカルフロー処理で算出される複数の速度ベクトルを取得するように構成される。
【0030】
ノイズ成分除去部６０は、必要に応じ、速度ベクトル取得部５８で取得された複数の速度ベクトルからノイズ成分を除去するように構成される。
【0031】
飛行状態取得部６２は、ＵＡＶ１０の飛行状態を示す少なくとも１つの変数の計測値を取得するように構成される。一実施形態では、飛行状態取得部６２は、ＵＡＶ１０の飛行状態を示す変数としてのＵＡＶ１０の飛行高さＨ及び／又はＵＡＶ１０の飛行速度ｖの計測値を取得するように構成される。
【0032】
ダウンウォッシュ速度取得部６４は、ＵＡＶ１０の飛行状態を示す少なくとも１つの変数に基づきＵＡＶ１０のダウンウォッシュ速度ｕを取得するように構成される。このダウンウォッシュ速度ｕは、地表におけるＵＡＶ１０の真下の位置Ｐｃの近傍における、地表Ｇに沿った方向のダウンウォッシュの流れの速度である。
【0033】
抽出部６６は、速度ベクトル取得部５８で取得された複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出するように構成される。
【0034】
表示画像生成部６８は、抽出部６６で抽出されたガス速度ベクトルに基づいて、監視対象のガスに関する情報（例えば、ガスの存在するガス領域や、ガス速度ベクトルの分布等に係る情報）を示す画像データ（監視用の画像データ）を生成するように構成される。表示画像生成部６８で生成された画像データは、画像を表示するための表示部７０に出力されるようになっていてもよい。
【0035】
また、ガス監視装置５０は、ＵＡＶ１０の飛行状態を示す変数とＵＡＶ１０のダウンウォッシュ速度との相関関係や、地表からの高さと風速との相関関係等を記憶するための記憶部６９を備えてもよい。記憶部６９は、ガス監視装置５０を構成する計算機の記憶装置（メモリデバイス；ＲＡＭ等）又は補助記憶装置を含んでもよく、あるいは、該計算機とネットワークを介して接続されるストレージデバイスを含んでもよい。
【0036】
（ガス監視方法のフロー）
以下、幾つかの実施形態に係るガス監視方法のフローについて説明する。なお、以下においては、上述のガス監視装置５０を用いて監視対象のガスであるＣＯ２ガスを監視する場合について説明するが、以下に説明する手順の一部又は全部を、他の装置を用いて、あるいは手動で行ってもよく、監視対象のガスはＣＯ２以外のガスであってもよい。
以下、幾つかの実施形態に係るガス監視方法のフローについて説明する。なお、以下においては、上述のガス監視装置５０を用いて監視対象のガスであるＣＯ２ガスを監視する場合について説明するが、以下に説明する手順の一部又は全部を、他の装置を用いて、あるいは手動で行ってもよく、監視対象のガスはＣＯ２以外のガスであってもよい。
【0037】
図３は、一実施形態に係るガス監視方法のフローチャートである。
【0038】
一実施形態では、まず、ＵＡＶ１０に搭載された撮像装置２０で複数の時刻にそれぞれ撮像する（Ｓ１００）。撮像装置２０は、ＵＡＶ１０が飛行している状態で、ＵＡＶ１０の下方のフィールドを時系列で連続的に撮影してもよい。図４は、撮像装置２０で撮影された画像の一例を示す図である。図４に示す画像Ｉ０は、飛行中のＵＡＶ１０の下方のフィールドを撮影したものであり、地表Ｇと、地表Ｇに設置されたＣＯ２ガスボンベＢ１，Ｂ２が写っている。本実施形態では、撮像装置２０として、監視対象のＣＯ２ガスが吸収する波長（４．３μｍ）の赤外線を選択的に透過させるフィルタを備えた赤外線カメラを用いる。
【0039】
画像取得部５２は、このように撮像装置２０で撮像された複数の画像のデータ（複数の画像データ）を取得する。
【0040】
次に、リフレーム処理部５４は、必要に応じ、ステップＳ１００にて画像取得部５２によって取得された複数の画像データについてリフレーム処理を行ってもよい（Ｓ２００）。リフレーム処理をして予めＵＡＶ１０のふらつきによるノイズを除去してから、後述のステップＳ３００でオプティカルフロー推定処理することができるため、後述のステップＳ８００でのガス速度ベクトルの抽出精度をより向上することができる。なお、ステップＳ２００のリフレーム処理については後述する。
【0041】
次に、オプティカルフロー処理部５６は、ステップＳ１００にて画像取得部５２によって取得された複数の画像データについて、あるいは、ステップＳ２００にてリフレーム処理が施された複数の画像データについて、オプティカルフロー推定処理を行う（Ｓ３００）。オプティカルフロー推定では、時間的に連続した時系列の２つの画像データを用いて、この２つの画像データ間での、画像内の複数の画素の動きをそれぞれ示す複数の速度ベクトルが算出される。オプティカルフロー推定処理は、例えば特許文献１に記載の方法で行うことができる。図５に示す画像Ｉ１は、ステップＳ３００のオプティカルフロー推定処理の対象となった画像（撮像装置で撮影されたもの）に、オプティカルフロー推定により算出された複数の速度ベクトルを重ねた画像の一例である。
【0042】
なお、本実施形態では、ＣＯ２ガスが吸収する波長を選択的にフィルタリングしてイメージングする赤外線カメラを用いるため、撮影視野の流体（気体）についてはＣＯ２ガスのみが画像データにイメージングされるとともに、オプティカルフロー推定の演算対象となる。
【0043】
次に、速度ベクトル取得部５８は、ステップＳ３００でのオプティカルフロー推定処理で算出された複数の速度ベクトルを取得する（Ｓ４００）。
【0044】
次に、ノイズ成分除去部６０は、必要に応じ、ステップＳ４００で取得された複数の速度ベクトルからノイズ成分を除去する（Ｓ５００）。これにより、後述のステップＳ８００でのガス速度ベクトルの抽出精度を向上することができる。なお、ステップＳ５００のノイズ成分除去については後述する。
【0045】
次に、飛行状態取得部６２は、ＵＡＶ１０の飛行状態を示す少なくとも１つの変数の計測値を取得する（Ｓ６００）。ステップＳ６００では、ＵＡＶ１０の飛行状態を示す変数としてのＵＡＶ１０の飛行高さＨ及び／又はＵＡＶ１０の飛行速度ｖの計測値を取得してもよい。
【0046】
次に、ダウンウォッシュ速度取得部６４は、ステップＳ６００で取得されたＵＡＶ１０の飛行状態を示す少なくとも１つの変数に基づきＵＡＶ１０のダウンウォッシュ速度ｕを取得する（Ｓ７００）。ＵＡＶ１０の飛行状態を示す上述の少なくとも１つの変数と、ＵＡＶ１０のダウンウォッシュ速度ｕとの相関関係は、予め取得され記憶部６９に記憶されていてもよい。ステップＳ７００において、ダウンウォッシュ速度取得部６４は、ステップＳ６００で取得されたＵＡＶ１０の飛行状態を示す少なくとも１つの変数に係る計測値と、記憶部６９から取得される上述の相関関係とに基づいて、該計測値に対応するダウンウォッシュ速度ｕを取得してもよい。
【0047】
上述のＵＡＶ１０の飛行状態を示す少なくとも１つの変数は、ＵＡＶ１０の飛行高さＨ及びＵＡＶ１０の飛行速度ｖを含んでもよい。ここで、図６及び図７は、それぞれ、ＵＡＶ１０の飛行高さＨ及びＵＡＶ１０の飛行速度ｖと、ＵＡＶ１０のダウンウォッシュ速度ｕとの相関関係を示すグラフ又はチャートの一例である。
【0048】
一実施形態では、図６に示す相関関係に基づいてダウンウォッシュ速度ｕを取得してもよい。すなわち、ＵＡＶ１０の飛行高さＨ及びＵＡＶ１０の飛行速度ｖとＵＡＶ１０のダウンウォッシュ速度ｕとは、例えば図６に示すような所定の相関関係を有する。この相関関係によれば、同じ飛行速度ｖでは、飛行高さが高いほどダウンウォッシュ速度ｕは大きく、また、飛行速度が大きいほどダウンウォッシュ速度ｕは大きい（ただし、図６において飛行速度ｖはｖ０＜ｖ１＜ｖ２である）。この相関関係は、ＵＡＶ１０を用いた試験により予め取得し、記憶部６９に記憶させておく。そして、ダウンウォッシュ速度取得部６４は、ステップＳ６００で取得されたＵＡＶ１０の飛行高さＨ及び飛行速度ｖの計測値を、記憶部６９から取得される相関関係に当てはめることで、ＵＡＶ１０のダウンウォッシュ速度ｕを取得することができる。
【0049】
一実施形態では、図７に示す相関関係に基づいて位置Ｐ１におけるダウンウォッシュ速度ｕを取得してもよい。なお、位置Ｐ１は、ＵＡＶ１０の真下の位置Ｐｃの近傍の位置である。すなわち、例えば図７に示すように、ダウンウォッシュ速度ｕは、ＵＡＶ１０の飛行高さＨ及びＵＡＶ１０の飛行速度ｖに加え、地表におけるＵＡＶ１０真下の位置ＰｃからＰ１までの距離ｒ、及び、地表において位置Ｐｃを中心としたときの位置Ｐ１の基準方位からの角度θ（図８参照）と所定の相関関係を有する。この相関関係は、ＵＡＶ１０を用いた試験により予め取得し、記憶部６９に記憶させておく。そして、ダウンウォッシュ速度取得部６４は、ステップＳ６００で取得されたＵＡＶ１０の飛行高さＨ及び飛行速度ｖの計測値、及び、上述の距離ｒ及び角度θの計測値を、記憶部６９から取得される相関関係に当てはめることで、ＵＡＶ１０のダウンウォッシュ速度ｕを取得することができる。なお、上述の距離ｒ及び角度θの計測値は、測定者が入力装置（キーボード、マウス、又はタッチパネル等）を介してガス監視装置５０に入力してもよい。
【0050】
なお、図８は、ダウンウォッシュ速度ｕの取得手順を説明するための図であり、地表におけるＵＡＶ１０真下の位置ＰｃからＰ１までの距離ｒ、及び、地表において位置Ｐｃｗｏ中心としたときの位置Ｐ１の基準方位からの角度θを示す図である。
【0051】
幾つかの実施形態では、ステップＳ７００では、地表Ｇにおける風速Ｗｇを考慮したダウンウォッシュ速度である補正ダウンウォッシュ速度ｕｇを取得してもよい。補正ダウンウォッシュ速度ｕｇは、上述のようにして求めたダウンウォッシュ速度ｕ及び地表における風速Ｗｇに基づいて算出してもよい。
【0052】
補正ダウンウォッシュ速度ｕｇ（ベクトル）は、例えば、上述のようにして求めたダウンウォッシュ速度ｕ（ベクトル）と、地表における風速Ｗｇ（ベクトル）との合成ベクトルとして取得することができる（図９参照）。ここで、図９は、補正ダウンウォッシュ速度ｕｇの取得手順を説明するための図である。
【0053】
地表における風速Ｗｇ（ベクトル）は、ＵＡＶ１０の飛行高さにおける風速Ｗ１（ベクトル）と、地表からの高さと風速との相関関係とに基づいて取得してもよい。ここで、図１０は、地表からの高さとの風速との相関関係の一例を示すグラフである。地表からの高さとの風速とは、風速補正係数αを介して、例えば図１０に示すような相関関係を有する。すなわち、図１０に示す例では、風速補正係数α（ただし地表における風速補正係数αをゼロとする）は、ＵＡＶ１０の飛行高さに対して線形の関数で表現される。この関数から、地表での風速Ｗｇと、飛行高さＨ１における風速Ｗ１との関係は、飛行高さＨ１における風速補正係数α１とすれば、Ｗｇ＝α１×Ｗ１で表すことができる。地表からの高さと風速との上述の相関関係は、予め取得され、記憶部６９に記憶されていてもよい。
【0054】
地表からの高さと風速との上述の相関関係は、地表（フィールド）上における風向とＵＡＶ１０の飛行高さにおける風向は同じであり、かつ、高さ方向における風速変化の割合は地表（フィールド）上の場所によらず同じであるとの仮定の下、
フィールドの任意地点でＵＡＶ１０の飛行高さＨにおける風速Ｗを複数点取得することで、取得してもよい。
フィールドの任意地点でＵＡＶ１０の飛行高さＨにおける風速Ｗを複数点取得することで、取得してもよい。
【0055】
ＵＡＶ１０の飛行高さにおける風速Ｗ１（ベクトル）は、ＵＡＶ１０のモータ電流値（プロペラ１４を回転駆動するためのモータの電流値）から取得されてもよい。ＵＡＶ１０は、指定された飛行速度（及び方向）に対して、各プロペラ１４のモータ電流値を制御して、一定の飛行速度を保つことができ、周囲に風がある場合でも、各プロペラ１４のモータ負荷を制御して飛行速度一定に保つことができる。したがって、ＵＡＶ１０の各プロペラ１４のモータ電流値から風速（風速Ｗ１（ベクトル）の大きさ）及び風向（風速Ｗ１（ベクトル）の向き）を特定することができる。
【0056】
あるいは、ＵＡＶ１０の飛行高さにおける風速Ｗ１は、ＵＡＶ１０に設けられた風速及び風向を計測するように構成された風速計（不図示）及び風向計（不図示）の計測結果に基づき取得されてもよい。なお、風速計及び風向計の計測結果を示す信号がガス監視装置５０に送られるようになっていてもよい。
【0057】
すなわち、ダウンウォッシュ速度取得部６４は、ＵＡＶ１０のモータ電流値を取得し、該モータ電流値に基づいてＵＡＶ１０の飛行高さにおける風速Ｗ１（ベクトル）を算出してもよい。あるいは、ダウンウォッシュ速度取得部６４は、ＵＡＶ１０に設けられた風速計及び風向計の計測結果を取得し、該計測値に基づいてＵＡＶ１０の飛行高さにおける風速Ｗ１（ベクトル）を算出してもよい。そして、ダウンウォッシュ速度取得部６４は、記憶部６９から地表からの高さと風速との相関関係を取得し、該相関関係と、上述のＵＡＶ１０の飛行高さにおける風速Ｗ１とに基づいて、地表における風速Ｗｇを取得してもよい。ダウンウォッシュ速度取得部６４は、このように取得した地表における風速Ｗｇと、上述のダウンウォッシュ速度ｕとに基づいて、補正ダウンウォッシュ速度ｕｇを取得してもよい。
【0058】
次に、抽出部６６は、ステップＳ４００で取得した複数の速度ベクトル（オプティカルフロー推定処理の結果得られた複数の速度ベクトル；図５参照）から、ステップＳ７００で取得したダウンウォッシュ速度ｕ（又は補正ダウンウォッシュ速度ｕｇ）に基づいて決定される範囲内の大きさを有するガス速度ベクトルＶＧを抽出する（Ｓ８００）。なお、ステップＳ５００でノイズ成分を除去した場合は、ステップＳ８００では、ステップＳ５００にてノイズ成分が除去された複数の速度ベクトルから、ステップＳ７００で取得したダウンウォッシュ速度ｕ（又は補正ダウンウォッシュ速度ｕｇ）に基づいて決定される範囲内の大きさを有するガス速度ベクトルＶＧを抽出する。ここで、図１１に示す画像Ｉ２は、ステップＳ３００のオプティカルフロー推定処理の対象となった画像（撮像装置で撮影されたもの）に、ステップＳ８００で抽出されたガス速度ベクトルＶＧを重ねた画像の一例である。
【0059】
一実施形態では、ステップＳ８００では、例えば、ステップＳ４００で取得した複数の速度ベクトルのうち、ステップＳ７００で取得したダウンウォッシュ速度ｕ（又は補正ダウンウォッシュ速度ｕｇ）の位置における速度ベクトルについて、その速度の大きさがダウンウォッシュ速度ｕ（又は補正ダウンウォッシュ速度ｕｇ）の速度の大きさと等しいものと見做す。そして、複数の速度ベクトルのうち、ダウンウォッシュ速度ｕ（又は補正ダウンウォッシュ速度ｕｇ）に近い速度の大きさを有するものを抽出する。
【0060】
ガス速度ベクトルＶＧの抽出の仕方の一例として、ここでは、補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルＶＧを抽出する場合について説明する。一例においては、複数の速度ベクトルから、補正ダウンウォッシュ速度ｕｇと同じ大きさの速度ベクトル（ＶＧ＿ｍａｘ）と、補正ダウンウォッシュ速度ｕｇよりも小さい速度ｕｍｉｎと同じ大きさの速度ベクトル（ＶＧ＿ｍｉｎ）との間の大きさを有する速度ベクトルを、ガス速度ベクトルＶＧとして抽出する（図１２参照）。例えば、複数の速度ベクトルのうち、０．２５×ｕｇ以上、かつ、１×ｕｇ以下の大きさを有するものを抽出する。なお、図１２は、ガス速度ベクトルＶＧの抽出の手順を説明するためのグラフである。
【0061】
なお、図１１に示す画像Ｉ２においては、抽出されたガス速度ベクトルＶＧの分布、及び、抽出されたガス速度ベクトルＶＧが存在する領域であるガス領域ＡＧが示されている。また、ガス速度ベクトルＶＧは、ＣＯ２ガスボンベから噴出されるＣＯ２ガスの流れを示すベクトルである。
【0062】
ステップＳ８００で抽出されたガス速度ベクトルＶＧに基づいて、表示画像生成部６８は、抽出部６６で抽出されたガス速度ベクトルに基づいて、監視対象のガスに関する情報（例えば、ガスの存在するガス領域や、ガス速度ベクトルの分布等に係る情報）を示す画像データ（監視用の画像データ；例えば図１１に示す画像Ｉ２のデータ等）を生成してもよい。このように生成された画像データを、表示部７０に送って表示部７０にて出力する（表示する）ようにしてもよい（Ｓ９００）。
【0063】
ステップＳ３００のオプティカルフロー推定処理では、複数の速度ベクトルの向き及び相対的な大きさを算出することができるが、絶対的な大きさ（速度の絶対値）を算出することはできない。これは、ＵＡＶ１０に搭載される撮像装置２０と撮影対象の距離が一定でなく、さらに速度算出の基準となるスケールが画像内に存在するとは限らないためである。
この点、上述の実施形態では、ＵＡＶ１０の飛行状態を示す少なくとも１つの変数（例えば飛行高さＨや飛行速度ｖ）に基づき、ＵＡＶ１０のダウンウォッシュ速度ｕを取得することができる。そこで、オプティカルフロー推定で得られた未知の大きさを有する複数の速度ベクトルのうち、ＵＡＶ１０のダウンウォッシュの速度ｕを示すベクトルの大きさを基準にして、該複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づき決定される範囲内の大きさ（又は、ダウンウォッシュ速度ｕから求まる補正ダウンウォッシュ速度ｕｇに基づき決定される範囲内の大きさ）を有する速度ベクトルを監視対象のガスの移動速度及び方向を示すガス速度ベクトルＶＧ（ガス領域）として特定することができる。このように抽出されたガス速度ベクトルＶＧを用いて、ガスを適切に監視することができる。
この点、上述の実施形態では、ＵＡＶ１０の飛行状態を示す少なくとも１つの変数（例えば飛行高さＨや飛行速度ｖ）に基づき、ＵＡＶ１０のダウンウォッシュ速度ｕを取得することができる。そこで、オプティカルフロー推定で得られた未知の大きさを有する複数の速度ベクトルのうち、ＵＡＶ１０のダウンウォッシュの速度ｕを示すベクトルの大きさを基準にして、該複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づき決定される範囲内の大きさ（又は、ダウンウォッシュ速度ｕから求まる補正ダウンウォッシュ速度ｕｇに基づき決定される範囲内の大きさ）を有する速度ベクトルを監視対象のガスの移動速度及び方向を示すガス速度ベクトルＶＧ（ガス領域）として特定することができる。このように抽出されたガス速度ベクトルＶＧを用いて、ガスを適切に監視することができる。
【0064】
また、一実施形態では、上述したように、ＵＡＶ１０のダウンウォッシュ速度ｕは、ＵＡＶの飛行高さＨや飛行速度ｖ及び相関関係を有する。したがって、ＵＡＶ１０の飛行高さＨ及びＵＶＡの飛行速度ｖに基づいて、これらとの相関関係から、ＵＡＶ１０のダウンウォッシュ速度ｕを適切に取得することができる。
【0065】
また、地表においてダウンウォッシュとは別の風が吹いてる場合、地表におけるダウンウォッシュの流れは該風の影響を受け得る。この点、一実施形態では、上述したように、ＵＡＶ１０の飛行状態を示す変数（飛行高さＨや飛行速度ｖ）に基づき得られたダウンウォッシュ速度ｕ及び地表における風速Ｗｇに基づいて風速Ｗｇを考慮した補正ダウンウォッシュ速度ｕｇを算出する。そして、オプティカルフロー推定で得られた複数のベクトルから、該補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルＶＧを抽出する。したがって、監視対象のガスの移動速度及び方向を示すガス速度ベクトルをより適切に抽出することができる。
【0066】
（ノイズ成分除去について）
ステップＳ５００のノイズ成分を除去する処理について説明する。幾つかの実施形態では、ステップＳ５００では、ステップＳ４００で得られた複数のベクトル（ステップＳ３００のオプティカルフロー推定処理で算出される複数のベクトル）から、該複数の速度ベクトルの大きさの度数分布、又は、該複数の速度ベクトルの向きの度数分布に基づいて、ノイズ成分の速度ベクトルを除去してもよい。
ステップＳ５００のノイズ成分を除去する処理について説明する。幾つかの実施形態では、ステップＳ５００では、ステップＳ４００で得られた複数のベクトル（ステップＳ３００のオプティカルフロー推定処理で算出される複数のベクトル）から、該複数の速度ベクトルの大きさの度数分布、又は、該複数の速度ベクトルの向きの度数分布に基づいて、ノイズ成分の速度ベクトルを除去してもよい。
【0067】
ノイズ成分の速度ベクトルは、大きさ又は向きに関して所定のパターンを有する傾向がある。上述の実施形態では、オプティカルフロー推定により得られる複数の速度ベクトルの大きさ又は向きの度数分布に基づいて、該複数のベクトルからノイズ成分を除去することができる。また、ステップＳ８００では、ノイズ成分が除去された複数のベクトルからダウンウォッシュ速度ｕに基づきガス速度ベクトルを抽出するので、ガス速度ベクトルＶＧの抽出精度を向上することができる。
【0068】
ここで、図１３は、ステップＳ４００で得られた複数の速度ベクトルの大きさの度数分布の一例を示すヒストグラムであり、図１４は、ステップＳ４００で得られた複数の速度ベクトルの向きの度数分布の一例を示すヒストグラムである。
【0069】
一実施形態では、ステップＳ５００では、複数の速度ベクトルの大きさの度数分布（図１３参照）において、度数が第１閾値（図１３における閾値１）以上の大きさの範囲に属し、かつ、複数の速度ベクトルの向きの度数分布（図１４参照）において、度数が第２閾値（図１４における閾値２）以下の向きの範囲に属する速度ベクトルを、ノイズ成分としてステップＳ４００で得られた複数の速度ベクトルから除去する。
【0070】
ノイズ成分の速度ベクトル（例えば、ＵＡＶ１０のふらつきを示す速度ベクトル）は、ほぼ同等の大きさを持ち、かつ、向きがばらついているパターンを有する場合がある。上述の実施形態では、大きさの度数分布において度数が比較的大きい範囲に属し、かつ、向きの度数分布において度数が比較的小さい範囲に属するベクトル（即ち上述のパターンに合致するベクトル）を、ノイズ成分として適切に除去することができる。
【0071】
なお、上述の第１閾値又は第２閾値として、各ヒストグラムの縦軸の平均値又は中央値を用いてもよい。
【0072】
なお、複数の速度ベクトルの大きさの度数分布（図１３参照）及び向きの度数分布（図１４参照）の両方において、度数が閾値よりも大きい速度ベクトル（例えば、図１３においてＰ１で示すピーク、及び、図１４においてＰ２で示すピークの両方に含まれるベクトル）は、監視対象のガスの速度ベクトルである可能性が高い。
【0073】
（リフレーム処理について）
ステップＳ２００のリフレーム処理について説明する。ここで、図１５～図１７は、複数の画像データのリフレーム処理を説明するための図である。図１５～図１７に示す画像Ｉ１１～Ｉ１３は、ＵＡＶ１０に搭載された撮像装置２０で時系列に撮影された画像である。リフレーム処理では、時系列に撮影された各画像の画像データについて、基準点Ｐ０がフレームＦ（画像の表示領域）内の所定位置（例えばフレームＦの中心位置）に位置するように画像の縦横比（ａ：ｂ）を調節して（即ち、画像のうちフレームＦ外の領域をトリミングして）、フレームＦに対する画像の位置を調整する。図示する例では、ＣＯ２ガスボンベＢ１，Ｂ２の位置を基準点Ｐ０として、各画像Ｉ１１～Ｉ１３の各々について基準点Ｐ０がフレームの中心に位置するように、画像の縦横比が、それぞれａ１：ｂ１、ａ２：ｂ２、ａ３：ｂ３となるように調節して画像の位置を調整している。
ステップＳ２００のリフレーム処理について説明する。ここで、図１５～図１７は、複数の画像データのリフレーム処理を説明するための図である。図１５～図１７に示す画像Ｉ１１～Ｉ１３は、ＵＡＶ１０に搭載された撮像装置２０で時系列に撮影された画像である。リフレーム処理では、時系列に撮影された各画像の画像データについて、基準点Ｐ０がフレームＦ（画像の表示領域）内の所定位置（例えばフレームＦの中心位置）に位置するように画像の縦横比（ａ：ｂ）を調節して（即ち、画像のうちフレームＦ外の領域をトリミングして）、フレームＦに対する画像の位置を調整する。図示する例では、ＣＯ２ガスボンベＢ１，Ｂ２の位置を基準点Ｐ０として、各画像Ｉ１１～Ｉ１３の各々について基準点Ｐ０がフレームの中心に位置するように、画像の縦横比が、それぞれａ１：ｂ１、ａ２：ｂ２、ａ３：ｂ３となるように調節して画像の位置を調整している。
【0074】
撮像装置２０での撮影中にＵＡＶ１０は移動する（例えばふらつく）ため、ＵＡＶ１０に搭載された撮像装置２０で連続的に撮影される複数の画像データは、撮影位置がずれる場合がある。この点、上述の実施形態によれば、複数の画像データについてリフレーム処理をするため、複数の画像データに含まれる撮影対象物のフレームＦ中の位置を固定（いわゆるブレ補正）することができる。このように、リフレーム処理をすることで、予めＵＡＶ１０のふらつきによるノイズを除去してから、ステップＳ３００でオプティカルフロー推定処理することで、ステップＳ８００におけるガス速度ベクトルの抽出精度をより向上することができる。
【0075】
上記各実施形態に記載の内容は、例えば以下のように把握される。
【0076】
（１）本発明の少なくとも一実施形態に係るガス監視方法は、
ＵＡＶ（１０）に搭載された撮像装置（２０）で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップ（Ｓ４００）と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップ（Ｓ７００）と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップ（Ｓ８００）と、
を備える。
ＵＡＶ（１０）に搭載された撮像装置（２０）で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得するステップ（Ｓ４００）と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するステップ（Ｓ７００）と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出ステップ（Ｓ８００）と、
を備える。
【0077】
上記（１）の方法では、ＵＡＶに搭載された撮像装置で撮像された時系列の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルは、検出対象のガスの速度ベクトル（ガス速度ベクトル）に加え、他の速度ベクトル（検出対象のガス以外の物体の移動を示す速度ベクトルやノイズ等）も含まれる。また、ＵＡＶと撮影対象の距離が一定でなく、さらに速度算出の基準となるスケールが画像内に存在するとは限らないため、該複数の速度ベクトルの各々の大きさ（速度の絶対値）はオプティカルフロー推定では演算できない。一方、上記（１）の構成では、ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき、ＵＡＶのダウンウォッシュ速度ｕを取得することができる。そこで、オプティカルフロー推定で得られた未知の大きさを有する複数の速度ベクトルのうち、ＵＡＶのダウンウォッシュの速度を示すベクトルの大きさを基準にして、該複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づき決定される範囲内の大きさを有する速度ベクトルを監視対象のガスの移動速度及び方向を示すガス速度ベクトル（ガス領域）として特定することができる。このように抽出されたガス速度ベクトルを用いて、ガスを適切に監視することができる。
【0078】
（２）幾つかの実施形態では、上記（１）の方法において、
前記ＵＡＶの飛行状態を示す前記少なくとも１つの変数は、前記ＵＡＶの飛行高さＨ及び前記ＵＡＶの飛行速度ｖを含む。
前記ＵＡＶの飛行状態を示す前記少なくとも１つの変数は、前記ＵＡＶの飛行高さＨ及び前記ＵＡＶの飛行速度ｖを含む。
【0079】
ＵＡＶのダウンウォッシュ速度ｕは、ＵＡＶの飛行高さＨや飛行速度ｖ及び相関関係を有する。上記（２）の方法によれば、ＵＡＶの飛行高さＨ及びＵＶＡの飛行速度ｖに基づいて、これらとの相関関係から、ＵＡＶのダウンウォッシュ速度ｕを適切に取得することができる。
【0080】
（３）幾つかの実施形態では、上記（１）又は（２）の方法において、
前記ガス監視方法は、
地表における風速Ｗｇを取得するステップ（Ｓ７００）と、
前記ダウンウォッシュ速度ｕ及び前記風速Ｗｇに基づいて、前記風速Ｗｇを考慮したダウンウォッシュ速度である補正ダウンウォッシュ速度ｕｇを算出するステップ（Ｓ７００）と、を備え、
前記抽出ステップでは、前記複数の速度ベクトルから、前記補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する。
前記ガス監視方法は、
地表における風速Ｗｇを取得するステップ（Ｓ７００）と、
前記ダウンウォッシュ速度ｕ及び前記風速Ｗｇに基づいて、前記風速Ｗｇを考慮したダウンウォッシュ速度である補正ダウンウォッシュ速度ｕｇを算出するステップ（Ｓ７００）と、を備え、
前記抽出ステップでは、前記複数の速度ベクトルから、前記補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する。
【0081】
地表においてダウンウォッシュとは別の風が吹いてる場合、地表におけるダウンウォッシュの流れは該風の影響を受け得る。上記（３）の方法によれば、ＵＡＶの飛行状態を示す変数に基づき得られたダウンウォッシュ速度u及び地表における風速Ｗｇに基づいて風速Ｗｇを考慮した補正ダウンウォッシュ速度ｕｇを算出し、オプティカルフロー推定で得られた複数のベクトルから、該補正ダウンウォッシュ速度ｕｇに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する。したがって、監視対象のガスの移動速度及び方向を示すガス速度ベクトルをより適切に抽出することができる。
【0082】
（４）幾つかの実施形態では、上記（３）の方法において、
前記ダウンウォッシュ速度ｕのベクトルと前記地表における風速Ｗｇのベクトルとを合成することで前記補正ダウンウォッシュ速度ｕｇのベクトルを算出する。
前記ダウンウォッシュ速度ｕのベクトルと前記地表における風速Ｗｇのベクトルとを合成することで前記補正ダウンウォッシュ速度ｕｇのベクトルを算出する。
【0083】
上記（４）の方法によれば、ダウンウォッシュ速度ｕと風速Ｗｇの合成ベクトルとして、補正ダウンウォッシュ速度ｕｇを適切に算出することができる。
【0084】
（５）幾つかの実施形態では、上記（３）又は（４）の方法において、
前記ガス監視方法は、
地表からの高さと風速との相関関係を取得するステップ（Ｓ７００）と、
前記ＵＡＶの飛行高さにおける風速Ｗ１を取得するステップ（Ｓ７００）と、を備え、
前記地表における風速Ｗｇを取得するステップでは、前記風速Ｗ１及び前記相関関係に基づいて前記地表における風速Ｗｇを取得する。
前記ガス監視方法は、
地表からの高さと風速との相関関係を取得するステップ（Ｓ７００）と、
前記ＵＡＶの飛行高さにおける風速Ｗ１を取得するステップ（Ｓ７００）と、を備え、
前記地表における風速Ｗｇを取得するステップでは、前記風速Ｗ１及び前記相関関係に基づいて前記地表における風速Ｗｇを取得する。
【0085】
上記（５）の方法によれば、地表からの高さと風速との相関関係に基づいて、ＵＡＶの飛行高さにおける風速Ｗ１から、地表における風速Ｗｇを取得する。このように取得した地表での風速Ｗｇに基づいて、上述の補正ダウンウォッシュ速度ｕｇを適切に算出することができる。
【0086】
（６）幾つかの実施形態では、上記（５）の方法において、
前記風速Ｗ１を取得するステップでは、前記ＵＡＶのモータ電流値に基づいて前記ＵＡＶの飛行高さにおける前記風速Ｗ１を算出する。
前記風速Ｗ１を取得するステップでは、前記ＵＡＶのモータ電流値に基づいて前記ＵＡＶの飛行高さにおける前記風速Ｗ１を算出する。
【0087】
ＵＡＶの飛行高さにおける風速Ｗ１と、ＵＡＶのプロペラを駆動するモータの電流値とは、所定の相関関係を有する。上記（６）の方法によれば、例えば風速計等を用いなくても、ＵＡＶのモータ電流値に基づいて、ＵＡＶの飛行高さにおける風速Ｗ１を適切に取得することができる。
【0088】
（７）幾つかの実施形態では、上記（１）乃至（６）の何れかの方法において、
前記撮像装置は、監視対象のガスが吸収する波長の赤外線を選択的に透過させるフィルタを備えた赤外線カメラを含む。
前記撮像装置は、監視対象のガスが吸収する波長の赤外線を選択的に透過させるフィルタを備えた赤外線カメラを含む。
【0089】
上記（７）の方法によれば、撮像装置として、監視対象のガスが吸収する波長の赤外線を選択的に透過させるフィルタを備えた赤外線カメラを用いるので、監視対象のガスを撮影することができる。そして、これにより得られる時系列の画像データをオプティカルフロー処理することで、監視対象ガスの移動を示す速度ベクトルを取得することができる。
【0090】
（８）幾つかの実施形態では、上記（１）乃至（７）の何れかの方法において、
前記ガス監視方法は、
前記複数の画像データをオプティカルフロー推定処理して得られる前記複数の速度ベクトルから、前記複数の速度ベクトルの大きさの度数分布又は前記複数の速度ベクトルの向きの度数分布に基づいてノイズ成分を除去するステップ（Ｓ５００）を備え、
前記抽出ステップでは、前記ノイズ成分が除去された前記複数の速度ベクトルから前記ガス速度ベクトルを抽出する。
前記ガス監視方法は、
前記複数の画像データをオプティカルフロー推定処理して得られる前記複数の速度ベクトルから、前記複数の速度ベクトルの大きさの度数分布又は前記複数の速度ベクトルの向きの度数分布に基づいてノイズ成分を除去するステップ（Ｓ５００）を備え、
前記抽出ステップでは、前記ノイズ成分が除去された前記複数の速度ベクトルから前記ガス速度ベクトルを抽出する。
【0091】
ノイズ成分の速度ベクトルは、大きさ又は向きに関して所定のパターンを有する傾向がある。上記（８）の方法によれば、オプティカルフロー推定により得られる複数の速度ベクトルの大きさ又は向きの度数分布に基づいて、該複数のベクトルからノイズ成分を除去することができる。また、ノイズ成分が除去された複数のベクトルからダウンウォッシュ速度ｕに基づきガス速度ベクトルを抽出するので、ガス速度ベクトルの抽出精度を向上することができる。
【0092】
（９）幾つかの実施形態では、上記（８）の方法において、
前記複数の速度ベクトルの大きさの度数分布において、度数が第１閾値以上の大きさの範囲に属し、かつ、前記複数の速度ベクトルの向きの度数分布において、度数が第２閾値以下の向きの範囲に属する速度ベクトルを、前記複数の速度ベクトルから前記ノイズ成分として除去する。
前記複数の速度ベクトルの大きさの度数分布において、度数が第１閾値以上の大きさの範囲に属し、かつ、前記複数の速度ベクトルの向きの度数分布において、度数が第２閾値以下の向きの範囲に属する速度ベクトルを、前記複数の速度ベクトルから前記ノイズ成分として除去する。
【0093】
ノイズ成分の速度ベクトルは、ほぼ同等の大きさを持ち、かつ、向きがばらついているパターンを有する場合がある。上記（９）の方法によれば、大きさの度数分布において度数が比較的大きい範囲に属し、かつ、向きの度数分布において度数が比較的小さい範囲に属するベクトル（即ち上述のパターンに合致するベクトル）を、ノイズ成分として適切に除去することができる。
【0094】
（１０）幾つかの実施形態では、上記（１）乃至（９）の何れかの方法において、
前記ガス監視方法は、
前記ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データを取得するステップ（Ｓ１００）と、
前記複数の画像データについてリフレーム処理を施すステップ（Ｓ２００）と、
前記リフレーム処理が施された複数の画像データをオプティカルフロー推定処理して前記複数の速度ベクトルを取得するステップ（Ｓ３００，Ｓ４００）と、
を備える。
前記ガス監視方法は、
前記ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データを取得するステップ（Ｓ１００）と、
前記複数の画像データについてリフレーム処理を施すステップ（Ｓ２００）と、
前記リフレーム処理が施された複数の画像データをオプティカルフロー推定処理して前記複数の速度ベクトルを取得するステップ（Ｓ３００，Ｓ４００）と、
を備える。
【0095】
撮像装置での撮影中にＵＡＶは移動する（例えばふらつく）ため、ＵＡＶに搭載された撮像装置で連続的に撮影される複数の画像データは、撮影位置がずれる場合がある。この点、上記（１０）の方法によれば、複数の画像データについてリフレーム処理をするため、複数の画像データに含まれる撮影対象物のフレーム中の位置を固定（いわゆるブレ補正）することができる。このように、リフレーム処理をすることで、予めＵＡＶのふらつきによるノイズを除去してから、オプティカルフロー推定処理することで、ガス速度ベクトルの抽出精度をより向上することができる。
【0096】
（１１）本発明の少なくとも一実施形態に係るガス監視装置（５０）は、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する速度ベクトル取得部（５８）と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するダウンウォッシュ速度取得部（６４）と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出部（６６）と、
を備える。
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する速度ベクトル取得部（５８）と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得するダウンウォッシュ速度取得部（６４）と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する抽出部（６６）と、
を備える。
【0097】
上記（１１）の構成では、ＵＡＶに搭載された撮像装置で撮像された時系列の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルは、検出対象のガスの速度ベクトル（ガス速度ベクトル）に加え、他の速度ベクトル（検出対象のガス以外の物体の移動を示す速度ベクトルやノイズ等）も含まれる。また、ＵＡＶと撮影対象の距離が一定でなく、さらに速度算出の基準となるスケールが画像内に存在するとは限らないため、該複数の速度ベクトルの各々の大きさ（速度の絶対値）はオプティカルフロー推定では演算できない。一方、上記（１１）の構成では、ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき、ＵＡＶのダウンウォッシュ速度ｕを取得することができる。そこで、オプティカルフロー推定で得られた未知の大きさを有する複数の速度ベクトルのうち、ＵＡＶのダウンウォッシュの速度を示すベクトルの大きさを基準にして、該複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づき決定される範囲内の大きさを有する速度ベクトルを監視対象のガスの移動速度及び方向を示すガス速度ベクトル（ガス領域）として特定することができる。このように抽出されたガス速度ベクトルを用いて、ガスを適切に監視することができる。
【0098】
（１２）本発明の少なくとも一実施形態に係るガス監視システム（１００）は、
撮像装置（２０）を搭載したＵＡＶ（１０）と、
前記撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して複数の速度ベクトルを算出するように構成されたオプティカルフロー推定処理部（５６）と、
前記複数の速度ベクトルから前記ガス速度ベクトルを抽出するように構成された上記（１１）に記載のガス監視装置（５０）と、
を備える。
撮像装置（２０）を搭載したＵＡＶ（１０）と、
前記撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して複数の速度ベクトルを算出するように構成されたオプティカルフロー推定処理部（５６）と、
前記複数の速度ベクトルから前記ガス速度ベクトルを抽出するように構成された上記（１１）に記載のガス監視装置（５０）と、
を備える。
【0099】
上記（１２）の構成では、ＵＡＶに搭載された撮像装置で撮像された時系列の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルは、検出対象のガスの速度ベクトル（ガス速度ベクトル）に加え、他の速度ベクトル（検出対象のガス以外の物体の移動を示す速度ベクトルやノイズ等）も含まれる。また、ＵＡＶと撮影対象の距離が一定でなく、さらに速度算出の基準となるスケールが画像内に存在するとは限らないため、該複数の速度ベクトルの各々の大きさ（速度の絶対値）はオプティカルフロー推定では演算できない。一方、上記（１２）の構成では、ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき、ＵＡＶのダウンウォッシュ速度ｕを取得することができる。そこで、オプティカルフロー推定で得られた未知の大きさを有する複数の速度ベクトルのうち、ＵＡＶのダウンウォッシュの速度を示すベクトルの大きさを基準にして、該複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づき決定される範囲内の大きさを有する速度ベクトルを監視対象のガスの移動速度及び方向を示すガス速度ベクトル（ガス領域）として特定することができる。このように抽出されたガス速度ベクトルを用いて、ガスを適切に監視することができる。
【0100】
（１３）本発明の少なくとも一実施形態に係るガス監視プログラムは、
コンピュータに、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する手順と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得する手順と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する手順と、
実行させる。
コンピュータに、
ＵＡＶに搭載された撮像装置で複数の時刻にそれぞれ撮像して得られる複数の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルを取得する手順と、
前記ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき前記ＵＡＶのダウンウォッシュ速度ｕを取得する手順と、
前記複数の速度ベクトルから、前記ダウンウォッシュ速度ｕに基づいて決定される範囲内の大きさを有するガス速度ベクトルを抽出する手順と、
実行させる。
【0101】
上記（１３）の構成では、ＵＡＶに搭載された撮像装置で撮像された時系列の画像データをオプティカルフロー推定処理して得られる複数の速度ベクトルは、検出対象のガスの速度ベクトル（ガス速度ベクトル）に加え、他の速度ベクトル（検出対象のガス以外の物体の移動を示す速度ベクトルやノイズ等）も含まれる。また、ＵＡＶと撮影対象の距離が一定でなく、さらに速度算出の基準となるスケールが画像内に存在するとは限らないため、該複数の速度ベクトルの各々の大きさ（速度の絶対値）はオプティカルフロー推定では演算できない。一方、上記（１３）の構成では、ＵＡＶの飛行状態を示す少なくとも１つの変数に基づき、ＵＡＶのダウンウォッシュ速度ｕを取得することができる。そこで、オプティカルフロー推定で得られた未知の大きさを有する複数の速度ベクトルのうち、ＵＡＶのダウンウォッシュの速度を示すベクトルの大きさを基準にして、該複数の速度ベクトルから、ダウンウォッシュ速度ｕに基づき決定される範囲内の大きさを有する速度ベクトルを監視対象のガスの移動速度及び方向を示すガス速度ベクトル（ガス領域）として特定することができる。このように抽出されたガス速度ベクトルを用いて、ガスを適切に監視することができる。
【0102】
以上、本発明の実施形態について説明したが、本発明は上述した実施形態に限定されることはなく、上述した実施形態に変形を加えた形態や、これらの形態を適宜組み合わせた形態も含む。
【0103】
本明細書において、「ある方向に」、「ある方向に沿って」、「平行」、「直交」、「中心」、「同心」或いは「同軸」等の相対的或いは絶対的な配置を表す表現は、厳密にそのような配置を表すのみならず、公差、若しくは、同じ機能が得られる程度の角度や距離をもって相対的に変位している状態も表すものとする。
例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
また、本明細書において、四角形状や円筒形状等の形状を表す表現は、幾何学的に厳密な意味での四角形状や円筒形状等の形状を表すのみならず、同じ効果が得られる範囲で、凹凸部や面取り部等を含む形状も表すものとする。
また、本明細書において、一の構成要素を「備える」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。
例えば、「同一」、「等しい」及び「均質」等の物事が等しい状態であることを表す表現は、厳密に等しい状態を表すのみならず、公差、若しくは、同じ機能が得られる程度の差が存在している状態も表すものとする。
また、本明細書において、四角形状や円筒形状等の形状を表す表現は、幾何学的に厳密な意味での四角形状や円筒形状等の形状を表すのみならず、同じ効果が得られる範囲で、凹凸部や面取り部等を含む形状も表すものとする。
また、本明細書において、一の構成要素を「備える」、「含む」、又は、「有する」という表現は、他の構成要素の存在を除外する排他的な表現ではない。
【符号の説明】
【0104】
１０ ＵＡＶ
１２ ＵＡＶ本体
１４ プロペラ
２０ 撮像装置
５０ ガス監視装置
５２ 画像取得部
５４ リフレーム処理部
５６ オプティカルフロー処理部
５８ 速度ベクトル取得部
６０ ノイズ成分除去部
６２ 飛行状態取得部
６４ ダウンウォッシュ速度取得部
６６ 抽出部
６８ 表示画像生成部
６９ 記憶部
７０ 表示部
１００ ガス監視システム
ＡＧ ガス領域
Ｂ１ ＣＯ２ガスボンベ
Ｂ２ ＣＯ２ガスボンベ
Ｄ ダウンウォッシュ
Ｆ フレーム
Ｇ 地表
Ｐ０ 基準点
ＶＧ ガス速度ベクトル
１２ ＵＡＶ本体
１４ プロペラ
２０ 撮像装置
５０ ガス監視装置
５２ 画像取得部
５４ リフレーム処理部
５６ オプティカルフロー処理部
５８ 速度ベクトル取得部
６０ ノイズ成分除去部
６２ 飛行状態取得部
６４ ダウンウォッシュ速度取得部
６６ 抽出部
６８ 表示画像生成部
６９ 記憶部
７０ 表示部
１００ ガス監視システム
ＡＧ ガス領域
Ｂ１ ＣＯ２ガスボンベ
Ｂ２ ＣＯ２ガスボンベ
Ｄ ダウンウォッシュ
Ｆ フレーム
Ｇ 地表
Ｐ０ 基準点
ＶＧ ガス速度ベクトル
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】