GNSS-R satellite coverage calculation method, electronic device and storage medium

Abstract

The present application provides a method, electronic device and storage medium for calculating the coverage of GNSS-R satellites. The embodiments of the present application obtain GNSS-R satellite data; determine the earth surface grid map based on a preset unit distance, and determine the longitude and latitude range of each grid; map the GNSS-R satellite data to the earth surface grid map according to the longitude and latitude data corresponding to the longitude and latitude range and the mirror reflection point, and obtain the number of revisits and target longitude and latitude data corresponding to each grid; determine the longitude variation range according to the average distance of the earth, the average radius of the earth and the longitude variation function; determine the revisit number matrix according to the number of revisits; determine the average coverage corresponding to each latitude interval based on the longitude variation range, the number of grids that are not 0 in any row of the revisit number matrix, the number of grids included in the longitude variation range and the preset unit distance. The present application can calculate the coverage of GNSS-R satellites and evaluate the detection capability of GNSS-R satellites.

Description

Technical Field

[0001] The present application relates to the field of GNSS-R satellite technology, and in particular to a method for calculating GNSS-R satellite coverage, an electronic device, and a storage medium. Background Art

[0002] The Global Navigation Satellite System-Reflectometry (GNSS-R) technology uses existing GNSS-R satellites and does not require active signal reflection. It has the advantages of miniaturization, low cost, and low power consumption. It is very suitable for being carried by small satellites to form constellations and achieve high temporal and spatial resolution observation of the earth. The main difference between GNSS-R technology and traditional remote sensing technology is that the observation points of GNSS-R are strip-shaped, and the sampling points within the strip are relatively dense. The analysis of the temporal and spatial resolution of GNSS-R satellites is very important for the design of GNSS-R constellations and the evaluation of detection capabilities. Coverage is an important indicator for measuring the temporal and spatial resolution of satellite constellations. It is defined as the ratio of the coverage area of ​​the earth's surface by the constellation observation data within a certain time and at a certain spatial resolution to the total surface area of ​​the earth. The larger the coverage, the higher the temporal and spatial resolution of the constellation and the better the observation effect.

[0003] At present, due to the unique observation system of GNSS-R technology, the calculation of its coverage rate also needs to be specifically considered based on its technical characteristics. There are few methods for calculating coverage rate in the existing technology, and the calculated coverage rate fails to meet the preset standard. The calculation accuracy of GNSS-R satellite coverage rate is low, which is not conducive to the evaluation of the detection capability of GNSS-R satellites. Summary of the invention

[0004] The purpose of the embodiments of the present application is to provide a method for calculating the coverage rate of GNSS-R satellites, an electronic device and a storage medium, which can solve the problem of low calculation accuracy in the process of calculating the coverage rate of GNSS-R satellites in related technologies. Since the coverage rate is an important indicator for measuring the spatiotemporal resolution of GNSS-R satellites, the larger the coverage rate, the higher the spatiotemporal resolution of the GNSS-R satellites and the better the observation effect. When the calculation accuracy of the coverage rate is low, the spatiotemporal resolution of the GNSS-R satellites will be reduced.

[0005] In order to solve the above technical problems, the embodiments of the present application are implemented as follows:

[0006] In a first aspect, an embodiment of the present application provides a method for calculating GNSS-R satellite coverage, which is applied to a GNSS-R system. The method includes:

[0007] Obtain GNSS-R satellite data;

[0008] Determine a grid map of the earth's surface based on a preset unit distance, and determine the latitude and longitude range of each grid;

[0009] According to the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, the GNSS-R satellite data is mapped into a grid map on the earth's surface to obtain the number of revisits and target longitude and latitude data corresponding to each grid;

[0010] Determine the longitude change range according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data;

[0011] Determine a revisit number matrix according to the revisit number;

[0012] Based on the longitude variation range, the number of grids that are not 0 in any row in the revisit count matrix, the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, and the preset unit distance, the average coverage rate corresponding to each latitude interval is determined.

[0013] Optionally, the longitude variation range is expressed as:

[0014] ;

[0015] in, is the average distance from the Earth, is the average radius of the Earth, cos( ) is the longitude change function.

[0016] Optionally, the average coverage rate corresponding to each latitude interval is expressed as:

[0017] ;

[0018] in, is the preset unit distance, is the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, is the number of grids that are not 0 in any row of the revisit count matrix, It is the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data.

[0019] Optionally, the method further comprises:

[0020] Based on the longitude data and latitude data corresponding to each mirror reflection point, determine the coordinate data corresponding to each mirror reflection point in the revisit number matrix;

[0021] Based on the coordinate data and the single grid, the number of revisits corresponding to the single grid is determined.

[0022] Optionally, determining the number of revisits corresponding to a single grid based on the coordinate data and the single grid includes:

[0023] Determining the number of specular reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the specular reflection points;

[0024] When the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is greater than or equal to the preset time threshold, determining that the corresponding number of revisits in the single grid is increased by 1;

[0025] When the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is less than the preset time threshold, it is determined that the corresponding number of revisits in the single grid is increased by 0.

[0026] Optionally, the longitude data corresponding to a single grid map in the earth surface grid map and latitude data It is expressed as:

[0027] ;

[0028] ;

[0029] in, is the preset unit distance, , , , 360 / d+1.

[0030] In a second aspect, an embodiment of the present application provides a device for calculating GNSS-R satellite coverage, the device comprising:

[0031] Acquisition module, used to acquire GNSS-R satellite data;

[0032] A determination module, used to determine a grid map of the earth's surface based on a preset unit distance, and determine the latitude and longitude range of each grid;

[0033] An obtaining module is used to map the GNSS-R satellite data into a grid map on the earth's surface according to the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtain the number of revisits and target longitude and latitude data corresponding to each grid;

[0034] The determination module is also used to determine the longitude change range according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data; determine the revisit number matrix according to the revisit number; based on the longitude change range, the number of grids that are not 0 in any row in the revisit number matrix, the number of grids contained in the longitude change range corresponding to each latitude interval in the target longitude and longitude data and the preset unit distance, determine the average coverage rate corresponding to each latitude interval.

