CN111401632B

CN111401632B - Weather forecast method and related device

Info

Publication number: CN111401632B
Application number: CN202010174457.XA
Authority: CN
Inventors: 梁卉林; 袁兴德; 陈斌; 王铁强; 曾垂宽; 王志军; 葛文涛; 兰水泉; 刘明全
Original assignee: China Resource Power Technology Research Institute
Current assignee: China Resource Power Technology Research Institute
Priority date: 2020-03-13
Filing date: 2020-03-13
Publication date: 2023-07-21
Anticipated expiration: 2040-03-13
Also published as: CN111401632A

Abstract

The embodiment of the application discloses a weather forecast method, which is used for solving the problem that weather forecast data with higher resolution is difficult to obtain in limited time and limited computing resources as the time for waiting for completion of GFS initial field downloading is longer. The method comprises the following steps: acquiring a configuration file of weather forecast, wherein the configuration file comprises an aging starting point, an aging end point and a time step of weather forecast data, the time step separates a time period formed by connecting the aging starting point and the aging end point into a plurality of aging points, the aging points comprise the aging starting point and the aging end point, GFS initial field data corresponding to the aging points are acquired one by one according to the configuration file, preprocessing is carried out on the GFS initial field data to obtain target GFS initial field data, weather forecast integral parameters are obtained, integral calculation is carried out according to the weather forecast integral parameters and the target GFS initial field data to obtain weather forecast data, and weather elements are extracted from the weather forecast data.

Description

Weather forecast method and related device

Technical Field

The embodiment of the application relates to the field of meteorological data processing, in particular to a weather forecast method and a related device.

Background

In the current weather forecast products, there is a block that the development cost is high and the computing resources are limited, the common practice is to build a business numerical weather forecast product manufacturing platform based on a weather research and forecast model (WRF) of a physical model, and manufacture and send weather forecast data of a specific area at regular time every day by using a Global Forecast System (GFS) manufactured and issued by a national environmental forecast center (nationalcenters forenvironmentalprediction, NCEP).

However, the GFS initial field released by the national environment prediction center is thicker in resolution, and in the prior art, after the GFS initial field is downloaded, the GFS initial field needs to be calculated by a specific algorithm to obtain weather prediction data with higher resolution. The GFS initial field is global coarse grid data which is manufactured and released by the national environment forecasting center every day at regular time. Because the GFS initial field contains a plurality of data packets, each data packet represents the GFS initial field of a different time point, when the network transmission speed is not ideal, the progress of downloading the GFS initial field can be dragged to be long, and a great amount of time is required to wait for downloading the GFS initial field. And because the weather forecast data has strong timeliness, on the premise of limited computing resources, the longer the time spent waiting for completion of GFS initial field downloading is, the shorter the time left for obtaining the weather forecast data with higher resolution through specific algorithm computing is, and further the weather forecast data with higher resolution is difficult to obtain in limited time and limited computing resources.

Disclosure of Invention

The embodiment of the application provides a weather forecast method and a related device, which are used for solving the problem that the longer the waiting time for completion of GFS initial field downloading is, the more difficult the weather forecast data with higher resolution is obtained in limited time and limited computing resources.

A first aspect of the present application provides a weather forecast method, including:

acquiring a configuration file of weather forecast, wherein the configuration file comprises an aging starting point, an aging end point and a time step of weather forecast data, the time step separates a time period formed by connecting the aging starting point and the aging end point to form a plurality of aging points, and the aging points comprise the aging starting point and the aging end point;

acquiring GFS initial field data corresponding to the aging points one by one according to the configuration file;

preprocessing the GFS initial field data to obtain target GFS initial field data;

acquiring weather forecast integral parameters;

performing integral calculation according to the weather forecast integral parameters and the target GFS initial field data to obtain weather forecast data;

and extracting meteorological elements from the weather forecast data.

Optionally, the step of obtaining GFS initial field data corresponding to the aging points one by one according to the configuration file includes:

Determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

acquiring a URL link of GFS initial field data corresponding to the current aging point;

downloading GFS initial field data corresponding to the current aging point according to the URL link;

judging whether the downloading of the GFS initial field data is completed or not;

if the downloading is completed, the GFS initial field data is moved to a specified path, and GFS initial field data corresponding to the current aging point is obtained;

adding the time step to the current aging point to obtain a later aging point of the current aging point;

judging whether the latter aging point is larger than the aging end point or not;

and if the latter aging point is smaller than or equal to the aging end point, regarding the latter aging point as a current aging point, and triggering the step of acquiring the URL link of the GFS initial field data corresponding to the current aging point.

Optionally, the preprocessing the GFS initial field data to obtain target GFS initial field data includes:

obtaining topographic data of a weather forecast area;

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

