【GNSS】单点定位程序总结 二 时间系统模块

最新推荐文章于 2024-05-15 11:58:26 发布
最新推荐文章于 2024-05-15 11:58:26 发布
阅读量820 收藏 15
点赞数 3
分类专栏： C++ GNSS 文章标签： gnss
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_32109917/article/details/113689905
版权
GNSS 同时被 2 个专栏收录
42 篇文章 38 订阅
订阅专栏
C++
26 篇文章 7 订阅
订阅专栏

宝藏就在这里，就看还在上这门课的学弟学妹能不能发现了

时间系统模块

1.1儒略日

时间定义与表示

儒略日（Julian day，JD）是指由公元前4713年1月1日，协调世界时中午12时开始所经过的天数，多为天文学家采用，用以作为天文学的单一历法，把不同历法的年表统一起来。如果计算相隔若干年的两个日期之间间隔的天数，利用儒略日就比较方便。

在程序中，用结构体JDTIME表示儒略日时间，结构体定义如下

1

2

3

4

5

6

7

8

9

10

11

12

/* 儒略日; */

struct JDTIME

{

    int Days;

    double FracDay;

 

    JDTIME()

    {

        Days = 0;

        FracDay = 0.0;

    }

};

 

1.2简化儒略日

 

1

2

3

4

5

6

7

8

9

10

11

12

/* 简化儒略日; */

struct MJDTIME

{

    int Days;

    double FracDay;

 

    MJDTIME()

    {

        Days = 0;

        FracDay = 0.0;

    }

};

 

时间定义与表示

1.3通用时

时间定义与表示

1

2

3

4

5

6

7

8

9

10

/* 通用时间定义; */

struct COMMONTIME

{

    unsigned short Year;

    unsigned short Month;

    unsigned short Day;

    unsigned short Hour;

    unsigned short Minute;

    double Second;

};

 

1.4GPST

时间定义与表示

1

2

3

4

5

6

7

8

9

10

11

12

13

/* GPS时间定义 */

struct GPSTIME

{

    unsigned short Week;

    double

        SecOfWeek;

 

    GPSTIME()

    {

        Week = 0;

        SecOfWeek = 0.0;

    }

};

 

1.5年积日

无

时间定义与表示

1.3时间系统相互转换公式与代码

1.3.1儒略日与简化儒略日

1.儒略日到简化儒略日转换公式

假设简化儒略日表示为MJD ，儒略日表示为JD 。我们可以得到JD=MJD-2400000.5

3.转换程序代码

儒略日到简化儒略日

1

2

3

4

5

6

7

bool JDTime2MJDTime(const JDTIME JD, MJDTIME*MJD)

{

 

    MJD->Days = JD.Days - 2400000 - (int)(JD.FracDay - 0.5-1);

    MJD->FracDay = FRAC(JD.FracDay + 1 - 0.5);

    return true;

}

简化儒略日到儒略日

1

2

3

4

5

6

7

bool MJDTime2JDTime(const MJDTIME MJD, JDTIME*JD)

{

   

    JD->Days = MJD.Days + 2400000 + (int)(MJD.FracDay + 0.5);

    JD->FracDay = FRAC(MJD.FracDay+0.5);

    return true;

}

1.3.2儒略日与通用时

1.儒略日到通用时转换公式

2.通用时到儒略日转换公式

JD=INT[365.25y]+INT[30.6001(m+1)]+D+UT/24+1720981.5

MJD=JD-2400000.5

其中：

       如果M≤2 ，则y=Y-1,m=M+12

       如果M>2 ，则y=Y,m=M

       JD 为儒略日

       Y 为年，M 为月， D为日，UT 为世界时

       INT[] 表示取实数的整数部分

3.转换程序代码

儒略日到通用时

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

bool JDTime2CommonTime(const JDTIME JD, COMMONTIME*Common)

