Introduction Imagerie Spatiale

Imagerie spatiale
Introduction
Cours Qualité Image : Introduction 1/90 Adapté des slides CNES/ cours QI V5.4 © CNES 2018
Objectifs du cours
• Mesurer l’utilité de la télédétection spatiale
• Introduire les particularités de la prise de vue spatiale
• Comprendre le fonctionnement d’un système d’observation spatial
• Exposer la notion de qualité image, afin de savoir :
 Evaluer les performances d’un système d’observation
 Dimensionner un système d’observation qui répond au besoin
Cours Qualité Image : Introduction 2/90 V5.4 © CNES 2018

Utiliser une image brute ?

?
La prise de vue
n’est donc pas
4
« parfaite » ?
Le satellite imageur
est un système physique réel
On sait décrire son comportement 
5
Position, vitesse
attitude
Image brute
=
F(Paysage observé)
Optique, distorsion,
détection
Atmosphère
Eclairement
Paysage
observé
Conditions de prise de vue réelles
Acquisition
bord
Produit utilisateur
Paysage
observé image brute
Traitements
Bord/Sol
Modèle de prise de vue
132 211 453 514
350 530 382 361
sortie instrument
=
• Modèle physique qui relie le paysage observé à la
652
412 153
286 501
72
340
301 311
423
291 502
175 245
164 653
231 275 341 98

341 506 401 428
184 311 173 72

401 364 746 271
738 352 462 215
193 801 508 96
8
Vue d’ensemble
Utilisation de la télédétection spatiale

Qu’est-ce que la télédétection ?
Diversité spectrale
Diversité spatiale

L’observation de la Terre depuis l’espace
image Vis MeteoSat - Eumetsat image Radarsat - ASC
image thermique AVHRR - NOAA Vue 3D Pléiades de la

calotte du Tungnafelljokull
© CNES 2013, Distribution Airbus D&S
Synthèse 1998-2007 concentration en phytoplancton
Venise par Pleiades

Sent-2 - ESA © CNES 2013, Distribution Airbus D&S

Utilité de l’observation de la Terre
Suivi des catastrophes naturelles
Images Spot 4 avant et après la tempête Xynthia qui a frappé les côtes
vendéennes dans la nuit 27 au 28 février 2010. © CNES Distribution Spot Image
Haïti Cyclone Sandy
Cours Qualité Image : Introduction 12/90

© CNES 2012, distribution Astrium Services / Spot Image SA, all rights reserved
V5.4 © CNES 2018
© CNES 2012, distribution Astrium Services / Spot Image SA, all rights reserved
Renseignement
Cartographie
Lybie Benina airfield BDA, 1986

Mesure du relief
© CNES 2007/Distribution Spot Image

Traitement CNRS-LEGOS.
Qu’est-ce que la Télédétection ?
C’est l’ensemble des TECHNIQUES utilisées pour DETERMINER A

DISTANCE les PROPRIETES d’objets, naturels ou artificiels, à partir
des RAYONNEMENTS qu’ils émettent ou réfléchissent.
• Pierre angulaire : étude des interactions Matière-Rayonnement
Objet Rayonnement Capteur
Emission Transmission Réception

Ondes électromagnétiques
OPTIQUE
RADAR

Domaines spectraux en optique
froid
PHR rouge/vert/bleu PHR PIR/rouge/vert
SPOT 4 MIR/PIR/Vert
chaud
neige nuages
InfraRouge AVHRR nov 2012
Bandes spectrales
Bandes optiques : désignation courante des bandes
Code abrégé Bande spectrale Longueurs d’onde

PA ou PAN Panchromatique 0.45 – 0.70 à 0.90 µm
B0 Bleu 0.40 – 0.55 µm
B1 Vert 0.50 – 0.60 µm
B2 Rouge 0.60 – 0.70 µm
B3 ou PIR / anglais NIR Proche Infrarouge 0.75 – 1.30 µm
MIR / anglais SWIR Moyen Infrarouge 1.30 – 3.00 µm
IRT / anglais TIR Infrarouge thermique 3 – 50 µm
Bandes Radar : désignation par des lettres K, X, S, C, L, P

Exemples de réponses spectrales
Réponse paysages et instrumentales PHR

300
250
Mer
Luminances spectrales
200 Neige_a8
Vgt_herbe1
Sable_dune
150 B0
B1
B2
100
B3
PAN
50
0
0,4 0,5 0,6 0,7 0,8 0,9
microns

Optique ou Radar ?
OPTIQUE : PASSIVE RADAR : ACTIF

capteur
émetteur:
Illumination Écho réfléchi
naturelle
Impulsion
transmise

Imagerie Radar : système « tout temps »
Waterford, Ireland, 9/08/91
Surface : 50 km x 50 km
Landsat pass time:

