WO2017179344A1 - Dispositif de calcul de hauteur d'onde, dispositif radar et procédé de calcul de hauteur d'onde - Google Patents

Dispositif de calcul de hauteur d'onde, dispositif radar et procédé de calcul de hauteur d'onde Download PDF

Info

Publication number
WO2017179344A1
WO2017179344A1 PCT/JP2017/009496 JP2017009496W WO2017179344A1 WO 2017179344 A1 WO2017179344 A1 WO 2017179344A1 JP 2017009496 W JP2017009496 W JP 2017009496W WO 2017179344 A1 WO2017179344 A1 WO 2017179344A1
Authority
WO
WIPO (PCT)
Prior art keywords
wave
correction coefficient
spectrum power
area
wave height
Prior art date
Application number
PCT/JP2017/009496
Other languages
English (en)
Japanese (ja)
Inventor
亮祐 森垣
敏志 川浪
健介 井芹
ミン トラン
祐也 燒山
陵 中島
Original Assignee
古野電気株式会社
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 古野電気株式会社 filed Critical 古野電気株式会社
Priority to JP2018511930A priority Critical patent/JP6676151B2/ja
Publication of WO2017179344A1 publication Critical patent/WO2017179344A1/fr

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • 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

Definitions

  • the present invention relates to a wave height calculation device that calculates the height of a wave generated on a water surface, a radar device including the wave height calculation device, and a wave height calculation method.
  • m 0 is the zero-order moment of the wave spectrum
  • alpha is a fixed coefficient.
  • the wave height H 1/3 described above is a so-called significant wave height.
  • the significant wave height is the average of the significant wave heights.
  • a significant wave means that when a continuous wave is observed at a certain point, the height of the wave height is the upper third with the number of waves observed as a parameter. (For example, if 100 waves are observed in 20 minutes, the larger 33 waves).
  • the significant wave height defined as described above substantially coincides with the wave height monitored visually.
  • Kanazawa Local Meteorological Observatory “Kameme Mame Knowledge”, [online], [Search February 24, 2016], Internet ⁇ http://www.jma-net.go.jp/kanazawa/mame/nami/nami.html >
  • the way the echoes of the wavy line are reflected has an orientation dependency.
  • a transmission wave is transmitted in the traveling direction of the wave, that is, the traveling direction of the wave based on the direction from the analysis area including the wave whose wave height is to be calculated to the receiver.
  • the relative wave direction defined by is near 0 degrees
  • the wave spectrum power ⁇ m 0 tends to be large, and the echoes of the wavy line are easily reflected clearly.
  • the transmission wave is transmitted along the traveling direction of the wave, that is, when the relative wave direction is around 180 degrees
  • the wave spectrum power ⁇ m 0 is likely to increase, and the echo of the wave line is It is easy to see clearly.
  • the transmission wave is transmitted in a direction perpendicular to the traveling direction of the wave, that is, when the relative wave direction is around 90 degrees or around 270 degrees, the wave spectrum power ⁇ m 0 becomes small. It is easy and the echo of the wavy line is not easily reflected. As described above, since the reflection of the wavy line has direction dependency, the value of ⁇ m 0 varies depending on the direction of the analysis area even when the wave height is the same. That is, in the conventional method, a wave height measurement error due to the bearing occurs.
  • the present invention is for solving the above-described problems, and its purpose is to accurately calculate the wave height regardless of the relative wave direction of the wave to be calculated.
  • a wave height calculation apparatus is obtained from a receiver that receives an echo that is transmitted and reflected by a wave generated on a water surface in a detection area.
  • a wave height calculation device for calculating a wave height of a wave based on an echo signal, a wave spectrum power calculation unit for calculating a wave spectrum power in an analysis area included in the detection area, and the received wave from the analysis area
  • a wave spectrum power correction unit that corrects the wave spectrum power based on a relative wave direction that is a traveling direction of the waves in the analysis area with reference to the direction toward the vessel, and based on the corrected wave spectrum power
  • a wave height calculator for calculating the wave height of the waves in the analysis area.
  • the wave spectrum power correction unit calculates the wave spectrum power according to a first correction coefficient calculation formula that is a function of the relative wave direction that takes a minimum value when the relative wave direction is 0 degree.
  • the wave spectrum power is corrected using a first correction coefficient obtained by substituting the relative wave direction of the waves in the analysis area.
  • the first correction coefficient calculation formula takes a minimum value when the relative wave direction is 180 degrees.
  • the wave spectrum power correction unit multiplies the wave spectrum power by the first correction coefficient obtained by the first correction coefficient calculation expression represented by the following expression (1) to thereby calculate the wave spectrum power. Correct the spectral power.
  • ⁇ is the first correction coefficient
  • ⁇ is the relative wave direction
  • A, B, and C are correction coefficient parameters.
  • the wave spectrum power correction unit obtains the second correction coefficient obtained based on the first correction coefficient obtained by the first correction coefficient calculation formula represented by the following formula (2):
  • the wave spectrum power is corrected by multiplying the wave spectrum power.
  • i is a natural number given corresponding to the direction
  • ⁇ i is the first correction coefficient calculated corresponding to the direction of the wave
  • ⁇ i is the relative wave direction calculated corresponding to the direction of the wave.
  • A, B, and C are correction coefficient parameters, respectively.
  • the correction coefficient parameter is set according to the direction in which the waves arrive.
  • the correction coefficient parameter is set according to the wave period.
  • the wave height calculation device further includes a storage unit that stores a plurality of correction coefficient parameters.
  • the wave height calculation apparatus further includes a correction coefficient parameter calculation processing unit that calculates a plurality of correction coefficient parameters.
  • the correction coefficient parameter calculation processing unit calculates a wave spectrum for each area, which is a wave spectrum power included in each of the plurality of data acquisition areas set in the detection area.
  • a relative wave for each area that calculates a relative wave direction for each area, which is a traveling direction of waves in each of the data acquisition areas, based on a direction from each data acquisition area to the receiver.
  • the sample point specified by the direction calculation unit, the wave wave power per area and the relative wave direction per area corresponds to the first axis corresponding to the wave wave power per area and the relative wave direction per area.
  • a graph generating unit for generating a wave power spectrum by relative wave direction by plotting on coordinates having the second axis, and the wave by relative wave direction; It has a correction coefficient for parameter calculation unit for calculating a parameter for the correction coefficient based on Spectral power graph, a.
  • the plurality of wave spectrum powers per area obtained while the receiver rotates 360 degrees along the horizontal plane are each divided by the wave power spectrums of areas having the largest value among them.
  • the graph generation unit plots the normalized wave spectrum power for each area and the relative wave direction for each area on the coordinates to plot the wave spectrum power graph for each relative wave direction. Generate.
  • a radar apparatus is reflected by a transmitter that transmits a transmission wave, and the transmission wave is reflected by waves generated on a water surface in a detection area.
  • a wave height calculation method is obtained from a receiver that receives an echo that is transmitted and reflected by waves generated on a water surface in a detection area.
  • a wave height calculation method for calculating the wave height of a wave based on an echo signal, the step of calculating wave spectrum power in an analysis area included in the detection area, and a direction from the analysis area toward the receiver The wave spectrum power is corrected based on the relative wave direction that is the traveling direction of the waves in the analysis area with reference to the wave, and the waves in the analysis area are corrected based on the corrected wave spectrum power. Calculating a wave height.
  • the wave height can be accurately calculated regardless of the relative wave direction of the wave that is the wave height calculation target.
  • FIG. 1 is a block diagram of a radar apparatus according to an embodiment of the present invention. It is a block diagram which shows the structure of the wave height calculation process part shown in FIG. It is a figure which shows typically an example of the echo image produced
  • a wave height calculation processing unit 10 as a wave height calculation device according to an embodiment of the present invention a radar apparatus 1 including the wave height calculation processing unit 10, and a wave height calculation method will be described with reference to the drawings.
  • the present invention can be widely applied to a wave height calculation device that calculates the height of a wave generated on the water surface, a radar device including the wave height calculation device, and a wave height calculation method.
  • FIG. 1 is a block diagram of a radar apparatus 1 according to an embodiment of the present invention.
  • the radar apparatus 1 according to the present embodiment is provided in a ship such as a fishing boat. According to the radar apparatus 1, as will be described in detail below, the wave height that is the height of the waves generated on the water surface can be accurately calculated.
  • the radar device 1 includes an antenna unit 2, an anemometer 3, a wave height calculation processing unit 10 as a wave height calculation device, and a display 4.
  • the antenna unit 2 includes an antenna 5, a receiving unit 6, and an A / D conversion unit 7.
  • the antenna 5 is a radar antenna capable of transmitting a pulsed radio wave as a highly directional transmission wave.
  • the antenna 5 is configured to receive a reflected wave from a target (a wave in the case of the present embodiment). That is, the antenna 5 functions as a transmitter that transmits a transmission wave and a receiver that receives a reflected wave of the transmitted transmission wave as a reception wave.
  • the radar apparatus 1 measures the time from when a pulsed radio wave is transmitted until the reflected wave is received. Thereby, the radar apparatus 1 can detect the distance r to the target.
  • the antenna 5 is configured to be able to rotate 360 ° on a horizontal plane.
  • the antenna 5 is configured to repeatedly transmit and receive radio waves while changing the transmission direction of pulsed radio waves (for example, changing the antenna angle). With the above configuration, the radar apparatus 1 can detect a target on a plane around the ship over 360 °.
  • 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 360 ° while transmitting and receiving radio waves is called “scan”.
  • the receiving unit 6 detects and amplifies an echo signal obtained from an echo received by the antenna 5.
  • the reception unit 6 outputs the amplified echo signal to the A / D conversion unit 7.
  • the A / D converter 7 samples an analog echo signal and converts it to digital data composed of a plurality of bits. This digital data is echo data.
  • the echo data includes data specifying the intensity of the echo signal obtained from the reflected wave received by the antenna 5.
  • the A / D converter 7 outputs the echo data to the wave height calculation processor 10.
  • the anemometer 3 measures the wind speed at sea (sea wind speed, surface wind speed), and is installed in the ship.
  • the anemometer 3 outputs data relating to the measured wind speed to the wave height calculation processing unit 10.
  • the wave height calculation processing unit 10 calculates the wave height (wave height) based on the echo data output from the antenna unit 2.
  • the wave height calculation processing unit 10 outputs data relating to the calculated wave height to the display 4.
  • the wave height calculation processing unit 10 does not calculate the wave height when the wind speed obtained by the anemometer 3 is equal to or less than a predetermined value. This is because when the wind speed is low, the wave height tends to be small, and it is difficult to calculate an accurate wave height.
  • the configuration and operation of the wave height calculation processing unit 10 will be described later in detail.
  • the display 4 displays data relating to the wave height output from the wave height calculation processing unit 10 (for example, a numerical value of the wave height). Thereby, the user can know the wave height at sea.
  • FIG. 2 is a block diagram showing a configuration of the wave height calculation processing unit 10 shown in FIG.
  • the wave height calculation processing unit 10 includes an image generation unit 11, an analysis area setting unit 12, a frequency analysis unit 13, a wave spectrum power calculation unit 14, a relative wave direction calculation unit 15, a correction coefficient calculation unit 16, and a wave A spectrum power correction unit 17 and a wave height calculation unit 18 are provided.
  • the wave height calculation processing unit 10 includes a hardware processor 8 (for example, CPU, FPGA, etc.) and a device such as a nonvolatile memory.
  • a hardware processor 8 for example, CPU, FPGA, etc.
  • the CPU reads the program from the nonvolatile memory and executes it, thereby causing the wave height calculation processing unit 10 to function as the image generation unit 11, the analysis area setting unit 12, the frequency analysis unit 13, the wave spectrum power calculation unit 14, and the like. be able to.
  • FIG. 3 is a diagram schematically illustrating an example of the echo image P generated by the image generation unit 11.
  • the image generation unit 11 generates an echo image P based on the echo data output from the antenna unit 2.
  • the echo image P is generated every time the antenna 5 rotates 360 ° (that is, every scan). Note that the example shown in FIG. 3 shows an example in which the wavy line w is relatively clearly reflected.
  • the analysis area setting unit 12 sets an analysis area Z for the echo image P (see FIG. 3).
  • the analysis area Z may be set as needed by the user, or may be set in advance when the apparatus is shipped.
  • analysis area Z as shown in FIG. 3 is set in front of the ship S 0.
  • the region where the analysis area Z is set at a region other than the ship S 0 backward are preferred. This is because the ship behind the region is because accurate wave analyzed by the undertow of the ship S 0 becomes difficult.
  • the frequency analysis unit 13 performs Fourier transform on the echo image in the analysis area, which is an echo image in the analysis area Z, obtained for each scan, and calculates the frequency spectrum S (f). Further, the frequency analysis unit 13 calculates the wave direction of the waves based on the frequency spectrum S (f) obtained by Fourier transform. In addition, since the calculation method of a wave direction is known, the description is abbreviate
  • the wave spectrum power calculation unit 14 calculates the zero-order moment m 0 of the wave spectrum in the analysis area Z using the following equation (3) based on the frequency spectrum S (f) obtained from the echo image in the analysis area. To do.
  • the relative wave direction calculation unit 15 refers to FIG. 3, and the wave direction dw of the waves in the analysis area Z calculated by the frequency analysis unit 13 based on the direction ds from the analysis area Z toward the ship S 0 as a reference.
  • the relative wave direction ⁇ defined by is calculated.
  • the relative wave direction calculation unit 15 calculates the relative wave direction ⁇ with the counterclockwise direction in FIG. 3 as the positive direction.
  • the correction coefficient calculation unit 16 substitutes the relative wave direction ⁇ calculated by the relative wave direction calculation unit 15 into the following equation (4), thereby calculating the wave spectrum power ⁇ calculated by the wave spectrum power calculation unit 14.
  • a correction coefficient ⁇ (first correction coefficient) for correcting m 0 is calculated.
  • Wave spectrum power correcting unit 17 the wave spectrum power ⁇ m 0 calculated by ocean wave spectrum power calculation unit 14, the wave spectrum power ⁇ m 0 by multiplying a correction coefficient calculated by the correction coefficient calculating unit 16 beta After correction, the corrected wave spectrum power ⁇ m 0 — new is calculated.
  • the wave height calculation unit 18 calculates the wave height H 1/3 of the wave in the analysis area Z by multiplying the corrected wave spectrum power ⁇ m 0_new calculated by the wave spectrum power correction unit 17 by the fixed coefficient ⁇ . Specifically, the wave height calculation unit 18 calculates the wave height H 1/3 based on the following equation (5).
  • the relative wave direction ⁇ is 0 degree (specifically, referring to FIG. 3, when ds and dw are in the same direction), or when 180 degrees (ds and dw are in opposite directions) ) Is a case where the transmission wave from the radar apparatus 1 is transmitted vertically toward the wave line w.
  • the wave spectrum power becomes relatively strong. Then, even if the wave height is calculated using the wave spectrum power as in the prior art, the wave height can be calculated relatively accurately.
  • the case where the relative wave direction ⁇ is 90 degrees or 270 degrees is a case where the transmission wave from the radar apparatus 1 is transmitted to the side of the wave line w.
  • the transmission wave only hits a relatively narrow range in the wave, and the wave spectrum power is calculated to be weak.
  • the tendency that the wave height is calculated to be lower than when the relative wave direction is 0 degree or 180 degrees becomes higher, and the wave height cannot be calculated accurately.
  • the correction coefficient ⁇ when the relative wave direction ⁇ is around 0 degrees or around 180 degrees, the correction coefficient ⁇ is a relatively small value.
  • the correction coefficient ⁇ when the relative wave direction ⁇ is around 90 degrees or around 270 degrees, the correction coefficient ⁇ is larger than when the relative wave direction is around 0 degrees or around 180 degrees. That is, by using the equation (4), even if the wave spectrum power ⁇ m 0 is calculated to be weak, the wave spectrum power ⁇ m 0 is corrected by the correction coefficient ⁇ so as to increase.
  • FIG. 4 is a flowchart for explaining each step performed when a preliminary experiment for determining the correction coefficient parameters A, B, and C for calculating the correction coefficient ⁇ shown in Expression (4) is performed. is there.
  • each step of the preliminary experiment for determining the correction coefficient parameters A, B, and C will be described with reference to FIG.
  • FIG. 5 is a diagram showing a position relative to the ship S 0 of wave spectrum power ⁇ m 0 is ⁇ data acquisition area Z1 calculated Z7 in preliminary experiments.
  • step S1 frequency analysis of echo signals included in each of the data acquisition areas Z1 to Z7 is performed based on the echoes obtained from each of the plurality of data acquisition areas Z1 to Z7 shown in FIG.
  • step S1 the echo images included in each of the data acquisition areas Z1 to Z7 are Fourier transformed to convert each echo image into a frequency spectrum.
  • step S2 wave directions in the respective data acquisition areas Z1 to Z7 are calculated based on the frequency spectrum generated in step S1.
  • step S2 the direction toward the ship S 0 from the data acquisition area Z1 ⁇ Z7 as a reference, a wave of wave direction in each data acquisition areas Z1 ⁇ Z7, the data acquisition area Z1 ⁇ Calculated as the relative wave direction of waves in Z7.
  • step S3 the wave spectrum power ⁇ m 0 in each of the data acquisition areas Z1 to Z7 is calculated before or after step S2 or in parallel with step S2.
  • the calculation method of the wave spectrum power ⁇ m 0 is the same as that of the wave spectrum power calculation unit 14.
  • step S4 ocean wave spectrum power ⁇ m 0 in each data acquisition areas Z1 ⁇ Z7 is normalized. Specifically, the wave spectrum power ⁇ m 0 in each of the data acquisition areas Z1 to Z7 is divided by the wave spectrum power of the data acquisition area having the highest wave spectrum power among the seven data acquisition areas Z1 to Z7. Is done. As a result, the wave spectrum power of each of the data acquisition areas Z1 to Z7 is normalized so that the wave spectrum power of the area having the highest wave spectrum power among the data acquisition areas Z1 to Z7 is 1. In addition, the wave spectrum power normalized in this way is hereinafter referred to as normalized wave spectrum power.
  • step S5 seven normalized wave spectrum powers acquired by performing steps S1 to S4 are calculated at each of a plurality of timings (ie, a plurality of scans).
  • step S6 as shown in FIG. 6, each of a large number of sample points having the normalized wave spectrum power calculated in step S5 as information has the relative wave direction ⁇ as the x axis and the normalized wave spectrum. Plotted on Cartesian coordinates with power as y-axis. Thereby, the normalized wave spectrum power graph classified by relative wave direction is generated.
  • the normalized wave spectrum power graph classified by relative wave direction is simply referred to as a scatter diagram SP.
  • sample points indicated by circles in FIG. 6 are sample points obtained when the wind class of the wind speed obtained by the anemometer is 5.
  • the sample points indicated by square marks are sample points obtained when the wind speed class obtained by the anemometer is 6.
  • the sample points indicated by triangles are sample points obtained when the wind speed class obtained by the anemometer is 7.
  • the sample points indicated by x are sample points obtained when the wind speed of the wind speed obtained by the anemometer is 8 or more.
  • step S7 correction coefficient parameters A, B, and C are calculated based on the scatter diagram SP generated in step S6. Specifically, in step S6, the sum of squares of the residuals between the above-described expression (4), the right side denominator (A + B cos ⁇ + C cos 2 ⁇ ), and each sample point constituting the scatter diagram SP is minimized. Then, correction coefficient parameters A, B, and C are calculated. That is, in step S6, correction coefficient parameters A, B, and C at A + B cos ⁇ + C cos 2 ⁇ are calculated by the least square method. Thus, the correction coefficient parameter can be calculated based on the wave spectrum power in each data acquisition area actually obtained. In FIG. 6, a graph indicating a mathematical formula (A + B cos ⁇ + C cos 2 ⁇ ) in which the correction coefficient parameters A, B, and C are calculated by the least square method is superimposed on the scatter diagram SP.
  • the wave height calculation processing unit 10 of the radar apparatus 1 corrects the wave spectrum power ⁇ m 0 based on the relative wave direction ⁇ to calculate the corrected wave spectrum power ⁇ m 0_new ,
  • the wave height H 1/3 is calculated by multiplying the corrected wave spectrum power ⁇ m 0 — new by a fixed coefficient ⁇ .
  • the wave height H 1/3 can be accurately calculated by correcting the wave spectrum power ⁇ m 0 based on the relative wave direction ⁇ as in the wave height calculation processing unit 10 of the present embodiment.
  • the wave height calculation processing unit 10 the wave height can be accurately calculated regardless of the relative wave direction ⁇ of the wave to be wave height calculation target.
  • a correction coefficient calculation formula for correcting the wave spectrum power ⁇ m 0 used in the wave height calculation processing unit 10 a formula that takes a minimum value when the relative wave direction ⁇ is 0 degree is employed.
  • the correction coefficient ⁇ can be appropriately set by adopting a correction coefficient calculation formula that reduces the correction coefficient when the wave spectrum power ⁇ m 0 tends to increase.
  • a correction coefficient calculation formula for correcting the wave spectrum power ⁇ m 0 used in the wave height calculation processing unit 10 a formula that takes a minimum value when the relative wave direction ⁇ is 180 degrees is adopted.
  • the correction coefficient ⁇ can be appropriately set by adopting a correction coefficient calculation formula that reduces the correction coefficient when the wave spectrum power ⁇ m 0 tends to increase.
  • Equation (4) the correction coefficient ⁇ increases when the wave spectrum power ⁇ m 0 tends to be small, specifically when the relative wave direction ⁇ is around 90 degrees or around 270 degrees. Therefore, the wave height calculation processing unit 10 can set the correction coefficient ⁇ more appropriately.
  • the correction coefficient parameters A, B, and C for calculating the correction coefficient ⁇ are preliminarily determined by preliminary experiments before the wave height calculation unit 18 calculates the wave height H 1/3. These values are set and stored in the correction coefficient calculation unit 16. This eliminates the need to calculate the correction coefficient parameters A, B, and C in parallel with the wave height calculation by the wave height calculation unit 18, thereby reducing the calculation load on the wave height calculation processing unit 10.
  • the radar apparatus 1 it is possible to provide a radar apparatus including a wave height calculation processing unit that can accurately calculate the wave height regardless of the wave traveling direction in which the wave height H 1/3 is calculated.
  • FIG. 7 is a block diagram of a wave height calculation processing unit 10a of a radar apparatus according to a modification.
  • the correction coefficient ⁇ is calculated using the correction coefficient parameters A, B, and C calculated based on preliminary experiments performed in advance.
  • correction coefficient parameters A, B, and C are calculated as needed in parallel with the calculation of the wave height by the wave height calculation unit 18.
  • the wave height calculation processing unit 10a includes a correction coefficient parameter calculation processing unit 20 in addition to the components included in the wave height calculation processing unit 10 of the above-described embodiment.
  • FIG. 8 is a block diagram showing the correction coefficient parameter calculation processing unit 20 shown in FIG.
  • the correction coefficient parameter calculation processing unit 20 includes a data acquisition area setting unit 21, an area frequency analysis unit 22, an area wave spectrum power calculation unit 23, a normalization unit 24, and an area relative wave direction calculation unit. 25, a graph generation unit 26, and a correction coefficient parameter calculation unit 27.
  • the data acquisition area setting unit 21 sets a plurality of data acquisition areas that are areas in which the wave spectrum power is calculated. For example, the data acquisition area setting unit 21 sets each of the data acquisition areas Z1 to Z7 shown in FIG. 5 as a data acquisition area.
  • the area-by-area frequency analysis unit 22 performs Fourier transform on the echo images in the data acquisition areas Z1 to Z7 obtained for each scan, and calculates a frequency spectrum. Further, the frequency analysis unit 22 for each area calculates the wave direction of the waves in each of the data acquisition areas Z1 to Z7 based on the frequency spectrum obtained by the Fourier transform.
  • the area-specific wave spectrum power calculation unit 23 calculates the area-specific wave spectrum power ⁇ m 0 in each data acquisition area based on the frequency spectrum obtained corresponding to each data acquisition area.
  • the wave wave power per area ⁇ m 0 is calculated in the same manner as in the above embodiment.
  • the normalizing unit 24 normalizes the wave spectrum power ⁇ m 0 for each area. Specifically, the normalization unit 24, by a plurality of areas each wave spectral power ⁇ m 0 obtained in the first scan, respectively, divided by the area per wave spectral power ⁇ m 0 most large value of them Normalize the wave spectrum power ⁇ m 0 per area. Thereby, the wave spectrum power for each area is normalized so that the wave spectrum power of the area having the highest wave spectrum power among the plurality of areas becomes 1. The normalization unit 24 normalizes the wave spectrum power ⁇ m 0 for each area for each scan.
  • the relative wave direction calculation unit 25 for each area is based on the wave direction calculated for each data acquisition area and the direction from each data acquisition area toward the ship. Calculate the relative wave direction of the waves.
  • the calculation method of the relative wave direction is the same as that of the relative wave direction calculation unit 15 of the above embodiment.
  • the graph generation unit 26 plots the sample points specified by the normalized wave spectrum power per area and the relative wave direction for each area obtained in the latest in the scatter diagram stored in the graph generation unit 26. Update the scatter plot to generate a new scatter plot.
  • the scatter diagram generated by the graph generation unit 26 is the same as that shown in FIG.
  • the correction coefficient parameter calculation unit 27 calculates correction coefficient parameters A, B, and C each time a new scatter diagram is generated by the graph generation unit 26. Specifically, the correction coefficient parameter calculation unit 27 uses correction coefficient parameters so that the sum of squares of the residuals between the expression represented by A + B cos ⁇ + C cos 2 ⁇ and each sample point constituting the updated scatter diagram is minimized. Parameters A, B, and C are calculated.
  • the wave spectrum power correction unit 17 uses the correction coefficient parameters A, B, and C calculated after taking into account the normalized wave spectrum power obtained at the latest timing.
  • the wave spectrum power in the analysis area Z is corrected based on the obtained correction coefficient ⁇ .
  • the wave height calculation unit 18 calculates the wave height in the analysis area Z based on the corrected wave spectrum power.
  • the correction coefficient parameter A calculated in consideration of the normalized wave spectrum power obtained at the latest timing. , B, and C are used to correct the wave spectrum power in the analysis area Z based on the correction coefficient ⁇ , and the wave height H 1/3 in the analysis area Z is calculated based on the corrected wave spectrum power. Is done. That is, according to the wave height calculation processing unit 10a of the present modification, the wave height H 1/3 in the analysis area Z can be calculated based on the latest data, and therefore the wave height H 1/3 can be calculated more accurately. .
  • the wave height calculation processing unit 10a as the configuration requirements of the correction coefficient parameter calculation processing unit 20, the area-specific wave spectrum power calculation unit 23, the area-specific relative wave direction calculation unit 25, the graph generation unit 26, and the correction coefficient parameter A calculation unit 27 is provided. Thereby, according to the wave height calculation process part 10a, the specific structure for calculating the wave height H1 / 3 in the analysis area Z can be provided based on the newest data.
  • the wave height calculation processing unit 10a normalizes a plurality of wave spectrum powers obtained for each scan and calculates correction coefficient parameters A, B, and C based on the normalized wave spectrum powers. Generates a scatter plot.
  • the wave height calculation processing unit 10a of this modification in the case of the wave height calculation processing unit 10a of this modification, the plurality of wave spectrum powers obtained for each scan are normalized for each scan, and the correction coefficient parameters A, A scatter diagram for calculating B and C is generated. In this way, the magnitude of the wave spectrum power caused by the wind speed can be made uniform, so that it is not necessary to generate a scatter diagram for each wind speed. Thereby, it is possible to reduce the calculation load on the wave height calculation processing unit 10a for calculating the correction coefficient parameters A, B, and C.
  • FIG. 9 is a block diagram showing a configuration of the wave height calculation processing unit 10b of the radar apparatus according to the modification.
  • the wave height calculation processing unit 10 that calculates the wave height of waves arriving from one direction has been described as an example, but the present invention is not limited thereto.
  • a wave height calculation processing unit 10b that can accurately calculate the wave height of each wave coming from a plurality of directions can be configured.
  • the configurations and operations of the frequency analysis unit 13a, the relative wave direction calculation unit 15a, and the correction coefficient calculation unit 16a are different from those in the above embodiment. Below, a different part from the said embodiment is demonstrated and description is abbreviate
  • the frequency analysis unit 13a calculates the frequency spectrum S (f) in the same manner as in the above embodiment, and then calculates the wave direction of each wave arriving from a plurality of directions based on the frequency spectrum S (f).
  • the relative wave direction calculation unit 15a calculates the relative wave direction ⁇ i of waves arriving from each direction.
  • the correction coefficient calculation unit 16a calculates a wave-by-wave correction coefficient ⁇ i (first correction coefficient) for each wave from each direction. Specifically, the correction coefficient calculation unit 16a calculates a wave-by-wave correction coefficient ⁇ i corresponding to waves from each direction based on the following equation (6). Then, the correction coefficient calculating unit 16a calculates the correction coefficient ⁇ (second correction coefficient) by substituting the wave-by-wave correction coefficient ⁇ i obtained corresponding to each direction into the following equation (7). . That is, Expression (6) is a first correction coefficient calculation expression, and Expression (7) is a second correction coefficient calculation expression.
  • the wave spectrum power correction unit 17 multiplies the wave spectrum power ⁇ m 0 calculated by the wave spectrum power calculation unit 14 by the correction coefficient ⁇ calculated by the correction coefficient calculation unit 16a to thereby generate the wave spectrum power ⁇ m.
  • the corrected wave spectrum power ⁇ m 0 — new is calculated by correcting 0 .
  • the wave height calculation processing unit 10b can accurately calculate the wave height regardless of the relative wave direction ⁇ of the wave that is the wave height calculation target, as in the case of the above-described embodiment.
  • equation (8) is used instead of equation (7) to correct the wave size from each direction (for example, the peak value of the wave spectrum power).
  • the coefficient ⁇ can be calculated. Thereby, the wave height can be calculated more accurately.
  • ⁇ i is a weighting coefficient, which is a numerical value determined in accordance with the magnitude of the wave from each direction. A larger value is set as the wave is larger, and a smaller value is set as the wave is smaller. .
  • FIG. 10 is a block diagram showing a configuration of the wave height calculation processing unit 10c of the radar apparatus according to the modification.
  • the correction coefficient ⁇ is calculated using the correction coefficient parameters A, B, and C calculated using the same scatter diagram SP regardless of the arrival direction of the waves.
  • this is not restrictive.
  • correction coefficient parameters A i , B i , and C i calculated based on a scatter diagram SP i (not shown) generated for each direction in which waves arrive are stored in the correction coefficient calculation unit 16b.
  • the correction coefficient calculation unit 16b calculates a wave-by-wave correction coefficient ⁇ i (first correction coefficient) corresponding to the waves from each direction based on the following equation (9).
  • the correction coefficient calculation unit 16b calculates the correction coefficient ⁇ (second correction coefficient) using Expression (7), as in the case of the modification described with reference to FIG.
  • the wave spectrum power correction unit 17 multiplies the wave spectrum power ⁇ m 0 calculated by the wave spectrum power calculation unit 14 by the correction coefficient ⁇ calculated by the correction coefficient calculation unit 16b to thereby generate the wave spectrum power ⁇ m.
  • the corrected wave spectrum power ⁇ m 0 — new is calculated by correcting 0 .
  • the wave height calculation processing unit 10c can also accurately calculate the wave height regardless of the relative wave direction ⁇ of the wave that is the wave height calculation target, as in the case of the above embodiment.
  • the correction coefficient ⁇ can be calculated in consideration of the waves from these directions. Therefore, according to this modification, the wave height can be calculated more accurately. Moreover, according to the present modification, the correction coefficient ⁇ can be calculated using the correction coefficient parameters A i , B i , and C i calculated for each azimuth, so that the wave height can be calculated more accurately.
  • FIG. 11 is a block diagram showing a configuration of the wave height calculation processing unit 10d of the radar apparatus according to the modification.
  • the wave height is calculated using the correction coefficient ⁇ calculated based on the same correction coefficient parameters A, B, and C regardless of the wave period for which the wave height is to be calculated.
  • the present invention is not limited to this, and different correction coefficient parameters may be used depending on the wave period.
  • the correction coefficient calculation unit 16c has a wave having a period of 8 seconds or less and a period of 8 seconds or more. Different correction coefficient parameters may be used for the waves.
  • waves include wind and swell, and these characteristics are different from each other.
  • the period of wind is approximately 8 seconds or less
  • the period of undulation is approximately 8 seconds or more. That is, as described above, the correction coefficient ⁇ corresponding to each of the wind and swell having different characteristics can be calculated by using different correction coefficient parameters depending on the wave period. Accordingly, since an appropriate correction coefficient ⁇ can be calculated according to each of waves having different characteristics (for example, wind and swell), the wave height can be calculated more accurately according to the type of the waves.
  • the wave period that is a wave height calculation target is calculated by the frequency analysis unit 13, and the correction coefficient calculation unit 16 c uses any correction coefficient parameter depending on the wave period value calculated by the frequency analysis unit 13. Decide whether to use.
  • Wave height calculation processing unit (wave height calculation device) 14 Wave Spectrum Power Calculation Unit 17 Wave Spectrum Power Correction Unit 18 Wave Height Calculation Unit

