(57)【要約】  漏洩源特定支援装置は、漏洩流体領域（ガス領域）を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理部を備え、前記処理部は、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する。
2020110410.tif 000002
【特許請求の範囲】
【請求項１】
漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理部を備え、
前記処理部は、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する、漏洩源特定支援装置。
【請求項２】
前記処理部は、前記第１画素の動きの方向と逆方向を示す軸を基準にして、所定の角度範囲内にある二以上の前記画素を前記第２画素の候補に決定し、前記候補の中から前記第２画素を決定する、請求項１に記載の漏洩源特定支援装置。
【請求項３】
前記所定の角度範囲は、前記軸を基準にして±９０度の範囲である、請求項２に記載の漏洩源特定支援装置。
【請求項４】
前記処理部は、前記軸に対する角度が最も小さい位置にある前記候補を、前記第２画素と決定する、請求項２または３に記載の漏洩源特定支援装置。
【請求項５】
前記処理部は、前記候補のうち、前記漏洩流体領域外に位置する前記候補を前記候補から除外する、請求項２〜４のいずれか一項に記載の漏洩源特定支援装置。
【請求項６】
前記処理部は、前記候補の全てが前記漏洩流体領域外に位置するとき、前記候補の決定に用いられた前記第１画素を前記最終移動元画素と決定する、請求項５に記載の漏洩源特定支援装置。
【請求項７】
前記処理部は、前記第１画素の動きの量が所定のしきい値より小さいとき、当該第１画素を前記最終移動元画素と決定する、請求項１〜６のいずれか一項に記載の漏洩源特定支援装置。
【請求項８】
前記画像を構成する複数の前記画素とは、前記漏洩流体領域を構成する全画素である、請求項１〜７のいずれか一項に記載の漏洩源特定支援装置。
【請求項９】
前記画像を構成する複数の前記画素のそれぞれの前記最終移動元画素において、最頻の前記最終移動元画素を漏洩源の位置と推定する推定部をさらに備える、請求項１〜８のいずれか一項に記載の漏洩源特定支援装置。
【請求項１０】
前記処理部は、前記画像を構成する複数の前記画素のそれぞれについて、前記最終移動元画素を決定することについて、異なる時刻の前記画素の動きをそれぞれ示す複数の前記画像のそれぞれに対して実行し、
前記推定部は、複数の前記画像のそれぞれについて得られた前記最終移動元画素において、最頻の前記最終移動元画素を前記漏洩源の位置と推定する、請求項９に記載の漏洩源特定支援装置。
【請求項１１】
漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理ステップを備え、
前記処理ステップは、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する、漏洩源特定支援方法。
【請求項１２】
漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理ステップを、コンピュータに実行させ、
前記処理ステップは、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する、漏洩源特定支援プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、漏洩流体の漏洩源の特定を支援する技術に関する。
【背景技術】
【０００２】
流体（例えば、ガス）の漏洩が発生したとき、漏洩源を特定し、漏洩箇所の修理が必要となる。漏洩源を検出できる技術として、例えば、特許文献１に開示された流体漏洩箇所検出装置がある。この装置は、対象領域を撮像した時系列な複数の画像を取得する画像取得部と、前記画像取得部で取得した前記時系列な複数の画像に基づいて、流体を収容する流体収容部の外部に在る漏洩流体の移動軌跡を求める移動軌跡演算部と、前記移動軌跡演算部で求めた前記移動軌跡に基づいて、前記漏洩流体が前記流体収容部から漏れ出る漏洩箇所を求める漏洩箇所演算部と、を備える。【０００３】
本発明者は、特許文献１の技術とは異なる原理で、漏洩流体領域（例えば、ガス領域）を含む画像を基にして、比較的簡単に漏洩源の特定を支援できる技術を検討した。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】国際公開第２０１６／２０８３１７号
【発明の概要】
【０００５】
本発明は、漏洩流体領域を含む画像を基にして、比較的簡単に漏洩源の特定を支援できる漏洩源特定支援装置、漏洩源特定支援方法、及び、漏洩源特定支援プログラムを提供することを目的とする。
【０００６】
上述した目的を実現するために、本発明の一側面を反映した漏洩源特定支援装置は、処理部を有する。前記処理部は、漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する。そして、前記処理部は、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する。
【０００７】
発明の１又は複数の実施形態により与えられる利点及び特徴は以下に与えられる詳細な説明及び添付図面から十分に理解される。これら詳細な説明及び添付図面は、例としてのみ与えられるものであり本発明の限定の定義として意図されるものではない。
【図面の簡単な説明】
【０００８】
【図１Ａ】実施形態に係る漏洩源特定支援システムの構成を示すブロック図である。
【図１Ｂ】図１Ａに示す漏洩源特定支援装置のハードウェア構成を示すブロック図である。
【図２】実施形態に係る漏洩源特定支援装置の動作を説明するフローチャートである。
【図３】赤外画像の動画データを構成するフレーム群の中から選択された４つのフレームを時系列に並べた画像図である。
【図４】ガス領域を見やすくする画像処理がされた画像を時系列で示す画像図である。
【図５】最大画素値で構成される画像の一例を示す画像図である。
【図６】図５に示す画像に対してノイズ除去および二値化処理をした画像の一例を示す画像図である。
【図７】図４に示す画像のオプティカルフローを示す画像において、ガス領域の一部を拡大した模式図である。
【図８】移動元画素の追跡を説明するフローチャートである。
【図９】移動元画素の追跡を説明する第１説明図である。
【図１０】移動元画素の追跡を説明する第２説明図である。
【図１１】移動元画素の追跡を説明する第３説明図である。
【図１２】第１画素（画素Ｐ−１）からの移動ベクトルを説明する説明図である。
【図１３】第１画素（画素Ｐ−３）からの移動ベクトルを説明する説明図である。
【図１４】図７に示す画像おいて、画素Ｐ−１について、移動元画素を追跡した結果を説明する説明図である。
【図１５】図７に示す画像において、ガス領域内の各画素と最終移動元画素との関係を説明する説明図である。
【発明を実施するための形態】
【０００９】
以下、図面を参照して、本発明の１又は複数の実施形態が説明される。しかし、発明の範囲は、開示された実施形態に限定されない。
【００１０】
各図において、同一符号を付した構成は、同一の構成であることを示し、その構成について、既に説明している内容については、その説明を省略する。本明細書において、総称する場合にはハイフンを省略した参照符号で示し（例えば、画像Ｉｍ２）、個別の構成を指す場合にはハイフンを付した参照符号で示す（例えば、画像Ｉｍ２−１）。
【００１１】
図１Ａは、実施形態に係る漏洩源特定支援システム１の構成を示すブロック図である。漏洩源特定支援システム１は、赤外線カメラ２と漏洩源特定支援装置３とを備える。実施形態では、漏洩する流体としてガスを例にして説明するが、漏洩する流体は液体でもよい。
【００１２】
赤外線カメラ２は、ガス漏れの監視対象（例えば、ガス輸送管どうしが接続されている箇所）を含む被写体について、赤外画像の動画を撮影し、動画を示す動画データＭＤを生成する。時系列に撮像された複数の赤外画像であればよく、動画に限定されない。赤外線カメラ２は、光学系４、フィルター５、二次元イメージセンサー６及び信号処理部７を備える。
【００１３】
光学系４は、被写体の赤外画像を二次元イメージセンサー６上で結像させる。フィルター５は、光学系４と二次元イメージセンサー６との間に配置され、光学系４を通過した光のうち、特定波長の赤外線のみを通過させる。赤外の波長帯のうち、フィルター５を通過させる波長帯は、検知するガスの種類に依存する。例えばメタンの場合、３．２〜３．４μｍの波長帯を通過させるフィルター５が用いられる。二次元イメージセンサー６は、例えば、冷却型インジウムアンチモン（ＩｎＳｂ）イメージセンサーであり、フィルター５を通過した赤外線を受光する。信号処理部７は、二次元イメージセンサー６から出力されたアナログ信号を、デジタル信号に変換し、公知の画像処理をする。このデジタル信号が、動画データＭＤとなる。
【００１４】
漏洩源特定支援装置３は、パーソナルコンピュータ、スマートフォン、タブレット端末等であり、機能ブロックとして、画像データ入力部８、演算処理部９、表示制御部１０、ディスプレイ１１及び入力部１２を備える。
【００１５】
画像データ入力部８は、赤外線カメラ２の通信部（不図示）と通信する通信インターフェイスである。画像データ入力部８には、赤外線カメラ２の通信部から送られてきた動画データＭＤが入力される。画像データ入力部８は、動画データＭＤを演算処理部９へ送る。
【００１６】
演算処理部９は、動画データＭＤを用いて、ガスの漏洩源を特定する支援に必要な処理をする。演算処理部９は、記憶部９１と、領域抽出部９２と、ベクトル算出部９３と、画素決定処理部９４と、漏洩源推定部９５と、を備える。
【００１７】
記憶部９１は、動画データＭＤ、ガスの漏洩源を特定する支援に必要な各種のプログラム等を記憶する。領域抽出部９２は、赤外画像に対して画像処理をして、赤外画像からガス領域を抽出する。ガス領域は、漏洩流体領域の一例である。ベクトル算出部９３は、後で説明する画像Ｉｍ２（図４）を構成する各画素の移動ベクトルを算出する。
【００１８】
画素決定処理部９４（処理部）は、ガス領域を含む画像を構成する画素の１つである第１画素の動きを基にして、その画像を構成する画素の１つであり、第１画素の移動元画素となる第２画素を決定する処理をし、第２画素を新たな第１画素として、その処理を繰り返すことにより、最終移動元画素を決定する。画素決定処理部９４は、その画像を構成する複数の画素のそれぞれについて、その処理を繰り返すことにより、複数の画素のそれぞれについて、最終移動元画素を決定する。なお、第２画素は、好ましくは、第１画素の近傍に位置しており、例えば、第１画素と第２の画素との間に画素が存在しなくてもよいし、第１画素と第２画素との間に数画素が存在してもよい。
【００１９】
上記複数の画素は、例えば、ガス領域を含む画像の全画素でもよいし、その画像の全画素の中から選択された画素でもよいし、ガス領域の全画素でもよいし、ガス領域の全画素の中から選択された画素でもよい。最終移動元画素の数は、複数の画素の数以下である。例えば、複数の画素が３つの画素とする。３つの画素の最終移動元画素がそれぞれ異なることもあるし、同じこともあるし、３つのうち２つの画素の最終移動元画素が同じこともある。
【００２０】
ユーザは、複数の画素のそれぞれの最終移動元画素を参考にして、漏洩源を特定してもよいし、次に説明するように、漏洩源推定部９５が、複数の画素のそれぞれの最終移動元画素を基にして、漏洩源を推定してもよい。ユーザは、この推定を参考にして、漏洩源を特定することができる。
【００２１】
漏洩源推定部９５（推定部）は、画像を構成する複数の画素のそれぞれの最終移動元画素において、最頻の最終移動元画素をガスの漏洩源の位置と推定する。
【００２２】
表示制御部１０は、動画データＭＤで示される動画、演算処理部９で演算処理された画像（例えば、ガス領域の抽出処理がされた画像）をディスプレイ１１に表示させる。
【００２３】
入力部１２は、ガスの漏洩源を特定する支援に関連する各種入力がされる。実施形態に係る漏洩源特定支援装置３は、表示制御部１０、ディスプレイ１１及び入力部１２を備えるが、これらを備えない漏洩源特定支援装置３でもよい。
【００２４】
図１Ｂは、図１Ａに示す漏洩源特定支援装置３のハードウェア構成を示すブロック図である。漏洩源特定支援装置３は、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）３ａ、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）３ｂ、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）３ｃ、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）３ｄ、液晶ディスプレイ３ｅ、通信インターフェイス３ｆ、キーボード等３ｇ、及び、これらを接続するバス３ｈを備える。液晶ディスプレイ３ｅは、ディスプレイ１１を実現するハードウェアである。液晶ディスプレイ３ｅの替わりに、有機ＥＬディスプレイ（Ｏｒｇａｎｉｃ  Ｌｉｇｈｔ  Ｅｍｉｔｔｉｎｇ  Ｄｉｏｄｅ  ｄｉｓｐｌａｙ）、プラズマディスプレイ等でもよい。通信インターフェイス３ｆは、画像データ入力部８を実現するハードウェアである。キーボード等３ｇは、入力部１２を実現するハードウェアである。キーボードの替わりに、タッチパネルでもよい。
【００２５】
ＨＤＤ３ｄには、演算処理部９及び表示制御部１０について、これらの機能ブロックをそれぞれ実現するためのプログラム、及び、各種データ（例えば、動画データＭＤ）が格納されている。演算処理部９を実現するプログラムは、動画データＭＤを取得し、動画データＭＤを用いて、ガス漏洩源の特定を支援する演算処理プログラムである。表示制御部１０を実現するプログラムは、例えば、動画データＭＤで示される動画をディスプレイ１１に表示させたり、演算処理部９で演算処理された画像をディスプレイ１１に表示させたりする表示制御プログラムである。これらのプログラムは、ＨＤＤ３ｄに予め記憶されているが、これに限定されない。例えば、これらのプログラムを記録している記録媒体（例えば、磁気ディスク、光学ディスクのような外部記録媒体）が用意されており、この記録媒体に記憶されているプログラムがＨＤＤ３ｄに記憶されてもよい。また、これらのプログラムは、漏洩源特定支援装置３とネットワーク接続されたサーバに格納されており、ネットワークを介して、これらのプログラムがＨＤＤ３ｄに送られ、ＨＤＤ３ｄに記憶されてもよい。これらのプログラムは、ＨＤＤ３ｄの替わりにＲＯＭ３ｃに記憶してもよい。漏洩源特定支援装置３は、ＨＤＤ３ｄの替わりに、フラッシュメモリを備え、これらのプログラムはフラッシュメモリに記憶してもよい。
【００２６】
ＣＰＵ３ａは、これらのプログラムを、ＨＤＤ３ｄから読み出してＲＡＭ３ｂに展開させ、展開されたプログラムを実行することによって、演算処理部９及び表示制御部１０が実現される。但し、演算処理部９の機能及び表示制御部１０の機能について、各機能の一部又は全部は、ＣＰＵ３ａによる処理に替えて、又は、これと共に、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）による処理によって実現されてもよい。又、同様に、各機能の一部又は全部は、ソフトウェアによる処理に替えて、又は、これと共に、専用のハードウェア回路による処理によって実現されてもよい。
【００２７】
なお、演算処理部９は、図１Ａに示す複数の要素によって構成される。従って、ＨＤＤ３ｄには、これらの要素を実現するためのプログラムが格納されている。すなわち、ＨＤＤ３ｄには、領域抽出部９２、ベクトル算出部９３、画素決定処理部９４、漏洩源推定部９５のそれぞれを実現するためのプログラムが格納されている。これらのプログラムは、領域抽出プログラム、ベクトル算出プログラム、画素決定処理プログラム、漏洩源推定プログラムと表現される。
【００２８】
これらのプログラムは、要素の定義を用いて表現される。漏洩源推定部９５および漏洩源推定プログラムを例にして説明する。漏洩源推定部９５は、画像を構成する複数の画素のそれぞれの最終移動元画素において、最頻の最終移動元画素を漏洩源の位置と推定する。漏洩源推定プログラムは、画像を構成する複数の画素のそれぞれの最終移動元画素において、最頻の最終移動元画素を漏洩源の位置と推定するプログラムである。
【００２９】
ＣＰＵ３ａによって実行されるこれらのプログラム（領域抽出プログラム、ベクトル算出プログラム等）のフローチャートが、次に説明する図２である。
【００３０】
赤外線カメラ２を用いて屋外の試験場所を撮影して得られた赤外画像の動画データＭＤを例にして、漏洩源特定支援装置３の動作を説明する。図２は、この動作を説明するフローチャートである。図３は、赤外画像の動画データＭＤを構成するフレーム群の中から選択された４つのフレームを時系列に並べた画像図である。試験場所には、ガスを噴出させることができる地点ＳＰがある。地点ＳＰではガスが噴出している。画像Ｉｍ１−１は、時刻Ｔ１に撮影された赤外画像（フレーム）である。画像Ｉｍ１−２は、時刻Ｔ１から０．２秒後の時刻Ｔ２に撮影された赤外画像（フレーム）である。画像Ｉｍ１−３は、時刻Ｔ１から０．４秒後の時刻Ｔ３に撮影された赤外画像（フレーム）である。画像Ｉｍ１−４は、時刻Ｔ１から０．６秒後の時刻Ｔ４に撮影された赤外画像（フレーム）である。
【００３１】
画像Ｉｍ１（赤外画像）からは、ガス領域（ここでは、ガス領域とは、地点ＳＰから噴出し、試験場所に漂っているガスを示す領域である）が明確に分からない。そこで、領域抽出部９２（図１Ａ）は、画像Ｉｍ１に含まれるガス領域を抽出する（図２のステップＳ１）。詳しく説明する。領域抽出部９２は、各画像Ｉｍ１に含まれるガス領域を見やすくする画像処理をする。図４は、この画像処理がされた画像Ｉｍ２を時系列で示す画像図である。
【００３２】
図３及び図４を参照して、画像Ｉｍ２−１は、画像Ｉｍ１−１に含まれるガス領域を見やすくする画像処理がされた画像である。画像Ｉｍ２−２は、画像Ｉｍ１−２に含まれるガス領域を見やすくする画像処理がされた画像である。画像Ｉｍ２−３は、画像Ｉｍ１−３に含まれるガス領域を見やすくする画像処理がされた画像である。画像Ｉｍ２−４は、画像Ｉｍ１−４に含まれるガス領域を見やすくする画像処理がされた画像である。画像Ｉｍ２の中心部に拡がっている白い領域がガス領域である。ガスは、地点ＳＰを起点にして、右から左へ流れていることが分かる。画像Ｉｍ２において、右側が風上であり、左側が風下となる。ガス領域は、言い換えれば、ガス雲の像である。
【００３３】
ガス領域を見やすくする画像処理は、公知技術を適用することができる。例えば、この公知技術は、国際公開ＷＯ２０１７０７３４３０に開示されている。この国際公開された技術は、本発明者らが発明した技術であり、具体的には、ガス漏れの監視対象を複数の時刻で撮影した赤外画像に対して画像処理をするガス検知用画像処理装置であって、漏れたガスによる温度変化を示す第１の周波数成分データよりも周波数が低く、前記監視対象の背景の温度変化を示す第２の周波数成分データを、前記赤外画像を示す画像データから除く処理をする画像処理部を備えるガス検知用画像処理装置を開示している。
【００３４】
また、本発明者らが発明した技術ではないが、例えば、特開２０１２−５８０９３号公報は、検査対象領域におけるガス漏れを検出するガス漏れ検出装置において、検査対象領域を撮影する赤外線カメラと、赤外線カメラにより撮影された赤外線画像を処理する画像処理部と、を有し、画像処理部は、時系列に並べられた複数の赤外線画像からガス漏れによる動的なゆらぎを抽出するゆらぎ抽出部を有する、ことを特徴とするガス漏れ検出装置を開示している。
【００３５】
領域抽出部９２は、時系列に並ぶ画像Ｉｍ２の中から所定期間（例えば１秒間）の動画を構成する画像Ｉｍ２を選択する。領域抽出部９２は、例えば、時刻Ｔ１の画像Ｉｍ２−１、時刻Ｔ１から０．２秒経過した時刻Ｔ２の画像Ｉｍ２−２、時刻Ｔ１から０．４秒経過した時刻Ｔ３の画像Ｉｍ２−３、時刻Ｔ１から０．６秒経過した時刻Ｔ４の画像Ｉｍ２−４、時刻Ｔ１から０．８秒経過した時刻Ｔ５の画像Ｉｍ２−５（不図示）を選択する。
【００３６】
領域抽出部９２は、これらの画像Ｉｍ２（画像Ｉｍ２−１〜画像Ｉｍ２−５）を用いて、最大画素値で構成される画像Ｉｍ３を生成する。図５は、画像Ｉｍ３の一例を示す画像図である。画像Ｉｍ３の生成について詳しく説明する。領域抽出部９２は、これらの画像Ｉｍ２において、同じ順番に位置する画素の中から、画素が示す値の最大値を決定する。領域抽出部９２は、この最大値を、画像Ｉｍ３の上記順番に位置する画素の値とする。具体的に説明すると、領域抽出部９２は、これらの画像Ｉｍ２において、１番目の画素が示す値の最大値を決定し、この値を、画像Ｉｍ３の１番目の画素の値とする。領域抽出部９２は、これらの画像Ｉｍ２において、２番目の画素が示す値の最大値を決定し、この値を、画像Ｉｍ３の２番目の画素の値とする。領域抽出部９２は、３番目以降の画素についても同様の処理をする。
【００３７】
領域抽出部９２は、画像Ｉｍ３に対してノイズ除去および二値化処理をした画像Ｉｍ４を生成する。図６は、画像Ｉｍ４の一例を示す画像図である。画像Ｉｍ４は、ガス領域２１（漏洩流体領域の一例）が抽出された画像である。
【００３８】
なお、ガス領域２１の抽出方法は、上述した方法に限定されない。また、ユーザがガス領域２１を指定してもよい。これによれば、ガス領域２１の抽出が不要となる。
【００３９】
図４を参照して、ベクトル算出部９３（図１Ａ）は、画像Ｉｍ２を構成する各画素の移動ベクトルを算出する（図２のステップＳ２）。詳しく説明する。ベクトル算出部９３は、異なる時刻の２つの画像Ｉｍ２を用いて、オプティカルフローを生成する。図７は、図４に示す画像Ｉｍ２のオプティカルフローを示す画像Ｉｍ５において、ガス領域２１の一部を拡大した模式図である。ここでは、時刻Ｔ１の画像Ｉｍ２−１から時刻Ｔ４の画像Ｉｍ２−４へ進む場合のオプティカルフローが生成されたとする。各画素に含まれる矢印は、各画素の移動ベクトルを示す。移動ベクトルは、各画素の移動方向と移動量を示す。
【００４０】
なお、移動ベクトルの算出は、オプティカルフローに限定されない。移動ベクトルは、例えば、テンプレートマッチングを用いて算出されてもよい。テンプレートマッチングは、例えば、ＳＡＤ（Ｓｕｍ  ｏｆ  Ａｂｓｏｌｕｔｅ  Ｄｉｆｆｅｒｅｎｃｅ）、ＮＣＣ（Ｎｏｒｍａｌｉｚｅｄ  Ｃｒｏｓｓ−Ｃｏｒｒｅｌａｔｉｏｎ）がある。テンプレートマッチングを用いた移動ベクトルの算出について、図４に示す時刻Ｔ１の画像Ｉｍ２−１と時刻Ｔ２の画像Ｉｍ２−２を例にして詳しく説明する。ベクトル算出部９３は、時刻Ｔ１の画像Ｉｍ２−１に含まれる多数の特徴点（例えば、エッジ）について、時刻Ｔ２の画像Ｉｍ２−２を用いて、テンプレートマッチングし、これらの特徴点の移動ベクトルを算出する。特徴点は、１画素でなく複数の画素で構成される領域である（例えば、５画素×５画素）。従って、例えば、特徴点Ａの移動ベクトルが、特徴点Ａを構成する複数の画素のそれぞれの移動ベクトルとなる。以上のようにして、画像Ｉｍ２−１を構成する各画素の移動ベクトルが算出される。なお、エッジ等の特徴点を用いずに、図４の画素で構成される領域をそのまま処理してもよい。
【００４１】
画素決定処理部９４（図１Ａ）は、移動元画素を追跡する（図２のステップＳ３）。図８は、これを説明するフローチャートである。図９〜図１１は、これを説明する説明図であり、図９は第１説明図、図１０は第２説明図、図１１は第３説明図である。図９〜図１１は、オプティカルフローを示す画像Ｉｍ５（図７）の一部を拡大した図であり、図示の便宜上、画素Ｐと画素Ｐとの間にスペースが設けられている。画素Ｐは、行列状に配置されている。画素Ｐには、移動ベクトルＶが示されている。
【００４２】
画素決定処理部９４（図１Ａ）は、ガス領域２１（図６）を構成する全画素Ｐに対して、移動元画素を追跡する。図８および図９を参照して、ガス領域２１を構成する画素Ｐの一つである画素Ｐ−１を例にして説明する。画素決定処理部９４は、画素Ｐ−１を第１画素として、移動元画素の追跡を開始する。画素決定処理部９４は、第１画素（画素Ｐ−１）の移動ベクトルＶの量が所定のしきい値より小さいか否かを判断する（ステップＳ３０）。
【００４３】
第１画素の移動ベクトルの量（動きの量）が小さければ、この第１画素の付近に漏洩源が位置する可能性がある。そこで、画素決定処理部９４は、第１画素（画素Ｐ−１）の移動ベクトルＶの量が所定のしきい値より小さいと判断したとき（ステップＳ３０でＹｅｓ）、第１画素（画素Ｐ−１）を最終移動元画素と決定する（ステップＳ３１）。所定のしきい値は、例えば、ガス領域２１を構成する全画素Ｐの平均移動ベクトルの量の１０分の１である。
【００４４】
図８および図１０を参照して、画素決定処理部９４は、第１画素（画素Ｐ−１）の移動ベクトルＶの量が所定のしきい値以上と判断したとき（ステップＳ３０でＮｏ）、第１画素（画素Ｐ−１）の移動ベクトルＶ−１の逆ベクトルｖ−１を算出する。画素決定処理部９４は、第１画素（画素Ｐ−１）と隣接する８個の画素Ｐ（画素Ｐ−２〜画素Ｐ−９）の中で、逆ベクトルｖ−１の方向を示す軸を基準にして、±９０度の範囲内にある画素Ｐを、第２画素の候補と決定する（ステップＳ３２）。第１画素（画素Ｐ−１）から見て、第１画素の動きの方向側に位置する画素Ｐ（画素Ｐ−５、画素Ｐ−７、画素Ｐ−８、画素Ｐ−９）は、この第１画素の移動元画素（第２画素）として不適当である。ここでは、画素Ｐ−２、画素Ｐ−３、画素Ｐ−４、画素Ｐ−６が第２画素の候補と決定されたとする。
【００４５】
画素決定処理部９４は、候補（画素Ｐ−２、画素Ｐ−３、画素Ｐ−４、画素Ｐ−６）の中で、ガス領域２１外に位置する候補があれば、この候補を除外する（ステップＳ３３）。例えば、画素Ｐ−２がガス領域２１外に位置するとき、画素Ｐ−２は、候補から除外される。ガスの漏洩源の位置は、通常、ガス領域２１内である。画素決定処理部９４が、ガス領域２１外に位置する画素Ｐを第２画素として、移動元画素を追跡する処理をすれば、最終移動元画素は、ガスの漏洩源の位置から大きく離れる可能性が高い。そこで、画素決定処理部９４は、候補のうち、ガス領域２１外に位置する候補を候補から除外する。
【００４６】
第２画素の候補が１つしかない場合を考える。第２画素の候補が第２画素として不適当であれば（例えば、ガス領域２１外に位置する画素Ｐ）、移動元画素を追跡する処理は、ここで止まり、漏洩源の位置に近づくことができない。実施形態によれば、第２画素の候補が二以上あるので、別の候補を第２画素にして、移動元画素を追跡する処理を続けることができる。
【００４７】
画素決定処理部９４は、候補（画素Ｐ−２、画素Ｐ−３、画素Ｐ−４、画素Ｐ−６）の全てが除外されたか否かを判断する（ステップＳ３４）。画素決定処理部９４は、候補の全てが除外されたと判断したとき（ステップＳ３４でＹｅｓ）、候補の決定に用いられた第１画素（すなわち、候補の中から第２画素が決定されたならば、この第２画素とペアとなる第１画素）を、最終移動元画素と決定する（ステップＳ３５）。第１画素は、画素Ｐ−１なので、画素決定処理部９４は、画素Ｐ−１を最終移動元画素と決定する。
【００４８】
画素決定処理部９４は、候補の全てが除外されたと判断しないとき、すなわち、除外されていない候補があると判断したとき（ステップＳ３４でＮｏ）、除外されなかった候補について、第１画素（画素Ｐ−１）からの移動ベクトルＶを算出する（ステップＳ３６）。ここでは、除外された候補がないとする。すなわち、候補（画素Ｐ−２、画素Ｐ−３、画素Ｐ−４、画素Ｐ−６）の全てが、ガス領域２１内にあるとする。図１２を参照して、画素決定処理部９４は、画素Ｐ−１から画素Ｐ−２へ移動する移動ベクトルＶ−１−２、画素Ｐ−１から画素Ｐ−３へ移動する移動ベクトルＶ−１−３、画素Ｐ−１から画素Ｐ−４へ移動する移動ベクトルＶ−１−４、画素Ｐ−１から画素Ｐ−６へ移動する移動ベクトルＶ−１−６を算出する。図１２は、これらの移動ベクトルＶを説明する説明図である。
【００４９】
画素決定処理部９４は、移動ベクトルＶ−１−２と逆ベクトルｖ−１との角度、移動ベクトルＶ−１−３と逆ベクトルｖ−１との角度、移動ベクトルＶ−１−４と逆ベクトルｖ−１との角度、移動ベクトルＶ−１−６と逆ベクトルｖ−１との角度を算出する。算出された角度の中で、最小の角度を得た候補が、第１画素の移動元画素（第２画素）と見なすのが妥当である。そこで、画素決定処理部９４は、算出された角度の中で、最小の角度を得た候補を、第２画素と決定する（ステップＳ３７）。ここでは、画素Ｐ−３が第２画素と決定される。
【００５０】
このように、画素決定処理部９４は、移動ベクトルＶ−１−ｎ（ここでは、ｎが２、３、４、６）と逆ベクトルｖ−１との角度（軸に対する角度）が最も小さい位置にある候補でもガス領域２１外に位置すれば、候補から除外し（ステップＳ３３）、残された候補の中で、移動ベクトルＶ−１−ｎと逆ベクトルｖ−１との角度が最も小さい位置にある候補を、第１画素の移動元画素（第２画素）に決定する（ステップＳ３７）。よって、ガス領域２１のエッジで移動元画素の追跡が終了し、エッジに最終移動元画素が集まることを防止できる。
【００５１】
以上により、画素Ｐ−１（第１画素）と画素Ｐ−３（第２画素）のペアが形成される。画素決定処理部９４は、このペアの形成方法と同じ方法を用いて、新たなペアを形成する。詳しく説明する。図８および図１０を参照して、画素決定処理部９４は、第１画素（画素Ｐ−３）の移動ベクトルＶの量が所定のしきい値より小さいか否かを判断する（ステップＳ３０）。画素決定処理部９４は、第１画素（画素Ｐ−３）の移動ベクトルＶの量が所定のしきい値より小さいと判断したとき（ステップＳ３０でＹｅｓ）、第１画素（画素Ｐ−３）を最終移動元画素と決定する（ステップＳ３１）。
【００５２】
図８および図１１を参照して、画素決定処理部９４は、第１画素（画素Ｐ−３）の移動ベクトルＶの量が所定のしきい値以上と判断したとき（ステップＳ３０でＮｏ）、第１画素（画素Ｐ−３）の移動ベクトルＶ−３の逆ベクトルｖ−３を算出する。画素決定処理部９４は、第１画素（画素Ｐ−３）と隣接する８個の画素Ｐ（画素Ｐ−１０、画素Ｐ−１１、画素Ｐ−１２、画素Ｐ−２、画素Ｐ−４、画素Ｐ−５、画素Ｐ−１、画素Ｐ−６）の中で、逆ベクトルｖ−３の方向を示す軸を基準にして、±９０度の範囲内にある画素Ｐを、第２画素の候補と決定する（ステップＳ３２）。ここでは、画素Ｐ−１０、画素Ｐ−１１、画素Ｐ−１２、画素Ｐ−４が第２画素の候補と決定されたとする。
【００５３】
画素決定処理部９４は、候補（画素Ｐ−１０、画素Ｐ−１１、画素Ｐ−１２、画素Ｐ−４）の中で、ガス領域２１外に位置する候補があれば、この候補を除外する（ステップＳ３３）。
【００５４】
画素決定処理部９４は、候補（画素Ｐ−１０、画素Ｐ−１１、画素Ｐ−１２、画素Ｐ−４）の全てが除外されたか否かを判断する（ステップＳ３４）。画素決定処理部９４は、候補の全てが除外されたと判断したとき（ステップＳ３４でＹｅｓ）、候補の決定に用いられた第１画素を、最終移動元画素と決定する（ステップＳ３５）。第１画素は、画素Ｐ−３なので、画素決定処理部９４は、画素Ｐ−３を最終移動元画素と決定する。
【００５５】
画素決定処理部９４は、候補の全てが除外されたと判断しないとき、すなわち、除外されなかった候補があると判断したとき（ステップＳ３４でＮｏ）、除外されなかった候補について、第１画素（画素Ｐ−３）からの移動ベクトルＶを算出する（ステップＳ３６）。ここでは、除外された候補がないとする。すなわち、候補（画素Ｐ−１０、画素Ｐ−１１、画素Ｐ−１２、画素Ｐ−４）の全てが、ガス領域２１内にあるとする。図１３を参照して、画素決定処理部９４は、画素Ｐ−３から画素Ｐ−１０へ移動する移動ベクトルＶ−３−１０、画素Ｐ−３から画素Ｐ−１１へ移動する移動ベクトルＶ−３−１１、画素Ｐ−３から画素Ｐ−１２へ移動する移動ベクトルＶ−３−１２、画素Ｐ−３から画素Ｐ−４へ移動する移動ベクトルＶ−３−４を算出する。図１３は、これらの移動ベクトルＶを説明する説明図である。
【００５６】
画素決定処理部９４は、移動ベクトルＶ−３−１０と逆ベクトルｖ−３との角度、移動ベクトルＶ−３−１１と逆ベクトルｖ−３との角度、移動ベクトルＶ−３−１２と逆ベクトルｖ−３との角度、移動ベクトルＶ−３−４と逆ベクトルｖ−３との角度を算出する。画素決定処理部９４は、算出された角度の中で、最小の角度を得た候補を、第２画素と決定する（ステップＳ３７）。ここでは、画素Ｐ−１２が第２画素と決定される。
【００５７】
以上により、画素Ｐ−３（第１画素）と画素Ｐ−１２（第２画素）のペアが形成される。画素決定処理部９４は、最終移動元画素を決定するまで（ステップＳ３１、ステップＳ３５）、ペアの形成の処理を繰り返すことにより、移動元画素を追跡する。図１４は、図７に示す画像Ｉｍ５おいて、画素Ｐ−１について、移動元画素を追跡した結果を説明する説明図である。画素Ｐ−１の最終移動元画素は、画素Ｐ−１００である。
【００５８】
図８を参照して、画素決定処理部９４は、最終移動元画素を決定したとき（ステップＳ３１、ステップＳ３５）、ガス領域２１の全画素Ｐについて、最終移動元画素を決定したか否かを判断する（ステップＳ３８）、画素決定処理部９４は、ガス領域２１の全画素Ｐについて、最終移動元画素を決定していないと判断したとき（ステップＳ３８でＮｏ）、最終移動元画素の決定がされていない画素Ｐに対して、ステップＳ３０〜ステップＳ３７の処理をする。
【００５９】
画素決定処理部９４は、ガス領域２１の全画素Ｐについて、最終移動元画素を決定したと判断したとき（ステップＳ３８でＹｅｓ）、移動元画素の追跡（図２のステップＳ３）が終了する。漏洩源推定部９５（図１Ａ）は、ガスの漏洩源の位置を推定する（図２のステップＳ４）。これについて詳しく説明する。図１５は、図７に示す画像Ｉｍ５において、ガス領域２１内の各画素Ｐと最終移動元画素との関係を説明する説明図である。ガス領域２１内の各画素Ｐと最終移動元画素とが直線でつながれることにより、ガス領域２１内の各画素Ｐと最終移動元画素との関係を示している。ガス領域２１内の全画素Ｐについて直線が示されると、図面が見にくくなるので、ガス領域２１内のいくつかの画素Ｐについて直線が示されている。
【００６０】
最終移動元画素が一つでなく、複数の場合、最頻の最終移動元画素が漏洩源の位置の可能性が高い。そこで、漏洩源推定部９５は、ガス領域２１の全画素Ｐのそれぞれの最終移動元画素において、最頻の最終移動元画素をガスの漏洩源の位置と推定する。例えば、最終移動元画素が画素Ｐ−１００と画素Ｐ−１０１とする。漏洩源推定部９５は、画素Ｐ−１００が最終移動元画素となる画素Ｐの数と、画素Ｐ−１０１が最終移動元画素となる画素Ｐの数とを比較し、多い方を、ガスの漏洩源の位置と推定する。図１５の場合、画素Ｐ−１００がガスの漏洩源の位置と推定される。
【００６１】
実施形態の主な効果を説明する。図８および図１４を参照して、画素決定処理部９４は、ガス領域２１を含む画像Ｉｍ５において、第１画素の移動ベクトルＶ（第１画素の動き）を基にして、第１画素の移動元画素となる第２画素を決定する処理を繰り返すことにより（言い換えれば、第１画素と第２画素のペアを決定する処理を繰り返す）、移動元画素を追跡し、最終移動元画素を決定する。漏洩したガスの流れの開始地点は、ガスの漏洩源なので、移動元画素を追跡することにより、ガスの漏洩源の位置に近づく可能性があり、最終移動元画素は、ガスの漏洩源の位置を示す可能性がある。
【００６２】
漏洩源推定部９５は、ガス領域２１の全画素Ｐの最終移動元画素を基にして、ガスの漏洩源を推定する（図２のステップＳ４）。ユーザは、この推定を参考にして、ガスの漏洩源を特定することができる。従って、実施形態によれば、ガス領域２１を含む画像Ｉｍ５を基にして、比較的簡単にガスの漏洩源の特定を支援することができる。
【００６３】
実施形態は、ガス領域２１を含むフレームで構成される動画を逆再生し、この動画を画像処理するのではなく、ガス領域２１を含む画像Ｉｍ５（静止画）に対して画像処理をする。従って、実施形態によれば、画像処理の量を少なくすることができるので、画像処理の速度を向上させることができる。
【００６４】
実施形態の変形例を説明する。実施形態では、オプティカルフローを示す画像Ｉｍ５が１つである。変形例では、異なる時刻のオプティカルフローを示す複数の画像Ｉｍ５（異なる時刻の画素の動きをそれぞれ示す複数の画像）の場合である。例えば、画素決定処理部９４は、図４に示す時刻Ｔ１の画像Ｉｍ２−１から時刻Ｔ２の画像Ｉｍ２−２へ進む場合のオプティカルフローを示す画像Ｉｍ５（図７）に対して、図２に示すステップＳ２およびステップＳ３の処理をし、かつ、図４に示す時刻Ｔ２の画像Ｉｍ２−２から時刻Ｔ３の画像Ｉｍ２−３へ進む場合のオプティカルフローを示す画像Ｉｍ５−１（不図示）に対して、図２に示すステップＳ２およびステップＳ３の処理をする。すなわち、画素決定処理部９４は、ガス領域２１の全画素Ｐ（複数の画素Ｐ）のそれぞれについて、最終移動元画素を決定することについて、異なる時刻の画素Ｐの動きをそれぞれ示す複数の画像のそれぞれに対して実行する。
【００６５】
漏洩源推定部９５は、オプティカルフローを示す画像Ｉｍ５（図７）から得られた最終移動元画素、および、オプティカルフローを示す画像Ｉｍ５−１（不図示）から得られた最終移動元画素の中で、最頻の最終移動元画素をガスの漏洩源の位置と推定する。すなわち、漏洩源推定部９５は、複数の画像のそれぞれについて得られた最終移動元画素において、最頻の最終移動元画素を前記漏洩源の位置と推定する。
【００６６】
変形例によれば、異なる時刻の画素Ｐの動きをそれぞれ示す複数の画像を用いるので、１つの画像を用いる場合と比べて、ガスの漏洩源の位置を推定する精度を向上させることが可能となる。
【００６７】
（実施形態の纏め）
実施形態の第１局面に係る漏洩源特定支援装置は、漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理部を備え、前記処理部は、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する。
【００６８】
漏洩する流体は、ガスでもよいし、液体でもよい。第２画素は、好ましくは、第１画素の近傍に位置しており、例えば、第１画素と第２の画素との間に画素が存在しなくてもよいし、第１画素と第２画素との間に数画素が存在してもよい。
【００６９】
処理部は、漏洩流体領域を含む画像において、第１画素の動き（例えば、第１画素の移動ベクトル）を基にして、第１画素の移動元画素となる第２画素を決定する処理を繰り返すことにより（言い換えれば、第１画素と第２画素のペアを決定する処理を繰り返す）、移動元画素を追跡し、最終移動元画素を決定する。漏洩流体の流れの開始地点は、漏洩源なので、移動元画素を追跡することにより、漏洩源の位置に近づく可能性があり、最終移動元画素は、漏洩源の位置を示す可能性がある。
【００７０】
処理部は、漏洩流体領域を含む画像を構成する複数の画素（複数の画素は、例えば、その画像の全画素でもよいし、その画像の全画素の中から選択された画素でもよいし、漏洩流体領域の全画素でもよいし、漏洩流体領域の全画素の中から選択された画素でもよい）のそれぞれについて、最終移動元画素を決定する。最終移動元画素の数は、複数の画素の数以下である。例えば、複数の画素が３つの画素とする。３つの画素の最終移動元画素がそれぞれ異なることもあるし、同じこともあるし、３つのうち２つの画素の最終移動元画素が同じこともある。
【００７１】
ユーザは、複数の画素のそれぞれの最終移動元画素を参考にして、漏洩源を特定してもよいし、後で説明するように、漏洩源特定支援装置が、複数の画素のそれぞれの最終移動元画素を基にして、漏洩源を推定してもよい。ユーザは、この推定を参考にして、漏洩源を特定することができる。
【００７２】
以上説明したように、実施形態の第１局面に係る漏洩源特定支援装置によれば、漏洩流体領域を含む画像を基にして、比較的簡単に漏洩源の特定を支援することができる。
【００７３】
上記構成において、前記処理部は、前記第１画素の動きの方向と逆方向を示す軸を基準にして、所定の角度範囲内にある二以上の前記画素を前記第２画素の候補に決定し、前記候補の中から前記第２画素を決定する。
【００７４】
第２画素として不適当な画素がある。例えば、漏洩流体領域外に位置する画素である。第２画素の候補が１つしかない場合を考える。第２画素の候補が第２画素として不適当であれば、移動元画素を追跡する処理は、ここで止まり、漏洩源の位置に近づくことができない。この構成によれば、第２画素の候補が二以上あるので、別の候補を第２画素にして、移動元画素を追跡する処理を続けることができる。
【００７５】
第２画素の候補が、第１画素の動きの方向と逆方向を示す軸を基準にして、所定の角度範囲内にある二以上の画素とされている。この理由は、第１画素から見て、第１画素の動きの方向側に位置する画素は、この第１画素の移動元画素（第２画素）として不適当だからである。所定の角度範囲は、例えば、逆方向を示す軸を基準にして±９０度の範囲である。
【００７６】
上記構成において、前記処理部は、前記軸に対する角度が最も小さい位置にある前記候補を、前記第２画素と決定する。
【００７７】
第２画素の候補の中で、第１画素の動きの方向と逆方向を示す軸に対する角度が最も小さい位置にある候補が、第１画素の移動元画素（第２画素）と見なすのが妥当である。この構成では、その候補を第２画素と決定する。
【００７８】
上記構成において、前記処理部は、前記候補のうち、前記漏洩流体領域外に位置する前記候補を前記候補から除外する。
【００７９】
漏洩源の位置は、通常、漏洩流体領域内である。処理部が、漏洩流体領域外に位置する画素を第２画素として、移動元画素を追跡する処理をすれば、最終移動元画素は、漏洩源の位置から大きく離れる可能性が高い。そこで、この構成では、第２画素の候補のうち、漏洩流体領域外に位置する候補を候補から除外する。
【００８０】
この構成と、前記軸に対する角度が最も小さい位置にある前記候補を、前記第２画素と決定する構成との組み合わせによれば、次の作用効果が生じる。この組み合わせによれば、軸に対する角度が最も小さい位置にある候補でも漏洩流体領域外に位置すれば、候補から除外し、残された候補の中で、軸に対する角度が最も小さい位置にある候補を、第１画素の移動元画素（第２画素）に決定する。これにより、漏洩流体領域のエッジで移動元画素の追跡が終了し、エッジに最終移動元画素が集まることを防止できる。
【００８１】
上記構成において、前記処理部は、前記候補の全てが前記漏洩流体領域外に位置するとき、前記候補の決定に用いられた前記第１画素を前記最終移動元画素と決定する。
【００８２】
漏洩流体領域外に位置する候補を第２画素とすることは不適当である。そこで、候補の全てが漏洩流体領域外に位置するとき、処理部は、これらの候補の決定に用いられた第１画素（すなわち、候補の中から第２画素が決定されたならば、この第２画素とペアとなる第１画素）を最終移動元画素に決定する。
【００８３】
上記構成において、前記処理部は、前記第１画素の動きの量が所定のしきい値より小さいとき、当該第１画素を前記最終移動元画素と決定する。
【００８４】
第１画素の動きの量が小さければ、この第１画素の付近に漏洩源が位置する可能性がある。そこで、処理部は、その場合、当該第１画素を最終移動元画素と決定する。所定のしきい値は、例えば、漏洩流体領域を構成する全画素の平均移動ベクトルの量の１０分の１である。
【００８５】
上記構成において、前記画像を構成する複数の前記画素のそれぞれの前記最終移動元画素において、最頻の前記最終移動元画素を漏洩源の位置と推定する推定部をさらに備える。
【００８６】
最終移動元画素が一つでなく、複数の場合、最頻の最終移動元画素が漏洩源の位置の可能性が高い。そこで、推定部は、最頻の最終移動元画素を漏洩源の位置と推定する。
【００８７】
上記構成において、前記処理部は、前記画像を構成する複数の前記画素のそれぞれについて、前記最終移動元画素を決定することについて、異なる時刻の前記画素の動きをそれぞれ示す複数の前記画像のそれぞれに対して実行し、前記推定部は、複数の前記画像のそれぞれについて得られた前記最終移動元画素において、最頻の前記最終移動元画素を前記漏洩源の位置と推定する。
【００８８】
異なる時刻の画素の動きをそれぞれ示す複数の画像について具体例で説明する。例えば、時刻Ｔ１に撮影された画像と時刻Ｔ２に撮影された画像とで生成されたオプティカルフローを示す画像、および、時刻Ｔ２に撮影された画像と時刻Ｔ３に撮影された画像とで生成されたオプティカルフローを示す画像の場合、この２つのオプティカルフローを示す画像が、異なる時刻の画素の動きをそれぞれ示す複数の画像である。
【００８９】
この構成によれば、異なる時刻の画素の動きをそれぞれ示す複数の画像を用いるので、１つの画像を用いる場合と比べて、漏洩源の位置を推定する精度を向上させることが可能となる。
【００９０】
実施形態の第２局面に係る漏洩源特定支援方法は、漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理ステップを備え、前記処理ステップは、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する。
【００９１】
実施形態の第２局面に係る漏洩源特定支援方法は、実施形態の第１局面に係る漏洩源特定支援装置を方法の観点から規定しており、実施形態の第１局面に係る漏洩源特定支援装置と同様の作用効果を有する。
【００９２】
実施形態の第３局面に係る漏洩源特定支援プログラムは、漏洩流体領域を含む画像を構成する画素の１つである第１画素の動きを基にして、前記画像を構成する前記画素の１つであり、前記第１画素の移動元画素となる第２画素を決定する処理をし、前記第２画素を新たな前記第１画素として、前記処理を繰り返すことにより、最終移動元画素を決定する処理ステップを、コンピュータに実行させ、前記処理ステップは、前記画像を構成する複数の前記画素のそれぞれについて、前記処理を繰り返すことにより、複数の前記画素のそれぞれについて、前記最終移動元画素を決定する。
【００９３】
実施形態の第３局面に係る漏洩源特定支援プログラムは、実施形態の第１局面に係る漏洩源特定支援装置をプログラムの観点から規定しており、実施形態の第１局面に係る漏洩源特定支援装置と同様の作用効果を有する。
【００９４】
本発明の実施形態が詳細に図示され、かつ、説明されたが、それは単なる図例及び実例であって限定ではない。本発明の範囲は、添付されたクレームの文言によって解釈されるべきである。
【００９５】
２０１８年１１月２７日に提出された日本国特許出願特願２０１８−２２１２１５は、その全体の開示が、その全体において参照によりここに組み込まれる。
【産業上の利用可能性】
【００９６】
本発明によれば、漏洩源特定支援装置、漏洩源特定支援方法、及び、漏洩源特定支援プログラムを提供することができる。
【図１Ａ】
2020110410.tif 000003
【図１Ｂ】
2020110410.tif 000004
【図２】
2020110410.tif 000005
【図３】
2020110410.tif 000006
【図４】
2020110410.tif 000007
【図５】
2020110410.tif 000008
【図６】
2020110410.tif 000009
【図７】
2020110410.tif 000010
【図８】
2020110410.tif 000011
【図９】
2020110410.tif 000012
【図１０】
2020110410.tif 000013
【図１１】
2020110410.tif 000014
【図１２】
2020110410.tif 000015
【図１３】
2020110410.tif 000016
【図１４】
2020110410.tif 000017
【図１５】
2020110410.tif 000018
【国際調査報告】
2020110410.tif 000019
2020110410.tif 000020
2020110410.tif 000021
2020110410.tif 000022