[0035] Optionally, the device further comprises:

[0036] The determination submodule is used to determine the coordinate data corresponding to each mirror reflection point in the revisit number matrix based on the longitude data and latitude data corresponding to each mirror reflection point; and determine the revisit number corresponding to a single grid based on the coordinate data and the single grid.

[0037] Optionally, the device further comprises:

[0038] The first determination submodule is used to determine the number of mirror reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection points; when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is greater than or equal to the preset time threshold, determine that the corresponding number of revisits in the single grid is increased by 1; when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is less than the preset time threshold, determine that the corresponding number of revisits in the single grid is increased by 0.

[0039] In a third aspect, an embodiment of the present application provides an electronic device, comprising a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other through the communication bus; the memory is used to store executable instructions, and the executable instructions enable the processor to execute the GNSS-R satellite coverage calculation method as described in any one of the above items.

[0040] In a fourth aspect, an embodiment of the present application provides a readable storage medium, on which a program or instruction is stored. When the program or instruction is executed by a processor, the method for calculating the GNSS-R satellite coverage as described in any one of the above items is implemented.

[0041] In an embodiment of the present application, a method for calculating GNSS-R satellite coverage, an electronic device, and a storage medium are provided. The embodiment of the present application obtains GNSS-R satellite data; determines a grid map of the earth's surface based on a preset unit distance, and determines the longitude and latitude range of each grid; maps the GNSS-R satellite data to the grid map of the earth's surface according to the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtains the number of revisits and target longitude and latitude data corresponding to each grid; determines the longitude variation range according to the average distance of the earth, the average radius of the earth, and the longitude variation function corresponding to each latitude interval in the target longitude and longitude data; determines the revisit number matrix according to the revisit number; and calculates the revisit number matrix based on the longitude variation range, the number of grids that are not 0 in any row in the revisit number matrix, the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and longitude data, and the preset unit distance. , determine the average coverage rate corresponding to each latitude interval. In the process of calculating the GNSS-R satellite coverage rate, the GNSS-R system can first obtain the GNSS-R satellite data, and then map the GNSS-R satellite data to the earth's surface grid map to obtain the number of revisits and the target longitude and latitude data corresponding to each grid. In the subsequent calculation process, the number of revisits and the target longitude and latitude data corresponding to each grid can be used for calculation, and then the longitude change range is determined according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data; finally, based on the longitude change range, the number of grids that are not 0 in any row in the revisit number matrix, the number of grids contained in the longitude change range corresponding to each latitude interval in the target longitude and longitude data, and the preset unit distance, the average coverage rate corresponding to each latitude interval is determined, and then the detection capability of the GNSS-R satellite can be evaluated.

[0042] The above description is only an overview of the technical solution of the present application. In order to more clearly understand the technical means of the present application, it can be implemented in accordance with the contents of the specification. In order to make the above and other purposes, features and advantages of the present application more obvious and easy to understand, the specific implementation methods of the present application are listed below. BRIEF DESCRIPTION OF THE DRAWINGS

[0043] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings required for use in the description of the embodiments of the present application will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0044] Figure 1 A step flow of a method for calculating GNSS-R satellite coverage provided in an embodiment of the present application Figure 1 ;

[0045] Figure 2 A step flow of a method for calculating GNSS-R satellite coverage provided in an embodiment of the present application Figure 2 ;

[0046] Figure 3 A schematic diagram of the position and grid of a mirror reflection point provided in an embodiment of the present application;

[0047] Figure 4 The coverage rate of the GNSS-R satellites provided in the embodiment of the present application at different latitudes in different time periods under a 0.25 degree grid;

[0048] Figure 5 A logic block diagram of a GNSS-R satellite coverage calculation device provided in an embodiment of the present application;

[0049] Figure 6 A schematic diagram of the structure of an electronic device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0050] The following will be combined with the drawings in the embodiments of the present application to clearly and completely describe the technical solutions in the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, not all of the embodiments. Based on the embodiments in the present application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0051] The terms "first", "second", etc. in the specification and claims of the present application are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable under appropriate circumstances, so that the embodiments of the present application can be implemented in an order other than those illustrated or described here, and the objects distinguished by "first", "second", etc. are generally of one type, and the number of objects is not limited. For example, the first object can be one or more. In addition, "and / or" in the specification and claims represents at least one of the connected objects, and the character " / " generally indicates that the objects associated with each other are in an "or" relationship.

[0052] Method Embodiment

[0053] The following is a detailed description of the method for calculating the GNSS-R satellite coverage provided in the embodiment of the present application through specific embodiments and application scenarios in conjunction with the accompanying drawings.

[0054] Reference Figure 1 , showing the steps of a method for calculating GNSS-R satellite coverage provided in an embodiment of the present application Figure 1 ,like Figure 1As shown, the method specifically includes steps S101 to S105:

[0055] Step S101: Acquire GNSS-R satellite data.

[0056] In this step, during the process of calculating the GNSS-R satellite coverage, the GNSS-R system may first acquire GNSS-R satellite data.

[0057] In this step, the longitude data corresponding to the GNSS-R satellite data may be any longitude data corresponding to the GNSS-R satellite data composed of 22 polar-orbiting satellites (orbit altitude 520 km), which is not limited in the embodiments of the present application.

[0058] In this step, the grid map of the earth's surface can be determined according to a preset unit distance, and the preset unit distance can be 0.25°, 0.5°, 1°, etc., which is not limited in the embodiments of the present application.

[0059] In this step, the GNSS-R satellite data may include GNSS-R satellite data composed of 22 polar-orbiting satellites (orbit altitude 520 km), which is not limited in the embodiments of the present application.

[0060] Step S102: Determine a grid map of the earth's surface based on a preset unit distance, and determine the longitude and latitude range of each grid.

[0061] In this step, during the calculation method of the GNSS-R satellite coverage, after acquiring the GNSS-R satellite data, the GNSS-R system can determine a grid map of the earth's surface based on a preset unit distance, and determine the longitude and latitude range of each grid.

[0062] In this step, a single grid can be determined based on a preset unit distance. The single grid can be an average distance of about 100d=25km on the earth when the preset unit distance is 0.25°, and then the single grid can be 25km×25km, which is not limited in the embodiments of the present application.

