Multiband weather radar reflectivity factor networking method and device

Abstract

The invention discloses a multiband weather radar reflectivity factor networking method and device. The method comprises the following steps: acquiring S-band weather radar stream data with the time resolution of 1 minute and X-band weather radar base data with the time resolution of 3 minutes; reconstructing S-band weather radar base data; radar station parameters of 9 elevation angles of the S-band weather radar and radar station parameters of 9 elevation angles of the X-band weather radar are extracted respectively; generating equal-height surface reflectivity factor data of the S wave band according to the radar station parameters of the 9 elevation angles of the S wave band, and generating equal-height surface reflectivity factor data of the X wave band according to the radar station parameters of the 9 elevation angles of the X wave band; and carrying out fusion networking on the S-band equal-height surface reflectivity factor data and the X-band equal-height surface reflectivity factor data. According to the method, the fusion networking time resolution can be improved, and meanwhile, the equal-height surface reflectivity factors of the multi-band radar can be quickly checked in real time.

Description

Technical Field

[0001] The present invention belongs to the field of meteorological technology, and more specifically, relates to a multi-band weather radar reflectivity factor networking method and device. Background Art

[0002] S-band operational weather radar can quantitatively detect the intensity, radial velocity, spectral width and other information of echoes. Its high temporal and spatial resolution and accurate detection capability have become an extremely effective monitoring and early warning tool for mesoscale disastrous weather. However, the detection range of a single ground-fixed radar is limited and insufficient to cover complex and changing weather systems. In addition, my country's terrain conditions are very complex, and most operational radars are generally blocked to varying degrees at an elevation angle of 0.5°, which limits the S-band operational weather radar network's ability to detect low-level meteorological targets.

[0003] The X-band dual-polarization Doppler weather radar has a higher resolution than the S-band operational weather radar, and can fill in the blind spots of the S-band operational weather radar's low-altitude detection and static vertebral blind spots. Compared with the S-band operational weather radar, the X-band weather radar is relatively small in size, and the site construction and radar installation are more portable and low in cost. More radars can be built with a small amount of money, which can provide detailed observation data for mountainous areas and key monitoring areas, and enhance the local disaster weather monitoring and early warning capabilities. Therefore, in order to improve the low-altitude detection capabilities, the joint network observation of the X-band radar and the S-band operational weather radar can better reveal the laws of atmospheric movement during the weather process and effectively enhance the detection capabilities of disaster weather.

[0004] At present, domestic X-band dual-polarization Doppler weather radar and S-band operational weather radar are usually operated separately. The two have their own business systems and are not integrated into a whole. With the construction of the X-band radar network, it will become an inevitable trend to integrate it with the current S-band operational weather radar network. Existing research has proposed a fusion method for S-band weather radar and X-band weather radar. This method can fuse the radar reflectivity factors of the overlapping areas of radars of different bands, but it is based on the S-band weather radar volume scanning data with a time resolution of 6 minutes, which cannot meet the timeliness requirements of actual applications. Summary of the invention

[0005] In view of the above defects or improvement needs of the prior art, the present invention provides a multi-band weather radar reflectivity factor networking method and device, which can improve the fusion networking resolution and can quickly view the multi-band radar reflectivity factor in real time.

[0006] To achieve the above object, according to one aspect of the present invention, a multi-band weather radar reflectivity factor networking method is provided, comprising the steps of:

[0007] S1, obtains S-band weather radar streaming data with a time resolution of 1 minute and X-band weather radar base data with a time resolution of 3 minutes;

[0008] S2, matching the stream data closest to the current time in the S-band weather radar stream data, and taking the stream data closest to the current time as the end stream data, matching N stream data forward, N is a preset number, and constructing the S-band weather radar base data with a time resolution of 1 minute according to the standard format of weather radar base data from the matched N+1 stream data;

[0009] S3, parse the reconstructed S-band weather radar base data and X-band weather radar base data, extract radar site parameters of 9 elevation angles of S-band weather radar and radar site parameters of 9 elevation angles of X-band weather radar, respectively. Each set of radar site parameters includes radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor;

[0010] S4, generates the S-band contour surface reflectivity factor data according to the radar site parameters of the 9 elevation angles of the S-band, and generates the X-band contour surface reflectivity factor data according to the radar site parameters of the 9 elevation angles of the X-band;

[0011] S5, the S-band contour surface reflectivity factor data and the X-band contour surface reflectivity factor data are fused and networked.

[0012] Furthermore, the step of constructing S-band weather radar base data with a time resolution of 1 minute from the matched N+1 stream data in a standard format of weather radar base data includes the following steps:

[0013] Get the matching N+1 stream data;

[0014] The common data blocks of the N+1 stream data are parsed, and the common data blocks include common header data, site configuration data, task configuration data, and scan configuration data.

[0015] Parsing radial data blocks of N+1 stream data, where the radial data blocks include radial headers and radial data;

[0016] Generate base data based on virtual volume scanning technology, and write the common data blocks of N+1 stream data into the base data of virtual volume scanning technology in binary format;

[0017] Rewrite the scanning start time in the task configuration module in the base data based on the virtual body scanning technology, and the start time is based on the earliest time in the combination of stream data based on the virtual body scanning technology;

[0018] After sorting the elevation angle numbers in the radial headers of the radial data blocks of the N+1 stream data in ascending order, the radial header blocks and radial data blocks corresponding to each elevation angle are written into the base data based on the virtual volume scanning technology.

[0019] Furthermore, the S3 comprises the steps of:

[0020] The basic data of S-band weather radar were parsed, and the radar site parameters of S-band weather radar were extracted according to 11 layers of elevation angles, 360 azimuth angles, and 920 distance libraries. Finally, the radar site parameters of 9 elevation angles of 0.5°, 1.5°, 2.4°, 3.4°, 4.3°, 6.6°, 9.9°, 14.6°, and 19.5° were retained.

[0021] The X-band weather radar basic data were parsed, and the radar site parameters of the X-band weather radar were extracted according to 11 layers of elevation angles, 360 azimuth angles, and 1000 distance libraries. Finally, the radar site parameters of 9 elevation angles of 0.5°, 1.5°, 2.4°, 3.4°, 4.3°, 6.6°, 9.9°, 14.6°, and 19.5° were retained.

[0022] Furthermore, the S4 comprises the steps of:

[0023] Read the radar site parameters of 9 elevation angles of the S band and the radar site parameters of 9 elevation angles of the X band, add or subtract 2.3° to the longitude and latitude values ​​of the S band radar site as the longitude and latitude grid range of the radar reflectivity factor, add or subtract 0.7° to the longitude and latitude values ​​of the X band radar site as the longitude and latitude grid range of the radar reflectivity factor, and set the grid resolution to 0.0025°×0.0025°. Set 17 layers of height in the vertical direction, with a height range of 0-12km, and the vertical resolution below 5km. The vertical resolution is 0.5km, and the vertical resolution of 5-12km is 1km. The height grids are 0.5km, 1.0km, 1.5km, 2.0km, 2.5km, 3.0km, 3.5km, 4.0km, 4.5km, 5.0km, 6.0km, 7.0km, 8.0km, 9.0km, 10.0km, 11.0km, and 12.0km, forming a three-dimensional grid of the S band and a three-dimensional grid of the X band composed of longitude, latitude, and altitude;

[0024] Taking the three-dimensional grid as the Cartesian coordinate, the elevation angle, azimuth angle, and range spherical coordinates of the radar station corresponding to the S band of each three-dimensional grid are inversely calculated, which are recorded as (r, a, e), where r represents the range spherical coordinate corresponding to the S band, a represents the azimuth angle corresponding to the S band, and e represents the elevation angle corresponding to the S band; the elevation angle, azimuth angle, and range spherical coordinates of the radar station corresponding to the X band of each three-dimensional grid are inversely calculated, which are recorded as (r′, a′, e′), where r′ represents the range spherical coordinate corresponding to the X band, a′ represents the azimuth angle corresponding to the X band, and e′ represents the elevation angle corresponding to the X band;

[0025] The reflectivity factor of each three-dimensional grid in the S band is calculated to generate the reflectivity factor data of the contour surface in the S band. The reflectivity factor of each three-dimensional grid in the X band is calculated to generate the reflectivity factor data of the contour surface in the X band.

[0026] Furthermore, the calculation formula for inversely calculating the elevation angle, azimuth angle, and range spherical coordinates of each three-dimensional grid corresponding to the S-band radar site is:

[0027]

[0028] where lon c ,lat c is the longitude and latitude of the S-band radar site, lon, lat are the longitude and latitude of any grid cell in the S-band Cartesian three-dimensional grid, and A is calculated as follows:

[0029]

[0030] The calculation formula of S is:

[0031] S=cos -1 (sin(lat c )×sin(lat)+cos(lat c )×cos(lat)×cos(lon-lon c ))

[0032]

[0033] where h c is the height of the S-band radar, and h is the height of any grid cell in the S-band Cartesian three-dimensional grid;

[0034] The slope distance r is calculated as follows:

[0035]

[0036] The calculation formula for inversely calculating the elevation angle, azimuth angle, and range spherical coordinates of each three-dimensional grid corresponding to the X-band radar site is:

[0037]

[0038] where lon′ c , lat′ c is the longitude and latitude of the X-band radar site, lon′, lat′ are the longitude and latitude of any grid cell in the X-band Cartesian three-dimensional grid, and A′ is calculated as follows:

[0039]

[0040] The calculation formula of S′ is:

[0041] S′=cos -1 (sin(lat′ c )×sin(lat′)+cos(lat′ c )×cos(lat′)×cos(lon′-lon′ c ))

[0042]

[0043] where h′ c is the height of the radar, h′ is the height of any grid cell in the Cartesian three-dimensional grid;

[0044] The slope distance r′ is calculated as follows:

[0045]

[0046] Furthermore, the calculation formula for calculating the reflectivity factor of each three-dimensional grid in the S band is:

[0047]

[0048] Among them, f(r,a,e) is the reflectivity factor of the coordinate (r,a,e), e1 and e2 are the upper and lower elevation angles closest to the elevation angle e, f(r,a,e1) and f(r,a,e2) are the reflectivity factor values ​​corresponding to the upper and lower elevation angles closest to the azimuth a and distance r, respectively;

[0049] The calculation formula for calculating the reflectivity factor of each three-dimensional grid in the X band is:

[0050]

[0051] Wherein, f(r′,a′,e′) is the reflectivity factor of the coordinate (r′,a′,e′), e′1 and e′2 are the upper and lower elevation angles closest to the elevation angle e′, and f(r′,a′,e′1) and f(r′,a′,e′2) are the reflectivity factor values ​​corresponding to the upper and lower elevation angles closest to the azimuth a′ and distance r′, respectively.

[0052] Furthermore, the S5 comprises the steps of:

[0053] For the area covered by S-band weather radar detection and the area not covered by X-band weather radar detection, the S-band weather radar detection reflectivity factor is used as the reflectivity factor value after networking;

