生成日期向量，其中所有年份都被格式化为闰年R.

Question

如何使用以下信息获取正确的日期：

All years are formatted as leap years with 29 days for February and Julian day 60 are missing in the no-leap years？

我的rasterstack的长度（nlayers）= 23376，我想在1950-01-01到2013-12-31之间添加记录维度。所以我试着像这样生成一个日期向量：

dates=seq(as.POSIXct("1950-01-01 "), as.POSIXct("2013-12-31"), by="day")

但是所有年份都没有如上所述，length(dates) = 23375。

感谢您的建议。 AT

Answer 1

根据@thelatemail和@ 42的建议，我使用以下方式获得了正确的日期：

dates=seq(as.Date("1950-01-01"), as.Date("2013-12-31"), by="day")
> length(dates)
[1] 23376
> table(format(dates,"%Y"))

使用as.POSIXct没有给我完整的日期：

dates=seq(as.POSIXct("1950-01-01 "), as.POSIXct("2013-12-31"), by="day")