How to convert the number of seconds from an epoch into UTC using the datetime function.
152 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Matthew Casiano
il 30 Lug 2021
Commentato: James Tursa
il 26 Nov 2024 alle 18:43
I am struggling to convert the number of seconds from an epoch (1980-01-06 GPS epoch) into UTC using the datetime function. The following commands give me a time of 18-Mar-2021 20:37:31, however in UTC this should be 18-Mar-2021 20:37:13 since I would need to subtract out the number of leap seconds from the total (currently there are 18 leap seconds since the 1980-01-06 GPS epoch).
NoOfSeconds=1300135051;
datetime(NoOfSeconds,'convertfrom','epochtime','Epoch','1980-01-06')
Can the datetime function automatically account for leapseconds in this conversion? I tried using the 'TimeZone','UTCLeapSeconds' name/value pair, but this adds 18 seconds since it assumes the time is already in UTC.
0 Commenti
Risposta accettata
Peter Perkins
il 2 Ago 2021
Matthew, this is easier than you think. In MATLAB, 'TimeZone','UTC' is sort of a ficticious timeline that pretends that leaps seconds don't exist (because hardly anyone wants to actually account for leap seconds). 'UTCLeapSeconds' is the timeline that knows all about leap seconds. The GPS timeline (in the real world) ticks off the same leap seconds as 'UTCLeapSeconds' in MATLAB does, it just doesn't treat them specially. So for the purposes of arithmetic, i.e. elapsed time, 'UTCLeapSeconds' is what you want. Then all you have to do is change to 'UTC' (or not, it's up to you).
>> epoch = datetime(1980,1,6,'TimeZone','UTCLeapSeconds')
epoch =
datetime
1980-01-06T00:00:00.000Z
>> dtUTC = epoch + seconds(1300135051)
dtUTC =
datetime
2021-03-18T20:37:13.000Z
>> dtUTC.TimeZone = 'UTC'
dtUTC =
datetime
18-Mar-2021 20:37:13
Also, this Answer might be helpful:
If I may ask, I would be interested in hearing more about how you work with GPS timestamps. Do you always get "seconds since 1980" as your input? Do you need to work with "week number, seconds within week"? Do you need to work with "human readable timestamps" (e.g. 18-Mar-2021 20:37:31) that refer to the GPS timeline (as opposed to UTC)? Do you need to convert back to one of those GPS time formats? Thanks in advance for anything you are able to share.
8 Commenti
Peter Perkins
il 25 Nov 2024 alle 20:22
James, you are mixing apples and oranges. "UTCLeapSeconds" represents the real-world timeline that official UTC clocks follow, including counting leap seconds. "UTC" is a ficticious timeline that makes life less painful for the overwhelming majority of people who do not want to complicate their life with leap seconds.
2016-12-31T23:59:60.500Z does not exist on the ficticious "UTC" timeline. If I want to convert that to "UTC", should it map to 23:59:59.5 on that timeline, or to 00:00:0.5 on the next day? I would argue that the former creates fewer problems because it does not change the date, and I guess you would argue that since POSIX defines the POSIX time during a leap second as the 0th second of the next day, that's how the conversion should work. But the real answer is that there is no answer. It can't possibly be a lossless conversion.
This
dtutc2 = datetime(1970,1,1,'TimeZone','UTC') + seconds(posixtime(dtgps2))
adds an elapsed time that was calculated under "UTCLeapSeconds" to a "UTC" datetime. They should not be mixed, and it's going to lead to inconsistencies.
It seems like maybe the underlying issue is that UTCLeapSeconds only supports the ISO display format, and you'd like more flexibility, so you are partially using "UTC", and partially "UTCLeapSeconds".
James Tursa
il 26 Nov 2024 alle 18:43
@Peter Perkins I am not trying to get into an argument here about what MATLAB should do in this lossy conversion. I prefer consistency among other MATLAB functions such as posixtime, and you prefer to keep the date the same. OK, fine. But that's not my point. My point is the user should be aware of what datetime does during this conversion and offer a different formulation so they can choose the behavior they want if needed. That's all.
Più risposte (2)
Sarthak
il 13 Mar 2024
function date = gps2utc(g0, g1)
date = datevec(datetime(g0*604800+g1,'convertfrom','epochtime','Epoch','1980-01-06'));
0 Commenti
Vedere anche
Categorie
Scopri di più su Dates and Time in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!