Effective wave height fusion method and system based on multi-source satellite radar altimeter data

Abstract

The present invention relates to a method and system for effective wave height fusion based on multi-source satellite radar altimeter data, including: obtaining observation field effective wave height data and background field effective wave height data; calculating observation field variance and background field variance respectively according to observation field effective wave height data and background field effective wave height data; calculating observation operator according to preset influence radius and longitude and latitude distance, longitude and latitude distance is the distance between longitude and latitude of background field effective wave height grid point and longitude and latitude of altimeter observation position; calculating projected effective wave height according to background field effective wave height data and observation operator; calculating fused effective wave height by observation field effective wave height data, background field effective wave height data, observation field variance, background field variance, observation operator and projected effective wave height. The present invention constructs effective wave height analysis field by using variational fusion method for observation field effective wave height data and background field effective wave height data, thereby improving the accuracy of fused effective wave height field.

Description

Technical Field

[0001] The present invention relates to the field of ocean remote sensing and applications, and in particular to a method and system for fusing effective wave heights based on multi-source satellite radar altimeter data. Background Art

[0002] Significant wave height is one of the main parameters describing the characteristics of ocean waves and plays an important role in ocean research and marine environmental forecasting. Significant wave height is an important factor that must be considered for production activities such as marine transportation, marine aquaculture, maritime salvage and marine engineering.

[0003] Yang Jinsong et al. used the inverse distance weighting method, Kriging method and continuous correction method to fuse the effective wave height data of GFO, Jason-1 and ENVISAT satellite radar altimeters in the China Sea and adjacent waters. The inverse distance weighting method has a high operating efficiency and requires at least three satellites for effective wave height fusion. Wang Jiuke et al. used the wide-swath sea surface wind field of the satellite microwave scatterometer and the effective wave height field of the satellite altimeter to invent a fusion inversion method for the effective wave height field of waves, achieving a breakthrough in satellite remote sensing observation of waves from point to surface, greatly improving the spatial coverage and spatiotemporal resolution of wave observation. Based on different satellite combination schemes, Han Weixiao used the inverse distance weighting method to fuse T / P, Jason-1 / 2, ENVISAT / RA-2, Cryosat-2, HY-2A, and ENVISAT / ASAR data, with a temporal resolution of 1 day and a spatial resolution of 0.25°×0.25°. The results showed that adding ENVISATASAR wave spectrum data or HY-2A data can improve the accuracy of wave fusion products, which is better than the existing wave remote sensing fusion products. Yang Jungang et al. used inverse distance weighting and used HY-2A, CryoSat-2, SARAL and Sentinel-3 multi-satellite joint observations to obtain wave significant wave height remote sensing fusion products with temporal and spatial resolutions better than 1d and 0.25° in the Arctic Ocean (north of 65°N).

[0004] Foreign research institutions use multi-source satellite remote sensing data to develop global ocean wave remote sensing data set products with various spatiotemporal resolutions and long time series. Among them, the AVISO (Archiving, Validation and Interpretation of Satellite Oceanographic data) product data set uses a low-order equivalent system filtering algorithm and fuses the ocean wave data of multiple altimeters such as TOPEX / POSEIDON, Jason-1 and ERS / ENVISAT to provide a global gridded effective wave height fusion product with a resolution of 1°×1°. The European Medium- and Long-Term Numerical Weather Forecast Product (ECMWF) can provide an effective wave height fusion product with a spatial resolution of 0.125°×0.125°. The existing fusion technology has the problem that the effective wave height observation points of multi-source altimeters are sparse and the accuracy of the fused effective wave height field is low. To address this problem, the present invention proposes an effective wave height fusion method and system based on multi-source satellite altimeter data, which improves the accuracy of the fused effective wave height field. Summary of the invention

[0005] The purpose of the present invention is to provide a method and system for the fusion of effective wave heights based on multi-source satellite radar altimeter data, which uses the effective wave height data of multi-source altimeter data as the observation field and the effective wave height of the numerical model forecast field as the background field, and constructs the effective wave height analysis field through the variational fusion method, thereby improving the accuracy of the fused effective wave height field.

[0006] To achieve the above object, the present invention provides the following solutions:

