Suite

Quelle méthode DOS utilise pour convertir la radiance au capteur en valeurs à la surface dans GRASS (données landsat, module i.landsat.toar) ?

Quelle méthode DOS utilise pour convertir la radiance au capteur en valeurs à la surface dans GRASS (données landsat, module i.landsat.toar) ?


j'essaye de calculer réflectance de surface valeurs des données Landsat DN brutes téléchargées depuis GLOVIS. J'ai importé des données Landsat dans GRASS 6 avec en-tête par i.landsat.import et reprojeté la scène dans différents PCS. Ensuite, je peux convertir DN en valeurs de radiance/réflectance au capteur en utilisant i.landsat.toar.

Par la suite, je voudrais convertir valeurs au capteur aux valeurs à la surface (nécessaire pour plus i.topo.corr correction de minnaert en raison de la zone montagneuse). La méthode simple est implémentée par les méthodes DOS[1-4] dans i.landsat.toar. Cependant, il existe 4 méthodes DOS (http://grass.osgeo.org/grass65/manuals/i.landsat.toar.html) :

  • DOS1 : TAUv = 1,0, TAUz = 1,0 et Esky = 0,0
  • DOS2 : TAUv = 1,0, Esky = 0,0 et TAUz = sin(e) pour toutes les bandes avec une longueur d'onde maximale inférieure à 1. (c'est-à-dire les bandes 4-6 MSS, 1-4 TM et 1-4 ETM+) autres bandes TAUz = 1,0
  • DOS3 : TAUv = exp[-t/cos(sat_zenith)], TAUz = exp[-t/sin(e)], Esky = rayleigh
  • DOS4 : TAUv = exp[-t/cos(sat_zenith)], TAUz = exp[-t/sin(e)], Esky = PI * radiance_dark

Avez-vous des expériences, laquelle des méthodes DOS est la plus simple et fournit les résultats les plus fiables en zone montagneuse ? Mon domaine d'intérêt varie de 800-2500 a.s.l. en Europe centrale, surface d'environ 1050 km2 (seulement une fraction de l'ensemble de la scène Landsat).


J'ai trouvé Dos3 plus réaliste pour Landsat 5 dans les plaines pour les cultures vertes. Aucune expérience dans d'autres situations.

Landsat 8 réagira différemment aux différentes méthodes de correction, par rapport à L5.

Je n'ai pas utilisé GRASS depuis un moment, j'ai migré vers le plugin pour QGIS - voir http://fromgistors.blogspot.com/