基于CGCS2000的青藏高原地区城市平面坐标系建立方法

第9卷第3期导航定位学报 Vol.9，No. 3 2021年6月Journal of Navigation and Positioning Jun.，2021
引文格式：胡斌，陈可可，黄志伟，等. 基于CGCS2000的青藏高原地区城市平面坐标系建立方法[J]. 导航定位学报, 2021, 9(3): 80-85.（HU Bin, CHEN Keke, HUANG Zhiwei, et al. Establishment of urban plane coordinate system in Qinghai Tibet plateau based on CGCS2000[J]. Journal of Navigation and Positioning, 2021, 9(3): 80-85.）DOI:10.16547/ki.10-1096.20210313.
基于CGCS2000的青藏高原地区城市平面坐标系
建立方法
胡斌1，陈可可2，黄志伟2，李春华2
（1. 四川省交通勘察设计研究院有限公司，成都610017；2. 成都市勘察测绘研究院，成都610081）
摘要：针对青藏高原地区大部分城市具有东西跨度大、海拔高、地形要素简单、人口密度低、城市建设相对缓慢等特点，探讨了处理城市全球卫星导航系统（GNSS）控制网中的投影长度变形的控制方法。

以某城市C级GNSS控制网为实例，计算分析了城市平面坐标系的三种建立方法。

结果表明：该城市采用2000国家大地坐标系（CGCS2000）的椭球参数、任意带加平均高程面（主要城市建设发展区域）的投影方法建立城市平面控制网，能够满足城市控制网对长度变形的要求。

关键词：青藏高原；2000国家大地坐标系、城市平面坐标系；长度变形；任意带；平均高程面
中图分类号：P228文献标志码：A 文章编号：2095 4999(2021)03 0080 06
Establishment of urban plane coordinate system in Qinghai Tibet plateau
based on CGCS2000
HU Bin1, CHEN Keke2, HUANG Zhiwei2, LI Chunhua2
（1. Sichuan Communication Survering & Design Institute Co. Ltd., Chengdu 610017, China；
2. Chengdu City Institute of Surveying and Investigation, Chengdu 610081, China）
Abstract：In view of the characteristics of most cities in the Qinghai Tibet plateau, such as large east-west span, high altitude, simple terrain elements, low population density and relatively slow urban construction, this paper discusses in detail the control method of processing projection length deformation in urban Global Navigation Satellite System (GNSS) control network.Taking the C-level GNSS control network of a city as an example, three methods of establishing urban plane coordinate system are calculated and analyzed.The results show that: the ellipsoid parameters of China Geodetic Coordinate System 2000 (CGCS2000) and the projection method of arbitrary zone plus average elevation plane (main urban construction and development area) are used to establish the urban horizontal control network, which can meet the requirements of the urban control network for length deformation.
Keywords：Qinghai Tibet plateau; China geodetic coordinate system 2000; urban plane coordinate system; length deformation; arbitrary zone; average elevation plane
收稿日期：2020 12 21
第一作者简介：胡斌（1968—），男，四川阆中人，大学，高级工程师，注册测绘师，研究方向为公路精密工程测量与测绘基准应用等。

通信作者简介：李春华（1969—），男，四川达州人，博士，教授级高级工程师，注册测绘师，研究方向为高精度卫星导航定位、精密工程测量等。

第3期胡斌，等. 基于CGCS2000的青藏高原地区城市平面坐标系建立方法81
0引言
城市或者工程平面坐标系统的关键问题是长度变形的控制方法。

地面观测值投影到高斯平面的长度变形，包括高程归化长度变形与高斯投影长度变形，前者是将地面边长归算到参考椭球面上，后者是将参考椭球面上的长度投影到高斯平面上。

经过以上两次投影改正后，地面实际长度被改变了，这种高斯平面上的长度与地面真实长度之差，称为长度变形[1-2]。

按照《城市测量规范》的要求，城市平面控制测量坐标系统的选择要求其投影长度变形不大于2.5 cm/km，即相对长度变形值不大于25×10-6[3]，才可以满足城市1:500地形图测量、面积测绘、市政测量（如管线、道路）等城市日常测量的要求，这也是公路、铁路、地铁、水电站等绝大部分工程平面坐标系建立的基本要求。

为此，我国绝大部分大中小城市，均建立了相对独立的平面坐标系，如北京、上海、重庆、成都、青岛、拉萨等。

