集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法技术

本发明专利技术公开了一种集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其通过在测站点上，使天顶筒的光轴指向天顶方向对天顶区域进行拍摄，获取CCD恒星图像，并对其中恒星的图像坐标进行计算，通过GNSS测得的大地坐标和GNSS时间信号控制曝光历元，结合EGM2008重力场模型计算的测站垂线偏差，在适用星表中获取天顶区域恒星信息，实现星表中天球的天顶切平面区域恒星与CCD恒星图像中恒星的匹配识别，利用最小二乘法迭代计算天顶天文坐标，结合大地坐标可计算出观测点的垂线偏差。本发明专利技术相对于现有技术，具有操作简便、省时省力、以及测量精度高等特点，适于高精度垂线偏差的快速测量。

【技术实现步骤摘要】

本专利技术涉及一种垂线偏差的天文大地测量方法，尤其涉及一种集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法。
技术介绍
地面上一点的重力向量g和相应椭球面上的法线向量n之间的夹角定义为该点的垂线偏差。垂线偏差表征重力的方向，是大地测量归算、地球重力场和大地水准面模型精化、研究地球形状的必要地球物理量，在大地测量和空间技术等方面有着重要用途。垂线偏差是地球内部质量分布的表现，可用于反演大气反常折射研究，监测地球内部质量迁移和能量积累，监测地震、火山爆发等自然灾害，也可用于探测地下质量异常，进行资源勘探。在军事领域，垂线偏差是导弹、火箭发射等必需的基础数据。现有技术中，常用的垂线偏差观测方法，主要包括传统天文大地测量方法、重力测量方法、天文重力测量方法以及GPS水准测量方法。但是，上述这些方法普遍存在基础数据采集过程复杂，或者根本就无法有效获取、后期数据处理过程繁琐复杂、费时费力等缺点或不足：传统天文大地测量方法，通常采用T4、J05等天文经纬仪，配以一套守时、授时系统。其天文纬度的观测常用Telgte方法，天文经度的观测常用双星等高法并同时收录无线电时号。这种方法获得的垂线偏差精度可达到0.3\。但是，对于这种传统天文大地测量方法，通常，其为了使观测精度足够高，所使用的天文经纬仪的体积都很大、十分笨重；且仪器的操作对观测者经验要求过高、观测时间很长、效率极低。而且，其常用的守时、授时系统，不仅操作复杂，精度较低，且容易出错。重力测量方法，主要借助于大地水准面和地球椭球面上的重力异常，利用Stokes或VeningMeinesz公式计算垂线偏差。然而，这种方法的前提是假定大地水准面之外没有扰动物质，且其前提是，需要对全球重力异常数据已知，这十分困难、且根本就不现实。因此，该方法还没有得到独立的应用，并且该方法的垂线偏差的测量精度也偏低。天文重力测量方法，其实质是结合天文大地测量方法和重力测量方法来确定垂线偏差。该方法首先需要一定数量的天文大地点的垂线偏差值，然后利用重力测量方法加密重力异常数据，采用最小二乘配置或多面函数方法确定其他点的垂线偏差值。这种方法的主要缺点在于，求出的垂线偏差精度偏低。GPS水准测量方法，需要同时进行GPS测量和水准测量，获取基线长、大地方位角和高程异常之差，就可以求出垂线偏差。然而，这种方法的使用方面有条件限制，比如，对于某些地形起伏较大的区域进行精密水准测量的难度太大，基线也不能过长，且由于基线长和大地方位角的影响，导致其求出的垂线偏差精度偏低。
技术实现思路
本专利技术的目的是，提供一种具有普适特性，且高精度垂线偏差的数据采集过程省时省力、数据的处理方法相对简单的一种集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法。本专利技术为实现上述目的所采用的技术方案是，一种集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特征在于，包括以下步骤：第一步，在测站点上架设CCD天顶筒，使其光轴指向天顶方向，对天顶区域连续在0°和180°方向进行拍摄，获取CCD恒星图像，并由高精度电子水平仪获取天顶筒的倾斜数据，由GNSS获取观测时间信息和测站点的大地坐标；第二步，对CCD天顶筒拍摄的FITS格式的CCD恒星图像进行读取，再进行背景噪声消除；第三步，对噪声消除后的CCD恒星图像进行恒星星象区域的自动搜索；第四步，运用改进的修正二维矩方法，计算上述自动搜索出的恒星星象区域的能量中心的影像坐标；第五步，对适用星表进行处理，通过GNSS获取的时间信息以及测站点的大地坐标，结合EGM2008计算测站点粗略天文坐标。根据结果对星表中与CCD恒星图像相对应的区域进行截取，并计算曝光时刻天顶区域的恒星赤道坐标与视星等信息，并将恒星赤道坐标投影到切平面上转换为平面坐标；第六步，对CCD影像中的星象与星表投影到平面上的恒星进行识别、匹配；第七步，依据匹配成功的星对，利用最小二乘法计算摄影测量投影变换公式中的坐标转换系数，并迭代计算测站点的初始赤道坐标，再转换为测站点的初始天文坐标；第八步，将高精度电子水平仪获取的天顶筒倾斜改正数据和极移改正数据转换至子午方向和卯酉方向上，确定出测站点的精确天文坐标；第九步，结合GNSS测量数据解算出的测站点大地坐标，利用Helmert公式，即可计算出测站点的垂线偏差。上述技术方案直接带来的技术效果是，垂线偏差快速测量方法基于GNSS和CCD天顶筒，在满足测量精度和分辨率需要的基础上，利用改进的数据采集及处理方法，降低了垂线偏差获取过程所消耗的时间，提高了测量效率，能够满足高的测量精度和无地域限制等要求。优选为，上述背景噪声消除所采用的方法为3×3中值滤波法，具体步骤如下：步骤一，确定一个以某个像素为中心点的3×3邻域；步骤二，将该邻域中各个像素的灰度值(x1,x2,x3,...,xn)按从小到大进行排序：x1≤x2≤x3...≤xn；步骤三，按公式计算出各个像素的灰度值的中值Y：步骤四，以中值Y代替原来中心点的灰度值，进行中值滤波；步骤五，确定一个以另一个像素为中心点的3×3邻域，重复上述步骤二至步骤四，直至全部像素滤波完成。该优选技术方案直接带来的技术效果是，去噪前后恒星识别率提高了约60％，增加了观测量，大大削弱了随机噪声对于CCD恒星图像处理的影响。垂线偏差子午分量和卯酉分量的标准差相较于去噪前分别平均减小了0.016\和0.021\，内符精度有所提高。3×3中值滤波是一种基于排序统计理论的能去除噪声的非线性处理方法，基本原理是把数字图像中的灰度值用该点的3×3邻域中各点灰度值的中值代替。这里的邻域通常被称为窗口，当窗口在图像中上下左右进行移动后，利用中值滤波可以很好地对图像进行去噪处理。进一步优选，上述恒星星象区域的自动搜索是采用区域生长算法进行的。该优选技术方案直接带来的技术效果是，采用区域生长算法进行CCD恒星图像中的星象自动搜索，其方法简单，所需花费的时间较短。原因在于，区域生长算法的基本原理为设定一个最小灰度阈值，如果CCD图像中某区域由灰度值相近的若干相邻像素组成，且该区域内各像素灰度值均大于最小灰度阈值，则判定该区域为一颗恒星星象区域。进一步优选，上述改进的修正二维矩方法，是在以往的修正二维矩方法上，运用迭代计算方法对阈值进行改进，其步骤如下：第1步，提取CCD恒星图像中最大灰度值gmax和最小灰度值gmin，令T0=gmin+gmax2---(2);]]>第本文档来自技高网...

