Formation à l'Orfeo ToolBox: Introduction des

Travaux Pratiques
Manuel Grizonnet (CNES), Julien Michel (CNES), Victor Poughon (CNES)

2016

Slides formation OTB 2016 1 / 67

 Introduction

Objectifs de la formation
Objectifs
Présenter l'OTB et mettre en oeuvre les outils disponibles dans les OTB
applications pour l'exploitation, la manipulation des images de télédection
Capacités acquises en n de formation
Mise en ÷uvre de chaine de traitements d'images en utilisant les applications
OTB, incluant entre autre :
Ï Extraction de primitives

Ï Calibration

Ï Classication

Ï Segmentation

Ï Traitements basiques radar

Slides formation OTB 2016 2 / 67

 Introduction

Contenu du kit de formation

Répertoire Data
Ï Contient les données nécessaires à chaque TP, dans des sous-répertoire
séparés
Ï Le sous-répertoire nécessaire est indiqué au début de l'énoncé du TP

Répertoire Documentation
Software Guide Guide de l'API C++ (contient également des descriptions
algorithmiques)
CookBook Guide pour les nons dévelopeurs (contient l'API complète des
applications)
QGIS User Guide Manuel d'utilisation de QGIS

Slides formation OTB 2016 3 / 67

 Introduction

Contenu du kit de formation

Répertoire Guide
Ï Guide de TP
Ï Planches (cette présentation)
Ï Guide d'installation
Ï Questionnaire d'évaluation
Ï Solutions (distribuées à la n)

Slides formation OTB 2016 4 / 67

 Introduction

Origine des données

LandSat-8 plusieurs dates concaténées, dénuagés (gap-lling) avec des

données terrains de référence au format ESRI Shapele (pour le TP
classication),
Pléiades image PHR Bundle PRIMARY Niveau 1A issue de la RTU sur le
©
site OSR MiPy, acquise en novembre 2013 ( CNES (2013),
distribution Airbus DS/ Spot Image),
Sentinel-1 Produit SLC (complexe) SM (strip Map, les plus résolues),
bi-polarisation (HH et HV) sur le sud du lac Constance
(Allemagne).

Slides formation OTB 2016 5 / 67

 Introduction

Accès aux données

LandSat-8 produit de niveau 2A disponible sur le site THEIA

([Link]
Pléiades accès privilégié UIA (institutionnels français) :
[Link]
Sentinel-1 Données gratuites disponible sur ESA
Hub([Link] ou sur le site CNES
PEPS([Link]

Slides formation OTB 2016 6 / 67

 Introduction

Informations pratiques

Repas

Code Wi

Autre ?

Slides formation OTB 2016 7 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Sommaire
1 Introduction

2 TP1 : Généralités
Utiliser Monteverdi et QGIS
Le mécanisme des applications Orfeo ToolBox
Les mécanismes internes de l'Orfeo ToolBox

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

4 TP3 : Classication supervisée pour les séries multi-t

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Slides formation OTB 2016 8 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Objectifs et Données
Objectifs
Ï Savoir visualiser une image et régler son rendu dans Monteverdi,
Ï Savoir visualiser une pile de données dans Monteverdi,
Ï Savoir visualiser une image et régler son rendu dans QGIS,
Ï Savoir visualiser une données vecteur dans QGIS,
Ï Savoir visualiser une pile de données dans QGIS,
Ï Échanger des trucs et astuces.

Données
Les données utilisées se trouvent dans les répertoires suivant :
Ï Data/preprocessing (extrait Pléiades)

Ï Data/classification/Extract16bits/ (séries temporelle LandSat8)

Ï Data/classification/training/ (chier vecteur shp)

Slides formation OTB 2016 9 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Déroulement

1 Visualiser une image dans Monteverdi

2 Visualiser une pile d'images dans Monteverdi
3 Visualiser une image dans QGIS
4 Visualiser une données vecteur dans QGIS

Slides formation OTB 2016 10 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Introduction

Ï Plusieurs logiciels qui partagent des fonctionnalités communes

Ï Mais des philosophies parfois diérentes (Monteverdi orienté et centré
"Image", QGIS plateforme SIG)
Ï Plusieurs façon de faire la même opération dans ces logiciels

Slides formation OTB 2016 11 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Monteverdi

Ï A l'origine c'est un outil simple de démonstration du potentiel de la librairie

OTB
Ï Devenu un outil intégré de visualisation qui embarque toutes les
OTB-Applications
Ï Interface graphique qui utilise la bibliothèque Qt
Ï le moteur de rendu utilise la librairie Ice (OpenGL+OTB)
Ï Eets tirant partie des capacités des cartes graphiques actuelles (shader) pour
eectuer des rendus à la volée (contraste local, transparence, gradient. . . )

Slides formation OTB 2016 12 / 67

TP1 : Généralités Utiliser Monteverdi et QGIS

Monteverdi

Slides formation OTB 2016 13 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Monteverdi

Ï Fenêtre principale découpée en plusieurs parties :

Ï Main menu (ouvrir image, fermer, préférences. . . )
Ï Top toolbar (raccourcis)
Ï Image displaying (vue principale)
Ï Right side dock (composition colorée, quicklook. . . )
Ï Stack layer (pile d'images)
Ï Voir A brief tour of Monteverdi

Slides formation OTB 2016 14 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Quantum GIS

Ï système d'information géographique (SIG) : système d'information capable

d'organiser et de présenter des données alphanumériques spatialement
référencées
Ï Les principales fonctionnalités d'un SIG sont les suivantes :
Ï Achage ou production de cartes
Ï Acquisition des données
Ï Analyse des données pour créer une nouvelle information
Ï Abstraction ou représentation des éléments choisis
Ï Archivage
Ï Beaucoup de support de formation, tutoriels disponibles (par exemple en
français : Intro ENTE Aix - ENS - Licence ouverte)

Slides formation OTB 2016 15 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Quantum GIS processing (a partir

de la version 2.0)

Ï Plateforme de traitements permettant d'intéragir avec GRASS, SAGA-GIS, R

et l'OTB !
Ï Historique : intégration dans QGIS du projet sextante (développée par Victor
Olaya)
Ï Accès à toutes les applications OTB (certaines sont repackagées avec une
interface simpliée)
Ï Voir QGIS processing documentation

Slides formation OTB 2016 16 / 67

TP1 : Généralités Utiliser Monteverdi et QGIS

Quantum GIS

Slides formation OTB 2016 17 / 67

 TP1 : Généralités Utiliser Monteverdi et QGIS

Principe

Ï Découverte du jeu de données et des logiciels

Ï Monteverdi :
Ï Savoir visualiser une image et régler son rendu
Ï Savoir visualiser une pile de données
Ï Explorer les "Eets" disponible
Ï QGIS :
Ï Savoir visualiser une image et régler son rendu
Ï Savoir visualiser une données vecteur
Ï Savoir visualiser une pile de données
Ï Travail 15 minutes en binome et en aveugle
Ï Tour de table (30 à 45 minutes), chaque binome passe devant les autres et
propose une astuce ou pose une question

Slides formation OTB 2016 18 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Sommaire
1 Introduction

2 TP1 : Généralités
Utiliser Monteverdi et QGIS
Le mécanisme des applications Orfeo ToolBox
Les mécanismes internes de l'Orfeo ToolBox

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

4 TP3 : Classication supervisée pour les séries multi-t

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Slides formation OTB 2016 19 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Objectifs et données

Objectifs
Ï Savoir rechercher une application dans la liste des applications disponibles
Ï Savoir paramétrer une application
Ï Savoir trouver la documentation d'une application
Ï Savoir utiliser plusieurs applications classiques

Données
Les données se trouvent dans le répertoire Data/stegano/.

Slides formation OTB 2016 20 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Déroulement

Le TP consiste à décoder successivement 6 messages encodés par stéganographie

dans les images, en utilisant les applications OTB.

Slides formation OTB 2016 21 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Rapide intro aux OTB-Applications

Ï Orfeo ToolBox est une bibliothèque orant des fonctionnalités de
télédétection
Ï Souvent nécessaire de combinant plusieurs fonctions (ltres) OTB dans une
chaîne de traitements
Ï Besoin interface de plus haut niveau (paramètres, les données d'entrée et de
sortie, les logs. . . )
Ï Chaque application correspond à une fonctionnalité de "haut niveau"
(segmentation, orthorectication, segmentation, calculatrice sur les
bandes. . . )
Ï Correspond à une classe C++ (otb : :Application) → librairie
Ï Mécanisme de plugin
Ï Possibilité de développer et distribuer ces applications homemade

Slides formation OTB 2016 22 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Rapide intro aux OTB-Applications

Ï En ligne de commande → otbcli_Orthorectication

Ï Dans une interface graphique → otbgui_Orthorectication
Ï En Python :
import otbApplication
orth=[Link]("OrthoRectification")

Ï Intégré dans QGIS → Menu processing

Ï Intégré dans Monteverdi

Slides formation OTB 2016 23 / 67

TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

GUI

Slides formation OTB 2016 24 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Paramétrage

Ï Acher l'aide de l'application d'OrthoRectication (bon exemple)

Ï Notion de groupe de paramètres
Ï Notion de paramètre obligatoire
Ï Dépendance entre paramètre (dynamique)
Ï Valeur par défaut
Ï Type des images en sortie (par défaut ottant)

Slides formation OTB 2016 25 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Autre Exemple : BandMath

Ï Calculatrice raster
Ï Utilise la librairie MuParser
Ï Prend en entrée une liste d'images (paramètre -il)
Ï Codage des expressions mathématiques :
Ï imXbY :
Ï X : numéro de l'image dans la liste (de 1 à N)
Ï Y : numéro de la bande dans l'image X (de 1 à N)

Ï Liste des opérateurs mathématiques disponibles

Ï Exemples :
Ï *-exp "(im1b4 - im1b1) / (im1b4 - im1b1))"
Ï *-exp "acos((123*im1b1+265*im1b2+652*im1b3)
/ (sqrt(123*123+265*265+652*652)
*sqrt(im1b1*im1b1+im1b2*im1b2+im1b3*im1b3)))"
Ï *-exp "im1b1>0.5 ?255 :0"*

Slides formation OTB 2016 26 / 67

 TP1 : Généralités Le mécanisme des applications Orfeo ToolBox

Liens utiles

Ï Introduction aux OTB applications

Ï Documentation de référence des applications
Ï Exemples utilisation des OTB applications
Ï QGIS processing
Ï Aide d'une application : otbcli_XXX -help

Slides formation OTB 2016 27 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Sommaire
1 Introduction

2 TP1 : Généralités
Utiliser Monteverdi et QGIS
Le mécanisme des applications Orfeo ToolBox
Les mécanismes internes de l'Orfeo ToolBox

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

4 TP3 : Classication supervisée pour les séries multi-t

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Slides formation OTB 2016 28 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Objectifs et données

Objectifs
Ï Comprendre les mécanismes transparents de l'Orfeo ToolBox
Ï Savoir inuencer la manière dont l'Orfeo ToolBox exécute un traitement
Ï Savoir ou trouver des informations complémentaires

Données
Les données se trouvent dans le répertoire Data/internals/.

Slides formation OTB 2016 29 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Déroulement

Le déroulement du TP est le suivant :

1 Encodage des images
2 Fichiers geom
3 Noms de chiers étendus
4 Streaming
5 Multi-threading

Slides formation OTB 2016 30 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Encodage des images

type plage de valeurs nombre de bits

uint8 [0,255] 8 bits
int16 [-32 767, +32 767 ] 16 bits
uint16 [0, 65 535] 16 bits
int32 [-2 147 483 647, -2 147 483 647] 32 bits
uint32 [0, 4 294 967 294] 32 bits
oat Ö
[ -3.402823 1038,3.402823 1038 ] Ö 32 bits
double [ -10308 , 10308 ] 64 bits

Slides formation OTB 2016 31 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Fichiers geom
[...]
samp_num_coeff_15: -8.69402623737171e-06
samp_num_coeff_16: -2.52010136133467e-09
samp_num_coeff_17: -5.70277370040739e-07
samp_num_coeff_18: -2.67844954240191e-07
samp_num_coeff_19: -7.80920276666093e-09
samp_off: 19999
samp_scale: 19999.5
sensor: PHR 1A
[...]
support_data.image_date: 2012-11-15T[Link].4Z
support_data.image_id: 602631101-001
support_data.image_size: 38187 38890
support_data.instrument: PHR
support_data.instrument_index: 1A
support_data.line_period: 0.0735
[...]
Slides formation OTB 2016 32 / 67
 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Noms de chiers étendus

Syntaxe comprise par tout code OTB lisant ou écrivant des images
"[Link]?&geom=[Link]"
"[Link]?&gdal:co:TILED=yes&streaming:type=none"

Slides formation OTB 2016 33 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Streaming et multi-threading (1/2)

Slides formation OTB 2016 34 / 67

 TP1 : Généralités Les mécanismes internes de l'Orfeo ToolBox

Streaming et multi-threading (2/2)

Slides formation OTB 2016 35 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Pré-traitements de l'imagerie THR optique

Sommaire

1 Introduction

2 TP1 : Généralités

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

Pré-traitements de l'imagerie THR optique
Segmentation et export vers un SIG

4 TP3 : Classication supervisée pour les séries multi-t

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Slides formation OTB 2016 36 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Pré-traitements de l'imagerie THR optique

Objectifs et Données

Objectifs
Ï Savoir réaliser une calibration optique
Ï Savoir réaliser une fusion (pan-sharpening)
Ï Savoir réaliser une ortho-rectication

Données
Les données se trouvent dans le répertoire Data/preprocessing/. Les
sous-répertoires SRTM et Geoid sont également utilisés.

Slides formation OTB 2016 37 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Pré-traitements de l'imagerie THR optique

Déroulement

Le déroulement du TP est le suivant :

1 Corrections atmosphériques
2 Fusion P+XS
3 Ortho-rectication

Slides formation OTB 2016 38 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Pré-traitements de l'imagerie THR optique

Calibration radiométrique
Correction atmosphérique

Slides formation OTB 2016 39 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Pré-traitements de l'imagerie THR optique

Fusion (pan-sharpening)

Pourquoi le pansharpening ?
Ï La plupart des capteurs THR sont constitués de deux voix :
Ï Une voie panchromatique (une bande spectrale large) très résolue
Ï Une voie multi-spectrale (plusieurs bandes spectrales plus étroites) moins
résolue (habituellement 4 fois moins)
Ï Pansharpening = image de synthèse avec la résolution du panchro et la
couleur du multi-spectral

Principe :
1 Superposition ne des voies P et XS
2 Application d'un algorithme de fusion

Slides formation OTB 2016 40 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Pré-traitements de l'imagerie THR optique

Orthorectication

Bundle-block Homologous
Adjustement Points

Sensor Fine Map

Model Regis- Projec-
tration tion

Input series Geo-referenced Series Registered Series

DEM

Slides formation OTB 2016 41 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Sommaire

1 Introduction

2 TP1 : Généralités

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

Pré-traitements de l'imagerie THR optique
Segmentation et export vers un SIG

4 TP3 : Classication supervisée pour les séries multi-t

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Slides formation OTB 2016 42 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Objectifs et Données

Objectifs
Ï Connaître les étapes pour réaliser une segmentation
Ï Savoir optimiser les paramètres de la segmentation
Ï Savoir exporter la segmentation vers un logiciel SIG

Données
Les données pour cet exercice se trouvent dans le répertoire Data/segmentation.

Slides formation OTB 2016 43 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Déroulement

Le déroulement du TP est le suivant :

1 Lissage de l'image par l'algorithme MeanShift
2 Segmentation
3 Traitement des petites régions
4 Vectorisation
5 Filtrage des polygones dans QGIS

Slides formation OTB 2016 44 / 67

TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Algorithme MeanShift (1/4)

Recherche dans la fenêtre des pixels proches spectralement

Slides formation OTB 2016 45 / 67

TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Algorithme MeanShift (2/4)

Moyennes spatiale et spectrale des pixels localisés dans la fenêtre

Slides formation OTB 2016 46 / 67

 TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Algorithme MeanShift (3/4)

Déplacement du pixel vers le barycentre et aectation de la moyenne spectrale

Slides formation OTB 2016 47 / 67

TP2 : Imagerie THR optique, des pré-traitements au SIG Segmentation et export vers un SIG

Algorithme MeanShift (4/4)

Retour à l'étape 1 jusqu'à convergence

Slides formation OTB 2016 48 / 67

 TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Sommaire

1 Introduction

2 TP1 : Généralités

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

4 TP3 : Classication supervisée pour les séries multi-t

Classication supervisée pour les séries multi-temporelles

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Slides formation OTB 2016 49 / 67

 TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Objectifs et Données
Objectifs
Les objectifs sont les suivants :
Ï Connaître les diérentes applications constituant la procédure de

classication supervisée
Ï Utiliser diérents algorithmes pour l'apprentissage

Ï Savoir mesurer les performances de la classication

Ï Connaître les post-traitements applicables à une classication

Données
Les données sont disponibles dans le répertoire Data/classification, avec les
sous-répertoires suivants :
Ï Extract16bits contient la série multi-temporelle LandSat8,

Ï training contient la donnée d'apprentissage au format shp,

Ï testing contient la donnée de validation au format shp.

Slides formation OTB 2016 50 / 67

 TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Déroulement

les étapes de l'exercice sont les suivantes :

1 Introduction aux données landsat8
2 Classication mono-date
3 Classication multi-date
4 Classication avec prol de NDVI
5 Post-traitements de la classication

Slides formation OTB 2016 51 / 67

 TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Présentation des données LandSat8

Résolution spatiale : 30 mètres
Dates
Bandes
2014-03-09
2014-04-01 0 Coastal aerosol
2014-04-17 1 Blue
2014-05-28 2 Green
2014-06-20 3 Red
2014-07-31 4 Near Infrared (NIR)
2014-09-01 5 SWIR 1
2014-10-03 6 SWIR 2
2014-10-26

Slides formation OTB 2016 52 / 67

TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Présentation des données de

référence
Code Nom #polygones
11 Éte 7898
12 Hiver 8171
31 Foret feuilles caduques 867
32 Foret feuilles persistantes 125
34 Pelouses 45
36 Lande ligneuse 386
41 Bati 4719
51 Eau 1280
211 Prairie 5647
221 Verger 204
222 Vigne 559

Slides formation OTB 2016 53 / 67

TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Classication supervisée

Slides formation OTB 2016 54 / 67

TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Algorithme SVM

Slides formation OTB 2016 55 / 67

 TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Algorithme RF
Ensemble d'arbres de décision aléatoires
Apprentissage
1 Séparer le jeu d'apprentissage en k ensembles Sk aléatoires
2 Pour chaque Sk choisir aléatoirement Fk primitives
3 Construire un arbre de décision récursivement, pour chaque noeud :
1 Choisir f ∈ Fk et le seuil tk qui sépare l'ensemble restant en 2 parties les plus
pures
2 Arrêter quand l'ensemble restant devient trop petit

Décision
Vote majoritaire de tous les arbres aléatoires

Slides formation OTB 2016 56 / 67

 TP3 : Classication supervisée pour les séries multi-t Classication supervisée pour les séries multi-temporelles

Matrice de confusion

Préd. 1 Préd. 2 Préd. 3

Réf. 1 Vrais pos. 1
Réf. 2 Vrais pos. 2
Réf. 3 Vrais pos. 3

Ï precision = P VPi
pred .i
VPi
Ï rappel = T P ref .i
P
VPi
Ï Accuracy = Total
Accuracy −chance
Ï Kappa = 1−chance

Slides formation OTB 2016 57 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Sommaire

1 Introduction

2 TP1 : Généralités

3 TP2 : Imagerie THR optique, des pré-traitements au SIG

4 TP3 : Classication supervisée pour les séries multi-t

5 TP4 : Traitements SAR pour l'imagerie Sentinel 1

Introduction au traitements des images RSO

Slides formation OTB 2016 58 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Objectifs et Données

Objectifs
Ï Savoir manipuler des images
Ï Savoir réaliser une calibration radiométrique
Ï Savoir réaliser une orthorectication
Ï Savoir utiliser le ltrage du speckle
Ï Savoir réaliser une analyse polarimétrique simple
Ï Savoir réaliser une extraction de primitives simple

Données
Les données pour cet exercice se trouvent dans le répertoire Data/sar.

Slides formation OTB 2016 59 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Déroulement

1 Introduction à l'imagerie RSO

2 Calibration radiométrique
3 Calibration géométrique
4 Filtrage du speckle
5 Polarimétrie
6 Extraction d'information

Slides formation OTB 2016 60 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Diérence Optique/SAR

(source CNES : IMAGERIE SPATIALE Des principes d'acquisition au traitement

des images optiques pour l'observation de la Terre)

Slides formation OTB 2016 61 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Introduction SAR

Ï RSO : Radar à Synthèse d'Ouverture (SAR en anglais)

Ï Impulsions hyperfréquences et enregistrement des échos
Ï Quasi-transparence de l'atmosphère
Ï Acquisition de jour comme de nuit
Ï Cours en ligne ENSG

Slides formation OTB 2016 62 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Principe de fonctionnement

Ï 1 antenne réalise des impulsions micro-ondes dans le plan perpendiculaire à la

trajectoire du porteur
Ï Écho recueilli par la même antenne
Ï Réecteur sont localisés en fonction de leur éloignement (distance)
Ï Répétition des échos réalise la deuxième dimension (azimuth)
Ï Signal reçu par le radar est caractérisé par son amplitude et sa phase

Slides formation OTB 2016 63 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Calibration radiométrique

Ï Comportement de la réexion dépend des caractéristiques de la surface

Ï Notamment de sa rugosité, de l'humidité des sols
Ï Comme pour l'optique, la calibration permet de transformer des comptes
numériques (CN) en grandeur physique
Ï Permet d'inter-comparer des images radars provenant de capteurs diérents
ou de capteurs identiques utilisés dans diérents modes

Slides formation OTB 2016 64 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Correction géométrique

Ï Spatialiser l'image
Ï Associer une coordonnée géographique au pixel
Ï le radar mesure la distance des objets obliquement
Ï Resolutions
Ï azimuthResolution = L∗H ∗λ
cos(θ )
c ∗prf
Ï rangeResolution = 2∗sin(θ)
Ï Il existe distorsions géométriques dues au déplacement du relief
(foreshortening et layover)

Slides formation OTB 2016 65 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Speckle

Ï Images RSO fortement aectées par le chatoiement (speckle)

Ï Bruit très fort
Ï Eet multiplicatif !
Ï Plusieurs méthodes pour réduire ce bruit
Ï Atténuer ce bruit en conservant le maximum de détails
Ï Filtrage permet d'améliorer énormément la qualité des images et leur
interprétation

Slides formation OTB 2016 66 / 67

 TP4 : Traitements SAR pour l'imagerie Sentinel 1 Introduction au traitements des images RSO

Pour aller plus loin

Ï S1 ToolBox (SNAP)
Ï Polarimétrie : ESA PolSARPro

Slides formation OTB 2016 67 / 67