Brevet FR2950438A1 — publié le 20110325

Localisation tridimensionnelle de zone cible terrestre par fusion d'images prises par deux capteurs satellitaires La présente invention concerne d'une manière générale la télédétection spatiale d'images notamment pour l'observation de la terre. Plus particulièrement, elle concerne un procédé et un système de localisation tridimensionnelle de zone cible terrestre par fusion d'images prises par deux capteurs ayant des propriétés radiométriques différentes, comme un capteur radar et un capteur optique, embarqués dans des satellites distincts. 1 o Deux capteurs différents produisent des images différentes de la même scène terrestre observée. Si les images sont traitées séparément, les localisations précises des pixels dans les images ne sont pas possible dans le cas général. Il est nécessaire de recourir à une information supplémentaire consistant en un modèle numérique d'altitude de terrain MNT ("Digital 15 Elevation Model" DEM en anglais) qui est interpolé pour trouver la localisation de chaque petite zone cible pour chaque image en utilisant les caractéristiques techniques de chaque capteur. Un traitement de ces images ensemble après corrections de diverses distorsions produit automatiquement l'appariement de points homologues dans 20 les deux images qui servent à estimer un modèle de déformation analytique fondé sur un calcul de mesure de similarité et au moins une interpolation pour chaque paires de points homologues. Dans l'article de Jordi Inglada and Alain Giros, "On the Possibility of Automatic Multisensor Image Registration," IEEE TRANSACTIONS ON 25 GEOSCIENCE AND REMOTE SENSING, Vol. 42, No. 10, pages. 2104-2120, Oct. 2004, le chapitre VII présente l'estimation d'une topographie en utilisant un capteur optique et un capteur radar sur les plateformes de deux satellites différents pointés vers un point cible. La hauteur h du point cible sur la terre considérée comme plate est donnée par la formule suivante : 30 h = Ax (tg a,)/((tg a, tg a2) + 1)), où a, et a2 dénotent les angles d'incidence du point cible vu par les capteurs par rapport à la terre, et Ax le décalage entre les images du point cible fournies par les capteurs. Le décalage Ax est égal à la somme de la distance entre la projection du point cible sur la terre et le point d'intersection de l'axe de visée 35 du capteur optique avec la terre et de la distance entre la projection du point cible sur la terre et la projection du point cible sur la terre perpendiculairement à l'axe de visée du capteur radar, en supposant un front d'onde plan émis par le capteur radar. Une telle estimation n'est pas assez précise pour localiser tridimensionnellement chaque petite zone cible, appelée pixel, dans les images afin d'enregistrer une topographie tridimensionnelle de la scène. Le capteur radar mesure implicitement la distance entre le centre de phase de l'antenne et la zone cible. Pour cette raison, l'angle d'incidence associé au capteur radar ne peut pas être mesuré avec précision. Si le capteur optique et le capteur radar sont embarqués sur la plate-forme d'un unique satellite, soit a, = a2 = a, la formule précédente s'écrit : h = Ax (tg a)/((tg a)2 + 1)), soit h = Ax (sin 2a)/2. Pratiquement, les capteurs ne seront pas montés pour avoir un angle d'incidence identique. Il y aura toujours un biais dans le calcul de localisation. Chacune de ces formules signifie que la hauteur h du point cible à mesurer est proportionnelle au décalage mesuré entre les images du point cible et indépendante des altitudes des capteurs. Cette approche trigonométrique de la fusion de données issues d'un capteur optique et d'un capteur radar est aussi trop simpliste. Pour avoir une position tridimensionnelle, il faut prendre en compte d'autres transformations liées aux positions et orientations des capteurs. L'objectif de l"invention est de localiser tridimensionnellement avec une plus grande précision une zone cible terrestre par la fusion de données acquises issues d'un capteur optique et d'un capteur radar, en dépendance des positions des satellites supportant les capteurs. Pour atteindre cet objectif, un procédé pour localiser tridimensionnellement une zone cible terrestre au moyen d'un capteur optique et d'un capteur radar embarqués respectivement dans des satellites, comprenant des mesures des positions des satellites lorsque les capteurs pointent vers la zone cible, est caractérisé en ce qu'il comprend une mesure d'une direction de visée du capteur optique pointant vers la zone cible, une 3o mesure de la distance entre la zone cible et le capteur radar pointant vers la zone cible, et une détermination de la distance entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar et de la direction de visée du capteur optique. La distance entre la zone cible et le capteur optique est ainsi déterminée 35 en fonction de quelques données simples mesurées précisément notamment à bord des satellites co-localisés, par une approche de fusion de deux équations relatives aux capteurs, sans faire des approximations sur la position des satellites, ce qui permet de fournir avec une grande précision la position de la zone cible visée par les capteurs. L'invention peut être utilisée pour établir un modèle numérique d'altitude de terrain tridimensionnelle de la scène ou pour des projections cartographiques précises à partir de l'enregistrement des images fournies respectivement par les satellites en utilisant des critères mutuels pour harmoniser les disparités entres deux images d'une même zone cible prises à des instants a priori différents par le capteurs. En outre, l'usage actuel seulement de capteurs optiques, ou de capteurs radar, à très haute résolution et la nécessité de produire des modèles numériques d'altitude de qualité élevée justifient l'exploitation industrielle de l'invention qui, grâce à la combinaison de quelques données fournies par un capteur optique avec quelques données fournies par un capteur radar, est peu complexe. Comme on le verra dans la suite, la détermination de la distance entre la zone cible et le capteur optique comprend une résolution d'une équation quadratique de cette distance qui résulte de la combinaison de deux équations liées aux capteurs. La résolution de l'équation quadratique revient à déterminer les deux solutions de distance suivantes : S= ûoT(PR-Po)±([ûoT (PR-Po )12- I PR-po 12 -F dR2 )1'2, dans lesquelles po et pR dénotent les positions des satellites embarquant le capteur optique et le capteur radar, dR dénote la distance entre la zone cible et le capteur radar et ûo dénote la direction de visée du capteur optique vers la zone cible, et à sélectionner celle des deux solutions de distance qui est inférieure au rayon de la terre. L'invention a aussi pour objet un système et un dispositif de localisation tridimensionnelle d'une zone cible terrestre. Le système comprend un capteur optique et un capteur radar embarqués respectivement dans des satellites, un moyen pour mesurer une position du satellite embarquant le capteur optique lorsque le capteur optique est apte à pointer vers la zone cible, et un moyen pour mesurer une position du satellite embarquant le capteur radar lorsque le capteur radar est apte à pointer vers la zone cible, et est caractérisé en ce qu'il comprend un moyen pour mesurer une direction de visée du capteur optique lorsque le capteur optique est apte à pointer vers la zone cible, un moyen pour mesurer la distance entre la zone cible et le capteur radar lorsque le capteur radar est apte à pointer vers la zone cible, et un dispositif pour déterminer la distance entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar et de la direction de visée du capteur optique. Le dispositif est apte à localiser tridimensionnellement une zone cible terrestre vers laquelle sont aptes à pointer un capteur optique et un capteur radar embarqués respectivement dans des satellites ayant des positions mesurées, et est caractérisé en ce qu'il est apte à déterminer la distance entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar lorsque le capteur radar est 1 o apte à pointer vers la zone cible, et de la direction de visée du capteur optique lorsque le capteur optique est apte à pointer vers la zone cible. D'autres caractéristiques et avantages de la présente invention apparaîtront plus clairement à la lecture de la description suivante de plusieurs 15 réalisations de l'invention données à titre d'exemples non limitatifs, en référence aux dessins annexés correspondants dans lesquels : - la figure 1 montre schématiquement un système de localisation tridimensionnelle de zone cible terrestre selon l'invention ; et - la figure 2 est un algorithme d'un procédé de localisation 20 tridimensionnelle de zone cible terrestre selon l'invention. En référence à la figure 1, un système de localisation tridimensionnelle de zone cible terrestre comprend un capteur optique CO, un capteur radar CR et un dispositif de localisation tridimensionnelle de zone terrestre DL. Le capteur 25 optique CO est embarqué sur la plate-forme d'un premier satellite STO. Le capteur optique comporte un système d'imagerie interne ayant une matrice de photodétecteurs à transfert de charge CCD ("Charged Coupled Device" en anglais) situé dans le plan focal d'un système optique, dont l'axe de visée a une direction par rapport à la terre mesurable dans le capteur, et une partie 30 électronique analogique et numérique pour restituer et mémoriser des images bidimensionnelles. Le capteur radar CR est un radar imageur à synthèse d'ouverture SAR ("Synthetic Aperture Radar" en anglais) embarqué sur la plate-forme d'un deuxième satellite STR. Le dispositif DL est par exemple un ordinateur situé au sol et par exemple relié, ou intégré, à au moins une station 35 terrestre de poursuite multi-satellitaire ST en communication avec les satellites STO et STR et synchrone avec les horloges des satellites. Le capteur optique CO et le capteur radar CR transmettent des signaux so et SR relayés par la station ST vers le dispositif DL afin que la station acquière périodiquement diverses données numériques transmises par les satellites, comme les positions géographiques des satellites mesurées par ceux-ci avec des données d'horodatage. Les signaux so et sR peuvent supporter à des instants, différents a priori, respectivement une image optique IO et une image radar IR avec leurs dates respectives, couvrant chacune une zone cible terrestre ZC. Les instants de prise de l'image optique et de l'image radar peuvent différer de quelques heures, jours ou mois. Chaque capteur mémorise des images respectives et les retransmet classiquement avec des données du capteur et du satellite embarquant le capteur au dispositif DL via la station ST lorsque celle-ci est en exposition directe de l'antenne du satellite pour que le dispositif DL enregistre les données, puis les traite ultérieurement selon l'invention. La zone cible ZC est de très petite dimension par rapport à la surface de la terre E. Pour cette raison, la zone cible peut être considérée comme confondue avec un point et être appelée "pixel". Par exemple, la zone ZC a un diamètre de quelques mètres et est définie par son point central sur la surface de la terre. La zone ZC est à localiser par le dispositif DL par rapport à un référentiel tridimensionnel OXYZ qui est par exemple fixe et lié à la terre et dont l'origine O est au centre de la terre. Dans le référentiel OXYZ, la position recherchée de la zone cible ZC est définie par un vecteur de position pc (po(X), po(Y), po(Z)), la position du capteur optique satellitaire CO est définie par un vecteur de position po (po(X), po(Y), po(Z)) lorsque le capteur CO pointe vers la zone cible ZC, et la position du capteur radar satellitaire CR est définie par un vecteur de position pR (pR(X), pR(Y), pR(Z)), lorsque le capteur CR pointe vers la zone cible ZC. Par rapport au capteur optique CO, le vecteur de position pc de la zone cible terrestre ZC dépend du vecteur de position po du capteur optique CO, de coordonnées x et y de la zone cible ZC dans un référentiel situé sur le plan focal du système optique interne au capteur optique CO, de la distance focale f du système optique, d'une matrice de rotation R et d'un facteur d'échelle k, selon l'équation vectorielle suivante : p-c=p-o+kR qui peut s'écrire : x y -f x Y = R-«P, - po )lk (2). ùf La matrice de rotation R transforme les coordonnées d'un vecteur de direction dans le plan focal en les coordonnées de ce vecteur de direction dans le référentiel OXYZ lié à la terre E : _ r11 r21 r13 X12 r22 r23 r13 r23 r33 Les coefficients rä à r33 de la matrice de rotation R dépendent notamment de l'orientation de la plate-forme supportant le capteur optique CO dans le satellite STO, de la position et de l'orientation de miroirs tournants sur le chemin optique dans le système optique, et de l'orientation du capteur optique par rapport aux récepteurs dans la station de poursuite ST reliée au dispositif de détermination de localisation DL et donc par rapport au référentiel OXYZ fixe. La résolution des trois équations de colinéarité déduites de l'équation vectorielle (2) conduit aux coordonnées x et y de la zone cible ZC en fonction de la distance focale f et des différences entre les coordonnées de la position pc (pc(X), pc(Y), pc(Z)) de la zone cible ZC et de la position po (po(X), po(Y), po(Z)) du capteur optique CO, après avoir remplacé le facteur k déduit de la troisième équation de colinéarité dans les deux premières équations de colinéarité : x = -f '1(0 (X)1+r21[pc(Y)+po(Y)1+r31[pc(Z)ùpo(Z)] (3) 'i3[Pc(X)ùpo(X)] +r23[Pc(Y)ùpo(Y)]+r33[Pc(Z)ùpo(Z)] r12[PC(X)po(X)]+r22[Pc(Y)+po(Y)]+r32[Pc(Z)po(Z)] 'i3 [Pc (X) ù po (X)] + r23 [Pc (Y) ù po (Y) ] + r33 [PC (Z) ù po (Z) Pour une seule image optique, les équations (3) et (4), qui définissent les coordonnées x et y de la zone cible visée par le capteur optique CO dépendant des propriétés du capteur optique, ne peuvent être résolues sans recourir à un modèle d'altitude MNT de la surface de la terre E. L'équation vectorielle (1) peut s'écrire différemment selon l'équation suivante (5) dans laquelle la position de la zone cible ZC s'exprime en fonction de la distance 8 en mètre entre la zone cible ZC et le capteur optique CO et d'un vecteur unitaire ûo suivant la direction de visée depuis le capteur optique vers la zone cible dans le référentiel OXYZ : pc=po+5ûo y=-f Par rapport au capteur radar CR ayant un vecteur de position connu pR et situé à une distance dR de la zone cible terrestre ZC déduite du retour des échos radar reçus par le capteur CR, le vecteur de position pc de la zone cible est déduit d'une première équation géométrique dans le référentiel OXYZ : dR = 1 PR-PC I (6) et d'une deuxième équation déduite du cône de fréquence Doppler ayant pour sommet le capteur CR et relative à la fréquence Doppler fdop en hertz du satellite STR supportant le capteur radar CR et ayant la vitesse VR par rapport 1 o à la vitesse vc de la zone cible ZC dans le référentiel OXYZ, pour une longueur d'onde radar en mètre correspondant à la fréquence centrale des signaux radar émis par le capteur radar CR du satellite STR vers la zone cible : Z (P -15 Y' fop = X wR ûvc) (7), ~PR û Pc ) T désignant l'opérateur de transposition matricielle. 15 D'une manière similaire aux deux équations (3) et (4), les équations (6) et (7) sont ambiguës pour une seule image radar. Des données externes sous la forme d'un modèle numérique d'altitude MNT de la scène observée autour de la zone cible sont nécessaires pour résoudre complètement les équations (6) et (7) dans lesquelles les coordonnées x et y de la zone cible dépendent des 20 propriétés du capteur radar, afin de localiser la zone cible par ses coordonnées. La géo-localisation tridimensionnelle de la zone cible terrestre par chacun des deux systèmes d'équations (3-4) et (6-7) est ainsi relativement complexe 25 puisqu'elle recourt à un modèle numérique d'altitude MNT qui doit être interpolé pour localiser précisément la zone cible. Pour chacun de ces deux systèmes d'équations, ceci peut être évité en collectant deux images de la zone cible selon le même procédé, optique ou radar, depuis un satellite correspondant équipé d'un capteur optique ou d'un capteur radar. Dans le cas 3o d'un satellite avec capteur optique, la géo-localisation est appelée stéréogrammétrie. Dans le cas d'un satellite avec capteur radar, la géolocalisation est appelée radargrammétrie et revient à trouver l'intersection de deux arcs pour localiser tridimensionnellement la zone cible. Avec l'apparition des systèmes de géo-localisation par satellites GPS ("Global Positioning 35 System" en anglais) et pour le cas de satellites équipés de capteurs optiques, ces difficultés sont quelque peu atténuées par une modélisation de l'orbite et de l'orientation des capteurs optiques ce qui conduit à la possibilité d'un traitement automatique d'images sans points de contrôle externes. Selon l'invention, le dispositif de localisation DL fusionne des données du capteur optique CO dans le premier satellite STO et des données du capteur radar CR dans le deuxième satellite STR, en remplaçant le vecteur de position pc de la zone cible dans l'équation (6) par sa valeur déduite de l'équation (5) afin d'obtenir la distance recherchée 8 entre la zone cible ZC et le capteur optique CO : dR = I P,-Po - 6 ûo I (8) Cette équation représente l'intersection du vecteur pR - po - 8 ûo avec une sphère de rayon dR centrée sur le capteur radar CR. On obtient en élevant au carré les deux membres de l'équation (8) : dR2 = (PR-Po - b uo)T (PR-Po - 8 ûo) puis en regroupant les différences des vecteurs de position du capteur radar et du capteur optique : dR2=1% -po)T(PR-Po)ù28 ûoT (pR-po)+82 uoT uo, une équation du second degré en 8 : 82ù2 ûoT (pR-po)8+(PR-Po)T(PR-Po)-dR2=0, dont les deux solutions sont : =ûoT (pR-po)+([ûoT(PR-Po)]2- I PR-Po I2+dR2)112 82= ûoT (PR-Po)-(LûoT (PR-Po)l2- I PR-Po I2+dR2)1'2 (9). Parmi les solutions 81 et 62, celle positive et inférieure au rayon RE de la terre est la distance correcte 8 depuis le capteur optique CO à la zone cible 2 5 ZC. En se référant maintenant à la figure 2, le procédé de localisation tridimensionnelle de zone cible terrestre selon l'invention comprend des étapes El à E7 afin que le dispositif DL détermine la distance 8. 30 Initialement à l'étape El, la station de poursuite ST acquiert périodiquement diverses données numériques supportées par les signaux so et SR émis par les satellites et incluant notamment les positions géographiques po et pR des satellites STO et STR, confondues avec celles des capteurs CO et CR, mesurées par les satellites dans le référentiel fixe OXYZ. En variante, 35 les positions des satellites STO et STR sont mesurées par la station ST en liaison avec des balises sur la terre. A l'étape E2, une image IO prise par le capteur optique CO et, en association avec l'image prise 10, notamment la position po du satellite STO, l'orientation ûo du capteur optique et des données d'horodatage HO de la prise de l'image IO sont transmises via le signal so par le satellite STO à la station de poursuite ST qui retransmet l'image IO avec les données u0, po et HO associées à l'image IO au dispositif DL qui les mémorise dans une base de données. De même, à l'étape E3 qui peut être antérieure ou postérieure à l'étape E2, le satellite STR transmet via le signal SR une image IR prise par le capteur radar CR et, en association avec l'image prise IR, notamment la distance dR du capteur radar à la terre E et la position pR mesurées par le capteur radar et des données d'horodatage HR de la prise de l'image IR à la station de poursuite ST qui retransmet l'image IR et les données précédentes associées à l'image IR au dispositif DL qui les mémorise dans la base de données. Le dispositif DL classent les images mémorisées IO et IR avec leurs données mémorisées associées telles que ûo, po, HO et dR, pR, HR en fonction de leurs données d'horodatage. Pour deux images IO et IR pour des prises d'image d'une zone cible prédéterminée ZC sur la terre, le dispositif DL les apparie avec leurs données, à l'étape E4. Puis à l'étape E5, le dispositif DL détermine selon un programme d'ordinateur préenregistré les solutions de distance Si et b2 entre la zone cible prédéterminée ZC et le capteur optique CO selon les relations (9), en fonction des données associées aux deux images telles que les positions po et pR des satellites, la distance dR entre la zone cible et le capteur radar CR et la direction de visée ûo du capteur optique CO vers la zone cible ZC. A l'étape E6, le dispositif DL compare chacune des solutions de distance 81 et b2 au rayon RE de la terre et sélectionne la solution qui est positive et inférieure au rayon RE de la terre et qui constitue la distance correcte 6 entre le capteur optique CO et la zone cible ZC lors de la prise de l'image 10. Finalement, le dispositif DL déduit à l'étape E7 les coordonnées du vecteur de position pc = po + 8 ûo de la zone cible ZC dans le référentiel OXYZ et particulièrement son altitude en fonction de la distance 6, selon l'équation (5). La distance correcte 6 et la position pc sont déterminées directement selon l'invention, sans l'aide d'un modèle numérique d'altitude MNT. L'invention décrite ici concerne un procédé, un système et un dispositif pour localiser tridimensionnellement une zone cible terrestre. Selon une implémentation, les étapes du procédé de l'invention sont déterminées par les instructions d'un programme d'ordinateur incorporé dans le dispositif. Le programme d'ordinateur apte à être mis en oeuvre dans le dispositif de localisation de l'invention comporte des instructions de programme qui, lorsque ledit programme est exécuté dans le dispositif dont le fonctionnement est alors commandé par l'exécution du programme, réalisent une détermination de la distance entre une zone cible et un capteur optique conformément au procédé de l'invention. En conséquence, l'invention s'applique également à un programme d'ordinateur, notamment un programme d'ordinateur enregistré sur ou dans un support d'enregistrement lisible par un ordinateur et tout dispositif de traitement de données, adapté à mettre en oeuvre l'invention. Ce programme peut utiliser n'importe quel langage de programmation, et être sous la forme de code 1 o source, code objet, ou de code intermédiaire entre code source et code objet tel que dans une forme partiellement compilée, ou dans n'importe quelle autre forme souhaitable pour implémenter le procédé selon l'invention. Le programme peut être téléchargé dans le dispositif via un réseau de communication, comme internet. 15 Le support d'enregistrement peut être n'importe quelle entité ou dispositif capable de stocker le programme. Par exemple, le support peut comporter un moyen de stockage sur lequel est enregistré le programme d'ordinateur selon l'invention, tel qu'une ROM, par exemple un CD ROM ou une ROM de circuit microélectronique, ou encore une clé USB, ou un moyen d'enregistrement 20 magnétique, par exemple une disquette (floppy disc) ou un disque dur. REVENDICATIONS1 û Procédé pour localiser tridimensionnellement une zone cible terrestre (ZC) au moyen d'un capteur optique (CO) et d'un capteur radar (CR) embarqués respectivement dans des satellites (STO, STR), comprenant des mesures (El) des positions (po, pR) des satellites lorsque les capteurs pointent vers la zone cible, caractérisé en ce qu'il comprend une mesure (E2) d'une direction de visée (ûo) du capteur optique pointant vers la zone cible, une mesure (E3) de la distance (dR) entre la zone cible et le capteur radar 1 o pointant vers la zone cible, et une détermination (E5, E6) de la distance (8) entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar et de la direction de visée du capteur optique. 15 2 û Procédé conforme à la revendication 1, selon lequel la détermination de la distance (8) entre la zone cible (ZC) et le capteur optique (CO) comprend une résolution (E5) d'une équation quadratique de la distance (8) entre la zone cible et le capteur optique dépendant des positions (po, PR) des satellites, de la distance (dR) entre la zone cible et le capteur radar et de la direction de 20 visée (ûo) du capteur optique vers la zone cible en deux solutions de distance, et une sélection (E6) de l'une des deux solutions de distance qui est inférieure au rayon de la terre. 3 û Procédé conforme à la revendication 1 ou 2, selon lequel la 25 détermination de la distance (8) entre la zone cible (ZC) et le capteur optique (CO) comprend une détermination (E5) des deux solutions de distance suivantes : = uoT(PR-Po)±(~ûoT (PR-Po)12_ I PR-PO I2+dR2)1'2 dans lesquelles po et pR dénotent les positions des satellites (STO, STR) 3o embarquant le capteur optique (CO) et le capteur radar (CR), dR dénote la distance entre la zone cible et le capteur radar et ûo dénote la direction de visée du capteur optique vers la zone cible, et une sélection (E6) de l'une des deux solutions de distance qui est inférieure au rayon de la terre. 35 4 û Système pour localiser tridimensionnellement une zone cible terrestre (ZC), comprenant un capteur optique (CO) et un capteur radar (CR) embarqués respectivement dans des satellites (STO, STR), un moyen (STO) pour mesurerune position (po) du satellite (STO) embarquant le capteur optique lorsque le capteur optique est apte à pointer vers la zone cible, et un moyen (STR) pour mesurer une position (1%) du satellite (STR) embarquant le capteur radar lorsque le capteur radar est apte à pointer vers la zone cible, caractérisé en ce qu'il comprend un moyen (CO) pour mesurer une direction de visée (ù >o) du capteur optique lorsque le capteur optique est apte à pointer vers la zone cible, un moyen (CR) pour mesurer la distance (dR) entre la zone cible et le capteur radar lorsque le capteur radar est apte à pointer vers la zone cible, et un dispositif (DL) pour déterminer la distance (8) entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar et de la direction de visée du capteur optique. 5 ù Dispositif (DL) pour localiser tridimensionnellement une zone cible terrestre (ZC) vers laquelle sont aptes à pointer un capteur optique (CO) et un capteur radar (CR) embarqués respectivement dans des satellites (STO, STR) ayant des positions mesurées (po, pR ), caractérisé en ce qu'il est apte à déterminer la distance (8) entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar lorsque le capteur radar est apte à pointer vers la zone cible, et de la direction de visée du capteur optique lorsque le capteur optique est apte à pointer vers la zone cible. 6 ù Programme d'ordinateur apte à être mis en oeuvre dans un dispositif (DL) pour localiser tridimensionnellement une zone cible terrestre (ZC) vers laquelle sont aptes à pointer un capteur optique (CO) et un capteur radar (CR) embarqués respectivement dans des satellites (STO, STR) ayant des positions mesurées (po, pR ), ledit programme étant caractérisé en ce qu'il comprend des instructions qui, lorsque le programme est exécuté dans ledit dispositif, réalisent une détermination (E5, E6) de la distance (8) entre la zone cible et le capteur optique en fonction des positions des satellites, de la distance entre la zone cible et le capteur radar lorsque le capteur radar est apte à pointer vers la zone cible, et de la direction de visée du capteur optique lorsque le capteur optique est apte à pointer vers la zone cible.