CN106291756B - The construction method of near space air virtual environment resource - Google Patents

Abstract

A method for constructing atmospheric virtual environment resources in adjacent space relates to a method for constructing atmospheric virtual environment resources. The invention aims to fill the problem of lack of atmospheric virtual environment resources in the adjacent space in the domestic virtual test system. The present invention acquires the historical data of the atmospheric environment resources in the adjacent space, and processes the abnormal data; then normalizes the height data obtained by each event; then converts the coordinates and obtains the spatial range of each event; The two-dimensional array of each height layer is processed by two-dimensional smoothing; finally, the three-dimensional array of each resource data is constructed to obtain the cubic grid of the three-dimensional array of temperature, density and pressure, and finally the data format is converted to complete the approach conforming to the SEDRIS specification Construction of space atmosphere virtual environment resources. The invention is applicable to the construction field of atmospheric environment resources near space in virtual experiment.

Description

Method for constructing atmospheric virtual environment resources in near space

technical field

The invention relates to a method for constructing atmospheric virtual environment resources.

Background technique

Virtual natural environment resources are an important part of supporting virtual experiments, and their perfection and credibility directly affect the fidelity of virtual experiments. Atmospheric environment is an important part of the comprehensive natural environment, and the construction of atmospheric virtual environment resources in adjacent space is of great significance for improving the virtual test system.

Adjacent space refers to the airspace 20 to 100 kilometers above the ground, including the stratosphere, mesosphere and low thermosphere. As an important part of comprehensive natural environment data and the basis of comprehensive natural environment database development, the generation and improvement of atmospheric environmental resources are very important. At present, the development and utilization of adjacent space is still in the initial stage, and the influencing factors and degree of influence of aircraft flying in this space are still unclear. If a physical aircraft is directly used to conduct experiments in atmospheric space, many unknown consequences will occur, and even bring immeasurable loss. Since the "Twelfth Five-Year Plan", domestic atmospheric virtual environment resources below 20km have been gradually improved, while the construction method of atmospheric virtual environment resources in adjacent spaces is only in its infancy.

Contents of the invention

The invention aims to fill the problem of lack of atmospheric virtual environment resources in the adjacent space in the domestic virtual test system.

A method for constructing a near-space atmospheric virtual environment resource includes the following steps:

Step 1. Obtain the historical data of atmospheric environment resources in the adjacent space, including longitude, latitude, height data and corresponding temperature data;

Step 2, processing the abnormal data in the longitude, latitude, height and temperature data;

Step 3. Normalize the height data acquired by each event: For the height corresponding to each event, first round the height data in each event; The same altitude value is regarded as the altitude layer of the same layer. If there are multiple temperature data in one altitude layer, the average value of multiple temperature data is taken as the temperature value of the altitude layer in the final event;

Step 4, coordinate conversion: convert the latitude and longitude coordinates into Cartesian coordinates;

Step 5. Obtain the spatial range of each event: take the minimum value of the starting point and the maximum value of the latitude and longitude of each event after the coordinate conversion, and draw the area range of a single event and the area range represented by all events at each level, which is convenient for later analysis and data processing;

Step 6. Construct a two-dimensional array of each altitude layer: take longitude and latitude as dimensions, and normalize it to κ′ kilometer interval, construct a two-dimensional array of temperature for each altitude layer; the area range represented by each event If there is an overlapping part, when constructing a two-dimensional array, the temperature value of the overlapping part is calculated as the temperature value of the area by calculating the average temperature of the two events, and if there is no overlap, it will be stored in the original temperature value;

Step 7, two-dimensional smoothing processing: performing 2D smoothing processing on each two-dimensional array;

Step 8. Construct the three-dimensional array of the temperature of this year and this month: the height starts from 20km and ends at 100km. For the altitude layer with an incremental interval of κ km, write the smoothed two-dimensional array of temperature layer by layer, and finally get the following 3D array with longitude, latitude and altitude as dimensions;

