区域似大地水准面的二次精化及其软件实现.docVIP

区域似大地水准面的二次精化及其软件实现 　第３８卷第６期　２０１４年１２月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ 物　探　与　化　探 Ｖｏｌ．３８，Ｎｏ．６　Ｄｅｃ．，２０１４　 ｄｏｉ：１０．１１７２０／ｗｔｙｈｔ．２０１４．６．３５２０１４．６．３５ ３８（６）：１２９０－１２９２．ｈｔｔｐ：／／ｄｏｉ．ｏｒｇ／１０．１１７２０／ｗｔｙｈｔ．２０１４．６．３５ （１．南京地质调查中心，江苏南京　２１００１６；２．南京林业大学，江苏南京　２１００３７） 摘要：为了满足区域重力调查高精度测地工作的要求，实现不同坐标系之间的转换，通过区域似大地水准面的二次精化方法及开发的软件，将ＧＰＳ大地高转换到正常高。采用实测２６个水准点外检核精度，全区软件转换精度为关键词：似大地水准面；坐标转换；二次精化；ＧＰＳ 中图分类号：Ｐ６３１；Ｐ２２８．４ 文献标识码：Ａ 文章编号：１０００－８９１８（２０１４）０６－１２９０－０３±５．５７ｃｍ，完全满足重力测量中高程精度要求。 张宝松１，郑加柱２，滕龙１，邸兵叶１，朱红兵１，张俊１ 　　目前，１∶５万高精度重力勘查测地均使用高精度ＧＰＳ进行测量。ＧＰＳ测量的坐标系是世界大地坐标系即ＷＧＳ?８４坐标系。在重力调查中使用的三维坐标系为１９５４年北京坐标系与１９８５国家高程基准，因此对ＧＰＳ测量的三维定位成果必须进行换算［１－７］。目前常用的有七参数转换法与三参数转换法，在区域范围较小，转换精度能满足重力调查测地平面与高程精度要求。当重力调查区范围较大，直在重力调查中使用高程异常进行高程求取可以较好接参数转换将导致高程转换存在较大的误差［８］。地保证高程精度，从而满足重力调查测地精度要求。通过增加高等级水准点资料，实现区域似大地水准面精化与软件编制，并对二次精化后的成果进行了检核，从而保障了调查区内重力测点高程转换精度。 笔者在省级似大地水准面精化成果基础上［９］， 值，组成不符值序列； 标组成多项式拟合“观测方程”，其中未知参数为多项式系数； （３）按最小二乘原理求解拟合多项式系数； （２）由不符值序列和相应ＧＰＳ水准点的球面坐 标，对似大地水准面进行拟合纠正，即可求得适配于该区域的ＧＰＳ水准网的最终似大地水准面。 以三次多项式为例，其数学模型为［９］ Δζ＝α０＋α１（φ－φｍ）＋α２（λ－λｍ）＋ α５（λ－λｍ）２＋α６（φ－φｍ）３＋α７（φ－φｍ）３（λ－λｍ）＋ （４）由拟合多项式系数和区域格网中心点坐 α３（φ－φｍ）２＋α４（φ－φｍ）（λ－λｍ）＋α８（φ－φｍ）（λ－λｍ）２＋α９（λ－λｍ）３ １　区域似大地水准面的二次精化数据模型 水准点成果，采用多项式拟合［６－７，９］，进一步精化测区局域似大地水准面。其计算方法、步骤和数学模型如下： （１）用省级２．５′×２．５′似大地水准面格网内插利用省级似大地水准面成果和外业观测的ＧＰＳ 式中：α０，…，α９为多项式拟合系数；Δζ为ＧＰＳ水准似大地水准面与现有似大地水准面之差，即Δζ＝ ２　实例及软件开发 件［９］。 ζＧＰＳ－ζｇｅｏｉｄ，φｍ、λｍ分别为拟合区的中心纬度和经度。 选用江苏１∶５万某重力调查项目，说明区域似大地水准面精化在测地中的应用实例与开发的软２．１　已有资料 ＧＰＳ水准点上的似大地水准面高（ζｇｅｏｉｄ），并求解与 ＧＰＳ水准点上的实测似大地水准面高（ζＧＰＳ）的差 收稿日期：２０１３－１２－０６；修回日期：２０１４－０２－２４ 基金项目：中国地质调查局地质调查工作项目（１２１２０１１０８７０４１） （１）分辨率为２．５′×２．５′的江苏省似大地水准面 ６期张宝松等： 图１　ＧＰＳ水准点分布位置 资料。 ５个点，点位分布见图１。 水准点，其中江苏范围内共２１个点，安徽范围内共 图２　似大地水准面精化流程 （２）测区内部及周边范围均匀选取２６个ＧＰＳ ２．２　似大地水准面精化及软件包开发 ５′×２．５′的格网似大地水准面数据。小地区似大地水准面精化的基本原理，是拟合测区内的ＧＰＳ水准似大地水准面和已有的似大地水准面上的差值，通过最小二乘法将二者的差异平滑处理，进一步消除和减小二者的差异。其基本流程如图 ２。 高程ｍ１３．７２８１１．４９６１２．９６７１０．５６６１０．１８８１４．６０２９７．６２６８８．１