JP7623616B2 - 処理装置、処理方法およびプログラム - Google Patents

Description

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

衛星に搭載されるセンサは、地球を取り巻く大気、地表面、または海面等からの電磁波の放射または反射を計測する。センサは、対象物の物理量または状態を実測ではない。センサが測定した電磁波から、所望の値に変換するモデルが用いられる。気象予測の分野において、データ同化が発展している（非特許文献１）。データ同化は、観測と数値モデルの間を相補的につなぐ。

気象庁、"第３章 データ同化"、［online］、［2021年5月27日検索］、インターネット〈URL：https://www.jma.go.jp/jma/kishou/books/nwptext/45/1_chapter3.pdf〉

自然環境の変化によって生じるノイズにより、モデルの精度が十分に得られない場合がある。この問題を解決するため、モデル適用前のキャリブレーション等で利用するための実測データが必要とされている。

本発明は、上記事情に鑑みてなされたものであり、本発明の目的は、データ同化において、モデル適用前のキャリブレーションに適したデータの観測地点を効率的に選定可能な技術を提供することである。

本発明の一態様の処理装置は、所定時刻に撮影された衛星画像データから分割された、画素の特徴値が類似する複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出する抽出部を備える。

本発明の一態様の処理方法は、コンピュータが、所定時刻に撮影された衛星画像データから分割された、画素の特徴値が類似する複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出する。

本発明の一態様は、上記処理装置として、コンピュータを機能させるプログラムである。

本発明によれば、データ同化において、モデル適用前のキャリブレーションに適したデータの観測地点を効率的に選定可能な技術を提供することができる。

図１は、本発明の実施の形態に係る処理装置の機能ブロックを説明する図である。 図２は、衛星画像データの一例を説明する図である。 図３は、衛星画像データから分割されたセグメントの一例を説明する図である。 図４は、代表値の時間変化を説明する図である。 図５は、代表値の区分の一例を説明する図である。 図６は、プラットフォームの航路の一例を説明する図である。 図７は、本発明の実施の形態に係る処理装置による処理方法を説明するフローチャートである。 図８は、処理装置に用いられるコンピュータのハードウエア構成を説明する図である。

以下、図面を参照して、本発明の実施形態を説明する。図面の記載において同一部分には同一符号を付し説明を省略する。

（処理装置）
本発明の実施の形態に係る処理装置１は、衛星が撮影した衛星画像データ１１から、データ同化のモデル適用前のキャリブレーションで利用する実測データを取得する実測セグメントと、その実測セグメントにおける実測ポイントを特定する。また処理装置１は、シミュレーションで利用する実測データを取得する実測セグメントと、その実測セグメントにおける実測ポイントを特定する。実測データは、所定の拠点から実測ポイントまでドローン等のプラットフォームを飛行させ、実測ポイントにおいてプラットフォーム搭載された測定装置で測定される。

図１に示すように処理装置１は、衛星画像データ１１、セグメントデータ１２、実測セグメントデータ１３および実測ポイントデータ１４の各データと、分割部２１、抽出部２２、特定部２３および出力部２４の各機能を備える。各データは、メモリ９０２またはストレージ９０３等の記憶装置に記憶される。各機能は、CPU９０１に実装される。

衛星画像データ１１は、衛星から撮影された地球上の画像である。衛星画像データ１１は、各画素について、複数の電磁波のバンドの値を有する。具体的には、各画素は、赤（Ｒ）、緑（Ｇ）、および青（青）の可視光の値のみならず、近赤外線、熱赤外線等の各バンドの値を有する。画素が有する各バンドの値は、特徴値とも称される。

衛星画像データ１１は、図２に示すように、複数の時刻および複数の高度で撮影された、複数の画像データの集合であっても良い。衛星画像データ１１は、衛星の軌道によって所定時間間隔に撮影される。図２に示す例において、衛星画像データ１１は、時刻ｔ－１で５つの高度で撮影された各画像のデータ、時刻ｔで５種の高度で撮影された各画像のデータ、および時刻ｔ＋１で５つの高度で撮影された各画像のデータを有する。

セグメントデータ１２は、所定の時刻かつ高度で撮影された衛星画像データ１１から分割された複数のセグメントを特定する。セグメントデータ１２は、例えば、撮影時刻および高度のほか、セグメントの識別子と、そのセグメントに属する画素の識別子を対応づける。セグメントデータ１２は、分割部２１によって生成される。衛星画像データ１１は、図３に示すように、複数のセグメントに分割される。図３に示す例において、衛星画像データ１１は、５つのセグメントに分割されるが、セグメント数はこれに限るものではない。

実測セグメントデータ１３は、分割部２１が特定する複数のセグメントのうち、キャリブレーション用またはシミュレーション用に実測対象となるセグメントを特定する。実測セグメントデータ１３は、例えば、セグメントの識別子および高度などを含む。実測セグメントデータ１３は、抽出部２２によって生成される。

実測ポイントデータ１４は、実測対象のセグメントのうち、実測対象のポイントを特定する。実測ポイントデータ１４は、特定部２３によって生成される。実測ポイントデータ１４は、実測対象の画素の識別子、またはその画素に対応する緯度および経度を含む。また実測ポイントデータ１４は、高度等の位置情報を含んでも良い。

分割部２１は、所定時刻に撮影された衛星画像データ１１を複数のセグメントに分割して、セグメントデータ１２を生成する。分割部２１は、例えば、画素の特徴値が類似し、かつ隣接する画素を、セグメントとして特定する（例えば、「高分解能衛星画像のテクスチャ特徴量とスペクトル特徴量を用いたオブジェクト指向型林分タイプ分類」参照（ＵＲＬ：https://www.jstage.jst.go.jp/article/jsprs1975/46/2/46_2_27/_pdf/-char/ja））。

抽出部２２は、所定時刻に撮影された衛星画像データ１１から分割された複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出する。抽出されるセグメントは、測定対象以外から電磁波を除去するキャリブレーション、またはシミュレーションに用いる実測データを取得する対象である。実測データは、抽出されたセグメント内の実測ポイントで取得されても良い。ここで「代表値が異なる」とは、代表値の値が同一でない場合のみならず、代表値の値が所定値以上離れる場合を含んでも良い。抽出部２２は例えば、算出された代表値を所定順序で並べ、順序の離れる２つのセグメントを抽出する。抽出部２２は、抽出したセグメントの識別子を含む実測セグメントデータ１３を出力する。

まず、キャリブレーションに用いる複数のセグメントを抽出する処理について説明する。

抽出部２２は、複数のセグメントのそれぞれについて、そのセグメントに属する画素の特徴値から代表値を算出する。代表値は、セグメントに属する各画素の特徴値の平均値または中央値、セグメントの重心における画素の特徴値等であって、セグメントの属性を特徴づける値である。各画素の特徴値が、バンド毎の値などの複数の要素を含む場合、代表値は、その複数の要素で特定されるベクトルの大きさである。

抽出部２２は、各セグメントについて代表値を算出すると、その代表値を所定順序で並べる。抽出部２２は、代表値がもっとも大きいセグメントと、最も小さいセグメントなど、代表値が乖離する２つのセグメントを抽出する。抽出部２２は、衛星画像データ１１に含まれる複数のセグメントのうち、ドローンが移動可能なセグメントから、キャリブレーションに用いる複数のセグメントを抽出しても良い。

例えば抽出部２２は、代表値を複数の区分に分割し、最大値および最小値を有する各区分からそれぞれセグメントを抽出する。最大値を有する区分に属するセグメントが１つの場合、抽出部２２は、そのセグメントを抽出する。同様に、最小値を有する区分に属するセグメントが１つの場合、抽出部２２は、そのセグメントを抽出する。処理装置１は、代表値が乖離する２つのセグメントからキャリブレーション用に実測するポイントを選択することにより、安定したキャリブレーションを支援することができる。

最大値または最小値を有する区分に属するセグメントが複数の場合、抽出部２２は、その複数のセグメントから実測セグメントを１つ選択する。抽出部２２は、時間変化または距離に着目して、１つの実測セグメントを選択する。

