24小时热门版块排行榜

>论坛更新日志 (3202)
>虫友互识 (703)
>休闲灌水 (109)
>导师招生 (104)
>博后之家 (92)
>基金申请 (92)
>论文投稿 (87)
>文献求助 (77)
>教师之家 (47)
>硕博家园 (43)
>考博 (25)
>药学 (21)
>招聘信息布告栏 (15)
>医学 (11)
>公派出国 (11)
>考研 (11)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zzahkj

木虫 (著名写手)

龙

应助: 42 (小学生)
金币: 1929.7
散金: 2554
红花: 30
帖子: 1520
在线: 465.5小时
虫号: 964751
注册: 2010-03-08
性别: GG
专业: 地理信息系统

[交流] Arcengine怎么把西安80坐标转换为经纬度？已有1人参与

菜鸟求助，给点代码吧！谢谢！

回复此楼

» 猜你喜欢

曹建劲团队在各种介质中发现含金属微粒，发明微粒找矿技术，获得8项国家发明专利授权已经有11人回复
城乡规划学招收调剂 1名已经有0人回复
地球物理学和空间物理学论文润色/翻译怎么收费? 已经有70人回复
首次在隐伏矿体上方发现金、硫酸铅、硝酸铅、三氧化钨等含金属纳米微粒已经有8人回复
求助：有没有大神有流体包裹体计算软件已经有0人回复
【招聘】科研助理+五险一金已经有20人回复

» 本主题相关价值贴推荐，对您同样有帮助:

80金币求助专业翻译已经有1人回复
Vb.net+ArcEngine生成点状图层的速度太慢了该怎么办？已经有10人回复
【分享】“北京54坐标系”转“西安80坐标系”的转换方法和步骤已经有11人回复
【求助】重金求助vasp 算过渡态 80 金币已经有69人回复
审稿人让加几篇参考文献怎么办已经有29人回复
【求助】discover 里面怎么分析氢键的键角和键长分布已经有12人回复

不要迷恋哥，哥只是一个过客！

1楼 2011-05-28 13:40:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zzahkj

木虫 (著名写手)

龙

应助: 42 (小学生)
金币: 1929.7
散金: 2554
红花: 30
帖子: 1520
在线: 465.5小时
虫号: 964751
注册: 2010-03-08
性别: GG
专业: 地理信息系统

★
yanggis(金币+1): 谢谢分享。。。 2011-05-29 21:47:02

引用回帖:

Originally posted by 昆昆的小虫 at 2011-05-28 14:18:31:
代码不方便共享建议去看下书《应用大地测量学》中国矿业大学出版社

哦主要是我的原图是无带号的即X、Y是6、7为的，另外我的题目也问错了，应该是经纬度转换为西安坐标系，可是我只能把经纬度转换为有带号的即X、Y是8、7位的，我的代码如下能不能给我指点一下怎么转换为6、7位的

CODE:

    Dim ProjNo = 0, ZoneWide As Integer '  带宽 

        Dim longitude1, latitude1, longitude0, latitude0, X0, Y0, xval, yval As Double

        Dim a, f, e2, ee, NN, T, C, M, iPI, Aa As Double

        iPI = 0.0174532925199433 '  3.1415926535898/180.0; 

        ZoneWide = 6   ' 6度带宽 

        a = 6378140.0

        f = 1 / 298.257 '   80年西安坐标系参数 

        ProjNo = (UserInput_Lon / ZoneWide)

        longitude0 = ProjNo * ZoneWide + ZoneWide / 2

        longitude0 = longitude0 * iPI

        latitude0 = 0

        longitude1 = UserInput_Lon * iPI '  经度转换为弧度

        latitude1 = UserInput_Lat * iPI '    纬度转换为弧度

        e2 = 2 * f - f * f

        ee = e2 * (1.0 - e2)

        NN = a / Math.Sqrt(1.0 - e2 * Math.Sin(latitude1) * Math.Sin(latitude1))

        T = Math.Tan(latitude1) * Math.Tan(latitude1)

        C = ee * Math.Cos(latitude1) * Math.Cos(latitude1)

        Aa = (longitude1 - longitude0) * Math.Cos(latitude1)

        M = a * ((1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256) * latitude1 - (3 * e2 / 8 + 3 * e2 * e2 / 32 + 45 * e2 * e2 * e2 / 1024) * Math.Sin(2 * latitude1) + (15 * e2 * e2 / 256 + 45 * e2 * e2 * e2 / 1024) * Math.Sin(4 * latitude1) - (35 * e2 * e2 * e2 / 3072) * Math.Sin(6 * latitude1))

        xval = NN * (Aa + (1 - T + C) * Aa * Aa * Aa / 6 + (5 - 18 * T + T * T + 72 * C - 58 * ee) * Aa * Aa * Aa * Aa * Aa / 120)

        yval = M + NN * Math.Tan(latitude1) * (Aa * Aa / 2 + (5 - T + 9 * C + 4 * C * C) * Aa * Aa * Aa * Aa / 24 + (61 - 58 * T + T * T + 600 * C - 330 * ee) * Aa * Aa * Aa * Aa * Aa * Aa / 720)

        X0 = 1000000L * (ProjNo + 1) + 500000L

        Y0 = 0

        xval = xval + X0

        yval = yval + Y0

        Dim LocationX As Integer = xval

        Dim LocationY As Double = yval