[0063] Furthermore, the longitude data corresponding to a single grid map in the earth surface grid map and latitude data It is expressed as:

[0064] ;

[0065] ;

[0066] in, is the preset unit distance, , , , 360 / d+1.

[0067] In this step, a single grid map may be the smallest unit in the grid map of the earth's surface, which is not limited in the embodiments of the present application.

[0068] In this step, the longitude data corresponding to a single grid map can be obtained with the center of the earth as the origin, which is not limited in the embodiments of the present application.

[0069] Exemplarily, in this step, taking the center of the earth as the origin, the longitude data corresponding to a single grid map can be 116 degrees 20 minutes east longitude, or 45°, etc., which is not limited in the embodiments of the present application.

[0070] In this step, the latitude data corresponding to a single grid map can be obtained with the center of the earth as the origin, which is not limited in the embodiments of the present application.

[0071] Exemplarily, in this step, taking the center of the earth as the origin, the latitude data corresponding to a single grid map can be obtained as 39 degrees 54 minutes north latitude, or 45°, etc., which is not limited in the embodiments of the present application.

[0072] Step S103: Map the GNSS-R satellite data to a grid map of the earth's surface according to the latitude and longitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtain the number of revisits and target longitude and latitude data corresponding to each grid.

[0073] After determining the earth's surface grid map based on the preset unit distance and determining the longitude and latitude range of each grid, the GNSS-R system can map the GNSS-R satellite data to the earth's surface grid map according to the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtain the number of revisits and target longitude and latitude data corresponding to each grid.

[0074] In this step, the longitude data corresponding to each mirror reflection point can be 100 degrees 20 minutes east longitude, or 45°, etc., which is not limited in the embodiments of the present application.

[0075] Exemplarily, in this step, the latitude data corresponding to each mirror reflection point may be 40 degrees and 50 minutes north latitude, or 35°, etc., which is not limited in the embodiments of the present application.

[0076] In this step, the process of mapping the GNSS-R satellite data into the earth surface grid map may include mapping the GNSS-R satellite observation data for a period of time into the earth surface grid map, thereby obtaining the mirror reflection point position data.

[0077] Among them, the number of revisits corresponding to each grid can be any integer greater than or equal to 0, for example: the number of revisits corresponding to each grid can be 8 times, or the number of revisits corresponding to each grid can be 2 times, etc., which is not limited in the embodiments of the present application.

[0078] It should be noted that the number of revisits refers to the number of revisits to the same grid.

[0079] The target longitude and latitude data may include target longitude data and target latitude data, and the target longitude and latitude data are both longitude and latitude data obtained after mapping.

[0080] Step S104: Determine the longitude change range according to the average earth distance, the average earth radius and the longitude change function corresponding to each latitude interval in the target longitude and latitude data.

[0081] After mapping the GNSS-R satellite data to the earth's surface grid map based on the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection points in the GNSS-R satellite data, and obtaining the number of revisits and target longitude and latitude data corresponding to each grid, the GNSS-R system can determine the longitude change range based on the average distance to the earth, the average radius of the earth, and the longitude change function corresponding to each latitude interval in the target longitude and latitude data.

[0082] Furthermore, the longitude range is expressed as:

[0083] ;

[0084] in, is the average distance from the Earth, is the average radius of the Earth, cos( ) is the longitude change function.

[0085] Exemplarily, the longitude change function can characterize how the longitude changes with distance and moving direction when moving a certain distance on the surface of the earth, which is not limited in the embodiments of the present application.

[0086] Step S105: Determine a revisit frequency matrix according to the revisit frequency.

[0087] After determining the longitude variation range based on the average distance to the earth, the average radius of the earth, and the longitude variation function corresponding to each latitude interval in the target longitude and latitude data, the GNSS-R system can determine the revisit number matrix based on the revisit number.

[0088] In this step, after comparing the time threshold between two adjacent mirror reflection points with the preset time threshold, the GNSS-R system may determine a revisit number matrix based on the corresponding revisit number in a single grid.

[0089] In this step, after determining the revisit number matrix based on the corresponding revisit number in a single grid, the GNSS-R system may also determine the revisit time matrix based on the revisit number matrix and the time parameter.

[0090] Furthermore, the revisit time matrix is ​​expressed as:

[0091] ;

[0092] in, is the time parameter, is the revisit count matrix, , , , 360 / d+1.

[0093] Exemplarily, in this step, Can be Etc., are not limited in the embodiments of the present application.

[0094] Exemplarily, in this step, It can be 361, etc., which is not limited in the embodiments of the present application.

[0095] In this step, after determining the revisit number matrix based on the corresponding revisit number in a single grid, the GNSS-R system may also determine the revisit time matrix based on the revisit number matrix and the time parameter.

[0096] Furthermore, the revisit time matrix is ​​expressed as:

[0097] ;

[0098] in, is the time parameter, is the revisit count matrix, , , , 360 / d+1.

[0099] Exemplarily, in this step, Can be Etc., are not limited in the embodiments of the present application.

[0100] Step S106, based on the longitude variation range, the number of grids that are not 0 in any row in the revisit times matrix, the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, and the preset unit distance, determine the average coverage rate corresponding to each latitude interval.

[0101] After determining the revisit count matrix according to the number of revisits, the GNSS-R system can determine the average coverage corresponding to each latitude interval based on the longitude change range, the number of grids that are not 0 in any row in the revisit count matrix, the number of grids contained in the longitude change range corresponding to each latitude interval in the target longitude and latitude data, and the preset unit distance.

[0102] In this step, the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data may be 361, etc., which is not limited in the embodiments of the present application.

[0103] In this step, the longitude variation range can be determined based on the average distance of the earth, the average radius of the earth and the longitude variation function.

[0104] In this step, the preset unit distance may be 0.25°, 0.5°, 1°, etc., which is not limited in the embodiments of the present application.

[0105] In this step, in the process of determining the corresponding number of revisits in a single grid, the GNSS-R system may also determine the corresponding number of revisits in a single grid based on the number of continuous sampling points falling within the single grid of the target longitude and latitude data.

