如何解决在python中有效地为多边形分配点
我有一个与 Assign points to polygons effeciently 几乎相同的问题,但我使用的是 Python。
我已经将一个经纬度列表转换为一个匀称点列表,并且我有一个从 shapefile 中读入的多边形列表:
import shapely
import shapefile
points = [Point(lon,lat) for lon,lat in zip(lons,lats)]
shp = shapefile.Reader('file.shp')
polys = shp.shapes()
records = shp.records()
我想生成一个列表,其中包含每个点的多边形 ID(存储在 records[0]
中)。
到目前为止,我的方法是这样的:
polyIDs = [0] * len(points)
for i in range(len(points)):
count = 0
for j in polys:
if points[i].within(shape(j)):
polyIDs[i] = records[count][0]
count += 1
我不确定这些代码是否有效并且需要很长时间(我有 40,000 个点和 2,000 个多边形),或者是否出现了问题?无论哪种方式,有没有办法改进这个过程?有没有更好的方法来检查一系列点落在几个多边形中的哪个而不检查每个可能的关系?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。