RTKLIB学习--前向滤波

#前言
如果要详细了解RTKLIB或进行二次开发，了解obs指针所存储每个历元的卫星观测数据是必不可少的环节，此文对RTKLIB的（由于后处理和实时运行都要用到前向滤波）前向滤波（从文件头读取观测数据到obs结构体中，进行每个历元的解算）进行详细的分析。结合代码调试来学习效果更佳

procpos函数

procpos是开启历元循环所在的函数，也是循环读取每个历元观测数据所在的函数。

solstatic=sopt->solstatic&&//静态解(popt->mode==PMODE_STATIC||popt->mode==PMODE_PPP_STATIC);rtkinit(&rtk,popt);rtcm_path[0]='\0';//实时数据流

while ((nobs=inputobs(obs,rtk.sol.stat,popt))>=0) {

在这个函数中循环读取一个历元的观测数据，而广播星历则是一次读取完。

/* exclude satellites 排除没有参与定位的卫星*/for (i=n=0;i<nobs;i++) {if ((satsys(obs[i].sat,NULL)&popt->navsys)&&popt->exsats[obs[i].sat-1]!=1) obs[n++]=obs[i];}if (n<=0) continue;//n为一个历元有效观测卫星数

if (!rtkpos(&rtk,obs,n,&navs)) continue;

正式进入定位模块，后面的条件判断则是一个历元定位结束后根据滤波方向决定下一个读取的历元。

inputobs()

在inputobs中通过不同的滤波方式读取基准站，流动站观测数据；由于后向滤波(从文件末尾读)只适用于后处理，在这里只介绍前向滤波(从文件头读取数据)。

其中obss.data[] 存储了所有历元的卫星观测数据，在inputobs(）函数中通过流动站索引 iobsu，基准站索引 iobsr 从 obss.data 中获取 obsd_t *obs 其中存储了一个历元的所有卫星的观测数据，由于在spp,ppp,模式中不涉及到基准站，inputobs()函数对基准站的处理不易读懂，算法很巧妙，需要下一定的功夫理通

* input obs data, navigation messages and sbas correction -------------------*/
static int inputobs(obsd_t *obs, int solq, const prcopt_t *popt)
{gtime_t time={0};char path[1024];int i,nu,nr,n=0;trace(3,"infunc  : revs=%d iobsu=%d iobsr=%d isbs=%d\n",revs,iobsu,iobsr,isbs);if (0<=iobsu&&iobsu<obss.n) {settime((time=obss.data[iobsu].time));if (checkbrk("processing : %s Q=%d",time_str(time,0),solq)) {aborts=1; showmsg("aborted"); return -1;}}if (!revs) { /* input forward data 从第一个历元读取*/if ((nu=nextobsf(&obss,&iobsu,1))<=0) return -1;//获取流动站观测卫星数/,&iobsu是上一个历元的基准站观测值索引，需要更新到当前历元的流动站观测值索引if (popt->intpref) {//内插参考观测值时，将iobsr索引(上一个历元的基准站索引)定位到当前历元的流动站索引for (;(nr=nextobsf(&obss,&iobsr,2))>0;iobsr+=nr)//最终iobsr被定位到当前历元流动站索引if (timediff(obss.data[iobsr].time,obss.data[iobsu].time)>-DTTOL) break;//iobsr所在历元时间段与当前历元的流动站历元时间差大于阈值时，退出}else {//不内插参考观测值时for (i=iobsr;(nr=nextobsf(&obss,&i,2))>0;iobsr=i,i+=nr)//获取上一个历元基准站观测数；iobsr被定位到当前历元基准站索引if (timediff(obss.data[i].time,obss.data[iobsu].time)>DTTOL) break;//上一个历元基准站索引时间与当前历元流动站索引时间差大于阈值时}

注：

①当 popt->intpref 为1时，表明参考站观测值是内嵌到流动站观测中的，也就是说一个历元观测值为(流动站观测值+基准站观测值)；

​ &iobsr 初始为上一个历元的基准站索引，nr=上一历元基准站观测数，（假设历元时间段为30s）此时timediff= -30<阈值；iobsr+=nr成为了当前历元的流动站索引值，nr=当前历元的流动站数，timediff=0;结束循环

②popt->intpref 为0时，流动站观测值和基准站观测值是分开的。

​ &iobsr 初始为上一个历元的基准站索引，nr=上一历元基准站观测数，（假设历元时间段为30s）此时timediff= -30<阈值；iobsr=i不变，i +=nr,成为了当前历元基准站索引，timediff=0<阈值; iobsr=i 成为了当前历元基准站索引，i +=nr成为了下一个历元的基准站索引，timediff=30>阈值，结束循环

nr=nextobsf(&obss,&iobsr,2);//确保iobsr被定位到当前历元基准站索引，获取当前历元基准站观测数

如果 popt->intpref 为1时 iobsr从当前历元的流动站索引定位到基准站索引，然后获取观测卫星数，popt->intpref 为0时只获取观测卫星数。

如果是单点定位模式，则在第一次调用nextobsf() 时，由于观测数据的rcv无法匹配到给定参数为“2"的rcv,&iobsr 被遍历到最后一个历元，在之后调用nextobsf() 时，不会进入到if语句 的内部

for (i=0;i<nu&&n<MAXOBS*2;i++) obs[n++]=obss.data[iobsu+i];//复制流动站观测卫星数据
for (i=0;i<nr&&n<MAXOBS*2;i++) obs[n++]=obss.data[iobsr+i];//复制基准站观测卫星数据
iobsu+=nu;//流动站索引此时变为参考站索引 或下一个历元的流动站索引(spp,ppp)

nextobsf()

search next observation data index 前向滤波的方式--------------------------*/
static int nextobsf(const obs_t *obs, int *i, int rcv)
{double tt;int n;/*如果rcv为流动站，i在上一级函数中通过+=nu成为了上一个历元的基准站索引；如果是rtk模式一个历元段包括(rcv=1:流动站段 + rcv=2:基准站段)(参考观测值内插时)*/for (;*i<obs->n;(*i)++) if (obs->data[*i].rcv==rcv) break;//rtk模式:rcv=1(基准站索引->当前历元流动站站索引)；rcv=2(-->下一段索引，可能是基准站，也可能是流动站)for (n=0;*i+n<obs->n;n++) {//遍历当前段观测卫星，tt=timediff(obs->data[*i+n].time,obs->data[*i].time);//第n个观测卫星时间与索引观测卫星时间差if (obs->data[*i+n].rcv!=rcv||tt>DTTOL) break;//如果卫星rcv改变或者历元时差过大,结束循环}return n;//返回索引所在历元rcv观测卫星数