Reference station network encryption method and system considering atmosphere delay experience precision information

Abstract

The invention discloses a reference station network encryption method and system considering atmospheric delay experience precision information. The method comprises the steps that grid division is carried out in a reference station encryption area according to a certain interval on the basis of geographic latitude and longitude; taking the positions of the divided grid points as the positions of the simulated users, and for each simulated user, selecting an existing reference station closest to the simulated user in the encrypted area to construct an observation model; according to the relative position relation between the user and the reference station, atmospheric delay experience precision information is calculated, and the random error of the observation value is estimated; determining an observation value error matrix according to the random error of the observation value, and calculating a corresponding precision matrix; and taking the sum of diagonal elements related to the coordinate values in the precision matrix as a precision index, and taking a plurality of grid points with the maximum precision index as positions for constructing the encryption base station.

Description

Technical Field

[0001] The present invention belongs to the field of global navigation system and positioning measurement technology, and specifically relates to a base station network encryption method that takes into account the empirical accuracy information of atmospheric delay. The method is used to calculate several areas where user positioning accuracy is easily affected by atmospheric delay as encryption station construction sites, thereby maximizing the service accuracy of the entire network after the encryption stations are built. Background Art

[0002] Satellite navigation positioning base station (base station) refers to a ground fixed observation station that conducts long-term continuous observation of satellite navigation signals and transmits the observation data to the data center in real time or regularly through communication facilities. The base station is the core component of the regional ground-based augmentation system. By obtaining continuous and stable high-precision satellite navigation observation data, it can provide high-precision positioning services for various types of users in the region. The current base station encryption method is mainly based on geographical, topographic maps and geological information, and uses field surveys and tests to select encryption sites. It only considers the base station construction environment and the base station's own data quality. It does not consider whether the base station site selection can effectively enhance the atmospheric modeling accuracy of the entire network and accurately enhance the positioning effect of users in areas with weak base station network service accuracy. For example, building a base station in an area with relatively calm atmosphere will not significantly improve the base station network service effect, resulting in a relative waste of resources. Summary of the invention

[0003] In order to solve the problem that the atmospheric enhancement degree of the base station network caused by the newly built base station is ignored in the existing encryption station site selection method during the encryption process of the base station network, resulting in the problem that the service effect of the whole network is not significantly improved, the present invention provides a base station network encryption method and system taking into account the empirical accuracy information of atmospheric delay. By considering the regional atmospheric delay accuracy information in the process of selecting the encrypted base stations, several areas where the user positioning accuracy is easily affected by the atmospheric delay are calculated as the construction sites of the encryption stations, thereby maximizing the improvement of the service accuracy of the whole network after the encryption stations are built.

[0004] According to one aspect of the present invention, there is provided a reference station network encryption method taking into account atmospheric delay empirical accuracy information, comprising:

[0005] Grid division is performed in the base station densification area according to the geographical longitude and latitude at the set intervals;

[0006] The divided grid point locations are used as the locations of simulated users. For each simulated user, the nearest existing reference station in the encrypted area is selected to construct an observation model.

[0007] According to the relative position relationship between the user and the reference station, the empirical accuracy information of atmospheric delay is calculated and the random error of the observation value is estimated;

[0008] Determine the observation error matrix from the random errors of the observations, and calculate the corresponding precision matrix;

[0009] The sum of the diagonal elements related to the coordinate values ​​in the precision matrix is ​​taken as the precision index, and several grid points with the largest precision index are used as the locations for the construction of encrypted reference stations.

[0010] As a further technical solution, grid division is performed in the base station densification area according to the geographical longitude and latitude at the set intervals, including:

[0011] For the base station network coverage area where encrypted stations need to be built, the east-west direction is the X-axis and the north-south direction is the Y-axis. The regional grid is divided according to the set intervals. The grid point position coordinates are expressed as:

[0012]

[0013] in: Represents the coordinates of the grid point; (X0, Y0) represents the coordinates of the grid point in the southwest corner of the area; n, m represent the number of rows and columns, n, m = 1, 2, 3...; ΔEW represents the grid spacing in the east-west direction; ΔNS represents the grid spacing in the north-south direction.

[0014] As a further technical solution, the divided grid point locations are used as the locations of simulated users. For each simulated user, the nearest existing reference station in the encrypted area is selected to construct an observation model, including:

[0015] Construct a user observation equation for each simulated user;

[0016] An indirect adjustment model is constructed based on the user observation equation.

