FR 2897163 A1 20070810 FR 0601125 A 20060208 PROCEDE DE GEO-LOCALISATION D'UNE OU PLUSIEURS CIBLES Le domaine de l'invention est celui de la géo-localisation d'une ou plusieurs cibles immobiles et plus généralement du géo-référencement d'une zone au moyen d'un capteur optronique passif monté sur un aéronef. On entend par capteur optronique passif un capteur sans émetteur-récepteur optique tel qu'un télémètre laser. On rappelle en relation avec la figure 1, que la géo-localisation d'une cible immobile consiste à déterminer la position absolue de la cible au sol représentée par un point P, à partir des écarts d'abscisse X et d'ordonnée Y entre la cible et la projection au sol Co de l'aéronef aussi désigné porteur ; la position Co est fournie par le système de navigation de l'aéronef, en général avec une erreur. Les angles R et G sont respectivement les angles de roulis et de gisement de la ligne de visée de la cible dans un référentiel terrestre. Les angles G et R dépendent notamment de l'attitude de l'aéronef et de l'orientation de la ligne de visée par rapport à celle de l'aéronef ; ils sont fournis par le système de navigation de l'aéronef. La figure 2 donne une coupe verticale de la figure 1, lorsque Y est nul ; sauf mention expresse, la suite du texte est basée sur ce cas, sans nuire à la généralité ; en effet, l'angle G est généralement petit, dans des conditions opérationnelles habituelles. On a alors : X = h tg R (1) h qui est la hauteur de l'aéronef à l'aplomb de Co, est fournie par le système de navigation. Le résultat obtenu pour X est sensible aux erreurs sur h et R 25 surtout pour des incidences rasantes c'est-à-dire lorsque R est voisin de 90 . On a en effet au premier ordre : Ex = (ôX / 3h) Eh + (3X / 3R) ER Ex/X = Eh/h + ER (1 + tg2 R) / tg R, Pour des incidences rasantes, on a 1/tg R tg R et on obtient 30 finalement : Ex/X = Eh/h + ER tg R, Ex, Eh, ER étant respectivement les erreurs sur X, h et R. Ainsi par exemple, pour h = 1000 m, Eh = 0 et X = 30 km, une erreur de 1 mrad sur R conduit à une erreur de géo-localisation de 900 m, c'est-à-dire de 3% ce qui est inacceptable d'un point de vue opérationnel. En effet, dans le cas d'une application à un système de reconnaissance, la précision de la géo-localisation est une spécification fondamentale de ce système. Dans le cas d'une mission préparée, le système doit acquérir des données dans une zone du théâtre d'opération définie par ses coordonnées absolues et désignée zone cible ; une erreur de géo-localisation conduit à augmenter la taille de la zone à couvrir pour être sûr qu'elle contienne la zone cible. De plus, la connaissance de la distance entre le capteur et les points à imager est nécessaire pour garantir la qualité de l'image (flou de bougé). En outre, une précision sur R à 1 mrad près est très difficile à atteindre en particulier si le capteur est monté sur une nacelle ( pod en anglais) qui elle-même bouge par rapport à l'aéronef : dans ce cas, une erreur de plusieurs mrad est usuelle. De plus, lorsque le terrain est accidenté et que la différence de hauteur entre la cible et Co n'est pas connue, la formule (1) conduit à une erreur sur la mesure de X encore plus importante comme illustré sur la figure 3. On obtient X mesuré au lieu de X. En outre, la réfraction atmosphérique introduit une erreur supplémentaire sur la mesure de X, comme illustré sur la figure 4. On a X mesuré au lieu de X. Cette erreur peut être de l'ordre de 900 mètres. Le but de l'invention est de déterminer au moyen d'un capteur optronique passif, l'abscisse X avec une meilleure précision, de l'ordre de 1% par exemple, notamment lorsque la différence de hauteur entre la cible et Co n'est pas connue et/ou en présence d'une éventuelle réfraction atmosphérique. Pour atteindre ce but, l'invention propose un procédé de géolocalisation d'au moins une cible P immobile à partir d'un aéronef au moyen d'un capteur optronique passif. Il est principalement caractérisé en ce qu'il comprend les étapes suivantes : acquisition par le capteur d'au moins une image 11 contenant la 35 cible P à partir d'une position C1 de l'aéronef et d'une image 12 contenant la cible à partir d'une position C2 de l'aéronef, ces images 11 et 12 ayant une zone de recouvrement, identification dans ladite zone de recouvrement d'au moins une cible P commune aux deux images 11 et 12 et détermination de la position de chaque cible P dans chacune des deux images, calcul de la distance d entre chaque cible P et un point C, de coordonnées connues par rapport à C1 et C2, en fonction de l'angle (31 entre une direction de référence et la ligne de visée de l'image 11, l'angle (32 entre la même direction de référence et la ligne de visée de l'image 12, de la position de chaque cible P dans l'image 11 et dans l'image 12, P étant défini dans un référentiel terrestre par son abscisse X et son ordonnée Y par rapport à la projection au sol de l'aéronef C, calcul pour chaque cible de X et Y en fonction des distances d. La quantité X ainsi obtenue est indépendante de la hauteur h, ce qui permet de s'affranchir des erreurs sur h et de la méconnaissance de la différence de hauteur entre la cible et Co. Elle est aussi quasiment indépendante des effets de réfraction atmosphérique, contrairement à celle obtenue par le calcul décrit dans le préambule. Selon une caractéristique de l'invention, la position de chaque cible P dans les images 11 et 12 est calculée en nombre de pixels et convertie respectivement en un angle E1 et un angle a. Avantageusement, il comprend une étape d'acquisition de l'angle de roulis R et de l'angle de gisement G de la ligne de visée en C dans un référentiel terrestre et la position de chaque cible P est calculée aussi en fonction de ces angles de roulis R et de gisement G. Selon une autre caractéristique de l'invention, la distance d est sensiblement égale à (B cos G)/(r32- (31 + E2-E1), B étant la distance entre les 30 points Cl et C2. Les angles de roulis R et de gisement G étant déterminés, X est de préférence sensiblement égal à d cos G.sin R et Y sensiblement égal à d sin G. Pour des incidences rasantes et notamment lorsque G est petit (cos G 35 proche de 1), on a tg R 1 et donc l'erreur EX/X est quasiment égale à sd/d, ce qui permet de relâcher l'exigence de précision sur l'angle de roulis R. Or, on obtient une erreur sur d d'environ 1% et donc une erreur sur X d'environ 1% également. Selon une autre caractéristique de l'invention, la géo-localisation est 5 effectuée à bord de l'aéronef. L'invention concerne également un procédé de géo-référencement d'une zone, caractérisé en ce qu'il comprend la répétition du procédé de géolocalisation pour différentes cibles de la zone. L'invention a aussi pour objet un coffret de traitement comportant au 10 moins une carte électronique, caractérisé en ce qu'au moins une carte électronique comprend des moyens de mise en oeuvre du procédé précédemment décrit, ainsi qu'une nacelle de reconnaissance destinée à être installée sur un aéronef, caractérisée en ce qu'elle comprend un tel coffret de traitement. 15 D'autres caractéristiques et avantages de l'invention apparaîtront à la lecture de la description détaillée qui suit, faite à titre d'exemple non limitatif et en référence aux dessins annexés dans lesquels : la figure 1 déjà décrite présente la géométrie dans l'espace du 20 problème posé de la géo-localisation d'un point P, la figure 2 déjà décrite illustre schématiquement le calcul de la différence d'abscisse X entre une cible P et CO selon l'état de la technique, la figure 3 déjà décrite illustre schématiquement le calcul de la différence d'abscisse X mesuré entre une cible P et CO selon l'état de la 25 technique, lorsque le terrain est accidenté, la figure 4 déjà décrite illustre l'erreur supplémentaire induite par la réfraction atmosphérique, pour une géolocalisation selon l'état de la technique, la figure 5 illustre schématiquement un exemple de calcul d'une 30 distance d par triangulation selon l'invention, dans le cas particulier d'un angle de gisement nul, la figure 6 illustre schématiquement une étape de calcul d'un angle El à partir de la position d'une cible P dans une image 11, la figure 7 illustre schématiquement le calcul de la distance d'une cible en présence de réfraction atmosphérique, dans les 3 dimensions d'un référentiel terrestre, selon l'invention, la figure 8 illustre schématiquement un exemple de calcul d'une 5 distance d par triangulation selon l'invention, dans le cas d'un angle de gisement G non nul, la figure 9 représente schématiquement une nacelle de reconnaissance selon l'invention. 10 Selon l'invention, X est déterminé par la formule : X = d sin R. On a alors : Ex= (3X / 3d) Ed + (3X / 3R) ER 3X/3d=sinR et 3X/3R=dcosR Ex /X= Ed / d + ER / tg R. Pour des incidences rasantes, on a tg R 1 et donc : Ex/X = Ed/d. 15 Dans ce cas, l'erreur sur X provient alors seulement de l'erreur sur d ; elle est quasiment indépendante de l'erreur sur R. Pour une erreur relative sur d d'environ 1%, on obtient alors une erreur sur X de seulement 300 m, pour une cible située à 30 km. Selon l'invention, la distance d est obtenue par triangulation 20 comme illustré sur la figure 5, dans le cas d'un angle de gisement nul (G=0). Le plan de la figure est celui des lignes de visées de l'aéronef qui se déplace entre une position C, et une position C2. Durant le processus d'acquisition, l'aéronef est successivement localisé en C, et C2 et le capteur monté sur l'aéronef acquiert respectivement 25 des images Il et 12 présentant une zone de recouvrement. Pour simplifier, on assimile la position du capteur à celle de l'aéronef. Les images acquises dans la zone d'acquisition ne sont pas nécessairement des images successives. L'image 11 est par exemple acquise lors d'un balayage aller de la zone d'acquisition et l'image 12 lors d'un balayage retour. Le capteur voit 30 donc chaque cible P de cette zone deux fois, c'est-à-dire une fois par image. La distance d entre cette cible P et un point C dont les coordonnées sont connues par rapport à ci et C2, est donnée approximativement par la formule : d=B/(a2ùa1) B est la distance entre C1 et C2 mesurée par le système de navigation, al l'angle entre une direction de référence et C1P, a2 l'angle entre la même direction de référence et C2P. La direction de référence est par exemple la direction du Nord, ou encore celle du segment C1C2. Sur la figure C est au milieu de C1C2 mais il est situé plus généralement dans le voisinage de C1 et C2. on a: a2ùa1=(32-(31 +c2-c1 ces quantités étant affectées d'erreurs de mesure. 131 mesuré par un système inertiel embarqué tel que par exemple un gyromètre est l'angle de la ligne de visée de l'image I1, c'est-à-dire l'angle entre une direction de référence et C1M1, où M1 est le centre de l'image I1, idem pour 132 avec M2 et 12. El est l'angle entre C1M1 et C1P, déterminé à partir de la position de P par rapport au centre de h mesurée en nombre de pixels, idem pour E2 avec 15 C2M2 et 12. On a représenté figure 6 une image I. La cible P est repérée dans cette image par ses coordonnées établies en nombre de pixels par rapport au centre M1. Comme la largeur de l'image mesurée en pixels correspond à un angle déterminé, vu de l'avion, l'angle c1 est directement 20 déduit des coordonnées de P : et est obtenu par conversion d'un nombre de pixels en un angle. Lorsque par exemple la largeur de l'image comprend 1000 pixels et correspond à un angle de 1 , l'angle ci pour une cible P située à 400 pixels de M1 est de 0.4 . 25 De même, les valeurs de roulis R associées aux différentes cibles P sont obtenues, à l'instar des angles a, en faisant la somme des roulis de la ligne de visée (équivalents aux 13) et d'angles (équivalents aux c) obtenus par conversion des positions en roulis des cibles P dans l'image, positions comptées en pixels. 30 L'erreur obtenue pour d dépend de : l'erreur sur B mesurée par le système de navigation et qui est typiquement inférieure à 0.3 %, l'erreur sur a2 ù al dit angle relatif dont la précision est 35 principalement limitée par la dérive du système inertiel de l'ordre de 1 /heure et par la marche aléatoire de l'ordre de 0.05 /heurel12 ; cette erreur est inférieure à 100 Arad en quelques secondes. Pour une distance B de 1000 mètres et une distance d de 30 kilomètres, l'erreur relative introduite par ce terme est de l'ordre de 0.3 %. On obtient finalement une erreur sur d inférieure à 10/0 et donc une erreur sur X également inférieure à 1%, dans les conditions défavorables décrites précédemment (incidence rasante). Avec ce calcul, la distance X est obtenue sans qu'il soit nécessaire d'acquérir la hauteur h, ce qui présente aussi un avantage. Cependant, la mesure de hauteur peut être utilisée à condition qu'on puisse lui accorder une confiance suffisante, par exemple si elle est assez précise et si la différence de hauteur entre la cible et Co est connue. Ce calcul est de préférence réitéré pour plusieurs cibles P de la 15 zone de recouvrement. Ce calcul présente un autre avantage : la distance obtenue est quasiment indépendante des effets de réfraction atmosphérique, contrairement à celle obtenue par le calcul décrit dans le préambule. On a illustré figure 4 l'erreur provoquée par la courbure concave du rayon réfracté 20 issu de la cible, lorsque la géolocalisation est réalisée selon l'état de la technique. Lorsque la distance est calculée par triangulation selon l'invention, le point estimé E de la cible est quasiment à la verticale de la cible réelle P : E est à l'intersection de C1V1 et C2V2, V1 et V2 étant respectivement les points visés à partir des positions C1 et C2 de l'aéronef, comme illustré figure 25 7. La figure 4 est une vue de la figure 7 dans un plan contenant les points C1, C10 et V1 ou C2, C20 et V2. C10 et C20 sont les projections au sol des positions C1 et C2 de l'aéronef. L'erreur introduite par l'utilisation du point estimé E est négligeable. En reprenant l'exemple précédent avec une erreur due à la 30 courbure du rayon égale à 900 m, on obtient une erreur sur la distance estimée de ù 9 m. L'invention a été décrite pour un angle de gisement G nul ; elle s'applique également lorsque la ligne de visée est dépointée en gisement d'un angle G non nul comme illustré figure 8. Dans ce cas, la distance d est sensiblement égale à (B cos G) / (a2 û ai), soit (B cos G)/((32- (31 + 82-El). Et l'on a X=dcosGsinR Y = d sin G. On réalise aussi selon l'invention, un géo-référencement d'une 5 zone en appliquant ce procédé de géo-localisation à plusieurs cibles de la zone à référencer, et en interpolant pour les autres cibles de la zone. Selon un mode de réalisation, le procédé selon l'invention est mis en oeuvre dans une nacelle de reconnaissance montée sur un aéronef. La 10 nacelle 10 représentée figure 9 comprend un capteur optronique passif 1 et un coffret de traitement 2 qui comporte des cartes électroniques de traitement 3. De manière classique, le capteur est orientable et acquiert les images ; le procédé est mis en oeuvre par un logiciel implémenté sur une ou plusieurs des cartes électroniques. 15 REVENDICATIONS 1. Procédé de géo-localisation d'au moins une cible P immobile à partir d'un aéronef au moyen d'un capteur optronique passif, caractérisé en ce qu'il comprend les étapes suivantes : acquisition par le capteur d'au moins une image 11 contenant la cible P à partir d'une position C1 de l'aéronef et d'une image 12 contenant la cible P à partir d'une position C2 de l'aéronef, ces images 11 et 12 ayant une zone de recouvrement, identification dans ladite zone de recouvrement d'au moins une cible P commune aux deux images 11 et 12 et détermination de la 1 o position de chaque cible P dans chacune des deux images, calcul de la distance d entre chaque cible P et un point C de coordonnées connues par rapport à C1 et C2, en fonction de l'angle (31 entre une direction de référence et la ligne de visée de l'image 11, l'angle [32 entre la même direction de référence et la ligne de 15 visée de l'image 12, et de la position de chaque cible P dans l'image 11 et dans l'image 12, P étant défini dans un référentiel terrestre par son abscisse X et son ordonnée Y par rapport à la projection au sol de l'aéronef C, calcul pour chaque cible de X et Y en fonction des distances d. 20 2. Procédé de géo-localisation selon la revendication précédente, caractérisé en ce que la position de chaque cible P dans les images 11 et 12 est calculée en nombre de pixels et convertie respectivement en un angle El et un angle E2. 25 3. Procédé de géo-localisation selon l'une quelconque des revendications précédentes, caractérisé en ce que la ligne de visée en C présentant un angle de roulis R et un angle de gisement G dans un référentiel terrestre, il comprend une étape d'acquisition des angles de roulis 30 R et de gisement G et en ce que la position de chaque cible P est calculée en fonction de ces angles de roulis R et de gisement G. 4. Procédé de géo-localisation selon la revendication précédente, caractérisé en ce que la distance d est sensiblement égale à (B cos G)/(f32- R1 + E2-E1), B étant la distance entre les points Cl et C2. 5. Procédé de géo-localisation selon l'une quelconque des revendications 3 ou 4, caractérisé en ce que les angles de roulis R et de gisement G étant déterminés, la position de chaque cible P par rapport à la projection au sol de l'aéronef en C est donnée par X = d cos G sin R et Y = d sin G. 6. Procédé de géo-localisation selon l'une quelconque des revendications précédentes, caractérisé en ce que la géo-localisation est effectuée à bord de l'aéronef. 7. Procédé de géo-référencement d'une zone, caractérisé en ce qu'il comprend la répétition du procédé de géo-localisation selon l'une des revendications précédentes pour différentes cibles de la zone. 20 8. Coffret de traitement (2) comportant au moins une carte électronique (3), caractérisé en ce qu'au moins une carte électronique (3) comprend des moyens de mise en oeuvre du procédé selon l'une des revendications précédentes. 25 9. Nacelle de reconnaissance (10) destinée à être installée sur un aéronef, caractérisée en ce qu'elle comprend un coffret de traitement selon la revendication précédente.15