[0106] In this step, the number of continuous sampling points falling within a single grid of the target longitude and latitude data may be any value greater than 0, for example, the number of continuous sampling points may be 100, or 300, etc., which is not limited in the embodiments of the present application.

[0107] Illustratively, in this step, the corresponding number of revisits in a single grid may be 5, which is not limited in the embodiments of the present application.

[0108] In this step, when determining the corresponding number of revisits in a single grid based on the number of continuous sampling points falling within a single grid of target longitude and latitude data, the GNSS-R system may first obtain the longitude and latitude data corresponding to each mirror reflection point.

[0109] In this step, the longitude data corresponding to each mirror reflection point can be 100 degrees 20 minutes east longitude, or 45°, etc., which is not limited in the embodiments of the present application.

[0110] Exemplarily, in this step, the latitude data corresponding to each mirror reflection point may be 40 degrees and 50 minutes north latitude, or 35°, etc., which is not limited in the embodiments of the present application.

[0111] In this step, after obtaining the longitude data and latitude data corresponding to each mirror reflection point, the GNSS-R system can also determine the coordinate data corresponding to each mirror reflection point in the revisit count matrix based on the longitude data and latitude data.

[0112] In this step, the coordinate data corresponding to each mirror reflection point in the revisit frequency matrix may include abscissa data and ordinate data, the abscissa data corresponds to latitude data, and the ordinate data corresponds to longitude data.

[0113] Exemplarily, the coordinate data may be the coordinate data of the 175th row and the 306th column in the revisit frequency matrix, etc., which is not limited in the embodiments of the present application.

[0114] In this step, after determining the coordinate data corresponding to each mirror reflection point in the revisit number matrix based on the longitude data and the latitude data, the GNSS-R system can also determine the revisit number corresponding to a single grid based on the coordinate data and the single grid.

[0115] In this step, the time parameter may be any data greater than 0, for example, the time parameter may be 10 days, or 5 days, etc., which is not limited in the embodiments of the present application.

[0116] In this step, the revisit time matrix may be determined based on the revisit number matrix and the time parameter.

[0117] Exemplarily, the revisit time matrix may be a time matrix with 181 rows and 361 columns, and the dimensions of the revisit time matrix and the revisit number matrix are the same.

[0118] In this step, when determining the revisit time matrix based on the corresponding revisit times and time parameters in a single grid, the GNSS-R system may also determine the revisit times matrix based on the corresponding revisit times in a single grid.

[0119] In this step, the revisit frequency matrix may be a matrix of 181 rows and 361 columns, and the elements in the matrix are determined by comparing the time threshold between two adjacent mirror reflection points with a preset time threshold.

[0120] In this step, in the process of determining the number of grids that are not zero in any row of the revisit count matrix, the GNSS-R system can first determine the number of mirror reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection point; and then determine the corresponding number of revisits in a single grid based on the time threshold between two adjacent mirror reflection points and the size of the preset time threshold, and then determine the number of grids that are not zero in any row.

[0121] In this step, when determining the revisit number matrix based on the corresponding revisit number in a single grid, the GNSS-R system can also determine the number of mirror reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection point.

[0122] In this step, the diagonal length corresponding to a single grid can be determined according to the length and width of the single grid. For example, when the length and width of a single grid are both 25 km, the diagonal length corresponding to the single grid can be determined to be 25 km. , which is not limited in the embodiments of the present application.

[0123] Exemplarily, in this step, the moving speed of the mirror reflection point may be 6 km / s, which is not limited in the embodiments of the present application.

[0124] In this step, the number of specular reflection points in a single grid can be determined according to the diagonal length corresponding to the single grid and the moving speed of the specular reflection points. For example, the diagonal length corresponding to the single grid is divided by the moving speed of the specular reflection points to obtain the number of specular reflection points in the single grid. For example, when the diagonal length is approximately , when the moving speed of the mirror reflection point is 6km / s, the number of mirror reflection points in a single grid can be expressed as: This is not limited in the embodiments of the present application.

[0125] In this step, the time threshold between two adjacent mirror reflection points may be 5 seconds, or 7 seconds, etc., which is not limited in the embodiments of the present application.

[0126] In this step, the preset time threshold may be equivalent to the number of mirror reflection points in a single grid. For example, when the number of mirror reflection points in a single grid is 5, the preset time threshold may be 5s, etc., which is not limited in the embodiments of the present application.

[0127] In this step, when the time threshold between two adjacent mirror reflection points is greater than or equal to the preset time threshold, the corresponding number of revisits in a single grid is increased by 1. For example, when the time threshold between two adjacent mirror reflection points is 6s and the preset time threshold is 5s, the corresponding number of revisits in a single grid is increased by 1, which is not limited in the embodiments of the present application.

[0128] In this step, when the time threshold between two adjacent mirror reflection points is less than the preset time threshold, the corresponding number of revisits in a single grid is increased by 0. For example, when the time threshold between two adjacent mirror reflection points is 5s and the preset time threshold is 6s, the corresponding number of revisits in a single grid is increased by 0, which is not limited in the embodiments of the present application.

[0129] Furthermore, the average coverage rate corresponding to each latitude interval is expressed as:

[0130] ;

[0131] in, is the preset unit distance, is the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, is the number of grids that are not 0 in any row of the revisit count matrix, It is the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data.

[0132] In this step, The number of grids that are not 0 in any row of the revisit number matrix can be any number greater than 0. For example, the number of grids that are not 0 in any row of the revisit number matrix can be 50, or the number of grids that are not 0 in any row of the revisit number matrix can be 80, etc., which is not limited in the embodiments of the present application.

[0133] Furthermore, in the process of performing the method for calculating the GNSS-R satellite coverage rate, the method for calculating the GNSS-R satellite coverage rate may further include steps S107 to S108:

[0134] Step S107: Determine the coordinate data corresponding to each mirror reflection point in the revisit times matrix based on the longitude data and latitude data corresponding to each mirror reflection point.

[0135] In this step, during the calculation method of the GNSS-R satellite coverage, the GNSS-R system can also determine the coordinate data corresponding to each mirror reflection point in the revisit number matrix based on the longitude data and latitude data corresponding to each mirror reflection point.

