A method for complementing missing meteorological monitoring data

Abstract

The present invention relates to the technical field of meteorological monitoring data, and discloses a method for completing missing meteorological monitoring data, including: selecting a monitoring station with missing data as the station to be completed, finding the nearest intermediate station with complete data for simulated missing processing, gradually increasing the number of modeling stations according to the distances between the intermediate station and surrounding monitoring stations, constructing an observation matrix and calculating the root mean square error, determining the optimal number of modeling stations, and using the data of the optimal modeling stations to complete the data of the station to be completed. It takes into account both the temporal sequence and the spatial distribution characteristics, has high interpretability and portability, and only relies on the physical characteristics driven by Beidou meteorological data itself without relying on external training data. The calculation is simple, there is no problem of differences in spatio-temporal resolution and update cycle between different data sources, and it can adaptively adjust to different monitoring station densities to ensure that the completed data approaches the true value, effectively improving the stability and integrity of Beidou meteorological monitoring.

Description

Technical Field

[0001] The present invention relates to the technical field of meteorological monitoring data, and in particular to a method for completing missing meteorological monitoring data. Background Art

[0002] The timeliness, accuracy and completeness of meteorological monitoring data are crucial for disaster warning, emergency response, sustainable development of agriculture and environment, and national security. BeiDou Satellite Navigation System, as an advanced technology independently developed by China, not only provides positioning, navigation and timing services, but also monitors ground meteorological data in real time, showing its great potential in the field of meteorological monitoring.

[0003] However, in actual applications, Beidou meteorological monitoring data is often missing due to various factors. Given the key role of meteorological data in climate research, disaster prediction, agricultural production and public safety, how to effectively supplement these missing data has become a core issue that needs to be solved urgently.

[0004] Traditional methods for missing meteorological data, such as interpolation and regression analysis, are often designed for general meteorological data and lack specialized solutions for Beidou meteorological monitoring data. At the same time, these methods usually only consider the temporal characteristics of the data, while ignoring the temporal and spatial distribution characteristics of Beidou meteorological monitoring stations. In recent years, although multi-source data fusion strategies have become a research hotspot, machine learning-based methods have achieved certain results, but their models have poor interpretability and portability, are highly dependent on a large amount of training data, have high computational complexity, and have differences in spatiotemporal resolution and update cycles between different data sources, which also limits the widespread application of this method in the meteorological field. Summary of the invention

[0005] Based on this, it is necessary to propose a method to complete the missing meteorological monitoring data in response to the above problems. It takes into account both temporal and spatial distribution characteristics, and has high interpretability and portability. At the same time, it only relies on the physical characteristics driven by the Beidou meteorological data itself, does not rely on external training data, is simple to calculate, and does not have the problem of differences in spatiotemporal resolution and update cycle between different data sources. It can also be adaptively adjusted to adapt to different monitoring station densities, ensuring that the completed data is close to the true value, effectively improving the stability and integrity of Beidou meteorological monitoring.

[0006] To achieve the above object, the present invention provides a method for completing missing meteorological monitoring data in a first aspect, the method comprising:

[0007] Obtain meteorological monitoring data of all monitoring stations in the target area within a preset time series, and select one of the monitoring stations with missing meteorological monitoring data as the station to be supplemented among all the monitoring stations;

[0008] Take the monitoring station that is closest to the whole station to be supplemented and has no missing meteorological monitoring data as the intermediate station, and use the meteorological monitoring data of the emptied intermediate station as the missing meteorological monitoring data of the intermediate station;

[0009] Take the n monitoring stations closest to the intermediate station as modeling stations, where n is a positive integer and the initial value of n is 1;

[0010] Form an observation matrix by arranging the missing meteorological monitoring data of the intermediate station and the meteorological monitoring data of the n modeling stations by column, and based on the observation matrix, complete the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station;

[0011] Determine the root mean square error based on the completed meteorological monitoring data of the intermediate station and the meteorological monitoring data;

[0012] Let n = n + 1, and return to execute the step of taking the n monitoring stations closest to the intermediate station as modeling stations until n is equal to the maximum modeling threshold, obtaining multiple root mean square errors, and take the n corresponding to the minimum root mean square error as the target modeling number, and set the target modeling stations as m;

[0013] Take the m monitoring stations closest to the whole station to be supplemented as modeling stations;

[0014] Form a to-be-completed observation matrix by arranging the meteorological monitoring data of the whole station to be supplemented and the meteorological monitoring data of the m modeling stations by column, and based on the to-be-completed observation matrix, complete the meteorological monitoring data of the whole station to be supplemented to obtain the completed meteorological monitoring data of the whole station to be supplemented.

[0015] Optionally, before the step of based on the observation matrix, completing the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station, the method further includes:

[0016] Determine whether the row vectors of the observation matrix at any observation epoch all satisfy non-empty according to the observation matrix;

[0017] In the case where there is a row vector at any observation epoch of the observation matrix that satisfies non-empty, let n = n + 1, and return to execute the step of taking the n monitoring stations closest to the intermediate station as modeling stations until the row vectors of the observation matrix at any observation epoch all do not satisfy non-empty.

[0018] Optionally, the step of based on the observation matrix, completing the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station includes:

[0019] Determine the covariance matrix according to the observation matrix;

[0020] Determine the non - missing row vector and the missing row vector at the i - th observation epoch according to the observation matrix, where i is a positive integer and the initial value of i is 1;

[0021] Determine the non - missing mean estimate value according to the non - missing row vector at the i - th observation epoch, and determine the missing mean estimate value according to the missing row vector at the i - th observation epoch;

[0022] Determine the non - missing estimation matrix corresponding to the covariance matrix of the non - missing row vector at the i - th observation epoch according to the covariance matrix, and determine the pre - cross - estimation matrix corresponding to the cross - covariance matrix between the non - missing row vector and the missing row vector at the i - th observation epoch;

[0023] Determine the regression coefficient according to the non - missing estimation matrix and the pre - cross - estimation matrix;

[0024] Determine the completed estimation row vector at the i - th observation epoch according to the regression coefficient, the non - missing mean estimate value, the missing mean estimate value, and the non - missing row vector at the i - th observation epoch;

[0025] Complete the row vector of the observation matrix at the i - th observation epoch according to the completed estimation row vector at the i - th observation epoch to obtain an intermediate matrix;

[0026] Extract the completed filling elements from the intermediate matrix to obtain the completed filling column vector and the completed filling column label;

[0027] Determine the stagnation tolerance according to the completed filling column vector and the completed filling column label;

[0028] Determine whether the intermediate matrix is completed according to the stagnation tolerance;

[0029] In the case that the intermediate matrix is not completed, determine the missing estimation matrix corresponding to the missing covariance matrix of the missing row vector at the i - th observation epoch according to the covariance matrix, and determine the post - cross - estimation matrix corresponding to the cross - covariance matrix between the missing row vector and the non - missing row vector at the i - th observation epoch;

[0030] Determine the residual covariance matrix according to the regression coefficient, the non - missing estimation matrix, the missing estimation matrix, the pre - cross - estimation matrix, and the post - cross - estimation matrix;

[0031] Determine the intermediate covariance matrix according to the residual covariance matrix and the covariance matrix;

[0032] Use the intermediate covariance matrix as the covariance matrix, let \(i = i + 1\), and return to execute the step of determining the non-missing row vector and the missing row vector at the \(i\)-th observation epoch according to the observation matrix until the intermediate matrix is completed;

[0033] Use the completed intermediate matrix as the target observation matrix;

[0034] According to the target observation matrix, complete the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station.