{

    MJDTIME MJD;

    MJD.Days = JD.Days - 2400000 - (int)(JD.FracDay - 0.5 - 1);

    MJD.FracDay = FRAC(JD.FracDay + 1 - 0.5);

    int a, b, c, d, e = 0;

    double UTDay = 0;

    double UT = 0;

    a = (int)(MJD.Days + MJD.FracDay + 2400000.5 + 0.5);

    b = a + 1537;

    c = (int)((1.0*b - 122.1) / (365.25));

    d = (int)(365.25*c);

    e = (int)((b - d)*1.0 / 30.6001);

    UTDay = b - d - (int)(30.6001*e) + FRAC(MJD.Days + MJD.FracDay + 2400000.5 + 0.5);

    Common->Day = (int)(UTDay);

    Common->Month = e - 1 - 12 * (int)(e*1.0 / 14);

    Common->Year = c - 4715 - (int)((7 + Common->Month)*1.0 / 10.0);

    UT = (FRAC(UTDay) * 24.0);

    Common->Hour = int(MJD.FracDay * 24);

    Common->Minute = int((MJD.FracDay * 24 - Common->Hour) * 60);

    Common->Second = MJD.FracDay * 24 * 3600 - Common->Hour * 3600 - Common->Minute * 60;

    return true;

}

 

通用时到儒略日

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

bool CommonTime2JDTime(const COMMONTIME Common, JDTIME*JD)

{

    int y, m = 0;

    double TemptJD, UT = 0;

    if (Common.Month <= 2)

    {

        y = Common.Year - 1;

        m = Common.Month + 12;

    }

    if (Common.Month > 2)

    {

        y = Common.Year;

        m = Common.Month;

    }

    UT = Common.Hour + Common.Minute / 60.0 + Common.Second / 3600.0;

    TemptJD = (int)(365.25*y) + (int)(30.6001*(m + 1)) + Common.Day + UT / 24.0 + 1720981.5;

    JD->Days = (int)(TemptJD);

    JD->FracDay = (TemptJD)-JD->Days;

    return true;

}

1.3.3儒略日与GPST

1.儒略日到GPST转换公式

先将儒略日转换为简化儒略日

MJD=JD+2400000.5

再从简化儒略日转换GPST

GPSWEEK=INT[(MJD-44244)/7]
SECOFWEEK=(MJD-44244-GPSWEEK×7)×86400

2.GPST到儒略日转换公式

MJD=44244+GPSWEEK×7+SECOFWEEK/86400

3.转换程序代码

儒略日到GPST

1

2

3

4

5

6

bool JDTime2GPSTime(const JDTIME JD, GPSTIME*GPS)

{

    GPS->Week = (int)((JD.Days + JD.FracDay - 2400000.5 - 44244) / 7.0);

    GPS->SecOfWeek = (JD.Days + JD.FracDay - 2400000.5 - 44244 - GPS->Week * 7) * 86400;

    return true;

}

 

GPST到儒略日

1

2

3

4

5

6

7

8

bool GPSTime2JDTime(const GPSTIME GPS, JDTIME*JD)

{

    double temptJD;

    temptJD = 44244 + GPS.Week * 7 + GPS.SecOfWeek / 86400 + 2400000.5;

    JD->Days = (int)(temptJD);

    JD->FracDay = FRAC(temptJD);

    return true;

}

1.3.4简化儒略日与通用时

1.简化儒略日到通用时转换公式

JD=MJD+2400000.5,
a=INT[JD+0.5]
b=a+1537
c=INT[(b-122.1)/365.25]
d=INT[365.25×c]
e=INT[(b-d)/30.6001]
Day=b-d-INT[30.6001×e]+FRAC[JD+0.5]
Month=e-1-12×INT[e/14]
Year=c-4715-INT[(7+M)/10]
Number=mod{INT[JD+0.5],7}(Number=0,Monday⋯Number=7,Sunday)
Hour=INT{FRAC[MJD]×24}
Minute=INT{FRAC[MJD]×24×60-Hour×60}
Second=INT{FRAC[MJD]×24×3600-Hour×3600-Minute×60}

2.通用时到简化儒略日转换公式

MJD=INT[365.25×y]+INT[30.6001×(m+1)]+Day+UT/24-679019

其中：

       如果M≤2 ，则y=Year-1,m=Mouth+12

       如果M>2 ，则y=Year,m=Mouth

       UT表示世界时

3.转换程序代码

简化儒略日到通用时

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

bool MJDTime2CommTime(const MJDTIME MJD, COMMONTIME*Common)

