从R中的LANDSAT8图像名(扩展)中提取数据

Question

我想从LANDSAT8文件夹名中提取日期(YYDDD格式)，并在为多幅图像(52)计算的植被索引中使用该数据。

以下是我到目前为止所拥有的：

setwd('E:/Landsat8')
folders <- list.dirs(full.names=FALSE, recursive = FALSE)

for(scene in folders){
  NIR <- raster(paste(scene,'/',scene,'_B5.TIF', sep=''))
  SWIR1 <- raster(paste(scene,'/',scene,'_B6.TIF', sep=''))
  LSWI <- overlay(x=NIR, y=SWIR1, fun=calcIndex)     #calculates LSWI 
  writeRaster(LSWI, filename=paste(scene,'/',scene,'_LSWI.TIF', sep=''), format='GTiff', datatype='FLT4S', overwrite=TRUE) 
  }

LANDSAT8图像中文件名的格式如下: LC81730382014069LGN00

我希望将文件安全为2014_069_LSWI.TIF。目前它正在将其保存为LC81730382014069LGN00_LSWI.TIF，我无法找到解决方案，只能从名称字符串的中间提取日期。希望有人有个主意！

Answer 1

如果没有弄错，Landsat文件名中的年份子串和日子串总是位于相同的位置。因此，这应该做到：

#Sample image name
nm <- "LC81730382014069LGN00"

#Extract year
yr <- substr(nm, start=10, stop=13)

#Extract day of year
dy <- substr(nm, start=14, stop=16)

#Make name
nm2 <- paste(yr, dy, "LSWI.TIF", sep="_")

然后，您可以将nm2作为文件名在writeRaster期间传递，如下所示。注意paste0的用法;-)

writeRaster(LSWI, filename=paste0(scene,'/', nm2), format='GTiff', datatype='FLT4S', overwrite=T)