擅长:python、mysql、java
<p>我看到您回答了您的问题(并且似乎对使用交集很满意),但也注意到<code>shapely</code>(以及底层的<code>geos</code>库)已经准备了<em>几何体</em>,以便对某些谓词(包含、正确地包含、覆盖和交集)进行更有效的批处理操作。
见<a href="https://shapely.readthedocs.io/en/stable/manual.html#prepared-geometry-operations" rel="nofollow noreferrer">Prepared geometry operations</a></p>
<p>根据您问题中的代码改编,可以这样使用:</p>
<pre><code>from shapely.prepared import prep
# determine maximum edges
polygon = shape(geojson['features'][i]['geometry'])
latmin, lonmin, latmax, lonmax = polygon.bounds
# create prepared polygon
prep_polygon = prep(polygon)
# construct a rectangular mesh
points = []
for lat in np.arange(latmin, latmax, resolution):
for lon in np.arange(lonmin, lonmax, resolution):
points.append(Point((round(lat,4), round(lon,4))))
# validate if each point falls inside shape using
# the prepared polygon
valid_points.extend(filter(prep_polygon.contains, points))
</code></pre>