[0136] In this step, the coordinate data corresponding to each mirror reflection point in the revisit frequency matrix may include abscissa data and ordinate data, the abscissa data corresponds to latitude data, and the ordinate data corresponds to longitude data.

[0137] Exemplarily, the coordinate data may be the coordinate data of the 175th row and the 306th column in the revisit frequency matrix, etc., which is not limited in the embodiments of the present application.

[0138] Step S108: Determine the number of revisits corresponding to the single grid based on the coordinate data and the single grid.

[0139] In this step, after determining the coordinate data corresponding to each mirror reflection point in the revisit number matrix based on the longitude data and latitude data corresponding to each mirror reflection point, the GNSS-R system can also determine the revisit number corresponding to a single grid based on the coordinate data and the single grid.

[0140] Furthermore, in the process of determining the number of revisits corresponding to a single grid based on the coordinate data and the single grid, the method for GNSS-R satellite coverage may further include steps S1081 to S1083:

[0141] Step S1081: Determine the number of specular reflection points in a single grid based on the length of the diagonal corresponding to the single grid and the moving speed of the specular reflection points.

[0142] In this step, in the process of determining the number of revisits corresponding to a single grid based on the coordinate data and the single grid, the GNSS-R system can also determine the number of mirror reflection points in the single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection points.

[0143] In this step, the diagonal length corresponding to a single grid can be determined according to the length and width of the single grid. For example, when the length and width of a single grid are both 25 km, the diagonal length corresponding to the single grid can be determined to be 25 km. , which is not limited in the embodiments of the present application.

[0144] Exemplarily, in this step, the moving speed of the mirror reflection point may be 6 km / s, which is not limited in the embodiments of the present application.

[0145] In this step, the number of specular reflection points in a single grid can be determined according to the diagonal length corresponding to the single grid and the moving speed of the specular reflection points. For example, the diagonal length corresponding to the single grid is divided by the moving speed of the specular reflection points to obtain the number of specular reflection points in the single grid. For example, when the diagonal length is approximately , when the moving speed of the mirror reflection point is 6 km / s, the number of mirror reflection points in a single grid can be obtained, which is not limited in the embodiments of the present application.

[0146] Step S1082: when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is greater than or equal to the preset time threshold, determine the corresponding number of revisits in a single grid and add 1.

[0147] In this step, after determining the number of mirror reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection point, the GNSS-R system can also determine the corresponding number of revisits in the single grid plus 1 when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is greater than or equal to the preset time threshold.

[0148] In this step, the time threshold between two adjacent mirror reflection points may be 5 seconds, or 7 seconds, etc., which is not limited in the embodiments of the present application.

[0149] In this step, the preset time threshold may be equivalent to the number of mirror reflection points in a single grid. For example, when the number of mirror reflection points in a single grid is 5, the preset time threshold may be 5s, etc., which is not limited in the embodiments of the present application.

[0150] In this step, when the time threshold between two adjacent mirror reflection points is greater than or equal to the preset time threshold, the corresponding number of revisits in a single grid is increased by 1. For example, when the time threshold between two adjacent mirror reflection points is 6s and the preset time threshold is 5s, the corresponding number of revisits in a single grid is increased by 1, which is not limited in the embodiments of the present application.

[0151] Step S1083: when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is less than the preset time threshold, determine the corresponding number of revisits in a single grid and add 0.

[0152] In this step, after determining the number of mirror reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection point, the GNSS-R system can also determine that the corresponding number of revisits in the single grid is plus 0 when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is less than the preset time threshold.

[0153] In this step, when the time threshold between two adjacent mirror reflection points is less than the preset time threshold, the corresponding number of revisits in a single grid is increased by 0. For example, when the time threshold between two adjacent mirror reflection points is 5s and the preset time threshold is 6s, the corresponding number of revisits in a single grid is increased by 0, which is not limited in the embodiments of the present application.

[0154] In summary, an embodiment of the present application provides a method for calculating GNSS-R satellite coverage, an electronic device, and a storage medium. The embodiment of the present application determines a grid map of the earth's surface based on a preset unit distance, and determines the longitude and latitude range of each grid; obtains GNSS-R satellite data; maps the GNSS-R satellite data to the grid map of the earth's surface according to the longitude and latitude range and the longitude data and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtains data corresponding to each grid, wherein the data corresponding to each grid includes the number of revisits to each grid and the target longitude and latitude data corresponding to each grid; determines the longitude change range according to the average distance of the earth, the average radius of the earth, and the longitude change function corresponding to each latitude interval in the target longitude and longitude data; determines the longitude change range based on the longitude change range, the number of grids that are not 0 in any row in the revisit number matrix, the number of grids contained in the longitude change range corresponding to each latitude interval in the target longitude and longitude data, and the preset unit distance. , determine the average coverage rate corresponding to each latitude interval. In the process of calculating the GNSS-R satellite coverage rate, the GNSS-R system can first obtain the GNSS-R satellite data, and then map the GNSS-R satellite data to the earth surface grid map to obtain the data corresponding to each grid. The data corresponding to each grid includes the number of revisits to each grid and the target longitude and latitude data corresponding to each grid. In the subsequent calculation process, the target longitude and latitude data are used for calculation, and then the longitude change range is determined according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data; finally, based on the longitude change range, the number of grids that are not 0 in any row in the revisit number matrix, the number of grids contained in the longitude change range corresponding to each latitude interval in the target longitude and longitude data, and the preset unit distance, the average coverage rate corresponding to each latitude interval is determined, and then the detection capability of the GNSS-R satellite can be evaluated.

