高斯投影坐标正反算VB程序.docx
《高斯投影坐标正反算VB程序.docx》由会员分享,可在线阅读,更多相关《高斯投影坐标正反算VB程序.docx(7页珍藏版)》请在冰豆网上搜索。
![高斯投影坐标正反算VB程序.docx](https://file1.bdocx.com/fileroot1/2022-11/27/969976fb-629b-48df-b6e4-93be1903dce2/969976fb-629b-48df-b6e4-93be1903dce21.gif)
高斯投影坐标正反算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.
注意选择需要的椭球。