将时间向量转换为儒略日答案

【问题标题】：Converting time vectors into Julian days将时间向量转换为儒略日
【发布时间】：2016-04-19 13:46:27
【问题描述】：

我有一个 .txt 文件，其中第二列表示 367 天（测量结束日期为 16 年 4 月 1 日）期间的测量时间（以秒为单位）（测量开始日期为 15 年 4 月 1 日）。第三个等是测量变量。

在我的分析中，我需要在儒略日（在处理这种类型的日期时，我通常将其用作“原点”01/01/1992）来评估其他文件中已经存在的其他变量儒略时间的日期。如果有的话，我也可以使用日历日期！

文件示例：

（第一天）

091.00000   0.00    102.63  16.79   81.40   5.07    144.20  0.00
091.00069   60.00   102.63  16.80   81.30   5.86    136.00  0.00
091.00139   120.00  102.63  16.75   81.20   2.46    144.20  0.00

(...)

091.99931   86340.00    102.53  17.09   82.60   1.60    117.00  0.00

（第二天）

092.00000   0.00    102.53  16.96   82.60   1.56    115.00  0.00
092.00069   60.00   102.53  16.87   82.90   1.38    114.70  0.00
092.00139   120.00  102.53  17.01   82.90   1.13    110.30  0.00

(...)

092.99931   86340.00    102.44  17.38   87.10   1.75    116.00  0.00

在整个 367 天的所有日子里都会发生这种情况（对应于总共 526970 次测量），我将所有这些都放在一个文件中。可能会丢失数据记录，因此我可以在这些间隙中使用缺失值 (NA)。

我正在使用 R。

【问题讨论】：

你的问题不是很清楚。一些示例数据会更清楚。
我编辑了这个问题，希望现在更清楚了。谢谢。

标签： r time-series julian-date

【解决方案1】：

第二列似乎是自当天开始以来的秒数。 24 小时/天 * 3600 秒每小时等于每天 86400 秒。所以SecondColumn/86400 给出儒略日期的小数。第二列不包含有关日期本身的信息。

【讨论】：

好的，所以我拥有的第一列已经是儒略日期（0.9100000、0.9100069、...）。那么当我将此列作为 x 轴绘制数据时，我一定做错了。我基本上将 tm2_1min=dat2$V1 （第一列）用于 x 轴，将 atmos_press_1min=dat2$V3 （第三列）用于一个，例如变量然后 plot(tm2_1min,atmos_press_1min,type="l",xaxt = "n") 轴（side=1,at=datJ,lab=labJ）定义为 datJ=seq(from=julian(04,01, 2015, origin=c(1,1,1992)), to=julian(02,04,2016, origin=c(1,1,1992))) labJ=seq.dates("04/01/2015", "02/04/2016",by="days")
嘿，我已经设法解决了我的问题。我在文件中的 Julia 日子从 2015 年开始计算，而不是像我习惯的那样从 1992 年开始计算。此外，当计数达到 365 时，它又回到 091.000，当我想绘制数据时，R 对此造成了麻烦。所以我只是更正了计数，使其从 365 一直上升到 456。我现在可以很好地绘制数据。感谢您的提示。