C++ GEOS库判断点与多边形的关系代码示例（例子）

C++中的GEOS库（Geometry Engine Open Source）是一个强大的地理空间操作库，它提供了一系列用于处理几何对象的功能。如果你想在C++中判断一个点与多边形的关系，可以使用`geos::geom::Point`和`geos::geom::Polygon`这两个类，以及它们的相关函数。

这里有一个简单的例子，假设你已经有了一个`geos::geom::Point`类型的点和一个`geos::geom::Polygon`类型的多边形：

#include <geos/geometry.hpp>
#include <geos/operation/overlay/InteriorityOperator.hpp>

// 假设已经创建了Point p和Polygon poly
std::unique_ptr<geos::geom::Geometry> point = ...; // 创建点对象
std::unique_ptr<geos::geom::Polygon> polygon = ...; // 创建多边形对象

// 使用InteriorityOperator判断点和多边形的关系
geos::operation::overlay::InteriorityOperator overlayOp;
geos::geom::IntersectionMatrix matrix = overlayOp.apply(*point, *polygon);

// 点在多边形内则matrix结果为INSIDE, 在边界上为BOUNDARY，外部为EXTerior
if (matrix.getOnSegment()) {
    std::cout << "点在多边形边上" << std::endl;
} else if (matrix.isInterior()) {
    std::cout << "点在多边形内部" << std::endl;
} else {
    std::cout << "点在多边形外部" << std::endl;
}