Step 9. According to steps 1 to 8, the temperature data of the same period in several years are constructed and averaged, and a reliable three-dimensional array of cubic grids is established as the adjacent space temperature data resource of the month in the virtual test;

Step 10, according to step 1 to step 9, carry out the same processing to density, pressure, obtain the cubic grid of the three-dimensional array of density and pressure;

Step 11, write the constructed data including temperature, density and pressure into text file T;

Step 12D, then read the atmospheric environment data from the text file T of writing data, according to the DRM (Data Representation Model, data representation model) of SEDRIS (Synthetic Environment Data Representation and Interchange Specification, comprehensive environmental data representation and exchange specification), The three specifications of SRM (Spatial Reference Model, spatial reference model) and EDCS (Environment Data Coding Specification, environmental data coding specification) represent the data, and save the atmospheric environment data in the SEDRIS standard format.

Preferably, the method for constructing the atmospheric virtual environment resource in the adjacent space also includes the construction process of the atmospheric horizontal wind field, and the specific process is as follows:

Step 12A, according to the range of latitude, determine the corresponding atmospheric horizontal wind field calculation formula:

(1) In the range of latitude from 15° to 80°, the geostrophic wind is calculated according to formula (1):

Among them, P is the air pressure, and ρ is the atmospheric density; is called the geostrophic parameter, Ω is the geostrophic angular velocity, is the geographic latitude; x is the distance to the east, and y is the distance to the north; u _g and v _g are the horizontal latitudinal and meridional geostrophic winds within the range of 15°-80°, respectively;

Then calculate the gradient wind according to formula (2):

in, a is the radius of the earth; u _gr and v _gr are the horizontal latitudinal and meridional gradient winds respectively;

(2) Within the latitude range of 15°S to 15°N, the gradient wind is calculated according to formula (3):

Among them, u _e and v _e are the horizontal latitudinal and meridional gradient winds within 15°S～15°N, respectively;

Step 12B. Construct a uniform cubic grid of the three-dimensional array of the atmospheric horizontal wind field: Substitute the previously constructed three-dimensional array of density and pressure in the area into the wind field calculation formula to obtain the uniform three-dimensional array of the geostrophic wind and gradient wind at the corresponding position Cube grid, which generates the atmospheric horizontal wind field resources in the adjacent space of the region;

Step 12C, write the constructed wind field data into the text file T.

Preferably, the specific process of processing the abnormal data in the longitude, latitude, height and temperature data described in step 2 is as follows:

For the longitude, latitude, altitude and temperature data in each event, the original average resolution of the data is obtained, and then the abnormal data is eliminated, and the previous data of the abnormal data is used as the basis, and the resolution is added on this basis to obtain Recursive data equal to the number of abnormal data removed.

Preferably, the specific process of coordinate transformation described in step 4 is as follows:

The grid of the finally established coordinate system is based on distance. Therefore, according to the coordinate projection system in Matlab2014a, the latitude and longitude data obtained in each event are converted into distance data in Cartesian coordinates and rounded to the nearest integer.

Preferably, said κ is equal to said κ'.

Preferably, κ=κ′=1.

Preferably, the specific process of obtaining the historical data of the atmospheric environment resources in the adjacent space in step 1 is:

Select the NC file of near-space atmospheric environment data corresponding to the specific space range and time range of SABER; read the NC file to obtain the historical data of near-space atmospheric environment resources.

Beneficial effect:

(1) The empirical atmospheric model adopted by the present invention can make the construction results of the atmospheric environment resources more accurate with the improvement of the accuracy of the detection data, without considering the complexity of the adjacent space atmosphere and the selection of parameterization schemes for many processes.

(2) The construction of near-space atmospheric environment resources is based on satellite detection data. The data source is authentic and reliable, and the construction results have high credibility.

