Trajectory meteorological data acquisition method and system based on spatial interpolation of probe data

Abstract

The present application provides a method and system for obtaining ballistic meteorological data based on spatial interpolation of probe data, which comprises: calculating the trajectory coordinates of a projectile according to standard meteorological conditions, the firing angle of a cannon, the initial velocity of the projectile, and longitude and latitude coordinates; converting the trajectory coordinates of the projectile into geodetic coordinates; interpolating the data of multiple upper air meteorological probes to each height of the trajectory coordinates of the projectile by an interpolation method; horizontally interpolating the meteorological probe data of the same height to a grid; selecting the meteorological probe data of four adjacent grid points corresponding to the height of each geodetic coordinate of the trajectory of the projectile, and interpolating the meteorological probe data to the coordinate point of the projectile. The present application solves the technical problem of low accuracy of full trajectory meteorological data acquisition, which restricts the firing accuracy of the cannon.

Description

Technical Field

[0001] This invention relates to the field of meteorological observation, specifically to a method and system for acquiring ballistic meteorological data based on spatial interpolation of observation data. Background Technology

[0002] Projectiles fired from artillery and other weapon platforms are affected by meteorological factors such as temperature, air pressure, relative humidity, and wind during their flight through the air. If accurate meteorological data is not provided during trajectory correction, firing accuracy will inevitably be affected. For example, the existing invention patent application document "Auxiliary Firing Device and Method" with publication number CN114216363A includes: a handheld meteorological instrument for collecting meteorological data; a handheld laser rangefinder for obtaining the target's distance and elevation angle; and a mobile terminal for displaying a first interface containing weapon selection controls; in response to a first operation on the weapon selection controls, determining the weapon type, obtaining the corresponding projectile shape coefficient and drag coefficient, and displaying a second interface containing read and start controls; in response to a second operation on the read controls, reading the meteorological data collected by the handheld meteorological instrument; and in response to a third operation on the start controls, performing trajectory calculations based on the projectile shape coefficient and drag coefficient, meteorological data, and the target's distance and elevation angle, obtaining and outputting elevation and azimuth corrections. The existing technology has a simple method for acquiring meteorological data, but this existing solution only uses a handheld meteorological instrument to acquire meteorological data. In the application scenario of long-range shooting, it cannot accurately and comprehensively acquire ballistic meteorological data. Furthermore, the existing invention patent application document CN115544710A, entitled "A Method and Device for Firing Simulation Evaluation of a Direct-Fire Weapon Simulation Terminal," describes a method that includes: confirming the line-of-sight relationship between the direct-fire weapon simulation terminal and the target through laser emission; when the direct-fire weapon simulation terminal and the target are in line of sight, calculating the center position of the light spot illuminated by the direct-fire weapon simulation terminal on the target, and calculating the simulated bullet trajectory based on the center position of the light spot; performing collision deduction between the simulated bullet trajectory and the target's movement path to determine whether the target has been hit; when the target is hit, calculating the hit position, and matching the damage status of the target from a damage comparison table according to the bullet type and the hit position. This existing technology is applicable to direct-fire weapons. While it considers the effects of air resistance and weather conditions on bullet trajectory during ballistic simulation, it is only suitable for short-range direct-fire weapons and cannot solve the problem of acquiring meteorological data for long-range weapons. Furthermore, with the continuous development of artillery technology and the increasing range of artillery, the coverage radius of current single-station meteorological detection is insufficient to cover the entire projectile trajectory, making it impossible to accurately acquire full-trajectory meteorological data and thus affecting artillery firing accuracy.

[0003] In summary, existing technologies suffer from low accuracy in acquiring full-trajectory meteorological data, which limits the accuracy of artillery firing. Summary of the Invention

[0004] The technical problem to be solved by this invention is how to solve the problem that the low accuracy of full-trajectory meteorological data acquisition restricts the accuracy of artillery firing.

[0005] The present invention solves the above-mentioned technical problems by adopting the following technical solution: a ballistic meteorological data acquisition method based on spatial interpolation of probe data includes:

[0006] S1. Acquire artillery firing and meteorological data, and calculate the trajectory coordinates of the projectile based on the artillery firing and meteorological data. The artillery firing and meteorological data include: standard meteorological conditions, artillery firing angle, projectile initial velocity, and latitude and longitude coordinates.

[0007] S2. Convert the projectile's trajectory coordinates to ground coordinates.

[0008] S3. Using pre-set interpolation logic, interpolate the data from at least two high-altitude meteorological stations into the trajectory coordinates of the projectile at at least two pre-set altitudes.

[0009] S4. Horizontally interpolate the meteorological station data at the preset altitude to the preset grid.

[0010] S5. Based on the preset altitude corresponding to the geodetic coordinates of each projectile's flight trajectory, select the meteorological observation data of the grid points adjacent to the preset altitude in the preset grid, and use the bilinear interpolation method to interpolate the meteorological observation data of the grid points to the projectile coordinate points to obtain the ballistic meteorological data.

[0011] This invention utilizes high-altitude meteorological data from multiple stations to ensure that the coverage radius extends across the entire projectile flight trajectory. It accurately acquires meteorological data for the entire trajectory, rapidly and accurately obtaining meteorological elements at different locations along the projectile's trajectory. This meets the needs of artillery, air defense, and other fields for meteorological support in long-range strike applications for acquiring ballistic meteorological data and correcting the trajectory.

[0012] In a more specific technical solution, in step S1, the coordinates of the projectile's flight trajectory are calculated using the following logic:

[0013]

[0014] In the formula, c is the ballistic coefficient, and H τ (y) is the air density function, G(v) rτ ) is the air resistance function, v x v y v zThese are the components of velocity along the x, y, and z axes, respectively. x w z These are the ballistic longitudinal wind and the ballistic crosswind, respectively, where g is the acceleration due to gravity.

[0015] In a more specific technical solution, the gravitational acceleration g is obtained using the following logical processing:

[0016] g=9.80616(1-0.0026373cos2B+0.0000059cos 2 2B)

[0017] In the formula, B is the latitude of the upper-air meteorological sounding station.

[0018] In a more specific technical solution, step S2 includes:

[0019] S21. Based on the projectile's flight trajectory coordinates, process the projectile's ground coordinates and convert the projectile's ground coordinates into geocentric rectangular coordinates.

[0020] S22. Convert the geocentric rectangular coordinates to the geodetic coordinates of the projectile's flight trajectory.

[0021] In a more specific technical solution, in step S21, the following logic is used to convert the projectile's geodetic coordinates into the projectile's flight trajectory geodetic coordinates:

[0022]

[0023] In the formula, B, L, and H are the latitude, longitude, and elevation coordinates of the upper-air meteorological sounding station, respectively; N is the radius of curvature of the vertex-groove passing through the station; and e 2 It is the square of the first eccentricity.

[0024] In a more specific technical solution, in step S22, the following logic is used to convert the geocentric rectangular coordinates into the geodetic coordinates of the projectile's flight trajectory:

[0025]

[0026] In the formula, x, y, and z are the station center coordinates of the upper-air meteorological station relative to the artillery position with the origin as the origin, X, Y, and Z are the geocentric rectangular coordinates of the upper-air meteorological station, and X0, Y0, and Z0 are the geocentric rectangular coordinates of the artillery position.

[0027] In a more specific technical solution, in step S3, the data from the upper-air meteorological sounding station is interpolated into the coordinates of the projectile's flight trajectory using the following logic:

[0028]

[0029] In the formula, y ih is the meteorological element value interpolated to the trajectory height of the projectile. i For the height of the projectile trajectory, y i+1 y i-1 These represent the meteorological element values ​​of the layers adjacent to the projectile's trajectory height, h. i-1 h i+1 These represent the heights of the adjacent layers above and below the projectile's trajectory height.

[0030] This invention utilizes simultaneous detection data from multiple upper-air meteorological stations for spatial interpolation to obtain meteorological data at non-detection coordinates. This not only ensures the accuracy of data acquisition but also solves the problem of not being able to obtain full-trajectory meteorological data, thereby improving the accuracy of correcting ballistic meteorological conditions.

[0031] In a more specific technical solution, in step S4, the meteorological station data is horizontally interpolated to a preset grid using the Kriging interpolation method, based on the following logic:

[0032]

[0033] In the formula, λ i Z(x) i The weights obtained by assigning weights, Z(x) i (x) is a known position. i The detected value at Z * (x0) is the estimated value at position x0, and n is the estimated value of Z. * The number of actual probe values ​​used for the (x0) value.

[0034] In a more specific technical solution, in step S5, the following logic is used to interpolate the meteorological observation data of the grid points to the coordinate points of the projectile:

[0035]

[0036]

[0037]

[0038] Among them: Q 11 Q 12 Q 21 Q 22 These are four adjacent grid points, and (x1-y1), (x1-y2), (x2-y1), and (x2-y2) represent Q. 11 Q 12 Q 21 Q 22 The coordinates of f(Q) 11 f(Q) 12 f(Q) 21 f(Q) 22f(R1) and f(R2) are the detection values ​​of four adjacent grid points, respectively. f(R1) and f(R2) are the linear interpolation values ​​of the four grid points on the x-axis, respectively. f(P) is the result of bilinear interpolation.

[0039] In more specific technical solutions, ballistic meteorological data acquisition systems based on spatial interpolation of probe data include:

[0040] The flight trajectory coordinate calculation module is used to acquire artillery firing and meteorological data, and calculate the projectile flight trajectory coordinates based on the artillery firing and meteorological data. The artillery firing and meteorological data include: standard meteorological conditions, artillery firing angle, projectile initial velocity, and latitude and longitude coordinates.

[0041] The flight trajectory coordinate conversion module is used to convert the projectile's flight trajectory coordinates into the projectile's flight trajectory geodetic coordinates. The flight trajectory coordinate conversion module is connected to the flight trajectory coordinate calculation module.

[0042] The upper-air meteorological station data interpolation module is used to interpolate at least two upper-air meteorological station data into the projectile flight trajectory coordinates according to at least two preset heights using preset interpolation logic. The upper-air meteorological station data interpolation module is connected to the flight trajectory coordinate conversion module.

[0043] The grid interpolation module is used to horizontally interpolate meteorological station data at a preset altitude to a preset grid. The grid interpolation module is connected to the upper-air meteorological station data interpolation module.

[0044] The projectile coordinate point difference module is used to select meteorological observation data of grid points adjacent to the preset altitude in the preset grid according to the preset altitude corresponding to the geodetic coordinates of each projectile flight trajectory. The meteorological observation data of the grid points is interpolated to the projectile coordinate points using the bilinear interpolation method to obtain the ballistic meteorological data. The projectile coordinate point difference module is connected to the grid interpolation module.

[0045] Compared with the prior art, the present invention has the following advantages: The present invention utilizes high-altitude meteorological detection data from multiple stations, and the meteorological detection from multiple stations enables the protection radius to cover the entire trajectory of the projectile, accurately acquire meteorological data of the entire trajectory, and quickly and accurately acquire meteorological elements at different positions on the projectile trajectory, meeting the needs of artillery, air defense meteorological support and other fields for ballistic meteorological data acquisition and trajectory correction in long-range strike application scenarios.

[0046] This invention utilizes simultaneous data from multiple upper-air meteorological stations for spatial interpolation to obtain meteorological data at non-detection coordinates. This ensures data acquisition accuracy while solving the problem of not being able to obtain full-trajectory meteorological data, thus improving the accuracy of ballistic meteorological condition corrections. This invention addresses the technical problem in existing technologies where low accuracy in acquiring full-trajectory meteorological data restricts artillery firing accuracy. Attached Figure Description

[0047] Figure 1 This is a schematic diagram of the basic steps of the ballistic meteorological data acquisition method based on spatial interpolation of detection data according to Embodiment 1 of the present invention;

[0048] Figure 2 This is a schematic diagram of the ballistic meteorological data acquisition method based on spatial interpolation of detection data according to Embodiment 1 of the present invention. Detailed Implementation

[0049] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below in conjunction with the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0050] Example 1