Abstract

La présente invention a pour but de calculer précisément une hauteur d'onde indépendamment de la direction de progression de l'onde dont la hauteur d'onde est en cours de calcul. À cet effet, la présente invention concerne un dispositif de calcul de hauteur d'onde (10) qui comporte : une unité de calcul de puissance de spectre d'onde (14) qui calcule une puissance de spectre d'onde √m0 à l'intérieur d'une zone d'analyse incluse dans une zone de détection ; une unité de correction de puissance de spectre d'onde (17) qui corrige la puissance de spectre d'onde √m0 sur la base d'une direction d'onde relative θ, qui est la direction de progression des ondes dans la zone d'analyse par rapport à une direction allant de la zone d'analyse vers un récepteur d'onde ; et une unité de calcul de hauteur d'onde (18) qui calcule une hauteur d'onde H1/3 des ondes dans la zone d'analyse sur la base de la puissance de spectre d'onde corrigée √m0_nouvelle.
PCT/JP2017/009496 2016-04-11 2017-03-09 Dispositif de calcul de hauteur d'onde, dispositif radar et procédé de calcul de hauteur d'onde WO2017179344A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2018511930A JP6676151B2 (ja) 2016-04-11 2017-03-09 波高算出装置、レーダ装置、及び波高算出方法

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2016-079175 2016-04-11
JP2016079175 2016-04-11