(3) This method of constructing near-space atmospheric environment resources has universal applicability. For any region of any size in the world, the resource construction process is basically unchanged, and it can meet the requirements of near-space atmospheric environment in most cases in virtual experiments.

Description of drawings

Figure 1 is a schematic diagram of the regional scope of an event;

Figure 2 is a schematic diagram of the area range of 10 events at a certain altitude;

Figure 3 is a temperature distribution diagram after an event has been processed;

Figure 4 is a schematic diagram of the temperature distribution after smoothing the events shown in Figure 3;

FIG. 5 is a schematic flow chart of Embodiment 1;

FIG. 6 is a schematic flow chart of Embodiment 2.

Detailed ways

Specific implementation mode 1: This implementation mode is described in conjunction with FIG. 5 ;

A method for constructing a near-space atmospheric virtual environment resource includes the following steps:

Step 1. Obtain the historical data of atmospheric environment resources in the adjacent space, including longitude, latitude, height data and corresponding temperature data;

Step 2, processing the abnormal data in the longitude, latitude, height and temperature data; thereby obtaining a two-dimensional matrix without abnormal data, which is convenient for subsequent data processing operations such as normalization and averaging;

Step 3. Normalize the height data acquired by each event: For the height corresponding to each event, first round the height data in each event; The same altitude value is regarded as the altitude layer of the same layer. If there are multiple temperature data in one altitude layer, the average value of multiple temperature data is taken as the temperature value of the altitude layer in the final event;

例：针对某个事件中的高度数据(km)99.98444366、99.61489105、99.24510956、98.87528229、98.50534821、98.13528442、97.76511383、97.39484406、97.0244751、96.65398407，该事件每个高度对应的温度值为225.5196838、225.0498047、224.7708435、224.5986786 . This level corresponds to two temperature values of 225.5196838 and 225.0498047, and the average value of the two needs to be 225.28474425 as the temperature at the 100km level of the event. The processing methods for other levels are similar.

Step 4, coordinate conversion: convert the latitude and longitude coordinates into Cartesian coordinates;

Step 5. Obtain the spatial range of each event: take the minimum value of the starting point and the maximum value of the latitude and longitude of each event after the coordinate conversion, and draw the area range of a single event and the area range represented by all events at each level, which is convenient for later analysis and data processing;

Taking a height layer of 1 km as an example, the spatial range of temperature at the height layer is obtained, as shown in Figures 1 and 2. It can be seen from the figure that there is overlap between the ranges represented by each event.

Step 6. Construct a two-dimensional array of each altitude layer: take longitude and latitude as dimensions, and normalize it to κ′ kilometer interval, construct a two-dimensional array of temperature for each altitude layer; the area range represented by each event There are overlapping parts, so when constructing a two-dimensional array, the temperature value of the overlapping part is calculated as the temperature value of the area by calculating the average temperature of the two events, and the original temperature value is stored if there is no overlap;

As shown in Figure 3, it is a temperature distribution diagram after an event is processed;

Step 7, two-dimensional smoothing processing: perform 2D smoothing processing on each two-dimensional array, so that the transition between the overlapping part and other parts is more uniform;

As shown in Figure 4, Figure 4 is the temperature distribution after smoothing the event shown in Figure 3;

Step 8. Construct a three-dimensional array of the temperature of this year and this month: Since the downloaded NC file takes months as the time range, the time resolution of the constructed three-dimensional array is the monthly average. According to specific requirements, a certain spatial range can be intercepted; the height is from 20km Start and end at 100km, write the smoothed temperature two-dimensional array layer by layer for the height layer with increasing interval of κ km, and finally get a three-dimensional array with longitude, latitude and height as dimensions;

Step 9. According to steps 1 to 8, the temperature data of the same period in several years (for example, the first 4 years) are constructed and averaged to establish a more reliable three-dimensional array of cubic grids, which has universal applicability in virtual experiments. Nearby spatial temperature data resources for the month;

