1

我最近才开始使用 shapefile。我有一个 shapefile，其中每个对象都是一个多边形。我想生成一个新的 shapefile，其中每个多边形的几何形状都被它的质心替换。有我的代码。

import geopandas as gp
from shapely.wkt import loads as load_wkt

fname = '../data_raw/bg501c_starazagora.shp'
outfile = 'try.shp'
shp = gp.GeoDataFrame.from_file(fname)

centroids = list()
index = list()

df = gp.GeoDataFrame()

for i,r in shp.iterrows():
    index.append(i)
    centroid = load_wkt(str(r['geometry'])).centroid.wkt
    centroids.append(centroid)

df['geometry'] = centroids
df['INDEX'] = index

gp.GeoDataFrame.to_file(df,outfile)

当我运行脚本时，我最终raise ValueError("Geometry column cannot contain mutiple " ValueError: Geometry column cannot contain mutiple geometry types when writing to file. 无法理解出了什么问题。有什么帮助吗？

4

1 回答 1

1

问题是您正在使用几何的字符串表示填充几何字段，而不是形状匀称的几何对象。

无需转换为 wkt。您的循环可能是：

for i,r in shp.iterrows():
    index.append(i)
    centroid = r['geometry'].centroid
    centroids.append(centroid)

但是，根本不需要遍历地理数据框。您可以创建一个新的 shapefile 质心，如下所示：

df=gp.GeoDataFrame(data=shp, geometry=shp['geometry'].centroid)
df.to_file(outfile)
于 2016-12-13T17:52:48.510 回答