CN112558076A - Volume scanning mode calculation method based on networking weather radar coverage area and application - Google Patents

Abstract

The invention discloses a volume scanning mode calculation method based on a networking weather radar coverage area and application thereof, and relates to the technical field of networking radars. The method comprises the following steps: constructing an initial weather radar net three-dimensional lattice point field based on networking weather radar data in a target area, and identifying lattice point coverage information in the three-dimensional lattice point field; determining a blind area, and identifying the boundary of the blind area; collecting data of radars to be added, identifying and configuring volume scanning mode parameters of the radars for each radar to be added based on a blind area, and fusing minimum azimuth angles, maximum azimuth angles and scanning ranges among set heights through calculation to obtain final sector scanning parameters. The method and the device can identify the blind area based on the radar coverage area information, and perform self-adaptive configuration on the volume scanning mode parameters of the blind area by combining the position parameters and other information of the weather radar to be added, so as to provide decision basis for introducing the radar to perform blind area compensation data fusion, performing radar collaborative networking layout and the like for users.

Description

Volume scanning mode calculation method based on networking weather radar coverage area and application

Technical Field

The invention relates to the technical field of networking radar, in particular to a volume scanning mode calculation method and application based on a networking weather radar coverage area.

Background

The weather radar plays an important role in a modern meteorological comprehensive observation system, has good monitoring capability on a large and medium-scale weather system, and provides visual data for short-term weather forecast. With the continuous improvement of the performance of the weather radar, the number of the erected radar parts is gradually increased, and a weather radar net is formed.

The flexibility and variability of the geometric structure and the information fusion mode of the networking weather radar enable the networking weather radar to be adjusted and suitable for various weather radar scenes and applications, and the networking weather radar becomes one of the research hotspots in the field of the current weather radar. The networking weather radar can increase the measurement dimension and confidence coefficient, and improve the fault tolerance and robustness of the system, so that the system is widely applied. As an example, ground coverage is assessed, such as using weather radar in combination with geographic information; for example, in the effective coverage and terrain occlusion analysis of the new generation weather radar netting design, the coverage range of the effective observation area of the radar under different observation elevation angles is calculated by combining high-resolution terrain data; for example, radar net coverage is combined with specific weather forecast application, and the radar net coverage is applied to quantitative precipitation estimation, hail coverage detection and the like.

The application of the networking weather radar mainly focuses on the aspects that the radar is shielded by high and large terrains to cause detection blind areas and influence on quantitative application such as precipitation estimation and the like, and generally focuses on the lower layer of a troposphere. However, with the development of the demands of urban fine management, smart weather forecast and service, the vertical observation of the weather elements over the cities is more and more urgent. Particularly, with the development of super-large city tests, the research on the refined landing points of multiple radar coverage including blind areas on a three-dimensional lattice point field of a large city and a surrounding area is increasingly important, so that a model can be provided for the research on the real construction of the three-dimensional reflectivity field, the refined prediction of cities on a surface scale and the analysis of a vertical reflectivity profile under non-interpolation on a point scale, and the model can be used for guiding how to more scientifically deploy atmosphere vertical observation equipment and how to better apply the existing observation data. The position of the networking weather radar and the volume coverage pattern vcp (volume coverage pattern) affect the coverage and the falling area of the three-dimensional space scale, and affect the acquisition of the space atmosphere scanning information.

The radar volume scanning mode is a mode in which the radar operates according to configuration parameters such as a preset pulse width, a preset pulse repetition frequency, a preset number of elevation layers, a preset antenna rotation rate and the like, and a scanning mode is usually set manually in the prior art. Along with the implementation of the concept of 'urban fine management' and the implementation of the 'first defense line of meteorological disaster prevention and reduction', the requirements for urban aerial fine observation and weather radar coverage area assessment are increasingly highlighted, and the scanning strategy highly dependent on manual operation obviously influences the detection precision to a certain extent.

In summary, since the networking weather radar joint observation puts higher requirements on spatial layout, cooperative observation and vertical profile, how to intelligently set the volume scanning strategy of the networking cooperative radar according to the existing radar coverage information to fully exert the detection capability of the radar and construct a more detailed and balanced stereo detection network is a technical problem which needs to be solved urgently at present.

Disclosure of Invention

The invention aims to: the defects of the prior art are overcome, and a volume scanning mode calculation method and application based on a networking weather radar coverage area are provided. The volume scanning mode calculation method based on the networking weather radar coverage area can identify the blind area based on the radar coverage area information, and self-adaptive configuration is carried out on volume scanning mode parameters of the blind area by combining information such as position parameters of weather radars to be added, so that decision basis is provided for users to introduce radars for blind area compensation data fusion, radar collaborative networking layout and the like.

In order to achieve the above object, the present invention provides the following technical solutions:

a volume scanning mode calculation method based on a networking weather radar coverage area comprises the following steps:

constructing an initial weather radar mesh three-dimensional grid field based on networking weather radar data in a target area, and identifying grid point coverage information in the initial weather radar mesh three-dimensional grid field, wherein the grid point coverage information comprises blind point identification and radar coverage point identification;

obtaining blind spot information of marks in a three-dimensional grid spot field of an initial weather radar net, determining a blind area region, and identifying the boundary of the blind area region;

collecting data of radars to be added, identifying and configuring volume scanning mode parameters of the radars to be added based on a blind area for each radar to be added, and comprising the following steps of: acquiring longitude and latitude, altitude and radial distance parameters of the radar, and establishing a three-dimensional coordinate system of the radar; judging whether the blind point is in an effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as an effective point when the blind point is in the effective radial range; calculating the azimuth angle and the initial angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range at each height; and according to the set height range, obtaining the minimum azimuth angle, the maximum azimuth angle and the scanning range among the set heights, and then fusing to obtain the final fan-scanning parameters.

Further, after obtaining the final fan-scanning parameters, the following steps are also executed:

and acquiring a pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain a final pitch angle parameter.

Further, the step of constructing the initial weather radar mesh three-dimensional grid field includes,

acquiring position data and volume scanning mode parameter data of a plurality of networking weather radars in a target area, wherein the position data comprises radar longitude and latitude and feed source altitude;

for each networking weather radar, converting the polar coordinates into a lattice field under Cartesian coordinates according to the position data and the volume scanning mode parameter data of the networking weather radar, and calculating the longitude range and the latitude range detected by each networking weather radar according to the detection radius of each networking weather radar;

constructing an initial weather radar mesh three-dimensional lattice field by combining the lattice fields of all the networking weather radars, and determining the horizontal range, the horizontal resolution, the height range, the equidistant height hs and the spatial resolution hr of the three-dimensional lattice field; the horizontal range of the three-dimensional lattice field of the initial weather radar net is determined by respectively comparing the longitude range and the latitude range detected by all the networking weather radars, wherein the horizontal range comprises a longitude interval range and a latitude interval range, the longitude interval range is the sum of the distance value between the radar with the largest detection longitude and the radar with the smallest detection longitude, the detection radius of the radar with the largest detection longitude and the detection radius of the radar with the smallest detection longitude, and the latitude interval range is the sum of the distance value between the radar with the largest detection latitude and the radar with the smallest detection latitude, the detection radius of the radar with the largest detection latitude and the detection radius of the radar with the smallest detection latitude; the horizontal resolution, height range, equidistant height hs, and spatial resolution hr are set by the user or the system.

Further, grid point coverage information in the three-dimensional grid point field of the initial weather radar mesh is identified by analyzing coverage information on the equal altitude surface, the steps are as follows,

lattice point conversion: according to the set height range, all grid points on the equidistant height hs on the set height Z are obtained, and any grid point P in the grid points is_tObtaining a grid point P_tCoordinate (X) of_t,Y_t,Z_t) Wherein T =1,2, … …, T, T is the total number of lattice points, X_tIs latitude, Y_tIs longitude, Z_tIs the height; determination of grid point P by radar beam propagation and large circle geometry theory_tPolar position (r) in a radar polar coordinate system_t,θ_t,φ_t) Wherein r is_tIs the pitch, θ_tIs an azimuth angle phi_tIs a pitch angle;

identification: for the foregoingArbitrary lattice point P_tOne by one will (r)_t,θ_t,φ_t) Matching with azimuth, elevation, beam width and distance information of single radar, and determining (r)_t,θ_t,φ_t) When the grid point P is in the detection range_tIdentifying as a radar coverage point, otherwise, identifying the grid point as a blind point; after all radars in the traversing weather radar network identify the grid points according to the method, the radar coverage quantity of each grid point is obtained, when one or more radars cover the grid point, the grid point is identified as a radar coverage point, and when no radar covers the grid point, the grid point is identified as a blind point.

Further, grid point coverage information in the initial weather radar mesh three-dimensional grid point field is identified by combining coverage information among equal-altitude layers, and the method comprises the following steps:

based on a single radar, all grid points on the space resolution hr at the set height Z are obtained according to the set height range, and any grid point P in the grid points is_jObtaining a grid point P_jCoordinate (X) of_j,Y_j,Z_j) And coordinate conversion is carried out to obtain the corresponding polar coordinate position (r)_j,θ_j,φ_j) Wherein j =1,2, … …, m, m is the total number of grid points; will (r)_j,θ_j,φ_j) Matching the azimuth, elevation, beam width and distance information of the single radar, and determining (r)_j,θ_j,φ_j) When the grid point P is in the detection range_jIdentifying as a radar coverage point, otherwise, identifying the grid point as a blind point;

regarding the height interval between any two adjacent equidistant heights as a layer L_q，L_qRepresenting a q-th layer, wherein the value of q is an integer which is more than or equal to 1; according to the set height Z, the lattice point P on the equal distance height hs_tLattice point P with spatial resolution hr_jMatching, and marking the grid points with the same height as repeated grid points; and acquiring all layers at a set height Z, L for each layer_qCalculate the layer L_qThe corresponding number of lattice points is equal to the value of the layer after the duplication removal of the sum of the number of lattice points on the equidistant height hs corresponding to the layer and the number of lattice points on the corresponding spatial resolution hr, and the layer is judgedL_qWhether the corresponding grid points are all blind points or not, and if so, the layer L is divided into the blind points_qMarking as a blind spot, otherwise, only any grid point is a radar coverage point, and then, marking the layer L as a blind spot_qMarking as a radar coverage point;

and after all radars in the traversing weather radar network identify the layers according to the method, acquiring the radar coverage quantity of each layer, identifying the radar coverage points when one or more radars cover the layer, and identifying the blind points when no radar covers the layer.

Further, the blind spot is identified by using a number 0, and when a radar coverage point is identified, a corresponding number k is identified according to the radar coverage quantity on the grid point, wherein k =1,2, … …, N is the total number of radars, and the radar coverage point comprises a single radar coverage point, a double radar coverage punctuation point, … … and k radar coverage points.

Further, the radar to be added is an X-band radar.

The invention also provides a collaborative adaptive scanning method of the networking weather radar, which comprises the following steps:

acquiring radar data to be added, which can be added into a target area weather radar network, in a preset area;

according to the method, volume scanning mode parameters of the radar to be added are configured, and observation is carried out.

The invention also provides a volume scanning mode calculation device based on the networking weather radar coverage area, which comprises the following structures:

the grid point identification module is used for constructing an initial weather radar grid three-dimensional grid point field according to networking weather radar data in a target area, and identifying grid point coverage information in the initial weather radar grid three-dimensional grid point field, wherein the grid point coverage information comprises blind point identification and radar coverage point identification;

the blind area identification module is used for determining a blind area after acquiring the blind point information of the mark in the three-dimensional grid point field of the initial weather radar network and identifying the boundary of the blind area;

the volume scanning mode processing module is used for collecting data of the radar to be added, and configuring volume scanning mode parameters of the radar according to the following steps based on blind area identification for each radar to be added:

acquiring longitude and latitude, altitude and radial distance parameters of the radar, and establishing a three-dimensional coordinate system of the radar; judging whether the blind point is in an effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as an effective point when the blind point is in the effective radial range; calculating the azimuth angle and the initial angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range at each height; and according to the set height range, obtaining the minimum azimuth angle, the maximum azimuth angle and the scanning range among the set heights, and then fusing to obtain the final fan-scanning parameters.

Further, the volume scanning mode processing module is further configured to:

and after the final sector scanning parameters are obtained, acquiring the pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain the final pitch angle parameters.

Due to the adoption of the technical scheme, compared with the prior art, the invention has the following advantages and positive effects as examples:

on one hand, the blind area can be identified based on the radar coverage area information, the volume scanning mode parameters of the blind area can be configured in a self-adaptive mode by combining information such as position parameters of weather radars to be added, and decision basis is provided for users to introduce radars for blind supplement data fusion, radar collaborative networking layout and the like.

On the other hand, when grid point coverage information in a three-dimensional grid point field is identified, an equal-height interlayer coverage algorithm is provided, coverage information identification can be carried out between equal-height layers, information of grid points at equal heights and grid points at spatial resolution is considered when a coverage area is identified through layer information combination, the problem of blind area expansion caused by only considering the coverage information on the equal-height surface is avoided, and the radar simulation precision efficiency is improved.

On the other hand, when the three-dimensional lattice field is constructed, the horizontal range area of the three-dimensional lattice field can be intelligently set according to the position information, the detection range and other information of each networking weather radar, and the construction precision and efficiency of the lattice field are improved.

Drawings

Fig. 1 is a schematic flowchart of a method for calculating a volume scanning pattern based on a networked weather radar coverage area according to an embodiment of the present invention.

Fig. 2 is a flowchart illustrating a process of calculating a radar volume scanning pattern according to an embodiment of the present invention.

Fig. 3 is a flowchart illustrating a process of identifying grid coverage information according to an embodiment of the present invention.

Fig. 4 is a diagram illustrating an example of setting of a horizontal range of a three-dimensional lattice point field of a 3-dimensional lattice weather radar when the 3-dimensional lattice weather radar is set according to an embodiment of the present invention.

Fig. 5 is a diagram illustrating a calculation example for obtaining overlay information on an equal-height surface and overlay information between equal-height layers according to an embodiment of the present invention.

Detailed Description

The method and application of calculating the volume scanning mode based on the coverage area of the networked weather radar disclosed by the invention are further described in detail with reference to the accompanying drawings and specific embodiments. It should be noted that technical features or combinations of technical features described in the following embodiments should not be considered as being isolated, and they may be combined with each other to achieve better technical effects. In the drawings of the embodiments described below, the same reference numerals appearing in the various drawings denote the same features or elements, which may be applied to different embodiments. Thus, once an item is defined in one drawing, it need not be further discussed in subsequent drawings.

It should be noted that the structures, proportions, sizes, and other dimensions shown in the drawings and described in the specification are only for the purpose of understanding and reading the present disclosure, and are not intended to limit the scope of the invention, which is defined by the claims, and any modifications of the structures, changes in the proportions and adjustments of the sizes and other dimensions, should be construed as falling within the scope of the invention unless the function and objectives of the invention are affected. The scope of the preferred embodiments of the present invention includes additional implementations in which functions may be executed out of order from that described or discussed, including substantially concurrently or in reverse order, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present invention.

Techniques, methods, and apparatus known to those of ordinary skill in the relevant art may not be discussed in detail, but are intended to be part of the specification where appropriate. In all examples shown and discussed herein, any particular value should be construed as merely illustrative, and not limiting. Thus, other examples of the exemplary embodiments may have different values.

In addition, the blind spot and the radar coverage point referred to in the present invention may be a geometric point (without dimension, with uniqueness, such as a grid point) or a block (with dimension, with uniqueness, such as a layer with height). Similar to the prior art pixels, they are often referred to as pixels, but are actually a block.

Examples

The position and the volume scanning mode of the networking weather radar influence the coverage and the falling area of a three-dimensional space scale and influence the acquisition of space atmosphere scanning information. At present, the basic scanning modes of a radar VCP are mainly classified into a clear sky mode and a precipitation mode, and generally include 9 volume scanning modes, which are respectively: VCP32, VCP31, VCP21, VCP11, VCP12, VCP121, VCP211, VCP212 and VCP221, wherein the scanning mechanisms of VCP211, VCP221 and VCP212 correspond to VCP11, VCP21 and VCP12, respectively, except that a phase encoding algorithm is employed at two elevation angles of the lower layer to improve detection of unambiguous velocity.

For clear sky mode, there are mainly two kinds of VCP31 and VCP32, including 5 elevation angles, and the scanning time is about 10 minutes. VCP31 has a pulse width of 4.7us and a PRF of approximately 322 Hz; and the pulse width of VCP32 is 1.57us, the PRF is about 1013Hz, and the detection precision and the speed measurement range are higher. For the precipitation mode, the VCP11 and VCP21 are mainly used for detecting, tracking and analyzing the precipitation type and the characteristics of strong convection weather, and the VCP12 mode is introduced for meeting the field observation requirement in the strong convection weather at the later stage. The details of the VCPs in precipitation mode are shown in table 1.

Table 1 weather radar common body scanning mode parameter table

The radar volume scanning mode is a mode in which the radar operates according to configuration parameters such as a predetermined pulse width, a predetermined pulse repetition frequency, a predetermined number of elevation layers, and a predetermined antenna rotation rate. When meteorological detection is performed, volume scanning mode information such as a sector scanning parameter (for example, a sector scanning center, a sector scanning range and the like) and a pitch angle parameter (for example, the number of elevation angles, the angle of elevation and the like) of a radar needs to be set through a related console, and the console sends parameter setting adjustment to change an antenna scanning mode.

According to the technical scheme provided by the invention, the blind area can be identified based on the radar coverage area information, and partial key body scanning mode parameters can be automatically configured by combining the position parameters and other information of the weather radar to be added, so that the dependence on manual operation is reduced, and a decision basis is provided for users to introduce the radar for blind area compensation data fusion, radar collaborative networking layout and the like.

Referring to fig. 1, a volume scanning mode calculation method based on a networking weather radar coverage area is provided in the present invention. The method comprises the following steps:

s100, constructing a three-dimensional grid point field of the networking weather radar and marking grid point coverage information.

Firstly, constructing an initial weather radar mesh three-dimensional grid field based on networking weather radar data in a target area, and identifying grid coverage information in the initial weather radar mesh three-dimensional grid field, wherein the grid coverage information comprises blind spot identification and radar coverage point identification.

In this embodiment, blind spots and radar coverage points may be identified numerically. As a typical preference, the blind spot is identified by the number 0. When identifying a radar coverage point, a corresponding number k may be identified according to the number of radar coverage on the grid point, where k =1,2, … …, N is the total number of radars, and the radar coverage point includes a single radar coverage point, a double radar coverage punctuation point, … …, k radar coverage points.

By way of example and not limitation, for example, N =3, then k =1,2, 3, respectively corresponding to a single radar coverage point, a double radar coverage punctuation point, and a triple radar coverage point.

Determining a blind area and identifying the boundary of the blind area.

And then, obtaining the blind spot information identified in the three-dimensional grid spot field of the initial weather radar net, determining a blind area region, and identifying the boundary of the blind area region.

Referring to fig. 2, according to the identification numbers of the blind spot and the radar coverage point, the step S200 may specifically include the following steps:

s210, determining a blind area: and acquiring blind spot information marked as 0 in the three-dimensional grid spot field of the initial weather radar network, and connecting all the blind spots marked as 0 to form a blind area region.

And (3) identifying the boundary of the blind area: and identifying the boundary of the blind area region through a preset identification model.

And collecting data of the radars to be added, and identifying and configuring volume scanning mode parameters of the radars for each radar to be added based on the blind area.

The radar to be added can be one or more, for example, weather radar outside the aforementioned target area, and the user can select the radar according to the needs. In this embodiment, the radar to be added is preferably an X-band radar. The X-band dual-polarization short-distance network radar with the rapid adaptive scanning capability (the detection radius of the radar is 30-40km, and the average distance between the radars is 25km) can make up the defect of the low-layer detection capability of the conventional radar network.

Specifically, with continued reference to fig. 2, the step of configuring the volume scanning mode parameters to be added to the radar based on the blind area identification may be as follows:

s310, acquiring parameters such as longitude and latitude, altitude, radial distance and the like of the newly added X-band radar, and establishing a three-dimensional coordinate system of the radar.

And judging whether the blind point is in the effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as the effective point when the blind point is in the effective radial range. The effective radial range is set by a user or a system.

And calculating the azimuth angle and the starting angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range on each height.

And according to the set height range, obtaining the minimum azimuth angle, the maximum azimuth angle and the scanning range among the set heights, and then fusing to obtain the final fan-scanning parameters.

And acquiring a pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain a final pitch angle parameter.

And forming a final volume scanning mode of the radar according to the final fan-scanning parameter and the final pitch angle parameter. Then, observation is performed based on the aforementioned final volume scan pattern.

Referring to fig. 3, in this embodiment, a specific process of constructing the initial weather radar mesh three-dimensional lattice field in step S100 is as follows:

firstly, position data and volume scanning mode parameter data of a plurality of networking weather radars in a target area are obtained, wherein the position data comprises radar longitude and latitude and feed source altitude.

And then, for each networking weather radar, converting the polar coordinates into a lattice field under Cartesian coordinates according to the position data and the volume scanning mode parameter data of the networking weather radar, and calculating the longitude range and the latitude range detected by each networking weather radar according to the detection radius of each networking weather radar.

And finally, combining the grid point fields of the various networking weather radars to construct an initial weather radar grid three-dimensional grid point field, and determining the horizontal range, the horizontal resolution, the height range, the equidistant height hs and the spatial resolution hr of the three-dimensional grid point field.

The horizontal range of the three-dimensional grid point field of the initial weather radar net is determined by respectively comparing the longitude range and the latitude range detected by all the networking weather radars.

Specifically, the horizontal range includes a longitude range and a latitude range. The longitude distance range is the sum of the distance value between the radar with the largest detection longitude and the radar with the smallest detection longitude, the detection radius of the radar with the largest detection longitude and the detection radius of the radar with the smallest detection longitude. The latitude distance range is the sum of the distance value between the radar with the largest detection latitude and the radar with the smallest detection latitude, the detection radius of the radar with the largest detection latitude and the detection radius of the radar with the smallest detection latitude.

The horizontal resolution, height range, equidistant height hs, and spatial resolution hr are set by the user or the system.

Preferably, the horizontal resolution is set to 1 km; the spatial resolution hr is set to a multiple of 250 meters, such as 250m, 500m, 750m, 1km or 2km, and may be set as desired.

Referring to fig. 4, the following describes how to determine the horizontal range of the three-dimensional lattice field in detail, taking the example that the weather radar net includes 3 radars.

The weather radar net comprises a radar A, a radar B and a radar C. The longitude and latitude of the A radar are (LonA, LatA), and the detection radius is Ra; the longitude and latitude of the radar B are (LonB, LatB), and the detection radius is Rb; the longitude and latitude of the C radar are (LonC, LatC), and the detection radius is Rc.

According to the positional relationship of the radar a, the radar B, and the radar C shown in fig. 4, the radar C is located at the maximum detection longitude (Lon) (corresponding to the rightmost side of fig. 4), and the radar a is located at the minimum detection longitude (corresponding to the leftmost side of fig. 4), and then the range of the longitude interval of the three-dimensional lattice field is the sum of the distance value between the radar C with the maximum detection longitude and the radar a with the minimum detection longitude, the detection radius of the radar C with the maximum detection longitude, and the detection radius of the radar a with the minimum detection longitude. That is, the longitude range L1 of the three-dimensional lattice field is the detection radius Ra of the a radar and the distance d between the AC radars_ACAnd the detection radius Rc of the C radar, i.e. L1= Ra + d_AC+Rc。

According to the position relationship among the radar a, the radar B, and the radar C shown in fig. 4, the radar B is the largest detected latitude (Lat) (corresponding to the uppermost edge of fig. 4), and the radar B is the smallest detected latitude(corresponding to the lowest side of fig. 4) is radar C, the range of the latitude interval of the three-dimensional lattice field is the sum of the distance value between radar B with the largest detection latitude and radar C with the smallest detection latitude, the detection radius of radar B with the largest detection latitude and the detection radius of radar C with the smallest detection latitude. That is, the range of the latitude interval L2 of the three-dimensional lattice field is the distance d between the detection radii Rb and BC of the B radar_BCAnd the detection radius Rc of the C radar, i.e. L1= Rb + d_BC+Rc。

With continued reference to fig. 3, in step S100, the grid point coverage information may be identified by analyzing the coverage information on the equal-height surface and the coverage information between the equal-height layers.

Specifically, the step of identifying the grid point coverage information in the initial weather radar mesh three-dimensional grid point field through the coverage information on the equal altitude surface may be as follows:

and S121, a lattice point conversion step.

According to the set height range, all grid points on the equidistant height hs on the set height Z are obtained, and any grid point P in the grid points is_tObtaining a grid point P_tCoordinate (X) of_t,Y_t,Z_t) Wherein T =1,2, … …, T, T is the total number of lattice points, X_tIs latitude, Y_tIs longitude, Z_tIs the height; determination of grid point P by radar beam propagation and large circle geometry theory_tPolar position (r) in a radar polar coordinate system_t,θ_t,φ_t) Wherein r is_tIs the pitch, θ_tIs an azimuth angle phi_tIs a pitch angle.

And identifying.

For the aforementioned arbitrary lattice point P_tOne by one will (r)_t,θ_t,φ_t) Matching with azimuth, elevation, beam width and distance information of single radar, and determining (r)_t,θ_t,φ_t) When the grid point P is in the detection range_tAnd identifying the grid points as radar coverage points, otherwise, identifying the grid points as blind points.

After all radars in the traversing weather radar network identify the grid points according to the method, the radar coverage quantity of each grid point is obtained, when one or more radars cover the grid point, the grid point is identified as a radar coverage point, and when no radar covers the grid point, the grid point is identified as a blind point.

The step of identifying the grid point coverage information in the initial weather radar mesh three-dimensional grid point field by the coverage information between the equal altitude layers may be as follows:

s131, based on a single radar, obtaining all grid points on the set height Z with the spatial resolution hr according to the set height range, and obtaining any grid point P in the grid points_jObtaining a grid point P_jCoordinate (X) of_j,Y_j,Z_j) And coordinate conversion is carried out to obtain the corresponding polar coordinate position (r)_j,θ_j,φ_j) Wherein j =1,2, … …, m, m is the total number of grid points; will (r)_j,θ_j,φ_j) Matching the azimuth, elevation, beam width and distance information of the single radar, and determining (r)_j,θ_j,φ_j) When the grid point P is in the detection range_jAnd identifying the grid points as radar coverage points, otherwise, identifying the grid points as blind points.

The height interval between any two adjacent equidistant heights is regarded as a layer L_q，L_qRepresenting a q-th layer, wherein the value of q is an integer which is more than or equal to 1; according to the set height Z, the lattice point P on the equal distance height hs_tLattice point P with spatial resolution hr_jMatching, and marking the grid points with the same height as repeated grid points; and acquiring all layers at a set height Z, L for each layer_qCalculate the layer L_qThe corresponding number of lattice points, which is equal to the number of the lattice points on the equidistant height hs corresponding to the layer and the number of the lattice points on the corresponding spatial resolution hr after the duplication is removed, the layer L is judged_qWhether the corresponding grid points are all blind points or not, and if so, the layer L is divided into the blind points_qMarking as a blind spot, otherwise, if any grid point is a radar coverage point, the layer L is marked as a blind spot_qThe markers are radar coverage points.

And after all radars in the traversing weather radar network identify the layers according to the method, acquiring the radar coverage quantity of each layer, identifying the radar coverage points when one or more radars cover the layer, and identifying the blind points when no radar covers the layer.

The manner in which the coverage areas are identified by the overlay information on the equal-altitude surface and between the equal-altitude layers is described in detail below in conjunction with fig. 5.

The set height Z, equidistant height hs, and spatial resolution hr may be set by a user or system as desired. Taking the example shown in fig. 5, height Z is set to include equidistant heights hs1, hs2, hs3 and hs4, including spatial resolutions hr1, hr2, hr3, hr4, hr5, hr6, hr7 and hr8, hs1 includes hr1 and hr2, hs2 includes hr3 and hr4, hs3 includes hr5 and hr6, hs4 includes hr7 and hr8, i.e. equidistant heights hs =2 hr.

When the coverage information on the equal-height surface is analyzed, it is found that 5 grid points, i.e., T =5, are present on the equidistant height hs at the set height Z, corresponding to P in fig. 5_i，P_i+2，P_i+4，P_i+6，P_i+8。

For an arbitrary grid point P, taking a single radar as an example, such as radar A_tBy its corresponding polar position (r)_t,θ_t,φ_t) And judging whether the grid points are radar coverage points of the radar A, if so, identifying the grid points as 1, otherwise, identifying the grid points as 0, namely blind points. For 5 lattice points P_i，P_i+2，P_i+4，P_i+6，P_i+8Are identified according to the above method.

Then, traversing all radars in the weather radar network, including the radar A, the radar B and the radar C, identifying the grid points according to the method, and acquiring the grid point P_i，P_i+2，P_i+4，P_i+6，P_i+8K, of radar coverage. Specifically, k may be set to an initial value of 0, and k + + may be performed every time 1 radar coverage is determined.

A grid point is identified as a radar coverage point when it has one or more radar coverage, and a value of k is identified, and a grid point is identified as a 0 when it has no radar coverage, i.e., a blind spot.

When the contour interlayer coverage information is analyzed, it is known that there are 9 grid points in the spatial resolution hr at the set height Z, i.e., m =9, corresponding to P in fig. 5_i，P_i+1，P_i+2，P_i+3，P_i+4，P_i+5，P_i+6，P_i+7，P_i+8。

Taking a single radar, such as radar A, for any grid point P therein_jBy its corresponding polar position (r)_j,θ_j,φ_j) And judging whether the grid points are radar coverage points of the radar A, if so, identifying the grid points as 1, otherwise, identifying the grid points as 0, namely blind points. For 9 grid points P_i，P_i+1，P_i+2，P_i+3，P_i+4，P_i+5，P_i+6，P_i+7，P_i+8Are identified according to the above method.

Regarding the height interval between any two adjacent equidistant heights as a layer L_q，L_qRepresenting the q-th level, a total of 4 levels in fig. 5, i.e. q =4, corresponding to L in the figure_i+1，L _i+2，L _i+3，L _i+4. Lattice point P on the peer-to-peer distance height hs_tLattice point P with spatial resolution hr_jMatching is performed, and grid points with the same height are marked as repeated grid points, in FIG. 5, P_i，P_i+2，P_i+4，P_i+6，P_i+8Are repeating grid points.

For all layers L at a set height Z_i+1，L _i+2，L _i+3And L _i+4For each layer L_qCalculate the layer L_qThe corresponding number of lattice points, which is equal to the number of the lattice points on the equidistant height hs corresponding to the layer and the number of the lattice points on the corresponding spatial resolution hr after the duplication is removed, the layer L is judged_qWhether the corresponding grid points are all blind points or not, and if so, the layer L is divided into the blind points_qMarking as a blind spot, otherwise, only any grid point is a radar coverage point, and then, marking the layer L as a blind spot_qThe markers are radar coverage points. With L_i+1For example, the number of corresponding grid points is equal to the number of grid points on the equidistant height hs (i.e. P) corresponding to the layer_iAnd P_i+2Wherein P is_iIs a lower equidistant height grid point of the layer, P_i+2Equal distance height grid points of a layer) and a corresponding number of spatial resolution hr grid points (i.e., P)_i，P_i+1And P_i+2) The sum of (a) and (b) is the value after de-duplication, i.e. 2+3-2=3, L_i+1The corresponding 3 grid points are respectively P_i，P_i+1，P_i+2And (5) grid points. Fault judgment L_i+1Whether the corresponding 3 grid points are all blind points or not, if the 3 grid points are all blind points, the layer L is divided into three layers_i+1Marking as 0, namely blind spot, otherwise, if any grid point is a radar coverage point, the layer L is marked_qThe label is 1, i.e. the radar coverage point.

And by analogy, all the layers are identified according to the steps. Specifically, in FIG. 5, L_i+2The corresponding 3 grid points are respectively P_i+2，P_i+3，P_i+4Lattice points, L_i+3The corresponding 3 grid points are respectively P_i+4，P_i+5，P_i+6Lattice points, L_i+4The corresponding 3 grid points are respectively P_i+6，P_i+7，P_i+8And (5) grid points. By way of example, and not limitation, layer L in FIG. 5_i+1As radar coverage points, L _i+2As radar coverage points, L _i+3Is a blind spot, L _i+4Is a radar coverage point.

And traversing all radars in the weather radar network, including the radar A, the radar B and the radar C, and identifying the layers according to the method to obtain the radar coverage quantity of each layer. Similarly, k can be initialized to 0 and executed every 1 radar coverage is determined.

A radar coverage point is identified when a layer has one or more radar coverage and a value of k is identified, and a value of 0, i.e. a blind spot, is identified when a layer has no radar coverage.

After grid point coverage information (including coverage information on an equal-altitude surface and coverage information between equal-altitude layers) of the initial weather radar grid three-dimensional grid point field is determined, a blind area region can be determined according to the blind point (including grid points and layers) information of the marked 0. After the boundary of the blind area is identified, the volume scanning mode parameters of the radar are configured based on the blind area identification.

According to the scheme provided by the invention, the scanning mode parameters of the newly added radar volume can be configured in a self-adaptive manner based on the blind area, and a decision basis is provided for users to introduce radars for blind area compensation data fusion, radar collaborative networking layout and the like. Meanwhile, by means of an equal-height interlayer covering algorithm, information of grid points at equal heights and grid points at spatial resolution is considered when a coverage area is identified, the problem of blind area expansion caused by only considering covering information on an equal-height surface is solved, and the radar simulation precision efficiency is improved.

The present embodiment will be described in detail below by taking an example in which the networked weather radar in a certain urban area includes 3 weather radars.

A networking weather radar network in a certain urban area comprises 1S-band WSR-88D weather radar, 1S-band CINRAD/SA weather radar and a shared C-band WRK-200 weather radar. The spacing between the WSR-88D weather radar and the CINRAD/SA weather radar is set to 90km, and the longitude of the WRK-200 weather radar is set between the two radars. When an initial weather radar network three-dimensional lattice field is constructed, the earth curvature, the beam broadening and all elevation angles of radar volume scanning are considered, and the beams are assumed to be transmitted under standard atmosphere, and the radar is observed by adopting different volume scanning modes VCP (volume coverage pattern) to acquire atmosphere information. The parameters of the volume scanning mode can comprise detection distance, radial resolution, pitch angle, azimuth angle and beam width, and each parameter value can be set as required. Wherein the beam width (also called beam spread) is typically set to 1 °.

At present, the WSR-88D weather radar and the CINRAD/SA radar mainly observe three types of modes of precipitation, namely VCP21, VCP12 and VCP 11. The VCP21 mode is suitable for stable lamellar cloud precipitation and is also a precipitation mode uniformly used by the domestic weather radar; the VCP12 mode is suitable for encrypted observation at a low-layer elevation angle; the VCP11 mode is suitable for encrypted observation at high elevation angles.

Firstly, constructing an initial weather radar mesh three-dimensional lattice field based on the 3 weather radars.

The three radars are located as shown in fig. 4, and the constructed horizontal range includes a longitude range (E) and a latitude range (N). The field longitude range consists of three parts: a distance between a maximum longitude radar and a minimum longitude radar, a sum of a radar radius of maximum longitude and a radar radius of minimum longitude. The inter-latitude distance of the field is composed of three parts, the distance between the maximum latitude radar and the minimum latitude radar, and the sum of the radiation radius of the minimum latitude radar and the radar radius of the maximum latitude. In this example, the horizontal range of the field is 500km long by 400km wide.

The grid point number of the three-dimensional grid point field constructed by the three radars is related to the grid point resolution. In this example, the ground is divided at equal intervals to 20km in height, with two resolutions: simulations were performed with normal resolution, i.e. 1km x 1km in horizontal plane, and height 0.5km and with high resolution, i.e. 0.5km x 0.5km in horizontal plane, and height 0.25 km. In this example, simulation was performed at high resolution, and the two-dimensional horizon (X, Z) constructed was divided 976 × 1107, i.e., 1080432 grid points, and the three-dimensional field was 1080432 × 80, i.e., 86434560 grid points.

In this example, the grid point coverage information is identified by analyzing the coverage information on the iso-level surface and the coverage information between the iso-level layers. Specifically, the grid point coverage information of the three-dimensional grid point field under the combination of the modes when the WSR-88D weather radar and the CINRAD/SA radar simultaneously adopt the VCP21, VCP11 and VCP12 modes and operate together with the WRK-200 weather radar operation mode is analyzed. Blind spots can be divided into two categories: one type is a blind spot (called as an outer blind spot for short) outside the boundary of a multi-radar coverage area in the two-dimensional surface range of the weather radar of each networking, and the change of the coverage range of the weather radar of the networking can be reflected by the increase and decrease of the number of the outer blind spots; the other type is blind spots (inner blind spots for short) in the boundary range of a multi-radar coverage area in the two-dimensional surface range of the weather radar of each layer network, and the increase and decrease of the number of the inner blind spots can reflect the conversion of multi-radar blind single-double multi-coverage.

Then, a blind area analysis of the three-dimensional lattice field was performed.

Through simulation analysis, the location of the blind spot in VCP21 mode was identified. The area of each layer of blind area at the height of 5-7km is larger than a preset threshold value and just above the central urban area, so that the observation benefit is greatly influenced. Therefore, the peripheral X-band radar is called to perform blind spot compensation.

The periphery of the city comprises two X-band radars, the longitude and latitude of the first radar are (E1 degrees, N1 degrees), the altitude is H1m, the longitude and latitude of the second radar are (E2 degrees, N2 degrees), and the altitude is H2 m.

The radar volume scanning mode calculation method provided by the invention is used for carrying out volume scanning mode parameter configuration and calculating the radial distance, the azimuth angle and the pitch angle. By calculation, at 5-7km, the azimuth starting and stopping angle of the first radar blind patch is 230.38-344.44 degrees, the radial distance starting and stopping range is about 20-80 km, and the pitch angle starting and stopping range is 4.4-26.37 degrees. If mechanical antenna scanning is selected, the pitch angle can be configured to increase by 1 deg., requiring 23 rotations. If the antenna rotates at 18/s-36/s a total of 230s-460s is required. And the radial distance of the blind area of the second radar is less than 70km at the height of 7km, but the horizontal range is 360 degrees, and the position of the blind area covers the position of the second radar, so that the blind area cannot be compensated.

The invention further provides a cooperative adaptive scanning method of the networking weather radar. The method comprises the following steps:

step 1, acquiring radar data to be added, which can be added into a target area weather radar network, in a preset area.

And 2, configuring the body scanning mode parameters of the radar to be added, and carrying out observation.

In this embodiment, the step of configuring the volume scanning mode parameter to be added to the radar specifically includes:

firstly, constructing an initial weather radar mesh three-dimensional grid field based on networking weather radar data in a target area, and identifying grid coverage information in the initial weather radar mesh three-dimensional grid field, wherein the grid coverage information comprises blind spot identification and radar coverage point identification.

In this embodiment, blind spots and radar coverage points may be identified numerically. As a typical preference, the blind spot is identified by the number 0. When identifying a radar coverage point, a corresponding number k may be identified according to the number of radar coverage on the grid point, where k =1,2, … …, N is the total number of radars, and the radar coverage point includes a single radar coverage point, a double radar coverage punctuation point, … …, k radar coverage points.

By way of example and not limitation, for example, N =3, then k =1,2, 3, respectively corresponding to a single radar coverage point, a double radar coverage punctuation point, and a triple radar coverage point.

And then, obtaining the blind spot information identified in the three-dimensional grid spot field of the initial weather radar net, determining a blind area region, and identifying the boundary of the blind area region.

The method specifically comprises the following steps: determining a blind area: and acquiring blind spot information marked as 0 in the three-dimensional grid spot field of the initial weather radar network, and connecting all the blind spots marked as 0 to form a blind area region. Identification of the boundary of the blind area: and identifying the boundary of the blind area region through a preset identification model.

And finally, collecting data of the radars to be added, and identifying and configuring volume scanning mode parameters of the radars for each radar to be added based on the blind area.

The radar to be added can be one or more, for example, weather radar outside the aforementioned target area, and the user can select the radar according to the needs. In this embodiment, the radar to be added is preferably an X-band radar. The X-band dual-polarization short-distance network radar with the rapid adaptive scanning capability (the detection radius of the radar is 30-40km, and the average distance between the radars is 25km) can make up the defect of the low-layer detection capability of the conventional radar network.

Specifically, the step of configuring the volume scanning mode parameters to be added into the radar based on the blind area identification is as follows: and acquiring parameters such as longitude and latitude, altitude, radial distance and the like of the newly added X-band radar, and establishing a three-dimensional coordinate system of the radar. And judging whether the blind point is in the effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as an effective point when the blind point is in the effective radial range. The effective radial range is set by a user or a system. And calculating the azimuth angle and the initial angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range at each height. And according to the set height range, obtaining the minimum azimuth angle, the maximum azimuth angle and the scanning range among the set heights, and then fusing to obtain the final fan-scanning parameters.

And acquiring a pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain a final pitch angle parameter. And forming a final volume scanning mode of the radar according to the final fan-scanning parameter and the final pitch angle parameter. Then, observation is performed based on the aforementioned final volume scan pattern.

Other technical features are referred to in the previous embodiments and are not described herein.

The invention further provides a volume scanning mode calculation device based on the networking weather radar coverage area.

The device comprises a lattice point identification module, a blind area identification module and a volume scanning mode processing module.

The grid point identification module is used for constructing an initial weather radar grid three-dimensional grid point field according to networking weather radar data in a target area, and identifying grid point coverage information in the initial weather radar grid three-dimensional grid point field, wherein the grid point coverage information comprises blind point identification and radar coverage point identification;

the blind area identification module is used for determining a blind area after acquiring the blind point information of the mark in the three-dimensional grid point field of the initial weather radar net and identifying the boundary of the blind area;

the volume scanning mode processing module is used for collecting data of radars to be added, and configuring volume scanning mode parameters of the radars according to the following steps based on blind area identification for each radar to be added: acquiring longitude and latitude, altitude and radial distance parameters of the radar, and establishing a three-dimensional coordinate system of the radar; judging whether the blind point is in an effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as an effective point when the blind point is in the effective radial range; calculating the azimuth angle and the initial angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range at each height; acquiring a minimum azimuth angle, a maximum azimuth angle and a scanning range among set heights according to the set height range, and then fusing to obtain final fan-scanning parameters; and after the final sector sweep parameter is obtained, acquiring the pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain a final pitch angle parameter.

The lattice point identification module comprises a lattice point field construction unit, an information identification unit covered on an equal-height surface and an information identification unit covered between equal-height layers.

The lattice field construction unit is configured to perform the steps of:

acquiring position data and volume scanning mode parameter data of a plurality of networking weather radars in a target area, wherein the position data comprises radar longitude and latitude and feed source altitude;

for each networking weather radar, converting the polar coordinates into a lattice field under Cartesian coordinates according to the position data and the volume scanning mode parameter data of the networking weather radar, and calculating the longitude range and the latitude range detected by each networking weather radar according to the detection radius of each networking weather radar;

constructing an initial weather radar mesh three-dimensional lattice field by combining the lattice fields of all the networking weather radars, and determining the horizontal range, the horizontal resolution, the height range, the equidistant height hs and the spatial resolution hr of the three-dimensional lattice field; the horizontal range of the three-dimensional lattice field of the initial weather radar net is determined by respectively comparing the longitude range and the latitude range detected by all the networking weather radars, wherein the horizontal range comprises a longitude interval range and a latitude interval range, the longitude interval range is the sum of the distance value between the radar with the largest detection longitude and the radar with the smallest detection longitude, the detection radius of the radar with the largest detection longitude and the detection radius of the radar with the smallest detection longitude, and the latitude interval range is the sum of the distance value between the radar with the largest detection latitude and the radar with the smallest detection latitude, the detection radius of the radar with the largest detection latitude and the detection radius of the radar with the smallest detection latitude; the horizontal resolution, height range, equidistant height hs, and spatial resolution hr are set by the user or the system.

The overlay information identification unit on the equal altitude surface is configured to perform the steps of:

lattice point conversion: according to the set height range, all grid points on the equidistant height hs on the set height Z are obtained, and any grid point P in the grid points is_tObtaining a grid point P_tCoordinate (X) of_t,Y_t,Z_t) Wherein T =1,2, … …, T, T is the total number of lattice points, X_tIs latitude, Y_tIs longitude, Z_tIs the height; determination of grid point P by radar beam propagation and large circle geometry theory_tPolar position (r) in a radar polar coordinate system_t,θ_t,φ_t) Wherein r is_tIs the pitch, θ_tIs an azimuth angle phi_tIs a pitch angle;

identification: for the aforementioned arbitrary lattice point P_tOne by one will (r)_t,θ_t,φ_t) Matching with azimuth, elevation, beam width and distance information of single radar, and determining (r)_t,θ_t,φ_t) When the grid point P is in the detection range_tIdentifying as a radar coverage point, otherwise, identifying the grid point as a blind point; after all radars in the traversing weather radar network identify the grid points according to the method, the radar coverage quantity of each grid point is obtained, when one or more radars cover the grid point, the grid point is identified as a radar coverage point, and when no radar covers the grid point, the grid point is identified as a blind point.

The equal-height interlayer coverage information identification unit is configured to perform the steps of:

based on a single radar, all grid points on the space resolution hr at the set height Z are obtained according to the set height range, and any grid point P in the grid points is_jObtaining a grid point P_jCoordinate (X) of_j,Y_j,Z_j) And coordinate conversion is carried out to obtain the corresponding polar coordinate position (r)_j,θ_j,φ_j) Wherein j =1,2, … …, m, m is the total number of grid points; will (r)_j,θ_j,φ_j) Matching the azimuth, elevation, beam width and distance information of the single radar, and determining (r)_j,θ_j,φ_j) When the grid point P is in the detection range_jIdentifying as a radar coverage point, otherwise, identifying the grid point as a blind point;

regarding a height interval between any two adjacent equidistant heights as a layer Lq, wherein the Lq represents a q-th layer, and the value of q is an integer greater than or equal to 1; matching grid points Pt on the equal distance height hs with grid points Pj of the spatial resolution hr according to the set height Z, and marking the grid points with the same height as repeated grid points; acquiring all layers on a set height Z, calculating the number of lattice points corresponding to each layer Lq, wherein the number of lattice points is equal to the number of lattice points on an equidistant height hs corresponding to the layer and the number of lattice points on a corresponding spatial resolution hr after duplication removal, judging whether the lattice points corresponding to the layer Lq are blind points, if so, marking the layer Lq as a blind point, otherwise, marking the layer Lq as a radar coverage point as long as any one lattice point is a radar coverage point;

and after all radars in the traversing weather radar network identify the layers according to the method, acquiring the radar coverage quantity of each layer, identifying the radar coverage points when one or more radars cover the layer, and identifying the blind points when no radar covers the layer.

Other technical features are described in the previous embodiment and are not described in detail herein.

In the foregoing description, the disclosure of the present invention is not intended to limit itself to these aspects. Rather, the various components may be selectively and operatively combined in any number within the intended scope of the present disclosure. In addition, terms like "comprising," "including," and "having" should be interpreted as inclusive or open-ended, rather than exclusive or closed-ended, by default, unless explicitly defined to the contrary. All technical, scientific, or other terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs unless defined otherwise. Common terms found in dictionaries should not be interpreted too ideally or too realistically in the context of related art documents unless the present disclosure expressly limits them to that. Any changes and modifications of the present invention based on the above disclosure will be within the scope of the appended claims.

Claims (10)

1. A volume scanning mode calculation method based on a networking weather radar coverage area is characterized by comprising the following steps:

constructing an initial weather radar mesh three-dimensional grid field based on networking weather radar data in a target area, and identifying grid point coverage information in the initial weather radar mesh three-dimensional grid field, wherein the grid point coverage information comprises blind point identification and radar coverage point identification;

obtaining blind spot information of marks in a three-dimensional grid spot field of an initial weather radar net, determining a blind area region, and identifying the boundary of the blind area region;

collecting data of radars to be added, identifying and configuring volume scanning mode parameters of the radars to be added based on a blind area for each radar to be added, and comprising the following steps of: acquiring longitude and latitude, altitude and radial distance parameters of the radar, and establishing a three-dimensional coordinate system of the radar; judging whether the blind point is in an effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as an effective point when the blind point is in the effective radial range; calculating the azimuth angle and the initial angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range at each height; and according to the set height range, obtaining the minimum azimuth angle, the maximum azimuth angle and the scanning range among the set heights, and then fusing to obtain the final fan-scanning parameters.

2. The method of claim 1, wherein: after the final fan-scan parameters are obtained, the following steps are performed,

and acquiring a pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain a final pitch angle parameter.

3. The method of claim 2, wherein: the step of constructing the initial weather radar mesh three-dimensional grid field comprises,

acquiring position data and volume scanning mode parameter data of a plurality of networking weather radars in a target area, wherein the position data comprises radar longitude and latitude and feed source altitude;

for each networking weather radar, converting the polar coordinates into a lattice field under Cartesian coordinates according to the position data and the volume scanning mode parameter data of the networking weather radar, and calculating the longitude range and the latitude range detected by each networking weather radar according to the detection radius of each networking weather radar;

constructing an initial weather radar mesh three-dimensional lattice field by combining the lattice fields of all the networking weather radars, and determining the horizontal range, the horizontal resolution, the height range, the equidistant height hs and the spatial resolution hr of the three-dimensional lattice field; the horizontal range of the three-dimensional lattice field of the initial weather radar net is determined by respectively comparing the longitude range and the latitude range detected by all the networking weather radars, wherein the horizontal range comprises a longitude interval range and a latitude interval range, the longitude interval range is the sum of the distance value between the radar with the largest detection longitude and the radar with the smallest detection longitude, the detection radius of the radar with the largest detection longitude and the detection radius of the radar with the smallest detection longitude, and the latitude interval range is the sum of the distance value between the radar with the largest detection latitude and the radar with the smallest detection latitude, the detection radius of the radar with the largest detection latitude and the detection radius of the radar with the smallest detection latitude; the horizontal resolution, height range, equidistant height hs, and spatial resolution hr are set by the user or the system.

4. The method of claim 3, wherein: the grid point coverage information in the three-dimensional grid point field of the initial weather radar mesh is identified by analyzing the coverage information on the equal altitude surface, the steps are as follows,

lattice point conversion: according to the set height range, all grid points on the equidistant height hs on the set height Z are obtained, and any grid point P in the grid points is_tObtaining a grid point P_tCoordinate (X) of_t,Y_t,Z_t) Wherein T =1,2, … …, T, T is the total number of lattice points, X_tIs latitude, Y_tIs longitude, Z_tIs the height; determination of grid point P by radar beam propagation and large circle geometry theory_tPolar position (r) in a radar polar coordinate system_t,θ_t,φ_t) Wherein r is_tIs the pitch, θ_tIs an azimuth angle phi_tIs a pitch angle;

identification: for the aforementioned arbitrary lattice point P_tOne by one will (r)_t,θ_t,φ_t) Matching with azimuth, elevation, beam width and distance information of single radar, and determining (r)_t,θ_t,φ_t) When the grid point P is in the detection range_tIdentifying as a radar coverage point, otherwise, identifying the grid point as a blind point; after all radars in the traversing weather radar network identify the grid points according to the method, the radar coverage quantity of each grid point is obtained, when one or more radars cover the grid point, the grid point is identified as a radar coverage point, and when no radar covers the grid point, the grid point is identified as a blind point.

5. The method of claim 4, wherein: combining the coverage information among the equal altitude layers to identify the lattice point coverage information in the initial weather radar mesh three-dimensional lattice point field, and the method comprises the following steps:

based on a single radar, all grid points on the space resolution hr at the set height Z are obtained according to the set height range, and any grid point P in the grid points is_jObtaining a grid point P_jCoordinate (X) of_j,Y_j,Z_j) And coordinate conversion is carried out to obtain the corresponding polar coordinate position (r)_j,θ_j,φ_j) Wherein j =1,2, … …, m, m is the total number of grid points; will (r)_j,θ_j,φ_j) Matching the azimuth, elevation, beam width and distance information of the single radar, and determining (r)_j,θ_j,φ_j) When the grid point P is in the detection range_jIdentifying as a radar coverage point, otherwise, identifying the grid point as a blind point;

regarding the height interval between any two adjacent equidistant heights as a layer L_q，L_qRepresenting a q-th layer, wherein the value of q is an integer which is more than or equal to 1; according to the set height Z, the lattice point P on the equal distance height hs_tLattice point P with spatial resolution hr_jMatching, and marking the grid points with the same height as repeated grid points; and acquiring all layers at a set height Z, for each layerL_qCalculate the layer L_qThe corresponding number of lattice points, which is equal to the number of the lattice points on the equidistant height hs corresponding to the layer and the number of the lattice points on the corresponding spatial resolution hr after the duplication is removed, the layer L is judged_qWhether the corresponding grid points are all blind points or not, and if so, the layer L is divided into the blind points_qMarking as a blind spot, otherwise, if any grid point is a radar coverage point, the layer L is marked as a blind spot_qMarking as a radar coverage point;

and after all radars in the traversing weather radar network identify the layers according to the method, acquiring the radar coverage quantity of each layer, identifying the radar coverage points when one or more radars cover the layer, and identifying the blind points when no radar covers the layer.

6. The method of claim 5, wherein: and identifying blind points by using a number 0, and when identifying radar coverage points, identifying corresponding numbers k according to the radar coverage quantity on the grid points, wherein k =1,2, … …, N and N are the total number of radars, and the radar coverage points comprise single radar coverage points, double radar coverage punctuations, … … and k radar coverage points.

7. The method according to any one of claims 1-6, wherein: the radar to be added is an X-band radar.

8. A collaborative adaptive scanning method for a networking weather radar is characterized by comprising the following steps:

acquiring radar data to be added, which can be added into a target area weather radar network, in a preset area;

the method according to any of claims 1-7, configuring volume scan mode parameters to be added to the radar, developing observations.

9. A volume scanning pattern calculation apparatus based on a networked weather radar coverage area, comprising:

the grid point identification module is used for constructing an initial weather radar grid three-dimensional grid point field according to networking weather radar data in a target area, and identifying grid point coverage information in the initial weather radar grid three-dimensional grid point field, wherein the grid point coverage information comprises blind point identification and radar coverage point identification;

the blind area identification module is used for determining a blind area after acquiring the blind point information of the mark in the three-dimensional grid point field of the initial weather radar network and identifying the boundary of the blind area;

the volume scanning mode processing module is used for collecting data of the radar to be added, and configuring volume scanning mode parameters of the radar according to the following steps based on blind area identification for each radar to be added:

acquiring longitude and latitude, altitude and radial distance parameters of the radar, and establishing a three-dimensional coordinate system of the radar; judging whether the blind point is in an effective radial range according to the distance information between the blind point of the blind area boundary and the radar center, and identifying the blind point as an effective point when the blind point is in the effective radial range; calculating the azimuth angle and the initial angle of each effective point, and acquiring the minimum azimuth angle, the maximum azimuth angle and the scanning range at each height; and according to the set height range, obtaining the minimum azimuth angle, the maximum azimuth angle and the scanning range among the set heights, and then fusing to obtain the final fan-scanning parameters.

10. The apparatus of claim 9, wherein: the volume scan mode processing module is further configured to,

and after the final sector scanning parameters are obtained, acquiring the pitch angle of each effective point, sequencing the pitch angles of the effective points in the set height range to obtain a minimum pitch angle and a maximum pitch angle, and combining the minimum pitch angle, the maximum pitch angle and the set beam spread value to obtain the final pitch angle parameters.

Images