if the latter aging point is smaller than or equal to the aging end point, judging whether the downloading of the current GFS initial field corresponding to the current aging point and the latter GFS initial field corresponding to the latter aging point is completed;

if the downloading is completed, analyzing the current GFS initial field and the subsequent GFS initial field to obtain two analyzed GFS initial fields;

matching the two analyzed GFS initial fields with the topographic data to obtain two target GFS initial field data;

deleting the current GFS initial field;

taking the latter time effect point as a current time effect point;

triggering the step of adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.

Optionally, the performing integral calculation according to the weather forecast integral parameter and the target GFS initial field data, to obtain weather forecast data includes:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

if the latter aging point is smaller than or equal to the aging end point, judging whether current target GFS initial field data corresponding to the current aging point and latter target GFS initial field data corresponding to the latter aging point exist or not;

if the current aging points exist, judging whether the current aging points are aging starting points or not;

if the current aging point is an aging starting point, integrating calculation is carried out by using target GFS initial field data corresponding to the aging starting point, target GFS initial field data corresponding to the latter aging point and the weather forecast integral parameter, so as to obtain weather forecast data from the aging starting point to the latter aging point;

if the current aging point is not the aging starting point, integrating calculation is carried out on the end point value of the weather forecast data obtained by integrating calculation with the previous aging point of the current aging point, the target GFS initial field data corresponding to the next aging point and the weather forecast integration parameter, so as to obtain weather forecast data from the current aging point to the next aging point;

deleting the current target GFS initial field data;

taking the initial field data of the target GFS as the initial field data of the current target GFS;

Optionally, the extracting the weather element from the weather forecast data includes:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

judging whether weather forecast data corresponding to the current aging point exists or not;

if so, extracting preset meteorological elements from weather forecast data corresponding to the current aging point;

storing the meteorological elements to a database;

deleting weather forecast data corresponding to the current aging point;

if the latter aging point is smaller than or equal to the aging end point, the latter aging point is taken as the current aging point;

triggering and judging whether weather forecast data corresponding to the current aging point exist or not.

Optionally, the GFS initial field data is GFS initial field data of a certain time point in a global coarse grid initial field which is manufactured and released by the national environment prediction center every day at regular time.

A second aspect of the present application provides a weather forecast system comprising:

the weather forecast system comprises a first acquisition unit, a second acquisition unit and a weather forecast unit, wherein the first acquisition unit is used for acquiring a configuration file of weather forecast, the configuration file comprises an ageing starting point, an ageing end point and a time step of weather forecast data, the time step separates a time period formed by connecting the ageing starting point and the ageing end point into a plurality of ageing points, and the ageing points comprise the ageing starting point and the ageing end point;

the second acquisition unit is used for acquiring GFS initial field data corresponding to the aging points one by one according to the configuration file;

the preprocessing unit is used for preprocessing the GFS initial field data to obtain target GFS initial field data;

the third acquisition unit is used for acquiring weather forecast integral parameters;

the calculation unit is used for carrying out integral calculation according to the weather forecast integral parameter and the target GFS initial field data to obtain weather forecast data;

and the extraction unit is used for extracting meteorological elements from the weather forecast data.

Optionally, when the second obtaining unit obtains GFS initial field data corresponding to the aging point one by one according to the configuration file, the second obtaining unit is specifically configured to:

determining an aging starting point in the configuration file;

Taking the aging starting point as a current aging point;

Optionally, the preprocessing unit is configured to preprocess the GFS initial field, and when obtaining target GFS initial field data, specifically configured to:

obtaining topographic data of a weather forecast area;

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current GFS initial field;

taking the latter time effect point as a current time effect point;

Optionally, the calculating unit performs integral calculation on the target GFS initial field data according to the weather forecast integral parameter, and is specifically used for when obtaining weather forecast data:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current target GFS initial field data;

Optionally, the extracting unit extracts the weather element from the weather forecast data includes:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

storing the meteorological elements to a database;

deleting weather forecast data corresponding to the current aging point;

Optionally, the GFS initial field data is initial field data of a certain time point in a global coarse grid initial field which is manufactured and released by the national environment prediction center every day at regular time.

A third aspect of the present application provides a weather forecast system comprising:

a processor, a memory, a bus, and an input/output device;

the processor is connected with the memory and the input and output equipment;

the bus is respectively connected with the processor, the memory and the input and output equipment;

the processor performs the method of any of the preceding first aspects.

A fourth aspect of the present application provides a computer storage medium having instructions stored therein, which when executed on a computer, cause the computer to perform the method of any of the preceding first aspects.

A fifth aspect of the present application provides a computer program product which, when executed on a computer, causes the computer to perform the method of any of the preceding first aspects.

From the above technical solutions, the embodiments of the present application have the following advantages:

