Python地理座標系和投影座標系轉換

2021-08-21 16:15:28 字數 2441 閱讀 6393

地心地固座標系(earth-centered, earth-fixed,ecef),簡稱地心座標系。

地理座標系統(geographic coordinate system,gcs)1

,座標系是地心座標系,用經緯度表示球面上的點。

世界大地測量系統(world geodetic system, wgs),比如wgs84,是一種地理座標系統,用於全球定位系統(gps)。

投影座標系統(projection coordinate system,pcs)2

通用橫軸墨卡托投影(universal transverse mercator,utm),是一種投影方法。

關於地理座標系和投影座標系更詳細的解釋可以檢視這篇文章:你必須知道的地理座標系和投影座標系

地理座標系統有不同的基準和方法,比如:xian_1980,beijing_1954,wgs_1984等。投影座標系統也有不同的基準和方法,utm和ups等。每乙個地理座標系統(gcs)和投影座標系統(pcs)都有乙個獨特的epsg**,**可在 epsg **查詢。

有一篇介紹pyproj進行地理投影座標系轉換的文章3

,但不夠全面。其中提到arcgis**上查詢 地理座標系 和 投影座標系 的方法很實用但不全。

整理使用python的第三方庫 pypro4

轉換經緯度表示的地理座標系統到投影座標系統。

pypro模組共有兩個函式:

函式描述

test()

執行模組測試

transform(p1, p2, x, y, z=none, radians=false)

用法:x2, y2, z2 = transform(p1, p2, x1, y1, z1, radians=false),將在座標系統p1下的點(x1, y1, z1)轉換到p2座標系統下

轉換經緯度到的投影座標系統;轉換乙個投影座標系統到另乙個投影座標系統;反向轉換,把投影座標系統上的點轉換到地理座標系統:

>>> p1 = pyproj.proj(init='epsg:26915') # 乙個投影座標系統epsg code

>>> p2 = pyproj.proj(init='epsg:26715') # 另乙個投影座標系統epsg code

>>> x1, y1 = p1(-92.199881,38.56694) # 投影到epsg code為26915的投影座標系統

>>>> '%9.3f %11.3f' % (x1,y1)

'569704.566 4269024.671'

>>> x2, y2 = pyproj.transform(p1,p2,x1,y1) # 轉換乙個投影座標系統到另乙個投影座標系統

>>> '%9.3f %11.3f' % (x2,y2)

'569722.342 4268814.027'

>>> '%8.3f %5.3f' % p2(x2,y2,inverse=true) # 反向轉換

' -92.200 38.567'

用元組傳入多個點

>>> lats = (38.83,39.32,38.75) # 所有緯度組成的元組

>>> lons = (-92.22,-94.72,-90.37) # 所有精度組成的元組

>>> x1, y1 = p1(lons,lats) # 轉換經緯度到投影座標系統

>>> x2, y2 = pyproj.transform(p1,p2,x1,y1) # 轉換乙個投影座標系統到另乙個投影座標系統

>>> lons, lats = p2(x2,y2,inverse=true) # 反向轉換

除了使用epsg code之外,還可以顯示指定座標系統名稱

>>> p1 = pyproj.proj(proj='latlong',datum='wgs84') # wgs84,gps使用的地理座標系統,epsg code為4326

>>> x1 = -111.5; y1 = 45.25919444444

>>> p2 = pyproj.proj(proj="utm",zone=10,datum='nad27') # 投影座標系統nad27 / utm zone 10n,epsg code為26710

>>> x2, y2 = pyproj.transform(p1, p2, x1, y1)

>>> "%s %s" % (str(x2)[:9],str(y2)[:9])

'1402285.9 5076292.4'

如果需要計算地球上兩點的距離,可以使用大圓距離來計算,不能使用投影座標系中的兩點計算距離

↩︎#map_projection ↩︎

↩︎ ↩︎

地理座標系和投影座標系

用經緯度表示地面點位的球面座標系。在大地測量學中,對於地理座標系統中的經緯度有三種描述 天文經緯度 大地經緯度 地心經緯度。天文經緯度 表示地面點在大地水準面上的位置,用天文經度和天文緯度表示。大地經緯度 表示地面點在參考橢球面上的位置,用大地精度 大地緯度 大地高h表示 大地座標均以橢球面法線來定...

地理座標系 投影座標系

1.基本概念 平時開展gis開發 研究 應用工作,總會接觸到座標系,也會遇到座標轉換的問題,如地理座標系 投影座標系等。地理座標系是球面座標,參考平面是橢球面,座標單位是經緯度 投影座標系是平面座標系,參考平面是水平面,座標單位是公尺 千公尺等。地理座標系轉換到投影座標系的過程理解為投影,即將不規則...

大地座標系 地理座標系 投影座標系

大地座標系 是大地測量中以參考橢球面為基準面建立起來的座標系。地面點的位置用 大地經度 大地緯度和大地高度 表示。大地座標系的確立包括選擇乙個橢球 對橢球進行定位和確定大地起算資料。乙個形狀 大小和定位 定向都已確定的地球橢球叫參考橢球。參考橢球一旦確定,則標誌著大地座標系已經建立。大地座標系是一種...