(57)【要約】【課題】水面波の偽像による影響を抑制する。【解決手段】本発明の観測装置は、時系列のレーダ反射波の波形データをフーリエ変換によって複数の周波数成分からなる成分波形データに変換するスペクトル解析部と、水面波の波形を表す波モデルデータを予め格納する波モデルデータベースと、波モデルデータベースが記憶する波モデルデータに基づき、現在の水面波と合致する波モデル成分波形データを決定する波モデル成分波形データ決定部と、スペクトル解析部が変換した成分波形データにおいて波モデル成分波形データ決定部が決定した波モデル成分波形データが含む周波数成分を抑制することで波成分抑制成分波形データを生成する波成分抑制部と、波成分抑制部が生成した波成分抑制成分波形データを逆フーリエ変換することでレーダによる水面に浮かぶ対象物の観測波を生成する観測波生成部とを備える。【選択図】図１
2018146540.tif 000002
【特許請求の範囲】
【請求項１】
レーダによって検出した時系列のレーダ反射波の波形データを取得する反射波取得部と、
前記時系列のレーダ反射波の波形データをフーリエ変換によって複数の周波数成分を含む成分波形データに変換するスペクトル解析部と、
予め定められた一つ又は複数の水面波の波形を表す波モデルデータの中から、現在の水面波と合致する波モデル成分波形データを決定する波モデル成分波形データ決定部と、
前記スペクトル解析部が変換した前記成分波形データにおいて、前記波モデル成分波形データ決定部が決定した前記波モデル成分波形データに含まれる周波数成分を抑制することで波成分抑制成分波形データを生成する波成分抑制部と、
前記波成分抑制部が生成した前記波成分抑制成分波形データを逆フーリエ変換することで前記レーダによる水面に浮かぶ対象物の観測波を生成する観測波生成部と、
を備える観測システム。
【請求項２】
前記波モデル成分波形データ決定部が、前記一つ又は複数の水面波の波形を表す波モデルデータに基づき、現在の風速または波高に応じて前記現在の水面波と合致する前記波モデル成分波形データを決定する、
請求項１に記載の観測システム。
【請求項３】
前記一つ又は複数の水面波の波形を表す波モデルデータには、波候の情報と対応付けられた前記波モデルデータが含まれ、
前記波モデル成分波形データ決定部が、前記一つ又は複数の水面波の波形を表す波モデルデータに基づき、現在の波候に応じて、前記現在の水面波と合致する前記波モデル成分波形データを決定する、
請求項１または請求項２に記載の観測システム。
【請求項４】
レーダと、
請求項１から請求項３の何れか１項に記載の観測システムと、
水面波の波形を表す波モデルデータを予め格納する波モデルデータベースと、
を備える水上観測システム。
【請求項５】
レーダによって検出した時系列のレーダ反射波の波形データを取得する反射波取得ステップと、
前記時系列のレーダ反射波の波形データをフーリエ変換によって複数の周波数成分を含む成分波形データに変換するスペクトル解析ステップと、
予め定められた一つ又は複数の水面波の波形を表す波モデルデータの中から、現在の水面波と合致する波モデル成分波形データを決定する波モデル成分波形データ決定ステップと、
前記スペクトル解析ステップにて変換した前記成分波形データにおいて、前記波モデル成分波形データ決定ステップにて決定した前記波モデル成分波形データに含まれる周波数成分を抑制することで波成分抑制成分波形データを生成する波成分抑制ステップと、
前記波成分抑制ステップにて生成した前記波成分抑制成分波形データを逆フーリエ変換することで前記レーダによる水面に浮かぶ対象物の観測波を生成する観測波生成ステップと、
を有する観測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、観測システム、水上観測システムおよび観測方法に関する。
【背景技術】
【０００２】
船に搭載したレーダから、海上に向けてレーダを発射し、その反射波によって海上に存在する他の船などを検知したり、その位置を判別したりする技術が存在する。また、陸上の観測施設から、同様に海上に向けてレーダを発射し、海上の船などを監視する技術が存在する。
【０００３】
特許文献１には、レーダ等の電波放射機器から放射される電波を受信し、それらの電波の電波諸元を利用することで電波放射源を分類識別する電波源の分類識別方法が記載されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平５−０４５３９０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところで、レーダによって検出した時系列のレーダ反射波に基づき水面に浮かぶ対象物の観測波を生成する場合、水面の波で電波が反射し、偽像をつくるときがある。この偽像については、例えばレーダ観測者の経験に基づき偽像を判断することで偽像による影響が回避されている。しかし、経験に基づく判断のみによると、偽像対策が不十分となってしまう可能性があるという課題がある。
【０００６】
本発明は、上記事情に鑑みてなされたものであり、水面波の偽像による影響を抑制することができる観測システム、水上観測システムおよび観測方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記課題を解決するため、本発明の一態様は、レーダによって検出した時系列のレーダ反射波の波形データを取得する反射波取得部と、前記時系列のレーダ反射波の波形データをフーリエ変換によって複数の周波数成分を含む成分波形データに変換するスペクトル解析部と、予め定められた一つ又は複数の水面波の波形を表す波モデルデータの中から、現在の水面波と合致する波モデル成分波形データを決定する波モデル成分波形データ決定部と、前記スペクトル解析部が変換した前記成分波形データにおいて、前記波モデル成分波形データ決定部が決定した前記波モデル成分波形データに含まれる周波数成分を抑制することで波成分抑制成分波形データを生成する波成分抑制部と、前記波成分抑制部が生成した前記波成分抑制成分波形データを逆フーリエ変換することで前記レーダによる水面に浮かぶ対象物の観測波を生成する観測波生成部と、を備える観測システムである。
【０００８】
また、本発明の一態様は、上記観測システムであって、前記波モデル成分波形データ決定部が、前記一つ又は複数の水面波の波形を表す波モデルデータに基づき、現在の風速または波高に応じて前記現在の水面波と合致する前記波モデル成分波形データを決定する。
【０００９】
また、本発明の一態様は、上記観測システムであって、前記一つ又は複数の水面波の波形を表す波モデルデータには、波候の情報と対応付けられた前記波モデルデータが含まれ、前記波モデル成分波形データ決定部が、前記一つ又は複数の水面波の波形を表す波モデルデータに基づき、現在の波候に応じて、前記現在の水面波と合致する前記波モデル成分波形データを決定する。
【００１０】
また、本発明の一態様は、レーダと、上記の観測システムと、水面波の波形を表す波モデルデータを予め格納する波モデルデータベースと、を備える水上観測システムである。
【００１１】
また、本発明の一態様は、レーダによって検出した時系列のレーダ反射波の波形データを取得する反射波取得ステップと、前記時系列のレーダ反射波の波形データをフーリエ変換によって複数の周波数成分を含む成分波形データに変換するスペクトル解析ステップと、予め定められた一つ又は複数の水面波の波形を表す波モデルデータの中から、現在の水面波と合致する波モデル成分波形データを決定する波モデル成分波形データ決定ステップと、前記スペクトル解析ステップにて変換した前記成分波形データにおいて、前記波モデル成分波形データ決定ステップにて決定した前記波モデル成分波形データに含まれる周波数成分を抑制することで波成分抑制成分波形データを生成する波成分抑制ステップと、前記波成分抑制ステップにて生成した前記波成分抑制成分波形データを逆フーリエ変換することで前記レーダによる水面に浮かぶ対象物の観測波を生成する観測波生成ステップと、を有する観測方法である。
【発明の効果】
【００１２】
本発明によれば、レーダによる水面に浮かぶ対象物の観測波における水面波の偽像による影響を抑制することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の一実施形態に係る観測装置を説明するためのブロック図である。
【図２】図１に示す観測装置１の動作例を説明するための模式図である。
【図３】図１に示す観測装置１の動作例を示すフローチャートである。
【発明を実施するための形態】
【００１４】
＜実施形態＞
本発明の実施形態について図１〜図３を参照して説明する。
図１は、本発明の一実施形態に係る観測装置１を説明するためのブロック図である。また、図２は、図１に示す観測装置１で処理される情報の流れを示す模式図である。
【００１５】
図１に示す水上観測システム１００は、観測装置１とレーダ２を備え、水面４に浮かぶ船等の対象物３の観測波を生成する。なお、この場合、水面４には、波高Ｈおよび周期Ｔの水面波４１が発生している。なお、以下では水面波（水面の高低運動）を波あるいは波浪等とも呼ぶ。一方、レーダ２は、アンテナ２１を有し、レーダ電波５１の発射方向を変化させながら所定のレーダ電波５１を間欠的に送信し、対象物３で反射したレーダ反射波５２を検出し、時系列のレーダ反射波５２を表す波形データを出力する。なお、レーダ２は、陸上あるいは海上等の水上に設置されたものであってもよいし、船上等の移動体上に設置されたものであってもよい。
【００１６】
観測装置１は、コンピュータであり、ＣＰＵ（中央処理装置）、記憶装置、入出力装置、通信装置等を備える。観測装置１は、記憶装置に記憶した所定のプログラムをＣＰＵで実行することで、図１に示す反射波取得部１１、スペクトル解析部１２、波成分抑制部１３、観測波生成部１４、および波モデル成分波形データ決定部１５として機能するとともに、記憶装置内に波モデルデータベース１６を記憶する。
【００１７】
反射波取得部１１は、レーダ２が検出した時系列のレーダ反射波の波形データ２０１（図２）を取得する。図２に示すように、時系列のレーダ反射波の波形データ２０１は、横軸を時間［ｓ］、縦軸を受信電力［ｄＢｍ］として、複数のレーダ反射波５２の受信電力の時間変化を表すデータである。反射波取得部１１は、例えば受信電力が所定のしきい値を越えてから下回るまでを一連の時系列データとして取得する。
【００１８】
スペクトル解析部１２は、反射波取得部１１が取得した時系列のレーダ反射波の波形データ２０１をフーリエ変換によって複数の周波数成分からなる成分波形データ２０２に変換する。図２に示すように、成分波形データ２０２は、例えば、横軸を周波数［Ｈｚ］、縦軸を受信電力［ｄＢｍ］として、波形データ２０１のスペクトルを表すデータである。
【００１９】
一方、波モデルデータベース１６は、水面波４１の波形を表す１または複数の波モデルデータ１７を予め格納するファイルである。波モデルデータ１７は、波モデル成分波形データ決定部１５において波モデル成分波形データ２０３を決定するために用いられる情報である。ここで、波モデル成分波形データ２０３は、水面波４１の波形を表す複数の周波数成分からなる情報であり、例えば横軸を周波数［Ｈｚ］、縦軸を受信電力［ｄＢｍ］として波浪のスペクトルを表すデータである。本実施形態において波モデルデータ１７は、例えば、（形式１）波モデルを表す情報、（形式２）波モデルと風速または波高さ以外のパラメータとの組み合わせを表す情報、または（形式３）波モデルとパラメータを用いて算出された算出結果とパラメータとを対応付けた情報である。波モデルデータベース１６は、（形式１）〜（形式３）の１種類または２種類以上が混在するように波モデルデータ１７を含む。なお、後述するように（形式１）の場合、波モデルデータ１７とは、ピアソン・モスコビッツスペクトルを表す関数Ｓ（ｆ）であって、波モデル成分波形データ決定部１５は風速又は波高のパラメータを関数Ｓ（ｆ）に入力して、現在の水面波と合致する波モデル成分波形データ２０３を算出する。また、（形式２）の波モデルデータ１７とは、ＪＯＮＳＷＡＰスペクトルを表す関数Ｓ（ｆ）であって、波モデル成分波形データ決定部１５は風速又は波高さ及びそれら以外のパラメータを関数Ｓ（ｆ）入力して波モデル成分波形データ２０３を算出する。また、（形式３）の波モデルデータ１７とは、予め（形式１）や（形式２）に基づいて算出した波モデル成分波形データ２０３を、それらの算出に用いたパラメータと対応付けて波モデルデータベース１６に記録しておき、波モデル成分波形データ決定部１５が、現在の波の状態を示すパラメータを用いて、対応する波モデル成分波形データ２０３を選択する。図２に示す例において、波モデルデータ１７は、波モデルを用いて算出された算出結果であり、横軸を周波数［Ｈｚ］、縦軸を周波数スペクトル密度関数（単にスペクトルともいう）Ｓ（ｆ）の値［ｍ２ｓ］とし、波浪の周波数スペクトルを表すデータである。ここで、ｆは周波数である。
【００２０】
本実施形態における波モデルは、一例として、波浪を表現する正弦波の重ね合わせからなる数値モデルであり、周波数スペクトル密度関数Ｓ（ｆ）で表すことができる。周波数スペクトル密度関数Ｓ（ｆ）は、例えば、過去の水面観測により得られた水面波の周波数スペクトルである。波モデルとしては、ピアソン・モスコビッツ（Ｐｉｅｒｓｏｎ−Ｍｏｓｋｏｗｉｔｚ）スペクトル、ＪＯＮＳＷＡＰスペクトル等がある。ピアソン・モスコビッツスペクトルは、風速をパラメータとし、風波が十分に発達した状態を表すスペクトルである。ＪＯＮＳＷＡＰスペクトルは、風速のほか、スペクトルのピーク付近の増幅率γ、スペクトルピーク付近のバンド幅ρ等をパラメータとし、風波の発達過程の状態を表すスペクトルである。なお、ピアソン・モスコビッツスペクトルやＪＯＮＳＷＡＰスペクトルにおける風速のパラメータは、統計的に、波高Ｈ（おおよそ風速の二乗に比例）または波高（有義波高）Ｈと周期Ｔの組み合わせに置き換えることができる。
【００２１】
なお、波浪は、風波とうねりに分類される。風波は風からエネルギを受けている波である。うねりは風からエネルギを受けずに伝搬している波である。また、波浪は、上述した波高Ｈと周期Ｔのほか、波向、波長、周波数、波速、波形勾配、波令、波数、角周波数等の特性によって定義することができる。また、風波は、風速のほか、吹走距離と吹続時間（連吹時間）によって発達の度合いが変化する。また、水深や海底地形、海流の影響によって波長や波高は変化する。また、ある場所における統計的な波浪特性を波候という。すなわち、波候は、特定地点における波浪特性の月別、季節別、年別等の所定期間の統計解析の結果を表す。また、場所や期間が波候を特定する情報である。
【００２２】
本実施形態において波モデルデータ１７を波モデルを表す情報とする場合（形式１）、波モデルデータ１７は、例えばピアソン・モスコビッツスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）を示す情報とすることができる。この場合、波モデル成分波形データ決定部１５は、例えば、現在の風速、波高または波高と周期を用いて、波モデルデータ１７である周波数スペクトル密度関数Ｓ（ｆ）を算出することで、波モデル成分波形データ２０３を決定することができる。
【００２３】
また、波モデルデータ１７を波モデルとパラメータとの組み合わせを表す情報とする場合（形式２）、波モデルデータ１７は、例えばＪＯＮＳＷＡＰスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）を示す情報と、上述した増幅率γ、バンド幅ρ等のパラメータの値を表す情報を含むとともに、波候を特定する場所や期間を示す情報に対応づけた情報とすることができる。この場合、例えば、第１の波モデルデータ１７は、第１の場所に対応づけられていて、第１の場所における予め定められた適切な第１の増幅率γ１、第１のバンド幅ρ１等を表す情報と、ＪＯＮＳＷＡＰスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）を表す情報を含む。また、例えば、第２の波モデルデータ１７は、第２の場所に対応づけられていて、第２の場所における予め定められた適切な第２の増幅率γ２、第２のバンド幅ρ２等を表す情報と、ＪＯＮＳＷＡＰスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）を表す情報を含む。この場合、波モデル成分波形データ決定部１５は、例えば、波モデルデータベース１６から現在の場所に対応する波モデルデータ１７を選択した後、現在の風速、波高または波高と周期を用いて、さらに当該波モデルデータ１７が含む予め定められた現在の場所に対応する適切なパラメータを適用することで、当該波モデルデータ１７が含む周波数スペクトル密度関数Ｓ（ｆ）を算出することで、波モデル成分波形データ２０３を決定することができる。なお、場所ごとの適切な増幅率γ、バンド幅ρ等のパラメータが予め定められていない場合、増幅率γ等のパラメータを外部から与えてもよい。
【００２４】
また、波モデルデータ１７を波モデルとパラメータを用いて算出された算出結果を示す情報とする場合（形式３）、波モデルデータ１７は、（形式１）または（形式２）の波モデルデータ１７が含む周波数スペクトル密度関数Ｓ（ｆ）を所定のパラメータを用いて所定の周波数範囲にわたって算出した結果を、使用したパラメータや当該波モデルデータ１７に対応づけられた波候を特定する情報に対応づけた情報とすることができる。この場合、例えば第３の波モデルデータ１７は、第３の風速に対応づけられた所定の周波数範囲のピアソン・モスコビッツスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）の算出結果を含む。また、第４の波モデルデータ１７は、第４の風速に対応づけられた所定の周波数範囲のピアソン・モスコビッツスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）の算出結果を含む。あるいは、例えば第５の波モデルデータ１７は、第５の風速と第５の波候を特定する情報に対応づけられた所定の周波数範囲のＪＯＮＳＷＡＰスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）の算出結果を含む。また、第６の波モデルデータ１７は、第６の風速と第６の波候を特定する情報に対応づけられた所定の周波数範囲のＪＯＮＳＷＡＰスペクトルを表す周波数スペクトル密度関数Ｓ（ｆ）の算出結果を含む。この場合、波モデル成分波形データ決定部１５は、例えば、波モデルデータベース１６から現在の風速または現在の風速と場所や期間に対応する波モデルデータ１７を選択することで、波モデル成分波形データ２０３を決定することができる。
【００２５】
波モデル成分波形データ決定部１５は、波モデルデータベース１６が記憶する波モデルデータ１７に基づき、現在の水面波と合致する波モデル成分波形データ２０３を決定する。波モデル成分波形データ決定部１５は、例えば、現在の風速、波高または波高と周期を用いて、波モデルデータ１７が表す周波数スペクトル密度関数Ｓ（ｆ）を算出し、さらに単位換算等を行って、波モデル成分波形データ２０３を決定する（形式１）。その際、波モデル成分波形データ決定部１５は、周波数スペクトル密度関数Ｓ（ｆ）のエネルギがピークの周波数ｆ２が、成分波形データ２０２のエネルギがピークの周波数ｆ１と一致するように、周波数軸を伸縮させたり、オフセットさせたりする変換処理を行ってもよい。あるいは、波モデル成分波形データ決定部１５は、例えば、波モデルデータベース１６から現在の場所に対応する波モデルデータ１７を選択した後、現在の風速、波高または波高と周期を用いて、当該波モデルデータ１７が含む現在の場所に対応するパラメータにて当該波モデルデータ１７が含む周波数スペクトル密度関数Ｓ（ｆ）を算出し、さらに単位換算等を行って、波モデル成分波形データ２０３を決定する（形式２）。あるいは、波モデル成分波形データ決定部１５は、例えば、波モデルデータベース１６から現在の風速または現在の風速と場所や期間に対応する波モデルデータ１７を選択し、さらに単位換算等を行って、波モデル成分波形データ２０３を決定する（形式３）。
【００２６】
なお、波モデル成分波形データ決定部１５が波モデル成分波形データ２０３を決定する際に用いるパラメータ等の情報は、例えばレーダ観測者が手作業で観測装置１に対して入力してもよいし、風速計、撮像装置および画像認識装置、全地球測位システム（ＧＰＳ）等を用いて波モデル成分波形データ決定部１５が自動的に取得できるようにしてもよい。
また、波モデル成分波形データ決定部１５が波モデル成分波形データ２０３を決定する方法は、上記の（形式１）〜（形式３）に限られず、例えば、対象となる海域の波の成分波形データが分かっているような場合、外部から波の成分波形データを入力し、波モデル成分波形データ決定部１５は、その入力された波の成分波形データを、波モデル成分波形データ２０３として決定してもよい。
【００２７】
波成分抑制部１３は、スペクトル解析部１２が変換した成分波形データ２０２において波モデル成分波形データ決定部１５が決定した波モデル成分波形データ２０３が含む周波数成分を抑制（除去）することで波成分抑制成分波形データ２０４を生成する。波成分抑制成分波形データ２０４は、例えば、横軸を周波数［Ｈｚ］、縦軸を受信電力［ｄＢｍ］として、成分波形データ２０２から波の周波成分を抑制した後のスペクトルを表すデータである。
【００２８】
観測波生成部１４は、波成分抑制部１３が生成した波成分抑制成分波形データ２０４を逆フーリエ変換することでレーダ２による水面４に浮かぶ対象物３の観測波２０５を生成する。観測波２０５は、横軸を時間［ｓ］、縦軸を受信電力［ｄＢｍ］として、波の周波成分を抑制した後の複数のレーダ反射波５２の受信電力の時間変化を表すデータである。観測波生成部１４が生成した観測波２０５は、観測装置１が備える表示装置等の出力装置に出力される。
【００２９】
次に、図３を参照して、観測装置１の動作例について説明する。観測装置１において、反射波取得部１１が時系列のレーダ反射波の波形データ２０１を取得すると（ステップＳ１）、スペクトル解析部１２が、時系列のレーダ反射波の波形データ２０１をフーリエ変換によって、複数の周波数成分からなる成分波形データ２０２に変換する（ステップＳ２）。次に、波モデル成分波形データ決定部１５が、波モデルデータベース１６が記憶する波モデルデータ１７に基づき、現在の水面波４１と合致する波モデル成分波形データ２０３を決定する（ステップＳ３）。次に、波成分抑制部１３が、スペクトル解析部１２が変換した成分波形データ２０２において波モデル成分波形データ２０３が含む周波数成分を抑制することで波成分抑制成分波形データ２０４を生成する（ステップＳ４）。次に、観測波生成部１４が、波成分抑制成分波形データ２０４を逆フーリエ変換することでレーダ２による水面４に浮かぶ対象物３の観測波２０５を生成し、所定の出力装置に出力する（ステップＳ５）。
【００３０】
以上のように、本実施形態の観測装置１では、反射波取得部１１が、レーダ２によって検出した時系列のレーダ反射波の波形データ２０１を取得する。スペクトル解析部１２が、時系列のレーダ反射波の波形データ２０１をフーリエ変換によって複数の周波数成分からなる成分波形データ２０２に変換する。波モデルデータベース１６が、水面波４１の波形を表す波モデルデータ１７を予め格納する。波モデル成分波形データ決定部１５が、波モデルデータベース１６が記憶する波モデルデータ１７に基づき、現在の水面波４１と合致する波モデル成分波形データ２０３を決定する。波成分抑制部１３が、スペクトル解析部１２が変換した成分波形データ２０２において波モデル成分波形データ決定部１５が決定した波モデル成分波形データ２０３が含む周波数成分を抑制することで波成分抑制成分波形データ２０４を生成する。観測波生成部１４が、波成分抑制部１３が生成した波成分抑制成分波形データ２０４を逆フーリエ変換することでレーダ２による水面４に浮かぶ対象物３の観測波２０５を生成する。したがって、本実施形態によれば、レーダ２による水面４に浮かぶ対象物３の観測波２０５における水面波４１による偽像による影響を抑制することができる。
【００３１】
なお、本実施形態において、波モデルデータベース１６は、波モデルデータ１７を風速や波高と対応付けて予め複数記憶することができ、波モデル成分波形データ決定部１５は、波モデルデータベース１６が記憶する波モデルデータ１７に基づき、現在の風速および波高のうち少なくとも一方に応じて現在の水面波４１と合致する波モデル成分波形データ２０３を決定することができる。また、波モデルデータベース１６は、波モデルデータ１７を波候に対応付けて予め複数記憶することができ、波モデル成分波形データ決定部１５が、波モデルデータベース１６が記憶する波モデルデータ１７に基づき、（現在の波候を特定する情報を指定することで）現在の波候に応じて、現在の水面波４１と合致する波モデル成分波形データ２０３を決定することができる。
【００３２】
なお、本発明の実施形態は上記のものに限定されず、例えば、以下のような変形が可能である。すなわち、例えば図１に示す波モデルデータベース１６は、例えば観測装置１の外部の記憶装置に格納されていて通信回線等を介して観測装置１に接続されていてもよい。また、波モデルは、上記のものに限定されず、例えばＴＭＡスペクトル等であってもよい。
【００３３】
また、本実施形態に係る水上観測システム１００の適用例として、水上観測システム１００を船Ｆ（図示せず）に搭載することが可能である。船Ｆはレーダ２から水面上に対してレーダ電波５１を発射し、それに対するレーダ反射波５２を受信する。そして、船Ｆ上から水面上を航行する他の対象物３を観測する。また、海岸沿いの陸上に建設された観測塔Ｇ（図示せず）や海底に固定された海洋観測塔Ｈ（図示せず）に水上観測システム１００を搭載することが可能である。観測塔Ｇ等のレーダ２から水面上に対して発射したレーダ電波５１のレーダ反射波５２を受信して、水面上を航行する対象物３を観測する。本実施形態の水上観測システム１００によれば、何れの場合も水面の波浪による偽像の影響を受けることなく正確に水面上に存在する船等の対象物を観測することができる。
【００３４】
なお、観測装置１は観測システムの一例である。上述した観測装置１における各処理の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムを観測装置１のコンピュータが読み出して実行することによって、上記処理が行われる。ここでコンピュータ読み取り可能な記録媒体とは、磁気ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、半導体メモリ等をいう。また、このコンピュータプログラムを通信回線によってコンピュータに配信し、この配信を受けたコンピュータが当該プログラムを実行するようにしてもよい。
【００３５】
また、上記プログラムは、前述した機能の一部を実現するためのものであってもよい。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるもの、いわゆる差分ファイル（差分プログラム）であってもよい。
また、観測装置１は、１台のコンピュータで構成されていても良いし、通信可能に接続された複数のコンピュータで構成されていてもよい。
【００３６】
その他、本発明の趣旨を逸脱しない範囲で、上記した実施の形態における構成要素を周知の構成要素に置き換えることは適宜可能である。また、この発明の技術範囲は上記の実施形態に限られるものではなく、本発明の趣旨を逸脱しない範囲において種々の変更を加えることが可能である。
【符号の説明】
【００３７】
１００  水上観測システム
１  観測装置
１１  反射波取得部
１２  スペクトル解析部
１３  波成分抑制部
１４  観測波生成部
１５  波モデル成分波形データ決定部
１６  波モデルデータベース
１７  波モデルデータ
２  レーダ
３  対象物
４  水面
４１  水面波
【図１】
2018146540.tif 000003
【図２】
2018146540.tif 000004
【図３】
2018146540.tif 000005