随着西藏地区国民经济的持续发展以及城市化进程的不断推进，迫切需要建立统一的城市平面坐标系统与高程基准，既是城市基础测绘的重要内容，也是城市集中统一规划建设与管理的重要前提条件。

我国已于2008年7月1日正式启用了地心坐标系，即2000国家大地坐标系（China geodetic coordinate system 2000, CGCS2000）[4]。

本文在CGCS2000三维地心坐标的基础上，针对青藏高原地区城市地域分布范围广、东西跨度大、海拔高、高差大等城市地形地貌与城市布局特点，详细探讨了如何妥善处理城市全球卫星导航系统（global navigation satellite system, GNSS）控制网中的投影长度的控制方法，以满足《城市测量规范》的要求，为该地区城市平面控制网的建立提供参考，实现市级CGCS2000覆盖主要城市规划建设区域，为城市自然资源管理规划、国民经济建设等提供统一的高精度空间定位基准，保障多规合一的顺利实施。

1长度变形的控制方法
为便于数学计算与分析，国际上均采用参考椭球来描述地球。

实际工作中，首先需要将地面观测值归算到参考椭球面上，再按照一定的投影方式，将椭球面上的观测值投影到平面上，我国采用高斯-克吕格（Gauss-Kruger）投影的方式，将球面坐标转换为平面坐标，椭球面投影到高斯平面的数学模型[5-6]为
(
)
()
()
cos
cos
cos cos
cos
X X N t B l N t t
B l
Y N B l N t B l
N t t t B l
η
η
η
ηη
⎫=+⋅⋅⋅+⋅−++⎪
⎪
⋅⎪
⎪
⎬=⋅⋅+−+⋅+⎪
⎪
⎪
−++−⋅⎪
⎭
2
222
44
4
2233
2422255
11
59
224
4
1
1
6
1
5181458
120
（1）
式中：X、Y分别为高斯平面上的纵坐标与横坐标，单位为m；B为投影点的大地纬度，单位为rad；X0为自赤道量起的子午线弧长，单位为m；N为椭球面卯酉圈的曲率半径，单位为m；l为经差，即投影点纬度与中央子午线的差值，单位为rad；
tan
t B
=；cos
e B
η′
=，e′为椭球第二偏心率，详细的计算方法见文献[5]。

将地面长度（平距）观测值归算到参考椭球面
的距离改正系数
1
m的计算方法[1-2]为
1
()
1
H H
m
R
−
⎛⎞
=+
⎜⎟
⎝⎠
（2）
式中：
H为投影高程面，单位为m；H为两点间的平均高程，单位为m；R为平均曲率半径，单位为m。

式（2）表明，当地面高H小于投影高
H时，
改化到参考椭球面上的距离改化系数
1
m＞1，这意味着改化到参考椭球面上的长度较实际长度长；
当地面高H大于投影高
H时，改化到参考椭球面
上的距离改化系数
1
m＜1，这意味着改化到参考椭球面上的长度较实际长度短。

实际计算结果表明，每升高或者降低6～7m，每千米的长度改化值约为±1 mm。

参考椭球面上的长度投影到高斯平面的距离
改正系数
2
m计算公式[1-2]为
24
2
2224
1
22424
m m
Y Y
Y
m
R R R
Δ
=+++（3）
式中：
m
Y为两点之间的平均横坐标，单位为m；Y
Δ为两点间的横坐标差，单位为m。

实际计算中可近似地只取前面两项。

显然，高斯投影引起的距离改正大小呈非线性变化、且恒大于等于1。

距中央子午线愈远，长度变形越大。

实际计算结果表明，如需满足每公里长度变形不大于2.5 cm的要求，其最边缘距离投影中央子午线的距离应不大于45 km。

地面距离观测值（平距）投影到高斯平面上的
82
导航定位学报
2021年6月
综合改正系数m 为
12m m m =⋅
（4）
为控制球面到平面上的长度变形，国家一般按照3°带或者6°带的分带投影方式，以满足中小比例尺地形图，如1:1万、1:5万等地形图的测图需要，但是这种标准分带投影的方法基本上不能满足城市测量对长度变形的要求。

城市独立坐标系的建立方法概括起来有以下几种主要方式[7-10]：
1）抵偿高程面。

