タグ付けされた質問 「pyproj」

Shapelyオブジェクト（つまり、PolygonとMultiPolygon）をある投影から別の投影に変換する簡単な方法はありますか？ 実際、この時点でそれらがShapelyオブジェクトであるかどうかは気にしません。フィーチャと投影を渡し、再投影されたフィーチャセットを取得するだけです。 この種の機能は存在しますか、それとも手作業でコーディングする必要がありますか？

26 python  coordinate-system  shapely  pyproj 

与えられたgeopandas GeoDataFrameポリゴンのシリーズを含む、私は私のリストにある各機能のキロ平方の領域を取得したいと思います。 これはかなり一般的な問題であり、過去に通常の提案された解決策を使用するようになっているshapelyとpyproj（例えば、直接こことここ）。 これを純粋に行う方法はありgeopandasますか？

16 python  area  shapely  geopandas  pyproj 

初心者はここでGISと格闘しています。私は、郡のウェブサイトの郡のウェブサイトにあるシェープファイルを使用して、ミルウォーキー市の病棟を作成しようとしています。私はいくつかの成功を収めてここのスレッドをフォローしています。私のコードは与えます： from pyproj import Proj, transform # wisconsing EPSG:32054 # epsg:4326 is for the entire world, wgs 84...not obvious inProj = Proj(init='epsg:32054') outProj = Proj(init='epsg:4326') x1,y1 = 2560131.496875003, 406816.434375003 x2,y2 = transform(inProj,outProj,x1,y1) print(x2,y2) 出力付き、 -65.70220967836329 43.08590211722421 問題はこれが間違っていることです。ミルウォーキーの経度/緯度は-87.863984および42.920816です。 次に、シェープファイル全体に対してプログラムでこれを行うにはどうすればよいですか。これをベースマップにプロットしたいと思います。このスレッドをフォローしようとすると、エラーコードが表示されます： with fiona.open("ward2012/ward.shp") as shp: ori = Proj(init='epsg:32054' ), dest= Proj(init='EPSG:4326',preserve_units=True) with fiona.open('ward2012/MKE_wards_lat_lon.shp', …

10 python  shapefile  fiona  pyproj 

方位と範囲を指定したポイントから目的地の位置を（Pythonで）計算する方法を探しています。 subject（geopyとpyproj）の2つのライブラリからの結果の比較に基づいて、最終的な出力の違いが大きくなっていることに気付きました。たとえば、100 kmでは、およそデシメートルのオーダーです。これは私が意味することの最小限の例です： from __future__ import absolute_import, division, print_function long_1 = -1.729722 lat_1 = 53.320556 bearing = 96.021667 distance = 124.8 #in km # using geopy import geopy from geopy.distance import VincentyDistance origin = geopy.Point(lat_1, long_1) destination = VincentyDistance(kilometers=distance).destination(origin, bearing) gp_lat_2 = destination.latitude gp_long_2 = destination.longitude # using pyproj from …

8 python  geodesy  pyproj  geopy  vincenty-formulae