利用.shp文件判断一个坐标是否在某国家或者是洲内

最近在把一个Matlab的toolbox转到Python语言，Matlab中一个inpolygon函数难住了我。通过谷歌和百度都没有查询到很有效的方法。

因为我需要判断一个坐标是否在.shp file读取的经纬度（这个经纬度画出来的范围实际上是世界地图上面的一个洲）内，但是Python的很多包，比如：

1. 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]))  #判断(2.5, 2.5)是否在五角星内

图片显示如下：

但是输出结果却为False. 后来我发现是由于中间的五边形会被Py