大风江调水工程控制网的投影变形抵偿
叶长干1 叶达忠2
(1.钦州市国土资源局 广西钦州 535000;
2.广西水利电力勘测设计研究院 广西南宁 530023)
【摘 要】 本文结合大风江调水工程,应用TrimbleGPS全球定位系统施测GPS工程测量网,采用平差软件包CosaGPS进行数据处理,并就工程控制网的投影变形问题、抵偿计算和解决方案进行了讨论。
【关键词】 GPS 投影变形 高程抵偿
1 工程概况
为进行大风江调水工程技术设计阶段的水库淹没实物指标调查工作,需进行库区征地界桩测量工作以及1:2000专用地形图测绘工作,钦州市国土资源局受委托承担此项工作。测区跨度:东经108°37′~109°00′(原属3度带第36带,中央子午线经度108°);中心经度L约108°48′;纬度跨度:北纬21°50′~21°57′,中心纬度B约21°53′;淹没线高程2.5m,平均高程约2.5米。
工程需要库区、闸口分别布设平面控制网,各个控制网必须联成一体,保持统一的系统和整体的精度,同时要求满足测图及淹没界桩放样的要求。
观测仪器采用Trimble 5700双频GPS接收机。平差软件:武汉大学研制的GPS工程测量网通用平差软件包CosaGPS;Trimble的随机软件Trimble Geomatics Office Software。
2 作业依据
作业依据:1、《全球定位系统(GPS)测量规范》,(GB/T 18314-2001);2、《水利水电工程测量规范》(规划设计阶段),(SL 197-97);2、《水利水电工程测量规范》(规划设计阶段),SL 197-97,水利部、电力工业部。
根据规范的要求,投影变形允许值为5cm/km,但因为要同时满足测图及淹没界桩放样的要求,所以我们控制投影变形允许值为2.5cm/km。
3 投影原理及变形量计算
参考椭球面变形通常可以抵消一部分长度变形,将地面实测边长S归算到参考椭球面上,边长变化计算公式:
△s1=-(S×Hm)/R (5)
式中:S—归算边长;
Hm—归算边高出参考椭球面的平均高程;
R—归算边方向参考椭球法截弧的曲率半径。
将参考椭球面上的边长S0归算到高斯投影面上,边长变形影响计算公式:
△s2=S0×(Ym/Rm)2 /2 (6)
式中:S0—投影归算边长;
Ym—归算边两端点横坐标平均值;
Rm—参考椭球面平均曲率半径。
变形补偿计算:
欲使两种变形抵消,即要(5)+(6)=△s1+△s2=0mm,则 ,即: 。根据工程需要,当△s1+△s2剩余的变形量大于2.5cm/Km时,应当考虑采用高斯正形投影任意带平面直角坐标系统,或建立独立坐标系统。
表1 测区起算数据
数据 |
中心经度L |
原3度带中央子午线经度L0 |
与中央子午线经差l |
与中央子午线经差l |
中心纬度B |
中心纬度B |
距离中央子午线Ym(高斯正算) |
平均高程Hm |
投影边长S |
单位 |
/°′″ |
/°′ |
/° |
/弧度 |
/°′″ |
/弧度 |
/Km |
/m |
(m) |
值 |
108.4500 |
108 |
0.75 |
0.0131 |
25.2000 |
0.442150077 |
75.5 |
20 |
1000 |
表2 曲率半径计算
椭球 |
长半轴a |
短半轴b |
极点处子午线曲率半径c |
第2偏心率e |
辅助函数V |
参考椭球面平均曲率半径Rm |
参考椭球面卯酉圈曲率半径N |
参考椭球面子午圈曲率半径M |
/m |
/m |
/m |
无 |
无 |
/Km |
/Km |
/Km |
1954克氏椭球 |
6378245 |
6356863.01877304 |
6399698.9017827 |
0.082088522 |
1.00275E+00 |
6364663 |
6382157 |
6347217 |
1975国际椭球 |
6378140 |
6356755.28815752 |
6399596.6519880 |
0.082094469 |
1.00275E+00 |
6364556 |
6382052 |
6347108 |
表3 变形量计算(抵偿前)
边长高出参考椭球面的平均高程Hm |
归算边方向椭球法截弧曲率半径R |
参考椭球面变形△S1 |
边长两端横坐标平均值Ym |
平均曲率半径Rm |
投影带变形△S2 |
总变形量的绝对值 |
/m |
/Km |
/m |
/Km |
/Km |
/m |
/m |
20 |
6365 |
-0.003 |
75.5 |
6365 |
0.07037 |
0.067 |
由表1~表3看出,总变形量为6.7cm/km超出规范的2.5cm/km,必须进行抵偿改正。
4 抵偿计算
投影抵偿计算见表4、表5。
表4 投影抵偿计算
投影高程后△H/m |
归算边方向椭球法截弧曲率半径R/Km |
参考椭球面变形△S1/m |
测区中部离子午线距离Ym/Km |
平均曲率半径Rm/Km |
投影带变形△S2/m |
总变形量/m |
备注 |
270 |
6365 |
-0.042 |
75.5 |
6365 |
0.070 |
0.028 |
纯高程抵偿面方案 |
20 |
6365 |
-0.003 |
25.0 |
6365 |
0.008 |
0.005 |
综合方案 |
表5 高斯反算求任意带经度
高程投影面 |
测区中部离子午线距离Ym |
参考椭球面卯酉圈曲率半径N |
中心纬度余弦cosB |
经差 |
测区中心经度L |
任意带(减)经度 |
任意带(加)经度 |
/m |
/Km |
/Km |
/无 |
/弧度 |
/弧度 |
/°′ |
/°′ |
-250 |
75.5 |
6382 |
0.90383 |
0.0131 |
1.8980 |
108.00 |
109.30 |
0 |
25 |
6382 |
0.90383 |
0.0043 |
1.8980 |
108.30 |
109.00 |
5 抵偿方案
根据表5的数据,我们确定选择方案2,投影带选在109度,投影面0米,抵偿后变形量残差为5mm/km,符合规范要求,方案可行。方案1虽然不改变投影带,但需要投影到-250米的高程面,不符合习惯做法。
6 结束语
本文最后采用方案2进行投影计算,结果与全站仪边长比较,非常吻合,证明上述计算无误,达到预期目的,同时解决了投影变形的问题,提高了控制网精度,并方便了用常规仪器进行界桩放样。
【参考文献】
[1] 叶达忠.控制测量中的正形投影变形及其计算.广西水利水电.2005.03
|