[0155] GNSS-R technology uses existing GNSS-R satellites and does not require active signal reflection. It has the advantages of miniaturization, low cost and low power consumption. It is very suitable for being carried by small satellites to form constellations and achieve high temporal and spatial resolution observation of the earth. The main difference between GNSS-R technology and traditional remote sensing technology is that the observation points of GNSS-R are strip-shaped, and the sampling points in the strip are dense. The analysis of the temporal and spatial resolution of GNSS-R satellites is very important for the constellation design and detection capability evaluation of GNSS-R. Coverage is an important indicator to measure the temporal and spatial resolution of satellite constellations. It is defined as the ratio of the coverage area of ​​the earth's surface by the constellation observation data within a certain time and at a certain spatial resolution to the total surface area of ​​the earth. The larger the coverage, the higher the temporal and spatial resolution of the constellation and the better the observation effect. Due to the unique observation system of GNSS-R technology, the calculation of its coverage also needs to be specifically considered based on its technical characteristics.

[0156] In order to calculate the coverage of GNSS-R satellites, this solution will first divide the earth into grids of a certain size, and then use the mirror reflection point positions in the GNSS-R satellite observation data over a period of time to calculate the number of revisits within each grid point. Since GNSS-R sampling is dense, in order to better reflect the temporal and spatial resolution of the constellation, it is necessary to eliminate continuous observation data within the same grid.

[0157] After obtaining the number of revisits in each grid, the coverage rate needs to be calculated at different latitudes because the coverage of the earth by GNSS-R satellites varies greatly with latitude. The coverage rate is calculated by dividing the number of grids covered by GNSS-R data at each latitude during a certain observation time by the total number of grids corresponding to that dimension.

[0158] Using GNSS-R satellite data composed of 22 polar-orbiting satellites (orbit altitude 520km), calculate the temporal and spatial resolution under the grid size d, d is generally 0.25°, 0.5°, 1°. In this case, d is 0.25 degrees (the average distance of 0.25° on the earth is about 100d=25km).

[0159] Reference Figure 2 , showing the steps of a method for calculating GNSS-R satellite coverage provided in an embodiment of the present application Figure 2 ,like Figure 2 As shown, the method specifically includes steps S201 to S206:

[0160] Step S201: determine the grid size based on a preset unit distance.

[0161] In the process of calculating the coverage, the GNSS-R system can first determine the grid size based on a preset unit distance.

[0162] Step S202: Initialize the global grid.

[0163] After determining the grid size based on a preset unit distance, the GNSS-R system can initialize the global grid.

[0164] Step S203: Obtain the longitude and latitude of the mirror reflection point corresponding to the GNSS-R satellite.

[0165] During the coverage calculation process, the GNSS-R system can also obtain the longitude and latitude of the GNSS-R mirror reflection point.

[0166] Step S204: Calculate the revisit frequency matrix based on the longitude and latitude of the mirror reflection point corresponding to the GNSS-R satellite.

[0167] After obtaining the longitude and latitude of the GNSS-R mirror reflection point, the GNSS-R system can also calculate a revisit number matrix based on the longitude and latitude of the GNSS-R mirror reflection point.

[0168] Step S205: Calculate the revisit time matrix based on the revisit frequency matrix.

[0169] After calculating the revisit number matrix based on the longitude and latitude of the GNSS-R mirror reflection point, the GNSS-R system can also calculate the revisit time matrix based on the revisit number matrix.

[0170] Step S206: Calculate coverage at different latitudes based on the revisit time matrix.

[0171] After calculating the revisit time matrix based on the revisit number matrix, the GNSS-R system can also calculate the coverage at different latitudes based on the revisit time matrix.

[0172] Specifically, the process of calculating the number of revisits of a GNSS-R satellite within k hours is as follows:

[0173] Grid the Earth The latitude of each grid can be expressed as:

[0174] (1)

[0175] in, is the preset unit distance, , .

[0176] The longitude corresponding to each grid can be expressed as:

[0177] (2)

[0178] in, is the preset unit distance, , 360 / d+1.

[0179] Initialize the revisit count matrix , the dimension is , The elements in are initialized to 0.

[0180] For each mirror reflection point in the k-hour GNSS-R data, its latitude and longitude are extracted. The subscript corresponding to the mirror reflection point in the revisit number matrix can be expressed as:

[0181] (3)

[0182] (4)

[0183] The {} symbol indicates rounding.

[0184] In order to avoid repeated calculation of revisit times due to consecutive sampling points of GNSS-R data falling in the same grid, it is necessary to consider when calculating the revisit times matrix that the moving speed of the GNSS-R mirror reflection point on the earth's surface is about 6 km / s, and the average diagonal length of a grid is about d, so there are at most Continuous observation points. Traverse all GNSS-R data and calculate the subscript in the revisit count matrix C for each mirror reflection point , ,when When is 0, or the difference in observation time between two observation data in the same grid is greater than m, the corresponding number of revisits is increased by one, which can be expressed as:

[0185] (5)

[0186] Reference Figure 3 , shows a schematic diagram of the position and grid of a mirror reflection point provided in an embodiment of the present application, such as Figure 3 shown.

[0187] Specifically, the process of calculating coverage at different latitudes is as follows:

[0188] Traverse the latitude range , calculate the average coverage rate within k hours at each latitude. For the latitude, the corresponding longitude range within a distance of 100d can be expressed as the following formula (6):

[0189]

[0190] in, is the average distance from the Earth, is the average radius of the Earth, cos( ) is the longitude change function.

[0191] The average coverage at latitude can be expressed as follows:

[0192]

[0193] in, is the preset unit distance, is the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, is the number of grids that are not 0 in any row of the revisit count matrix, It is the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data.

[0194] Reference Figure 4 , shows the coverage of the GNSS-R satellite provided in the embodiment of the present application at different latitudes in different time periods under a 0.25 degree grid, such as Figure 4 shown.

[0195] In summary, when calculating the number of revisits, it is necessary to consider the GNSS-R data sampling interval and eliminate the continuous sampling in the same grid to reflect the actual situation; the accurate calculation method of coverage and analysis at different latitudes.

[0196] Device Embodiment

[0197] like Figure 5 As shown, Figure 5 A logic block diagram of a GNSS-R satellite coverage calculation device provided in an embodiment of the present application is shown, and the device includes:

[0198] An acquisition module 1001 is used to acquire GNSS-R satellite data;

[0199] A determination module 1002 is used to determine a grid map of the earth's surface based on a preset unit distance, and determine the latitude and longitude range of each grid;