[0051] like Figure 1 As shown, the ballistic meteorological data acquisition method based on spatial interpolation of probe data provided by the present invention includes the following basic steps:

[0052] Step S1: Calculate the trajectory coordinates of the projectile;

[0053] In this embodiment, the coordinates of the projectile's flight trajectory center are calculated based on data such as standard meteorological conditions, artillery firing angle, projectile initial velocity, and latitude and longitude coordinates.

[0054] In this embodiment, the specific method for calculating the trajectory coordinates of the projectile is as follows:

[0055]

[0056] Where c is the ballistic coefficient, H τ (y) is the air density function, G(v) rτ ) is the air resistance function, v x v y v z These are the components of velocity along the x, y, and z axes, respectively. x w z These represent the ballistic longitudinal wind and the ballistic crosswind, respectively, where g is the acceleration due to gravity. The initial conditions are:

[0057]

[0058] Where θ0 is the launch angle and v0 is the initial velocity of the projectile. The acceleration due to gravity g is calculated as follows:

[0059] g=9.80616(1-0.0026373cos2B+0.0000059cos 2 2B)

[0060] Where B represents the latitude of the upper-air meteorological station.

[0061] Step S2: Convert the trajectory coordinates to geodetic coordinates;

[0062] In this embodiment, the coordinates of each upper-air meteorological detection station are converted into the coordinates of the artillery position center.

[0063] In this embodiment, the geodetic coordinates of each upper-air meteorological sounding station are first converted to geocentric rectangular coordinates, and then the geocentric rectangular coordinates are converted to the center coordinates of the artillery position. The calculation method for converting geodetic coordinates to geocentric rectangular coordinates is as follows:

[0064]

[0065] Where B, L, and H are the latitude, longitude, and elevation coordinates of the upper-air meteorological sounding station, respectively, N is the radius of curvature of the zonal circle passing through the station, and e 2 e is the square of the first eccentricity. 2 =0.006694380022900788.

[0066] The method for converting geocentric rectangular coordinates to artillery position center coordinates is as follows:

[0067]

[0068] Where x, y, and z are the station-centered coordinates of the upper-air meteorological station relative to the artillery position with the origin as the origin, X, Y, and Z are the geocentric rectangular coordinates of the upper-air meteorological station, and X0, Y0, and Z0 are the geocentric rectangular coordinates of the artillery position.

[0069] Step S3: Interpolation of meteorological data at the projectile height level;

[0070] In this embodiment, data from multiple upper-air meteorological stations are interpolated into the various altitudes of the projectile's flight trajectory coordinates using an interpolation method.

[0071] In this embodiment, data from multiple upper-air meteorological stations are interpolated to various altitudes of the projectile's flight trajectory coordinates using linear interpolation. The calculation method is as follows:

[0072]

[0073] Where: y i h is the meteorological element value interpolated to the trajectory height of the projectile. i For the height of the projectile trajectory, y i+1 y i-1These represent the meteorological element values ​​of the layers adjacent to the projectile's trajectory height, h. i-1 h i+1 These represent the heights of the adjacent layers above and below the projectile's trajectory height.

[0074] Step S4: Interpolate station data to the grid;

[0075] In this embodiment, meteorological observation data at the same altitude are horizontally interpolated onto the grid;

[0076] In this embodiment, meteorological observation data at the same altitude are horizontally interpolated to the grid using Kriging interpolation, employing a power model semivariogram. The calculation method is as follows:

[0077]

[0078] Where: λ i Z(x) i The weights obtained by assigning weights, Z(x) i (x) is a known position. i The detected value at Z * (x0) is the estimated value at position x0, and n is the estimated value of Z. * The number of actual probe values ​​used for the (x0) value.

[0079] Step S5: Interpolate grid data to projectile coordinates;

[0080] In this embodiment, meteorological observation data from four grid points with adjacent altitudes corresponding to the geodetic coordinates of each projectile's flight trajectory are selected and interpolated to the projectile's coordinate points.

[0081] In this embodiment, bilinear interpolation is used to interpolate meteorological observation data from four grid points adjacent to each projectile's flight trajectory geodetic coordinate altitude. The calculation method is as follows:

