JP2010042774A - 位置誤差補正装置及び位置誤差補正プログラム - Google Patents

Abstract

【課題】衛星画像に関する光学マルチセンサデータのシステム補正処理から幾何モデルパラメータ調整までの一連の処理を再帰的に自動で行う装置を提供する。
【解決手段】画像処理部１０は、光学マルチセンサデータ、人工衛星の軌道・姿勢データ、幾何モデルパラメータに基づき、複数帯域の衛星画像データを作成する。バンド間位置誤差算出部２０は、所定の帯域の衛星画像データを基準画像データとし、基準画像データ以外の各帯域の衛星画像データを参照画像データとし、基準画像データと各参照画像データとの間の位置誤差をバンド間位置誤差として参照画像データごとに算出する。幾何モデルパラメータ算出部４０は、各バンド間位置誤差に基づき、新たな幾何モデルパラメータを算出する。画像処理部１０は、前記光学マルチセンサデータ、前記軌道・姿勢衛星データ、新たに算出された幾何モデルパラメータとから新たに複数帯域の衛星画像データを作成する。
【選択図】図３

Description

この発明は、衛星画像におけるバンド間の位置合わせの補正処理に関わり、特に、光学マルチセンサによって得られた画像データのバンド間の位置ずれ補正に関するものである。

光学センサの打ち上げ初期の衛星画像データにおいて、地球観測衛星を打ち上げ時の振動による、撮像素子の位置ずれ、あるいは光学センサの取り付け位置の微小な位置ずれが要因となり、各バンド（帯域）のフィルタをもつ光学センサで取得した画像データを重ね合わせた場合、位置ずれが発生することが知られている。

上記のような画像データの位置ずれを補正する方法として、地上での画像処理による補正は行わず、観測衛星のオンボードにて、撮像素子の光学系指向方向のずれを補正する従来技術が知られている（例えば、特許文献１）。
特開平６−８８９７号公報

従来は、光学マルチセンサデータの補正処理、バンド間の位置誤差算出、位置誤差の精度評価、幾何モデルパラメータ調整を、作業者の手作業で実施していたため、位置誤差の目標精度を満たすまで、何度も繰り返し実施する必要があり、作業者への負担が大きく、かつ時間がかかるという課題があった。

また、特許文献１の場合、オンボード装置に障害が発生した場合、観測衛星のオンボード装置を修理することができず、光学センサ（撮像素子）の位置ずれを補正できないという課題がある。その場合、地上において位置ずれの修正を行わなくてはならない。

この発明は、バンド間の位置誤差精度を満たす幾何モデルパラメータを作成するため、上記に示す光学マルチセンサデータのシステム補正処理から幾何モデルパラメータ調整までの処理を再帰的に行う装置を提供し、作業者の負担を大幅に軽減することを目的とする。

この発明の位置誤差補正装置は、
所定の帯域のフィルタを持つ撮像素子である複数の光学センサから構成される光学マルチセンサによって地球の領域を対象として撮像された光学マルチセンサデータと、前記光学マルチセンサが搭載された人工衛星の軌道データと姿勢データとを含む衛星データと、前記衛星データに含まれる前記軌道データと前記姿勢データとから地球に対する前記光学マルチセンサの向きの特定に使用するパラメータである幾何モデルパラメータとを取得し、前記光学マルチセンサデータと、前記衛星データと、前記幾何モデルパラメータとに基づき、地球の同一領域における複数の帯域の衛星画像データを作成する画像処理部と、
前記画像処理部によって作成された前記複数の帯域の衛星画像データを取り込み、所定の帯域の衛星画像データを基準画像データとし、前記基準画像データ以外のそれぞれの帯域の前記衛星画像データを参照画像データとし、前記基準画像データと各参照画像データとの間の位置誤差をバンド間位置誤差として前記参照画像データごとに算出するバンド間位置誤差算出部と、
前記バンド間位置誤差算出部によって算出されたそれぞれの前記バンド間位置誤差に基づき、前記幾何モデルパラメータに代替する新たな幾何モデルパラメータを算出する幾何モデルパラメータ算出部と
を備え、
前記画像処理部は、
前記光学マルチセンサデータ、前記衛星データ、及び新たに算出された前記幾何モデルパラメータとから、新たに前記複数の帯域の衛星画像データを作成することを特徴とする。