Publications (1)

Publication Number Publication Date
WO2017179344A1 true WO2017179344A1 (fr) 2017-10-19

Family

ID=60041620

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/009496 WO2017179344A1 (fr) 2016-04-11 2017-03-09 Dispositif de calcul de hauteur d'onde, dispositif radar et procédé de calcul de hauteur d'onde

Country Status (2)

Country Link
JP (1) JP6676151B2 (fr)
WO (1) WO2017179344A1 (fr)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109141376A (zh) * 2018-08-06 2019-01-04 上海海事大学 一种基于单目视觉的浪向检测方法
CN110823190A (zh) * 2019-09-30 2020-02-21 广州地理研究所 基于随机森林的岛礁浅海水深预测方法
CN111965628A (zh) * 2020-08-11 2020-11-20 中国人民解放军91550部队 一种垂直出水航行体瞬时波浪参数估计方法
CN112197749A (zh) * 2020-09-30 2021-01-08 国家海洋环境预报中心 一种波浪浮标有效波高的交叉校准方法及装置
WO2021100402A1 (fr) 2019-11-21 2021-05-27 古野電気株式会社 Dispositif radar à semi-conducteurs

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003185743A (ja) * 2001-12-19 2003-07-03 Yokohama Tlo Co Ltd 波高算出装置、波高算出方法、記録媒体、及び船舶
JP2003307563A (ja) * 2002-04-16 2003-10-31 Michio Nobeoka レーダ波高計測装置、レーダ波高補正係数作成方法、及びレーダ波高補正係数作成装置
US20100315284A1 (en) * 2009-09-02 2010-12-16 Trizna Dennis B Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents
JP2015014472A (ja) * 2013-07-03 2015-01-22 古野電気株式会社 エコー信号処理装置、波浪レーダ装置、エコー信号処理方法、及びエコー信号処理プログラム

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003185743A (ja) * 2001-12-19 2003-07-03 Yokohama Tlo Co Ltd 波高算出装置、波高算出方法、記録媒体、及び船舶
JP2003307563A (ja) * 2002-04-16 2003-10-31 Michio Nobeoka レーダ波高計測装置、レーダ波高補正係数作成方法、及びレーダ波高補正係数作成装置
US20100315284A1 (en) * 2009-09-02 2010-12-16 Trizna Dennis B Method and apparatus for coherent marine radar measurements of properties of ocean waves and currents
JP2015014472A (ja) * 2013-07-03 2015-01-22 古野電気株式会社 エコー信号処理装置、波浪レーダ装置、エコー信号処理方法、及びエコー信号処理プログラム

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YUKIHARU HISAKI ET AL.: "Observations of sea wave by HF Ocean Radar", PROCEEDINGS OF THE 12TH JAPANESE CONFERENCE ON REMOTE SENSING, 7 May 1992 (1992-05-07), pages 63 - 66, ISSN: 09163611 *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109141376A (zh) * 2018-08-06 2019-01-04 上海海事大学 一种基于单目视觉的浪向检测方法
CN109141376B (zh) * 2018-08-06 2021-02-26 上海海事大学 一种基于单目视觉的浪向检测方法
CN110823190A (zh) * 2019-09-30 2020-02-21 广州地理研究所 基于随机森林的岛礁浅海水深预测方法
CN110823190B (zh) * 2019-09-30 2020-12-08 广州地理研究所 基于随机森林的岛礁浅海水深预测方法
WO2021100402A1 (fr) 2019-11-21 2021-05-27 古野電気株式会社 Dispositif radar à semi-conducteurs
CN111965628A (zh) * 2020-08-11 2020-11-20 中国人民解放军91550部队 一种垂直出水航行体瞬时波浪参数估计方法
CN111965628B (zh) * 2020-08-11 2023-09-12 中国人民解放军91550部队 一种垂直出水航行体瞬时波浪参数估计方法
CN112197749A (zh) * 2020-09-30 2021-01-08 国家海洋环境预报中心 一种波浪浮标有效波高的交叉校准方法及装置
CN112197749B (zh) * 2020-09-30 2021-06-29 国家海洋环境预报中心 一种波浪浮标有效波高的交叉校准方法及装置