Step 10, according to step 1 to step 9, carry out the same processing to density, pressure, obtain the cubic grid of the three-dimensional array of density and pressure;

Step 11, write the constructed data including temperature, density and pressure into text file T;

Step 12D, then read the atmospheric environment data from the text file T of writing data, according to the DRM (Data Representation Model, data representation model) of SEDRIS (Synthetic Environment Data Representation and Interchange Specification, comprehensive environmental data representation and exchange specification), The three specifications of SRM (Spatial Reference Model, spatial reference model) and EDCS (Environment Data Coding Specification, environmental data coding specification) represent the data, and save the atmospheric environment data in the SEDRIS standard format.

Specific implementation mode two: this implementation mode is described in conjunction with FIG. 6 ;

A method for constructing a near-space atmospheric virtual environment resource includes the following steps:

Step 1. Obtain the historical data of atmospheric environment resources in the adjacent space, including longitude, latitude, height data and corresponding temperature data;

Step 2, processing the abnormal data in the longitude, latitude, height and temperature data; thereby obtaining a two-dimensional matrix without abnormal data, which is convenient for subsequent data processing operations such as normalization and averaging;

Step 3. Normalize the height data acquired by each event: For the height corresponding to each event, first round the height data in each event; The same altitude value is regarded as the altitude layer of the same layer. If there are multiple temperature data in one altitude layer, the average value of multiple temperature data is taken as the temperature value of the altitude layer in the final event;

例：针对某个事件中的高度数据(km)99.98444366、99.61489105、99.24510956、98.87528229、98.50534821、98.13528442、97.76511383、97.39484406、97.0244751、96.65398407，该事件每个高度对应的温度值为225.5196838、225.0498047、224.7708435、224.5986786 . This level corresponds to two temperature values of 225.5196838 and 225.0498047, and the average value of the two needs to be 225.28474425 as the temperature at the 100km level of the event. The processing methods for other levels are similar.

Step 4, coordinate conversion: convert the latitude and longitude coordinates into Cartesian coordinates;

Step 5. Obtain the spatial range of each event: take the minimum value of the starting point and the maximum value of the latitude and longitude of each event after the coordinate conversion, and draw the area range of a single event and the area range represented by all events at each level, which is convenient for later analysis and data processing;

Taking a height layer of 1 km as an example, the spatial range of temperature at the height layer is obtained, as shown in Figures 1 and 2. It can be seen from the figure that there is overlap between the ranges represented by each event.

Step 6. Construct a two-dimensional array of each altitude layer: take longitude and latitude as dimensions, and normalize it to κ′ kilometer interval, construct a two-dimensional array of temperature for each altitude layer; the area range represented by each event There are overlapping parts, so when constructing a two-dimensional array, the temperature value of the overlapping part is calculated as the temperature value of the area by calculating the average temperature of the two events, and the original temperature value is stored if there is no overlap;

As shown in Figure 3, it is a temperature distribution diagram after an event is processed;

Step 7, two-dimensional smoothing processing: perform 2D smoothing processing on each two-dimensional array, so that the transition between the overlapping part and other parts is more uniform;

As shown in Figure 4, Figure 4 is the temperature distribution after smoothing the event shown in Figure 3;

Step 8. Construct a three-dimensional array of the temperature of this year and this month: Since the downloaded NC file takes months as the time range, the time resolution of the constructed three-dimensional array is the monthly average. According to specific requirements, a certain spatial range can be intercepted; the height is from 20km Start and end at 100km, write the smoothed temperature two-dimensional array layer by layer for the height layer with increasing interval of κ km, and finally get a three-dimensional array with longitude, latitude and height as dimensions;

Step 9. According to steps 1 to 8, the temperature data of the same period in several years (for example, the first 4 years) are constructed and averaged to establish a more reliable three-dimensional array of cubic grids, which has universal applicability in virtual experiments. Nearby spatial temperature data resources for the month;

