0% ont trouvé ce document utile (0 vote)
56 vues167 pages

Imagerie Radar : Techniques et Applications

Transféré par

Sahbi Ouali
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
56 vues167 pages

Imagerie Radar : Techniques et Applications

Transféré par

Sahbi Ouali
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Imagerie Radar

PPMD

Sébastien Giordano

Univ. Paris-Est, LASTIG MATIS, IGN, ENSG, F-94160 Saint-Mande, France

vendredi 30 novembre 2018


Introduction

Géométrie
des images
radar

Radiométrie
1 Introduction 3 Radiométrie
Radargrammétrie
et Radarcli-
nométrie
2 Géométrie des images radar 4 Radargrammétrie et
Interférommétrie
Radarclinométrie
polarimetrie

Sébastien Giordano  30-11-2018 2/2 Radar PPMD


IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
Introduction
polarimetrie

Sébastien Giordano  30-11-2018 3/3 Radar PPMD


Télédétection Active

RADAR vs. LIDAR


RADAR : RAdio Detection And Ranging utilise les ondes radio
IIntroduction
(micro-ondes ou hyperfréquences)
Géométrie
des images LIDAR : LIght Detection And Ranging utilise les ondes du visible (vert)
radar
ou du proche infrarouge (lumière)
Radiométrie

Radargrammétrie
et Radarcli-
nométrie
Radar : système actif
Interférommétrie Systèmes de surveillance des avions et des bateaux
polarimetrie
Radar altimètre : (aéronautique : aide au pilotage ; océanographie)
Systèmes imageurs (radars à ouverture réelle ou à ouverture synthétique)

Dans ce cours on s’intéresse aux systèmes radar imageurs (à ouverture


réelle ou à ouverture synthétique).

1
1
[Link] électromagné[Link]
Sébastien Giordano  30-11-2018 3/3 Radar PPMD
Bandes Radar et Systèmes Radar

Bandes radar :

IIntroduction Bande ∆λ Exemples d’application


Géométrie
des images VHF 1m-6m pénétration couvert forestier, sol
radar
P-Band 30 cm - 1 m humidité du sol, biomasse
Radiométrie
L-Band 15 cm - 30 cm foresterie, humidité du sol
Radargrammétrie
et Radarcli- C-Band 3.75 cm - 7.5 cm ocean, agriculture
nométrie
X-Band 2.5 cm - 3.75 cm ocean, agriculture, urbain
Interférommétrie
Ku-Band 1.67 cm - 2.5 cm glaciologie (couvert neigeux)
polarimetrie
Ka-Band 0.75 cm - 1.11 cm haute résolution radar

Radars Imageurs → monofréquence

Satellite λ
Terrasar-X 3.1 cm
Radarsat-2 5.5 cm
ALOS-PALSAR-2 24 cm
Biomass (2020) 69 cm

Sébastien Giordano  30-11-2018 4/4 Radar PPMD


Radar actif : système de surveillance

IIntroduction déterminer la position ainsi que la vitesse d’objets tels que les avions,
Géométrie les bateaux, ou la pluie,
des images
radar ondes envoyées par l’émetteur sont réfléchies par la cible, distance est
Radiométrie obtenue grâce au temps aller/retour du signal,
Radargrammétrie
et Radarcli-
direction grâce à la position angulaire de l’antenne et décalage de
nométrie fréquence du signal de retour généré selon l’effet Doppler,
Interférommétrie

polarimetrie

Ecran radar bateau Radar dual pol bande X - ENPC

Sébastien Giordano  30-11-2018 5/5 Radar PPMD


Radar actif : Radar altimètre
déterminer avec une très grande précision la distance entre le satellite et
la surface d’intérêt, afin de déduire la hauteur (ou la variation de hauteur)
IIntroduction
mesurer la hauteur des vagues qui permet d’estimer la force du vent
Géométrie Exemples: ERS-1 (08/1991 - 06/1996), Topex/Poseidon (09/1992 -
des images
radar
10/2005), Jason-2 (07/2008- ), SRAL (Sar Radar Altimeter) sur Sentinel
Radiométrie
3 bi-fréquence bande Ku et en bande C.
Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Sea level trend from 21 years of radar altimeter data (1993 to 2014) derived
from retracked Topex, Jason-1 and Jason-2 data.
Sébastien Giordano  30-11-2018 6/6 Radar PPMD
Système passif

Système passif
Radiomètre micro-onde : contenu en vapeur d’eau et eau liquide de
IIntroduction
l’atmosphère. Exemple : SMOS (Soil Moisture and Ocean Salinity),
Géométrie
des images radiomètre interféromètre en bande L, programme commun d’observation
radar
de la Terre ESA, CNES et CDTI.
Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Soil moisture, November 2010 and 20113

3
[Link] moisture November 2010 and 2011
Sébastien Giordano  30-11-2018 7/7 Radar PPMD
Satellites SAR de première génération

Satellites SAR de première génération

IIntroduction
SEASAT (1978) : premier satellite SAR : Radar à Synthèse d’Ouverture
Géométrie ERS (1991 - ) : ESA. ERS-2 de 1995 à 2000 sur la même orbite que
des images
radar ERS-1, bande C
Radiométrie JERS (1992 - 1998) : Japanese Earth Ressource Satellite, bande L
Radargrammétrie
et Radarcli- Radarsat (1995 - ) : bande C, Agence Spatiale Canadienne
nométrie
ENVISAT (2002 - 2012) : ESA, 10 instruments dont SAR en bande C
Interférommétrie

polarimetrie

Utilisation de la navette spatiale


Navette endeavour du 11 au 22 février 2000, mission SRTM (Shuttle Radar
Topography Mission, bande C, MNT par interférométrie.

Sébastien Giordano  30-11-2018 8/8 Radar PPMD


Satellites SAR de seconde génération

Années 90 Depuis 2006

IIntroduction
radars SAR décamétriques amélioration de la résolution spatiale
Géométrie
en bande C et L bande X, C, L et bientôt P
des images polarisation simple ou double polarimétrie
radar

Radiométrie

Radargrammétrie Exemples de satellites SAR de seconde génération


et Radarcli-
nométrie ALOS-PALSAR-2, SAR en bande L, satellite japonais ALOS-PALSAR-1
Interférommétrie (2006-2011):
polarimetrie Mode spotlight, single pol, 0.625m/1.25m, 25km×25km
Mode Stripmap 2.5m/5.0m single or dual, Pol 55km×70km
Full (Quad.) Polarimetry, full polar, 3.125m/6.25m, 50km×70km
Radarsat-2 (2007 - ) : Agence Spatiale Canadienne, ultra-fine
3m × 3m, full polar 11m × 9m
Terrasar-x (2007 - ) EADS et DLR, bande X polarimétrique, résolution
métrique, Tandem-X
Cosmoskymed : Agence spatiale Italienne ASI, mode spotlight mieux
que 1m, possibilité full polar
Sentinel-1 (2014(A)/2016(B)), bande C
Mode interférométrique : 250 km, résolution est de 5m×20m
Sébastien Giordano  30-11-2018 9/9 Radar PPMD
TerraSAR-X : mode Staring SpotLight 25 cm, λ=3,1cm

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  10 / 10 Radar PPMD


TerraSAR-X, mode StripMap 6m VV/HH, Sanchagang, Chine

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  11 / 11 Radar PPMD


COSMO-Skymed SCS B product 4

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

4
COSMO-SkyMed Image ASI
c (year of acquisition). All rights reserved
Sébastien Giordano  30-11-2018  12 / 12 Radar PPMD
Radarsat2 - full polar HH, HV, VV - bande C

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Forcalquier
Sébastien Giordano  30-11-2018  13 / 13 Radar PPMD
ALOS/PALSAR - full polar HH, HV, VV - bande L

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Fontainebleau
Sébastien Giordano  30-11-2018  14 / 14 Radar PPMD
Satellite Radar : Bilan

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  15 / 15 Radar PPMD


TerraSAR-X et Tandem-X 5

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie (Loading)

5
[Link]
Sébastien Giordano  30-11-2018  16 / 16 Radar PPMD
Global 25m Resolution PALSAR-2/PALSAR Mosaic and Forest/N
Forest Map

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

6
[Link] fnf/fnf [Link]
Sébastien Giordano  30-11-2018  17 / 17 Radar PPMD
WorldDEM (Tandem-X)

IIntroduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  18 / 18 Radar PPMD


Radar aéroporté

IIntroduction
SETHI
Géométrie ONERA : SETHI (Système Expérimental de Télédétection
des images
radar Hyperfréquence Imageur)
Radiométrie radars en bandes P, L et X ; interférométrie et polarimétrie complète
Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  19 / 19 Radar PPMD


Radar : Amplitude et Phase

Radar : Amplitude et Phase par canal de polarisation

IIntroduction SXY = |SXY | e jφXY (1)


Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Amplitude |SVV | Figure: Amplitude φVV

Sébastien Giordano  30-11-2018  20 / 20 Radar PPMD


Avantage des techniques microondes

IIntroduction

Géométrie
des images
radar Par rapport aux capteurs optiques, l’intérêt de l’emploi des radars imageurs est
Radiométrie lié à :
Radargrammétrie la faible sensibilité aux conditions atmosphériques ; indépendance de
et Radarcli-
nométrie la radiation solaire permettant l’obtention d’images de jour, de nuit et à
Interférommétrie travers 100% de couvert nuageux ;
polarimetrie
la possibilité de contrôle de la radiation électromagnétique émise :
puissance, fréquence, polarisation
l’accès à des paramètres physiques différents de ceux détectés dans le
visible et l’infrarouge

Sébastien Giordano  30-11-2018  21 / 21 Radar PPMD


Systèmes radar actif

Principe : envoyer vers la surface terrestre un rayonnement électromagnétique


IIntroduction puis mesurer la quantité de rayonnement électromagnétique renvoyé .
Géométrie
des images
radar
Caractéristiques de l’onde électromagnétique
Radiométrie direction de propagation
Radargrammétrie
et Radarcli-
amplitude
nométrie
longueur d’onde
Interférommétrie

polarimetrie
polarisation
phase

Sébastien Giordano  30-11-2018  22 / 22 Radar PPMD


Principe de base du radar

IIntroduction Transmission et réception d’impulsions dans un faisceau étroit des bandes cen-
Géométrie timétriques.
des images
radar

Radiométrie
Les signaux de retour sont enregistrés :
Radargrammétrie force relative
et Radarcli-
nométrie temps mis par l’impulsion pour atteindre l’objet et en revenir
Interférommétrie
phase
polarimetrie

Sébastien Giordano  30-11-2018  23 / 23 Radar PPMD


Définitions

IIntroduction

Géométrie
des images
radar Systèmes radar de deux types :
Radiométrie Radar Bistatique : émetteur et le récepteur séparés.
Radargrammétrie
et Radarcli- Radar Monostatique : émetteur et le récepteur au même endroit
nométrie

Interférommétrie
Radar Imageurs :
polarimetrie
Radar à antenne réelle RAR (Real Aperture Radar ) ou SLAR Side
Looking Airbone Radar
Radars à antenne synthétique SAR (Synthetic Aperture Radar )

Sébastien Giordano  30-11-2018  24 / 24 Radar PPMD


Introduction

IGéométrie
des images
radar
Résolution
radiale
Résolution
azimutale
Distorsions
géométriques Géométrie des images radar
Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  25 / 25 Radar PPMD


Ambiguités de distance
Traitement radar basé sur le temps de parcours des impulsions
Visée verticale : ambiguı̈tés en temps et par conséquent en distance
Introduction
Ambiguı̈tés
IGéométrie
des images
radar AB = AB 0 : l’onde rétrodiffusée par B arrive à l’antenne A en même temps
Résolution
radiale
que B 0
Résolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Ambiguı̈tés en distances provoquées par l’observation verticale (source :


[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991])
Sébastien Giordano  30-11-2018  25 / 25 Radar PPMD
Visée oblique
Utilisation d’une visée oblique pour minimiser les ambiguı̈tés
Ambiguı̈tés
Introduction
Point A : Near Range, Point B : Far Range
IGéométrie
des images
radar
Résolution
radiale
Résolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Détermination des distances obliques (source :


[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991])
Sébastien Giordano  30-11-2018  26 / 26 Radar PPMD
Résolutions des images radar

Résolution spatiale radiale (across track, range resolution)


Résolution spatiale azimutale (along-track resolution)
Introduction
Cas du SLAR et du SAR
IGéométrie
des images
radar
Résolution
radiale
Résolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

7
[Link]
Sébastien Giordano  30-11-2018  27 / 27 Radar PPMD
SLAR : Résolution radiale

Résolution radiale
Résolution radiale déterminée par la durée de l’impulsion
Introduction

Géométrie
des images
radar
IRésolution
radiale
Résolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Durée de l’impulsion radar (source :


[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991])

Sébastien Giordano  30-11-2018  28 / 28 Radar PPMD


SLAR : Résolution radiale

Résolution radiale
Introduction Distance oblique entre
Géométrie A et B = L2
des images
radar
PA rétrodiffusion de P
IRésolution
radiale par A
Résolution
azimutale
PB rétrodiffusion de P
Distorsions
géométriques par B
Radiométrie
superposition de PA et
Radargrammétrie
et Radarcli- PB = confondus dans
nométrie
1 cellule de résolution
Interférommétrie

polarimetrie
Si
distanceoblique > L2
= cellules de
résolution différentes

Figure: Effet de la durée d’impulsion émise sur la


résolution radiale (source :
[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991])

Sébastien Giordano  30-11-2018  29 / 29 Radar PPMD


Résolution radiale

Résolution radiale oblique


L cτ
Résolution radiale oblique (slant range) : δr = 2
= 2
Introduction
cτ cτ
Géométrie Résolution radiale sur le terrain (ground range) : δr = 2 cos(α)
= 2 sin(θ)
des images
radar
IRésolution
avec :
radiale
Résolution c : vitesse de la lumière
azimutale
Distorsions τ : durée de l’impulsion
géométriques

Radiométrie α : angle de dépression


Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
Remarques sur la résolution
polarimetrie
radiale (× 2)
Comment améliorer la résolution
radiale ?

Sébastien Giordano  30-11-2018  30 / 30 Radar PPMD


Résolution radiale

Introduction

Géométrie La résolution radiale:


des images
radar
IRésolution
est indépendante de la distance entre l’antenne et la cible pour un α
radiale donné,
Résolution
azimutale
Distorsions
varie à travers la fauchée en fonction de α.
géométriques
améliorer la résolution radiale : réduire les impulsions radar (et donc
Radiométrie
augmenter la puissance)
Radargrammétrie
et Radarcli-
nométrie
difficulté d’émettre de courtes impulsions et de fortes énergies
Interférommétrie

polarimetrie
Solution : la compression d’impulsion
émission d’un signal fortement modulé en fréquence : chirp

Sébastien Giordano  30-11-2018  31 / 31 Radar PPMD


Chirp

Introduction

Géométrie
des images
radar
IRésolution
radiale
Résolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  32 / 32 Radar PPMD


Résolution radiale : Chirp

filtrage adapté : auto-corrélation du signal reçu avec le signal émis

Introduction

Géométrie
des images
radar
IRésolution
radiale
Résolution
azimutale
Distorsions
géométriques

Radiométrie
8
Radargrammétrie
et Radarcli-
avant filtrage adapté après filtrage adapté
nométrie

Interférommétrie

polarimetrie

8
[Link] d’impulsion
Sébastien Giordano  30-11-2018  33 / 33 Radar PPMD
Modulation de fréquence : Chirp

but : impulsion longue (énergie) et bonne résolution radiale ?


Chirp : émission d’un signal modulé en fréquence
Introduction 0
1
après compression d’impulsion... τ = ∆f avec ∆f largeur de la bande de
Géométrie
des images modulation autour de la fréquence centrale (our porteuse), largeur du pic
radar
IRésolution
d’autocorrélation d’auto-corrélation.
compression d’impulsion our que ττ0 inférieur à 1, de l’ordre de 20 à 30.
radiale
Résolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Résolution radiale avec compression d’impulsions


c
δr = (2)
2∆f

Sébastien Giordano  30-11-2018  34 / 34 Radar PPMD


Résolution azimutale

Résolution azimutale
Introduction Capacité de discriminer deux points voisins
Géométrie alignés parallèlement à la trajectoire de
des images
radar l’antenne
Résolution
radiale
IRésolution
Pour un radar à ouverture réelle la résolution azimutale dépend
azimutale
Distorsions
uniquement de la largeur du faisceau (A) et donc des dimensions de
géométriques
l’antenne.
Radiométrie
Pour une longueur d’onde donnée λ l’angle d’ouverture du faisceau β
Radargrammétrie
et Radarcli-
dépend de la dimension effective D de l’antenne β = Dλ .
nométrie
Les cibles 1 et 2 sont discernables mais les cibles 3 et 4 plus éloignées ne
Interférommétrie
le sont pas.
polarimetrie

Résolution azimutale
λR
δaz = R.β = D
λ : longueur d’onde du radar
D : dimension de l’antenne dans la direction de son déplacement
R : distance oblique

Sébastien Giordano  30-11-2018  35 / 35 Radar PPMD


Principe de la synthèse d’ouverture

Introduction

Géométrie SLAR sur une plateforme satellite :


des images
radar SLAR sur une plateforme satellite : trop faible résolution azimutale
Résolution
radiale grande taille antenne physiquement pas réalisable
IRésolution
azimutale
Distorsions
Exemple ERS (orbite 900km): tâche au sol 5km along-track. Pour
géométriques obtenir une résolution décamétrique il faudrait une antenne de plusieurs
Radiométrie kilomètres de long.
Radargrammétrie
et Radarcli-
nométrie Principe
Interférommétrie
La compression d’impulsions permet de simuler une antenne très longue en
polarimetrie
utilisant une antenne physique relativement courte. Ce traitement est rendu
possible par l’Effet Doppler.

Sébastien Giordano  30-11-2018  36 / 36 Radar PPMD


Effet Doppler

Effet Doppler : décalage de fréquence d’une onde entre la mesure à l’émission


et la mesure à la réception lorsque la distance entre l’émetteur et le récepteur
Introduction
varie au cours du temps
Géométrie
des images
radar
Résolution
radiale
IRésolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie (Loading)
Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  37 / 37 Radar PPMD


Effet Doppler
Antenne synthétique : chaque cible est illuminée par beaucoup d’impulsions
et produira un grand nombre d’échos
Introduction

Géométrie
des images
radar
Résolution
radiale
IRésolution
azimutale
Distorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Fréquence Doppler

9
9
souce [FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991]
Sébastien Giordano  30-11-2018  38 / 38 Radar PPMD
Effet Doppler

Introduction

Géométrie
des images
radar Effet Doppler
Résolution
radiale
IRésolution
f0 porteuse, fréquence centrale
azimutale
Distorsions
point terrain commence à être illuminé : en avance du radar, donc f > f0
géométriques

Radiométrie
fréquence Doppler : écart entre f et la porteuse
Radargrammétrie impulsion après impulsion la fréquence doppler diminue
et Radarcli-
nométrie nulle lorsque le point se trouve dans la direction perpendiculaire à la
Interférommétrie trajectoire
polarimetrie
devient négative quand le point se trouve en arrière du radar

Sébastien Giordano  30-11-2018  39 / 39 Radar PPMD


Synthèse d’ouverture

Présence d’un chirp


Traitement de compression
Introduction
d’impulsions
Géométrie
des images
radar
Analogie Compression en azimut
Résolution
radiale
/ Compression en distance
IRésolution
azimutale
Résolution en azimut dépend de
Distorsions
géométriques
la largeur de bande ∆f
Radiométrie

Radargrammétrie résolution azimutale, synthèse


et Radarcli-
nométrie d’ouverture
D
Interférommétrie
δaz = 2
polarimetrie

10
10
[Link] IIN/RESSOURCES/Cours/radar/[Link]
Sébastien Giordano  30-11-2018  40 / 40 Radar PPMD
Antenne synthétique

Introduction

Géométrie
des images
radar
Résolution
radiale
IRésolution
Antenne Synthétique
azimutale
Distorsions La synthèse d’ouverture permet de simuler une antenne de très grande
géométriques
taille
Radiométrie

Radargrammétrie dans le cas du SAR, l’antenne doit être la plus petite possible
et Radarcli-
nométrie mais antenne trop petite = mauvais rapport signal-sur-bruit
Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  41 / 41 Radar PPMD


Les effets du relief

Introduction

Géométrie
des images
radar
Résolution
radiale
Résolution
azimutale
Distorsions
IDistorsions
géométriques Distorsions géométriques et radiométriques présentes sur les images radar :
Radiométrie Distorsions dues à l’élévation du terrain (Elevation distortion)
Radargrammétrie
et Radarcli- Distorsions dues à la distance (Range distorsion)
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  42 / 42 Radar PPMD


Rétrodiffusion et angle d’incidence local

Une pente faisant face au faisceau aura une rétrodiffusion plus élevée qu’une
Introduction pente située à l’opposé et apparaı̂tra donc sur l’image avec une tonalité plus
Géométrie claire.
des images
radar
Résolution
radiale
Résolution
azimutale
IDistorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: Relation entre intensité et angle d’incidence (source


[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991])

Sébastien Giordano  30-11-2018  43 / 43 Radar PPMD


Les Ombres

Introduction

Géométrie
des images
radar
Résolution
radiale
Résolution
azimutale
IDistorsions
Figure: Relation entre ombre et distance oblique (source
géométriques [FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991]
Radiométrie

Radargrammétrie
et Radarcli- Des pentes ayant un angle plus abrupt que l’angle de dépression sont sujettes à
nométrie
des ombres portées. Les ombres sont plus importantes pour les points éloignées
Interférommétrie
(far range) que pour les points rapprochés (near range). Aucune information
polarimetrie
dans les ombres.

Figure: Ombre portée (source


[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991])

Sébastien Giordano  30-11-2018  44 / 44 Radar PPMD


Les Ombres

Introduction

Géométrie Exploitation altimétrique des ombres


des images
radar
Résolution
En géométrie slant-range : H = ∆Rslant cos(θ)
radiale
Résolution
En géométrie ground-range : H = ∆Rslant cos(θ) sin(θ)
azimutale
IDistorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: source : [?]

Sébastien Giordano  30-11-2018  45 / 45 Radar PPMD


Raccourcissement et Basculement

Introduction

Géométrie
des images
radar
Résolution
radiale
Résolution
azimutale
IDistorsions
géométriques

Radiométrie Figure: source :


Radargrammétrie [FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991]
et Radarcli-
nométrie

Interférommétrie

polarimetrie Raccourcissement et Basculement


Raccourcissement des distances - Foreshortening : A0 B 0 < AB 00 .
Raccourcissement des pentes inclinées vers le radar.
Basculement du relief - Layover. Dans des cas extrêmes, la position
des points peut être inversée. Distance oblique à B plus courte qu’à A.

Sébastien Giordano  30-11-2018  46 / 46 Radar PPMD


Réflexions multiples

Introduction

Géométrie
des images
radar
Résolution
radiale
Résolution
azimutale
IDistorsions
géométriques

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Figure: source : [?]

Sébastien Giordano  30-11-2018  47 / 47 Radar PPMD


Introduction

Géométrie
des images
radar

IRadiométrie
Facteurs
Speckle
Réflectivité
radar d’un
canal de
Radiométrie
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  48 / 48 Radar PPMD


Quantité mesurée par le radar

Introduction

Géométrie Equation radar


des images
radar
G .λ2 .q
IRadiométrie PR = PT σ0 (3)
Facteurs
(4πR)3
Speckle
Réflectivité
radar d’un avec :
canal de
polarisation
SXY PR = Puissance moyenne reçue
Radargrammétrie PT = Puissance transmise
et Radarcli-
nométrie
G = Gain d’antenne
Interférommétrie

polarimetrie
λ = Longueur d’onde émise
q = Constante tenant compte des gains et pertes diverses
R = Distance oblique entre l’antenne et la surface imagée

Sébastien Giordano  30-11-2018  48 / 48 Radar PPMD


Calibration radiométrique

Introduction

Géométrie Plusieurs modes de calibrations, à choisir selon l’application.


des images
radar
Beta Naught β0 : mesure du coefficient de rétrodiffusion radar en slant
IRadiométrie
Facteurs
range, sans dimension. Cette mesure du coefficient de rétrodiffusion
Speckle dépend de la position dans l’image (angle d’incidence, ...) et de la
Réflectivité
radar d’un géométrie du terrain (pente) en plus des propriétés de la cible.
canal de
polarisation
SXY
Sigma Naught σ0 sans dimension, normalisé en unité de surface sur le
Radargrammétrie terrain ( [m2/m2], exprimé généralement en decibel (dB)). besoin d’un
et Radarcli-
nométrie
MNT pour passer de Beta Nought au Sigma Nought. dépend de l’angle
Interférommétrie
d’incidcence local.
polarimetrie Gamma Naught γ0 : Sigma Nought normalisé par l’angle d’incidence.
plus de dépendance de la visée en distance.

Sébastien Giordano  30-11-2018  49 / 49 Radar PPMD


Coefficient de rétrodiffusion radar

Introduction
σ0 : section efficace de diffusion radar par unité de surface de la cible, coef-
Géométrie
des images ficient de rétrodiffusion radar . Il est sans unité.
radar

IRadiométrie Coefficient de rétrodiffusion radar


Facteurs
Speckle Généralement exprimé en décibels (de l’ordre de 105 )
Réflectivité
radar d’un
canal de
polarisation
σ0dB = 10 log(σ0 ) (4)
SXY

Radargrammétrie Le coefficient de rétrodiffusion est fonction de :


et Radarcli-
nométrie
de la longueur d’onde ,
Interférommétrie

polarimetrie
de la polarisation,
de l’angle d’incidence
propriétés géométriques et géophysiques de la cible

Sébastien Giordano  30-11-2018  50 / 50 Radar PPMD


Influence de longueur d’onde

Introduction La longueur d’onde λ entre en compte


Géométrie dans :
des images
radar
la profondeur de pénétration de
Radiométrie
l’onde dans le milieu observé
Facteurs
Longueur
d’onde
le concept de rugosité de la
Angle
d’incidence
surface observée
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
11
nométrie

Interférommétrie

polarimetrie

11
[Link] [Link]
Sébastien Giordano  30-11-2018  51 / 51 Radar PPMD
Influence de longueur d’onde

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

12
Sébastien Giordano  30-11-2018  52 / 52 Radar PPMD
12
Influence de l’angle d’incidence

Angle d’incidence
Introduction

Géométrie L’angle d’incidence est formé par la perpendiculaire de la surface observée


des images
radar avec la direction de l’onde électromagnétique incidente.
Radiométrie
Facteurs
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie
Figure: Angle d’incidence d’après [?]
Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  53 / 53 Radar PPMD


Rugosité de surface

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

13
13
[FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture, 1991]
Sébastien Giordano  30-11-2018  54 / 54 Radar PPMD
Rugosité de surface

Introduction

Géométrie
des images La rétrodiffusion dépend de l’aspect de la superficie considérée.
radar

Radiométrie
Radar : une surface sera dite rugueuse si sa structure ou superficie ont
Facteurs des dimensions correspondants à une fraction importante de la longueur
Longueur
d’onde
d’onde incidente.
Angle
d’incidence une surface donnée peut être lisse pour λ = ( 25)cm et rugueuse pour
Rugosité de
surface λ = ( 3)cm
Humidité
Polarisation Exemple : rugosité de surface pour cartographie du trait de côte
Speckle
Réflectivité eau : très légèrement rugueuse, légère rétrodiffusion vers le radar
radar d’un
canal de
polarisation
surfaces terrestre plus rugueuses
SXY

Radargrammétrie
les zones d’eau apparaissent sombres
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  55 / 55 Radar PPMD


Rugosité de surface

Introduction

Géométrie
des images
radar
Critère de Rayleigh
Radiométrie
λ
Facteurs lisse si h < 8 cos(θ)
Longueur
d’onde λ
Angle
rugueuse si h > 8 cos(θ)
d’incidence
Rugosité de
surface
avec :
Humidité
Polarisation
h : hauteur moyenne des irrégularités
Speckle
λ : longueur d’onde
Réflectivité
radar d’un
canal de θ : angle d’incidence
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  56 / 56 Radar PPMD


Humidité

Introduction

Géométrie
des images l’humidité influe sur les propriétés électriques du sol et de la végétation
radar
car la constante diélectrique de l’eau est très élevée,
Radiométrie
Facteurs constante diélectrique ou permittivité : la réponse d’un milieu donné à
Longueur
d’onde un champ électrique appliqué,
Angle
d’incidence une surface humide fournira un écho plus fort qu’une surface plus sèche
Rugosité de
surface
Humidité
la pénétration de l’onde électromagnétique dans le milieu est une
Polarisation fonction inverse de la teneur en eau,
Speckle
Réflectivité Cas de la végétation : la profondeur de pénétration dépend de l’humidité,
radar d’un
canal de de la densité et de la structure géométriques des plantes. (feuilles,
polarisation
SXY branches).
Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  57 / 57 Radar PPMD


Influence de la polarisation

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs Polarisation
Longueur
d’onde


Angle
Direction du vecteur champ electrique E
d’incidence
Rugosité de
surface
Humidité
Direction perpendiculaire à la direction de propagation de l’onde
Polarisation
Speckle
Plan d’incidence : plan perpendiculaire à la surface observée
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  58 / 58 Radar PPMD


Influence de la polarisation

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie 14 15

Interférommétrie

polarimetrie

14
[Link] [Link]
15
[Link] [Link]
Sébastien Giordano  30-11-2018  59 / 59 Radar PPMD
Onde plane monochromatique

Introduction Onde plane : solution aux equations de Maxwell dans le vide


Géométrie
des images A une grande distance de l’antenne on peut considérer que le front
radar
d’onde est plan
Radiométrie
Facteurs
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY


Radargrammétrie E résultat d’une composante horizontale et verticale
et Radarcli- −

nométrie E prévisible : onde complètement polarisée
Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  60 / 60 Radar PPMD


Ellipse de polarisation

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Ellipse de polarisation
Réflectivité
radar d’un A : amplitude
canal de
polarisation
SXY
α : phase absolue
φ : angle d’orientation, φ ∈ − π2 ; π2
 
Radargrammétrie
et Radarcli-
τ : ellipticité (ou excentricité), τ ∈ 0; π4
 
nométrie

Interférommétrie

polarimetrie 16

16
illustrations : [?]
Sébastien Giordano  30-11-2018  61 / 61 Radar PPMD
Émission système radar

Introduction
antenne conçue pour transmettre et capter des ondes dans une
Géométrie
polarisation définie
des images
radar émission d’ondes de différentes polarisations séparément en utilisant un
Radiométrie commutateur
Facteurs
Longueur
réception simultanée de différentes polarisations de l’onde
d’onde
Angle polarisations linéaires horizontale, H, et verticale, V les plus utilisées
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
17

polarimetrie

17
illustration : [?]
Sébastien Giordano  30-11-2018  62 / 62 Radar PPMD
Polarisation

Configurations de polarisation des systèmes radar imageurs:


Introduction HH: émission horizontale , réception horizontale,
Géométrie
des images
HV : émission horizontale , réception verticale,
radar
VV : émission verticale, réception verticale,
Radiométrie
Facteurs VH : émission verticale , réception horizontale.
Longueur
d’onde
Angle
d’incidence
Rugosité de
surface
Humidité
Polarisation
Speckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
La polarisation peut apporter des informations sur les différentes couches du
polarimetrie
milieu observé : pénétration selon la polarisation choisie, renseigne sur la forme
et l’orientation des éléments diffuseurs 18
18
illustration : [?]
Sébastien Giordano  30-11-2018  63 / 63 Radar PPMD
Speckle

Principal inconvénient dans l’interprétation visuelle ou automatique des images


Introduction SAR : speckle (ou châtoiement)
Géométrie
des images
radar
Speckle
Radiométrie dû à la cohérence de l’onde radar : interférence des champs rétrodiffusés par
Facteurs
ISpeckle
les objets à l’intérieur de chaque cellule de résolution
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

19

19
schéma : [?]
Sébastien Giordano  30-11-2018  64 / 64 Radar PPMD
Marche aléatoire dans le plan complexe

Introduction

Géométrie
des images
radar

Radiométrie speckle : lié au déphasage de chaque diffuseur élémentaire φn .


Facteurs
ISpeckle problème déterministe (si les N diffuseurs et l’antenne radar sont placés
Réflectivité
radar d’un
exactement à la même position, on obtient la même valeur de SXY ), mais
canal de
polarisation en pratique impossible à reproduire (micro-configuration des
SXY
diffuseurs).
Radargrammétrie
et Radarcli-
nométrie
somme cohérente des réponses des N cibles ponctuelles est vue comme
Interférommétrie
une marche aléatoire dans le plan complexe [Goodman, 1975].
polarimetrie

Sébastien Giordano  30-11-2018  65 / 65 Radar PPMD


Speckle peinement développé

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs
ISpeckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

20

20
cours Irena Hajnsek, DLR
Sébastien Giordano  30-11-2018  66 / 66 Radar PPMD
Marche aléatoire

=(SXY )
Introduction

Géométrie
|SnXY |
des images
radar φnXY
Radiométrie
Facteurs
ISpeckle
Réflectivité
radar d’un
canal de
polarisation
φXY
SXY

Radargrammétrie
et Radarcli- <(SXY )
nométrie

Interférommétrie [3]
polarimetrie |SXY |

Sébastien Giordano  30-11-2018  67 / 67 Radar PPMD


Speckle

Introduction
On modélise le speckle comme un bruit multiplicatif mais ce n’est pas un
Géométrie
des images
bruit.
radar
Description statistique du speckle : modèle de Goodman du speckle
Radiométrie
Facteurs
(1976)
ISpeckle
Réflectivité
L’énergie reçue au niveau de l’antenne est la somme de contributions
radar d’un
canal de
indépendantes et aléatoirement phasées de diffuseurs à l’intérieur de la
polarisation
SXY même cellule de résolution
Radargrammétrie
et Radarcli-
les ondes réfléchies doivent être en assez grand nombre, uniformément
nométrie distribuées et aléatoirement phasées
Interférommétrie

polarimetrie
Hypothèse
Toutes ces hypothèses assurent que les parties réelles et imaginaires du speckle
sont des variables aléatoires qui obéissent à une distribution gaussienne.

Sébastien Giordano  30-11-2018  68 / 68 Radar PPMD


Zone

Introduction

Géométrie
des images
radar

Radiométrie
Facteurs
ISpeckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  69 / 69 Radar PPMD


Speckle : introduction

Problème : forte variabilité de la mesure radar pour des zones homogènes en


Introduction
occupation du sol.
Géométrie
des images
radar

Radiométrie
Facteurs
ISpeckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie
Extrait type 2 Extrait type 4 Extrait type 8
Figure: Trois décompositions de Pauli 1-look sur des zones homogènes en
occupation du sol. Ces figures illustrent le phénomène du speckle.

Sébastien Giordano  30-11-2018  70 / 70 Radar PPMD


Statistiques du Speckle d’après [Goodman, 1976]

Introduction

Géométrie Somme cohérente de diffuseurs élémentaires


des images
radar N
X
Radiométrie SXY = |SnXY | e jφnXY (5)
Facteurs
n=1
ISpeckle
Réflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

(a) Diffuseurs (b) Image Radar (c) Amplitude d’une


élémentaires d’une une 1-look zone homogène
cellule de résolution

Sébastien Giordano  30-11-2018  71 / 71 Radar PPMD


Somme cohérente

Introduction
On considère qu’une cellule de résolution radar est constituée de N diffuseurs
Géométrie
élémentaires. Ces N diffuseurs élémentaires sont répartis aléatoirement dans
des images
radar
la cellule de rétrodiffusion. Chacun de ces diffuseurs élémentaires possède une
Radiométrie
amplitude élémentaire |SnXY | et une phase élémentaire φnXY pour un canal de
Facteurs polarisation XY . Chaque canal de polarisation est constitué par la somme
ISpeckle
Réflectivité
complexe :
radar d’un
canal de
N

polarisation
|SnHH | e jφnHH
P
SXY


 SHH =
n=1

Radargrammétrie 

et Radarcli-
 N
|SnHV | e jφnHV
 P
 SHV =
nométrie


n=1
Interférommétrie
N (6)
|SnVH | e jφnVH
P
SVH =
polarimetrie





 n=1
N


|SnVV | e jφnVV
 P
 SVV =

n=1

Sébastien Giordano  30-11-2018  72 / 72 Radar PPMD


Hypothèses de Goodman du speckle pleinement développé

Introduction

Géométrie
des images Pour déterminer les lois de distribution de l’amplitude et de l’intensité, [Goodman, 19
radar
fait les hypothèses suivantes (hypothèses du speckle pleinement développé):
Radiométrie
Facteurs 1 Le nombre de diffuseurs élémentaires N par cellule de résolution est
ISpeckle
Réflectivité suffisamment grand;
radar d’un
canal de
polarisation 2 L’amplitude et la phase (|SnXY |, φnXY ) d’un diffuseur n sont
SXY
indépendantes de celles des autres diffuseurs;
Radargrammétrie
et Radarcli-
nométrie
3 Pour chaque diffuseur n l’amplitude |SnXY | et la phase φnXY sont
Interférommétrie indépendantes entre elles;
polarimetrie 4 Les phases φnXY sont uniformément distribuées entre [−π, π].

Sébastien Giordano  30-11-2018  73 / 73 Radar PPMD


Statistiques
Pour un canal de polarisation SXY , la somme cohérente (6) s’écrit en partie
réelle et imaginaire avec l’équation (7).
Introduction
 N
Géométrie P
 iXY = |SnXY | cos(φnXY )
des images


radar n=1
N (7)
Radiométrie P
 qXY = |SnXY | sin(φnXY )


Facteurs
Speckle n=1
IRéflectivité
radar d’un
canal de
En prenant l’espérance mathématique E () on obtient :
polarisation
SXY
 N
P
 E (iXY ) = E ( |SnXY | cos(φnXY ))

Radargrammétrie 
et Radarcli- n=1
nométrie N (8)
P
 E (qXY ) = E ( |SnXY | sin(φnXY ))


Interférommétrie
n=1
polarimetrie

Avec l’hypothèse (2) de [Goodman, 1975] on peut écrire :


 N
P
 E (iXY ) = E (|SnXY | cos(φnXY ))


n=1
N (9)
P
 E (qXY ) = E (|SnXY | sin(φnXY ))


n=1

Sébastien Giordano  30-11-2018  74 / 74 Radar PPMD


Statistiques

Introduction
Avec l’hypothèse (3) de [Goodman, 1975] on a :
Géométrie  N
des images P
 E (iXY ) = E (|SnXY |).E (cos(φnXY ))

radar

n=1
Radiométrie N (10)
P
Facteurs
 E (qXY ) = E (|SnXY |).E (sin(φnXY ))


Speckle
n=1
IRéflectivité
radar d’un
canal de
polarisation L’hypothèse (4) de [Goodman, 1975] assure que :
SXY

Radargrammétrie E (cos(φnXY ) = 0
et Radarcli- (11)
nométrie E (sin(φnXY ) = 0
Interférommétrie

polarimetrie Et par conséquent :



E (iXY ) = 0
(12)
E (qXY ) = 0

Sébastien Giordano  30-11-2018  75 / 75 Radar PPMD


Statistiques du speckle

L’espérance mathématique de la partie réelle et de la partie imaginaire d’un


canal de polarisation est donc nulle. On s’intéresse maintenant à l’espérance
Introduction
de l’amplitude et de l’intensité d’un canal de polarisation E (AXY ) et E (IXY ).
Géométrie
des images
En considérant les hypothèses (2) et (3) de [Goodman, 1975] on a :
radar
N P N

Radiométrie 2 P
Facteurs


 E (iXY ) = E (|SkXY | . |SmXY |).E (cos(φkXY ). cos(φmXY ))
k=1 m=1


Speckle 
IRéflectivité
 N P N
2 P
radar d’un
canal de
E (qXY )= E (|SkXY | . |SmXY |).E (sin(φkXY ). sin(φmXY )) (13)
polarisation

 k=1 m=1
SXY

 PN P N
E (|SkXY | . |SmXY |).E (cos(φkXY ). sin(φmXY ))

Radargrammétrie  E (iXY .qXY ) =

et Radarcli- k=1 m=1
nométrie

Interférommétrie L’hypothèse (4) de [Goodman, 1975] sur les phases uniformément distribuées
polarimetrie donne :
  1
2
si k = m
E (cos(φkXY ). cos(φmXY )) = E (sin(φkXY ). sin(φmXY )) =

0 si k 6= m
E (cos(φkXY ). sin(φmXY )) = 0 ∀k, m

(14)

Sébastien Giordano  30-11-2018  76 / 76 Radar PPMD


Statistiques

Introduction

Géométrie
On obtient alors :
des images 
radar N
2
) = 21 E ( |SnXY |2 )
P
E (iXY


Radiométrie


 n=1
Facteurs
N (15)
E (qXY ) = 2 E ( |SnXY |2 )
2 1
Speckle
P
IRéflectivité


radar d’un

 n=1
E (iXY .qXY ) = 0
canal de

polarisation
SXY

Radargrammétrie On note souvent R la réflectivité radar (16). La réflectivité radar augmente


et Radarcli-
nométrie avec le nombre de diffuseurs élémentaires contenus dans la cellule.
Interférommétrie
N
X
2 2
polarimetrie
RXY = E (IXY ) = E (iXY ) + E (qXY ) = E( |SnXY |2 ) (16)
n=1

Sébastien Giordano  30-11-2018  77 / 77 Radar PPMD


Lois de distribution

Introduction

Géométrie
des images
Pour calculer les lois de distribution de l’amplitude AXY et de l’intensité IXY ,
radar [Goodman, 1975] fait l’hypothèse que le nombre de diffuseurs élémentaires N
Radiométrie est très grand (hypothèse (1)). La partie réelle iXY et la partie imaginaire qXY
Facteurs
Speckle
correspondent à une somme de variables aléatoires indépendantes. Avec N très
IRéflectivité
radar d’un
grand et en faisant aussi l’hypothèse d’une loi à variance finie, le théorème
canal de
polarisation
central limite montre sous ces conditions que iXY et qXY convergent vers des
SXY
lois normales de moyenne nulle et de même variance (17).
Radargrammétrie
et Radarcli-

i2
nométrie  p(i |R) = √1 e − XY

R
XY πR
Interférommétrie
2
qXY
(17)
polarimetrie  p(qXY |R) = √1 e − R

πR

Sébastien Giordano  30-11-2018  78 / 78 Radar PPMD


p
2 2
Introduction
Sachant que AXY = iXY + qXY (44) et que p(iXY |R) et p(qXY |R) suivent des
Géométrie
des images
lois normales, on sait par définition que p(AXY |R) suit une loi de Rayleigh ce
radar qui est traduit par l’équation (17). IXY est distribuée selon une loi exponentielle
Radiométrie décroissante (19) d’espérance E (IXY ) = R et de variance var (IXY ) = R 2 .
Facteurs
Speckle
2A − AR2
IRéflectivité
radar d’un p(AXY |R) = e A≥0 (18)
canal de R
polarisation
SXY

Radargrammétrie 1 − RI
et Radarcli- p(IXY |R) = e I ≥0 (19)
nométrie R
Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  79 / 79 Radar PPMD


lois

Ces lois de distribution sont données en figures 20, 21 et 22 pour les types
Introduction d’occupation du sol définis enCHAPITRE
table 1. [Link]
CHAPITRE RAPPELS
RAPPELSSUR LA
LAPOLARIMÉTRIE
SURd’occupation
type POLARIMÉTRIERADAR
RADAR
du sol est
Géométrie
des images
caractérisé par sa propre réflectivité radar RXY .
radar

Radiométrie ��������� �����


��������� ������� ��������� ����� �� ��
��������� �����

Facteurs � ������
������� ����� �� ��
�����
����� ����� ����� ������ �������
� ����� ����� ����� �������� ��� � �
����� ����� �����
����� ������� ����� ����� �����
���� ������ ���
������ �������� ����� ������ �
������� ������ �������� ����� �������� ���
� �
������ �������� �����
����� �������
�����
�� �� ������ �������� �����
���� ������ ���
Speckle � ������
������� � ����� � �
������� ����� ������ ����
����� ������� ������� ����� ������ ���
������� ����� ������ ������� ������� ����� ������ ������� ���
������ ���
���� �������� ����� ������ �
������� ���� �������� ����� �������� ���
� �
���� �������� �����
����� �������
�� �� ���� �������� �����
���� ������ ���
IRéflectivité ������ ������� ���������
������� �����
��������� ����� ������
����� �������

������� �����
����� ������� ���������
������� �����
��������� ����� �������� ���
������ �
���

� ������� ����
����� ����� ������ �
������� ����� ����� �������� ���
� �
����� �����
����� ������� ����� �����
���� ������ ���
radar d’un ������� ������
������� �����
������ ����� ������
����� �������

������� ����� �� ��
����� �������
������������� �����
������ ����� ��������
���� ���
������ � �
���
� ������
�������
canal de
������������
������������

������������
������������
���� �� ��
� �����
polarisation � ������
�������

SXY ����� �� ��
�����

� ������
�������

���� �� ��
� �����
Radargrammétrie � ������
�������
���� �� ��
�����
et Radarcli- �

������
�������
nométrie �
����� �� ��
�����

�� �� � ��

�� �� ����
����� ����
����� ����
����� ����
����� �����
������ �� �� ������
������� � � � � � � � � � � � � � �
Interférommétrie � � � � � � �����
����� ������
� ������� �����
����� �������
������� �����
����� �������
������� �����
�����

������ ������

polarimetrie
(A
PP(A
P(A HH
HH
HH |R)
|R)
|R) P(I(IHH
PP(I HH |R)
|R)
HH|R)
3.3.9– –Lois
Figure3.3.9
Figure Loisdededistribution
distributionthéoriques
théoriquesdedeAAHH etetIIHH destypes
HHdes typesd’occupation
d’occupationdu dusol
sol
Figure: Lois
dusite
de distribution
sited’Oraison
d’Oraisonsachant
théoriques
sachantleurs
de AHH
leursréflectivités
etHH
réflectivitésradar
IHH des
radar(paramètres
(paramètres
types d’occupation
depauli
paulide
delalatable
du
table3.1)
sol
3.1)
du de
du site d’Oraison sachant leurs réflectivités radar (paramètres de Pauli de la table 1)

���������
��������� ���
����� � ��� ���������
��������� � ��
������
�����

�������
������ �� ��
�������
� ������

����������
���������� ��������
����� ������� � �
������ ����� ����� ������� ������ � �
���
�� ������ ����� ����� ����� ���������� ���
������
������ ��������
�������� ��������
����� ������� � �
������ ������ �������� ������� ������ � �
���
� �������
������ �� ������ �� ��
����
�����

������ �������� ����� ���������� ���
�������
������� ���������� ��������
����� ������� � �
������ ������� ����� ������� ������ � �
���
�� ������ ������� ����� ����� ���������� ���

�������� ��������
�������� ��������
����� ������� � �
������ ���� �������� ������� ������ � �
���
�������
�� ������ �� ��
���� �������� ����� ���������� ���
� ������ ������� �������� �
��������������������������
� � � �� ������������
� � � ����� � � �������
������ ������� �������� �
�������������������������
� � � �� � � � ����
���������� � �
�� �

���������� ��������
����� ������� � �
������ ����� ������� ������ � �
���
Sébastien Giordano������ �� ������  80 /��������
80 Radar PPMD ����� ����� ���������� ���
������ � �
 30-11-2018
�������
������� ������
������ ��������
����� �������
�� ������
� �
������ �� �� � � �������
������� ������
������ �������
����� ������
���������� ���
� �
���
�����
�� � �
� � �

�� ����� ����� ����� ����� � ����� �
� ������� ����� � ������� � ����� � ������� � ����� � ������� � ������
� � ����� ����� ����� ����� ������ �� �������

�����

�������

�����

�������

�����

�������

�����
��� ���
��� ���

Lois PP(A PP(I


(AHH
HH |R)
|R) (IHH
HH|R)
|R)
Figure3.3.9
Figure 3.3.9––Lois
Loisdededistribution
distributionthéoriques
théoriquesde deAAHH
HH et
et II
HH
HH des
des types
types d’occupation
d’occupation du du sol
sol
du
dusite
sited’Oraison
d’Oraisonsachant
sachantleurs
leursréflectivités
réflectivitésradar
radar(paramètres
(paramètres dede pauli
pauli de
de la
la table
table 3.1)
3.1)

Introduction

Géométrie
des images ��������� ������
��������� ������ ��������� �
������ �
��������� ������

� ������
������� ������ �� ��
�������
radar ����� �����
����� �����
����� �������
�������� ������

�������

�����
����� ����� ������
����� ����
��

���������� �
������������
������ �������� ����� ������� �
������� ������ �������� ������

���������� �
� ������
�������
������ �������� �������� ������ ������ �� ������ �������� ����
��������������
� � ����� �
�������
������������ �����
����� �������
�������� �������
������ �������
������� ����� ������
����� ����
�� ���������� �
������������
� �
Radiométrie � ������
�������
����
����
�������
��������

���������
�����
��������
�����
�������
��������
�������
�������
������

������� �������
�� ��
������
����
����
�������
��������
��������
���������
������
����
������
�� ���������� �
������������

���������� �
������� ��������� �������� ������ � ������� ��������� ����
��������������
����� � �� � ���� �� � � ����� ������

���������� �
����������������������
� � � � � ����� ����
��������������
Facteurs ������
������� �������
������������� �����
������ ����� �������
������

������� ������ �� �������
������� ������ ������
������ ����

���������� �
������������

����� � ����� ��
������������
������������

������������
Speckle

������������
� �����
������ � ������ �� ��
�������

IRéflectivité � ������
������� � ������ ��
�����

radar d’un
� ������
������� � ������ �� ��
�������
canal de
polarisation � ������
������� � ������ ��
�����

SXY
� ������
������� � ������ ��
�����

Radargrammétrie � � ��
�� �� ������� ���� ���� ���� ���� ���� ���� ����
� � ��
�� �� ������ ������ ������ ������ ��
� � ����� � ����� � ����� � ����� � ����� � ����� � ����� ������
� � ������
� � ������
� � ������
� � ����
����
� � ������� �
������
� � ��������� �
������
et Radarcli- ���
��� ���
���
nométrie

Interférommétrie
(AHV
PP(A
P(A HV |R)
|R)
HV|R) P (IHV
P(I
P(I HV |R)
|R)
HV|R)
3.3.10––Lois
Figure3.3.10
Figure Loisde
dedistribution
distributionthéoriques
théoriquesde deAAHV etIIHV
et des types
HV des types d’occupation
d’occupation du du sol
sol
polarimetrie Figure: Lois
dusite
de distribution
sited’Oraison
d’Oraisonsachant
théoriques
sachantleurs
de AHV
leursréflectivités
réflectivitésradar
etHV
IHV des
radar(paramètres
(paramètres
types d’occupation
de pauli
pauli de
de la
la table
du sol
table 3.1)
3.1)
du de
du site d’Oraison sachant leurs réflectivités radar (paramètres de Pauli de la table 1)

���������
��������� �������
����� ���������
��������� �����
���
���
���
������
����� �� ��
�������

� � ������ �
�����
����� ���������������������
����� ����������������� �� �����
����� ����������������������
����� �����������������

������
�������
� ������
������ �������������
������������������ �����������
������
���
���
�� ������
������ �������������������������
��������������� ����������
��
�������
������� ����������
���������������
������
������� ���
���
��

�� ��
�������
������ �������
������� ����������������������
������������ ����������
��


�������
������ �������� �������������
������������������
������
������� ���
���
�� ����
���� �������������������������
��������������� ����������
��
�������
������� ��������������
�������������������
������
������� ���
���
�� �������
������� ��������������������������
���������������� ����������
��
������� �����
����� �����
����������
������
������� ���
���
��
������ ��
����� �����
����� �����������������
������� ����������
��
������ �
��� ��

�������
������� ������������������������
����������� ������������
�������
������� �����������
����������� ��
�����������
�����������
�� ��
�������
������
������������

������������

�� ��

������������

�����
������������

� ����
�������
� ������
���� ��
�� �
������
Sébastien Giordano������ �

 30-11-2018  81 /����
81 �
Radar PPMD
� ������ �����
�� �

� � �

�� ���� ����� ����� ����
� ����� ����� ����� � ���� �
� ������
� ������� ������� ������� ����� � ������� � ������� �
� � ���� ����� ����� ����� ����� ����� ����� ����� �� ������� ������� ������� �������

�����

�������

�������
��� ���
��� ���

VV PP(A PP(I
(AHV
HV|R)
|R) (IHV
HV|R)
|R)
Figure3.3.10
Figure 3.3.10––Lois
Loisde
dedistribution
distributionthéoriques
théoriquesde deAAHV
HV et
et II
HV
HV des
des types
types d’occupation
d’occupation du du sol
sol
du
dusite
sited’Oraison
d’Oraisonsachant
sachantleurs
leursréflectivités
réflectivitésradar
radar(paramètres
(paramètres de
de pauli
pauli de
de la
la table
table 3.1)
3.1)

Introduction

Géométrie
des images ��������� ������
��������� ������ ��������� �
������ �
��������� ������

� �����
������ ������ �� ��
�������
radar ����� �����
����� �����
����� ����� �������
����� ������

�������

�����
����� ����� ������
����� ����

���������� �
������������
��
�������������� ����� ����������� �
������� ������ �������� ������

���������� �
� ������
������� ������ �������� �����
����� ������ ������ �������� ���� ������������
��
������� ����� ����� ������� �
������� �� �� ������� ����� ������

���������� �
������� ����� �����
����� ������ � ������
������� ������� ����� ���� ������������
��
� �
Radiométrie � ������
�������
����
�������
��������
����
���������
�����
�������� �����
�����
�������
�����
�������
�������
������

�������
����
����
�������
��������
��������
���������
������
����
������
���������� �
������������
��

���������� �
������� ��������� �����
����� ������ ������� ��������� ���� ������������
��
����� ����� ������� �
������� �� �� ����� ������

���������� �
� ������
������� ����� �����
����� ������ �����
����� ����� ���� ������������
��
Facteurs ������� ������
������� �����
������ ����� �����������
������

������� �������
������� ������ ������
������ ����

����������� �
����� �������������
��
� ������
�������
������������
������������

������������
Speckle �� ��

������������
�����
�����

� ������
�������
IRéflectivité �� ��
������ �����
�����
�������
radar d’un �

canal de � ������
�������
������� ��
�����

polarisation � ������
�������

SXY ������� ��
�����
� ������
�������

Radargrammétrie �� ��
�� �� ���� ���� ���� ���� ���� ����
� ��

�� �� ������ � � �� �� �� ��
� ����� � ����� � ����� � ����� � ����� � ����� � ������� ����
����
� � ������
������
� � ����
����
� � ������
������
� � ����
����
� � ��������� �
������ ��������
����
et Radarcli- ������ ���
���
nométrie

Interférommétrie
P(A(AVVV
PP(A |R)
|R)
VVV|R) (IVVVV
P(I
P(I
P |R)
|R)
VV|R)

3.3.11––Lois
Figure3.3.11
Figure Loisde
dedistribution
distributionthéoriques
théoriquesdedeAA et etIIVVVV des
des types
types d’occupation
d’occupation du du sol
sol
polarimetrie Figure: Lois de distribution théoriques de AVV etVVVIVVV des types d’occupation du sol
dusite
du sited’Oraison
d’Oraisonsachant
sachantleurs
leursréflectivités
réflectivitésradar
radar(paramètres
(paramètres de de pauli
pauli de
de la
la table
table 3.1)
3.1)
du site d’Oraison sachant leurs réflectivités radar (paramètres de Pauli de la table 1)

versiondu
version du30
30septembre
septembre2015
2015 55
55

Sébastien Giordano  30-11-2018  82 / 82 Radar PPMD


Filtrage
En radar une observation seule ne suffit pas. Il faut donc, à partir de plusieurs
observations, pouvoir estimer E (IXY ) ou E (AXY ) ou dans le cas polarimétrique
Introduction
E ([C ]) ou E ([T ]). Deux stratégies sont envisageables :
Géométrie le filtrage temporel : utilisation d’une série temporelle
des images
radar filtrage spatial : une seule image est utilisée et plusieurs pixels (par
Radiométrie exemple un voisinage local) contribuent à l’estimation.
Facteurs
Speckle Approche la plus simple BoxCar Filter : moyenne locale dans un voisinage
IRéflectivité
radar d’un
N × M pixels.
canal de
polarisation
N X M
SXY
1 X
Radargrammétrie < IXY >N×M = IXY (i, j) (20)
et Radarcli- N × M i=1 j=1
nométrie

Interférommétrie

polarimetrie N X M
1 X
< C >N×M = C (i, j) (21)
N × M i=1 j=1

Plus on augmente le nombre de looks (N×M) et plus l’estimation de l’espérance


est correcte. Cependant ce moyennage se fait au détriment de la résolution
spatiale et la probabilité d’avoir des types d’occupations du sol différents dans
l’échantillon des observations utilisées pour le calcul de l’espérance augmente.
Sébastien Giordano  30-11-2018  83 / 83 Radar PPMD
Filtrage

Introduction

Géométrie Idée : sélectionner de façon adaptative les pixels radar qui contribuent au
des images
radar multilooking plutôt que d’utiliser un simple voisinage local
Radiométrie Quatre catégories proposées par [Deledalle et al., 2015] : Oriented Windows,
Facteurs
Region Growing, Point-Wise et Patch-Wise.
Speckle
IRéflectivité
radar d’un Voisinage orienté : Lee Refined [Lee, 1981] sélection parmi huit fenêtres
canal de
polarisation orientées.
SXY

Radargrammétrie Croissance de régions : [Vasile et al., 2006] proposent une méthode


et Radarcli-
nométrie
d’estimation à partir d’un voisinage donné par une croissance de région.
Interférommétrie Sélection des pixels dans un voisinage : sélectionner seulement
polarimetrie certains pixels dans un voisinage local. C’est l’idée du très utilisé Sigma
filter proposé par [?].

Sébastien Giordano  30-11-2018  84 / 84 Radar PPMD


Lois de distribution pour les données multilookées
Dans le cas de données radar N-look, les lois de distribution de l’intensité (22)
et de l’amplitude (23) sont données en équations (24) et (25) avec Γ() la loi
Introduction
gamma d’Euler.
Géométrie
des images N
radar 1 X
INXY = In (22)
Radiométrie
N n=1 XY
Facteurs
Speckle
IRéflectivité
radar d’un
canal de v
polarisation
SXY
u
u1 X N
Radargrammétrie ANXY =t In (23)
et Radarcli- N n=1 XY
nométrie

Interférommétrie

polarimetrie N N INN−1 NIN


XY
p(INXY |R) = XY
N
e− R (24)
Γ(N)R

NA2
2NA2N−1
NXY
NXY
p(ANXY |R) = e− R (25)
Γ(N)R N

Sébastien Giordano  30-11-2018  85 / 85 Radar PPMD


CHAPITRE 3. RAPPELS SUR LA POLARIMÉTRIE RADAR

Lois de distribution pour les données multilookées


La figure 3.4.13 donne la loi de distribution de l’intensité N-look INXY pour le type d’occupation
Lasol
du figure 23 donne
4 (culture d’orge) la
et loi
pourde distribution
plusieurs valeurs de l’intensité
de N . N-look INXY pour le type
d’occupation du sol 4 (culture d’orge) et pour plusieurs valeurs de N.
Introduction ����������������������������������������������

Géométrie � ������
���
des images
���
radar
���
Radiométrie
� ������ ���
Facteurs
���
����
Speckle
IRéflectivité
� ������
������������

radar d’un
canal de
polarisation
SXY � ������
Radargrammétrie
et Radarcli-
nométrie � ������
Interférommétrie

polarimetrie � ������

��
�� �����
� �����
� �����
� ������ ������ ������ ������ ������ ������
����

Loi 3.4.13
Figure
Figure: – Loi de distribution
de distribution théoriquethéorique de l’intensité
de l’intensité N-lookN-look pour l’image
INHH l’image
INHH pour
RS2-FORCA-003 pour le type d’occupation du sol culture orge
RS2-FORCA-003 pour le type d’occupation du sol culture orgedéfini en table 3.1
défini en table 1
Sébastien Giordano  30-11-2018  86 / 86 Radar PPMD
Images radar polarimétrique
résolution spatiale

Slant Range
Introduction
[4, 73 m × 4, 74 m]
Géométrie
des images Ground Range
radar
[7, 48 m × 4, 74 m]
Radiométrie
Facteurs
Speckle
IRéflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  87 / 87 Radar PPMD


Images radar polarimétrique
résolution spatiale, speckle

Ground Range
Introduction
[7, 48 m × 4, 74 m]
Géométrie
des images 1-look
radar

Radiométrie
Facteurs
Speckle
IRéflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  87 / 87 Radar PPMD


Images radar polarimétrique
résolution spatiale, speckle

Introduction Ground Range


Géométrie [7, 48 m × 4, 74 m]
des images
radar 1-look
Radiométrie
Facteurs
Speckle
IRéflectivité
radar d’un
canal de
polarisation
SXY

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
81-looks
polarimetrie

Lee Sigma Filter

Sébastien Giordano  30-11-2018  87 / 87 Radar PPMD


Images radar polarimétrique
résolution spatiale, speckle, mécanismes de rétrodiffusion

Introduction Ground Range


Géométrie [7, 48 m × 4, 74 m]
des images
radar Nomenclature
Radiométrie
Facteurs 1 plan d’eau
Speckle 2 jeune mais
IRéflectivité
radar d’un
3 galets
canal de
polarisation
4 orge
SXY 5 urbain
Radargrammétrie 6 sainfoin
et Radarcli- 7 foret
nométrie
8 colza
Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  87 / 87 Radar PPMD


Exercice Matlab

Théorème central limite


Introduction Créer un vecteur de N échantillons pouvant décrire la variation
Géométrie temporelle d’un même pixel dans une série d’images radar, ou N pixels
des images
radar radar de la même image radar : observant dans les deux cas le même
Radiométrie milieu homogène.
Facteurs
Speckle
Pour chaque échantillon n de N, considérer que le résultat est la somme
IRéflectivité
radar d’un
de nbdiff (exemple nbdiff = 1000) diffuseurs répartis dans la cellule de
canal de
polarisation
diffusion (pixel).
SXY
Créer une variable aléatoire φ, qui pour chaque n est le résultat d’une
Radargrammétrie
et Radarcli- tirage aléatoire entre 0 et 2π.
nométrie

Interférommétrie
Observer la distribution des valeurs de φ, par exemple pour la dernière
polarimetrie
itération ( hist(Phi,100))
Faire varier le nombre N d’échantillons (exemple : progressivement de
100 à 100 000).
Calculer pour chaque itération la moyenne des φ(n) et observer la
distribution des valeurs de la moyenne hist(MoyPhi,100). Vers quoi tend
cette distribution en augmentant le nombre d’échantillons N ?

Sébastien Giordano  30-11-2018  88 / 88 Radar PPMD


Exercice Matlab

Introduction

Géométrie
des images
radar
Statistiques de l’Amplitude et l’Intensité Radar
Radiométrie
Facteurs Pour une observation radar SXY = iXY + [Link] = |SXY | e j.φXY
Speckle
IRéflectivité on prendra pour chaque diffuseur élémentaire une amplitude donnée
radar d’un
canal de
polarisation
(exemple : an = 1000)
SXY
Ajouter à votre code le calcul de la distribution de :
Radargrammétrie
et Radarcli-
1 la moyenne de la partie réelle iXY ,
nométrie 2 la moyenne de la partie imaginaire qXY ,
Interférommétrie
3 la moyenne de l’amplitude,
polarimetrie
4 la moyenne de l’intensité.

Sébastien Giordano  30-11-2018  89 / 89 Radar PPMD


Exercice Matlab

Introduction

Géométrie
des images Speckle Pattern
radar

Radiométrie
Créer un image de 200 × 200pixels...
Facteurs
Speckle
Pour les 4 blocs de 100 × 100 pixels prendre 4 valeurs de an , exemple
IRéflectivité an : 100, 300, 500, 100 et affecter au pixel la moyenne de l’amplitude avec
radar d’un
canal de
polarisation
la formule de l’étape précédente.
SXY
Normaliser pour que valeurs varient entre 0 et 1.
Radargrammétrie
et Radarcli-
nométrie
Enregistrer et observer l’image.
Interférommétrie Observer la forme de la loi de distribution et calculer l’écart-type pour
polarimetrie chacun des 4 blocs. Conclure sur la modélisation de ce phénomène
comme un bruit multiplicatif.

Sébastien Giordano  30-11-2018  90 / 90 Radar PPMD


Introduction

Géométrie
des images
radar

Radiométrie

IRadargrammétrie
et Radarcli-
nométrie
Radargrammétrie
Radargrammétrie et Radarclinométrie
Radarclinométrie

Interférommétrie

polarimetrie

Sébastien Giordano  30-11-2018  91 / 91 Radar PPMD


Radargrammétrie

Introduction

Géométrie
des images
radar Fonction de localisation : associer à un pixel (i, j) et étant donnée une
Radiométrie altitude z, les coordonnées géographiques (λ, φ) du point visé.
Radargrammétrie
et Radarcli- Intersection d’un modèle de terre (plan, ellipsoı̈de, MNT)
nométrie
IRadargrammétrie
Radarclinométrie Orientation interne
Interférommétrie
transformer (i, j) en grandeurs physiques (date, distance)
polarimetrie
R = R0 + (i − 1)∆R
t = t0 + (j − 1)∆t

Sébastien Giordano  30-11-2018  91 / 91 Radar PPMD


Radargrammetrie

Introduction
position de l’antenne
Géométrie
des images fréquence doppler centrale fd (généralement nulle)
radar

Radiométrie
longueur d’onde λ
Radargrammétrie Position P du point terrain
et Radarcli-
nométrie
IRadargrammétrie
Radarclinométrie
Equations de localisation
Interférommétrie
(p − s)2 = R 2
polarimetrie λfd R
( ds
dt
)(p − s) = 2

intersection d’une sphère et d’un cône


fréquence Doppler le plus souvent nulle (le cône est un plan)
( ds
dt
)(p − s) = 0

Sébastien Giordano  30-11-2018  92 / 92 Radar PPMD


Radargrammetrie

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie
IRadargrammétrie
Radarclinométrie

Interférommétrie

polarimetrie

ht!
Figure: Intersection sphère et cône Doppler

Sébastien Giordano  30-11-2018  93 / 93 Radar PPMD


Radargrammetrie

Géométrie Radar
Géométrie radar localement équivalente à une géométrie optique dont la
Introduction
direction de visée est la tangente au cercle distance Doppler.
Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie
IRadargrammétrie
Radarclinométrie

Interférommétrie

polarimetrie

ht!
Figure: Analogie avec PVA optique

forte distorsions = incidence plongeante (petit angle de visée radar) =


grand angle de visée PVA optique
cas extrême = radar proche du nadir, visée optique à l’horizon
Sébastien Giordano  30-11-2018  94 / 94 Radar PPMD
Vision stéréoscopique

Introduction

Géométrie
des images
radar

Radiométrie
Observation stéréoscopique SAR difficile
Radargrammétrie
et Radarcli-
nométrie
Speckle décorrélé entre images, filtrage par l’oeil humain
IRadargrammétrie
Radarclinométrie
fort B/H implique fortes différences radiométriques
Interférommétrie exemple : visées opposées, radiométrie corrélées négativement
polarimetrie
difficulté majeure : décorrélation des radiométries

Sébastien Giordano  30-11-2018  95 / 95 Radar PPMD


Corrélation automatique

Introduction

Géométrie
des images
radar Problèmes pour la corrélation automatique :
Radiométrie
optique : similitudes radiométriques (même éclairement), différent radar
Radargrammétrie
et Radarcli-
capteur actif
nométrie
IRadargrammétrie
speckle différent dans les images : hautes fréquences décorrélées
Radarclinométrie
difficile de restituer les microreliefs
Interférommétrie

polarimetrie
utiliser des objets de haut niveau
utilisation de plus de 2 images
zones d’ombres

Sébastien Giordano  30-11-2018  96 / 96 Radar PPMD


Radarclinométrie

Introduction

Géométrie
Radarclinométrie
des images
radar utilisation de la radiométrie pour la reconstruction du relief
Radiométrie
orientation locale du terrain à partir de la valeur radiométrique du pixel
Radargrammétrie
et Radarcli- changement de pente changement d’intensité : modèle de retrodiffusion
nométrie
Radargrammétrie intensité de l’image en fonction de l’angle d’incidence :angle entre le
IRadarclinométrie
direction antenne terrain et la direction perpendiculaire au terrain
Interférommétrie

polarimetrie
radarclinométrie : inversion du modèle

Difficultés :
intensité aussi liée à la réflectivité du terrain
à un angle d’incidence donné correspond non pas une unique orientation
de la surface mais une infinité d’orientations possibles

Sébastien Giordano  30-11-2018  97 / 97 Radar PPMD


Radarclinométrie

Exemple inversion de modèle pour un terrain de réflectance lambertienne,


Introduction l’équation de la radarclinométrie peur s’écrire de la façon suivante :
Géométrie
des images
radar
Equation de la radarclinométrie
2
Radiométrie
Ir = K σ0 δd δa sin(θ) cos(θ−α)
sin(θ−α)
cos(β)
Radargrammétrie
et Radarcli-
nométrie avec :
Radargrammétrie
IRadarclinométrie K constante qui exprime tous les paramètres liés au radar
Interférommétrie σ0 paramètres de réflectance liés au terrain
polarimetrie
Ir intensité reçue par la portion de terrain dans la cellule de résolution
radar
δd et δa résolution en distance et en azimut
α et β : deux angles définissant l’orientation de la parcelle par rapport
aux axes Ox et Oy
β angle de visée

Sébastien Giordano  30-11-2018  98 / 98 Radar PPMD


Radarclinométrie

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie
Radargrammétrie
IRadarclinométrie

Interférommétrie

polarimetrie

Figure: Thèse Pascal Daniel, Peut-on extraire le relief d’une seule image ?

Sébastien Giordano  30-11-2018  99 / 99 Radar PPMD


Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

IInterférommétrie
Interférommétrie
Mesurer
l’élévation
Mesurer le
mouvement

polarimetrie

Sébastien Giordano  30-11-2018  100 / 100 Radar PPMD


Interférométrie : introduction

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

IInterférommétrie
Mesurer
l’élévation
Mesurer le
mouvement

polarimetrie

Sébastien Giordano  30-11-2018  100 / 100 Radar PPMD


Introduction

Introduction

Géométrie
des images
radar
Interférométrie
Radiométrie L’interférométrie radar produit deux types d’information
Radargrammétrie
et Radarcli-
information topographique : mesurer l’élévation
nométrie
interférométrie de déplacement : mesurer le mouvement
IInterférommétrie
Mesurer
l’élévation
Mesurer le Principe
mouvement

polarimetrie Interférométrie radar depuis un couple d’images acquises d’un avion,


d’une navette spatiale, d’un satellite.
On calcule une différence de phase entre les deux images.

Sébastien Giordano  30-11-2018  101 / 101 Radar PPMD


Principe de l’interférométrie radar

Principe
Introduction
Point P sur le terrain vue de 2 images radar de deux points de vues
Géométrie
des images différents (M et S)
radar

Radiométrie
Phase d’un point : phase géométrique + phase propre
Radargrammétrie
et Radarcli-
nométrie

IInterférommétrie
Phases φM et φS
Mesurer
l’élévation
Mesurer le
mouvement φM = φgeom−M + φscattM
polarimetrie 4πMP
= + φM−scatt (26)
λ

φS = φgeom−S + φscatt−S
4πSP
PixelM = |SM | e jφM et = + φscatt−S (27)
λ
PixelS = |SS | e jφS

Sébastien Giordano  30-11-2018  102 / 102 Radar PPMD


Phase interférométrique

La technique de l’intérférométrie exploite la différence de phase ∆φint = φS −


Introduction φM : phase interférométrique ou interférogramme. ∆φint liée (entre autres) à
Géométrie l’altitude du point → possibilité de générer des modéles numériques d’élévation.
des images
radar
Phase interférométrique
Radiométrie

Radargrammétrie 4π(SP − MP)


et Radarcli- ∆φint = φS − φM = + (φscatt−S − φscatt−M ) (28)
nométrie λ
Interférommétrie
IMesurer Si l’état de surface n’a pas trop changé, et ligne de base pas trop importante :
l’élévation
Mesurer le pas de décorrélation spatiale ou temporelle
mouvement

polarimetrie Propriété

φscatt−S − φscatt−M = 0 (29)

Cette propriété mieux valable sur des surfaces stables (milieu urbain, roche,...)
moins bien en présence de végétation, champ labouré...
21

21
[Crosetto et al., 2016]
Sébastien Giordano  30-11-2018  103 / 103 Radar PPMD
Conditions de faisabilité

Introduction

Géométrie
des images
Pour que ∆φM soit petit il faut que :
radar
les deux phases propres de la cible correspondant aux deux échos
Radiométrie
présentent un fort coefficient de corrélation entre elles
Radargrammétrie
et Radarcli-
nométrie
les caractéristiques physiques de la cible n’aient pas changées entre les
Interférommétrie
deux prises de vues.
IMesurer
l’élévation
Mesurer le
mouvement Décorrélation
polarimetrie
On distingue dans la décorrélation deux origines principales :
la décorrélation spatiale
la décorrélation temporelle

Sébastien Giordano  30-11-2018  104 / 104 Radar PPMD


Topographie

Introduction

Géométrie
des images
radar
En première approche on a :
Radiométrie
∆φint est appelé interférogramme
Radargrammétrie
et Radarcli-
nométrie 4π(SP − MP)
∆φint = φS − φM = (30)
Interférommétrie λ
IMesurer
l’élévation
Mesurer le
mouvement
Mesure précise mais connue modulo 2π
polarimetrie l’information topographique créée des franges interférométriques
techniques de déroulement de phases

Sébastien Giordano  30-11-2018  105 / 105 Radar PPMD


Information topographique

information topographique

Introduction 4πd 4πhδi


Géométrie
∆φint = = (31)
des images
λ λ sin(i)
radar
λsin(i) Borth
Radiométrie h= ∆φ avec : δi = (32)
4πδi Dist
Radargrammétrie
et Radarcli- (33)
nométrie

Interférommétrie
IMesurer
l’élévation
Mesurer le
mouvement

polarimetrie

Sébastien Giordano  30-11-2018  106 / 106 Radar PPMD


Topographie

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
IMesurer
l’élévation
Mesurer le
mouvement

polarimetrie
Mont Etna, radar embarqué sur une Navette
Comparaison des phases de deux images prises à un jour d’intervalle depuis
des positions orbitales très proches. Une frange correspond à une différence
d’altitude de 500 m.
23

23
Massonnet et Inglada
Sébastien Giordano  30-11-2018  107 / 107 Radar PPMD
Altitude d’ambiguı̈té

Introduction

Géométrie
des images
radar A quelle variation d’altitude correspond une frange ?
Radiométrie λ sin(i) λDist sin(i)
∆φ = 2π implique h = 2δi
= 2Borth
Radargrammétrie
et Radarcli-
nométrie
Altitude d’ambiguı̈té
Interférommétrie
λDist sin(i)
IMesurer
l’élévation
Ea = 2Borth
Mesurer le
mouvement

polarimetrie
Exemple numérique : B orth : 100 m, Dist=890 km, i=23 degree, λ=5.6 cm
Ea=97 m
Précision : dépend de l’altitude d’ambiguı̈té (de 10m à 100m en fonction de la
base SM, et plus spécifiquement de la projection orthogonale Borth )

Sébastien Giordano  30-11-2018  108 / 108 Radar PPMD


Influence de l’altitude d’ambiguı̈té sur les franges

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
IMesurer
l’élévation
Mesurer le
mouvement

polarimetrie

24
24
Massonnet et Inglada
Sébastien Giordano  30-11-2018  109 / 109 Radar PPMD
Interférométrie différentielle

Introduction
Principe : suppression des franges topographiques
Géométrie Comment : utilisation d’un MNE ou d’une 3eme image SAR
des images
radar Résultat : franges résiduelles → mouvement
Radiométrie
Précision : inférieure au cm
Radargrammétrie
et Radarcli-
nométrie Point P vu d’une antenne en M (mesure de phase φM ). Le point P bouge en
0
Interférommétrie P , seconde acquisition depuis une antenne en S (mesure de phase φS ).
Mesurer
l’élévation
IMesurer le
0
mouvement mouvement P → P
polarimetrie 0
4πSP
φS = φgeom−S + φscatt−S = + φscatt−S (34)
λ
4πMP
φM = φgeom−M + φscattM = + φM−scatt (35)
λ
25

25
[Crosetto et al., 2016]
Sébastien Giordano  30-11-2018  110 / 110 Radar PPMD
Interférométrie différentielle

Dans ce cas le mesure de phase interférométrique est donnée par :

Introduction Phase interférométrique


Géométrie
des images 4π(SP 0 − MP)
radar ∆φint = φS − φM = + (φscatt−S − φscatt−M ) (36)
λ
Radiométrie
4πSP
Radargrammétrie
et Radarcli-
En ajoutant et retranchant le terme λ
l’équation devient:
nométrie

Interférommétrie Phase interférométrique


Mesurer
l’élévation
IMesurer le 0
mouvement SP − MP SP − SP
polarimetrie
∆φint = φS − φM = λ
+ λ
+ (φscatt−S − φscatt−M ) (37)
4pi 4pi

SP−MP
λ composante topographique φTopo

0
SP −SP
λ composante de déplacement φDispl

26

26
[Crosetto et al., 2016]
Sébastien Giordano  30-11-2018  111 / 111 Radar PPMD
Interférométrie Différentielle

Introduction

Géométrie
des images
radar

Radiométrie
En faisant l’hypothèse que φscatt−S − φscatt−M = 0. Avec un MNE on peut
Radargrammétrie
simuler φtopo . On obtient ∆φD−int , DInSAR phase.
et Radarcli-
nométrie
DInSAR phase
Interférommétrie
Mesurer
l’élévation ∆φD−int = ∆φInt − φToposimul = φDispl (38)
IMesurer le
mouvement
27
polarimetrie

27
[Crosetto et al., 2016]
Sébastien Giordano  30-11-2018  112 / 112 Radar PPMD
Franges résiduelles, mouvement

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
Mesurer
l’élévation
IMesurer le
mouvement

polarimetrie

Mise en évidence d’un soulèvement au sud de la gare St-Lazare à Paris. Une


frange (un cycle de couleur) représente 3 mm de déplacement vertical (Le
Mouélic et al., GRL 2002) 28
28
[Link]
Sébastien Giordano  30-11-2018  113 / 113 Radar PPMD
Franges résiduelles, mouvement

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
Mesurer
l’élévation
IMesurer le
mouvement

polarimetrie

Les déformations du Piton de la fournaise (Réunion) sont mises en évidence en


combinant la phase de deux images du satellite Radarsat. Une frange représente
un déplacement vertical de 3 cm 29
29
[Link]
Sébastien Giordano  30-11-2018  114 / 114 Radar PPMD
Interférométrie Différentielle
Jusqu’à présent équations simplifiées, en réalité :
Eq. (7) represents a simplified DInSAR observation equation : erreurs liées à
Introduction l’atmosphère, orbitographie, topographie
Géométrie
des images Phase interférométrique différentielle
radar

Radiométrie

Radargrammétrie ∆φD−Int = ∆φInt − φTopos imul


et Radarcli-
nométrie = φDispl + φTopores + φAtm−S − φAtm−M
Interférommétrie
Mesurer + φOrb−S − φOrb−M + φNoise + 2kπ (39)
l’élévation
IMesurer le
mouvement

polarimetrie
φTopores : residual topographic error, erreur précision MNT, φTopos imu
φOrb composante de la phase liée aux erreurs d’orbites (erreurs dans la
détermination de la position de M et S) retards atmosphériques
φAtm composante atmosphérique, date d’acquisition de chaque image
φNoise perte de cohérence : décorrélation spatiale et temporelle
2kπ k ambiguı̈té de phase, on connait le déphasage [2π]
30
30
[Crosetto et al., 2016]
Sébastien Giordano  30-11-2018  115 / 115 Radar PPMD
Mesure de déplacement

Introduction

Géométrie
des images
radar
Objectif
Radiométrie

Radargrammétrie Séparer φDispl des autres composantes.


et Radarcli-
nométrie

Interférommétrie
Limitations :
Mesurer
l’élévation
décorrelations temporels et géométrqiues → φnoise
IMesurer le
mouvement item the atmospheric component
polarimetrie déroulement de phase k
[Crosetto et al., 2016]

Sébastien Giordano  30-11-2018  116 / 116 Radar PPMD


Solution : Coregistration

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
La calcul de l’interférogramme est toujours précédé d’un calcul de coregistration
nométrie par corrélation et en préservant la phase: recalage de l’image slave sur master.
Interférommétrie Meilleur est la superpostion des images (par corrélation) plus faible est le rang
Mesurer
l’élévation d ’ambiguı̈té.
IMesurer le
mouvement

polarimetrie

Sébastien Giordano  30-11-2018  117 / 117 Radar PPMD


Solution : Permanent Scatterer Interferometry

Introduction

Géométrie
des images
radar
Utiliser des pixels caractérisés par un faible φNoise
Radiométrie
Permanent Scatterers PS
Radargrammétrie
et Radarcli-
nométrie
mécanisme de rétrodiffusion du signal radar dominé par un diffuseur
Interférommétrie
élémentaire unique.
Mesurer
l’élévation
phase interférométrique peu affectée par la décorrélation géométrique. Si
IMesurer le
mouvement
pas de changement dans le temps décorrélation temporelle négligeable.
polarimetrie bâtiments, structures métalliques, éléments géologiques naturels ..

[Crosetto et al., 2016]

Sébastien Giordano  30-11-2018  118 / 118 Radar PPMD


Solution : systèmes d’acquistion

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli- système d’acquisition : SRTM / Tandem-X
nométrie

Interférommétrie bases très courtes pour la mesure de déformation : Sentinel-1


Mesurer
l’élévation image Tandem-X image navette SRTM
IMesurer le
mouvement

polarimetrie

Sébastien Giordano  30-11-2018  119 / 119 Radar PPMD


Carte de cohérence

Introduction

Géométrie
des images S1 = A1 e jφ1 et S2 = A2 e jφ2 Après que les images sont recalées géométriquement
radar
: L’interférogramme complexe correspondant au couple est donné par le produit
Radiométrie
de S1 et du conjugué de S2 .
Radargrammétrie
et Radarcli-
nométrie Interférogramme
Interférommétrie PN
Mesurer (< s1 >n . < s2 ∗ >n )
l’élévation γ = p PNi=1 (40)
( i=1 (< s1 >n . N
P
i=1 (< s2 >n )
IMesurer le
mouvement

polarimetrie
N : nombre de pixels de la fenêtre de voisinage
|γ|2 : coefficient de corrélation appelé cohérence

Sébastien Giordano  30-11-2018  120 / 120 Radar PPMD


Carte de cohérence

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie
Mesurer
l’élévation
IMesurer le
mouvement

polarimetrie

(a) InSAR and (b) MAI interferograms; (c) coherence map and (d) the expected
MAI displacement accuracy map generated from ERS interferometric pair 31
31
[Link]
Sébastien Giordano  30-11-2018  121 / 121 Radar PPMD
Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie
polarimetrie
Interférommétrie

Ipolarimetrie

Sébastien Giordano  30-11-2018  122 / 122 Radar PPMD


Onde plane monochromatique

Onde plane
Introduction
Onde plane : solution aux equations de Maxwell dans le vide + à une grande
Géométrie
des images distance de l’antenne on peut considérer que le front d’onde est plan
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie



E résultat d’une composante horizontale et verticale


E prévisible : onde complètement polarisée

Sébastien Giordano  30-11-2018  122 / 122 Radar PPMD


Ellipse de polarisation

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

Ellipse de polarisation
A : amplitude
α : phase absolue
φ angle d’orientation, φ ∈ [− π2 ; π2 ]
τ ellipticité (ou execentricité), τ ∈ [0; π4 ]

Sébastien Giordano  30-11-2018  123 / 123 Radar PPMD


Polarisations remarquables

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie composantes Ex et Ey en
Interférommétrie phase : polarisation linéaire
Ipolarimetrie
∆φ = π2 ou 3π
2
: polarisation
circulaire

32 33
32
[Link]
33
Sébastien Giordanofiles/images/resource/tutor/polarim/chapter1/images/fig1
 30-11-2018  124 / 124 [Link] Radar PPMD
Émission système radar

antenne conçue pour transmettre et capter des ondes dans une


Introduction
polarisation définie
Géométrie
des images émission d’ondes de différentes polarisations séparément en utilisant un
radar
commutateur
Radiométrie

Radargrammétrie réception simultanée de différentes polarisations de l’onde


et Radarcli-
nométrie polarisations linéaires horizontale, H, et verticale, V les plus utilisées
Interférommétrie

Ipolarimetrie

[?]

Sébastien Giordano  30-11-2018  125 / 125 Radar PPMD


Polarisation

Configurations de polarisation des systèmes radar imageurs:


HH: émission horizontale , réception horizontale,
Introduction

Géométrie
HV : émission horizontale , réception verticale,
des images
radar VV : émission verticale, réception verticale,
Radiométrie VH : émission verticale , réception horizontale.
Radargrammétrie
et Radarcli-
nométrie Co-polarisation et Cross-polarisation
Interférommétrie
canaux copolarisés : canaux HH et VV .
Ipolarimetrie
canaux de cross-polarisation ou en polarisation croisée pour HV et VH.

34
34
illustration : [?]
Sébastien Giordano  30-11-2018  126 / 126 Radar PPMD
Matrice de diffusion [S]

Introduction Les propriétés de diffusion d’une cible sont entièrement contenues dans la ma-
Géométrie trice de diffusion [S].
des images
radar
Matrice de diffusion [S]
Radiométrie
 s
  EHi EHi
    
Radargrammétrie EH S SHV
et Radarcli- = S = HH (41)
nométrie EVs EVi SVH SVV EVi
Interférommétrie

Ipolarimetrie
< i > : onde incidente
< s > : onde diffusée (scattered)
chaque élément : amplitude et phase mesurées par les canaux du radar
polarimétrique
radar monostatique : SHV = SVH , [S] symétrique
contient l’intégralité de l’information polarimétrique.

Sébastien Giordano  30-11-2018  127 / 127 Radar PPMD


à déplacer

Chaque élément SXY (XY = {HH, HV , VH, VV }) est un nombre complexe.


Notation en amplitude |SXY | et phase φXY , et la notation en partie réelle iXY
Introduction
et partie imaginaire qXY .
Géométrie
des images
radar SXY = |SXY | e jφXY
Radiométrie
(42)
SXY =iXY + [Link]
Radargrammétrie
et Radarcli-
nométrie Pour chaque canal de polarisation SXY , IXY intensité et AXY amplitude du
Interférommétrie canal SXY .
Ipolarimetrie

IXY = |SXY |2 = iXY


2 2
+ qXY (43)
q
2 2
AXY = |SXY | = iXY + qXY (44)
Quand la même antenne radar réalise l’émission et la réception des ondes, le
radar est dit monostatique. Si deux antennes différentes sont utilisées pour
l’émission et la réception, le radar est dit bistatique. Dans le cas monostatique
(cas majoritaire dans le cas du radar satellite) on a SHV = SVH et par conséquent
la matrice [S] est symétrique.

Sébastien Giordano  30-11-2018  128 / 128 Radar PPMD


Matrice de Covariance [C ]

Introduction

Géométrie
des images
radar La matrice de covariance [C ] :
Radiométrie est une représentation en puissance de l’information polarimétrique,
Radargrammétrie
et Radarcli- est construite à partir du vecteur de diffusion kc .
nométrie

Interférommétrie
Vecteur de diffusion, cas monostatique
Ipolarimetrie
 
√SHH
kc =  2SHV  (45)
SVV

Sébastien Giordano  30-11-2018  129 / 129 Radar PPMD


Matrice de Covariance

[C ] une matrice hermitienne semi-définie positive


Introduction

Géométrie Matrice de covariance dans le cas monostatique


des images
radar
 ∗
√ ∗ ∗ 
√SHH .SHH∗ 2SHH .SHV SHH .SVV
Radiométrie

Radargrammétrie
et Radarcli-
C = kc .kc∗t =  2SHV .SHH √2SHV .SHV ∗
2SHV .SVV∗ 
∗ ∗ ∗
nométrie SVV .SHH 2SVV .SHV SVV .SVV
Interférommétrie  2
√ ∗ ∗ 
Ipolarimetrie √ |SHH | ∗ 2SHH SHV
2
√SHH SVV∗
=  2SHV SHH √2 |SHV | 2SHV SVV  (46)

SVV SHH ∗
2SVV SHV |SVV |2
(47)

On utilise traditionnellement les notations suivantes :

 
C11 C12 C13
C = C21 C22 C23  (48)
C31 C32 C33

Sébastien Giordano  30-11-2018  130 / 130 Radar PPMD


Matrice de cohérence

Introduction

Géométrie
des images La matrice de cohérence [C ] :
radar

Radiométrie
est une représentation en puissance de l’information polarimétrique,
Radargrammétrie est construite à partir du vecteur de Pauli kp
et Radarcli-
nométrie

Interférommétrie Vecteur de Pauli dans le cas monostatique


Ipolarimetrie  
S + SVV
1  HH
kp = √ SHH − SVV  (49)
2 2S HV

le vecteur de Pauli kp traduit mieux les types de rétrodiffusion.

Sébastien Giordano  30-11-2018  131 / 131 Radar PPMD


Matrice de cohérence

Matrice de cohérence dans le cas monostatique


Introduction
∗t
Géométrie [T ] = kp .kp
des images

(SHH + SVV )(SHH − SVV )∗ ∗ 

radar 1 (SHH + SVV )(SHH + SVV ) 2(SHH + SVV )SHV
= (SHH − SVV )(SHH + SVV )∗ (SHH − SVV )(SHH − SVV )∗ ∗ 
2(SHH − SVV )SHV

Radiométrie 2 2SHV (SHH + SVV ) 2SHV (SHH − SVV )∗ ∗
4SHV SHV (50)
Radargrammétrie
|SHH + SVV |2 (SHH + SVV )(SHH − SVV )∗ ∗ 

et Radarcli- 1 2(SHH + SVV )SHV
nométrie = (SHH − SVV )(SHH + SVV )∗ |SHH − SVV |2 ∗ 
2(SHH − SVV )SHV
2
2SHV (SHH + SVV )∗ 2SHV (SHH − SVV )∗ 4 |SHV |2
Interférommétrie

Ipolarimetrie

[C ] une matrice hermitienne semi-définie positive


On peut exprimer chaque Tij comme une combinaison linéaire de Cij .
[C ] et [T ] ont les mêmes valeurs propres;
structure des vecteurs propres de [T ] est plus appropriée pour l’analyse
polarimétrique
les termes diagonaux sont réels et sont liés à des mécanismes
simples de rétrodiffusion

Sébastien Giordano  30-11-2018  132 / 132 Radar PPMD


Estimation locale de [C ] ou [T ]

Introduction

Géométrie
des images
radar
Estimation ces matrices sur un voisinage local (spatial ou temporel).
Radiométrie
notation < [C ] >M×M et de même < Cij >M×M désigne un élément Cij
Radargrammétrie
de < [C ] >M×M .
et Radarcli-
nométrie Le nombre d’observations entrant en compte dans le calcul de ces
Interférommétrie paramètres moyens est appelé nombre de look ou nombres de vues
Ipolarimetrie Dans le cas de l’estimation locale sur un voisinage fixe, on utilise
généralement une taille [11 × 11] pour des données acquises en 1-look et
des tailles [7 × 7] pour des images 4-look [?]
voisinage locaux les plus communes allant de [7 × 7] à [11 × 11]
.

Sébastien Giordano  30-11-2018  133 / 133 Radar PPMD


Décomposition polarimétrique

Introduction

Géométrie
Théorèmes de décomposition polarimétrique
des images
radar Les théorèmes de décomposition polarimétrique ont pour but de décomposer
Radiométrie les représentations de l’information polarimétrique afin d’en retrouver les
Radargrammétrie mécanismes de rétrodiffusion.
et Radarcli-
nométrie

Interférommétrie Les théorèmes peuvent s’appliquer à la matrice de diffusion, à la matrice


Ipolarimetrie de covariance, à la matrice de cohérence [T ] ou à d’autres
représentations de l’information polarimétrique
Il existe de nombreux théorèmes de décomposition. État de l’art dans :
[Cloude and Pottier, 1996], [Lee, 2009] et [?].
Dans ce cours on présente 2 décompositions :
la décomposition de Paui
la décomposition de Cloude et Pottier

Sébastien Giordano  30-11-2018  134 / 134 Radar PPMD


Décomposition polarimétrique

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

Sébastien Giordano  30-11-2018  135 / 135 Radar PPMD


Décomposition de Pauli

Introduction La décomposition de Pauli s’appuie sur :


Géométrie
des images la décomposition de la matrice de diffusion [S] (décomposition cohérente)
radar
ou la matrice de cohérence [T ] (décomposition incohérente) sur une base
Radiométrie
simple de mécanismes de rétrodiffusion.
Radargrammétrie
et Radarcli-
nométrie
Les mécanismes de rétrodiffusion considérés sont la réflexion de surface,
Interférommétrie
la réflexion paire (ou de doubles rebonds) et la réflexion de volume.
Ipolarimetrie Les matrices de Pauli sont utilisées pour la décomposition de la matrice
de diffusion.
Matrices de Pauli
       
1 0 1 0 0 1 0 −j
σ0 = σ1 = σ2 = σ3 = (51)
0 1 0 −1 1 0 j 0

Sébastien Giordano  30-11-2018  136 / 136 Radar PPMD


Décomposition de Pauli

Evi
Introduction

Géométrie Ehi
des images
radar

Radiométrie

Radargrammétrie Ehs Ev
et Radarcli-
nométrie Evs Eh
Interférommétrie

Ipolarimetrie Ehs
Evs

Evi
Ehi

Figure: Réflexion de surface Figure: Réflexion paire, doubles


rebonds

Sébastien Giordano  30-11-2018  137 / 137 Radar PPMD


Décomposition de Pauli

En utilisant la base (σ0 , σ1 , σ2 ) des matrices de Pauli, on peut décomposer la


Introduction
matrice de diffusion.
Géométrie
des images
radar Décomposition de la matrice de diffusion avec les bases de Pauli
Radiométrie
       
Radargrammétrie S SHV SHH + SVV 1 0 SHH − SVV 1 0 0 1
et Radarcli- S = HH = + + SHV
nométrie SVH SHV 2 0 1 2 0 −1 1 0
Interférommétrie (52)
Ipolarimetrie

On retrouve les intensités des composantes du vecteur de Pauli dans la définition


de la matrice de cohérence [T ].

Matrice de cohérence
|SHH + SVV |2 (SHH + SVV )(SHH − SVV )∗ ∗ 

1 2(SHH + SVV )SHV
T = (SHH − SVV )(SHH + SVV )∗ |SHH − SVV |2 ∗ 
2(SHH − SVV )SHV (53)
2
2SHV (SHH + SVV )∗ 2SHV (SHH − SVV )∗ 4 |SHV |2

Sébastien Giordano  30-11-2018  138 / 138 Radar PPMD


Visualisation

Introduction

Géométrie
des images
radar

Radiométrie La décomposition de Pauli est communément représentée sous la forme d’une


Radargrammétrie composition colorée avec les conventions suivantes :
et Radarcli-
nométrie 2
Interférommétrie Canal Bleu : diffusion de surface ( |SHH +S
2
VV |
)
Ipolarimetrie Canal Vert : diffusion de volume (2|SHV |2 )
2
Canal Rouge diffusion de doubles rebonds ( |SHH −S
2
VV |
)

Sébastien Giordano  30-11-2018  139 / 139 Radar PPMD


Mécanismes simples de rétrodiffusion

Diffusion de surface (σ0 )


Introduction
représente une interaction simple entre l’onde électromagnétique incidente et
Géométrie une surface plane. composantes horizontales et verticales subissent un
des images
radar
déphasage de π ce qui se traduit par EHs = −EHi et , EVs = −EVi . La forme
Radiométrie
canonique de la matrice de diffusion [S] est donnée par [S] = −σ0 .
Radargrammétrie
et Radarcli-
nométrie Diffusion de doubles rebonds (σ1 )
Interférommétrie
représente l’interaction d’une onde électromagnétique avec un dièdre. Dans le
Ipolarimetrie
cas d’une réflexion paire (deux interactions entre l’onde et l’objet) la
composante horizontale est conservée EHs = EHi alors que seule la composante
verticale subit un déphasage de π : EVs = −EVi . La forme canonique
[S] = −σ1 de la matrice de diffusion décrit ce comportement.

Diffusion de volume (σ2 )


La forme canonique [S] = −σ2 illustre la diffusion de volume. Il s’agit de
multiples interactions entre l’onde incidente et l’objet. Dans ce cas beaucoup
d’énergie revient dans les canaux de polarisation croisée.

Sébastien Giordano  30-11-2018  140 / 140 Radar PPMD


Décomposition de Pauli

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

Sébastien Giordano  30-11-2018  141 / 141 Radar PPMD


Décomposition de Pauli

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

Sébastien Giordano  30-11-2018  141 / 141 Radar PPMD


Introduction

Géométrie
des images
radar
1 : plan d’eau 2 : jeune maı̈s 3 : galets Durance
Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

4 : culture orge 5 : zone urbaine 6 : prairie sainfoin

7 : forêt 8 : culture colza


Sébastien Giordano  30-11-2018  142 / 142 Radar PPMD
Décomposition de Pauli

Introduction

Géométrie
des images Table: Paramètres moyens de la décomposition de Pauli (d.n.) de l’image
radar RS2-FORCA-003 pour chacune des huit zones du site d’Oraison
Radiométrie

Radargrammétrie
et Radarcli-
nométrie Num. Type E (< T11 >7×7 ) E (< T22 >7×7 ) E (< T33 >7×7 )
Interférommétrie
1 plan d’eau 0, 16 × 106 0, 10 × 106 0, 07 × 106
Ipolarimetrie 2 jeune maı̈s 2, 08 × 106 0, 31 × 106 0, 26 × 106
3 galets 5, 92 × 106 1, 04 × 106 0, 97 × 106
4 culture d’orge 1, 91 × 106 1, 83 × 106 0, 46 × 106
5 zone urbaine 1, 14 × 107 1, 71 × 107 0, 98 × 106
6 prairie de sainfoin 1, 72 × 106 1, 15 × 106 0, 89 × 106
7 forêt 2, 77 × 106 1, 09 × 106 1, 09 × 106
8 culture colza 4, 55 × 106 3, 36 × 106 3, 12 × 106

Sébastien Giordano  30-11-2018  143 / 143 Radar PPMD


Décomposition de Cloude et Pottier

Introduction
La décomposition incohérente de Cloude et Potter [Cloude and Pottier, 1997]
Géométrie
des images s’appuie sur une décomposition en éléments propres de la matrice de cohérence
radar
[T ].
Radiométrie

Radargrammétrie Décomposition incohérente de Cloude et Pottier


et Radarcli-
nométrie

Interférommétrie
   ∗t
λ1 0 0
Ipolarimetrie
[T ] = u1 u2 u3   0 λ2 0  u1 u2 u3 
0 0 λ3
u1 , u2 , u3 base de vecteurs propres orthogonaux (54)
jδi jγi T
ui de la forme ui = [cos(αi ), sin(αi ) cos(βi )e , sin(αi ) sin(βi )e ]
λ1 , λ2 , λ3 valeurs propres réelles

Sébastien Giordano  30-11-2018  144 / 144 Radar PPMD


Décomposition de Cloude et Pottier : Entropie Polarimétrique

À partir de la décomposition en éléments propres trois paramètres sont définis


: l’entropie polarimétrique, l’anisotropie polarimétrique et l’angle α.
Introduction

Géométrie Entropie polarimétrique H


des images
radar
Mesure du degré de désordre des ondes réfléchies : H ∈ [0, 1]. Si H = 0.
Radiométrie
Dans ce cas un seul mécanisme de rétrodiffusion est à l’œuvre dans la cellule
Radargrammétrie
et Radarcli- de résolution et seulement λ1 est non nul si λ1 > λ2 > λ3 . A l’inverse si
nométrie
H = 1 les ondes reçues ont été entièrement dépolarisées et λ1 = λ2 = λ3 = 13 .
Interférommétrie
L’entropie peut aussi être vue comme une mesure de la qualité des
Ipolarimetrie
interprétations qui peuvent être réalisées à partir d’une décomposition
polarimétrique. Il est communément admis que quand H > 0, 7 il est difficile
d’identifier les mécanismes de rétrodiffusion.

Entropie polarimétrique

3
X λi
H=− Pi . log3 (Pi ) avec Pi = P3 (55)
k=1 k=1 (λk )

Sébastien Giordano  30-11-2018  145 / 145 Radar PPMD


Décomposition de Cloude et Pottier : Anisotropie Polarimétrique

Introduction
Anisotropie polarimétrique A
Géométrie
des images
radar
L’anisotropie polarimétrique A (56) est comme l’entropie polarimétrique
Radiométrie
définie à partir des valeurs propres λi de la décomposition en éléments propres
Radargrammétrie
de [T ] et par définition A ∈ [0, 1]. L’anisotropie polarimétrique compare la
et Radarcli-
nométrie
deuxième et la troisième valeur propre. L’anisotropie est une information
Interférommétrie
complémentaire à l’entropie. Particulièrement quand H ∈ [0, 7, 0, 9], A → 1
Ipolarimetrie
indique qu’un deuxième mécanisme de diffusion important intervient. Si
A → 0 il n’y a pas de comportement secondaire (λ2 ≈ λ3 ).

Anisotropie polarimétrique A

λ2 − λ3
A= (56)
λ2 + λ3

Sébastien Giordano  30-11-2018  146 / 146 Radar PPMD


Décomposition de Cloude et Pottier : Angle α

Introduction

Géométrie
des images Angle α
radar

Radiométrie
Angle α L’angle α est déterminé à partir des valeurs propres et des vecteurs
Radargrammétrie
propres de la décomposition en éléments propres de [T ]. L’angle α est un
et Radarcli-
nométrie
angle ∈ [0, 90◦ ] qui traduit le type de mécanisme de rétrodiffusion dominant
Interférommétrie
dans la cellule de résolution. α = 0 caractérise la diffusion par une surface,
Ipolarimetrie
α = 45◦ la diffusion par un volume, α = 90◦ la diffusion par un dièdre.

Angle α

α = P1 α1 + P2 α2 + P3 α3 (57)

Sébastien Giordano  30-11-2018  147 / 147 Radar PPMD


Décomposition incohérente de Cloude et Pottier

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

Entropie H Anisotropie A Angle α

0 Entropie 1
0 Anisotropie 1
0 Angle α 90◦

Figure: Paramètres de [Cloude and Pottier, 1997] : Entropie H, Anisotropie A et


Angle α, estimés sur un voisinage 7 × 7

Sébastien Giordano  30-11-2018  148 / 148 Radar PPMD


Paramètres polarimétriques de Cloude et Pottier

Introduction

Géométrie
des images Table: Moyenne des paramètres (H, A, α) de [Cloude and Pottier, 1997] pour les 8
radar zones du site d’Oraison
Radiométrie

Radargrammétrie Num. Type H̄ Ā ᾱ


et Radarcli-
nométrie
1 plan d’eau 0,91 0,26 49,43
Interférommétrie
2 jeune maı̈s 0,57 0,24 23,57
Ipolarimetrie
3 galets Durance 0,64 0,27 27,18
4 culture orge 0,68 0,34 48,19
5 zone urbaine 0,58 0,67 52,26
6 prairie sainfoin 0,88 0,28 48,62
7 forêt 0,86 0,24 42,26
8 culture colza 0,92 0,23 51,07

Sébastien Giordano  30-11-2018  149 / 149 Radar PPMD


Segmentation du plan (H, α)

Introduction
Table: Types Z 1 à Z 9 (définitions et codes couleurs) proposés par
Géométrie
des images
[Cloude and Pottier, 1997] correspondant à des régions du plan (H, α)
radar

Radiométrie

Radargrammétrie
et Radarcli- Z1 High Entropy Multiple Scattering (branch/crown structure)
nométrie
Z2 High Entropy Vegetation Scattering (cloud of anisotropic needles)
Interférommétrie

Ipolarimetrie
Z3 High Entropy Surface Scatter (no feasible region)
Z4 Medium Entropy Multiple Scattering (forestry double bounce)
Z5 Medium Entropy Vegetation Scattering (vegetation)
Z6 Medium Entropy Surface Scatte (surface roughness)
Z7 Low Entropy Multiple Scattering Events (dihedral scatterer)
Z8 Low Entropy Dipole Scattering (dipole)
Z9 Low Entropy Surface Scatter (Bragg surface)

Sébastien Giordano  30-11-2018  150 / 150 Radar PPMD


Segmentation du plan (H, α)

���������������������������������������������������������������������������

Introduction
���
��������������
Géométrie ���
des images
radar
��� �� ��
Radiométrie ��
���
Radargrammétrie
et Radarcli-
�� � �
nométrie �
� �
�� �� ��

Interférommétrie
��� �
Ipolarimetrie
��� ��


��� �� ��

���

��
�� ���� ���
� ���� ���� ���� ���� ���� ���� � ��� ��

Figure: Plan (H, α), régions empiriques du plan et représentation des valeurs de
(H, α) pour chaque type d’occupation du sol du site d’Oraison

Sébastien Giordano  30-11-2018  151 / 151 Radar PPMD


PolInSAR

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

Sébastien Giordano  30-11-2018  152 / 152 Radar PPMD


Bibliographie

Introduction

Géométrie
des images
radar

Radiométrie

Radargrammétrie
et Radarcli-
nométrie

Interférommétrie

Ipolarimetrie

[Massonnet, 2008], [Maı̂tre, 2001]

Sébastien Giordano  30-11-2018  153 / 153 Radar PPMD


Merci pour votre attention

Question(s) ?

[Link]@[Link]
Bibliographie I

Cloude, S. R. and Pottier, E. (1996).


Introduction A review of target decomposition theorems in radar polarimetry.
Géométrie Geoscience and Remote Sensing, IEEE Transactions on, 34(2):498–518.
des images
radar

Radiométrie
Cloude, S. R. and Pottier, E. (1997).
Radargrammétrie
An entropy based classification scheme for land applications of
et Radarcli- polarimetric sar.
nométrie
Geoscience and Remote Sensing, IEEE Transactions on, 35(1):68–78.
Interférommétrie

Ipolarimetrie Crosetto, M., Monserrat, O., Cuevas-González, M., Devanthéry, N., and
Crippa, B. (2016).
Persistent scatterer interferometry: A review.
ISPRS Journal of Photogrammetry and Remote Sensing, 115:78–89.
Deledalle, C.-A., Denis, L., Tupin, F., Reigber, A., and Jager, M. (2015).
Nl-sar: A unified nonlocal framework for resolution-preserving (pol)(in)
sar denoising.
Geoscience and Remote Sensing, IEEE Transactions on,
53(4):2021–2038.

Sébastien Giordano  30-11-2018  155 / 155 Radar PPMD


Bibliographie II

Introduction
FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture,
Géométrie
. (1991).
des images Principes d’imagerie radar.
radar

Radiométrie Goodman, J. W. (1975).


Radargrammétrie Statistical properties of laser speckle patterns.
et Radarcli-
nométrie Laser speckle and related phenomena, Vol. 9:9–75.
Interférommétrie
Goodman, J. W. (1976).
Ipolarimetrie
Some fundamental properties of speckle.
J. Opt. Soc. Am., 66(11):1145–1150.
Lee, J.-S. (1981).
Refined filtering of image noise using local statistics.
Computer graphics and image processing, 15(4):380–389.
Lee, P. (2009).
Polarimetrics Radar Imaging, From Basic to Applications.

Sébastien Giordano  30-11-2018  156 / 156 Radar PPMD


Bibliographie III

Introduction

Géométrie
des images
radar Maı̂tre, H., editor (2001).
Radiométrie Traitement des images de radar à synthèse d’ouverture.
Radargrammétrie Traite IC2. Hermès science publications, Paris.
et Radarcli-
nométrie
Massonnet, S. (2008).
Interférommétrie
Imaging with Synthetic Aperture Radar.
Ipolarimetrie

Vasile, G., Trouvé, E., Lee, J., and Buzuloiu, V. (2006).


Intensity-driven adaptive-neighborhood technique for polarimetric and
interferometric sar parameters estimation.
IEEE Transactions on Geoscience and Remote Sensing, 44(6):1609–1621.

Sébastien Giordano  30-11-2018  157 / 157 Radar PPMD

Vous aimerez peut-être aussi