[0017] As a further technical solution, calculate the empirical accuracy information of atmospheric delay and estimate the random error of the observation value for:

[0018]

[0019] Among them, D obs is the error of pseudorange or carrier observation; is the inter-station ionospheric error, is the inter-station tropospheric error, and The accuracy information together constitutes the empirical accuracy information of atmospheric delay.

[0020] As a further technical solution, the observation error matrix is ​​determined by the random error of the observation value, and the corresponding accuracy matrix is ​​calculated as follows

[0021]

[0022] Where A=(B T PB) -1 B T P, D L is the observation error matrix.

[0023] According to one aspect of the present invention, there is provided a reference station network encryption system taking into account empirical accuracy information of atmospheric delay, comprising:

[0024] The first processing module is used to perform grid division in the base station encryption area according to the geographical longitude and latitude at the set interval;

[0025] The second processing module is used to use the divided grid point location as the simulated user location, and for each simulated user, select the existing reference station closest to it in the encrypted area to build an observation model;

[0026] The third processing module is used to calculate the atmospheric delay empirical accuracy information and estimate the random error of the observation value according to the relative position relationship between the user and the reference station;

[0027] A fourth processing module is used to determine an observation value error matrix from the random errors of the observation values ​​and calculate a corresponding precision matrix;

[0028] The fifth processing module is used to take the sum of the diagonal elements related to the coordinate values ​​in the precision matrix as the precision index, and several grid points with the largest precision index as the locations for constructing encrypted reference stations.

[0029] According to one aspect of the present invention, there is provided a reference station network encryption device taking into account the empirical precision information of atmospheric delay, comprising a memory and a processor, wherein the memory stores program instructions executed by the processor, and the processor calls the program instructions to execute the steps of the reference station network encryption method taking into account the empirical precision information of atmospheric delay.

[0030] According to one aspect of the present specification, a non-transitory computer-readable storage medium is provided, wherein the non-transitory computer-readable storage medium stores computer instructions, wherein the computer instructions enable the computer to execute the steps of the reference station network encryption method taking into account the empirical accuracy information of atmospheric delay.

[0031] Compared with the prior art, the present invention has the following beneficial effects:

[0032] The present invention provides a base station network encryption method that takes into account the empirical accuracy information of atmospheric delay, so as to solve the problem that the existing encryption station site selection method in the encryption process of the base station network ignores the atmospheric enhancement degree of the new base station to the base station network, thereby resulting in the problem that the service effect of the whole network is not significantly improved. Specifically, the present invention performs grid division in the base station encryption area according to the geographical longitude and latitude at a certain interval; the divided grid point position is used as the position of the simulated user, and for each simulated user, the existing base station closest to it in the encryption area is selected to build an observation model; according to the relative position relationship between the user and the base station, the empirical accuracy information of atmospheric delay is calculated and the random error of the observation value is estimated; the observation value error matrix is ​​determined by the random error of the observation value, and the corresponding accuracy matrix is ​​calculated; the sum of the diagonal elements related to the coordinate value in the accuracy matrix is ​​used as the accuracy index, and several grid points with the largest accuracy index are used as the location of the encrypted base station construction. The present invention considers the regional atmospheric delay accuracy information in the process of selecting the encrypted base station, and calculates several areas where the user positioning accuracy is easily affected by the atmospheric delay as the construction site of the encrypted station, thereby maximizing the improvement of the service accuracy of the whole network after the encrypted station is built. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, a brief introduction is given below to the drawings used in the embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0034] Figure 1 A flowchart of a reference station network encryption method taking into account empirical accuracy information of atmospheric delay provided by an embodiment of the present invention.

[0035] Figure 2 A schematic diagram of grid division within a region provided by an embodiment of the present invention.

[0036] Figure 3 A schematic diagram of a reference station network encryption system taking into account empirical accuracy information of atmospheric delay provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0037] It should be noted that:

[0038] According to the relative position relationship between the reference station and the user, a random model of the user-side observation equation can be constructed, and the observation error can be introduced into the positioning model, so that the user-side positioning accuracy can be estimated in advance. Based on this, the present invention proposes a reference station network encryption method that takes into account the empirical accuracy information of atmospheric delay. In the process of selecting the encrypted station site, based on the existing reference station distribution data, based on the ionosphere, troposphere and observation value error empirical model, a user observation model covering the service range of the reference station network is established, and the user accuracy information throughout the network is extracted, and the location of the encrypted station construction is selected based on this.