[0082]

[0083]

[0084]

[0085] Among them: Q 11 Q 12 Q 21 Q 22 These are four adjacent grid points, and (x1-y1), (x1-y2), (x2-y1), and (x2-y2) represent Q. 11 Q 12 Q 21 Q 22 The coordinates of f(Q) 11 f(Q)12 f(Q) 21 f(Q) 22 f(R1) and f(R2) are the detection values ​​of four adjacent grid points, respectively. f(R1) and f(R2) are the linear interpolation values ​​of the four grid points on the x-axis, respectively. f(P) is the result of bilinear interpolation.

[0086] like Figure 2 As shown, in this embodiment, firstly, the trajectory is calculated based on standard meteorological conditions, artillery firing angle, projectile initial velocity, latitude and longitude coordinates, etc., to obtain a series of coordinates of the projectile's flight trajectory; then, the projectile's flight trajectory coordinates are converted into a series of geodetic coordinates; meteorological data obtained from multiple upper-air meteorological stations are interpolated to the height h of the projectile at a certain moment of flight through meteorological data heights h1 and h2, respectively. n Up; Next, through h n Meteorological data from four upper-air meteorological stations, p1 to p4, were horizontally interpolated to the grid using Kriging interpolation; finally, h... n Meteorological data at coordinate p of the projectile were obtained by bilinear interpolation of meteorological observation data from four adjacent grid points on the horizontal coordinate of the projectile.

[0087] In summary, this invention utilizes high-altitude meteorological data from multiple stations to ensure that the coverage radius covers the entire trajectory of the projectile, accurately acquires meteorological data for the entire trajectory, and quickly and accurately obtains meteorological elements at different locations on the projectile's trajectory. This meets the needs of artillery, air defense, and other fields for meteorological support in long-range strike applications for acquiring ballistic meteorological data and correcting the trajectory.

[0088] This invention utilizes simultaneous data from multiple upper-air meteorological stations for spatial interpolation to obtain meteorological data at non-detection coordinates. This ensures data acquisition accuracy while solving the problem of not being able to obtain full-trajectory meteorological data, thus improving the accuracy of ballistic meteorological condition corrections. This invention addresses the technical problem in existing technologies where low accuracy in acquiring full-trajectory meteorological data restricts artillery firing accuracy.

[0089] The above embodiments are only used to illustrate the technical solutions of the present invention, and are not intended to limit it. Although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for acquiring ballistic meteorological data based on spatial interpolation of probe data, characterized in that, The method includes: S1. Acquire artillery firing and meteorological data, and calculate the trajectory coordinates of the projectile based on the artillery firing and meteorological data. The artillery firing and meteorological data include: standard meteorological conditions, artillery firing angle, projectile initial velocity, and latitude and longitude coordinates. S2. Convert the projectile's trajectory coordinates to ground coordinates. S3. According to no less than two preset altitudes, interpolate the data from no less than two high-altitude meteorological stations into the coordinates of the projectile's flight trajectory; In S3, the following logic is used to interpolate the upper-air meteorological station data into the projectile's flight trajectory coordinates: In the formula, These are meteorological element values ​​interpolated to the height of the projectile trajectory. The height of the projectile trajectory. , These represent meteorological element values ​​for the layers adjacent to the projectile's trajectory height. , These represent the heights of the adjacent layers above and below the projectile's trajectory height; S4. Horizontally interpolate the meteorological station data at the preset altitude to the preset grid. In S4, the meteorological station data is horizontally interpolated to a preset grid using the Kriging interpolation method, based on the following logic: In the formula, for The assigned weights, For known location The detected value at that location, For position The estimated value at that location, n For estimation The number of actual detection values ​​used; S5. Based on the preset altitude corresponding to the geodetic coordinates of each projectile's flight trajectory, select the meteorological observation data of the grid points adjacent to the preset altitude in the preset grid, and use the bilinear interpolation method to interpolate the meteorological observation data of the grid points to the projectile coordinate points to obtain the ballistic meteorological data.

