Method, device, equipment and medium for determining atmospheric density based on near-ground high altitude

Abstract

The present application discloses a method, device, equipment and medium for determining atmospheric density based on the near-ground high layer. The method includes: determining the first target period in which the target moment is located, and the first period before the first target period is the second target period; determining the radiation level based on the pixel value of the pixel point in the solar image of the second target period; determining the radiation index M10.7 at the target moment based on the sum DOY at the target moment; determining the second target period in which the target moment is located, and obtaining the ap index of the first N second periods of the second target period; based on the ap index of the first N second periods, evaluating the geomagnetic disturbance intensity at the target moment in real time, and obtaining the Ap* index at the target moment; replacing the solar activity index in the input parameters of the atmospheric density estimation model from F10.7 to M10.7, replacing the geomagnetic activity index from Ap index to Ap* index, and determining the atmospheric density of the near-ground high layer at the target moment based on M10.7 and Ap* index. According to the present application, the accuracy of atmospheric density estimation can be improved.

Description

Technical Field

[0001] The present application relates to the field of space environment detection technology, in particular to the field of atmospheric density detection technology, and in particular to a method, device, equipment and medium for determining atmospheric density based on the near-earth high layer. Background Art

[0002] The precise calculation of the density of the upper atmosphere near the Earth is crucial for the precise calculation of the orbits of low-orbit and ultra-low-orbit near-Earth satellites. As my country's ultra-low-orbit satellite constellation progresses, the importance of this technology is becoming increasingly prominent.

[0003] At present, the scheme for estimating the density of the upper atmosphere near the Earth mainly relies on the atmospheric density estimation model, which estimates the atmospheric density through historical data and observation data. Specifically, by inputting various input parameters such as solar activity indicators and geomagnetic activity indicators into the model, the atmospheric density of the upper atmosphere near the Earth can be estimated. However, the intense solar activity will cause significant changes in the density of the upper atmosphere, which makes the existing scheme for estimating the density of the upper atmosphere near the Earth prone to atmospheric density estimation errors during the estimation. Summary of the invention

[0004] In view of this, the embodiments of the present application provide a method, device, equipment, medium and product for determining atmospheric density based on the near-ground high layer, which can reduce the error of atmospheric density estimation and improve its estimation accuracy.

[0005] In a first aspect, an embodiment of the present application provides a method for determining atmospheric density based on the near-ground upper layer, the method comprising: determining a first target period in which a target moment is located from a plurality of first periods, and obtaining M solar images of the first M first periods of the first target period, wherein the solar image is updated once every first period; for each solar image, determining the radiation level of the sun in the 304 angstrom band based on the pixel value of a pixel point in the solar image , and get the corresponding solar image ; Based on the DOY corresponding to the target time and the corresponding , determine the radiation index M10.7 corresponding to each solar image, and determine the target M10.7 corresponding to the target time by combining the M10.7 corresponding to the M solar images, wherein M10.7 is used to characterize the radiation degree of the sun in the far ultraviolet FUV band and the extreme ultraviolet EUV band; determine the second target period in which the target time is located from multiple second periods, and obtain the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period; based on the ap index of the first N second periods, evaluate the geomagnetic disturbance intensity at the target time in real time, and obtain the sliding geomagnetic Ap* index corresponding to the target time, wherein the ap index and the Ap* index are updated every second period, and the second period is less than or equal to 3h; replace the solar activity index in the input parameters of the atmospheric density estimation model from the solar radio flux F10.7 to the target M10.7, and replace the geomagnetic activity index from the geomagnetic Ap index to the Ap* index, and determine the atmospheric density of the near-earth upper layer at the target time based on the target M10.7 and the Ap* index.

[0006] In some implementations of the first aspect, the radiation level of the sun in the 304 angstrom band is determined based on the pixel value of the pixel point in the solar image. , including: converting the RGB three-channel solar image into a single-channel grayscale image; calculating the average value of the pixel matrix corresponding to the grayscale image to obtain the solar radiation level in the 304 angstrom band .

[0007] In some implementation methods of the first aspect, before converting the RGB three-channel solar image into a single-channel grayscale image, it also includes: setting the pixel values ​​of pixel points within a target range in the solar image to zero; wherein the target range is used to display identification information, and the identification information includes satellite information, band information, and time information.

[0008] In some implementations of the first aspect, based on the DOY corresponding to the target time and the , determine the radiation index M10.7 corresponding to each solar image, and combine the M10.7 corresponding to the M solar images to determine the target M10.7 corresponding to the target time, including: based on the DOY corresponding to the target time and the M10.7 corresponding to each solar image , determine the amount of adjustment of the solar radiation level in the 304 angstrom band ; Based on the amount of adjustment of the solar radiation level in the 304 angstrom band , determine the radiation degree of the sun in the FUV band and the EUV band, and obtain the radiation index M10.7 corresponding to each solar image; determine that the average of M M10.7 corresponding to M solar images is the target M10.7.

[0009] In some implementations of the first aspect, based on the DOY corresponding to the target time and the , determine the amount of adjustment of the solar radiation level in the 304 angstrom band , including: If DOY<183, then ; If DOY>=183, then .

[0010] In some implementations of the first aspect, the amount of adjustment is based on the solar radiation level in the 304 angstrom band. , determine the degree of solar radiation in the FUV band and EUV band, and obtain the radiation index M10.7 corresponding to each solar image, including: ,but ;like ,but .

[0011] In some implementable methods of the first aspect, based on the ap index of the first N second periods, the geomagnetic disturbance intensity at the target time is evaluated in real time to obtain the sliding geomagnetic Ap* index corresponding to the target time, including: calculating the average of the ap index of the first N second periods to obtain the Ap* index corresponding to the target time.

