VBϸß˹×ø±ê±ä»»µÄʵÏÖ
ÔøÊ¥  ³Âΰ

¸ß˹-¿ËÂÀ¸ñ(Gauss-Kruger)ͶӰ£¬ÊÇÒ»ÖÖ¡°µÈ½ÇºáÇÐÔ²ÖùͶӰ¡±¡£µÂ¹úÊýѧ¼Ò¡¢ÎïÀíѧ¼Ò¡¢ÌìÎÄѧ¼Ò¸ß˹£¨Carl Friedrich £Çauss£¬1777Ò» 1855£©ÓÚÊ®¾ÅÊÀ¼Í¶þÊ®Äê´úÄⶨ£¬ºó¾­µÂ¹ú´óµØ²âÁ¿Ñ§¼Ò¿ËÂÀ¸ñ£¨Johannes Kruger£¬1857¡«1928£©ÓÚ 1912Äê¶ÔͶӰ¹«Ê½¼ÓÒÔ²¹³ä£¬¹ÊÃû¡£ÉèÏëÓÃÒ»¸öÔ²ÖùºáÇÐÓÚÇòÃæÉÏͶӰ´øµÄÖÐÑë¾­Ïߣ¬°´ÕÕͶӰ´øÖÐÑë¾­ÏßͶӰΪֱÏßÇÒ³¤¶È²»±äºÍ³àµÀͶӰΪֱÏßµÄÌõ¼þ£¬½«ÖÐÑë¾­ÏßÁ½²àÒ»¶¨¾­²î·¶Î§ÄÚµÄÇòÃæÕýÐÎͶӰÓÚÔ²ÖùÃæ.
¸ß˹-¿ËÂÀ¸ñ(Gauss-Kruger)ͶӰÊǺáÖáÄ«¿¨ÍÐͶӰµÄ±äÖÖ£¬¸ß˹-¿ËÂÀ¸ñͶӰÊÇ¡°µÈ½ÇºáÇÐÔ²ÖùͶӰ¡±£¬Í¶Ó°ºóÖÐÑë¾­Ïß±£³Ö³¤¶È²»±ä£¬¼´±ÈÀýϵÊýΪ1¡£¸ß˹-¿ËÂÀ¸ñͶӰµÄ¶«Î±Æ«ÒÆÊÇ500¹«ÀͶӰ±±Î±Æ«ÒÆΪÁã¡£
¸ß˹-¿ËÂÀ¸ñͶӰÕý½â¹«Ê½£ºÒÑÖª£¨B,L£©ÖµÇó½â(X,Y)£¬£¨Ô­µãγ¶È 0£¬ÖÐÑë¾­¶ÈL0£©

ÉÏÃ湫ʽÖж«Î³Æ«ÒÆFE = 500000Ã× £¨±¾³ÌÐòÖÐÖÐÖ»Éè¼Æ¼Ó500000Ã׵ij£Êý£¬È绹Ҫ´øºÅµÄ»°FE = 500000+ ´øºÅ * 1000000Ã×£©£»¸ß˹-¿ËÂÀ¸ñͶӰ±ÈÀýÒò×Ók0 = 1
¸ß˹-¿ËÂÀ¸ñͶӰ·´½â¹«Ê½£ºÒÑÖª(X,Y)Çó½â£¨B,L£©£¬£¨Ô­µãγ¶È 0£¬ÖÐÑë¾­¶ÈL0£©       

a -- ÍÖÇòÌ峤°ëÖá
b -- ÍÖÇòÌå¶Ì°ëÖá
f -- ±âÂÊ
e -- µÚһƫÐÄÂÊ
e¡¯ -- µÚ¶þÆ«ÐÄÂÊ 
N --  îÓÏȦÇúÂʰ뾶 
R -- ×ÓÎçȦÇúÂʰ뾶
B -- γ¶È£¬L -- ¾­¶È£¬µ¥Î»»¡¶È(RAD)
X  -- ×ÝÖ±½Ç×ø±ê, Y -- ºáÖ±½Ç×ø±ê£¬µ¥Î»Ã×(M)
¹ØÓÚÍÖÇòÌå²ÎÊý£¬ÎÒ¹ú³£ÓõÄ3¸öÍÖÇòÌå²ÎÊýÈçÏ£º

