Corrections radiométriques des images optiques satellitales (Travail pratique 2)

1. Corrections radiométriques des images optiques satellitales (Travail pratique 2) Pour le moment dans notre cours: • Rayonnement solaire réfléchi  trouver la réflectance au sol des mesures des capteurs • Rayonnement émis  trouver la température au sol des mesures des capteurs

2. Rayonnement solaire réfléchi (RSR) • Étape 1: Conversion des valeurs numériques en luminances au capteur • Étape 2: Calcul des luminances au sol par élimination des effets atmosphériques • Étape 3: Conversion des luminances au sol en réflectances au sol

3. Rayonnement émis (IRT) • Étape 1: Conversion des valeurs numériques en luminances au capteur • Étape 2: Calcul des luminances au sol par élimination des effets atmosphériques • Étape 3: Conversion des luminances apparentes au sol en luminances vraies au sol en éliminant les effets de l’émissivité • Étape 4: Calcul des températures au sol par les luminances vraies

4. Étape 1 (RSR et IRT) • Étalonnage du capteur  Metadonnées Exemple: Landsat-7

5. Étape 2 (RSR et IRT) •  transmittance atmosphérique surface-capteur •  luminance atmosphérique «parasite» Question: comment connaître les paramètres atmosphériques au moment du passage du satellite?

6. Étape 3 (RSR) c  conversion W/m2 en W/m2/sr Question: comment connaître l’éclairement au sol au moment du passage du satellite?

7. Étape 3 (IRT) c  conversion W/m2 en W/m2/sr Question: comment connaître l’éclairement au sol au moment du passage du satellite; l’émissivité et la réflectance de surface dans l’IRT?

8. Étape 4 (IRT) • Calcul de la température par le Lsol (vraie)  approximation de la loi de Plank

9. Réponses aux questions • Plusieurs façons d’aborder le problème  ici modélisation physique du transfert radiatif dans le système surface-atmosphère (code atmosphérique)

10. RSR (Transfert radiatif)

11. Qu’est-ce qu’un code peut faire? Un code simule la quantité du rayonnement solaire (éclairement, luminance) mesurable à un point quelconque dans le système Terre-atmosphère et à une longueur d’onde quelconque. Pour ce faire : • Solution approximative de l’équation intégro- différentiale du transfert du rayonnement solaire dans le système Terre-atmosphère, pour cela il faut aussi: • Fixer les conditions aux limites du système • Fixer la structure de l’atmosphère • Décrire les deux phénomènes de base: absorption + diffusion

12. Exemple: le CODE 6S (Second Simulation of a Satellite Signal in the Solar Spectrum) • Solution: méthode de qualité, relativement rapide à exécuter • Conditions aux limites du système: plusieurs possibilités pour la surface terrestre, possibilités de simuler diverses conditions géométriques d’illumination et d’observation • Structure de l’atmosphère: atmosphère multi-couches, introduction de l’altitude du terrain, plusieurs possibilités quant aux composants (gaz, aérosol) • Phénomènes: absorption, diffusion, possibilité de mieux approcher les effets de la vapeur d’eau.

13. Structure de l’atmosphère: chaque couche est peuplée d’un nombre de molécules des gaz et de particules en suspension (aérosol)

14. Les conditions aux limites: L’éclairement solaire aux confins de l’atmosphère (supposé connu  observations)

15. Conditions géométriques: illumination-observation (cas satellite)

16. Passage par l’atmosphère

17. Diffusion et absorption La constitution de l’atmosphère en molécules des gaz pas mal connue et stable à l’exception de la vapeur d’eau et dans une moindre mesure de l’ozone stratosphérique; Les propriétés de diffusion et d’absorption des molécules des gaz pas mal connues; Le grand inconnu les particules de l’aérosol;

18. ProfilAtmophérique : zone tempérée : été Température [K] Log(O3 [gm-3]) P H2O T O3 Altitude [km] Log(H2O [gm-3]) 2-3 Log(P [mb]) -32 Pour les calculs de l’absorption et de la diffusion par les molécules des gaz nécessaire de connaître profil vertical température/pression ainsi que le profil vertical ozone + vapeur d’eau: modèles par défaut

19. Pour les aérosols on suppose la composition connue (rural, urbain, maritime,…) et on cherche à estimer la charge par les données satellitales mêmes ou d’autres approches (plus loin). Pour la distribution verticale autres hypothèses…. Également le plus souvent on suppose l’indépendance des phénomènes d’absorption et de diffusion

20. Le code alors calcule l’éclairement direct incident à un élément au sol Éclairement direct du soleil

21. Le code alors calcule l’éclairement du ciel incident à un élément au sol Éclairement diffus du ciel

22. Comment réagit la surface face à l’éclairement solaire direct? La réflectance bidirectionnelle

23. Comment réagit la surface face à l’éclairement diffus du ciel? La réflectance hémisphérique-directionnelle Pour toute direction de provenance on trouve la réflectance bidirectionnelle et l’on fait la somme

24. Attention!!! LES COURBES DE SIGNATURES SPECTRALES  RÉFLECTANCE HÉMISPHÉRIQUE-HÉMISPHÉRIQUE

25. Réflectance bidirectionnelle • Le code 6S propose une série de modèles adaptées à différents types de surface (surtout couvert végétal) • Dans la pratique difficultés d’appliquer ces modèles car on ne connaît pas le type de surface

27. Le plus souvent (malgré la validité questionnable de l’hypothèse) on suppose que la surface est isotrope alors: La réflectance bidirectionnelle = réflectance hémisphérique-directionnelle; alors Luminance quittant la surface = Lsol = Etot*ρ/ 

28. Le code calcule la transmittance Luminance de la surface est atténuée par l’atmosphère

29. Le code calcule la luminance parasite À cette luminance s’ajoute la luminance propre à l’atmosphère Lp

30. Rayonnement dans le visible quittant l’atmosphère selon le type de l’aérosol Maritime Rural Urbain

31. Réponse spectrale [%] Réponse spectrale [%] Longueur d’onde [mm] Longueur d’onde [mm] La luminance par longueur d’onde qui arrive au capteur est finalement: Lsat = Lsol Tr + Lp Le capteur filtre le rayonnement reçu selon la bande spectrale: la luminance totale est la somme de toute luminance spectrale qui passe par le filtre pondérée par la sensibilité du détecteur de cette bande SENSIBILITÉ SPECTRALE DES BANDES TM3,TM4 et TM5 TM 3 TM 4

32. 1. Capteurs spéciaux2. Cibles obscures3. Visibilité Estimation de la charge totale d’aérosols  Profondeur optique d’aérosols

33. 1. Capteurs spéciaux : MISR Résolution basse par de système disponible équivalent pour les résolutions moyennes et hautes

34. 2. La méthode des cibles obscures Lsat = Lsol Tr + Lp Lsat = Lp Application du code d’une manière itérative  trouve la profondeur optique qui donne une luminance parasite équivalente à celle observée

35. 3. La visibilité Si l’on dispose d’une mesure de la visibilité (exemple aéroport) 6S a des modèles de calcul de la profondeur optique (méthode très approximative)

36. Un exemple de correction

37. Les choses se complexifient lorsque le terrain est accidenté Le modèle de réflectance isotrope simplifie les calculs mais des erreurs importantes.

38. Les capteurs aéroportés: effets importants de réflexion bidirectionnelle

39. IRT (Transfert radiatif)

40. Formulation • Approximation : Comment faire pour enlever l’effet de l’émissivité?

41. Solutions? • 1 seul canal (ex. Landsat) • 2 et plus canaux (ex. ASTER)

42. 1 canal Créer des cartes approximatives d’émissivité Ex. type d’occupation du sol  valeurs par défaut

43. Plusieurs canaux Ex. Algorithme TES (temperature-emissivityseparation) Idée de base: peu importe le canal la température de surface reste la même; Alors processus itératif: quelles sont les valeurs de l’émissivité pour que les valeurs de Lsol (vraie) calculées pour une température, qui reste constante pour l’ensemble de bandes, sont équivalentes à Lsol (vraie) observées simultanément pour tous les canaux