according to the method and the device, the GFS initial field data of the required aging points are downloaded one by one through the indication in the preset configuration file, after the GFS initial field data of the corresponding aging points are downloaded, the GFS initial field data of the aging points are subjected to preprocessing, integration and other calculation steps immediately, and the preprocessing, integration and other calculation steps are not needed after the GFS initial field data of all the aging points are downloaded like the traditional GFS initial field data are downloaded, so that the problem that the longer the waiting time for the GFS initial field data to be downloaded is, the more difficult to obtain weather forecast data with higher resolution is solved.

Drawings

FIG. 1 is a schematic flow chart of an weather forecast method according to an embodiment of the present application;

fig. 2 is a schematic flow chart of a weather forecast method related to GFS initial field data downloading in the embodiment of the present application;

FIG. 3 is a schematic flow chart of an weather forecast method related to a pretreatment process according to an embodiment of the present application;

FIG. 4 is a schematic flow chart of the method for forecasting weather about the integral calculation process in the embodiment of the present application;

FIG. 5 is a schematic flow chart of a weather forecast method for extracting meteorological elements according to an embodiment of the present application;

FIG. 6 is a schematic diagram of a weather forecast system according to an embodiment of the present application;

fig. 7 is another schematic structural diagram of an weather forecast system according to an embodiment of the present application.

Detailed Description

The embodiment of the application provides a weather forecast method, which is used for solving the problem that weather forecast data with higher resolution is difficult to obtain in limited time and limited computing resources as the time for waiting for completion of GFS initial field downloading is longer.

