readMeta fails to read Landsat 8 MSS metadata (*MTL.txt) files

[jkoellin]

Hello bleutner,

thanks for your work with the Toolbox, it helped me in the past.

I have a problem with my new project. I try to read landsat8 data with the readMeta function and I get an error message.

This is my code:

install.packages("RStoolbox")
library(devtools)
install_github("bleutner/RStoolbox")
library(RStoolbox)
metaData_prefire <- readMeta("C:\\Users\\le-ob-dual-1\\Desktop\\Studium\\Ulluoa\\Semesterprojekt\\R_Program\\LC08_L1TP_193024_20180824_20200831_02_T1_richtig\\LC08_L1TP_193024_20180824_20200831_02_T1_MTL.txt")

and this is the error message I get:

Fehler in CRS(paste0(c("+proj=", "+zone=", "+units=m +datum="), pars, collapse = " ")) : 
  No spaces permitted in PROJ4 argument-value pairs: +proj= +zone= +units=m +datum=

my file was downloaded from the earthexplorer and I Included it:
LC08_L1TP_193024_20180824_20200831_02_T1_MTL.txt
Do you have any suggestions, how I could read Landsat8 files like that or maybe have a idea, what causes the error?

Thanks,
Jonas

[RichardLemoine]

Dear @bleutner,

adittionally to this, I'm experiencing another error when trying to read the metadata of a Landsat 8 collection 2 image.
My code:

library(RStoolbox)
meta <- readMeta("E:/Landsat/LC08_L1TP_025046_20190412_20190422_02_T1_MTL.txt")

This is the error:

Error in as.POSIXct.default(if (!legacy) meta$PRODUCT_METADATA["DATE_ACQUIRED", : c("do not know how to convert 'if (!legacy) meta$PRODUCT_METADATA["DATE_ACQUIRED", ] else meta$PRODUCT_METADATA["ACQUISITION_DATE", ' to class “POSIXct”", "do not know how to convert ' ]' to class “POSIXct”")

I think that the error should be due to a change in the structure of the metada in collection 2 in comparison to collection 1. I tested with the same scene from collection 1 and it works as expected. I hope this is an eassy to solve issue.

Regards,
Richard

[jgomezd]

Dear Bleutner,

I am also having the problem that the previous people report when reading a Landsat 8 Image Level 2 of Collection 2:

library(RStoolbox)
M <- readMeta("imagen_satelite/LC08_L2SP_025046_20210127_20210305_02_T1_MTL.txt")

I have the following message:
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=

R version 4.0.4 (2021-02-15) -- "Lost Library Book"
RStoolbox version 0.2.6
Windows 10, x64

[RikFerreira]

Dear @bleutner,

I'm experiencing exactly the same issue related previously with the MTL relative to Landsat 5 TM image LT05_L1TP_218072_20100801_20161015_01_T1 downloaded from USGS.

For the record, the MTL file from Landsat 5 TM image LT05_L1TP_218072_20100716_20161014_01_T1 works just as expected.

• Working file
• Bad file
• Diff

Does anyone who went through this problem found any solution?

Thank you all for your time,
Rik

[FelipeTC91]

Dear @bleutner

I have the same problem with the MTL file for landsat 8

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

I really appreciate your help with this error

Thanks

[jjarosch]

The projection error arises here:

RStoolbox/R/readMeta.R

Line 69 in 7b85a19

pars <- meta$PROJECTION_PARAMETERS[c("MAP_PROJECTION","UTM_ZONE","DATUM"),]

The group/slot is now called LEVEL1_PROJECTION_PARAMETERS. A number of groups have LEVEL1_ prefixed in Collection 2.

[bleutner]

Thanks everyone for your reports. Fixed now.