模型，并根据所述步骤13中拟合 得到的各所述GNSS测站对应的所述周期常数项心、年周期项的振幅与初始相位&和di、以 及半年周期项的振幅与初始相位A2、d2，通过反距离加权内插方法计算得到所述经炜度格网 函数模型中的每个格网点对应的周期常数项&、年周期项的振幅与初始相位&和di、以及 半年周期项的振幅与初始相位八2和d2。
2. 根据权利要求1所述的天顶对流层延迟建模方法，其特征在于，所述经炜度格网函 数模型在经炜方向的分辨率为rxr至5°X5°。
3. -种天顶对流层延迟建模装置，其特征在于，包括： 接收模块，其用于接收分布在中国大陆区域的陆态网GNSS测站分别实测的对流层天 顶延迟量； 连接至所述接收模块的高程改正模型建立模块，其用于计算各所述GNSS测站对应的 对流层天顶延迟量的单日值ZTDh(doy)和年均值ZTDh，并根据所述年均值ZTDh和式（1) 拟合得到各所述GNSS测站对应的高程改正系数0和归化到地球椭球面上的年均延迟量 ZTD。，然后根据各GNSS测站对应的所述高程改正系数0，由式（1'）将各GNSS测站对应的 单日值ZTDh(doy)改正到地球椭球面上，以获得高程改正后的单日延迟量ZTDjdoy): ZTDh=ZTDoXe^ (1)， ZTDh(doy) =ZTD〇(doy)XePh (1，） 在式（1)中，e为自然常数，h表示各所述GNSS测站的实际高程，doy表示年积日，其 中，0随着炜度范围不同而不同； 连接至所述接收模块和所述高程改正模型建立模块的周期改正模型建立模块，其根据 各所述GNSS测站对应的所述高程改正后的单日延迟量ZTDjdoy)和式（2)拟合得到各所 述GNSS测站对应的周期常数项\、年周期项的振幅与初始相位心和di、以及半年周期项的 振幅与初始相位A2、d2:
以及 连接至所述周期改正模型建立模块的格网函数模型建立模块，其用于建立基于中国大 陆区域的经炜度格网函数模型，并根据所述步骤13中拟合出的各所述GNSS测站对应的所 述周期常数项心、年周期项的振幅与初始相位&和di、以及半年周期项的振幅与初始相位 A2、d2，通过反距离加权内插方法计算得到所述经炜度格网函数模型中的每个格网点对应的 所述周期常数项心、年周期项的振幅与初始相位&和di、以及半年周期项的振幅与初始相 位八2和d2。
4. 根据权利要求3所述的天顶对流层延迟建模方法，其特征在于，所述经炜度格网函 数模型在经炜方向的分辨率为rxr至5°X5°。
5. -种天顶对流层延迟测量方法，以基于权利要求3所述的天顶对流层延迟建模装置 测量用户站在预定年积日doy对应的天顶对流层延迟量，其特征在于，该方法包括以下步 骤： 步骤S21，根据所述用户站所在位置的经炜度，查询所述用户站在所述经炜度格网函数 模型中对应的格网，并获取该格网的四个格网点的位置及其对应的所述周期常数项心、年 周期项的振幅与初始相位4和di、以及半年周期项的振幅与初始相位^和d2; 步骤S22,根据所述用户站所在格网的所述四个格网点对应的所述周期常数项&、年周 期项的振幅与初始相位4和di、以及半年周期项的振幅与初始相位^和d2,通过反距离加 权内插方法计算所述用户站对应的所述周期常数项&、所述年周期项的振幅与初始相位Ai 和屯、以及所述半年周期项的振幅与初始相位八2和d2; 步骤S23,将所述用户站对应的所述周期常数项\、年周期项的振幅与初始相位心和dl、和半年周期项的振幅与初始相位^和d2引入式（2)，以计算所述用户站在所述年积日 doy对应的归化到地球椭球面上的对流层天顶延迟量，记为ZTDd(doy); 步骤S24,根据所述用户站的炜度查询其对应的所述高程改正系数0，并根据所述式 (1'）计算所述用户站在所述年积日doy对应的实际对流层天顶延迟量ZTDh (doy)。
6. -种天顶对流层延迟测量装置，以基于权利要求5所述的天顶对流层延迟测量方法 测量用户站在预定年积日对应的天顶对流层延迟量，其特征在于，包括： 连接至的所述天顶对流层延迟建模装置的格网查询模块，其根据所述用户站所在位置 的经炜度，查询所述用户站在所述经炜度格网函数模型中对应的格网，并获取该格网的四 个格网点的位置及其对应的所述周期常数项心、年周期项的振幅与初始相位&和di、以及 半年周期项的振幅与初始相位八2和d2; 连接至的所述格网查询模块的用户站参数获取模块，其根据所述用户站所在格网的所 述四个格网点对应的所述周期常数项心、年周期项的振幅与初始相位&和di、以及半年周 期项的振幅与初始相位4和(12，通过反距离加权内插方法计算所述用户站对应的所述周期 常数项心、年周期项的振幅与初始相位^和七、以及半年周期项的振幅与初始相位4和(12; 连接至的所述用户站参数获取模块的周期修正模块，其将所述用户站对应的所述周期 常数项Ac!、年周期项的振幅与初始相位AJPdi、和半年周期项的振幅与初始相位八2和d2弓丨 入式（2)，以计算所述用户站在所述年积日对应的归化到地球椭球面上的对流层天顶延迟 量，记为ZTDQ(doy); 连接至的所述周期修正模块的高程修正模块，其根据所述用户站的炜度查询其对应的 所述高程改正系数0，并根据所述式（1'）计算所述用户站在所述年积日对应的实际对流 层天顶延迟量ZTDh(doy)。
【专利摘要】本发明提供一种天顶对流层延迟建模方法，包括：S11，获取各测站实测的对流层天顶延迟量；S12，计算各测站对应的对流层天顶延迟量的年均值ZTDh，并根据ZTDh拟合得到各测站对应的高程改正系数β和归化到椭球面上的年均延迟量ZTD0，再根据β获取各测站对应的高程改正后的单日延迟量ZTD0(doy)；S13，根据ZTD0(doy)拟合得到各测站对应的周期常数项A0、年周期项的振幅与初始相位A1和d1及半年周期项的振幅与初始相位A2、d2；S14，建立基于中国大陆区域的经纬度格网函数模型，并根据各测站对应的参数A0、A1、d1、A2和d2计算得到格网函数模型中的各格网点对应的A0、A1、d1、A2、d2。本发明适用于中国大陆区域，精度明显提高，可提升中国大陆区域卫星导航定位系统的整体服务水平。
【IPC分类】G01S19-07
【公开号】CN104777488
【申请号】CN201510112093
【发明人】王君刚, 陈俊平, 王解先, 杨赛男
【申请人】中国科学院上海天文台, 同济大学
【公开日】2015年7月15日
【申请日】2015年3月13日