pointinpoly.cxx

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>
#include <iostream>

using namespace boost::geometry;
using namespace std;

int main()
{
    typedef model::d2::point_xy<double> pointtype;
    pointtype point1(3, 3);
    pointtype point2(2, 12.1);

    model::polygon<pointtype> poly;
    pointtype points[] = {pointtype(0, 0), pointtype(0, 10), pointtype(5, 15),
                          pointtype(10, 10), pointtype(10, 0), pointtype(0, 0)};
    append(poly, points);

    cout << wkt(point1) << " is "
         << (within(point1, poly) ? "inside" : "outside of") << " the "
         << wkt(poly) << endl;
    cout << wkt(point2) << " is "
         << (within(point2, poly) ? "inside" : "outside of") << " the "
         << wkt(poly) << endl;
    return 0;
}