10H43 a.m.
ERS-1 pass time:

11H25 a.m.
ERS (radar C band) Landsat TM
Le radar actif permet d’imager la Terre :

• de jour comme de nuit
• quelles que soit les conditions météorologiques (sans nuages ! )

Diversité spectrale : SPOT (2/5)
B1 ~vert B2 ~rouge B3 ~proche IR
Composition colorée :
Bleu = B1
Vert = B2
Rouge = B3

Diversité spectrale : Sentinel-2
SWIR
VNIR
Visible
Aerosols Water vapor Cirrus
60 m Snow/ice/cloud
B9
B1 B10
Vegetation
20 m red edge
B5 B7
B6 B8a B11 B12
10 m
B2 B3 B4 B8
400 600 800 1000 1200 1400 1600 1800 2000 2200 2400
nm nm nm nm nm nm nm nm nm nm nm

Transfert radiatif : modélisation du signal observé
Réflectance spectrale de la végétation chlorophyllienne
Visible Proche Moyen Infra-Rouge

Infra-Rouge
REFLECTION
pigment Structure
de la feuille Teneur en eau
cellulaire
ABSORPTION
Chlorophylle Eau
Longueur d’onde (mm)

Transfert radiatif : modélisation du signal observé
Pics d'absorption propres à chaque type de molécules présent
dans l'atmosphère  notion de fenêtres atmosphériques
1,0
0,9 O2 H2O
Transmission spectrale.
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0,0
0,4 0,9 1,4 1,9 2,4
Longueur d'onde (µm)
1,0
0,9
CO2 O3
Transmission spectrale.
0,8
0,7
0,6
0,5
0,4
0,3
0,2
0,1
0,0
3 4 5 6 7 8 9 10 11 12
Longueur d'onde (µm)

Imagerie  spectroscopie
• Résolution spatiale: • Résolution spectrale:

– 1m – 1000 bandes entre 3 et 15 µm
• Echantillonnage fin dans le • Echantillonnage fin dans le

domaine spatial domaine spectral
• Pourquoi ? • Pourquoi ?
– Cartographier – Remonter à des propriétés
– Identifier physiques
Imagerie
hyperspectrale
PLEIADES Spectro-imageur IASI
Diversité spatiale
• Grande variété des échelles d ’observation :

 pixel ~ 1 à 100 km : bilan radiatif, météorologie, atmosphère
 échelle globale ou continentale
 pixel ~ 100 à 1000 m : agriculture, forêt, océanographie,
environnement
 échelle régionale ou continentale
 pixel ~ 10 à 100 m : agriculture, cartographie, géologie, gestion des
risques
 échelle régionale ou locale
 pixel ~ 1 à 10 m : cartographie de précision, urbanisme, forêts
 échelle locale
 pixel ~ 1m à 0.20 m : haute et très haute résolution, renseignement
 Identification d’objets

Image Météo
résolution 4 km
champ 12000 km

Image Végétation
résolution 1 km
champ 2000 km
composition
colorée
B3 B2 B0

Image SPOT4
résolution 10 m
champ 60 km

Image SPOT5
résolution 2.5 m
champ 60 km

Image Pléiades
résolution 70 cm
CoursCours
Qualité
QI ©Image
CNES –:Introduction
Introduction 32/90
- 32/88 V5.4 © CNES 2018
Image WorldView2
résolution 50 cm

Image aéroportée
résolution 10 cm

Tibesti vu par PHR1B
CoursCours
Qualité
QI ©Image
Introduction 35/90
- 35/88 V5.4 © CNES 2018
Les pyramides vues par PHR1B
CoursCours
Qualité
QI ©Image
Introduction 36/90
- 36/88 V5.4 © CNES 2018
Manhattan vu par PHR1B
CoursCours
Qualité
QI ©Image
Introduction 37/90
- 37/88 V5.4 © CNES 2018
Le système spatial
de télédétection optique
Les orbites
Acquisition des images
Intérêt de l’imagerie par satellite
Conception d’un système d’observation

Acquisition des images : Illustration Pléiades

Observer : où, pourquoi et comment ?
• Accès mondial ou régional : quelles zones veut-on observer ?
• Couverture : quelle zone veut-on acquérir depuis une orbite donnée ?
• Répétitivité : combien de temps entre 2 images de la même zone ?
• Réactivité : combien de temps entre la demande et la livraison ?
• Résolution spatiale : quels détails sont visibles ?
• Résolution radiométrique : quelle précision sur la mesure physique ?
• Bandes spectrales : quelles couleurs/ longueurs d’onde/ thématiques ?
• Localisation : quelle précision de positionnement des pixels ?

Généralités sur les orbites satellitaires
• Trajectoires képlériennes (terre sphérique et homogène) :
 elles sont planes
 ce sont des coniques de foyer le centre attracteur (centre Terre)
 leur loi horaire est décrite par la loi des aires :
2 2
T 4
avec m, constante de gravitation terrestre : m = 398600 km3/s2 3 
T est la période et a le demi grand axe de l’orbite a m
• Orbites utilisées pour l’observation de la Terre :
 circulaires (ou faiblement elliptiques) => altitude h et vitesse
inertielle constantes
130
En orbite basse (h~800 km) :
Période orbitale en min.
120
Période T ~ 1 h 40
110
Vitesse orbitale ~ 7 km /s
100
90
a = req + h
80 req = 6378 km
0 200 400 600 800 1000 1200 1400 1600 1800 2000
altitude en km

Les orbites pour l’observation de la Terre
• Satellite Polaire
• Satellite Géostationnaire
 Couverture quasi-mondiale
 Accès permanent sur hémisphère
 Accessibilité limitée (selon
 Altitude élevée (36000 km) période, débattement, latitude)
 Altitudes faibles possibles
1 tour ~1h40
1 tour = 24h
1 tour/jour ~15 tours/jour

Fixe % Terre Vsubsat ~7 km/s
Résumé : choix d’une orbite
• Durée de vie, stabilité
 altitude > 300 km
• Haute résolution géométrique
 basse altitude
• Rapidité d’acquisition (visibilités station et site visé)
 haute altitude
• Échelle constante
 orbite circulaire
• Minimisation des variations d’éclairement
 orbite héliosynchrone
• Maximisation de la surface couverte
 orbite quasi polaire (i ~ 90°)
• Répétitivité
 phasage de l’orbite (notion de cycle)
• Couverture complète
 nombre important d’orbites (cycle assez long) ou large fauchée
Avion ou satellite ?
Règlementation aérienne / vulnérabilité Accès global / immunité

Choix de la trajectoire % site Trajectoire orbitale figée
Très haute résolution accessible Résolution limitée physiquement
Qualité géométrique variable Géométrie stable

Principes d’acquisition
1 pixel 1 ligne 1 image
scanners push-broom matrices
Contribution NASA
http://earthobservatory.nasa.gov/Features/EO1/

Organisation générale système
Demandes et distribution de produits
Données images
Information-demande de mise Station principale

en œuvre du système de Kiruna
(Suède)
Centre de contrôle
et programmation
Utilisateurs
Station de
Utilisateurs réception directe
Utilisateurs
Centre utilisateur de Toulouse

Panorama de quelques missions
d’observation de la Terre

RADARSAT II (Canada, 2001)
Mode pleine
polarimétrie
Bande C (5,3 GHz)

Période : 24 jours / Enregistrement embarqué
Résolution au sol : 3 m-100m Fauchée : 10 km-500 km
Pleine polarimétrie : res. (distance sol): 25 m x 28 m
LANDSAT 8
Masse : 2600 kg
Tir : 11/02/2013
Satellite civil américain (NASA/NOAA/USGS)

Acquisition par balayage push-broom, pas de dépointage
Bandes : PA 15 m, VIS-PIR-MIR (9 bandes) 30m, IRT (2 bandes) 100m Fauchée : 185 km
Orbite : cycle 16 jours Heure locale : 10:30 Altitude: 705 km
Codage 12 bits
Stockage des données: mémoire de masse
SENTINEL-2
Masse : 1200 kg
Tir S2A : 23/06/2015
Tir S2B : 06/03/2017
2 Satellites civils européens (ESA)

Acquisition par balayage push-broom,
13 bandes : 4 VNIR 10m, 6 PIR+SWIR 20m, 3 Atmosphère 60m Fauchée : 290 km
Orbite : cycle 10 jours (5j S2A+S2B) Heure locale : 10:30 Altitude: 786 km
Codage 12 bits
Stockage des données: mémoire de masse
Masse : 1000 kg
PLEIADES
Tir : 1A 17/12/2011,
1B 2/12/2012
Satellite civil/militaire (CNES) Lancement : 2011 et 2012

Bandes spectrales et résolution PAN 0,7m B0-B1-B2-B3 2,8 m
Orbite : cycle 26 jours Heure locale : 10:15 Altitude: 698 km
Fauchée : 20 km Codage 12 bits Masse : 1000 kg
Instrument : IPV : Focale 12,90 m, Diamètre 65 cm, PAN TDI, XS CCD quadrichrome
Stockage des mémoire à état solide (450 Gb)
données :
GEOEYE-2
Mass: 2050 kg
Stocké pour lancement 2016
US civil satellite, TDI Push-broom acquisition

Bands: PE 0.34m or 0.25m, B0-B1-B2-B3 1.36m or 1.00m Swath: 14.5km or 10.5km (Nadir)
Local time: 10:30 Altitude: 681 km or 500km
11-bit encoding
Data storage: solid-state memory

WORLDVIEW 3
Mass: 2800 kg
Launch: 13/08/2014
US civil satellite (DigitalGlobe), TDI Push-broom acquisition

Bands: PE 0.31 m XS : 8bands at 1.24 SWIR : 20 bands 3.7m Swath: 13.1 km (Nadir)
Orbit: revisit interval of 1 days <1m Local time: 13:30 Altitude: 617 km
11-bit encoding & 14 bit SWIR
Data storage: solid-state memory

La Qualité des images

Où est le pixel?
Géométrie
Quels détails sont visibles ?

Résolution
Que signifie la valeur du pixel ?

Radiométrie

Fil conducteur de la Qualité Image :
Comprendre et mesurer les défauts
de la prise de vue pour pouvoir les corriger
erreur de
localisation rayures
flou
bruit
QI QI
distorsions Couleurs altérées

La chaine image bord/sol
Prise de vue à bord
Géométrie Radiométrie Traitements bord
• orbite • FTM • prétraitements Stockage
Paysage • attitude • détection • compression
• visées • bruit • formattage
• relief • quantification • chiffrement
Transmission
Qualité
Géométrie Radiométrie Pré-Traitements

• corrections • égalisation • Déchiffrement
• déconvolution • Déformattage
Produit capteur
• projection • débruitage • Décompression
image • Effet du relief • fusion
Corrections au sol

1er Atelier Utilisateur 30/31 Mai 2017
Infrastructure technique de
l’Observatoire de la
reconstruction
www.recovery-observatory.org

59

60

61

62

63

64

65

66

67

68

69

70

71

Introduction Imagerie Spatiale

Transféré par

Droits d'auteur :

Formats disponibles

Introduction Imagerie Spatiale

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Introduction Imagerie Spatiale

Transféré par

Droits d'auteur :

Formats disponibles

Imagerie spatiale

• Introduire les particularités de la prise de vue spatiale

• Comprendre le fonctionnement d’un système d’observation spatial

• Exposer la notion de qualité image, afin de savoir :

 Evaluer les performances d’un système d’observation

 Dimensionner un système d’observation qui répond au besoin

Cours Qualité Image : Introduction 2/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 3/90 V5.4 © CNES 2018

350 530 382 361

231 275 341 98

184 311 173 72

738 352 462 215

193 801 508 96

Utilisation de la télédétection spatiale

Cours Qualité Image : Introduction 9/90 V5.4 © CNES 2018

image thermique AVHRR - NOAA Vue 3D Pléiades de la

Synthèse 1998-2007 concentration en phytoplancton

Venise par Pleiades

Cours Qualité Image : Introduction 10/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 12/90

V5.4 © CNES 2018

Lybie Benina airfield BDA, 1986

© CNES 2007/Distribution Spot Image

C’est l’ensemble des TECHNIQUES utilisées pour DETERMINER A

• Pierre angulaire : étude des interactions Matière-Rayonnement

Objet Rayonnement Capteur

Emission Transmission Réception

Cours Qualité Image : Introduction 16/90 V5.4 © CNES 2018

Bandes optiques : désignation courante des bandes

Code abrégé Bande spectrale Longueurs d’onde

Bandes Radar : désignation par des lettres K, X, S, C, L, P

Cours Qualité Image : Introduction 18/90 V5.4 © CNES 2018

Réponse paysages et instrumentales PHR

Cours Qualité Image : Introduction 19/90 V5.4 © CNES 2018

OPTIQUE : PASSIVE RADAR : ACTIF

Cours Qualité Image : Introduction 20/90 V5.4 © CNES 2018

Landsat pass time:

ERS-1 pass time:

ERS (radar C band) Landsat TM

Le radar actif permet d’imager la Terre :

Cours Qualité Image : Introduction 21/90 V5.4 © CNES 2018

B1 ~vert B2 ~rouge B3 ~proche IR

Cours Qualité Image : Introduction 22/90 V5.4 © CNES 2018

Aerosols Water vapor Cirrus

Cours Qualité Image : Introduction 23/90 V5.4 © CNES 2018

Réflectance spectrale de la végétation chlorophyllienne

Visible Proche Moyen Infra-Rouge

Longueur d’onde (mm)

Cours Qualité Image : Introduction 25/90 V5.4 © CNES 2018

• Résolution spatiale: • Résolution spectrale:

• Echantillonnage fin dans le • Echantillonnage fin dans le

• Grande variété des échelles d ’observation :

Cours Qualité Image : Introduction 27/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 28/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 29/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 30/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 31/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 33/90 V5.4 © CNES 2018

Cours Qualité Image : Introduction 34/90 V5.4 © CNES 2018