[0007] A method for effective wave height fusion based on multi-source satellite radar altimeter data, comprising:

[0008] Obtain the effective wave height data of the observation field based on the data collected by the multi-source satellite radar altimeter;

[0009] Obtain the background field significant wave height data based on the numerical model forecast field significant wave height data;

[0010] Calculate the observation field variance according to the observation field effective wave height data, and calculate the background field variance according to the background field effective wave height data;

[0011] Calculate the observation operator according to the preset influence radius and the longitude and latitude distance, wherein the longitude and latitude distance is the distance between the longitude and latitude of the effective wave height grid point of the background field and the longitude and latitude of the altimeter observation position;

[0012] Calculating a projected effective wave height according to the background field effective wave height data and the observation operator, wherein the projected effective wave height is the effective wave height of the background field projected to the altimeter observation position;

[0013] The fused effective wave height is calculated according to the observed field effective wave height data, the background field effective wave height data, the observed field variance, the background field variance, the observation operator and the projected effective wave height.

[0014] The present invention also provides an effective wave height fusion system based on multi-source satellite radar altimeter data, comprising:

[0015] The first data acquisition module is used to acquire the effective wave height data of the observation field according to the data collected by the multi-source satellite radar altimeter;

[0016] The second data acquisition module is used to obtain the background field effective wave height data according to the numerical model forecast field wave height data;

[0017] A variance calculation module, used to calculate the variance of the observation field according to the effective wave height data of the observation field, and to calculate the variance of the background field according to the effective wave height data of the background field;

[0018] An observation operator calculation module, used to calculate the observation operator according to a preset influence radius and a longitude and latitude distance, wherein the longitude and latitude distance is the distance between the longitude and latitude of the effective wave height grid point of the background field and the longitude and latitude of the altimeter observation position;

[0019] A projected effective wave height calculation module, used to calculate the projected effective wave height according to the background field effective wave height data and the observation operator, wherein the projected effective wave height is the effective wave height of the background field projected to the altimeter observation position;

[0020] The fused effective wave height calculation module is used to calculate the fused effective wave height according to the observed field effective wave height data, the background field effective wave height data, the observed field variance, the background field variance, the observation operator and the projected effective wave height.

[0021] According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:

[0022] The present invention provides a method and system for effective wave height fusion based on multi-source satellite radar altimeter data, including: obtaining observation field effective wave height data and background field effective wave height data; calculating observation field variance and background field variance respectively according to observation field effective wave height data and background field effective wave height data; calculating observation operator according to preset influence radius and longitude and latitude distance, longitude and latitude distance is the distance between longitude and latitude of background field effective wave height grid point and longitude and latitude of altimeter observation position; calculating projected effective wave height according to background field effective wave height data and observation operator; calculating fused effective wave height using observation field effective wave height data, background field effective wave height data, observation field variance, background field variance, observation operator and projected effective wave height. The present invention uses multi-source altimeter data effective wave height data as observation field, uses numerical model forecast field effective wave height as background field, constructs effective wave height analysis field through variational fusion method, and improves the accuracy of fused effective wave height field. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0024] Figure 1 A flow chart of a method for fusion of effective wave height based on multi-source satellite radar altimeter data provided in Example 1 of the present invention;

[0025] Figure 2 A sea area division map provided in Example 1 of the present invention;

[0026] Figure 3 The global significant wave height background field map provided in Example 1 of the present invention;

[0027] Figure 4 The altimeter significant wave height data diagram provided in Example 1 of the present invention;

[0028] Figure 5 The fused significant wave height data graph provided in Example 1 of the present invention;

[0029] Figure 6 A schematic diagram of the difference in data before and after global significant wave height fusion provided in Example 1 of the present invention;

[0030] Figure 7 A block diagram of an effective wave height fusion system based on multi-source satellite radar altimeter data provided in Example 2 of the present invention. DETAILED DESCRIPTION

[0031] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0032] The purpose of the present invention is to provide a method and system for the fusion of effective wave heights based on multi-source satellite radar altimeter data, which uses the effective wave height data of multi-source altimeter data as the observation field and the effective wave height of the numerical model forecast field as the background field, and constructs the effective wave height analysis field through the variational fusion method, thereby improving the accuracy of the fused effective wave height field.

