一、概念
DCB(Differential Code Bias 差分码偏差)是全球卫星导航系统(GNSS)中,通过不同信号得到的观测值之间存在的系统性偏差。DCB是由卫星和接收机硬件延迟的影响造成的。一般来说接收机端的DCB可以跟接收机钟差一起解算,被接收机钟差所吸收;而卫星端的DCB在精密定位中必须得到补偿。
DCB主要有两种
相同频率不同码之间存在的偏差(如P1-C1、P2-C2等);
不同频率之间存在的偏差(如P1-P2);
二、数据来源
目前DCB数据来源主要有两个组织:IGS 和 CODE。
IGS(International GNSS Service)组织是国际最权威的GNSS精密应用服务组织之一。随着GPS的现代化以及BDS和Galileo等导航系统的建设,IGS于2011年开始建立多模GNSS实验跟踪网(MGEX, Multi-GNSS Experiment),用于多模GNSS导航信号监测及相关技术研究。基于MGEX多系统观测数据,德国宇航中心(DLR)自2013年起开始向IGS组织提交包括BDS、GPS、GLONASS及Galileo在内的多系统DCB产品。自2015年起,武汉的中国科学院测量与地球物理研究所(IGG of CAS)成为全球第二家向IGS提交多系统DCB产品的机构。这些数据可以在 ftp://cddis.nasa.gov/gnss/products/bias/ 下载。瑞士伯尼尔大学的欧洲定轨中心(CODE:Centre for Orbit Determination in Europe)也提供DCB产品下载,其历史更加悠久。下面是其FTP下载地址:ftp://ftp.aiub.unibe.ch/CODE/ 。
三、应用实现
对于双频接收机,在电离层延时中我们得到经过电离层校正后的伪距为
RTKLIB 中的单点定位就是用这种方式进行计算的。可以参考如下代码:
gamma=SQR(lam[j])/SQR(lam[i]); /* f12/f22 */
P1=obs->P[i];
P2=obs->P[j];
P1_P2=nav->cbias[obs->sat-1][0];
P1_C1=nav->cbias[obs->sat-1][1];
P2_C2=nav->cbias[obs->sat-1][2];
/* if no P1-P2 DCB, use TGD instead /
if (P1_P20.0&&(sys&(SYS_GPS|SYS_GAL|SYS_QZS))) {
P1_P2=(1.0-gamma)*gettgd(obs->sat,nav);
}
if (opt->ionooptIONOOPT_IFLC) { / dual-frequency */
if (P1==0.0||P2==0.0) return 0.0;
if (obs->code[i]==CODE_L1C) P1+=P1_C1; /* C1->P1 */
if (obs->code[j]==CODE_L2C) P2+=P2_C2; /* C2->P2 */
/* iono-free combination */
PC=(gamma*P1-P2)/(gamma-1.0);
}
else { /* single-frequency */
if (P1==0.0) return 0.0;
if (obs->code[i]==CODE_L1C) P1+=P1_C1; /* C1->P1 */
PC=P1-P1_P2/(1.0-gamma);
}