阅读:4665回复:8
谁有WGS84转北京54的程序啊,搞了好几天啦,没成果
谁有WGS84转北京54的程序啊,搞了好几天啦,没成果
|
|
1楼#
发布于:2005-06-03 14:24
<P>我有VB的代码,你自己参考一下吧:</P>
<P>'经纬度BL换算到高斯平面直角坐标XY(高斯投影正算)<BR>Private Function Deg2XY(ByRef F2 As Double, ByRef E2 As Double, _<BR> ByRef s2 As Double, ByRef t2 As Double) As Boolean<BR> 'A2 输入中央子午线,以度.分秒形式输入,如115度30分则输入115.30; 起算数据L0<BR> 'F2 以度小数形式输入经度值, L<BR> 'E2 以度小数形式输入纬度值,B<BR> 'S2 计算结果,横坐标Y,,请特别注意,这里生成的是高斯坐标<BR> 'T2 计算结果,纵坐标X<BR> '投影带号计算 n=[L/6]+1 如:兰州测得经度103.XXXX,故n=[103.X/6]+1=17+1=18<BR> '中央经线经度 L0 = n*6-3 = [L/6]*6+3<BR> <BR> 'ByRef A2 As Double,<BR> Dim A2 As Double<BR> Dim B2 As Double<BR> 'Dim G2 As Double<BR> Dim H2 As Double<BR> Dim i2 As Double<BR> Dim J2 As Double<BR> Dim K2 As Double<BR> Dim l2 As Double<BR> Dim M2 As Double<BR> Dim n2 As Double<BR> Dim O2 As Double<BR> Dim P2 As Double<BR> Dim Q2 As Double<BR> Dim R2 As Double<BR> Dim NN As Integer '投影代号<BR>' A2 = CInt(F2 / 6) * 6 - 3<BR> <BR>' B2 = Int(A2) + (Int(A2 * 100) - Int(A2) * 100) / 60 + (A2 * 10000 - Int(A2 * 100) * 100) / 3600<BR> '把L0化成度(A2)<BR> NN = CInt(F2 / 6) + 1 '投影代号n=[L/6]+1 '<BR> B2 = NN * 6 - 3 '中央经线<BR> 'G2 = F2 - B2 ' L -L0<BR> 'H2 = G2 / 57.2957795130823 '化作弧度<BR> H2 = (F2 - B2) / 57.2957795130823 '将经差的单位化为弧度<BR> i2 = Tan(E2 / 57.2957795130823) 'Tan (B)<BR> J2 = Cos(E2 / 57.2957795130823) ' Cos (B)<BR> K2 = 0.006738525415 * J2 * J2<BR> l2 = i2 * i2<BR> M2 = 1 + K2<BR> n2 = 6399698.9018 / Sqr(M2)<BR> O2 = H2 * H2 * J2 * J2<BR> P2 = i2 * J2<BR> Q2 = P2 * P2<BR> R2 = (32005.78006 + Q2 * (133.92133 + Q2 * 0.7031))<BR> s2 = ((((l2 - 18) * l2 - (58 * l2 - 14) * K2 + 5) * O2 / 20 + M2 - l2) * O2 / 6 + 1) * n2 * (H2 * J2)<BR> <BR> '计算结果,横坐标Y,这里生成的是高斯坐标<BR> s2 = s2 + NN * 1000000 + 500000 '18500000 '在计算的基础上加上了“带号”(18)和“东移”(500KM)<BR> '计算结果,纵坐标X<BR> t2 = 6367558.49686 * E2 / 57.29577951308 - P2 * J2 * R2 + ((((l2 - 58) * l2 + 61) * _<BR> O2 / 30 + (4 * K2 + 5) * M2 - l2) * O2 / 12 + 1) * n2 * i2 * O2 / 2<BR> <BR> Deg2XY = True<BR>End Function</P> |
|
|
2楼#
发布于:2005-06-04 09:18
谢谢老大
|
|
3楼#
发布于:2007-07-07 22:31
感谢。。。
|
|
4楼#
发布于:2008-05-12 09:36
<img src="images/post/smile/dvbbs/em02.gif" /><img src="images/post/smile/dvbbs/em01.gif" />
|
|
5楼#
发布于:2009-03-11 17:18
下个COORD软件,看帮助就OK了。
|
|
6楼#
发布于:2009-04-04 01:47
请教一下,WGS84到54的转换不需要转换参数吗?
|
|
7楼#
发布于:2009-04-07 11:58
<P>肯定要转换参数,如中央经线等</P>
|
|
8楼#
发布于:2009-05-26 15:11
<img src="images/post/smile/dvbbs/em02.gif" />shou le a xie xie
|
|