[0012] In a second aspect, an embodiment of the present application provides an atmospheric density determination device based on the near-ground upper layer, the device comprising: an acquisition module, used to determine the first target period in which the target moment is located from multiple first periods, and to obtain M solar images of the first M first periods of the first target period, wherein the solar image is updated once every first period; a determination module, used to determine the radiation level of the sun in the 304 angstrom band for each solar image based on the pixel value of the pixel point in the solar image , and get the corresponding solar image ; Determination module, also used based on the DOY corresponding to the target time and each solar image corresponding , determine the radiation index M10.7 corresponding to each solar image, and combine the M10.7 corresponding to the M solar images to determine the target M10.7 corresponding to the target time, wherein M10.7 is used to characterize the radiation degree of the sun in the far ultraviolet FUV band and the extreme ultraviolet EUV band; the acquisition module is also used to determine the second target period where the target time is located from multiple second periods, and obtain the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period; the evaluation module is used to evaluate the ap index based on the first N second periods , real-time evaluation of the geomagnetic disturbance intensity at the target time, and obtaining the sliding geomagnetic Ap* index corresponding to the target time, wherein the ap index and the Ap* index are updated every second period, and the second period is less than or equal to 3h; the determination module is also used to replace the solar activity index in the input parameters of the atmospheric density estimation model from the solar radio flux F10.7 to the target M10.7, and replace the geomagnetic activity index from the geomagnetic Ap index to the Ap* index, and determine the atmospheric density of the near-earth upper layer at the target time based on the target M10.7 and the Ap* index.

[0013] In some implementations of the second aspect, the determination module is specifically used to include: converting the RGB three-channel solar image into a single-channel grayscale image; calculating the average value of the pixel matrix corresponding to the grayscale image to obtain the solar radiation level in the 304 angstrom band .

[0014] In some implementations of the second aspect, it also includes: a pixel adjustment module, which is used to set the pixel values ​​of pixel points within a target range in the solar image to zero before converting the RGB three-channel solar image into a single-channel grayscale image; wherein the target range is used to display identification information, and the identification information includes satellite information, band information and time information.

[0015] In some implementations of the second aspect, the determining module includes: a first determining unit for determining the DOY corresponding to the target time and the DOY corresponding to each solar image based on the DOY corresponding to the target time. , determine the amount of adjustment of the solar radiation level in the 304 angstrom band ; A second determination unit for adjusting the amount based on the solar radiation level in the 304 angstrom band , determine the radiation degree of the sun in the FUV band and the EUV band, and obtain the radiation index M10.7 corresponding to each solar image; the third determination unit is used to determine that the average of M M10.7 corresponding to M solar images is the target M10.7.

[0016] In some implementations of the second aspect, the first determining unit is specifically configured to: if DOY<183, then ; If DOY>=183, then .

[0017] In some implementations of the second aspect, the second determining unit is specifically configured to: ,but ;like ,but .

[0018] In some implementations of the second aspect, the evaluation module is specifically used to: calculate the average of the ap indexes of the first N second periods to obtain the Ap* index corresponding to the target time.

[0019] In a third aspect, an embodiment of the present application provides an electronic device, comprising: a processor and a memory storing computer program instructions; when the processor executes the computer program instructions, the steps of the method for determining atmospheric density based on the near-ground high layer as in the first aspect are implemented.

[0020] In a fourth aspect, an embodiment of the present application provides a computer-readable storage medium having computer program instructions stored thereon. When the computer program instructions are executed by a processor, the steps of the method for determining atmospheric density based on the near-ground upper layer as in the first aspect are implemented.

[0021] In a fifth aspect, an embodiment of the present application provides a computer program product, which is stored in a non-volatile storage medium. When the computer program product is executed by a processor, it implements the steps of the method for determining atmospheric density based on the near-ground high layer as in the first aspect.

[0022] In a sixth aspect, an embodiment of the present application provides a chip, which includes a processor and a communication interface, wherein the communication interface and the processor are coupled, and the processor is used to run programs or instructions to implement the steps of the method for determining atmospheric density based on near-ground high altitudes as in the first aspect.

[0023] The present application provides a method, device, equipment, medium and product for determining atmospheric density based on the upper layer near the ground, determining the first target period in which the target time is located, and obtaining M solar images of the first M first periods before the first target period, and determining the radiation level of the sun in the 304 angstrom band based on the pixel value of each pixel point in the solar image ; Based on the corresponding The DOY at the target time is used to determine the degree of solar radiation in the far ultraviolet FUV band and the extreme ultraviolet EUV band, and the radiation index M10.7 of each solar image is obtained. The M10.7 of M solar images is combined for comprehensive evaluation to obtain the target M10.7 corresponding to the target time. The M10.7 index represents the radiation situation in the EUV band. The atmosphere at an altitude of about 400km has a higher absorption rate for electromagnetic waves in this band, so using M10.7 as a solar radiation parameter to predict the atmospheric density at an altitude of about 400km is more accurate than using F10.7. Determine the second target period where the target time is located, and obtain the ap index of the first N second periods of the second target period; Based on the ap index of the first N second periods, evaluate the geomagnetic disturbance intensity at the target time in real time, and obtain the sliding geomagnetic Ap* index corresponding to the target time; Replace the solar activity index in the input parameters of the atmospheric density estimation model from F10.7 to the target M10.7, and replace the geomagnetic activity index from the Ap index to the Ap* index, and determine the atmospheric density of the near-ground high layer at the target time based on the target M10.7 and Ap* index. In this way, compared with the daily frequency F10.7 and Ap, the time granularity of the sliding geomagnetic index Ap* and the high-frequency M10.7 index is finer, making the model more sensitive to changes in solar activity and more accurately reflecting the real-time status of solar activity. Therefore, using the target M10.7 and Ap* index can improve the accuracy of the estimation of the atmospheric density of the near-ground high layer at the target time. BRIEF DESCRIPTION OF THE DRAWINGS

[0024] In order to more clearly illustrate the technical solution of the embodiments of the present application, the following briefly introduces the drawings in the embodiments of the present application.

[0025] Figure 1 It is a flowchart of a method for determining atmospheric density based on a near-ground high altitude provided in an embodiment of the present application;

[0026] Figure 2 is an exemplary schematic diagram of a solar image provided by an embodiment of the present application;

