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é ?

Annonce

Suite à un problème technique intervenu entre le 22 et le 23 mars, nous avons du procéder dans la soirée du 25 mars, à la restauration de la base de données du 24 mars (matinée).

En clair, nous avons perdu vos contributions et inscriptions du dimanche 24 et du lundi 25 mars.
Nous vous prions de nous excuser.

#1 Tue 12 October 2021 10:38

othx
Participant occasionnel
Date d'inscription: 20 Apr 2020
Messages: 10

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
Participant occasionnel
Date d'inscription: 28 Sep 2009
Messages: 11

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
Participant occasionnel
Date d'inscription: 20 Apr 2020
Messages: 10

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
Participant occasionnel
Date d'inscription: 28 Sep 2009
Messages: 11

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
Participant occasionnel
Date d'inscription: 20 Apr 2020
Messages: 10

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