具有抵偿面的任意带高斯投影直角坐标系选取方法探讨.docVIP

具有抵偿面的任意带高斯投影直角坐标系选取方法探讨 　　【摘 要】本文通过对正形投影长度变换基本理论的详细研究，给出一种选取具有高程抵偿投影面的任意带坐标系的可行方案。并将该方法成功应用到实际工程当中，本文的理论推导和实例验证均表明该方案能够较好地抑制实测边长投影至高斯平面过程中产生的长度变形。 　　【关键词】投影变形；中央子午线；抵偿投影面；坐标转换；高斯平面 　　引言 　　在工程测量实践中，为保证测量结果的精确性，在限制测区内长度的投影变形时经常会遇到如何选取投影带和投影面的问题。如现行行业标准《火力发电厂工程测量技术规程》DL/T5001-2004与《电力工程施工测量技术规范》DL/T 5445-2010，均规定主测区内投影长度变形不应大于2.5cm/km，再如国家标准《330kV～750kV架空输电线路勘测规范》GB 50548-2010与《1000kV架空输电线路勘测规范》GB 50741-2012，则明确规定中央子午线的选择应考虑投影变形，投影变形值的允许偏差应不大于5cm/km要求。考虑到地面实测边长经过高斯正形投影变换到平面直角坐标系内所产生的投影变形的本征属性，本文通过同时选取恰当的高斯投影中央子午线平移参数及合适的投影抵偿面高程，以保证测区内投影长度变形，即“具有抵偿高程面的任意带高斯正形投影”法。 　　1 投影变形基本理论 　　外业测量所测得的数据的参考面为大地水准面，基准线为铅垂线，而工程图纸所使用的坐标系都为高斯平面直角坐标系。因此要保证外业测量数据的准确性必须考虑投影变形的影响。解决投影变形问题主要是平面控制网的投影面和投影带的选取问题。而投影变形主要包括长度的高程归化与投影面上的长度变形这两个问题。 　　（1）将地面上观测的长度归算至椭球面上，称为长度的高程归化，其值为ΔS1。 　　ΔS1=- HmRm D （1） 　　式中，Hm为归算边高出参考椭球面的平均高程；D为实测边长；Rm为归算边方向参考椭球法截弧的曲率半径。 　　（2）投影带对长度归算的影响，既参考椭球面的边长归算到高斯投影平面引起的长度变形量，其值为ΔS2。 　　ΔS2=y2m2R2m S （2） 　　式中，ym为归算边在相应高斯平面直角坐标系下两端点的平均横坐标值；S为椭球面上归算边对应的大地线长度。 　　由此可知长度投影变形值表达式为：ΔS=ΔS1+ΔS2。考虑到实际应用中，DΔS1， 且有S=D+ΔS1。因此我们可以认为长度投影变形值ΔS可以近似表达为： 　　ΔS =（y2m2R2m -HmRm ）D （3） 　　当取Rm≈6371000m时，根据（3）式可计算的每公里投影变形随ym和Hm变化的情况，并可归纳投影变形的主要特征如下： 　　1）地面实测长度归算至参考椭球体表面上时，会发生缩短。并且地面点与参考椭球面的高差越大、变形越大。 　　2）椭球面上长度归算至高斯面上，总是相应增长。且距离中央子午线愈远变形愈大。 　　3）由于高程归化投影变形与高斯投影变形符号相反，所以在一定的区域内，两种变形可以相互抵偿，即ΔS。 　　2 具有抵偿面的任意带高斯投影直角坐标系的选取方法探讨 　　根据对投影变形原因的分析，为了更好的抑制地面控制网长度投影变形，选取应该合适，满足“实测边长投影到该面上的测区内可能取得的最大的单位长度尽可能的小”。该约束条件可以用数学表达式表示为： 　　max{?OΔSi?O}=min （4） 　　其中，ΔSi =ΔS1i+ΔS2i=[ （ym-y0）22R2m - H0Rm ]Si，i=1，2，3，….n； 　　式中存在着两个变量， 而又不存在具有实际意义的自由极值， 要得到max{ΔSi}=min比较困难。为了简化问题， 我们可以先选取合适的中央子午线使得ΔS2最小.然后选取合适的抵偿投影面使得max{ΔSi’}=min（此后用′符号代表平移中央子午线后的坐标值）。 　　若使ΔS2i最小，对于整个测区而言，很明显y0=ymin+ymax2 。在选定中央子午线后，max{ΔSi’}用相对误差形式表述如下： 　　max{?OΔSi’?OSi }=max{?O（ymax-ymin）23R2 - H0R ?O， H0R }， （5） 　　此时，问题已经转化为了如何求取平移中央子午线后坐标转换的最佳高程抵偿面。根据最佳抵偿面选取的思想，先求出某点的y坐标y0使得在y0点ΔSi’=0且整个测区中max{?OΔSi’?O}=min。 　　由ΔSi’=0且根据（3）式得： 　　H0=y0’22R （6） 　　此时max={?OΔSi’?O }={?Oy0’2-y’max2?O， y0’2 }，且只有在 　　?Oy0’2-y’max2?O= y0’2时才有最小值。