SMAC prend en charge de nouveaux satellites

=>De nouveaux coefficients ont été ajoutés au site du CESBIO . Les nouveaux satellites pris en compte sont :

  • Landsat 8, RapidEye
  • Quickbird, Worldview2, Ikonos
  • Pléiades1A (PHR1A)
  • DMC-DEIMOS1

 Le Simplifié Modèle d’Atmosphérique Correction (SMAC) est parfaitement adapté à l’implémentation rapide et approchée de corrections atmosphériques. Il s’agit de fonctions analytiques dérivées du modèle 5S. Les 49 coefficients de ce modèle sont ajustés à partir de simulations de transfert radiatif obtenues avec le modèle 6S (l’ancienne version, pas la récente version vectorielle). SMAC n’est pas un modèle très précis (beaucoup moins que MACCS), et il faut lui fournir des données auxiliaires pour l’épaisseur optique des aérosols ou pour les contenus atmosphériques en ozone et vapeur d’eau. Quand ces données sont précisément connues, la précision des simulations est en général meilleure que deux à trois pour cent, sauf parfois pour les grands angles (au dessus de 70°) ou dans de fortes bandes d’absorption et si on ne prend pas en compte les effets d’environnement et les effets de pente. SMAC a été conçu pour être facile à utiliser :

#lecture des 49 cofficientsnom_smac ='COEFS/coef_FORMOSAT2_B1_CONT.dat'coefs=coeff(nom_smac)#Lire la réflectance TOA de la bande à traiter#(celà va dépendre du format de l'image)#Lire les valeurs des angles dans les métadonnées (ici, on les fixe)theta_s=30phi_s=180theta_v=0phi_v=0# calculer la pression atmosphérique à l'altitude du pixelpressure=PdeZ(1300)#trouver les valeurs des variables atmsophériques (ici, on les fixe arbitrairement)AOT550=0.1UO3=0.3UH2O=3#calculer la correction atmosphériquer_surf=smac_inv(r_toa,theta_s,phi_s,theta_v,phi_v,pressure,AOT,UO3,UH2O,coefs)

Dans la dernière ligne ci-dessus :

  • theta_s, phi_s sont resp. l’angle solaire zenithal et azimuthal
  • theta_v, phi_v sont resp. l’angle de visée zenithal et azimutha
  • AOT est l’épaisseur optique à 550 nm qui peut provenir d’une station Aeronet, ou fixée au jugé, ou égale à 0.1 (si on se contente d’une correction très approchée).
  • UO3 est le contenu en ozone, en cm.atm (0.3 est souvent suffisamment précis)
  • UH2O est le contenu en vapeur d’eau, en kg/m². J’utilise souvent une valeur égale à 3, quand je ne cherche pas une grande précision, dans des bandes ou l’absorption est faible
References

[1] Rahman, H., & Dedieu, G. (1994). SMAC: a simplified method for the atmospheric correction of satellite measurements in the solar spectrum. REMOTE SENSING, 15(1), 123-143.« [2] »Tanré, D., Deroo, C., Duhaut, P., Herman, M., Morcrette, J. J., Perbos, J., & Deschamps, P. Y. (1990). Technical note Description of a computer code to simulate the satellite signal in the solar spectrum: the 5S code. International Journal of Remote Sensing, 11(4), 659-668.« [3] »Vermote, E. F., Tanré, D., Deuze, J. L., Herman, M., & Morcette, J. J. (1997). Second simulation of the satellite signal in the solar spectrum, 6S: An overview. Geoscience and Remote Sensing, IEEE Transactions on, 35(3), 675-686.>« [4] »Kotchenova, S. Y., Vermote, E. F., Matarrese, R., & Klemm Jr, F. J. (2006). Validation of a vector version of the 6S radiative transfer code for atmospheric correction of satellite data. Part I: Path radiance. Applied Optics, 45(26), 6762-6774.« [5] »Kotchenova, S. Y., & Vermote, E. F. (2007). Validation of a vector version of the 6S radiative transfer code for atmospheric correction of satellite data. Part II. Homogeneous Lambertian and anisotropic surfaces. Applied Optics, 46(20), 4455-4464. 

Plus d'actualités

Dates de fauche en France

En 2024, des collègues du Cesbio ont publié un article* sur la cartographie de la date de fauche en France en 2022 à partir des données satellitaires Sentinel-2.   Leur magnifique figure 10 a attiré mon attention. La version à haute résolution fournie avec l’article (2861 × 2911 pixels) correspond à une image dont les pixels font […]

Évolution du jour de déneigement dans les Alpes françaises et les Pyrénées

Les socio-écosystèmes des Alpes et des Pyrénées dépendent étroitement des fluctuations annuelles du manteau neigeux. En particulier, le moment de l’année où la neige disparait détermine le début de la saison de croissance de la végétation de montagne et donc la période des estives. Le changement climatique est en train de bouleverser ce rythme saisonnier. […]

Evolution de l’altitude de la ligne de neige au cours des 41 dernières années dans le bassin versant du Vénéon (Oisans)

Pour contribuer à caractériser les conditions hydrométéorologiques lors de la crue torrentielle qui a frappé la Bérarde en juin, j’ai analysé une nouvelle série de cartes d’enneigement qui couvre la période 1984-2024 [1]. Grâce à la profondeur temporelle de cette série, on constate que l’altitude de la ligne de neige dans le bassin versant du […]

Rechercher