JP7222689B2 - Sea level measurement device and sea level measurement program - Google Patents
Sea level measurement device and sea level measurement program Download PDFInfo
- Publication number
- JP7222689B2 JP7222689B2 JP2018234949A JP2018234949A JP7222689B2 JP 7222689 B2 JP7222689 B2 JP 7222689B2 JP 2018234949 A JP2018234949 A JP 2018234949A JP 2018234949 A JP2018234949 A JP 2018234949A JP 7222689 B2 JP7222689 B2 JP 7222689B2
- Authority
- JP
- Japan
- Prior art keywords
- wave
- sea
- image
- sea surface
- divided image
- 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
Images
Classifications
-
- 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
- Radar Systems Or Details Thereof (AREA)
Description
本開示は、海面の波浪を計測する技術に関する。 The present disclosure relates to technology for measuring sea surface waves.
海面の波浪を計測する技術が、特許文献1等に開示されている。特許文献1では、複数の浮標を海面に配置しカメラで撮影することにより、海面の波浪を計測する。
A technique for measuring waves on the sea surface is disclosed in
特許文献1では、複数の浮標を保守管理する必要があるが、保守管理は手間がかかり、複数の浮標の動作を画像検出する必要があるが、画像検出は困難であった。
In
そこで、前記課題を解決するために、本開示は、海面の波浪を計測するにあたり、保守管理を軽減するとともに、画像検出を容易にすることを目的とする。 Therefore, in order to solve the above problems, an object of the present disclosure is to reduce maintenance and management and facilitate image detection when measuring waves on the sea surface.
前記課題を解決するために、レーダ画像の海面クラッタを計測することにより、海面の波浪を計測することとした。具体的には、レーダ画像を分割した各分割画像での波数空間解析又は実空間解析に基づいて、各分割画像での海面の擾乱を検出することとした。 In order to solve the above problems, we decided to measure sea surface waves by measuring sea surface clutter in radar images. Specifically, based on the wavenumber space analysis or real space analysis of each divided image obtained by dividing the radar image, the sea surface disturbance is detected in each divided image.
具体的には、本開示は、レーダ画像を分割して各分割画像を生成する画像分割部と、各分割画像での波数空間解析又は実空間解析に基づいて、各分割画像での海面の波長、波向き、波速及び波高の少なくともいずれかを計測する海面計測部と、一の分割画像が表す評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、前記一の分割画像の周囲分割画像が表す周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、の比較に基づいて、前記評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出する擾乱検出部と、を備えることを特徴とする海面計測装置である。 Specifically, the present disclosure includes an image dividing unit that divides a radar image to generate each divided image, and based on wave number space analysis or real space analysis in each divided image, the wavelength of the sea surface in each divided image , a sea surface measurement unit that measures at least one of wave direction, wave speed and wave height; at least one of sea surface wavelength, wave direction, wave speed and wave height in an evaluation sea area represented by one divided image; Based on the comparison with at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in the surrounding sea area represented by the surrounding divided image, at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in the evaluation sea area and a disturbance detection unit that detects the presence of disturbance.
また、本開示は、レーダ画像を分割して各分割画像を生成する画像分割ステップと、各分割画像での波数空間解析又は実空間解析に基づいて、各分割画像での海面の波長、波向き、波速及び波高の少なくともいずれかを計測する海面計測ステップと、一の分割画像が表す評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、前記一の分割画像の周囲分割画像が表す周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、の比較に基づいて、前記評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出する擾乱検出ステップと、を順にコンピュータに実行させるための海面計測プログラムである。 Further, the present disclosure is based on an image division step of dividing a radar image to generate each divided image, and wave number space analysis or real space analysis on each divided image. , a sea surface measurement step of measuring at least one of wave speed and wave height; at least one of wavelength, wave direction, wave speed and wave height of sea surface in an evaluation sea area represented by one divided image; and surrounding divided images of said one divided image Presence of disturbances in the wavelength, direction, speed and/or height of the sea surface in the assessment sea area based on a comparison of the wavelength, direction, speed and/or height of the sea surface in the surrounding sea area represented by and a sea level measurement program for causing a computer to sequentially execute a disturbance detection step for detecting .
これらの構成によれば、海面の波浪を計測するにあたり、レーダシステムの保守管理を行うのみでよく、海面クラッタの画像検出を行うのみでよい。そして、各分割画像が表す各海域において、海面の擾乱を個別に検出することができる。 According to these configurations, when measuring waves on the sea surface, it is sufficient to perform only maintenance of the radar system and to perform image detection of sea surface clutter. Then, in each sea area represented by each divided image, sea surface disturbance can be detected individually.
また、本開示は、前記擾乱検出部は、前記評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、前記周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかの前記周囲海域にわたる平均値と、の間の差分に基づいて、前記評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出することを特徴とする海面計測装置である。 Further, in the present disclosure, the disturbance detection unit detects at least one of the sea surface wavelength, wave direction, wave speed, and wave height in the evaluation sea area, and at least any of the sea surface wavelength, wave direction, wave speed, and wave height in the surrounding sea area. The sea surface characterized by detecting the presence of disturbance of at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in the evaluation sea area based on the difference between the average value over the surrounding sea area and It is a measuring device.
この構成によれば、各分割画像が表す各海域において、上記差分が大きい(小さい)ときには、海面の擾乱が存在する(存在しない)ことを個別に検出することができる。 According to this configuration, in each sea area represented by each divided image, when the difference is large (small), it is possible to individually detect the presence (absence) of sea surface disturbance.
また、本開示は、前記画像分割部が、各分割画像の大きさを可変に設定することにより、前記擾乱検出部は、前記一の分割画像が表す前記評価海域の大きさを可変に設定することを特徴とする海面計測装置である。 Further, in the present disclosure, the image dividing unit variably sets the size of each divided image, so that the disturbance detection unit variably sets the size of the evaluation sea area represented by the one divided image. This is a sea level measuring device characterized by:
この構成によれば、小さい分割画像での局所的な海面計測結果と、大きい分割画像での大局的かつ高精度な海面計測結果と、のいずれかを選択することができる。 According to this configuration, it is possible to select either a local sea level measurement result with a small divided image or a global and highly accurate sea level measurement result with a large divided image.
また、本開示は、前記画像分割部が、各分割画像の形状を可変に設定することにより、前記擾乱検出部は、前記一の分割画像が表す前記評価海域の形状を可変に設定することを特徴とする海面計測装置である。 Further, in the present disclosure, the image dividing unit variably sets the shape of each divided image, so that the disturbance detection unit variably sets the shape of the evaluation sea area represented by the one divided image. It is a sea level measuring device characterized by
この構成によれば、海面の擾乱の異方性が海岸等で大きい(沖合等で小さい)ときには、各分割画像の形状に異方性を持たせる(持たせない)ことができる。 According to this configuration, when the anisotropy of the sea surface disturbance is large at the coast or the like (small at the offshore or the like), the shape of each divided image can be made anisotropic (or not made to have).
このように、本開示は、海面の波浪を計測するにあたり、保守管理を軽減するとともに、画像検出を容易にすることができる。 Thus, the present disclosure can reduce maintenance and facilitate image detection when measuring waves on the sea surface.
添付の図面を参照して本開示の実施形態を説明する。以下に説明する実施形態は本開示の実施の例であり、本開示は以下の実施形態に制限されるものではない。 Embodiments of the present disclosure will be described with reference to the accompanying drawings. The embodiments described below are examples of implementing the present disclosure, and the present disclosure is not limited to the following embodiments.
本開示のレーダ目標検出システムの構成を図1に示す。本開示の海面計測処理の手順を図2に示す。レーダ目標検出システムRは、レーダ送信部1、レーダ受信部2、海面計測装置3及びレーダ映像表示部4から構成される。海面計測装置3は、画像分割部31、海面計測部32及び擾乱検出部33から構成され、図2に示した海面計測プログラムをコンピュータにインストールすることにより実現することができる。
FIG. 1 shows the configuration of the radar target detection system of the present disclosure. FIG. 2 shows the procedure of sea level measurement processing according to the present disclosure. The radar target detection system R is composed of a
レーダ送信部1は、海上に向けてレーダビームを照射する。レーダ受信部2は、海上で反射されたレーダビームを受信する。海面計測装置3は、レーダ画像の海面クラッタを計測することにより、海面の波浪を計測する。レーダ映像表示部4は、海面計測装置3が取得、処理及び作成したレーダ映像のデータを映像化して表示する。
A
本開示の第1の海面計測処理の手順を図3に示す。図3の左上欄では、画像分割部31は、レーダ映像を取得する。図3の右上欄では、画像分割部31は、レーダ画像を分割して各分割画像を生成する(ステップS1)。例えば、レーダ画像は、x方向に4分割され、y方向に4分割され、全体として4×4=16分割されている。
FIG. 3 shows the procedure of the first sea level measurement process of the present disclosure. In the upper left column of FIG. 3, the
図3の左右下欄では、海面計測部32は、各分割画像での波数空間解析又は実空間解析に基づいて、各分割画像での海面の波長、波向き、波速及び波高の少なくともいずれかを計測する(ステップS2)。例えば、海面の波長は、図3の左右下欄の矢印の長さにより表され、海面の波向きは、図3の左右下欄の矢印の方向で表され、海面の波速は、図3の左右下欄の数値で表され、海面の波高は、図3の左右下欄の階調で表される。
In the left and right lower columns of FIG. 3, the sea
図3の左右下欄では、擾乱検出部33は、一の分割画像が表す評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、一の分割画像の周囲分割画像が表す周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、の比較に基づいて、評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出する(ステップS3)。ここで、16個の各分割画像のうち、上から3行目かつ左から1列目の分割画像において、定置網が設置されており、海面の擾乱が存在しうる。
In the left and right lower columns of FIG. 3, the
図3の左下欄では、16個の各分割画像のうち、上から2行目かつ左から3列目の分割画像について、評価海域とされている。そして、16個の各分割画像のうち、上記の分割画像の周囲の分割画像について、周囲海域とされている。さらに、上記の評価海域での海面の波向き(図3の左下欄の矢印の方向)は、上記の周囲海域での海面の波向きと比べて、大差はない。よって、上記の評価海域での海面の波向きの擾乱は存在しない。 In the lower left column of FIG. 3, among the 16 divided images, the divided image on the second row from the top and the third column from the left is the evaluation sea area. Among the 16 divided images, the divided images surrounding the above divided images are regarded as the surrounding sea area. Furthermore, the direction of waves on the sea surface in the evaluation sea area (the direction of the arrow in the lower left column of FIG. 3) does not differ greatly from the direction of waves on the sea surface in the surrounding sea area. Therefore, there is no wave direction disturbance on the sea surface in the above evaluation sea area.
図3の右下欄では、16個の各分割画像のうち、上から3行目かつ左から1列目の分割画像について、評価海域とされている。そして、16個の各分割画像のうち、上記の分割画像の周囲の分割画像について、周囲海域とされている。さらに、上記の評価海域での海面の波向き(図3の右下欄の矢印の方向)は、上記の周囲海域での海面の波向きと比べて、大差がある。よって、上記の評価海域での海面の波向きの擾乱が存在する。 In the lower right column of FIG. 3, among the 16 divided images, the divided image on the third row from the top and the first column from the left is regarded as the evaluation sea area. Among the 16 divided images, the divided images surrounding the above divided images are regarded as the surrounding sea area. Furthermore, the direction of waves on the sea surface in the evaluation sea area (the direction of the arrow in the lower right column of FIG. 3) differs greatly from the direction of waves on the sea surface in the surrounding sea area. Therefore, there is a disturbance in the wave direction of the sea surface in the above evaluation sea area.
このように、海面の波浪を計測するにあたり、レーダ目標検出システムRの保守管理を行うのみでよく、海面クラッタの画像検出を行うのみでよい。そして、各分割画像が表す各海域において、海面の擾乱を個別に検出することができる。 In this way, when measuring waves on the sea surface, it is only necessary to perform maintenance and management of the radar target detection system R, and it is only necessary to perform image detection of sea surface clutter. Then, in each sea area represented by each divided image, sea surface disturbance can be detected individually.
本開示の第2の海面計測処理の手順を図4に示す。擾乱検出部33は、評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかの周囲海域にわたる平均値と、の間の差分に基づいて、評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出する(ステップS3)。
FIG. 4 shows the procedure of the second sea level measurement process of the present disclosure. The
図4の左上欄では、16個の各分割画像のうち、上から2行目かつ左から3列目の分割画像について、評価海域とされている。そして、16個の各分割画像のうち、上記の分割画像の周囲の分割画像について、周囲海域とされている。さらに、上記の評価海域について、海面の波浪のパラメータがeとされ、上記の周囲海域について、海面の波浪のパラメータがp1、p2、p3、p4、p5、p6、p7、p8とされている。 In the upper left column of FIG. 4, among the 16 divided images, the divided image on the second row from the top and the third column from the left is the evaluation sea area. Among the 16 divided images, the divided images surrounding the above divided images are regarded as the surrounding sea area. Further, for the evaluation sea area, the sea surface wave parameter is set to e, and for the surrounding sea area, the sea surface wave parameters are p 1 , p 2 , p 3 , p 4 , p 5 , p 6 , p 7 , p8 .
図4の右上欄及び左下欄では、上記の周囲海域について、海面の波浪のパラメータの平均値がp=Σpi/8とされ、上記の評価海域について、海面の波浪のパラメータと上記の平均値との間の差分がe’=e-pとされ、上記の周囲海域について、海面の波浪のパラメータと上記の平均値との間の差分がpi’=pi-p(i=1~8)とされている。 In the upper right column and the lower left column of FIG. 4, the average value of the sea surface wave parameter is p=Σp i /8 for the above surrounding sea area, and the sea surface wave parameter and the above average value for the above evaluation sea area is given as e′=ep, and for the above surrounding sea area, the difference between the sea surface wave parameter and the above average value is p i '=p i -p (i=1 to 8).
図4の右下欄では、上記の評価海域について、上記の差分e’がある程度は0と一致するときには、上記の周囲海域と比べて、海面の波浪のパラメータの擾乱が存在しない。一方で、上記の評価海域について、上記の差分e’がある程度も0と一致しないときには、上記の周囲海域と比べて、海面の波浪のパラメータの擾乱が存在する。 In the lower right column of FIG. 4, there is no perturbation of the sea surface wave parameters for the above evaluation sea area compared to the above surrounding sea area when the above difference e' is to some extent equal to 0. On the other hand, for the above evaluation sea area, when the above difference e' is not equal to 0 to some extent, there is a perturbation of the sea surface wave parameters compared to the above surrounding sea area.
このように、各分割画像が表す各海域において、上記の差分e’が大きい(小さい)ときには、海面の擾乱が存在する(存在しない)ことを個別に検出することができる。 In this way, in each sea area represented by each divided image, when the difference e' is large (small), it is possible to individually detect the presence (absence) of sea surface disturbance.
本開示の第3の海面計測処理の手順を図5に示す。画像分割部31が、各分割画像の大きさを可変に設定することにより(ステップS1)、擾乱検出部33は、一の分割画像が表す評価海域の大きさを可変に設定する(ステップS3)。
FIG. 5 shows the procedure of the third sea level measurement process of the present disclosure. The
図5の左欄では、レーダ画像は、x方向に4分割され、y方向に4分割され、全体として4×4=16分割されている。そして、16個の各分割画像について、評価海域とされている。さらに、16個の各評価海域について、海面の波浪のパラメータの擾乱の有無が検出されている。よって、大局的かつ高精度な海面計測を実行することができる。 In the left column of FIG. 5, the radar image is divided into 4 parts in the x direction and 4 parts in the y direction, for a total of 4×4=16 parts. Each of the 16 divided images is set as an evaluation sea area. Further, for each of the 16 evaluation sea areas, the presence or absence of perturbation of sea surface wave parameters is detected. Therefore, global and highly accurate sea level measurement can be performed.
図5の右欄では、レーダ画像は、x方向に8分割され、y方向に8分割され、全体として8×8=64分割されている。そして、64個の各分割画像について、評価海域とされている。さらに、64個の各評価海域について、海面の波浪のパラメータの擾乱の有無が検出されている。よって、局所的な海面計測を実行することができる。 In the right column of FIG. 5, the radar image is divided into 8 divisions in the x direction and 8 divisions in the y direction, for a total of 8×8=64 divisions. Then, each of the 64 divided images is set as an evaluation sea area. Further, for each of the 64 evaluation sea areas, the presence or absence of perturbation of sea surface wave parameters is detected. Thus, local sea level measurements can be performed.
このように、小さい分割画像での局所的な海面計測結果と、大きい分割画像での大局的かつ高精度な海面計測結果と、のいずれかを選択することができる。なお、レーダ画像のうち、ある分割画像を小さくしてもよく、他の分割画像を大きくしてもよい。 In this way, it is possible to select either a local sea level measurement result with a small divided image or a global and highly accurate sea level measurement result with a large divided image. Note that, among the radar images, a certain divided image may be made smaller, and other divided images may be made larger.
本開示の第4の海面計測処理の手順を図6に示す。画像分割部31が、各分割画像の形状を可変に設定することにより(ステップS1)、擾乱検出部33は、一の分割画像が表す評価海域の形状を可変に設定する(ステップS3)。
FIG. 6 shows the procedure of the fourth sea level measurement process of the present disclosure. The
図6の左欄では、レーダ画像は、x方向に4分割され、y方向に4分割され、全体として4×4=16分割されている。これは、海面の擾乱のx方向(例えば、東西方向。)とy方向(例えば、南北方向。)との間の異方性が、沖合等で小さいからである。そして、16個の各分割画像について、評価海域とされている。さらに、16個の各評価海域について、海面の波浪のパラメータの擾乱の有無が検出されている。 In the left column of FIG. 6, the radar image is divided into 4 parts in the x direction, 4 parts in the y direction, and 4×4=16 parts in total. This is because the anisotropy between the x-direction (for example, east-west direction) and the y-direction (for example, north-south direction) of the sea surface disturbance is small offshore. Each of the 16 divided images is set as an evaluation sea area. Further, for each of the 16 evaluation sea areas, the presence or absence of perturbation of sea surface wave parameters is detected.
図6の右欄では、レーダ画像は、x方向に1分割され、y方向に4分割され、全体として1×4=4分割されている。これは、海面の擾乱のx方向(例えば、海岸線と平行な方向。)とy方向(例えば、海岸から沖合への方向。)との間の異方性が、海岸等で大きいからである。そして、4個の各分割画像について、評価海域とされている。さらに、4個の各評価海域について、海面の波浪のパラメータの擾乱の有無が検出されている。 In the right column of FIG. 6, the radar image is divided into 1 part in the x direction, 4 parts in the y direction, and 1×4=4 parts as a whole. This is because the anisotropy between the x-direction (for example, the direction parallel to the coastline) and the y-direction (for example, the direction from the coast to the offshore) of the sea surface disturbance is large at the coast. Each of the four divided images is set as an evaluation sea area. Further, for each of the four evaluation sea areas, the presence or absence of perturbation of sea surface wave parameters is detected.
このように、海面の擾乱の異方性が海岸等で大きい(沖合等で小さい)ときには、各分割画像の形状に異方性を持たせる(持たせない)ことができる。なお、レーダ画像のうち、ある分割画像を長方形としてもよく、他の分割画像を正方形としてもよい。 In this way, when the anisotropy of the sea surface disturbance is large at the coast or the like (small at the offshore or the like), the shape of each divided image can be made anisotropic (or not). Of the radar images, a certain divided image may be rectangular, and the other divided images may be square.
まず、海面計測部32が、各分割画像での「波数空間解析」に基づいて、各分割画像での海面の波長、波向き、波速及び波高を計測する方法を説明する。
First, a method for measuring the wavelength, wave direction, wave speed, and wave height of the sea surface in each divided image by the sea
本開示の海面の波長及び波向きの計測方法を図7に示す。海面計測部32は、各フーリエスペクトルの波数空間でのピーク位置に基づいて、各分割画像での海面の波長及び波向きの少なくともいずれかを計測する(ステップS2)。
FIG. 7 shows the method of measuring sea surface wavelength and wave direction of the present disclosure. The sea
具体的には、Bスコープ画像をXY座標画像に座標変換し、STC処理を実行し、円形画像から方形画像のみ切り出し、海面計測対象のレーダ画像とする。次に、前回スキャン及び今回スキャンの各分割画像について、それぞれ、前回スキャン及び今回スキャンの各スペクトルを生成する。次に、前回スキャン及び今回スキャンの各スペクトルを乗算し、前回スキャンと今回スキャンとの間の各クロススペクトルCSを生成する。次に、前回スキャンと今回スキャンとの間の各クロススペクトルCSを、kx及びkyを座標軸とする波数座標からλx=1/kx及びλy=1/kyを座標軸とする波長座標へと座標変換する。 Specifically, the B-scope image is coordinate-transformed into an XY coordinate image, STC processing is performed, and only a square image is cut out from the circular image to be used as a radar image for sea level measurement. Next, each spectrum of the previous scan and the current scan is generated for each divided image of the previous scan and the current scan. Next, each spectrum of the previous scan and the current scan is multiplied to generate each cross spectrum CS between the previous scan and the current scan. Next, each cross spectrum CS between the previous scan and the current scan is converted from wave number coordinates with k x and k y as coordinate axes to wavelengths with λ x = 1/k x and λ y = 1/k y as coordinate axes. Coordinate transformation to coordinates.
ここで、前回スキャンと今回スキャンとの間の各クロススペクトルにおいて、ピーク位置は(λxp、λyp)である。そこで、波長λ及び波向きθを数式1、2により算出する。
このように、波数空間解析により実空間解析と比べて、海面の波長及び波向きの少なくともいずれかを容易に計測することができる。そして、各分割画像が表す各海域において、海面の波長及び波向きの少なくともいずれかを個別に計測することができる。なお、海面の波長及び波向きの少なくともいずれかを計測するにあたり、乗算後の各クロススペクトルを利用してもよいが、乗算前の各スペクトルを利用してもよい。また、乗算後の各クロススペクトル又は乗算前の各スペクトルを算出するにあたり、過去から現在に渡る複数個の乗算後の各クロススペクトル又は乗算前の各スペクトルを平均してもよい。 Thus, wave number space analysis makes it possible to easily measure at least one of the wavelength and wave direction of the sea surface compared to real space analysis. Then, at least one of the wavelength and wave direction of the sea surface can be individually measured in each sea area represented by each divided image. In measuring at least one of the wavelength and wave direction of the sea surface, each cross spectrum after multiplication may be used, but each spectrum before multiplication may be used. Further, in calculating each cross spectrum after multiplication or each spectrum before multiplication, a plurality of cross spectra after multiplication or each spectrum before multiplication from the past to the present may be averaged.
本開示の海面の波速の計測方法を図8に示す。海面計測部32は、先の時刻での各フーリエスペクトルのピーク位相と、後の時刻での各フーリエスペクトルのピーク位相と、の間の差分に基づいて、各分割画像での海面の波速を計測する(ステップS2)。
FIG. 8 shows the method of measuring the wave velocity on the sea surface of the present disclosure. The sea
具体的には、Bスコープ画像をXY座標画像に座標変換し、STC処理を実行し、円形画像から方形画像のみ切り出し、海面計測対象のレーダ画像とする。次に、前回スキャン及び今回スキャンの各分割画像について、それぞれ、前回スキャン及び今回スキャンの各スペクトルを生成する。次に、前回スキャン及び今回スキャンの各スペクトルを乗算し、前回スキャンと今回スキャンとの間の各クロススペクトルCSを生成する。 Specifically, the B-scope image is coordinate-transformed into an XY coordinate image, STC processing is performed, and only a square image is cut out from the circular image to be used as a radar image for sea level measurement. Next, each spectrum of the previous scan and the current scan is generated for each divided image of the previous scan and the current scan. Next, each spectrum of the previous scan and the current scan is multiplied to generate each cross spectrum CS between the previous scan and the current scan.
ここで、前回スキャンと今回スキャンとの間の各クロススペクトルにおいて、ピーク位置は(kxp、kyp)である。そこで、各クロススペクトルのピーク位相CSθ(kxp、kyp)を数式3により算出し、波速vを数式4により算出する。ただし、τは、前回スキャンと今回スキャンとの間の時間間隔(例えば、アンテナの回転周期。)である。
このように、波数空間解析により実空間解析と比べて、海面の波速を容易に計測することができる。そして、各分割画像が表す各海域において、海面の波速を個別に計測することができる。なお、ナイキスト条件2πv/λ<π/τを満たさないときには、海面の波速を正しく計測することができないことに留意すべきである。また、乗算後の各クロススペクトル又は乗算前の各スペクトルを算出するにあたり、過去から現在に渡る複数個の乗算後の各クロススペクトル又は乗算前の各スペクトルを平均してもよい。 Thus, the wave velocity on the sea surface can be easily measured by the wave number space analysis compared to the real space analysis. Then, the wave velocity on the sea surface can be individually measured in each sea area represented by each divided image. It should be noted that when the Nyquist condition 2πv/λ<π/τ is not satisfied, the wave velocity on the sea surface cannot be measured correctly. Further, in calculating each cross spectrum after multiplication or each spectrum before multiplication, a plurality of cross spectra after multiplication or each spectrum before multiplication from the past to the present may be averaged.
本開示の海面の波高の計測方法を図9に示す。海面計測部32は、(1)先の時刻での各フーリエスペクトルの各波数での位相と、後の時刻での各フーリエスペクトルの各波数での位相と、の間の差分に基づいて、各フーリエスペクトルの各波数に対応する計測波速を算出し、(2)重力波の分散関係及び各フーリエスペクトルの各波数に対応する波長に基づいて、各フーリエスペクトルの各波数に対応する理論波速を算出し、(3)各フーリエスペクトルの各波数に対応する計測波速と理論波速との間の一致程度に基づいて、各フーリエスペクトルの各波数での振幅を信号成分と雑音成分とに分類し、(4)有義波高の経験式及び各フーリエスペクトルの全波数にわたる信号対雑音比に基づいて、各分割画像での海面の波高を計測する(ステップS2)。 FIG. 9 shows the method of measuring the wave height of the sea surface of the present disclosure. (1) Based on the difference between the phase at each wavenumber of each Fourier spectrum at the earlier time and the phase at each wavenumber of each Fourier spectrum at the later time, each Calculate the measured wave velocity corresponding to each wavenumber of the Fourier spectrum, and (2) Calculate the theoretical wave velocity corresponding to each wavenumber of each Fourier spectrum based on the dispersion relationship of gravitational waves and the wavelength corresponding to each wavenumber of each Fourier spectrum. and (3) classifying the amplitude at each wavenumber of each Fourier spectrum into a signal component and a noise component based on the degree of agreement between the measured wave velocity and the theoretical wave velocity corresponding to each wavenumber of each Fourier spectrum, ( 4) Based on the empirical formula of significant wave height and the signal-to-noise ratio over all wavenumbers of each Fourier spectrum, measure the wave height of the sea surface in each split image (step S2).
具体的には、Bスコープ画像をXY座標画像に座標変換し、STC処理を実行し、円形画像から方形画像のみ切り出し、海面計測対象のレーダ画像とする。次に、前回スキャン及び今回スキャンの各分割画像について、それぞれ、前回スキャン及び今回スキャンの各スペクトルを生成する。次に、前回スキャン及び今回スキャンの各スペクトルを乗算し、前回スキャンと今回スキャンとの間の各クロススペクトルCSを生成する。 Specifically, the B-scope image is coordinate-transformed into an XY coordinate image, STC processing is performed, and only a square image is cut out from the circular image to be used as a radar image for sea level measurement. Next, each spectrum of the previous scan and the current scan is generated for each divided image of the previous scan and the current scan. Next, each spectrum of the previous scan and the current scan is multiplied to generate each cross spectrum CS between the previous scan and the current scan.
そして、各クロススペクトルの各波数での位相CSθ(kx、ky)を数式5により算出し、各クロススペクトルの各波数での振幅CSp(kx、ky)を数式6により算出し、各クロススペクトルの各波数に対応する計測波速vm(kx、ky)を数式7により算出し、各クロススペクトルの各波数に対応する理論波速vt(kx、ky)を数式8により算出する。ここで、gは、重力加速度(=9.8m/s2)であり、重力波の分散関係において、沖合では水深は波高より十分に深いことを考慮している。
そして、計測波速vm(kx、ky)と理論波速vt(kx、ky)とがある程度は一致するときには、振幅CSp(kx、ky)を信号成分S(kx、ky)に分類する。一方で、計測波速vm(kx、ky)と理論波速vt(kx、ky)とがある程度も一致しないときには、振幅CSp(kx、ky)を雑音成分N(kx、ky)に分類する。さらに、各クロススペクトルの全波数にわたる信号対雑音比SNRを数式9により算出し、波高Hを数式10により算出する。ここで、a及びbは、有義波高の実測結果に基づく値である。
このように、波数空間解析により実空間解析と比べて、海面の波高を容易に計測することができる。そして、各分割画像が表す各海域において、海面の波高を個別に計測することができる。なお、海面の波高を計測するにあたり、乗算後の各クロススペクトルを利用してもよいが、乗算前の各スペクトルを利用してもよい。また、乗算後の各クロススペクトル又は乗算前の各スペクトルを算出するにあたり、過去から現在に渡る複数個の乗算後の各クロススペクトル又は乗算前の各スペクトルを平均してもよい。 Thus, the wave height of the sea surface can be easily measured by the wave number space analysis compared to the real space analysis. Then, the wave height of the sea surface can be individually measured in each sea area represented by each divided image. In measuring the wave height of the sea surface, each cross spectrum after multiplication may be used, but each spectrum before multiplication may be used. Further, in calculating each cross spectrum after multiplication or each spectrum before multiplication, a plurality of cross spectra after multiplication or each spectrum before multiplication from the past to the present may be averaged.
次に、海面計測部32が、各分割画像での「実空間解析」に基づいて、各分割画像での海面の波長、波向き、波速及び波高を計測する方法を説明する。
Next, a method for measuring the wavelength, wave direction, wave speed, and wave height of the sea surface in each split image by the sea
海面計測部32は、(1)レーダ反射強度が高い海面クラッタを、海面の波頭として観測し、(2)海面の波頭の間隔に基づいて、海面の波長を計測し、(3)海面の波頭の位置の変化に基づいて、海面の波向き及び波速を計測し、(4)海面クラッタのレーダ反射強度のコントラストに基づいて、海面の波高を計測する(ステップS2)。 The sea surface measurement unit 32 (1) observes sea surface clutter with high radar reflection intensity as wave crests on the sea surface, (2) measures the wavelength of the sea surface based on the interval between the wave crests on the sea surface, and (3) wave crests on the sea surface. (4) measure the wave height on the sea surface based on the contrast of the radar reflection intensity of the sea surface clutter (step S2).
本開示の海面計測装置及び海面計測プログラムは、例えば、海面の波浪が激しい海域を特定することができ、船舶又は航空機が近づけない海域を特定することができる。 The sea level measurement device and the sea level measurement program of the present disclosure can, for example, identify sea areas where sea surface waves are severe, and can identify sea areas that ships or aircraft cannot approach.
R:レーダ目標検出システム
1:レーダ送信部
2:レーダ受信部
3:海面計測装置
4:レーダ映像表示部
31:画像分割部
32:海面計測部
33:擾乱検出部
R: Radar target detection system 1: Radar transmission unit 2: Radar reception unit 3: Sea level measurement device 4: Radar image display unit 31: Image division unit 32: Sea level measurement unit 33: Disturbance detection unit
Claims (5)
各分割画像での波数空間解析又は実空間解析に基づいて、各分割画像での海面の波長、波向き、波速及び波高の少なくともいずれかを計測する海面計測部と、
一の分割画像が表す評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、前記一の分割画像の周囲分割画像が表す周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、の比較に基づいて、前記評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出する擾乱検出部と、
を備え、前記画像分割部が、各分割画像の形状の異方性を可変に設定することにより、前記擾乱検出部は、前記一の分割画像が表す前記評価海域の海面の擾乱の異方性に応じて、前記一の分割画像が表す前記評価海域の形状の異方性を可変に設定する
ことを特徴とする海面計測装置。 an image dividing unit that divides a radar image to generate each divided image;
a sea surface measurement unit that measures at least one of the wavelength, wave direction, wave speed, and wave height of the sea surface in each divided image based on wave number space analysis or real space analysis in each divided image;
At least one of the wavelength, wave direction, wave speed, and wave height of the sea surface in the evaluation sea area represented by one divided image, and the wavelength, wave direction, wave speed, and wave height of the sea surface in the surrounding sea area represented by the surrounding divided images of the one divided image a disturbance detection unit that detects the presence of disturbance of at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in the evaluation sea area based on a comparison of at least one of
wherein the image dividing unit variably sets the anisotropy of the shape of each divided image, so that the disturbance detection unit detects the anisotropy of the sea surface disturbance in the evaluation sea area represented by the one divided image variably set the anisotropy of the shape of the evaluation sea area represented by the one divided image according to
A sea level measuring device characterized by:
ことを特徴とする、請求項1に記載の海面計測装置。 The disturbance detection unit detects at least one of the sea surface wavelength, wave direction, wave speed, and wave height in the evaluation sea area, and at least one of the sea surface wavelength, wave direction, wave speed, and wave height in the surrounding sea area. Detecting the presence of disturbance of at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in the evaluation sea area based on the difference between the average value and Sea level measurement device.
ことを特徴とする、請求項1又は2に記載の海面計測装置。 The image dividing unit variably sets the size of each divided image, so that the disturbance detection unit variably sets the size of the evaluation sea area represented by the one divided image. The sea level measuring device according to claim 1 or 2.
ことを特徴とする、請求項1から3のいずれかに記載の海面計測装置。The sea level measuring device according to any one of claims 1 to 3, characterized in that:
各分割画像での波数空間解析又は実空間解析に基づいて、各分割画像での海面の波長、波向き、波速及び波高の少なくともいずれかを計測する海面計測ステップと、
一の分割画像が表す評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、前記一の分割画像の周囲分割画像が表す周囲海域での海面の波長、波向き、波速及び波高の少なくともいずれかと、の比較に基づいて、前記評価海域での海面の波長、波向き、波速及び波高の少なくともいずれかの擾乱の存在を検出する擾乱検出ステップと、
を順にコンピュータに実行させ、前記画像分割ステップが、各分割画像の形状の異方性を可変に設定することにより、前記擾乱検出ステップは、前記一の分割画像が表す前記評価海域の海面の擾乱の異方性に応じて、前記一の分割画像が表す前記評価海域の形状の異方性を可変に設定する
ことを特徴とする海面計測プログラム。 an image segmentation step of segmenting a radar image to generate each segmented image;
a sea surface measurement step of measuring at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in each divided image based on the wave number space analysis or real space analysis in each divided image;
At least one of the wavelength, wave direction, wave speed, and wave height of the sea surface in the evaluation sea area represented by one divided image, and the wavelength, wave direction, wave speed, and wave height of the sea surface in the surrounding sea area represented by the surrounding divided images of the one divided image a disturbance detection step of detecting the presence of a disturbance of at least one of the wavelength, wave direction, wave speed and wave height of the sea surface in the evaluation sea area based on a comparison of at least one of
are sequentially executed by the computer, and the image division step variably sets the anisotropy of the shape of each divided image, so that the disturbance detection step detects the disturbance of the sea surface in the evaluation sea area represented by the one divided image variably set the anisotropy of the shape of the evaluation sea area represented by the one divided image according to the anisotropy of
A sea level measurement program characterized by :
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018234949A JP7222689B2 (en) | 2018-12-14 | 2018-12-14 | Sea level measurement device and sea level measurement program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018234949A JP7222689B2 (en) | 2018-12-14 | 2018-12-14 | Sea level measurement device and sea level measurement program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2020094997A JP2020094997A (en) | 2020-06-18 |
JP7222689B2 true JP7222689B2 (en) | 2023-02-15 |
Family
ID=71084225
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2018234949A Active JP7222689B2 (en) | 2018-12-14 | 2018-12-14 | Sea level measurement device and sea level measurement program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP7222689B2 (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005156192A (en) | 2003-11-20 | 2005-06-16 | Japan Radio Co Ltd | Wave observing apparatus |
CN103364784A (en) | 2013-07-22 | 2013-10-23 | 武汉大学 | Looking-around synthetic aperture imaging radar |
JP2016085206A (en) | 2014-10-22 | 2016-05-19 | 三菱電機株式会社 | Tsunami monitoring system |
-
2018
- 2018-12-14 JP JP2018234949A patent/JP7222689B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005156192A (en) | 2003-11-20 | 2005-06-16 | Japan Radio Co Ltd | Wave observing apparatus |
CN103364784A (en) | 2013-07-22 | 2013-10-23 | 武汉大学 | Looking-around synthetic aperture imaging radar |
JP2016085206A (en) | 2014-10-22 | 2016-05-19 | 三菱電機株式会社 | Tsunami monitoring system |
Also Published As
Publication number | Publication date |
---|---|
JP2020094997A (en) | 2020-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106569196B (en) | A kind of compressed sensing based ground radar multi-target detection method | |
CN102508219A (en) | Turbulent current target detection method of wind profiler radar | |
JP2008292449A (en) | Automatic target identifying system for detecting and classifying object in water | |
CN109116352A (en) | A kind of circle sweeps ISAR mode ship super-resolution imaging method | |
JP2007248293A (en) | Ocean radar device | |
US20150051883A1 (en) | Process to identify and classify oil seep areas at the seabed through inverse modeling | |
US7714766B2 (en) | Method of processing a radar image | |
Sediono | Method of measuring Doppler shift of moving targets using FMCW maritime radar | |
CN109581384B (en) | Clear sky vertical wind profile detection method and system based on Doppler weather radar | |
JPWO2017179343A1 (en) | Signal processing apparatus and radar apparatus | |
JP7222689B2 (en) | Sea level measurement device and sea level measurement program | |
JP7219603B2 (en) | Sea level measurement device and sea level measurement program | |
CN113419238A (en) | Mountain landslide monitoring method based on millimeter wave radar, electronic device and storage medium | |
CN113514833A (en) | Sea surface arbitrary point wave direction inversion method based on sea wave image | |
CN102721411B (en) | Wave scale monitoring method based on water wave image | |
CN115546526B (en) | Three-dimensional point cloud clustering method, device and storage medium | |
KR101834063B1 (en) | Apparatus of cross-range scaling for inverse synthetic aperture radar image using principal component analysis and method thereof | |
JP2012255657A (en) | Radar image processing device and radar image processing method | |
Wengrove et al. | Monitoring morphology of the Sand Engine leeside using Argus’ cBathy | |
Wang et al. | A novel algorithm in estimating signal-to-noise ratio for ocean wave height inversion from X-band radar images | |
CN115201816A (en) | Three-dimensional ground penetrating radar quantitative imaging method based on improved backward projection | |
US20150109163A1 (en) | Target Detection System and Method | |
JP7164424B2 (en) | Water depth measurement device and water depth measurement program | |
Voulgaris et al. | 2-D inner-shelf current observations from a single VHF WEllen RAdar (WERA) station | |
KR101990078B1 (en) | Simulation Apparatus for Radar Signal Processing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20211209 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20221028 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20221122 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230105 |
|
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: 20230131 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20230203 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7222689 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |