主要是由postpos.c中的这段代码实现。
这里的obs为由inputobs()函数提取出的某一个历元的所有观测值,nobs为这一历元的观测值数量。在读取观测值和广播星历后就按照卫星的PRN号和时间进行了重新排列与重复数据的删除,所以无论renix文件中观测值数据是如何排列的,在rtklib中都变成了按照每一历元,卫星PRN号从小到大的原则排列。
stasys()函数的作用是根据卫星的PRN号范围卫星所属的系统,注意这里卫星系统在rtklib.h中采用宏定义定义成了16进制整数,因此可以用位运算符(&)很方便的判断一个卫星是否属于所选系统。
这里的for循环遍历了当前历元的每一个数据,将需要用到的卫星数据往前放,同时n记录了需要用到的数据个数。这样对obs重排后取前n个数就是需要用到的观测值数据。
/* 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;