まず時間変化に着目して実測セグメントを選択する処理を説明する。最大値または最小値を有する区分に属するセグメントが複数の場合、抽出部２２は、代表値が、所定時刻とは異なる時刻に撮影された衛星画像データにおける代表値との差分の小さいセグメントを、抽出しても良い。最大値を有する区分に属するセグメントが複数の場合、抽出部２２は、その複数のセグメントのうち代表値の時間変化の少ないセグメントを、実測セグメントとして抽出する。同様に、最大値を有する区分に属するセグメントが複数の場合、抽出部２２は、その複数のセグメントのうち代表値の時間変化の少ないセグメントを、実測セグメントとして抽出する。処理装置１は、既に選択されたセグメントと代表値が乖離するセグメントの候補が複数ある場合、時間変化のより小さいセグメントを抽出する。処理装置１は、時間変化のより小さいセグメントからキャリブレーション用に実測するポイントを選択することにより、安定したキャリブレーションを支援することができる。なお抽出部２２は、撮影時刻が異なる２つの衛星画像データから代表値の時間変化を特定しても良いし、撮影時刻が異なる３つ以上の衛星画像データから代表値の時間変化を特定しても良い。

ここで、図４を参照して時間変化を説明する。本発明の実施の形態において、衛星画像データ１１は、時間毎の同一の緯度経度かつ高度の衛星画像のデータを有する。各セグメントは、その条件によって、特徴値が時間毎に変化する場合がある。図４は、２つのセグメントの時間毎の代表値の変化を示す。図４において実線は時間毎の変化が小さく、時間変化が小さいことを示す。代表値の時間変化の小さいセグメントは、安定的な補正が可能になるので、キャリブレーションの実測データの測定に好適である。一方、破線は、時間毎の代表値の変化が大きく、時間変化が大きいことを示す。代表値の時間変化の小さいセグメントは、多種多様なデータをシミュレーションに反映させることができるので、シミュレーションの実測データの測定に好適である。

例えば図５に示すように、抽出部２２は、代表値の幅を、１０などの複数の区分に分割し、各区分に、その区分に属する代表値を有するセグメントの識別子を対応づける。図５に示す例において、最大値を有する区分１と最小値を有する区分１０のそれぞれから、セグメントが抽出される。最大値を有する区分１に属するセグメント数は１なので、抽出部２２は、区分１について、セグメントＡを実測セグメントとして抽出する。

一方最小値を有する区分１０に属するセグメント数は２である。そこで抽出部２２は、他の時刻に撮影された衛星画像データを参照して、セグメントＢとセグメントＣのそれぞれにおける代表値を算出する。抽出部２２は、セグメントＢとセグメントＣのうち、衛星画像データ１１における代表値と他の時刻に撮影された衛星画像データにおける代表値との差分の小さい方のセグメントを抽出する。

次に距離に着目して実測セグメントを選択する処理を説明する。最大値または最小値を有する区分に属するセグメントが複数の場合、もう一方の区分から選択されるセグメントと距離的に近いセグメントを抽出しても良い。例えば、最大値を有する区分するセグメントが１つで、最小値を有する区分に属するセグメントが複数の場合、抽出部２２は、最小値を有する区分に属する複数のセグメントのうち最大値を有する区分に属するセグメントに近いセグメントを、実測セグメントとして抽出する。これにより、実測におけるドローンの飛行距離を抑制することができる。

次に、シミュレーションに用いる複数のセグメントを抽出する処理について説明する。

抽出部２２は、代表値が、所定時刻とは異なる時刻に撮影された衛星画像データ１１における代表値との差分の大きいセグメントを、シミュレーションのための実測セグメントとして抽出する。抽出部２２は、衛星画像データ１１に含まれる複数のセグメントのうち、ドローンが移動可能なセグメントから、シミュレーションに用いる複数のセグメントを抽出しても良い。抽出部２２は、抽出したセグメントの識別子を含む実測セグメントデータ１３を出力する。時間変化の大きいセグメントをシミュレーションのための実測セグメントとすることにより、処理装置１は、多様な実測値を、シミュレーションに反映させることができる。

また抽出部２２は、キャリブレーション用に選択されたセグメントのうちの１つを、シミュレーション用のセグメントとしても良い。これにより、ドローンの移動距離を抑制して、効率的に実測値を収集することができる。

特定部２３は、実測セグメントデータ１３で特定されたセグメントのうち、測定器による実測ポイントを特定する。特定部２３は、特定したポイントを実測ポイントデータ１４として出力する。特定部２３は、衛星画像データ１１における画素を特定し、その画素そのものを、実測ポイントとして出力しても良いし、その画素に対応する緯度経度を、実測ポイントとして出力しても良い。

まず、キャリブレーションに用いる実測ポイントを特定する処理について説明する。