[0027] Figure 3 is an exemplary schematic diagram of an algorithm interface of an MSIS20 model provided in an embodiment of the present application;

[0028] Figure 4 It is a structural schematic diagram of an atmospheric density determination device based on the near-ground high layer provided in an embodiment of the present application;

[0029] Figure 5 It is a schematic diagram of the hardware structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0030] The principles and spirit of the present application will be described below with reference to several exemplary embodiments. It should be understood that the purpose of providing these embodiments is to make the principles and spirit of the present application clearer and more thorough, so that those skilled in the art can better understand and implement the principles and spirit of the present application. The exemplary embodiments provided herein are only part of the embodiments of the present application, not all of the embodiments. Based on the embodiments herein, all other embodiments obtained by ordinary technicians of the art without creative work are within the scope of protection of this application.

[0031] In this document, terms such as first, second, third, etc. are only used to distinguish one entity (or operation) from another entity (or operation), but not to require or imply any order or relationship between these entities (or operations).

[0032] The following is a brief description of the concepts and technical terms that may be involved in the embodiments of the present application.

[0033] Solar activity index: refers to the index used to describe solar activity and reflect the intensity of solar activity, such as F10.7 (radiation) index and M10.7 (radiation) index.

[0034] F10.7: is an index used to describe the intensity of solar radiation. It represents the solar radio flux at a wavelength of 10.7 cm, and therefore can reflect the intensity of solar activity. F10.7 is used as a good representative of solar ultraviolet radiation. The measurement of F10.7 is the average value of the intensity of solar radio radiation measured within one hour on the ground at a 100 MHz frequency band centered at 2800 MHz, and a value is updated daily. F10.7 is closely related to the activity on the surface of the sun, and also has a long and short periodic variation law of 11 years and 27 days. It is an important index to characterize solar activity, and there is a strong correlation between geomagnetic field disturbances and the activity level of the sun.

[0035] M10.7: is an index used to describe solar activity. It is calculated from the ratio of the core to the wing of the Mg II 280 nm band. This ratio reflects the emission of solar active regions in the chromosphere and is theoretically independent of instrument sensitivity changes. The M10.7 index is used as a good proxy for solar far ultraviolet (FUV) and extreme ultraviolet (EUV) radiation, which have a significant impact on the Earth's atmosphere, especially the thermosphere (located between about 85 km and 600 km above the Earth's surface).

[0036] SDO: Solar Dynamics Observatory, launched in 2010, is designed to monitor the activity and changes of the sun.

[0037] 304Å Solar Image: A solar image taken by the SDO satellite in the 304Å band, mainly showing the specific structure and activity of the Sun in the extreme ultraviolet (EUV) band.

[0038] Geomagnetic activity indicators: such as the Kp index and the (geomagnetic) Ap index, which reflect the extent to which the Earth's magnetic field is affected by the solar wind.

[0039] ap index: also known as the equivalent planetary amplitude of geomagnetic disturbance, the unit is 2nT, a new ap value is generated every 3 hours, and 8 ap indexes are generated every day in chronological order. According to relevant standards, an ap value greater than 39 corresponds to a magnetic storm period. The maximum value of the ap index is 400.

[0040] Ap index: It is an index of the global daily geomagnetic disturbance intensity, also known as the planetary equivalent daily amplitude. It is an indicator to measure the level of geomagnetic activity. A new Ap value is generated every day, which is the average of 8 AP values ​​every day. Therefore, the Ap index can only reflect the overall situation of geomagnetic disturbance throughout the day.

[0041] Ap* index: It is an index of the intensity of periodic geomagnetic disturbances and a real-time indicator for measuring the level of geomagnetic activity. Compared with the Ap index, it is more real-time and can therefore more accurately assess the overall situation of geomagnetic disturbances at a certain moment.

[0042] The MSIS (Mass Spectrometer and Incoherent Scatter Radar Extended) model is one of the widely used empirical models, which can provide atmospheric physical properties such as neutral atmospheric density and temperature from the ground to the thermosphere (0~1000km).

[0043] In related technologies, the schemes for predicting the density of the upper atmosphere near the ground mainly include empirical models and theoretical models. The empirical models predict the atmospheric density through historical data and observational data, while the theoretical models rely on the understanding of atmospheric physical processes.

[0044] The input parameters of the existing atmospheric density prediction model include: longitude and latitude, DOY, F10.7 index, and Ap index. When the MSIS series model calculates atmospheric density based on the above input parameters, the average error during the solar calm period is 15%, and it can reach 25%-48% during the outburst period.

[0045] Some of the reasons for this large discrepancy are:

[0046] The factor that has the greatest impact on the accurate prediction of atmospheric density is the activity of the sun. In the existing prediction model, the parameters corresponding to MSIS and solar activity are F10.7 and Ap, which correspond to the intensity of solar radiation and the intensity of geomagnetic disturbances caused by the solar wind, respectively, and both have a daily value. In reality, both solar radiation and geomagnetic disturbances may change dramatically during the day, which means that F10.7 and Ap can fluctuate greatly during the day. Therefore, the coarse time granularity of the input parameters leads to large errors in the prediction results, which in turn causes large errors in the atmospheric density prediction.

[0047] At present, most models use the F10.7 parameter to represent the solar radiation intensity, but in reality, F10.7 cannot fully represent the solar radiation intensity. First, solar radiation includes electromagnetic waves of various bands. The electromagnetic wave fluxes of various bands radiated by the sun are not necessarily proportional at the same time, and there may be a lag relationship between the radiation of different bands; secondly, the absorption rate of electromagnetic waves of different bands at different altitudes is completely different, which leads to completely different heating effects of electromagnetic waves of different bands of solar radiation on the atmosphere at different altitudes. In summary, the atmospheric density prediction model only uses F10.7 as the input parameter of solar radiation intensity. The error caused by using F10.7 to fully represent the radiation intensity will inevitably cause density prediction errors at different atmospheric altitudes.