[0200] The obtaining module 1003 is used to map the GNSS-R satellite data into a grid map of the earth's surface according to the latitude and longitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtain the number of revisits and target longitude and latitude data corresponding to each grid;

[0201] The determination module is also used to determine the longitude change range according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data; determine the revisit number matrix according to the revisit number; based on the longitude change range, the number of grids that are not 0 in any row in the revisit number matrix, the number of grids contained in the longitude change range corresponding to each latitude interval in the target longitude and longitude data and the preset unit distance, determine the average coverage rate corresponding to each latitude interval.

[0202] Optionally, the device further comprises:

[0203] The determination submodule is used to determine the coordinate data corresponding to each mirror reflection point in the revisit number matrix based on the longitude data and latitude data corresponding to each mirror reflection point; and determine the revisit number corresponding to a single grid based on the coordinate data and the single grid.

[0204] Optionally, the device further comprises:

[0205] The first determination submodule is used to determine the number of mirror reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the mirror reflection points; when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is greater than or equal to the preset time threshold, determine that the corresponding number of revisits in the single grid is increased by 1; when the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is less than the preset time threshold, determine that the corresponding number of revisits in the single grid is increased by 0.

[0206] Optionally, the longitude data corresponding to a single grid map in the earth surface grid map and latitude data It is expressed as:

[0207] ;

[0208] ;

[0209] in, is the preset unit distance, , , , 360 / d+1.

[0210] Optionally, the average coverage rate corresponding to each latitude interval is expressed as:

[0211] ;

[0212] in, is the preset unit distance, is the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, is the number of grids that are not 0 in any row of the revisit count matrix, It is the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data.

[0213] Optionally, the longitude variation range is expressed as:

[0214] ;

[0215] in, is the average distance from the Earth, is the average radius of the Earth, cos( ) is the longitude change function.

[0216] As for the device embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiment.

[0217] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0218] The calculation device of the GNSS-R satellite coverage in the embodiment of the present application can be an electronic device, or a component in the electronic device, such as a calculation circuit or chip of the GNSS-R satellite coverage. The electronic device can be a terminal, or it can be other devices other than the terminal. Exemplarily, the electronic device can be a GPU BOX, a mobile phone, a tablet computer, a laptop computer, a PDA, a vehicle-mounted electronic device, a mobile Internet device (Mobile Internet Device, MID), an augmented reality (augmented reality, AR) / virtual reality (virtual reality, VR) device, a robot, a wearable device, an ultra-mobile personal computer (ultra-mobile personal computer, UMPC), a netbook or a personal digital assistant (personal digital assistant, PDA), etc. It can also be a server, a network attached storage (Network Attached Storage, NAS), a personal computer (personal computer, PC), a television (television, TV), a teller machine or a self-service machine, etc., which is not specifically limited in the embodiment of the present application.

[0219] The GNSS-R satellite coverage calculation device in the embodiment of the present application may be a device having an operating system. The operating system may be an Android operating system, a Linux operating system, a Windows operating system, or other possible operating systems, which are not specifically limited in the embodiment of the present application.

[0220] The GNSS-R satellite coverage calculation device provided in the embodiment of the present application can achieve Figure 1 To avoid repetition, the various processes implemented by the method embodiment are not described here.

[0221] Alternatively, if Figure 6 As shown, an embodiment of the present application further provides an electronic device 1100, including a processor 1101 and a memory 1102, wherein the memory 1102 stores a program or instruction that can be executed on the processor 1101, and when the program or instruction is executed by the processor 1101, the various steps of the above-mentioned GNSS-R satellite coverage calculation method embodiment are implemented, and the same technical effect can be achieved. To avoid repetition, it will not be repeated here.

[0222] In an embodiment of the present application, the memory 1102 may be used to store software programs and various data. The memory 1102 may mainly include a first storage area for storing programs or instructions and a second storage area for storing data, wherein the first storage area may store an operating system, an application program or instructions required for at least one function (such as a sound playback function, an image playback function, etc.), etc. In addition, the memory 1102 may include a volatile memory or a non-volatile memory, or the memory 1102 may include both volatile and non-volatile memories. Among them, the non-volatile memory may be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory may be a random access memory (RAM), a static random access memory (SRAM), a dynamic random access memory (DRAM), a synchronous dynamic random access memory (SDRAM), a double data rate synchronous dynamic random access memory (DDRSDRAM), an enhanced synchronous dynamic random access memory (ESDRAM), a synchronous link dynamic random access memory (SLDRAM) and a direct memory bus random access memory (DRRAM). The memory 1102 in the embodiment of the present application includes but is not limited to these and any other suitable types of memory.

[0223] The processor 1101 may include one or more processing units; optionally, the processor 1101 includes an application processor and a modem processor, wherein the application processor mainly processes operations involving an operating system, a user interface, and application programs, and the modem processor mainly processes wireless communication signals, such as a baseband processor. It is understandable that the above-mentioned modem processor may also be incorporated into the processor 1101.

[0224] An embodiment of the present application also provides an electronic device, which includes a processor, a memory, a communication interface and a communication bus, wherein the processor, the memory and the communication interface communicate with each other through the communication bus; the memory is used to store executable instructions, and the executable instructions enable the processor to execute the various processes of the embodiment of the GNSS-R satellite coverage calculation method as described above, and can achieve the same technical effect. To avoid repetition, it will not be repeated here.

[0225] An embodiment of the present application also provides a readable storage medium, on which a program or instruction is stored. When the program or instruction is executed by a processor, each process of the above-mentioned GNSS-R satellite coverage calculation method embodiment is implemented, and the same technical effect can be achieved. To avoid repetition, it will not be repeated here.

[0226] The processor is the processor in the electronic device described in the above embodiment. The readable storage medium includes a computer readable storage medium, such as a computer read-only memory ROM, a random access memory RAM, a magnetic disk or an optical disk.

[0227] An embodiment of the present application further provides a chip, which includes a processor and a communication interface, wherein the communication interface is coupled to the processor, and the processor is used to run programs or instructions to implement the various processes of the above-mentioned GNSS-R satellite coverage calculation method embodiment, and can achieve the same technical effect. To avoid repetition, it will not be repeated here.