[0035] Optionally, the determining the regression coefficient according to the non-missing estimation matrix and the pre-cross estimation matrix includes:

[0036] Use the formula to determine the regression coefficient;

[0037] where, is the regression coefficient, is the diagonal matrix of the non-missing estimation matrix, is the non-missing estimation matrix, is the ridge parameter of ridge regression, is the identity matrix, is the pre-cross estimation matrix.

[0038] Optionally, the determining the completed estimation row vector at the \(i\)-th observation epoch according to the regression coefficient, the non-missing mean estimation value, the missing mean estimation value, and the non-missing row vector at the \(i\)-th observation epoch includes:

[0039] Use the formula to determine the completed estimation row vector at the \(i\)-th observation epoch;

[0040] where, is the completed estimation row vector at the \(i\)-th observation epoch, is the missing mean estimation value, is the non-missing row vector at the \(i\)-th observation epoch, is the non-missing mean estimation value, is the regression coefficient.

[0041] Optionally, the determining the residual covariance matrix according to the regression coefficient, the non-missing estimation matrix, the missing estimation matrix, the pre-cross estimation matrix, and the post-cross estimation matrix includes:

[0042] Use the formula to determine the residual covariance matrix;

[0043] where, is the residual covariance matrix, is the missing estimation matrix, is the regression coefficient, is the transpose symbol, is the non-missing estimation matrix, is the forward cross-estimation matrix, is the backward cross-estimation matrix.

[0044] Optionally, determining the covariance matrix according to the observation matrix includes:

[0045] Determining the column mean row vector according to the observation matrix;

[0046] Determining the centering matrix according to the column mean row vector and the observation matrix;

[0047] Determining the covariance matrix according to the centering matrix;

[0048] Determining the non-missing row vector and the missing row vector at the i-th observation epoch according to the observation matrix includes:

[0049] Determining the non-missing row vector and the missing row vector at the i-th observation epoch according to the centering matrix;

[0050] Completing the row vector of the observation matrix at the i-th observation epoch according to the completed estimated row vector at the i-th observation epoch to obtain an intermediate matrix, including:

[0051] Completing the row vector of the centering matrix at the i-th observation epoch according to the completed estimated row vector at the i-th observation epoch to obtain the intermediate matrix;

[0052] Determining the stagnation tolerance according to the completed filling column vector and the completed filling column label includes:

[0053] Determining the filled label mean row vector according to the completed filling column label and the column mean row vector;

[0054] Determining the stagnation tolerance according to the filled label mean row vector and the completed filling column vector;

[0055] Taking the completed intermediate matrix as the target observation matrix includes:

[0056] Determining an intermediate observation matrix according to the completed intermediate matrix and the column mean row vector;

[0057] Taking the intermediate observation matrix as the target observation matrix.

[0058] Optionally, determining the stagnation tolerance according to the filled label mean row vector and the completed filling column vector includes:

[0059] Using the formula to determine the stagnation tolerance;

[0060] wherein, is the matrix norm function, is the filling column vector, is the total number of missing elements of the observation matrix or the centralized matrix, is the filling mean row vector, is the transpose symbol, is the stagnation tolerance.

[0061] Optionally, determining whether the intermediate matrix is completed according to the stagnation tolerance includes:

[0062] When the stagnation tolerance is greater than or equal to the stagnation tolerance threshold, it is determined that the intermediate matrix is not completed;

[0063] When the stagnation tolerance is less than the stagnation tolerance threshold, it is determined that the intermediate matrix is completed.

[0064] Optionally, complementing the meteorological monitoring data of the station to be complemented according to the observation matrix to be complemented to obtain the complemented meteorological monitoring data of the station to be complemented includes:

[0065] Taking the observation matrix to be complemented as the observation matrix, returning to execute the step of determining the covariance matrix according to the observation matrix to obtain the target observation matrix, and taking the target observation matrix as the target complemented observation matrix;

[0066] Complementing the meteorological monitoring data of the station to be complemented according to the target complemented observation matrix to obtain the complemented meteorological monitoring data of the station to be complemented.

[0067] To achieve the above object, the present invention provides a device for complementing missing meteorological monitoring data in a second aspect, and the device includes:

[0068] An acquisition and selection module, configured to acquire the meteorological monitoring data of all monitoring stations in the target area within a preset time series, and select one of the monitoring stations with missing meteorological monitoring data as the station to be complemented among all monitoring stations;

[0069] A determination and nulling module, configured to use a monitoring station that is closest to the station to be complemented and has no missing meteorological monitoring data as the intermediate station, and use the meteorological monitoring data of the nulled intermediate station as the missing meteorological monitoring data of the intermediate station;

[0070] An iterative initial module for taking the n monitoring stations closest to the intermediate station as modeling stations, where n is a positive integer and the initial value of n is 1;

[0071] A first component completion module for forming an observation matrix by column with the missing meteorological monitoring data of the intermediate station and the meteorological monitoring data of the n modeling stations, and complementing the missing meteorological monitoring data of the intermediate station according to the observation matrix to obtain the complemented meteorological monitoring data of the intermediate station;

[0072] A root mean square error determination module for determining the root mean square error according to the complemented meteorological monitoring data of the intermediate station and the meteorological monitoring data;

[0073] An iterative module for setting n = n + 1, returning to execute the step of taking the n monitoring stations closest to the intermediate station as modeling stations until n is equal to the maximum modeling threshold, obtaining multiple root mean square errors, and taking the n corresponding to the minimum root mean square error as the target modeling number, and setting the target modeling stations as m;

[0074] A modeling determination module for taking the m monitoring stations closest to the station to be complemented as modeling stations;

[0075] A second component completion module for forming a to-be-complemented observation matrix by column with the meteorological monitoring data of the to-be-complemented station and the meteorological monitoring data of the m modeling stations, and complementing the meteorological monitoring data of the to-be-complemented station according to the to-be-complemented observation matrix to obtain the complemented meteorological monitoring data of the to-be-complemented station.

[0076] To achieve the above object, in a third aspect of the present invention, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by a processor, the processor is caused to execute the method according to any one of the first aspect.

[0077] To achieve the above object, in a fourth aspect of the present invention, a computer device is provided, including a memory and a processor, the memory storing a computer program, and when the computer program is executed by the processor, the processor is caused to execute the method according to any one of the first aspect.

[0078] Adopting the embodiment of the present invention has the following beneficial effects: By obtaining the meteorological monitoring data of all monitoring stations in the target area within a preset time series, and among all the monitoring stations, selecting one monitoring station with missing meteorological monitoring data as the station to be completely supplemented; taking the monitoring station that is closest to the station to be completely supplemented and has no missing meteorological monitoring data as the intermediate station, and taking the meteorological monitoring data of the emptied intermediate station as the missing meteorological monitoring data of the intermediate station; taking the n monitoring stations closest to the intermediate station as the modeling stations, where n is a positive integer and the initial value of n is 1; forming an observation matrix by arranging the missing meteorological monitoring data of the intermediate station and the meteorological monitoring data of the n modeling stations in columns, and based on the observation matrix, complementing the missing meteorological monitoring data of the intermediate station to obtain the complemented meteorological monitoring data of the intermediate station; determining the root mean square error according to the complemented meteorological monitoring data of the intermediate station and the meteorological monitoring data; letting n = n + 1, and returning to execute the step of taking the n monitoring stations closest to the intermediate station as the modeling stations until n is equal to the maximum modeling threshold, obtaining multiple root mean square errors, and taking the n corresponding to the minimum root mean square error as the target modeling number, and setting the target modeling stations as m; taking the m monitoring stations closest to the station to be completely supplemented as the modeling stations; forming a to-be-complemented observation matrix by arranging the meteorological monitoring data of the station to be completely supplemented and the meteorological monitoring data of the m modeling stations in columns, and based on the to-be-complemented observation matrix, complementing the meteorological monitoring data of the station to be completely supplemented to obtain the complemented meteorological monitoring data of the station to be completely supplemented. That is, by selecting the monitoring station with missing data as the station to be completely supplemented, finding the closest intermediate station with complete data for simulating the missing data, then gradually increasing the number of modeling stations according to the distance between the intermediate station and the surrounding monitoring stations, determining the optimal number of modeling stations by constructing an observation matrix and calculating the root mean square error, and finally using the data of the optimal modeling stations to complement the data of the station to be completely supplemented. This method not only considers the time series but also takes into account the spatial distribution characteristics, and has high interpretability and portability. At the same time, this method only depends on the physical characteristics driven by the Beidou meteorological data itself, does not rely on external training data, has simple calculations, does not have the problem of differences in spatio-temporal resolution and update cycle between different data sources, and can adaptively adjust to different monitoring station densities to ensure that the complemented data approaches the true value, effectively improving the stability and integrity of Beidou meteorological monitoring. BRIEF DESCRIPTION OF THE DRAWINGS

[0079] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the following drawings are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0080] Wherein:

[0081] Figure 1Schematic diagram of a method for completing missing meteorological monitoring data in an embodiment of the present application;

[0082] Figure 2 Schematic diagram of a device for completing missing meteorological monitoring data in an embodiment of the present application;

[0083] Figure 3 Internal structure diagram of a computer device in some embodiments. Detailed implementation manners

[0084] Next, the technical solutions in the embodiments of the present invention will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0085] The timeliness, accuracy, and integrity of meteorological monitoring data are crucial for aspects such as disaster warning, emergency response, sustainable development of agriculture and the environment, and national security. The Beidou satellite navigation system, as an advanced technology independently developed in China, not only provides positioning, navigation, and timing services but also can monitor ground meteorological data in real time, demonstrating its great potential in the field of meteorological monitoring.

[0086] However, in actual applications, Beidou meteorological monitoring data is often missing due to various factors. Given the key role of meteorological data in fields such as climate research, disaster prediction, agricultural production, and public safety, how to effectively complete these missing data has become a core problem to be solved urgently.

[0087] Traditional methods for completing missing meteorological data, such as interpolation and regression analysis, are often designed for general meteorological data and lack a dedicated solution for Beidou meteorological monitoring data. At the same time, these methods usually only consider the time characteristics of the data and ignore the temporal and spatial distribution characteristics simultaneously possessed by Beidou meteorological monitoring stations. In recent years, although the multi-source data fusion strategy has become a research hotspot, the methods based on machine learning have achieved certain results, but the interpretability and portability of their models are poor, and they highly rely on a large amount of training data, with high computational complexity. The differences in spatio-temporal resolution and update cycle between different data sources also limit the wide application of this method in the meteorological field.

[0088] In view of the above problems, the present application proposes a method for filling in missing meteorological monitoring data, which takes into account both temporal sequence and spatial distribution characteristics, has high interpretability and portability, and only relies on the physical characteristics driven by Beidou meteorological data itself without relying on external training data. It has simple calculations, does not have the problem of differences in spatio-temporal resolution and update cycle between different data sources, and can adaptively adjust to different monitoring station densities to ensure that the filled data approaches the true value, effectively improving the stability and integrity of Beidou meteorological monitoring. The specific implementation principle will be described in detail in the following embodiments.

[0089] The present application provides a method for filling in missing meteorological monitoring data in a first aspect.

[0090] Please refer to Figure 1 , which is a schematic diagram of a method for filling in missing meteorological monitoring data in an embodiment of the present application. The method includes:

[0091] Step 110: Obtain the meteorological monitoring data of all monitoring stations in the target area within a preset time series, and select one of the monitoring stations with missing meteorological monitoring data as the station to be filled in among all monitoring stations.

[0092] Among them, both the target area and the preset time series can be set by the operator based on a large amount of experience, experiments or statistics. Of course, they can also be set by the operator according to actual needs.

[0093] In some embodiments, the target area can be a research area determined by the operator; the preset time series can be a complete observation time series determined by the operator based on a large amount of experience, experiments or statistics.

[0094] Step 120: Select a monitoring station that is closest to the station to be filled in and has no missing meteorological monitoring data as the intermediate station, and use the meteorological monitoring data of the emptied intermediate station as the missing meteorological monitoring data of the intermediate station.

[0095] It should be noted that, in order to determine the optimal number of modeling stations for filling in data at the station to be filled in, the present application first selects a monitoring station that is closest to the station to be filled in and has no missing meteorological monitoring data as the intermediate station for modeling simulation, and then determines the optimal number of modeling.

[0096] It can be understood that since the intermediate station is close to the station to be filled in, it can be regarded as having similar meteorological conditions. Therefore, there is a strong spatio-temporal correlation between the meteorological monitoring data of the intermediate station and the station to be filled in.

[0097] Step 130: Select the n monitoring stations closest to the intermediate station as modeling stations, where n is a positive integer and the initial value of n is 1.

[0098] It should be noted that, in order to comprehensively consider the temporal and spatial characteristics of meteorological monitoring data, this application uses the n monitoring stations closest to the intermediate station as the modeling stations and performs iteration, aiming to simulate the most suitable number of modeling stations for the station to be supplemented.

[0099] It can be understood that there is a strong spatio-temporal correlation between the meteorological monitoring data between the intermediate station and the station to be supplemented, and there is also a strong spatio-temporal correlation between the meteorological monitoring data of the n monitoring stations and the intermediate station. Moreover, the terrain, meteorological data, monitoring station density, and missing degree in different regions are different, and the most suitable number of modeling stations is also different. Therefore, in order to determine the most suitable number of modeling stations for the station to be supplemented, this application uses spatio-temporal correlation to simulate the most suitable number of modeling stations for the station to be supplemented.

[0100] In some embodiments, the n monitoring stations closest to the intermediate station may or may not include the station to be supplemented. It can be understood that since the station to be supplemented is a monitoring station with missing meteorological monitoring data and subsequent iteration is required to determine the root mean square error and determine the optimal number of modeling stations through the root mean square error, whether the n monitoring stations closest to the intermediate station include the station to be supplemented will not affect the optimal number of modeling stations.

[0101] Step 140: Combine the missing meteorological monitoring data of the intermediate station and the meteorological monitoring data of the n modeling stations by column to form an observation matrix, and based on the observation matrix, complete the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station.

[0102] In some embodiments, the data of the intermediate station and the n modeling stations can be combined according to one column of data for each station to obtain an observation matrix. For example, when the data of each station in the preset time series is 100 data, the observation matrix is a matrix with 100 rows and n + 1 columns.

[0103] In some embodiments, a matrix completion algorithm (such as singular value decomposition or iterative optimization) can be used to complete the missing meteorological monitoring data of the intermediate station according to the observation matrix to obtain the completed meteorological monitoring data of the intermediate station.

[0104] Step 150: Determine the root mean square error based on the completed meteorological monitoring data of the intermediate station and the meteorological monitoring data.

[0105] In some embodiments, the formula can be used to determine the root mean square error; where is the root mean square error, is the total number of data of the completed meteorological monitoring data or the meteorological monitoring data of the intermediate station, is the r-th data in the completed meteorological monitoring data of the intermediate station, It is the r-th data in the meteorological monitoring data of the intermediate station.

[0106] Step 160: Let n = n + 1, and return to execute the step of taking the n monitoring stations closest to the intermediate station as modeling stations until n is equal to the maximum modeling threshold, obtaining multiple root mean square errors, and taking the n corresponding to the minimum root mean square error as the target modeling number. Let the target modeling stations be m.

[0107] Among them, the maximum modeling threshold can be set by the operator according to a large amount of experience, experiments or statistics. Of course, it can also be set by the operator according to actual needs.

[0108] It should be noted that the target modeling number is the number of models most suitable for the station to be complemented simulated by this application by utilizing spatio-temporal correlation.

[0109] Step 170: Take the m monitoring stations closest to the station to be complemented as modeling stations.

[0110] In some embodiments, the m monitoring stations closest to the station to be complemented may include the intermediate station.

[0111] Step 180: Compose the meteorological monitoring data of the station to be complemented and the meteorological monitoring data of the m modeling stations into a to-be-complemented observation matrix by columns, and according to the to-be-complemented observation matrix, complement the meteorological monitoring data of the station to be complemented to obtain the complemented meteorological monitoring data of the station to be complemented.

[0112] For the detailed description of this embodiment, reference can be made to the relevant content of Step 140 above, and it will not be elaborated here.

[0113] In the embodiments of this application, by selecting the monitoring station with missing data as the station to be complemented, finding the closest intermediate station with complete data for simulating the missing data, and then gradually increasing the number of modeling stations according to the distance between the intermediate station and the surrounding monitoring stations, determining the optimal number of modeling stations by constructing an observation matrix and calculating the root mean square error, and finally complementing the data of the station to be complemented with the data of the optimal modeling stations. This method not only considers the time series but also takes into account the spatial distribution characteristics, and has high interpretability and portability. At the same time, this method only depends on the physical characteristics driven by the Beidou meteorological data itself, does not rely on external training data, has simple calculations, does not have the problem of differences in spatio-temporal resolution and update period between different data sources, and can adaptively adjust to adapt to different monitoring station densities, ensuring that the complemented data approaches the true value, effectively improving the stability and integrity of Beidou meteorological monitoring.

[0114] In a feasible implementation manner, before step 140 in the above embodiment, that is, before completing the missing meteorological monitoring data of the intermediate station according to the observation matrix to obtain the completed meteorological monitoring data of the intermediate station, the method further includes: determining whether the row vectors of the observation matrix at any observation epoch satisfy non-emptiness according to the observation matrix; in the case that there exists a row vector at any observation epoch of the observation matrix that satisfies non-emptiness, let n = n + 1, and return to execute the step of taking the n monitoring stations closest to the intermediate station as the modeling stations until the row vectors of the observation matrix at any observation epoch do not satisfy non-emptiness.

[0115] For meteorological monitoring data, the observation epoch refers to the time points in the preset time series, and for the observation matrix, the observation epoch refers to the row label in the observation matrix.

[0116] It should be noted that in the case that there exists a row vector at any observation epoch of the observation matrix that satisfies non-emptiness, there will be a situation where it cannot be completely completed. At this time, the root mean square error determined by iteration must be relatively large. Therefore, the n in this iteration is definitely not the optimal number of modeling stations. In order to reduce the calculation amount and improve the calculation efficiency, the present application needs to perform an iterative judgment on whether it satisfies non-emptiness before the above step 140.

[0117] In the embodiment of the present application, by detecting and judging the non-emptiness of the observation matrix through pre-iteration, the modeling process is optimized, the calculation efficiency is improved, the feasibility and accuracy of the completion are ensured, and at the same time, the robustness of the method is enhanced, which has important practical significance and application value for the completion work of missing meteorological monitoring data.

[0118] In a feasible implementation manner, in step 140 of the above embodiments, the missing meteorological monitoring data of the intermediate station is complemented according to the observation matrix to obtain the complemented meteorological monitoring data of the intermediate station, including: determining the covariance matrix according to the observation matrix; determining the non-missing row vector and the missing row vector at the i-th observation epoch according to the observation matrix, where i is a positive integer and the initial value of i is 1; determining the non-missing mean estimate value according to the non-missing row vector at the i-th observation epoch, and determining the missing mean estimate value according to the missing row vector at the i-th observation epoch; determining the non-missing estimation matrix corresponding to the covariance matrix of the non-missing row vector at the i-th observation epoch according to the covariance matrix, and the pre-cross estimation matrix corresponding to the cross-covariance matrix between the non-missing row vector and the missing row vector at the i-th observation epoch; determining the regression coefficient according to the non-missing estimation matrix and the pre-cross estimation matrix; determining the complemented estimation row vector at the i-th observation epoch according to the regression coefficient, the non-missing mean estimate value, the missing mean estimate value, and the non-missing row vector at the i-th observation epoch; complementing the row vector of the observation matrix at the i-th observation epoch according to the complemented estimation row vector at the i-th observation epoch to obtain an intermediate matrix; extracting the completed filling elements from the intermediate matrix to obtain the completed filling column vector and the completed filling column label; determining the stagnation tolerance according to the completed filling column vector and the completed filling column label; determining whether the intermediate matrix is complemented according to the stagnation tolerance; in the case where the intermediate matrix is not complemented, determining the missing estimation matrix corresponding to the missing covariance matrix of the missing row vector at the i-th observation epoch according to the covariance matrix, and the post-cross estimation matrix corresponding to the cross-covariance matrix between the missing row vector and the non-missing row vector at the i-th observation epoch; determining the residual covariance matrix according to the regression coefficient, the non-missing estimation matrix, the missing estimation matrix, the pre-cross estimation matrix, and the post-cross estimation matrix; determining the intermediate covariance matrix according to the residual covariance matrix and the covariance matrix; taking the intermediate covariance matrix as the covariance matrix, setting i = i + 1, and returning to execute the step of determining the non-missing row vector and the missing row vector at the i-th observation epoch according to the observation matrix until the intermediate matrix is complemented; taking the complemented intermediate matrix as the target observation matrix; and complementing the missing meteorological monitoring data of the intermediate station according to the target observation matrix to obtain the complemented meteorological monitoring data of the intermediate station.

[0119] In determining the covariance matrix according to the observation matrix, in some embodiments, it is also necessary to set the missing elements in the observation matrix to zero; it can be understood that since the missing elements are missing values and the observation matrix cannot be calculated, it is necessary to set the missing elements in the observation matrix to zero.

[0120] For the determination methods of non-missing row vectors and missing row vectors at the i-th observation epoch, in some embodiments, the row vectors at the i-th observation epoch can be first determined according to the observation matrix, and then the missing and non-missing row vectors at the i-th observation epoch can be separated to obtain the non-missing row vectors and missing row vectors at the i-th observation epoch; for example, when the row vectors at the i-th observation epoch are [0, 4.01, 3.99, 0, 4.03], after separating the missing and non-missing row vectors, the non-missing row vectors at the i-th observation epoch obtained are [4.01, 3.99, 4.03], and the missing row vectors at the i-th observation epoch are [0, 0].

[0121] For the determination methods of the mean estimated value and the covariance matrix, existing methods can be used for calculation, which will not be elaborated here.

[0122] For the determination method of whether to complete the intermediate matrix, in some embodiments, it can be determined whether to complete the intermediate matrix according to the comparison result between the stagnation tolerance and the threshold.

[0123] For the determination method of the estimated matrix, in some embodiments, the matrix block division method can be used for calculation, that is, by performing block processing on the covariance matrix to obtain a non-missing estimated matrix, a pre-cross estimated matrix, a missing estimated matrix, and a post-cross estimated matrix.

