高斯投影坐标正反算VB程序.docx
《高斯投影坐标正反算VB程序.docx》由会员分享，可在线阅读，更多相关《高斯投影坐标正反算VB程序.docx（7页珍藏版）》请在冰豆网上搜索。
高斯投影坐标正反算VB程序
高斯投影坐标正反算
院:
级:
口
号:
名:
课程名称:
指导老师:
实验目的：
1.了解高斯投影坐标正反算的基本思想；
2.学会编写高斯正反算程序，加深了解。
实验原理：
高斯投影正算公式中应满足的三个条件：
1.中央子午线投影后为直线；
2.中央子午线投影后长度不变；
3.投影具有正形性质，即正形投影条件。
高斯投影反算公式中应满足的三个条件：
1.x坐标轴投影成中央子午线，是投影的对称轴；
2.x轴上的长度投影保持不变；
3.正形投影条件，即高斯面上的角度投影到椭球面上后角度没
有变形，仍然相等。
操作工具：
计算机中的VB6.0
代码：
DimaAsDouble,bAsDouble,xAsDouble,yAsDouble,y_#
Diml_AsDouble,b_AsDouble,a0#,a2#,a4#,a6#,a8#,m2#,m4#,m6#,m8#,m0#,l0#,e#,e1#
Dimdeg1AsDouble,min1AsDouble,sec1AsDouble,deg2AsDouble,min2AsDouble,sec2AsDouble
PrivateSubCommand1_Click（）
Dimx_AsDouble,t#,eta#,N#,W#,k1#,k2#,ik1%,ik2%,dh%deg1=Val（Text1.Text）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-3
Else
Ifdh=3Then
l0=3*ik2
Else
MsgBox"error",48,"error":
ExitSub
EndIf
EndIfl=l_-l0*3.14159/180e=Sqr（a*a-b*b）/a
m0=a*（1-e*e）a0=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*Sin（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））eta=Sqr（e1*e1*Cos（b）*Cos（b））
W=Sqr（1-e*e*Sin（b_）*Sin（b_））x=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（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/120
Text7=x
Ifdh=6Theny_=y+500000+1000000*ik1
Else
Ifdh=3Theny_=y+500000+1000000*ik2
Else
MsgBox"error",48,"error":
ExitSub
EndIf
EndIf
Text8=y_
EndSub
PrivateSubCommand2_Click（）
Dimbf#,j%,Wf#,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）
m6=m4*
a0=m0+m2/2+m4*3/8+m6*5/16+m8*35/128
a2=m2/2+m4/2+m6*15/32+m8*7/16
a4=m4/8+m6*3/16+m8*7/32
a6=m6/32+m8/16
a8=m8/128
bf=x/a0
Forj=1To10
bf=（x+a2/2*Sin（2*bf）-a4*Sin（4*bf）/4+a6*Sin（6*bf）/6-a8*Sin（8*bf）/8）/a0
Nextj
e1=Sqr（a*a-b*b）/b
Vf=Sqr（1+e1*e1*Cos（bf）*Cos（bf））
Wf=Sqr（1-e*e*Sin（bf）*Sin（bf））
Nf=a/Wf
Mf=c/Vf八3
tf=Tan（bf）
e1=Sqr（（a*a-b*b）/（b*b））
etaf=Sqr（e1*e1*Cos（bf）*Cos（bf））
ik=Int（y_/1000000）
y=y_-ik*1000000-500000
b_=bf-tf*y*y/（2*Mf*Nf）+tf*（5+3*tf*tf+etaf*etaf-9*etaf*etaf*tf*tf）*y*y*y*y/（24*Mf*Nf*Nf*Nf）+tf*（61+90*tf*tf+45*tf*tf*tf*tf）*y*y*y*y*y*y/（720*Mf*Nf*Nf*Nf*Nf*Nf）
l=y/（Nf*Cos（bf））-（1+2*tf*tf+etaf*etaf）*y*y*y/（6*
Nf*Nf*Nf*Cos（bf））+（5+28*tf*tf+24*tf*tf*tf*tf+6*etaf*etaf+8*etaf*etaf*tf*tf）*y*y*y*y*y/（120*Nf*Nf*Nf*
Nf*Nf*Cos（bf））dh=Val（Text9.Text）
Ifdh=6Thenl0=6*ik-3
Else
Ifdh=3Thenl0=3*ik
Else
MsgBox"error",48,"error":
ExitSub
EndIf
EndIfl_=l0*3.14159/180+lsec1=l_*206265
deg1=Int（sec1/3600）
min1=Int（（sec1-deg1*3600）/60）
sec1=sec1-deg1*3600-min1*60
sec2=b_*206265
deg2=Int（sec2/3600）
min2=Int（（sec2-deg2*3600）/60）
sec2=sec2-deg2*3600-min2*60
Text1=deg1
Text2=min1
Text3=Round（sec1,5）
Text4=deg2
Text5=min2
Text6=Round（sec2,5）
EndSub
PrivateSubCommand3_Click（）
End
EndSub
PrivateSubCommand4_Click（）
Text1.Text=
HH
Text2.Text=
HH
Text3.Text=
HH
Text4.Text=
HH
Text5.Text=
HH
Text6.Text=
HH
Text7.Text=
HH
Text8.Text=
HH
Text9.Text=
HH
Text10.Text=""
EndSub
PrivateSubOption1_Click（）a=6378245b=6356863.01877305
EndSub
PrivateSubOption2_Click（）
a=6378140b=6356755.28815753
EndSub
PrivateSubOption3_Click（）a=6378137b=6356752.3142
EndSub
PrivateSubOption4_Click（）a=6378137b=6356752.3141
EndSub
PrivateSubOption5_Click（）a=Val（InputBox（"a","plsaseinput"））b=Val（InputBox（"b","pleaseinput"））
EndSub
界面截图如下
经度L：
X:
rri
纬度肌
Y：
m.
度号;
清选择椭球：
r克拉索夬
C船S-34
BLTZOiOO
r自定文
坐标正®
坐标反尊
1.
使用方法:
如下所示，高斯投影坐标正算时，在经度纬度相应的文本框里输
入经纬度，再在度号对应的文本框里输入是六度带还是三度带,
2.
最后按坐标正算按钮即可，答案会显示在X，丫相应的文本框里;
如果要进行坐标反算，则输入X，Y与度号，最后按坐标反算按
钮即可得到所需的大地经纬度;
3.
注意选择需要的椭球。
升级会员 