いくつかの(必ずしも凸ではない)形状の境界を定義する点の大きなリストがあります。次に、クエリポイントがあり、ポイントの境界で定義された領域内にある(x, y)
かどうかを確認したいと思います。(x, y)
だから、十分に単純な質問。クエリポイントが境界ポイントによって形成される形状の内側にあるかどうかを判断するにはどうすればよいですか?そして、これのための素晴らしいブーストモジュールはありますか?調べてboost::geometry
いますが、まだ何も見つかりません。
いくつかの(必ずしも凸ではない)形状の境界を定義する点の大きなリストがあります。次に、クエリポイントがあり、ポイントの境界で定義された領域内にある(x, y)
かどうかを確認したいと思います。(x, y)
だから、十分に単純な質問。クエリポイントが境界ポイントによって形成される形状の内側にあるかどうかを判断するにはどうすればよいですか?そして、これのための素晴らしいブーストモジュールはありますか?調べてboost::geometry
いますが、まだ何も見つかりません。
あなたが探しているようですwithin
、違いますか?
http://www.boost.org/libs/geometry/doc/html/geometry/reference/algorithms/within/within_2.html
彼らがページに示している例は、実際には、ポイントインポリゴンです。
#include <iostream>
#include <list>
#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>
#include <boost/geometry/domains/gis/io/wkt/wkt.hpp>
int main()
{
typedef boost::geometry::model::d2::point_xy<double> point_type;
typedef boost::geometry::model::polygon<point_type> polygon_type;
polygon_type poly;
boost::geometry::read_wkt(
"POLYGON((2 1.3,2.4 1.7,2.8 1.8,3.4 1.2,3.7 1.6,3.4 2,4.1 3,5.3 2.6,5.4 1.2,4.9 0.8,2.9 0.7,2 1.3)"
"(4.0 2.0, 4.2 1.4, 4.8 1.9, 4.4 2.2, 4.0 2.0))", poly);
point_type p(4, 1);
std::cout << "within: " << (boost::geometry::within(p, poly) ? "yes" : "no") << std::endl;
return 0;
}
更新: @ildjarnが指摘しているようcovered_by
に、ポリゴンエッジ自体にあるポイントをカウントしたい場合は、次のように使用できます。
wrtエッジの動作はwithin
「依存」するため、ドキュメントでそのニュアンスに注意してください。