[0039] The terms "including" and "having" and any variations thereof in the specification and claims of the present invention and the above-mentioned drawings are intended to cover non-exclusive inclusions, for example, a process, method, system, product or apparatus comprising a series of steps or units is not necessarily limited to the steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to these processes, methods, products or apparatuses.

[0040] In order to make the purpose, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments in the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention. In addition, the technical features in the various embodiments or single embodiments provided by the present invention are arbitrarily combined with each other to form a new technical solution. This combination is not restricted by the sequence of steps and / or the structural composition mode, but must be based on the ability of ordinary technicians in this field to achieve. When the combination of technical solutions is contradictory or cannot be achieved, it should be considered that this combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0041] An embodiment of the present invention provides a base station network encryption method that takes into account the empirical accuracy information of atmospheric delay. First, grid division is performed in the base station encryption area according to the geographical longitude and latitude at a certain interval; then, the divided grid point position is used as the position of the simulated user, and for each simulated user, the existing base station closest to it in the encryption area is selected to build an observation model; then, based on the relative position relationship between the user and the base station, the atmospheric delay empirical accuracy information is calculated and the random error of the observation value is estimated; then, the observation value error matrix is ​​determined by the random error of the observation value, and the corresponding accuracy matrix is ​​calculated; finally, the sum of the diagonal elements related to the coordinate value in the accuracy matrix is ​​used as the accuracy index, and several grid points with the largest accuracy index are used as the locations for the construction of the encrypted base stations.

[0042] The embodiment of the present invention considers the regional atmospheric delay accuracy information in the process of selecting encrypted reference stations, calculates several areas where the user positioning accuracy is easily affected by atmospheric delay as the construction sites of encrypted stations, thereby maximizing the service accuracy of the entire network after the encrypted stations are built. The embodiment of the present invention solves the problem that in the existing encryption station site selection method during the encryption of the reference station network, the degree of atmospheric enhancement of the reference station network by the newly built reference station is ignored, resulting in the problem that the service effect of the entire network is not significantly improved.

[0043] See also Figure 1 The reference station network encryption method taking into account the empirical accuracy information of atmospheric delay provided by the embodiment of the present invention includes the following specific steps:

[0044] Step 1: Grid division is performed in the base station densification area according to the geographic longitude and latitude at a certain interval.

[0045] For the base station network coverage area where encrypted stations need to be built, the east-west direction is the X-axis and the north-south direction is the Y-axis. The regional grid is divided at a certain interval. The grid point position coordinates can be expressed as:

[0046]

[0047] in, represents the coordinates of the grid point; (X0, Y0) represents the coordinates of the grid point in the southwest corner of the region; n, m represent the number of rows and columns, n, m = 1, 2, 3...; ΔEW represents the grid spacing in the east-west direction; ΔNS represents the grid spacing in the north-south direction. The schematic diagram of the grid division in the region is as follows Figure 2 shown.

[0048] Step 2: Use the divided grid point locations as the locations of simulated users. For each simulated user, select the existing reference station closest to it in the encrypted area to build an observation model.

[0049] First, the user observation equation can be expressed as:

[0050]

[0051] in:

[0052] i and j represent satellite numbers, r and u represent receiver numbers, and λ represents the wavelength of satellite signals;

[0053] is the double-difference pseudorange observation between ij satellites and ru receivers;

[0054] is the double-difference carrier observation between the ij satellites and the ru receivers;

[0055] is the double difference geometric distance between ij satellites and ru receivers;

[0056] is the tropospheric error between ij satellites and ru receivers, which is generally independent of the satellite number when calculated;

[0057] is the ionospheric error between the ij satellites and the ru receivers;

[0058] is the double difference ambiguity between the ij satellites and the ru receivers.

[0059] Next, an indirect adjustment model is constructed based on the user observation equation:

[0060]

[0061] in:

[0062] v is the residual vector of the equation;

[0063] It is the geometric double difference distance calculated from the a priori site coordinates;

[0064] is the sight vector of unit length;

[0065] Δx ru is the coordinate increment relative to the prior position.

[0066] The model can be simplified as:

[0067]

[0068] The P matrix can be calculated according to the altitude model or the signal-to-noise ratio model. B and L correspond to the corresponding parts in the indirect adjustment model, B is the design matrix, P is the weight matrix, and L is the observation vector.

[0069] Step 3: Calculate the empirical accuracy of atmospheric delay and estimate the random error of the observation value based on the relative position relationship between the user and the reference station.

[0070]