この発明の位置誤差補正装置により、補正処理から幾何モデルパラメータ調整までの処理を再帰的に行うことにより、バンド間の位置ずれの補正を自動処理し、効率的に行うことができる。

実施の形態１．
図１〜図５を参照して実施の形態１を説明する。図１は、コンピュータで実現される衛星画像評価装置１００（位置誤差補正装置）の外観の一例を示す図である。図１において、衛星画像評価装置１００は、システムユニット８３０、ＣＲＴ（Ｃａｔｈｏｄｅ・Ｒａｙ・Ｔｕｂｅ）やＬＣＤ（液晶）の表示画面を有する表示装置８１３、キーボード８１４（Ｋｅｙ・Ｂｏａｒｄ：Ｋ／Ｂ）、マウス８１５、ＦＤＤ８１７（Ｆｌｅｘｉｂｌｅ・Ｄｉｓｋ・ Ｄｒｉｖｅ）、コンパクトディスク装置８１８（ＣＤＤ：Ｃｏｍｐａｃｔ Ｄｉｓｋ Ｄｒｉｖｅ）、プリンタ装置８１９などのハードウェア資源を備え、これらはケーブルや信号線で接続されている。

図２は、コンピュータで実現される衛星画像評価装置１００のハードウェア資源の一例を示す図である。図２において、衛星画像評価装置１００は、プログラムを実行するＣＰＵ８１０（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ：中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、プロセッサともいう）を備えている。ＣＰＵ８１０は、バス８２５を介してＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）８１１、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）８１２、表示装置８１３、キーボード８１４、マウス８１５、通信ボード８１６、ＦＤＤ８１７、ＣＤＤ８１８、プリンタ装置８１９、磁気ディスク装置８２０と接続され、これらのハードウェアデバイスを制御する。磁気ディスク装置８２０の代わりに、光ディスク装置、フラッシュメモリなどの記憶装置でもよい。

ＲＡＭ８１２は、揮発性メモリの一例である。ＲＯＭ８１１、ＦＤＤ８１７、ＣＤＤ８１８、磁気ディスク装置８２０等の記憶媒体は、不揮発性メモリの一例である。これらは、記憶装置あるいは記憶部、格納部、保存部、バッファの一例である。通信ボード８１６、キーボード８１４、ＦＤＤ８１７、ＣＤＤ８１８などは、取得部、入力部、入力装置の一例である。また、通信ボード８１６、表示装置８１３、プリンタ装置８１９などは、出力部、出力装置の一例である。

通信ボード８１６は、ネットワーク（例えばＬＡＮ：Ｌｏｃａｌ Ａｒｅａ Ｎｅｔｗｏｒｋ）に接続されている。通信ボード８１６は、ＬＡＮに限らず、インターネット、ＩＳＤＮ等のＷＡＮ（ワイドエリアネットワーク）などに接続されていても構わない。

磁気ディスク装置８２０には、オペレーティングシステム８２１（ＯＳ）、ウィンドウシステム８２２、プログラム群８２３、ファイル群８２４が記憶されている。プログラム群８２３のプログラムは、ＣＰＵ８１０、オペレーティングシステム８２１、ウィンドウシステム８２２により実行される。

上記プログラム群８２３には、以下に述べる実施の形態の説明において「〜部」として説明する機能を実行するプログラムが記憶されている。プログラムは、ＣＰＵ８１０により読み出され実行される。