2. The method for acquiring ballistic meteorological data based on spatial interpolation of probe data according to claim 1, characterized in that, In step S1, the coordinates of the projectile's flight trajectory are calculated using the following logic: In the formula, c is the ballistic coefficient. It is a function of air density. It is a function of air resistance. , , The speeds are respectively at x , y , z Components on the axis, , These are the ballistic longitudinal winds and the ballistic crosswinds, respectively. g This is the acceleration due to gravity.

3. The method for acquiring ballistic meteorological data based on spatial interpolation of probe data according to claim 2, characterized in that, The gravitational acceleration is obtained using the following logic processing. g : In the formula, The latitude of the upper-air meteorological sounding station.

4. The method for acquiring ballistic meteorological data based on spatial interpolation of probe data according to claim 1, characterized in that, Step S2 includes: S21. Based on the projectile's flight trajectory coordinates, process the projectile's ground coordinates and convert the projectile's ground coordinates into geocentric rectangular coordinates. S22. Convert the geocentric rectangular coordinates into the geodetic coordinates of the projectile's flight trajectory.

5. The method for acquiring ballistic meteorological data based on spatial interpolation of probe data according to claim 4, characterized in that, In step S21, the projectile's geodetic coordinates are converted into the projectile's flight trajectory geodetic coordinates using the following logic: In the formula, , and These are the latitude, longitude, and elevation coordinates of the upper-air meteorological sounding station. Let the radius of curvature of the zonal loop passing through the station be denoted as . It is the square of the first eccentricity.

6. The method for acquiring ballistic meteorological data based on spatial interpolation of probe data according to claim 4, characterized in that, In step S22, the geocentric rectangular coordinates are converted into the geodetic coordinates of the projectile's flight trajectory using the following logic: In the formula, , and The coordinates of the upper-air meteorological station relative to the artillery position are the station center coordinates. , and The coordinates of the upper-air meteorological sounding station are in geocentric rectangular coordinates. , and Let be the geocentric rectangular coordinates of the artillery position.

7. The method for acquiring ballistic meteorological data based on spatial interpolation of probe data according to claim 1, characterized in that, In step S5, the meteorological data from the grid points is interpolated to the coordinates of the projectile using the following logic: in: , , , Each of the four adjacent grid points is a separate grid. , , and They are respectively , , , coordinates , , , These are the detection values ​​of four adjacent grid points. , These are the linear interpolations of the detection values ​​from the four grid points on the x-axis. This is the result of bilinear interpolation.

8. A ballistic meteorological data acquisition system based on spatial interpolation of probe data, used to execute the ballistic meteorological data acquisition method based on spatial interpolation of probe data as described in any one of claims 1 to 7, characterized in that, The system includes: The flight trajectory coordinate calculation module is used to acquire artillery firing and meteorological data, and calculate the projectile flight trajectory coordinates based on the artillery firing and meteorological data, wherein the artillery firing and meteorological data include: standard meteorological conditions, artillery firing angle, projectile initial velocity, and latitude and longitude coordinates; The flight trajectory coordinate conversion module is used to convert the projectile's flight trajectory coordinates into the projectile's flight trajectory geodetic coordinates. The flight trajectory coordinate conversion module is connected to the flight trajectory coordinate calculation module. The upper-air meteorological station data interpolation module is used to interpolate at least two upper-air meteorological station data into the trajectory coordinates of the projectile at at least two preset altitudes using preset interpolation logic. The upper-air meteorological station data interpolation module is connected to the trajectory coordinate conversion module. A grid interpolation module is used to horizontally interpolate the meteorological station data at the preset altitude to a preset grid. The grid interpolation module is connected to the upper-air meteorological station data interpolation module. The projectile coordinate point difference module is used to select meteorological observation data of grid points adjacent to the preset altitude in the preset grid according to the preset altitude corresponding to the geodetic coordinates of each projectile flight trajectory, and interpolate the meteorological observation data of the grid points to the projectile coordinate points using a bilinear interpolation method to obtain ballistic meteorological data. The projectile coordinate point difference module is connected to the grid interpolation module.

Images