Python的经纬度与xy坐标系相互转换方式

目录

• 经纬度与xy坐标系相互转换
• 注意事项
• 奉上代码及运行结果
• 总结

经纬度与xy坐标系相互转换

代码是已经运行过的，可以直接输入参数使用

注意事项

• 经纬度参数不可以写颠倒，否则会由于math.log参数不能为负数和0的原因报错
• 由于地球半径的不确定性，所以(x,y)与(lon, lat)相互转换会有误差

奉上代码及运行结果

• 经纬度转换成xy坐标

实例经维度参数 (114.52,36.4049)

def millerToXY (lon, lat): xy_coordinate = [] #地球周长 L = *math.pi*2 #平面展开，将周长视为X轴 W = L #Y轴约等于周长一般 H = L/2 #米勒投影中的一个常数，范围大约在正负2.3之间 mill = 2.3 #将经度从度数转换为弧度 x = lon*math.pi/180 # 将纬度从度数转换为弧度 y = lat*math.pi/180 #这里是米勒投影的转换 y = 1.25*math.log(math.tan(0.25*math.pi+0.4*y)) # 这里将弧度转为实际距离 ，转换结果的单位是公里 x = (W/2)+(W/(2*math.pi))*x y = (H/2)-(H/(2*mill))*y xy_coordinate.append((int(round(x)),int(round(y)))) return xy_coordinate

• xy坐标转换成经纬度

def xy_to_coor(x, y): lonlat_coordinate = [] L =  * math.pi*2 W = L H = L/2 mill = 2.3 lat = ((H/2-y)*2*mill)/(1.25*H) lat = ((math.atan(math.exp(lat))-0.25*math.pi)*180)/(0.4*math.pi) lon = (x-W/2)*360/W # TODO 最终需要确认经纬度保留小数点后几位 lonlat_coordinate.append((round(lon,15),round(lat,15))) return lonlat_coordinate

• 调用两个函数

import math print(millerToXY(114.52,36.4049)) print(xy_to_coor(, ))

打印结果

总结

以上为个人经验，希望能给大家一个参考，也希望大家多多支持本网站。

您可能感兴趣的文章:

• python逐像素获取栅格经纬度分别保存在两个矩阵中(代码收藏)
• python库geopy计算多组经纬度距离的实现方式
• Python如何批量处理经纬度数据并生成位置信息
• Python利用Redis计算经纬度距离案例
• Python pyecharts实时画图自定义可视化经纬度热力图