ファイル群８２４には、以下に述べる実施の形態の説明において、「バンド間位置誤差」、「評価データ」などとして説明する情報や、「〜の判定結果」、「〜の算出結果」、「〜の抽出結果」、「〜の生成結果」、「〜の処理結果」として説明する情報や、データや信号値や変数値やパラメータなどが、「〜ファイル」や「〜データベース」の各項目として記憶されている。「〜ファイル」や「〜データベース」は、ディスクやメモリなどの記録媒体に記憶される。ディスクやメモリなどの記憶媒体に記憶された情報やデータや信号値や変数値やパラメータは、読み書き回路を介してＣＰＵ８１０によりメインメモリやキャッシュメモリに読み出され、抽出・検索・参照・比較・演算・計算・処理・出力・印刷・表示などのＣＰＵの動作に用いられる。抽出・検索・参照・比較・演算・計算・処理・出力・印刷・表示のＣＰＵの動作の間、情報やデータや信号値や変数値やパラメータは、メインメモリやキャッシュメモリやバッファメモリに一時的に記憶される。

また、以下に述べる実施の形態の説明においては、データや信号値は、ＲＡＭ８１２のメモリ、ＦＤＤ８１７のフレキシブルディスク、ＣＤＤ８１８のコンパクトディスク、磁気ディスク装置８２０の磁気ディスク、その他光ディスク、ミニディスク、ＤＶＤ（Ｄｉｇｉｔａｌ・Ｖｅｒｓａｔｉｌｅ・Ｄｉｓｋ）等の記録媒体に記録される。また、データや信号は、バス８２５や信号線やケーブルその他の伝送媒体によりオンライン伝送される。

また、以下に述べる実施の形態の説明において、「〜部」として説明するものは、「手段」、「〜回路」、「〜機器」であってもよく、また、「〜ステップ」、「〜手順」、「〜処理」であってもよい。すなわち、「〜部」として説明するものは、ＲＯＭ８１１に記憶されたファームウェアで実現されていても構わない。或いは、ソフトウェアのみ、或いは、素子・デバイス・基板・配線などのハードウェアのみ、或いは、ソフトウェアとハードウェアとの組み合わせ、さらには、ファームウェアとの組み合わせで実施されても構わない。ファームウェアとソフトウェアは、プログラムとして、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、ＤＶＤ等の記録媒体に記憶される。プログラムはＣＰＵ８１０により読み出され、ＣＰＵ８１０により実行される。すなわち、プログラムは、以下に述べる「〜部」としてコンピュータを機能させるものである。あるいは、以下に述べる「〜部」の手順や方法をコンピュータに実行させるものである。

図３は、実施の形態１における光学マルチセンサ幾何モデル評価システム１の構成図である。光学マルチセンサ幾何モデル評価システム１は、衛星画像利用システム２と、衛星画像評価装置１００とからなる。衛星画像評価装置１００は、画像処理部１０と、バンド間位置誤差算出部２０と、バンド間位置誤差評価部３０と、幾何モデルパラメータ算出部４０とを備える。また、衛星画像評価装置１００は、ネットワークを介して、衛星画像利用システム２と接続している。

図４は、画像処理部１０（ステップＳ１１〜Ｓ１４）、バンド間位置誤差算出部２０（Ｓ２１〜Ｓ２３）の動作を示すフローチャートである。図５は、バンド間位置誤差評価部３０（Ｓ３１〜Ｓ３５）、幾何モデルパラメータ算出部４０（Ｓ４１〜Ｓ４４）の動作を示すフローチャートである。図４、図５を参照して動作を説明する。

（画像処理部１０の動作）
まず、画像処理部１０の動作を説明する。図３に示すように、画像処理部１０は、画像処理用データ取得部１１と、画像処理用データ保存部１２と、システム補正部１３と、画像データ出力部１４とを備える。

