Временной ряд данных MODIS в пикселях с датой на оси x

Question

У меня есть 65 изображений MOD11A2. Я установил дату, используя as.Date функция, а затем я создал стек временных рядов, используя timeStackMODIS функция bfastSpatial пакет. Я хотел создать график временных рядов пикселей после этого поста. Проблема в том, что мне не удалось установить дату по оси х. Какие-нибудь мысли? Вот пример моего кода.

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)

Полученный график

Чтобы установить дату, я использовал этот код:

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])
}