一、引言随着2000国家大地坐标系CGCS2000在全国范围内的全面推广各地原有的城市独立坐标系面临着向CGCS2000统一过渡的任务。广州市作为国内率先完成坐标系升级的城市之一建立了基于CGCS2000椭球的新一代城市坐标系——“广州2000坐标系”。该坐标系在技术指标上实现了“100%点位误差控制在4.66厘米以内平面点位中误差1.21厘米”-。然而在实际工作中测绘、规划、国土等领域的从业人员经常面临“广州2000坐标系”与“国家2000坐标系CGCS2000”之间数据互转的需求。本文将从转换模型、参数计算、实际操作等角度系统介绍两套坐标系之间的转换方法。二、坐标系基本概念2.1 国家2000大地坐标系CGCS20002000国家大地坐标系是我国当前法定的国家大地坐标系其原点位于地球质量中心采用地心坐标系定义-。CGCS2000是一个三维空间直角坐标系所有坐标值均基于地球质心。2.2 广州2000坐标系“广州2000坐标系”是广州市基于CGCS2000椭球建立的城市独立坐标系-。它与传统“广州坐标系”基于1954北京椭球有本质区别——前者与CGCS2000“同根同源”三维空间基准完全一致区别仅在于投影参数的优化-。项目类别详细定义与规格系统名称广州2000坐标系执行启用时间2016年7月1日全市全面启用旧坐标系成果逐步转换系统性质基于CGCS2000建立的相对独立的城市平面坐标系统空间基准椭球采用2000国家大地坐标系CGCS2000的参考椭球长半轴a6378137m扁率f1/298.257222101空间基准完全一致投影模型高斯-克吕格投影横轴椭圆柱等角投影投影参数针对广州地区优化定制包含特定的中央子午线及投影面高程旨在最大限度减少区域长度变形核心精度指标- 平面点位中误差1.21 cm- 点位误差控制100%的点位误差在4.66 cm以内法律地位法定标准依据《广州市测绘地理信息管理办法》第八条明确为本市法定坐标系统与CGCS2000的关系三维同基准通过严密转换模型如布尔莎七参数可高精度互转理论上无精度损失生产实现方式通过广州市连续运行卫星定位服务系统GZCORS直接播发用户可实时获取该坐标成果建立使命终结了广州市历史上54、80及旧广州坐标系并存的混乱局面统一全市测绘基准保障跨部门空间数据共享与融合补充说明执行时间的依据为广州市原国土资源和规划委员会发布的《关于全面使用2000国家大地坐标系和广州2000坐标系的通告》过渡期内允许旧成果逐步转换但新项目必须使用新坐标系。如需查阅通告原文可访问广州市规划和自然资源局官网。三、转换模型布尔莎七参数模型广州2000坐标系与国家2000坐标系之间的转换采用的是布尔莎Bursa七参数模型-。该模型是三维空间坐标转换的经典方法包含七个独立参数-3个平移参数ΔX、ΔY、ΔZ单位米3个旋转参数εX、εY、εZ单位弧度或秒1个尺度参数k无量纲通常以ppm表示七参数模型能够完整描述两个空间直角坐标系之间存在的微小差异从而支撑厘米级乃至毫米级的空间一致性保障-。3.1 布尔莎模型数学公式布尔莎坐标转换模型的基本形式如下-其中(Xs,Ys,Zs)(Xs,Ys,Zs)为源坐标系中的空间直角坐标(Xt,Yt,Zt)(Xt,Yt,Zt)为目标坐标系中的空间直角坐标-。3.2 转换流程使用七参数进行坐标转换一般遵循以下步骤-高斯反算将源坐标系的平面投影坐标x, y转换为大地坐标B, L, H大地坐标转空间直角坐标将大地坐标转换为空间直角坐标X, Y, Z七参数转换代入七参数模型计算目标坐标系的空间直角坐标空间直角坐标转大地坐标将转换后的空间直角坐标转为大地坐标高斯正算将大地坐标投影为目标坐标系的平面坐标四、广州2000 ⇄ 国家2000 转换七参数根据勘测技术资料广州2000坐标系与国家2000坐标系之间的转换参数如下国家2000转广州2000残差计算4.1 广州2000 → 国家2000正向转换参数类型参数符号假定参数值平移参数XΔX-*.4405 m平移参数YΔY***.4741 m平移参数ZΔZ***.6464 m旋转参数XεX***秒旋转参数YεY-***秒旋转参数ZεZ-***秒尺度参数k-0***4.2 国家2000 → 广州2000反向转换反向转换可以通过对上述七参数求逆得到或直接使用以下参数参数类型参数符号参数值平移参数XΔX*** m平移参数YΔY*** m平移参数ZΔZ-*** m旋转参数XεX*** 秒旋转参数YεY***秒旋转参数ZεZ***秒尺度参数k0.***⚠️ 重要提示以上参数来源于公开技术资料在实际生产应用中建议通过已知控制点进行检核验证后再使用-1。不同转换软件的实现方式可能存在差异需注意参数的符号约定和单位统一。五、参数计算原理5.1 七参数的求解方法七参数通常通过公共点同时拥有两套坐标系坐标的控制点来求解-。具体步骤如下选取公共点选择至少3个实际应用中建议6个以上均匀分布的控制点这些点同时具有广州2000坐标系和CGCS2000坐标系的空间直角坐标。建立误差方程将每个公共点的两套坐标代入布尔莎模型建立误差方程组。最小二乘平差采用最小二乘原理求解七个参数的最优估计值-。精度评定计算转换残差评估转换精度。5.2 精度控制要点研究表明在坐标转换中平移参数的影响最大其次是旋转参数尺度参数的影响最小但尺度参数不容忽略-。因此在实际参数计算中应注意公共点的数量和分布直接影响参数求解的可靠性单一的点数过多或过少都不能提高转换精度-应进行粗差剔除保证公共点坐标成果的质量-六、二维四参数模型的补充说明除了三维七参数模型外在实际工程应用中有时也采用二维四参数模型进行平面坐标转换-1。四参数模型适用于同一椭球基准下、不同投影参数的两套平面坐标之间的转换。四参数转换公式为-1其中a,da,d为平移参数b,cb,c包含旋转和尺度信息。广州市在“三规合一”等项目中就采用了分别计算、分别转换的方法来处理不同部门之间的坐标统一问题-。七、精度验证7.1 理论精度“广州2000坐标系”在技术指标上实现了平面点位中误差1.21厘米-。计算对比实践精度高7.2 实践验证研究表明基于CGCS2000椭球建立的广州2000坐标系与原广州坐标系基于1954北京椭球的坐标成果较差在±3厘米以内实现了新旧坐标成果的平稳过渡-。城市2000坐标系参数反演算法也能够实现新建2000独立坐标系与原坐标系同名点坐标较差在±5cm以内-。八、实际操作建议8.1 参数使用注意事项确认坐标系定义在进行转换前务必确认源数据和目标数据各自的坐标系定义包括椭球参数、投影方式、中央子午线等。注意参数符号不同软件对旋转参数的正负号约定可能不同使用时需仔细核对。进行检核验证使用已知控制点对转换参数进行验证确保转换精度满足项目要求-。选择合适的模型三维空间数据使用七参数模型二维平面数据可使用四参数模型。8.2 推荐工具COORD软件支持七参数计算与坐标转换-ArcGIS支持自定义七参数坐标转换-MATLAB可编程实现七参数解算-九、总结广州2000坐标系与国家2000坐标系CGCS2000之间的转换核心是采用布尔莎七参数模型。本文提供的七参数可作为广州地区坐标转换的参考依据。在实际应用中建议通过官方渠道获取权威的转换参数或使用已知控制点对参数进行检核验证以确保转换成果的准确性和可靠性-1。参考文献*程晓晖.广州坐标系与CGCS2000坐标系转换关系的研究[J].城市勘测,2016(6):100-102*-广州2000坐标转换参数CSDN文库-12000独立坐标系参数反演算法研究及应用
