Le portail francophone de la géomatique

bonjour
merci pour les documents je l'ai suivi mais toujours j'ai problème
1-    La première étape :
J’ai calibré les trois séparément  suivant basics tools+proccessing +calibrations utilities +land sat clibration
2-    Deuxième étape
J’ai ajusté les unités sur band math+expression b1/10 de chaque bande spectrale
3-    Je l’ai mis dans un méta file (file+ Save as+ENVI meta) et j’ai nommé le fichier contenant  les trois bandes « scène »
4-    Basics Tools +convert data +j’ai sélectionné « scène » puis ok le fichier sortant nommer « scène BIL»
5-    Après spéctral +FLAASH les menus apparaissent. En cliquant sur « input radiance img » j’ai sélectionné  « scène BIL.hdr» +ok (1) le menu dialogue (plus haut est écris –Enter ASCII file name containg waveleghths and Fwnm-) apparaît  j’ai sélectionné « scène BIL» en cliquant sur ouvrir un message d’avertissement  apparait « no valid data   encountred in this file »
Ma question est :
1-    Asque il ya un manque d’un fichier que je dois le produire pour inclure en place de (1) ?
2-    En convertissant les unités asque pour chaque bande ou le méta file ?
 
merci d'avance

Bonjour amele:

La démarche des étapes est dans le bon sens sauf 02 interventions que je dois signalé our pouvoir calibré ta scène. Pour les étapes que je signale:

2- l'ajustement (b1/10) est pour le méta file (radiance.img) que t'as produire avec le suivi du path: Basic tools===Band Math et tu saisi l'expression: b1/10.0 et tu fais ok , ensuite il t'afficheras une fènetre pour indiquerles variable band , tu click sur:Variable Map to input File"et tu choisis ton méta file et tu le saveguarde avec la barre qui s'affiche au dessous.

5- Après convertir ton méta file (scaled : b1/10) en "scène BIL", son ficher header nécessite une information sur la longeur d'onde (wavelengths) de chaque band de la scène; ceci soit tu le trouve indiqué dans un fichier bloc notes (*.txt) ou word, s'il est absent, tu as besoin de céer un fichier ASCII (.txt) qui porte les valeurs par défaut des longueurs d'onde pour une scène de type "Landsat" de chaque band et tu l'ajoute dans fichier header ".hdr" (les wavelengths par défauts sont un intervalle en micromètres comme par exemple (cas d'une scène spot):
(Band1 XS: 0.50 - 0.59 , Band 2 XS: 0.61 - 0.72, Band3 XS: 0.74 - 0.94), alors que les longueurs d'ondes sont une information indispensable pour que l'input " radiance BIL" soit accepté. Bon courage

N.B: si l'étape des wavelengths t'apparait compliqué, je t'orienterai vers une autre méthode pour indiquer les wavelengths/fwhmn.

Bonsoir:

    Puisque la création d'un fichier ASCII te semble compliquée , je t'indique cette méthode simplifiée pour introduire les longueurs d'onde, en répondant sur ta première question oui l'opération se fait pour chaque band séparablement : aprés ouvrir l'image concernée ( Radiance pour ton cas), à partir de la barre menu de l'ENVI tu prend le chemin: File===>Edit ENVI Header , tu choisi l'image, une fenêtre s'affichera , tu click sur "edit attributes , un message s'affiche porte toute une information que tu peux ajouter sur ta scène, tu choisis "wavelengths", une nouvelle fenêtre s'affiche (Edit wavelengths value) avec les trois bands affichée comme suivant :

Band 1: 1.00000
Band 2: 2.00000
Band 3: 3.00000

tu click sur chaque band séparablement, elle s'affichera au dessous 'edit selected item' dans le carré de 1.0000 tu tape la valeur (0.525) , ainsi de suite tu choisis band 2 et tu écris (0.630) et pour la band 3 aussi tu tape la valeur (0.760) dont (0.525, 0.630, 0.760) sont les valeurs moyennes de lambda à chaque band pour une image de type landsat TM (en 03 band RGB) tes bands auront leurs longueurs d'ondes moyenne comme suivant:
-band1: 0.525
-band2: 0.630
-band3: 0.760

il te reste de indiqué l'unité de longueur qui est le micromètre(µm), toujours dans la fenêtre "edit wavelengths value" au dessous tu as : wavelengths/FWHM units: unknown , tu choisis micromètre et tu click ok , et puis ok l'image se ferme mais dans le menu "available band list" tu trouveras les longueurs d'onde affiché et comme sa tu as indiqué la valeurs des wavelengths.

Important: tu indique les wavelengths pour l'image radiance avant de faire la conversion (b1/10.0).

Boncourage.

Dernière modification par ibra_tld10 (Wed 20 October 2010 21:59)

.......

Bonjour,

Le support technique d'ENVI précise que les valeurs de radiance après l'application du "scale factor" de 10 doivent être comprises entre 2 et 30. Vérifies bien, après cette étape, que tes valeurs sont comprises dans cet ordre de grandeur.
Si c'est le cas, il se peut, par exemple, que le modèle atmosphérique choisi ne soit pas adapté à la situation de l'atmosphère lors de l'acquisition de l'image...Il se peut également que le fait d'avoir des pixels négatifs ne soit pas complétement "anormal" (voir les liens postés).

Dans l'autre cas (radiances non comprises entre 2 et 30), il y a surement un problème lors de l'étape de calibration des comptes numériques en radiance (luminance en français). Peut être que les valeurs de gain et de biais appliquées ne sont pas adaptées pour la calibration de ton image.

Voici une liste de liens venant du forum des utilisateurs d'ENVI qui pourraient t'aider à mieux comprendre les subtilités de FLAASH et éventuellement à résoudre ton problème.

http://www.ittvis.com/UserCommunity/Use … fault.aspx
http://www.ittvis.com/UserCommunity/Use … fault.aspx
http://www.ittvis.com/UserCommunity/Use … fault.aspx

Pour vérifier la précision de la calibration en réflectance au sol, je ne connais pas de moyens très fiables quand on ne dispose pas de données de terrain. Personnellement, je compare à l'allure générale des spectres présentés dans les manuels de télédétection ou aux librairies spectrales fournies avec ENVI (auparavant, il faut rééchantillonner les longueurs d'ondes des spectres selon les bandes spectrales de LANDSAT).

Ranfosic

Bonjour,
je en suis pas un expert en ENVi, mais je pense que c'est lié à l'affichage.
les valeurs théoriques du NDVI vont entre - 1 et 1 ... donc visuellement si ton affichage se fait sur une gamme plus large, par ex de  0 à 255, tu ne verra pas grand chose. La plupart des logiciels permettent "d'étaler l'histogramme des valeurs" en modifiant les seuils d'affichage. Dans ton cas, de -1 à 1 par exemple.

Cordialement

Julien

Bonjour,
pour ce qui est de l'ajustement des valeur de rélflectances. il faut que tu vérifie l'unité de radiance (dans le fichier .dim). Si celle si est en "W/m2*sr*µm", il faut dans tous les cas diviser par 10 pour l'apporter en "µW/cm2*sr*nm", unité utilisée pour la correction FLAASH.
Une autre facon de faire c'est ne pas utiliser Band MAth, mais lorsque tu va chercher ton fichier .bil dans "input" de Flaash, tu utilise using single scale factor --> 10
Je rajouterai qu'avec flaash, les résultats des valeurs de réflectances par défauts sont multipliée par 10.000 (Advancing setting -->output reflectance scale factor--> 10000. il suffit juste que tu divise par 10000.
Autre chose: la réflectance est soit comprise entre 0 (pixel noir) et 1 (pixel blanc) soit entre 0 et 100%. 2 et 32 est à titre d'information selon l'échelle que tu veux utilisée.
Pour ton traitement NDVI, il faut également modifier le contraste (linéaire , gaussian, arbitraire.... ) selon si c'est sur scrool ou image ou zoom que tu veux modifier ta gamme de NDVI, et généralement il faut mettre le minimum à 0 et pour le maximum cela dépend de ce que tu veux faire ressortir comme végétation

les valeur de réflectance allant de 780 à 1260 (soit  7,8% à 12% DONC TES VALEURS SONT BIEN COMPRISES ENTRE 2 ET 32)  semblent correct (tout dépend de la densité végétale de la zone);

Non il n'y a pas de valeurs fixe de NDVI car cela dépend de la densité végétal (biomasse), plus celle ci est forte plus le NDVI sera élevé.

Si tu veux envoie moi ton image traité en NDVI, je verrai ou est le probleme.

bonsoir
je vous remercie infiniment pour vos compléments 
concernant les valeurs de réfléctances obtenus après la correction sont les suivants je les vérifier en utilisant   Quick stats :
Stats             Min     Max      Mean             Stdev
     Band 1  -1890   5149   562.461126  684.955043
     Band 2   -674   6017   851.482967  966.754390
     Band 3      0   4862  1095.160520  791.390399
asque c'est bon jusque là ?

Bonjour,

A priori, l'ordre de grandeur pour ton image multispectrale peut être plausible (par exemple 0 à 4862 pour la bande Bleue).
Pour ton NDVI, le min et le max présentent des valeurs aberrantes ( théoriquement le ndvi est compris entre -1 et +1) mais ta valeur moyenne peut être possible.
Il faudrait que tu précises sur quelle région tu travailles. Sur ton image, il y a des grandes étendues d'eau, non?

Le plus simple serait que tu mettes en pièces jointes des images de profils spectraux de différents milieux naturels pour que se fasse une idée de la justesse de ta calibration.

Sinon, pour afficher les classes de valeurs du NDVI: Tools ->  Color Mapping -> Density Slice

Ranfosic

Dernière modification par ranfosic (Wed 08 December 2010 08:41)

Bonjour,

Tu ne dois pas masquer le lac avant de faire ta correction.  Je pensais que ton image comportait de grandes étendues d'eau car tu as des valeurs de réflectance plutôt faibles (comme celles de l'eau).
Je ne comprends pas vraiment pourquoi tu as ces valeurs (peut être un problème si tu prends en compte le "background", c'est à dire le fond noir autour de la scène?). As tu effectué la correction atmosphérique sur toute la scène ou seulement sur une partie de l'image?

Je pense qu'il faudrait que tu postes des images du Z profile de quelques milieux. Par exemple sol nu et forêt de chêne vert.

Ranfosic

bonsoir

pour but d'identifier les types des sols sur ma région d'étude .j'ai voulus utilisé indice de rougeur des sols IR=R^2/V^2

as qu'il y a un manuel adéquat qui nous permettre d'identifier les grandes unités des types des sols suivant la couleur et les valeurs obtenus aprés un traitement par IR?

je vous remercie;amele

bonsoir Ranfosic

mon but d'utilisé cet indice c'est pour donné la nature des sols para port aux types d'occupations exemple :forêt de chêne vert   ou de pin dense sur un sol rouge brun et des maquis (stade de dégradation des formation forestières) les sols devient calcaire et parcours ou végétation herbacées très dégradé sur des sols squelettiques dans cette dernière la couleur devient pale (diminution dioxyde de fer qui donne la couleur rouge des sols ) ......
j'ai essayé un indice de salinité sol-végétation IS et j'ai des résultats raisonnables
pensez vous que l'utilisation IR et de couleur pour complété mes résultats ,est une méthode adéquat pour justifier plus ma classification des unités d'occupations des sols ?

merci ;amele