Step 10, according to step 1 to step 9, carry out the same processing to density, pressure, obtain the cubic grid of the three-dimensional array of density and pressure;

Step 11, write the constructed data including temperature, density and pressure into text file T;

Step 12A, determine the calculation formula of the atmospheric horizontal wind field: the atmospheric horizontal wind field refers to the horizontal gradient wind field, including latitudinal wind and meridional wind; since the characteristics of the wind field at different latitudes are also different, its calculation method is different in latitude In the range of 15° to 80° and above the equator, the gradient wind needs to be calculated according to different formulas; while the wind field between 15°S and 15°N can be obtained by linear interpolation;

According to the latitude range, determine the corresponding wind field calculation formula:

(1) In the range of latitude from 15° to 80°, the geostrophic wind is calculated according to formula (1):

Among them, P is the air pressure, and ρ is the atmospheric density; is called the geostrophic parameter, Ω is the geostrophic angular velocity, is the geographic latitude; x is the distance to the east, and y is the distance to the north; u _g and v _g are the horizontal latitudinal and meridional geostrophic winds within the range of 15°-80°, respectively;

Then calculate the gradient wind according to formula (2):

in, a is the radius of the earth; u _gr and v _gr are the horizontal latitudinal and meridional gradient winds respectively;

(2) The sky over the equator requires a special solution. Within the latitude range of 15°S to 15°N, the gradient wind is calculated according to formula (3):

Among them, u _e and v _e are the horizontal latitudinal and meridional gradient winds within 15°S～15°N, respectively;

Step 12B, constructing a uniform cubic grid of three-dimensional arrays of the atmospheric horizontal wind field: the atmospheric horizontal wind field is mainly calculated using the atmospheric density and pressure, the radius of the earth's orbit, and the angular velocity, etc., so the previously constructed density and pressure of the region The three-dimensional array is substituted into the calculation formula of the wind field to obtain the uniform cubic grid of the three-dimensional array of the geostrophic wind and the gradient wind at the corresponding position, thereby generating the atmospheric horizontal wind field resources in the adjacent space of the region;

Step 12C, writing the constructed wind field data into the text file T;

Step 12D, then read the atmospheric environment data from the text file T where the data is written, represent the data according to the three specifications of SEDRIS DRM, SRM and EDCS, and save the atmospheric environment data in the SEDRIS standard format.

So far, the construction of near-space atmospheric environment data resources has been completed, and the database including temperature, density, pressure and wind field has been constructed. The data of each atmospheric element is saved in the form of a three-dimensional uniform grid. In order to facilitate subsequent data conversion, it is written into a txt file.

Specific implementation mode three:

The specific process of processing the abnormal data in the longitude, latitude, height and temperature data described in step 2 of this embodiment is as follows:

For the longitude, latitude, altitude and temperature data in each event, the original average resolution of the data is obtained, and then the abnormal data is eliminated, and the previous data of the abnormal data is used as the basis, and the resolution is added on this basis to obtain Recursive data equal to the number of abnormal data removed.

Take the abnormal value processing of temperature data as an example: the abnormal maximum value is generally distributed in the last two data of the temperature data in each event. Taking the last 10 temperature data of the temperature data as an example, the original average resolution is obtained according to the temperature data in advance.率，温度数据的后10个温度数据：225.5196838、225.0498047、224.7708435、224.5986786、224.3938904、224.054245、223.5292511、222.8206482、221.9805908、221.1124725、220.2646179、219.3830872、218.3756866、217.1251678、215.4690247、213.3180847、210.8113708、208.0403748、9.96921E+ 36. 9.96921E+36, it can be seen that the last two data are abnormal data. Firstly, the original average resolution is obtained according to the temperature data, and all the temperature data before the outlier are taken as the adjacent temperature difference, and then these differences are averaged, and the average value is used as the original temperature resolution of the event. The last two data in the data are abnormal data. Based on the penultimate data, the last two data are recursively obtained by adding the original resolution on the basis of the penultimate data.

