用x轴上的日期绘制MODIS数据的像素时间序列

Question

我有65个MOD11A2图像。我使用as.Date函数设置日期，然后使用timeStackMODIS包的bfastSpatial函数创建了一个时间序列堆栈。我想在此post之后创建一个像素的时间序列图。问题是我无法在x轴上设置日期。有什么想法吗？这是我的代码示例。

library(bfastSpatial)
library(raster)
library(rgdal)

setwd("my_dir")

#create list
list = list.files(path = "my_dir", full.names = FALSE, pattern = "\\.tif$" )

#time series layer stack of the list using the bfastspatial
m8 = timeStackMODIS(list, pattern = glob2rx("\\.tif$"))

plot_ts <- function(x) {
  plot(x[[1]])
  z <- c(click(x, n=1, show=FALSE))
  plot(z, type='l', lwd=2, ylab='y', xlab='time', las=1)
  z
}

z <- plot_ts(m8)

结果图是这个one

设置我使用此代码的日期：

li<-as.data.frame(list.files(pattern = ".tif|.TIF"))
li$nn<-paste0(substr(li[,1],1,9),format(as.Date(substr(li[,1],10,16), "%Y%j"),"%Y.%m.%d"),substr(li[,1],17,75))
for(i in 1:nrow(li)){
  file.rename(as.character(li[i,1]),li[i,2])
}