Example: math_polygon_intersection

Example for the C++ function mrpt::math::intersect(const TPolygon2D& subject, const TPolygon2D& clipping) for polygon clipping or intersection using the Sutherland-Hodgman algorithm.

C++ example source code:

/* +------------------------------------------------------------------------+
   |                     Mobile Robot Programming Toolkit (MRPT)            |
   |                          https://www.mrpt.org/                         |
   |                                                                        |
   | Copyright (c) 2005-2024, Individual contributors, see AUTHORS file     |
   | See: https://www.mrpt.org/Authors - All rights reserved.               |
   | Released under BSD License. See: https://www.mrpt.org/License          |
   +------------------------------------------------------------------------+ */

#include <mrpt/gui/CDisplayWindowPlots.h>
#include <mrpt/math/TObject2D.h>
#include <mrpt/math/TPolygon2D.h>
#include <mrpt/math/geometry.h>

#include <iostream>

static void polygonClippingExample()
{
    // Define the polygons:
    const mrpt::math::TPolygon2D subject = {{
        {0.0, 0.0},
        {5.0, 0.0},
        {7.0, 3.0},
        {3.0, 6.0},
        {-4.0, 4.0},
        {-1.0, -1.0}
        //
    }};

    const mrpt::math::TPolygon2D clipping = {{
        {-6.0, 0.5}, {8.0, 2.0}, {10.0, 4.0}, {-7.0, 3.0}
        //
    }};

    // Compute intersection:
    mrpt::math::TObject2D clippedObj;
    bool doIntersect = mrpt::math::intersect(subject, clipping, clippedObj);
    ASSERT_(doIntersect);

    mrpt::math::TPolygon2D clippedPoly;
    bool isPoly = clippedObj.getPolygon(clippedPoly);
    ASSERT_(isPoly);

    // GUI:
    mrpt::gui::CDisplayWindowPlots winIn(
        "Inputs (red:subject, blue:clipping)", 500, 500);
    winIn.setPos(20, 50);

    const auto [x1, y1] = subject.getPlotData();
    winIn.plot(x1, y1, "-r3");
    winIn.hold_on();

    const auto [x2, y2] = clipping.getPlotData();
    winIn.plot(x2, y2, ":b2");

    mrpt::gui::CDisplayWindowPlots winOut("Result", 500, 500);
    winOut.setPos(600, 50);

    winOut.plot(x2, y2, ":b2");
    winOut.hold_on();

    const auto [xp, yp] = clippedPoly.getPlotData();
    winOut.plot(xp, yp, "-r3");

    winIn.axis_equal();
    winOut.axis_equal();

    std::cout << "Press any key on input window or close it to end program.\n";
    winIn.waitForKey();
}

int main(int argc, char** argv)
{
    try
    {
        polygonClippingExample();
        return 0;
    }
    catch (const std::exception& e)
    {
        std::cout << e.what() << std::endl;
        return 1;
    }
}