基于Python的经纬度与xy坐标系相互转换

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

注意事项：

• 经纬度参数不可以写颠倒，否则会由于math.log参数不能为负数和0的原因报错

• 由于地球半径的不确定性，所以(x,y)与(lon, lat)相互转换会有误差

奉上代码及运行结果：

• 经纬度转换成xy坐标

实例经维度参数 (114.5456282282352,36.631262731204049)

def millerToXY (lon, lat):    xy_coordinate = []    #地球周长    L = 6381372*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 = 6381372 * 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 mathprint(millerToXY(114.5456282282352,36.631262731204049))print(xy_to_coor(32805300, 7107399))

打印结果

[(32805300, 7107399)][(114.54562994611808, 36.63126574363226)]

来源地址：https://blog.csdn.net/m0_58625397/article/details/129747071