七宝大人-CSDN博客

原创 BDS和GPS、电离层相关SSR数据解码

网上找不到RTCM最新版协议，RTCM2.3、RTCM3.2、RTCM3.3协议中，主要给出了GPS轨道、钟差等SSR信息结构（1057、1068、1060），并没有BDS和电离层SSR的信息结构。在BNC软件中，SSR数据流信息整理的比较完善，按照bnchelp.html，截图如下：对应解码部分可以从源码中移植，对于BDS轨道、钟差SSR，主要比GPS多出TOEMOD信息（BNC：#define G_BDS_TOEMOD(a) GETBITSFACTOR(a...

2021-09-27 14:41:23 1044

原创龙格库塔4阶方法

在groops源码数值积分时看到orbitPropagatorRungeKutta4.h，源码如下：//4级龙格-库塔方法，求解位置，速度，加速度，并输出inline OrbitArc OrbitPropagatorRungeKutta4::integrateArc(OrbitEpoch startEpoch, Time sampling, UInt posCount, ForcesPtr forces,

2021-07-12 18:15:08 2776

原创 groops软件与源码阅读

工作需要，对GROOPS源码进行了一段时间的阅读，参考源码说明，进行了简要的整理工作，后期会更新更详细更细节的源码阅读记录文件。简介GROOPS是由于奥地利格拉茨大学（TU Graz）、德国波恩大学等团队学者开发的开源的面向对象的重力场恢复的程序系统，主要代码为C++ 语言（+部分Fortran代码），功能包括：地球重力场恢复、精密单点定位（PPP）、低轨卫星定轨、局域椭球确定等。GROOPS安装（windows）2.1 获取GROOPS源代码开源代码下载地址：https://github.co

2021-05-10 13:11:51 1976 3

原创 RTCM-SSR信息读取与解码

标准RTCM格式可查看RTCM格式说明书读取流程：1.选择合适的挂载点，接收相应的实时数据流。2.读取、解码实时数据流实时数据流结构如下图所示，解码判断引导字后再解码判断信息类型，按照RTCM对应的信息结构读取解码相应信息。3.改正数信息参考点说明SSR实时轨道改正信息中包含相对于广播星历的卫星位置改正值及速度改正值，用户接收到NTRIP播发的实时改正后，结合对应时刻的广播星历轨道参数，可以反算得到高精度的实时轨道，进行实时定位。其播发的改正值的参考点一般有两种：天线相位中心(APC)和卫星

2021-03-25 12:34:56 2249

原创 20210325BERNESE软件中LEO_POD模块说明

1.流程图如下图所示：2.主要步骤为：下载观测数据、精密星历、精密钟差、ERP文件、DCB文件、ATT文件等。由GPS精密星历建立BERNESE自定义格式的卫星标准轨道。将GPS观测数据转化为BERNESE格式；生成先验轨道，利用伪距数据求解概略的轨道坐标，并将其转化为低轨卫星精密星历SP3C格式，接着再以低轨卫星星历来建立卫星标准轨道；时钟同步；数据预处理，进行周跳探测与修复，检测相位观测质量并标记粗差；参数估计；轨道积分，求解简化动力学轨道。将步骤6的轨道成果文件利用BERNESE

2021-03-25 11:33:03 777

原创 2021-03-10 太阳代码

