MTL file error

[abenson-ut]

Hello,

I am trying to do some atmospheric corrections on Landsat 8 data in Rstudio. Following previous code that I have seen I have run this so far.

install.packages("raster")
install.packages("rgdal")
install.packages("ggplot2")
install.packages("devtools")
install.packages("sf")

Define the multiple mutiple packages in a vector

PkgToLoad <- c("raster", "rgdal", "ggplot2", "devtools", "sf")

Load the packages

lapply(PkgToLoad, library, character.only = TRUE)

library(devtools)

The package 'devtools' is important in order to be able to import packages directly from the Github repository

Had to update all packages before this code would work

install_github("bleutner/RStoolbox")

library("RStoolbox")

#perform the radiometric correction
#Import meta-data and bands based on MTL file

readMeta("C:/Users/Allyson/University of Texas at El Paso/Microsoft forest SD project - General/ABP/PR_Trial/PuertoRicoL8USGS/LC08_L2SP_005047_20180616_20200831_02_T1_MTL.txt") #sep = '=', stringsAsFactors = F package = "RStoolbox")

-when I try to read the MTL file I get this error

Error in .rowNamesDF<-(x, value = value) :
duplicate 'row.names' are not allowed
In addition: Warning message:
non-unique value when setting 'row.names': ‘QA_NA’

I attached the MTL text file. Can anyone give me input on how I can fix this?

LC08_L2SP_005047_20180616_20200831_02_T1_MTL.txt

[Rapsodia86]

Try raw=TRUE, so it will read all info.
x <- readMeta("LC08_L2SP_005047_20180616_20200831_02_T1_MTL.txt",raw = TRUE)

looks like readMeta() has no LEVEL2 product included when you go default with raw=FALSE:
lines 81-84:

prodid  <- .getParFromAnyGroup("LANDSAT_PRODUCT_ID", "LEVEL1") 
level   <- .getParFromAnyGroup("PROCESSING_LEVEL", "LEVEL1") 
colNum  <- .getParFromAnyGroup("COLLECTION_NUMBER", "LEVEL1") 
colTier <- .getParFromAnyGroup("COLLECTION_CATEGORY", "LEVEL1") 

Then problem in line 100
bands <- regmatches(files, regexpr("(?<=[\\dQA]_)[^_]*(?=\\.TIF)|B6_VCID.*(?=\\.TIF)", files, perl = TRUE))
finds only:
[1] "AEROSOL" "PIXEL" "RADSAT"
while changing to:
bands <- regmatches(files, regexpr("(?<=[\\dQA_])[^_]*(?=\\.TIF)|B6_VCID.*(?=\\.TIF)", files, perl = TRUE))
gives:

[1] "B1"      "B2"      "B3"      "B4"      "B5"      "B6"      "B7"      "B10"     "TRAD"    "URAD"    "DRAD"    "ATRAN"   "EMIS"    "EMSD"    "CDIST"   "AEROSOL" "QA"      "PIXEL"   "RADSAT" 

Next thing, In level 2, there is no ANGLE tifs, so the MTL's PRODUCT_CONTENTS is different than those in the test folder.
So, line 104 returns all FALSE obviously.

@bleutner including LEVEL 2 would be great!:)

Btw. Landsat L2 SP is already atmospherically corrected (https://www.usgs.gov/landsat-missions/landsat-collection-2-surface-reflectance)

[andrewmaclachlan]

I also had a similar issue when using Landsat 9 data, as it is only distributed as a level 2 product