关于从GNSS信号中相位平滑伪距提取STEC中的DCB

m0_65087756

于 2024-08-11 16:48:54 发布

阅读量1.1k

点赞数 45

文章标签：算法

本文链接：https://blog.csdn.net/m0_65087756/article/details/141106788

版权

最近在从GNSS信号中提取STEC，在这步之前的工作先略过。但是这个DCB处理的公式很迷惑。找了好多文献才找到靠谱的解释。

DCB可以从IGS官方发布的文件中获取，问题在于，这个文件并不包括全球所有测站，所以，自己解算DCB是比较普遍的做法。

先来看公式：：

$\begin{aligned}&P_{r,f}^{s}=\rho_{r}^{s}+c(\delta t_{r}-\delta t^{s})+T_{r}^{s}+I_{r,f}^{s}+B_{r,f}+B_{f}^{s}+\varepsilon_{r,f}^{s}\\&L_{r,f}^{s}=\rho_{r}^{s}+c(\delta t_{r}-\delta t^{s})+T_{r}^{s}-I_{r,f}^{s}+b_{r,f}+b_{f}^{s}+\lambda_{f}N_{r,f}^{s}+\xi_{r,f}^{s}\end{aligned}$

公式就不多解释了，r,s,f分别代表接收机，卫星，信号频率，P是伪距，L是载波相位测量。

B是伪距硬件延迟，b是相位硬件延迟。

通过差分f1,f2信号观测数据，也就是无几何距离组合，可以消除频率无关项，提取电离层。

也就是 P1-P2， L1-L2，这两种方法都能提取出STEC。

$\begin{aligned}&P_{r,4}^{s}=P_{r,f_{1}}^{s}-P_{r,f_{2}}^{s}=I_{r,f_{1}}^{s}-I_{r,f_{2}}^{s}+B_{r,f_{1}}+B_{f_{1}}^{s}-(B_{r,f_{2}}+B_{f_{2}}^{s})+\varepsilon_{r,4}^{s}\\&L_{r,4}^{s}=L_{r,f_{1}}^{s}-L_{r,f_{2}}^{s}=-I_{r,f_{1}}^{s}+I_{r,f_{2}}^{s}+b_{r,f_{1}}+b_{f_{1}}^{s}-(b_{r,f_{2}}+b_{f_{2}}^{s})+N_{r,4}^{s}+\xi_{r,4}^{s}\end{aligned}$

其中：

$\begin{aligned}&N_{r,4}^{s}=\lambda_{f_{1}}N_{r,f_{1}}^{s}-\lambda_{f_{2}}N_{r,f_{2}}^{s}\\&\varepsilon_{r,4}^{s}=\varepsilon_{r,f_{1}}^{s}-\varepsilon_{r,f_{2}}^{s}\\&\xi_{r,4}^{s}=\xi_{r,f_{1}}^{s}-\xi_{r,f_{2}}^{s}\end{aligned}$

$\begin{aligned}&SATEC_{r}^{s}=\alpha P_{r,4}^{s}=T_{S,r}^{s}+DCB_{r,f_{1}f_{2}}+DCB_{f_{1}f_{2}}^{s}+\alpha\varepsilon_{r,4}^{s}\\&SRTEC_{r}^{s}=-\alpha L_{r,4}^{s}=T_{S,r}^{s}-DPB_{r,f_{1}f_{2}}-DPB_{f_{1}f_{2}}^{s}-\alpha N_{r,4}^{s}-\alpha\xi_{r,4}^{s}\end{aligned}$

SATEC为斜向绝对电子密度总数，由伪距差分提取得到的TEC观测数据；
SRTEC为斜向相对电子密度总数，由相位差分而来的TEC观测数据

另外，式中其它项表示：
$\begin{cases}\alpha=\frac{f_1^2f_2^2}{40.3(f_2^2-f_1^2)}\\DCB_{r,f_1f_2}=\alpha(B_{r,f_1}-B_{r,f_2})\\DCB_{f_1f_2}^s=\alpha(B_{f_1}^s-B_{f_2}^s)\\DPB_{r,f_1f_2}=\alpha(b_{r,f_1}-b_{r,f_2})\\DPB_{f_1f_2}^s=\alpha(b_{f_1}^s-b_{f_2}^s)\end{cases}$

注意，这里DCB、DPB以TECu为单位，除以(a*c)即可转换至时间单位，其中 c为光速。从上式可以看到，伪距及相位经过双频差分后均可以提取出包含误差的TEC数据

我的疑惑在于，在大多数的文章中，DPB这一项（有时用b，或者d表示），在推导中，推着推着就没了，但是，这一项，是合并进了DCB还是消除了，很多文章并没有很明确的解释。

继续：

相位差分提取TEC:

对于相位差分提取的TEC观测数据即SRTEC，误差包括相位对应的硬件延迟DPB（包含卫星以及接收机的DPB）、差分后的整周模糊度对应项, $\alpha N_{r,4}^s$ 以及噪声及多径造成的误差。其中DPB比较稳定， $\alpha N_{r,4}^s$ 在相位观测值未发生周跳的情况下即一个连续的观测弧段下，被认为是固定值。令
$dpb_{r,f_1f_2}^s=-DPB_{r,f_1f_2}-DPB_{f_1f_2}^s-\alpha N_{r,4}^s ,$ ，则 $dpb_{r,f_1f_2}^s$ 在每一个连续观测弧段对应一个固定值，将其改写为 $dpb_{arc}$ 下标arc表示观测弧段编号.