boost::geometry::difference

#include <iostream>

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

#include <boost/foreach.hpp>

int main()
{
    typedef boost::geometry::model::d2::point_xy<double> POINT;
    typedef boost::geometry::model::linestring<POINT> LINE;
    typedef boost::geometry::model::polygon<POINT> POLYGON;

    LINE line;
    POLYGON polygon;

    //                     Y
    //                     |
    //               -------------
    //               |     |     |
    //               |     |     |
    //   ------A>>>>>C>>>>>>>>>>>D>>>>>B-----> X
    //       start   |     |     |    end
    //       point   |     |     |   point
    //               ------|------
    //                     |
    //

    // Create line from A to B.
    boost::geometry::read_wkt("linestring(-2 0, 2 0)", line);
    boost::geometry::correct(line);

    boost::geometry::read_wkt("POLYGON((1 1, -1 1, -1 -1, 1 -1))", polygon);
    boost::geometry::correct(polygon);

    std::vector<LINE> output;
    boost::geometry::difference(line, polygon, output);

    // Result will be 2 lines.
    // A to C : LINESTRING(-2 0,-1 0)
    // D to B : LINESTRING(1 0,2 0)
    BOOST_FOREACH(LINE const& l, output)
    {
        std::cout << boost::geometry::wkt(l) << std::endl;
    }

    return 0;
}