我有一个列表,其中包含我的所有要点:[(lon, lat), ...]
(红色,绿色和蓝色。)现在我想检查哪些点是我的多边形的一部分,我的结果集应包含所有蓝色和所有绿色节点。
我发现http://geospatialpython.com/2011/08/point-in-polygon-2-on-line.html,形状和matplotlib。第一个根本不工作。 Shapely不会返回我需要的所有节点。 Matplotlib似乎也考虑了外部边界上的点。
获得所有(蓝色和绿色)节点最便宜的方式是什么?
答案 0 :(得分:0)
如果您可以与C ++接口,我建议添加Luke Simonson和Gyuszi Suto的BOOST库。
论文:http://www.boost.org/doc/libs/1_46_1/libs/polygon/doc/GTL_boostcon2009.pdf
图书馆:http://www.boost.org/doc/libs/1_55_0/libs/polygon/doc/index.htm
基于作者开发的优雅的多边形微积分,这些C ++的加法以"荒谬的速度"运行。