Initialisation d'un polygone dans boost :: geometry

Question

Je suis nouveau dans la bibliothèque de géométrie générique qu'il est proposé d'inclure avec boost:

http://geometrylibrary.geodan.nl/

J'ai deux vecteurs vector<int> Xb, Yb à partir desquels je tente de créer un polygone. J'essaie d'obtenir quelque chose qui ressemble à l'extrait de code suivant:

 polygon_2d P;

 vector<double>::const_iterator xi;
 vector<double>::const_iterator yi;

 for (xi=Xb.begin(), yi=Yb.begin(); xi!=Xb.end(); ++xi, ++yi)
  P.push_back (make<point_2d>(*xi, *yi));

Le code ci-dessus ne fonctionne pas, en se plaignant que P n'a pas de fonction membre push_back. Comment initialiser le polygone à partir de points dont les coordonnées sont vector<int> Xb,vector<int> Yb?

Answer 1

append(P, make<point_2d>(*xi, *yi));

Answer 2

Voici un exemple d'extension de la question initiale que vous avez posée en commentaire sous la réponse de Kirill: Des intersections entre polygones sont-elles possibles?

Oui, les Boost.Geometry (aka GGL) prennent en charge les intersections polygone-polygone.

#include <iostream>
#include <vector>
#include <boost/geometry/geometry.hpp>
#include <boost/geometry/geometries/cartesian2d.hpp>
#include <boost/geometry/geometries/adapted/c_array_cartesian.hpp>

using namespace boost::geometry;

int main(void)
{
    // Define a polygons and fill the outer rings.
    polygon_2d a;
    {
        const double c[][2] = {
            {160, 330}, {60, 260}, {20, 150}, {60, 40}, {190, 20}, {270, 130}, {260, 250}, {160, 330}
        };
        assign(a, c);
    }
    correct(a);
    std::cout << "A: " << dsv(a) << std::endl;

    polygon_2d b;
    {
        const double c[][3] = {
            {300, 330}, {190, 270}, {150, 170}, {150, 110}, {250, 30}, {380, 50}, {380, 250}, {300, 330}
        };
        assign(b, c);
    }
    correct(b);
    std::cout << "B: " << dsv(b) << std::endl;

    // Calculate interesection
    typedef std::vector<polygon_2d > polygon_list;
    polygon_list v;

    intersection_inserter<polygon_2d>(a, b, std::back_inserter(v));
    std::cout << "Intersection of polygons A and B" << std::endl;
    for (polygon_list::const_iterator it = v.begin(); it != v.end(); ++it)
    {
        std::cout << dsv(*it) << std::endl;
    }

    return 0;
}

Ici correspond le résultat (l'intersection du polygone est déplacée vers le sud pendant meilleure visibilité):

J'espère que cela fonctionnera pour vous.

Answer 3

Vous pouvez également utiliser un tuple pour initialiser le polygone

#include <boost/geometry/geometries/adapted/boost_tuple.hpp>

et

boost::geometry::assign_points(
    polygon, boost::assign::tuple_list_of
        (300, 330) (190, 270) (150, 170) (150, 110) (250, 30) (380, 50)
        (380, 250) (300, 330)
);