【技术保护点】
一种集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特征在于，包括以下步骤：第一步，在测站点上架设CCD天顶筒，使其光轴指向天顶方向，对天顶区域连续在0°和180°方向进行拍摄，获取CCD恒星图像，并由高精度电子水平仪获取天顶筒的倾斜数据，由GNSS获取观测时间信息和测站点的大地坐标；第二步，对CCD天顶筒拍摄的FITS格式的CCD恒星图像进行读取，再进行背景噪声消除；第三步，对噪声消除后的CCD恒星图像进行恒星星象区域的自动搜索；第四步，运用改进的修正二维矩方法，计算上述自动搜索出的恒星星象区域的能量中心的影像坐标；第五步，对适用星表进行处理，通过GNSS获取的时间信息以及测站点的大地坐标，结合EGM2008计算测站点的粗略天文坐标；根据计算的测站点粗略天文坐标对星表中与CCD恒星图像相对应的区域进行截取，并计算曝光时刻天顶区域的恒星赤道坐标与视星等信息，并将恒星赤道坐标投影到切平面上转换为平面坐标；第六步，对CCD影像中的星象与星表投影到平面上的恒星进行识别、匹配；第七步，依据匹配成功的星对，利用最小二乘法计算摄影测量投影变换公式中的坐标转换系数，并迭代计算测站点的初始赤道坐标，再转换为测站点的初始天文坐标；第八步，将高精度电子水平仪获取的天顶筒倾斜改正数据和极移改正数据转换至子午方向和卯酉方向上，确定出测站点的精确天文坐标；第九步，结合GNSS测量数据解算出的测站点大地坐标，并根据Helmert公式，即可计算出测站点的垂线偏差。...

