地形测绘成果控制系统一转换方法探究.docVIP

地形测绘成果控制系统一转换方法探究摘要：通过结合实际任务探讨电子地形图平面、高程系统转换方法，解决利用已有资源快速高效完成市政工程设计所需测量任务。 关键词：禅西大道电子地形图平面、高程系统转换 0前言 在全国市政建设的大潮中，有很多大、中、小城市平面、高程系统不统一，各家建设方、设计方采用不同的坐标系、高程系，为解决不必要的测绘资源浪费，快速高效利用已有地形测量成果资源完成各项工程建设所需的统一控制系测量任务，满足工程设计需要测量资料。下面以广东省佛山市禅西大道测量为例地形测绘成果控制系统一转换方法研究探讨。 在广东省佛山市禅西大道二期测量任务（5.5公里长、600米宽带状电子地形图）中，利用当地有关单位已有资源（1：500电子地形图），其中3.0公里是1954北京坐标系、黄海高程系，另2.5公里是南海坐标系、黄海高程系。而设计需要的是1：500佛山坐标系、珠江高程系电子地形图满足投影变形规范要求。因此需要对已有电子图资料进行平面、高程统一系统转换。1954年北京、南海、佛山坐标系3度带中央子午线经度L3中分别为：114度、113度、114度，佛山坐标系是采用高斯正形投影3度带平面直角坐标系、投影面采用抵偿高程面。两坐标系间的转换参数一般是利用重合点的两套坐标值通过一定的数学模型进行计算，然后利用所求得的转换参数对另一部分待联测点进行转换。一般来说，在小范围测区内，有三个公共点就可以进行两坐标系间网眼式的坐标转换，坐标系间转换参数有平面系统的平移、旋转、尺度因子等；高程系统采用联测三个以上已知?程点，求高程系统差平均值的方法进行高程系统转换。 1 解算控制系统间转换参数 1.1 控制系统间平面系统转换参数 三种坐标系：O-XY为1954年北京坐标系；O-X’Y’为南海坐标系；o-xy为佛山坐标系。平移值以A、B表示，旋转值以a表示，尺度因子以λ表示。 以下平面计算均以此三种坐标系间关系为依据。 已知条件由电子图提供单位提供如下表： 点名 坐 标 系 1.1.11954年北京坐标系和佛山坐标系间平面系统转换参数 根据已知条件算得下表： 方向边 1954年北京坐标系 佛山坐标系 较值 距离（米） 距离（米） 距离差值/54北京距离（ppm） 坐标方位角（度、分、秒） 坐标方位角（度、分、秒） 坐标方位角a（度、分、秒） 1.1.1.1解算旋转值a、尺度因子λ： 考虑现场三个已知点和设计道路中心线几何关系，J点离道路中心线较X和R点远,故边R-J权取1/3，边R-X权取2/3。故按照加权平均值解算方法得：a=-00’11.38’*1/3+(-0’0’12.76’*2/3)=-0’0’12.3’，故按照加权平均值解算方法得尺度因子： λ=1-（110.204*1/3+118.214*2/3）ppm=1-115.544ppm=0.999884456 1.1.1.2解算平移值A、B值： 根据1954年北京坐标系和佛山两坐标系间关系可得出： AXcosa-sinax =-λ ** BYsina cosay 计算过程用CASIOfx-4500PA计算器编小程序输入1954年北京坐标系各点坐标X、Y得出输入各点的考虑旋转和尺度因子的坐标值，再分别减输入各点的佛山坐标系坐标x、y值，即得到各个点算得的平移值A、B，并比较取A、B加权平均值。 程序：XY:Fixm:Pol(X,Y:Rec(0.999884456v,w+00’12.3’?w 以旋转值a=-00’12.3’、尺度因子λ=0.999884456计算结果表： a值 1.1.1.3据以上解算得到的转换参数A、B、a、λ，将1954年北京坐标系的各已知点坐标转换为佛山坐标系坐标 根据两坐标系间关系可得出： xcosa sinaX-A =(1/ λ)* * y -sina cosaY-B 计算过程用CASIOfx-4500PA计算器编小程序输入1954年北京坐标系各点坐标X、Y直接得出所求各点的佛山坐标系坐标x、y：程序：XY:Fixm:Pol(X,Y:441.6353+Rec(0.999884456v,w+0’0’12.3’?w-288.6528 根据程序算得三个已知点的佛山坐标系坐标比较结果表： 点名 坐 标 系 佛山x、y 差值（m） 1954北京坐标系 已知佛山坐标系坐标 算得佛山坐标系坐标 通过上表四结果X、R点转换后佛山坐标差值比J点小，是因为J点离设计主线远，加权时J点权取最小，利用三个公共点转换引起的加权位置中误差计算得到是±1.03cm，所以确定