{

    int a, b, c, d, e = 0;

    double UTDay = 0;

    double UT = 0;

    a = (int)(MJD.Days+MJD.FracDay + 2400000.5+0.5);

    b = a + 1537;

    c = (int)((1.0*b - 122.1) / (365.25));

    d = (int)(365.25*c);

    e = (int)((b - d)*1.0 / 30.6001);

    UTDay = b - d - (int)(30.6001*e) + FRAC(MJD.Days + MJD.FracDay + 2400000.5 + 0.5);

    Common->Day = (int)(UTDay);

    Common->Month = e - 1 - 12 * (int)(e*1.0 / 14);

    Common->Year = c - 4715 - (int)((7 + Common->Month)*1.0 / 10.0);

    UT = (FRAC(UTDay) * 24.0);

    Common->Hour = int(MJD.FracDay * 24);

    Common->Minute = int((MJD.FracDay * 24 - Common->Hour) * 60);

    Common->Second = MJD.FracDay * 24 * 3600 - Common->Hour * 3600 - Common->Minute * 60;

    return true;

}

通用时到简化儒略日

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

bool CommonTime2MJDTime(const COMMONTIME Common, MJDTIME*MJD)

{

    int y, m = 0;

    double JD, UT = 0;

    JDTIME  jd;

 

    if (Common.Month<=2)

    {

        y = Common.Year - 1;

        m = Common.Month + 12;

    }

    if (Common.Month>2)

    {

        y = Common.Year;

        m = Common.Month;

    }

    UT = Common.Hour + Common.Minute / 60.0 + Common.Second / 3600.0;

    JD = (int)(365.25*y) + (int)(30.6001*(m + 1)) + Common.Day + UT / 24.0 + 1720981.5;

    MJD->Days = (int)(365.25*y) + (int)(30.6001*(m + 1)) + Common.Day + UT / 24.0 + 1720981.5-2400000.5;

    MJD->FracDay = (UT / 24.0 ) - (int)(UT / 24.0 );

 

    return true;

}

 

 

