Python 使用pysal映射坐标的位置

我正在尝试绘制我所拥有的一组坐标所在的县。我希望代码遍历每个坐标，检查它所在的县，然后打印相应的县。但是，当前在空闲状态下运行时，光标只是闪烁，似乎什么也没有发生。对于我的代码可能出错的任何帮助，我们将不胜感激。我不确定代码是否正确地检查了县，或者是否可以从shapefile中获取县的名称

import pysal as ps
from ast import literal_eval

COUNTIES = '/Users/JoshuaHawley/Desktop/counties/English Ceremonial Counties.shp'
shps = ps.open(COUNTIES)

counties = shps.read()

f = open('/Users/JoshuaHawley/sundayCoordinates.txt').read()
seq = literal_eval(f.replace("\n", ","))

coordinates = seq[0:]
for line in coordinates:
    points = coordinates[:2]

NotFound = True

while NotFound == True:
    if coordinates in counties:
        NotFound = False
        name_of_county = counties.getNamePolygon()
        print(name_of_county)

---

此代码不正确。Countries是多边形列表，不支持contains（即Countries中的坐标）。您需要首先创建一个定位器，假设县是多边形，那么您将使用PolygonLocator（）getNamePolygon也不是一个真正的函数，您需要创建自己的逻辑来根据多边形在shapefile中的偏移在DBF中查找多边形的名称。