CN113075752A - Method and device for judging correctness of three-dimensional space position of meteorological observation station - Google Patents

Method and device for judging correctness of three-dimensional space position of meteorological observation station Download PDF

Info

Publication number
CN113075752A
CN113075752A CN202110413289.XA CN202110413289A CN113075752A CN 113075752 A CN113075752 A CN 113075752A CN 202110413289 A CN202110413289 A CN 202110413289A CN 113075752 A CN113075752 A CN 113075752A
Authority
CN
China
Prior art keywords
data
statistical analysis
air pressure
observation station
analysis index
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110413289.XA
Other languages
Chinese (zh)
Other versions
CN113075752B (en
Inventor
崔喜爱
赵培涛
施丽娟
李翠娜
林雪娇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CMA Meteorological Observation Centre
Original Assignee
CMA Meteorological Observation Centre
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CMA Meteorological Observation Centre filed Critical CMA Meteorological Observation Centre
Priority to CN202110413289.XA priority Critical patent/CN113075752B/en
Publication of CN113075752A publication Critical patent/CN113075752A/en
Application granted granted Critical
Publication of CN113075752B publication Critical patent/CN113075752B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/10Devices for predicting weather conditions
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W1/18Testing or calibrating meteorological apparatus
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Environmental & Geological Engineering (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Atmospheric Sciences (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention relates to a method and a device for judging the correctness of a three-dimensional space position of a meteorological observation station, wherein the method comprises the following steps: acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method; carrying out statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, and further determining the correctness of longitude data and latitude data of the meteorological observation station; for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actual measurement altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode; and judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.

Description

Method and device for judging correctness of three-dimensional space position of meteorological observation station
Technical Field
The disclosure relates to the technical field of meteorological observation, in particular to a method and a device for judging correctness of a three-dimensional space position of a meteorological observation station.
Background
The weather detection information is an important basic resource for national economic and social construction and weather service and scientific research. The longitude, the latitude and the altitude of the meteorological observation station are the most important meteorological detection metadata information of the meteorological observation station, the longitude and the latitude identify the specific geographic position of the observation station, and the altitude identifies the vertical spatial position of the observation station, so that the longitude, the latitude and the altitude are the most important parameters for developing and analyzing meteorological observation data of the meteorological observation station.
The judgment of the correctness of the longitude, the latitude and the altitude of the meteorological observation station is also a key link for ensuring the high quality of observation data. In the actual service operation process, the weather observation station always generates the situation of station address migration, the longitude, the latitude and the altitude before and after the station address migration all change, and the weather elements detected by the corresponding weather sensors also change in the whole trend. Sometimes, because an error occurs in a certain node such as manual measurement, data entry, database update, etc., an erroneous value may occur in the longitude, latitude, or altitude of the weather observation station, which is a great hidden danger for quality control, product processing, and subsequent service application of subsequent observation data.
Disclosure of Invention
In order to overcome the problems in the related art, the present disclosure provides a method and an apparatus for judging the correctness of a three-dimensional spatial position of a weather-observing station, so as to judge whether the longitude, the latitude and the altitude of the weather-observing station are correct in time, and correct the incorrect longitude, the latitude and the altitude of the weather-observing station in time.
According to a first aspect of the embodiments of the present disclosure, there is provided a method for determining correctness of a three-dimensional space position of a weather observation station, the method including:
acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
determining correctness of longitude data and latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actually-measured altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
and judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
For the meteorological observation station with incorrect longitude data, latitude data or altitude data, retest reminders can be output to the meteorological observation station, and the approval process of retest is started, so that retest can be ensured to be carried out.
In one embodiment, preferably, the performing a statistical analysis index calculation according to the air pressure data of the weather observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index includes:
according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, the first statistical analysis index, the second statistical analysis index and the third statistical analysis index are calculated by adopting the following first formula:
Figure BDA0003024797490000031
wherein Index1(PRE) represents the first statistical analysis indicator, PRE (i) represents the air pressure data of the ith weather observation station, SimulatedPRE (i) represents the simulated air pressure data of the ith corresponding air numerical prediction mode, NUMPER represents the number of air pressure data of the weather observation stations, Index2(PRE) represents the second statistical analysis indicator, Index3(PRE) represents the third statistical analysis indicator, wherein the first statistical analysis index is used for representing the coincidence degree between the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the first statistical analysis index is used for representing the standard deviation of a difference value sequence of an air pressure data sequence of a meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period; and the third statistical analysis index is used for representing the root mean square value of a difference value sequence of the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in a preset time period.
In one embodiment, preferably, the determining the correctness of the longitude data and the latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index includes:
when the first statistical analysis index, the second statistical analysis index and the third statistical analysis index meet the following conditions, determining that longitude data and latitude data of the meteorological observation station are abnormal data, and otherwise, determining that the longitude data and the latitude data are normal data;
i Index1(PRE) | > 4 or
Index2(PRE) > 6 or
Index3(PRE)>6
In one embodiment, the predicted altitude data is preferably calculated using the following second formula,
Figure BDA0003024797490000041
wherein H _ B represents simulated altitude data of an atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES _ B represents simulated barometric pressure data of the atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES represents barometric pressure data of the target weather-observing station, TEM1 represents air temperature data of the target weather-observing station, Alti represents measured altitude data of the target weather-observing station, and log _10 is a logarithmic function with a base 10.
In one embodiment, preferably, the determining the correctness of the measured altitude data of the target weather observation station according to the difference between the predicted altitude data and the measured altitude data includes:
and when the difference value between the predicted altitude data and the actually measured altitude data is larger than a preset threshold value, judging that the actually measured altitude data of the target meteorological observation station is incorrect.
According to a second aspect of the embodiments of the present disclosure, there is provided a three-dimensional spatial position correctness determination apparatus for a weather observation station, the apparatus including:
the system comprises a preprocessing module, a data acquisition module and a data processing module, wherein the preprocessing module is used for acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
the first calculation module is used for performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
the determining module is used for determining the correctness of the longitude data and the latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
the second calculation module is used for calculating predicted altitude data of a target meteorological observation station with correct longitude data and latitude data according to the air pressure data, the air temperature data, the actually-measured altitude data of the target meteorological observation station, and the simulated altitude data and the simulated air pressure data of the corresponding atmospheric numerical prediction mode;
and the judging module is used for judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
In one embodiment, preferably, the first calculation module is configured to:
according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, the first statistical analysis index, the second statistical analysis index and the third statistical analysis index are calculated by adopting the following first formula:
Figure BDA0003024797490000051
wherein Index1(PRE) represents the first statistical analysis indicator, PRE (i) represents the air pressure data of the ith weather observation station, SimulatedPRE (i) represents the simulated air pressure data of the ith corresponding air numerical prediction mode, NUMPER represents the number of air pressure data of the weather observation stations, Index2(PRE) represents the second statistical analysis indicator, Index3(PRE) represents the third statistical analysis indicator, wherein the first statistical analysis index is used for representing the coincidence degree between the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the first statistical analysis index is used for representing the standard deviation of a difference value sequence of an air pressure data sequence of a meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period; the third statistical analysis index is used for representing the root mean square value of a difference value sequence of an air pressure data sequence of the meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period;
the determination module is to:
when the first statistical analysis index, the second statistical analysis index and the third statistical analysis index meet the following conditions, determining that longitude data and latitude data of the meteorological observation station are abnormal data, and otherwise, determining that the longitude data and the latitude data are normal data;
i Index1(PRE) | > 4 or
Index2(PRE) > 6 or
Index3(PRE)>6
In one embodiment, the predicted altitude data is preferably calculated using the following second formula,
Figure BDA0003024797490000061
wherein H _ B represents simulated altitude data of an atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES _ B represents simulated air pressure data of the atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES represents air pressure data of the target weather-observing station, TEM1 represents air temperature data of the target weather-observing station, Alti represents measured altitude data of the target weather-observing station, and log _10 is a logarithmic function with a base 10;
the judging module is used for:
and when the difference value between the predicted altitude data and the actually measured altitude data is larger than a preset threshold value, judging that the actually measured altitude data of the target meteorological observation station is incorrect.
According to a third aspect of the embodiments of the present disclosure, there is provided a three-dimensional spatial position correctness determination apparatus for a weather observation station, the apparatus including:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to:
acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
determining correctness of longitude data and latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actually-measured altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
and judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
According to a fourth aspect of embodiments of the present disclosure, there is provided a computer-readable storage medium having stored thereon computer instructions which, when executed by a processor, implement the steps of the method of any one of the first aspects.
The technical scheme provided by the embodiment of the disclosure can have the following beneficial effects:
in the embodiment of the invention, the air pressure data and the air temperature data of the weather observation station which has been operated for more than the preset time and the corresponding simulated air pressure data and simulated altitude data of the atmospheric numerical prediction mode are obtained, so that whether the longitude, the latitude and the altitude of the weather observation station are correct or not can be automatically judged, and the incorrect longitude, the latitude and the altitude of the weather observation station can be corrected in time.
It is to be understood that both the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure.
Drawings
The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the present disclosure and together with the description, serve to explain the principles of the disclosure.
FIG. 1 is a flow chart illustrating a method for determining the correctness of the three-dimensional spatial position of a weather observation station according to an exemplary embodiment.
Fig. 2a and 2b are schematic diagrams showing incorrect longitude and latitude data of a certain weather observation station 1 according to an exemplary embodiment.
Fig. 3a and 3b are schematic diagrams illustrating incorrect longitude and latitude data for a certain weather-observing station 2 according to an exemplary embodiment.
Fig. 4a and 4b are schematic diagrams showing incorrect data of altitude for a certain weather-observing station 3 according to an exemplary embodiment.
Fig. 5a and 5b are schematic diagrams showing incorrect data for altitude for a certain weather-observing station 3, according to an exemplary embodiment.
FIG. 6 is a block diagram illustrating a three-dimensional spatial position correctness determination apparatus for a weather observation station according to an exemplary embodiment.
Detailed Description
Reference will now be made in detail to the exemplary embodiments, examples of which are illustrated in the accompanying drawings. When the following description refers to the accompanying drawings, like numbers in different drawings represent the same or similar elements unless otherwise indicated. The implementations described in the exemplary embodiments below are not intended to represent all implementations consistent with the present disclosure. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present disclosure, as detailed in the appended claims.
FIG. 1 is a flowchart illustrating a method for determining correctness of three-dimensional spatial positions of a weather observation station, according to an exemplary embodiment, as shown in FIG. 1, the method includes:
step S101, acquiring air pressure data and air temperature data of a weather observation station which runs for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
the method comprises the following steps of preprocessing by using a threshold limiting method to remove obviously abnormal data, and retaining the data according to the following principle: the air pressure data value should be within the range of 300 to 1100 hPa; the temperature data value should be in the range of-55 deg.C to 55 deg.C.
Step S102, carrying out statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
in one embodiment, preferably, the performing a statistical analysis index calculation according to the air pressure data of the weather observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index includes:
according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, the first statistical analysis index, the second statistical analysis index and the third statistical analysis index are calculated by adopting the following first formula:
Figure BDA0003024797490000091
wherein Index1(PRE) represents the first statistical analysis indicator, PRE (i) represents the air pressure data of the ith weather observation station, SimulatedPRE (i) represents the simulated air pressure data of the ith corresponding air numerical prediction mode, NUMPER represents the number of air pressure data of the weather observation stations, Index2(PRE) represents the second statistical analysis indicator, Index3(PRE) represents the third statistical analysis indicator, wherein the first statistical analysis index is used for representing the coincidence degree between the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the numerical value of the first statistical analysis index is inversely related to the coincidence degree, that is, a smaller value indicates a better matching degree, and if the longitude and latitude data are abnormal, the matching degree is poor, that is, the value of Index1 is high. The second statistical analysis Index is used for representing the standard deviation of the difference value sequence of the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the smaller the value is, the smaller the deviation of the difference value sequence from the average value is, if the longitude and latitude data are abnormal, the larger the deviation of the difference value sequence from the average value is, namely, the value of Index2 is high. The third statistical analysis Index is used for representing the root mean square value of a difference value sequence of the air pressure data sequence of the weather observation station in the preset time period and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode, and is also an Index reflecting the dispersion degree of the difference value sequence, the smaller the numerical value is, the smaller the dispersion degree of the difference value sequence is, if the longitude and latitude data are abnormal, the larger the root mean square value of the difference value sequence is, namely, the numerical value of Index3 is very high.
Step S103, determining the correctness of longitude data and latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
in one embodiment, preferably, the determining the correctness of the longitude data and the latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index includes:
when the first statistical analysis index, the second statistical analysis index and the third statistical analysis index meet the following conditions, determining that longitude data and latitude data of the meteorological observation station are abnormal data, and otherwise, determining that the longitude data and the latitude data are normal data;
i Index1(PRE) | > 4 or
Index2(PRE) > 6 or
Index3(PRE)>6
Step S104, for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actual measurement altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
in one embodiment, the predicted altitude data is preferably calculated using the following second formula,
Figure BDA0003024797490000101
wherein H _ B represents simulated altitude data of an atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES _ B represents simulated barometric pressure data of the atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES represents barometric pressure data of the target weather-observing station, TEM1 represents air temperature data of the target weather-observing station, Alti represents measured altitude data of the target weather-observing station, and log _10 is a logarithmic function with a base 10.
And step S105, judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
In one embodiment, preferably, the determining the correctness of the measured altitude data of the target weather observation station according to the difference between the predicted altitude data and the measured altitude data includes:
and when the difference value between the predicted altitude data and the actually measured altitude data is larger than a preset threshold value, judging that the actually measured altitude data of the target meteorological observation station is incorrect.
For example, if the preset threshold is 20m, the difference between the predicted altitude data and the measured altitude data is greater than 20m, that is, it is determined that the measured altitude data of the target weather-observing station is incorrect.
For the meteorological observation station with incorrect longitude, latitude and altitude data, the surrounding detection environment, the working state of the sensor and other multi-aspect information (such as experience data) of the station can be combined, a conclusion can be given after comprehensive judgment, and the longitude, the latitude and the altitude of the meteorological observation station with incorrect data can be corrected in time.
In the above embodiment, by acquiring the air pressure data and the air temperature data of the weather observation station that has been operated for more than the preset time, and the simulated air pressure data and the simulated altitude data of the atmospheric numerical prediction mode corresponding thereto, it is possible to automatically determine whether the longitude, the latitude, and the altitude of the weather observation station are correct, thereby correcting the longitude, the latitude, and the altitude of the weather observation station that are incorrect in time.
The above technical solution of the present invention is described in detail in a specific case.
Through the data processing and analyzing experiments of the correctness judgment test of the longitude, the latitude and the altitude of the national-level automatic station for more than one year, the invention is found that the metadata information of the observation station with wrong longitude, latitude or altitude can be effectively judged. The following are several cases:
(1) case 1 for determining correctness of latitude and longitude
The weather observation station longitude, latitude and altitude correctness judgment test analysis shows that the longitude and latitude information of a certain observation station 1 is abnormal as shown in fig. 2a and fig. 2b, and the longitude and latitude correctness judgment of the observation station is not passed. Fig. 2b shows the data of the station, wherein the longitude and latitude of the scribed lines are 46.4628 and 133.83, respectively, which are the error data determined by the method, and the normal longitude and latitude should be 41.8983 and 119.7131, respectively.
(2) Case 2 for determining correctness of latitude and longitude
The correctness of the longitude, the latitude and the altitude of the country-level automatic station is judged, the correctness of the longitude and the latitude of an observation station 2 is not judged, and the longitude and the latitude information of the station are displayed to be abnormal, which is shown in figure 3. The right side of fig. 3 shows the data of the station, wherein the longitude and latitude of the drawing are 33.4333 and 101.4833, respectively, which are the error data determined by the method, and the normal longitude and latitude should be 32.0667 and 121.6, respectively.
(3) Case 1 of altitude correctness determination
The longitude, latitude and altitude correctness judgment test analysis of the country-level automatic station finds that a certain observation station 3 does not pass through an altitude correctness judgment module of the observation station, and displays that the altitude information of the station is abnormal, as shown in fig. 4. The station altitude 1041.5m is shown on the left side of FIG. 4, which is the erroneous data identified by the method, and after repeated measurements, the altitude should be 1003.2 m.
(4) Case 2 for judging altitude correctness
The longitude, latitude and altitude correctness judgment test analysis of the country-level automatic station finds that a certain observation station 4 does not pass through an altitude correctness judgment module of the observation station, and displays that the altitude information of the station is abnormal, as shown in fig. 5. The left side of fig. 5 shows that the station has an altitude of 247.4m, which is the erroneous data identified by the method, and after repeated measurements, the altitude is found to be 219.7 m.
FIG. 6 is a block diagram illustrating a three-dimensional spatial position correctness determination apparatus for a weather observation station according to an exemplary embodiment.
As shown in fig. 6, according to a second aspect of the embodiments of the present disclosure, there is provided a three-dimensional spatial position correctness determination apparatus for a weather observation station, the apparatus including:
the preprocessing module 61 is used for acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
the first calculation module 62 is configured to perform statistical analysis index calculation according to the air pressure data of the weather observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index, and a third statistical analysis index;
a determining module 63, configured to determine correctness of longitude data and latitude data of the weather observation station according to the first statistical analysis index, the second statistical analysis index, and the third statistical analysis index;
a second calculation module 64, configured to calculate, for a target weather observation station with correct longitude data and latitude data, predicted altitude data according to air pressure data, air temperature data, actual measurement altitude data of the target weather observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
and the judging module 65 is configured to judge the correctness of the actual measurement altitude data of the target weather observation station according to the difference between the predicted altitude data and the actual measurement altitude data.
In one embodiment, the first calculation module 62 is preferably configured to:
according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, the first statistical analysis index, the second statistical analysis index and the third statistical analysis index are calculated by adopting the following first formula:
Figure BDA0003024797490000131
wherein Index1(PRE) represents the first statistical analysis indicator, PRE (i) represents the air pressure data of the ith weather observation station, SimulatedPRE (i) represents the simulated air pressure data of the ith corresponding air numerical prediction mode, NUMPER represents the number of air pressure data of the weather observation stations, Index2(PRE) represents the second statistical analysis indicator, Index3(PRE) represents the third statistical analysis indicator, wherein the first statistical analysis index is used for representing the coincidence degree between the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the first statistical analysis index is used for representing the standard deviation of a difference value sequence of an air pressure data sequence of a meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period; the third statistical analysis index is used for representing the root mean square value of a difference value sequence of an air pressure data sequence of the meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period;
the determining module 63 is configured to:
when the first statistical analysis index, the second statistical analysis index and the third statistical analysis index meet the following conditions, determining that longitude data and latitude data of the meteorological observation station are abnormal data, and otherwise, determining that the longitude data and the latitude data are normal data;
i Index1(PRE) | > 4 or
Index2(PRE) > 6 or
Index3(PRE)>6
In one embodiment, the predicted altitude data is preferably calculated using the following second formula,
Figure BDA0003024797490000141
wherein H _ B represents simulated altitude data of an atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES _ B represents simulated air pressure data of the atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES represents air pressure data of the target weather-observing station, TEM1 represents air temperature data of the target weather-observing station, Alti represents measured altitude data of the target weather-observing station, and log _10 is a logarithmic function with a base 10;
the judging module 65 is configured to:
and when the difference value between the predicted altitude data and the actually measured altitude data is larger than a preset threshold value, judging that the actually measured altitude data of the target meteorological observation station is incorrect.
According to a third aspect of the embodiments of the present disclosure, there is provided a three-dimensional spatial position correctness determination apparatus for a weather observation station, the apparatus including:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to:
acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
determining correctness of longitude data and latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actually-measured altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
and judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
According to a fourth aspect of embodiments of the present disclosure, there is provided a computer-readable storage medium having stored thereon computer instructions which, when executed by a processor, implement the steps of the method of any one of the first aspects.
It is further understood that the use of "a plurality" in this disclosure means two or more, as other terms are analogous. "and/or" describes the association relationship of the associated objects, meaning that there may be three relationships, e.g., a and/or B, which may mean: a exists alone, A and B exist simultaneously, and B exists alone. The character "/" generally indicates that the former and latter associated objects are in an "or" relationship. The singular forms "a", "an", and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It will be further understood that the terms "first," "second," and the like are used to describe various information and that such information should not be limited by these terms. These terms are only used to distinguish one type of information from another and do not denote a particular order or importance. Indeed, the terms "first," "second," and the like are fully interchangeable. For example, first information may also be referred to as second information, and similarly, second information may also be referred to as first information, without departing from the scope of the present disclosure.
It is further to be understood that while operations are depicted in the drawings in a particular order, this is not to be understood as requiring that such operations be performed in the particular order shown or in serial order, or that all illustrated operations be performed, to achieve desirable results. In certain environments, multitasking and parallel processing may be advantageous.
Other embodiments of the disclosure will be apparent to those skilled in the art from consideration of the specification and practice of the disclosure disclosed herein. This application is intended to cover any variations, uses, or adaptations of the disclosure following, in general, the principles of the disclosure and including such departures from the present disclosure as come within known or customary practice within the art to which the disclosure pertains. It is intended that the specification and examples be considered as exemplary only, with a true scope and spirit of the disclosure being indicated by the following claims.
It will be understood that the present disclosure is not limited to the precise arrangements described above and shown in the drawings and that various modifications and changes may be made without departing from the scope thereof. The scope of the present disclosure is limited only by the appended claims.

Claims (10)

1. A method for judging the correctness of a three-dimensional space position of a meteorological observation station is characterized by comprising the following steps:
acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
determining correctness of longitude data and latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actually-measured altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
and judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
2. The method of claim 1, wherein performing a statistical analysis indicator calculation based on the air pressure data of the weather observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis indicator, a second statistical analysis indicator, and a third statistical analysis indicator comprises:
according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, the first statistical analysis index, the second statistical analysis index and the third statistical analysis index are calculated by adopting the following first formula:
Figure FDA0003024797480000021
wherein Index1(PRE) represents the first statistical analysis indicator, PRE (i) represents the air pressure data of the ith weather observation station, SimulatedPRE (i) represents the simulated air pressure data of the ith corresponding air numerical prediction mode, NUMPER represents the number of air pressure data of the weather observation stations, Index2(PRE) represents the second statistical analysis indicator, Index3(PRE) represents the third statistical analysis indicator, wherein the first statistical analysis index is used for representing the coincidence degree between the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the first statistical analysis index is used for representing the standard deviation of a difference value sequence of an air pressure data sequence of a meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period; and the third statistical analysis index is used for representing the root mean square value of a difference value sequence of the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in a preset time period.
3. The method of claim 1, wherein determining the correctness of the longitude data and the latitude data of the weather observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index comprises:
when the first statistical analysis index, the second statistical analysis index and the third statistical analysis index meet the following conditions, determining that longitude data and latitude data of the meteorological observation station are abnormal data, and otherwise, determining that the longitude data and the latitude data are normal data;
i Index1(PRE) | > 4 or
Index2(PRE) > 6 or
Index3(PRE)>6。
4. The method of claim 1, wherein said predicted altitude data is calculated using a second formula,
Figure FDA0003024797480000031
wherein H _ B represents simulated altitude data of an atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES _ B represents simulated barometric pressure data of the atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES represents barometric pressure data of the target weather-observing station, TEM1 represents air temperature data of the target weather-observing station, Alti represents measured altitude data of the target weather-observing station, and log _10 is a logarithmic function with a base 10.
5. The method of claim 1, wherein said determining the correctness of the measured altitude data of the target weather observation station based on the difference between the predicted altitude data and the measured altitude data comprises:
and when the difference value between the predicted altitude data and the actually measured altitude data is larger than a preset threshold value, judging that the actually measured altitude data of the target meteorological observation station is incorrect.
6. An apparatus for determining correctness of a three-dimensional spatial position of a weather observation station, the apparatus comprising:
the system comprises a preprocessing module, a data acquisition module and a data processing module, wherein the preprocessing module is used for acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
the first calculation module is used for performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
the determining module is used for determining the correctness of the longitude data and the latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
the second calculation module is used for calculating predicted altitude data of a target meteorological observation station with correct longitude data and latitude data according to the air pressure data, the air temperature data, the actually-measured altitude data of the target meteorological observation station, and the simulated altitude data and the simulated air pressure data of the corresponding atmospheric numerical prediction mode;
and the judging module is used for judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
7. The apparatus of claim 6, wherein the first computing module is configured to:
according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode, the first statistical analysis index, the second statistical analysis index and the third statistical analysis index are calculated by adopting the following first formula:
Figure FDA0003024797480000041
wherein Index1(PRE) represents the first statistical analysis indicator, PRE (i) represents the air pressure data of the ith weather observation station, SimulatedPRE (i) represents the simulated air pressure data of the ith corresponding air numerical prediction mode, NUMPER represents the number of air pressure data of the weather observation stations, Index2(PRE) represents the second statistical analysis indicator, Index3(PRE) represents the third statistical analysis indicator, wherein the first statistical analysis index is used for representing the coincidence degree between the air pressure data sequence of the meteorological observation station and the simulated air pressure data sequence of the corresponding atmospheric numerical prediction mode in the preset time period, the first statistical analysis index is used for representing the standard deviation of a difference value sequence of an air pressure data sequence of a meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period; the third statistical analysis index is used for representing the root mean square value of a difference value sequence of an air pressure data sequence of the meteorological observation station and a simulated air pressure data sequence of a corresponding atmospheric numerical prediction mode in a preset time period;
the determination module is to:
when the first statistical analysis index, the second statistical analysis index and the third statistical analysis index meet the following conditions, determining that longitude data and latitude data of the meteorological observation station are abnormal data, and otherwise, determining that the longitude data and the latitude data are normal data;
i Index1(PRE) | > 4 or
Index2(PRE) > 6 or
Index3(PRE)>6。
8. The apparatus of claim 6, wherein the predicted altitude data is calculated using a second formula,
Figure FDA0003024797480000051
wherein H _ B represents simulated altitude data of an atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES _ B represents simulated air pressure data of the atmospheric numerical prediction mode corresponding to the target weather-observing station, PRES represents air pressure data of the target weather-observing station, TEM1 represents air temperature data of the target weather-observing station, Alti represents measured altitude data of the target weather-observing station, and log _10 is a logarithmic function with a base 10;
the judging module is used for:
and when the difference value between the predicted altitude data and the actually measured altitude data is larger than a preset threshold value, judging that the actually measured altitude data of the target meteorological observation station is incorrect.
9. An apparatus for determining correctness of a three-dimensional spatial position of a weather observation station, the apparatus comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to:
acquiring air pressure data and air temperature data of a weather observation station which operates for more than preset time, and preprocessing the air pressure data and the air temperature data by using a threshold limiting method;
performing statistical analysis index calculation according to the air pressure data of the meteorological observation station and the simulated air pressure data of the corresponding atmospheric numerical prediction mode to determine a first statistical analysis index, a second statistical analysis index and a third statistical analysis index;
determining correctness of longitude data and latitude data of the meteorological observation station according to the first statistical analysis index, the second statistical analysis index and the third statistical analysis index;
for a target meteorological observation station with correct longitude data and latitude data, calculating predicted altitude data according to air pressure data, air temperature data, actually-measured altitude data of the target meteorological observation station, and simulated altitude data and simulated air pressure data of a corresponding atmospheric numerical prediction mode;
and judging the correctness of the actually measured altitude data of the target meteorological observation station according to the difference value between the predicted altitude data and the actually measured altitude data.
10. A computer-readable storage medium having stored thereon computer instructions, which, when executed by a processor, carry out the steps of the method according to any one of claims 1 to 5.
CN202110413289.XA 2021-04-16 2021-04-16 Method and device for judging correctness of three-dimensional space position of meteorological observation station Active CN113075752B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110413289.XA CN113075752B (en) 2021-04-16 2021-04-16 Method and device for judging correctness of three-dimensional space position of meteorological observation station

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110413289.XA CN113075752B (en) 2021-04-16 2021-04-16 Method and device for judging correctness of three-dimensional space position of meteorological observation station

Publications (2)

Publication Number Publication Date
CN113075752A true CN113075752A (en) 2021-07-06
CN113075752B CN113075752B (en) 2022-05-17

Family

ID=76617823

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110413289.XA Active CN113075752B (en) 2021-04-16 2021-04-16 Method and device for judging correctness of three-dimensional space position of meteorological observation station

Country Status (1)

Country Link
CN (1) CN113075752B (en)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635281A (en) * 2015-02-17 2015-05-20 南京信息工程大学 Method for controlling data quality of automatic meteorological station based on strong weather process correction
CN104778248A (en) * 2015-04-14 2015-07-15 中国气象局气象探测中心 Network information management system for comprehensive meteorological observation operation monitoring system
CN107037439A (en) * 2017-03-28 2017-08-11 武汉大学 For the laser ceilometer atmosphere delay range error modification method of land target
CN108227041A (en) * 2017-12-27 2018-06-29 中国海洋大学 Horizontal visibility forecasting procedure based on website measured data and model results
US20190368872A1 (en) * 2018-06-04 2019-12-05 Nextnav, Llc Systems and methods for determining which reference-level pressures are used when estimating an altitude of a mobile device
CN111913948A (en) * 2020-09-29 2020-11-10 江西省气象信息中心(江西省气象培训中心、江西省农村经济信息中心) Processing method and system for collecting meteorological observation data and computer equipment
CN112285807A (en) * 2019-11-12 2021-01-29 北京嘉韵楷达气象科技有限公司 Meteorological information prediction method and device
CN112380499A (en) * 2020-11-11 2021-02-19 成都市气象局 Regional automatic weather station objective clustering method
CN112612781A (en) * 2019-12-11 2021-04-06 北京金风慧能技术有限公司 Data correction method, device, equipment and medium

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104635281A (en) * 2015-02-17 2015-05-20 南京信息工程大学 Method for controlling data quality of automatic meteorological station based on strong weather process correction
CN104778248A (en) * 2015-04-14 2015-07-15 中国气象局气象探测中心 Network information management system for comprehensive meteorological observation operation monitoring system
CN107037439A (en) * 2017-03-28 2017-08-11 武汉大学 For the laser ceilometer atmosphere delay range error modification method of land target
CN108227041A (en) * 2017-12-27 2018-06-29 中国海洋大学 Horizontal visibility forecasting procedure based on website measured data and model results
US20190368872A1 (en) * 2018-06-04 2019-12-05 Nextnav, Llc Systems and methods for determining which reference-level pressures are used when estimating an altitude of a mobile device
CN112285807A (en) * 2019-11-12 2021-01-29 北京嘉韵楷达气象科技有限公司 Meteorological information prediction method and device
CN112612781A (en) * 2019-12-11 2021-04-06 北京金风慧能技术有限公司 Data correction method, device, equipment and medium
CN111913948A (en) * 2020-09-29 2020-11-10 江西省气象信息中心(江西省气象培训中心、江西省农村经济信息中心) Processing method and system for collecting meteorological observation data and computer equipment
CN112380499A (en) * 2020-11-11 2021-02-19 成都市气象局 Regional automatic weather station objective clustering method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
闫荞荞 等: "观测数据中气象站位置信息检验及应对策略研究", 《气象科学》 *

Also Published As

Publication number Publication date
CN113075752B (en) 2022-05-17

Similar Documents

Publication Publication Date Title
GB2613296A (en) Greenhouse gas emission monitoring systems and methods
CN111896914A (en) Cooperative positioning method, device, equipment and storage medium
CN109635053B (en) Map quality inspection method, device, system and storage medium
CN112596545B (en) Multispectral-based water pollution source head unmanned aerial vehicle troubleshooting method and system and storage medium
CN111814740B (en) Pointer instrument reading identification method, device, computer equipment and storage medium
CN113218537B (en) Training method, training device, training equipment and training storage medium for temperature anomaly detection model
CN110619497B (en) Address checking method, device, electronic equipment and storage medium
CN113065538A (en) Pressure sensor detection method, device and equipment based on image recognition
CN115994248A (en) Data detection method and system for valve faults
CN106599367A (en) Method for detecting abnormal state of spacecraft
CN112686322A (en) Part difference identification method, device, equipment and storage medium
CN114492629A (en) Abnormality detection method, abnormality detection device, electronic apparatus, and storage medium
CN109523129B (en) Method for fusing information of multiple sensors of unmanned vehicle in real time
CN117671303B (en) Stereoscopic ranging hidden danger early warning method and system based on twin network
CN113075752B (en) Method and device for judging correctness of three-dimensional space position of meteorological observation station
CN112541161B (en) Regional multi-source precipitation data quality control method and system
CN115294373A (en) Equipment parameter information updating method and updating system
CN111679298B (en) Integrity monitoring method and device of navigation system and electronic equipment
CN115222691A (en) Image defect detection method, system and related device
CN114841255A (en) Detection model training method, device, equipment, storage medium and program product
CN114192583A (en) Scada platform-based strip steel rolling process quality monitoring method and system
CN114488247A (en) Method for analyzing mobility of equipment based on high-precision Beidou differential positioning
CN113624530A (en) Analog measuring point association method, device, equipment and storage medium
CN112243193A (en) Indoor positioning method and device, computer equipment and readable storage medium
CN112528417A (en) Aircraft semi-physical simulation evaluation method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant