Nous utilisons des cookies pour vous garantir la meilleure expérience sur notre site. Si vous continuez à utiliser ce dernier, nous considèrerons que vous acceptez l'utilisation des cookies. J'ai compris ! ou En savoir plus !.
banniere

Le portail francophone de la géomatique


Toujours pas inscrit ? Mot de passe oublié ?
Nom d'utilisateur    Mot de passe              Toujours pas inscrit ?   Mot de passe oublié ?

#1 Tue 12 October 2021 10:38

othx
Membre
Date d'inscription: 20 Apr 2020
Messages: 7

Corrections radiométriques avec R

Bonjour,

Je tente de passer des fichiers de bandes spectrales en comptes numériques (Landsat 8) en fichiers de réflectance de surface via R.

Voici le code que je mets en place :

metaData  <- readMeta("LC08_L1TP_183033_20210623_20210630_02_T1_MTL.txt")
lsat <- stackMeta("LC08_L1TP_183033_20210623_20210630_02_T1_MTL.txt")
lsat_sref <- radCor(lsat, metaData, method = "dos")

J'obtiens cette erreur :
Error in CRS(paste0(c("+proj=", "+zone=", "+units=m +datum="), pars, collapse = " ")) :
  No spaces permitted in PROJ4 argument-value pairs: +proj= +zone= +units=m +datum=

Pourriez-vous m'aider svp ?
Merci
Olivier

Hors ligne

 

#2 Tue 12 October 2021 20:15

Marc Gauthier
Membre
Date d'inscription: 28 Sep 2009
Messages: 9

Re: Corrections radiométriques avec R

Bonsoir Olivier,
Un petit traceback() dans RStudio/Rgui permettra d'identifier la fonction générant l'erreur.
Marc

Hors ligne

 

#3 Tue 12 October 2021 20:35

othx
Membre
Date d'inscription: 20 Apr 2020
Messages: 7

Re: Corrections radiométriques avec R

> traceback()
2: stop("metaData must be a path to the MTL file or an ImageMetaData object (see readMeta)")
1: radCor(lsat, metaData, method = "dos")

C'est à dire ? En faisant cela par exemple ?

Hors ligne

 

#4 Wed 13 October 2021 19:30

Marc Gauthier
Membre
Date d'inscription: 28 Sep 2009
Messages: 9

Re: Corrections radiométriques avec R

un summary(metaData) répond quoi ?

J'ai joué l'exemple du package sans erreur :

Code:

library(RStoolbox)
library(rgdal)
## Example metadata file (MTL)
mtlFile  <- system.file("external/landsat/LT52240631988227CUB02_MTL.txt", package="RStoolbox")

## Read metadata
metaData <- readMeta(mtlFile)
summary(metaData)

## Load rasters based on metadata file
lsat <- stackMeta(mtlFile)

lsat_sref <- radCor(lsat, metaData, method = "dos")

Dernière modification par Marc Gauthier (Wed 13 October 2021 19:32)

Hors ligne

 

#5 Wed 13 October 2021 21:50

othx
Membre
Date d'inscription: 20 Apr 2020
Messages: 7

Re: Corrections radiométriques avec R

J'ai trouvé ma réponse ici https://stackoverflow.com/questions/695 … ons-with-r
Merci pour vos réponses.

Hors ligne

 

Pied de page des forums

Powered by FluxBB

Partagez  |