特定部２３は、キャリブレーション用に抽出された２つのセグメントのそれぞれから、実測ポイントを特定する。特定部２３は、抽出されたセグメントにおいて、隣接するセグメントから所定の画素数離れた各画素のうち、所定時刻とは異なる時刻に撮影された衛星画像データにおける特徴値との差分の小さい画素を特定する。特定部２３は、特定された画素の位置を、キャリブレーションのための実測ポイントとして特定する。

キャリブレーション用の実測ポイントの特定において特定部２３は、まず、セグメントの境界から数画素以内の画素を、実測ポイントの候補から除外する。セグメントの境界に近い画素は、隣接するセグメントの特徴値の影響を受ける可能性がある。このような画素を実測ポイントから除外することにより、特定部２３は、セグメントの特徴をより反映した画素から、実測ポイントを選択することができる。

特定部２３は、セグメントの境界に近い画素を除外した各画素のうち、時間変化の小さい画素を、キャリブレーション用に測定するポイントとして特定する。時間変化の小さい画素を、キャリブレーション用の実測ポイントとすることにより、処理装置１は、安定したキャリブレーションを支援することができる。

次に、シミュレーションに用いる実測ポイントを特定する処理について説明する。

特定部２３は、キャリブレーション用に抽出されたセグメント、またはシミュレーション用に抽出されたセグメントから、シミュレーションのための実測ポイントを特定する。特定部２３は、抽出されたセグメントにおいて、隣接するセグメントから所定の画素数離れた各画素のうち、所定時刻とは異なる時刻に撮影された衛星画像データ１１における特徴値との差分の大きい画素を特定する。特定部２３は、特定された画素の位置を、シミュレーションのための実測ポイントとして特定する。

シミュレーション用の実測ポイントの特定において特定部２３は、まず、セグメントの境界に近い画素を、実測ポイントの候補から除外する。セグメントの境界に近い画素は、隣接するセグメントの特徴値の影響を受ける可能性がある。このような画素を実測ポイントから除外することにより、特定部２３は、セグメントの特徴をより反映した画素から、実測ポイントを選択することができる。

特定部２３は、セグメントの境界に近い画素を除外した各画素のうち、時間変化の大きい画素を、シミュレーション用に測定するポイントとして特定する。時間変化の大きい画素を、シミュレーション用の実測ポイントとすることにより、処理装置１は、多様な実測値を用いたシミュレーションを支援することができる。

出力部２４は、実測セグメントデータ１３および実測ポイントデータ１４を出力する。

出力部２４は、抽出部２２によって、キャリブレーション用に抽出された２つのセグメントを、測定対象以外からの電磁波を除去するキャリブレーションのための実測セグメントとして出力する。出力部２４はさらに、特定部２３によって、キャリブレーション用に抽出された２つのセグメントのそれぞれから特定された実測ポイントを出力する。

出力部２４は、抽出部２２によって、シミュレーション用に抽出されたセグメントを、シミュレーションのための実測セグメントとして出力する。出力部２４はさらに、特定部２３によって、シミュレーション用に抽出された２つのセグメントのそれぞれから特定された実測ポイントを出力する。

出力部２４は、実測セグメントデータ１３に含まれるセグメントとして、少なくとも、キャリブレーション用の実測セグメント２つと、シミュレーション用の実測セグメント１つを出力する。シミュレーションようの実測セグメントは、キャリブレーション用の実測セグメント２つのうちの１つと同じでも良いし、異なっても良い。

出力部２４は、実測ポイントデータ１４に含まれるポイントとして、キャリブレーション用の実測セグメント２つのそれぞれから特定された２つの実測ポイントと、シミュレーション用の実測セグメントから特定された１つの実測ポイントを出力する。なお、本発明の実施の形態においてキャリブレーション用のセグメントとシミュレーション用のセグメントが共通する場合でも、キャリブレーション用の実測ポイントとシミュレーション用の実測ポイントは異なる。キャリブレーション用の実測ポイントは時間変化が小さく、シミュレーション用の実測ポイントは時間変化が大きい。