Also Published As

Publication number Publication date
JPWO2017179344A1 (ja) 2019-02-14
JP6676151B2 (ja) 2020-04-08

Similar Documents

Publication Publication Date Title
WO2017179344A1 (fr) Dispositif de calcul de hauteur d'onde, dispositif radar et procédé de calcul de hauteur d'onde
US10324190B2 (en) Wind measuring apparatus
JP2012017995A (ja) 閾値設定方法、物標探知方法、閾値設定装置、物標探知装置、閾値設定プログラム、および物標探知プログラム
JP2011053028A (ja) ドップラーレーダ装置及びドップラー速度算出方法
US11249185B2 (en) Signal processing device and radar apparatus
JP6404212B2 (ja) 表層潮流推定装置、レーダ装置、表層潮流推定方法、及び表層潮流推定プログラム
JP5767002B2 (ja) 超音波送受信装置、および魚量検出方法
JP5664869B2 (ja) 測定装置、測定システム、測定方法、及びプログラム
JP6095899B2 (ja) 物標運動推定装置、物標運動推定方法、およびレーダ装置
JP6420136B2 (ja) 波高算出装置
JP2013246022A (ja) 物標探知装置、レーダ装置、物標探知方法、および物標探知プログラム
JP6154219B2 (ja) エコー信号処理装置、波浪レーダ装置、エコー信号処理方法、及びエコー信号処理プログラム
JP6043083B2 (ja) 物標運動推定装置、物標運動推定方法、およびレーダ装置
JP6289672B2 (ja) 合成開口レーダ信号処理装置及び合成開口レーダ信号処理プログラム
JP5196959B2 (ja) レーダ装置
JPH0228116B2 (fr)
JP6362523B2 (ja) 風計測装置
JP2015232509A (ja) 信号処理装置、レーダ装置、信号処理方法、及び信号処理プログラム
JP3777576B2 (ja) 位置推定方法および装置
JP6339074B2 (ja) 海象検出装置、レーダ装置、海象検出方法、および、プログラム
JP7269784B2 (ja) 水中探知装置、水中探知方法およびプログラム
JP2000338238A (ja) レーダ装置
CN116819534A (zh) 岸基超高频雷达风速反演方法及系统
JP2006071363A (ja) レーダ装置
JP2013088329A (ja) 超音波距離画像生成装置

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 2018511930

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17782162

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 17782162

Country of ref document: EP

Kind code of ref document: A1