我有一个包含我的所有点的列表：[(lon, lat), ...]（红色，绿色和蓝色的），现在我要检查哪些点是我的多边形的一部分，我的结果集应该包含所有蓝色和所有绿色节点。 我发现http://geospatialpython.com/2011/08/point-in-polygon-2-on-line.html，身材匀称和matplotlib。第一个根本不工作。匀称不返回我需要的所有节点。 Mat
我有一个数据列表如下： from shapely.geometry import box
data = [box(1,2,3,4), box(5,6,7,8), box(1,2,3,4)]
codes = ['A','B','C']
名单“数据”有以下元素： A = box(1,2,3,4)
B = box(5,6,7,8)
C = box(1,2,3,4)
我要检查是否一个元素
我正在处理2D地理数据。我有一长串轮廓路径。现在我想为我的域中的每个点确定它所在的轮廓（即我想计算由轮廓表示的特征的空间频率分布）。 为了说明什么，我想做的事，这里的第一个很天真的实现： import numpy as np
from shapely.geometry import Polygon, Point
def comp_frequency(paths,lonlat):
"
我在Shapely中有一个Polygon和一个MultiLineString。我想扩展没有到达多边形边界的LineString段，以便它确实到达多边形边界。没关系，它可以延伸过去，因为之后我可以轻松地将它夹到边界。 理想情况下，它会以相同的角度继续，但我想这比直接延伸到边界要困难得多。 有没有人有任何建议，我该怎么做呢？ 我产生的几何形状使用下面的代码（如什么，我真的需要做一个简单的例子）： im