创建新 Shapefile 时使用 PyShp 的多部分形状问题

Question

我有一个 shapefile (blockgroups.shp)，我从以下网址下载： https://github.com/GeospatialPython/pyshp/tree/master/shapefiles

我想使用 PyShp 创建一个只有四个属性（bkg_key、pop1990、white 和 black）的新 Shapefile。

我已尝试使用此代码：

import shapefile
file= "blockgroups.shp"
outFile = "newblockgroup3000"
sf = shapefile.Reader(file)
shapeRec= sf.shapeRecords()
w = shapefile.Writer(shapefile.POLYGON)
w.field('BKG_KEY', 'C', 40)
w.field('POP1990', 'N', 40)
w.field('NWHITES', 'N', 40)
w.field('NBLACKS', 'N', 40)
for row in shapeRec:
    bkg_key = row.record[1]
    pop1990 = row.record[2]
    white = row.record[7]
    black = row.record[8]
    points = row.shape.points
    parts = row.shape.parts
    w.parts = parts
    w.poly([points])
    w.record(bkg_key,pop1990,white,black)
    w.save(outFile)

它适用于除一个以外的所有形状。

有一个记录有多个部分。包含多个部分的记录是'BKG_KEY = 060750601001'和'POP = 4531'。在新的 shapefile 中，这条记录的形状很奇怪，因为 pyShp 会自动连接来自特征不同部分的第一个和最后一个顶点。

如果我只选择 'POP1990 4531' 的记录（不包括提到的记录）它可以工作，所以只有当有多个部分的记录时才会出现问题。

当我创建新的 shapefile 时，有什么方法可以保留原始 shapefile 的部分数量吗？我该如何处理这个问题。

我将不胜感激。 谢谢

Answer 1

我在寻找使用 pyshp 将多个形状保存在一个 shapefile 中的解决方案时发现了这个问题。你的代码帮助我解决了我的问题，所以我会为此提供一个解决方案。

import shapefile

file= "blockgroups.shp"
outFile = "newblockgroup3000"
sf = shapefile.Reader(file)
shapeRec= sf.shapeRecords()
w = shapefile.Writer(shapefile.POLYGON)
w.field('BKG_KEY', 'C', 40)
w.field('POP1990', 'N', 40)
w.field('NWHITES', 'N', 40)
w.field('NBLACKS', 'N', 40)

for row in shapeRec:
    bkg_key = row.record[1]
    pop1990 = row.record[2]
    white = row.record[7]
    black = row.record[8]
    points = row.shape.points
    parts = row.shape.parts
    geom = []

    # handle parts here
    for i in range(len(parts)):
        xy = []
        pt = None
        if i < len(parts) - 1:
            pt = points[parts[i]:parts[i + 1]]
        else:
            pt = points[parts[i]:]
        for x, y in pt:
            xy.append([x, y])
        geom.append(xy)

    w.poly(geom)
    w.record(bkg_key,pop1990,white,black)
w.save()

编辑前的结果：

编辑后的结果：