[0054] For the area not covered by S-band weather radar detection and the area covered by X-band weather radar detection, the reflectivity factor of X-band weather radar detection is used as the reflectivity factor value after networking;

[0055] For areas not covered by both S-band weather radar and X-band weather radar detection, there is no reflectivity factor value after networking;

[0056] For areas covered by both S-band weather radar and X-band weather radar, the larger reflectivity factor of the two is used as the reflectivity factor value after networking.

[0057] Furthermore, N=10.

[0058] According to another aspect of the present invention, the present invention provides a multi-band weather radar reflectivity factor networking device, comprising:

[0059] A data acquisition module is used to acquire S-band weather radar stream data with a time resolution of 1 minute and X-band weather radar base data with a time resolution of 3 minutes;

[0060] The base data construction module is used to match the stream data closest to the current time in the S-band weather radar stream data, and use the stream data closest to the current time as the end stream data, match N stream data forward, N is a preset number, and construct the S-band weather radar base data with a time resolution of 1 minute according to the standard format of weather radar base data from the matched N+1 stream data;

[0061] The radar site parameter extraction module is used to parse the reconstructed S-band weather radar base data and X-band weather radar base data, and extract the radar site parameters of the S-band weather radar at 9 elevation angles and the radar site parameters of the X-band weather radar at 9 elevation angles, respectively. Each set of radar site parameters includes the radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor.

[0062] A contour surface reflectivity factor data generation module is used to generate the contour surface reflectivity factor data of the S band according to the radar site parameters of the 9 elevation angles of the S band, and to generate the contour surface reflectivity factor data of the X band according to the radar site parameters of the 9 elevation angles of the X band;

[0063] The fusion networking module is used to fuse and network the S-band contour surface reflectivity factor data and the X-band contour surface reflectivity factor data.

[0064] Further, the radar site parameter extraction module includes radar site parameter extraction submodules distributed on each radar single station, and the radar site parameter extraction submodule on each radar single station is used to extract radar site parameters of 9 elevation angles in sequence, and after extracting the radar site parameters of a single elevation angle, while extracting the radar site parameters of the next elevation angle, the radar site parameters of the single elevation angle are saved;

[0065] The contour surface reflectivity factor data generation module includes a contour surface reflectivity factor data generation submodule distributed on each radar station. The contour surface reflectivity factor data generation submodule on each radar station is used to generate the contour surface reflectivity factor data of each altitude in sequence, and after extracting the contour surface reflectivity factor data of a single altitude, extract the contour surface reflectivity factor data of the next altitude, and save the contour surface reflectivity factor data of the single altitude at the same time.

[0066] In general, the above technical solutions conceived by the present invention have beneficial effects compared with the prior art:

[0067] (1) The original reflectivity factor fusion product with a time resolution of 6 minutes can be upgraded to a time resolution of 1 minute.

[0068] (2) It can automatically and quickly match the base data involved in the fusion and automatically generate reflectivity factor fusion products for users to quickly access in real time. BRIEF DESCRIPTION OF THE DRAWINGS

[0069] Figure 1 is a flow chart of a multi-band weather radar reflectivity factor networking method according to an embodiment of the present invention;

[0070] Figure 2 It is a schematic diagram of radar site parameter extraction of a multi-band weather radar reflectivity factor networking device according to an embodiment of the present invention;

[0071] Figure 3 It is a schematic diagram of generating contour surface reflectivity factor data of a multi-band weather radar reflectivity factor networking device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0072] In order to make the purpose, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.

[0073] In the description of the embodiments of the present application, the terms "first" and "second" are used for descriptive purposes only and are not to be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one of the features. "Multiple" means at least two, such as two, three, etc., unless otherwise clearly and specifically defined.

[0074] Unless otherwise specified, "plurality" means two or more.

[0075] The terms "including" and "having" and any variations thereof in the embodiments of the present invention are intended to cover non-exclusive inclusions. For example, a process, method, apparatus, product or equipment comprising a series of steps or modules is not necessarily limited to those steps or modules clearly listed, but may include other steps or modules not clearly listed or inherent to these processes, methods, products or equipment.

[0076] The naming or numbering of the steps in the embodiments of the present invention does not mean that the steps in the method flow must be executed in the time / logical sequence indicated by the naming or numbering. The execution order of the named or numbered process steps can be changed according to the technical purpose to be achieved, as long as the same or similar technical effects can be achieved.

[0077] Reference to "embodiments" herein means that a particular feature, structure, or characteristic described in conjunction with the embodiments may be included in at least one embodiment of the present invention. The appearance of the phrase in various places in the specification does not necessarily refer to the same embodiment, nor is it an independent or alternative embodiment that is mutually exclusive with other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0078] The embodiments of the present invention provide a multi-band weather radar reflectivity factor networking method and system, which are described below respectively.

[0079] like Figure 1 As shown, an embodiment of the present invention provides a multi-band weather radar reflectivity factor networking method, comprising the steps of:

[0080] S1, obtains S-band weather radar streaming data with a time resolution of 1 minute and X-band weather radar base data with a time resolution of 3 minutes.

[0081] S2, matches the stream data closest to the current time in the S-band weather radar stream data, and takes the stream data closest to the current time as the end stream data, matches N stream data forward, N is a preset number, and constructs the S-band weather radar base data with the matched N+1 stream data according to the standard format of weather radar base data.

[0082] The standard format of weather radar basic data is a standard data format developed to solve the problem of diversity and inconsistency of meteorological data formats.

[0083] Further, S2 comprises sub-steps:

[0084] Step 2.1 matches the stream data closest to the current time in the S-band weather radar stream data.

[0085] Specifically, monitor the radar station stream data folder, read the current time, sort the files according to the file modification date, and match the stream data closest to the current time. The stream data file name naming rules are as follows:

[0086] Z_RADA_I_{IIii}i_{YYYYMMDDhhmmss}_O_DOR-CUT_SAD_CAP_{elevation angle}_{serial number}_FMT.bin.bz2. The {} represent the variables in the naming rules, and the definitions are as follows:

[0087] Z_RADAR_I: unified prefix for radar stream data, Z is a fixed code indicating business data; RADA is a fixed code indicating radar data; I is a fixed code indicating that the following IIiii field is the zone station number of the radar station;

[0088] IIiii: indicates the radar zone station number, such as Z9270;

[0089] YYYYMMDDhhmmss: Universal time, scan start time, in the format of “year month day hour minute second”;

[0090] O: fixed code, indicating basic data;

[0091] DOR-CUT: fixed code, indicating Doppler radar stream data;

[0092] SAD: radar model, SAD means SA dual polarization radar;

[0093] Elevation: indicates the elevation angle, including 5, 15, 25, 34, 43, 60, 99, 146, 195;

[0094] Serial number: indicates the elevation layer number, including 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11;

[0095] FMT: Fixed code, indicating data in a standard format;

[0096] bin: indicates that the file format is binary data;

[0097] bz2: An optional field indicating a file compressed using bzip2.

[0098] Step 2.2 takes the stream data closest to the current time as the end stream data, and matches N stream data forward, where N is a preset number.

[0099] In one embodiment, N is 10, that is, the latest stream data is used as the end stream data, 10 stream data are matched forward, and these 11 stream data are combined into a stream data combination based on the virtual body scanning technology;

[0100] For example, if the current time is 11:00 on July 1, 2024, the file name closest to the current time in the radar station stream data folder is:

[0101] "Z_RADR_I_Z9710_20240626030018_O_DOR-CUT_SAD_CAP_5_1_FMT.bin.bz2"

[0102] The time in the file name is Universal Time.

[0103] Then match 10 stream data, the file names are:

[0104] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_5_2_FMT.bin.bz2"

[0105] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_15_3_FMT.bin.bz2"

[0106] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_24_5_FMT.bin.bz2"

[0107] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_34_6_FMT.bin.bz2"

[0108] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_43_7_FMT.bin.bz2"

[0109] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_60_8_FMT.bin.bz2"

[0110] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_99_9_FMT.bin.bz2"

[0111] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_146_10_FMT.bin.bz2"

[0112] "Z_RADR_I_Z9710_20240626025441_O_DOR-CUT_SAD_CAP_195_11_FMT.bin.bz2"

[0113] Step 2.3 constructs S-band weather radar base data from the matched N+1 stream data according to the standard format of weather radar base data.

[0114] The data structure of the base data is shown in Table 1.

[0115] Table 1 Weather radar binary base data structure

[0116]

[0117]

[0118] The common data blocks of the stream data and the base data are the same, except that the stream data has only one radial data block, while the base data has M radial data blocks.

[0119] Further, taking N=10 as an example, S203 includes sub-steps:

[0120] (1) Batch read the above 11 binary weather radar stream data.

[0121] (2) Parse the common data blocks in the stream data, which include a common header (32 bytes); site configuration (128 bytes); task configuration (256) bytes; and scan configuration (P*256 bytes).

[0122] (3) Parse the radial data blocks in the stream data. Each radial data block includes a radial header (64 bytes), K radial data headers (32×K bytes), and K radial data (1×K bytes).

[0123] (4) Generate the initialized virtual volume scanning technology base data with standard format, and write the common data blocks in the stream data read in step (2) including the common header (32 bytes); site configuration (128 bytes); task configuration (256) bytes; scanning configuration (P*256 bytes) into the base data based on virtual volume scanning technology in binary format. The common data blocks of the 11 stream data are the same, and any one of them can be selected and written in.

[0124] (5) Rewrite the scan start time in the task configuration module in the base data based on the virtual body scan technology in step (4), and the start time is based on the earliest time in the stream data combination based on the virtual body scan technology. In other words, replace the "scan start time in the initialization base data" with "the earliest time in the stream data combination of the virtual body scan technology".

[0125] For the 10 streams matched in the above example, except for the first one "20240626030018", the time of the following 10 data are all "20240626025441". "20240626030018" is later than "20240626025441", so the earliest time is "20240626025441", and the earliest time "20240626025441" is used to replace "the scan start time in the initialization base data".

[0126] (6) After sorting the elevation angle numbers in the radial header blocks of the 11 stream data radial data blocks read in step (3) in ascending order, the radial header blocks and radial data blocks corresponding to each elevation angle are written into the base data based on the virtual volume scanning technology in step (4).

[0127] Each time the current time is updated, steps 2.1-2.3 need to be executed again. The time resolution is 1 minute, that is, every one minute, steps 2.1-2.3 need to be executed again to construct the S-band weather radar basic data with a time resolution of 1 minute.

[0128] S3, parse the S-band weather radar base data and the X-band weather radar base data, extract the radar site parameters of the 9 elevation angles of the S-band weather radar and the radar site parameters of the 9 elevation angles of the X-band weather radar, respectively. Each set of radar site parameters includes the radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor.

[0129] The specific method for analyzing the S-band weather radar basic data and the X-band weather radar basic data based on the virtual volume scanning technology includes the following steps:

[0130] Step 3.1 monitors the S-band weather radar base data folder and the X-band weather radar base data folder based on the virtual volume scanning technology.

[0131] Step 3.2 Sort the base data in the S / X band monitoring folder by update time, select the latest base data for analysis, and write the analyzed file into the log. This is because the program runs in real time, so each update of data requires analysis of the new base data. At the same time, the processed base data should be written into the log to avoid repeated processing of the processed base data during operation, which would take up time.

[0132] Step 3.3 In the S / X band monitoring folder, compare the previous ten basic data with the analysis records in the log. If the previous ten basic data have been analyzed, no more analysis is required. If the previous ten basic data contain unanalyzed data, the unanalyzed basic data still needs to be analyzed. This is because sometimes various reasons (such as network problems) will cause data upload delays, and a lot of data will not be uploaded in the middle. When a large amount of data is uploaded at a certain time, it is impossible to process all the data, which will take up time. Therefore, compare the log file and select the 10 most recently uploaded basic data for processing.

[0133] Step 3.4 follows the following rules for base data analysis: Use Python programming language and its third-party libraries Pycwr and Pycinrad to extract the values ​​corresponding to the radar site longitude, radar site latitude, elevation, azimuth, distance library, and reflectivity factor bytes in the base data. However, since the azimuth and distance library sizes of each layer of elevation are not fixed, they need to be regularized. After regularization, the elevation, azimuth, and distance library sizes of the S-band radar are all fixed arrays of 11*360*920; after regularization, the elevation, azimuth, and distance library sizes of the X-band radar are all fixed arrays of 11*360*1000.

[0134] In step 3.5, due to the radar volume scanning strategy, only the reflectivity factor data of the 9 elevation angles in step 2.4 are retained and entered into CAPPI processing. The 9 elevation angles are 0.5°, 1.5°, 2.4°, 3.4°, 4.3°, 6.6°, 9.9°, 14.6°, and 19.5°, and the corresponding elevation angle numbers are 1, 3, 5, 6, 7, 8, 9, 10, and 11. The radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor of the 9 elevation angles are stored in the specified folder in NC file format. NC file refers to Network Common Data Format, also known as NetCDF file, which is the network common data format, mainly used for the storage of meteorological data.

[0135] S4, generates the S-band contour surface reflectivity factor data according to the radar site parameters of the 9 elevation angles of the S-band, and generates the X-band contour surface reflectivity factor data according to the radar site parameters of the 9 elevation angles of the X-band.

[0136] Step 4.1 Monitor 9 elevation reflectivity factor NC folders, read 9 elevation reflectivity factor NC files, extract radar site longitude, radar site latitude, azimuth, distance library, reflectivity factor, add or subtract 2.3° from the longitude and latitude values ​​of the S-band radar site as the longitude and latitude grid range of the radar reflectivity factor, add or subtract 0.7° from the longitude and latitude values ​​of the X-band radar site as the longitude and latitude grid range of the radar reflectivity factor, and set the grid resolution to 0.0025°×0.0025°, approximately 250m×250m. A total of 17 layers of height are set in the vertical direction, with a height range of 0-12km. The vertical resolution below 5km is 0.5km, and the vertical resolution between 5-12km is 1km, that is, the height grids are (0.5km, 1.0km, 1.5km, 2.0km, 2.5km, 3.0km, 3.5km, 4.0km, 4.5km, 5.0km, 6.0km, 7.0km, 8.0km, 9.0km, 10.0km, 11.0km, 12.0km), forming a three-dimensional grid of the S band and a three-dimensional grid of the X band composed of longitude, latitude and height;

[0137] Step 3.2: Using the three-dimensional grid as the Cartesian coordinate, inversely calculate the elevation angle, azimuth angle, and range spherical coordinates of the radar site corresponding to the S band for each three-dimensional grid, denoted as (r, a, e), where r represents the range spherical coordinate corresponding to the S band, a represents the azimuth angle corresponding to the S band, and e represents the elevation angle corresponding to the S band; inversely calculate the elevation angle, azimuth angle, and range spherical coordinates of the radar site corresponding to the X band for each three-dimensional grid, denoted as (r′, a′, e′), where r′ represents the range spherical coordinate corresponding to the X band, a′ represents the azimuth angle corresponding to the X band, and e′ represents the elevation angle corresponding to the X band.

[0138] The azimuth angle a is calculated as follows:

[0139]

[0140] where lon c ,lat c is the latitude and longitude of the S-band radar site, lon, lat are the latitude and longitude of any grid cell in the Cartesian three-dimensional grid, and A is calculated as follows:

[0141]

[0142] The calculation method of S is as follows:

[0143] S=cos -1 (sin(lat c )×sin(lat)+cos(lat c )×cos(lat)×cos(lon-lon c ))

[0144] The elevation angle e is calculated as follows:

[0145]

[0146] where h c is the height of the S-band radar, and h is the height of any grid cell in the Cartesian three-dimensional grid.

[0147] The slope distance r is calculated as follows:

[0148]

[0149] The azimuth angle a′ is calculated as follows:

[0150]

[0151] where lon′ c , lat′ c is the longitude and latitude of the X-band radar site, lon′, lat′ are the longitude and latitude of any grid cell in the Cartesian three-dimensional grid, and A′ is calculated as follows:

[0152]

[0153] The calculation formula of S′ is:

[0154] S′=cos -1 (sin(lat′ c )×sin(lat′)+cos(lat′ c )×cos(lat′)×cos(lon′-lon′ c ))

[0155]

[0156] where h′ c is the height of the X-band radar, and h′ is the height of any grid cell in the Cartesian three-dimensional grid.

[0157] The slope distance r′ is calculated as follows:

[0158]

[0159] Step 3.3 Calculate the reflectivity factor of each three-dimensional grid in the S band, generate the contour reflectivity factor data of the S band, calculate the reflectivity factor of each three-dimensional grid in the X band, and generate the contour reflectivity factor data of the X band:

[0160]

[0161] Where e is the elevation angle on the Cartesian coordinate grid point, e1 and e2 are the upper and lower elevation angles closest to the elevation angle e, and f(r,a,e1) and f(r,a,e2) are the reflectivity factor values ​​corresponding to the upper and lower elevation angles closest to the azimuth a and distance r.

[0162] The calculation formula for calculating the reflectivity factor of each three-dimensional grid in the X band is:

[0163]

[0164] Wherein, f(r′,a′,e′) is the reflectivity factor of the coordinate (r′,a′,e′), e′1 and e′2 are the upper and lower elevation angles closest to the elevation angle e′, and f(r′,a′,e′1) and f(r′,a′,e′2) are the reflectivity factor values ​​corresponding to the upper and lower elevation angles closest to the azimuth a′ and distance r′, respectively.

[0165] According to the above method, CAPPI products are generated for each radar in real time, and the longitude, latitude, altitude and reflectivity factor are written into the NC file in the specified folder for subsequent puzzle calls.

[0166] S5, the S-band contour surface reflectivity factor data and the X-band contour surface reflectivity factor data are fused and networked.

[0167] Step 5.1 monitor the CAPPI product NC folder of each radar, and select the CAPPI NC file closest to the current time every 1 minute and input it into the jigsaw algorithm;

[0168] Step 5.2: Fusion of all radar and other high reflectivity factor products. The fusion follows the following four principles:

[0169] For the area covered by S-band weather radar detection and the area not covered by X-band weather radar detection, the S-band weather radar detection reflectivity factor is used as the reflectivity factor value after networking;

[0170] For the area not covered by S-band weather radar detection and the area covered by X-band weather radar detection, the reflectivity factor of X-band weather radar detection is used as the reflectivity factor value after networking;

[0171] For areas not covered by both S-band weather radar and X-band weather radar detection, there is no reflectivity factor value after networking;

[0172] For areas covered by both S-band weather radar and X-band weather radar, the maximum value method is used to determine the network reflectivity factor value, that is, the larger reflectivity factor of the two is used as the reflectivity factor value after networking.

[0173] A multi-band weather radar reflectivity factor networking device according to an embodiment of the present invention includes:

[0174] A data acquisition module is used to acquire S-band weather radar stream data with a time resolution of 1 minute and X-band weather radar base data with a time resolution of 3 minutes;

[0175] The base data construction module is used to match the stream data closest to the current time in the S-band weather radar stream data, and use the stream data closest to the current time as the end stream data, match N stream data forward, N is a preset number, and construct the S-band weather radar base data with a time resolution of 1 minute according to the standard format of weather radar base data from the matched N+1 stream data;

[0176] The radar site parameter extraction module is used to parse the S-band weather radar base data and the X-band weather radar base data, and respectively extract the radar site parameters of the 9 elevation angles of the S-band weather radar and the radar site parameters of the 9 elevation angles of the X-band weather radar. Each set of radar site parameters includes the radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor.

[0177] A contour surface reflectivity factor data generation module is used to generate the contour surface reflectivity factor data of the S band according to the radar site parameters of the 9 elevation angles of the S band, and to generate the contour surface reflectivity factor data of the X band according to the radar site parameters of the 9 elevation angles of the X band;

[0178] The fusion networking module is used to fuse and network the S-band contour surface reflectivity factor data and the X-band contour surface reflectivity factor data.

[0179] Further, the radar site parameter extraction module includes radar site parameter extraction submodules distributed on each radar single station, and the radar site parameter extraction submodule on each radar single station is used to extract radar site parameters of 9 elevation angles in sequence, and after extracting the radar site parameters of a single elevation angle, while extracting the radar site parameters of the next elevation angle, the radar site parameters of the single elevation angle are saved;

[0180] The contour surface reflectivity factor data generation module includes a contour surface reflectivity factor data generation submodule distributed on each radar station. The contour surface reflectivity factor data generation submodule on each radar station is used to generate the contour surface reflectivity factor data of each altitude in sequence, and after extracting the contour surface reflectivity factor data of a single altitude, extract the contour surface reflectivity factor data of the next altitude, and save the contour surface reflectivity factor data of the single altitude at the same time.

[0181] The implementation principle and technical effect of the multi-band weather radar reflectivity factor networking device are the same as the above-mentioned multi-band weather radar reflectivity factor networking method, which will not be repeated here.

[0182] Specifically, the working diagram of the multi-band weather radar reflectivity factor networking device is as follows: Figure 2 and Figure 3 shown.

[0183] When generating a single-station radar CAPPI product, the data of a single radar station must be executed on one CPU core (within one process). Therefore, the smallest parallel execution unit designed is the single radar station. One radar station data starts a process and occupies one CPU core.

[0184] like Figure 2 As shown, there are multiple elevation angles in the internal data of a single radar station. The data between the elevation angles are logically related, but the data of each elevation angle can be calculated logically independently. The generated regularized single-station elevation reflectivity factor data must be written as a file in NC format. The reading and writing of files is a slow IO process, so the processing of the single-station elevation reflectivity factor data and the saving of its NC file are designed to be executed concurrently as a coroutine to make full use of CPU resources to speed up the processing process.

[0185] like Figure 3 As shown in the figure, after the single-station radar data is parsed, a coroutine is started after each elevation data calculation is completed to save the NC file of the elevation data. IO operation is a time-consuming slow waiting operation. During this waiting process for IO, the CPU can continue to process the next elevation data. By analogy, the total processing time of a single station data is almost the same as the saving time of all elevation NC files, and the time taken for data processing can be almost ignored. Based on the above conclusions, multiple radar station data are designed to be executed in parallel on multiple CPUs using a multi-process method. Then the overall time consumption is the parallel saving time of all elevation data NC files of all radar stations (the bottleneck is in the storage system), and the time consumption of the calculation process is reduced to a minimum.

[0186] It will be easily understood by those skilled in the art that the above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions and improvements made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims

1. A multi-band weather radar reflectivity factor networking method, characterized in that: Includes steps: S1, obtains S-band weather radar streaming data with a time resolution of 1 minute and X-band weather radar base data with a time resolution of 3 minutes; S2, matching the stream data closest to the current time in the S-band weather radar stream data, and taking the stream data closest to the current time as the end stream data, matching N stream data forward, N is a preset number, and constructing the S-band weather radar base data with a time resolution of 1 minute according to the standard format of weather radar base data from the matched N+1 stream data; S3, parse the reconstructed S-band weather radar base data and X-band weather radar base data, extract radar site parameters of 9 elevation angles of S-band weather radar and radar site parameters of 9 elevation angles of X-band weather radar, respectively. Each set of radar site parameters includes radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor; S4, generates the S-band contour surface reflectivity factor data according to the radar site parameters of the 9 elevation angles of the S-band, and generates the X-band contour surface reflectivity factor data according to the radar site parameters of the 9 elevation angles of the X-band; S5, the S-band contour surface reflectivity factor data and the X-band contour surface reflectivity factor data are fused and networked.

2. A multi-band weather radar reflectivity factor networking method as claimed in claim 1, characterized in that: The step of constructing S-band weather radar base data with a time resolution of 1 minute according to the standard format of weather radar base data by using the matched N+1 stream data includes the following steps: Get the matching N+1 stream data; The common data blocks of the N+1 stream data are parsed, and the common data blocks include common header data, site configuration data, task configuration data, and scan configuration data. Parsing radial data blocks of N+1 stream data, where the radial data blocks include radial headers and radial data; Generate base data based on virtual volume scanning technology, and write the common data blocks of N+1 stream data into the base data of virtual volume scanning technology in binary format; Rewrite the scanning start time in the task configuration module in the base data based on the virtual body scanning technology, and the start time is based on the earliest time in the combination of stream data based on the virtual body scanning technology; After sorting the elevation angle numbers in the radial headers of the radial data blocks of the N+1 stream data in ascending order, the radial header blocks and radial data blocks corresponding to each elevation angle are written into the base data based on the virtual volume scanning technology.

3. A multi-band weather radar reflectivity factor networking method as claimed in claim 1, characterized in that: The S3 comprises the steps of: The basic data of S-band weather radar were parsed, and the radar site parameters of S-band weather radar were extracted according to 11 layers of elevation angles, 360 azimuth angles, and 920 distance libraries. Finally, the radar site parameters of 9 elevation angles of 0.5°, 1.5°, 2.4°, 3.4°, 4.3°, 6.6°, 9.9°, 14.6°, and 19.5° were retained. The X-band weather radar basic data were parsed, and the radar site parameters of the X-band weather radar were extracted according to 11 layers of elevation angles, 360 azimuth angles, and 1000 distance libraries. Finally, the radar site parameters of 9 elevation angles of 0.5°, 1.5°, 2.4°, 3.4°, 4.3°, 6.6°, 9.9°, 14.6°, and 19.5° were retained.

4. A multi-band weather radar reflectivity factor networking method as claimed in claim 1, characterized in that: The S4 comprises the steps of: Read the radar site parameters of 9 elevation angles of the S band and the radar site parameters of 9 elevation angles of the X band, add or subtract 2.3° to the longitude and latitude values ​​of the S band radar site as the longitude and latitude grid range of the radar reflectivity factor, add or subtract 0.7° to the longitude and latitude values ​​of the X band radar site as the longitude and latitude grid range of the radar reflectivity factor, and set the grid resolution to 0.0025°×0.0025°. Set 17 layers of height in the vertical direction, with a height range of 0-12km, and the vertical resolution below 5km. The vertical resolution is 0.5km, and the vertical resolution of 5-12km is 1km. The height grids are 0.5km, 1.0km, 1.5km, 2.0km, 2.5km, 3.0km, 3.5km, 4.0km, 4.5km, 5.0km, 6.0km, 7.0km, 8.0km, 9.0km, 10.0km, 11.0km, and 12.0km, forming a three-dimensional grid of the S band and a three-dimensional grid of the X band composed of longitude, latitude, and altitude; Taking the three-dimensional grid as the Cartesian coordinate, the elevation angle, azimuth angle, and range spherical coordinates of the radar station corresponding to the S band of each three-dimensional grid are inversely calculated, which are recorded as (r, a, e), where r represents the range spherical coordinate corresponding to the S band, a represents the azimuth angle corresponding to the S band, and e represents the elevation angle corresponding to the S band; the elevation angle, azimuth angle, and range spherical coordinates of the radar station corresponding to the X band of each three-dimensional grid are inversely calculated, which are recorded as (r′, a′, e′), where r′ represents the range spherical coordinate corresponding to the X band, a′ represents the azimuth angle corresponding to the X band, and e′ represents the elevation angle corresponding to the X band; The reflectivity factor of each three-dimensional grid in the S band is calculated to generate the reflectivity factor data of the contour surface in the S band. The reflectivity factor of each three-dimensional grid in the X band is calculated to generate the reflectivity factor data of the contour surface in the X band.

