geopandas：如何“規范化”地理序列

Question

我有一個帶有一些多邊形的GeoDataFrame（來自某些County形狀文件）。

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,...

它們是相同的狀態，曲線圖也相同。 我是shapefile的新手，但是認為如果可以“標准化”這些多邊形（將x維度中的所有點除以它們的最小值，與y維度相同），那么兩個DataFrame中的多邊形應該是可比的。

我試過了：

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，python 2.7（我知道...該繼續）了，OSX。

Answer 1

您需要將其中一種幾何形狀重新投影或轉換為通用坐標參考系統。

請參閱GeoSeries.to_crs(crs=None, epsg=None)文檔以執行此操作。

第一個GeoDataFrame看起來像EPSG：4326（這是猜測，還有其他地理坐標參考系統）。 第二個可能是UTM或狀態平面投影。