pyshp генерирует случайные точки в шейп-файле

Question

Я хочу сгенерировать 1000 случайных точек в указанном c шейп-файле Табличной области почтового индекса с использованием pyshp. Мой код:

import shapefile

zctashape = shapefile.Reader('C:/mypath/tl_2019_us_zcta510.shp')

shapefile_len = len(zctashape.shapes())

#identify the index for zcta 84049 and designate number of points.
zcta_to_use = '84049'
pointcount = 1000

for i in range(0,shapefile_len):
    if zctashape.record(i)[1] == zcta_to_use:
        record_i=i
        break

record_i

Как мне перейти отсюда? Я прошу прощения, если это основа c, я в основном пользователь R, и на этом языке это очень легко. Шейп-файл загружен с https://www.census.gov/cgi-bin/geo/shapefiles/index.php?year=2019&layergroup=ZIP+Code+Tabulation+Areas

Answer 1

Я использую geopandas, так как я думаю, что это облегчает. Я не знаю, есть ли простой способ достичь этого исключительно с помощью pyshp:

import numpy as np
import random
from shapely.geometry import Point
import geopandas as gpd
import matplotlib.pyplot as plt

filename = 'tl_2019_us_zcta510.shp'
gdf = gpd.read_file(filename)

#identify the index for zcta 84049 and designate number of points.
zcta_to_use = '84049'
pointcount = 100

aoi = gdf[gdf['ZCTA5CE10'] == zcta_to_use]
aoi_geom = aoi.unary_union

# find area bounds
bounds = aoi_geom.bounds
xmin, ymin, xmax, ymax = bounds

xext = xmax - xmin
yext = ymax - ymin

points = []
while len(points) < pointcount:
    # generate a random x and y
    x = xmin + random.random() * xext
    y = ymin + random.random() * yext
    p = Point(x, y)
    if aoi_geom.contains(p):  # check if point is inside geometry
        points.append(p)

gs = gpd.GeoSeries(points)

fig, ax = plt.subplots()

aoi.plot(ax=ax, facecolor='none', edgecolor='steelblue')
gs.plot(ax=ax, color='r')