альтернатива boost::geometry::read_wkt?

Я пытаюсь проверить, находится ли точка внутри многоугольника или нет. Для этого я хочу использовать библиотеку boost. Мой вопрос заключается в том, как я могу изменить пример из boost, чтобы вместо чтения точек с помощью read_wkt читать их из стандартного вектора.

вот пример кода из boost:

#include <iostream>
#include <list>

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


int main()
{
    typedef boost::geometry::model::d2::point_xy<double> point_type;
    typedef boost::geometry::model::polygon<point_type> polygon_type;

    polygon_type poly;
    boost::geometry::read_wkt(
        "POLYGON((2 1.3,2.4 1.7,2.8 1.8,3.4 1.2,3.7 1.6,3.4 2,4.1 3,5.3 2.6,5.4 1.2,4.9 0.8,2.9 0.7,2 1.3)"
            "(4.0 2.0, 4.2 1.4, 4.8 1.9, 4.4 2.2, 4.0 2.0))", poly);

    point_type p(4, 1);

    std::cout << "within: " << (boost::geometry::within(p, poly) ? "yes" : "no") << std::endl;

    return 0;
}

и вот мои vecPoints:

struct PNTS{
int x;
int y;
}
std::vector<PNTS> vecPoints;

ОБРАТИТЕ ВНИМАНИЕ: мой вопрос НЕ связан с алгоритмом точки в многоугольнике.


person H'H    schedule 07.11.2014    source источник
comment
Почему бы не перебрать ваш вектор и не построить строку WKT? Будет ли это работать?   -  person Anas Alkhatib    schedule 07.11.2014
comment
Не работает, наверное...   -  person H'H    schedule 07.11.2014
comment
@AnasAlkhatib добавил этот подход в качестве моего второго ответа   -  person sehe    schedule 08.11.2014


Ответы (3)


если вы хотите добавить свой векторный элемент внутри многоугольника в boost:

using boost::geometry::append;
using boost::geometry::make;
using boost::geometry::correct;
std::vector<PNTS>::iterator iter;
for(iter = vecPoints.begin(); iter != vecPoints.end(); iter++)
{
  PNTS ver = *iter;
  append( poly, make<boost2dPoint>(ver.x, ver.y) );
}
// you have to close polygon by inserting first element as the last again
PNTS last = vecPoints[vecPoints.size()-1];
append( poly, make<boost2dPoint>(last.x, last.y) );
// you can also correct the polygon orientation
correct(poly);
person user2090491    schedule 09.11.2014
comment
как это может работать нормально, он дважды добавляет последнюю точку в конце.. - person Pa_; 27.02.2018

Я не очень разбираюсь в концепциях геометрии. Тем не менее, похоже, что это должно быть близко или, по крайней мере, дать вам некоторое вдохновение о том, как что-то делать.

Обновление добавило сериализацию в WKT, чтобы мы могли подтвердить, что полион идентичен.

Жить на Coliru

#include <iostream>
#include <list>

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

int main()
{
    typedef boost::geometry::model::d2::point_xy<double> point_type;
    typedef boost::geometry::model::polygon<point_type> polygon_type;

    polygon_type poly;
    poly.outer().assign({
        point_type {     2,   1.3   },
        point_type {   2.4,   1.7   },
        point_type {   2.8,   1.8   },
        point_type {   3.4,   1.2   },
        point_type {   3.7,   1.6   },
        point_type {   3.4,     2   },
        point_type {   4.1,     3   },
        point_type {   5.3,   2.6   },
        point_type {   5.4,   1.2   },
        point_type {   4.9,   0.8   },
        point_type {   2.9,   0.7   },
        point_type {     2,   1.3   },
    });

    poly.inners().emplace_back();
    poly.inners().back().assign({
        {   4.0,   2.0   },
        {   4.2,   1.4   },
        {   4.8,   1.9   },
        {   4.4,   2.2   },
        {   4.0,   2.0   },
    });

    point_type p(4, 1);

    std::cout << "within: " << (boost::geometry::within(p, poly) ? "yes" : "no") << std::endl;

    std::cout << boost::geometry::wkt(poly) << "\n";
}

Отпечатки

within: yes
POLYGON((2 1.3,2.4 1.7,2.8 1.8,3.4 1.2,3.7 1.6,3.4 2,4.1 3,5.3 2.6,5.4 1.2,4.9 0.8,2.9 0.7,2 1.3),(4 2,4.2 1.4,4.8 1.9,4.4 2.2,4 2))
person sehe    schedule 07.11.2014
comment
спасибо, есть ошибка, с которой я сталкиваюсь при использовании boost 1.56: ошибка: нет соответствующей функции для вызова 'boost::geometry::model::ring‹boost::geometry::model::d2::point_xy‹double›, true, true, std::vector, std::allocator›::assign(‹список инициализаторов, заключенный в фигурные скобки›)' - person H'H; 09.11.2014
comment
Ну что ж. Очевидно, что ваш компилятор немного озадачен С++ 11. Для удобства я бы инициализировал временный вектор по диапазону итератора: paste.ubuntu.com/8901103 (Coliru на данный момент не работает) - person sehe; 09.11.2014

В случае, если, как сказал комментатор, вы считаете, что проще создать WKT из векторов вручную, вот один из способов сделать это в ~ 2 строках кода, не используя ничего из Boost Geometry:

using namespace boost::spirit::karma;
std::cout << format_delimited("POLYGON(" << *('(' << auto_%',' << ')') << ")\n", ' ', rings);

Полная демонстрация Live On Coliru

#include <iostream>
#include <boost/fusion/adapted/boost_tuple.hpp>
#include <boost/spirit/include/karma.hpp>

int main()
{
    using Ring = std::vector<boost::tuple<double,double>>;

    std::vector<Ring> rings = {
        {   {   2, 1.3 }, { 2.4, 1.7 }, { 2.8, 1.8 }, { 3.4, 1.2 },
            { 3.7, 1.6 }, { 3.4,   2 }, { 4.1,   3 }, { 5.3, 2.6 },
            { 5.4, 1.2 }, { 4.9, 0.8 }, { 2.9, 0.7 }, {   2, 1.3 },
        },
        {   { 4.0, 2.0 }, { 4.2, 1.4 }, { 4.8, 1.9 }, { 4.4, 2.2 }, { 4.0, 2.0   },
        }
    };

    using namespace boost::spirit::karma;
    std::cout << format_delimited("POLYGON(" << *('(' << auto_%',' << ')') << ")\n", ' ', rings);
}

Отпечатки:

POLYGON( ( 2.0 1.3 , 2.4 1.7 , 2.8 1.8 , 3.4 1.2 , 3.7 1.6 , 3.4 2.0 , 4.1 3.0 , 5.3 2.6 , 5.4 1.2 , 4.9 0.8 , 2.9 0.7 , 2.0 1.3 ) ( 4.0 2.0 , 4.2 1.4 , 4.8 1.9 , 4.4 2.2 , 4.0 2.0 ) )

person sehe    schedule 07.11.2014