太阳喝月球void Sun ( const MJDTIME* Mjd_TT, double SunPos[] ){ double eps, T; double L, M, r; double R[9], SunOrb[3]; T = ( Mjd_TT->Days + Mjd_TT->FracDay - MJD_J2000 ) / 36525.0; /* Obliquity of J2000 ecliptic */ eps = 23.43929111 * Rad -( 4

2021-03-10 17:16:36 415

原创 2021-03-08动力学方程

运动方程惯性系下，由卫星的摄动理论可知，卫星的运动方程为：等式右边第一项为中心引力，第二项为摄动力，表示卫星所受的各种摄动力产生的加速度总和。Pd为卫星的动力学模型参数矢量（如光压参数、大气阻力参数等），认为其为常矢量。卫星所受力可以分为保守力和非保守力。图1给出了卫星上所受各种力产生的加速度大小随轨道高度的变化；图2给出了三类卫星各摄动力的量级（LEO、MEO、GEO）；可以看出：N体运动、太阳光压随轨道高度增加而增大；地球中心引力及镶滚摄动力随轨道高度高度增加而减小。MEO和GEO卫星受到的

2021-03-08 17:18:29 1631 4

原创 2021-03-05天线相位中心偏差和天线相位中心变化

天线相位中心偏差GNSS测量是从卫星发射天线的相位中心至接收机天线相位中心间的距离。而IGS精密星历给出的是卫星质心的坐标，两者之间并不一致。

2021-03-05 16:03:32 4551 1

原创 2021-03-04低轨卫星定轨流程整理学习04

各类坐标系及坐标系之间的转换为了减少牵连运动所引起的附加速度，使卫星的运动方程相对简化，一般在惯性坐标系中建立和求解卫星的运动方程。而地球对卫星的引力（地球重力场模型）是在地固坐标系中定义的。因此要建立卫星卫星运动方程和参数估计的观测方程，必须掌握这些坐标系统的定义及其相互转换关系。①2000.0惯性坐标系原点位于地球质心，基准面为2000.0地球平赤道面，X轴在基准面内由地球质心指向2000.0的平春分点。Z轴为基准面的法向。Y轴与X、Z轴成右手系。②地固坐标系...

2021-03-04 17:15:07 1371 1

原创 2021-03-03BERNESE软件定轨批处理详细教程

BERNESE软件定轨批处理详细教程非专业人士，工作需要又处理一天的grace卫星轨道数据，感觉自己之前写的教程自己都有些看不懂。遂在此细写一份grace卫星数据下载、linux下数据格式转换、Bernese安装、定轨批处理所需文件、文件位置与格式要求、结果输出等说明。①grace卫星数据下载建议下载使用软件FILAZILLA进行下载（网页下载容易中断，造成数据不完整）下载地址：ftp://isdcftp.gfz-potsdam.de/grace/Level-1B/JPL/INSTRUMENT/R

2021-03-03 15:09:07 1851 2

原创 2021-03-02低轨卫星定轨流程整理学习02

时间系统及各时间系统间转换代码、软件、网站整理1 时间系统及转换①国际原子时（TAI）原子时是基于原子的能量跃迁产生的电磁振荡定义的时间。国际原子时（TAI）的时间单位定义为国际单位秒（SI），规定 1958 年 1 月 1 日世界时（UT）0 时为时间起算点，使TAI与UT符合，国际原子时是由多台世界天文台的原子钟来共同维持的，是一个连续并均匀变化的时间系统。②恒星时（ST）恒星时是通过春分点在当地的时角来定义的，通常使用格林尼治恒星时，格林尼治恒星时分为格林尼治平恒星时（GMST）和格林尼治真

2021-03-02 14:29:50 1213 1

原创 2021-03-01低轨卫星定轨流程整理学习01

2021-03-01低轨卫星定轨流程整理学习01（简要介绍）一、什么是定轨？低轨卫星精密定轨数据处理是通过低轨卫星所能够提供的观测资料以及已知的低轨卫星运动学规律通过统计的方法来确定低轨卫星在指定时空基准下的位置和速度。#二、定轨方法假如已知低轨卫星在某一历元t0 的状态向量X0，并给定低轨卫星的运动微分方程，则利用数值积分的方法可以求得低轨卫星在t0后任一历元 t 的状态向量。然而在低轨卫星实际运动过程中，其初始状态向量X0无法提前知，并且运动微分方程中的物理常量和动力学模型以及观测模型不可避免地

2021-03-01 16:33:18 2789

qq_35519270的博客