转换日期（YYYY MM DD HH：MM：SS）到十进制日（YYYY MM DD.ddddd）

Question

我制作了一个脚本，它在望远镜上观察观察的日期和时间，并将其转换为带小数日的日期。此脚本还占用了一半的曝光时间并将其添加到十进制日，以便我可以获得最准确的天体测量计算。然后将其提交给小星球中心，该中心仅​​接受年月日的观察时间.XXXXX这样做的最小准确方法是5小数点是什么？这是我目前使用的方式，是的它非常凌乱，但它确实让我得到十进制日。偶尔它会在时间转换中关闭一秒钟。

expmpc = float.("10")
utcstartmpc = "05:45:19.03"
datempc = "2015-02-14"


ddexposure = expmpc/(60.0*60.0*24.0)
ddexposure = round(ddexposure, 5)
seconds = int(utcstartmpc[6:8]) / 60.0
minutes = (int(utcstartmpc[3:5]) + seconds) / 60
hours = (int(utcstartmpc[:2]) + minutes) / 24
hours = round(hours, 5)
expadd = str(hours + ddexposure)
datempc = datempc.replace("-", " ")
utcmpc = "C%s%s" % (datempc, expadd[1:])
utcmpc = utcmpc.ljust(17, "0")

正如你所看到的那样，这非常混乱，它涉及到大量数据，我相信我在舍入时会失去准确性。代码的完成结果留下了这样的时间： C2015 02 14.23986

是否有更好的模块？

感谢您的帮助。

Answer 1

这里有一些使用更多内置时间/日期模块的东西（我认为它可以做你想要的）：

import time
import datetime

def my_date_converter(date_in):
    parse_str, _ = date_in.split(".")
    date_str, time_str = parse_str.split(" ")
    parsed_time = time.strptime(time_str, "%H:%M:%S")
    total_seconds = datetime.timedelta(hours=parsed_time.tm_hour,minutes=parsed_time.tm_min,seconds=parsed_time.tm_sec).total_seconds()
    seconds_in_day = 86400.0
    decimal_seconds = total_seconds / seconds_in_day
    utcmpc = "C%s%s" % (date_str.replace("-", " "), str(round(decimal_seconds, 5))[1:])
    utcmpc = utcmpc.ljust(17, "0")
    return utcmpc


def main():
    to_convert = "2015-02-14 05:45:19.03"
    converted = my_date_converter(to_convert)
    print "%s => %s" % (to_convert, converted)

if __name__ == '__main__':
    main()

示例输出：2015-02-14 05:45:19.03 => C2015 02 14.23980

Answer 2

您是否尝试过优秀的星座库？他们有一个处理时间和转换的软件包：astropy.time

Answer 3

您可以将问题分成两个任务：

• 从输入datetime，expmpc，utcstartmpc
获取表示UTC时间的datempc个对象
• 将datetime对象转换为Year Month Day.XXXXX格式的字符串。

从输入datetime，expmpc，utcstartmpc

获取datempc个对象

from datetime import datetime, timedelta

expmpc = "10"
utcstartmpc = "05:45:19.03"
datempc = "2015-02-14"
dt = datetime.strptime(datempc + " " + utcstartmpc, "%Y-%m-%d %H:%M:%S.%f")
dt += timedelta(seconds=int(expmpc))

将datetime对象转换为Year Month Day.XXXXX格式

的字符串

from datetime import datetime, time, timedelta

s = dt.strftime("C%Y %m %d.")
day_fraction =  dt - datetime.combine(dt, time.min)
s += ("%.0f" % (100000*day_fraction / timedelta(days=1)))
# -> C2015 02 14.23992

注意：结果与您的问题略有不同（最后为92与86）。

在Python 2中，您需要将td / timedelta(1)替换为td.total_seconds() / 86400。

这是另一种方式：

>>> from __future__ import division
>>> assert day_fraction.days == 0
>>> '%.0f' % ((day_fraction.seconds*10**6 + day_fraction.microseconds) / 864000)
'23992'

真正除法之前的所有算术运算都以无限精度执行。

另一种方式也产生92：

>>> from datetime import timezone
>>> ('%.5f' % ((dt.replace(tzinfo=timezone.utc).timestamp() % 86400) / 86400))[2:]
'23992'