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。