从tzdata中提取历史性的闰秒
有没有一种方法可以从分布在大多数linux发行版的时区数据库中提取历史性闰秒的时刻？我正在寻找python中的解决方案，但是在命令行上工作的任何东西都可以。
我的用例是在GPS时间(基本上是自1980年第一颗GPS卫星启用以来的秒数)和UTC或当地时间之间进行转换。UTC不时对闰秒进行调整，而gps时间则呈线性增长.这相当于UTC和泰语之间的转换。TAI也忽略了闰秒，所以TAI和gps时间应该总是以相同的偏移量演化.在工作中，我们使用gps-时间作为时间标准来同步世界各地的天文观测.
我有在gps-time和UTC之间转换的工作函数，但是我不得不对一个闰秒表进行硬编码，得到这里 (文件tzdata2013xx.tar.gz包含一个名为leapseconds的文件)。我必须手动更新这个文件，每隔几年，当一个新的飞跃被宣布。我更愿意从标准的tzdata中获取这些信息，该数据每年通过系统更新几次自动更新。
我很确定这些信息隐藏在/usr/share/zoneinfo/中的一些二进制文件中。我已经能够使用struct.unpack提取其中的一些内容(man tzfile给出了一些有关格式的信息)，但我从未让它完全工作。是否有任何标准包可以访问此信息？我知道比兹，它似乎从同一个数据库获取标准的DST信息，但它不允许访问闰秒。我也找到了tai64n，但是看看它的源代码，它只包含一个硬编码的表。
编辑
在steveha的回答和pytz/tzfile.py中的一些代码的启发下，我终于得到了一个工作解决方案(在py2.5和py2.7上进行了测试)：
from struct import unpack, calcsize
from datetime import datetime
def print_leap(tzfile = '/usr/share/zoneinfo/right/UTC'):
with open(tzfile, 'rb') as f:
# read header
fmt = '>4s c 15x 6l'
(magic, format, ttisgmtcnt, ttisstdcnt,leapcnt, timecnt,
typecnt, charcnt) = unpack(fmt, f.read(calcsize(fmt)))
assert magic == 'TZif'.encode('US-ASCII'), 'Not a timezone file'
print 'Found %i leapseconds:' % leapcnt
# skip over some uninteresting data
fmt = '>%(timecnt)dl %(timecnt)dB %(ttinfo)s %(charcnt)ds' % dict(
timecnt=timecnt, ttinfo='lBB'*typecnt, charcnt=charcnt)
f.read(calcsize(fmt))
#read leap-seconds
fmt = '>2l'
for i in xrange(leapcnt):
tleap, nleap = unpack(fmt, f.read(calcsize(fmt)))
print datetime.utcfromtimestamp(tleap-nleap+1)有结果
In [2]: print_leap()
Found 25 leapseconds:
1972-07-01 00:00:00
1973-01-01 00:00:00
1974-01-01 00:00:00
...
2006-01-01 00:00:00
2009-01-01 00:00:00
2012-07-01 00:00:00虽然这确实解决了我的问题，但我可能不会赞成这个解决办法。相反，我将按照Matt的建议，在代码中包括leap-seconds.list。这似乎是用作tzdata源的权威列表，并且可能每年由NIST更新两次。这意味着我必须手工进行更新，但是这个文件很容易解析，并且包含一个过期日期( tzdata似乎缺少这个日期)。
回答 2
Stack Overflow用户
发布于 2013-10-29 01:02:55
我只是做了man 5 tzfile，并计算了一个偏移量，可以找到闰秒信息，然后读取闰秒信息。
您可以取消对"DEBUG:“print语句的注释，以查看它在文件中找到的更多内容。
编辑:程序更新到现在是正确的。它现在使用文件/usr/share/zoneinfo/right/UTC，现在它找到了要打印的闰秒。
最初的程序没有跳过说明页中记录的timezeone缩略语字符，但有点隐藏("...and tt_abbrind用作文件中遵循ttinfo结构的时区缩写字符数组的索引。“)
import datetime
import struct
TZFILE_MAGIC = 'TZif'.encode('US-ASCII')
def leap_seconds(f):
"""
Return a list of tuples of this format: (timestamp, number_of_seconds)
timestamp: a 32-bit timestamp, seconds since the UNIX epoch
number_of_seconds: how many leap-seconds occur at timestamp
"""
fmt = ">4s c 15x 6l"
size = struct.calcsize(fmt)
(tzfile_magic, tzfile_format, ttisgmtcnt, ttisstdcnt, leapcnt, timecnt,
typecnt, charcnt) = struct.unpack(fmt, f.read(size))
#print("DEBUG: tzfile_magic: {} tzfile_format: {} ttisgmtcnt: {} ttisstdcnt: {} leapcnt: {} timecnt: {} typecnt: {} charcnt: {}".format(tzfile_magic, tzfile_format, ttisgmtcnt, ttisstdcnt, leapcnt, timecnt, typecnt, charcnt))
# Make sure it is a tzfile(5) file
assert tzfile_magic == TZFILE_MAGIC, (
"Not a tzfile; file magic was: '{}'".format(tzfile_magic))
# comments below show struct codes such as "l" for 32-bit long integer
offset = (timecnt*4 # transition times, each "l"
+ timecnt*1 # indices tying transition time to ttinfo values, each "B"
+ typecnt*6 # ttinfo structs, each stored as "lBB"
+ charcnt*1) # timezone abbreviation chars, each "c"
f.seek(offset, 1) # seek offset bytes from current position
fmt = '>{}l'.format(leapcnt*2)
#print("DEBUG: leapcnt: {} fmt: '{}'".format(leapcnt, fmt))
size = struct.calcsize(fmt)
data = struct.unpack(fmt, f.read(size))
lst = [(data[i], data[i+1]) for i in range(0, len(data), 2)]
assert all(lst[i][0] < lst[i+1][0] for i in range(len(lst)-1))
assert all(lst[i][1] == lst[i+1][1]-1 for i in range(len(lst)-1))
return lst
def print_leaps(leap_lst):
# leap_lst is tuples: (timestamp, num_leap_seconds)
for ts, num_secs in leap_lst:
print(datetime.datetime.utcfromtimestamp(ts - num_secs+1))
if __name__ == '__main__':
import os
zoneinfo_fname = '/usr/share/zoneinfo/right/UTC'
with open(zoneinfo_fname, 'rb') as f:
leap_lst = leap_seconds(f)
print_leaps(leap_lst)Stack Overflow用户
发布于 2013-10-28 23:05:07
PyEphem有一个delta_t函数，它返回陆地时间和世界时间(秒)之间的差值。您可以从其中减去32.184以获得飞跃秒(参考)。
import ephem, datetime
ephem.delta_t(datetime.datetime.now()) - 32.184
Out[2]: 35.01972996360122https://stackoverflow.com/questions/19332902
复制相似问题
腾讯云开发者
