排版设计模板网站,wordpress如何打开,河北省住房和城身建设厅网站,西安网站制作定制本文主要记录对RTKLIB源码中postpos、execses_b、execses_r 函数的源码解读#xff0c;不涉及其中的天线、星历等文件读取的内容#xff0c;且为个人理解#xff0c;如果有误#xff0c;欢迎交流讨论。 一、postpos 函数部分
/rxn2rtkp函数 → postpos函数传递参数#x… 本文主要记录对RTKLIB源码中postpos、execses_b、execses_r 函数的源码解读不涉及其中的天线、星历等文件读取的内容且为个人理解如果有误欢迎交流讨论。 一、postpos 函数部分
/rxn2rtkp函数 → postpos函数传递参数 postpos函数接收参数 各参数含义
• gtime_t ts: 处理开始时间。ts.time 0 表示没有时间限制。
• gtime_t te: 处理结束时间。te.time 0 表示没有时间限制。
• double ti: 处理时间间隔秒。ti 0 表示处理整个时间段。
• double tu: 处理段时间秒。tu 0 表示处理所有时间。
• const prcopt_t *popt: 处理选项包含数据处理的各种参数配置。
• const solopt_t *sopt: 解算选项定义输出格式等参数。
• const filopt_t *fopt: 文件选项包含输入输出文件的路径等。
• char **infile: 输入文件列表包含观测文件、导航文件等。
• int n: 输入文件数量。
• char *outfile: 输出文件路径可为空字符串表示输出到标准输出。
• const char *rov: 移动站rover的ID列表。(用空格分隔)
• const char *base: 基准站base station的ID列表。(用空格分隔)
• 备注: 输入文件应包含观测数据、导航数据、精密星历/钟差数据可选、SBAS 日志文件可选、SSR
• 消息日志文件可选和 TEC 网格文件可选。只有第一个输入文件中的观测数据被视为流动站数据。
• 输入文件的类型通过文件扩展名识别如下
• .sp3,.SP3,.eph*,.EPH*: 精密星历 (sp3c)
• .sbs,.SBS,.ems,.EMS : SBAS 消息日志文件 (rtklib 或 ems)
• .rtcm3,.RTCM3 : SSR 消息日志文件 (rtcm3)
• .*i,.*I : TEC 网格文件 (ionex)
• 其他 : RINEX 观测、导航、GNAV、HNAV、QNAV 或钟差文件
• 输入文件可以包含通配符*:如果文件包含通配符将会使用展开后的多个文件。
• 输入文件可以包含关键词:如果文件包含关键词关键词将会被日期、时间、流动站 ID 和基准站 ID 替换进
• 行多次会话分析。请参阅 reppath() 了解关键词。
• 输出文件同样可以包含关键词如果输出文件不包含关键词所有多次会话分析的结果将输出到一个单一的输
• 出文件中。
• SSR 改正仅对前向估计有效。第一个 if 语句如果开始时间、结束时间 ≠0处理单元时间 0 reppath函数将文件路径中的关键字替换为 - 日期、时间、基准站流动站ID * args : char *path I 文件路径
* char *rpath O 存储替换关键字后的文件路径
* gtime_t time I time (gpst) (time.time0: not replaced)
* char *rov I 流动站ID (: not replaced)
* char *base I 基准站ID (: not replaced) 1. extern int reppath(const char *path, char *rpath, gtime_t time, const char *rov,
2. const char *base)
3. { //初始化参数
4. double ep[6],ep0[6]{2000,1,1,0,0,0}; //用于存储时间的各部分年、月、日、小时、分钟、秒
5. int week,dow,doy,stat0; //GPS周、周内天0-6、年积日一年中的第几天、记录替换状态
6. char rep[64];
7.
8. strcpy(rpath,path); //将path主函数传参过来的infile或outfile中的路径复制到rpathpostpos函数内部定义的用于后续传参的ifile或ofile 9. if (!strstr(rpath,%)) return 0; //检查路径中是否有%关键字没有直接返回0
10. if (*rov ) stat|repstr(rpath,%r,rov ); //检查rov、base是否非空替换路径中%r、%b关键字
11. if (*base) stat|repstr(rpath,%b,base); //rov、base为流动站和基准站的ID列表
12. // | 位或操作符其实就是判断返回值是否为非0就成功
13.
14. if (time.time!0) { //如果传入时间不为0执行后续所有代码
15. time2epoch(time,ep); //将传入gps时间转换为历元时间年月日时分秒并存储在 ep 数组中
16. ep0[0]ep[0]; //更新初始参数时间为当前时间的年份
17. dow(int)floor(time2gpst(time,week)/86400.0); //计算 GPS 周内天数
0表示星期天1 表示星期一以此类推
18. doy(int)floor(timediff(time,epoch2time(ep0))/86400.0)1; //计算当年第几天 //根据路径中不同关键字对应格式进行替换这里要注意传入的时间代表什么
19. sprintf(rep,%02d, ((int)ep[3]/3)*3); stat|repstr(rpath,%ha,rep);
20. sprintf(rep,%02d, ((int)ep[3]/6)*6); stat|repstr(rpath,%hb,rep);
21. sprintf(rep,%02d, ((int)ep[3]/12)*12); stat|repstr(rpath,%hc,rep);
22. sprintf(rep,%04.0f,ep[0]); stat|repstr(rpath,%Y,rep);
23. sprintf(rep,%02.0f,fmod(ep[0],100.0)); stat|repstr(rpath,%y,rep);
24. sprintf(rep,%02.0f,ep[1]); stat|repstr(rpath,%m,rep);
25. sprintf(rep,%02.0f,ep[2]); stat|repstr(rpath,%d,rep);
26. sprintf(rep,%02.0f,ep[3]); stat|repstr(rpath,%h,rep);
27. sprintf(rep,%02.0f,ep[4]); stat|repstr(rpath,%M,rep);
28. sprintf(rep,%02.0f,floor(ep[5])); stat|repstr(rpath,%S,rep);
29. sprintf(rep,%03d, doy); stat|repstr(rpath,%n,rep);
30. sprintf(rep,%04d, week); stat|repstr(rpath,%W,rep);
31. sprintf(rep,%d, dow); stat|repstr(rpath,%D,rep);
32. sprintf(rep,%c, a(int)ep[3]); stat|repstr(rpath,%H,rep);
33. sprintf(rep,%02d, ((int)ep[4]/15)*15); stat|repstr(rpath,%t,rep);
34. }
35. //如果路径中存在以下时间关键字但没有给出有效时间则返回-1 2
36. else if (strstr(rpath,%ha)||strstr(rpath,%hb)||strstr(rpath,%hc)||
37. strstr(rpath,%Y )||strstr(rpath,%y )||strstr(rpath,%m )||
38. strstr(rpath,%d )||strstr(rpath,%h )||strstr(rpath,%M )||
39. strstr(rpath,%S )||strstr(rpath,%n )||strstr(rpath,%W )||
40. strstr(rpath,%D )||strstr(rpath,%H )||strstr(rpath,%t )) {
41. return -1; /* no valid time */
42. }
43. return stat; //这个返回状态是只要有路径中关键字被替换就会返回非0值
44. }
* notes : 路径中的以下关键字将被日期、时间和名称替换
* %Y - yyyy : year (4 digits) (1900-2099)
* %y - yy : year (2 digits) (00-99)
* %m - mm : month (01-12)
* %d - dd : day of month (01-31)
* %h - hh : hours (00-23)
* %M - mm : minutes (00-59)
* %S - ss : seconds (00-59)
* %n - ddd : day of year (001-366)
* %W - wwww : gps week (0001-9999)
* %D - d : day of gps week (0-6)
* %H - h : hour code (a0,b1,c2,...,x23)
* %ha- hh : 3 hours (00,03,06,...,21)
* %hb- hh : 6 hours (00,06,12,18)
* %hc- hh : 12 hours (00,12)
* %t - mm : 15 minutes (00,15,30,45)
* %r - rrrr : rover id
* %b - bbbb : base station id 二、 execses_b 函数
/postpos函数 → execses_b函数传递参数
在postpos函数中根据是否输入了开始时间、结束时间和段时间对输入文件进行处理但是不管如何处理将输入文件分成多少个传递到exceses_b函数都是对某段时间进行处理只需要知道开始时间、结束时间和采样频率就可以。
该函数其实就是
if (ts.time!0te.time!0tu0.0) else if (ts.time!0) else execses_b函数接收参数 三、 execses_r 函数 四、总体分析
