Imagerie Radar : Techniques et Applications
Imagerie Radar : Techniques et Applications
PPMD
Sébastien Giordano
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
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Introduction
polarimetrie
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)
1
1
[Link] électromagné[Link]
Sébastien Giordano 30-11-2018 3/3 Radar PPMD
Bandes Radar et Systèmes Radar
Bandes radar :
Satellite λ
Terrasar-X 3.1 cm
Radarsat-2 5.5 cm
ALOS-PALSAR-2 24 cm
Biomass (2020) 69 cm
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
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
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
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
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
IIntroduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
IIntroduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
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
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
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
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
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
polarimetrie
polarisation
phase
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
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 )
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
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
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
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
Interférommétrie
Remarques sur la résolution
polarimetrie
radiale (× 2)
Comment améliorer la résolution
radiale ?
Introduction
polarimetrie
Solution : la compression d’impulsion
émission d’un signal fortement modulé en fréquence : 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
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
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
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
Introduction
Radiométrie
Radargrammétrie
et Radarcli-
nométrie (Loading)
Interférommétrie
polarimetrie
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
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
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
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
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
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.
Introduction
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
Introduction
Géométrie
des images
radar
Résolution
radiale
Résolution
azimutale
IDistorsions
géométriques
Interférommétrie
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
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
Introduction
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
Introduction
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
polarimetrie
de la polarisation,
de l’angle d’incidence
propriétés géométriques et géophysiques de la cible
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
Radargrammétrie
et Radarcli-
nométrie
Figure: Angle d’incidence d’après [?]
Interférommétrie
polarimetrie
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
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
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
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
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
polarimetrie
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
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
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
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 |
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.
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
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.
Introduction
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
polarimetrie
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
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 [−π, π].
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
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)
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
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
Radargrammétrie 1 − RI
et Radarcli- p(IXY |R) = e I ≥0 (19)
nométrie R
Interférommétrie
polarimetrie
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
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
�������
������� ������
������ ��������
����� �������
�� ������
� �
������ �� �� � � �������
������� ������
������ �������
����� ������
���������� ���
� �
���
�����
�� � �
� � �
�
�� ����� ����� ����� ����� � ����� �
� ������� ����� � ������� � ����� � ������� � ����� � ������� � ������
� � ����� ����� ����� ����� ������ �� �������
�
�����
�
�������
�
�����
�
�������
�
�����
�
�������
�
�����
��� ���
��� ���
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
Interférommétrie
polarimetrie N X M
1 X
< C >N×M = C (i, j) (21)
N × M i=1 j=1
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
Interférommétrie
NA2
2NA2N−1
NXY
NXY
p(ANXY |R) = e− R (25)
Γ(N)R N
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
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
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
81-looks
polarimetrie
polarimetrie
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 ?
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é.
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.
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
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
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
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
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
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
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
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
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 ?
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
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
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
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
polarimetrie Propriété
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
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
information topographique
Interférommétrie
IMesurer
l’élévation
Mesurer le
mouvement
polarimetrie
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 )
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
SP−MP
λ composante topographique φTopo
4π
0
SP −SP
λ composante de déplacement φDispl
4π
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
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
Radiométrie
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
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]
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
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 ..
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli- système d’acquisition : SRTM / Tandem-X
nométrie
polarimetrie
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
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
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
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 ]
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
Ipolarimetrie
[?]
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.
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
C11 C12 C13
C = C21 C22 C23 (48)
C31 C32 C33
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
Ipolarimetrie
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]
.
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
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Ipolarimetrie
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
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
Introduction
Géométrie
des images
radar
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Ipolarimetrie
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Ipolarimetrie
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
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
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
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
Entropie polarimétrique
3
X λi
H=− Pi . log3 (Pi ) avec Pi = P3 (55)
k=1 k=1 (λk )
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
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)
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Ipolarimetrie
0 Entropie 1
0 Anisotropie 1
0 Angle α 90◦
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
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)
���������������������������������������������������������������������������
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
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Ipolarimetrie
Introduction
Géométrie
des images
radar
Radiométrie
Radargrammétrie
et Radarcli-
nométrie
Interférommétrie
Ipolarimetrie
Question(s) ?
[Link]@[Link]
Bibliographie I
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.
Introduction
FAO(Organisation des Nations Unies pour l’Alimentation et l’Agriculture,
Géométrie
. (1991).
des images Principes d’imagerie radar.
radar
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