Other steps and parameters are the same as those in Embodiment 1 or Embodiment 2.

Specific implementation mode four:

The specific process of the coordinate conversion described in step 4 of the present embodiment is as follows:

Because the grid of the final coordinate system is based on distance, according to the coordinate projection system in Matlab2014a, the latitude and longitude data obtained in each event are converted into distance data in Cartesian coordinates and rounded to the nearest integer.

Other steps and parameters are the same as those in Embodiments 1 to 3.

Specific implementation mode five:

κ described in this embodiment is equal to κ′ described above. When κ=κ', the grid formed by longitude, latitude and height is a cube grid.

Other steps and parameters are the same as in one of the specific embodiments 1 to 4.

Specific implementation method six:

In this embodiment, κ=κ′=1.

Other steps and parameters are the same as one of the specific embodiments 1 to 5.

Specific implementation mode seven:

The specific process of obtaining the historical data of the atmospheric environment resources in the adjacent space in step 1 described in this embodiment is as follows:

Select the NC file of near-space atmospheric environment data corresponding to the specific space range and time range of SABER; read the NC file to obtain the historical data of near-space atmospheric environment resources.

Take the four-year data of a certain region as an example, 48 NC files in 48 months from 2012 to 2015, each file contains NC file information, stored data such as longitude, latitude, altitude, temperature, etc., using The NetCDF4Excel plug-in and matlab read NC files to obtain historical data of atmospheric environment resources in adjacent spaces.

Other steps and parameters are the same as one of the specific embodiments 1 to 6.

Claims (8)

1. The method for constructing the near space atmosphere virtual environment resource is characterized by comprising the following steps of:

step 1, acquiring historical data of atmospheric environment resources in an adjacent space, wherein the historical data comprises longitude data, latitude data, altitude data and corresponding temperature data;

step 2, processing abnormal data in the longitude, latitude, altitude and temperature data;

step 3, normalizing the height data acquired by each event: for the height corresponding to each event, rounding the height data in each event; then, layering is carried out at intervals of k kilometers, the same height value is regarded as the height layer of the same layer, and if a plurality of temperature data exist in one height layer, the average value of the plurality of temperature data is taken as the temperature value of the height layer in the final event;

step 4, coordinate conversion: converting the longitude and latitude coordinates into rectangular coordinates;

step 5, obtaining the spatial range of each event: taking the minimum value of the starting point and the maximum value of the ending point of the latitude in each event after coordinate conversion, and drawing the area range of a single event and the area range represented by all events in each height layer;

step 6, constructing a two-dimensional array of each height layer: taking longitude and latitude as dimensions, normalizing the dimensions into k' kilometer intervals, and constructing a temperature two-dimensional array of each height layer; the area range represented by each event has an overlapping part, when a two-dimensional array is constructed, the temperature value of the overlapping part is used as the temperature value of the area by solving the mean value of the two event temperatures, and if the two event temperatures are not overlapped, the temperature value is stored into the original temperature value;

and 7, two-dimensional smoothing: 2D smoothing is carried out on each two-dimensional array;

step 8, constructing a three-dimensional array of the temperature in this month of the year: starting from 20km and ending at 100km, writing the smoothed two-dimensional temperature array layer by layer aiming at an altitude layer with kappa kilometers as incremental intervals, and finally obtaining a three-dimensional array with longitude, latitude and altitude as dimensions;

step 9, constructing temperature data of the same period in a plurality of years according to the steps 1 to 8, averaging the temperature data, and establishing a reliable cubic grid of a three-dimensional array as a temperature data resource of an adjacent space of the month in a virtual test;

step 10, performing the same processing on the density and the pressure according to the steps 1 to 9 to obtain a cubic grid of a three-dimensional array of the density and the pressure;

