Преобразование имен файлов в дату в год + недели возвращает Ошибка в charToDate (x): строка символов не в стандартном однозначном формате - PullRequest
Новая
       19

Преобразование имен файлов в дату в год + недели возвращает Ошибка в charToDate (x): строка символов не в стандартном однозначном формате

1 голос
/ 05 апреля 2019

Для анализа временных рядов более 1000 растров в стопке растров мне нужна дата. Данные почти еженедельно в структуре файлов "... 1981036 .... TIF" Ноль разделяет год и неделю Мне нужно что-то вроде: "1981-36"

но всегда получаю ошибку Ошибка в charToDate (x): символьная строка не имеет стандартного однозначного формата 

library(sp)
library(lubridate)
library(raster)
library(Zoo)

raster_path <- ".../AVHRR_All"
all_raster <- list.files(raster_path,full.names = TRUE,pattern = ".tif$")
all_raster

приносит мне: all_raster 

".../VHP.G04.C07.NC.P1981036.SM.SMN.Andes.tif"
".../VHP.G04.C07.NC.P1981037.SM.SMN.Andes.tif"
".../VHP.G04.C07.NC.P1981038.SM.SMN.Andes.tif"
…

Чтобы узнать год и соответствующую неделю, я использовал следующий код: 

timeline <- data.frame(
  year= as.numeric(substr(basename(all_raster), start = 17, stop = 17+3)),
  week= as.numeric(substr(basename(all_raster), 21, 21+2))
)
timeline

приносит мне: хронология 

     year week
1    1981   35
2    1981   36
3    1981   37
4    1981   38
…

Но мне нужно что-то вроде = "1981-35", чтобы я мог построить свой временной ряд позже

Я пробовал это: 

timeline$week <- as.Date(paste0(timeline$year, "%Y")) + week(timeline$week -1, "%U")

и получите ошибку: Ошибка в charToDate (x): символьная строка не в стандартном однозначном формате

или я попробовал это 

fileDates <- as.POSIXct(substr((all_raster),17,23), format="%y0%U")

и получите ту же ошибку

Ответы [ 2 ]

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

Вы можете сделать это, если укажете, что понедельник является рабочим днем ​​1 с %u: 

w <- c(35,36,37,38)
y <- c(1981,1981,1981,1981)
s <- c(1,1,1,1)
df <- data.frame(y,w,s)
df$d <- paste(as.character(df$y), as.character(df$w),as.character(df$s), sep=".")
df$date <- as.Date(df$d, "%Y.%U.%u")

# So here we have variable date as date if you need that for later. 
class(df$date)
#[1] "Date"

# If you want it to look like Y-W, you can do the final formatting:
df$date <- format(df$date, "%Y-%U")

#     y  w s         d    date
# 1 1981 35 1 1981.35.1 1981-35
# 2 1981 36 1 1981.36.1 1981-36
# 3 1981 37 1 1981.37.1 1981-37
# 4 1981 38 1 1981.38.1 1981-38

# NB: though it looks correct, the resulting df$date is actually a character: 
class(df$date)
#[1] "character"

В качестве альтернативы, вы можете сделать то же самое, установив воскресенье на 0 с помощью %w.

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

до тех пор, пока кто-нибудь не опубликует лучший способ сделать это, вы можете попробовать: 

x <- c(".../VHP.G04.C07.NC.P1981036.SM.SMN.Andes.tif", ".../VHP.G04.C07.NC.P1981037.SM.SMN.Andes.tif",
       ".../VHP.G04.C07.NC.P1981038.SM.SMN.Andes.tif")

xx <- substr(x, 21, 27)


library(lubridate)


dates <- strsplit(xx,"0")
dates <- sapply(dates,function(x) {
  year_week <- unlist(x)
  year <- year_week[1]
  week <- year_week[2]
  start_date <- as.Date(paste0(year,'-01-01'))
  date <- start_date+weeks(week)
  #note here: OP asked for beginning of week.  
  #There's some ambiguity here, the above is end-of-week; 
  #uncommment here for beginning of week, just subtracted 6 days.  
  #I think this might yield inconsistent results, especially year-boundaries
  #hence suggestion to use end of week.  See below for possible solution
  #date <- start_date+weeks(week)-days(6)

  return (as.character(date))
})


newdates <- as.POSIXct(dates)
format(newdates, "%Y-%W")

Спасибо @Soren, который разместил этот ответ здесь: Получите месяц недели года

...