Annonce
L'appel à participation est ouvert jusqu'au 19 janvier 2025!
#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: 15
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: 15
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
