PROPOSITION D’UNE FONCTION PERMETTANT DE SUIVRE,
INDEPENDAMMENT DES CONTRASTES OMBRE LUMIERE,
po LES STRUCTURES GÉOLOGIQUES DANS LES RÉGIONS MONTAGNEUSES
ie
Catherine MERING*
et
Jean-François PARROT**
RÉSUMÉ
En région montagneuse, l'analyse numérique peut permettre de suivre dans les zones d'ombre les
traits structuraux que dessinent dans les zones éclairées les unités lithologiques. Pour ce faire, nous
avons choisi de rendre comparables les rapports des valeurs des pixels qui limitent deux unités litho-
logiques dans une zone d'ombre aux rapports des valeurs des pixels qui limitent ces deux mémes
unités dans une zone éclairée. Pour atteindre ce but, on applique aux valeurs des pixels consécutifs la
fonction suivante, soit dans la direction des lignes, soit dans celle des colonnes :
Log (sup (ID(K), ID (K«1)) * M4) « M,
| g[ID (K),ID (K--1)] - rM;
gen- Log (inf (ID (K), ID (K*1))* M4)
ID (K) et ID (K*1) correspondent aux valeurs brutes de deux pixels K et K*1 consécutifs. M4 et M9
on sont des constantes positives. Mo, fixée à 500, permet de se déplacer dans le domaine [0,255],
l'information élémentaire étant codée sur un octet. M4 est une valeur empirique destinée à
ces modifier l'échelle des valeurs résultantes selon la problématique envisagée. Les différents essais
re. effectués sur une région test du Ladakh ont conduit à lui attribuer une valeur égale à 20 appli-
cable à l'image de cette région.
En fait, les structures n'apparaissent vraiment que si l'on utilise une fonction f définie à partir
de la fonction g et ne prenant en compte que les valeurs de réflectance dont le gradient en ligne
est positif. Cette fonction qui filtre très sévèrement les valeurs, répond à la formulation suivante :
(a) si ID(K) 2 ID(K-*1)
11ite Log (ID(K) * M4) * M5
f [ID(K), ID(K«1)] = ——————— —£ — 2
ing Log (ID(K+1) + My)
(b) si ID(K) € ID(K-1)
mage f [ID(K), ID(K+1)] = O
A l'issue de ce traitement numérique fait, soit de gauche à droite dans la direction que définissent
les lignes de l'image, soit de haut en bas selon les colonnes, on obtient un document cartographique
dont toutes les ombres ont été supprimées (partant, le relief) et où ne subsistent que les traits
structuraux et le réseau hydrographique. || est alors possible de mettre en évidence, ainsi que
l'illustre l'étude d'une zone test centrée sur les nappes de la Shillakong dans le Ladakh (Inde),
non seulement les directions structurales dessinées par les unités lithologiques, mais encore les
plis qui affectent cet ensemble.
Une étude comparative des résultats obtenus par d'autres méthodes telles que le filtre gradient,
le filtre de Sobel ou le ratioing, montre que la fonction proposée répond mieux au probléme
que pose la reconnaissance dans les zones d'ombre d'une région montagneuse des unités reconnues
dans les zones éclairées, et qu'elle se révéle un bon indicateur des directions structurales dans des
régions de ce type. m
s Bureau de Télédétection, SSC. ORSTOM, 70 route d'Aulnay - F. 93140 BONDY
** Laboratoire de Géologie, SSC. ORSTOM, 70 route d'Aulnay - F. 93140 BONDY
123
I - 7 9757
e. a a I A sand E A —