出力部２４は、例えば、キャリブレーション用の実測ポイント２つと、シミュレーション用の実測ポイント１つに移動するための航路を出力しても良い。その場合、航路は、キャリブレーション用の１つの実測ポイント、シミュレーション用の実測ポイント、およびキャリブレーション用のもう１つの実測ポイントの順となることが好ましい。キャリブレーション用の２つのポイントに時間差を設けることにより、より安定したデータを取得することが可能になる。またキャリブレーション用に２つ目の実測ポイントを最後に測定するように航路を設定しても良い。キャリブレーション用およびシミュレーション用にそれぞれ１ポイントずつ測定することで、最小限のデータを取得して、測定の目的を達成することができる。またドローンなどのプラットフォームの動力等に余裕がある場合にキャリブレーション用に２つめのポイントを測定することで、より精度の高いキャリブレーションを支援することができる。

図６に示すように、１つの実測ポイントに対し、八の字を描くように航路が設定される。図６において航路は、（１）から（８）の順序で設定される。これにより、プラットフォームにかかる力学的な隔たりを軽減することができる。具体的には、実測ポイントＰを中心とした四角形を平均的に観測できるように、四角形内において対角線を横切り、かつ、外周も順次測定する。

出力部２４は、衛星が新たに撮影する衛星画像データと同期できる時間内に、プラットフォームによる測定が完了するように、プラットフォームの航路、移動速度、台数等を設定する。ここで、「衛星画像データと同期できる」とは、衛星が衛星画像データの新たな撮影に要する時間内に、プラットフォームによる測定が終了することを意味する。実測ポイントＰの実測データと衛星によって新たに撮影された衛星画像データとの連携により、キャリブレーションまたはシミュレーション等が可能になる。

処理装置１は、出力部２４が出力した航路に従ってプラットフォームが移動しているか否かを監視する制御部（図示せず）を備えても良い。制御部は、プラットフォームの位置を追跡し、航路から所定距離離れると、航路から離れた旨の警告を出力して、航路に戻るよう促しても良い。

図７を参照して、本発明の実施の形態に係る処理方法を説明する。

まずステップＳ１において処理装置１は、所定時刻の所定高度の衛星画像データを、複数のセグメントに分割する。

ステップＳ２において処理装置１は、各セグメントの代表値を基準に、キャリブレーション用の２つ以上のセグメントを抽出する。ステップＳ３において処理装置１は、ステップＳ２において抽出した各セグメントから、キャリブレーション用の実測ポイントを特定する。

ステップＳ４において処理装置１は、各セグメントの代表値を基準に、シミュレーション用のセグメントを抽出する。ステップＳ５において処理装置１は、ステップＳ４において抽出したセグメントから、シミュレーション用の実測ポイントを特定する。

ステップＳ６において処理装置１は、ステップＳ３およびステップＳ５で特定した実測ポイントを通過する航路のデータを出力する。

複数の高度について実測ポイントを特定する場合、処理装置１は、各高度の衛星画像データについて、ステップＳ１ないしステップＳ５の処理を繰り返す。

本発明の実施の形態に係る処理装置１は、衛星などの常設の観測システムで取得さあれる相対的な広域の衛星画像データを用いて。キャリブレーション等の所定の目的に合致したデータを局所的にかつ動的に取得することが可能になる、これにより、観測コストを低減し、かつ所望の衛星画像データに好適な実データの測定が可能になる。

また処理装置１は、衛星画像データにおいて特徴値が類似する複数のセグメントに分割し、複数のセグメントから、測定の目的に合致したセグメントを抽出する。処理装置１は、抽出したセグメントから、測定に好適な実測ポイントを特定すること。これにより処理装置１は、所望の実測ポイントを効率的に選択することが可能になる。

本発明の実施の形態に係る処理装置１は、データ同化において、モデル適用前のキャリブレーションに適したデータの観測地点を効率的に選定することができる。

上記説明した本実施形態の処理装置１は、例えば、CPU（Central Processing Unit、プロセッサ）９０１と、メモリ９０２と、ストレージ９０３（HDD：Hard Disk Drive、SSD：Solid State Drive）と、通信装置９０４と、入力装置９０５と、出力装置９０６とを備える汎用的なコンピュータシステムが用いられる。このコンピュータシステムにおいて、CPU９０１がメモリ９０２上にロードされたプログラムを実行することにより、処理装置１の各機能が実現される。

なお、処理装置１は、１つのコンピュータで実装されてもよく、あるいは複数のコンピュータで実装されても良い。また処理装置１は、コンピュータに実装される仮想マシンであっても良い。

