cartographie
La génération automatique de spatiocartes
Simon BAILLARIN et Marc BERNARD
L’instrument HRS, qui permet d’acquérir des couples
d’images stéréoscopiques quasi-simultanées
sur de très vastes étendues, est l’une
des innovations majeures de SPOT 5.
L’exploitation de ses images (assurée conjointement
par IGN et Spot Image au travers de la coédition
du produit Reference3D) a permis l’avènement
d’un outil entièrement
automatique de production mots clés
d’orthoimages d’une précision SPOT, HRS, ortho-
rectification, processus
géométrique unique sur le
automatique Figure 1 : Emprises des différents objets géographiques
marché: ANDORRE.
mis en jeu dans le processus
L
a fonction d’ANDORRE est d’assurer l’orthorectification mosaiquées de façon à obtenir une seule orthoimage et un
automatique de toute scène SPOT, en s’appuyant sur seul MNT. Le plus souvent, cela met en jeu 1 ou 4 dalles
Reference3D : Reference3D, puisque leur taille avoisine couramment les
• L’orthoimage HRS incluse dans Reference3D, aussi précise 10,000 km2 (1° de longitude par 1° de latitude).
qu’une carte régulière (précision horizontale de localisation
Ensuite, le modèle physique de la scène SPOT et le MNT issu
absolue de 15 m CE90), est utilisée comme source de calage.
de Reference3D sont utilisés afin de créer l’image “simulée”
L’idée est de retrouver automatiquement des points homo-
dans le plan focal. L’image ainsi créée prend exactement en
logues dans l’orthoimage Reference3D et dans la scène
compte les variations du relief. Ainsi le décalage mesuré entre
SPOT à traiter, et d’en déduire un modèle géomètrique de
l’image “simulée” et la scène SPOT réelle est uniquement dû
cette dernière.
à l’imprécision de localisation de son modèle physique (cal-
• le MNE DTED2 de Reference3D intervient ensuite pour l’or-
culé à partir des données auxiliaires).
thorectification de la scène SPOT.
La Figure 2 ci-dessous représente un extrait d’une scène
Le cœur algorithmique d’ANDORRE, dont l’idée est née des
SPOT4 HRV (IMAGE_N1A) et l’image simulée correspon-
travaux IGN sur la mission HRS, a été développé sous maî-
dante, le rectangle rouge délimitant une zone ayant les
trise d’ouvrage CNES. La première version de ces algorithmes
sera intégrée fin 2004 par Spot Image dans son atelier
[ANDORRE] de production d’orthoimages. Les adaptations
d’ANDORRE à des sources images autres que SPOT sont envi-
sagées.
Principe de la méthode automatique
d’orthorectification
La méthode d’ortho rectification automatique se décompose
en trois étapes principales :
• Génération de l’image ramenée au plan focal de SPOT
• Calcul de la correction de localisation de la scène SPOT par
correlation multi-échelle Figure 2 : Extraits d’une
• Ortho-rectification de l’image SPOT scène HRV (haut gauche) et
de son image simulée (haut
Génération de l’image “simulée”
droite) Le rectangle rouge
On nomme ainsi l’orthoimage Reference3D ramenée dans la (repris dans les zooms)
géométrie du capteur CCD au plan focal de SPOT. représente la même zone
Les dalles Reference3D touchant la scène SPOT à corriger sont de l’image
92 Revue XYZ • N°100 – 3e trimestre 2004
mêmes coordonnées dans les deux images. L’important déca- Facteur_initial = Precision_Loc / (Ech_image * (Taille_Fenetre/2) )
lage Est-Ouest parfaitement visible matérialise l’imprécision • Ech_image est le pas d’échantillonnage “pleine résolution”
de localisation du modèle de cette image HRV. de l’image SPOT à traiter (20 m, 10 m, 5 m). La dernière ité-
ration se fait à la résolution de Reference3D, soit environ 5 m.
Détermination de la correction à appliquer au
• Precision_Loc est la précision à 3 (du modèle physique de
modèle de localisation de la scène
localisation pour le satellite considéré (150 m pour SPOT 5,
La correction est calculée via un processus de corrélation 1 500 m pour SPOT2 à 4)
multi-échelle. A chaque étape on détermine un grand nombre •Taille_Fenetre est la taille de la fenêtre de recherche de la cor-
de points homologues entre la scène SPOT et son image rélation (5 pixels).
“simulée”. On en déduit le modèle de correction à appliquer
Par exemple, dans le cas d’une image SPOT4 panchroma-
pour ramener la scène SPOT dans la géométrie de l’image
tique, un niveau de dézoom initial de 1/64 est nécessaire (soit
simulée. Le modèle de correction est affiné à chaque niveau
7 niveaux de dézoom).
de la pyramide de résolution.
Pour chaque facteur de dézoom, la corrélation calcule les déca-
Les deux images initiales sont d’abord sous-échantillonées puis lages géométriques entre les deux images. Le résultat est filtré
ramenées progressivement à la pleine résolution. Le facteur ini- afin d’extraire un nombre de points homologues équi-répartis
tial de sous-échantillonnage est calculé pour conserver une
fenêtre de recherche d’une taille raisonnable (5 pixels):
sur l’ensemble de l’image et de modéliser le décalage mesuré.
Ce modèle est ensuite utilisé comme prédicteur de décalage ...
Figure 3 : Schématisation du processus multi-échelle utilisé dans ANDORRE
Revue XYZ • N°100 – 3e trimestre 2004 93
cartographie
Figure 4 : Scène corrigée
(gauche) et orthoimage de
référence (droite) Figure 5 : Orthoimage SPOT 5m (ci-dessus) calculée sur
La parfaite superposition le MNT Reference3D (ci-dessous)
des zones imagées dans le
rectangle (zooms) illustre la
précision du processus.
pour la corrélation au niveau d’échantillonnage supérieur.
Ce processus de corrélations multi-échelle converge vers un
modèle de recalage. Ce modèle est affiné en pleine résolution
(plusieurs itérations sont réalisées à F=1 afin de s’assurer de
la convergence du modèle).
Les points homologues obtenus à l’issue de cette étape sont
conservés dans le cas où une image identique (ou issu du
même segment de prise de vue) devrait être produite plus tard.
Ortho-rectification de l’image
La grille de rectification pour la scène SPOT est calculée en
combinant le modèle précédemment calculé avec la grille de
Contact
localisation associé au modèle de prise de vue inverse de Simon BAILLARIN
l’image SPOT, dans le référentiel WGS84, éventuellement CNES - Toulouse
combinée avec la grille de conversion entre WGS84 et la pro-
Marc BERNARD
jection cartographique finale souhaitée.
Spot Image
Enfin, la dernière opération consiste à orthorectifier la scène [Link]@[Link]
SPOT en utilisant cette grille.
La Figure 4 ci-dessous reprend la même zone illustrée Figure
2, une fois mise en projection UTM et recalée sur la dalle
Reference3D(tm). Le rectangle rouge a les mêmes coordon- ABSTRACT
nées dans les 2 images. On constate que les deux images sont Following the SPOT5 launch, Spot Image and the French
parfaitement superposables. National Geographic Institute (IGN) have decided to
La parfaite superposition des zones imagées dans le rectangle design a worldwide accurate database called
(zooms) illustre la précision du processus. Reference3D™ using data from the High Resolution
Stereoscopic SPOT5 instrument (HRS). Spot Image now
commercializes a system called ANDORRE to
Conclusion automatically produce ortho-rectified images thanks to
Reference3D™ data. ANDORRE has been designed to take
Par son fonctionnement entièrement automatique, ANDORRE advantage of Reference3D™ planimetric and altimetric
ramène les coût et délai de production de l’orthoimage au accuracy to automatically register and rectify any image
niveau de ceux pratiqués jusque là pour les images brutes. from SPOT satellites. In this framework, CNES is acting as
C’est donc une petite révolution : grâce à ANDORRE, l’orthoi- prime contractor to realize and industrialize the
mage, directement utilisable par les SIG des utilisateurs mais algorithms following a preliminary study undertaken by
IGN. This paper focuses on the algorithms and the
handicapée par un différentiel de prix, peut enfin réellement
architecture of the ANDORRE process.
devenir le produit “image spatiale” de base. ●
94 Revue XYZ • N°100 – 3e trimestre 2004