[0228] It should be understood that the chip involved 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.

[0229] An embodiment of the present application provides a computer program product, which is stored in a storage medium. The program product is executed by at least one processor to implement the various processes of the above-mentioned GNSS-R satellite coverage calculation method embodiment, and can achieve the same technical effect. To avoid repetition, it will not be repeated here.

[0230] It should be noted that, in this article, the terms "comprise", "include" or any other variant thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of further restrictions, an element defined by the sentence "comprises one..." does not exclude the presence of other identical elements in the process, method, article or device including the element. In addition, it should be pointed out that the scope of the method and device in the embodiment of the present application is not limited to performing functions in the order shown or discussed, and may also include performing functions in a substantially simultaneous manner or in reverse order according to the functions involved, for example, the described method may be performed in an order different from that described, and various steps may also be added, omitted or combined. In addition, the features described with reference to certain examples may be combined in other examples.

[0231] Through the description of the above implementation methods, those skilled in the art can clearly understand that the above-mentioned embodiment methods can be implemented by means of software plus a necessary general hardware platform, and of course by hardware, but in many cases the former is a better implementation method. Based on such an understanding, the technical solution of the present application, or the part that contributes to the relevant technology, can be embodied in the form of a computer software product, which is stored in a storage medium (such as ROM / RAM, disk, CD), and includes a number of instructions for a terminal (which can be a mobile phone, computer, server, or network device, etc.) to execute the methods described in each embodiment of the present application.

[0232] The embodiments of the present application are described above in conjunction with the accompanying drawings, but the present application is not limited to the above-mentioned specific implementation methods. The above-mentioned specific implementation methods are merely illustrative and not restrictive. Under the guidance of the present application, ordinary technicians in this field can also make many forms without departing from the purpose of the present application and the scope of protection of the claims, all of which are within the protection of the present application.

Claims

1. A method for calculating GNSS-R satellite coverage, characterized in that: Applied to a GNSS-R system, the method comprises: Obtain GNSS-R satellite data; Determine a grid map of the earth's surface based on a preset unit distance, and determine the latitude and longitude range of each grid; According to the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, the GNSS-R satellite data is mapped into a grid map on the earth's surface to obtain the number of revisits and target longitude and latitude data corresponding to each grid; Determine the longitude change range according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data; Determine a revisit number matrix according to the revisit number; Determine the average coverage rate corresponding to each latitude interval based on the longitude variation range, the number of grids that are not 0 in any row of the revisit count matrix, the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, and the preset unit distance; Wherein, the longitude variation range is expressed as: ； in, is the average distance from the Earth, is the average radius of the Earth, cos( ) is the longitude change function; The average coverage rate corresponding to each latitude interval is expressed as: ； in, is the preset unit distance, is the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, is the number of grids that are not 0 in any row of the revisit count matrix, It is the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data.

2. The method according to claim 1, characterized in that: The method further comprises: Based on the longitude data and latitude data corresponding to each mirror reflection point, determine the coordinate data corresponding to each mirror reflection point in the revisit number matrix; Based on the coordinate data and the single grid, the number of revisits corresponding to the single grid is determined.

3. The method according to claim 2, characterized in that The step of determining the number of revisits corresponding to a single grid based on the coordinate data and the single grid includes: Determining the number of specular reflection points in a single grid based on the diagonal length corresponding to the single grid and the moving speed of the specular reflection points; When the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is greater than or equal to the preset time threshold, determining that the corresponding number of revisits in the single grid is increased by 1; When the time threshold between the coordinate data corresponding to two adjacent mirror reflection points is less than the preset time threshold, it is determined that the corresponding number of revisits in the single grid is increased by 0.

4. The method according to claim 1, characterized in that Longitude data corresponding to a single grid map in the earth surface grid map and latitude data It is expressed as: ； ； in, is the preset unit distance, , , , 360 / d+1.

5. A GNSS-R satellite coverage calculation device, characterized in that: The device comprises: Acquisition module, used to acquire GNSS-R satellite data; A determination module, used to determine a grid map of the earth's surface based on a preset unit distance, and determine the latitude and longitude range of each grid; An obtaining module is used to map the GNSS-R satellite data into a grid map on the earth's surface according to the longitude and latitude range and the longitude and latitude data corresponding to the mirror reflection point in the GNSS-R satellite data, and obtain the number of revisits and target longitude and latitude data corresponding to each grid; The determination module is further used to determine the longitude change range according to the average distance of the earth, the average radius of the earth and the longitude change function corresponding to each latitude interval in the target longitude and latitude data; determine the revisit number matrix according to the revisit number; based on the longitude change range, the number of grids that are not 0 in any row of the revisit number matrix, the number of grids included in the longitude change range corresponding to each latitude interval in the target longitude and longitude data, and the preset unit distance, determine the average coverage rate corresponding to each latitude interval; Wherein, the longitude variation range is expressed as: ； in, is the average distance from the Earth, is the average radius of the Earth, cos( ) is the longitude change function; The average coverage rate corresponding to each latitude interval is expressed as: ； in, is the preset unit distance, is the longitude variation range corresponding to each latitude interval in the target longitude and latitude data, is the number of grids that are not 0 in any row of the revisit count matrix, It is the number of grids contained in the longitude variation range corresponding to each latitude interval in the target longitude and latitude data.

6. The device according to claim 5, characterized in that The device also includes: The determination submodule is used to determine the coordinate data corresponding to each mirror reflection point in the revisit number matrix based on the longitude data and latitude data corresponding to each mirror reflection point; and determine the revisit number corresponding to a single grid based on the coordinate data and the single grid.

7. An electronic device, characterized in that: The electronic device comprises a processor, a memory, a communication interface and a communication bus, and the processor, the memory and the communication interface communicate with each other via the communication bus; The memory is used to store executable instructions, and the executable instructions enable the processor to execute the method for calculating the GNSS-R satellite coverage rate according to any one of claims 1 to 4.

8. A readable storage medium, characterized in that: The readable storage medium stores a program or instruction, and when the program or instruction is executed by the processor, the method for calculating the GNSS-R satellite coverage rate according to any one of claims 1 to 4 is implemented.

Images