（画像処理用データ取得部１１）
画像処理用データ取得部１１は、バンド間位置誤差評価部３０からの処理指示を受信し、処理指示に記載された条件を満たす光学マルチセンサデータ１０１と、軌道・姿勢データ１０２（衛星データ）とを、衛星画像利用システム２から取得し、また、幾何モデルパラメータ１０３を、幾何モデルパラメータ算出部４０から取得する（Ｓ１１）。ここで「光学マルチセンサデータ１０１」とは、所定の帯域のフィルタを持つ撮像素子である複数の光学センサから構成される光学マルチセンサによって地球の領域を対象として撮像されたデータである。また「軌道・姿勢データ１０２」（衛星データ）とは、光学マルチセンサが搭載された人工衛星の軌道データと姿勢データとを含むデータである。また「幾何モデルパラメータ１０３」とは、軌道データと姿勢データとから地球に対する光学マルチセンサの向きの特定に使用するパラメータである。後述のように、これら光学マルチセンサデータと、軌道・姿勢データと、幾何モデルパラメータとに基づき、画像処理部１０によって、地球の同一領域における複数の帯域の衛星画像データが作成される。画像処理用データ取得部１１は、取得した光学マルチセンサデータ１０１、軌道・姿勢データ１０２、幾何モデルパラメータ１０３を、画像処理用データ保存部１２に保管する。そして、光学マルチセンサデータ１０１、軌道・姿勢データ１０２、幾何モデルパラメータ１０３は、画像処理用データ保存部１２によって、システム補正部１３に出力される（Ｓ１２）。

（システム補正部１３）
システム補正部１３は、光学マルチセンサデータ１０１と、軌道・姿勢データ１０２と、幾何モデルパラメータ１０３を読み込み、バンド毎に放射量歪み、及び幾何的歪みを補正した補正済み画像データ１１０を作成し、画像データ出力部１４に出力する（Ｓ１３）。

（画像データ出力部１４）
画像データ出力部１４は、各バンドの補正済み画像データ１１０ａ〜１１０ｄを、バンド毎にそれぞれファイルに出力する（Ｓ１４）。

補正済み画像データ１１０は、バンドａ〜バンドｄの４バンドからなる画像データであり、補正済み画像データ１１０ａは、バンドａの画像データに対応し、
補正済み画像データ１１０ｂは、バンドｂの画像データに対応し、
補正済み画像データ１１０ｃは、バンドｃの画像データに対応し、
補正済み画像データ１１０ｄは、バンドｄの画像データに対応する。
なお、位置誤差を評価するための基準となるバンドは、バンドａとする。

（バンド間位置誤差算出部２０の動作）
次に、バンド間位置誤差算出部２０の動作を説明する。図３に示すように、バンド間位置誤差算出部２０は、画像データ取得部２１と、画像相関処理部２２と、位置誤差算出部２３とを備える。

（画像データ取得部２１）
画像データ取得部２１は、画像処理部１０にて処理した補正済み画像データ１１０ａ〜１１０ｄを取り込み、基準画像を補正済み画像データ１１０ａ、参照画像を幾何補正済み画像データ１１０ｂ〜１１０ｄとして、それぞれ画像相関処理部２２に出力する（Ｓ２１）。

（画像相関処理部２２）
画像相関処理部２２は、入力された基準画像である補正済み画像データ１１０ａと、参照画像である補正済み画像データ１１０ｂとの間で、例えば、面積相関法により、相関係数を求め、相関係数分布を作成する（Ｓ２２）。

（位置誤差算出部２３）
位置誤差算出部２３は、相関係数分布から、最大となる相関係数値を求め、あらかじめ決められた閾値以上の場合、バンド間位置誤差を算出する（Ｓ２３）。バンド間位置誤差は、それぞれクロストラック（Ｃｒｏｓｓ Ｔｒａｃｋ）方向成分と、アロングトラック（Ａｌｏｎｇ Ｔｒａｃｋ）方向成分から構成される。ここで、「クロストラック方向」とは、観測衛星（人工衛星）の進行方向に対して垂直な方向であり、「アロングトラック方向」とは、観測衛星の進行方向に対して並行方向である。

