2

看来我在正确使用 boost 几何库时遇到了麻烦（Boost::geometry ：计算多边形的质心）。我非常感谢对我上一个问题的帮助，并想问一些关于 boost::geometry::within 方法的问题，回答一个几何图形是否包含在另一个几何图形中。

我在我的代码中使用它来检查一个点是否包含在多边形中，并且我遇到了奇怪的结果，其中一个点绝对不应该在很远的多边形内，但该方法True在调用时仍然返回。

我正在考虑我在声明我的多边形时缺少的一个微妙之处，我真的很想解决这个问题。虽然来回查看我的代码，但我在调试时缺乏想法，感觉就像一个隧道视野。这就是为什么我想对这个特定的例子有所提示：

我的观点有坐标：221.703 , 256

多边形是坐标，逐点：

266.158 256 266.447 256.5 267.024 256.5 267.313 257 267.024 257.5 267.313 258

这显然不应该包含上面给出的观点。

我很抱歉问这个迂腐的东西，但我真的很感谢任何愿意把鼻子放在这个问题上的人

我的代码：

#include <iostream>
#include <boost/geometry.hpp>

using namespace std;
namespace bg = boost::geometry;
typedef bg::model::point<float, 2, bg::cs::cartesian> point;
typedef bg::model::box<point> box;
typedef bg::model::polygon<point, false, true> polygon;


int main(int argc, char * argv[]){
polygon pol;
pol.outer().push_back(point(266.158,256));
pol.outer().push_back(point(266.447,256.5));
pol.outer().push_back(point(267.024,256.5));
pol.outer().push_back(point(267.313,257));
pol.outer().push_back(point(267.024,257.5));
pol.outer().push_back(point(267.313,258));

double x = atof(argv[1]);
double y = atof(argv[2]);
cout << "Is inside: " << ((bg::within(point(x,y),pol)) ? "yes" : "no") << endl;
return 0;
}
4

1 回答 1

3

当您使用时，bg::model::polygon<point, false, true>您正在定义一个多边形，该多边形point用作其点类型，其点按逆时针顺序排列并且是闭合的（意味着它的最后一个点等于它的第一个点）。如果您“关闭”多边形或使用开放多边形，则 的行为bg::within似乎是您所期望的：

在 Coliru 上运行

#include <iostream>
#include <boost/geometry.hpp>

using std::cout;
using std::endl;
namespace bg = boost::geometry;
typedef bg::model::point<float, 2, bg::cs::cartesian> point;
typedef bg::model::box<point> box;
typedef bg::model::polygon<point, false, true> closed_polygon;
typedef bg::model::polygon<point, false, false> open_polygon;



int main(int argc, char * argv[])
{
    {
        closed_polygon pol;
        pol.outer().push_back(point(266.158,256));
        pol.outer().push_back(point(266.447,256.5));
        pol.outer().push_back(point(267.024,256.5));
        pol.outer().push_back(point(267.313,257));
        pol.outer().push_back(point(267.024,257.5));
        pol.outer().push_back(point(267.313,258));
        pol.outer().push_back(point(266.158,256));//you need to close the polygon

        double x = 222;
        double y = 257;
        cout << "Is " << bg::wkt<point>(point(x,y)) << " inside: " << ((bg::within(point(x,y),pol)) ? "yes" : "no") << endl;

        x = 267;
        y = 257;
        cout << "Is " << bg::wkt<point>(point(x,y)) << " inside: " << ((bg::within(point(x,y),pol)) ? "yes" : "no") << endl;
    }

    {
        open_polygon pol;
        pol.outer().push_back(point(266.158,256));
        pol.outer().push_back(point(266.447,256.5));
        pol.outer().push_back(point(267.024,256.5));
        pol.outer().push_back(point(267.313,257));
        pol.outer().push_back(point(267.024,257.5));
        pol.outer().push_back(point(267.313,258));

        double x = 222;
        double y = 257;
        cout << "Is " << bg::wkt<point>(point(x,y)) << " inside: " << ((bg::within(point(x,y),pol)) ? "yes" : "no") << endl;

        x = 267;
        y = 257;
        cout << "Is " << bg::wkt<point>(point(x,y)) << " inside: " << ((bg::within(point(x,y),pol)) ? "yes" : "no") << endl;
    }

    return 0;
}
于 2014-05-07T14:00:33.110 回答