常用的几种时间系统GMT、UTC等以及使用Python代码GPS时间转UTC时间

内容转载自我的博客

时间标准

• 格林尼治标准时间GMT（Greenwich Mean Time）：它是指位于英国伦敦郊区的皇家格林尼治天文台当地的平太阳时，因为本初子午线被定义为通过那里的经线。自1924年2月5日开始，格林尼治天文台负责每隔一小时向全世界发放调时信息。格林尼治标准时间的正午是指当平太阳横穿格林尼治子午线时（也就是在格林尼治上空最高点时）的时间。由于地球每天的自转是有些不规则的，而且正在缓慢减速，因此格林尼治平时基于天文观测本身的缺陷，已经被原子钟报时的协调世界时（UTC）所取代
• 世界时UT（Universal Time）：它是一种以格林威治子夜起算的平太阳时。世界时是以地球自转为基准得到的时间尺度，其精度受到地球自转不均匀变化和极移的影响，为了解决这种影响，1955年国际天文联合会定义了UT0、UT1和UT2三个系统：UTO表示未经改正的世界时，UT1表示经过极移改正的世界时，UT2表示进一步经过地球自转速度的季节性改正后的世界时
• 国际原子时TAI（International Atomic Time）：英语简称为IAT，法语简称为TAI，最终使用法语简称TAI。它是根据秒的定义（1秒为铯-133原子基态两个超精细能级间跃迁辐射振荡9192631770周所持续的时间）的一种国际参照时标，是均匀的时间尺度，与地球的空间位置不关联，起点为1958年1月1日0时0分0秒。随着时间的迁延，TAI和UT1两种时间尺度的时间差越来越大。为此，在UT1和TAI之间进行协调，这就产生了协调世界时UTC
• 协调世界时UTC（Coordinated Universal Time）：它是最主要的世界时间标准，基于国际原子时，并通过不规则的加入闰秒来抵消地球自转变慢的影响。它又称为世界统一时间、世界标准时间、国际协调时间。由于英文（CUT）和法文（TUC）的缩写不同，作为妥协，简称UTC。协调世界时（UTC）是世界上调节时钟和时间的主要时间标准，是最接近格林威治标准时间（GMT）的几个替代时间系统之一。对于大多数用途来说，UTC时间被认为能与GMT时间互换。人们对该时间系统进行过数次调整，直到1972年引入了闰秒机制，调整工作得以简化。闰秒在必要的时候会被插入到UTC中，以保证协调世界时（UTC）与世界时（UT1）相差不超过0.9秒
• Unix或POSIX时间戳：它是UNIX或类UNIX系统使用的时间表示方式。一般定义为从协调世界时(UTC时间)1970年1月1日0时0分0秒起至现在的总秒数（10位是精确到秒，13位是精确到毫秒）。考虑到闰秒的话，更精确的定义为从协调世界时(UTC时间)1970年1月1日0时0分0秒起至现在经过闰秒调整之后的总秒数
• GPS时间：它是GPS原子时，它的时间基准是1980年1月6日0时0分0秒与世界协调时刻UTC相一致，以后按原子时TAI秒长累积计时。GPS时间跟UTC时间之差为秒的整倍数（因为UTC比TAI慢，而GPS是按照TAI来计时的，所以UTC也比TAI慢）。如1989年为5秒，1996年为11秒，2002年为13秒，到现在2020年08月为止为18秒。

闰秒

它是在协调世界时（UTC）中增加或减少一秒，使它与平太阳时贴近所做调整。需要闰秒的部分原因是因为平均太阳日（mean solar day）的长度正以非常缓慢的速度增加中，另一个原因是原子钟赋予秒固定的时间长度。而当两者结合时，就已经比当时的太阳时的秒短少了一点点。时间现在是以稳定的原子钟来测量（TAI或国际原子时），因为地球自转有着许多的变数，所以以前的秒定义（地球绕着轴自转和绕太阳的公转，以平均太阳日的1/86400来定义）被废除。当要增加正闰秒时，这一秒是增加在第二天的00:00:00之前，效果是延缓UTC第二天的开始。当天23:59:59的下一秒被记为23:59:60，然后才是第二天的00:00:00。如果是负闰秒的话，23:59:58的下一秒就是第二天的00:00:00了，但目前还没有负闰秒调整的需求
Leap_Second.dat可以查看每次闰秒的时间和第几次闰秒，主要内容如下：

#    MJD        Date        TAI-UTC (s)
#           day month year
#    ---    --------------   ------   
#
    41317.0    1  1 1972       10
    41499.0    1  7 1972       11
    41683.0    1  1 1973       12
    42048.0    1  1 1974       13
    42413.0    1  1 1975       14
    42778.0    1  1 1976       15
    43144.0    1  1 1977       16
    43509.0    1  1 1978       17
    43874.0    1  1 1979       18
    44239.0    1  1 1980       19
    44786.0    1  7 1981       20
    45151.0    1  7 1982       21
    45516.0    1  7 1983       22
    46247.0    1  7 1985       23
    47161.0    1  1 1988       24
    47892.0    1  1 1990       25
    48257.0    1  1 1991       26
    48804.0    1  7 1992       27
    49169.0    1  7 1993       28
    49534.0    1  7 1994       29
    50083.0    1  1 1996       30
    50630.0    1  7 1997       31
    51179.0    1  1 1999       32
    53736.0    1  1 2006       33
    54832.0    1  1 2009       34
    56109.0    1  7 2012       35
    57204.0    1  7 2015       36
    57754.0    1  1 2017       37

TAI、GPST、UTC换算