このようにして、残りの組み合わせである、補正済み画像データ１１０ａと補正済み画像データ１１０ｃ、および補正済み画像データ１１０ａと補正済み画像データ１１０ｄについても同様に、バンド位置誤差を算出する。また、図３では４バンドの光学マルチセンサデータを使用した動作を示しているが、４バンドに限定するものではない。４バンド以外の複数のバンドであっても構わない。

（バンド間位置誤差評価部３０の動作）
図３に示すように、バンド間位置誤差評価部３０は、位置誤差取得部３１、位置誤差量変換部３２、評価データ保存部３３、位置誤差評価部３４を備える。

（位置誤差取得部３１）
位置誤差取得部３１は、バンド間位置誤差算出部２０からバンド間位置誤差を取り込み、位置誤差量変換部３２に出力する（Ｓ３１）。

（位置誤差量変換部３２）
位置誤差量変換部３２は、位置誤差をピクセル単位の大きさから、センサ画素単位の大きさに単位変換する。この変換を行うことにより、直下視の位置誤差とポインティング時の位置誤差とを同じ大きさとして扱うことが可能である。変換した位置誤差は、評価データ保存部３３に保存される（Ｓ３２）。

（評価データ保存部３３）
評価データ保存部３３は、位置誤差量変換部３２で作成された変換後の位置誤差、及び位置誤差評価部３４によって後述の評価結果（評価データ）が作成された場合に、この評価データを取り込み、保存する（Ｓ３３）。

（位置誤差評価部３４）
位置誤差評価部３４は、作業者３の指示により、評価データ保存部３３に保存された変換後の位置誤差を取得し、クロストラック方向成分、及びアロングトラック方向成分のそれぞれに対して回帰分析を行い、バンド間位置誤差が、予め設定されている評価基準を満たしているか判定し、前記評価基準を満たすバンド間位置誤差からなるデータである評価データを作成すると共に、評価結果を作業者３に通知する（Ｓ３４）。なお、「評価基準」は、あらかじめ作業者３が位置誤差評価部３４に設定しておく。

また、位置誤差評価部３４は、作業者３の要求により、画像処理部１０から保存された画像処理用データのデータ一覧を取得し、作業者３に通知する。作業者３は、通知されたデータ一覧から、位置誤差を評価するための画像処理用データを選択し、位置誤差評価部３４を介して画像処理部１０に処理指示を通知する（Ｓ３５）。

（幾何モデルパラメータ算出部４０の動作）
図３に示すように、幾何モデルパラメータ算出部４０は、評価データ取得部４１、位置誤差補正量算出部４２、幾何モデルパラメータ作成部４３、幾何モデルパラメータ保存部４４、幾何モデルパラメータ管理部４５を備えている。

（評価データ取得部４１）
評価データ取得部４１は、バンド間位置誤差評価部３０からの作成要求を受信し、あるいは、作業者３が指示した作成要求を受信し、評価データ保存部３３から、評価データを取得し、位置誤差補正量算出部４２に出力する（Ｓ４１）。

（位置誤差補正量算出部４２）
位置誤差補正量算出部４２は、入力された評価データから、クロストラック方向の位置誤差と、アロングトラック方向の位置誤差を抽出し、それぞれ位置誤差から幾何歪みの補正量を算出し、算出結果を幾何モデルパラメータ作成部４３に出力する（Ｓ４２）。

（幾何モデルパラメータ作成部４３）
幾何モデルパラメータ作成部４３は、算出した補正量から新たな幾何モデルパラメータ１０４を作成し、幾何モデルパラメータ保存部４４に出力する（Ｓ４３）。

（幾何モデルパラメータ管理部４５）
幾何モデルパラメータ管理部４５は、画像処理部１０からの要求により、幾何モデルパラメータ保存部４４から新たな幾何モデルパラメータ１０４を取得し、画像処理部１０に出力する（Ｓ４４）。また、作業者３の要求により、幾何モデルパラメータ保存部４４から新たな幾何モデルパラメータ１０４を取得し、作業者３に出力する。