[0048] The following is a detailed description of the method for determining the atmospheric density based on the near-ground high layer provided in the embodiment of the present application through specific embodiments and their application scenarios in conjunction with the accompanying drawings.

[0049] Figure 1 It is a flow chart of a method for determining atmospheric density based on a high altitude near the ground provided in one embodiment of the present application. The executor of the method for determining atmospheric density based on a high altitude near the ground may be an atmospheric density prediction device.

[0050] The following takes the atmospheric density prediction device as an example to illustrate the atmospheric density determination method based on the near-ground high-level. It should be noted that the above-mentioned execution subject and application scenario do not constitute a limitation on the present application.

[0051] like Figure 1 As shown, the method for determining atmospheric density based on the near-ground high layer provided in the embodiment of the present application may include steps 110 to 160.

[0052] Step 110, determining a first target period in which the target time is located from a plurality of first periods, and acquiring M solar images of the first M first periods of the first target period;

[0053] Step 120, for each solar image, determine the solar radiation level in the 304 angstrom band based on the pixel values ​​of the pixels in the solar image. , and get the corresponding solar image ;

[0054] Step 130, based on the DOY corresponding to the target time and the , determine the radiation index M10.7 corresponding to each solar image, and combine the M10.7 corresponding to the M solar images to determine the target M10.7 corresponding to the target time;

[0055] Step 140, determining the second target period in which the target time is located from the plurality of second periods, and obtaining the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period;

[0056] Step 150, based on the ap indexes of the first N second periods, the geomagnetic disturbance intensity at the target time is evaluated in real time to obtain the sliding geomagnetic Ap* index corresponding to the target time;

[0057] Step 160, replace the solar activity index in the input parameters of the atmospheric density estimation model from the solar radio flux F10.7 to the target M10.7, and replace the geomagnetic activity index from the geomagnetic Ap index to the Ap* index, and determine the atmospheric density of the upper layer near the ground at the target time based on the target M10.7 and the Ap* index.

[0058] The method for determining atmospheric density based on the upper layer near the ground provided in the embodiment of the present application determines the first target period in which the target time is located, and obtains M solar images of the first M first periods before the first target period, and determines the radiation level of the sun in the 304 angstrom band based on the pixel value of the pixel point in each solar image. ; Based on the corresponding The DOY at the target time is used to determine the degree of solar radiation in the far ultraviolet FUV band and the extreme ultraviolet EUV band, and the radiation index M10.7 of each solar image is obtained. The M10.7 of M solar images is combined for comprehensive evaluation to obtain the target M10.7 corresponding to the target time. The M10.7 index represents the radiation situation in the EUV band. The atmosphere at an altitude of about 400km has a higher absorption rate for electromagnetic waves in this band, so using M10.7 as a solar radiation parameter to predict the atmospheric density at an altitude of about 400km is more accurate than using F10.7. Determine the second target period where the target time is located, and obtain the ap index of the first N second periods of the second target period; Based on the ap index of the first N second periods, evaluate the geomagnetic disturbance intensity at the target time in real time, and obtain the sliding geomagnetic Ap* index corresponding to the target time; Replace the solar activity index in the input parameters of the atmospheric density estimation model from F10.7 to the target M10.7, and replace the geomagnetic activity index from the Ap index to the Ap* index, and determine the atmospheric density of the near-ground high layer at the target time based on the target M10.7 and Ap* index. In this way, compared with the daily frequency F10.7 and Ap, the time granularity of the sliding geomagnetic index Ap* and the high-frequency M10.7 index is finer, making the model more sensitive to changes in solar activity and more accurately reflecting the real-time status of solar activity. Therefore, using the target M10.7 and Ap* index can improve the accuracy of the estimation of the atmospheric density of the near-ground high layer at the target time.

[0059] The specific implementation of the above steps will be described in detail below in conjunction with specific embodiments.

[0060] Involving step 110, determining the first target period in which the target moment is located from multiple first periods, and acquiring M solar images of the first M first periods of the first target period.

[0061] In step 110, the first cycle is an update cycle of the solar image, which is an image of the sun in the 304 angstrom band taken by the SDO satellite. The first target cycle is the first cycle where the target time is located, and M is a positive integer.

[0062] In step 120, for each solar image, based on the pixel value of the pixel point in the solar image, the radiation level of the sun in the 304 angstrom band is determined. , and get the corresponding solar image .

[0063] In step 120, the solar image is updated every first period, and the first period is less than 12 hours, such as 3 hours, 5 hours, etc. The solar image is a 512×512 image with RGB three channels, such as Figure 2 As shown, the lower left area of ​​the image contains satellite, corresponding wave and time information.

[0064] In some embodiments of the present application, the above step 120 may specifically include:

[0065] Convert the RGB three-channel sun image to a single-channel grayscale image;

[0066] Calculate the average value of the pixel matrix corresponding to the grayscale image to obtain the solar radiation level in the 304 angstrom band .

[0067] Specifically, the RGB three-channel solar image is converted into a single-channel grayscale image. At this time, the size of the pixel matrix is ​​converted from 512×512×3 to 512×512×1. The average value of the pixel matrix of the single-channel grayscale image is calculated based on formula (1) as the solar radiation level at 304Å. .

[0068] (1)

[0069] in, Represents the value of the i-th pixel from the left and the j-th pixel from the top in the pixel matrix of the grayscale image, where i and j both start from 0.

[0070] In this way, the pixel values ​​of the pixels in the single-channel grayscale image of the solar image can be used to accurately calculate the solar radiation level in the 304 angstrom band. .

[0071] In some embodiments of the present application, in order to further improve the radiation level Before converting the RGB three-channel sun image into a single-channel grayscale image, the following steps may be included:

[0072] Set the pixel values ​​of the pixels within the target range in the sun image to zero.

[0073] The target range is used to display identification information, which includes satellite information, band information, and time information. The target range corresponds to 24 pixels extending upward from the lower left corner of the image and 256 pixels extending to the right. The number of pixels in the target range is 256×24×3.

