Method and apparatus for providing integrity information for checking atmospheric correction parameters for correcting atmospheric disturbances for satellite navigation for a vehicle
11194051 · 2021-12-07
Assignee
Inventors
Cpc classification
G01S19/07
PHYSICS
G01S19/27
PHYSICS
G01S19/08
PHYSICS
International classification
G01S19/08
PHYSICS
G01S19/27
PHYSICS
Abstract
A method for providing integrity information for checking atmospheric correction parameters for the correction of atmospheric disturbances for satellite navigation for a vehicle includes reading state signals relating to a state of an atmosphere between at least one satellite receiver and at least one satellite of the at least one satellite receiver. Each state signal represents certain state data that are transmitted between a satellite and a satellite receiver. The method further includes using at least one satellite signal and that are dependent on a state of the atmosphere between the satellite and the satellite receiver. The method further includes determining the integrity information using the state data. A variation of the state data against time is analyzed.
Claims
1. A method for providing integrity information for checking atmospheric correction parameters for the correction of atmospheric disturbances for satellite navigation for a vehicle, the method comprising: reading state signals, relating to a state of an atmosphere between at least one satellite receiver and at least one satellite, using the at least one satellite receiver, wherein each state signal represents certain state data that are transmitted between the at least one satellite and the at least one satellite receiver using at least one satellite signal, and that are dependent on a state of the atmosphere between the at least one satellite and the at least one satellite receiver; and determining the integrity information, wherein the integrity information is based upon a variation of the state data over a predetermined period of time.
2. The method according to claim 1, further comprising: determining the integrity information using at least one of at least one static reference satellite receiver and a determining device that is external to the vehicle, and a satellite receiver of the vehicle.
3. The method according to claim 1, wherein the at least one satellite receiver is a static reference receiver, the method further comprising: using position information to determine the integrity information of the static reference satellite receiver.
4. The method according to claim 1, wherein the at least one satellite receiver is a satellite receiver of a vehicle, the method further comprising: transmitting at least one of the state signals and the integrity information to the satellite receiver of the vehicle.
5. The method according to claim 4, further comprising: using a result of a threshold value decision to implement the transmission of the state signals, wherein a comparison of the state signals or the integrity information with a threshold value related to the atmospheric disturbances is carried out when making the threshold value decision.
6. The method according to claim 5, further comprising: transmitting at least one of the state signals of selected static reference satellite receivers, and the integrity information based on the state signals of the selected static reference satellite receivers, wherein the selected static reference satellite receivers are selected depending on a geographic position of the satellite receiver of the vehicle.
7. The method according to claim 1, wherein: the state data represents at least one of a total electron content as a parameter of Earth's atmosphere, and quality indicators of the at least one satellite receiver; and determining the integrity information includes at least one of determining a gradient of the total electron content over a specified time interval, determining a standard deviation of a gradient of the total electron content over a specified time interval, and determining a change against time of the quality indicators over a specified time interval, wherein the total electron content is defined as a product of electron density and distance, measured in electrons per square meter, and wherein the quality indicators include at least one of standard deviations, and variances relating to a signal transmission using the at least one satellite receiver.
8. The method according to claim 1, wherein an apparatus is configured to at least one of carry out, and actuate steps of the method in corresponding units.
9. The method according to claim 1, wherein a computer program is configured to at least one of carry out, and actuate the method.
10. The method according to claim 9, wherein the computer program is stored on a non-transitory machine-readable memory medium.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1) Exemplary embodiments of the approach proposed here are represented in the drawings and are described in detail in the following description. In the figures:
(2)
(3)
(4)
(5)
(6)
(7)
DETAILED DESCRIPTION
(8) In the following description of favorable exemplary embodiments of the disclosure, the same or similar reference characters are used for the similarly acting elements represented in the different figures, wherein a repeated description of said elements is not provided.
(9)
(10) The static reference satellite receiver 130 and the user-side satellite receiver 160 are embodied to receive signals 115, in particular two-frequency GNSS signals 115, from the GNSS satellite 110. Furthermore, the static satellite receiver 130 is embodied to transmit state signals 135 or measurement data of the static reference satellite receiver 130 for network balancing, i.e. for example a parameter estimate including estimating an electron density and quality indicators, to the data processing device 140. The data processing device 140 is embodied to transmit correction data 145 or atmospheric correction parameters 145 to the user-side satellite receiver 160 by means of the transmission station 150 and the geostationary transmission satellite 120.
(11) In other words, in
(12) At the GNSS reference stations 130 in the correction data network, in particular the electron density or the so-called total electron content (TEC) for each signal path between satellites and receivers of the respective GNSS reference station 130 is determined using a common method. The TEC values are used as a basis for an ionospheric model by the correction service or correction data provider for the calculation of a global vertical total electron content (VTEC), i.e. a VTEC value can be calculated for each geographic length, geographic width and time, wherein interpolation estimates are used between the GNSS reference stations 130. For example, the use of polynomials, spherical surface functions, kriging interpolations, voxels or B-splines is common. Each of said methods results in particular in smoothing of the data between measurement base points, however. I.e. between measurement points the accuracy of the determined TEC depends on the interpolation method used. It can be assumed from this that the value is most accurate in the region of the individual GNSS reference stations 130.
(13) A region in the high atmosphere extending approximately between 50 kilometers and 1000 kilometers of altitude and transitioning above into the plasmasphere is described as the ionosphere. The ionosphere is characterized by an increased concentration of charged particles that exist as a result of ionization processes depending on the solar radiation. Electromagnetic signals are refracted in the ionosphere depending on the frequency and intensity of the ionization. As a result, signal delays sometimes occur in the observations of global satellite navigation systems, such as for example the GNSS, which must be corrected for positioning and navigation applications. With the use of geodetic two-frequency receivers, although signal combinations can be used in order to eliminate a first-order ionospheric effect, however the ionospheric state should be known in order to obtain very short convergence times for a highly accurate position solution. Single-frequency receivers resort to atmospheric models in order to correct the signals and to achieve accuracies of less than 1 meter. According to exemplary embodiments, as described with reference to the following figures, in particular such an ionospheric effect can be corrected reliably and accurately for positioning and navigation applications, wherein indications regarding the integrity of correction data 145 are also enabled.
(14) As a basis, for example the integrated electron density or the so-called total electron content (TEC) can be used as a correction parameter, because the same can be directly extracted from two-frequency GNSS observations, for example. In addition to other GNSS correction data, such as high-precision satellite track data, satellite clocks, tropospheric corrections or diverse biases, for example ionospheric correction data of so-called GNSS correction services can also be supplied in order to enable a centimeter-accurate position and short convergence times for a GNSS user. GNSS correction services calculate such correction data from regional or global station networks with GNSS receivers or GNSS reference stations 130. Whereas track, clock, tropospheric and bias parameters can be determined reliably based on the characteristics thereof, this is traditionally not the case or is only the case to a limited extent for ionospheric corrections. Because of physical conditions, in view of which a station distribution is inadequate in order to detect ionospheric disturbances at the location of the GNSS user, therefore traditionally no assurances regarding probabilities of error occurrence are provided by correction services. Instead of this, the responsibility for integrity assurance and thereby the detection of erroneous ionospheric corrections lies with the user. This concerns in particular small-scale temporal or spatial, dynamic and irregularly occurring phenomena, such as so-called traveling ionospheric disturbances (TIDs) differentiated as small-scale, medium-scale and large-scale TIDs, sudden ionospheric disturbances (SIDs), spread F-layer, sporadic E-layer, ionospheric storms and scintillations. According to exemplary embodiments however, as described with reference to the following figures, indications regarding the integrity of such correction data 145 may also be enabled.
(15) So-called ROT maps or ROTI maps (ROT=Rate of electron density, rate of TEC; TEC=Total Electron Content, electron density or total electron content−ROTI=rate of change of electron density index; rate of change of TEC index) can be generated, i.e. grid-based ROT values or ROTI values, which reflect ionospheric disturbances, as far as an underlying station distribution allows. Based on wired transmission paths, the transmission and monitoring of the ROTI would be conceivable, wherein for vehicle applications wireless links, such as the broadcast method by means of geostationary transmission satellites 120, are mainly used for transmission in this case and the correction data 145 are transmitted via L-Band. It should be noted here that correction parameters can be transmitted with the quality indicators or QI values thereof for a region that is visible from or covered by the satellite 110. In order to provide efficient data transfer at low bandwidth, in this case for example model coefficients are transmitted, from which the corresponding parameters can be extensively reconstructed.
(16) According to exemplary embodiments, it is described herein that in particular regional monitoring of ROT or ROTI and of QI values and the dynamics thereof in the surroundings of the user may be enabled while taking into account a transmission bandwidth and the costs.
(17)
(18)
(19)
(20) The static satellite receivers 130 are embodied to receive the satellite signals 115 from the satellites 110. The data processing device 140 is embodied to receive the state signals 135 from the static satellite receiver 130. The user-side satellite receiver 160 is disposed in a vehicle 470. The user-side satellite receiver 160 is embodied to receive the atmospheric correction parameters 145 for the correction of atmospheric disturbances during satellite navigation of the vehicle 470 from the data processing device 140. The satellite receiver 160 of the vehicle 470 is also embodied to receive satellite signals 115 from the satellite 110. Thus, on the user side ROT, ROTI or ROQ information can be determined in order to obtain an estimate of a current atmospheric or ionospheric activity and/or to check the plausibility of incoming correction data for atmospheric errors or ionospheric errors. ROQ stands for rate of quality in this case and will be described in detail further below in connection with
(21) According to the exemplary embodiment represented in
(22) The reader 441 is embodied to read the state signals 135 from the static satellite receiver 130. The state signals 135 are signals relating to a state of the atmosphere between the satellite receiver 130 and the satellite 110. Each state signal 135 represents certain state data that is dependent on a state of the atmosphere between a satellite 110 and a static satellite receiver 130 using at least one satellite signal 115 transmitted between the satellite 110 and the static satellite receiver 130.
(23) The determining device 442 is embodied to determine the integrity information 445 using the state data of the state signals 135 read by means of the reader 441. In this case, the determining device 442 is embodied to analyze a variation of the state data against time.
(24) According to one exemplary embodiment, the reader 441 is embodied to read state signals 135, the state data of which represent a total electron content as a parameter of the Earth's atmosphere and/or quality indicators of at least one satellite receiver 130. The total electron content is defined for this as the product of electron density and distance, measured in electrons per square meter. The quality indicators have standard deviations and/or variances relating to a signal transmission by means of the at least one satellite receiver. In this case, the determining device 442 is embodied to determine a gradient of the total electron content over a specified time interval, a standard deviation of a gradient of the total electron content over a specified time interval and/or a change against time of the quality indicators over a specified time interval as integrity information 445.
(25) According to the exemplary embodiment represented in
(26) According to one exemplary embodiment, the reader 441 is embodied to read position information of the static reference satellite receiver 130 as the source of the state signals 135, in particular together with the state signals 135. In this case, the determining device 442 is embodied to determine the integrity information 445 using the position information.
(27) According to the exemplary embodiment represented in
(28) The satellite receiver 160 of the vehicle 470 is embodied to check the atmospheric correction parameters 145 using the integrity information 445 that is provided by means of the apparatus 440 in order to produce atmospheric correction parameters that have been checked for integrity. Thus, the user-side satellite receiver 160 can be embodied to determine an own position of the user-side satellite receiver 160 or the vehicle 470 using the checked atmospheric correction parameters 145.
(29) According to the exemplary embodiment represented in
(30)
(31)
(32) In a reading step 610, state signals that are to be provided with the method 600 and that are related to a state of the atmosphere between at least one satellite receiver and at least one satellite of the at least one satellite receiver are read. Each state signal represents certain state data that are dependent on a state of the atmosphere between a satellite and a satellite receiver using at least one satellite signal that is transmitted between the satellite and the satellite receiver.
(33) Subsequently, the integrity information is determined using the state data in a determination step 620. In this case, a variation of the state data against time is analyzed.
(34) According to an exemplary embodiment, the determination step 610 is implemented by means of at least one static reference satellite receiver, by means of a determining device that is external to the vehicle and/or by means of a satellite receiver of the vehicle.
(35) According to a further exemplary embodiment, the providing method 600 also comprises a transmission step 630. In the transmission step 630, the state signals or the integrity information is or are transmitted to a satellite receiver of the vehicle. The transmission step 630 can be carried out before or after the determination step 620.
(36) With reference to the figures described above, in particular the
(37) First, the integrity monitoring based on ROT or ROTI by means of the apparatus 440 will be discussed.
(38) It can be assumed from this that interpolation errors of the TEC vertically above the static reference satellite receiver 130 are small because typically a plurality of satellites 110 that provide information for determining the TEC can be observed simultaneously in the immediate surroundings. It is therefore proposed to determine the rate of the TEC (ROT) at the positions of the static reference satellite receiver 130 as a quotient
(39)
wherein ΔVTEC is the TEC gradient against time at the zenith or the vertical TEC and Δt is the sampling against time. Furthermore, in this case Φ.sub.GF stands for a geometry-free linear combination of two-frequency phase measurements at frequencies f.sub.1 and f.sub.2. The use of an oblique TEC (STEC) between the static reference satellite receiver 130 and the satellite 110 is ignored here because spatial TEC variations (due to satellite movement) and temporal TEC variations would otherwise mix together in the calculation of the ROT and would provide the ROT satellite-specifically. It should be noted that calculations of the VTEC or STEC are based on common methods and are therefore not described in detail.
(40) In general, ROT thus represents a TEC gradient against time for a given time interval and is used as a basis for the calculation of ROTI:
(41)
(42) ROTI is defined as the standard deviation of ROT over a time interval or over a defined base number N. With ROT and ROTI there are two measures for the detection of ionospheric or atmospheric disturbances 202. Depending on a measurement frequency of the respective GNSS receiver 130 or static reference satellite receiver 130, scintillations can thus also be detected that can normally only be detected by high-cost scintillation receivers in conventional geodetic applications. In the course of developments for highly automated driving, for which high measurement frequencies and high accuracies using expensive sensors are required, for example with regard to motorway travel, it is thus conceivable that even inexpensive satellite receivers 130 or 160 could be used. Proposed interval lengths for the high-resolution temporal detection of ionospheric or atmospheric 202 disturbances are for example time steps of 1/30 s for the calculation of ROT and 5 minute time window lengths for determining ROTI.
(43) According to the exemplary embodiment represented in
(44) According to the exemplary embodiment represented in
(45) In both cases, ROT values or ROTI values for station positions of the correction networks are thus available as integrity information 445 for the user-side satellite receiver 160, being at least representative of the surroundings of a geographic position of the user-side satellite receiver 160, i.e. in both cases position information or station coordinates of the static reference satellite receivers 130 can therefore be provided. In particular, individual stationary reference satellite receivers 130 can thus be evaluated in the surroundings of the user-side satellite receiver 160 in order to determine whether discrepancies indicate local disturbances in the surroundings. The occurrence of increased ROT values or ROTI values can imply ionospheric or atmospheric disturbances 202, which could possibly result in erroneous correction data or atmospheric correction parameters 145.
(46) As a further measure, in the case of a two-frequency GNSS receiver as the user-side satellite receiver 160 it is possible to determine ROT or ROTI for the own position thereof without using correction data or atmospheric correction parameters 145 and comparing with values of the correction service arising from the surroundings. In this case, the VTEC at the geographic position of the user-side satellite receiver 160 can be derived from STEC measurements for a plurality of satellites 110. The following optional measures are proposed: ROT or ROTI should be derived from measurements on satellites 110 at higher elevation or close to the zenith in order to minimize error influences by means of near-horizon measurements. In principle, the ΔTEC is also determined here at satellites 110 in the surroundings of the user-side satellite receiver 160, i.e. an interpolation, for example a polynomial interpolation, can be carried out on the user side, by which the ΔVTEC can be determined at a geographic position of the user-side satellite receiver 160. Discrepancies between ROT or ROTI of the correction service or correction data provider and the self-determined integrity information 445 can imply that either there is a local ionospheric or atmospheric disturbance 202 in the area of the user-side satellite receiver 160 that is unrecognized by the correction service, for example a scintillation, or there is an atmospheric disturbance 202 at a static reference satellite receiver 130 that may reach the user-side satellite receiver 160 in the future, for example travelling ionospheric disturbances (TID). The latter can indicate that the correction data or atmospheric correction parameters 145 may already be corrupted.
(47) The integrity monitoring by means of the apparatus 440 based on TEC quality indicators will now be discussed.
(48) In contrast to ROT parameters or ROTI parameters, which can be used for the detection of receiver-specific disturbances, quality indicators (QI), typically standard deviations or variances that result from network balancing in the correction network and that are output together with the atmospheric correction parameters 145 or the TEC corrections, can be used in addition or alternatively for integrity monitoring. On the one hand, QI values can be used directly for this and on the other hand a change of the QI values at the static reference satellite receivers 130 (similarly to for ROT) can be monitored. Said parameter based on the change is referred to below as ROQ (Rate of Quality). The QI are not directly receiver-specific but are dependent on the network balancing and an underlying mathematical model. The standard deviations or variances that result from the network balancing and that are contained in a covariance matrix are often used as QI. For the case in which the QI is transmitted as the standard deviation, the following relationship results for example, similarly to for ROT:
(49)
wherein Δσ.sub.VTEC represents a change against time of the standard deviation that is determined from the network of the correction service or the correction data provider over the time window Δt. It is proposed to use the methods described for ROT or ROTI for quality indicators or ROQ. There are therefore the following corresponding cases:
(50) According to the exemplary embodiment represented in
(51) According to the exemplary embodiment represented in
(52) ROQ represents uncertainties in the determination of the correction data or atmospheric correction parameters 145, so that in this case deficiencies in the determination of the atmospheric correction parameters 145 because of ionospheric or atmospheric disturbances 202 at the individual static reference satellite receivers 130 can be revealed. This indication expands ROTI for a static reference satellite receiver 130 in the sense that the TEC gradients, which actually turn out to be small owing to smoothing in an interpolation or by down-weighting a disturbed static reference satellite receiver 130 in the adjustment and do not result in suspect ROTI values, in which ROQ values can be revealed, because ionospheric or atmospheric disturbances 202 are reflected in larger uncertainties of parameter estimation, i.e. larger standard deviations.
(53) Furthermore, it is also possible to calculate the ROQ for the geographic position of the user-side satellite receiver 160 so that here too there is a possibility of checking the plausibility of uncertainties. However, it should be noted that in contrast to ROT, ROQ at the position of the user-side satellite receiver 160 can only be determined from the atmospheric correction parameters 145, because there is no parameter estimation of the TEC here or no measurement redundancy for such a compensation and hence there are normally no specially calculated quality indicators.
(54) If an exemplary embodiment comprises an “and/or” connection between a first feature and a second feature, then this is to be read as that the exemplary embodiment comprises both the first feature and the second feature according to one embodiment and comprises either only the first feature or only the second feature according to a further embodiment.