【技术实现步骤摘要】
一种基于去趋势互相关分析的GNSS影像生成方法
本专利技术涉及GNSS数据精密处理
，尤其涉及一种基于去趋势互相关分析的GNSS影像生成方法。
技术介绍
近年来，国内外建立了各种GNSS监测网络，测站数目日益增多，极大地扩展了GNSS测站的覆盖范围，产生了大量的观测数据，为监测地壳形变提供了基础；因此，利用GNSS坐标时间序列生成地壳形变影像已成为现实(Hammondetal.,2016)。然而，地壳形变在不同空间尺度下呈现出既有广泛分布，又有局部突变的空间分布特征；若不能有效地描述地壳形变的空间特征，会严重影响地壳形变影像结果的可靠性与地学解释。以往，测站之间的相关性主要利用任意两个测站组成的测站对的残差时间序列进行计算，这样，一方面会引入由于速度模型建模导致的误差，另一方面，会大大减少原始时间序列的有效信息。本专利技术提供了一种采用原始时间序列，基于去趋势互相关分析描述测站间速度相关性，利用临近测站的速度对研究区域进行空间插值，最终生成地壳形变影像(即GNSS影像)的方法。
技术实现思路
本专利...
【技术保护点】
1.一种基于去趋势互相关分析的GNSS影像生成方法，其特征在于，该方法包括以下步骤：/nS1、获取GNSS测站坐标时间序列的观测值，并获取每个GNSS测站的坐标(B
【技术特征摘要】
1.一种基于去趋势互相关分析的GNSS影像生成方法，其特征在于，该方法包括以下步骤：
S1、获取GNSS测站坐标时间序列的观测值，并获取每个GNSS测站的坐标(Bi，Li)，其中，i＝1，，…，n，n为总测站数；Bi，Li为第i个测站在大地坐标系下的纬度与经度；
S2、提取GNSS测站坐标时间序列观测值中共同跨度的部分，记为其中，mi为第i个测站共同跨度部分的历元个数；
S3、计算GNSS测站坐标时间序列共同跨度部分的速度及其不确定性，速度记为v1，…，vn，不确定性记为un1，…，unn；
S4、选择第i个GNSS测站作为参考站，第j个GNSS测站作为对比站，对Yi和Yj进行去趋势互相关分析，计算去趋势互相关系数DCCACCij，其中，i＝1，…，n；j＝1，…，n；i≠j；
S5、逐一对所有GNSS测站进行空间滤波，记滤波后速度为v′1，…，v′n；
S6、将研究区域格网化，逐一对格网点进行空间插值，最终生成GNSS影像。
2.根据权利要求1所述的基于去趋势互相关分析的GNSS影像生成方法，其特征在于，步骤S4中计算去趋势互相关系数的具体方法为：
S41、获得原始序列Yi和Yj的累加序列Yi′和Yj′：
其中，mean(*)表示取均值；
S42、将Yi′和Yj′分别划分为L个相互重叠的数据窗口，每个窗口内有w个历元，其中：
其中，min(*)表示取最小值，w为窗口宽度；
S43、基于最小二乘准则，建立各数据窗口的一次多项式模型，形成累加序列Yi′和Yj′的函数模型，记为LSi，l和LSj，l，l＝1，2，…，min(mi，mj)；
S44、计算窗口为w1时的去趋势互相关系数：
其中：
S45、改变窗口宽度，计算最终的去趋势互相关系数：
DCCACCij＝median(ρDCCA(w1)，…，ρDCCA(wN))；
其中，median(*)表示取中位数，w1，…，wN表示不同的窗口宽度，窗口宽度的取值范围为用户自定。
3.根据权利要求1所述的基于去趋势互相关分析的GNSS影像生成方法，其特征在于，步骤S5中对GNSS测站进行空间滤波的具体方法为：
S51、基于测站坐标构建狄洛尼三角网，选取与滤波测站相连的测站作为备选测站；
S52、去除其中与滤波测站去趋势互相关系数为负的测站，记滤波测站速度为ve，剩余的备选测站速...
【专利技术属性】
技术研发人员：周晓慧，杨艺林，范文蓝，
申请(专利权)人：武汉大学，
类型：发明
国别省市：湖北;42
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。