在朱利安日期蟒蛇中创建日期范围

Question

我必须在两个日期之间创建一个范围，日期之间的间隔为几分钟，在儒略日期中，我创建了一个代码，但花费了很多时间(15 分钟，例如。

我的代码是：

from astropy.time import Time
import pandas as pd    
timedelta = "600s"
start = "2018-01-01" 
end = "2018-06-30"
dateslist = pd.date_range(start,end, freq =timedelta  ).tolist()
dates =  pd.DataFrame({'col':dateslist})
dates["col2"] =""
for i in range(len(dateslist)):
#print(i," / ", len(dateslist))
dates["col2"][i] = (Time(str(dateslist[i]).replace(" ", "T"), format="fits").jd)

我尝试在没有for的情况下使用时间，但出现错误

time = str(list(dates['col'])).replace("[Timestamp('","").replace(" Timestamp('","").replace("')","").replace(" ","T").split(",")
time
Time(time, format="fits")

值

错误：输入值与格式类匹配

有没有办法快速做到这一点？

谢谢现在，

Answer 1

使用DatetimeIndex.to_julian_date：

dates["col2"] = pd.date_range(start,end, freq = timedelta).to_julian_date()

Answer 2

astropy 中的等效方式是：

from astropy.time import Time
import astropy.units as u
timedelta = 600 * u.s
start = "2018-01-01" 
end = "2018-06-30"
dates["col2"] = np.arange(Time(start).jd, Time(end).jd, timedelta.to_value('day'))

另一种(也许在占星术中更惯用的方式(是：

start = Time("2018-01-01")
end = Time("2018-06-30")
timedelta = 600 * u.s
dates = start + timedelta * np.arange((end - start) / timedelta)

这为您提供了一个矢量Time对象，您可以通过jd属性将其转换为 JD。