Python geopandas:如何;“正常化”;地质系列

Python geopandas:如何;“正常化”;地质系列,python,gis,shapely,geopandas,Python,Gis,Shapely,Geopandas,我有一个带有一些多边形的GeoDataFrame(来自一些郡形状文件) 我正试图在另一个GeoDataFrame(相同状态)中获取这些多边形与一些多边形的交点。然而,多边形的坐标是完全不同的 0 POLYGON ((552346.2856999999 380222.8998000007,... 1 POLYGON ((529754.7249999996 409135.9135999996,... 2 POLYGON ((552333.0957000004 37

我有一个带有一些多边形的GeoDataFrame(来自一些郡形状文件)

我正试图在另一个GeoDataFrame(相同状态)中获取这些多边形与一些多边形的交点。然而,多边形的坐标是完全不同的

0       POLYGON ((552346.2856999999 380222.8998000007,...
1       POLYGON ((529754.7249999996 409135.9135999996,...
2       POLYGON ((552333.0957000004 371726.9315000009,...
3       POLYGON ((542734.0181999998 376530.5831000004,...
4       POLYGON ((552291.6204000004 362052.4251000006,...
他们是同一个州,情节看起来也一样。我不熟悉形状文件,但我想如果我能“规范化”这些多边形(将x维度中的所有点除以它们的最小值,y维度也是如此),那么两个数据帧中的多边形应该是可比较的

我试过:

df.scale(origin=(1,1))
而且

df.scale(xfact=1.0/df.total_bounds[0], yfact=1.0/df.total_bounds[1])
没有任何运气。它们只返回其原始值(或接近,值应介于1和2之间)。我可以在这些多边形中循环,将它们的点解包并进行分割,但我希望有一种更优雅的解决方法


运行Geopandas0.1.0、Python2.7的最新开发版本(我知道……是时候继续了),OSX

您需要将其中一个几何体重新投影或转换为公共坐标系

请务必做到这一点

第一个地理数据框看起来像EPSG:4326(这是一个猜测,还有其他地理坐标参考系)。第二种可能是UTM或状态平面投影

df.scale(xfact=1.0/df.total_bounds[0], yfact=1.0/df.total_bounds[1])