根据文件Leap_Second.dat可以得知，截止2020年08月：
TAI = UTC + 37
由于GPST从1980年1月6日0时0分0秒开始计时，所以1980年1月1日及以前的闰秒不考虑，则：
GPST = UTC + 18
另外，leapsecond可以查看实时的UTC、GPST、TAI时间

时区划分

时区是指地球上的某一个区域使用同一个时间定义。GMT时间或者UT时间，都是表示地球自转速率的一种形式。从太阳升起到太阳落下，时刻从0到24变化。这样，不同经度的地方时间自然会不相同。为了解决这个问题，人们把地球按经度划分为不同的区域，每个区域内使用同一个时间定义，相邻的区域时间差为1个小时。时区又分为理论时区和法定时区

• 理论时区：按经度，每15°为一个时区，将地球划分为24个时区，以本初子午线为中心，向东西两侧各延伸7.5°的区域为0时区
• 法定时区：法定时区是在理论时区的基础上，根据某些地区的国界线做了调整之后的时区。为实际使用的时区。例如中国横跨东五区到东九区五个时区，但统一使用东八区时间（北京时间）
• 时差：某个地方的时刻与0时区的时刻差称为时差，时差东正西负。以本初子午线为中心，每向东一跨过一个时区，时刻增加一个小时，每向西跨过一个时区，时刻减少一个小时
• 国际日期变更线：大体以180度经线为日界线。当自西向东穿过日期变更线时，日期需要减少一天，反之，日期增加一天

UTC时间转本地时间（东八区）

根据此代码，修改最后一行fromtimestamp(timestamp, timezone(timedelta(hours=8)))即可实现UTC转不同时区

from datetime import datetime, timedelta, timezone

# UTC时间转本地时间（北京）时间
# 1. 把utc的str转为datetime（无时区信息）
# 2. 添加时区信息为utc时区
# 3. datetime转为时间戳
# 4. 从时间戳得到本地时间datetime
# 输入格式为：'2020-08-05 02:03:03.815650'
# 输出格式为：datetime.datetime(2020, 8, 5, 10, 3, 3, 815650)
def utc_to_local(utc_time):
    datetimeformat = "%Y-%m-%d %H:%M:%S.%f"
    # 得到不包含时区的datetime
    dt_no_tz = datetime.strptime(utc_time, datetimeformat)
    # 设置时区为UTC
    # timezone.utc与timezone(timedelta(hours=0))一样
    utc_datetime = dt_no_tz.replace(tzinfo=timezone(timedelta(hours=0)))
    t = utc_datetime.timestamp()
    # 根据时间戳得到UTC时间
    # datetime.utcfromtimestamp(t)
    # 如果要将时间戳转化为东八区datetime
    # fromtimestamp(timestamp, timezone(timedelta(hours=8)))
    # 根据时间戳得到本地时间fromtimestamp(t, tz=None)
    return datetime.fromtimestamp(t)

本地时间（东八区）转UTC

根据此代码，修改replace(tzinfo=timezone(timedelta(hours=8)))即可实现某个时区转UTC时间

from datetime import datetime, timedelta, timezone

# 本地时间转UTC时间
# 输入格式为：'2020-08-05 10:03:03.815650'
# 输出格式为：datetime.datetime(2020, 8, 5, 2, 3, 3, 815650)
def local_to_utc(local_time):
    datetimeformat = "%Y-%m-%d %H:%M:%S.%f"
    # 得到不包含时区的datetime
    dt_no_tz = datetime.strptime(local_time, datetimeformat)
    # 设置时区为本地时区（北京，东八区）
    # timezone.utc与timezone(timedelta(hours=0))一样
    local_datetime = dt_no_tz.replace(tzinfo=timezone(timedelta(hours=8)))
    t = local_datetime.timestamp()
    # 根据时间戳得到UTC时间
    return datetime.utcfromtimestamp(t)

GPS时间转UTC时间

将GPS时间转换为UTC时间

from datetime import datetime, timedelta

# 闰秒
LEAP_SECONDS = 18

# 输入：GPS周、GPS周内秒、闰秒（可选，gps时间不同，闰秒值也不同，由Leap_Second.dat文件决定）
# 输出：UTC时间（格林尼治时间）
# 输入示例： gps_week_seconds_to_utc(2119, 214365.000)
# 输出示例： '2020-08-18 11:32:27.000000'
def gps_week_seconds_to_utc(gpsweek, gpsseconds, leapseconds=LEAP_SECONDS):
    datetimeformat = "%Y-%m-%d %H:%M:%S.%f"
    epoch = datetime.strptime("1980-01-06 00:00:00.000", datetimeformat)
    # timedelta函数会处理seconds为负数的情况
    elapsed = timedelta(days=(gpsweek*7), seconds=(gpsseconds-leapseconds))
    return datetime.strftime(epoch+elapsed, datetimeformat)

UTC时间转GPS时间

将UTC时间转换为GPS时间

from datetime import datetime, timedelta

# 闰秒
LEAP_SECONDS = 18

# 输入：UTC时间（datetime类型）
# 输出：GPS周、周内日、周内秒、毫秒
def utc_to_gps_week_seconds(utc, leapseconds=LEAP_SECONDS):
    datetimeformat = "%Y-%m-%d %H:%M:%S.%f"
    epoch = datetime.strptime("1980-01-06 00:00:00.000", datetimeformat)
    tdiff = utc - epoch + timedelta(seconds=leapseconds)
    gpsweek = tdiff.days // 7
    gpsdays = tdiff.days - 7*gpsweek
    gpsseconds = tdiff.seconds + 86400*(tdiff.days -7*gpsweek)
    return gpsweek, gpsdays, gpsseconds, tdiff.microseconds