画像処理部１０では、新たな幾何モデルパラメータ１０４と、先に用いた「光学マルチセンサデータ」及び軌道・姿勢データとから、新たに補正済み画像データを生成する。衛星画像評価装置１００は、以上の画像処理部１０〜幾何モデルパラメータ算出部４０の処理を位置誤差が目標精度に達するまで繰り返えす。この処理により、バンド間の位置ずれを自動的に補正することが可能となり、従来は手作業で補正していた作業者の補正作業の負担を軽減することができる。

以上の実施の形態１では、装置としての衛星画像評価装置１００を説明したが、衛星画像評価装置１００において、画像処理部１０、バンド間位置誤差算出部２０、バンド間位置誤差評価部３０、幾何モデルパラメータ算出部４０の一連の動作は互いに関連しており、これら一連の動作は、コンピュータに実施させるプログラムによる処理に置き換えることができる。また、この「プログラム」を、コンピュータ読み取り可能な記録媒体に記録させることで、プログラムを記録したコンピュータ読み取り可能な記録媒体の実施形態とすることができる。また、画像処理部１０、バンド間位置誤差算出部２０、バンド間位置誤差評価部３０、幾何モデルパラメータ算出部４０の一連の動作をステップととらえることで、衛星画像評価方法、あるいは位置誤差補正方法に置き換えることができる。

以上の実施の形態では、
光学マルチセンサデータ、軌道・姿勢データ、幾何モデルパラメータに基づき放射量歪み及び幾何的歪みを補正した補正済み画像データを作成し、各バンド（帯域）の補正済み画像データを出力する画像処理部１０と、
各バンド（帯域）の補正済み画像データを取り込み、各バンドの画像データをもとに求めた相関係数より位置誤差を算出するバンド間位置誤差算出部２０と、
位置誤差を取り込み、位置誤差が評価基準を満たしているか否かの判定を行い、画像処理部１０へ画像処理指示を出力行うとともに、幾何モデルパラメータ作成要求を出力するバンド間位置誤差評価部３０と、
バンド間位置誤差評価部３０より幾何モデルパラメータ作成要求を受信した場合、位置誤差に基づき幾何モデルパラメータ作成するとともに、幾何モデルパラメータを画像処理部１０に出力する幾何モデルパラメータ算出部４０と
を備え、
画像処理部１０、バンド間位置誤差算出部２０、バンド間位置誤差評価部３０、幾何モデルパラメータ算出部４０の処理を再帰的に実行する衛星画像評価装置１００を説明した。

実施の形態１における衛星画像評価装置１００の外観の一例を示す図。 実施の形態１における衛星画像評価装置１００のハードウェア資源の一例を示す図。 実施の形態１における衛星画像評価装置１００のブロック図。 実施の形態１における画像処理部１０、バンド間位置誤差算出部２０の動作を示すフローチャート。 実施の形態１におけるバンド間位置誤差評価部３０、幾何モデルパラメータ算出部４０の動作を示すフローチャート。

符号の説明

１ 光学マルチセンサ幾何モデル評価システム、２ 衛星画像利用システム、３ 作業者、１０ 画像処理部、１１ 画像処理用データ取得部、１２ 画像処理用データ保存部、１３ システム補正部、１４ 画像データ出力部、２０ バンド間位置誤差算出部、２１ 画像データ取得部、２２ 画像相関処理部、２３ 位置誤差算出部、３０ バンド間位置誤差評価部、３１ 位置誤差取得部、３２ 位置誤差量変換部、３３ 評価データ保存部、３４ 位置誤差評価部、４０ 幾何モデルパラメータ算出部、４１ 評価データ取得部、４２ 位置誤差補正量算出部、４３ 幾何モデルパラメータ作成部、４４ 幾何モデルパラメータ保存部、４５ 幾何モデルパラメータ管理部、１００ 衛星画像評価装置、１０１ 光学マルチセンサデータ、１０２ 軌道・姿勢データ、１０３ 幾何モデルパラメータ、１０４ 新たな幾何モデルパラメータ、１１０ａ，１１０ｂ，１１０ｃ，１１０ｄ 補正済み画像データ。

