高斯投影正反算

《高斯投影正反算》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、高斯投影正反算学院：资源与环境工程工程学院专业：测绘工程学号：X51414012姓名：孙超一、高斯投影概述想象有一个椭圆柱面横套在地球椭球体外面，并与某一条了午线相切，椭I员I柱的中心轴通过椭球体的中心，然后用一定投影方法，将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上，再将此柱面展开即成为投影面。高斯投影由于是正形投影，故保证了投影的角度不变性，图形的相似性以及在某点各方向上长度比的同一性。由丁采用了同样法则的分带投彩，这即限制了长度变形，又保证了在不同投影带中采用相同的简便公式和数表进行变形引起的各项改止的计算，并且带与带间的互相换算也能用相同的公式和方法进

2、行。高斯投影的这些优点必将使它得到广泛的推广和具有国际意义。二、高斯投影坐标正算公式1•高斯投影必须满足以下三个条件1）屮央子午线投影后为直线2）中央子午线投影后t度不变3）投影具冇正形性质，即正形投影条件2•高斯止算公式推导1）由第一个条件可知，由于地球椭球体是一个旋转椭球体，所以高斯投影必然有这样一个性质，即屮央子午线东西两侧的投影必然对称于屮央子午线。2）由于高斯投影是换带投影，在每带内经丼1是不大的，丄是一个微小量，所以P可以将X=X（1,q）,Y=Y（1,q）展开为经差为1的幕级数，它可写成如下的形式X=mo+m2l2+iiiil4+*e*Y=mil+m3l+

3、m5l°+••-式中m0,ml,m2,…是待定系数，他们都是纬度B的函数。3）由笫三个条件：竽¥和务-学,将上式分别对1和q求偏导odqddqx=m0+mJ+m2l2++m4/4+••…y=n()+q/+n2l2+n3l3+n4/4+-ax一两办一刃空力空刃厂丨——Idmnrdm.小dm，dm>“dm«——+/——+/2——+1’——L+r——+dqdqdqdqdq叭+2附2?+3mJ，+AmJ2+....如“也1_+厂如*尸虹+/4如+dqdqdqdqdqnr+2nd+3n,/2+4nJ34-....可得到下式dm..1dm.1dg1dm.-—，吗=g=，坷=，…

4、rdq~2dq33dq44dqVdiu1dn}1dn91dn.mi=_/774=--—<••dq2dq3dq4dq经过计算可以得出兀=X+仝―sinBcosB-/2+卫一sinBcos3B(5-t2+9772+474)/4224+-^sinBcos5B(61-58/2+Z4)/6720y=NcosB•/4-cos3B(1-Z2+72)/3・6+cos53(5-18r2+r4+14/72-58t?2t2)l5120三、高斯投影坐标反算公式推导1•思路：级数展开，应用高斯投影三个条件，待定系数法求解。2.投影公式在底点处展开q=f；(x,y)1=(兀,y)展开为fff9f3

5、/4q=m（）+77?!y+m2+y+m4y+l=n；+n[y+nf2y2+n^y3+ny4+.....3•引入高斯投影条件之一：正形条件—-m[+2加；丿+3my2+4/M*y3+....dldnLdn[zdn3dn4dn—=―+y—4-y—+y—+y—+dxdxdxchdxdxQJ=〃；+2n2y+3zi；y2+4ny34-....°y4•曲于可得到儿=0,带入上式可得到q=傀++加；）』+••…/=n[y+〃；y3+n5y5+5•引入高斯投影条件Z三：中央了午线投影后长度不变傀=Qffdm'secBfn.=——-=dxN$卅二1讥f严Bf22dx2N；

6、,sec/r9禺=_6N；（1+2匚+77；）,ffSecBf9°%=242：（5+6f；+»wcBf（5+2勺；+6“；）5]20N；f，f四、高斯投影的特点1•当I等于常数时，随着B的增加x的值增大，y的值减小，无论B值为止或为负，y值不变。这就是说，椭球面上除中央子午线外，其它子午线投影后，均向屮央子午线弯曲，并向两极收敛，同吋还对称于屮央子午线和赤道。2.当B等于常数时，随着I的增加，x值和y值都增大。所以在椭球面上对称于赤道的纬圈，投影后仍成为对称的曲线，同时与子午线的投影曲线互相垂直凹向两极。3.据屮央子午线越远的子午线，投影后弯曲越厉害，长度变形越大。五、

7、MATLAB编程实现坐标正反算1編写main函数functionmaindispf欢迎使用高斯投影正反算及相邻带的坐标换算程序J;disp（'l：高斯正算2：高斯反算3：换带计算'）；K=0;while（K<1

8、

9、K>3）K二input（'请根据上列选择计算类型K='）；switchKcase1GSZS;case2GSFS;case3HDJS;otherwisedisp（'K值无效（1-3）'）；enddispf程序作者：亚里士多墩J;dispC指导老师：亚里士多德J;end2.编写高斯正算GSZS函数functionGSZS%GSZS是