Les meilleurs résultats ont été obtenus avec le cofficient de corrélation
classique
C = Covariance //vax . (1) x var, (2)
3.2. Recherche du maximum de corrélation
Le calcul donne, pour chaque valeur de j, n valeurs C(j'). Une pro-
cédure simple recherche le maximum de la fonction C(x), prenant les valeurs
C(j') pour x = j'. Si x max est l'abcisse du maximum et C max sa valeur, la
parallaxe est égale a XT j et le poids Q qui lui est attribué est le produit
de C max par la courbure en x, de la fonction C (x) . S'il n'y a pas de. maxi-
mum, ou plusieurs maximums dont aucun n'est prépondérant, ou si C max est né-
gatif, la parallaxe est prise égale à O0 avec un poids nul.
4. Filtrage et lissage de la fonction parallaxe
Le modèle numérique P (i,j) de la fonction parallaxe brute ainsi ob-
tenue est filtré et lissé par la procédure de la grille élastique [1]. Aux
équations d'observation P(i,j) = Km (i,j) - j affectées du poids Q (i,j) on
ajoute les équations "de régularité" :
p (1,4-1) = 29 (4,5) + Pti,s + 1) =D
et PUM IY Un IP, AY Te Piet, Y= 0
affectées du poids 1.
Le systéme est résolu par la méthode des gradients conjugués au moyen
d'un A.P. (Processeur de tableaux). L'A.P. utilisé ne pouvant traiter que des
blocs de 10 000 pixels au maximum, le modéle numérique est décomposé automati-
quement en blocs de 100 x 100, avec des recouvrements assurant la continuité.
Un premier passage sert à éliminer les mesures aberrantes, auxquelles est at-
tribué un poids nul. Le deuxiéme passage donne les résultats définitifs.
5. Ré-échantillonnage de l'image (2) avec correction de la parallaxe et
harmonisation des densités
L'image (2) peut étre corrigée ligne par ligne en fonction de la pa-
rallaxe en affectant au pixel j de l'image corrigée la densité 6 interpolée
pour j + P (i,j). Il est toutefois préférable, pour éviter une cascade de ré-
échantillonnages nuisibles à la qualité de l'image, d'opérer la correction à
partir de l'image d'origine. Pour une ligne i quelconque, l'interpolation
dans le modéle de déformation donne les coordonnées-pixel p, q de l'image
d'origine correspondant ou pixel j + P (i,j) et la densité interpolée pour
P, q est affectéé en pixel j de l'image corrigée.
On peut profiter du ré-échantillonnage pour harmoniser les densités
de l'image (1) avec celles de l'image (2) en leur appliquant une correction
de densité moyenne 4^ (i, j) et une correction de contraste K (i, j) selon
la formule
8 tu, jy A (1,5) -5.K ti, 3)
Les fonctions À et K étant définies par deux modèles numériques calculés
par moindres carrés et pouvant être construits sur la même grille que le mo-
dèle de déformation géométrique.
Après correction, les parallaxes subsistant entre les deux images (fig.
ne représentent plus que les imperfections de la fonction parallaxe utilisée.
Au stade suivant, le nombre de pas d'exploration pourra être considérable-
ment réduit. En outre la corrélation sera bien meilleure, méme avec des fené-
tres relativement grandes, et le pic du maximum de corrélation beaucoup plus
accusé.
La procédure d'approximations successives permet de tenir un compte
exact des déformations dues au relief, sans qu'il soit nécessaire de modéli-
d E» LP
X
A+” A^
x NS "
i N A U
A Eh P
3 e eo
RA 32s D ier
PE
WS ^6 Tit
RP RE 2
Yat 5002425
TES.
A" S "iR A3
VS QM vl
4 » y Gr A
NO