step 11, writing the constructed data including temperature, density and pressure into a text file T;

step 12D, reading the atmospheric environment data from the text file T written with the data, expressing the data according to three specifications of DRM, SRM and EDCS of the SEDRIS, and storing the atmospheric environment data into a SEDRIS standard format; the SEDRIS is a comprehensive environment data representation and exchange specification; DRM is a data representation model; SRM is a space reference model; the EDCS is an environmental data coding specification.

2. The method for constructing the virtual environment resource in the atmosphere near the space according to claim 1, further comprising a construction process of an atmospheric horizontal wind field, which comprises the following specific processes:

step 12A, determining a corresponding atmospheric horizontal wind field calculation formula according to the latitude range:

firstly, calculating the wind-turning according to the formula (1) in a latitude range of 15-80 degrees:

wherein, P is air pressure, and rho is atmospheric density;referred to as the turn parameter, omega is the turn angular velocity,is the geographic latitude; x is the eastward distance and y is the northward distance; u. of_g、v_gWind is respectively blown in the horizontal latitudinal direction and the longitudinal direction within the range of 15-80 degrees;

and then calculating gradient wind according to the formula (2):

wherein,a is the earth radius; u. of_gr、v_grRespectively horizontal latitudinal direction gradient wind and longitudinal direction gradient wind;

(II) calculating gradient wind according to a formula (3) within a latitude range of 15 degrees S-15 degrees N:

wherein u is_e、v_eThe wind is respectively gradient wind in the horizontal latitudinal direction and the longitudinal direction at an angle of 15 degrees S-15 degrees N;

step 12B, constructing a uniform cubic grid of a three-dimensional array of the atmospheric horizontal wind field: substituting the previously constructed density and pressure three-dimensional array of the area range into a wind field calculation formula to obtain a uniform cubic grid of the three-dimensional array of the wind-turning and gradient wind at the corresponding position, thereby generating atmospheric horizontal wind field resources of the space adjacent to the area;

and 12C, writing the constructed wind field data into the text file T.

3. The method for constructing the virtual environment resource in the near space atmosphere according to claim 1 or 2, wherein the specific process of processing the abnormal data in the longitude, latitude, altitude and temperature data in the step 2 is as follows:

and acquiring the original average resolution of the data aiming at longitude, latitude, altitude and temperature data in each event, then eliminating abnormal data, taking the previous data of the abnormal data as a basis, and adding recursion data which is obtained by recursion by utilizing the resolution on the basis and has the same quantity with the eliminated abnormal data.

4. The method for constructing the virtual environment resource in the atmosphere near the space according to claim 3, wherein the specific process of the coordinate transformation in the step 4 is as follows:

the grid of the finally established coordinate system is in distance units, so the longitude and latitude data acquired in each event are converted into distance data in rectangular coordinates according to the coordinate projection system in Matlab2014a, and are rounded up nearby.

5. The method as claimed in claim 4, wherein k is equal to k'.

6. The method for constructing the virtual environment resource in the near space atmosphere according to claim 5, wherein the specific process of acquiring the historical data of the atmospheric environment resource in the near space in the step 1 is as follows:

selecting an NC file of the atmospheric environment data of the adjacent space corresponding to the specific spatial range and the time range of the SABER; and reading the NC file, and acquiring historical data of the atmospheric environment resources in the adjacent space.

7. The method for constructing the virtual environment resource in the near space atmosphere as claimed in claim 5, wherein k ═ k' ═ 1.

8. The method for constructing the virtual environment resource in the near space atmosphere according to claim 7, wherein the specific process of acquiring the historical data of the atmospheric environment resource in the near space in the step 1 is as follows:

selecting an NC file of the atmospheric environment data of the adjacent space corresponding to the specific spatial range and the time range of the SABER; and reading the NC file, and acquiring historical data of the atmospheric environment resources in the adjacent space.