Claims (3)

1. 所定の帯域のフィルタを持つ撮像素子である複数の光学センサから構成される光学マルチセンサによって地球の領域を対象として撮像された光学マルチセンサデータと、前記光学マルチセンサが搭載された人工衛星の軌道データと姿勢データとを含む衛星データと、前記衛星データに含まれる前記軌道データと前記姿勢データとから地球に対する前記光学マルチセンサの向きの特定に使用するパラメータである幾何モデルパラメータとを取得し、前記光学マルチセンサデータと、前記衛星データと、前記幾何モデルパラメータとに基づき、地球の同一領域における複数の帯域の衛星画像データを作成する画像処理部と、
   前記画像処理部によって作成された前記複数の帯域の衛星画像データを取り込み、所定の帯域の衛星画像データを基準画像データとし、前記基準画像データ以外のそれぞれの帯域の前記衛星画像データを参照画像データとし、前記基準画像データと各参照画像データとの間の位置誤差をバンド間位置誤差として前記参照画像データごとに算出するバンド間位置誤差算出部と、
   前記バンド間位置誤差算出部によって算出されたそれぞれの前記バンド間位置誤差に基づき、前記幾何モデルパラメータに代替する新たな幾何モデルパラメータを算出する幾何モデルパラメータ算出部と
   を備え、
   前記画像処理部は、
   前記光学マルチセンサデータ、前記衛星データ、及び新たに算出された前記幾何モデルパラメータとから、新たに前記複数の帯域の衛星画像データを作成することを特徴とする位置誤差補正装置。
2. 前記位置誤差補正装置は、さらに、
   前記バンド間位置誤差算出部によって算出された前記参照画像データごとの前記バンド間位置誤差を取り込み、取り込まれた各バンド間位置誤差が所定の評価基準を満たすかどうかを判定し、前記評価基準を満たすバンド間位置誤差からなるデータである評価データを作成するバンド間位置誤差評価部を備え、
   前記幾何モデルパラメータ算出部は、
   前記バンド間位置誤差評価部によって作成された前記評価データに基づいて、前記幾何モデルパラメータに代替する新たな幾何モデルパラメータを算出することを特徴とする請求項２記載の位置誤差補正装置。
3. （１）所定の帯域のフィルタを持つ撮像素子である複数の光学センサから構成される光学マルチセンサによって地球の領域を対象として撮像された光学マルチセンサデータと、前記光学マルチセンサが搭載された人工衛星の軌道データと姿勢データとを含む衛星データと、前記衛星データに含まれる前記軌道データと前記姿勢データとから地球に対する前記光学マルチセンサの向きの特定に使用するパラメータである幾何モデルパラメータとを取得し、前記光学マルチセンサデータと、前記衛星データと、前記幾何モデルパラメータとに基づき、地球の同一領域における複数の帯域の衛星画像データを作成する処理、
   （２）作成された前記複数の帯域の衛星画像データを取り込み、所定の帯域の衛星画像データを基準画像データとし、前記基準画像データ以外のそれぞれの帯域の前記衛星画像データを参照画像データとし、前記基準画像データと各参照画像データとの間の位置誤差をバンド間位置誤差として前記参照画像データごとに算出する処理、
   （３）算出されたそれぞれの前記バンド間位置誤差に基づき、前記幾何モデルパラメータに代替する新たな幾何モデルパラメータを算出する処理、
   （４）前記光学マルチセンサデータ、前記衛星データ、及び新たに算出された前記幾何モデルパラメータとから、新たに前記複数の帯域の衛星画像データを作成する処理、
   をコンピュータに実行させることを特徴とする位置誤差補正プログラム。

Images