Cette
ir de
POT
al par
ation
non-
2? 46'
rabie)
ngers
lacée
B.P.,
ns la
a duré
à l'axe
Letea,
ur en
delta,
e une
lation,
), d'où
Chilia
etea et
0,5 %
repère.
ordons
randes
le bras
corghe
ASS de
V (XS)
onnées
Tableau 1. Caractéristiques des données satellitaires
1975 1981 1984 1986 1993
24 juillet 5 juin 28 juin 27 mai 27 avril
Satellite Landsat 2 Landsat 2 Landsat 5 SPOT 1 SPOT 1
Intrument MSS MSS TM HRV 2 HRV 2
Acquisition < multibande >
Bandes 1.2.3.4 1234 1,2,3.4,5,6,7 1,23 123
K/J 195/029 195/029 180/029 98/259 98/259
Résolution (métres) 79 x79 79 x 79 30 x30 20 x 20 20 x 20
Rectification (pixels) +0475 +0,877 +0,813 +0,391 +0,275
Les traitements numériques de ces images ont été effectués référence de HRV (XS).
dans un environnement micro-ordinateur Unitek EISA 486-
DX-50 avec 32 mégaoctets de RAM et 5,35 gigaoctets de
mémoire sur les deux disques durs du systéme. La partie des
traitements — concernant l'application des — méthodes
photogrammétriques numériques a été effectuée avec le logiciel
Desktop Mapping System (DMS) version 3.1. La partie des
traitements de base des images satellitaires a été effectuée avec
le logiciel PCI, version 5.2 EASUPACE et, finalement, pour les
sorties sur l'imprimante couleurs DeskJet 560 C de Hewlett-
Packard on a utilisé le logiciel Idrisi pour Windows.
Les corrections géométriques ont été appliquées sur la base d'un
réseau de points de contróle déterminé à partir des mesures
d'aérotriangulation effectuées sur le restituteur analytique Wild
Aviolyt BC2 ( Barbalata, 1988, 1993).
Une campagne de vérité-terrain (1991-1993), effectué par des
chercheurs roumains a permis de recueillir les données
nécessaires pour dresser une carte thématique et de valider la
plupart des interprétations.
3.2. Méthode
La détection des changements à l'aide de données multibandes
implique des traitements préliminaires, des classifications et des
procédures spécifiques de type post-classification (Allum and
Dreisinger,1987; Hill and Sturm, 1991; Jensen ef al., 1995).
3.2.1. Traitements préliminaires. Pour assurer la comparabilité
des données radiométriques multidates, il a fallu ramener les
données MSS et TM à la résolution spatiale de l'image de
55
À cause des résolutions différentes de ces capteurs, 80 m pour
MSS, 30 m pour TM et 20 m pour HRV (XS) les données ont
été soumises à un suréchantillonnage de 4 x 4 pour MSS et de
1,5 x 1,5 pour TM, ce qui a conduit à l'amélioration des
contrastes et à la sélection de points d'appui. Les 35 points qui
représentent des objets invariants ont servi au calcul des
fonctions de transformation géométrique representées par des
polynomes de 3-éme degré. Ces transformations correspondent
à l'ajustement des images sur-échantillonnées sur l'image de
référence HRV (XS) d'avril 1993. Etant donné que les
coordonnées obtenues aprés les transformations ne
correspondent pas à un pixel précis, il a été nécessaire
d'appliquer une convolution cubique. Cette fonction
d'échantillonnage, qui utilise une approximation bicubique de
16 voisins les plus proches du pixel concerné, permet de
modéliser l'image localement par une surface polynomiale
(Defourny, 1994). À la fin on obtient des plans spectraux qui
vont coïncider parfaitement au niveau de chaque pixel, avec
l’image de référence. La précision de ces transformations varie
selon les images utilisées. Les erreurs moyennes quadratiques
pour chaque scène rectifiée sont présentées par le tableau 1.
Dans le cas de la rectification des quatre images on constate que
les précisions varient de +0,275 à + 0,877 pixels (tableau 1) et
que les précisions observées par exemple pour la rectification
de l'image TM- Landsat (juin 1984) sur l'image HRV-SPOT
(avril 1993) en utilisant 35 points de contróle, varient de + 0,119
à +1 .207 pixels, avec un écart maxima de l'ordre de 20 m.
Des compositions-couleurs infrarouge ont été réalisées en
applicant les traitements suivants:
International Archives of Photogrammetry and Remote Sensing. Vol. XXXI, Part B7. Vienna 1996
