欧美bbbwbbbw肥妇,免费乱码人妻系列日韩,一级黄片

Python的經(jīng)緯度與xy坐標(biāo)系相互轉(zhuǎn)換方式

 更新時(shí)間：2024年02月01日 09:43:42   作者：kalani呀  
這篇文章主要介紹了Python的經(jīng)緯度與xy坐標(biāo)系相互轉(zhuǎn)換方式,具有很好的參考價(jià)值,希望對(duì)大家有所幫助,如有錯(cuò)誤或未考慮完全的地方,望不吝賜教

經(jīng)緯度與xy坐標(biāo)系相互轉(zhuǎn)換

代碼是已經(jīng)運(yùn)行過(guò)的，可以直接輸入?yún)?shù)使用

注意事項(xiàng)

  • 經(jīng)緯度參數(shù)不可以寫(xiě)顛倒，否則會(huì)由于math.log參數(shù)不能為負(fù)數(shù)和0的原因報(bào)錯(cuò)
  • 由于地球半徑的不確定性，所以(x,y)與(lon, lat)相互轉(zhuǎn)換會(huì)有誤差

奉上代碼及運(yùn)行結(jié)果

  • 經(jīng)緯度轉(zhuǎn)換成xy坐標(biāo)

實(shí)例經(jīng)維度參數(shù) (114.5456282282352,36.631262731204049)

def millerToXY (lon, lat):
    xy_coordinate = []
    #地球周長(zhǎng)
    L = 6381372*math.pi*2 
    #平面展開(kāi)，將周長(zhǎng)視為X軸
    W = L 
    #Y軸約等于周長(zhǎng)一般
    H = L/2 
    #米勒投影中的一個(gè)常數(shù)，范圍大約在正負(fù)2.3之間
    mill = 2.3 
    #將經(jīng)度從度數(shù)轉(zhuǎn)換為弧度
    x = lon*math.pi/180 
    # 將緯度從度數(shù)轉(zhuǎn)換為弧度
    y = lat*math.pi/180 
    #這里是米勒投影的轉(zhuǎn)換
    y = 1.25*math.log(math.tan(0.25*math.pi+0.4*y)) 
    # 這里將弧度轉(zhuǎn)為實(shí)際距離 ，轉(zhuǎn)換結(jié)果的單位是公里
    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坐標(biāo)轉(zhuǎn)換成經(jīng)緯度
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 最終需要確認(rèn)經(jīng)緯度保留小數(shù)點(diǎn)后幾位
    lonlat_coordinate.append((round(lon,15),round(lat,15)))
    return lonlat_coordinate
  • 調(diào)用兩個(gè)函數(shù)
import math
print(millerToXY(114.5456282282352,36.631262731204049))
print(xy_to_coor(32805300, 7107399))

打印結(jié)果

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

總結(jié)

以上為個(gè)人經(jīng)驗(yàn)，希望能給大家一個(gè)參考，也希望大家多多支持腳本之家。

相關(guān)文章

最新評(píng)論