It should be noted that, for ease of understanding and implementation of the embodiments of the present application, the 5 executable files geogrid.exe, ungrib.exe, metagrid.exe, real.exe, wrf.exe mentioned in the present application are compiled of WRF installation packages that can be downloaded through an official website (https:// www2.Mmm.ucar. Edu/WRF/users/download/get_sources. Html#wps).

Referring to fig. 1, an embodiment of a weather forecast method of the present application includes:

101. the method comprises the steps of obtaining a configuration file of weather forecast, wherein the configuration file comprises an aging starting point, an aging end point and a time step of weather forecast data, and the time step separates a time period formed by connecting the aging starting point and the aging end point to form a plurality of aging points, and the aging points comprise the starting point and the aging end point.

The configuration file of the application refers to a configuration file which is edited and manufactured by self according to the weather forecast requirement, and the file stores information such as forecast starting time, start-stop forecast aging, time step, file storage path, working path during calculation and the like of a weather forecast system, wherein the configuration file is a parameter configuration file of the whole weather forecast system.

Since the GFS initial field issued by the national environment prediction center contains a plurality of data packets, the GFS initial field is GFS initial field data of global coarse grid (0.5 degree x 0.5 degree) issued by NCEP which is manufactured at 4 times per day (UTC time 00, 06, 12, 18, beijing time 08, 14, 22, next day 02) and issued, the forecast aging is 0-384 hours, and the time interval is 3 hours (namely, 129 data packets are included in the initial field data issued by each manufacturing, and the suffix names are f000, f003, … … and f384 respectively). Of the 129 packets, only the first packet (suffix f 000) is the actual GFS initial field measurement data, and the last 128 packets are calculated GFS initial field measurement data at different future times based on the first packet. The quantity of data packets possibly needed by different weather forecast requirements is different, GFS initial field data needed by the weather forecast is needed to be obtained according to related configuration files, the configuration files comprise an aging starting point, an aging end point and a time step of the weather forecast data, the time step separates a time period formed by connecting the aging starting point and the aging end point to form a plurality of aging points, and the aging points comprise the aging starting point and the aging end point.

If the global coarse grid (0.5 degree x 0.5 degree) GFS initial field data issued by the NCEP (corresponding to beijing time 08 at UTC time 00) is forecast aged for 0-384 hours, the time interval is 3 hours (namely 129 data packets are generated for the issued GFS initial field data, and the suffix names are f000, f003, … … and f384 respectively). Just when the configuration file requires that the ageing starting point of the acquired GFS initial field data is UTC time 00, when the ageing starting point corresponds to Beijing time 08, the ageing end point is 384 hours after the Beijing time 08, the time step is 3 hours, and at the moment, the ageing starting point and the ageing end point are separated by the time step to form a plurality of ageing points which are in one-to-one correspondence with NCEP (NCEP) to manufacture 129 issued global coarse grid (0.5 degree multiplied by 0.5 degree) GFS initial field data. It will be appreciated that the configuration file may be generated from GFS initial field data of 129 global coarse grids (0.5 ° ×0.5 °) issued by the NCEP generation, so that a subsequent step obtains GFS initial field data of part or all of the global coarse grids (0.5 ° ×0.5 °), without limitation.

In addition, data such as integration start-stop time and the like required by the preprocessing process need to be set and written into a text file namelist.wps for calling by the preprocessing process (geogrid.exe, ungrib.exe and metgrid.exe).

102. And acquiring GFS initial field data corresponding to the aging points one by one according to the configuration file.

After knowing the data packet of the GFS initial field data corresponding to the aging point required for the weather forecast in step 101, the data packet corresponding to the aging point may be obtained one by one.

103. Preprocessing the GFS initial field to obtain target GFS initial field data.

Since the GFS initial field obtained in step 102 is not necessarily suitable for the weather forecast system in format or form, the GFS initial field needs to be preprocessed to obtain target GFS initial field data meeting the usage standard. This step is performed after GFS initial field data corresponding to one aging point is obtained in step 102.

104. And obtaining weather forecast integral parameters.

In order to obtain weather forecast data with higher resolution, besides obtaining the target GFS initial field data meeting the use standard in step 103, weather forecast integral parameters are required to be obtained, where the weather forecast integral parameters are generally integral conditions set according to the actual needs of weather forecast. The integration condition is to set parameters such as integration start-stop time, boundary condition of integration area, spatial resolution, parameterization scheme, IO setting and the like required by the integration calculation process, and write a text file namellist.

105. And carrying out integral calculation on the target GFS initial field data according to the weather forecast integral parameters to obtain weather forecast data.

After obtaining the target GFS initial field in step 103 and the weather forecast integral parameter in step 104, integral calculation can be performed on the target GFS initial field data according to the weather forecast integral parameter to obtain weather forecast data with higher resolution.

106. And extracting meteorological elements from the weather forecast data.

After integrating and calculating the weather forecast data with higher resolution in the step 105, the weather elements can be extracted from the weather forecast data, and the weather elements are the information such as wind direction, wind speed, air temperature, air pressure, relative humidity and the like required by weather forecast in the application.

Referring to fig. 2, an embodiment of the weather forecast method of the present application regarding GFS initial field downloading includes:

201. an age start point in the configuration file is determined.

The aging start point of GFS initial field data to be downloaded is determined in the configuration file of fig. 1, and is generally determined according to the requirements of the weather forecast system. In practice, GFS initial fields are global coarse grid (0.5 ° ×0.5 °) initial field data which are regularly produced 4 times a day by NCEP (UTC times 00, 06, 12, 18, corresponding to beijing times 08, 14, 22, 02 next day), forecast aging is 0-384 hours, time intervals are 3 hours (i.e. there are 129 data packets for initial field data which are produced each time, suffix names are f000, f003, … …, f384, respectively), and in this embodiment, the aging start point of GFS initial field data which is the latest at the current time point is generally selected. For example, when the current time is Beijing time 09, the latest aging start point is Beijing time 08, and so on.

202. And taking the aging starting point as the current aging point.

After determining the ageing starting point in step 201, the ageing starting point is taken as the current ageing point.

203. And acquiring a URL link of GFS initial field data corresponding to the current time point.

URL links refer to Uniform Resource Locators (URLs), which are Uniform Resource Locators (URL) of the world wide web, that is, network addresses. Specifically, this step is used to obtain the network address of the GFS initial field data that downloads the current aging point. The network address of the GFS initial field data of the starting point is typically obtained in a configuration file, and the network address of the subsequent GFS initial field data associated with the starting point can be obtained by correlating the network address of the GFS initial field data of the starting point.

204. And downloading GFS initial field data corresponding to the current aging point according to the URL link.

After obtaining the network address of the GFS initial field data of the current aging point in step 203, the GFS initial field data corresponding to the current aging point may be downloaded through the URL link.

205. Whether the GFS initial field data is downloaded is determined, if so, step 206 is executed, and if not, the download is continued.

Since the data amount of GFS initial field data produced by the NCEP at 4 times per day is relatively large, it takes a considerable time for the complete download to be completed in terms of current network download speed. According to the method and the device, the data packets of GFS initial field data manufactured by NCEP are downloaded one by one as required according to the requirements of the configuration file, and the downloading state of the data packets is judged in real time. Generally, by using the Vtable of the g2print data, compared with the Vtable of the WRF, if the Vtable is consistent, the GFS initial field data of the current point of time is determined to be downloaded, and if the Vtable is inconsistent, the GFS initial field data of the current point of time is determined to be not downloaded.

206. And moving the downloaded GFS initial field data to a specified path to obtain GFS initial field data corresponding to the current aging point.

After determining that the GFS initial field downloading of the current aging point is completed in step 205, moving the downloaded GFS initial field data to a specified path to obtain GFS initial field data corresponding to the current aging point.

207. And adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.

After the current aging point is obtained in step 206, the GFS initial field data of the next aging point of the current aging point needs to be downloaded, and at this time, information of the next aging point needs to be obtained, 129 data packets are made and issued according to the NCEP initial field data, and the suffix names are respectively (f 000, f003, … …, f 384) rules, that is, the time step is 3, and the next aging point of the current aging point can be obtained according to the suffix name of the current aging point and the time step.

208. If the subsequent aging point is determined to be greater than the aging end point, step 209 is executed, and if the subsequent aging point is greater than the aging end point, step 210 is executed.

Since the last aging point is calculated by the accumulating method in step 207, and the initial field data made and issued by the NCEP has 129 data packets, the suffix names (f 000, f003, … …, f 384) are limited, so as to avoid errors when accumulating to the aging end point (suffix name f 384) and further accumulating, the judgment of this step is specially increased. Of course, this step is also applicable to the case where the user's profile requires only downloading the number of packets less than the initial field data created and issued by the NCEP, and thus, no more packets need to be downloaded, and the waste of network resources is reduced, for example, if the profile specifies that the aging end point is the suffix f003, then only GFS initial field data with the suffix f003 is downloaded, although there are 129 packets of initial field data created and issued by the NCEP.

209. The latter aging point is regarded as the current aging point, triggering execution of step 203.

After the next aging point satisfying the condition is obtained in step 208, the next aging point is regarded as the current aging point, and step 203 is triggered to be executed to obtain the URL link of the GFS initial field corresponding to the current aging point.

210. The download is stopped.

When it is determined in step 208 that the later aging point is greater than the aging end point, it is verified that the GFS initial field data corresponding to all the aging points has been downloaded according to the requirement, and the downloading of the GFS initial field data should be stopped.

The weather forecast method in fig. 2 can operate independently and modularly about GFS initial field downloading, and only depends on the requirements of configuration files, so that GFS initial field data at different timepoints can be downloaded one by one.

Referring to fig. 3, an embodiment of the weather forecast method related to preprocessing in the embodiment of the present application includes:

301. and obtaining the topographic data of the weather forecast area.

And (3) interpolating and manufacturing the topographic data geo_em by using other accessible topographic data of the officially issued WRF topographic data set through the geogrid.exe execution file, wherein generally, after the setting of the follow-up integrating area, if the weather forecast demand is unchanged and the WRF integrating area is unchanged, the topographic data geo_em does not need to be manufactured again during daily business forecast. The so-called terrain data geo em reflects the terrain conditions of the area where the weather forecast is required.

302. An age start point in the configuration file is determined.

The execution of this step is similar to the execution of step 201 in fig. 2 and will not be described again here.

303. And taking the aging starting point as the current aging point.

The execution of this step is similar to the execution of step 202 in fig. 2 and will not be described again here.

304. And adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.

After determining the current aging point in step 303, the GFS initial field data of the next aging point of the current aging point needs to be preprocessed, information of the next aging point needs to be acquired at this time, 129 data packets are made and issued according to the NCEP initial field data, the suffix names are respectively (f 000, f003, … …, f 384) rules, that is, the time step is 3, and the next aging point of the current aging point can be obtained according to the suffix name of the current aging point and the time step.

305. Whether the subsequent aging point is greater than the aging end point is determined, if so, step 306 is executed, and if so, step 311 is executed.

The execution of this step is similar to the execution of step 208 of fig. 2, and will not be described in detail here.

306. Judging whether the downloading of the current GFS initial field data corresponding to the current time point and the subsequent GFS initial field data corresponding to the subsequent time point is completed, if so, executing step 307, and if not, continuing to download.

Specifically, it is determined whether the current GFS initial field data corresponding to the current time point and the subsequent GFS initial field data corresponding to the subsequent time point exist in the designated path in step 206 of fig. 2.

307. Analyzing the current GFS initial field data and the subsequent GFS initial field data to obtain two analyzed GFS initial fields.

Specifically, the current GFS initial field data and the subsequent GFS initial field data downloaded in step 306 are in a grib format, and the current GFS initial field data and the subsequent GFS initial field data in the grib format are parsed by using an ungrab. Exe execution FILE, so as to obtain current FILE data corresponding to the current GFS initial field data and subsequent FILE data corresponding to the subsequent GFS initial field data.

308. And matching the two analyzed GFS initial fields with the topographic data to obtain two target GFS initial field data.

Specifically, the current FILE data corresponding to the current GFS initial field data obtained in step 307 and the subsequent FILE data corresponding to the subsequent GFS initial field data are respectively interpolated into the horizontal forecast grid of the topographic data in step 301 by using the meta executable FILE, so as to obtain corresponding meta_em data. And storing the met_em data corresponding to the current time point and the later time point to the designated path required by the configuration file.

309. The current GFS initial field data is deleted.

After the met_em data corresponding to the current aging point is stored in the designated path required by the configuration file in step 308, in order to facilitate the subsequent step, confusion is not easily generated when the subsequent aging point is used as the current aging point, and the current GFS initial field data needs to be deleted.

310. The next time point is taken as the current time point, and the step 304 is triggered and executed.

After deleting the current GFS initial field data in step 309, the next time point is taken as the current time point, and step 304 is triggered and executed to add a time step to the current time point, so as to obtain the next time point of the current time point.

311. The pretreatment is stopped.

When step 305 determines that the latter aging point is greater than the aging end point, it proves that the preprocessing of the GFS initial field data corresponding to all the aging points has been completed according to the preprocessing requirement, and at this time, the preprocessing of the GFS initial field data should be stopped.

The preprocessing of the weather forecast method with respect to GFS initial field data in fig. 3 may be modularized and may operate independently, and only depends on the requirements in the configuration file and the GFS initial field data of the current aging point after downloading, so that the preprocessing may be immediately started when the GFS initial field data corresponding to the current aging point is obtained in step 206 in fig. 2.

Referring to fig. 4, an embodiment of an weather forecast method related to an integral calculation process in the embodiment of the present application includes:

401. an age start point in the configuration file is determined.

402. And taking the aging starting point as the current aging point.

403. And adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.

After the acquisition of the current aging point is completed in step 402, since the process of integral calculation needs not only the target GFS initial field data of the current aging point but also the GFS initial field data of the subsequent aging point, the GFS initial field data of the subsequent aging point is used for defining the target GFS initial field data of the current aging point and is calculated until the GFS initial field data of the subsequent aging point is integrated according to a preset integral rule. This means that the integral calculation is to perform the integral calculation of weather forecast data between two aging points based on the target GFS initial field data of the two aging points as a boundary condition.

404. It is determined whether the subsequent aging point is greater than the aging end point, if it is less than or equal to the aging end point, step 405 is executed, and if it is greater than, step 411 is executed.

405. Judging whether the current target GFS initial field data corresponding to the current aging point and the target GFS initial field data of the next aging point exist or not, if so, executing step 406, and if not, continuing executing the step.

Specifically, it is determined whether there is met_em data corresponding to the current aging point and met_em data corresponding to the next aging point in the designated path in step 308 in fig. 3.

406. Whether the current aging point is an aging start point is determined, if yes, step 407 is executed, and if no, step 408 is executed.

From the foregoing description, it can be seen that, of the 129 packets that are made at the NCEP timing, only the first packet (suffix name f 000) is the actual GFS initial field measurement data, and the last 128 packets are the GFS initial field measurement data calculated from the first packet at different future times. Different integral calculation strategies are employed for this case.

407. And integrating calculation is carried out by using target GFS initial field data corresponding to the aging starting point, the next target GFS initial field data corresponding to the next aging point and the weather forecast integral parameters, so as to obtain weather forecast data from the aging starting point to the next aging point.

After determining in step 406 that the obtained target GFS initial field data of the current aging point is the target GFS initial field data of the aging start point, specifically, interpolating met_em data of the current aging point to the vertical prediction height layer of the topographic data in step 301 through a real. Exe executable file to obtain wrfbdy, wrfinput data of the aging start point, wherein wrfbdy is a lower boundary (or under-pad) condition, wrfbput is a side boundary condition, and an integral region of the in-situ topographic data is boundary in four directions of southwest and northwest. Similarly, the meta_em data of the next time point is interpolated into the vertical prediction height layer of the topographic data in step 301 by the real. Exe executable file to obtain wrfbdy, wrfinput data of the next time point, wherein wrfbdy is a lower boundary (or under-pad surface) condition, wrfput is a side boundary condition, and an integral area of the in-situ topographic data is bounded in four directions of southeast and northwest.

Combining the text file namellist.input of the integration condition in the step 104, wherein the text file namellist.input of the integration condition records parameters such as start-stop time, time resolution and the like required by integration; reading a lower boundary condition and a side boundary condition from wrfbdy, wrfinput data in a wrf.exe executable file, starting with target GFS initial field data of an aging starting point and ending with target GFS initial field of a later aging point of the aging starting point, performing integral calculation, so as to calculate weather forecast data between two aging points, and storing the weather forecast data between the aging starting point and the later aging point into a designated path in a configuration file. It will be appreciated that in the integration calculation using two time points as the time start point and the time end point, if weather forecast data with higher time resolution accuracy (two time points with smaller time interval) is desired, the more times the integration calculation is needed, the longer the time for the integration calculation is needed. The weather forecast data after integral calculation can output wrforut data according to IO setting.

408. And carrying out integral calculation by using the end point value of the weather forecast data obtained by the integral calculation of the previous time point of the current time point, the target GFS initial field data corresponding to the next time point and the weather forecast integral parameter to obtain the weather forecast data from the current time point to the next time point.

This step is performed as in 407, except that the target GFS initial field data of the current aging point in step 407 is replaced with the end point value of the weather forecast data obtained by participating in the integration calculation of the previous aging point of the current aging point. And storing weather forecast data from the previous aging point to the next aging point to a designated path in the configuration file.

409. The current target GFS initial field data is deleted.

After obtaining the weather forecast data from the current aging point to the next aging point in step 408 or step 409, it may be necessary to perform integral calculation with the target initial field data of the next aging point as the time starting point (it may not be necessary to use the target initial field of the next aging point as the time starting point as shown in step 408), and the target initial field data of the next aging point as the time ending point, so that the current target GFS initial field needs to be deleted to avoid confusion.

410. The next target GFS initial field data is taken as the current target GFS initial field data, triggering execution of step 403.

After deleting the current target GFS initial field data in step 409, the next target GFS initial field data may be used as the current target GFS initial field data, and step 403 is triggered to be executed to add a time step to the current aging point, so as to obtain the next aging point of the current aging point.

412. The integration calculation is stopped.

When it is determined in step 404 that the latter aging point is greater than the aging end point, it is verified that the integral calculation of the target GFS initial field data corresponding to all the aging points has been completed according to the integral calculation requirement, and at this time, the integral calculation of the target GFS initial field data should be stopped.

Referring to fig. 5, an embodiment of an weather forecast method related to extracting meteorological elements in an embodiment of the present application includes:

501. an age start point in the configuration file is determined.

502. And taking the aging starting point as the current aging point.

503. Judging whether weather forecast data corresponding to the current aging point exists or not, if so, executing step 504, and if not, continuing executing the step.

Specifically, this step mainly determines whether the weather forecast data corresponding to the current aging point exists in the designated path in step 406 in fig. 4.

504. And extracting preset meteorological elements from weather forecast data corresponding to the current aging point.

Specifically, after the weather forecast data between two time points with higher resolution is obtained by integral calculation in step 406 in fig. 4, a certain time point of the weather forecast data can be linearly interpolated to a height layer required by the topographic data, and then the required weather elements are directly extracted or obtained by calculation, wherein the weather elements are the forecast results of wind direction, wind speed, air temperature, air pressure, relative humidity and the like of the weather forecast in the application, and the weather element values on a required single time point can be extracted according to the requirement and output to a database for power forecast and call. In practical application, nc data is generated by running a post-processing script that is edited by itself, interpolating to a desired height layer and calculating desired weather elements.

505. The meteorological elements are stored to a database.

Specifically, the nc data generated in step 504 is stored in a designated database.

506. And deleting the weather forecast data of the current aging point.

In step 505, after the weather elements of the weather forecast data of the current aging point are extracted and stored in the database, the weather elements of the weather forecast data of the current aging point need to be extracted as the weather forecast data of the current aging point in the next time, and in order to avoid confusion, the weather forecast data of the current aging point need to be deleted after step 505.

507. And adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.

Since the weather forecast data between two aging points is obtained in step 406 in fig. 4, the weather forecast data between two aging points also needs to be known, that is, the current aging point is added with a time step to obtain the latter aging point of the current aging point on the premise of knowing the current aging point.

508. Whether the subsequent aging point is greater than the aging end point is determined, if so, step 509 is executed, and if so, step 510 is executed.

509. The execution of step 503 is triggered by taking the latter aging point as the current aging point.

After deleting the weather forecast data of the current aging point in step 508, the next aging point may be used as the current aging point, and step 503 is triggered to be executed to add a time step to the current aging point, so as to obtain the next aging point of the current aging point.

510. The extraction process is stopped.

When it is confirmed in step 508 that the latter aging point is greater than the aging end point, it is verified that the extraction of the weather elements of the weather forecast data between all the two aging points has been completed according to the requirements of the configuration file, and at this time, the extraction of the weather elements of the weather forecast data between all the two aging points should be stopped.

The weather forecast method in fig. 5 can operate independently and modularly with respect to the process of extracting the weather elements, and only depends on the requirement of configuration query price and the current weather forecast data in fig. 4, so that the weather forecast data can be extracted immediately after the weather forecast data is obtained in step 406 in fig. 4.

According to the modularized independent operation mode, the linear flow of the traditional operation weather forecast system is changed, all nodes are modularized, the time for waiting for the former step in the subsequent step is reduced, more integral calculation time is reserved for the weather forecast system, and weather forecast data with higher time resolution is obtained under the condition of limited calculation resources of a weather forecast system server and limited integral calculation time.

The weather forecast method of the present application is described above, and the weather forecast system of the present application is described below, referring to fig. 6, including:

A first obtaining unit 601, configured to obtain a configuration file of weather forecast, where the configuration file includes an aging start point, an aging end point, and a time step of weather forecast data, where the time step separates a time period formed by connecting the aging start point and the aging end point to form a plurality of aging points, and the aging points include the aging start point and the aging end point;

a second obtaining unit 602, configured to obtain GFS initial field data corresponding to the aging points one by one according to the configuration file;

a preprocessing unit 603, configured to preprocess the GFS initial field data to obtain target GFS initial field data;

a third obtaining unit 604, configured to obtain weather forecast integral parameters;

the calculating unit 605 is configured to perform integral calculation according to the weather forecast integral parameter and the target GFS initial field data, so as to obtain weather forecast data;

an extracting unit 606 is configured to extract weather elements from the weather forecast data.

Optionally, when the second obtaining unit 602 obtains GFS initial field data corresponding to the aging point one by one according to the configuration file, the second obtaining unit is specifically configured to:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

Optionally, the preprocessing unit 603 performs preprocessing on the GFS initial field, and is specifically configured to:

obtaining topographic data of a weather forecast area;

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current GFS initial field;

taking the latter time effect point as a current time effect point;

Optionally, the calculating unit 605 performs integral calculation on the target GFS initial field data according to the weather forecast integral parameter, and is specifically configured to:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current target GFS initial field data;

Optionally, the extracting unit 606 extracts weather elements from weather forecast data includes:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

if so, performing post-processing operations such as spatial interpolation, physical conversion and the like on weather forecast data corresponding to the current aging point to obtain the required meteorological elements

Storing the meteorological elements to a database;

deleting weather forecast data corresponding to the current aging point;

In this embodiment, the operations performed by the weather forecast system are similar to those described in fig. 1 to 5, and are not described here again.

The weather forecast system of the present application is described below with reference to fig. 7, including:

the weather forecast system 700 may include one or more central processors (centralprocessing units, CPU) 701 and memory 705, with one or more applications or data stored in the memory 705. Wherein the memory 705 is persistent storage. The program stored in the memory 705 may include one or more modules, each of which may include a series of instruction operations in the weather forecast system. Still further, the central processor 701 may be configured to communicate with the memory 705 and execute a series of instruction operations in the memory 705 on the weather forecast system 700. The weather forecast system 700 may also include one or more power supplies 702, one or more wired or wireless network interfaces 703, one or more input/output interfaces 704, and an operating system linux, among others. The cpu 701 may perform the operations in the embodiments shown in fig. 1 to 5, and detailed descriptions thereof are omitted herein.

It will be clear to those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described systems, apparatuses and units may refer to corresponding procedures in the foregoing method embodiments, which are not repeated herein.

In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods may be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative, e.g., the division of the units is merely a logical function division, and there may be additional divisions when actually implemented, e.g., multiple units or components may be combined or integrated into another system, or some features may be omitted or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be an indirect coupling or communication connection via some interfaces, devices or units, which may be in electrical, mechanical or other form.

The units described as separate units may or may not be physically separate, and units shown as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.

In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may exist alone physically, or two or more units may be integrated in one unit. The integrated units may be implemented in hardware or in software functional units.

The integrated units, if implemented in the form of software functional units and sold or used as stand-alone products, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be embodied in essence or a part contributing to the prior art or all or part of the technical solution in the form of a software product stored in a storage medium, including several instructions to cause a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the methods described in the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a read-only memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk, or the like, which can store program codes.

Claims

1. A method of weather forecast, comprising:

acquiring weather forecast integral parameters;

extracting meteorological elements from the weather forecast data;

the step of obtaining the GFS initial field data corresponding to the aging points one by one according to the configuration file comprises the following steps:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

if the latter aging point is smaller than or equal to the aging end point, the latter aging point is regarded as a current aging point, and the step of acquiring URL links of GFS initial field data corresponding to the current aging point is triggered;

the step of carrying out integral calculation according to the weather forecast integral parameter and the target GFS initial field data to obtain weather forecast data comprises the following steps:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current target GFS initial field data;

2. The weather forecast method of claim 1, wherein the preprocessing the GFS initial field data to obtain target GFS initial field data includes:

Obtaining topographic data of a weather forecast area;

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current GFS initial field;

taking the latter aging point as the current aging point;

3. The method of claim 1, wherein the extracting weather elements from the weather forecast data comprises:

Determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

storing the meteorological elements to a database;

deleting weather forecast data corresponding to the current aging point;

4. A weather forecast method according to any one of claims 1 to 3, characterized in that the GFS initial field data is GFS initial field data of a certain aging point in a global coarse grid initial field made and issued by the national environment forecast center at daily timing.

5. A weather forecast system, comprising:

an extracting unit for extracting weather elements from the weather forecast data;

the second obtaining unit is specifically configured to:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

the computing unit is specifically configured to:

determining an aging starting point in the configuration file;

taking the aging starting point as a current aging point;

deleting the current target GFS initial field data;

6. A weather forecast system, comprising:

a processor, a memory, a bus, and an input/output device;

the processor is connected with the memory and the input and output equipment;

the processor performs the method of any one of claims 1 to 4.

7. A computer storage medium having instructions stored therein, which when executed on a computer, cause the computer to perform the method of any of claims 1 to 4.