地质公园:如何“规范化”地质学

Question

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

0       POLYGON ((-118.71455 46.822508, -118.67662 46....
1       POLYGON ((-118.663169 47.043673, -118.66334 47...
2       POLYGON ((-118.484262 46.73667, -118.495995 46...
3       POLYGON ((-118.483051 46.839378, -118.482689 4...
4       POLYGON ((-118.82874 46.77891, -118.822177 46....
5       POLYGON ((-118.606476 46.96391999999999, -118....

我试着把这些多边形和另一个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,...

他们是同一个州，情节看起来是一样的。我对shapefiles还不熟悉，但是如果我能够“规范化”这些多边形(用x维中的所有点除以最小值，y维也是一样的)，那么两个DataFrames中的多边形应该是相似的。

我试过：

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

而且还

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

没有任何运气。他们只是返回他们的原始值(或者接近，值应该在1到2之间)。我可以循环这些多边形，并打开他们的点，并划分他们，但我希望有一个更优雅的解决办法。

运行最新的geopandas 0.1.0，python2.7(我know...time继续前进)，OSX的开发版本。

Answer 1

您需要将其中一个几何图形重新规划或转换为一个公共的坐标参考系统。

请参阅GeoSeries.to_crs(...)来完成此操作。

第一个GeoDataFrame看起来像EPSG:4326 (这是一个猜测，还有其他的地理坐标参考系统)。第二种可能是UTM或国家平面投影。