[0124] For the determination methods of the regression coefficient and the residual covariance matrix, in some embodiments, the ridge regression algorithm can be used for calculation.

[0125] For the determination method of the intermediate covariance matrix, in some embodiments, the sum matrix between the residual covariance matrix and the covariance matrix is used as the intermediate covariance matrix.

[0126] For the determination method of the complemented meteorological monitoring data of the intermediate station, in some embodiments, the column elements corresponding to the intermediate station in the target observation matrix can be used as the complemented meteorological monitoring data of the intermediate station.

[0127] In the embodiments of the present application, by comprehensively considering strategies such as temporal and spatial characteristics, iterative optimization and dynamic adjustment, self-driven completion based on physical properties, high interpretability and portability, and computational efficiency and optimization, the efficient and accurate completion of missing Beidou meteorological monitoring data is realized, which is of great significance and application value for improving the stability and integrity of meteorological monitoring, supporting fields such as climate research, disaster prediction, agricultural production, and public safety.

[0128] It is understandable that by comprehensively considering temporal and spatial features: by constructing an observation matrix and combining the calculation of the covariance matrix, this method not only takes into account the temporality of meteorological monitoring data but also fully utilizes the spatial correlation between monitoring stations. This comprehensive consideration makes the completion result closer to the real meteorological conditions and improves the accuracy of completion; iterative optimization and dynamic adjustment: an iterative optimization strategy is adopted in this method. By gradually increasing the number of modeling stations and calculating the root mean square error, the optimal number of modeling stations is determined. This dynamic adjustment process can adaptively adapt to the density and missing degree of different monitoring stations, ensuring the flexibility and robustness of the completion method; self-driven completion based on physical characteristics: this method only relies on the physical characteristics driven by Beidou meteorological data itself and does not rely on external training data. This avoids the problems of differences in spatio-temporal resolution and update cycle that may exist between different data sources, making the completion result more reliable and consistent; high interpretability and portability: each step and algorithm design in this method are based on clear mathematical principles and statistical methods, with high interpretability. At the same time, since it does not rely on specific external data or models, this method is easy to be transplanted and applied in different regions and meteorological monitoring systems; computational efficiency and optimization: by detecting and judging the non-emptiness of the observation matrix through pre-iterative, the modeling process is optimized, reducing unnecessary computational amounts. At the same time, algorithms such as matrix block division and ridge regression are used for calculation, further improving the computational efficiency and accuracy.

[0129] In a feasible implementation manner, determining the regression coefficient according to the non-missing estimation matrix and the pre-cross estimation matrix in the above embodiment includes:

[0130] Using the formula to determine the regression coefficient;

[0131] where, is the regression coefficient, is the diagonal matrix of the non-missing estimation matrix, is the non-missing estimation matrix, is the ridge parameter of ridge regression, is the identity matrix, is the pre-cross estimation matrix.

[0132] It should be noted that is the ridge parameter in the ridge regression algorithm and can be preset by the operator.

[0133] In the embodiments of the present application, by using the formula to determine the regression coefficient, it has the advantages of improving the completion accuracy, enhancing stability, having clear mathematical principles and interpretability, strong adaptability, and high computational efficiency in the completion of missing meteorological monitoring data, making this method have broad application prospects and practical value in the field of meteorological monitoring.

[0134] It is understandable that improving the accuracy of completion: determining the regression coefficient through the ridge regression algorithm can more accurately reflect the relationship between non-missing data and missing data. Ridge regression performs well in dealing with multicollinearity and data instability, which helps to provide more robust estimates in missing data completion for meteorological monitoring; enhancing stability: ridge regression can control its complexity to a certain extent and prevent overfitting by introducing ridge parameters. In meteorological monitoring data, since the data may contain noise and outliers, this feature of ridge regression helps to enhance stability and make the completion results more reliable; clear mathematical principles and interpretability: using formulas to determine regression coefficients, based on clear mathematical principles and statistical methods, with high interpretability, which makes the design and implementation of the method more transparent, easy to understand and verify; strong adaptability: the ridge regression algorithm is applicable to various types of data sets, including meteorological monitoring data, which may have temporal and spatial correlation. By adjusting the ridge parameters, it can adapt to different data characteristics and missing situations and provide flexible completion solutions; high computational efficiency: formulated methods usually have high computational efficiency because they can directly obtain results through mathematical operations without complex iterations or search processes. This efficiency advantage is particularly important in the missing completion of meteorological monitoring data because of the need to process large amounts of data and real-time monitoring needs.

[0135] In a feasible implementation, the above embodiment determines the completed estimated row vector at the i-th observation epoch according to the regression coefficient, the non-missing mean estimate and the missing mean estimate, and the non-missing row vector at the i-th observation epoch, including:

[0136] Using the formula Determine the row vector of the completed estimate at the i-th observation epoch;

[0137] in, is the row vector of the completed estimate at the i-th observation epoch, is the missing mean estimate, is the non-missing row vector at the i-th observation epoch, is the non-missing mean estimate, is the regression coefficient.

[0138] In the embodiment of the present application, by using the formula to determine the completed estimated row vector at the i-th observation epoch, not only the accuracy and stability of the completion are improved, but also the interpretability and adaptability of the method are enhanced, while the computational efficiency is improved. These advantages make this method have broad application prospects and practical value in the field of missing completion of meteorological monitoring data.

[0139] In a feasible implementation manner, determining the residual covariance matrix according to the regression coefficient, the non-missing estimation matrix, the missing estimation matrix, the forward cross-estimation matrix, and the backward cross-estimation matrix in the above embodiments includes:

[0140] Using the formula to determine the residual covariance matrix;

[0141] where is the residual covariance matrix, is the missing estimation matrix, is the regression coefficient, is the transpose symbol, is the non-missing estimation matrix, is the forward cross-estimation matrix, is the backward cross-estimation matrix.

[0142] In the embodiments of the present application, determining the residual covariance matrix by using the formula plays an important role in the missing data completion of meteorological monitoring data. It not only improves the accuracy and stability of the completion, but also enhances the interpretability and adaptability of the method, providing strong support for the accurate and complete acquisition of meteorological monitoring data.

[0143] It can be understood that improving the completion accuracy: the calculation of the residual covariance matrix takes into account the complex relationship between the missing data and the non-missing data, and quantifies this relationship through an accurate mathematical formula, which helps to more accurately estimate the missing values during the iterative process, thereby improving the accuracy of the completed data; enhancing the robustness: in meteorological monitoring data, due to the influence of various factors (such as equipment failures, weather interference, etc.), the data may contain noise and outliers. By calculating the residual covariance matrix in a formulaic manner, the influence of these outliers on the completion result can be suppressed to a certain extent, enhancing the robustness of the method; supporting iterative optimization: the residual covariance matrix is a key input in the iterative optimization process, providing information about the difference between the current estimated value and the true value. By continuously updating the residual covariance matrix, the iterative process can be guided to converge to a better solution; clear mathematical principles and interpretability: using the formula to determine the residual covariance matrix is based on clear mathematical principles and statistical methods, with high interpretability, which makes the design and implementation of the method more transparent, facilitating understanding and verification.

[0144] In a feasible implementation manner, determining the covariance matrix according to the observation matrix in the above embodiments includes: determining the column mean row vector according to the observation matrix; determining the centering matrix according to the column mean row vector and the observation matrix; determining the covariance matrix according to the centering matrix.

