CN111401632B - Weather forecast method and related device - Google Patents
Weather forecast method and related device Download PDFInfo
- Publication number
- CN111401632B CN111401632B CN202010174457.XA CN202010174457A CN111401632B CN 111401632 B CN111401632 B CN 111401632B CN 202010174457 A CN202010174457 A CN 202010174457A CN 111401632 B CN111401632 B CN 111401632B
- Authority
- CN
- China
- Prior art keywords
- aging
- point
- gfs
- current
- initial field
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000032683 aging Effects 0.000 claims abstract description 604
- 238000004364 calculation method Methods 0.000 claims abstract description 60
- 238000007781 pre-processing Methods 0.000 claims abstract description 26
- 230000010354 integration Effects 0.000 claims description 22
- 230000001960 triggered effect Effects 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 7
- 230000001550 time effect Effects 0.000 description 6
- 230000003203 everyday effect Effects 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000002354 daily effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/95—Retrieval from the web
- G06F16/955—Retrieval from the web using information identifiers, e.g. uniform resource locators [URL]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- General Business, Economics & Management (AREA)
- Databases & Information Systems (AREA)
- Marketing (AREA)
- Development Economics (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Entrepreneurship & Innovation (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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
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;
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;
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;
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;
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;
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 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;
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;
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;
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 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;
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;
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 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;
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;
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;
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 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;
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;
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;
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 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.
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.
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.
The execution of this step is similar to the execution of step 201 in fig. 2 and will not be described again here.
402. 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.
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.
The execution of this step is similar to the execution of step 208 of fig. 2, and will not be described in detail here.
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.
The execution of this step is similar to the execution of step 201 in fig. 2 and will not be described again here.
502. 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.
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.
The execution of this step is similar to the execution of step 208 of fig. 2, and will not be described in detail here.
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;
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 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;
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;
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 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;
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;
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;
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 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;
judging whether weather forecast data corresponding to the current aging point exists or not;
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;
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;
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 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.
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 (7)
1. A method of weather forecast, comprising:
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;
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;
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;
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;
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;
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;
triggering the step of adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.
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;
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;
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 aging point as the current aging 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.
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;
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;
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;
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.
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:
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;
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;
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;
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 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;
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;
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;
triggering the step of adding the time step to the current aging point to obtain the subsequent aging point of the current aging point.
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 bus is respectively connected with the processor, 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.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010174457.XA CN111401632B (en) | 2020-03-13 | 2020-03-13 | Weather forecast method and related device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010174457.XA CN111401632B (en) | 2020-03-13 | 2020-03-13 | Weather forecast method and related device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111401632A CN111401632A (en) | 2020-07-10 |
CN111401632B true CN111401632B (en) | 2023-07-21 |
Family
ID=71428747
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010174457.XA Active CN111401632B (en) | 2020-03-13 | 2020-03-13 | Weather forecast method and related device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111401632B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112565398B (en) * | 2020-12-03 | 2021-12-10 | 中科三清科技有限公司 | GFS data downloading method and device, electronic equipment and storage medium |
CN114090266B (en) * | 2021-12-01 | 2022-12-09 | 中科三清科技有限公司 | Air quality forecast generation method and device |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443833A (en) * | 2015-08-06 | 2017-02-22 | 中国电力科学研究院 | Numerical weather prediction method based on cloud computation |
CN110852523A (en) * | 2019-11-19 | 2020-02-28 | 上海眼控科技股份有限公司 | Weather forecasting method, device, equipment and storage medium based on numerical mode |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9952352B2 (en) * | 2015-09-25 | 2018-04-24 | The United States Of America As Represented By The Secretary Of The Army | Device and system for generating a weather forecast for a selected zone |
-
2020
- 2020-03-13 CN CN202010174457.XA patent/CN111401632B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106443833A (en) * | 2015-08-06 | 2017-02-22 | 中国电力科学研究院 | Numerical weather prediction method based on cloud computation |
CN110852523A (en) * | 2019-11-19 | 2020-02-28 | 上海眼控科技股份有限公司 | Weather forecasting method, device, equipment and storage medium based on numerical mode |
Non-Patent Citations (1)
Title |
---|
基于SMS的华中区域中尺度数值天气预报系统设计和实现;赖安伟等;《气象科技》(第01期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN111401632A (en) | 2020-07-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111414416B (en) | Data processing method, device, equipment and storage medium | |
CN111401632B (en) | Weather forecast method and related device | |
CN110633277A (en) | Time sequence data storage method and device, computer equipment and storage medium | |
US8370303B1 (en) | Generating snapshots of data tables | |
CN110569222B (en) | Link tracking method and device, computer equipment and readable storage medium | |
CN112910723B (en) | Edge terminal management method, device, equipment and storage medium | |
CN103064670A (en) | Method and system for innovation platform data management based on place net | |
CN112651693A (en) | Car system BOM structure data processing method for production material distribution | |
CN110852523A (en) | Weather forecasting method, device, equipment and storage medium based on numerical mode | |
CN111629001A (en) | Meteorological message generation method and device, computer equipment and storage medium | |
CN112651692A (en) | Production material distribution method based on vehicle BOM tree structure analysis | |
CN112488745A (en) | Intelligent charge control management method, device, equipment and storage medium | |
CN111105076A (en) | Weather forecast data generation method and device, computer equipment and storage medium | |
CN109766488B (en) | Data acquisition method based on Scapy | |
CN111158999A (en) | Method and device for calculating formatting time | |
CN110968929A (en) | Wind power plant wind speed prediction method and device and electronic equipment | |
CN113760983A (en) | Data processing method, system and non-transitory computer readable storage medium | |
CN116302455B (en) | Timing task processing method and device, storage medium and computer equipment | |
CN111221517A (en) | Model creating method and device, computer equipment and readable storage medium | |
CN111415200A (en) | Data processing method and device | |
CN116629015A (en) | Digital twin task processing method, device, equipment and storage medium | |
CN110222877A (en) | A kind of load prediction system and load forecasting method based on customized neural network | |
CN106570193B (en) | Method for loading time sequence big data | |
CN115907234A (en) | Method and device for predicting plant diseases and insect pests, electronic equipment and storage medium | |
CN111078748A (en) | Weather forecast data generation method and device, computer equipment and storage medium |
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 |