从文本文件获取坐标行并检查点是否在多边形内

w

我在文本文件中以以下顺序(xy)拥有数百个坐标(控制点):

5480.000 -4880.000
5480.000 -4860.000
5480.000 -4840.000
5480.000 -4820.000
5480.000 -4800.000
5480.000 -4780.000
5480.000 -4760.000
5480.000 -4740.000
5480.000 -4720.000
5480.000 -4700.000

和多边形的坐标(12):

259.59 229.87
329.07 284.19
397.76 262.62
480.03 397.60
336.26 545.37
325.88 400.00
142.17 440.74
92.65 275.40
158.95 265.02
180.51 196.33
265.97 119.65
259.59 229.87

我要做的是检查控制点是否在多边形内。通常,如果我的脚本中只有几个坐标,就可以这样做,但是我不知道如何读取这样的元组。我有一个暂时不起作用的脚本。

它返回

File "C:/......../test.py", line 16, in <module> xc.append(row[0]) IndexError: list index out of range

最好的解决方案是在终端内是否出现多边形内部的坐标。

from shapely.geometry import Point, Polygon
import csv

poly = open(r"polygon.txt","r")
control = open(r"controlPoints.txt","r")

xp, yp = [], []
for l in poly:
    row = l.split()
    xp.append(row[0])
    yp.append(row[1])

xc,yc = [], []
for l in control:
    row = l.split()
    xc.append(row[0])
    yc.append(row[1])
print(xc)
print(yc)

# Create Point objects

p = Point(xc, yc)

# Create a Polygon 

coords = [xp, yp]
poly1 = Polygon(coords)
print(poly1)

# CHECK 

p.within(poly1)
print(p)

check=poly1.contains(p)
print(check)
莫里斯·迈耶(Maurice Meyer)

您可以结合使用列表理解和map()来创建形状所需的嵌套列表:

poly = Polygon([tuple(map(float, row.split())) for row in open(r"polygon.txt","r") if len(row.split()) == 2])
rawPoints = [tuple(map(float, row.split())) for row in open(r"controlPoints.txt","r") if len(row.split()) == 2]

for pt in rawPoints:
    p = Point(*pt)
    print (p, p.within(poly), poly.contains(p))

打印:

POINT (5480 -4880) False False
POINT (5480 -4860) False False
POINT (5480 -4840) False False
POINT (5480 -4820) False False
POINT (5480 -4800) False False
POINT (5480 -4780) False False
POINT (5480 -4760) False False
POINT (5480 -4740) False False
POINT (5480 -4720) False False
POINT (5480 -4700) False False

本文收集自互联网,转载请注明来源。

如有侵权,请联系 [email protected] 删除。

编辑于
0

我来说两句

0 条评论
登录 后参与评论

相关文章