小马甲的新马甲
关注
  • 3
    点赞
  • 踩
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
gps时间转换，包括年积儒略日，简化儒略日，gps周，年月日等
11-06
gps时间转换，包括年积儒略日，简化儒略日，gps周gps时间转换，包括年积儒略日，简化儒略日，gps周，年月日等
MATLAB GNSS单点定位程序
04-02
MATLAB GNSS单点定位程序
GPS时、北斗时、儒略日时间转换.zip
03-20
C++ 代码：GPS时、北斗时、儒略日转换
空间大地测量与GPS导航定位时间系统相互转换,格里高利时通用时儒略日,GPS时,年积日相互转换
FreeBird的专栏
08-26 5789
进行空间大地测量所需时间系统的相互转换                                      1.时间标示法 从古至今，世界各国曾出现过多种历法，如罗马历、儒略日、格里高利以及我国的农历等，目前世界上广泛采用的历法是格里高利历。 1.1.格里高利时（通用时） 格里高利历也称公历，现被世界各国广泛采用。格里高利历是一个由146097天所组成的400年周期为基础，1年的平
时钟服务器之GNSS_gnss时间系统(1)
最新发布
2301_79055239的博客
05-15 400
全球导航卫星系统（Global Navigation Satellite System，GNSS），又称全球卫星导航系统，是能在地球表面或近地空间的任何地点为用户提供全天候的3维坐标和速度以及时间信息的空基无线电导航定位系统。其包括一个或多个卫星星座及其支持特定工作所需的增强系统。全球卫星导航系统国际委员会公布的全球4大卫星导航系统供应商，包括中国的北斗卫星导航系统（BDS）、美国的全球定位系统（GPS）、俄罗斯的格洛纳斯卫星导航系统（GLONASS）和欧盟的伽利略卫星导航系统（GALILEO）。
【GNSS】GNSS处理中的时间系统
小马甲的新马甲
01-30 3027
精密GNSS数据处理中主要涉及三种时间类型：以地球自转为基础建立的世界时系统；基于物质内部的原子运动特性的原子时系统；基于天体动力学理论的动力学时系统。 概念介绍 世界时（UT） 世界时是指通过观测恒星的视运动表达的时间系统，如恒星时（ST）指以春分点为参考连续两次经过地方子午圈的时间间隔。格林尼治平恒星时（GMST）消去了章动影响。 同理，以太阳为参考称为太阳时。以真太阳时的平均角速度在天球赤道上作周年视运动的 平太阳作为参考的为平太阳时。若平太阳时起点为平子夜，称为世界时（UT0）。UT0 进行 极移改
GPS涉及到的各种时间转换（年月日，年积日，儒略日，GPS周及周内日或周内秒，星期几）python
热门推荐
qq_40185784的博客
12-08 1万+
因为写一些程序或脚本经常要涉及GPS的时间转换，比如下载各种产品遍历的时候可能需要用年积日，但是如果日期有跨年的话年积日也不方便，可能就需要先转儒略日再转对应的年月日，而有些产品命名又是GPS周相关，所以就写了一个转换时间的代码，方便每次的调用。
GNSS单点定位.rar
09-19
%读取RINEX V3.04电文文件和观测文件 %可分别对GPS Galileo BD GLONASS进行单点定位（C1C） %GLONASS处理数据为 2014年以后的数据（UTC跳秒为16 或 18）
GNSS精密单点定位代码_单点定位_精密单点定位_GNSS_精密单点定位（PPP）_
10-02
给出了精密单点定位的代码，方便使用者在此基础上实现其他的改进方法。
双系统双频伪距单点定位程序
07-23
该程序采用C/C++（一开始是C++，后来为了做嵌入式，改成了C）语言编写而成，支持RENIX3.04文件和OEM7二进制数据文件的SPP解算，可自由选择GPS、BDS系统，可自由选择单频、双频、双频无电离层组合SPP解算。该程序简单...
GNSS+精密单点定位+IGS产品使用说明+Kouba
04-15
【2015-Kouba】A GUIDE TO USING INTERNATIONAL GNSS SERVICE (IGS) PRODUCTS 原文，可能有的朋友找不到，特此分享一下。 文章摘要如下： Since 1994, the International GNSS Service (IGS) has provided precise ...
公历与简化儒略日转换(GPS)
12-10
运用C#窗体编写的一个小应用程序，仅供娱乐用途
TimeConvert_儒略日_儒略日utc_timeconvert_时间转换_GPST与UTC_
09-29
GPST、UTC、儒略日等之间的时间转换。
格里历、儒略日、GPS时、Glonass时、北斗时、年积日互转
11-25
C# 时间转换 代码
GNSS精密单点定位技术.pptx
06-26
GNSS精密单点定位技术概述PPT，详细介绍了PPP技术的原理、研究进展、开源软件等相关信息，对PPP技术的未来发展方向进行了分析总结。
卫星定位授时系统的时间表示主要有以下几种方法：
我们编程吧
03-14 4679
卫星定位授时系统的时间表示主要有以下几种方法： 1、通用时间表示方法： 日历表示法：年、月、日、时、分、秒 2、儒略日（JULIAN DAY）：是指从-4712年1月1日（即公元前4713年1月1日）正午开始的天数。 由J.J.Scaliger在1583年提出。它的特点是连续，利于数学表达，但是不直观。 3、新儒略日（Modified Julian Day –MJD）:从儒略日中减去2,400,000.5天来得到，给出的是从1858年11月17日子夜开始的天数。 特点是数值比儒略日小。 4、年
GNSS时间转化（GPS时-儒略日-年月日时分秒）
weixin_45858729的博客
03-06 2350
使用python进行时间转换代码如下： 代码引用自：https://blog.csdn.net/qq_40185784/article/details/104956594 import math def cal2mjd(cal): # cal2jd 将公历年月日时分秒转换到简化儒略日。 # 输入公历时间列表，返回儒略日 if (len(cal) < 6): for i in range(len(cal), 6): cal.app..
年月日、年积日、简化儒略日、GPS周之间相互转换（C++）
qq_41696018的博客
08-12 7656
在学习GNSS过程中，不管是数据下载又或是不同文件名中时间的匹配，经常会遇到需要在不同时间之间转换，笔者在此归纳了几种经常会遇到的时间转换： 1、年月日转年积日 3、年积日转简化儒略日 5、简化儒略日转GPS周 2、年积日转年月日 4、简化儒略日转年积日 6、GPS周转简化儒略日 笔者在此同时分享一个小软件，用于帮助大家检验通过代码转换的结果：链接：h...
gnss单点定位算法代码
07-14
GNSS（全球导航卫星系统）单点定位算法是一种利用卫星信号进行定位的方法，主要通过接收多颗卫星的信号进行测量和计算来确定用户的位置。下面是一个示例的GNSS单点定位算法的代码： ```python import numpy as np ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值