JP6689961B2 - Signal processing device, radar device, and signal processing method - Google Patents
Signal processing device, radar device, and signal processing method Download PDFInfo
- Publication number
- JP6689961B2 JP6689961B2 JP2018511928A JP2018511928A JP6689961B2 JP 6689961 B2 JP6689961 B2 JP 6689961B2 JP 2018511928 A JP2018511928 A JP 2018511928A JP 2018511928 A JP2018511928 A JP 2018511928A JP 6689961 B2 JP6689961 B2 JP 6689961B2
- Authority
- JP
- Japan
- Prior art keywords
- wave
- spectrum
- evaluation
- reliability
- unit
- 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.)
- Active
Links
- 238000003672 processing method Methods 0.000 title claims description 7
- 238000001228 spectrum Methods 0.000 claims description 332
- 238000011156 evaluation Methods 0.000 claims description 314
- 238000004458 analytical method Methods 0.000 claims description 56
- 238000004364 calculation method Methods 0.000 claims description 34
- 238000000034 method Methods 0.000 claims description 32
- 238000000605 extraction Methods 0.000 claims description 20
- 239000006185 dispersion Substances 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 7
- 238000006243 chemical reaction Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 19
- 101100365087 Arabidopsis thaliana SCRA gene Proteins 0.000 description 16
- 238000001514 detection method Methods 0.000 description 9
- 230000010354 integration Effects 0.000 description 9
- 238000011045 prefiltration Methods 0.000 description 7
- 238000012986 modification Methods 0.000 description 6
- 230000004048 modification Effects 0.000 description 6
- 238000013500 data storage Methods 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 2
- 230000012447 hatching Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
Images
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
- 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
本発明は、波浪からのエコーを処理して波向などの波浪情報を算出する信号処理装置、この信号処理装置を備えたレーダ装置、および、信号処理方法に関する。 The present invention relates to a signal processing device that processes an echo from a wave to calculate wave information such as a wave direction, a radar device including the signal processing device, and a signal processing method.
従来、波浪に関する情報としての波浪情報を得るための装置として、特許文献１に開示される波浪観測レーダ（すなわち、レーダ装置）が知られている。このレーダ装置では、特許文献１の図２を参照して、海上における所定範囲からのエコーから得られたエコー信号が２次元ＦＦＴ処理されることにより、特許文献１の図３を参照して、２次元フーリエ変換信号Ｓｆが導出される。そして、このレーダ装置では、その２次元フーリエ変換信号Ｓｆに基づき、波浪情報（波向、波長など）が算出される。 BACKGROUND ART Conventionally, a wave observation radar (that is, a radar device) disclosed in
ところで、波浪観測レーダからユーザに示される波浪情報について、実際の波浪状態を正確に示しているか否か、ユーザが判断し難い場合がある。すなわち、ユーザなどにおいて、波浪情報の信頼性を判断し難い場合がある。 By the way, there are cases where it is difficult for the user to determine whether or not the actual wave state is accurately indicated for the wave information shown to the user by the wave observation radar. That is, it may be difficult for the user or the like to judge the reliability of the wave information.
本発明は、上記課題に鑑みてなされており、その目的の１つは、波浪情報を信頼できるか否かについて、より確実に評価するための構成を有する信号処理装置、レーダ装置、および、信号処理方法を提供することにある。 The present invention has been made in view of the above problems, and one of the objects thereof is a signal processing device, a radar device, and a signal having a configuration for more reliably evaluating whether or not wave information is reliable. It is to provide a processing method.
上記課題を解決するため、本発明のある局面に係る信号処理装置は、所定の観測エリアにおける波浪を観測することで得られた観測データに基づいて、波数と周波数との関係を含む波数データを取得する波数データ取得部と、前記波数データに所定のフィルタ処理を施すことで前記波数データに含まれる波浪データを検出する波浪情報取得部と、前記観測エリアの一部である所定の解析領域について、前記波浪データで特定される波浪情報の信頼性を評価する、信頼性評価部と、を備えている。 To solve the above problems, the signal processing device according to an aspect of the present invention, based on the observation data obtained by observing the waves in a predetermined observation area, the wave number data including the relationship between the wave number and the frequency. Wave number data acquisition unit to obtain, the wave information acquisition unit to detect the wave data included in the wave number data by performing a predetermined filtering process on the wave number data, and a predetermined analysis region that is part of the observation area And a reliability evaluation unit that evaluates the reliability of the wave information specified by the wave data.
また、上記課題を解決するため、本発明のある局面に係る信号処理方法は、所定の観測エリアにおける波浪を観測することで得られた観測データに基づいて、波数と周波数との関係を含む波数データを取得し、前記波数データに所定のフィルタ処理を施すことで前記波数データに含まれる波浪データを検出し、前記観測エリアの一部である所定の解析領域について、前記波浪データで特定される波浪情報の信頼性を評価する。 Further, in order to solve the above problems, a signal processing method according to an aspect of the present invention, based on the observation data obtained by observing the waves in a predetermined observation area, the wave number including the relationship between the wave number and the frequency. Data is acquired, wave data included in the wave number data is detected by performing a predetermined filtering process on the wave number data, and a predetermined analysis region that is a part of the observation area is specified by the wave data. Evaluate the reliability of wave information.
本発明によれば、波浪情報を信頼できるか否かについて、より確実に評価するための構成を実現できる。 According to the present invention, it is possible to realize a configuration for more reliably evaluating whether or not wave information is reliable.
以下、本発明に係る信号処理装置、およびこの信号処理装置を備えたレーダ装置の実施形態について図面を参照しつつ説明する。本発明は、波浪情報を処理するための信号処理装置、この信号処理装置を備えたレーダ装置、および、波浪情報を処理するための信号処理方法に広く適用することができる。 Hereinafter, embodiments of a signal processing device according to the present invention and a radar device including the signal processing device will be described with reference to the drawings. INDUSTRIAL APPLICABILITY The present invention can be widely applied to a signal processing device for processing wave information, a radar device provided with this signal processing device, and a signal processing method for processing wave information.
図１は、本発明の実施形態に係るレーダ装置１の信号処理装置２０の模式的な構成を示すブロック図である。図２は、レーダ装置１の模式的な構成を示すブロック図である。図３は、信号処理装置２０の第１信頼性評価ユニット２１の模式的な構成を示すブロック図である。図４は、レーダ装置１に備えられる信号処理装置２０の第２信頼性評価ユニット２２の模式的な構成を示すブロック図である。 FIG. 1 is a block diagram showing a schematic configuration of a
図２を参照して、本実施形態に係るレーダ装置１は、送波した送信波が波浪に反射して帰来するエコーに基づいて、波浪に関する情報である波浪情報（具体的には、波向など）を算出する。本実施形態のレーダ装置１は、たとえば、漁船などの船舶としての自船に備えられている。なお、レーダ装置１は、自船の周囲に存在する船舶を探知可能に構成されていてもよい。 With reference to FIG. 2, the
レーダ装置１は、アンテナユニット２と、信号処理装置２０と、表示器３と、を備えている。 The
アンテナユニット２は、送波器および受波器の双方として機能するアンテナ４と、受信部５と、Ａ／Ｄ変換部６と、を含んでいる。 The
アンテナ４は、指向性の強い送信波としてのパルス状電波を送波可能なレーダアンテナである。また、アンテナ４は、当該アンテナ４から送波されたパルス状電波が物標（本実施形態の場合、波浪）に反射して帰来する反射波（すなわち、エコー）を受波するように構成されている。レーダ装置１は、パルス状電波を送波してから反射波を受波するまでの時間を測定する。これにより、レーダ装置１は、物標までの距離ｒを検出することができる。アンテナ４は、水平面上で３６０°回転可能に構成されている。アンテナ４は、パルス状電波の送波方向を変えながら（たとえば、アンテナ角度を変えながら）、電波の送受波を繰り返し行うように構成されている。以上の構成で、レーダ装置１は、自船周囲の平面上の物標を、３６０°にわたり探知することができる。 The
なお、以下の説明では、パルス状電波を送波してから次のパルス状電波を送波するまでの動作を「スイープ」という。また、電波の送受信を行いながらアンテナを３６０°回転させる動作を「スキャン」と呼ぶ。 In the following description, the operation from the transmission of a pulsed radio wave to the transmission of the next pulsed radio wave is referred to as “sweep”. The operation of rotating the antenna by 360 ° while transmitting and receiving radio waves is called "scan".
受信部５は、アンテナ４で受波したエコーから得られるエコー信号を検波して増幅する。受信部５は、増幅したエコー信号を、Ａ／Ｄ変換部６へ出力する。Ａ／Ｄ変換部６は、アナログ形式のエコー信号をサンプリングし、複数ビットからなるデジタルデータに変換する。このデジタルデータは、エコーデータである。エコーデータは、アンテナ４が受波した反射波から得られたエコー信号の強度を特定するデータを含んでいる。Ａ／Ｄ変換部６は、エコーデータを、信号処理装置２０へ出力する。 The
図１を参照して、信号処理装置２０は、周波数領域スペクトル生成部８と、信頼性評価部９と、表示用データ生成部１０と、を有している。 With reference to FIG. 1, the
信号処理装置２０は、ハードウェア・プロセッサ７（たとえば、ＣＰＵ、ＦＰＧＡなど）および不揮発性メモリなどのデバイスで構成される。たとえば、ＣＰＵが不揮発性メモリからプログラムを読み出して実行することにより、信号処理装置２０を、周波数領域スペクトル生成部８、信頼性評価部９、および、表示用データ生成部１０として機能させることができる。 The
信号処理装置２０は、アンテナ４で受波されたエコーを処理する。 The
図５は、自船１００の周囲における波浪の一例を示す海上の模式的な平面図である。図５では、自船１００の周囲において自船１００に向かう２種類の波１０１，１０２が存在している様子を示している。図５では、たとえば、自船１００の右舷前方における波１０１が、自船１００の周囲で最も強い波である状態を示している。また、図５では、自船１００の右舷後方における波１０２が、波１０１の次に強い波である状態を示している。図５では、波１０１，１０２について、明瞭な波峰線ｗが表示されている。なお、自船１００の周囲には、波１０１，１０２以外にも、比較的波高が低い波浪が存在しているけれども、それらの波浪については図示を省略している。 FIG. 5 is a schematic plan view on the sea showing an example of waves around the
図１および図５を参照して、周波数領域スペクトル生成部８は、自船１００の周囲の所定の観測エリアとしての自船周囲領域Ｚ１における所定の解析領域Ｚ１０内に含まれる波浪から得られたエコーを周波数解析する。自船周囲領域Ｚ１は、レーダ装置１によって探知可能な領域である。周波数領域スペクトル生成部８は、解析領域Ｚ１０における方向周波数スペクトル（周波数領域スペクトル）を生成する。なお、本実施形態では、図５に示すように、解析領域Ｚ１０が自船の右舷前方に設定されている。 With reference to FIGS. 1 and 5, the frequency domain
周波数領域スペクトル生成部８は、データ蓄積部１２と、周波数解析部１３と、波浪情報取得部１４と、分散関係算出部１５と、ＭＴＦ処理部１６と、ωθ変換部１７と、２次元スペクトル変換部１８と、を有している。 The frequency domain
データ蓄積部１２は、自船周囲領域Ｚ１内に含まれるエコーデータを蓄積する。データ蓄積部１２は、１回のスキャンで得られる自船周囲領域Ｚ１における解析領域Ｚ１０内のエコー画像に含まれるエコーデータを、スキャン毎に抽出する。これにより、データ蓄積部１２は解析領域Ｚ１０について、複数スキャン分（たとえば３２枚分）のエコーデータを蓄積する。 The
周波数解析部１３は、本発明の「波数データ取得部」の一例である。周波数解析部１３は、解析領域Ｚ１０における波浪を観測することで得られた観測データ（すなわち、エコーデータ）に基づいて、波数と当該波数における周波数との関係を含む波数データＤ１を取得する。 The
本実施形態では、周波数解析部１３は、解析領域Ｚ１０について、複数スキャン分のエコーデータを用いて周波数解析する。具体的には、周波数解析部１３は、解析領域Ｚ１０について、複数スキャン分のエコーデータを用いて３次元の高速フーリエ変換（３ＤＦＦＴ、Three Dimensional Fast Fourier Transform）処理を行う。これにより、解析領域Ｚ１０について、図６に示すように、ｘ軸およびｙ軸の単位をそれぞれｒａｄ／ｍ、ｚ軸の単位をｒａｄ／ｓｅｃとした３次元データである波数データＤ１を生成する。 In the present embodiment, the
図６は、エコーデータを３ＤＦＦＴ処理した結果を示すグラフである。この３次元の波数データＤ１の処理結果を示すグラフにおけるｘ軸は、たとえば、東西方向における波数ｋｘであり、ｙ軸は、南北方向における波数ｋｙであり、ｚ軸は、角周波数ωである。３次元グラフを特定する波数データＤ１において、同一のエコーを示す結果が、周波数の正方向および負方向に対称的に表れる。 FIG. 6 is a graph showing the result of 3DFFT processing of echo data. The x-axis in the graph showing the processing result of the three-dimensional wave number data D1 is, for example, the wave number kx in the east-west direction, the y-axis is the wave number ky in the north-south direction, and the z-axis is the angular frequency ω. In the wave number data D1 that specifies the three-dimensional graph, the results showing the same echo appear symmetrically in the positive and negative directions of the frequency.
図１を参照して、周波数解析部１３は、得られた３次元の波数データＤ１を、波浪情報取得部１４および信頼性評価部９に出力する。波浪情報取得部１４は、本発明の「波浪情報取得部」の一例である。波浪情報取得部１４は、周波数解析部１３によって得られた３次元の波数データＤ１の中から、波浪に起因する成分である波浪成分を抽出することで、波浪データＤ２を取得する。具体的には、波浪情報取得部１４は、３次元の波数データＤ１に対して、分散関係算出部１５で算出された波浪の分散関係を用いるフィルタ処理を施す。これにより、波浪情報取得部１４は、波浪成分、すなわち、波数データＤ１に含まれる波浪成分データを検出する。分散関係算出部１５は、波浪成分に関する分散関係を算出する部分として設けられている。分散関係算出部１５は、たとえば、以下の（１）式で表わされる分散関係式を算出する。 Referring to FIG. 1,
ω２＝ｇｋ…（１）ω 2 = gk (1)
但し、ωは角周波数、ｋは波数、ｇは重力加速度である。この分散関係式の算出結果は、重力波を示す内容となる。 Here, ω is the angular frequency, k is the wave number, and g is the gravitational acceleration. The calculation result of this dispersion relational expression has the content indicating a gravitational wave.
波浪情報取得部１４は、周波数解析部１３によって得られた３次元データのうち、上記式（１）に近接するスペクトルのデータを抽出するというフィルタ処理を施すことにより、波浪成分を抽出する。これにより、周波数解析部１３は、図７に示す、２次元のスペクトルデータである波浪データを生成する。図７は、周波数解析部１３の処理結果を示すグラフである。この２次元データにおけるｘ軸およびｙ軸は、図６の３次元データの次元と同じである。 The wave
図７において一例として示すｋｘｋｙスペクトルは、座標の各地点を構成するサンプル点のスペクトル強度を、情報として有している。なお、図７から後述する図１２までの各図においては、各サンプル点のスペクトル強度を、ハッチングの濃さに対応させて図示している。すなわち、ハッチングが濃い部分を構成するサンプル点のエコー強度は、ハッチングが薄い部分を構成するサンプル点のエコー強度よりも高い。波浪情報取得部１４は、波浪成分を特定する波浪データＤ２を、信頼性評価部９およびＭＴＦ処理部１６に出力する。 The kxky spectrum shown as an example in FIG. 7 has, as information, the spectral intensities of sample points forming each point of coordinates. In each of FIGS. 7 to 12 to be described later, the spectral intensity of each sample point is shown in correspondence with the darkness of the hatching. That is, the echo intensity of the sample points forming the darkly hatched portion is higher than the echo intensity of the sample points forming the lightly hatched portion. The wave
ＭＴＦ処理部１６は、波数ｋに応じた所定の係数を、波浪データＤ２で特定される値に乗じることで、波浪データＤ３を生成する。 The
ωθ変換部１７は、波浪情報取得部１４によって抽出されＭＴＦ処理部１６で係数を乗じられた波浪データＤ３で特定されるスペクトルを、極座標（ωθ座標）に変換することで、ωθスペクトルを生成する。このωθ座標は、横軸が自船１００を基準とした波浪の方向θ（方位方向）に対応し且つ縦軸が角周波数ωに対応する座標である。ωθスペクトルは、波浪データＤ３で特定されるスペクトルの、極座標における、方位方向θと周波数方向ωの２次元スペクトルである。このωθスペクトルは、方位方向θのスペクトル成分としての方位スペクトルと、周波数方向ωのスペクトル成分としての周波数スペクトルと、を含んでいる。なお、ωθスペクトルは、スムージング処理を施されてもよい。 The
図８は、ωθ変換部１７によって生成されたωθスペクトルの一例を示す図である。図１、図５および図８を参照して、一例として示すωθスペクトルは、極座標の各地点を構成するサンプル点のスペクトル強度を、情報として有している。図８では、一例として、解析領域Ｚ１０に含まれる波１０１に対応する第１ピーク領域３１と、解析領域Ｚ１０に含まれる波１０２に対応する第２ピーク領域３２と、が示されている。本実施形態では、第１ピーク領域３１におけるスペクトルのピーク値Ｐ１は、第２ピーク領域３２におけるスペクトルのピーク値Ｐ２よりも大きい。 FIG. 8 is a diagram showing an example of the ωθ spectrum generated by the
ωθ変換部１７は、ωθスペクトルを特定するデータとしてのωθスペクトルデータＤ４を、２次元スペクトル変換部１８および信頼性評価部９へ出力する。 The
２次元スペクトル変換部１８は、ωθスペクトルをｘ−ｙ直交座標に変換することで、図９に示す波浪スペクトルを特定するデータとしての波浪スペクトルデータＤ５を生成する。図９は、波浪スペクトルの一例を示す模式的な図である。なお、この直交座標において、ｘ軸は、たとえば、東西方向を示し、ｙ軸は、たとえば、南北方向を示す。 The two-dimensional
ωｘ−ωｙ直交座標において、原点からの向きは、波向を示している。また、原点からの波浪スペクトルの距離は、波周期を示している。また、各サンプル点のスペクトル強度を、ハッチングの濃さに対応させて図示している。この波浪スペクトルは、解析領域Ｚ１０内に最も大きな波１０１が存在するとともに、解析領域Ｚ１０内に２番目に大きな波１０２が存在する状態を示している。２次元スペクトル変換部１８は、波浪スペクトルデータＤ５を、表示用データ生成部１０へ出力する。 In the ωx-ωy orthogonal coordinates, the direction from the origin indicates the wave direction. The distance of the wave spectrum from the origin indicates the wave period. In addition, the spectrum intensity of each sample point is shown in correspondence with the density of hatching. This wave spectrum shows a state in which the
次に、信頼性評価部９の構成を、より具体的に説明する。図１、図３および図４を参照して、信頼性評価部９は、波浪スペクトルで特定される波浪情報について、実際の波浪の状態を示しているか否かの信頼性（換言すれば、波浪情報の信頼度）を評価するために設けられている。 Next, the configuration of the reliability evaluation unit 9 will be described more specifically. With reference to FIG. 1, FIG. 3, and FIG. 4, the reliability evaluation unit 9 determines whether or not the wave information specified by the wave spectrum indicates the actual wave state (in other words, the wave state). It is provided to evaluate the reliability of information).
信頼性評価部９は、第１信頼性評価ユニット２１と、第２信頼性評価ユニット２２と、を有している。 The reliability evaluation section 9 includes a first
第１信頼性評価ユニット２１は、フィルタ前パワー算出部２３と、フィルタ後パワー算出部２４と、パワー比算出部２５と、指標判定部２６と、を有している。 The first
フィルタ前パワー算出部２３は、周波数解析部１３で生成された３次元データ（すなわち、波数データＤ１）で特定される３次元の各成分のパワーを、フィルタ前パワーＰＷ１として算出する。具体的には、３次元データで特定される３次元座標での配列（３次元配列）における全ての要素の和が、フィルタ前パワーＰＷ１として算出される。なお、３次元配列における全ての要素の「要素」とは、各波数、周波数におけるエコー強度を意味する。フィルタ前パワー算出部２３は、算出したフィルタ前パワーＰＷ１の値を、パワー比算出部２５へ出力する。 The pre-filter
フィルタ後パワー算出部２４は、波浪データＤ２で特定される２次元の各成分のパワーを、フィルタ後パワーＰＷ２として算出する。具体的には、波浪情報取得部１４のおける分散関係フィルタリング（換言すれば、重力波曲面上にある成分が抽出された）後の２次元データである、波浪成分データ（２次元配列）における全て要素の和を２倍した値を、フィルタ後パワーＰＷ２として算出する。なお、２次元配列における全ての要素の「要素」とは、各波数におけるエコー強度を意味する。 The post-filter
また、波浪成分データ（２次元配列）における全て要素の和を２倍している理由は、分散関係は周波数の正方向および負方向のうちの正方向のみをみているため、フィルタリング後の要素を２倍する必要があるからである。フィルタ後パワー算出部２４は、算出したフィルタ後パワーＰＷ２の値を、パワー比算出部２５へ出力する。 The reason for doubling the sum of all the elements in the wave component data (two-dimensional array) is that the dispersion relationship is only seen in the positive direction of the positive and negative directions of the frequency. This is because it needs to be doubled. The post-filter
パワー比算出部２５は、フィルタ後パワーＰＷ２とフィルタ前パワーＰＷ１との比ＰＷ２／ＰＷ１を算出し、算出したパワー比ＰＷ２／ＰＷ１の値を、パワー比指標判定部２６に出力する。 The
パワー比指標判定部２６は、パワー比ＰＷ２／ＰＷ１が所定のしきい値Ｐｗｔｈ以上であるか否かを判定する。しきい値Ｐｗｔｈは、所定値に設定されている。なお、エコーデータに含まれる波浪成分が多いほど、パワー比ＰＷ２／ＰＷ１は大きくなる。一方、エコーデータにおいて、波浪成分以外の成分（たとえば、陸地からのエコーによる成分）が大きければ、パワー比ＰＷ２／ＰＷ１は小さくなる。 The power ratio
つまり、パワー比ＰＷ２／ＰＷ１が大きければ大きいほど、エコーデータに含まれる波浪成分が多いこととなり、エコーデータについて、実際の波浪をより正確に示すデータであるといえる。パワー比指標判定部２６は、パワー比ＰＷ２／ＰＷ１がしきい値Ｐｗｔｈ以上であれば、エコーデータの信頼性が高いと判定する。一方、パワー比指標判定部２６は、パワー比ＰＷ２／ＰＷ１がしきい値Ｐｗｔｈ未満であれば、エコーデータの信頼性が低いと判定する。パワー比指標判定部２６は、判定結果を示すデータとしてのパワー比判定結果データＤ６を、２次元スペクトル変換部１８へ出力する。 That is, the larger the power ratio PW2 / PW1, the more wave components are included in the echo data, and it can be said that the echo data is data that more accurately indicates the actual wave. The power ratio
次に、第２信頼性評価ユニット２２の構成を、より具体的に説明する。 Next, the configuration of the second
図１、図４、図５および図８を参照して、信頼性評価部９は、後述する第１評価用エリアＡＲ１０および第２評価用エリアＡＲ２０について、波浪データＤ３で特定される波浪情報の信頼性を評価する。 With reference to FIGS. 1, 4, 5, and 8, the reliability evaluation unit 9 determines the wave information specified by the wave data D3 for a first evaluation area AR10 and a second evaluation area AR20 described later. Evaluate reliability.
信頼性評価部９は、第１評価用スペクトル算出部４１と、第２評価用スペクトル算出部４２と、第１スペクトル評価部４３と、第２スペクトル評価部４４と、ピーク比算出部４５と、ピーク比評価部４６と、を有している。 The reliability evaluation unit 9 includes a first evaluation
第１評価用スペクトル算出部４１は、ωθスペクトルにおけるピーク値が最大である領域を含む領域としての第１評価用エリアＡＲ１０に関して、波浪データＤ２に基づいて、波浪に関するスペクトルである第１評価用スペクトルＳＰ１０の算出処理を行う。 The first evaluation
第１評価用スペクトル算出部４１は、第１ピーク検出部４７と、第１評価用エリア設定部４８と、第１評価用スペクトル抽出部４９と、を有している。 The first evaluation
第１ピーク検出部４７は、解析領域Ｚ１０に含まれるエコー画像から得られたωθスペクトルのうち、ピーク値（エコー強度）が最も高い点としての第１ピーク点ＰＴ１を検出する。なお、ωθスペクトルは、解析領域Ｚ１０内に含まれる波浪から得られたエコーに基づいて算出された、波浪の角周波数と波浪の到来方向とで特定されるωθ座標の各位置に該波浪のエコー強度が表示されたωθスペクトルであるといえる。本実施形態では、第１ピーク検出部４７は、たとえば、図８のうち、色の最も濃い点、すなわち、ピーク値Ｐ１を有する点としての第１ピーク点Ｐ１を検出する。 The first
図１０は、第１評価用エリアＡＲ１０について説明するための概念図である。図４および図１０を参照して、第１評価用エリア設定部４８は、第１評価用エリアＡＲ１０を設定する。第１評価用エリア設定部４８は、第１ピーク点ＰＴ１を含む、所定の第１評価用エリアＡＲ１０を設定する。第１評価用エリアＡＲ１０は、第１ピーク点ＰＴ１を中心として、方位方向θおよび周波数方向ωのそれぞれに、所定の範囲に亘って設定される。なお、本明細書では、第１評価用エリアＡＲ１０および第２評価用エリアＡＲ２０は、方位方向θおよび周波数方向ωのそれぞれについて、誇張して示している。 FIG. 10 is a conceptual diagram for explaining the first evaluation area AR10. Referring to FIGS. 4 and 10, the first evaluation
なお、第１ピーク点ＰＴ１がθ＝０ｄｅｇ付近、または、３６０ｄｅｇ付近に存在する場合、第１評価用エリア設定部４８は、ωθスペクトルのθ＝０ｄｅｇの部分とθ＝３６０ｄｅｇの部分とが連続していると扱うことで、ピーク値が不連続に変化する不自然な第１評価用エリアＡＲ１０の生成を抑制できる。同様に、周波数方向ωにおいて、第１ピーク点ＰＴ１が下限（たとえば、０ｒａｄ／ｓ）または上限付近に存在する場合、第１評価用エリア設定部４８は、ωθスペクトルのうち、周波数方向ωにおいて、ωθスペクトルの外側の領域について、ωθスペクトルにおける最小の信号強度と同じ強度のスペクトルを補完する。これにより、ピーク値が不連続に変化する不自然な第１評価用エリアＡＲ１０の生成を抑制できる。 When the first peak point PT1 is near θ = 0 deg or around 360 deg, the first evaluation
第１評価用スペクトル抽出部４９は、第１評価用スペクトルＳＰ１０を抽出する。第１評価用スペクトルＳＰ１０は、第１評価用エリアＡＲ１０におけるスペクトルで特定される波浪情報が実際の波浪状態を正確に示しているか否かについての、信頼性を評価するために用いられるスペクトルである。第１評価用スペクトルＳＰ１０は、第１評価用エリアＡＲ１０において最も大きいピーク値Ｐ１を有するスペクトルである。 The first evaluation
第１評価用スペクトル抽出部４９が抽出する第１評価用スペクトルＳＰ１０の成分は、第１方位スペクトルＳＰ１１と、第１周波数スペクトルＳＰ１２と、を含んでいる。第１方位スペクトルＳＰ１１は、第１評価用エリアＡＲ１０における、方位方向θに沿う強度のスペクトルを示す。第１方位スペクトルＳＰ１１の抽出方法として、２種類例示することができる。より具体的には、（１）積分法と、（２）中心抽出法と、を例示することができる。 The components of the first evaluation spectrum SP10 extracted by the first evaluation
上記（１）の積分法では、方位方向θの任意の１点において、方位方向θの座標が同じスペクトルの強度を周波数方向ωに積分した積分値が、当該１点における第１方位スペクトルＳＰ１の強度となる。そして、この積分値を方位方向θに沿ってプロットする処理を行うことで、第１方位スペクトルＳＰ１１を抽出することができる。なお、図１０では、一例としてこの積分法で抽出された第１方位スペクトルＳＰ１１を示している。 In the integration method of (1) above, at any one point in the azimuth direction θ, the integrated value obtained by integrating the intensity of the spectrum having the same coordinate in the azimuth direction θ in the frequency direction ω is the first azimuth spectrum SP1 at the one point. It becomes strength. Then, the first azimuth spectrum SP11 can be extracted by performing a process of plotting this integrated value along the azimuth direction θ. Note that FIG. 10 shows the first azimuth spectrum SP11 extracted by this integration method as an example.
上記（２）の中心抽出法では、ωθスペクトルにおいて、第１ピーク点ＰＴ１を通過し且つ方位方向θに延びる基準線Ｌ１上での方位θと強度との関係を示すスペクトル（すなわち、１次元配列スペクトル）が、第１方位スペクトルＳＰ１１として抽出される。 In the center extraction method of (2) above, in the ωθ spectrum, a spectrum showing the relationship between the azimuth θ and the intensity on the reference line L1 that passes through the first peak point PT1 and extends in the azimuth direction θ (that is, a one-dimensional array). Spectrum) is extracted as the first azimuth spectrum SP11.
第１周波数スペクトルＳＰ１２は、第１評価用エリアＡＲ１０における、周波数方向ωに沿う強度のスペクトルを示す。第１周波数スペクトルＳＰ１２の抽出方法として、上記した２種類を例示することができる。 The first frequency spectrum SP12 shows a spectrum of intensity along the frequency direction ω in the first evaluation area AR10. As the extraction method of the first frequency spectrum SP12, the above two types can be exemplified.
上記（１）の積分法では、周波数方向ωの任意の１点において、周波数方向ωの座標が同じスペクトルの強度を方位方向θに積分した積分値が、当該１点における第１周波数スペクトルＳＰ１２の強度となる。そして、この積分値を周波数方向ωに沿ってプロットする処理を行うことで、第１周波数スペクトルＳＰ１２を抽出することができる。なお、図１０では、一例としてこの積分法で抽出された第１周波数スペクトルＳＰ１２を示している。 In the integration method of (1) above, at any one point in the frequency direction ω, the integrated value obtained by integrating the intensity of the spectrum having the same coordinate in the frequency direction ω in the azimuth direction θ is the first frequency spectrum SP12 at the one point. It becomes strength. Then, the first frequency spectrum SP12 can be extracted by performing a process of plotting this integrated value along the frequency direction ω. Note that FIG. 10 shows the first frequency spectrum SP12 extracted by this integration method as an example.
上記（２）の中心抽出法では、ωθスペクトルにおいて、第１ピーク点ＰＴ１を通過し且つ周波数方向ωに延びる基準線Ｌ２上での周波数ωと強度との関係を示すスペクトル（すなわち、１次元配列スペクトル）が、第１周波数スペクトルＳＰ１２として抽出される。 In the center extraction method of (2) above, in the ωθ spectrum, a spectrum showing the relationship between the frequency ω and the intensity on the reference line L2 that passes through the first peak point PT1 and extends in the frequency direction ω (that is, a one-dimensional array). Spectrum) is extracted as the first frequency spectrum SP12.
図４および図８を参照して、第２評価用スペクトル算出部４２は、図８に示すωθ座標に含まれる領域である評価対象エリアＡＲ１のうち、第１評価用エリアＡＲ１０を除くエリアについて、ωθスペクトルにおけるピーク値（エコー強度）が最も大きい領域を含む領域としての第２評価用エリアＡＲ２０に関して、波浪データＤ２に基づいて波浪に関するスペクトルである第２評価用スペクトルＳＰ２０の算出処理を行う。 With reference to FIG. 4 and FIG. 8, the second evaluation
第２評価用スペクトル算出部４２は、第１評価用エリア置換部５６と、第２ピーク検出部５７と、第２評価用エリア設定部５８と、第２評価用スペクトル抽出部５９と、を有している。 The second evaluation
第１評価用エリア置換部５６は、ωθスペクトルのうち、第１評価用エリアＡＲ１０の各座標における第１評価用スペクトルＳＰ１０の強さを低減させる処理（すなわち、置換処理）を行う。具体的には、第１評価用エリア置換部５６は、ωθスペクトルのうち、第１評価用エリアＡＲ１０の各座標における信号強度を、ωθスペクトルにおける最小強度またはゼロに置換する。これにより、図１１に示すように、第１ピーク点ＰＴ１を含む第１評価用エリアＡＲ１０の信号強度は、ノイズ程度の小さな値となる。図１１は、ωθスペクトルにおける第１評価用エリアＡＲ１０が置換処理された状態を示す模式図である。 The first evaluation
再び図４および図８を参照して、第２ピーク検出部５７は、評価対象エリアＡＲ１のうち、第１評価用エリアＡＲ１０を除くエリアにおいて、ωθスペクトルのうち、ピーク値（エコー強度）が最も高い点としての第２ピーク点ＰＴ２を抽出する。本実施形態では、第２ピーク検出部５７は、たとえば、図１１におけるωθスペクトルのうち、色の最も濃い点を検出する。 Referring again to FIGS. 4 and 8, the second
図１２は、第２評価エリアＡＲ２０について説明するための概念図である。図４および図１２を参照して、第２評価用エリア設定部５８は、図８に示すωθスペクトルにおいて、第２ピーク点ＰＴ２を含む、所定の第２評価用エリアＡＲ２０を設定する。第２評価用エリア設定部５８による第２評価エリアＡＲ２０の設定方法は、第１評価用エリア設定部４８による第１評価用エリアＡＲ１０の設定方法と同様である。 FIG. 12 is a conceptual diagram for explaining the second evaluation area AR20. Referring to FIGS. 4 and 12, the second evaluation
具体的には、第２評価用エリアＡＲ２０は、評価対象エリアＡＲ１のうち、第１評価用エリアＡＲ１０を除くエリアにおいて最も強いエコー強度が観測された第２ピーク点ＰＴ２を中心として、方位方向θおよび周波数方向ωのそれぞれに、所定の範囲に亘って設定される。 Specifically, the second evaluation area AR20 has the azimuth direction θ around the second peak point PT2 where the strongest echo intensity is observed in the areas other than the first evaluation area AR10 in the evaluation target area AR1. And in the frequency direction ω, respectively, over a predetermined range.
前述したように、第２評価エリアＡＲ２０の設定に用いられるωθスペクトルは、第１評価用エリアＡＲ１０が置換処理される前の当該ωθスペクトルであり、図８に示すωθスペクトルである。なお、第１評価用エリアＡＲ１０が置換処理された後の当該ωθスペクトル（図１１に示すωθスペクトル）を基に第２評価用エリアＡＲ２０が設定される場合、第１評価用エリアＡＲ１０と第２評価用エリアとの重複部分において、強度に不自然なエッジが生じてしまうことが考えられる。よって、第２評価エリアＡＲ２０の設定に用いられるωθスペクトルを、第１評価用エリアＡＲ１０が置換処理される前の当該ωθスペクトル（図８に示すωθスペクトル）とすることで、このような不自然なエッジを抑制できる。 As described above, the ωθ spectrum used for setting the second evaluation area AR20 is the ωθ spectrum before the first evaluation area AR10 is replaced, and is the ωθ spectrum shown in FIG. 8. When the second evaluation area AR20 is set on the basis of the ωθ spectrum (ωθ spectrum shown in FIG. 11) after the first evaluation area AR10 has been replaced, the first evaluation area AR10 and the second evaluation area AR10 are set. It is conceivable that an unnatural edge is generated in the strength in the overlapping portion with the evaluation area. Therefore, by using the ωθ spectrum used for setting the second evaluation area AR20 as the ωθ spectrum before the first evaluation area AR10 is subjected to the replacement process (ωθ spectrum shown in FIG. 8), such an unnatural It is possible to suppress various edges.
また、第２ピーク点ＰＴ２がθ＝０ｄｅｇ付近、または、θ＝３６０ｄｅｇに存在する場合、第２評価用エリア設定部５８は、ωθスペクトルのθ＝０ｄｅｇの部分とθ＝３６０ｄｅｇの部分とが連続していると扱うことで、ピーク値が不連続に変化する不自然な第２評価用エリアＡＲ２０の生成を抑制できる。同様に、周波数方向ωにおいて、第２ピーク点ＰＴ２が下限（たとえば、０ｒａｄ／ｓ）または上限付近に存在する場合、第２評価用エリア設定部５８は、ωθスペクトルのうち、周波数方向ωにおいて、ωθスペクトルの外側の領域について、ωθスペクトルにおける最小の信号強度と同じ強度のスペクトルを補完する。これにより、ピーク値が不連続に変化する不自然な第２評価用エリアＡＲ２０の生成を抑制できる。 When the second peak point PT2 exists near θ = 0 deg or at θ = 360 deg, the second evaluation
第２評価用スペクトル抽出部５９は、第２評価用スペクトルＳＰ２０を抽出する。第２評価用スペクトルＳＰ２０は、第２評価用エリアＡＲ２０におけるスペクトルで特定される波浪情報が実際の波浪状態を正確に示しているか否かについての、信頼性を評価するために用いられるスペクトルである。第２評価用スペクトルＳＰ２０は、第２評価用エリアＡＲ２０の中央において最も大きいピーク値Ｐ２を有するスペクトルである。 The second evaluation
第２評価用スペクトル抽出部５９が抽出する第２評価用スペクトルＳＰ２０の成分は、第２方位スペクトルＳＰ２１と、第２周波数スペクトルＳＰ２２と、を含んでいる。第２方位スペクトルＳＰ２１は、第２評価用エリアＡＲ２０における、方位方向θに沿う強度のスペクトルを示す。第２方位スペクトルＳＰ２１の抽出方法として、２種類例示することができる。より具体的には、前述した（１）積分法と、（２）中心抽出法と、を例示することができる。 The components of the second evaluation spectrum SP20 extracted by the second evaluation
よって、第２評価用スペクトルＳＰ２０についての積分法、および、中心抽出法についての詳細な説明は、省略する。なお、図１２では、一例として積分法で抽出された第２方位スペクトルＳＰ２１を示している。 Therefore, detailed description of the integration method and the center extraction method for the second evaluation spectrum SP20 is omitted. Note that FIG. 12 shows the second azimuth spectrum SP21 extracted by the integration method as an example.
第２周波数スペクトルＳＰ２２は、第２評価用エリアＡＲ２０における、周波数方向ωに沿う強度のスペクトルを示す。第２周波数スペクトルＳＰ２２の抽出方法として、上記した２種類を例示することができる。なお、図１２では、一例として積分法で抽出された第２周波数スペクトルＳＰ２２を示している。 The second frequency spectrum SP22 shows a spectrum of intensity along the frequency direction ω in the second evaluation area AR20. As the method of extracting the second frequency spectrum SP22, the above two types can be exemplified. Note that FIG. 12 shows the second frequency spectrum SP22 extracted by the integration method as an example.
図４および図１０を参照して、第１スペクトル評価部４３は、周波数解析によって生成されたωθスペクトルにおける第１評価用エリアＡＲ１０で特定される波浪情報が、実際の波浪状態を示しているか否かの信頼性を、第１評価用スペクトルＳＰ１０を用いて判定する。 With reference to FIGS. 4 and 10, the first
具体的には、第１スペクトル評価部４３は、第１方位スペクトルＳＰ１１と第１周波数スペクトルＳＰ１２のそれぞれについて、強度の平均値ＰＡ１１，ＰＡ１２を算出する。そして、第１スペクトル評価部４３は、第１方位スペクトルＳＰ１１の基準値としての中央ピーク値Ｐ１１を平均値ＰＡ１１で除した値Ｐ１１／ＰＡ１１を、信頼性判定指標Ｂ１１として算出する。同様に、第１スペクトル評価部４３は、第１周波数スペクトルＳＰ１２の中央ピーク値Ｐ１２を平均値ＰＡ１２で除した値Ｐ１２／ＰＡ１２を、信頼性判定指標Ｂ１２として算出する。なお、中央ピーク値とは、対応する方位方向θまたは周波数方向ωにおける中央位置での信号強度をいう。 Specifically, the first
第１スペクトル評価部４３は、各信頼性判定指標Ｂ１１，Ｂ１２が、所定のしきい値Ｂ１ｔｈ以上であるか否かを判定する。これは第１方位スペクトルＳＰ１１および第１周波数スペクトルＳＰ１２のそれぞれにおいて、対応する信頼性判定指標Ｂ１１，Ｂ１２が大きいと、スペクトルのピーク形状が鋭いこととなり、波浪成分が明確であることにより信頼性が高いと判定できるからである。 The 1st
第１スペクトル評価部４３は、各信頼性判定指標Ｂ１１，Ｂ１２が何れもしきい値Ｂ１ｔｈ以上である場合、第１評価用スペクトルＳＰ１０で特定される波浪情報の信頼性が高いと判定する。一方、第１スペクトル評価部４３は、各信頼性判定指標Ｂ１１，Ｂ１２の少なくとも一方がしきい値Ｂ１ｔｈ未満である場合、第１評価用スペクトルＳＰ１０で特定される波浪情報の信頼性が低いと判定する。 The first
図４および図１２を参照して、第２スペクトル評価部４４は、周波数解析によって生成されたωθスペクトルにおける第２評価用エリアＡＲ２０で特定される波浪情報が、実際の波浪状態を示しているか否かの信頼性を、第２評価用スペクトルＳＰ２０を用いて判定する。本実施形態では、第２スペクトル評価部４４における処理は、第１スペクトル評価部４３における処理と同様である。 With reference to FIGS. 4 and 12, the second
具体的には、第２スペクトル評価部４４は、第２方位スペクトルＳＰ２１と第２周波数スペクトルＳＰ２２のそれぞれについて、強度の平均値ＰＡ２１，ＰＡ２２を算出する。そして、第２スペクトル評価部４４は、第２方位スペクトルＳＰ２１の基準値としての中央ピーク値Ｐ２１を平均値ＰＡ１２で除した値Ｐ２１／ＰＡ２１を、信頼性判定指標Ｂ２１として算出する。同様に、第２スペクトル評価部４４は、第２周波数スペクトルＳＰ２２の中央ピーク値Ｐ２２を平均値ＰＡ２２で除した値Ｐ２２／ＰＡ２２を、信頼性判定指標Ｂ２２として算出する。 Specifically, the second
第２スペクトル評価部４４は、各信頼性判定指標Ｂ２１，Ｂ２２が、所定のしきい値Ｂ２ｔｈ（たとえば、１．５）以上であるか否かを判定する。第２スペクトル評価部４４は、各信頼性判定指標Ｂ２１，Ｂ２２が何れもしきい値Ｂ２ｔｈ以上である場合、第２評価用スペクトルＳＰ２０で特定される波浪情報の信頼性が高いと判定する。一方、第２スペクトル評価部４４は、各信頼性判定指標Ｂ２１，Ｂ２２の少なくとも一方がしきい値Ｂ２ｔｈ未満である場合、第２評価用スペクトルＳＰ２０で特定される波浪情報の信頼性が低いと判定する。 The second
図４、図１０および図１２を参照して、ピーク比算出部４５は、第１評価用スペクトルＳＰ１０と第２評価用スペクトルＳＰ２０の比を算出する。具体的には、ピーク比算出部４５は、第２方位スペクトルＳＰ２１の中央ピーク値Ｐ２１を第１方位スペクトルＳＰ１１の中央ピーク値Ｐ１１で除した値Ｐ２１／Ｐ１１を、信頼性判定指標Ｂ３として算出する。同様に、ピーク比算出部４５は、第２周波数スペクトルＳＰ２２の中央ピーク値Ｐ２２を第１周波数スペクトルＳＰ１２の中央ピーク値Ｐ１２で除した値Ｐ２２／Ｐ１２を、信頼性判定指標Ｂ４として算出する。 With reference to FIGS. 4, 10 and 12, the peak
ピーク比評価部４６は、各信頼性判定指標Ｂ３，Ｂ４が、所定のしきい値Ｂ３４ｔｈ以上であるか否かを判定する。ここで、第２評価用スペクトルＳＰ２０の中央ピーク値Ｐ２１，Ｐ２２が第１評価用スペクトルＳＰ１０の対応するピーク値Ｐ１１，Ｐ１２よりも十分に小さければ、第２評価用スペクトルＳＰ２０で特定される波浪は、第１評価用スペクトルＳＰ１０で特定される波浪に対してノイズ程度の存在となり、信頼度が低いか、または、無視できると考えることができる。 The peak
ピーク比評価部４６は、各信頼性判定指標Ｂ３，Ｂ４が何れもしきい値Ｂ３４ｔｈ以上である場合、第２評価用スペクトルＳＰ２０で特定される波浪情報の信頼性が高いと判定する。一方、ピーク比評価部４６は、各信頼性判定指標Ｂ３，Ｂ４の少なくとも一方がしきい値Ｂ３４ｔｈ未満である場合、第２評価用スペクトルＳＰ２０で特定される波浪情報の信頼性が低いと判定する。 The peak
なお、信頼性評価部９の第１スペクトル評価部４３、第２スペクトル評価部４４、および、ピーク比評価部４６のぞれぞれの評価結果を示すデータＤ７，Ｄ８，Ｄ９は、図１に示す２次元スペクトル変換部１８に出力される。２次元スペクトル変換部１８は、各評価用エリアＡＲ１０，ＡＲ２０の座標について、ωθ座標からωｘ−ωｙ座標に変換する。２次元スペクトル変換部１８で処理されたデータは、表示用データ生成部１０へ出力される。なお、本実施形態では、波周期，波向など、波浪情報のうちの一部のみが出力される構成となっている。よって、２次元スペクトル変換部１８を省略し、ωθ変換部１７から出力されるデータＤ４、および、信頼性評価部９から出力されるデータＤ６〜Ｄ９は、直接表示用データ生成部１０に出力されてもよい。これにより、信号処理装置２０の構成をより簡素にすることが可能である。 The data D7, D8, and D9 indicating the evaluation results of the first
表示用データ生成部１０は、２次元スペクトル変換部１８から出力された２次元の波浪スペクトルデータＤ５と、信頼性評価データＤ６〜Ｄ９と、を用いて、表示器３に表示させる画像のためのデータを生成する。 The display
図１３は、表示器３の表示画面の一例を示す図である。図１、図５および図１３を参照して、表示器３は、たとえば、液晶ディスプレイ装置であり、表示用データに基づいて画像を表示する。図１３では、一例として、波１０１の検出結果についての信頼性が高く、波１０２の検出結果についての信頼性が低いと判定された場合を示している。 FIG. 13 is a diagram showing an example of a display screen of the
表示用データ生成部１０は、たとえば、波浪スペクトルデータＤ５で特定されるスペクトルのうち、第１評価用エリアＡＲ１０に相当する領域の座標および信号強度などから、自船周囲領域ＡＲ１における最も強い波１０１の波浪情報を検出する。この場合の波浪情報は、たとえば、波高、方位、および、波長である。さらに、表示用データ生成部１０は、第１スペクトル評価部４３の出力結果に基づいて、当該波１０１の検出結果が十分な信頼性を有しているか否かを示すデータを生成する。 For example, the display
表示用データ生成部１０は、たとえば、波１０１の波浪情報が十分な信頼性を有している場合、波１０１について、波浪情報を表示させることを許可するデータを生成する。一方、表示用データ生成部１０は、たとえば、波１０１の波浪情報が十分な信頼性を有していない場合、波１０１について、波浪情報を表示することを許可しないデータを生成する。 For example, when the wave information of the
上記と同様に、表示用データ生成部１０は、たとえば、波浪スペクトルデータＤ５で特定されるスペクトルのうち、第２評価用エリアＡＲ２０に相当する領域の座標および信号強度などから、自船周囲領域ＡＲ１において２番目に強い波１０２の波浪情報を検出する。この場合の波浪情報は、たとえば、波高、方位、および、波長である。さらに、表示用データ生成部１０は、第２スペクトル評価部４４の出力結果に基づいて、当該波１０２の波浪情報が十分な信頼性を有しているか否かを示すデータを生成する。 Similar to the above, for example, the display
表示用データ生成部１０は、たとえば、波１０２の波浪情報が十分な信頼性を有している場合、波１０２について、波浪情報を表示させることを許可するデータを生成する。一方、表示用データ生成部１０は、たとえば、波１０２の波浪情報が十分な信頼性を有していない場合、波１０２について、波浪情報を表示することを許可しないデータを生成する。 For example, when the wave information of the
上記の処理によって生成された表示用データは、表示用データ生成部１０から表示器３へ出力される。 The display data generated by the above processing is output from the display
表示器３は、たとえば、波１０１，１０２のそれぞれについて、表示用データによって特定される波浪情報を表示可能である。表示器３は、第１波用表示領域６１と、第２波用表示領域６２とを有している。 The
波１０１に関する波浪情報の信頼性が高い場合、図１３に示すように、第１波用表示領域６１には、波１０１に関する波浪情報（本実施形態では、波高、方位、波長）がそれぞれ表示される。一方、波１０１に関する波浪情報の信頼性が低い場合、第１波用表示領域６１には、波１０１に関する波浪情報は表示されず、たとえば、横棒が表示される。 When the reliability of the wave information about the
上記と同様に、波１０２に関する波浪情報の信頼性が高い場合、第２波用表示領域６２には、波１０２に関する波浪情報（本実施形態では、波高、方位、波長）がそれぞれ表示される。一方、波１０２に関する波浪情報の信頼性が低い場合、図１３に示すように、第２波用表示領域６２には、波１０２に関する波浪情報は表示されず、たとえば、横棒が表示される。 Similarly to the above, when the reliability of the wave information about the
次に、信号処理装置２０における処理の流れの一例を説明する。図１４は、信号処理装置２０における処理の流れの一例を説明するためのフローチャートである。なお、以下では、フローチャートを参照して説明する場合、フローチャート以外の図も適宜参照しながら説明する。 Next, an example of the flow of processing in the
信号処理装置２０では、まず、データ蓄積部１２が、受信部５から出力されるエコーデータを、複数スキャン（本実施形態では、３２スキャン）分、蓄積する（ステップＳ１１）。 In the
次に、周波数解析部１３が、複数スキャン分のエコーデータを用いて、周波数解析をすることで、波数データＤ１を取得する（ステップＳ１２）。すなわち、所定の観測エリアとしての自船周囲領域ＡＲ１のうち解析領域Ｚ１０における波浪を観測することで得られた観測データに基づいて、波数と周波数との関係を含む３次元の波数データを取得する。 Next, the
次に、波浪情報取得部１４は、３次元の波数データＤ１に対して、分散関係算出部１５で算出された波浪の分散関係を用いるフィルタ処理を施す。これにより、波浪情報取得部１４は、波浪成分、すなわち、波数データＤ１に含まれる波浪データＤ２を検出する（ステップＳ１３）。 Next, the wave
次に、ＭＴＦ処理部１６は、波浪データＤ２に所定の係数を乗じる処理（ステップＳ１４）を行うことで、波浪データＤ３を生成する。 Next, the
その後、ωθ変換部１７は、波浪データＤ３の座標をωθの極座標に変換する（ステップＳ１５）。 After that, the
その後、信頼性評価部９は、周波数解析部１３から出力された波数データＤ１と、波浪データＤ２と、ωθ変換部１７から出力されたωθスペクトルデータＤ４と、を用いて、波浪データＤ３の信頼性を評価する（ステップＳ１６）。この際、信頼性評価部９は、ωθスペクトルデータＤ４を用いて、評価対象エリアＡＲ１の一部である第１評価用エリアＡＲ１０および、第２評価用エリアＡＲ２０について、波浪データＤ３で特定される波浪情報の信頼性を評価する。 After that, the reliability evaluation unit 9 uses the wave number data D1 output from the
２次元スペクトル変換部１８は、ωθスペクトルから出力されたωθスペクトルデータＤ４と、信頼性評価部９から出力された、各評価用エリアＡＲ１０，ＡＲ２０の信頼性評価結果（本実施形態では、評価データＤ６〜Ｄ９）と、について、ｋｙｋｙの２次元座標のデータに変換する（ステップＳ１７）。 The two-dimensional
次に、表示用データ生成部１０は、２次元スペクトル変換部１８から出力されたωθスペクトルデータＤ５を用いて、最も大きな波１０１と２番目に大きな波１０２のそれぞれの波向、波高、周波数について、信頼性評価に合格した波の情報を表示するためのデータを生成する。このデータが用いられることで、表示器３に信頼性の高い波浪データが表示される。 Next, the display
次に、信号処理装置２０の信頼性評価部９における処理（ステップＳ１６）の流れの一例を説明する。図１５は、信号処理装置２０の信頼性評価部９における処理の流れの一例を説明するためのフローチャートである。 Next, an example of the flow of processing (step S16) in the reliability evaluation unit 9 of the
信頼性評価部９では、まず、第１信頼性評価ユニット２１において、パワー比信頼性を判定する。第１信頼性評価ユニット２１の指標判定部２６は、パワー比ＰＷ２／ＰＷが前述のしきい値Ｐｗｔｈ以上である場合、信号処理装置２０において十分な大きさの波浪成分が検出されていると判定し（ステップＳ２１でＯＫ）、信頼性評価部９による信頼性評価（ステップＳ２２〜Ｓ２８）を許可する。 In the reliability evaluation unit 9, first, the first
一方、第１信頼性評価ユニット２１の指標判定部２６は、パワー比ＰＷ２／ＰＷが前述のしきい値Ｐｗｔｈ未満である場合、信号処理装置２０において十分な大きさの波浪成分が検出されていないと判定する（ステップＳ２１でＮＧ）。この場合、第１信頼性評価ユニット２１の指標判定部２６は、波１０１，１０２のそれぞれについての波浪情報の信頼性が低いことを示すデータＤ６を、２次元スペクトル変換部１８に出力する（ステップＳ２３）。 On the other hand, when the power ratio PW2 / PW is less than the threshold value Pwth, the
第１信頼性評価ユニット２１が、第２信頼性評価ユニット２２による信頼性評価を許可した場合（ステップＳ２１でＯＫ）、第２信頼性評価ユニット２２の第１スペクトル評価部４３は、第１評価用スペクトルＳＰ１０が信頼性を有する値を示しているか否かを判定する（ステップＳ２２）。 When the first
第１方位スペクトルＳＰ１１および第１周波数スペクトルＳＰ１２の少なくとも一方において、対応する判定指標Ｂ１１，Ｂ１２がしきい値Ｂ１ｔｈ未満である場合（ステップＳ２２でＮＧ）、第１スペクトル評価部４３は、第１評価用エリアＡＲ１０における波１０１、すなわち、自船周囲領域ＡＲ１において最も強い波１０１に関して、波浪情報の信頼度が低いと判定する（ステップＳ２３）。 In at least one of the first azimuth spectrum SP11 and the first frequency spectrum SP12, when the corresponding determination indicators B11 and B12 are less than the threshold value B1th (NG in step S22), the first
なお、波１０１に関して波浪情報の信頼度が低い場合、波１０２についての波浪情報の信頼度も当然に低いと考えられる。よって、この場合、第１スペクトル評価部４３は、第２評価用スペクトルＳＰ２０で特定される波１０２についても、波浪情報の信頼度が低いと判定する（ステップＳ２３）。この場合、第１スペクトル評価部４３は、この判定結果を示すデータＤ７を２次元スペクトル変換部１８へ出力する。 When the reliability of the wave information on the
一方、第１評価用スペクトルＳＰ１０の第１方位スペクトルＳＰ１１および第１周波数スペクトルＳＰ１２の何れにおいても、対応する判定指標Ｂ１１，Ｂ１２がしきい値Ｂ１ｔｈ以上である場合（ステップＳ２２でＯＫ）、第１スペクトル評価部４３は、第１評価用エリアＡＲ１０における波１０１、すなわち、自船周囲領域ＡＲ１の解析エリアＺ１０において最も強い波１０１に関して、波浪情報の信頼度が高いと判定する（ステップＳ２４）。この場合、第１スペクトル評価部４３は、この判定結果を示すデータＤ７を２次元スペクトル変換部１８へ出力する。 On the other hand, in any of the first azimuth spectrum SP11 and the first frequency spectrum SP12 of the first evaluation spectrum SP10, when the corresponding determination indicators B11 and B12 are equal to or more than the threshold value B1th (OK in step S22), the first The
次に、波１０１に関する波浪情報の信頼性が高いと判定された場合、ピーク比評価部４６は、ピーク比Ｐ２１／Ｐ１１，Ｐ２２／Ｐ１２のそれぞれについての信頼性判定指標Ｂ３，Ｂ４がしきい値Ｂ３４ｔｈ以上であるか否かを判定する（ステップＳ２５）。 Next, when it is determined that the reliability of the wave information regarding the
信頼性判定指標Ｂ３，Ｂ４の少なくとも一方がしきい値Ｂ３４ｔｈ未満である場合（ステップＳ２５でＮＧ）、ピーク比評価部４６は、第２評価用スペクトルＳＰ２０で特定される波１０２の波浪情報について、信頼度が低いと判定する（ステップＳ２７）。ピーク比評価部４６は、この判定結果を示すデータＤ９を２次元スペクトル変換部１８へ出力する。 When at least one of the reliability determination indexes B3 and B4 is less than the threshold value B34th (NG in step S25), the peak
一方、信頼性判定指標Ｂ３，Ｂ４の双方がしきい値Ｂ３４ｔｈ以上である場合（ステップＳ２５でＯＫ）、第２スペクトル評価部４４による、第２評価用スペクトルＳＰ２０の評価が行われる。より具体的には、第２スペクトル評価部４４は、第２評価用スペクトルＳＰ２０が信頼性を有する値を示しているか否かを判定する（ステップＳ２６）。 On the other hand, when both the reliability determination indexes B3 and B4 are equal to or greater than the threshold value B34th (OK in step S25), the second
第２方位スペクトルＳＰ２１および第２周波数スペクトルＳＰ２２の少なくとも一方において、対応する判定指標Ｂ２１，Ｂ２２がしきい値Ｂ２ｔｈ未満である場合（ステップＳ２６でＮＧ）、第２スペクトル評価部４４は、解析領域Ｚ１０における波１０２に関して、波浪情報の信頼度が低いと判定する（ステップＳ２７）。そして、第２スペクトル評価部４４は、この判定結果を示すデータＤ８を２次元スペクトル変換部１８へ出力する。 In at least one of the second azimuth spectrum SP21 and the second frequency spectrum SP22, when the corresponding determination indicators B21 and B22 are less than the threshold value B2th (NG in step S26), the second
一方、第２評価用スペクトルＳＰ２０の第２方位スペクトルＳＰ２１および第２周波数スペクトルＳＰ２２の何れにおいても、対応する判定指標Ｂ２１，Ｂ２２がしきい値Ｂ２ｔｈ以上である場合（ステップＳ２６でＯＫ）、第２スペクトル評価部４４は、解析領域Ａ１０（第２評価用エリアＡＲ２０）における波１０２に関して、波浪情報の信頼度が高いと判定する（ステップＳ２８）。この場合、第２評価用スペクトル評価部４４は、この判定結果を示すデータＤ８を２次元スペクトル変換部１８へ出力する。 On the other hand, in any of the second azimuth spectrum SP21 and the second frequency spectrum SP22 of the second evaluation spectrum SP20, when the corresponding determination indicators B21 and B22 are equal to or more than the threshold value B2th (OK in step S26), the second The
以上の処理の流れによって、表示用データ生成部１０には、波１０１についての波浪情報が信頼性を有しているか否かについてのデータと、波１０２についての波浪情報が信頼性を有しているか否かについてのデータと、が与えられる。表示用データ生成部１０は、これらのデータ（すなわち、信頼性データ）に基づいて、表示器３に画像を表示させるためのデータを生成する。 According to the above processing flow, the display
以上説明したように、本実施形態によると、信頼性評価部９は、自船周囲領域Ｚ１の一部である解析領域Ｚ１０について、波浪情報の信頼性を評価する。この構成によると、信頼性評価部９は、波浪データのうち、解析領域Ｚ１０において比較的大きな波が生じていることを示すデータに着目して、波浪情報の信頼性を評価することができる。このような構成により、信頼性評価部９は、波浪情報を信頼できるか否かについて、より確実に判定することが可能である。 As described above, according to the present embodiment, the reliability evaluation unit 9 evaluates the reliability of the wave information for the analysis area Z10 that is a part of the ship surrounding area Z1. According to this configuration, the reliability evaluation unit 9 can evaluate the reliability of the wave information by paying attention to the data indicating that a relatively large wave is generated in the analysis region Z10 among the wave data. With such a configuration, the reliability evaluation unit 9 can more reliably determine whether or not the wave information is reliable.
また、本実施形態によると、信頼性評価部９は、解析領域Ｚ１０に関して、波浪データに基づいて波浪に関するスペクトルＳＰ１０，ＳＰ２０を算出する評価用スペクトル算出部４１，４２と、スペクトルＳＰ１０，ＳＰ２０に基づいて信頼性を評価するスペクトル評価部４３，４４と、を含んでいる。この構成によると、信頼性評価部９は、波浪に関するスペクトルＳＰ１０，ＳＰ２０を用いて、波浪情報の信頼性を評価できる。 Further, according to the present embodiment, the reliability evaluation unit 9 calculates the spectrums SP10 and SP20 relating to the waves based on the wave data with respect to the analysis region Z10, and the evaluation
また、本実施形態によると、スペクトル評価部４１，４２は、対応するスペクトルＳＰ１０，ＳＰ２０の強さの基準値としての中央ピーク値Ｐ１１，Ｐ１２，Ｐ２１，Ｐ２２と対応する平均値ＰＡ１１，ＰＡ１２，ＰＡ１２，ＰＡ２２との関係から信頼性の判定指標Ｂ１１，Ｂ１２，Ｂ２１，Ｂ２２を算出する。また、この判定指標Ｂ１１，Ｂ１２，Ｂ２１，Ｂ２２は、中央ピーク値Ｐ１１，Ｐ１２，Ｐ２１，Ｐ２２のそれぞれを対応する平均値ＰＡ１１，ＰＡ１２，ＰＡ２１，ＰＡ２２で除した値である。また、スペクトル評価部４３，４４は、判定指標Ｂ１１，Ｂ１２，Ｂ２１，Ｂ２２が対応する所定のしきい値Ｂ１ｔｈ，Ｂ２ｔｈ以上であるか否かに基づいて、信頼性の有無を評価する。この構成によると、各成分スペクトルＳＰ１１，ＳＰ１２，ＳＰ２１，ＳＰ２２の形状の急峻さ（急峻度合）に基づいて、波浪情報の信頼性をより正確に評価できる。 Further, according to the present embodiment, the
また、本実施形態によると、信頼性評価部９の第１スペクトル抽出部４９は、第１評価用エリアＡＲ１０において、最も大きいピーク値を有する第１評価用スペクトルＳＰ１０を抽出する。そして、第１スペクトル評価部４３は、第１評価用スペクトルＳＰ１０に基づいて波浪情報の信頼性を判定する。この構成によると、信頼性評価部９は、解析領域Ｚ１０内の波に関するスペクトル形状を最も明確に評価し易い当該スペクトルを用いて、波浪情報の信頼性を評価できる。 Further, according to the present embodiment, the first
また、本実施形態によると、信頼性評価部９の第２スペクトル抽出部５９は、評価対象エリアＡＲ１のうち第１評価用エリアＡＲ１０を除くエリアにおいて、最も強いエコー強度が観測された箇所を中央に含む第２評価用エリアＡＲ２０であって、ピーク値Ｐ２を有する第２評価用スペクトルＳＰ２０を、抽出する。そして、第２スペクトル評価部４４は、第２評価用スペクトルＳＰ２０に基づいて波浪情報の信頼性を判定する。この構成によると、信頼性評価部９は、解析領域Ｚ１０内の波に関するスペクトル形状を最も明確に評価し易い第１評価用スペクトルＳＰ１０に加えて、自船周囲領域ＡＲ１の波に関するスペクトル形状を２番目に評価し易い第２評価用スペクトルＳＰ２０を用いて、波浪情報の信頼性を評価できる。 Further, according to the present embodiment, the second
また、本実施形態によると、第２評価用エリア設定部５８は、第１評価用エリアＡＲ１０における第１評価用スペクトルＳＰ１０の強さを低減させた状態で、第２評価用エリアＡＲ２０を設定する。この構成によると、第２評価用エリア設定部５８は、第１評価用エリアＡＲ１０におけるスペクトル強度に影響されることを抑制された状態で、より正確に第２評価用エリアＡＲ２０を設定できる。 Further, according to the present embodiment, the second evaluation
また、本実施形態によると、信頼性評価部９のピーク比算出部４５は、第１評価用スペクトルＳＰ１０の中央ピーク値Ｐ１１，Ｐ１２と第２評価用スペクトルＳＰ２０の対応する中央ピーク値Ｐ２１，Ｐ２２と、の比に基づいて、信頼性の判定指標Ｂ３，Ｂ４を算出する。そして、ピーク比評価部４６は、判定指標Ｂ３，Ｂ４が所定のしきい値Ｂ３４ｔｈ以上であるか否かに基づいて、波浪情報の信頼性を評価する。この構成によると、第２評価用スペクトルＳＰ２０の中央ピーク値Ｐ２が第１評価用スペクトルＳＰ１０の中央ピーク値Ｐ１と比べて十分に小さい場合、信頼性評価部９は、第２評価用スペクトルＳＰ２０で特定される、解析領域Ｚ１０内における２番目に大きな波について、ノイズ程度の小さなスペクトルデータであると判定することができる。 Further, according to the present embodiment, the peak
また、第１評価用スペクトルＳＰ１０および第２評価用スペクトルＳＰ２０は、それぞれ、極座標における方位方向θと周波数方向ωの２次元スペクトルである。そして、第１評価用スペクトルＳＰ１０および第２評価用スペクトルＳＰ２０は、それぞれ、方位方向θのスペクトル成分としての方位スペクトルＳＰ１１，ＳＰ２１および、周波数方向ωのスペクトル成分としての周波数スペクトルＳＰ１２，ＳＰ２２を含んでいる。この構成によると、信頼性評価部９は、スペクトルの成分毎に、波浪情報の信頼性評価を行うことができる。これにより、信頼性評価部９は、より正確に、波浪情報の信頼性を評価できる。 The first evaluation spectrum SP10 and the second evaluation spectrum SP20 are two-dimensional spectra in the azimuth direction θ and the frequency direction ω in polar coordinates, respectively. The first evaluation spectrum SP10 and the second evaluation spectrum SP20 include azimuth spectra SP11 and SP21 as spectrum components in the azimuth direction θ and frequency spectra SP12 and SP22 as spectrum components in the frequency direction ω, respectively. There is. With this configuration, the reliability evaluation unit 9 can perform reliability evaluation of wave information for each spectrum component. Thereby, the reliability evaluation unit 9 can evaluate the reliability of the wave information more accurately.
また、本実施形態によると、表示器３には、信頼性評価部９での評価結果に応じて、波浪情報が表示されたり、非表示となったりする。具体的には、信頼性評価部９において信頼性が高いと評価された波に関する波浪情報については、表示器３に表示される。一方、信頼性評価部９において信頼性が高いと評価された波に関する波浪情報については、表示器３に表示されない。これにより、ユーザに対して、信頼性が高い波浪情報のみを提供することができる。 Further, according to the present embodiment, the wave information is displayed or hidden on the
以上、本発明の実施形態について説明したけれども、本発明はこれらに限定されず、本発明の趣旨を逸脱しない限りにおいて種々の変更が可能である。 Although the embodiments of the present invention have been described above, the present invention is not limited to these, and various modifications can be made without departing from the spirit of the present invention.
（１）上記実施形態では、図１５を参照して、第１スペクトルの信頼度が低い場合、具体的には、第１方位スペクトルＳＰ１１および第１周波数スペクトルＳＰ１２の少なくとも一方において、対応する判定指標Ｂ１１，Ｂ１２がしきい値Ｂ１ｔｈ未満である場合（ステップＳ２２でＮＧ）、第２スペクトルの信頼度の判定が行われないが、これに限らない。例えば、第１スペクトルの信頼度が低い場合であっても、第２スペクトルの信頼度の判定が行われてもよい。 (1) In the above embodiment, with reference to FIG. 15, when the reliability of the first spectrum is low, specifically, in at least one of the first azimuth spectrum SP11 and the first frequency spectrum SP12, the corresponding determination index. When B11 and B12 are less than the threshold value B1th (NG in step S22), the reliability of the second spectrum is not determined, but not limited to this. For example, even when the reliability of the first spectrum is low, the reliability of the second spectrum may be determined.
（２）上記実施形態では、２つの波１０１，１０２に関する波浪情報の信頼度を判定する例を挙げて説明したが、これに限らず、３つ以上のｎ個（ｎは自然数）の波に関する波浪情報の信頼度が判定されてもよい。これにより、多数の波の波浪情報について、信頼性を評価できる。 (2) In the above-described embodiment, an example in which the reliability of the wave information regarding the two
（３）図１６は、変形例に係るレーダ装置の表示器の表示画面の一例を示す図である。上述した実施形態では、信頼性が高い波浪情報は表示器３に表示される一方、信頼性が低い波浪情報については表示器３に表示されないが、これに限らない。具体的には、例えば各波について算出された波浪情報が信頼性の高低に関わらず全て表示器に表示され、それとともに、各波に関する波浪情報に対応する位置にその波に関する信頼性の度合を示す指標が表示されてもよい。本変形例では、第１波用表示領域６１及び第２波用表示領域６２のそれぞれに、各波の信頼性の度合を示す指標として、インジケータ６３が表示される。図１６で示すインジケータ６３では、信頼性が高くなるにつれて、矩形状の指示部６４の長さが右側へ向かって延びて長くなる一方、信頼性が低くなるにつれて、指示部６４の長さが短くなる。 (3) FIG. 16 is a diagram showing an example of a display screen of a display unit of the radar device according to the modification. In the above-described embodiment, the wave information having high reliability is displayed on the
以上のように、本変形例によれば、波浪情報が信頼性の度合（本変形例の場合、インジケータ６３）とともに表示されるため、ユーザは、その信頼性の度合を見ることにより、波浪情報として表示されている数値の確からしさを把握できる。 As described above, according to the present modification, the wave information is displayed together with the degree of reliability (in the case of the present modification, the indicator 63). Therefore, the user can see the degree of reliability to determine the wave information. You can grasp the certainty of the value displayed as.
なお、本変形例では、インジケータとして、図１６に示すような棒グラフ状のインジケータを例に挙げたが、これに限らず、インジケータは、その他の表示方式であってもよい。例えば一例として、信頼性が低い場合には赤が表示され、信頼性が高くなるに連れて、橙、黄、緑、青と、徐々に色調が変化するインジケータを設けてもよい。この場合、波浪情報の文字色等をそのインジケータの色調に連動させてもよい。 Note that, in the present modified example, the bar graph-shaped indicator as shown in FIG. 16 is taken as an example of the indicator, but the present invention is not limited to this, and the indicator may be another display system. For example, an indicator may be provided in which red is displayed when the reliability is low, and the color tone gradually changes from orange to yellow, green, and blue as the reliability increases. In this case, the character color of the wave information may be linked to the color tone of the indicator.
１ レーダ装置
４ アンテナ（送波器、受波器）
９ 信頼性評価部
１３ 周波数解析部（波数データ取得部）
１４ 波浪情報取得部
２０ 信号処理装置
４１ 第１スペクトル算出部
４２ 第２スペクトル算出部
４３ 第１スペクトル評価部
４４ 第２スペクトル評価部
４５ ピーク比算出部
４６ ピーク比評価部
４８ 第１評価用エリア設定部
４９ 第１スペクトル抽出部
５８ 第２評価用エリア設定部
５９ 第２スペクトル抽出部
ＡＲ１ 評価対象エリア
ＡＲ１０，ＡＲ２０ 評価用エリア
Ｂ１１，Ｂ１２，Ｂ２１，Ｂ２２，Ｂ３，Ｂ４ 信頼性の判定指標
Ｂ１ｔｈ，Ｂ２ｔｈ，Ｂ３４ｔｈ しきい値
Ｄ１ 波数データ
Ｄ２ 波浪データ
Ｐ１１，Ｐ１２，Ｐ２１，Ｐ２２ スペクトルの強さの中央ピーク値
ＰＡ１１，ＰＡ１２，ＰＡ１２，ＰＡ２２ スペクトルの強さの平均値
ＳＰ１０ 第１評価用スペクトル
ＳＰ１１，ＳＰ２１ 方位スペクトル
ＳＰ１２，ＳＰ２２ 周波数スペクトル
ＳＰ２０ 第２評価用スペクトル
Ｚ１ 自船周囲領域（観測エリア）
Ｚ１０ 解析領域1
9
14 wave
Z10 analysis area
Claims (8)
前記波数データに所定のフィルタ処理を施すことで前記波数データに含まれる波浪データを検出する波浪情報取得部と、
前記観測エリアの一部である所定の解析領域について、前記波浪データで特定される波浪情報の信頼性を評価する、信頼性評価部と、
を備え、
前記信頼性評価部は、
前記解析領域に関して、前記波浪データに基づいて前記波浪に関するスペクトルを算出するスペクトル算出部と、
前記解析領域内に含まれる波浪から得られたエコーに基づいて算出された、該波浪の角周波数と波浪の到来方向とで特定されるωθ座標の各位置に該波浪のエコー強度が表示されたωθスペクトル、において、最も強いエコー強度が観測された箇所を含む所定の第１評価用エリアを前記評価用エリアとして設定する第１評価用エリア設定部と、
前記第１評価用エリアにおいて最も大きいピーク値を有する所定の第１評価用スペクトルを抽出する第１評価用スペクトル抽出部と、
前記ωθ座標に含まれる領域である評価対象エリアのうち前記第１評価用エリアを除くエリアにおいて、最も強いエコー強度が観測された箇所を中心に含む所定の第２評価用エリアを前記評価用エリアとして設定する第２評価用エリア設定部と、
前記第２評価用エリアの中央においてピーク値を有する所定の第２評価用スペクトルを抽出する第２評価用スペクトル抽出部と、
前記第１評価用スペクトルの前記ピーク値と前記第２評価用スペクトルの前記中央におけるピーク値と、の比に基づいて、前記信頼性の判定指標を算出するピーク比算出部と、
を備えることを特徴とする信号処理装置。 Based on the observation data obtained by observing the waves in a predetermined observation area, a wave number data acquisition unit that acquires wave number data including the relationship between the wave number and the frequency,
A wave information acquisition unit that detects wave data included in the wave number data by performing a predetermined filtering process on the wave number data,
For a predetermined analysis region that is a part of the observation area, to evaluate the reliability of the wave information specified by the wave data, a reliability evaluation unit,
Equipped with
The reliability evaluation unit,
With respect to the analysis region, a spectrum calculation unit that calculates a spectrum regarding the wave based on the wave data,
The echo intensity of the wave was displayed at each position of the ωθ coordinate specified by the angular frequency of the wave and the arrival direction of the wave, which was calculated based on the echo obtained from the wave included in the analysis region. In the ωθ spectrum, a first evaluation area setting unit that sets a predetermined first evaluation area including a portion where the strongest echo intensity is observed as the evaluation area,
A first evaluation spectrum extraction unit that extracts a predetermined first evaluation spectrum having the largest peak value in the first evaluation area;
In the evaluation target area that is the area included in the ωθ coordinate, excluding the first evaluation area, a predetermined second evaluation area including a location where the strongest echo intensity is observed is the evaluation area. A second evaluation area setting unit that is set as
A second evaluation spectrum extraction unit for extracting a predetermined second evaluation spectrum having a peak value in the center of the second evaluation area;
A peak ratio calculator that calculates the reliability determination index based on a ratio of the peak value of the first evaluation spectrum and the peak value of the second evaluation spectrum at the center;
A signal processing device comprising:
前記信頼性評価部は、ピーク比評価部を含み、
前記ピーク比評価部は、前記判定指標が所定のしきい値以上であるか否かに基づいて、
前記信頼性を評価することを特徴とする、信号処理装置。 The signal processing device according to claim 1, wherein
The reliability evaluation unit includes a peak ratio evaluation unit,
The peak ratio evaluation unit, based on whether the determination index is a predetermined threshold or more,
A signal processing device, characterized in that the reliability is evaluated.
前記第２評価用エリア設定部は、前記第１評価用エリアにおける前記第１評価用スペクトルの強さを低減させた状態で、前記第２評価用エリアを設定することを特徴とする、信号処理装置。 The signal processing device according to claim 1 or 2, wherein
The second evaluation area setting unit sets the second evaluation area in a state in which the strength of the first evaluation spectrum in the first evaluation area is reduced, and signal processing. apparatus.
前記波数データに分散関係に近接するスペクトルのデータを抽出するフィルタ処理を施すことで前記波数データに含まれる波浪データを検出する波浪情報取得部と、
前記観測エリアの一部である所定の解析領域について、前記波浪データで特定される波浪情報の信頼性を評価する、信頼性評価部と、
を備え、
前記信頼性評価部は、
前記解析領域に関して、前記波浪データに基づいて前記波浪に関するスペクトルを算出するスペクトル算出部と、
前記解析領域内に含まれる波浪から得られたエコーに基づいて算出された、該波浪の角周波数と波浪の到来方向とで特定されるωθ座標の各位置に該波浪のエコー強度が表示されたωθスペクトル、において、最も強いエコー強度が観測された箇所を含む所定の第１評価用エリアを前記評価用エリアとして設定する第１評価用エリア設定部と、
前記第１評価用エリアにおいて最も大きいピーク値を有する所定の第１評価用スペクトルを抽出する第１評価用スペクトル抽出部と、
を含み、
前記信頼性評価部は、前記第１評価用スペクトルの強さの基準値を、前記第１評価用スペクトルの強さの平均値で除した値に基づいて信頼性を評価することを特徴とする、信号処理装置。 Based on the observation data obtained by observing the waves in a predetermined observation area, a wave number data acquisition unit that acquires wave number data including the relationship between the wave number and the frequency,
A wave information acquisition unit that detects wave data included in the wave number data by performing a filtering process to extract spectrum data that is close to a dispersion relation to the wave number data,
For a predetermined analysis region that is a part of the observation area, to evaluate the reliability of the wave information specified by the wave data, a reliability evaluation unit,
Equipped with
The reliability evaluation unit,
With respect to the analysis region, a spectrum calculation unit that calculates a spectrum regarding the wave based on the wave data,
The echo intensity of the wave was displayed at each position of the ωθ coordinate specified by the angular frequency of the wave and the arrival direction of the wave, which was calculated based on the echo obtained from the wave included in the analysis region. In the ωθ spectrum, a first evaluation area setting unit that sets a predetermined first evaluation area including a portion where the strongest echo intensity is observed as the evaluation area,
A first evaluation spectrum extraction unit that extracts a predetermined first evaluation spectrum having the largest peak value in the first evaluation area;
Including,
The reliability evaluation unit evaluates reliability based on a value obtained by dividing a reference value of the strength of the first evaluation spectrum by an average value of the strength of the first evaluation spectrum. , Signal processing device.
前記スペクトルは、極座標における方位方向と周波数方向の２次元スペクトルであり、
前記方位方向のスペクトル成分としての方位スペクトル、および、前記周波数方向のスペクトル成分としての周波数スペクトルを含んでいることを特徴とする、信号処理装置。 The signal processing device according to any one of claims 1 to 4,
The spectrum is a two-dimensional spectrum in the azimuth direction and the frequency direction in polar coordinates,
A signal processing device comprising an azimuth spectrum as a spectrum component in the azimuth direction and a frequency spectrum as a spectrum component in the frequency direction.
前記送波器から送波された送信波が波浪に反射して帰来するエコーを受波する受波器と、
前記受波器で受波されたエコーを処理する請求項１から請求項５の何れか１項に記載の信号処理装置と、
を備えていることを特徴とする、レーダ装置。 A transmitter that transmits a transmission wave,
A wave receiver that receives an echo that is returned from the wave transmitted by the wave transmitter and reflected in waves,
The signal processing device according to any one of claims 1 to 5, which processes an echo received by the wave receiver,
A radar device comprising:
前記信頼性評価部での評価結果に応じて前記波浪情報の表示と非表示とが切り替えられる、又は前記信頼性評価部での評価結果を示す指標を表示する表示器、を更に備えていることを特徴とする、レーダ装置。 The radar device according to claim 6, wherein
The display further comprises a display for switching between display and non-display of the wave information according to an evaluation result of the reliability evaluation unit, or an indicator showing an evaluation result of the reliability evaluation unit. A radar device.
前記波数データに分散関係に近接するスペクトルのデータを抽出するフィルタ処理を施すことで前記波数データに含まれる波浪データを検出するステップと、
前記観測エリアの一部である所定の解析領域について、前記波浪データで特定される波浪情報の信頼性を評価するステップと、を含むことを特徴とする、信号処理方法であって、
前記信頼性を評価するステップは、
前記解析領域に関して、前記波浪データに基づいて前記波浪に関するスペクトルを算出するステップと、
前記解析領域内に含まれる波浪から得られたエコーに基づいて算出された、該波浪の角周波数と波浪の到来方向とで特定されるωθ座標の各位置に該波浪のエコー強度が表示されたωθスペクトル、において、最も強いエコー強度が観測された箇所を含む所定の第１評価用エリアを前記評価用エリアとして設定するステップと、
前記第１評価用エリアにおいて最も大きいピーク値を有する所定の第１評価用スペクトルを抽出するステップと、
を含み、
前記信頼性を評価するステップは、前記第１評価用スペクトルの強さの基準値を、前記第１評価用スペクトルの強さの平均値で除した値に基づいて信頼性を評価することを特徴とする、信号処理方法。 Based on the observation data obtained by observing the waves in a predetermined observation area, a step of acquiring wave number data including the relationship between the wave number and the frequency,
A step of detecting the wave data included in the wave number data by performing a filtering process for extracting the spectrum data close to the wave number data in a dispersion relation ;
A predetermined analysis region that is a part of the observation area, a step of evaluating the reliability of the wave information specified by the wave data, a signal processing method, comprising:
The step of evaluating the reliability includes:
Calculating a spectrum relating to the wave based on the wave data with respect to the analysis region;
The echo intensity of the wave was displayed at each position of the ωθ coordinate specified by the angular frequency of the wave and the arrival direction of the wave, which was calculated based on the echo obtained from the wave included in the analysis region. In the ωθ spectrum, a step of setting a predetermined first evaluation area including a portion where the strongest echo intensity is observed as the evaluation area,
Extracting a predetermined first evaluation spectrum having the largest peak value in the first evaluation area;
Including,
In the step of evaluating the reliability, the reliability is evaluated based on a value obtained by dividing a reference value of the strength of the first evaluation spectrum by an average value of the strength of the first evaluation spectrum. And the signal processing method.
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2016079168 | 2016-04-11 | ||
| JP2016079168 | 2016-04-11 | ||
| PCT/JP2017/009494 WO2017179342A1 (en) | 2016-04-11 | 2017-03-09 | Signal processing device, radar device, and signal processing method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2017179342A1 JPWO2017179342A1 (en) | 2019-02-14 |
| JP6689961B2 true JP6689961B2 (en) | 2020-04-28 |
Family
ID=60042091
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018511928A Active JP6689961B2 (en) | 2016-04-11 | 2017-03-09 | Signal processing device, radar device, and signal processing method |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP6689961B2 (en) |
| WO (1) | WO2017179342A1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP7219603B2 (en) * | 2018-12-14 | 2023-02-08 | 日本無線株式会社 | Sea level measurement device and sea level measurement program |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002008786A1 (en) * | 2000-07-21 | 2002-01-31 | Gkss-Forschungszentrum Geesthacht Gmbh | Method for determining hydrographic parameters, which describe a sea swell field in situ, using a radar device |
| JP3888671B2 (en) * | 2001-12-19 | 2007-03-07 | よこはまティーエルオー株式会社 | Wave height calculation device, wave height calculation method, recording medium, and ship |
| JP3992101B2 (en) * | 2003-05-16 | 2007-10-17 | よこはまティーエルオー株式会社 | Individual wave prediction / warning system |
| JP3786205B2 (en) * | 2003-06-13 | 2006-06-14 | 日本無線株式会社 | Radar wave observation device |
| JP5996259B2 (en) * | 2012-05-09 | 2016-09-21 | 古野電気株式会社 | Tsunami detection device, tsunami detection method, and tsunami detection program |
| JP6154219B2 (en) * | 2013-07-03 | 2017-06-28 | 古野電気株式会社 | Echo signal processing device, wave radar device, echo signal processing method, and echo signal processing program |
-
2017
- 2017-03-09 JP JP2018511928A patent/JP6689961B2/en active Active
- 2017-03-09 WO PCT/JP2017/009494 patent/WO2017179342A1/en active Application Filing
Also Published As
| Publication number | Publication date |
|---|---|
| JPWO2017179342A1 (en) | 2019-02-14 |
| WO2017179342A1 (en) | 2017-10-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN202083394U (en) | Intelligent navigation system for ship during navigating | |
| JP6712313B2 (en) | Signal processing device and radar device | |
| CN103969639B (en) | The signal processing system of five wave beam fish detectors and signal processing method thereof | |
| CN107110967B (en) | Tracking processing device and tracking processing method | |
| JP6714261B2 (en) | Underwater detection device, underwater detection method, and underwater detection program | |
| CN106483520B (en) | A kind of Ship Radiated-Noise index of modulation estimation method | |
| JP2016206153A (en) | Signal processor and radar device | |
| JP6689961B2 (en) | Signal processing device, radar device, and signal processing method | |
| GB2553404A (en) | Underwater detection system | |
| JP5650420B2 (en) | Signal processing apparatus, radar apparatus, signal processing method, and signal processing program | |
| EP3144697B1 (en) | Radar apparatus | |
| KR101293840B1 (en) | The ship safety navigation system by real-time sharing of wave information | |
| JP5423713B2 (en) | Active sonar device and signal processing method thereof | |
| JP2016206152A (en) | Wave radar device | |
| CN202126500U (en) | Ship navigation monitoring control system | |
| JP2011112373A (en) | Radar signal processing apparatus | |
| EP3754371A1 (en) | Underwater detection apparatus and underwater detection method | |
| JP3905880B2 (en) | Wave observation device | |
| EP3144698A1 (en) | Radar apparatus | |
| JP2011002426A (en) | Radar system | |
| JP2012141189A (en) | Representative position detection device, display device, radar device, representative position detection method and representative position detection program | |
| JP2016121970A (en) | Echo signal processing device, radar device, echo signal processing method, and program | |
| JP2013205268A (en) | Radar signal processor | |
| JP5936857B2 (en) | Ultrasonic diagnostic apparatus and image processing program thereof | |
| JP2014020906A (en) | Signal processor, underwater searching device, signal processing method, and program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20181024 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20190903 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20191105 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200407 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200408 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6689961 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |