Antiaircraft gun hail suppression operation point firing range diagram generation method and system and medium

Abstract

The invention relates to an antiaircraft gun hail suppression operation point firing range map generation method and system and a medium, and the method comprises the steps: obtaining multi-element geographic data through aerial survey of an unmanned aerial vehicle, and generating a data base map; determining an air density parameter based on the altitude of the operation point extracted from the multi-element geographic data, and correcting the trajectory motion equation of the antiaircraft gun by using the air density parameter so as to solve the range and the height of the antiaircraft gun at each shooting angle; constructing a working period wind field dynamic model, and carrying out ballistic trajectory offset correction on an antiaircraft gun ballistic motion equation by utilizing a wind field vector simulated by the wind field dynamic model so as to solve ranges and azimuth angles in different working periods; and fusing the firing range and the firing height under each firing angle and the firing range and the azimuth angle of each operation period, generating an airspace firing boundary on the data base map, and outputting a two-dimensional safe firing boundary map. Therefore, through geographic data fusion and dynamic trajectory correction, the accuracy and safety of the antiaircraft gun hail suppression work firing range map are remarkably improved.

Description

Technical Field

[0001] This invention relates to the field of meteorological safety operation technology, and in particular to a method, system and medium for generating a hail suppression point boundary map for anti-aircraft artillery. Background Technology

[0002] Artificial hail suppression is a key means of agricultural meteorological disaster prevention and mitigation, and is crucial for ensuring food security. Safety management is the core element for the healthy development and risk prevention of weather modification. In the safety management system for artificial anti-aircraft gun hail suppression operations, the safety firing boundary map occupies a fundamental and critical position. As a "safety navigation map" for the operation site, it accurately delineates the safety firing boundary zone, which is the most effective technical measure to avoid operational safety accidents, reflecting the progress of safety management from extensive to intensive. Initially, the safety firing boundary map was drawn based on the principle requirements of the "Regulations on the Management of Weather Modification." With the release of the meteorological industry standard QX / T 256-2015 "Specification for Drawing Safety Firing Boundary Maps of 37mm Anti-aircraft Gun Weather Modification Operation Sites," it gained preliminary technical support. In recent years, it has been further upgraded to the national standard GB / T 39782-2021 "Specification for Drawing Safety Firing Boundary Maps of Rocket Weather Modification Operation Sites," marking the establishment of a unified national-level technical specification in this field.

[0003] However, despite the aforementioned standards and functions, the current creation and application of artificial hail suppression safety boundary maps still suffer from a series of shortcomings, limiting their safety assurance effectiveness. At the data and technology level, because many operational sites rely on outdated maps or low-resolution satellite imagery as base maps, information on terrain features (such as newly built residential areas) within the hail suppression boundary is inaccurate. Furthermore, the widespread adoption of automated mapping systems means that many areas still heavily rely on traditional manual drawing methods, resulting in low efficiency, susceptibility to human error, and difficulty in ensuring the standardization and scientific accuracy of the maps. At the standardization and application level, this is manifested in the ineffective implementation of existing national and industry standards. Some hail suppression maps do not strictly adhere to the prescribed safe hail suppression zone boundaries, and the degree of standardization varies considerably. More importantly, because hail suppression maps have not been effectively integrated into modern operational command systems (such as integrated meteorological monitoring, real-time airspace reporting, and automated operational command systems), they cannot provide dynamic, real-time safety boundary guidance during rapidly changing operational decision-making processes, thus limiting their practical value. At the dynamic management level, since safety boundary maps are usually regarded as static administrative approval documents and lack a mandatory periodic review and dynamic update mechanism, they are difficult to respond to the rapid changes in the surrounding environment in a timely manner. Summary of the Invention

[0004] (a) Technical problems to be solved

[0005] In view of the above-mentioned shortcomings and deficiencies of the prior art, the present invention provides a method, system and medium for generating a hail suppression operation point firing boundary map for anti-aircraft guns, which solves the technical problems of insufficient base map accuracy, reliance on manual drawing, standard execution deviation, disconnection from the operation command system and lack of real-time updates in the existing manual hail suppression safety firing boundary map.

[0006] (II) Technical Solution

[0007] To achieve the above objectives, the main technical solutions adopted by the present invention include:

[0008] In a first aspect, embodiments of the present invention provide a method for generating a point firing boundary map for anti-aircraft artillery hail suppression operations, comprising:

[0009] The aerial survey was conducted using drones to obtain multi-element geographic data of anti-hail gun operation points and preset areas, and a data base map was generated.

[0010] The air density parameter is determined based on the altitude of the work site extracted from multi-element geographic data, and the air density parameter is used to perform altitude-adaptive correction on the pre-constructed anti-aircraft gun trajectory equation in order to calculate the range and altitude at each firing angle.

[0011] Based on the acquired historical meteorological observation data and meteorological reanalysis data, a wind field dynamic model for the operation period is constructed. The wind field vector simulated by the wind field dynamic model is used to correct the trajectory deviation of the anti-aircraft gun ballistic motion equation, so as to calculate the range and azimuth of different operation periods.

[0012] By integrating the range and altitude at various firing angles and the range and azimuth at various operational periods, an airspace firing boundary is generated on the data base map, and a safe firing boundary map for anti-aircraft artillery hail suppression operations is output.

[0013] Optionally, aerial surveying using drones can be used to acquire multi-element geographic data of the anti-aircraft gun hail suppression operation point and a preset range, and a data base map can be generated, including:

[0014] Based on the coordinates of the currently pre-deployed and historically deployed hail suppression operation points, the drone equipped with lidar, photography module and multispectral sensor is controlled to first perform lidar scanning to generate an elevation reference layer according to the preset radial route, and then switch to the preset spiral progressive route to perform aerial surveying tasks.

[0015] During aerial surveying, the elevation data of the lidar is compared with the elevation reference layer in real time. When a terrain change zone with an elevation difference exceeding a preset threshold is detected, the photography module is triggered to perform multi-angle image verification of the terrain change zone. When the verification is confirmed, the scanning frequency of the lidar is adjusted in coordination to generate a geo-registered enhanced elevation layer.

[0016] Based on the spectral data of ground features collected by multispectral sensors and the distribution characteristics of the enhanced elevation layer, vegetation mask and artificial building outlines are extracted through a pre-constructed intelligent segmentation model. Combined with historical topographic benchmark data, dynamic ground feature interference is eliminated through comparison and iteration, and a clear topographic map containing static ground feature vectors and elevation information is output.

[0017] Spatially fuse the topographic map of the clear area with the acquired imagery, combine it with artificial buildings including houses, villages and factories based on the outline of artificial buildings, and inject plane coordinate system parameters and time labels to generate a two-dimensional data base map including topographic elevation, static feature distribution and image base.

[0018] Optionally, air density parameters are determined based on the altitude of the operational site extracted from multi-element geographic data, and the pre-constructed anti-aircraft gun trajectory equations are modified for altitude adaptability using air density parameters to calculate the range and altitude at various firing angles, including:

[0019] The altitude of the work site is extracted from multi-element geographic data. Combined with the obtained vertical lapse rate of atmospheric temperature and predefined sea level reference parameters, and based on the correlation between gas molecule properties and thermodynamic constants, the air density corresponding to the altitude is calculated.

[0020] The nonlinear relationship of the drag coefficient is determined by the ratio of the instantaneous velocity of the projectile to the speed of sound in the air medium. The reciprocal of the ballistic coefficient is calculated by combining the inherent properties of the projectile and the air density.

[0021] Based on the horizontal influence term that couples the projectile velocity and horizontal velocity components using the reciprocal of the ballistic coefficient as a scaling factor, the rate of change of the horizontal velocity in the trajectory equation of the anti-aircraft gun is dynamically adjusted.

[0022] Based on the vertical velocity attenuation benchmark caused by gravitational acceleration, a product term of projectile velocity and vertical velocity components is introduced, with the reciprocal of the ballistic coefficient as a scaling factor, to construct a vertical influence term of drag and gravity. Based on the vertical influence term, the rate of change of vertical velocity in the trajectory equation of the anti-aircraft gun is dynamically adjusted.

[0023] The time-step numerical iterative calculation is applied to the first revised anti-aircraft gun ballistic equation to update the instantaneous motion state of the projectile in the horizontal and vertical directions one by one until the conditions for the projectile to hit the ground are met, and the calculation process is terminated. Finally, the range and altitude corresponding to each firing angle are output.

[0024] Among them, the ballistic motion equation is a dynamic model composed of a set of velocity differential equations in the horizontal and vertical directions and a set of position update equations. It constructs a complete description of the projectile's trajectory by using the inverse of the ballistic coefficient to quantify the nonlinear attenuation effect of air resistance on the velocity component and superimposing the effect of gravity.

[0025] Optionally, the expression for the first revised anti-aircraft gun trajectory equation is:

[0026] ;

[0027] In the formula, The reciprocal of the ballistic coefficient, ρ is the air density, and c is the air density. d Where A is the drag coefficient, v is the projectile cross-sectional area, m is the projectile velocity, and v is the projectile mass. x For the horizontal velocity component, v y denoted as the vertical velocity component, and g as the acceleration due to gravity.

[0028] Optionally, a dynamic wind field model for the operational period is constructed based on historical meteorological observation data and meteorological reanalysis data. The wind field vector simulated by the dynamic wind field model is used to correct the trajectory deviation of the anti-aircraft gun's ballistic motion equation, so as to calculate the range and azimuth angle for different operational periods, including:

[0029] By integrating historical meteorological observation data and meteorological reanalysis data of the acquired work site and surrounding preset range, a dynamic wind field model for the work period is constructed.

[0030] Based on the wind field dynamic model, the temporal changes of the wind field at the work site are numerically simulated under various meteorological modes. The wind field vector is obtained and decomposed to the ballistic coordinate system to obtain the longitudinal wind component along the projectile range and the crosswind component perpendicular to the firing direction.

[0031] In the longitudinal and transverse wind dimensions, the velocity vector of the projectile relative to the air is quantified based on the wind speed vector and the measured velocity of the projectile relative to the ground.

[0032] In the first revised anti-aircraft gun trajectory equation, the projectile velocity is replaced by the velocity vector of the projectile relative to the air, and the rate of velocity change in the range direction caused by the longitudinal wind and the rate of velocity change in the perpendicular to the range direction caused by the crosswind are calculated respectively.

[0033] The displacement correction in the range direction is calculated by time-series integration along the trajectory based on the rate of change of velocity in the range direction, and the drift displacement perpendicular to the range direction is solved by time-series integration along the trajectory based on the rate of change of velocity perpendicular to the range direction.

[0034] The range displacement correction and the lateral drift displacement are spatially vector-superimposed to generate corrected ballistic trajectory parameters that take into account the wind field coupling effect.

[0035] Based on the modified ballistic trajectory parameters that take into account the wind field coupling effect, the range and azimuth angle corresponding to different operational periods are calculated.

[0036] Optionally, by integrating the range and elevation at various firing angles and the range and azimuth at various operational periods, an airspace firing boundary is generated on the data base map, outputting a safe firing boundary map for anti-aircraft gun hail suppression operations, including:

[0037] By integrating the range and altitude at various firing angles and the range and azimuth at various operational phases, a ballistic trajectory parameter set and a parameter space composed of altitude, range, and azimuth are obtained.

[0038] The ballistic parameter set and the data base map are registered in coordinate system, and airspace areas that conflict with terrain obstacles are dynamically eliminated based on the maximum firing height threshold of the projectile constrained by ballistic performance.

[0039] The accessibility analysis of the airspace region after culling operation is performed on an azimuth-by-azimuth basis within the parameter space composed of firing height, range and azimuth angle. Discrete firing envelopes are generated. Multi-angle firing envelope interpolation is performed on the discrete firing envelopes to generate a continuous and smooth airspace firing boundary surface.

[0040] By overlaying the airspace firing boundary surface onto the data base map and combining it with the elevation data for layered rendering, a two-dimensional safe firing boundary map for anti-aircraft artillery hail suppression operations is generated.

[0041] Secondly, embodiments of the present invention provide a system for generating a point firing boundary map for anti-aircraft gun hail suppression operations, comprising:

[0042] The basic database subsystem is used to store operational site information data, terrain data, satellite remote sensing data, UAV aerial photography data, meteorological observation data, and meteorological reanalysis data.

[0043] The data analysis subsystem is used to execute the methods described above;

[0044] The radiation boundary map creation subsystem is used to draw a safety radiation boundary map with a preset accuracy based on the data output by the basic database subsystem and the data analysis subsystem.

[0045] The boundary map management subsystem is used to respond to commands or detect status changes, and to provide boundary map review reminders and update alarms.

[0046] Optionally, the data analysis subsystem includes:

[0047] The data acquisition module is used to acquire multi-element geographic data of anti-aircraft gun hail suppression operation points and preset ranges through UAV aerial surveying, and generate a data base map;

[0048] The first analysis module is used to determine the air density parameter based on the altitude of the work site extracted from multi-element geographic data, and to use the air density parameter to perform altitude-adaptive correction on the pre-constructed anti-aircraft gun trajectory equation in order to calculate the range and altitude at each firing angle.

[0049] The second analysis module is used to construct a wind field dynamic model for the operation period based on the acquired historical meteorological observation data and meteorological reanalysis data, and to use the wind field vector simulated by the wind field dynamic model to correct the trajectory deviation of the anti-aircraft gun ballistic motion equation, so as to calculate the range and azimuth of different operation periods.

[0050] The firing range map output module is used to integrate the range and altitude at various firing angles and the range and azimuth at various operational periods to generate airspace firing range boundaries on the data base map and output a safe firing range map for anti-aircraft artillery hail suppression operations.

[0051] Thirdly, embodiments of the present invention provide a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the method described above.

[0052] (III) Beneficial Effects

[0053] The beneficial effects of this invention are:

[0054] First, the high-resolution data base map generated by the present invention based on drone aerial photography can effectively identify the surrounding environment of the anti-aircraft gun operation point, capture high-definition environmental base map data, avoid the firing boundary error caused by insufficient resolution of traditional satellite data, and ensure the authenticity of the geographic information of the artificial hail suppression operation point and its surrounding environment.

[0055] Secondly, the air density correction mechanism based on the actual altitude of the operation point breaks through the limitation of existing technologies that do not consider altitude differences. It specifically corrects the air resistance and ballistic coefficient in the ballistic equation, overcomes the problem of range calculation deviation in high-altitude areas caused by the use of general formulas in existing methods, and ensures that the range and altitude calculations at different firing angles conform to the actual terrain constraints.

[0056] Meanwhile, the ballistic deviation correction driven by the wind field dynamic model uses historical meteorological data and reanalysis data to simulate the impact of the wind field on the ballistic trajectory during the operation period. It incorporates the combined effects of longitudinal and crosswinds on the trajectory into the equation correction, realizing real-time deviation compensation for range and azimuth angle during different operation periods. This overcomes the risk of outdated safe firing range caused by ignoring wind field changes in traditional methods.

[0057] Finally, the corrected ballistic parameters are fused with the geographic base map to output a digital safe firing range map. This avoids subjective errors from manual drawing and is compatible with the data interface of the operational command system, supporting real-time updates and dynamic adjustments. This closed-loop technology systematically solves the multi-dimensional shortcomings of traditional methods, such as base map accuracy, reliance on manual labor, standard execution, system collaboration, and update timeliness. Attached Figure Description

[0058] Figure 1This is a schematic diagram of the overall process of the method provided in the embodiments of the present invention;

[0059] Figure 2 This is a schematic diagram illustrating the specific process of step S1 of the method provided in this embodiment of the invention;

[0060] Figure 3 This is a detailed flowchart illustrating step S2 of the method provided in this embodiment of the invention;

[0061] Figure 4 This is a detailed flowchart illustrating step S3 of the method provided in this embodiment of the invention;

[0062] Figure 5 This is a detailed flowchart illustrating step S4 of the method provided in this embodiment of the invention;

[0063] Figure 6 This is a schematic diagram of the system provided in an embodiment of the present invention. Detailed Implementation

[0064] To better explain and facilitate understanding of the present invention, the present invention will be described in detail below with reference to the accompanying drawings and specific embodiments.

[0065] like Figure 1 As shown in the figure, an embodiment of the present invention proposes a method for generating a firing boundary map for anti-aircraft artillery hail suppression operations. The method includes: acquiring multi-element geographic data of the anti-aircraft artillery hail suppression operation point and a preset range through UAV aerial surveying, and generating a data base map; determining air density parameters based on the altitude of the operation point extracted from the multi-element geographic data, and using the air density parameters to perform altitude-adaptive correction on the pre-constructed anti-aircraft artillery ballistic motion equation to calculate the range and altitude at various firing angles; constructing a wind field dynamic model for the operation period based on acquired historical meteorological observation data and meteorological reanalysis data, and using the wind field vector simulated by the wind field dynamic model to perform trajectory offset correction on the anti-aircraft artillery ballistic motion equation to calculate the range and azimuth at different operation periods; integrating the range and altitude at various firing angles and the range and azimuth at various operation periods to generate an airspace firing boundary on the data base map, and outputting a safe firing boundary map for anti-aircraft artillery hail suppression operations.

[0066] First, the high-resolution data base map generated by the present invention based on drone aerial photography can effectively identify the surrounding environment of the anti-aircraft gun operation point, capture high-definition environmental base map data, avoid the firing boundary error caused by insufficient resolution of traditional satellite data, and ensure the authenticity of the geographic information of the artificial hail suppression operation point and its surrounding environment.

[0067] Secondly, the air density correction mechanism based on the actual altitude of the operation point breaks through the limitation of existing technologies that do not consider altitude differences. It specifically corrects the air resistance and ballistic coefficient in the ballistic equation, overcomes the problem of range calculation deviation in high-altitude areas caused by the use of general formulas in existing methods, and ensures that the range and altitude calculations at different firing angles conform to the actual terrain constraints.

[0068] Meanwhile, the ballistic deviation correction driven by the wind field dynamic model uses historical meteorological data and reanalysis data to simulate the impact of the wind field on the ballistic trajectory during the operation period. It incorporates the combined effects of longitudinal and crosswinds on the trajectory into the equation correction, realizing real-time deviation compensation for range and azimuth angle during different operation periods. This overcomes the risk of outdated safe firing range caused by ignoring wind field changes in traditional methods.

[0069] Finally, the corrected ballistic parameters are fused with the geographic base map to output a digital safe firing range map. This avoids subjective errors from manual drawing and is compatible with the data interface of the operational command system, supporting real-time updates and dynamic adjustments. This closed-loop technology systematically solves the multi-dimensional shortcomings of traditional methods, such as base map accuracy, reliance on manual labor, standard execution, system collaboration, and update timeliness.

[0070] To better understand the above technical solutions, exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. Although exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that the present invention can be understood more clearly and thoroughly, and that the scope of the present invention can be fully conveyed to those skilled in the art.

[0071] Specifically, embodiments of the present invention provide a method for generating a point firing boundary map for anti-aircraft artillery hail suppression operations, comprising:

[0072] S1. Obtain multi-element geographic data of the anti-aircraft artillery hail suppression operation point and its preset range through UAV aerial surveying, and generate a data base map. This invention is based on UAV aerial photography technology, using a civilian multi-rotor UAV to complete aerial photography of the anti-aircraft artillery hail suppression operation point and its surroundings. Combined with deployed UAV data processing software, high-definition map drawing of the area around the operation point is achieved, with a map resolution of less than 10 centimeters.

[0073] Furthermore, such as Figure 2 As shown, step S1 includes:

[0074] S11. Based on the coordinates of the currently pre-deployed and historically deployed hail suppression operation points, control the UAV equipped with lidar, photography module and multispectral sensor to first perform lidar scanning to generate an elevation reference layer according to the preset radial route, and then switch to the preset spiral progressive route to perform aerial surveying tasks.

[0075] Specifically, the UAV is controlled to fly along a pre-defined radial flight path (referring to a fan-shaped path radiating outwards from the work point) based on the coordinates of pre-deployed work points (not only pre-deployed but also already deployed work points. As required, the work point design map is updated periodically (generally every 5 years) or irregularly (when there are significant changes around the work point)). It utilizes lidar to rapidly scan the surface, generating a two-dimensional initial elevation reference layer (DEM) containing terrain elevation and slope distribution characteristics for aerial survey parameter calibration. After parameter calibration, the UAV switches to a spiral progressive flight path (referring to the UAV starting from the outer boundary of the work point coverage area, spiraling inwards along a continuous path, gradually reducing the spiral radius until the central area is covered), performing the main aerial survey task. Simultaneously, it collects lidar point clouds, multi-angle images, and multispectral data, ensuring a balance between high-resolution data acquisition and full-area coverage efficiency for the target terrain area.

[0076] S12. During aerial surveying, the elevation data of the lidar is compared with the elevation reference layer in real time. When a terrain abrupt change area with an elevation difference exceeding a preset threshold is detected, the photography module is triggered to perform multi-angle image verification of the terrain abrupt change area. When the verification is confirmed, the scanning frequency of the lidar is adjusted collaboratively to generate a geo-registered enhanced elevation layer. In this step, during aerial surveying, the current lidar elevation data is compared with the pre-generated elevation reference layer in real time to identify terrain abrupt change areas with elevation differences exceeding a preset threshold (e.g., ±0.5 meters). For abnormal areas, the local high-resolution multi-angle image re-measurement of the photography module is triggered. The authenticity of the terrain abrupt change (e.g., geological collapse or temporary obstacles) is verified by image texture analysis and comparison with historical benchmarks. If the abrupt change is verified, the local scanning frequency of the lidar is dynamically increased (e.g., from 100Hz to 200Hz) to generate a geo-registered enhanced DEM layer (resolution increased from 10cm to 5cm).

[0077] S13. Based on the spectral data of ground features acquired by multispectral sensors and the distribution characteristics of the enhanced elevation layer, a pre-constructed intelligent segmentation model is used to extract vegetation masks and artificial building outlines. This is combined with historical topographic benchmark data for comparison and iterative elimination of dynamic ground feature interference, outputting a clear topographic map containing static ground feature vectors and elevation information. In the data processing stage, the spectral data of ground features acquired by multispectral sensors and the topographic distribution characteristics of the enhanced elevation layer are first spatially aligned and feature-level fused to form a multi-channel input dataset. A deep learning-based U-Net segmentation model is pre-trained for intelligent ground feature segmentation. This model extracts spectral-topographic joint features through an encoder, while the decoder gradually restores spatial details and generates pixel-level classification results: the Normalized Difference Vegetation Index (NDVI) is calculated using the near-infrared and red bands in the ground feature spectral data, and a dynamic threshold is set to segment vegetation mask areas; simultaneously, based on the abrupt changes in reflectance in the short-wave infrared bands and the detection of local elevation extrema in the ground feature spectral data, the outlines of regular buildings such as houses and factories are extracted. Subsequently, historical topographic baseline data was retrieved for spatiotemporal difference analysis. A sliding window was used to traverse and compare the spatial consistency between the current building outline and historical vector data. If a building feature exists in the current data but has no corresponding entry in the historical record, and the duration does not exceed a preset threshold (e.g., 45 days), it is determined to be a temporary structure and removed. After multiple rounds of iterative verification and manual outline optimization, the final output is a clear topographic map containing permanent building vector boundaries, vegetation distribution range, and high-precision topographic elevation information.

[0078] S14. Spatially fuse the topographic map with the acquired images, combine the artificial buildings including houses, villages and factories based on the outline annotation of artificial buildings, and inject plane coordinate system parameters and time labels to generate a two-dimensional data base map including topographic elevation, static feature distribution and image base.

[0079] S2. Based on the altitude of the work site extracted from multi-element geographic data, the air density parameter is determined, and the air density parameter is used to perform altitude-adaptive correction on the pre-constructed anti-aircraft gun trajectory equation in order to calculate the range and altitude at each firing angle.

[0080] Furthermore, such as Figure 3 As shown, step S2 includes:

[0081] S21. Extract the altitude of the work site from the multi-element geographic data, combine the obtained vertical lapse rate of atmospheric temperature and the predefined sea level reference parameters, and calculate the air density corresponding to the altitude based on the correlation between gas molecule properties and thermodynamic constants.

[0082] S22. Determine the nonlinear relationship of the drag coefficient based on the ratio of the instantaneous velocity of the projectile to the speed of sound in the air medium. Combine the inherent properties of the projectile with the air density to calculate the reciprocal of the ballistic coefficient.

[0083] S23. Based on the horizontal influence term of the projectile velocity and horizontal velocity components coupled with the reciprocal of the ballistic coefficient as a scaling factor, the rate of change of the horizontal velocity in the trajectory equation of the anti-aircraft gun is dynamically adjusted.

[0084] S24. Based on the vertical velocity attenuation benchmark caused by gravitational acceleration, a product term of projectile velocity and vertical velocity components is introduced, with the reciprocal of the ballistic coefficient as a scaling factor, to construct a vertical influence term of drag and gravity. Based on the vertical influence term, the rate of change of vertical velocity in the trajectory equation of the anti-aircraft gun is dynamically adjusted.

[0085] S25. Apply time-step numerical iterative calculation to the first revised anti-aircraft gun trajectory equation, and update the instantaneous motion state of the projectile in the horizontal and vertical directions one by one until the conditions for the projectile to hit the ground are met, and terminate the calculation process. Finally, output the range and altitude corresponding to each firing angle.

[0086] In this step, after the initial correction of the anti-aircraft gun trajectory equations, the numerical iterative solution stage begins. First, a time-step numerical calculation method is used to discretize the continuous anti-aircraft gun trajectory equations into a sequence of state updates within minute time intervals. Within each time step, based on the current horizontal and vertical velocity components of the projectile, its instantaneous velocity and spatial position are updated synchronously.

[0087] Subsequently, the flight status of the projectile is progressively advanced through a loop iteration mechanism. After each iteration, it is checked in real time whether the vertical coordinate of the projectile has reached the ground contact threshold (e.g., the projectile height is lower than the ground reference). If the condition is met, the calculation process is terminated immediately; if not, the parameters for the next time step are loaded and recalculated.

[0088] During this process, parameters such as air density and the reciprocal of the ballistic coefficient are dynamically updated according to the projectile's position, ensuring real-time feedback on the effects of drag. Finally, when the projectile hits the ground, its horizontal displacement is recorded as the range, the peak height of its flight trajectory is recorded as the altitude, and the calculation results for different firing angles are summarized and output.

[0089] It is important to understand that the ballistic motion equation is a dynamic model composed of a set of velocity differential equations in the horizontal and vertical directions and a set of position update equations. It constructs a complete description of the projectile's trajectory by using the inverse of the ballistic coefficient to quantify the nonlinear attenuation effect of air resistance on the velocity component and superimposing the effect of gravity.

[0090] In one specific embodiment, the ballistic calculation process is implemented through a set of differential equations driven by dynamic coupling of air density and the reciprocal of the ballistic coefficients, as well as a set of position update equations.

[0091] First, based on the actual altitude h of the anti-aircraft gun emplacement, the altitude parameters of the operational point are extracted from multi-element geographic data. Combined with the atmospheric temperature lapse rate L (0.0065 K / m, in the troposphere), sea-level standard density ρ0 (1.225 kg / m³), sea-level standard temperature T0 (288.15 K), air molar mass M (0.0289644 kg / mol), universal gas constant R (8.31446 J / (mol•K)), and gravitational acceleration g, the real-time air density is calculated using the following modified atmospheric model:

[0092] ;

[0093] This formula reflects the characteristic that air density decreases exponentially with increasing altitude.

[0094] Secondly, substitute density into the air resistance formula:

[0095] ;

[0096] In the formula, F d Let be the air resistance (N), ρ be the air density (kg / m³), and ρ be a function of altitude h. d The drag coefficient is a function of the Mach number Ma = v / c (velocity divided by the local speed of sound). A is the cross-sectional area of ​​the projectile (m²), and v is the instantaneous velocity of the projectile (m / s).

[0097] Subsequently, the motion equations of the anti-aircraft gun trajectory were modeled, and the force analysis based on Newton's second law was performed. In the horizontal direction, the component of air resistance in the velocity direction can be decomposed into... ,in This is the angle between the instantaneous velocity vector and the horizontal plane (i.e., the trajectory inclination angle). Because... (v) x (where the horizontal velocity component is the equation of motion in the horizontal direction) can be expressed as:

[0098] ;

[0099] For forces acting in the vertical direction, both gravity and air resistance must be considered, specifically:

[0100] ;

[0101] Combination (v) y Given the relationship between the vertical velocity components, the vertical acceleration can be simplified to:

[0102] ;

[0103] Furthermore, the gravity term g dominates the decrease in vertical velocity, while the air resistance term further exacerbates the decrease in velocity.

[0104] The position update equation is directly defined by the velocity components, that is:

[0105] ;

[0106] This pair of differential equations describes how the displacement of the projectile in the horizontal and vertical directions changes over time.

[0107] In summary, the system of equations of motion can be integrated as follows:

[0108] ;

[0109] In the formula, It is the reciprocal of the ballistic coefficient, and its specific calculation process is as follows: Look up the drag coefficient c based on the ratio of the projectile velocity v to the local speed of sound c (Mach number Ma = v / c). d The nonlinear variation of the projectile, combined with the projectile mass m, cross-sectional area A, and current air density ρ, is used to calculate the reciprocal of the ballistic coefficient, which is physically represented as the proportionality factor of air resistance per unit mass.

[0110] S3. Based on the acquired historical meteorological observation data and meteorological reanalysis data, a dynamic wind field model for the operational period is constructed. The wind field vector simulated by the dynamic wind field model is used to correct the trajectory deviation of the anti-aircraft gun's trajectory equation, so as to calculate the range and azimuth for different operational periods. This invention is based on meteorological observation data (surface wind speed and direction, wind speed and direction from radiosonde observations) of the operational site and its surroundings over many years (more than 10 years), combined with meteorological reanalysis data (EAR5, etc.), and uses meteorological models to calculate the wind field changes at the operational site during key operational periods (wind field changes over time; generally, the peak operational period for hail suppression is from May to October each year). Based on the actual wind field data and the trajectory equation, the range and azimuth of the shells are corrected for different operational periods.

[0111] Furthermore, such as Figure 4 As shown, step S3 includes:

[0112] S31. Integrate historical meteorological observation data and meteorological reanalysis data of the acquired operation site and surrounding preset range to construct a wind field dynamic model during the operation period; wherein, the wind field dynamic model is a spatiotemporal evolution model constructed based on multi-source meteorological data fusion and numerical simulation technology, used to characterize the spatiotemporal variation characteristics of the wind field in the hail suppression operation site and surrounding airspace during the operation period.

[0113] To quantify the impact of wind field in ballistic correction, historical meteorological observation data (such as ground station wind speed and radiosonde vertical wind profiles) and meteorological reanalysis data (such as ERA5 and NCEP global reanalysis data) of the operation site and surrounding pre-defined areas are first integrated. Based on data assimilation technology, multi-source heterogeneous data are fused into a unified spatiotemporal grid to construct a spatiotemporal evolution model of the wind field during the operation period. This model is deployed on a high-performance computing cluster, using the MPI parallel computing framework to improve computational efficiency, and outputs wind field vector data to characterize the spatiotemporal variation of wind speed and direction within the hail suppression operation airspace.

[0114] S32. Based on the wind field dynamic model, numerical simulations are performed on the temporal changes of the wind field at the work site under various meteorological modes to obtain the wind field vector and decompose it to the ballistic coordinate system, thereby obtaining the longitudinal wind component along the projectile range direction and the crosswind component perpendicular to the firing direction.

[0115] In this step, batch numerical simulation tasks are initiated for various typical meteorological models (such as calm winds, strong convection, and stratospheric disturbances). A Python automated script calls the wind field model kernel to generate a dataset of wind field temporal changes for the operational site over the next 6 hours. To adapt to the requirements of ballistic calculations, coordinate transformation is used to convert the wind field vector in the geographic coordinate system to the ballistic coordinate system, obtaining the longitudinal wind component (along the projectile's range direction, directly affecting the projectile's axial drag and range deviation) and the crosswind component (horizontal wind perpendicular to the projectile's range direction, causing lateral drift and attitude disturbances in the trajectory).

[0116] S33. In the longitudinal and transverse wind dimensions, the velocity vector of the projectile relative to the air is quantified based on the wind speed vector and the measured velocity of the projectile relative to the ground.

[0117] S34. In the first revised anti-aircraft gun trajectory equation, the projectile velocity is replaced by the velocity vector of the projectile relative to the air, and the velocity change rate in the range direction caused by the longitudinal wind and the velocity change rate perpendicular to the range direction caused by the crosswind are calculated respectively.

[0118] S35. Calculate the displacement correction in the range direction by integrating along the trajectory time series based on the velocity change rate in the range direction, and solve for the drift displacement perpendicular to the range direction by integrating along the trajectory time series based on the velocity change rate perpendicular to the range direction.

[0119] S36. The range displacement correction and the lateral drift displacement are spatially superimposed to generate corrected ballistic trajectory parameters that take into account the wind field coupling effect.

[0120] S37. Based on the modified ballistic trajectory parameters considering the wind field coupling effect, calculate the range and azimuth for different operational periods.

[0121] In one specific embodiment, the anti-aircraft gun trajectory wind field coupling correction is achieved through the following steps: First, based on historical meteorological data and reanalysis data of the operation point and preset area, a spatiotemporal evolution wind field dynamic model is constructed. This model accurately characterizes the spatiotemporal variation characteristics of the wind field in the horizontal and vertical directions during the operation period through multi-source data fusion and numerical simulation technology. Next, the wind field model is used to numerically simulate the temporal evolution of the wind field under typical meteorological patterns, and the simulated wind field vector is decomposed into the ballistic coordinate system to obtain multiple velocity vectors.

[0122] Next, the projectile velocity in the air resistance formula needs to be replaced with the projectile's velocity vector relative to the air, that is:

[0123] ;

[0124] In the formula, Let V be the velocity vector of the projectile relative to the ground. Wind speed vector (w) x w y w z In this embodiment, only longitudinal and transverse winds are typically considered, while vertical winds are not considered.

[0125] Its mold length is:

[0126] ;

[0127] Since the invention ultimately generates a two-dimensional ray map, vertical wind is not considered in practical applications. However, in order to ensure the accuracy of the overall calculation values, the effects of longitudinal wind, crosswind, and vertical wind must be handled together in the calculation, but their values ​​are ignored.

[0128] The specific corrections to the equations of motion are reflected in the following directional dynamic modeling:

[0129] (1) The wind causes the relative velocity between the projectile and the air to change The formula for the influence of longitudinal wind is:

[0130] ;

[0131] (2) In the lateral drift direction, lateral resistance is introduced, and the formula for the effect of crosswind is:

[0132] ;

[0133] (3) The formula for the influence of vertical wind is:

[0134] ;

[0135] In addition, the position update equation is synchronously corrected as follows:

[0136] ;

[0137] In summary, the second-corrected trajectory equations for anti-aircraft guns are obtained:

[0138] .

[0139] Finally, by iteratively solving the modified set of motion equations, the ballistic range and azimuth under different wind conditions during different operational periods can be obtained, and after ignoring vertical wind, it provides a basis for subsequent two-dimensional firing range maps.

[0140] S4. By integrating the range and altitude at various firing angles and the range and azimuth at various operational periods, an airspace firing boundary is generated on the data base map, and a safe firing boundary map for anti-aircraft artillery hail suppression operations is output.

[0141] Furthermore, such as Figure 5 As shown, step S4 includes:

[0142] S41. By fusing the range and altitude at various firing angles and the range and azimuth at various operational phases, a ballistic trajectory parameter set and a parameter space composed of altitude, range, and azimuth are obtained. Specifically, by fusing the range, altitude, and azimuth parameters at different firing angles and during operational phases, a parameter space is generated with altitude H, range χ, and azimuth θ as axes. Within this space, each discrete point corresponds to the limiting ballistic performance (such as maximum range and altitude) at a specific angle.

[0143] S42. Register the ballistic parameter set with the base map using coordinate systems, and dynamically eliminate airspace regions conflicting with terrain obstacles based on the maximum shell elevation threshold constrained by ballistic performance. In this step, based on the base map, affine transformations are used to align the ballistic parameter space with the geographic coordinate system. Spatial algebraic operations are performed to calculate the geometric intersection of the ballistic trajectory and terrain obstacles in real time: if the maximum shell elevation is lower than the top elevation of the obstacle (i.e., H... 弹道 <H 障碍 If a region is identified as a spatial conflict region, it will be dynamically removed from the parameter space.

[0144] S43. In the parameter space composed of altitude, range, and azimuth, perform reachability analysis on the airspace region after culling operations on a per-directional angle basis, generating discrete radii envelopes. Then, perform multi-angle radii envelope interpolation on these discrete radii envelopes to generate a continuous and smooth airspace radii boundary surface. In the reachability analysis for each azimuth angle θ, first extract the set of out-of-field boundary points {χ(θ), H(θ)} of the non-conflict region. Subsequently, use a non-uniform rational B-spline interpolation algorithm to fit the discrete boundary points along the azimuth dimension, generating a smooth and closed airspace radii envelope surface.

[0145] S44. Overlay the airspace firing boundary surface onto the data base map, and combine it with elevation data for layered rendering to generate a two-dimensional safe firing boundary map for anti-aircraft artillery hail suppression operations. Through the above process, the airspace safety boundary for anti-aircraft artillery hail suppression operations is transformed from abstract parameters into an intuitive and visual firing boundary map. This firing boundary map is used for paper-based deployment at work points through planar projection conversion, and is also embedded into the command system as a digital layer, providing a spatial carrier for terrain registration of airspace safety parameters and providing visual decision support for adjusting firing elevation angles, predicting ammunition trajectories, and providing real-time risk warnings under complex terrain conditions.

[0146] Additionally, this invention, based on computer and GIS technology, and in accordance with GB / T 39782-2021 "Specifications for Drawing Safety Fire Boundary Maps of Artificial Weather Modification Rocket Operation Sites," develops a system for generating fire boundary maps of anti-aircraft gun hail suppression operation sites, referencing... Figure 6 ,include:

[0147] The basic database subsystem stores operational site information, terrain data, satellite remote sensing data, UAV aerial photography data, meteorological observation data, and meteorological reanalysis data. It also strictly adheres to the requirements of GB / T 39782-2021 to build a standardized data warehouse, supporting spatiotemporal indexing and multi-resolution data fusion.

[0148] The data analysis subsystem executes the methods described above; the firing range map production subsystem, according to GB / T39782-2021 "Specifications for Drawing Safety Firing Range Maps of Artificial Weather Modification Rocket Operation Sites," draws a safety firing range map of preset accuracy based on the data output by the basic database subsystem and the data analysis subsystem. The data analysis subsystem includes: a data acquisition module, used to acquire multi-element geographic data of the anti-aircraft gun hail suppression operation site and its preset range via UAV aerial surveying, and generate a data base map; and a first analysis module, used to determine air density parameters based on the altitude of the operation site extracted from the multi-element geographic data, and utilize air... The density parameter is used to perform altitude-adaptive correction on the pre-constructed anti-aircraft gun trajectory equations to calculate the range and altitude at various firing angles. The second analysis module is used to construct a dynamic wind field model for the operational period based on the acquired historical meteorological observation data and meteorological reanalysis data. The wind field vector simulated by the dynamic wind field model is used to correct the trajectory deviation of the anti-aircraft gun trajectory equations to calculate the range and azimuth for different operational periods. The firing range map output module is used to integrate the range and altitude at various firing angles and the range and azimuth for different operational periods, generate the airspace firing range boundary on the data base map, and output a safe firing range map for anti-aircraft gun hail suppression operations.

[0149] The boundary map management subsystem is used to respond to instructions or detect status changes, and to provide boundary map review reminders and update alarms. Specifically, it adds functions for regular boundary map review, regular update reminders and early warnings, and realizes reminder and alarm functions based on relevant requirements and time limits, reminding work units to review and update the boundary map on time, realizing dynamic updates of the boundary map, and improving safety management capabilities.

[0150] The aforementioned system adopts a modular architecture, connecting various subsystems through microservice interfaces. It supports cloud deployment and multi-terminal collaborative operation, meeting the intelligent firing range management needs of anti-aircraft artillery hail suppression operations in complex terrain conditions. This invention integrates electronic firing range maps into the anti-aircraft artillery operation system, using technical means to forcibly restrict operational parameters, forming an "electronic fence" to effectively eliminate the risks of human error. Simultaneously, the system incorporates a built-in function for periodic review and update reminders, establishing a dynamic management mechanism for firing range maps to ensure continuous compliance with safety requirements. Through a technical approach of "data precision, automated production, electronic supervision, and dynamic management," this invention constructs a safer and more reliable firing range map management system, providing a solid technical guarantee for artificial hail suppression operations.

[0151] Furthermore, embodiments of the present invention provide a computer-readable storage medium storing computer-executable instructions, which, when executed by a processor, implement the method for generating the anti-aircraft gun hail suppression point firing boundary map as described above.

[0152] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0153] This invention is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each block of the flowchart illustrations and / or block diagrams, as well as combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions.

[0154] It should be noted that any reference numerals placed between parentheses in the claims should not be construed as limiting the claims. The word "comprising" does not exclude the presence of components or steps not listed in the claims. The word "a" or "an" preceding a component does not exclude the presence of a plurality of such components. The invention can be implemented by means of hardware comprising several different components and by means of a suitably programmed computer. In claims that enumerate several means, several of these means may be embodied by the same hardware. The use of the terms first, second, third, etc., is merely for convenience of expression and does not indicate any order. These terms can be understood as part of the component names.

[0155] Furthermore, it should be noted that in the description of this specification, the terms "one embodiment," "some embodiments," "embodiment," "example," "specific example," or "some examples," etc., refer to specific features, structures, materials, or characteristics described in connection with that embodiment or example, which are included in at least one embodiment or example of the present invention. In this specification, the illustrative expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in one or more embodiments or examples. Furthermore, without contradiction, those skilled in the art can combine and integrate the different embodiments or examples described in this specification, as well as the features of different embodiments or examples.

[0156] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the claims should be interpreted to include both the preferred embodiments and all changes and modifications falling within the scope of the invention.

[0157] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, then this invention should also include these modifications and variations.

Claims

1. A method for generating a point firing boundary map for anti-aircraft gun hail suppression operations, characterized in that, include: The aerial survey was conducted using drones to obtain multi-element geographic data of anti-hail gun operation points and preset areas, and a data base map was generated. The air density parameter is determined based on the altitude of the work site extracted from multi-element geographic data, and the air density parameter is used to perform altitude-adaptive correction on the pre-constructed anti-aircraft gun trajectory equation in order to calculate the range and altitude at each firing angle. Based on the acquired historical meteorological observation data and meteorological reanalysis data, a wind field dynamic model for the operation period is constructed. The wind field vector simulated by the wind field dynamic model is used to correct the trajectory deviation of the anti-aircraft gun ballistic motion equation, so as to calculate the range and azimuth of different operation periods. By integrating the range and altitude at various firing angles and the range and azimuth at various operational periods, an airspace firing boundary is generated on the data base map, and a safe firing boundary map for anti-aircraft artillery hail suppression operations is output.

2. The method for generating a firing boundary map for anti-aircraft gun hail suppression operations as described in claim 1, characterized in that, Aerial surveys using drones were conducted to acquire multi-element geographic data of anti-aircraft artillery hail suppression operation points and preset areas, and a data base map was generated, including: Based on the coordinates of the currently pre-deployed and historically deployed hail suppression operation points, the drone equipped with lidar, photography module and multispectral sensor is controlled to first perform lidar scanning to generate an elevation reference layer according to the preset radial route, and then switch to the preset spiral progressive route to perform aerial surveying tasks. During aerial surveying, the elevation data of the lidar is compared with the elevation reference layer in real time. When a terrain change zone with an elevation difference exceeding a preset threshold is detected, the photography module is triggered to perform multi-angle image verification of the terrain change zone. When the verification is confirmed, the scanning frequency of the lidar is adjusted in coordination to generate a geo-registered enhanced elevation layer. Based on the spectral data of ground features collected by multispectral sensors and the distribution characteristics of the enhanced elevation layer, vegetation mask and artificial building outlines are extracted through a pre-built intelligent segmentation model. Combined with historical topographic benchmark data, dynamic ground feature interference is eliminated through comparison and iteration, and a clear topographic map containing static ground feature vectors and elevation information is output. Spatially fuse the topographic map with the acquired imagery, combine it with artificial buildings including houses, villages and factories based on the outline annotation of artificial buildings, and inject plane coordinate system parameters and time labels to generate a two-dimensional data base map including topographic elevation, static feature distribution and image base.

3. The method for generating a firing boundary map for anti-aircraft gun hail suppression operations as described in claim 1, characterized in that, The air density parameter is determined based on the altitude of the operational site extracted from multi-element geographic data. This air density parameter is then used to perform altitude-adaptive corrections on the pre-constructed anti-aircraft gun trajectory equations to calculate the range and altitude at various firing angles, including: The altitude of the work site is extracted from multi-element geographic data. Combined with the obtained vertical lapse rate of atmospheric temperature and predefined sea level reference parameters, and based on the correlation between gas molecule properties and thermodynamic constants, the air density corresponding to the altitude is calculated. The nonlinear relationship of the drag coefficient is determined by the ratio of the instantaneous velocity of the projectile to the speed of sound in the air medium. The reciprocal of the ballistic coefficient is calculated by combining the inherent properties of the projectile and the air density. Based on the horizontal influence term that couples the projectile velocity and horizontal velocity components using the reciprocal of the ballistic coefficient as a scaling factor, the rate of change of the horizontal velocity in the trajectory equation of the anti-aircraft gun is dynamically adjusted. Based on the vertical velocity attenuation benchmark caused by gravitational acceleration, a product term of projectile velocity and vertical velocity components is introduced, with the reciprocal of the ballistic coefficient as a scaling factor, to construct a vertical influence term of drag and gravity. Based on the vertical influence term, the rate of change of vertical velocity in the trajectory equation of the anti-aircraft gun is dynamically adjusted. The time-step numerical iterative calculation is applied to the first revised anti-aircraft gun ballistic equation to update the instantaneous motion state of the projectile in the horizontal and vertical directions one by one until the conditions for the projectile to hit the ground are met, and the calculation process is terminated. Finally, the range and altitude corresponding to each firing angle are output. Among them, the ballistic motion equation is a dynamic model composed of a set of velocity differential equations in the horizontal and vertical directions and a set of position update equations. It constructs a complete description of the projectile's trajectory by using the inverse of the ballistic coefficient to quantify the nonlinear attenuation effect of air resistance on the velocity component and superimposing the effect of gravity.

4. The method for generating a point firing boundary map for anti-aircraft gun hail suppression operations as described in claim 1, characterized in that, The first revised expression for the trajectory equation of the anti-aircraft gun is as follows: ； In the formula, The reciprocal of the ballistic coefficient, ρ is the air density, and c is the air density. d Where A is the drag coefficient, v is the projectile cross-sectional area, m is the projectile velocity, and v is the projectile mass. x For the horizontal velocity component, v y denoted as the vertical velocity component, and g as the acceleration due to gravity.

5. The method for generating a point firing boundary map for anti-aircraft gun hail suppression operations as described in claim 4, characterized in that, Based on historical meteorological observation data and meteorological reanalysis data, a dynamic wind field model for the operational period was constructed. The wind field vector simulated by the dynamic wind field model was used to correct the trajectory deviation of the anti-aircraft gun's ballistic motion equation, so as to calculate the range and azimuth angle for different operational periods, including: By integrating historical meteorological observation data and meteorological reanalysis data of the acquired work site and surrounding preset range, a dynamic wind field model for the work period is constructed. Based on the wind field dynamic model, the temporal changes of the wind field at the work site are numerically simulated under various meteorological modes. The wind field vector is obtained and decomposed to the ballistic coordinate system to obtain the longitudinal wind component along the projectile range and the crosswind component perpendicular to the firing direction. In the longitudinal and transverse wind dimensions, the velocity vector of the projectile relative to the air is quantified based on the wind speed vector and the measured velocity of the projectile relative to the ground. In the first revised anti-aircraft gun trajectory equation, the projectile velocity is replaced by the velocity vector of the projectile relative to the air, and the rate of velocity change in the range direction caused by the longitudinal wind and the rate of velocity change in the perpendicular to the range direction caused by the crosswind are calculated respectively. The displacement correction in the range direction is calculated by time-series integration along the trajectory based on the rate of change of velocity in the range direction, and the drift displacement perpendicular to the range direction is solved by time-series integration along the trajectory based on the rate of change of velocity perpendicular to the range direction. The range displacement correction and the lateral drift displacement are spatially vector-superimposed to generate corrected ballistic trajectory parameters that take into account the wind field coupling effect. Based on the modified ballistic trajectory parameters that take into account the wind field coupling effect, the range and azimuth angle corresponding to different operational periods are calculated.

6. The method for generating a point firing boundary map for anti-aircraft gun hail suppression operations as described in any one of claims 1-5, characterized in that, By integrating the range and elevation at various firing angles and the range and azimuth at various operational phases, an airspace firing boundary is generated on the data base map, outputting a safe firing boundary map for anti-aircraft artillery hail suppression operations, including: By integrating the range and altitude at various firing angles and the range and azimuth at various operational phases, a ballistic trajectory parameter set and a parameter space composed of altitude, range, and azimuth are obtained. The ballistic parameter set and the data base map are registered in coordinate system, and airspace areas that conflict with terrain obstacles are dynamically eliminated based on the maximum firing height threshold of the projectile constrained by ballistic performance. The accessibility analysis of the airspace region after culling operation is performed on an azimuth-by-azimuth basis within the parameter space composed of firing height, range and azimuth angle. Discrete firing envelopes are generated. Multi-angle firing envelope interpolation is performed on the discrete firing envelopes to generate a continuous and smooth airspace firing boundary surface. By overlaying the airspace firing boundary surface onto the data base map and combining it with the elevation data for layered rendering, a two-dimensional safe firing boundary map for anti-aircraft artillery hail suppression operations is generated.

7. A system for generating a point firing boundary map for anti-aircraft gun hail suppression operations, characterized in that, include: The basic database subsystem is used to store operational site information data, terrain data, satellite remote sensing data, UAV aerial photography data, meteorological observation data, and meteorological reanalysis data. A data analysis subsystem is configured to perform the method as described in any one of claims 1-6; The radiation boundary map creation subsystem is used to draw a safety radiation boundary map with a preset accuracy based on the data output by the basic database subsystem and the data analysis subsystem. The boundary map management subsystem is used to respond to commands or detect status changes, and to provide boundary map review reminders and update alarms.

8. The anti-aircraft gun hail suppression operation point firing boundary map generation system as described in claim 7, characterized in that, The data analysis subsystem includes: The data acquisition module is used to acquire multi-element geographic data of anti-aircraft gun hail suppression operation points and preset ranges through UAV aerial surveying, and generate a data base map; The first analysis module is used to determine the air density parameter based on the altitude of the work site extracted from multi-element geographic data, and to use the air density parameter to perform altitude-adaptive correction on the pre-constructed anti-aircraft gun trajectory equation in order to calculate the range and altitude at each firing angle. The second analysis module is used to construct a wind field dynamic model for the operation period based on the acquired historical meteorological observation data and meteorological reanalysis data, and to use the wind field vector simulated by the wind field dynamic model to correct the trajectory deviation of the anti-aircraft gun ballistic motion equation, so as to calculate the range and azimuth of different operation periods. The firing range map output module is used to integrate the range and elevation at various firing angles and the range and azimuth at various operational periods to generate airspace firing range boundaries on the data base map and output a safe firing range map for anti-aircraft artillery hail suppression operations.

9. A computer-readable storage medium storing computer-executable instructions thereon, characterized in that, When the executable instructions are executed by the processor, they implement the method as described in any one of claims 1-6.

Images