高斯投影坐标正反算VB程序文档格式.docx
《高斯投影坐标正反算VB程序文档格式.docx》由会员分享，可在线阅读，更多相关《高斯投影坐标正反算VB程序文档格式.docx（9页珍藏版）》请在冰点文库上搜索。
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）/206265
b_=（deg2*3600+min2*60+sec2）/206265
dh=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=6Then
l0=6*ik1-3
Else
Ifdh=3Then
l0=3*ik2
MsgBox"
error"
48,"
:
ExitSub
EndIf
l=l_-l0*3.14159/180
e=Sqr（a*a-b*b）/a
m0=a*（1-e*e）
m2=e*e*m0*3/2
m4=e*e*m2*5/4
m6=m4*e*e*7/6
m8=e*e*m6*9/8
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
x_=a0*b_-a2*Sin（2*b_）/2+a4*Sin（4*b_）/4-Sin（6*b_）*a6/6+Sin（8*b_）*a8/8
t=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_））
N=a/W
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/720
y=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
y_=y+500000+1000000*ik1
y_=y+500000+1000000*ik2
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）
m2=e*e*m0*3/2
m4=e*e*m2*5/4
m6=m4*e*e*7/6
m8=e*e*m6*9/8
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
c=a*a/b
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））
l0=6*ik-3
l0=3*ik
l_=l0*3.14159/180+l
sec1=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）
PrivateSubCommand3_Click（）
End
PrivateSubCommand4_Click（）
Text1.Text="
"
Text2.Text="
Text3.Text="
Text4.Text="
Text5.Text="
Text6.Text="
Text7.Text="
Text8.Text="
Text9.Text="
Text10.Text="
PrivateSubOption1_Click（）
a=6378245
b=6356863.01877305
PrivateSubOption2_Click（）
a=6378140
b=6356755.28815753
PrivateSubOption3_Click（）
a=6378137
b=6356752.3142
PrivateSubOption4_Click（）
b=6356752.3141
PrivateSubOption5_Click（）
a=Val（InputBox（"
a"
"
plsaseinput"
））
b=Val（InputBox（"
b"
pleaseinput"
界面截图如下：
使用方法：
1.如下所示，高斯投影坐标正算时，在经度纬度相应的文本框里输入经纬度，再在度号对应的文本框里输入是六度带还是三度带，最后按坐标正算按钮即可，答案会显示在X，Y相应的文本框里；
2.如果要进行坐标反算，则输入X，Y与度号，最后按坐标反算按钮即可得到所需的大地经纬度；
3.注意选择需要的椭球。
升级会员 