【技术特征摘要】
1.一种集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特征在于，包括以下
步骤：
第一步，在测站点上架设CCD天顶筒，使其光轴指向天顶方向，对天顶区域连续在0°和
180°方向进行拍摄，获取CCD恒星图像，并由高精度电子水平仪获取天顶筒的倾斜数据，由
GNSS获取观测时间信息和测站点的大地坐标；
第二步，对CCD天顶筒拍摄的FITS格式的CCD恒星图像进行读取，再进行背景噪声
消除；
第三步，对噪声消除后的CCD恒星图像进行恒星星象区域的自动搜索；
第四步，运用改进的修正二维矩方法，计算上述自动搜索出的恒星星象区域的能量中心
的影像坐标；
第五步，对适用星表进行处理，通过GNSS获取的时间信息以及测站点的大地坐标，结
合EGM2008计算测站点的粗略天文坐标；
根据计算的测站点粗略天文坐标对星表中与CCD恒星图像相对应的区域进行截取，并计
算曝光时刻天顶区域的恒星赤道坐标与视星等信息，并将恒星赤道坐标投影到切平面上转换
为平面坐标；
第六步，对CCD影像中的星象与星表投影到平面上的恒星进行识别、匹配；
第七步，依据匹配成功的星对，利用最小二乘法计算摄影测量投影变换公式中的坐标转
换系数，并迭代计算测站点的初始赤道坐标，再转换为测站点的初始天文坐标；
第八步，将高精度电子水平仪获取的天顶筒倾斜改正数据和极移改正数据转换至子午方
向和卯酉方向上，确定出测站点的精确天文坐标；
第九步，结合GNSS测量数据解算出的测站点大地坐标，并根据Helmert公式，即可计
算出测站点的垂线偏差。
2.根据权利要求1所述的集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特
征在于，所述背景噪声消除所采用的方法为3×3中值滤波法，具体步骤如下：
步骤一，确定一个以某个像素为中心点的3×3邻域；
步骤二，将该邻域中各个像素的灰度值(x1,x2,x3,...,xn)按从小到大进行排序：
x1≤x2≤x3...≤xn；
步骤三，按公式计算出各个像素的灰度值的中值Y：
步骤四，以中值Y代替原来中心点的灰度值，进行中值滤波；
步骤五，确定一个以另一个像素为中心点的3×3邻域，重复上述步骤二至步骤四，直至
全部像素滤波完成。
3.根据权利要求1所述的集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特
征在于，所述恒星星象区域的自动搜索是采用区域生长算法进行的。
4.根据权利要求1所述的集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特
征在于，所述改进的修正二维矩方法，是在以往的修正二维矩方法上，运用迭代计算方法对
阈值进行改进，其步骤如下：
第1步，提取CCD恒星图像中最大灰度值gmax和最小灰度值gmin，令
T0=gmin+gmax2---(2);]]>第2步，根据阈值Tk把CCD恒星图像区分为星象区域和背景区域，其中，星象区域位
置(xi,yj)的灰度值gk+1(xi,yj)＞0，而背景区域的灰度值为0，如公式所示：
则，按公式(4)求得星象区域的平均灰度值：
gk+1=Σgk(xi,yj)>Tk(gk(xi,yi)-Tk)Σgk(xi,yi)>Tk---(4)]]>令将式(4)代入，可得：
Tk+1=Σgk(xi,yj)>Tk(gk(xi,yi)-Tk)Σgk(xi,yi)>Tk---(5)]]>第3步，设置一个正数μ，如果|Tk-Tk+1|＜μ，则终止迭代计算过程，按公式(6)求得
最优阈值：
T＝Tk+B   (6)
如果|Tk-Tk+1|≥μ，则令Tk＝Tk+1，继续进行第二步的迭代计算；
第4步，按公式(7)和式(8)，求出CCD恒星图像中星象区域中心坐标：
x0=ΣΣxi(g(xi,yj)-T)ΣΣ(g(xi,yj)-T)---(7)]]>y0=ΣΣyi(g(xi,yi)-T)ΣΣ(g(xi,yj)-T)---(8)]]>5.根据权利要求1所述的集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特
征在于，所述适用星表的处理是通过调用AURIGA中的NOVAS-F程序，并结合GNSS获取
的时间信息，对观测历元观测到的恒星位置与适用星表历元恒星平位置存在差异的各影响要
素进行改正，以获取星表中恒星的视位置；
然后，对星表中与CCD恒星星像相对应的区域进行截取，截取过程分为两步：
第一步，在集成CCD和GNSS数字天顶筒进行观测之前，首先对星表中本观测时段可能
通过测站点天顶区域的恒星信息进行截取，星表截取范围以将观测时可能出现在天顶区域视
场内的恒星全部覆盖为原则，其测站点赤道坐标(α2,δ2)由如下公式(9)求得：
δ2＝Φ2α2＝Λ2+GAST   (9)
上式(9)中，(Φ2,Λ2)为测站点的初始天文坐标，该初始天文坐标是采用GNSS大地测
量结果和EGM2008重力场模型计算出测站点的垂线偏差，再推得测站点的天文坐标，取近
似值得到的；
所述星表截取范围大于数字天顶筒视场；
在截取范围确定后，运用恒星视位置计算模型对观测时段内可能出现在该范围内的恒星
信息进行截取，并将这些信息进行存储；
第二步，在进行恒星匹配工作前，需对以上星表区域进行进一步截取，使该区域能够完
全与CCD恒星图像区域对应，其截取条件如下：
αmin＜α1＜αmaxδmin＜δ1＜δmax   (10)
上式(10)中：
(α1,δ1)为天顶区域视场内恒星的赤道坐标；
αmin、αmax、δmin、δmax取决于测站点的赤道坐标(α2,δ2)和数字天顶筒的视场Fi，其计算
公式如下：
αmin＝α2-Fi/2
αmax＝α2+Fi/2
δmin＝δ2-Fi/2   (11)
δmax＝δ2+Fi/2
所述将恒星赤道坐标投影到平面上转换为平面坐标是按投影变换公式(12)和(13)计
算得出的：
l=tan(α1-α2)cos(q-δ2)---(12)]]>m＝tan(q-δ2)   (13)
上式(12)和(13)中：
cotq＝cotδ1cos(α1-α2)；
(l,m)为恒星平面坐标；
(α1,δ1)为天顶区域视场内恒星的赤道坐标。
6.根据权利要求1所述的集成GNSS和CCD天顶筒的高精度垂线偏差快速测量方法，其特
征在于，所述恒星的识别、匹配是采用四边形算法进行的，其匹配原则...

【专利技术属性】
技术研发人员：郭金运，沈毅，王建波，代杰，刘新，孔巧丽，刘智敏，
申请(专利权)人：山东科技大学，
类型：发明
国别省市：山东;37