[0074] In this way, by setting the pixel values ​​of the pixels within the target range to zero and the identification information to black, the pixel values ​​of the pixels in the single-channel grayscale image of the sun image are calculated. When the identification information is eliminated, the influence of the calculation results is eliminated, and the radiation level of the 304 angstrom band is further improved. The calculation accuracy.

[0075] Involving step 130, based on the DOY corresponding to the target time and the corresponding , determine the radiation index M10.7 corresponding to each solar image, and combine the M10.7 corresponding to the M solar images to determine the target M10.7 corresponding to the target time.

[0076] In step 130, the DOY (day of year) of the target time is the day of the year at which the target time is located. Based on the DOY, an M10.7 can be calculated in combination with each solar image. M solar images can obtain M M10.7s. The M10.7s are averaged to obtain the target M10.7 at the target time. The M10.7 corresponding to each solar image is the M10.7 corresponding to each first period.

[0077] In some embodiments of the present application, the above step 130 may specifically include the following steps:

[0078] Based on the DOY corresponding to the target time and the corresponding , determine the amount of adjustment of the solar radiation level in the 304 angstrom band ;

[0079] Adjustment based on the solar radiation level at 304 angstroms , determine the degree of solar radiation in the FUV band and EUV band, and obtain the radiation index M10.7 corresponding to each solar image;

[0080] Determine the mean of the M M10.7s corresponding to the M solar images as the target M10.7.

[0081] In the embodiment of the present application, according to each sun image corresponding to and the DOY at the target time, the amount of adjustment of the solar radiation level in the 304 angstrom band can be accurately determined Based on this, the amount of radiation level adjustment is based on the solar image , the radiation index M10.7 corresponding to the solar image can be calculated. The M10.7 index represents the radiation situation in the EUV band. The atmosphere at an altitude of about 400km has a higher absorption rate for electromagnetic waves in this band. Therefore, using M10.7 as a solar radiation parameter to predict the atmospheric density at an altitude of about 400km is more accurate than using F10.7. Since M10.7 has a lag in the assessment of atmospheric density, the M10.7 of the first cycle before the target time is combined, and the average of the M historical M10.7 is input as an input parameter to the atmospheric density estimation model, which provides a reliable basis for the model's estimation and improves the accuracy of atmospheric density estimation.

[0082] In some embodiments of the present application, the above-mentioned DOY corresponding to the target time and each solar image corresponding to , determine the amount of adjustment of the solar radiation level in the 304 angstrom band , which may specifically include:

[0083] If DOY<183, then ;

[0084] If DOY>=183, then .

[0085] In some embodiments of the present application, the amount of adjustment is based on the solar radiation level in the 304 angstrom band. , determine the degree of solar radiation in the FUV band and EUV band, and obtain the radiation index M10.7 corresponding to each solar image, which can specifically include:

[0086] like ,but ;

[0087] like ,but .

[0088] In the embodiment of the present application, the time granularity of M10.7 can be refined based on the update cycle of the solar image. For example, if the solar image provided by SDO is updated every 3 hours, the M10.7 index calculated based on the 304Å solar image can achieve a time granularity of 3 hours. Similarly, the target M10.7 calculated based on M10.7 can achieve a time granularity of 3 hours. In this way, replacing the daily F10.7 index with the target M10.7 index with a higher update frequency can more accurately reflect the real-time status of solar activity and improve the accuracy of the calculation of the atmospheric density of the near-earth high-level at the target time.

[0089] Involving step 140, determining the second target period in which the target moment is located from multiple second periods, and obtaining the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period.

[0090] In step 140, the second period is an update period of the AP index, the second period is less than or equal to 3 hours, and the second target period is the second period where the target time is located.

[0091] Involving step 150, based on the ap indexes of the first N second periods, the geomagnetic disturbance intensity at the target time is evaluated in real time to obtain the sliding geomagnetic Ap* index corresponding to the target time.

[0092] In step 150, N×the duration of the second cycle≤24h, and the value of N can be specifically set according to the duration of the second cycle, which is not specifically limited in this application. The Ap* index corresponding to the target time is the Ap* index of the second cycle where the target time is located.

[0093] Since the Ap* index at any moment is obtained based on the first N second periods of the second period at that moment, the Ap* index can also be updated periodically. Both the Ap* index and the Ap* index are updated every second period, and the second period is less than or equal to 3h.

[0094] In some embodiments of the present application, the above step 150 evaluates the geomagnetic disturbance intensity at the target time in real time based on the ap index of the first N second periods to obtain the sliding geomagnetic Ap* index corresponding to the target time, which may specifically include:

[0095] Calculate the average of the ap indexes of the first N second periods to obtain the Ap* index corresponding to the target time.

[0096] Specifically, the Ap* index corresponding to the target time t can be calculated using formula (2):

[0097] (2)

[0098] in, represents the sliding Ap index at the target time t, Represents the AP index of the second cycle (i.e., the second target cycle) at time t, It represents the ap index of the second period of the i-th period before time t. If the update frequency of the ap index is 3h, then Ap* will be updated every 3h.

[0099] In one example, the second cycle has a duration of 1 hour, N=24, then .

[0100] In another example, the second cycle has a cycle length of 3 hours, N=8, then .

[0101] In the embodiment of the present application, compared with the daily updated geomagnetic index Ap, the sliding geomagnetic index Ap* is updated faster and therefore more real-time, and can better reflect the real-time geomagnetic disturbance. Therefore, using Ap* to replace Ap for atmospheric density estimation can more accurately reflect the real-time state of solar activity and improve the accuracy of the estimated atmospheric density in the near-earth high layer at the target time.

[0102] Involving step 160, the solar activity index in the input parameters of the atmospheric density estimation model is replaced by the solar radio flux F10.7 with the target M10.7, and the geomagnetic activity index is replaced by the geomagnetic Ap index with the Ap* index. Based on the target M10.7 and the Ap* index, the atmospheric density of the upper layer near the ground at the target time is determined.

