批量经纬度转成空间直角坐标系

最新推荐文章于 2022-12-05 18:18:18 发布

兰小莫

最新推荐文章于 2022-12-05 18:18:18 发布

阅读量4.2k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： Python 、GDAL

本文链接：https://blog.youkuaiyun.com/qq_14906811/article/details/73724719

Python 、GDAL 专栏收录该内容

11 篇文章

订阅专栏

本代码示例展示了如何从文本文件中读取经纬度数据，并将其转换为三维空间直角坐标系中的X, Y, Z坐标。该过程考虑到了地球椭球体的特性，如扁率和长半轴等参数。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

直接上代码。。

import numpy
#从键盘自定义输入扁率，长半轴，椭球高（区域中心高）等
#f=1/(298.257)
f=float(input("f:"))
#a=6378140.0
a=float(input("a:"))
L=0.0
e=2*f-f*f
B=0.0
H=float(input("H:"))
f1=open(r'C:\Users\suns\Desktop\POI1.txt','r')
f2=open(r'C:\Users\suns\Desktop\new.txt','w')
#从文本中读取经纬度转成空间直角坐标X,Y,Z
for i in f1:
     st=i.split(',')
     L=float(st[0])
     B=float(st[1])
     n1=e*numpy.math.sin(B/57.29578)
     n2=n1*n1
     n3=1-n2
     n4=numpy.math.sqrt(n3)  
     n5=a/n4
     X=(H+n5)*(numpy.math.cos(B/57.29578)*numpy.math.cos(L/57.29578))
     Y=(H+n5)*(numpy.math.cos(B/57.29578)*numpy.math.sin(L/57.29578))
     Z=(n5*(1-e*e)+H)*(numpy.math.sin(B/57.29578))
     f2.write(str(X)+","+str(Y)+","+str(Z)+"\n")
f1.close()
f2.close()

其实应该要和坐标转换结合起来，有待改进........