[0033] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, the present invention is further described in detail below with reference to the accompanying drawings and specific embodiments.

[0034] Example 1

[0035] See also Figure 1 This embodiment provides a method for fusing effective wave height based on multi-source satellite radar altimeter data, including:

[0036] Step S1: Obtaining effective wave height data of the observation field based on data collected by a multi-source satellite radar altimeter;

[0037] Step S2: Obtaining background field effective wave height data according to the numerical model forecast field wave height data;

[0038] Step S3: calculating the observation field variance according to the observation field effective wave height data, and calculating the background field variance according to the background field effective wave height data;

[0039] The calculating the observation field variance according to the observation field effective wave height data specifically includes:

[0040] Calculate the first average value of the effective wave height data of the observation field; the expression of the first average value is: Among them, X i represents the effective wave height data of the observation field; h represents the number of effective wave height data of the observation field;

[0041] The observation field variance is calculated according to the first average value; the expression of the observation field variance is:

[0042] Calculating the background field variance according to the background field effective wave height data specifically includes:

[0043] Calculate the second average value of the background field effective wave height data; the expression of the second average value is: Among them, Y i represents the effective wave height data of the background field; g represents the number of effective wave height data of the background field;

[0044] The background field variance is calculated according to the second average value; the expression of the background field variance is:

[0045] Step S4: Calculate the observation operator according to the preset influence radius and the longitude and latitude distance, wherein the longitude and latitude distance is the distance between the longitude and latitude of the effective wave height grid point of the background field and the longitude and latitude of the altimeter observation position, or the distance between the longitude and latitude of the fused effective wave height field grid point and the longitude and latitude of the altimeter observation position;

[0046] The expression for calculating the observation operator according to the preset influence radius and longitude and latitude distance is:

[0047]

[0048] Wherein, W represents the observation operator; R represents the preset influence radius; and r represents the longitude and latitude distance.

[0049] Step S5: calculating the projected effective wave height according to the background field effective wave height data and the observation operator, wherein the projected effective wave height is the effective wave height of the background field projected to the altimeter observation position;

[0050] The expression for calculating the projected effective wave height according to the background field effective wave height data and the observation operator is:

[0051]

[0052] Among them, H b Represents the effective wave height data of the background field; ∑W represents the sum of all observation operators within the preset influence radius.

[0053] Step S6: Calculate the fused effective wave height according to the observation field effective wave height data, the background field effective wave height data, the observation field variance, the background field variance, the observation operator and the projected effective wave height.

[0054] The expression for calculating the fused effective wave height is:

[0055]

[0056] Among them, S o Represents the effective wave height data of the observation field.

[0057] The derivation process of the fused effective wave height expression is given below:

[0058] Let S ois the effective wave height of the altimeter in the observation field, E is the variance of the observation field, B is the variance of the background field, and the cost function J(H) is established.

[0059] J(H)=[S o -S(H)] T ·E -1 ·[S o -S(H)]+(H b -H) T ·B -1 ·(H b -H) T =min

[0060] Taking the derivative of J(H), we can get

[0061]

[0062] For S(H) in H b In the first-order Taylor expansion,

[0063] S(H)=S(H b )+S′(H b )·(HH b )

[0064] Substitute S(H) into Formula

[0065]

[0066] Further

[0067] S′ T (H) E -1 ·[S o -S(H b )-S′(H b )·(HH b )]+B -1 ·(H b -H)=0

[0068] Move Item

[0069] -[S′ T (H) E -1 ·S′(H b )+B -1 ]·(HH b )+S′ T (H) E -1 ·[S o -S(H b )]=0

[0070] It can be concluded that the effective wave height field to be determined is

[0071] H=H b +[S′ T (H) E -1 ·S′(H b )+B -1 ] -1 ·S′ T (H) E -1 ·[S o -S(H b )]

[0072] Approximately, S′ T (H)≈S′ T (H b ), that is,

[0073] H=H b +[S′ T (H b )·E -1 ·S′(H b )+B -1 ] -1 ·S′ T (H b )·E -1 ·[S o -S(H b )]