[0103] In step 160, the atmospheric density estimation model may be an MSIS series model such as the MSIS20 model or a DTM model, etc., and this application does not make any specific limitation on this.

[0104] As a specific example, the atmospheric density estimation model may be the MSIS20 model. The input parameters and output parameters of the MSIS20 model are as follows: Figure 3 As shown, the algorithm interface calc obtains various input parameters by calling instructions. The present application replaces the Ap index in the input parameter with the Ap* index and inputs it into the interface Ap(1), and replaces the daily F10.7 index with the target M10.7 index with a higher update frequency and inputs it into the interface SFLUX to output atmospheric density and temperature with higher accuracy.

[0105] In addition, experimental data show that when the first period is set to 12 minutes, the second period is set to 3 hours, and N is set to 8, the same model MSIS20 is used in this application. Only by optimizing the input parameters, replacing F10.7 with the target M10.7, and replacing the Ap index with the Ap* index, the prediction error of the 400km high calm period is reduced from 15% to 9%, which effectively improves the accuracy of the estimated atmospheric density.

[0106] Corresponding to the method embodiment of the present application, the present application also provides an atmospheric density determination device based on the near-ground high altitude.

[0107] Figure 4 Schematic diagram of the structure of an atmospheric density determination device based on the near-ground high-altitude provided in an embodiment of the present application. Figure 4 As shown, the atmospheric density determination device 400 based on the near-ground upper layer may include: an acquisition module 410 , a determination module 420 , and an evaluation module 430 .

[0108] The acquisition module 410 is used to determine the first target period in which the target time is located from multiple first periods, and obtain M solar images of the first M first periods of the first target period, wherein the solar image is updated once every first period; the determination module 420 is used to determine the radiation level of the sun in the 304 angstrom band for each solar image based on the pixel value of the pixel point in the solar image. , and get the corresponding solar image ; Determine module 420, also used based on the target time corresponding to the DOY and each solar image corresponding , determine the radiation index M10.7 corresponding to each solar image, and determine the target M10.7 corresponding to the target time by combining the M10.7 corresponding to the M solar images, wherein M10.7 is used to characterize the radiation degree of the sun in the far ultraviolet FUV band and the extreme ultraviolet EUV band; the acquisition module 410 is further used to determine the second target period in which the target time is located from multiple second periods, and obtain the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period; the evaluation module 430 is used to evaluate the ap index of the first N second periods based on the ap index of the first N second periods. The module 420 is used to evaluate the intensity of the geomagnetic disturbance at the target time in real time, and obtain the sliding geomagnetic Ap* index corresponding to the target time, wherein the Ap index and the Ap* index are updated every second period, and the second period is less than or equal to 3h; the module 420 is also used to replace the solar activity index in the input parameters of the atmospheric density estimation model from the solar radio flux F10.7 to the target M10.7, and replace the geomagnetic activity index from the geomagnetic Ap index to the Ap* index, and determine the atmospheric density of the near-earth upper layer at the target time based on the target M10.7 and the Ap* index.

[0109] The atmospheric density determination device based on the near-ground upper layer provided in the embodiment of the present application determines the first target period in which the target time is located, and obtains M solar images of the first M first periods before the first target period, and determines the radiation level of the sun in the 304 angstrom band based on the pixel value of the pixel point in each solar image. ; Based on the corresponding The DOY at the target time is used to determine the degree of solar radiation in the far ultraviolet FUV band and the extreme ultraviolet EUV band, and the radiation index M10.7 of each solar image is obtained. The M10.7 of M solar images is combined for comprehensive evaluation to obtain the target M10.7 corresponding to the target time. The M10.7 index represents the radiation situation in the EUV band. The atmosphere at an altitude of about 400km has a higher absorption rate for electromagnetic waves in this band, so using M10.7 as a solar radiation parameter to predict the atmospheric density at an altitude of about 400km is more accurate than using F10.7. Determine the second target period where the target time is located, and obtain the ap index of the first N second periods of the second target period; Based on the ap index of the first N second periods, evaluate the geomagnetic disturbance intensity at the target time in real time, and obtain the sliding geomagnetic Ap* index corresponding to the target time; Replace the solar activity index in the input parameters of the atmospheric density estimation model from F10.7 to the target M10.7, and replace the geomagnetic activity index from the Ap index to the Ap* index, and determine the atmospheric density of the near-ground high layer at the target time based on the target M10.7 and Ap* index. In this way, compared with the daily frequency F10.7 and Ap, the time granularity of the sliding geomagnetic index Ap* and the high-frequency M10.7 index is finer, making the model more sensitive to changes in solar activity and more accurately reflecting the real-time status of solar activity. Therefore, using the target M10.7 and Ap* index can improve the accuracy of the estimation of the atmospheric density of the near-ground high layer at the target time.

[0110] In some embodiments of the present application, the determination module is specifically used to include: converting the RGB three-channel solar image into a single-channel grayscale image; calculating the average value of the pixel matrix corresponding to the grayscale image to obtain the solar radiation level in the 304 angstrom band .

[0111] In some embodiments of the present application, it also includes: a pixel adjustment module, which is used to set the pixel values ​​of pixel points within a target range in the solar image to zero before converting the RGB three-channel solar image into a single-channel grayscale image; wherein the target range is used to display identification information, and the identification information includes satellite information, band information and time information.

[0112] In some embodiments of the present application, the determination module includes: a first determination unit for determining the DOY corresponding to the target time and the DOY corresponding to each solar image. , determine the amount of adjustment of the solar radiation level in the 304 angstrom band ; A second determination unit for adjusting the amount based on the solar radiation level in the 304 angstrom band , determine the radiation degree of the sun in the FUV band and the EUV band, and obtain the radiation index M10.7 corresponding to each solar image; the third determination unit is used to determine that the average of M M10.7 corresponding to M solar images is the target M10.7.

[0113] In some embodiments of the present application, the first determining unit is specifically configured to: if DOY<183, then ; If DOY>=183, then .

[0114] In some embodiments of the present application, the second determining unit is specifically configured to: ,but ;like ,but .