该方法仍然是采用3°带中央子午线（不在城市中间），加上虚拟的高程面，用来抵偿高斯投影引起的长度变化。

该方法计算较为简便，但控制面积有限。

2）任意带投影。

以某一高等级控制点的经度
或者城市范围内的平均经度0L 进行高斯投影。

该方法主要适用于地势起伏较小、海拔较低的东南沿海平坦地区。

3）平均中央子午线与平均高程面结合使用。

即采用城市平均中央子午线和城市平均高程面进行高斯投影。

该方法更适合绝大部分城市的实际情况，尤其海拔较高的内陆或者高原地区。

目前，绝大部分城市均采用这种方法建立城市基本控制网。

2 实例分析
西藏高原某市辖区位于88°25′E ～89°42′E 、28°56′N ～29°39′N 之间，东西跨度约120 km ，南北跨度约80 km ，面积约3 700 km 2，平均海拔约为
4 000 m ，主要建设区域的海拔约为3 900 m 。

在该地
区主要建设区域内，布设了21个C 级GNSS 点作为该地区的首级三维控制网，其GNSS 点的CGCS 2000坐标介于88°21′E ～89°39′E ，29°09′N ～29°26′N 之间，平均高程约为3 850 m ，GNSS 点大地高介于3 780～3 900 m 之间，地貌特征与点位分布见图1。

图1 某市地形起伏与控制点分布示意
【审图号GS （2021）3659号】
该区域位于3°带的第29度与第30度（中央子午线分别为87°、90°）范围内，但绝大部分区域位于第30度带范围内，最东端距离该中央子午线以西约35.5 km ，下面只讨论第30度带的长度改化值（或者长度变形值）。

以C 级GNSS 点的Y 坐标和大地高为边长（1 km 边长）的平均Y 坐标与平均大地高，对其3°带（第30度带）投影边长进行长度变形计算，其结果见表1，边长改化大小统计结果见图2。

表1 3°带（第30度带）投影长度变形计算表
序号 平均Y 坐标/m
距离L 0的距离/km
平均高度/m
平均纬度
每千米长度改化/mm
高程归化
高斯投影 综合 1 344 479 -155.5 3 873 29°19′42″N -608.3 298.3 -310.2 2 340 098 -159.9 3 885 29°12′52″N -610.1 315.4 -294.9 3 361 765 -138.2 3 962 29°09′37″N -622.2 235.7 -386.7 4 359 738 -140.3 3 857 29°21′34″N -605.7 242.7 -363.2 5 372 662 -127.3 3 826 29°21′30″N -601.0 200.0 -401.1 6 378 208 -121.8 3 854 29°10′34″N -605.3 183.0 -422.4 7 386 801 -113.2 3 843 29°21′54″N -603.5 158.1 -445.6 8 400 318 -99.7 3 799 29°20′10″N -596.6 122.6 -474.1 9 412 549 -87.5 3 975 29°25′55″N -624.4 94.3 -530.1 10 422 499 -77.5 3 878 29°20′19″N -609.1 74.1 -535.1 11 434 083 -65.9 3 799 29°20′22″N -596.7 53.6 -543.2 12 451 265 -48.7 3 793 29°20′07″N -595.7 29.3 -566.4 13
398 354
-101.6
3 983
29°10′26″N
-625.5
127.4
-498.2
第3期胡斌，等. 基于CGCS2000的青藏高原地区城市平面坐标系建立方法83
（续）序号平均Y坐标/m 距离L0的距离/km 平均高度/m平均纬度
每千米长度改化/mm
高程归化高斯投影综合14422773 -77.2390829°14′21″N -613.8 73.6-540.2 15438019 -62.0386029°16′44″N -606.2 47.4-558.9 16456030 -44.0399829°16′43″N -628.0 23.8-604.2 17464522 -35.5378429°20′38″N -594.2 15.5-578.7 18392022-108.0381829°16′04″N -599.6143.8-455.9 19406637 -93.4386229°09′43″N -606.6107.5-499.1 20419640 -80.4384329°10′29″N -603.6 79.7-524.0 21392518-107.5380629°15′00″N -597.8142.5-455.4
表1中：
L为高斯投影中央子午线。

图2CGCS2000坐标3°带投影长度改化分布
表1表明，每千米高程归化、高斯投影以及综合长度改化分别介于-629~-594、15.5~315.4、-604.2~-294.9mm之间，远远无法满足《城市测量规范》每千米不大于2.5 cm的要求。

