资源描述：
《高斯投影坐标正反算vb程序》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、高斯投影坐标正反算高斯投影坐标正反算学院：班级：学号：姓名：课程名称:指导老师：-12-高斯投影坐标正反算实验目的：1.了解高斯投影坐标正反算的基本思想；2.学会编写高斯正反算程序，加深了解。实验原理：高斯投影正算公式中应满足的三个条件：1.中央子午线投影后为直线；2.中央子午线投影后长度不变；3.投影具有正形性质，即正形投影条件。高斯投影反算公式中应满足的三个条件：1.x坐标轴投影成中央子午线，是投影的对称轴；2.x轴上的长度投影保持不变；3.正形投影条件，即高斯面上的角度投影到椭球面上后角度没有变形，仍然相等。操作工具：计算机中的VB6.0代码：DimaAsDouble,bA
2、sDouble,xAsDouble,yAsDouble,y_#Diml_AsDouble,b_AsDouble,a0#,a2#,a4#,a6#,a8#,-12-高斯投影坐标正反算m2#,m4#,m6#,m8#,m0#,l0#,e#,e1#Dimdeg1AsDouble,min1AsDouble,sec1AsDouble,deg2AsDouble,min2AsDouble,sec2AsDoublePrivateSubCommand1_Click()Dimx_AsDouble,t#,eta#,N#,W#,k1#,k2#,ik1%,ik2%,dh%deg1=Val(Text1.Text
3、)min1=Val(Text2.Text)sec1=Val(Text3.Text)deg2=Val(Text4.Text)min2=Val(Text5.Text)sec2=Val(Text6.Text)l_=(deg1*3600+min1*60+sec1)/206265b_=(deg2*3600+min2*60+sec2)/206265dh=Val(Text9.Text)k1=((l_*180/3.14159+3)/6)k2=(l_*180/3.14159/3)ik1=Round(k1,0)ik2=Round(k2,0)Ifdh=6Thenl0=6*ik1-3Else-12-高斯
4、投影坐标正反算Ifdh=3Thenl0=3*ik2ElseMsgBox"error",48,"error":ExitSubEndIfEndIfl=l_-l0*3.14159/180e=Sqr(a*a-b*b)/am0=a*(1-e*e)m2=e*e*m0*3/2m4=e*e*m2*5/4m6=m4*e*e*7/6m8=e*e*m6*9/8a0=m0+m2/2+m4*3/8+m6*5/16+m8*35/128a2=m2/2+m4/2+m6*15/32+m8*7/16a4=m4/8+m6*3/16+m8*7/32a6=m6/32+m8/16a8=m8/128x_=a0*b_-a2*S
5、in(2*b_)/2+a4*Sin(4*b_)/4-Sin(6*b_)*a6/6+Sin(8*b_)*a8/8t=Tan(b_)e1=Sqr((a*a-b*b)/(b*b))-12-高斯投影坐标正反算eta=Sqr(e1*e1*Cos(b)*Cos(b))W=Sqr(1-e*e*Sin(b_)*Sin(b_))N=a/Wx=x_+N*Sin(b_)*Cos(b_)*l*l/2+N*Sin(b_)*Cos(b_)^3*(5-t*t+9*eta*eta+4*eta^4)*l^4/24+N*Sin(b_)*Cos(b_)^5*(61-58*t*t+t^4)*l^6/720y=N*Cos
6、(b_)*l+N*Cos(b_)^3*(1-t*t+eta*eta)*l*l*l/6+N*Cos(b_)^5*(5-18*t*t+t^4+14*eta*eta-58*eta*eta*t*t)*l^5/120Text7=xIfdh=6Theny_=y+500000+1000000*ik1ElseIfdh=3Theny_=y+500000+1000000*ik2ElseMsgBox"error",48,"error":ExitSubEndIfEndIfText8=y_-12-高斯投影坐标正反算EndSubPrivateSubCommand2_Click()Dimbf#,j%,Wf#,
7、Vf#,Nf#,Mf#,c#,tf#,etaf#,dh%,ik%x=Val(Text7.Text)y_=Val(Text8.Text)e=Sqr((a*a-b*b)/(a*a))m0=a*(1-e*e)m2=e*e*m0*3/2m4=e*e*m2*5/4m6=m4*e*e*7/6m8=e*e*m6*9/8a0=m0+m2/2+m4*3/8+m6*5/16+m8*35/128a2=m2/2+m4/2+m6*15/32+m8*7/16a4=m4/8+m6*3/16+m8*7/32a6