[0115] In some embodiments of the present application, the evaluation module is specifically used to: calculate the average of the ap indexes of the first N second periods to obtain the Ap* index corresponding to the target time.

[0116] The atmospheric density determination device based on the near-ground high-altitude provided in the embodiment of the present application can achieve Figure 1-3 The various processes implemented by the service platform in the method embodiment can achieve the same technical effect. To avoid repetition, they will not be described here.

[0117] Figure 5 It is a schematic diagram of the hardware structure of an electronic device provided in an embodiment of the present application.

[0118] like Figure 5 As shown, the electronic device 500 includes a memory 501 , a processor 502 , and a computer program stored in the memory 501 and executable on the processor 502 .

[0119] In an example, the processor 502 may include a central processing unit (CPU), or an application specific integrated circuit (ASIC), or may be configured to implement one or more integrated circuits of the embodiments of the present application.

[0120] The memory 501 may include a read-only memory (ROM), a random access memory (RAM), a disk storage medium device, an optical storage medium device, a flash memory device, an electrical, optical or other physical / tangible memory storage device. Therefore, generally, the memory includes one or more tangible (non-transitory) computer-readable storage media (e.g., a memory device) encoded with software including computer executable instructions, and when the software is executed (e.g., by one or more processors), it can be operated to perform the operations described with reference to the card opening method in the embodiment according to the first aspect of the present application.

[0121] The processor 502 runs a computer program corresponding to the executable program code by reading the executable program code stored in the memory 501, so as to implement the card activation method in the embodiment of the first aspect.

[0122] In some examples, the electronic device 500 may further include a communication interface 503 and a bus 510. Figure 5 As shown, the memory 501, the processor 502, and the communication interface 503 are connected via a bus 510 and communicate with each other.

[0123] The communication interface 503 is mainly used to realize the communication between the modules, devices, units and / or equipment in the embodiment of the present application. The communication interface 503 can also be used to access input devices and / or output devices.

[0124] The bus 510 includes hardware, software, or both, coupling the components of the electronic device 500 to each other. By way of example and not limitation, the bus 510 may include an Accelerated Graphics Port (AGP) or other graphics bus, an Enhanced Industry Standard Architecture (EISA) bus, a Front Side Bus (FSB), a Hyper Transport (HT) interconnect, an Industry Standard Architecture (ISA) bus, an InfiniBand interconnect, a Low Pin Count (LPC) bus, a Memory bus, a Micro Channel Architecture (MCA) bus, a Peripheral Component Interconnect (PCI) bus, a PCI-Express (PCI-E) bus, a Serial Advanced Technology Attachment (SATA) bus, a Video Electronics Standards Association Local Bus (VLB) bus, or other suitable buses or a combination of two or more of these. Where appropriate, the bus 510 may include one or more buses. Although embodiments of the present application describe and illustrate a particular bus, the present application contemplates any suitable bus or interconnect.

[0125] The electronic device provided in the embodiment of the present application can realize Figure 1-3 The various processes implemented by the electronic device in the method embodiment can achieve the same technical effect, and to avoid repetition, they will not be described here.

[0126] In combination with the method for determining atmospheric density based on the upper layer near the ground in the above-mentioned embodiment, the embodiment of the present application can provide a computer storage medium for implementation. The computer storage medium stores computer program instructions; when the computer program instructions are executed by the processor, the steps of any one of the methods for determining atmospheric density based on the upper layer near the ground in the above-mentioned embodiment are implemented.

[0127] In combination with the method for determining atmospheric density based on a high altitude near the ground in the above-mentioned embodiment, the embodiment of the present application may provide a computer program product for implementation. The (computer) program product is stored in a non-volatile storage medium, and when the program product is executed by at least one processor, the steps of any one of the methods for determining atmospheric density based on a high altitude near the ground in the above-mentioned embodiment are implemented.

[0128] An embodiment of the present application further provides a chip, which includes a processor and a communication interface, wherein the communication interface and the processor are coupled, and the processor is used to run programs or instructions to implement the various processes of the above-mentioned embodiment of the method for determining atmospheric density based on the near-ground high layer, and can achieve the same technical effect. To avoid repetition, it will not be repeated here.

[0129] It should be understood that the chip mentioned in the embodiments of the present application can also be called a system-level chip, a system chip, a chip system or a system-on-chip chip, etc.

[0130] It should be clear that the present application is not limited to the specific configuration and processing described above and shown in the figures. For the sake of simplicity, a detailed description of the known method is omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of the present application is not limited to the specific steps described and shown, and those skilled in the art can make various changes, modifications and additions, or change the order between the steps after understanding the spirit of the present application.

[0131] The functional blocks shown in the structural block diagram described above can be implemented as hardware, software, firmware or a combination thereof. When implemented in hardware, it can be, for example, an electronic circuit, an application-specific integrated circuit (ASIC), appropriate firmware, a plug-in, a function card, etc. When implemented in software, the elements of the present application are programs or code segments used to perform the required tasks. The program or code segment can be stored in a machine-readable medium, or transmitted on a transmission medium or a communication link via a data signal carried in a carrier. "Machine-readable medium" may include any medium capable of storing or transmitting information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, ROMs, flash memories, erasable ROMs (EROMs), floppy disks, CD-ROMs, optical disks, hard disks, optical fiber media, radio frequency (RF) links, etc. The code segment can be downloaded via a computer network such as the Internet, an intranet, etc.

[0132] It should also be noted that the exemplary embodiments mentioned in this application describe some methods or systems based on a series of steps or devices. However, this application is not limited to the order of the above steps, that is, the steps can be performed in the order mentioned in the embodiment, or in a different order from the embodiment, or several steps can be performed simultaneously.