[0074] In order to determine the accuracy of the fused effective wave height, error analysis can be performed:

[0075] Obtain a variety of data released by ECWMF (European Centre for Medium-Range Weather Forecasts), including ECWMF numerical model forecast field data and ECWMF reanalysis data.

[0076] ECWMF numerical model forecast field data is a forecast field at a certain moment in the future obtained by numerical model calculation based on multi-source observation data and background field data. The present invention uses a 24-hour forecast field. The forecast field has a relatively low accuracy and can be used as a background field.

[0077] ECWMF reanalysis data is a complete set of reanalysis data obtained by quality controlling and assimilating observation data from various sources (ground, ships, radiosondes, wind balloons, aircraft, satellites, etc.) at a certain moment. It has high accuracy and can be used as a real field to verify the accuracy of fusion products, so ECWMF reanalysis data can be regarded as the actual effective wave height data observed at the observation station.

[0078] For example: during the verification, the 24-hour ECWMF numerical model forecast field of significant wave height at 00:00 on December 31, 2016 is used as the background field at 00:00 on January 1, 2017. The multi-source altimeter significant wave height data is integrated to obtain the required significant wave height field (fused significant wave height), which is compared with the ECWMF reanalysis data at 00:00 on January 1, 2017, and the root mean square error is calculated.

[0079] The significant wave height fusion product is matched with the ECWMF (European Centre for Medium-Range Weather Forecasts) reanalysis data in time and space, and error analysis is carried out. Therefore, after calculating the fused significant wave height, the following is also included:

[0080] Obtain the actual significant wave height data observed at the observation site;

[0081] Calculating a root mean square error based on the actual significant wave height data and the fused significant wave height data;

[0082] The calculation formula of the root mean square error is:

[0083]

[0084] in, Indicates the actual effective wave height data; (x j ,y j ) represents the latitude and longitude of the observation site; j = 1, 2, 3, ..., n, where n is the number of observation sites; represents the fused effective wave height data, which can also be called the fused effective wave height field position (x i ,y i ), i.e., the estimated value of the fused significant wave height;

[0085] The accuracy of the fused effective wave height data is judged according to the root mean square error.

[0086] The data fusion processing technology of the present invention is used to perform fusion processing on the global significant wave height data, and then the six sea areas are segmented.

[0087] The spatial range of the fused dataset (see 6 regions) Figure 2 ):

[0088] Northwest Pacific (0°~50°N, 100°E~180°);

[0089] Northeast Pacific (0°~40°N, 180°~120°W);

[0090] Southwest Pacific (0°~50°S, 120°E~160°W);

[0091] Southeast Pacific (0°~50°S, 180°~100°W);

[0092] Indian Ocean (50°S~30°N, 40°E~120°E).

[0093] The method used to prepare the dataset is to first fuse the effective wave height data, and then cut it according to six sea areas to obtain six sea area data products.

[0094] The present invention uses the ECWMF numerical model forecast field as the background field and the satellite observation significant wave height data as the observation field to perform global data fusion and then perform data clipping. The global significant wave height background field on January 1, 2017 is shown in Figure 2. Figure 3 , the effective wave height data of the altimeter on that day can be found at Figure 4 After the variational fusion of the present invention, the effective wave height data product is shown in Figure 5 The difference between the global significant wave height before and after fusion is +0.5m. Figure 6 , we can see the modulation effect of the effective wave height of the observation field. Figure 6 It can also be seen that the variational fusion method can make the sea area with a large difference between the effective wave height of the observation field and the effective wave height of the background field slowly transition, and can correct the effective wave height of the background field. In the sea area with a small difference between the effective wave height of the observation field and the effective wave height of the background field, there is basically no difference between the fused effective wave height and the effective wave height of the background field.

[0095] The spatial range of the fused dataset is divided into the China Coastal Sea (15°N~40°N, 105°E~130°E), Northwest Pacific (0°~50°N, 100°E~180°), Northeast Pacific (0°~40°N, 180°~120°W), Southwest Pacific (0°~50°S, 120°E~160°W), South Pacific (0°~50°S, 180°~100°W), and Indian Ocean (50°S~30°N, 40°E~120°E), so the error analysis of the fused effective wave height is carried out separately from the above six sea areas.