[0145] Determining the non-missing row vector and the missing row vector at the i-th observation epoch according to the observation matrix in the above embodiments includes: determining the non-missing row vector and the missing row vector at the i-th observation epoch according to the centering matrix.

[0146] Completing the row vectors of the observation matrix at the i-th observation epoch according to the estimated completed row vector at the i-th observation epoch to obtain an intermediate matrix, including: completing the row vectors of the centering matrix at the i-th observation epoch according to the estimated completed row vector at the i-th observation epoch to obtain an intermediate matrix.

[0147] Determining the stagnation tolerance according to the completed filling column vector and the completed filling column label in the above embodiments includes: determining the mean row vector of filling labels according to the completed filling column label and the column mean row vector; determining the stagnation tolerance according to the mean row vector of filling labels and the completed filling column vector.

[0148] Regarding the completed intermediate matrix as the target observation matrix in the above embodiments includes: determining the intermediate observation matrix according to the completed intermediate matrix and the column mean row vector; regarding the intermediate observation matrix as the target observation matrix.

[0149] In determining the covariance matrix according to the observation matrix, in some embodiments, before determining the covariance matrix according to the centering matrix, it is also necessary to set the missing elements in the centering matrix to zero.

[0150] Regarding the determination method of the column mean row vector, in some embodiments, its missing elements can be ignored to determine the column mean row vector according to the valid value matrix of the observation matrix.

[0151] Regarding the determination method of the centering matrix, in some embodiments, the difference matrix between the observation matrix and the column mean row vector can be used as the centering matrix.

[0152] Regarding the determination method of the covariance matrix, in some embodiments, the formula can be used to determine the covariance matrix; where is the covariance matrix, is the centering matrix, is the transpose symbol, is the total number of rows of the centering matrix.

[0153] Regarding the determination method of the non-missing row vector and the missing row vector at the i-th observation epoch, in some embodiments, the row vector at the i-th observation epoch can be first determined according to the centering matrix, and then the missing and non-missing separation of the row vector at the i-th observation epoch can be performed to obtain the non-missing row vector and the missing row vector at the i-th observation epoch.

[0154] For the determination method of the filling label mean row vector, in some embodiments, according to the completed filling column labels, elements with the same column labels as the completed filling column labels can be extracted from the column mean row vector, and then these are combined into a row vector, and this row vector is used as the filling label mean row vector.

[0155] For the determination method of the intermediate observation matrix, in some embodiments, the sum matrix between the completed intermediate matrix and the column mean row vector can be used as the intermediate observation matrix.

[0156] In the embodiments of the present application, by introducing centering and refining relevant steps, significant effects have been achieved in the missing data completion of meteorological monitoring data, improving the accuracy, robustness, computational efficiency, interpretability, and portability of the completion. At the same time, temporal and spatial features are comprehensively considered, providing strong support for the accurate and complete acquisition of meteorological monitoring data.

[0157] It can be understood that the accuracy of completion is improved: by calculating the column mean row vector and centering the observation matrix, the mean shift in the data is eliminated, making subsequent calculations more accurate. The use of the centering matrix makes the calculation of the covariance matrix more precise, reflecting the true correlation between the data. During the completion process, the centering matrix is used to determine non-missing and missing row vectors, as well as calculate the completion estimation row vector, improving the accuracy of the completed data; the robustness is enhanced: the missing elements in the centering matrix are set to zero to avoid the interference of missing values on the calculation, enhancing the robustness of the method; through iterative optimization and dynamic adjustment strategies, combined with the judgment of the stagnation tolerance, the stability and convergence of the completion process are ensured; the computational efficiency is improved: by pre-calculating the column mean row vector and the centering matrix, subsequent calculation steps are simplified, improving the computational efficiency; using matrix operations and formula-based methods, such as the calculation of the covariance matrix and the determination of regression coefficients, further improves the calculation speed; the interpretability and portability are enhanced: each step and algorithm design in the method is based on clear mathematical principles and statistical methods, with high interpretability. Since it does not depend on specific external data or models, this method is easy to be transplanted and applied in different regions and meteorological monitoring systems; temporal and spatial features are comprehensively considered: by constructing the observation matrix and the centering matrix, and combining the calculation of the covariance matrix, this method not only considers the temporality of meteorological monitoring data but also fully utilizes the spatial correlation between monitoring stations.

[0158] In a feasible implementation manner, the determination of the stagnation tolerance according to the filling label mean row vector and the completed filling column vector in the above embodiments includes:

[0159] Using the formula to determine the stagnation tolerance;

[0160] where, is the matrix norm function, is the filled column vector, is the total number of missing elements in the observation matrix or the centering matrix, is the filled mean row vector, is the transpose symbol, is the stagnation tolerance.

[0161] In the embodiments of the present application, by using the formula to determine the stagnation tolerance, it has significant advantages in the missing data completion of meteorological monitoring data. It not only improves the accuracy and robustness of judgment, but also supports the iterative optimization process, improves the calculation efficiency, and enhances the interpretability and portability of the method. These advantages make the method have broad application prospects and practical value in the field of meteorological monitoring.

[0162] It can be understood that improving the judgment accuracy: by calculating the difference between the filled column vector and the filled mean row vector through the matrix norm function, the degree of change in the filling process can be quantified. This quantification method is more accurate than simple comparison or counting methods because it considers the overall differences of all elements; enhancing robustness: the calculation of the stagnation tolerance considers the completion of all missing elements, rather than relying solely on the changes of individual elements. This makes the method more robust to outliers or noise because large changes in individual elements will not overly affect the overall stagnation judgment; supporting iterative optimization: the stagnation tolerance is an important indicator in the iterative optimization process, used to judge whether the filling process has converged or reached a stable state. By calculating the stagnation tolerance in a formulaic way, the progress of the iterative process can be more objectively evaluated, thereby guiding the termination or continuation of the iterative process; improving the calculation efficiency: the calculation of the matrix norm function usually has high efficiency, especially with the support of modern computers and numerical calculation libraries. This makes the calculation of the stagnation tolerance not become the bottleneck of the entire filling process, thus ensuring the calculation efficiency of the overall method; enhancing interpretability and portability: using the formula to determine the stagnation tolerance, based on clear mathematical principles and statistical methods, has high interpretability. This method does not rely on specific external data or models and is easy to be transplanted and applied in different regions and meteorological monitoring systems.

[0163] In a feasible implementation manner, determining whether the intermediate matrix is completed according to the stagnation tolerance in the above embodiments includes: when the stagnation tolerance is greater than or equal to the stagnation tolerance threshold, determining that the intermediate matrix is not completed; when the stagnation tolerance is less than the stagnation tolerance threshold, determining that the intermediate matrix is completed.

[0164] Among them, the stagnation tolerance threshold can be set by the operator based on a large amount of experience, experiments or statistics. Of course, it can also be set by the operator according to actual needs.

[0165] In the embodiments of the present application, by using the stagnation tolerance and threshold to determine whether the intermediate matrix is completed, it has significant advantages in the missing data completion of meteorological monitoring data. It not only improves the objectivity, robustness and calculation efficiency of the determination, but also supports the iterative optimization process, enhances the interpretability and portability of the method. These advantages make the method have broad application prospects and practical value in the field of meteorological monitoring.

[0166] In a feasible implementation manner, step 180 in the above embodiment, according to the observation matrix to be completed, completes the meteorological monitoring data of the station to be completed, and obtains the completed meteorological monitoring data of the station to be completed, including: taking the observation matrix to be completed as the observation matrix, and returning to execute the step of determining the covariance matrix according to the observation matrix to obtain the target observation matrix, and taking the target observation matrix as the target completed observation matrix; according to the target completed observation matrix, completing the meteorological monitoring data of the station to be completed to obtain the completed meteorological monitoring data of the station to be completed.

