Questions tagged «pyproj»

3
变换形状多边形和MultiPolygon对象
是否有一种简单的方法可以将Shapely对象(即Polygons和MultiPolygons)从一个投影转换为另一个投影,而无需手动挖掘和提取坐标? 实际上,我现在甚至都不在乎它们是否是Shapely对象,我只想传递要素和投影,并获得重新投影的要素集。 这种功能是否存在,或者必须手动编码?

2
使用geopandas获取多边形区域?
给定一个geopandas GeoDataFrame包含一系列多边形的区域,我想获取列表中每个要素的面积(平方公里)。 这是一个很常见的问题,并且在过去通常建议的解决方案是使用shapely和pyproj直接(如这里和这里)。 有没有办法做到这一点geopandas?

1
使用Python从投影坐标转换​​shapefile
新手在这里苦苦于GIS。我正在尝试使用在其县网站(县)网站上找到的shapefile来为密尔瓦基市(Milwuakee)绘制地图。我在这里取得了一些成功。我的代码是给: 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 问题是这是错误的。密尔沃基的lon / lat为-87.863984和42.920816。 其次,我如何以编程方式对整个shapefile执行此操作。我想将其绘制到底图中。当我尝试遵循此线程时,出现错误代码: with fiona.open("ward2012/ward.shp") as shp: ori = Proj(init='epsg:32054' ), dest= Proj(init='EPSG:4326',preserve_units=True) …
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.