5. A multi-band weather radar reflectivity factor networking method as claimed in claim 4, characterized in that: The calculation formula for inversely calculating the elevation angle, azimuth angle, and range spherical coordinates of each three-dimensional grid corresponding to the S-band radar site is: where lon c ,lat c is the longitude and latitude of the S-band radar site, lon, lat are the longitude and latitude of any grid cell in the S-band Cartesian three-dimensional grid, and A is calculated as follows: The calculation formula of S is: S=cos -1 (son(years c )×sin(lat)+cos(lat c )×something(lat)×something(lon-lon c )) where h c is the height of the S-band radar, and h is the height of any grid cell in the S-band Cartesian three-dimensional grid; The slope distance r is calculated as follows: The calculation formula for inversely calculating the elevation angle, azimuth angle, and range spherical coordinates of each three-dimensional grid corresponding to the X-band radar site is: where lon′ c , lat′ c is the longitude and latitude of the X-band radar site, lon′, lat′ are the longitude and latitude of any grid cell in the X-band Cartesian three-dimensional grid, and A′ is calculated as follows: The calculation formula of S′ is: S′=something -1 (sin(lat′ c )×sin(lat′)+cos(lat′ c )×cos(lat′)×cos(lon′-lon′ c )) where h′ c is the height of the radar, h′ is the height of any grid cell in the Cartesian three-dimensional grid; The slope distance r′ is calculated as follows:

6. A multi-band weather radar reflectivity factor networking method as claimed in claim 4, characterized in that: The calculation formula for calculating the reflectivity factor of each three-dimensional grid in the S band is: Among them, f(r,a,e) is the reflectivity factor of the coordinate (r,a,e), e1 and e2 are the upper and lower elevation angles closest to the elevation angle e, f(r,a,e1) and f(r,a,e2) are the reflectivity factor values ​​corresponding to the upper and lower elevation angles closest to the azimuth a and distance r, respectively; The calculation formula for calculating the reflectivity factor of each three-dimensional grid in the X band is: Wherein, f(r′,a′,e′) is the reflectivity factor of the coordinate (r′,a′,e′), e′1 and e′2 are the upper and lower elevation angles closest to the elevation angle e′, and f(r′,a′,e′1) and f(r′,a′,e′2) are the reflectivity factor values ​​corresponding to the upper and lower elevation angles closest to the azimuth a′ and distance r′, respectively.

7. A multi-band weather radar reflectivity factor networking method as claimed in claim 4, characterized in that: The S5 comprises the steps of: For the area covered by S-band weather radar detection and the area not covered by X-band weather radar detection, the S-band weather radar detection reflectivity factor is used as the reflectivity factor value after networking; For the area not covered by S-band weather radar detection and the area covered by X-band weather radar detection, the reflectivity factor of X-band weather radar detection is used as the reflectivity factor value after networking; For areas not covered by both S-band weather radar and X-band weather radar detection, there is no reflectivity factor value after networking; For areas covered by both S-band weather radar and X-band weather radar, the larger reflectivity factor of the two is used as the reflectivity factor value after networking.

8. A multi-band weather radar reflectivity factor networking method as claimed in claim 1, characterized in that: N＝10。 9. A multi-band weather radar reflectivity factor networking device, characterized in that: include: A data acquisition module is used to acquire S-band weather radar stream data with a time resolution of 1 minute and X-band weather radar base data with a time resolution of 3 minutes; The base data construction module is used to match the stream data closest to the current time in the S-band weather radar stream data, and use the stream data closest to the current time as the end stream data, match N stream data forward, N is a preset number, and construct the S-band weather radar base data with a time resolution of 1 minute according to the standard format of weather radar base data from the matched N+1 stream data; The radar site parameter extraction module is used to parse the reconstructed S-band weather radar base data and X-band weather radar base data, and extract the radar site parameters of the S-band weather radar at 9 elevation angles and the radar site parameters of the X-band weather radar at 9 elevation angles, respectively. Each set of radar site parameters includes the radar site longitude, radar site latitude, azimuth, distance library, and reflectivity factor. A contour surface reflectivity factor data generation module is used to generate the contour surface reflectivity factor data of the S band according to the radar site parameters of the 9 elevation angles of the S band, and to generate the contour surface reflectivity factor data of the X band according to the radar site parameters of the 9 elevation angles of the X band; The fusion networking module is used to fuse and network the S-band contour surface reflectivity factor data and the X-band contour surface reflectivity factor data.

10. A multi-band weather radar reflectivity factor networking device as claimed in claim 9, characterized in that: The radar site parameter extraction module includes radar site parameter extraction submodules distributed on each radar single station, and the radar site parameter extraction submodule on each radar single station is used to extract radar site parameters of 9 elevation angles in sequence, and after extracting the radar site parameters of a single elevation angle, while extracting the radar site parameters of the next elevation angle, the radar site parameters of the single elevation angle are saved; The contour surface reflectivity factor data generation module includes a contour surface reflectivity factor data generation submodule distributed on each radar station. The contour surface reflectivity factor data generation submodule on each radar station is used to generate the contour surface reflectivity factor data of each altitude in sequence, and after extracting the contour surface reflectivity factor data of a single altitude, extract the contour surface reflectivity factor data of the next altitude, and save the contour surface reflectivity factor data of the single altitude at the same time.

Images