在点和多边形上使用内部方法的麻烦

Boost::geometry Trouble in using within method on point and polygons?

本文关键字：内部 方法 麻烦 多边形      更新时间：2023-10-16

似乎我在正确使用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;
}

当您使用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;
}
相关文章：