大风江调水工程控制网的投影变形抵偿

叶长干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

地址:广西南宁市建政路5号  邮编:530023  Tel:0771-5606397  Email:webmaster@digitalgx.com
广西基础地理信息中心版权所有 2005-2010 广西基础地理信息中心制作