³ÌÐòʵÏÖ:
'¸ß˹ÕýËãÇóX
Public Function X(ByVal B#, ByVal L#, ByVal L0#) As Double
Dim n#, T#, T2#, m#, m2#, ng2#
Dim S#, C#
X = A1 * B + A2 * Sin(2 * B) + A3 * Sin(4 * B) + A4 * Sin(6 * B)? '×ÓÎçÏß»¡³¤
S = Sin(B)
C = Cos(B)
T = Tan(B)
T2 = T * T
n = a / Sqr(1 - e12 * S * S)'îÓÏȦÇúÂʰ뾶
m = C * (L - L0)
m2 = m * m
ng2 = C * C * e12 / (1 - e12)????
X = X + n * T * ((0.5 + ((5 - T2 + 9 * ng2 + 4 * ng2 * ng2) / 24# + (61 - 58 * T2 + T2 * T2) * m2 / 720#) * m2) * m2)?
End Function
'¸ß˹ÕýËãÇóY
Public Function Y(ByVal B#, ByVal L#, ByVal L0#) As Double
Dim n#, T#, T2#, m#, m2#, ng2#
Dim S#, C#
S = Sin(B)
C = Cos(B)
T = Tan(B)
T2 = T * T
n = a / Sqr(1 - e12 * S * S)
m = C * (L - L0)
m2 = m * m
ng2 = C * C * e12 / (1 - e12)
Y = n * m * (1 + m2 * ((1 - T2 + ng2) / 6# + m2 * (5 - 18 * T2 + T2 * T2 + 14 * ng2 - 58 * ng2 * T2) / 120#))
Y = Y + Y0
End Function
'¸ß˹·´ËãÇóB£¨Î³¶È)
Public Function B(ByVal X#, ByVal Y#) As Double
Dim S#, C#, T#, T2#, n#, ng2#, V#, yN#
Dim preB0#, B0#
Dim eta#
Y = Y - Y0
B0 = X / A1
Do
preB0 = B0
B0 = (X - (A2 * Sin(2 * B0) + A3 * Sin(4 * B0) + A4 * Sin(6 * B0))) / A1
If Abs(B0 - preB0) < 0.000000001 Then Exit Do
Loop
S = Sin(B0)
C = Cos(B0)
T = Tan(B0)
T2 = T * T
n = a / Sqr(1 - e12 * S * S)
ng2 = C * C * e12 / (1 - e12)
V = Sqr(1 + ng2)
yN = Y / n
B = B0 - (yN * yN - (5 + 3 * T2 + ng2 - 9 * ng2 * T2) * yN * yN * yN * yN / 12# + (61 + 90 * T2 + 45 * T2 * T2) * yN * yN * yN * yN * yN * yN / 360#) * V * V * T / 2#
End Function

³ÌÐòÔËÐнçÃæ:

¼ì²âÊý¾Ý:
£¨1£©¡¢ÒÑÖªÔÚ±±¾©×ø±êϵÏÂ
ÖÐÑë×ÓÎçÏߣº117¶È£¬Î³¶ÈB£º28¶È32·Ö14.5Ã룬¾­¶ÈL£º116¶È54·Ö12.3Ãë
ÕýËãÇó½â±±·½ÏòX¡¢¶«·½ÏòYÖµ¡£µÃ:
X= 3158054.118Ã×
Y= 490547.265Ã×
£¨2£©¡¢ÒÑÖªÔÚ±±¾©×ø±êϵÏÂ
ÖÐÑë×ÓÎçÏߣº114¶È£¬±±·½ÏòXÖµ£º2834561.381Ãס¢¶«·½ÏòYÖµ£º583726.735Ã×
·´ËãÇó½âB¡¢LÖµ¡£µÃ:
B= 25¶È36·Ö56.1Ãë
L= 114¶È50·Ö00.8Ãë
£¨3£©¡¢ÒÑÖªÔÚ±±¾©×ø±êϵÏÂ
»»´øÇ°ÖÐÑë×ÓÎçÏߣº117¶È£¬»»´øºóÖÐÑë×ÓÎçÏߣº116¶È20·Ö£¬±±·½ÏòXÖµ£º 3158038.731Ã×£¬¶«·½ÏòYÖµ£º490539.097Ã×
»»´ø¼ÆËãÖÐÑë×ÓÎçÏß116¶È20·ÖÏÂеÄX¡¢YÖµ¡£µÃ
X= 3158167.489Ã×
Y= 555787.193Ã×

½¨ÒéʹÓÃ1024*768·Ö±æÂÊ IE5.5ÒÔÉÏ°æ±¾ä¯ÀÀÆ÷
°æȨËùÓР ®½­Î÷ʡúÌïµØÖʾֲâ»æ´ó¶Ó  ±¸°¸ºÅ:¸ÓICP±¸05009195
E-mail£ºjxsmtch@vip.sina.com