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