[0096] Substitute the ECWMF numerical model forecast field and ECWMF reanalysis field at the corresponding time into the root mean square error formula to calculate the root mean square error of the background field effective wave height. Substitute the fused effective wave height field and ECWMF reanalysis field at the corresponding time into the root mean square error formula to fuse the root mean square error of the effective wave height. The overall situation of the root mean square error of the effective wave height of the background field and the fused effective wave height in the six sea areas is obtained, as shown in Table 1.

[0097] Table 1 Overall situation of the root mean square error of the background field effective wave height and the fused effective wave height in the six sea areas

[0098]

[0099] In general, the root mean square error of the fused effective wave height is less than 0.5 m. Furthermore, the root mean square error of the fused effective wave height is lower than the root mean square error of the background effective wave height, indicating that the effective wave height of the observation field plays a positive role.

[0100] In this embodiment, the effective wave height data of the observation field and the effective wave height data of the background field are obtained, and the effective wave height analysis field is constructed by using the variational fusion method to fuse the effective wave height data of the observation field and the background field, thereby improving the accuracy of the fused effective wave height field.

[0101] Example 2

[0102] See also Figure 7 This embodiment provides an effective wave height fusion system based on multi-source satellite radar altimeter data, including:

[0103] The first data acquisition module L1 is used to obtain the effective wave height data of the observation field according to the data collected by the multi-source satellite radar altimeter;

[0104] The second data acquisition module L2 is used to obtain the background field effective wave height data according to the numerical model forecast field wave height data;

[0105] A variance calculation module L3, used to calculate the variance of the observation field according to the effective wave height data of the observation field, and to calculate the variance of the background field according to the effective wave height data of the background field;

[0106] An observation operator calculation module L4 is used to calculate the observation operator according to a preset influence radius and a longitude and latitude distance, wherein the longitude and latitude distance is the distance between the longitude and latitude of the effective wave height grid point of the background field and the longitude and latitude of the altimeter observation position;

[0107] The projected effective wave height calculation module L5 is used to calculate the projected effective wave height according to the background field effective wave height data and the observation operator, and the projected effective wave height is the effective wave height of the background field effective wave height projected to the altimeter observation position;

[0108] The fused effective wave height calculation module L6 is used to calculate the fused effective wave height according to the observed field effective wave height data, the background field effective wave height data, the observed field variance, the background field variance, the observation operator and the projected effective wave height.

[0109] In order to judge the accuracy of the fused effective wave height, error analysis can be performed, so the accuracy judgment module L7 of the fused effective wave height data can also be set, which specifically includes:

[0110] Obtain the actual significant wave height data observed at the observation site;

[0111] Calculating a root mean square error based on the actual significant wave height data and the fused significant wave height data;

[0112] The accuracy of the fused effective wave height data is judged according to the root mean square error.

[0113] As for the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the description of the method part.

[0114] The principles and implementation methods of the present invention are described in this article using specific examples. The description of the above embodiments is only used to help understand the method and core idea of ​​the present invention. At the same time, for those skilled in the art, according to the idea of ​​the present invention, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as limiting the present invention.

Claims

