pyproj -

Question

私は緯度と経度するX、Yシステムを変換しようとしています、そして多くの人々が彼らのサンプルコードではpyproj

をお勧めしているようですが、彼らはp1 = Proj(init='epsg:26915')を使用するすべてのUTMゾーンの初期点を決定する方法と、

p1 = Proj(init='epsg:26910') 
p2 = Proj(init='epsg:26710') 
x,y = 492893.53, 5458938.38 
x2,y2 = transform(p1, p2, x, y) 

：p2 = Proj(init='epsg:26715')

しかし、それはここではUTM Zone15 のためですが、私のサンプルのx = 492893.53、Y = 5458938.38 とUTM Zone10Nのための私のサンプルコードです出力X2、Y2は、私がUTMゾーン10NのためのEPSG値をチェックしてみました右の緯度と経度

はありませんが、オンラインの3つのバージョンがあります：私はそれらのすべてを試したversion 1、version 2、version 3 は、それらのどれも動作しません。

p1、p2は実際に何を意味し、どのUTMゾーンのp1、p2の値を見つけることができますか？

Answer 1

lat-longに変換する場合は、lat-long座標系の投影が必要です。ほとんどの場合、これはWGS84 LAT-ロングEPSGコードであるGPS、4326

import pyproj 
ll = pyproj.Proj(init="epsg:4326") 

で使用されるように、あなたはあなたのデータポイントからの投影を必要としています。この場合、あなたはUTM10と言っていますが、あなたはどちらのバリエーションを言っているのかは分かりませんが、座標から告げる方法はありません。最初のは、WGS84バージョン試すことができます：

>>> utm10 = pyproj.Proj(init="epsg:32610") 
>>> pyproj.transform(utm10, ll, 492893.53, 5458938.38) 
(-123.0977143330712, 49.28315708900033) 

をこれがUTM10の有効領域のために（実際には（-123、49）が長い緯度である）可能LAT-長いように見えます。だから私はあなたの質問に「彼らのどれも働かない」ことが何を意味するのか分かりません。

他の2つの投影コードは、ほぼ同じ結果をもたらします。あなたが使用しているのはどれですかは数字だけで決めることはできません。数字を入力した人にはどの投影があるのか​​尋ねる必要があります。

x座標とy座標を指定して、 xとyの数字は各ゾーンの経度バンドの各ゾーンで重複しているため、地球上のどのゾーンでもかまいません。