下面分别就抵偿面投影、任意带投影以及任意带（平均中央子午线）加平均高程面三种投影方式控制长度变形的方法进行详细讨论。

2.1抵偿面投影
根据表1每千米综合长度改化的结果分析，该地区21条边的每千米的平均长度改正约为-475.6 mm，对应的抵偿高程面约为3030 m（大地高），在3°带的基础上，采用该抵偿面进行高斯投影，投影后各长度改化量见图3。

图3抵偿高程面投影长度改化分布
采用抵偿高程面投影后，虽然该地区21条边长的每千米平均综合长度变形量不足1 mm，但不大于2.5 cm/km的仅有5条边，仅占24%，大部分地区仍然不能满足要求，进一步证实了抵偿高程面控制长度变形的局限性，且该抵偿面与实际平均高程面不一致，没有实际意义。

2.2任意带投影
以该地区平均中央子午线（约89°）进行任意带高斯投影，计算结果表明，高斯投影每千米的距离改化介于0~48.5 mm之间，其中仅有5条边的每千米的长度改化大于25 mm，而大于35 mm的距离改化仅剩3条边，距中央子午线分别为59、62、63 km。

综合距离改化量介于-559.7~-625.3 mm 之间，其大小主要为高程归化（即海拔高）的贡献。

详细距离改化分布见图4。

图4任意带高斯投影长度改化分布
2.3平均中央子午线加平均高程面投影
以上分析表明，选取任意带投影基本上可以满足绝大部分地区高斯投影长度变形的要求，同时，鉴于该主要建设区域地势起伏不大，拟选定任意带中央子午线加平均高程面的投影方式建立城市独立坐标系。

该城市范围东西跨度约120 km，按照式（3）及3.3节中的计算结果说明，宜选取两条中央子午线进行任意带投影，才能满足长度变形不大于2.5 cm/km
的要求。

但考虑该地区地形地貌特征与城市建设
84导航定位学报2021年6月
发展的特点、以及多带投影，即两个或两个以上的城市独立平面坐标系不便于日常使用与管理的基本原则，拟采用单一任意带（城市平均中央子午线）投影方式加城市主要建设地区的平均高程面的高斯投影方法建立城市统一的独立坐标系。

通过C级GNSS点以及其它资料分析后，拟选定的投影参数如下：CGCS2000参考椭球、城市独立平面坐标系高斯投影的中央子午线为89°、投影面大地高为3860 m。

投影后相应的边长投影变形见表2，边长改化大小见图5。

表2城市独立坐标系下长度变形计算表
序号平均Y坐标/m 距离L0的距离/km 平均高度/m平均纬度
每千米长度改化/mm
高程归化高斯投影综合1441583-58.4387329°19′42″N -2.142.1 40.0 2437308-62.7388529°12′52″N -3.948.5 44.6 3459036-41.0396229°09′37″N -16.020.7 4.7 4456819-43.2385729°21′34″N 0.523.0 23.5 5469750-30.2382629°21′30″N 5.311.3 16.6 6475472-24.5385429°10′34″N 1.0 7.4 8.4 7483890-16.1384329°21′54″N 2.7 3.2 5.9 8497442 -2.6379929°20′10″N 9.6 0.1 9.7 9509588 9.6397529°25′55″N -18.1 1.1-17.0 10519632 19.6387829°20′19″N -2.8 4.8 1.9 11531222 31.2379929°20′22″N 9.512.0 21.5 12548418 48.4379329°20′07″N 10.628.9 39.5 13495630 -4.4398329°10′26″N -19.3 0.2-19.0 14520000 20.0390829°14′21″N -7.5 4.9 -2.6 15535218 35.2386029°16′44″N 0.015.3 15.3 16553239 53.2399829°16′43″N -21.835.0 13.2 17561675 61.7378429°20′38″N 12.046.9 58.9 18489206 -10.8381829°16′04″N 6.6 1.4 8.1 19503929 3.9386229°09′43″N -0.3 0.2 -0.1 20516926 16.9384329°10′29″N 2.6 3.5 6.2 21489719-10.3380629°15′00″N 8.4 1.3 9.7
图5城市独立坐标系下长度改化分布
与表1数据比较表明，该方法东西两边距离任意带中央子午线的距离约62 km，每千米长度变形略大于25 mm，其余大部分地区均能满足要求。

