FR3141254A1 - Detection and correction of drift of a navigation device by star sighting - Google Patents
Detection and correction of drift of a navigation device by star sighting Download PDFInfo
- Publication number
- FR3141254A1 FR3141254A1 FR2210995A FR2210995A FR3141254A1 FR 3141254 A1 FR3141254 A1 FR 3141254A1 FR 2210995 A FR2210995 A FR 2210995A FR 2210995 A FR2210995 A FR 2210995A FR 3141254 A1 FR3141254 A1 FR 3141254A1
- Authority
- FR
- France
- Prior art keywords
- star
- coordinates
- stars
- navigation device
- calculation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000012937 correction Methods 0.000 title claims abstract description 10
- 238000001514 detection method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 39
- 238000004364 calculation method Methods 0.000 claims abstract description 33
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 238000012545 processing Methods 0.000 claims abstract description 6
- 230000035945 sensitivity Effects 0.000 claims description 16
- 235000019892 Stellar Nutrition 0.000 claims description 5
- 101000582320 Homo sapiens Neurogenic differentiation factor 6 Proteins 0.000 claims description 3
- 102100030589 Neurogenic differentiation factor 6 Human genes 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims description 2
- 238000012795 verification Methods 0.000 abstract description 3
- 238000010187 selection method Methods 0.000 description 6
- 230000004807 localization Effects 0.000 description 3
- 238000005070 sampling Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/16—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using electromagnetic waves other than radio waves
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
Détection et correction de dérive d’un dispositif de navigation par visée stellaire Procédé de correction de dérive d’un dispositif de navigation d’un véhicule à l’aide d’un moyen de géolocalisation par visée stellaire, le procédé comprenant :- une étape de pré-sélection d’étoiles observables à utiliser, puis- une étape de mesures comprenant un relevé des coordonnées géographiques du viseur par le dispositif de navigation, puis, pour chacune des étoiles sélectionnées, une vérification de l’accessibilité de l’étoile et une mesure des coordonnées horizontales, puis- une étape de traitement comportant, pour toutes les étoiles pour lesquelles des coordonnées horizontales mesurées ont été enregistrées un calcul des coordonnées géographiques du viseur stellaire, et une estimation de l’erreur entre les coordonnées géographiques calculées et les coordonnées géographiques relevées par le dispositif de navigation lors de la visée de l’étoile, puis- une étape d’estimation de l’erreur de dérive du dispositif de navigation, et- une étape de correction de la dérive du dispositif de navigation. Figure pour l’abrégé : Fig. 8Detection and correction of drift of a navigation device by star sighting Method of correcting drift of a navigation device of a vehicle using a means of geolocation by star sighting, the method comprising: - a step pre-selection of observable stars to be used, then- a measurement step comprising a reading of the geographical coordinates of the viewfinder by the navigation device, then, for each of the selected stars, a verification of the accessibility of the star and a measurement of the horizontal coordinates, then- a processing step comprising, for all the stars for which measured horizontal coordinates have been recorded, a calculation of the geographical coordinates of the star finder, and an estimation of the error between the calculated geographical coordinates and the geographical coordinates noted by the navigation device when aiming at the star, then - a step of estimating the drift error of the navigation device, and - a step of correcting the drift of the navigation device. Figure for abstract: Fig. 8
Description
L’invention concerne la géolocalisation de navires de surfaces disposant d’une centrale de navigation à leur bord, et plus particulièrement un système et un procédé de détection et de correction d’erreur d’un système de géolocalisation embarqué à bord d’un navire de surface, notamment lors de son évolution en pleine mer.The invention relates to the geolocation of surface ships having a navigation center on board, and more particularly a system and a method for detecting and correcting errors in a geolocation system on board a ship. surface, particularly during its evolution in the open sea.
Une centrale de navigation est un système de capteurs inertiels fournissant une estimation de la position et de l’altitude du porteur sur Terre. L’atout principal d’une centrale de navigation est qu’elle est autonome, c’est-à-dire que son estimation de position ne dépend d’aucun autre système externe. Cependant, son principal inconvénient est que l’estimation de position se détériore petit à petit dans le temps, de l’ordre de quelques kilomètres par jour. On dit alors que l’estimation de position dérive dans le temps.A navigation center is a system of inertial sensors providing an estimate of the wearer's position and altitude on Earth. The main advantage of a navigation center is that it is autonomous, that is to say that its position estimation does not depend on any other external system. However, its main disadvantage is that the position estimation deteriorates little by little over time, on the order of a few kilometers per day. We then say that the position estimate drifts over time.
Pour estimer cette erreur de position, il est nécessaire de disposer épisodiquement d’une source d’information externe fiable. Généralement, la référence externe utilisée est un système de positionnement par satellite du fait de leur précision et de leur disponibilité. Ces systèmes de navigation globale par satellite sont désignés par le sigle GNSS pour Global Navigation Satellite Systems en anglais. L’exemple le plus connu de GNSS est le système Américain GPS (Global Positioning System). Cependant, il est possible que le GNSS soit inaccessible ou que l’on reçoive des données GNSS leurrées : dans ce cas, le système de localisation du navire peut être défaillant, et, dans le cas d’un navire de combat, le système de combat peut être défaillant.To estimate this position error, it is necessary to occasionally have a reliable external source of information. Generally, the external reference used is a satellite positioning system due to their precision and availability. These global navigation satellite systems are designated by the acronym GNSS for Global Navigation Satellite Systems in English. The best-known example of GNSS is the American GPS (Global Positioning System). However, it is possible that the GNSS is inaccessible or that deceptive GNSS data is received: in this case, the ship's localization system may be faulty, and, in the case of a combat ship, the tracking system combat may fail.
Plusieurs solutions existent pour retrouver la position d’un navire sans référence de type GNSS disponible.Several solutions exist to find the position of a vessel without a GNSS type reference available.
La première solution connue nécessite de disposer d’un amer géo-référencé, c’est-à-dire un point de repère fixe et identifiable sans ambiguïté tel qu’un phare par exemple. Cette solution est précise, mais applicable qu’en des circonstances particulières en présence d’un amer en navigation côtière avec utilisation par exemple de l'almanach du marin Breton (L'almanach du marin Breton, 2022).The first known solution requires having a geo-referenced landmark, that is to say a fixed and unambiguously identifiable landmark such as a lighthouse for example. This solution is precise, but only applicable in particular circumstances in the presence of a landmark in coastal navigation with use for example of the Breton sailor's almanac (L'almanac du marin Breton, 2022).
La seconde solution connue nécessite de disposer d’éphémérides du soleil, d’un garde-temps et d’un sextant. Cette solution est peu précise (de l’ordre de 3 nautiques), et n’est applicable que lorsque le soleil est visible et au zénith. Le principe de l’estimation de position par une mesure au sextant est expliqué, par exemple, dans document d’Umland, de 2019 intitulé « A Short Guide to Celestial Navigation ».The second known solution requires having solar ephemeris, a timepiece and a sextant. This solution is not very precise (of the order of 3 nautical), and is only applicable when the sun is visible and at the zenith. The principle of position estimation by sextant measurement is explained, for example, in Umland's 2019 document entitled “A Short Guide to Celestial Navigation”.
La troisième solution connue nécessite de disposer d’éphémérides d’étoiles, d’un garde-temps, d’un système de visée précis et d’un inclinomètre (système d’accéléromètres donnant la verticale du lieu). C’est ce que l’on appelle la navigation stellaire. Le principe de cette méthode est décrit, par exemple, dans (Wei, Cui, Wang, & Wan, 2019). La compensation des erreurs de visée utilise la méthode QUEST (QUaternion ESTimator). Cependant, il n’y a pas de garantie sur la précision de l’estimateur de position.The third known solution requires having star ephemeris, a timepiece, a precise aiming system and an inclinometer (accelerometer system giving the vertical of the place). This is called stellar navigation. The principle of this method is described, for example, in (Wei, Cui, Wang, & Wan, 2019). Aiming error compensation uses the QUEST method (QUaternion ESTimator). However, there is no guarantee on the accuracy of the position estimator.
Utiliser les étoiles pour retrouver sa position sur Terre est bien connu. En revanche, les méthodes existantes peuvent être imprécises car les méthodes de discrimination des étoiles à observer ne sont pas assez développées. En effet, on peut remarquer que certaines visées d’étoiles retournent des précisions de localisation peu fiables : soit parce qu’elles sont situées dans des régions singulières de l’espace, c’est-à-dire intrinsèquement non-observables, soit parce qu’elles sont dans des régions de l’espace où elles vont amplifier les défauts du viseur.Using the stars to find your position on Earth is well known. On the other hand, existing methods can be imprecise because the methods for discriminating which stars to observe are not sufficiently developed. Indeed, we can notice that certain star sightings return unreliable location precisions: either because they are located in singular regions of space, that is to say intrinsically unobservable, or because that they are in regions of space where they will amplify the defects of the viewfinder.
L'invention vise à pallier les inconvénients mentionnés ci-dessus et à fournir une solution pour détecter ce type de dérive d’un dispositif de navigation rapidement, en pleine mer par exemple, ou retrouver sa position en l’absence de signal du dispositif de navigation.The invention aims to overcome the drawbacks mentioned above and to provide a solution for detecting this type of drift of a navigation device quickly, in the open sea for example, or finding its position in the absence of a signal from the navigation device. navigation.
La présente invention permet en outre d’identifier des régions de l’espace défavorables à la géolocalisation par visée stellaire, ainsi que des régions plus favorables pour la géolocalisation, moins sujettes aux erreurs de visée.The present invention also makes it possible to identify regions of space unfavorable for geolocation by star sighting, as well as regions more favorable for geolocation, less subject to sighting errors.
Ce but est atteint tout d’abord grâce à un procédé de sélection d’étoiles observables à utiliser pour une géolocalisation par un véhicule, notamment un navire en mer, à partir d’un moyen de géolocalisation par visée stellaire associé à au moins un dispositif de navigation, la sélection d’étoiles étant réalisée à partir d’étoiles recensées dans un catalogue où chaque étoile est associée à des coordonnées équatoriales à une date spécifique, les coordonnées équatoriales comportant une ascension droite, α, et une déclinaison, δ.This goal is achieved first of all thanks to a method of selecting observable stars to be used for geolocation by a vehicle, in particular a ship at sea, from a means of geolocation by stellar sighting associated with at least one device navigation, the selection of stars being made from stars listed in a catalog where each star is associated with equatorial coordinates on a specific date, the equatorial coordinates comprising a right ascension, α, and a declination, δ.
Le procédé comprend un relevé des coordonnées géographiques d’un viseur stellaire à un instant de départ par le dispositif de navigation, les coordonnées géographiques du viseur stellaire comprenant une longitude, L, et une latitude, ϕ.The method comprises a reading of the geographical coordinates of a star finder at a starting time by the navigation device, the geographical coordinates of the star finder comprising a longitude, L, and a latitude, ϕ.
Puis il comprend, pour chaque étoile du catalogue, un calcul des coordonnées de visée de l’étoile, une détermination de l’accessibilité de l’étoile à partir des coordonnées de visée calculées, et, si l’étoile est accessible, un calcul d’un paramètre de sensibilité d’observation de l’étoile, les coordonnées de visée de l’étoile correspondant aux coordonnées horizontales de l’étoile par rapport à la position géographique du viseur stellaire, les coordonnées horizontales comportant un azimut, A, et une élévation, h, puis un stockage des coordonnées de visée et du paramètre de sensibilité d’observation calculés pour chaque étoile, et un classement des étoiles ainsi stockée par ordre croissant de paramètre de sensibilité de l’étoile, les étoiles dont le paramètre de sensibilité est nul étant exclues du classement, et une sélection d’un nombre déterminé d’étoiles observables, les étoiles étant sélectionnées selon l’ordre de classement des étoiles stockées.Then it includes, for each star in the catalog, a calculation of the sighting coordinates of the star, a determination of the accessibility of the star from the calculated sighting coordinates, and, if the star is accessible, a calculation of a sensitivity parameter for observing the star, the sighting coordinates of the star corresponding to the horizontal coordinates of the star relative to the geographical position of the star finder, the horizontal coordinates comprising an azimuth, A, and an elevation, h, then storage of the sighting coordinates and the observation sensitivity parameter calculated for each star, and a classification of the stars thus stored in ascending order of sensitivity parameter of the star, the stars whose sensitivity is zero being excluded from the classification, and a selection of a determined number of observable stars, the stars being selected according to the order of classification of the stored stars.
Selon un premier mode de mise en œuvre de l’invention du procédé de sélection, le calcul des coordonnées de visée de chaque étoile peut comprendre un calcul des coordonnées horizontales de l’étoile à partir des coordonnées équatoriales de l’étoile et de la position géographique du viseur, les coordonnées équatoriales comprenant un azimut, A, de l’étoile et une élévation, h, de l’étoile, et
la détermination de l’accessibilité de l’étoile peut comprendre une comparaison de l’élévation calculée pour l’étoile à un seuil d’élévation minimale, l’étoile étant inaccessible si son élévation calculée est inférieure au seuil d’élévation minimale.According to a first mode of implementation of the invention of the selection method, the calculation of the sighting coordinates of each star can include a calculation of the horizontal coordinates of the star from the equatorial coordinates of the star and the position geographical position of the viewfinder, the equatorial coordinates including an azimuth, A, of the star and an elevation, h, of the star, and
determining the accessibility of the star may include comparing the calculated elevation for the star to a minimum elevation threshold, the star being inaccessible if its calculated elevation is less than the minimum elevation threshold.
De préférence, les coordonnées géographiques du viseur utilisées pour calculer les coordonnées horizontales de l’étoile dépendent de l’instant auquel le calcul est réalisé, les coordonnées géographiques du viseur étant estimées à chaque instant à partir des coordonnées géographiques du viseur relevées à l’instant de départ et d’un calcul tenant compte du cap et de la vitesse du véhicule à l’instant du relevé.Preferably, the geographic coordinates of the viewfinder used to calculate the horizontal coordinates of the star depend on the instant at which the calculation is carried out, the geographic coordinates of the viewfinder being estimated at each instant from the geographic coordinates of the viewfinder noted at departure time and a calculation taking into account the heading and speed of the vehicle at the time of the reading.
Dans un deuxième mode de mise en œuvre du procédé de sélection, le calcul des coordonnées de visée de l’étoile à partir des coordonnées géographiques du viseur stellaire et des coordonnées équatoriales de l’étoile considérée peut être réalisé à partir du système d’équations (S) :In a second mode of implementation of the selection method, the calculation of the sighting coordinates of the star from the geographical coordinates of the star finder and the equatorial coordinates of the star considered can be carried out from the system of equations (S):
[Math 1] S
Avec GAST désignant le temps sidéral apparent de Greenwich et LAST désignant le temps sidéral apparent local.With GAST denoting Greenwich apparent sidereal time and LAST denoting local apparent sidereal time.
Dans un troisième mode de mise en œuvre du procédé de sélection, le paramètre, Xi, de sensibilité d’observation de l’étoile, Sti, calculé dépend de préférence uniquement de l’azimut, A, de l’élévation, h, et de la déclinaison, δ, de l’étoile considérée.In a third mode of implementation of the selection method, the parameter, X i , of observation sensitivity of the star, St i , calculated preferably depends only on the azimuth, A, of the elevation, h , and the declination, δ, of the star considered.
Dans un quatrième mode de mise en œuvre du procédé de sélection, le paramètre, Xi, de sensibilité d’observation de l’étoile, Sti, est calculé via l’équation :In a fourth mode of implementation of the selection method, the parameter, X i , of observation sensitivity of the star, St i , is calculated via the equation:
[Math 2]
avec [Math 3]
with [Math 3]
[Math 4]
[Math 5]
[Math 6]
avec
Le but de l’invention est atteint dans un second temps grâce à un procédé de correction de dérive d’un dispositif de navigation d’un véhicule à l’aide d’un moyen de géolocalisation par visée stellaire.The aim of the invention is achieved in a second step thanks to a method of correcting drift of a navigation device of a vehicle using a means of geolocation by star sighting.
Le procédé de correction comprend une étape de pré-sélection d’étoiles comportant une sélection d’étoiles à utiliser par le dispositif de géolocalisation par visée stellaire à partir du procédé de sélection d’étoiles observables à utiliser telle que définie ci-dessus.The correction method includes a star pre-selection step comprising a selection of stars to be used by the geolocation device by star sighting from the method of selecting observable stars to use as defined above.
Puis le procédé comprend une étape de mesures comprenant, tout d’abord, un relevé des coordonnées géographiques du viseur stellaire par le dispositif de navigation, les coordonnées géographiques du viseur stellaire comprenant une longitude (L) et une latitude (ϕ), puis, pour chacune des étoiles sélectionnées, une vérification de l’accessibilité de l’étoile et une mesure des coordonnées horizontales de l’étoile, les coordonnées horizontales mesurées étant enregistrées en association avec les coordonnées équatoriales correspondantes pour l’étoile et l’instant de mesure, puis une étape de traitement comportant, pour toutes les étoiles pour lesquelles des coordonnées horizontales mesurées ont été enregistrées un calcul des coordonnées géographiques du viseur stellaire à partir des coordonnées horizontales mesurées enregistrées pour l’étoile, de l’instant de mesure enregistré, des coordonnées équatoriales enregistrées et du système d’équations (S) défini ci-dessus selon [Math 1], une estimation de l’erreur entre les coordonnées géographiques calculées et les coordonnées géographiques relevées par le dispositif de navigation lors de la visée de l’étoile, puis une étape d’estimation de l’erreur de dérive du dispositif de navigation à partir des erreurs estimées sur les coordonnées géographiques, et une étape de correction de la dérive du dispositif de navigation à partir de l’erreur de dérive estimée.Then the method comprises a measurement step comprising, first of all, a reading of the geographical coordinates of the star finder by the navigation device, the geographical coordinates of the star finder including a longitude (L) and a latitude (ϕ), then, for each of the selected stars, a check of the accessibility of the star and a measurement of the horizontal coordinates of the star, the measured horizontal coordinates being recorded in association with the corresponding equatorial coordinates for the star and the time of measurement , then a processing step comprising, for all the stars for which measured horizontal coordinates have been recorded a calculation of the geographical coordinates of the star finder from the measured horizontal coordinates recorded for the star, the recorded measurement instant, recorded equatorial coordinates and the system of equations (S) defined above according to [Math 1], an estimate of the error between the calculated geographic coordinates and the geographic coordinates noted by the navigation device when aiming the star, then a step of estimating the drift error of the navigation device from the errors estimated on the geographic coordinates, and a step of correcting the drift of the navigation device from the estimated drift error.
Dans un mode de mise en œuvre du procédé de correction, l’étape de mesures comprend un calcul des coordonnées horizontales de l’étoile sélectionnée à partir de l’instant de calcul, des coordonnées géographiques du viseur stellaire à l’instant de calcul, et des coordonnées équatoriales indiquées dans un catalogue où chaque étoile est associée à des coordonnées équatoriales à une date spécifique appelée par exemple J2000, les coordonnées équatoriales comportant une ascension droite (α) et une déclinaison (δ) et les coordonnées horizontales comportant un azimut (A) et une élévation (h), puis une orientation d’un viseur stellaire selon les coordonnées horizontales calculées, une vérification de l’accessibilité de l’étoile sélectionnée, si l’étoile sélectionnée est inaccessible, une nouvelle réalisation des trois étapes précédentes avec une nouvelle étoile de la sélection d’étoiles, si l’étoile est accessible, un recalage de l’étoile sélectionnée au centre du viseur stellaire, une mesure par le dispositif de navigation des coordonnées horizontales de l’étoile visée par le viseur stellaire après recalage, et un enregistrement de l’instant de mesure, des coordonnées horizontales recalées ainsi mesurées, des coordonnées équatoriales associées et des coordonnées géographiques du viseur à l’instant de mesure.In one mode of implementing the correction method, the measurement step includes a calculation of the horizontal coordinates of the selected star from the calculation instant, of the geographical coordinates of the star finder at the calculation instant, and equatorial coordinates indicated in a catalog where each star is associated with equatorial coordinates on a specific date called for example J2000, the equatorial coordinates comprising a right ascension (α) and a declination (δ) and the horizontal coordinates comprising an azimuth ( A) and an elevation (h), then an orientation of a star finder according to the calculated horizontal coordinates, a verification of the accessibility of the selected star, if the selected star is inaccessible, a new realization of the three previous steps with a new star from the selection of stars, if the star is accessible, a resetting of the selected star to the center of the star finder, a measurement by the navigation device of the horizontal coordinates of the star targeted by the star finder after registration, and a recording of the measurement instant, the adjusted horizontal coordinates thus measured, the associated equatorial coordinates and the geographical coordinates of the viewfinder at the measurement instant.
De préférence, la mesure par le dispositif de navigation des coordonnées horizontales de l’étoile visée par le viseur stellaire après recalage est réalisée pendant 5 secondes.Preferably, the measurement by the navigation device of the horizontal coordinates of the star targeted by the star finder after adjustment is carried out for 5 seconds.
Avantageusement, le dispositif de navigation peut comprendre un système de positionnement par satellites.Advantageously, the navigation device may include a satellite positioning system.
Avantageusement, le dispositif de navigation peut comprendre une centrale inertielle.Advantageously, the navigation device can include an inertial unit.
Naviguer à l’aide des étoiles nécessite de travailler avec trois systèmes de coordonnées différentes : les coordonnées géographiques, les coordonnées équatoriales et les coordonnées horizontales.Navigating using the stars requires working with three different coordinate systems: geographic coordinates, equatorial coordinates and horizontal coordinates.
Comme cela est illustré sur la
Elles permettent de repérer un point sur la surface de la Terre T, la surface de la Terre étant assimilée à un géoïde de centre O. Sur la figure est représenté un trièdre x, y et z se coupant au centre O du géoïde terrestre T. La longitude L est un angle mesuré par rapport au centre O du géoïde. La longitude est un premier angle qui a pour origine le méridien de Greenwich mget est comptée positivement en tournant autour de l’axe z vers l’Est E et négativement vers l’ouest W. La latitude Φ est un second angle mesuré par rapport à l’équateur péet est compté positivement en tournant autour de l’axe y vers le pôle Nord N, négativement vers le pôle Sud S.They make it possible to locate a point on the surface of the Earth T, the surface of the Earth being assimilated to a geoid with center O. The figure shows a trihedron x, y and z intersecting at the center O of the terrestrial geoid T. Longitude L is an angle measured relative to the center O of the geoid. Longitude is a first angle which has its origin on the Greenwich meridian m g and is counted positively by rotating around the z axis towards the East E and negatively towards the west W. The latitude Φ is a second angle measured by relative to the equator p é and is counted positively by rotating around the y axis towards the North Pole N, negatively towards the South Pole S.
Un point M sur le géoïde terrestre T a une longitude L correspondant à l’angle entre le méridien de Greenwich mget le méridien mmde M, et une latitude Φ correspondant à l’angle entre l’équateur péet la parallèle pmà l’équateur passant par M.A point M on the terrestrial geoid T has a longitude L corresponding to the angle between the Greenwich meridian m g and the meridian m m of M, and a latitude Φ corresponding to the angle between the equator p é and the parallel p m to the equator passing through M.
Comme cela est illustré sur la
Cette fois-ci, la Terre T se trouve au centre d’une Sphère céleste Cs sur lequel passe un cercle horaire Ch qui passe par le pôle nord céleste NC et le pôle sud céleste SC. La représentation comprend en outre le cercle de l’équateur céleste Cc, prolongation de l’équateur sur la sphère céleste Cs, et le cercle de l’écliptique Ce, les deux cercles étant sécant au point vernal Pv. L’arc de cercle s’étendant entre l’étoile St et l’équateur céleste Cc le long du cercle horaire Ch correspond à la déclinaison δ. L’arc de cercle s’étendant le long de l’équateur céleste Ce entre le point vernal Pv et le Cercle horaire Ch correspond à l’ascension droite α.This time, Earth T is at the center of a Celestial Sphere Cs over which passes a hour circle Ch which passes through the north celestial pole NC and the south celestial pole SC. The representation also includes the circle of the celestial equator Cc, extension of the equator on the celestial sphere Cs, and the circle of the ecliptic Ce, the two circles intersecting at the vernal point Pv. The arc of a circle extending between the star St and the celestial equator Cc along the hour circle Ch corresponds to the declination δ. The arc of a circle extending along the celestial equator Ce between the vernal point Pv and the hourly circle Ch corresponds to the right ascension α.
Les coordonnées équatoriales permettent ainsi de repérer une étoile St sur la sphère céleste, une extension virtuelle de la Terre à l’infini. L’ascension droite α correspond, en un certain sens, à la longitude L sur cette sphère céleste et a pour origine le point vernal Pv, nom donné à l’un des deux points d’intersection du plan de l’équateur céleste Cc et du plan de l’écliptique Ce (plan contenant la trajectoire apparente du Soleil vue depuis la Terre) sur la sphère céleste. La déclinaison δ correspond, en un certain sens, à la latitude sur la sphère céleste, et a pour origine l’équateur céleste Cc.The equatorial coordinates thus make it possible to locate a star St on the celestial sphere, a virtual extension of the Earth to infinity. The right ascension α corresponds, in a certain sense, to the longitude L on this celestial sphere and has its origin at the vernal point Pv, name given to one of the two points of intersection of the plane of the celestial equator Cc and of the plane of the ecliptic Ce (plane containing the apparent trajectory of the Sun seen from the Earth) on the celestial sphere. The declination δ corresponds, in a certain sense, to the latitude on the celestial sphere, and originates from the celestial equator Cc.
Comme cela est illustré sur la
Elles permettent de mesurer la direction d’un objet dans le ciel, telle qu’une étoile St, depuis un point sur Terre. L’azimut A est défini comme l’angle entre le Nord géographique N et l’objet visé St. Il est compté positivement vers l’Est. La hauteur h est définie comme l’angle entre le plan tangent à la surface de la Terre et l’objet visé St. Elle est comptée positivement vers le haut. On parle aussi de repère géographique local, noté [g], associé à ce système de coordonnées : l’axe x est dirigé vers le Nord, l’axe z vers le zénith du lieu (perpendiculaire au plan tangent au géoïde) et l’axe Y formant un trièdre direct (non visible sur cette figure).They make it possible to measure the direction of an object in the sky, such as a star St, from a point on Earth. Azimuth A is defined as the angle between geographic North N and the target object St. It is counted positively towards the East. The height h is defined as the angle between the plane tangent to the surface of the Earth and the target object St. It is counted positively upwards. We also speak of a local geographical reference, denoted [g], associated with this coordinate system: the x axis is directed towards the North, the z axis towards the zenith of the place (perpendicular to the plane tangent to the geoid) and the Y axis forming a direct trihedron (not visible in this figure).
Le passage d’un système de coordonnées à l’autre s’effectue par des considérations géométriques et temporelles.The transition from one coordinate system to another is carried out by geometric and temporal considerations.
La
On remarque sur la
[Math 1] S=
avec GAST désignant le temps sidéral apparent de Greenwich et LAST désignant le temps sidéral apparent local.with GAST denoting Greenwich apparent sidereal time and LAST denoting local apparent sidereal time.
Les coordonnées équatoriales (α,δ) d’un nombre important d’étoiles sont données par des catalogues (ou éphémérides) dans un repère céleste J2000 par exemple à une date spécifique appelée J2000 : le 1erJanvier 2000 à 12h, Temps Terrestre. Notons que le Temps Terrestre diffère du temps universel coordonnée (UTC pour « coordinated universal time ») utilisé dans la vie de tous les jours par un certain nombre de secondes entières appelées secondes intercalaires, évoluant chaque année. Le nombre d’étoiles renseigné varie suivant le type de catalogue ; il peut atteindre plusieurs millions. Le passage des coordonnées équatoriales d’une étoile à la date J2000 à une autre date et heure d’observation se fait par des corrections physiques du mouvement de la Terre (précession, nutation) et des étoiles elles-mêmes (mouvement propre) qui sont détaillées notamment dans l’article de Meeus, J. (1998) intitulé « Astronomical Algorithms » aux éditions « Willmann-Bell ». Ces corrections physiques ne sont pas détaillées ici.The equatorial coordinates (α, δ) of a large number of stars are given by catalogs (or ephemeris) in a celestial reference frame J2000 for example on a specific date called J2000: January 1, 2000 at 12 p.m., Terrestrial Time. Note that Earth Time differs from Coordinated Universal Time (UTC) used in everyday life by a certain number of whole seconds called leap seconds, changing each year. The number of stars indicated varies depending on the type of catalog; it can reach several millions. The transition from the equatorial coordinates of a star on date J2000 to another date and time of observation is done by physical corrections of the movement of the Earth (precession, nutation) and of the stars themselves (proper movement) which are detailed in particular in the article by Meeus, J. (1998) entitled “Astronomical Algorithms” published by “Willmann-Bell”. These physical corrections are not detailed here.
La base de données VizieR (https://vizier.cds.unistra.fr/viz-bin/VizieR) de l’université de Strasbourg permet d’accéder à un certain nombre de catalogues.The VizieR database (https://vizier.cds.unistra.fr/viz-bin/VizieR) of the University of Strasbourg provides access to a certain number of catalogs.
Le calcul du temps sidéral moyen à Greenwich à 0h Temps Universel (UT), noté GMST_0 dans le système (S), se déduit à l’aide de la formule suivante (une version similaire à cette formule, exprimée en secondes horaires, est donnée dans le document de Meeus de 1998 déjà mentionné) :The calculation of the mean sidereal time at Greenwich at 0h Universal Time (UT), denoted GMST_0 in the system (S), is deduced using the following formula (a version similar to this formula, expressed in hourly seconds, is given in the 1998 Meeus document already mentioned):
Où T est le nombre de siècles Juliens écoulés depuis J2000 (pour une définition des jours Juliens, voir le document de Meeus de 1998 déjà mentionné). De
De GMST, nous déduisons le temps sidéral apparent à Greenwich GAST :From GMST, we deduce the apparent sidereal time at Greenwich GAST:
Où εγdésigne l’équation des équinoxes, dont le document de Meeus de 1998 déjà mentionné donne une méthode de calcul. Disposant des coordonnées équatoriales (α,δ) d’une étoile à observer et du temps sidéral apparent de Greenwich GAST, il reste à déterminer l’azimut et l’élévation (A,h) de l’étoile depuis le lieu d’observation. Pour cela, il est nécessaire d’avoir un moyen de mesure précis, permettant d’obtenir les coordonnées horizontales (A,h).Where ε γ designates the equation of the equinoxes, for which the already mentioned 1998 Meeus document gives a calculation method. Having the equatorial coordinates (α, δ) of a star to be observed and the apparent sidereal Greenwich time GAST, it remains to determine the azimuth and elevation (A, h) of the star from the observation location . For this, it is necessary to have a precise measuring means, making it possible to obtain the horizontal coordinates (A,h).
Le procédé selon l’invention propose alors d’utiliser un viseur 10 deux axes, c’est-à-dire un ensemble de caméras monté sur un système de cardans commandable suivant deux axes θ et ψ, site et gisement, comme cela est illustré sur la
Ce viseur 10 est destiné à être utilisé par exemple sur un navire 12 naviguant en pleine mer 14. Le viseur 10 doit comporter une centrale inertielle permettant d’estimer le nord géographique et la verticale du lieu d’observation ; ainsi, le viseur est pilotable en Azimut (A) et Elévation (h). La centrale inertielle permet aussi la stabilisation inertielle de la Ligne de Visée (c’est-à-dire l’axe optique d’une caméra), afin que l’image restituée par la caméra soit nette quelles que soient les mouvements du navire.This viewfinder 10 is intended to be used for example on a ship 12 sailing on the open sea 14. The viewfinder 10 must include an inertial unit making it possible to estimate the geographic north and the vertical of the observation location; thus, the viewfinder can be controlled in Azimuth (A) and Elevation (h). The inertial unit also allows the inertial stabilization of the Line of Sight (i.e. the optical axis of a camera), so that the image rendered by the camera is clear whatever the movements of the ship.
Toutes les étoiles ne sont cependant pas adéquates pour réaliser une fonction de géolocalisation, et à contrario certaines sont à privilégier.However, not all stars are suitable for carrying out a geolocation function, and on the contrary some are to be preferred.
En effet, il existe des lieux singuliers et des lieux d’insensibilité aux erreurs du viseur en Azimut et Elévation pour le calcul de la localisation. Ces lieux singuliers et d’insensibilité sont intrinsèques au système d’équations [Math 1] et sont à compléter avec d’autres précautions : la visée proche de l’horizon est déconseillée à cause de la réfraction atmosphérique (difficile à estimer), la visée proche du zénith est parfois irréalisable suivant le type de système de visée utilisé, etc.Indeed, there are singular places and places of insensitivity to viewfinder errors in Azimuth and Elevation for the calculation of the location. These singular places of insensitivity are intrinsic to the system of equations [Math 1] and must be completed with other precautions: aiming close to the horizon is not recommended because of atmospheric refraction (difficult to estimate), the aiming close to the zenith is sometimes unachievable depending on the type of aiming system used, etc.
Pour calculer les lieux singuliers propres au problème de localisation par visée stellaire, c’est-à-dire propres au système d’équations [Math 1], il est nécessaire de calculer les différentielles de la longitude L et de la latitude Φ.To calculate the singular locations specific to the problem of location by star sighting, that is to say specific to the system of equations [Math 1], it is necessary to calculate the differentials of longitude L and latitude Φ.
Supposons que les coordonnées équatoriales (α,δ) peuvent être considérées comme constantes par rapport à l’échelle de temps considérée (une journée d’observation). Alors, on montre que les différentielles de la longitude L et de la latitude Φ s’écrivent, pour
Avec l1, l2, p1, p2, ΔL, Δϕ tels que définis ci-dessus dans les équations [Math 3] à [Math 6].With l 1 , l 2 , p 1 , p 2 , ΔL, Δϕ as defined above in equations [Math 3] to [Math 6].
On voit alors apparaître des lieux singuliers ΔL=0 et ΔΦ=0, des lieux d’insensibilité en Azimut l1=0 et p1=0, et des lieux d’insensibilité en élévation l2=0 et p2=0. Ces lieux sont représentés sur le graphique de la
Les courbes en pointillé de la
A contrario, la courbe en trait plein et la courbe en trait mixte décrivent des zones à privilégier pour ce même objectif. La courbe en trait plein représente une zone d’insensibilité de la localisation à une erreur d’élévation h, et la courbe en trait mixte représente une zone d’insensibilité de la localisation à une erreur d’azimut A.Conversely, the solid line curve and the dotted line curve describe areas to be favored for this same objective. The solid line curve represents a zone of insensitivity of the location to an elevation error h, and the dotted line curve represents a zone of insensitivity of the location to an azimuth error A.
La courbe en trait plein d’insensibilité de la localisation à une erreur d’élévation h est définie par l’équation suivante :The solid line curve of insensitivity of location to an elevation error h is defined by the following equation:
La courbe en trait mixte d’insensibilité de la localisation à une erreur d’azimut A est définie par l’équation suivante :The dotted line curve of insensitivity of location to an azimuth error A is defined by the following equation:
En outre, aux points de localisation A=π, on a une double insensibilité de la localisation à la fois à une erreur d’élévation h et à une erreur d’azimut A, comme cela est illustré par les courbes verticales en trait (un trait long, et deux points).Furthermore, at the location points A=π, there is a double insensitivity of the location to both an elevation error h and an azimuth error A, as illustrated by the vertical curves in line (a long line, and two dots).
En outre, sur la représentation graphique, les étoiles en points représentent des exemples d’étoiles insensibles à une erreur de visée en azimut A, les étoiles en trait plein représentent des exemples d’étoiles insensibles à une erreur de visée en élévation h, et les étoiles en pleines (grisées) représentent des exemples d’étoiles exploitables pour la géolocalisation.Furthermore, on the graphical representation, the stars in dots represent examples of stars insensitive to a sighting error in azimuth A, the stars in solid lines represent examples of stars insensitive to a sighting error in elevation h, and the solid stars (grayed) represent examples of stars usable for geolocation.
De plus, il est possible d’obtenir des expressions analytiques de ces courbes et de les exploiter pour discriminer des étoiles en temps réel lors d’observations, et ainsi d’assurer la convergence et l’optimalité de l’estimation de position réalisée.In addition, it is possible to obtain analytical expressions of these curves and exploit them to discriminate stars in real time during observations, and thus ensure the convergence and optimality of the position estimation carried out.
D’un point de vue pratique, l’estimation de la position et de l’attitude du système est réalisée à l’aide de la centrale inertielle : celle-ci est un signal discret, échantillonné à une fréquence fe(généralement de 100 Hz), et donc mis-à-jour régulièrement. La visée d’étoile permet alors d’affiner ou de corriger cette estime.From a practical point of view, the estimation of the position and attitude of the system is carried out using the inertial unit: this is a discrete signal, sampled at a frequency f e (generally 100 Hz), and therefore updated regularly. Star sighting then makes it possible to refine or correct this estimate.
Le logigramme de la
Soit
A l’instant UTC t0, on sauvegarde la position estimée par la centrale inertielle, (L0, ϕ0).At time UTC t 0 , we save the position estimated by the inertial unit, (L 0 , ϕ 0 ).
A partir d’un catalogue d’étoiles et du temps UTC t0, l’algorithme calcule l’ensemble des étoiles observables depuis (L0, ϕ0) à l’aide du système d’équation [Math 1].From a catalog of stars and the UTC time t 0 , the algorithm calculates all the stars observable from (L 0 , ϕ 0 ) using the equation system [Math 1].
L’algorithme choisi alors M étoiles observables (par exemple, M=10 afin de disposer de suffisamment d’étoiles pour un traitement statistique), en éliminant celles à hauteur trop faible h<hinf(par exemple, h < 45° afin d’éviter les problèmes de réfraction atmosphérique) ainsi que celles trop proches des lieux singuliers ΔL=0 et Δϕ=0 et en privilégiant celles proches des lieux d’insensibilité l1=0, p1=0, l2=0 ou p2=0. Ce choix est réalisé en calculant le paramètre [Math 2]
Plus explicitement, la sélection d’étoiles est réalisée à partir d’étoiles recensées dans un catalogue où chaque étoile est associée à des coordonnées équatoriales à une date spécifique appelée J2000, les coordonnées équatoriales comportant une ascension droite, α, et une déclinaison, δ.More explicitly, the selection of stars is carried out from stars listed in a catalog where each star is associated with equatorial coordinates on a specific date called J2000, the equatorial coordinates comprising a right ascension, α, and a declination, δ .
Le procédé de sélection d’étoiles comprend tout d’abord, dans une première étape 700, un relevé des coordonnées géographiques d’un viseur stellaire à un instant de départ par le dispositif de navigation, les coordonnées géographiques du viseur stellaire comprenant une longitude, L, et une latitude, ϕ.The star selection method firstly comprises, in a first step 700, a reading of the geographical coordinates of a star finder at a starting time by the navigation device, the geographical coordinates of the star finder including a longitude, L, and a latitude, ϕ.
Puis, dans une deuxième étape 710, on réalise, pour chaque étoile du catalogue, un traitement de fond. Le traitement de fond pour chaque étoile comprend les sous-étapes suivantes.Then, in a second step 710, a background treatment is carried out for each star in the catalog. The background processing for each star includes the following substeps.
Un calcul 712 des coordonnées de visée de l’étoile comportant un calcul, à partir du système d’équations [Math 1], des coordonnées horizontales de l’étoile à partir des coordonnées équatoriales de l’étoile et de la position géographique du viseur, les coordonnées équatoriales comprenant un azimut, A, de l’étoile et une élévation, h, de l’étoile.A calculation 712 of the sighting coordinates of the star comprising a calculation, from the system of equations [Math 1], of the horizontal coordinates of the star from the equatorial coordinates of the star and the geographical position of the viewfinder , the equatorial coordinates including an azimuth, A, of the star and an elevation, h, of the star.
Une détermination de l’accessibilité de l’étoile à partir des coordonnées de visée calculées, la détermination d’accessibilité comportant une comparaison 714 de l’élévation calculée pour l’étoile à un seuil d’élévation minimale, l’étoile étant inaccessible si son élévation calculée est inférieure au seuil d’élévation minimale.A determination of the accessibility of the star from the calculated sighting coordinates, the determination of accessibility comprising a comparison 714 of the elevation calculated for the star to a minimum elevation threshold, the star being inaccessible if its calculated elevation is less than the minimum elevation threshold.
Et, si l’étoile est accessible, un calcul 716 d’un paramètre de sensibilité d’observation de l’étoile, les coordonnées de visée de l’étoile correspondant aux coordonnées horizontales de l’étoile par rapport à la position géographique du viseur stellaire, les coordonnées horizontales comportant un azimut, A, et une élévation, h.And, if the star is accessible, a calculation 716 of an observation sensitivity parameter of the star, the sighting coordinates of the star corresponding to the horizontal coordinates of the star relative to the geographical position of the viewfinder stellar, the horizontal coordinates including an azimuth, A, and an elevation, h.
Puis, dans une troisième étape 720, on réalise un stockage des coordonnées de visée et du paramètre de sensibilité d’observation calculés pour chaque étoile, et un classement des étoiles ainsi stockée par ordre croissant de paramètre de sensibilité de l’étoile, les étoiles dont le paramètre de sensibilité est nul étant exclue du classement.Then, in a third step 720, the aiming coordinates and the observation sensitivity parameter calculated for each star are stored, and the stars thus stored are classified in increasing order of the sensitivity parameter of the star, the stars whose sensitivity parameter is zero being excluded from the classification.
Enfin dans une quatrième étape 730, on réalise une sélection d’un nombre déterminé d’étoiles observables, les étoiles étant sélectionnées selon l’ordre de classement des étoiles stockées.Finally, in a fourth step 730, a selection of a determined number of observable stars is made, the stars being selected according to the order of classification of the stored stars.
Pour réaliser une géolocalisation fiable, par exemple sur un navire en pleine mer, et donc pour s’affranchir d’une dérive du dispositif de navigation, l’invention propose de mettre en œuvre le procédé de correction de dérive illustré sur le logigramme de la
Le procédé de correction de dérive d’un dispositif de navigation d’un véhicule à l’aide d’un moyen de géolocalisation par visée stellaire présenté sur le logigramme de la
Puis, le procédé de correction comprend une étape suivante 810 de mesures comprenant, tout d’abord, un relevé des coordonnées géographiques du viseur stellaire par le dispositif de navigation, les coordonnées géographiques du viseur stellaire comprenant une longitude L et une latitude ϕ, puis, pour chacune des étoiles sélectionnées, une vérification de l’accessibilité de l’étoile et une mesure des coordonnées horizontales de l’étoile, les coordonnées horizontales mesurées étant enregistrées en association avec les coordonnées équatoriales correspondantes pour l’étoile et l’instant de mesure.Then, the correction method comprises a following step 810 of measurements comprising, first of all, a reading of the geographical coordinates of the star finder by the navigation device, the geographical coordinates of the star finder comprising a longitude L and a latitude ϕ, then , for each of the selected stars, a check of the accessibility of the star and a measurement of the horizontal coordinates of the star, the measured horizontal coordinates being recorded in association with the corresponding equatorial coordinates for the star and the instant of measure.
L’étape de mesures 810 comprend un calcul 811 des coordonnées horizontales de l’étoile sélectionnée à partir de l’instant de calcul, des coordonnées géographiques du viseur stellaire à l’instant de calcul, et des coordonnées équatoriales indiquées dans un catalogue où chaque étoile est associée à des coordonnées équatoriales à une date spécifique appelée J2000, les coordonnées équatoriales comportant une ascension droite α et une déclinaison δ et les coordonnées horizontales comportant un azimut A et une élévation h. Elle comprend ensuite une orientation 812 d’un viseur stellaire selon les coordonnées horizontales calculées, et une vérification 813 de l’accessibilité de l’étoile sélectionnée.The measurement step 810 includes a calculation 811 of the horizontal coordinates of the star selected from the calculation instant, of the geographical coordinates of the star finder at the calculation instant, and of the equatorial coordinates indicated in a catalog where each star is associated with equatorial coordinates on a specific date called J2000, the equatorial coordinates having a right ascension α and a declination δ and the horizontal coordinates having an azimuth A and an elevation h. It then includes an orientation 812 of a star finder according to the calculated horizontal coordinates, and a verification 813 of the accessibility of the selected star.
Si l’étoile sélectionnée est inaccessible, elle comprend une nouvelle réalisation des trois étapes précédentes avec une nouvelle étoile de la sélection d’étoiles, alors que si l’étoile est accessible, elle comprend un recalage 814 de l’étoile sélectionnée au centre du viseur stellaire, une mesure 815 pendant 5 secondes par le dispositif de navigation des coordonnées horizontales de l’étoile visée par le viseur stellaire après recalage, et un enregistrement 816 de l’instant de mesure, des coordonnées horizontales recalées ainsi mesurées, des coordonnées équatoriales associées et des coordonnées géographiques du viseur à l’instant de mesure.If the selected star is inaccessible, it includes a new realization of the three previous steps with a new star from the selection of stars, whereas if the star is accessible, it includes an adjustment 814 of the selected star to the center of the star finder, a measurement 815 for 5 seconds by the navigation device of the horizontal coordinates of the star targeted by the star finder after adjustment, and a recording 816 of the instant of measurement, of the readjusted horizontal coordinates thus measured, of the equatorial coordinates associated and geographical coordinates of the viewfinder at the time of measurement.
Puis, le procédé de correction comprend une étape suivante 820, de traitement comportant, pour toutes les étoiles pour lesquelles des coordonnées horizontales mesurées ont été enregistrées, un calcul 822 des coordonnées géographiques du viseur stellaire à partir des coordonnées horizontales mesurées enregistrées pour l’étoile, de l’instant de mesure enregistré, des coordonnées équatoriales enregistrées et du système d’équations [Math 1], et une estimation 824 de l’erreur entre les coordonnées géographiques calculées et les coordonnées géographiques relevées par le dispositif de navigation lors de la visée de l’étoile.Then, the correction method comprises a following processing step 820 comprising, for all the stars for which measured horizontal coordinates have been recorded, a calculation 822 of the geographical coordinates of the star finder from the measured horizontal coordinates recorded for the star , the recorded measurement instant, the recorded equatorial coordinates and the system of equations [Math 1], and an estimate 824 of the error between the calculated geographic coordinates and the geographic coordinates recorded by the navigation device during the sight of the star.
Puis, il comprend une étape suivante 830 d’estimation de l’erreur de dérive du dispositif de navigation à partir des erreurs estimées sur les coordonnées géographiques.Then, it includes a following step 830 of estimating the drift error of the navigation device from the errors estimated on the geographic coordinates.
Enfin, il comprend une étape 840 de correction de la dérive du dispositif de navigation à partir de l’erreur de dérive estimée.Finally, it includes a step 840 of correcting the drift of the navigation device based on the estimated drift error.
En d’autres termes, si on prend l’étoile choisie parmi les N étoiles présélectionnées, t le temps UTC échantillonné à fe, et (Li(t), ϕi(t)) l’estimation de position renvoyée par la centrale inertielle et/ou le système de navigation par satellite du dispositif de navigation, i.e. un signal échantillonné à fe.In other words, if we take the star chosen from the N preselected stars, t the UTC time sampled at f e , and (L i (t), ϕ i (t)) the position estimate returned by the inertial unit and/or the satellite navigation system of the navigation device, ie a signal sampled at f e .
Dynamiquement, c’est-à-dire à chaque période d’échantillonnage, le procédé calcule précisément (Ai(t), hi(t)) de l’étoile Stià partir du catalogue donnant les coordonnées équatoriales (αi(t), δi(t)) de Stià l’instant UTC t, de (Li(t), ϕi(t)) et du système d’équations [Math 1].Dynamically, that is to say at each sampling period, the process precisely calculates (A i (t), h i (t)) of the star St i from the catalog giving the equatorial coordinates (α i (t), δ i (t)) of St i at UTC time t, of (L i (t), ϕ i (t)) and the system of equations [Math 1].
Si (Ai(t), hi(t)) est trop proche d’une singularité, alors l’étoile
Sinon, on pointe le viseur dans la direction (Ai(t), hi(t)) et, si l’étoile n’est pas visible, alors celle-ci est ignorée et on passe à l’étoile suivante. Sinon, on mesure
On sauvegarde les signaux
A la suite de la séquence de mesure, pour chaque étoile mesurée
On calcule alors l’erreur entre la position retournée par le dispositif de navigation et la position estimée par visée stellaire :We then calculate the error between the position returned by the navigation device and the position estimated by star sighting:
Cette erreur représente une estime de l’erreur de position de la centrale. Celle-ci doit avoir une dynamique lente, et donc être de pente nulle sur l’étendue de la mesure (5 secondes). Le calcul de la moyenne de ce signal renvoie une estime de l’erreur de position de la centrale pour la visée stellaire numéro i.This error represents an estimate of the position error of the unit. This must have a slow dynamic, and therefore be of zero slope over the extent of the measurement (5 seconds). Calculating the average of this signal returns an estimate of the position error of the center for star sight number i.
On sauvegarde les moyennes
Claims (10)
le procédé comprenant :
- un relevé des coordonnées géographiques d’un viseur stellaire à un instant de départ (t0) par le dispositif de navigation, les coordonnées géographiques du viseur stellaire comprenant une longitude (L) et une latitude (ϕ), puis
- pour chaque étoile du catalogue, un calcul des coordonnées de visée de l’étoile, une détermination de l’accessibilité de l’étoile à partir des coordonnées de visée calculées, et, si l’étoile est accessible, un calcul d’un paramètre (X) de sensibilité d’observation de l’étoile, les coordonnées de visée de l’étoile correspondant aux coordonnées horizontales de l’étoile par rapport à la position géographique du viseur stellaire, les coordonnées horizontales comportant un azimut (A) et une élévation (h), puis
- un stockage des coordonnées de visée et du paramètre de sensibilité d’observation calculés pour chaque étoile, et un classement des étoiles ainsi stockée par ordre croissant de paramètre de sensibilité de l’étoile, les étoiles dont le paramètre (X) de sensibilité est nul étant exclue du classement, et
- une sélection d’un nombre déterminé d’étoiles observables, les étoiles étant sélectionnées selon l’ordre de classement des étoiles stockées.Method for selecting observable stars to be used for geolocation by a vehicle, in particular a ship at sea, from a means of geolocation by star sighting associated with at least one navigation device, the selection of stars being carried out from stars listed in a catalog where each star is associated with equatorial coordinates on a specific date, the equatorial coordinates comprising a right ascension (α) and a declination (δ),
the process comprising:
- a reading of the geographical coordinates of a star finder at a starting time (t 0 ) by the navigation device, the geographical coordinates of the star finder including a longitude (L) and a latitude (ϕ), then
- for each star in the catalog, a calculation of the sighting coordinates of the star, a determination of the accessibility of the star from the calculated sighting coordinates, and, if the star is accessible, a calculation of a parameter (X) of observation sensitivity of the star, the sighting coordinates of the star corresponding to the horizontal coordinates of the star relative to the geographical position of the star finder, the horizontal coordinates comprising an azimuth (A) and an elevation (h), then
- storage of the sighting coordinates and the observation sensitivity parameter calculated for each star, and a classification of the stars thus stored in ascending order of sensitivity parameter of the star, the stars whose sensitivity parameter (X) is null being excluded from the classification, and
- a selection of a determined number of observable stars, the stars being selected according to the order of classification of the stored stars.
dans lequel la détermination de l’accessibilité de l’étoile comprend une comparaison de l’élévation calculée pour l’étoile à un seuil d’élévation minimale, l’étoile étant inaccessible si son élévation calculée est inférieure au seuil d’élévation minimale.Method for selecting observable stars according to claim 1, in which the calculation of the sighting coordinates of each star comprises a calculation of the horizontal coordinates of the star from the equatorial coordinates of the star and the geographical position of the viewfinder, the equatorial coordinates including an azimuth (A) of the star and an elevation (h) of the star, and
wherein determining the accessibility of the star includes comparing the calculated elevation for the star to a minimum elevation threshold, the star being inaccessible if its calculated elevation is less than the minimum elevation threshold.
[Math 1]
avec GAST désignant le temps sidéral apparent de Greenwich et LAST désignant le temps sidéral apparent local.Method for selecting observable stars according to one of claims 2 or 3, in which the calculation of the sighting coordinates of the star from the geographical coordinates of the star finder and the equatorial coordinates of the star considered is carried out from of the system of equations (S):
[Math 1]
with GAST denoting Greenwich apparent sidereal time and LAST denoting local apparent sidereal time.
[Math 2]
avec [Math 3]
[Math 4]
[Math 5]
[Math 6]
avec
[Math 2]
with [Math 3]
[Math 4]
[Math 5]
[Math 6]
with
- une étape de pré-sélection d’étoiles comportant une sélection d’étoiles à utiliser par le dispositif de géolocalisation par visée stellaire à partir du procédé de sélection d’étoiles observables à utiliser selon l’une des revendications 1 à 6, puis
- une étape de mesures comprenant, tout d’abord, un relevé des coordonnées géographiques du viseur stellaire par le dispositif de navigation, les coordonnées géographiques du viseur stellaire comprenant une longitude (L) et une latitude (ϕ), puis, pour chacune des étoiles sélectionnées, une vérification de l’accessibilité de l’étoile et une mesure des coordonnées horizontales de l’étoile, les coordonnées horizontales mesurées étant enregistrées en association avec les coordonnées équatoriales correspondantes pour l’étoile et l’instant de mesure, puis
- une étape de traitement comportant, pour toutes les étoiles pour lesquelles des coordonnées horizontales mesurées ont été enregistrées :
- un calcul des coordonnées géographiques du viseur stellaire à partir des coordonnées horizontales mesurées enregistrées pour l’étoile, de l’instant de mesure enregistré, des coordonnées équatoriales enregistrées et du système d’équations (S) :
[Math 1]
avec GAST désignant le temps sidéral apparent de Greenwich et LAST désignant le temps sidéral apparent local. - une estimation de l’erreur entre les coordonnées géographiques calculées et les coordonnées géographiques relevées par le dispositif de navigation lors de la visée de l’étoile, puis
- une étape de correction de la dérive du dispositif de navigation à partir de l’erreur de dérive estimée.Method for correcting drift of a navigation device of a vehicle using a means of geolocation by star sighting, the correction method comprising:
- a star pre-selection step comprising a selection of stars to be used by the geolocation device by stellar sighting from the method of selecting observable stars to be used according to one of claims 1 to 6, then
- a measurement step comprising, first of all, a reading of the geographical coordinates of the star finder by the navigation device, the geographical coordinates of the star finder including a longitude (L) and a latitude (ϕ), then, for each of the selected stars, a check of the accessibility of the star and a measurement of the horizontal coordinates of the star, the measured horizontal coordinates being recorded in association with the corresponding equatorial coordinates for the star and the time of measurement, then
- a processing step comprising, for all the stars for which measured horizontal coordinates have been recorded:
- a calculation of the geographical coordinates of the star finder from the measured horizontal coordinates recorded for the star, the recorded measurement time, the recorded equatorial coordinates and the system of equations (S):
[Math 1]
with GAST denoting Greenwich apparent sidereal time and LAST denoting local apparent sidereal time. - an estimate of the error between the calculated geographic coordinates and the geographic coordinates recorded by the navigation device when aiming at the star, then
- a step of correcting the drift of the navigation device based on the estimated drift error.
- un calcul des coordonnées horizontales de l’étoile sélectionnée à partir de l’instant de calcul, des coordonnées géographiques du viseur stellaire à l’instant de calcul, et des coordonnées équatoriales indiquées dans un catalogue où chaque étoile est associée à des coordonnées équatoriales à une date spécifique, les coordonnées équatoriales comportant une ascension droite (α) et une déclinaison (δ) et les coordonnées horizontales comportant un azimut (A) et une élévation (h), puis
- une orientation d’un viseur stellaire selon les coordonnées horizontales calculées,
- une vérification de l’accessibilité de l’étoile sélectionnée,
- si l’étoile sélectionnée est inaccessible, une nouvelle réalisation des trois étapes précédentes avec une nouvelle étoile de la sélection d’étoiles,
- si l’étoile est accessible, un recalage de l’étoile sélectionnée au centre du viseur stellaire, une mesure par le dispositif de navigation des coordonnées horizontales de l’étoile visée par le viseur stellaire après recalage, et un enregistrement de l’instant de mesure, des coordonnées horizontales recalées ainsi mesurées, des coordonnées équatoriales associées et des coordonnées géographiques du viseur à l’instant de mesure.Method for correcting drift of a vehicle navigation device according to claim 7, in which the measurement step comprises:
- a calculation of the horizontal coordinates of the selected star from the calculation instant, the geographical coordinates of the star finder at the calculation instant, and the equatorial coordinates indicated in a catalog where each star is associated with equatorial coordinates on a specific date, the equatorial coordinates comprising a right ascension (α) and a declination (δ) and the horizontal coordinates comprising an azimuth (A) and an elevation (h), then
- an orientation of a star finder according to the calculated horizontal coordinates,
- a check of the accessibility of the selected star,
- if the selected star is inaccessible, a new execution of the three previous steps with a new star from the star selection,
- if the star is accessible, a resetting of the selected star to the center of the star finder, a measurement by the navigation device of the horizontal coordinates of the star targeted by the star finder after resetting, and a recording of the instant measurement, the reset horizontal coordinates thus measured, the associated equatorial coordinates and the geographical coordinates of the viewfinder at the measurement instant.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2210995A FR3141254A1 (en) | 2022-10-24 | 2022-10-24 | Detection and correction of drift of a navigation device by star sighting |
PCT/FR2023/051667 WO2024089353A2 (en) | 2022-10-24 | 2023-10-24 | Detection and correction of drift in a star tracking navigation device |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR2210995 | 2022-10-24 | ||
FR2210995A FR3141254A1 (en) | 2022-10-24 | 2022-10-24 | Detection and correction of drift of a navigation device by star sighting |
Publications (1)
Publication Number | Publication Date |
---|---|
FR3141254A1 true FR3141254A1 (en) | 2024-04-26 |
Family
ID=85685649
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
FR2210995A Pending FR3141254A1 (en) | 2022-10-24 | 2022-10-24 | Detection and correction of drift of a navigation device by star sighting |
Country Status (2)
Country | Link |
---|---|
FR (1) | FR3141254A1 (en) |
WO (1) | WO2024089353A2 (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080269965A1 (en) * | 2006-12-28 | 2008-10-30 | Luk-Paszyc Jerzy W | Calibration of ship attitude reference |
CN114001756A (en) * | 2021-11-03 | 2022-02-01 | 中国科学院光电技术研究所 | Outfield ground star finding method of small-view-field star sensor |
-
2022
- 2022-10-24 FR FR2210995A patent/FR3141254A1/en active Pending
-
2023
- 2023-10-24 WO PCT/FR2023/051667 patent/WO2024089353A2/en unknown
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080269965A1 (en) * | 2006-12-28 | 2008-10-30 | Luk-Paszyc Jerzy W | Calibration of ship attitude reference |
CN114001756A (en) * | 2021-11-03 | 2022-02-01 | 中国科学院光电技术研究所 | Outfield ground star finding method of small-view-field star sensor |
Non-Patent Citations (2)
Title |
---|
J.E. MCINTYRE ET AL: "A star scan/attitude determination experiment conducted on the geostationary meteorological satellite", ACTA ASTRONAUTICA, vol. 7, no. 2, 1 February 1980 (1980-02-01), GB, pages 137 - 154, XP055538366, ISSN: 0094-5765, DOI: 10.1016/0094-5765(80)90056-9 * |
MEEUS, J.: "Astronomical Al-gorithms", 1998 |
Also Published As
Publication number | Publication date |
---|---|
WO2024089353A2 (en) | 2024-05-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP2353024B1 (en) | Method for geolocating an object by multitelemetry | |
US7447591B2 (en) | Daytime stellar imager for attitude determination | |
US7349803B2 (en) | Daytime stellar imager | |
US7349804B2 (en) | Daytime stellar imager | |
US9791278B2 (en) | Navigating with star tracking sensors | |
US20110004405A1 (en) | Earth horizon sensor | |
US10621741B2 (en) | System and method for a digital image sextant | |
EP2490042B1 (en) | Method and system for determining navigation parameters of an aircraft | |
FR3018383A1 (en) | METHOD AND DEVICE FOR DETERMINING NAVIGATION PARAMETERS OF AN AIRCRAFT DURING A LANDING PHASE | |
US8767072B1 (en) | Geoposition determination by starlight refraction measurement | |
CN102575960A (en) | Apparatus and method for navigation | |
EP2385346B1 (en) | Method for creating a navigation phase in a navigation system involving terrain correlation | |
EP2478334B1 (en) | Three-dimensional location of target land area by merging images captured by two satellite-based sensors | |
WO2014199073A1 (en) | Electronic sextant with a triaxial inertial system and method for position determination | |
WO2014090878A1 (en) | Method for accurately geolocating an image sensor installed on board an aircraft | |
EP1281093B1 (en) | Method and device for instantaneous determination of orientation, based on satellite positioning signals | |
Zhan et al. | High‐accuracy absolute positioning for the stationary planetary rover by integrating the star sensor and inclinometer | |
US11073397B2 (en) | Magnetic-inertial global positioning system | |
EP4272010A1 (en) | Method for determining, using an optronic system, positions in a scene, and associated optronic system | |
FR3141254A1 (en) | Detection and correction of drift of a navigation device by star sighting | |
WO2011039457A1 (en) | Method and device for dynamically locating a moving body | |
US20190104492A1 (en) | Cell Phone-Based Land Navigation Methods and Systems | |
FR3071624B1 (en) | DISPLAY SYSTEM, DISPLAY METHOD, AND COMPUTER PROGRAM | |
EP3980720B1 (en) | Method and device for resetting an inertial unit of a transport means on the basis of information delivered by a viewfinder of the transport means | |
EP3374736B1 (en) | Method for designing a navigation path and method for orienting a sighting member from said navigation path |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PLFP | Fee payment |
Year of fee payment: 2 |
|
PLSC | Publication of the preliminary search report |
Effective date: 20240426 |