[0133] Aspects of the present disclosure are described above with reference to the flowcharts and / or block diagrams of the methods, devices (systems) and computer program products according to the embodiments of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device to produce a machine so that these instructions executed by the processor of the computer or other programmable data processing device enable the implementation of the functions / actions specified in one or more boxes of the flowchart and / or block diagram. Such a processor can be, but is not limited to, a general-purpose processor, a special-purpose processor, a special application processor, or a field programmable logic circuit. It can also be understood that each box in the block diagram and / or flowchart and the combination of boxes in the block diagram and / or flowchart can also be implemented by dedicated hardware that performs the specified function or action, or can be implemented by a combination of dedicated hardware and computer instructions.

[0134] The above is only a specific implementation of the present application. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the systems, modules and units described above can refer to the corresponding processes in the aforementioned method embodiments, and will not be repeated here. It should be understood that the protection scope of the present application is not limited to this. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed in this application, and these modifications or replacements should be included in the protection scope of this application.

Claims

1. A method for determining atmospheric density based on the upper layer near the ground, characterized in that: include: Determine a first target period in which a target moment is located from a plurality of first periods, and obtain M solar images of the first M first periods before the first target period, wherein the solar images are updated every first period; For each solar image, the solar radiation level in the 304 angstrom band is determined based on the pixel values ​​of the pixels in the solar image. , get the corresponding ; Based on the DOY corresponding to the target time and the DOY corresponding to each solar image , determine the radiation index M10.7 corresponding to each solar image, and determine the target M10.7 corresponding to the target time by combining the M10.7 corresponding to the M solar images, wherein M10.7 is used to characterize the radiation degree of the sun in the far ultraviolet FUV band and the extreme ultraviolet EUV band; Determine a second target period in which the target time is located from a plurality of second periods, and obtain the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period; Based on the ap indexes of the first N second periods, the geomagnetic disturbance intensity at the target time is evaluated in real time to obtain the sliding geomagnetic Ap* index corresponding to the target time, wherein the ap index and the Ap* index are updated every second period, and the second period is less than or equal to 3 hours; The solar activity index in the input parameters of the atmospheric density estimation model is replaced by the target M10.7 from the solar radio flux F10.7, and the geomagnetic activity index is replaced by the Ap* index from the geomagnetic Ap index. Based on the target M10.7 and the Ap* index, the atmospheric density of the high altitude near the ground at the target time is determined.

2. The method according to claim 1, characterized in that Based on the pixel values ​​of the pixels in the solar image, the solar radiation level F in the 304 angstrom band is determined. 304 ,include: Convert the RGB three-channel sun image to a single-channel grayscale image; Calculate the average value of the pixel matrix corresponding to the grayscale image to obtain the solar radiation level in the 304 angstrom band .

3. The method according to claim 2, characterized in that Before converting the RGB three-channel sun image into a single-channel grayscale image, the method further includes: Setting the pixel values ​​of the pixels within the target range in the solar image to zero; The target range is used to display identification information, and the identification information includes satellite information, band information and time information.

4. The method according to claim 1, characterized in that Based on the DOY corresponding to the target time and the DOY corresponding to each solar image , determining the radiation index M10.7 corresponding to each solar image, and combining the M10.7 corresponding to the M solar images to determine the target M10.7 corresponding to the target time, including: Based on the DOY corresponding to the target time and the DOY corresponding to each solar image , determine the amount of adjustment of the solar radiation level in the 304 angstrom band ; Adjustment based on the solar radiation level at 304 angstroms , determine the radiation degree of the sun in the FUV band and the EUV band, and obtain the radiation index M10.7 corresponding to each solar image; Determine the average of the M M10.7s corresponding to the M solar images as the target M10.

7.

5. The method according to claim 4, characterized in that The DOY corresponding to the target time and each sun image corresponding to , determine the amount of adjustment of the solar radiation level in the 304 angstrom band ,include: If DOY<183, then ; If DOY>=183, then .

6. The method according to claim 4, characterized in that The amount of adjustment based on the solar radiation level in the 304 angstrom band , determine the degree of solar radiation in the FUV band and the EUV band, and obtain the radiation index M10.7 corresponding to each solar image, including: like ,but ; like ,but .

7. The method according to claim 1, characterized in that The step of evaluating the geomagnetic disturbance intensity at the target time in real time based on the ap indexes of the first N second periods to obtain the sliding geomagnetic Ap* index corresponding to the target time includes: Calculate the average of the ap indexes of the first N second periods to obtain the Ap* index corresponding to the target time.

8. A device for determining atmospheric density based on the upper layer near the ground, characterized in that: include: an acquisition module, used to determine a first target period in which a target moment is located from a plurality of first periods, and acquire M solar images of the first M first periods before the first target period, wherein the solar images are updated once every first period; A determination module is used to determine the radiation level of the sun in the 304 angstrom band for each solar image based on the pixel values ​​of the pixels in the solar image. , get the corresponding ; The determination module is further used to determine the radiation index M10.7 corresponding to each solar image based on the DOY corresponding to the target time and each solar image, and to determine the target M10.7 corresponding to the target time in combination with the M10.7 corresponding to the M solar images, wherein M10.7 is used to characterize the radiation degree of the sun in the far ultraviolet FUV band and the extreme ultraviolet EUV band; The acquisition module is further used to determine the second target period in which the target moment is located from the multiple second periods, and to acquire the geomagnetic disturbance equivalent planetary amplitude ap index of the first N second periods of the second target period; An evaluation module is used to evaluate the geomagnetic disturbance intensity at the target time in real time based on the ap index of the first N second periods, and obtain the sliding geomagnetic Ap* index corresponding to the target time, wherein the ap index and the Ap* index are updated every second period, and the second period is less than or equal to 3 hours; The determination module is also used to replace the solar activity index in the input parameters of the atmospheric density calculation model from the solar radio flux F10.7 to the target M10.7, and to replace the geomagnetic activity index from the geomagnetic Ap index to the Ap* index, and to determine the atmospheric density of the near-earth high layer at the target time based on the target M10.7 and the Ap* index.

9. An electronic device, characterized in that: The electronic device comprises: a processor and a memory storing computer program instructions; when the electronic device executes the computer program instructions, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer program instructions, and when the computer program instructions are executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images