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

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

0 голосов
/ 25 апреля 2019

У меня 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)

Полученный график такой: 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])
}
...