RTKlib中的GPS秒转历法格里历的表示方法

RTKlib中的GPS秒转历法格里历的表示方法

在RTKlib里面有个time2epoch方法，按其注释就是convert gtime_t struct to calendar day/time。

/* time to calendar day/time ---------------------------------------------------
* convert gtime_t struct to calendar day/time
* args   : gtime_t t        I   gtime_t struct
*          double *ep       O   day/time {year,month,day,hour,min,sec}
* return : none
* notes  : proper in 1970-2037 or 1970-2099 (64bit time_t)
*-----------------------------------------------------------------------------*/
extern void time2epoch(gtime_t t, double *ep)
{
    const int mday[]={ /* # of days in a month */
        31,28,31,30,31,30,31,31,30,31,30,31,31,28,31,30,31,30,31,31,30,31,30,31,
        31,29,31,30,31,30,31,31,30,31,30,31,31,28,31,30,31,30,31,31,30,31,30,31
    };
    int days,sec,mon,day;
    
    /* leap year if year%4==0 in 1901-2099 */
    days=(int)(t.time/86400);
    sec=(int)(t.time-(time_t)days*86400);
    for (day=days%1461,mon=0;mon<48;mon++) {
        if (day>=mday[mon]) day-=mday[mon]; else break;
    }
    ep[0]=1970+days/1461*4+mon/12; ep[1]=mon%12+1; ep[2]=day+1;
    ep[3]=sec/3600; ep[4]=sec%3600/60; ep[5]=sec%60+t.sec;
}

首先，它定义了一个数组mday，仔细看看一共有48个元素，对应着4年的48个月，从1970年算起，1970、1971、1972、1973，四年一个循环，为什么是四个月呢？我们要考虑闰年，一般来说非世纪年只要能被4整除就是闰年，世纪年要被400整除才是闰年，1970年至当前2020年，只有一个世纪年2000年，它也是闰年，所以不用考虑世纪年是不是闰年的问题，只用考虑只要被4整除就是闰年了。

从1970年开始，第一个闰年是1972年，也就是mday数组里面的第三年，它的二月有29天，所以说从1972年开始每四年就是闰年。
定义了四个int变量days,sec,mon,day：

days=(int)(t.time/86400);

days是天数，即从GPS秒除86400（即一天的24h）取整，即为天数。

sec=(int)(t.time-(time_t)days*86400);

sec是秒数，即天内秒。

for (day=days%1461,mon=0;mon<48;mon++) {
        if (day>=mday[mon]) day-=mday[mon]; else break;
    }

mon是月数，这里的月数是取0-47。

day是模1461后的剩余天数，这里的1461即为mday数组里四年48个月的总天数。用day去和mday数组里每个月的天数比较，大于或等于的话，就减去这个月的天数，去和下个月比较，直到day小于当前元素对应的天数，那么时间就是这个月的。

ep[0]=1970+days/1461*4+mon/12; //年数
ep[1]=mon%12+1; //月
ep[2]=day+1;//日
ep[3]=sec/3600; //时
ep[4]=sec%3600/60; //分
ep[5]=sec%60+t.sec;//秒

mon因为是0-47，所以模12后要加1。
day因为是通过days求得的，days是取整天数后的值，所以应该是第二天，所以要加1。