1. A method for effective wave height fusion based on multi-source satellite radar altimeter data, characterized in that: include: Obtain the effective wave height data of the observation field based on the data collected by the multi-source satellite radar altimeter; Obtain the background field effective wave height data based on the wave height data predicted by the numerical model; Calculate the observation field variance according to the observation field effective wave height data, and calculate the background field variance according to the background field effective wave height data; Calculate the observation operator according to the preset influence radius and the longitude and latitude distance, wherein the longitude and latitude distance is the distance between the longitude and latitude of the effective wave height grid point of the background field and the longitude and latitude of the altimeter observation position; Calculating a projected effective wave height according to the background field effective wave height data and the observation operator, wherein the projected effective wave height is the effective wave height of the background field projected to the altimeter observation position; Calculating a fused effective wave height according to the observation field effective wave height data, the background field effective wave height data, the observation field variance, the background field variance, the observation operator and the projected effective wave height; The step of calculating the observation field variance according to the observation field effective wave height data specifically includes: Calculate the first average value of the effective wave height data of the observation field; the expression of the first average value is: In the formula, X i represents the effective wave height data of the observation field; h represents the number of effective wave height data of the observation field; The observation field variance is calculated according to the first average value; the expression of the observation field variance is: The background field variance is calculated according to the background field effective wave height data, specifically including: Calculate the second average value of the background field effective wave height data; the expression of the second average value is: Where Y i represents the effective wave height data of the background field; g represents the number of effective wave height data of the background field; The background field variance is calculated according to the second average value; the expression of the background field variance is: The expression for calculating the observation operator according to the preset influence radius and longitude and latitude distance is: Where W represents the observation operator; R represents the preset influence radius; r represents the longitude and latitude distance; Among them, the expression for calculating the projected effective wave height according to the background field effective wave height data and the observation operator is: In the formula, H b Represents the effective wave height data of the background field; ΣW represents the sum of all observation operators within the preset influence radius; The expression for calculating the fused effective wave height is: In the formula, S o Represents the effective wave height data of the observation field.

2. The method according to claim 1, characterized in that After calculating the fused effective wave height, it also includes: Obtain the actual significant wave height data observed at the observation site; Calculating a root mean square error based on the actual significant wave height data and the fused significant wave height data; The accuracy of the fused effective wave height data is judged according to the root mean square error.

3. The method according to claim 2, characterized in that The calculation formula of the root mean square error is: in, Indicates the actual effective wave height data; (x j ,y j ) represents the latitude and longitude of the observation site; j = 1, 2, 3, ..., n, where n is the number of observation sites; Indicates fused significant wave height data.

4. An effective wave height fusion system based on multi-source satellite radar altimeter data, characterized in that: include: The first data acquisition module is used to acquire the effective wave height data of the observation field according to the data collected by the multi-source satellite radar altimeter; The second data acquisition module is used to obtain the background field effective wave height data according to the numerical model forecast field wave height data; A variance calculation module, used to calculate the variance of the observation field according to the effective wave height data of the observation field, and to calculate the variance of the background field according to the effective wave height data of the background field; The step of calculating the observation field variance according to the observation field effective wave height data specifically includes: Calculate the first average value of the effective wave height data of the observation field; the expression of the first average value is: In the formula, X i represents the effective wave height data of the observation field; h represents the number of effective wave height data of the observation field; The observation field variance is calculated according to the first average value; the expression of the observation field variance is: The background field variance is calculated according to the background field effective wave height data, specifically including: Calculate the second average value of the background field effective wave height data; the expression of the second average value is: Where Y i represents the effective wave height data of the background field; g represents the number of effective wave height data of the background field; The background field variance is calculated according to the second average value; the expression of the background field variance is: An observation operator calculation module, used to calculate the observation operator according to a preset influence radius and a longitude and latitude distance, wherein the longitude and latitude distance is the distance between the longitude and latitude of the effective wave height grid point of the background field and the longitude and latitude of the altimeter observation position; The expression for calculating the observation operator according to the preset influence radius and longitude and latitude distance is: Where W represents the observation operator; R represents the preset influence radius; r represents the longitude and latitude distance; A projected effective wave height calculation module, used to calculate the projected effective wave height according to the background field effective wave height data and the observation operator, wherein the projected effective wave height is the effective wave height of the background field projected to the altimeter observation position; Among them, the expression for calculating the projected effective wave height according to the background field effective wave height data and the observation operator is: In the formula, H b Represents the effective wave height data of the background field; ΣW represents the sum of all observation operators within the preset influence radius; A fused effective wave height calculation module, used for calculating a fused effective wave height according to the observed field effective wave height data, the background field effective wave height data, the observed field variance, the background field variance, the observation operator and the projected effective wave height; The expression for calculating the fused effective wave height is: In the formula, S o Represents the effective wave height data of the observation field.

5. The system according to claim 4, characterized in that It also includes an accuracy judgment module that integrates effective wave height data, including: Obtain the actual significant wave height data observed at the observation site; Calculating a root mean square error based on the actual significant wave height data and the fused significant wave height data; The accuracy of the fused effective wave height data is judged according to the root mean square error.

Images