(57)【要約】【課題】  少ない計算量で精度の高いレーダ画像からの波浪方向スペクトル逆推定方法及びシステムを提供すること。【解決手段】  レーダで電波を海面に照射し、その後方散乱の強度を処理して得られた画像から、波浪の方向スペクトルを逆推定する方法において、波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用し、レーダの方位と周波数の座標で展開し、ベイズモデルを導入して後方散乱波の強度画像から波浪の方向スペクトルを求める。積分方程式に対して、対数をとり、波浪の波向を半円方向に限定し、積分を離散化することで、積分方程式を非線形代数方程式に近似してもよい。【選択図】  図４
2005043333.tif 000002
【特許請求の範囲】
【請求項１】
レーダで電波を海面に照射し、その後方散乱の強度を処理して得られた画像から、波浪の方向スペクトルを逆推定する方法であって、
波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、
これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用し、
レーダの方位と周波数の座標で展開し、ベイズモデルを導入して後方散乱波の強度画像から波浪の方向スペクトルを求める
ことを特徴とするレーダ画像からの波浪方向スペクトル逆推定方法。
【請求項２】
波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用した結果得られる積分方程式に対して、
対数をとり、波浪の波向を半円方向に限定し、積分を離散化することで、
積分方程式を非線形代数方程式に近似する
請求項１に記載のレーダ画像からの波浪方向スペクトル逆推定方法。
【請求項３】
滑らかな分布をする波浪の方向スペクトルの先験条件を、ベイズモデルの事前分布に組み込み、
その先験条件の重み係数を、ＡＢＩＣ（赤池のベイズ型情報基準）の最小化によって決定し、
ベイズの最尤法によって、波浪の方向スペクトルを逐次近似で求める
請求項１または２に記載のレーダ画像からの波浪方向スペクトル逆推定方法。
【請求項４】
レーダで電波を海面に照射し、その後方散乱の強度を処理して得られた画像から、波浪の方向スペクトルを逆推定するシステムであって、
波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、
これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用し、
レーダの方位と周波数の座標で展開し、ベイズモデルを導入して後方散乱波の強度画像から波浪の方向スペクトルを求める演算部を備える
ことを特徴とするレーダ画像からの波浪方向スペクトル逆推定システム。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、航空機や人工衛星等の飛翔体に搭載された映像化レーダ（合成開口レーダ、実開口レーダ等）で観測された海表面からの後方散乱波の強度を表す画像から、波浪の方向スペクトルを逆推定する方法とシステムに関するものである。
【背景技術】
【０００２】  航空機や人工衛星に搭載されている合成開口レーダ（ＳＡＲ：Ｓｙｎｔｈｅｔｉｃ  Ａｐｅｒｔｕｒｅ  Ｒａｄａｒ）は、マイクロ波を地表面や海面に照射して、そこからの反射や後方散乱を受信する構成になっている。
マイクロ波の照射方向は、通常、航空機や人工衛星の進行方向に対して横向きでかつ斜め下向きである。このため、進行方向に対する反射・後方散乱の左右の曖昧さを避けることができる。マイクロ波を直下に照射するセンサではこの左右の曖昧さのため、左右どちらからの反射・後方散乱かを特定することはできない。
また、地表面及び海面の映像化は、一般的なレーダのようにアンテナを動かして掃引するのではなく、鉛直方向に幅広く水平方向に幅狭いビームを使用し、アンテナを固定したままで航空機や人工衛星の移動とともに地表面や海面を帯状に掃引する。【０００３】  ＳＡＲは、太陽光を必要とせずに地表面や海面を観測することができる能動型センサである。このため、季節や昼夜に関係なく観測を行うことができる。また、ＳＡＲで用いるマイクロ波は、大気における透過率が高いので、雲による影響を受けることなく定期的な観測を行うことができる。これにより、従来の可視・近赤外のセンサでは観測することができなかった台風や低気圧下における海面の様子を観測することができる。【０００４】  可視・近赤外センサは、観測高度が高くなると地表面の空間分解能は低下する。通常、この種のセンサは観測高度が２倍になると空間分解能は４分の１になる。
ＳＡＲの空間分解能は、観測に用いているマイクロ波パルスの長さ、パルスの帯域幅、アンテナの大きさなどから決まるため，観測高度や用いているマイクロ波の波長には無関係である。そのため、人工衛星の観測高度からも高い空間分解能で観測することができる。また、ＳＡＲは、パルス圧縮や合成開口処理によって、他のマイクロ波センサと比較して極めて高い空間分解能を有している。【０００５】  ＳＡＲは、マイクロ波を斜め下方向に照射して、海面からの反射・後方散乱波の強度を測定している。マイクロ波の反射・散乱現象は海面の粗さと傾斜による影響を強く受ける。特に、ＳＡＲで海面を観測すると、海面に存在する波浪によってマイクロ波の入射角の大きさが連続的に変化する様子を観測することができる。このマイクロ波の入射角の変化を解析することにより、波浪の情報（波高、波向、周期、方向スペクトル）を得ることができる。【０００６】  ＳＡＲ画像をスペクトル解析することによって得られる２次元波数スペクトルと、波浪の方向スペクトルとの関係については、これまで様々な研究が行われ、両者を関係づける定式化がなされてきた。例えば、Ｈａｓｅｌｌｍａｎ（１９９１）は、準線形理論に基づいて合成開口処理の際に生じる非線形の効果を考慮して定式化を行っている。
ＳＡＲ画像の２次元波数スペクトルから波浪の方向スペクトルを求めるためには、これらの関係式の逆解析を行わなければならない。これまで、この逆問題を解析する方法には以下のようなものがある。【０００７】
【非特許文献１】  Ｈａｓｓｅｌｍａｎｎ，Ｋ．ａｎｄ  Ｓ．Ｈａｓｓｅｌｍａｎｎ，“Ｏｎ  ｔｈｅ  Ｎｏｎｌｉｎｅａｒ  Ｍａｐｐｉｎｇ  ｏｆ  ａｎ  Ｏｃｅａｎ  Ｗａｖｅ  Ｓｐｅｃｔｒａ  ｉｎｔｏ  ａ  Ｓｙｎｔｈｅｔｉｃ  Ａｐｔｕｒｅ  Ｒａｄａｒ  Ｉｍａｇｅ  Ｓｐｅｃｔｒｕｍ  ａｎｄ  ｉｔｓ  Ｉｎｖｅｒｓｉｏｎ”，Ｊ．Ｇｅｏｐｈｙｓ．Ｒｅｓ．，Ｖｏｌ．９６，Ｎｏ．Ｃ６，１９９１，ｐｐ．１０７１３−１０７９９．
【非特許文献２】  Ｅｎｇｅｎ，Ｇ．，Ｈ．Ｊｏｈｎｓｅｎ，Ｈ．Ｅ．Ｋｒｏｇｓｔａｄ，Ｓ．Ｆ．Ｂａｒｓｔｏｗ，”Ｄｉｒｅｃｔｉｏｎａｌ  ｗａｖｅ  ｓｐｅｃｔｒａ  ｂｙ  ｉｎｖｅｒｓｉｏｎ  ｏｆ  ＥＲＳ−１  ｓｙｎｔｈｅｔｉｃ  ａｐｅｒｔｕｒｅ  ｒａｄａｒ  ｏｃｅａｎ  ｉｍａｇｅｒｙ”，ＩＥＥＥ  Ｔｒａｎｓ．  Ｏｎ  Ｇｅｏｓｃｉｅｎｃｅ  ａｎｄ  Ｒｅｍｏｔｅ  ｓｅｎｓｉｎｇ，Ｖｏｌ．３２，Ｎｏ．２，１９９４，ｐｐ．３４０−３５２
【非特許文献３】  泉宮尊司，居場博之，”合成開口レーダによる海面画像の順線形理論に基づいた逆解析解の存在とその特性”，海岸工学論文集，第４７巻，２０００，ｐｐ．１３２６−１３３０
【０００８】  非特許文献１の逆解析法は、波浪推算で与えた初期条件に依存するような解析法なので、初期条件として与える波浪推算が存在しない場合には解析することができない。いわゆる初期値に依存した解析法（初期値問題がある解析法）になっている。
非特許文献２及び３の逆解析法は、解を一意に決めるために人為的な処理を行っていて、完全に逆問題を解いた形になっていないために、必ずしも精度良く波浪の方向スペクトルを推定できない。また、ＳＡＲ画像中に存在するノイズの大きさに依存して、その精度が大きく変化してしまう。
【発明の開示】
【発明が解決しようとする課題】
【０００９】  そこで、本発明は、合成開口レーダに代表されるレーダ画像から、統計学と情報学に基づいて、初期値や画像中のノイズに依存することなく高精度で、波浪の方向スペクトルを逆推定する方法及びシステムを提供することを課題とする。
【課題を解決するための手段】
【００１０】  上記課題を解決するために、本発明のレーダ画像からの波浪方向スペクトル逆推定方法は、次の構成を備える。
すなわち、レーダで電波を海面に照射し、その後方散乱の強度を処理して得られた画像から、波浪の方向スペクトルを逆推定する方法において、波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用し、レーダの方位と周波数の座標で展開し、ベイズモデルを導入して後方散乱波の強度画像から波浪の方向スペクトルを求めることを特徴とする。【００１１】  ここで、波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用した結果得られる積分方程式に対して、対数をとり、波浪の波向を半円方向に限定し、積分を離散化することで、積分方程式を非線形代数方程式に近似して、画像の２次元波数スペクトルの変動を滑らかにし、波浪方向スペクトルの推定精度の向上に寄与させてもよい。【００１２】  滑らかな分布をする波浪の方向スペクトルの先験条件を、ベイズモデルの事前分布に組み込み、その先験条件の重み係数を、ＡＢＩＣ（赤池のベイズ型情報基準）の最小化によって決定し、ベイズの最尤法によって、波浪の方向スペクトルを逐次近似で求めて、過剰な先験条件と先験条件に対する重み係数の客観的な決定に対応させてもよい。【００１３】  また、本発明のレーダ画像からの波浪方向スペクトル逆推定システムは、レーダで電波を海面に照射し、その後方散乱の強度を処理して得られた画像から、波浪の方向スペクトルを逆推定するシステムであって、波浪の方向スペクトルを、指数部に離散的一定値関数を有する指数関数で表し、これを、準線形理論に基づく波浪の方向スペクトルとレーダ画像の２次元波数スペクトルとの関係式に適用し、レーダの方位と周波数の座標で展開し、ベイズモデルを導入して後方散乱波の強度画像から波浪の方向スペクトルを求める演算部を備えることを特徴とする。
【発明の効果】
【００１４】  本発明によると、少ない計算量で精度高く、レーダ画像から波浪方向スペクトルを逆推定することができる。特に、実際の観測で画像中に生じるスペックルノイズの影響を最小限に抑えて、推定精度を向上させることができる。
【発明を実施するための最良の形態】
【００１５】  ＳＡＲ等のレーダ画像をスペクトル解析することによって得られる２次元波数スペクトルと、波浪の方向スペクトルとの関係については、これまで様々な研究が行われ、両者を関係づける定式化がなされている。
本発明では、Ｈａｓｅｌｌｍａｎ（１９９１）が提案した準線形理論に基づいて定式化された以下の式を用いて、２次元波数スペクトルから波浪の方向スペクトルを逆推定する。【００１６】
【数１】
2005043333.tif 000003
【００１７】  ここで、ｋ＝（ｋａ，ｋｒ）は波浪の波数ベクトル、ｋ（＝｜ｋ｜）は波数、ｋａはアジマス方向の波数、ｋｒはレンジ方向の波数、Ｒは人工衛星と観測地点との距離、Ｖは人工衛星の移動速度、ＳＳＡＲ（ｋ）はＳＡＲ画像の２次元波数スペクトル、Ｓｗ（ｋ）は波浪の方向スペクトル、Ｔｖ（ｋ）とＴｓ（ｋ）はそれぞれ変調関数を表し、以下のように定式化されている。【００１８】
【数２】
2005043333.tif 000004
【００１９】
【数３】
2005043333.tif 000005
【００２０】  ここで、ωは波浪の周波数、ｋｒａｄａｒはＳＡＲで使用しているマイクロ波の波数、θはマイクロ波の海面への入射角、Ｍｔｉｔｌ（ｋ）、Ｍｈｙｄｒｏ（ｋ）、Ｍｖｂ（ｋ）はそれぞれ傾き変調関数、水理学的変調関数、速度バンチング変調関数を表している。
傾き変調関数Ｍｔｉｔｌ（ｋ）、水理学的変調関数Ｍｈｙｄｒｏ（ｋ）、速度バンチング変調関数Ｍｖｂ（ｋ）はそれぞれ以下に示すように定式化されている。【００２１】
【数４】
2005043333.tif 000006
【００２２】
【数５】
2005043333.tif 000007
【００２３】
【数６】
2005043333.tif 000008
【００２４】
2005043333.tif 000009
する。これにより、方向スペクトルＳ（ｆ，θ）は、指数部に離散的一定値関数を有する指数関数を用いて以下のように定義する。【００２５】
【数７】
2005043333.tif 000010
【００２６】  ここで、Ｉは周波数分割数、Ｊは方向分割数で、δｉ，ｊ（ｆ，θ）はデルタ関数をそれぞれ表している。デルタ関数δｉ，ｊ（ｆ，θ）は以下のような関係を満足している。【００２７】
【数８】
2005043333.tif 000011
【００２８】  また、式７における周波数ｆと方向角θは、それぞれ次式によって離散化する。【００２９】
【数９】
2005043333.tif 000012
【００３０】  式７を式１に代入すれば、式１は未知数Ｘ＝（ｘ１，１，・・・，ｘＩ，Ｊ）ｔを含む積分方程式になる。
しかしながら、式１において指数部に積分が存在するため、容易には式１から波浪の方向スペクトルを逆推定することはできない。
そこで、本発明では、式１の両辺の対数を取り、以下のように式変形を行った。【００３１】
【数１０】
2005043333.tif 000013
【００３２】  上記のように変形することにより、観測より得られるＳＡＲなどのレーダ画像の２次元波数スペクトルの変動は滑らかになり、推定する方向スペクトルの推定精度が向上する。
特に、実際の観測では、画像中にはスペックルノイズと呼ばれる時空間での変動が激しいノイズ成分が混入し、画像の２次元波数スペクトルの変動を大きくしている。そのため、この大きな変動ために推定される方向スペクトルの精度が著しく低下する。本発明では、対数を取ることにより、このノイズの影響を最小限にし、推定精度を向上させた。【００３３】  画像から特定される波浪の波向には１８０度のあいまいさが存在する。
この１８０度のあいまいさを表現するために式１０の右辺第２項は、｛｜ＴＳ（ｋ）｜２ＳＷ（ｋ）＋｜ＴＳ（−ｋ）｜２ＳＷ（−ｋ）｝の形を取っている。
本発明では、推定する方向スペクトルにおいて、波浪の波向を０〜１８０度に限定することにより式１０の第２項中の｜ＴＳ（−ｋ）｜２ＳＷ（−ｋ）を省略することができるようにした。つまり、式１０は以下のように変形することができる。【００３４】
【数１１】
2005043333.tif 000014
【００３５】  式１１の積分を離散化することにより、積分方程式は非線形代数方程式で近似することができる。
本発明では、式１１の離散化した積分を計算するために、求める波浪の方向スペクトル値をその４近傍の格子点上の方向スペクトルを用いて、以下のようにして内挿することにより算定する。【００３６】
【数１２】
2005043333.tif 000015
【００３７】
2005043333.tif 000016
きるようにした。
さらに、画像の２次元波数スペクトルの誤差Ｅεｋ４を考慮し、最終的には、式１の積分方程式は、次式で表される未知数Ｘ＝（ｘ１，１，・・・，ｘＩ，Ｊ）ｔを含む非線形代数方程式で近似される。【００３８】
【数１３】
2005043333.tif 000017
【００３９】  ここで、Ｆｋ（Ｘ）は波数ごとに積分方程式を離散化したベクトル関数、サフィック
2005043333.tif 000018
は、次式で与えられる。【００４０】
【数１４】
2005043333.tif 000019
【００４１】  波浪の方向スペクトルは離散的一定値関数として近似したが、各区間のエネルギーの相関については考慮していない。波の線形性からは各微小区間毎のエネルギー分布はそれぞれ独立であると見なされるが、方向スペクトルが不連続なエネルギー分布をしているとは考えにくい。また、一般に方向スペクトルは滑らかな連続関数であると見なされている。
したがって、ここでは方向スペクトルＳ（ｆ，θ）が滑らかであるという仮定の表現として、８方向の微分オペレータを用いて以下のように先験条件を規定する。【００４２】
【数１５】
2005043333.tif 000020
【００４３】  また，解析対象周波数の上限（ｉ＝Ｉ）及び下限（ｉ＝１）では、５方向の微分オペレータを用いて以下のように先験条件を規定する。【００４４】
【数１６】
2005043333.tif 000021
【００４５】  以上、式１５と式１６を用いると先験条件は、次のように表せる。【００４６】
【数１７】
2005043333.tif 000022
【００４７】  ただし、ｘｉ，０＝ｘｉ，Ｊ、ｘｉ，−１＝ｘｉ，Ｊ−１とする。式１３は作用行列Ｄを導入することにより、次のように行列表現できる。【００４８】
【数１８】
2005043333.tif 000023
【００４９】  式１８は、その値が小さいほど、方向スペクトルの推定値が滑らかになることを表している。
したがって、方向スペクトルＳ（ｆ，θ）の推定値としては、式１８をあまり大きくしない範囲で尤度（式１４）の大きいものが望ましい。
これを定式化すると、適当なパラメータｕ２（超パラメータ、重み係数）を用いて、【００５０】
【数１９】
2005043333.tif 000024
【００５１】を最大化するＸ＝（ｘ１，１，・・・，ｘＩ，Ｊ）ｔを求めればよい。
これはベイズの推論方法において、事後分布ｐＰＯＳＴ（  ）の【００５２】
【数２０】
2005043333.tif 000025
【００５３】の関係式において、Ｘ＝（ｘ１，１，・・・，ｘＩ，Ｊ）ｔの事前分布として、次式を想定したことに相当する。【００５４】
【数２１】
2005043333.tif 000026
【００５５】  ｕ２を与えれば、式２１を最大化するＸは、λ２に無関係に決まり、【００５６】
【数２２】
2005043333.tif 000027
を最小化することにより得られる。【００５７】  つまり、画像の２次元波数スペクトルから方向スペクトルを逆推定する逆問題を、式２２で表される最小値検索問題に帰着したことに相当する。
ここで、式２２の第１項は客観的な誤差エネルギーを表し、第２項は主観的な誤差エネルギーをそれぞれ表している。
客観的な誤差エネルギーは、観測値と推定値の一致度を表し、この値が小さいほど観測値と推定値が良く一致していることを表している。これに対して、主観的な誤差エネルギーは、推定された方向スペクトルが先験条件をどの程度満たしているのかを表している。【００５８】  超パラメータｕは主観的な誤差エネルギーの重みを与える定数で、その値が小さくなれば推定される方向スペクトルは画像の２次元波数スペクトルに一致するように推定され、その値が０のときは事前確率分布（先験条件）を考慮しない最小２乗法によって推定される方向スペクトルと同じ結果を得る。これに対して、超パラメータｕの値が大きくなると、推定される方向スペクトルは主観（先験条件）にかなったものになり、画像の２次元波数スペクトルの値が推定に反映されなくなる。
超パラメータｕとλ２は、次式で表されるＡＢＩＣ（赤池のベイズ型情報量基準、赤池（１９８０））を最小化することで、解の確からしさと滑らかさの両方の観点から望ましいｕが決められる。【００５９】
【数２３】
2005043333.tif 000028
【００６０】  方向スペクトルを推定する場合、式１７の最小化及び式１８の積分と最小化を実行しなければならない。しかしながら、いまの場合それらを解析的に行うことは不可能である。
そこで、テイラーの微分補正法を用いて式２２を線形化し、繰り返し計算によ
2005043333.tif 000029
てＸ０のまわりでＦｋ（Ｘ）をＴａｙｌｏｒ展開すると、【００６１】
【数２４】
2005043333.tif 000030
【００６２】
【数２５】
2005043333.tif 000031
【００６３】となる。
式２４を用いて式１３を行列表示すると、次のようになる。【００６４】
【数２６】
2005043333.tif 000032
【００６５】
【数２７】
2005043333.tif 000033
【００６６】
2005043333.tif 000034
まず、適当なｕ２及びＸの初期値Ｘ０を与えて、式２８で示される最小２乗法によ
2005043333.tif 000035
【００６７】
【数２８】
2005043333.tif 000036
【００６８】  ここで、最初にＸ０を与えて式２８に最小２乗法を適用し、第１近似解Ｘ（１）を算出する。次に式２６のＸ０をＸ（１）に置き換えて、式２８に最小２乗法を適用し、Ｘ（２）を
2005043333.tif 000037
【００６９】  得られた結果を用いて、次式により、与えられたｕ２に対するＡＢＩＣを算出する。【００７０】
【数２９】
2005043333.tif 000038
【００７１】
【数３０】
2005043333.tif 000039
【００７２】次いで、ｕ２を種々に変えて上記の計算を繰り返す。
2005043333.tif 000040
【００７３】
2005043333.tif 000041
Ｓ（ｆ，θ）が得られる。
なお、本発明では、適用する際の数値計算の利便性を考え、全ての計算ケースで初期値Ｘ０を０として計算し、数値計算の安定性を確認している。また、式２８の繰り返し計算の収束基準としては、（ｋ＋１）ステップ目における未知数の変化の大きさ‖Ｘ（ｋ＋１）−Ｘ（ｋ）‖と（ｋ）ステップ目における未知数の大きさ‖Ｘ（ｋ）‖との比が１０−２以下に
2005043333.tif 000042
また、超パラメータｕ２の設定に際しては、一般には下式を用い、格子探索法による繰り返し計算を行っている。【００７４】
【数３１】
2005043333.tif 000043
【実施例】
【００７５】  本発明の推定精度及び問題点を明らかにするために、数値シミュレーションによる検証を行った。
この検証では、モデル波浪の方向スペクトルからＳＡＲ画像の２次元波数スペクトルを算定し、本発明による波浪の方向スペクトル逆推定法によって波浪の方向スペクトルを逆推定した。推定結果と計算の初期条件で与えたモデル波浪の方向スペクトルとの比較を行うことにより、本発明で提案した推定法の精度及び問題点を検証した。
計算の初期条件で与えるモデル波浪の方向スペクトルＳＳ（ｆ，θ）のは以下に示すＢｒｅｔｓｃｈｎｅｉｄｅｒ−Ｍｉｔｓｕｙａｓｕ型周波数スペクトルＳｆ（ｆ）と光易型型方向分布関数ＧＳ（θ）のを掛け合わせたものを用いた。【００７６】
【数３２】
2005043333.tif 000044
【００７７】
【数３３】
2005043333.tif 000045
【００７８】
【数３４】
2005043333.tif 000046
【００７９】
【数３５】
2005043333.tif 000047
【００８０】  ここで、ｆは周波数、θは波向、Ｈ１／３は有義波高、Ｔ１／３は有義周期、θ０は主波向、
2005043333.tif 000048
をそれぞれ表している。【００８１】  図１に、モデル波浪の方向スペクトルの一例を示す。主波向０度、有義波高２［ｍ］、有義周期６［ｓ］、集中度パラメータ８０の場合の方向スペクトルを示している。
図１に示される波浪の方向スペクトルを式６に代入して計算したＳＡＲ画像の２次元波数スペクトルを図２に示す。図２は、ＳＡＲによって海表面を観測して得られる画像を２次元のスペクトル解析をすることによって得られる。
図３にＳＡＲ画像の２次元波数スペクトル（図２）から、波浪の方向スペクトルを逆推定した結果を示す。
図１と図３を比較すると、波浪の方向スペクトルの形状が非常に類似していることがわかる。また、両者の有義波高と有義周期がほぼ一致し、ＳＡＲ画像の２次元波数スペクトルから波浪の方向スペクトルを精度良く推定できていることがわかる。
図４〜７に主波向を変化（９０度、１８０度、２２５度、２７０度）させたときの逆推定結果を示す。これらの図より、波向によらず、ＳＡＲ画像の２次元波数スペクトルから波浪の方向スペクトルを精度良く逆推定することができることがわかる。
【産業上の利用可能性】
【００８２】  本発明によると、実際の観測で画像中に生じるスペックルノイズの影響を抑えて、レーダ画像から波浪方向スペクトルを実用レベルの高精度で逆推定することができる。これによって、気象海象予報や海難防止の他、海洋学や、海岸工学、水産業等の分野でも有用である。
【図面の簡単な説明】
【００８３】
【図１】  モデル波浪の方向スペクトルの例を示すグラフ
【図２】  ＳＡＲ画像の２次元波数スペクトルのグラフ
【図３】  ＳＡＲ画像２次元波数スペクトルから、波浪方向スペクトルを逆推定した結果のグラフ
【図４】  主波向を変化（９０度）させたときの真値とその逆推定結果のグラフ
【図５】  主波向を変化（１８０度）させたときの真値とその逆推定結果のグラフ
【図６】  主波向を変化（２２５度）させたときの真値とその逆推定結果のグラフ
【図７】  主波向を変化（２７０度）させたときの真値とその逆推定結果のグラフ
【図１】
2005043333.tif 000049
【図２】
2005043333.tif 000050
【図３】
2005043333.tif 000051
【図４】
2005043333.tif 000052
【図５】
2005043333.tif 000053
【図６】
2005043333.tif 000054
【図７】
2005043333.tif 000055