JPH0228116B2 - - Google Patents
Info
- Publication number
- JPH0228116B2 JPH0228116B2 JP58131339A JP13133983A JPH0228116B2 JP H0228116 B2 JPH0228116 B2 JP H0228116B2 JP 58131339 A JP58131339 A JP 58131339A JP 13133983 A JP13133983 A JP 13133983A JP H0228116 B2 JPH0228116 B2 JP H0228116B2
- Authority
- JP
- Japan
- Prior art keywords
- wave height
- output
- intensity
- circuit
- sea surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000004364 calculation method Methods 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 9
- 238000000034 method Methods 0.000 claims description 8
- 230000000737 periodic effect Effects 0.000 claims description 7
- 238000010183 spectrum analysis Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 description 9
- 230000003321 amplification Effects 0.000 description 6
- 238000003199 nucleic acid amplification method Methods 0.000 description 6
- 238000005311 autocorrelation function Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000007781 pre-processing Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/95—Radar or analogous systems specially adapted for specific applications for meteorological use
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Description
【発明の詳細な説明】
本発明はマイクロ波パルスを送受信して、受信
波の対数増幅検波出力を情報として利用する船舶
用レーダ装置に接続して、海面反射信号を計算・
処理し、海面のうねりの高さの測定を可能とす
る、波高レーダ観測方法及び装置に関する。
波の対数増幅検波出力を情報として利用する船舶
用レーダ装置に接続して、海面反射信号を計算・
処理し、海面のうねりの高さの測定を可能とす
る、波高レーダ観測方法及び装置に関する。
海上を航行中の船舶が、暗夜や荒天で視界の悪
い時にも、その周囲のうねりの高さを即時に自動
的に測定できる実用的な装置があれば、船舶の安
全運航に役立たせることができる。しかし、従来
そのような装置は存在しなかつた。
い時にも、その周囲のうねりの高さを即時に自動
的に測定できる実用的な装置があれば、船舶の安
全運航に役立たせることができる。しかし、従来
そのような装置は存在しなかつた。
一般に船舶用レーダでは、波浪が強くなると、
周囲の海面からの反射信号が観測できるが、それ
は主に、風によつて海面に発生したさざ波（以
下、風浪と称する）がレーダ電波を反射している
状態であり、反射信号の強さとうねりの高さには
直接の関係がないとされていた。
周囲の海面からの反射信号が観測できるが、それ
は主に、風によつて海面に発生したさざ波（以
下、風浪と称する）がレーダ電波を反射している
状態であり、反射信号の強さとうねりの高さには
直接の関係がないとされていた。
以下、本発明の原理について説明する。うねり
の高さが数ｍ以上で、その表面に風浪が重畳して
いる時の海面の状況は第１図のようになつてい
る。このような海面における船舶用レーダのPPI
画像の例を第２図に示す、うねりの波峰線と関連
する周期的パターンが表われている。これは、第
１図の左側から、水平面に対して数度の角度でレ
ーダ電波が入射しているとすると、波峰の左側で
は反射が強く、波峰の右側では反射が弱くなるた
めである。従つて、レーダで観測できる海面反射
信号は風浪が激しいほど、うねりが高いほど強く
なる。このような関係を次式で表わしてみる。
の高さが数ｍ以上で、その表面に風浪が重畳して
いる時の海面の状況は第１図のようになつてい
る。このような海面における船舶用レーダのPPI
画像の例を第２図に示す、うねりの波峰線と関連
する周期的パターンが表われている。これは、第
１図の左側から、水平面に対して数度の角度でレ
ーダ電波が入射しているとすると、波峰の左側で
は反射が強く、波峰の右側では反射が弱くなるた
めである。従つて、レーダで観測できる海面反射
信号は風浪が激しいほど、うねりが高いほど強く
なる。このような関係を次式で表わしてみる。
E(t)＝Ew(t)×Es(t) ……(1)
ここに、E(t)：レーダのパルスを送信後、時間
ｔにおける海面反射強度。レーダからの距離ｘ＝
0.5ct（ｃは光速）の関係を有する。
ｔにおける海面反射強度。レーダからの距離ｘ＝
0.5ct（ｃは光速）の関係を有する。
Ew(t)：風速による反射成分で不規則に速く変動
する。
する。
Es(t)：うねりによる反射成分で、周期的にゆつ
くり変動する。
くり変動する。
時間領域で(1)式の関係を持つと、周波数領域で
は、 S(f)＝∫+∞ -∞Sw(〓)Ss（ｆ−λ）dλ……(2) ここに、S(f)：海面反射強度のスペクトルで、
E(t)のフーリエ変換値、すなわち、S(f)＝∫+∞ -∞E(t)
exp（−2πjft）dt ……(2)′ Sw(f)：風浪反射成分のスペクトルで、Ew(t)のフ
ーリエ変換値 Ss(f)：うねりによる反射成分のスペクトルで、
Es(t)のフーリエ変換値 (2)式の右辺はたたみ込み積分である。
は、 S(f)＝∫+∞ -∞Sw(〓)Ss（ｆ−λ）dλ……(2) ここに、S(f)：海面反射強度のスペクトルで、
E(t)のフーリエ変換値、すなわち、S(f)＝∫+∞ -∞E(t)
exp（−2πjft）dt ……(2)′ Sw(f)：風浪反射成分のスペクトルで、Ew(t)のフ
ーリエ変換値 Ss(f)：うねりによる反射成分のスペクトルで、
Es(t)のフーリエ変換値 (2)式の右辺はたたみ込み積分である。
風浪による反射成分Ew(t)は、時間領域で不規則
に速く変動するので、そのスペクトル分布を一様
分布と仮定して、 Sw(f)＝Ｗ（一定） ……(3) とすると、(2)式より Ss(f)＝S(f)／Ｗ ……(4) となつて、うねりによる反射のスペクトルSs(f)が
求まる。
に速く変動するので、そのスペクトル分布を一様
分布と仮定して、 Sw(f)＝Ｗ（一定） ……(3) とすると、(2)式より Ss(f)＝S(f)／Ｗ ……(4) となつて、うねりによる反射のスペクトルSs(f)が
求まる。
以下、(4)式について具体的に説明する。第２図
のレーダ画像で、点線にそつた海面反射強度のス
ペクトルS(f)は第３図のようになる。第３図で横
軸の周波数がＦのスペクトルS(F)は、うねりの高
さと関連する、周波数Ｆから離れた領域のスペク
トルは、風浪の激しさと関連する（図中のＷ）。
この場合、うねりの高さはS(F)／Ｗと比例関係を
有する。
のレーダ画像で、点線にそつた海面反射強度のス
ペクトルS(f)は第３図のようになる。第３図で横
軸の周波数がＦのスペクトルS(F)は、うねりの高
さと関連する、周波数Ｆから離れた領域のスペク
トルは、風浪の激しさと関連する（図中のＷ）。
この場合、うねりの高さはS(F)／Ｗと比例関係を
有する。
次に、うねりの高さを自動的に計算する方法に
ついて説明する。スペクトルの最大値S(F)は容易
に計算できる。風浪強度Ｗは次の近似式で算出す
る。
ついて説明する。スペクトルの最大値S(F)は容易
に計算できる。風浪強度Ｗは次の近似式で算出す
る。
Ｗ≒｛１／Fmax∫Fmax 0｜S(f)｜2df｝1/2 ……(5)
この場合、(5)式の右辺の積分値には、うねりの
成分S(F)の影響は少ないものと仮定している。パ
ーシブルの定理によりＷは時間領域でも容易に求
めることができる、すなわち、 Ｗ≒｛１／Tmax∫Tmax 0｜E(t)｜2dt｝1/2 ……(6) (6)式のＷは海面反射強度の、ある距離範囲での
実効値であり、従来、風浪との相関が強いと言わ
れて来た。
成分S(F)の影響は少ないものと仮定している。パ
ーシブルの定理によりＷは時間領域でも容易に求
めることができる、すなわち、 Ｗ≒｛１／Tmax∫Tmax 0｜E(t)｜2dt｝1/2 ……(6) (6)式のＷは海面反射強度の、ある距離範囲での
実効値であり、従来、風浪との相関が強いと言わ
れて来た。
以下、実測例について説明する。第４図の縦軸
は風速、横軸は観測時刻である。図中の実線は風
速計による測定値、点線は(6)式に従つて求めた海
面反射強度の実効値Ｗを、風速に直線回帰させた
結果である。このとき、風速計と海面反射強度の
実効値から求めた風速との相関係数は0.77で、強
い関連を示している。これは今まで言われて来た
関係を再確認したものである。従来の理論はこの
段階までであつた。
は風速、横軸は観測時刻である。図中の実線は風
速計による測定値、点線は(6)式に従つて求めた海
面反射強度の実効値Ｗを、風速に直線回帰させた
結果である。このとき、風速計と海面反射強度の
実効値から求めた風速との相関係数は0.77で、強
い関連を示している。これは今まで言われて来た
関係を再確認したものである。従来の理論はこの
段階までであつた。
第５図の縦軸は平均波高、横軸は確測時刻であ
る。図中の実線は波高計による測定値、点線は(4)
式および(6)式によつて求めたS(F)／Ｗを波高に直
線回帰させた結果である。このとき、波高計とレ
ーダのデータから算出した波高との相関係数は
0.84、実効誤差は0.46mであつた。そして、第４
図の風速と第５図の波高の関連はほとんどない。
これは、(1)式〜(6)式によつて説明した本発明の波
高測定の原理が正しい事を実際の波の観測で証明
するものである。
る。図中の実線は波高計による測定値、点線は(4)
式および(6)式によつて求めたS(F)／Ｗを波高に直
線回帰させた結果である。このとき、波高計とレ
ーダのデータから算出した波高との相関係数は
0.84、実効誤差は0.46mであつた。そして、第４
図の風速と第５図の波高の関連はほとんどない。
これは、(1)式〜(6)式によつて説明した本発明の波
高測定の原理が正しい事を実際の波の観測で証明
するものである。
ここまでの原理の概要は、海面反射強度のスペ
クトル解析をして、周波数領域でうねりの周期成
分の強度を求め、それを風浪と関連する反射強度
の実効値で除算して、うねりの高さを求めるとい
うものである。これと同じ結果を時間領域の計算
で求める方法も容易に考える事ができる。海面反
射強度E(t)の自己相関関数φEE(〓)は次のように定義
できる。
クトル解析をして、周波数領域でうねりの周期成
分の強度を求め、それを風浪と関連する反射強度
の実効値で除算して、うねりの高さを求めるとい
うものである。これと同じ結果を時間領域の計算
で求める方法も容易に考える事ができる。海面反
射強度E(t)の自己相関関数φEE(〓)は次のように定義
できる。
φEE(〓)＝１／Ｔ∫T 0E(t)E(t-〓)dτ ……(7)
ここに、τ：時間差
Ｔ：観測時間
第２図のレーダ画像で、点線にそつた海面反射
強度の自己相関関数φEE(〓)は第６図のようになる。
第６図のφEE(〓)をフーリエ変換すると第３図のス
ペクトラムが求まる。第６図からも海面反射強度
の実効値Ｗ、うねりの周期成分の強度S(F)を求め
て、うねりの波高を算出することができる。以
下、スペクトル解析による実施例について詳細に
説明する。
強度の自己相関関数φEE(〓)は第６図のようになる。
第６図のφEE(〓)をフーリエ変換すると第３図のス
ペクトラムが求まる。第６図からも海面反射強度
の実効値Ｗ、うねりの周期成分の強度S(F)を求め
て、うねりの波高を算出することができる。以
下、スペクトル解析による実施例について詳細に
説明する。
第７図は本発明の一実施例である。１は送信
機、２は送受切換器、３はアンテナ、４はミク
サ、５は対数増幅検波器、６はA/D変換器、７
は距離補正回路、８は逆対数変換回路、９はウイ
ンドウ処理回路、１０は実効値計算回路、１１は
フーリエ変換回路、１２は最大値検出回路、１３
は除算器、１４は波高計算回路、１５は波高出力
端子である。送信機１、送受切換器２、アンテナ
３、ミクサ４は従来の一般的な船舶用レーダと同
一のものでよい。
機、２は送受切換器、３はアンテナ、４はミク
サ、５は対数増幅検波器、６はA/D変換器、７
は距離補正回路、８は逆対数変換回路、９はウイ
ンドウ処理回路、１０は実効値計算回路、１１は
フーリエ変換回路、１２は最大値検出回路、１３
は除算器、１４は波高計算回路、１５は波高出力
端子である。送信機１、送受切換器２、アンテナ
３、ミクサ４は従来の一般的な船舶用レーダと同
一のものでよい。
対数増幅検波器５は、ミクサ４の出力である中
間周波数信号を増幅検波するもので、その入力と
出力が対数関数の関係を有するものである。これ
により広範囲の振幅を持つ海面反射信号を飽和す
ることなく増幅検波できる。
間周波数信号を増幅検波するもので、その入力と
出力が対数関数の関係を有するものである。これ
により広範囲の振幅を持つ海面反射信号を飽和す
ることなく増幅検波できる。
A/D変換器６は対数増幅検波器のアナログ出
力をデイジタル変換する。船舶の安全運航のため
のおこなう波高観測のうねりの最短波長は30m位
であり、これを測定するためのレーダの距離分解
能は15m以下（送信パルス幅0.1μS以下）である。
従つてA/D変換器６が反射信号をサンプリング
する周期も0.1μS以下にする必要がある。デイジ
タル化の精度は６〜８ビツト、距離方向のサンプ
リング数は64〜128個でよい。
力をデイジタル変換する。船舶の安全運航のため
のおこなう波高観測のうねりの最短波長は30m位
であり、これを測定するためのレーダの距離分解
能は15m以下（送信パルス幅0.1μS以下）である。
従つてA/D変換器６が反射信号をサンプリング
する周期も0.1μS以下にする必要がある。デイジ
タル化の精度は６〜８ビツト、距離方向のサンプ
リング数は64〜128個でよい。
海面の状態が一様であつても、レーダからの距
離が近いほど海面反射信号は強くなる。アンテナ
入力において海面反射強度は距離の３乗にほぼ逆
比例して変化する。従つて、距離補正は対数領域
の方が容易におこなえる。距離補正回路７は、
Ａ／Ｄ変換器６の出力に対して、距離方向の強度
補正をおこなう。
離が近いほど海面反射信号は強くなる。アンテナ
入力において海面反射強度は距離の３乗にほぼ逆
比例して変化する。従つて、距離補正は対数領域
の方が容易におこなえる。距離補正回路７は、
Ａ／Ｄ変換器６の出力に対して、距離方向の強度
補正をおこなう。
逆対数変換器８は、距離補正回路７の出力を逆
対数変換するもので、その出力を受信波の包絡線
の振幅に比例したものにする（直線検波出力）。
これは、フーリエ変換のための前処理の一つであ
る。
対数変換するもので、その出力を受信波の包絡線
の振幅に比例したものにする（直線検波出力）。
これは、フーリエ変換のための前処理の一つであ
る。
ウインドウ処理回路９はフーリエ変換の際の周
波数応答を最適化するための重み付け処理で、
Hammingウインドウ関数などを用いる。
波数応答を最適化するための重み付け処理で、
Hammingウインドウ関数などを用いる。
実効値計算回路１０は(6)式に従つて、海面反射
強度の実効値Ｗを算出する。
強度の実効値Ｗを算出する。
フーリエ変換回路１１はウインドウ処理回路９
の出力をフーリエ変換する。フーリエ変換は基本
的には(2)′式で定義されるが、実際にはFFT（高
速フーリエ変換）計算法を使用する。
の出力をフーリエ変換する。フーリエ変換は基本
的には(2)′式で定義されるが、実際にはFFT（高
速フーリエ変換）計算法を使用する。
最大値検出回路１２は、フーリエ変換器１１の
出力から、うねりの周波数Ｆの成分の最大値S(F)
を検出する。
出力から、うねりの周波数Ｆの成分の最大値S(F)
を検出する。
除算器１３は、最大値検出回路１２の出力S(F)
を、実効値計算回路１０の出力Ｗで除算して、
S(F)／Ｗを出力する。
を、実効値計算回路１０の出力Ｗで除算して、
S(F)／Ｗを出力する。
波高計算回路１４は、除算器１３の出力Ｘを波
高Ｈに換算するための計算（Ｈ＝AX＋Ｂ）をお
こなう。Ａ、Ｂは回帰係数で、波高計のデータに
よりあらかじめ校正・算出しておく。このように
して、レーダによるうねりの波高の測定値Ｈは端
子１５に出力される。
高Ｈに換算するための計算（Ｈ＝AX＋Ｂ）をお
こなう。Ａ、Ｂは回帰係数で、波高計のデータに
よりあらかじめ校正・算出しておく。このように
して、レーダによるうねりの波高の測定値Ｈは端
子１５に出力される。
以上説明した第７図の実施例において、対数増
幅検波器５のかわりに直線増幅検波器も使用でき
る。その場合は８の逆対数変換器は不要となる
が、増幅器が飽和しないように利得を調節する必
要がある。距離補正回路７、逆対数変換器８、ウ
インドウ処理回路９は、フーリエ変換のための前
処理をおこなうものであるが、前処理の中に海面
反射信号の平均値（直流分）を除去する処理を含
めてもよい。また簡略化のために、その全部また
は一部を省略してもよい。また、処理の順番を変
更しても同様の効果を得ることができる。波高計
算回路１４の計算式（Ｈ＝AX＋Ｂ）は、実際の
波高測定精度を改良するために他の計算式（たと
えば、Ｈ＝AX2＋BX＋Ｃ）を使用してもよい。
７〜１４の演算装置は、ハードウエアまたはソフ
トウエアおよび両方の組合せによつて、容易に実
現できる。また、第７図は１次元での実施例であ
るが、これを２次元での処理に応用することは容
易である。(7)式で説明した自己相関関数による方
法は演算時間が長いので現状では第７図の実施例
の方法が実用的である。
幅検波器５のかわりに直線増幅検波器も使用でき
る。その場合は８の逆対数変換器は不要となる
が、増幅器が飽和しないように利得を調節する必
要がある。距離補正回路７、逆対数変換器８、ウ
インドウ処理回路９は、フーリエ変換のための前
処理をおこなうものであるが、前処理の中に海面
反射信号の平均値（直流分）を除去する処理を含
めてもよい。また簡略化のために、その全部また
は一部を省略してもよい。また、処理の順番を変
更しても同様の効果を得ることができる。波高計
算回路１４の計算式（Ｈ＝AX＋Ｂ）は、実際の
波高測定精度を改良するために他の計算式（たと
えば、Ｈ＝AX2＋BX＋Ｃ）を使用してもよい。
７〜１４の演算装置は、ハードウエアまたはソフ
トウエアおよび両方の組合せによつて、容易に実
現できる。また、第７図は１次元での実施例であ
るが、これを２次元での処理に応用することは容
易である。(7)式で説明した自己相関関数による方
法は演算時間が長いので現状では第７図の実施例
の方法が実用的である。
以上説明したように、本発明は船舶用レーダに
おける海面反射強度のスペクトル解析をして、周
波数領域でうねりの周期成分の強度を求め、それ
を風浪の激しさと関連する反射強度の実効値で除
算してうねりの高さを求めるもので、従来、不可
能であつた船舶用レーダによるうねりの高さの測
定を可能とするものである。そのための装置は容
易に構成できて、小形・低価格・高信頼性であ
り、船舶の安全運航・海難防止に役立たせること
ができる。
おける海面反射強度のスペクトル解析をして、周
波数領域でうねりの周期成分の強度を求め、それ
を風浪の激しさと関連する反射強度の実効値で除
算してうねりの高さを求めるもので、従来、不可
能であつた船舶用レーダによるうねりの高さの測
定を可能とするものである。そのための装置は容
易に構成できて、小形・低価格・高信頼性であ
り、船舶の安全運航・海難防止に役立たせること
ができる。
第１図はうねりと風浪が重畳しているときの海
面の状況を示す図、第２図はうねりの波峰線が表
われているレーダPPI画像の図、第３図は第２図
の点線の方向の海面反射信号のスペクトルの例を
示す図、第４図は風速計で観測した風速とレーダ
のデータから算出した風速を示す図、第５図は波
高計で観測した波高と本発明の原理に従つてレー
ダデータから算出した波高を示す図、第６図は第
２図の点線の方向の海面反射信号の自己相関関数
の例を説明する図、第７図は本発明の一実施例を
示す図である。 １……送信機、２……送受切換器、３……アン
テナ、４……ミクサ、５……対数増幅検波器、６
……A/D変換器、７……距離補正回路、８……
逆対数変換回路、９……ウインドウ処理回路、１
０……実効値計算回路、１１……フーリエ変換回
路、１２……最大値検出回路、１３……除算器、
１４……波高計算回路。
面の状況を示す図、第２図はうねりの波峰線が表
われているレーダPPI画像の図、第３図は第２図
の点線の方向の海面反射信号のスペクトルの例を
示す図、第４図は風速計で観測した風速とレーダ
のデータから算出した風速を示す図、第５図は波
高計で観測した波高と本発明の原理に従つてレー
ダデータから算出した波高を示す図、第６図は第
２図の点線の方向の海面反射信号の自己相関関数
の例を説明する図、第７図は本発明の一実施例を
示す図である。 １……送信機、２……送受切換器、３……アン
テナ、４……ミクサ、５……対数増幅検波器、６
……A/D変換器、７……距離補正回路、８……
逆対数変換回路、９……ウインドウ処理回路、１
０……実効値計算回路、１１……フーリエ変換回
路、１２……最大値検出回路、１３……除算器、
１４……波高計算回路。
Claims (1)
- 【特許請求の範囲】 １ 船舶用パルスレーダにおいて、受信した海面
反射信号を増幅検波した信号から海面反射強度の
スペクトル解析をし、周波数領域で最大値を求め
てうねりの周期成分の強度を算出し、他方におい
て前記検波した信号から風浪の強さに関連する反
射強度の実効値を求め、前記うねりの周期成分の
強度を前記反射強度の実効値で除算し、除算した
値から予め求めておいた回帰係数に基づく波高値
計算式に従つてうねりの波高を求めることを特徴
とする波高レーダ観測方法。 ２ 船舶用パルスレーダにおいて、受信した海面
反射信号を増幅検波した信号をデイジタル化する
Ａ／Ｄ変換器と、該A/D変換器の出力に対して距
離方向の強度補正をおこなう距離補正回路と、該
距離補正回路の出力に対して周波数応答を最適化
するための重み付けをするウインドウ処理回路
と、該ウインドウ処理回路の出力をフーリエ変換
するフーリエ変換器と、該フーリエ変換器の出力
からうねりの周期成分の強度を算出する最大値検
出回路と、他方において前記距離補正回路の出力
から風浪の強さに関連する反射強度の実効値を算
出する実効値計算回路と、前記最大値検出回路の
出力を前記実効値計算回路の出力で除算する除算
器と、該除算器の出力から波高を算出する波高計
算回路とを備えてうねりの波高を得ることを特徴
とする波高レーダ観測装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58131339A JPS6022681A (ja) | 1983-07-19 | 1983-07-19 | 波高レ−ダ観測方法及び装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP58131339A JPS6022681A (ja) | 1983-07-19 | 1983-07-19 | 波高レ−ダ観測方法及び装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS6022681A JPS6022681A (ja) | 1985-02-05 |
| JPH0228116B2 true JPH0228116B2 (ja) | 1990-06-21 |
Family
ID=15055625
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP58131339A Granted JPS6022681A (ja) | 1983-07-19 | 1983-07-19 | 波高レ−ダ観測方法及び装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPS6022681A (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6249214A (ja) * | 1985-08-29 | 1987-03-03 | Mitsui Eng & Shipbuild Co Ltd | 海上におけるうねりの観測方法 |
| JPH02118482A (ja) * | 1988-10-28 | 1990-05-02 | Japan Radio Co Ltd | レーダ式波高測定装置 |
| JP3533260B2 (ja) * | 1995-04-06 | 2004-05-31 | 日本無線株式会社 | 波浪周期検出装置 |
| JP3466908B2 (ja) * | 1998-02-20 | 2003-11-17 | 長野日本無線株式会社 | 短波または超短波を用いた海洋レーダ |
| JP4719373B2 (ja) * | 2001-04-23 | 2011-07-06 | 次清 平山 | 波高算出装置、波高算出方法及びコンピュータ可読媒体 |
| WO2019199326A1 (en) * | 2018-04-13 | 2019-10-17 | Hewlett-Packard Development Company, L.P. | Antenna direction weightings |
| CN111007572B (zh) * | 2019-11-22 | 2020-11-10 | 北京中科蓝图科技有限公司 | 道路地下空洞的自动识别方法、装置及系统 |
-
1983
- 1983-07-19 JP JP58131339A patent/JPS6022681A/ja active Granted
Also Published As
| Publication number | Publication date |
|---|---|
| JPS6022681A (ja) | 1985-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7317659B2 (en) | Measurement of air characteristics in the lower atmosphere | |
| US4318100A (en) | Automatic ground clutter rejection in weather pulse radar system | |
| JPH08506894A (ja) | 冶金工程における二以上の表面位置の同時測定方法 | |
| JPH0228116B2 (ja) | ||
| US4715707A (en) | Laser doppler velocimeter with laser backscatter discriminator | |
| US4139847A (en) | Automatic ground-clutter rejection in weather pulse radar system | |
| KR100661748B1 (ko) | Ｆｍｃｗ 레이더의 누설신호 제거장치 | |
| RU2585401C1 (ru) | Устройство обнаружения объектов в водной среде | |
| JPS60263880A (ja) | 地下埋設物探査方法 | |
| US3745517A (en) | Distance measuring device and method | |
| RU2510040C2 (ru) | Устройство для определения состояния морской поверхности | |
| JPH0341796B2 (ja) | ||
| US6229761B1 (en) | Estimating ship velocity through the water and over the ground | |
| JP2951045B2 (ja) | 超音波反射強度測定装置 | |
| JPH02115782A (ja) | 座礁予測距離装置 | |
| RU2019855C1 (ru) | Параметрический эхоледомер | |
| JP2653747B2 (ja) | シークラッタ抑圧方法 | |
| JPS6128108B2 (ja) | ||
| JPS6093317A (ja) | レ−ダ波高測定装置 | |
| JPH06273198A (ja) | 波浪特性抽出装置 | |
| JP2839406B2 (ja) | 方位測定装置 | |
| JPH0228115B2 (ja) | ||
| JPS61241685A (ja) | 超音波距離計 | |
| JPH045033Y2 (ja) | ||
| JPH0395477A (ja) | 超音波探知装置 |