考虑到东西边缘的城市建设活动不频繁，为便于应用管理，实际工作中可根据需要适当放宽长度变形的要求，特殊工程可单独进行技术处理。

3结束语
青藏高原地区大部分城市具有东西跨度大、海拔高、地形要素简单、人口密度低、城市建设相对缓慢等特点，在充分估计城市建设活跃频繁区域的基础上，采用CGCS2000的椭球参数、任意带加平均高程面（主要城市建设发展区域）的投影方法建立城市平面控制网，能够满足城市控制网对长度变形的要求，极大地避免了多投影带对城市建设管理带来的不便。

同时，可以通过一定的技术手段，如联测国际GNSS服务组织（International GNSS Service, IGS）站点或地震监测与省市级连续运行参考站（continuously operating reference stations, CORS）站点，采用精密星历进行数据处理，可实现三维基准的动态维护与青藏高原板块运动变形监测。

第3期胡斌，等. 基于CGCS2000的青藏高原地区城市平面坐标系建立方法85
参考文献
[1]李春华, 张瑛, 罗贤茂, 等. 城市控制网改造中的关键问题[J]. 测绘科学, 2008, 33(1): 77-82.
[2]李春华. 基于网络GPS和精化大地水准面的区域实时三维定位理论与应用[D]. 成都: 西南交通大学, 2010.
[3]中华人民共和国住房与城乡建设部. 城市测量规范: CJJ/T 8-2011[S]. 北京: 中国建筑工业出版社, 2011.
[4]欧海平, 秦亮军, 喻永平, 等. 广州2000坐标系建立技术与应用[J]. 工程勘察, 2018 , 46(11): 55-58.
[5]孔祥元, 郭际明, 刘宗泉. 大地测量学基础[M]. 武汉: 武汉大学出版社, 2010.
[6]徐欢, 张棋, 魏猛. 基于CGCS2000椭球的大地测量实用公式[J]. 导航定位学报, 20153(3): 105-109.
[7]李江卫, 解斌, 蔡国兴, 等. 基于CGCS2000的城市平面坐标系最佳选取[J]. 城市勘测, 2011(2): 118-123.
[8]陈顺宝, 任建春, 亓月, 等. 抵偿任意带高斯投影平面坐标系选择的研究[J]. 测绘通报, 2005(7): 21-24.
[9]谯勇, 陈达玉, 刘俊, 等. 城市GNSS控制网建立的基本方法与精度分析[J] . 测绘科学, 2011, 36(增刊1): 252-
253.
[10]周飞. 拉萨城市区域控制网加密及投影转换[J]. 城市勘测, 2019(6): 117-119.
（上接第59页）
[6]YOUSSEF M, AGRAWALA A. The Horus WLAN location determination system[J]. Wireless Networks, 2005,
14(3): 357-374.
[7]BAHL P, PADAMANABHAN V N. RADAR: an in-building rf-based user location and tracking system[EB/OL].
[2020-07-28].https://sci-hub.se/10.1109/INFCOM.2000.832252.
[8]NASRABADI N M. Pattern recognition and machine learning[J]. Journal of Electronic Imaging, 2017, 16(4):
049901-049923.
[9]李欢, 刘晓东, 赵翔, 等. 混响室随机多径衰落电磁环境的实验研究[J].无线电工程, 2019, 49(8): 715-719.
[10] RIZK H, TORKI M, YOUSSEF M. CellinDeep:robust and accurate cellular-based indoor localization via deep
learning[J]. IEEE Sensors Journal, 2018, 34(6): 34-43.
[11]张永顺, 朱卫纲, 贾鑫, 等. 基于块稀疏贝叶斯学习的直扩通信窄带干扰检测与参数估计[J]. 系统工程与电
子技术, 2019, 41(4): 889-897.
[12]YOUSSEF M, AGRAWALA A. Continuous space estimation for WLAN location determination systems[EB/OL].
[2020-07-28]./~moustafa/papers/ic3n04.pdf.
[13]WANG X, GAO L, MAO S, et al. DeepFi: Deep learning for indoor fingerprinting using channel state information
[EB/OL].[2020-07-28].https://sci-hub.se/10.1109/WCNC.2015.7127718.。