[0167] For the determination method of the completed meteorological monitoring data of the station to be completed, in some embodiments, the column elements corresponding to the station to be completed in the target completed observation matrix can be used as the completed meteorological monitoring data of the station to be completed.

[0168] In the embodiments of the present application, by reprocessing the observation matrix to be completed based on the most appropriate number of models, and completing the meteorological monitoring data of the station to be completed according to the target completed observation matrix, this method shows significant advantages in the missing data completion of meteorological monitoring data, including improving the completion accuracy, enhancing the robustness, supporting dynamic adjustment, improving the calculation efficiency, and maintaining the consistency and portability of the method. These advantages make the method have broad application prospects and practical value in the field of meteorological monitoring.

[0169] The present application provides a device for completing missing meteorological monitoring data in a second aspect.

[0170] Please refer to Figure 2 , which is a schematic diagram of a device for completing missing meteorological monitoring data in the embodiments of the present application. The device 210 includes:

[0171] An acquisition and selection module 211, configured to acquire the meteorological monitoring data of all monitoring stations in the target area within a preset time series, and select one of the monitoring stations with missing meteorological monitoring data as the station to be completed among all monitoring stations;

[0172] A determination and nulling module 212, configured to use the monitoring station closest to the station to be completed and without missing meteorological monitoring data as the intermediate station, and use the meteorological monitoring data of the nulled intermediate station as the missing meteorological monitoring data of the intermediate station;

[0173] The iterative initial module 213 is used to take the n monitoring stations closest to the intermediate station as the modeling stations, where n is a positive integer and the initial value of n is 1;

[0174] The first component completion module 214 is used to form an observation matrix by column with the missing meteorological monitoring data of the intermediate station and the meteorological monitoring data of the n modeling stations, and based on the observation matrix, complete the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station;

[0175] The root mean square error determination module 215 is used to determine the root mean square error according to the completed meteorological monitoring data of the intermediate station and the meteorological monitoring data;

[0176] The iterative module 216 is used to set n = n + 1, return to execute the step of taking the n monitoring stations closest to the intermediate station as the modeling stations until n is equal to the maximum modeling threshold, obtain multiple root mean square errors, and take the n corresponding to the minimum root mean square error as the target modeling number, and set the target modeling stations as m;

[0177] The modeling determination module 217 is used to take the m monitoring stations closest to the station to be completed as the modeling stations;

[0178] The second component completion module 218 is used to form a to-be-completed observation matrix by column with the meteorological monitoring data of the station to be completed and the meteorological monitoring data of the m modeling stations, and based on the to-be-completed observation matrix, complete the meteorological monitoring data of the station to be completed to obtain the completed meteorological monitoring data of the station to be completed.

[0179] In the embodiments of the present application, the relevant contents of the above-mentioned acquisition and selection module 211, determination and blanking module 212, iterative initial module 213, first component completion module 214, root mean square error determination module 215, iterative module 216, modeling determination module 217 and second component completion module 218 can refer to Figure 1 the content in the shown embodiments, which will not be elaborated here.

[0180] It should be noted that the device 210 of the present application further includes some other modules. It can be understood that there is a one-to-one correspondence between the method of the present application and the device 210. Therefore, some other modules of the device 210 of the present application are the corresponding contents of the method of the present application in the above embodiments.

[0181] In the embodiment of the present application, by selecting the monitoring station with missing data as the station to be supplemented, finding the intermediate station with the closest distance and complete data for simulating the missing data, then gradually increasing the number of modeling stations according to the distances between the intermediate station and the surrounding monitoring stations, determining the optimal number of modeling stations by constructing an observation matrix and calculating the root mean square error, and finally using the data of the optimal modeling stations to complete the data of the station to be supplemented. This method not only considers the temporal sequence but also takes into account the spatial distribution characteristics, and has high interpretability and portability. At the same time, this method only depends on the physical characteristics driven by the Beidou meteorological data itself, does not rely on external training data, has simple calculations, does not have the problem of differences in spatio-temporal resolution and update cycle between different data sources, and can adaptively adjust to different monitoring station densities to ensure that the completed data approaches the true value, effectively improving the stability and integrity of Beidou meteorological monitoring.

[0182] In a third aspect, the present application also provides a computer-readable storage medium storing a computer program, which when executed by a processor causes the processor to execute a method for completing missing meteorological monitoring data in the above method embodiment.

[0183] In a fourth aspect, the present application also provides a computer device including a memory and a processor. The memory stores a computer program, which when executed by the processor causes the processor to execute a method for completing missing meteorological monitoring data in the above method embodiment.

[0184] Figure 3 The internal structure diagram of the computer device in some embodiments is shown. The computer device can specifically be a terminal, a server, or a gateway. As Figure 3 shown, the computer device includes a processor, a memory, and a network interface connected through a system bus.

[0185] Among them, the memory includes a non-volatile storage medium and an internal memory. The non-volatile storage medium of the computer device stores an operating system and can also store a computer program, which when executed by the processor can cause the processor to implement the various steps in the above method embodiment. The internal memory can also store a computer program, which when executed by the processor can cause the processor to execute the various steps in the above method embodiment. Those skilled in the art can understand that Figure 3 the structure shown in

[0186] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The program can be stored in a non-volatile computer-readable storage medium. When the program is executed, it can include the processes of the embodiments of the above methods.

[0187] Among them, any reference to a memory, storage, database, or other medium used in the embodiments provided in the present application may include non-volatile and / or volatile memories. Non-volatile memory may include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), or flash memory. Volatile memory may include random access memory (RAM) or an external cache memory. By way of illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDR SDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), Rambus direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and Rambus dynamic RAM (RDRAM), etc.

[0188] The technical features of the above embodiments can be combined arbitrarily. For the sake of brevity of description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered as the scope described in this specification.

[0189] The above-described embodiments merely represent several implementation manners of the present application. The description is relatively specific and detailed, but it should not be construed as a limitation on the patent scope of the present application. It should be noted that for those of ordinary skill in the art, without departing from the concept of the present application, several modifications and improvements can still be made, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the appended claims.

Claims

1. A method for complementing missing meteorological monitoring data, characterized in that, The method includes: Obtain the meteorological monitoring data of all monitoring stations in the target area within a preset time series, and among all the monitoring stations, select one monitoring station with missing meteorological monitoring data as the station to be fully supplemented; Take the monitoring station that is closest to the station to be fully supplemented and has no missing meteorological monitoring data as the intermediate station, and use the meteorological monitoring data of the emptied intermediate station as the missing meteorological monitoring data of the intermediate station; Take the n monitoring stations closest to the intermediate station as the modeling stations, where n is a positive integer and the initial value of n is 1; Form an observation matrix by arranging the missing meteorological monitoring data of the intermediate station and the meteorological monitoring data of the n modeling stations in columns, and based on the observation matrix, complete the missing meteorological monitoring data of the intermediate station to obtain the completed meteorological monitoring data of the intermediate station; Determine the root mean square error according to the completed meteorological monitoring data of the intermediate station and the meteorological monitoring data; Let n = n + 1, and return to execute the step of taking the n monitoring stations closest to the intermediate station as the modeling stations until n is equal to the maximum modeling threshold, obtaining multiple root mean square errors, and take the n corresponding to the minimum root mean square error as the target modeling number, and set the target modeling number as m; Take the m monitoring stations closest to the station to be fully supplemented as the modeling stations; Form a to-be-completed observation matrix by arranging the meteorological monitoring data of the station to be fully supplemented and the meteorological monitoring data of the m modeling stations in columns, and based on the to-be-completed observation matrix, complete the meteorological monitoring data of the station to be fully supplemented to obtain the completed meteorological monitoring data of the station to be fully supplemented.