処理装置１のプログラムは、HDD、SSD、USB（Universal Serial Bus）メモリ、CD (Compact Disc)、DVD (Digital Versatile Disc)などのコンピュータ読取り可能な記録媒体に記憶することも、ネットワークを介して配信することもできる。

なお、本発明は上記実施形態に限定されるものではなく、その要旨の範囲内で数々の変形が可能である。

１ 処理装置
１１ 衛星画像データ
１２ セグメントデータ
１３ 実測セグメントデータ
１４ 実測ポイントデータ
２１ 分割部
２２ 抽出部
２３ 特定部
２４ 出力部
９０１ CPU
９０２ メモリ
９０３ ストレージ
９０４ 通信装置
９０５ 入力装置
９０６ 出力装置

Claims (8)

1. 所定時刻に撮影された衛星画像データから分割された、画素の特徴値が類似する複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出する抽出部を備え、
   前記２つのセグメントは、測定対象以外から電磁波を除去するキャリブレーションのための実測セグメントであって、
   前記抽出部は、
   前記代表値を複数の区分に分割し、最大値および最小値を有する各区分からそれぞれセグメントを抽出し、
   前記最大値または最小値を有する区分に属するセグメントが複数の場合、前記代表値が、前記所定時刻とは異なる時刻に撮影された衛星画像データにおける代表値との差分の小さいセグメントを、抽出する
   処理装置。
2. 所定時刻に撮影された衛星画像データから分割された、画素の特徴値が類似する複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出する抽出部を備え、
   前記２つのセグメントは、測定対象以外から電磁波を除去するキャリブレーションのための実測セグメントであって、
   前記抽出部は、
   前記代表値を複数の区分に分割し、最大値および最小値を有する各区分からそれぞれセグメントを抽出し、
   前記最大値または最小値を有する区分に属するセグメントが複数の場合、もう一方の区分から選択されるセグメントと距離的に近いセグメントを抽出する
   処理装置。
3. 前記２つのセグメントは、測定対象以外から電磁波を除去するキャリブレーションのための実測セグメントであって、
   抽出されたセグメントにおいて、隣接するセグメントから所定の画素数離れた各画素のうち、前記所定時刻とは異なる時刻に撮影された衛星画像データにおける特徴値との差分の小さい画素を特定し、特定された画素の位置を、実測ポイントとして特定する特定部
   をさらに備える請求項１または２に記載の処理装置。
4. 抽出されたセグメントにおいて、隣接するセグメントから所定の画素数離れた各画素のうち、前記所定時刻とは異なる時刻に撮影された衛星画像データにおける特徴値との差分の大きい画素を特定し、特定された画素の位置を、シミュレーションのための実測ポイントとして特定する特定部
   をさらに備える請求項１または２に記載の処理装置。
5. 前記抽出部は、さらに、
   前記代表値が、前記所定時刻とは異なる時刻に撮影された衛星画像データにおける代表値との差分の大きいセグメントを、シミュレーションのための実測セグメントとして抽出する
   請求項１または２に記載の処理装置。
6. コンピュータが、所定時刻に撮影された衛星画像データから分割された、画素の特徴値が類似する複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出し、
   前記２つのセグメントは、測定対象以外から電磁波を除去するキャリブレーションのための実測セグメントであって、
   前記コンピュータが、
   前記代表値を複数の区分に分割し、最大値および最小値を有する各区分からそれぞれセグメントを抽出し、
   前記最大値または最小値を有する区分に属するセグメントが複数の場合、前記代表値が、前記所定時刻とは異なる時刻に撮影された衛星画像データにおける代表値との差分の小さいセグメントを、抽出する
   処理方法。
7. コンピュータが、所定時刻に撮影された衛星画像データから分割された、画素の特徴値が類似する複数のセグメントのそれぞれについて代表値を算出し、算出された代表値が異なる２つのセグメントを抽出し、
   前記２つのセグメントは、測定対象以外から電磁波を除去するキャリブレーションのための実測セグメントであって、
   前記コンピュータが、
   前記代表値を複数の区分に分割し、最大値および最小値を有する各区分からそれぞれセグメントを抽出し、
   前記最大値または最小値を有する区分に属するセグメントが複数の場合、もう一方の区分から選択されるセグメントと距離的に近いセグメントを抽出する
   処理方法。
8. コンピュータを、請求項１または請求項２に記載の処理装置として機能させるためのプログラム。

Images