[0071] Where D obs is the error of pseudorange or carrier observation; is the inter-station ionospheric error; It should be noted that the empirical accuracy information of atmospheric delay described in the present invention refers to the accuracy information of the inter-station ionospheric error and the inter-station tropospheric error.

[0072] The corresponding error is calculated as follows:

[0073] When the observation type is pseudorange, D obs The value of is 1. When the observation type is carrier, D obs The value of is 0.001.

[0074] Inter-station ionospheric error The calculation is performed according to the following empirical model:

[0075]

[0076] Where exp(·) represents the logarithmic function, h ru represents the distance between the grid point user and the reference station, h0 is 200km, E i is the altitude angle of satellite i.

[0077] Tropospheric error between stations Estimate as follows:

[0078]

[0079] Among them, std(·) represents the standard deviation calculation, ZTD u,Δt It represents the ZTD (zenith tropospheric delay) sequence within the time period Δt calculated at the grid point u according to the empirical model (such as UNB model, GZTD model, GPT model).

[0080] Step 4: Determine the observation error matrix based on the random errors of the observations and calculate the corresponding precision matrix.

[0081] The observation error matrix is:

[0082]

[0083] n is the total number of pseudorange and carrier observation equations.

[0084] The corresponding precision matrix Calculate as follows:

[0085]

[0086] in:

[0087] A=(B T PB) -1 B T P.

[0088] Step 5: The sum of the diagonal elements related to the coordinate values ​​in the precision matrix is ​​used as the precision index, and several grid points with the largest precision index are used as the locations for the construction of encrypted reference stations.

[0089] It can be expressed as:

[0090]

[0091] Where D Coor is the element matrix related to the coordinate value, D Coor The size is 3×3, D Others is the element matrix related to other quantities to be determined.

[0092]

[0093] D Coor Each element in is related to the coordinate value. Here we only consider the diagonal elements. xx , D yy , D zz To a certain extent, it can better reflect the accuracy.

[0094] The precision index R of the corresponding grid point grid Calculate as follows:

[0095] R grid =D xx +D yy +D zz

[0096] According to R grid Arrange the grid points from large to small, take the largest grid points as the probabilistic locations for the construction of encryption stations, and complete the encryption of the base station network.

[0097] The implementation basis of each embodiment of the present invention is to implement programmed processing through a device with a processor function. Therefore, in engineering practice, the technical solutions and functions of each embodiment of the present invention are encapsulated into various modules. Based on this reality, on the basis of the above embodiments, an embodiment of the present invention provides a reference station network encryption system that takes into account the empirical accuracy information of atmospheric delay, and the system is used to execute the reference station network encryption method that takes into account the empirical accuracy information of atmospheric delay in the above method embodiment.

[0098] See also Figure 3The system includes: a first processing module, which is used to perform grid division in the base station encryption area according to the set interval based on the geographical longitude and latitude; a second processing module, which is used to use the divided grid point position as the simulated user's position, and for each simulated user, select the existing base station closest to it in the encryption area to build an observation model; a third processing module, which is used to calculate the atmospheric delay empirical accuracy information and estimate the random error of the observation value based on the relative position relationship between the user and the base station; a fourth processing module, which is used to determine the observation value error matrix from the random error of the observation value, and calculate the corresponding accuracy matrix; a fifth processing module, which is used to use the sum of the diagonal elements related to the coordinate value in the accuracy matrix as the accuracy index, and a number of grid points with the largest accuracy index as the location of the encrypted base station construction.

[0099] The embodiment of the present invention provides a reference station network encryption system taking into account the empirical accuracy information of atmospheric delay. The system aims to solve the problem that the existing encryption station site selection method in the reference station network encryption process ignores the atmospheric enhancement degree of the newly built reference station to the reference station network, thereby resulting in an insignificant improvement in the service effect of the entire network. Figure 3 Several modules in it consider the regional atmospheric delay accuracy information in the process of selecting encrypted base stations, calculate several areas where user positioning accuracy is easily affected by atmospheric delay as the construction sites of encrypted stations, so as to maximize the service accuracy of the entire network after the encrypted stations are built.

[0100] It should be noted that the system embodiment provided by the present invention is used to implement the methods in the above method embodiment as well as the methods in other method embodiments provided by the present invention. The only difference is that the corresponding functional modules are set. The principle is basically the same as the principle of the above system embodiment provided by the present invention. As long as the technical personnel in this field refer to the specific technical solutions in other method embodiments on the basis of the above system embodiment, obtain the corresponding technical means and the technical solutions composed of these technical means by combining the technical features, and on the premise of ensuring the practicality of the technical solutions, improve the modules in the above system embodiment to obtain the corresponding system class embodiments, which are used to implement the methods in other method class embodiments. For example:

[0101] Based on the content of the above system embodiment, as a preferred embodiment, in the reference station network encryption system taking into account the empirical accuracy information of atmospheric delay provided in the embodiment of the present invention, the first processing module is also used to execute the following instructions:

[0102] For the base station network coverage area where encrypted stations need to be built, the east-west direction is the X-axis and the north-south direction is the Y-axis. The regional grid is divided according to the set intervals. The grid point position coordinates are expressed as:

[0103]

[0104] in: Represents the coordinates of the grid point; (X0, Y0) represents the coordinates of the grid point in the southwest corner of the area; n, m represent the number of rows and columns, n, m = 1, 2, 3...; ΔEW represents the grid spacing in the east-west direction; ΔNS represents the grid spacing in the north-south direction.

[0105] Based on the content of the above system embodiment, as a preferred embodiment, a reference station network encryption system that takes into account the empirical accuracy information of atmospheric delay is provided in the embodiment of the present invention, and the second processing module is also used to execute the following instructions: construct a user observation equation for each simulated user; and construct an indirect adjustment model based on the user observation equation.

[0106] Based on the content of the above system embodiment, as a preferred embodiment, a reference station network encryption system taking into account the empirical accuracy information of atmospheric delay is provided in the embodiment of the present invention, and the third processing module is also used to execute the following instructions: calculate the empirical accuracy information of atmospheric delay and estimate the random error of the observation value for:

[0107]

[0108] Among them, D obs is the error of pseudorange or carrier observation; is the inter-station ionospheric error, is the inter-station tropospheric error, and The accuracy information together constitutes the empirical accuracy information of atmospheric delay.

[0109] Based on the content of the above system embodiment, as a preferred embodiment, a reference station network encryption system taking into account the empirical accuracy information of atmospheric delay is provided in the embodiment of the present invention, wherein the fourth processing module is further used to execute the following instructions: determine the observation value error matrix from the random error of the observation value, and calculate the corresponding accuracy matrix in the following manner

[0110]

[0111] Where A=(B T PB) -1 B T P, D L is the observation error matrix.

[0112] Based on the same inventive concept as the above-mentioned embodiment, an embodiment of the present invention provides a reference station network encryption device taking into account the empirical accuracy information of atmospheric delay, including a memory and a processor, wherein the memory stores program instructions executed by the processor, and the processor calls the program instructions to execute the steps of the reference station network encryption method taking into account the empirical accuracy information of atmospheric delay.

[0113] In an embodiment of the present invention, the memory may be a non-volatile memory, such as a hard disk drive (HDD) or a solid-state drive (SSD), etc., or a volatile memory (volatile memory), such as a random-access memory (RAM). The memory is any other medium that can be used to carry or store the desired program code in the form of instructions or data structures and can be accessed by a computer, but is not limited thereto. The memory in an embodiment of the present invention may also be a circuit or any other device that can implement a storage function, for storing program instructions and / or data.

[0114] In the embodiments of the present invention, the processor may be a general-purpose processor, a digital signal processor, an application-specific integrated circuit, a field programmable gate array or other programmable logic device, a discrete gate or transistor logic device, or a discrete hardware component, and may implement or execute the methods, steps, and logic block diagrams disclosed in the embodiments of the present invention. A general-purpose processor may be a microprocessor or any conventional processor, etc. The steps of the method disclosed in the embodiments of the present invention may be directly embodied as being executed by a hardware processor, or may be executed by a combination of hardware and software modules in the processor.

[0115] Based on the same inventive concept as the above-mentioned embodiment, an embodiment of the present invention provides a non-transitory computer-readable storage medium, wherein the non-transitory computer-readable storage medium stores computer instructions, and the computer instructions enable the computer to execute the steps of the reference station network encryption method taking into account the empirical accuracy information of atmospheric delay.

