最近在把一个Matlab的toolbox转到Python语言,Matlab中一个inpolygon函数难住了我。通过谷歌和百度都没有查询到很有效的方法。
因为我需要判断一个坐标是否在.shp file读取的经纬度(这个经纬度画出来的范围实际上是世界地图上面的一个洲)内,但是Python的很多包,比如:
- matplotlib 里面包含了Path这个类,可以通过Path产生一个图形,下面是一个简单的例子:
route = Path([[1,0],[0,1],[1,1],[1,0]])
plt.plot([1,0,1,1],[0,1,1,0])
plt.show()
print(route.contains_point([0.8,0.6])) # 判断(0.8, 0.6)是否在三角形内
图片显示是这样的:
输出结果是True
. 但是,一旦把图片的形状改为,五角形状,代码如下:
route = Path([[1,3],[4,3],[1.5,1.5],[2.5,4],[3.5,1.5],[1,3]])
plt.plot([1,4,1.5,2.5,3.5,1],[3,3,1.5,4,1.5,3])
plt.show()
print(route.contains_point([2.5,2.5]))