【手続補正書】
【提出日】平成16年9月22日(2004.9.22)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】特許請求の範囲
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
処理部は、観測対象エリア内の各画素に対応して、衛星による植生の季節変動についての時系列の第１の観測データを記憶した第１の観測データファイルから、画素毎の第１の観測データを読み込むステップと、
処理部は、各画素毎に、ある時期又は時点に対応して、該時期又は時点の前後所定期間の局所最大値を抽出することにより時系列の第２の観測データを求め、それを各画素に対応して第２の観測データファイルに記憶するステップと、
処理部は、第２の観測データファイルに記憶された第２の観測データを、長期的非周期項及び季節変化項を含む次式の関数を用い、各第１の係数ｃｉを最小二乗法又は他の推定法により決定することでフィッティングし、各第１の係数ｃｉを第１のパラメータファイルに記憶するステップと、
【数１】
2003296702.tif 000001
［  ここで、
ｃｉ：求める第１の係数
ｔ：時期又は時間
Ｎ：周期関数のペア数（この場合ｓｉｎとｃｏｓのペア数）
ｌ：周期関数の仮の番号
ｋｌ：１２ヶ月が回帰の整数倍になる数列（例：１，２，３，４，６，１２）
Ｍ：１年間のデータ数（シーン数）  ］
処理部は、第１の係数ｃｉで定められた前記関数に従い、第２の観測データから雑音データを除去した第１の推定観測データを求め、それを第１の推定観測データファイルに記憶する、又は、出力部に出力若しくは表示するステップと
を含む地球観測衛星データのノイズ除去方法。
【請求項２】
処理部は、第２の観測データについて次式の時系列の状態空間モデルを適用し、該状態空間モデルの各第２の係数ｃｉ（ｔ）を推定する処理を実行し、各第２の係数ｃｉ（ｔ）を第２のパラメータファイルに記憶するステップと、
【数２】
2003296702.tif 000002
［  ここで、
ｃｉ（ｔ）：求める第２の係数
ｗｋ（ｔ）：システム雑音
ｙ（ｔ）：観測される値
ｖ（ｔ）：平均０、分散がσｖの白色雑音（観測雑音）
【数３】
2003296702.tif 000003
]
処理部は、第２の係数ｃｉ（ｔ）で定められた前記関数に従い、第２の観測データから雑音データを除去した第２の推定観測データを求め、それを第２の推定観測データファイルに記憶する、又は出力部に出力、若しくは表示するステップ
をさらに含む請求項１に記載の地球観測衛星データのノイズ除去方法。
【請求項３】
時刻ｋに対応して、時期ｋ−ｎから時期ｋまでのデータについて最大値を抽出するステップと、
時刻ｋに対応して、時期ｋから時期ｋ＋ｎのデータについて最大値を抽出するステップと、
時刻ｋに対応して、得られたこれら二つの最大値のうち、値の小さいものを局所的最大値とするステップと、
を含む請求項１又は２に記載の地球観測衛星データのノイズ除去方法。
【請求項４】
時期ｋに、推定値と実際の画素値のずれが、設定したしきい値よりも大きい場合には雲やシステムノイズの影響とみなしてその時期の観測データは採用しないようにして、採用されたデータだけを用いて再度、係数及び／又は推定観測データを求めることを特徴とする請求項１又は２に記載の地球観測衛星データのノイズ除去方法。
【請求項５】
前記第１の観測データは、植生指数データ、水分指数、中間赤外線バンド、熱赤外線輝度バンド、温度データのいずれか又は複数であることを特徴とする請求項１又は２に記載の地球観測衛星データのノイズ除去方法。
【請求項６】
前記雑音データは、雲、霧、雨、ヘイズ等の大気状態に起因するもの、及び／又は、センサ感度等のシステムに起因するものを含むことを特徴とする請求項１又は２に記載の地球観測衛星データのノイズ除去方法。
【請求項７】
処理部は、観測データおよび推定データの複数を組み合わせて出力部に出力又は表示するステップをさらに含む請求項１又は２に記載の地球観測衛星データのノイズ除去方法。
【請求項８】
処理部は、観測対象エリア内の各画素に対応して、衛星による植生の季節変動についての時系列の第１の観測データを記憶した第１の観測データファイルから、画素毎の第１の観測データを読み込むステップと、
処理部は、各画素毎に、ある時期又は時点に対応して、該時期又は時点の前後所定期間の局所最大値を抽出することにより時系列の第２の観測データを求め、それを各画素に対応して第２の観測データファイルに記憶するステップと、
処理部は、第２の観測データファイルに記憶された第２の観測データを、長期的非周期項及び季節変化項を含む次式の関数を用い、各第１の係数ｃｉを最小二乗法又は他の推定法により決定することでフィッティングし、各第１の係数ｃｉを第１のパラメータファイルに記憶するステップと、
【数４】
2003296702.tif 000004
［  ここで、
ｃｉ：求める第１の係数
ｔ：時期又は時間
Ｎ：周期関数のペア数（この場合ｓｉｎとｃｏｓのペア数）：
ｌ：周期関数の仮の番号
ｋｌ：１２ヶ月が回帰の整数倍になる数列（例：１，２，３，４，６，１２）
Ｍ：１年間のデータ数（シーン数）  ］
処理部は、第１の係数ｃｉで定められた前記関数に従い、第２の観測データから雑音データを除去した第１の推定観測データを求め、それを第１の推定観測データファイルに記憶する、又は出力部に出力若しくは表示するステップ
をコンピュータに実行させるための地球観測衛星データのノイズ除去処理プログラム。
【請求項９】
処理部は、第２の観測データについて次式の時系列の状態空間モデルを適用し、該状態空間モデルの各第２の係数ｃｉ（ｔ）を推定する処理を実行し、各第２の係数ｃｉ（ｔ）を第２のパラメータファイルに記憶するステップと、
【数５】
2003296702.tif 000005
［  ここで、
ｃｉ（ｔ）：求める第２の係数
ｗｋ（ｔ）：システム雑音
ｙ（ｔ）：観測される値
ｖ（ｔ）：平均０、分散がσｖの白色雑音（観測雑音）
【数６】
2003296702.tif 000006
］
処理部は、第２の係数ｃｉ（ｔ）で定められた前記関数に従い、第２の観測データから雑音データを除去した第２の推定観測データを求め、それを第２の推定観測データファイルに記憶する、又は出力部に出力若しくは表示するステップ
をさらにコンピュータに実行させるための請求項８に記載の地球観測衛星データのノイズ除去処理プログラム。
【請求項１０】
処理部は、観測対象エリア内の各画素に対応して、衛星による植生の季節変動についての時系列の第１の観測データを記憶した第１の観測データファイルから、画素毎の第１の観測データを読み込むステップと、
処理部は、各画素毎に、ある時期又は時点に対応して、該時期又は時点の前後所定期間の局所最大値を抽出することにより時系列の第２の観測データを求め、それを各画素に対応して第２の観測データファイルに記憶するステップと、
処理部は、第２の観測データファイルに記憶された第２の観測データを、長期的非周期項及び季節変化項を含む次式の関数を用い、各第１の係数ｃｉを最小二乗法又は他の推定法により決定することでフィッティングし、各第１の係数ｃｉを第１のパラメータファイルに記憶するステップと、
【数７】
2003296702.tif 000007
［  ここで、
ｃｉ：求める第１の係数
ｔ：時期又は時間
Ｎ：周期関数のペア数（この場合ｓｉｎとｃｏｓのペア数）：
ｌ：周期関数の仮の番号
ｋｌ：１２ヶ月が回帰の整数倍になる数列（例：１，２，３，４，６，１２）
Ｍ：１年間のデータ数（シーン数）  ］
処理部は、第１の係数ｃｉで定められた前記関数に従い、第２の観測データから雑音データを除去した第１の推定観測データを求め、それを第１の推定観測データファイルに記憶する、又は出力部に出力、若しくは表示するステップ
をコンピュータに実行させるための地球観測衛星データのノイズ除去処理プログラムを記録したコンピュータ読み取り可能な記録媒体。
【請求項１１】
処理部は、第２の観測データについて次式の時系列の状態空間モデルを適用し、該状態空間モデルの各第２の係数ｃｉ（ｔ）を推定する処理を実行し、各第２の係数ｃｉ（ｔ）を第２のパラメータファイルに記憶するステップと、
【数８】
2003296702.tif 000008
［  ここで、
ｃｉ（ｔ）：求める第２の係数
Ｗｋ（ｔ）：システム雑音
ｙ（ｔ）：観測される値
ｖ（ｔ）：平均０、分散がσｖの白色雑音（観測雑音）
【数９】
2003296702.tif 000009
］
処理部は、第２の係数ｃｉ（ｔ）で定められた前記関数に従い、第２の観測データから雑音データを除去した第２の推定観測データを求め、それを第２の推定観測データファイルに記憶する、又は出力部に出力、若しくは表示するステップ
をさらにコンピュータに実行させるための地球観測衛星データのノイズ除去処理プログラムを記録した請求項１０に記載のコンピュータ読み取り可能な記録媒体。