[0116] In summary of the above embodiments, the present invention discloses a base station network encryption method that takes into account the empirical accuracy information of atmospheric delay, including grid division in the base station encryption area according to the geographic longitude and latitude at a certain interval; using the divided grid point position as the simulated user's position, for each simulated user, selecting the existing base station closest to it in the encryption area to build an observation model; calculating the empirical accuracy information of atmospheric delay and estimating the random error of the observation value based on the relative position relationship between the user and the base station; determining the observation value error matrix from the random error of the observation value, and calculating the corresponding accuracy matrix; taking the sum of the diagonal elements related to the coordinate value in the accuracy matrix as the accuracy index, and using several grid points with the largest accuracy index as the location of the encrypted base station construction. The present invention takes into account the regional atmospheric delay accuracy information in the process of selecting the encrypted base station, and calculates several areas where the user positioning accuracy is easily affected by the atmospheric delay as the construction sites of the encrypted stations, thereby maximizing the service accuracy of the entire network after the encrypted stations are built.

[0117] It should be understood that parts not elaborated in detail in this specification belong to the prior art.

[0118] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or replace some or all of the technical features therein with equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the technical solutions of the embodiments of the present invention.

Claims

1. A reference station network encryption method taking into account the empirical accuracy information of atmospheric delay, characterized in that: include: Grid division is performed in the base station densification area according to the geographical longitude and latitude at the set intervals; The divided grid point locations are used as the locations of simulated users. For each simulated user, the nearest existing reference station in the encrypted area is selected to build an observation model. According to the relative position relationship between the user and the reference station, the empirical accuracy information of atmospheric delay is calculated and the random error of the observation value is estimated; Determine the observation error matrix from the random errors of the observations, and calculate the corresponding precision matrix; The sum of the diagonal elements related to the coordinate values ​​in the precision matrix is ​​taken as the precision index, and several grid points with the largest precision index are used as the locations for the construction of encrypted reference stations.

2. The reference station network encryption method according to claim 1, wherein: Grid division is performed in the base station densification area according to the set intervals based on the geographic longitude and latitude, including: For the base station network coverage area where encrypted stations need to be built, the east-west direction is the X-axis and the north-south direction is the Y-axis. The regional grid is divided according to the set intervals. The grid point position coordinates are expressed as: in: Represents the coordinates of the grid point; (X0, Y0) represents the coordinates of the grid point in the southwest corner of the area; n, m represent the number of rows and columns, n, m = 1, 2, 3...; ΔEW represents the grid spacing in the east-west direction; ΔNS represents the grid spacing in the north-south direction.

3. The reference station network encryption method according to claim 1, characterized in that: The divided grid point locations are used as the locations of simulated users. For each simulated user, the nearest existing reference station in the encrypted area is selected to construct an observation model, including: Construct a user observation equation for each simulated user; An indirect adjustment model is constructed based on the user observation equation.

4. The reference station network encryption method according to claim 1, wherein: Calculate empirical accuracy information for atmospheric delays and estimate random errors in observations for: Among them, D obs is the error of pseudorange or carrier observation; is the inter-station ionospheric error, is the inter-station tropospheric error, and The accuracy information together constitutes the empirical accuracy information of atmospheric delay.

5. The reference station network encryption method according to claim 1, characterized in that: The observation error matrix is ​​determined by the random error of the observation value, and the corresponding accuracy matrix is ​​calculated as follows Where A=(B T PB) -1 B T P, D L is the observation error matrix.

6. A base station network encryption system taking into account the empirical accuracy information of atmospheric delay, characterized in that: include: The first processing module is used to perform grid division in the base station encryption area according to the geographical longitude and latitude at the set interval; The second processing module is used to use the divided grid point location as the simulated user location, and for each simulated user, select the existing reference station closest to it in the encrypted area to build an observation model; The third processing module is used to calculate the atmospheric delay empirical accuracy information and estimate the random error of the observation value according to the relative position relationship between the user and the reference station; A fourth processing module is used to determine an observation value error matrix from the random errors of the observation values ​​and calculate a corresponding precision matrix; The fifth processing module is used to take the sum of the diagonal elements related to the coordinate values ​​in the precision matrix as the precision index, and several grid points with the largest precision index as the locations for constructing encrypted reference stations.

7. A reference station network encryption device taking into account the empirical accuracy information of atmospheric delay, characterized in that: The invention comprises a memory and a processor, wherein the memory stores program instructions executed by the processor, and the processor calls the program instructions to execute the steps of the reference station network encryption method taking into account the empirical accuracy information of atmospheric delay as described in any one of claims 1 to 5.

8. A non-transitory computer-readable storage medium, characterized in that: The non-transitory computer-readable storage medium stores computer instructions, and the computer instructions enable the computer to execute the steps of the reference station network encryption method taking into account the atmospheric delay empirical accuracy information as described in any one of claims 1 to 5.

Images