关键词:
GNSS-R
原始中频数据
高分辨率DDM
河流流速检测
湖泊边界检测
摘要:
地球反射全球导航卫星系统(Global Navigation Satellite System,GNSS)信号可由地基、机载或专用轨道星载接收机接收,用于遥感和地球观测目的。近年来,GNSS-R(Global navigation satellite system reflectometry technology)实验任务大多只能提供GNSS-R Level-1(L1)数据即信号反射功率的延迟多普勒图(Delay Doppler Maps,DDM)作为主要数据产品,由于其时空分辨率较低且固定,已不能满足目前高精度的反演需求。除了DDM产品外,其中一些GNSS-R任务实验还收集了大量的原始中频(Intermediate Frequency,IF)数据,原始中频数据是指在模数转换器(Analog to Digital Converter,ADC)之后和任何数字处理之前记录的原始卫星信号样本的数据流。这些原始中频数据的未处理性质为探索GNSS反射测量(GNSS-R)技术在先进地球物理应用提供了独特的机会。相较于GNSS-R L1 DDM产品数据,使用软件定义的接收机对GNSS-R的原始中频数据进行处理,可以提取出更多的潜在信息(如干涉载波相位观测值、残差多普勒观测值等),挖掘更多的反射面物理特性信息,获得更高分辨率的跟踪相关结果,反演结果具有更高的精度和更高的时空分辨率。本研究以GNSS-R原始中频数据为基础进一步获取高时空分辨率、高精度的信号特征量进而实现GNSS-R高精度反演,以目前本课题组实施的岸基GNSS-R河流流速项目实验数据和公开的星载GNSS-R原始中频数据为依托,以岸基平台下河流流速检测与星载平台下湖泊检测为科学目标,旨在提取GNSS-R反射信号中的物理反演特征、提升GNSS-R反演精度、挖掘更多GNSS-R应用可行性。本文主要研究工作和贡献如下:(1)给出了岸基实验条件下GNSS直射信号和反射信号捕获跟踪方法以及直射/反射信号之间码片延迟和载波相位观测值延迟的精确提取算法。(2)针对GNSS-R星载遥感卫星原始反射信号高动态性、低信噪比的特点,将二阶锁频环辅助三阶锁相环的闭环跟踪方法运用于GNSS-R原始反射中频数据的处理上,成功对星载实验条件下(以CYGNSS(cyclone global navigation satellite system),TDS-1(Tech Demo Sat-1)为例)的原始中频数据实现了跟踪。对于CYGNSS,分别给出了使用CYGNSS原始中频数据对海洋和陆地表面的反射事件进行相关后生成DDM的结果,使用皮尔逊相关系数定量地对CYGNSS原始中频数据的相关结果与对应的L1 DDM进行验证。其中,海洋反射事件的皮尔逊相关系数为0.976,陆地反射事件的皮尔逊相关系数为0.807,均表现出高相关性。对TDS-1原始数据相关结果分析,分别给出了其海面反射事件和海冰反射事件的相关结果,定性分析相关结果符合菲涅尔反射原理。CYGNSS和TDS-1的结果证明了使用本文设计的跟踪环路,实施对星载实验条件下的GNSS-R原始中频数据解调的泛用性,使用该方法可以生成更高分辨率的DDM,为后续星载实验条件下的高分辨率反演奠定了基础。(3)在实现对岸基平台下GNSS-R原始中频数据解调与跟踪的基础上,本研究使用基于载波相位观测值的河流流速反演模型。干涉相位是多普勒频率的积分。通过第一章介绍的闭环跟踪方法对GNSS-R原始中频数据集进行跟踪处理,得到河流流动产生的多普勒频移观测值,进而实现速度反演。实验数据为本课题组于2021年4月22日,在江苏省南京市大胜关长江南岸进行了近两个小时的岸基河流实验数据。在使用GPS L1频段GNSS-R原始中频数据实现河流流速反演后,进一步实现了BDS B1I GEO卫星与IGSO卫星的联合反演,验证了利用BDS反射信号反演河流流速的可行性。与实际河流流速相比,GPS L1 PRN4(第一阶段)反演精度可达0.028 m/s(平均绝对误差,MAE)和0.036 m/s(均方根误差,RMSE)。BDS GEO 2反演精度可达0.048 m/s(MAE)和0.063 m/s(RMSE),BDS IGSO 10反演精度可达0.061 m/s(MAE)和0.073 m/s(RMSE)。进一步分析实验环境表明,卫星高度角变化率和镜面点与流量表的距离对河流流速反演精度有一定的负面影响。镜面点被障碍物阻挡或距离测速仪太远,都可能导致MAE和RMSE的不确定误差。无论是卫星高度角还是信号强度都与反演精度没有明显的相关性,这与理论原理是一致的。(4)在实现对星载平台下CYGNSS、TDS-1原始中频数据解调与跟踪的基础上,本文首次使用CYGNSS原始中频数据,使用第一章对星载平台下GNSS-