CORRELATION NUMERIQUE D'IMAGES QUELCONQUES SELON LES LIGNES 
QUASI-EPIPOLAIRES,PAR APPROXIMATIONS SUCCESSIVES. 
par G. de MASSON d'AUTUME 
Ingénieur général géographe honoraire 
France 
$ 
Commission III, Groupe de travail III/4 
1. Introduction 
Bien que conçue et mise au point en vue des traitements des images 
du satellite SPOT, dont le lancement est prévu pour 1985, la procédure décrite 
est tout à fait générale ; ses caractéristiques principales sont : 
- un ré-échantillonnage préalable des deux images réduisant au minimum leurs 
différences et tel que l'homologue d'un point quelconque de l'image 1, se 
trouve sur la ligne de méme numéro de l'image 2 ; 
- une procédure d'approximations successives au cours de laquelle les deux ima- 
ges subissent des corrections géométriques et radiométriques, jusqu'à deve- 
nir pratiquement identiques ; 
- la possibilité de traiter des images ayant des propriétés radiométriques to- 
talement différentes, par réduction aux contours ; 
- l'utilisation systématique de l'interpolation par grille élastique pour le 
filtrage des mesures aberrantes et la reconstitution des portions de surface 
dépourvues d'accidents. 
2 - Ré-échantillonnage selon les lignes quasi-épipolaires [2] 
2.1 - Théorie générale (fig.1) 
Les deux images sont ré-échantillonnées de telle façon que, si un 
point quelconque de l'espace-objet a pour coordonnées-pixel ii, A et ior jo: 
la quantité (i,-i,) est du second ordre par rapport à la parallaxe P - (j,-3j,). 
Un tel ré-échantillonnage est toujours possible, méme en l'absence de 
points connus. Soit en effet un jeu de valeurs des paramétres de prise de vues 
tel que les rayons homologues s'intersectent ; c'est le cas des valeurs exac- 
tes des paramétres, généralement inacccessibles, mais tout jeu satisfaisant la 
condition d'intersection fait aussi bien l'affaire. Le lieu des points d'inter- 
section est une surface S. Soit Sy une surface régulière voisine de S ; selon 
l'étendue de la zone à traiter et la for- 
me de S, So peut étre un plan, un cylin- 
dre, une sphére, etc ... 
Soit un systéme de coordonnées i, j 
défini sur S,, tel que les lignes i = ct, 
soient tangentes en tout point m à la 
droite mt intersection du plan tangent à 
S, en M avec le plan mt, t, des tangentes 
aux deux rayons passent par M. Les va- 
leurs entières de i, j définissent une 
grille de ré-échantillonnage. Pour chaque 
noeud de cette grille, la loi de forma- 
tion d'image permet de calculer les coor- 
données-pixel p, q de l'image d'origine 
et la densité Gij est interpolée pour p, 
q (on peut aussi prendre la densité du plus proche pixel). 
  
    
Fig. 1 
Soit alors M un point-objet quelconque, my (i,, J,) et m, (i,, jo) les 
points d'impact sur So des rayons (1) et (2) passant par M. Si M se déplace 
sur le rayon (1) en tendant vers m,, le rayon (2) engendre une surface I dont