2. The method for complementing missing meteorological monitoring data according to claim 1, wherein Before the step of completing the missing meteorological monitoring data of the intermediate station according to the observation matrix to obtain the completed meteorological monitoring data of the intermediate station, the method further includes: Determine whether the row vectors of the observation matrix at any observation epoch are all non-empty according to the observation matrix; In the case where there is a row vector at any observation epoch of the observation matrix that is non-empty, let n = n + 1, and return to execute the step of taking the n monitoring stations closest to the intermediate station as the modeling stations until the row vectors of the observation matrix at any observation epoch are all not non-empty.

3. The method for completing missing meteorological monitoring data according to claim 1, characterized in that, The step of completing the missing meteorological monitoring data of the intermediate station according to the observation matrix to obtain the completed meteorological monitoring data of the intermediate station includes: Determine the covariance matrix according to the observation matrix; Determine the non-missing row vector and the missing row vector at the i-th observation epoch according to the observation matrix, where i is a positive integer and the initial value of i is 1; Determine the non-missing mean estimate value according to the non-missing row vector at the i-th observation epoch, and determine the missing mean estimate value according to the missing row vector at the i-th observation epoch; Determine the non-missing estimation matrix corresponding to the covariance matrix of the non-missing row vector at the i-th observation epoch and the pre-cross estimation matrix corresponding to the cross-covariance matrix between the non-missing row vector and the missing row vector at the i-th observation epoch according to the covariance matrix; Determine the regression coefficient according to the non-missing estimation matrix and the pre-cross estimation matrix; Determine a completed estimated row vector at the $i$-th observation epoch according to the regression coefficient, the non-missing mean estimate value, the missing mean estimate value, and the non-missing row vector at the $i$-th observation epoch; Complete the row vector of the observation matrix at the $i$-th observation epoch according to the completed estimated row vector at the $i$-th observation epoch to obtain an intermediate matrix; Extract the completed filling elements from the intermediate matrix to obtain a completed filling column vector and completed filling column labels; Determine a stagnation tolerance according to the completed filling column vector and the completed filling column labels; Determine whether the intermediate matrix is completed according to the stagnation tolerance; In the case where the intermediate matrix is not completed, determine a missing estimation matrix corresponding to the missing covariance matrix of the missing row vector at the $i$-th observation epoch according to the covariance matrix, and a post-cross estimation matrix corresponding to the cross-covariance matrix between the missing row vector and the non-missing row vector at the $i$-th observation epoch; Determine a residual covariance matrix according to the regression coefficient, the non-missing estimation matrix, the missing estimation matrix, the pre-cross estimation matrix, and the post-cross estimation matrix; Determine an intermediate covariance matrix according to the residual covariance matrix and the covariance matrix; Take the intermediate covariance matrix as the covariance matrix, let $i = i + 1$, and return to execute the step of determining the non-missing row vector and the missing row vector at the $i$-th observation epoch according to the observation matrix until the intermediate matrix is completed; Take the completed intermediate matrix as the target observation matrix; Complete the missing meteorological monitoring data of the intermediate station according to the target observation matrix to obtain the completed meteorological monitoring data of the intermediate station.

4. The method for completing missing meteorological monitoring data according to claim 3, characterized in that The determining the regression coefficient according to the non-missing estimation matrix and the pre-cross estimation matrix includes: Using the formula to determine the regression coefficient; wherein, is the regression coefficient, is the diagonal matrix of the non-missing estimation matrix, is the non-missing estimation matrix, is the ridge parameter of ridge regression, is the identity matrix, is the forward cross-estimation matrix.

5. The method for completing missing meteorological monitoring data according to claim 3, characterized in that, The determining the completed estimated row vector at the $i$-th observation epoch according to the regression coefficient, the non-missing mean estimate value, the missing mean estimate value, and the non-missing row vector at the $i$-th observation epoch includes: Using the formula to determine the completed estimation row vector at the i-th observation epoch; wherein, is the completed estimation row vector at the i-th observation epoch, is the missing mean estimation value, is the non-missing row vector at the i-th observation epoch, is the non-missing mean estimation value, is the regression coefficient.

6. The method for complementing missing meteorological monitoring data according to claim 3, characterized in that, The determining the residual covariance matrix according to the regression coefficient, the non-missing estimation matrix, the missing estimation matrix, the pre-cross estimation matrix, and the post-cross estimation matrix includes: Using the formula to determine the residual covariance matrix; Among them, is the residual covariance matrix, is the missing estimation matrix, is the regression coefficient, is the transpose symbol, is the non-missing estimation matrix, is the front cross-estimation matrix, is the back cross-estimation matrix.

7. The method for complementing missing meteorological monitoring data according to claim 3, characterized in that The determining the covariance matrix according to the observation matrix includes: Determine a column mean row vector according to the observation matrix; Determine a centering matrix according to the column mean row vector and the observation matrix; Determine the covariance matrix according to the centering matrix; The determining the non-missing row vector and the missing row vector at the $i$-th observation epoch according to the observation matrix includes: Determine the non-missing row vector and the missing row vector at the $i$-th observation epoch according to the centering matrix; The completing the row vector of the observation matrix at the $i$-th observation epoch according to the completed estimated row vector at the $i$-th observation epoch to obtain an intermediate matrix includes: Complete the row vector of the centering matrix at the $i$-th observation epoch according to the completed estimated row vector at the $i$-th observation epoch to obtain the intermediate matrix; Determining the stagnation tolerance according to the completed filling column vector and the completed filling column label includes: Determining the mean row vector of filling labels according to the completed filling column label and the column mean row vector; Determining the stagnation tolerance according to the mean row vector of filling labels and the completed filling column vector; Regarding the completed intermediate matrix as the target observation matrix includes: Determining the intermediate observation matrix according to the completed intermediate matrix and the column mean row vector; Regarding the intermediate observation matrix as the target observation matrix.

8. The method for completing missing meteorological monitoring data according to claim 7, characterized in that Determining the stagnation tolerance according to the mean row vector of filling labels and the completed filling column vector includes: Use the formula to determine the stalling tolerance; Among them, is the matrix norm function, is the filling column vector, is the total number of missing elements of the observation matrix or the centering matrix, is the mean row vector of the filling labels, is the transpose symbol, is the stagnation tolerance.

9. The method for complementing missing meteorological monitoring data according to claim 3, characterized in that, Determining whether the intermediate matrix is completed according to the stagnation tolerance includes: When the stagnation tolerance is greater than or equal to the stagnation tolerance threshold, determining that the intermediate matrix is not completed; When the stagnation tolerance is less than the stagnation tolerance threshold, determining that the intermediate matrix is completed.

10. The method for complementing missing meteorological monitoring data according to any one of claims 3 to 9, characterized in that, Completing the meteorological monitoring data of the station to be completed according to the observation matrix to be completed to obtain the completed meteorological monitoring data of the station to be completed includes: Regarding the observation matrix to be completed as the observation matrix, returning to execute the step of determining the covariance matrix according to the observation matrix to obtain the target observation matrix, and regarding the target observation matrix as the target completed observation matrix; Completing the meteorological monitoring data of the station to be completed according to the target completed observation matrix to obtain the completed meteorological monitoring data of the station to be completed.

Images