Question

В Boost есть отличная библиотека для геометрии. Это позволяет также рисовать изображения SVG. Я хочу использовать его в каком-то проекте, но он работает для меня очень странно (см. Изображение ниже).

Итак, у нас есть 3 точки пикселя, представленные в виде квадратных многоугольников в 2d пространстве

 1 1
 0 1

enter image description here pic 1

мы хотим получить от них объединение и упростить его, чтобы при масштабировании мы получили треугольник, подобный

1 1 1 1 1 1
1 1 1 1 1 1  
1 1 1 1 1 1
0 1 1 1 1 1 
0 0 1 1 1 1 
0 0 0 1 1 1

enter image description here pic 2

но мы получаем это:

enter image description here

, где желтая пунктирная линия - это объединение, а зеленая - упрощение.

Исходник:

#include <iostream>
#include <fstream>
#include <boost/assign.hpp>

#include <boost/algorithm/string.hpp>
#include <boost/geometry/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/multi/geometries/multi_polygon.hpp>
#include <boost/geometry/algorithms/envelope.hpp>

#include <boost/geometry/extensions/io/svg/svg_mapper.hpp>

template <typename Geometry1, typename Geometry2>
void create_svg(std::string const& filename, Geometry1 const& a, Geometry2 const& b)
{
    typedef typename boost::geometry::point_type<Geometry1>::type point_type;
    std::ofstream svg(filename.c_str());

    boost::geometry::svg_mapper<point_type> mapper(svg, 400, 400);
    mapper.add(a);
    mapper.add(b);

    mapper.map(a, "fill-opacity:0.5;fill:rgb(153,204,0);stroke:rgb(153,204,0);stroke-width:2");
    mapper.map(b, "opacity:0.8;fill:none;stroke:rgb(255,128,0);stroke-width:4;stroke-dasharray:1,7;stroke-linecap:round");
}

int main()
{

    // create points (each point == square poligon)
    boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > one, two, three;

    boost::geometry::read_wkt(
        "POLYGON((1 1, 1 0, 0 0, 0 1))", one);

    boost::geometry::read_wkt(
        "POLYGON((2 2, 2 1, 1 1, 1 2))", two);

    boost::geometry::read_wkt(
        "POLYGON((1 1, 1 2, 0 2, 0 1))", three);

    // create a container for joined points structure
    boost::geometry::model::multi_polygon< boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > > output, simpl;

    // join points one by one (because one day we would have many=))
    boost::geometry::union_(one, two, output);
    boost::geometry::union_( output , three, output);

    // simplify joined structure
    boost::geometry::simplify(output, simpl, 0.5);

    // create an svg image
    create_svg("make_envelope.svg", simpl, output );
}

требуется как минимум повышение 1.47.0 и 3 файла из boost / geometry / extensions / io / svg /

Итак, как сделать так, чтобы я упростил, как я хочу, чтобы смысл обрел форму, подобную pic 2?

Обновление

Создан новый код, работает корректно, достаточно проверено:

#include <iostream>
#include <fstream>
#include <boost/assign.hpp>

//Boost
#include <boost/algorithm/string.hpp>
#include <boost/geometry/geometry.hpp>
#include <boost/geometry/geometries/geometries.hpp>
#include <boost/geometry/multi/geometries/multi_polygon.hpp>
#include <boost/geometry/geometries/adapted/boost_tuple.hpp>

#include <boost/foreach.hpp>

//and this is why we use Boost Geometry from Boost trunk 
#include <boost/geometry/extensions/io/svg/svg_mapper.hpp>

BOOST_GEOMETRY_REGISTER_BOOST_TUPLE_CS(cs::cartesian)

template <typename Geometry1, typename Geometry2>
void create_svg(std::string const& filename, Geometry1 const& a, Geometry2 const& b)
{
    typedef typename boost::geometry::point_type<Geometry1>::type point_type;
    std::ofstream svg(filename.c_str());

    boost::geometry::svg_mapper<point_type> mapper(svg, 400, 400);
    mapper.add(a);
    mapper.add(b);

    mapper.map(a, "fill-rule:nonzero;fill-opacity:0.5;fill:rgb(153,204,0);stroke:rgb(153,204,0);stroke-width:2;");
    mapper.map(b, "opacity:0.8;fill:none;stroke:rgb(255,128,0);stroke-width:4;stroke-dasharray:1,7;stroke-linecap:round");
}


void make_point(int x, int y,  boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > & ring)
{
    using namespace boost::assign;

    boost::geometry::append(  ring,     boost::geometry::model::d2::point_xy<double>(x-1, y-1));
    boost::geometry::append(  ring,     boost::geometry::model::d2::point_xy<double>(x, y-1));
    boost::geometry::append(  ring,      boost::geometry::model::d2::point_xy<double>(x, y));
    boost::geometry::append(  ring,      boost::geometry::model::d2::point_xy<double>(x-1, y));
    boost::geometry::append(  ring,     boost::geometry::model::d2::point_xy<double>(x-1, y-1));
    boost::geometry::correct(ring);
}

void create_point(int x, int y, boost::geometry::model::multi_polygon< boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > > & mp)
{
    boost::geometry::model::multi_polygon< boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > > temp;
    boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > ring;
    make_point(x, y, ring);
    boost::geometry::union_(mp, ring, temp);
    boost::geometry::correct(temp);
    mp=temp;
}

int main()
{
    using namespace boost::assign;

    typedef boost::geometry::model::polygon
        <
        boost::geometry::model::d2::point_xy<double>
        > polygon;

    typedef boost::geometry::model::multi_polygon<polygon> mp;

    polygon ring;

    mp pol, simpl;
    polygon exring;

    create_point(1,1, pol);
    create_point(2, 1, pol);
    create_point(3, 1, pol);
    create_point(4,1, pol);
    create_point(5, 1, pol);

    create_point(1,2, pol);
    create_point(2, 2, pol);
    create_point(3, 2, pol);
    create_point(4,2, pol);
    create_point(5, 2, pol);

    create_point(2, 3, pol);
    create_point(3, 3, pol);
    create_point(5, 3, pol);

    create_point(3, 4, pol);

    create_point(5, 3, pol);

    create_point(5, 5, pol);

    //boost::geometry::dissolve(ring, pol); // Baad
    boost::geometry::simplify(pol, simpl, 0.5); // Good

    create_svg("make_envelope.svg",pol,  simpl );
}

И этот код создает такое изображение:

enter image description here

И за 3 балла он возвращает изображения одинаково @ J. Calleja ответ:

enter image description here

J. Calleja · Answer 1 · 09 ноября 2011

Я думаю, что есть несколько проблем с кодом:

Полигоны, которые вы определяете:

1 1
1 0

То есть:

three two
one    -

Таким образом, ожидаемый результат отличается от pic2.

Полигоны должны быть замкнуты и направленыпо часовой стрелке.

Отсутствует точка закрытия, а третий многоугольник не направлен по часовой стрелке.Взгляните на правильный метод.В этом примере вы должны вызывать его для каждого определенного полигона.

Нельзя использовать один и тот же аргумент для ввода и вывода при использовании _union .

Вы должны использовать временную переменную:

  boost::geometry::union_(one, two, outputTmp);    
  boost::geometry::union_( outputTmp, three, output);

Ваш ожидаемый результат может не быть результатом алгоритма.

После выполнения исправленного кода результат будет:

simplify result

Это может быть допустимым упрощением вашего многоугольника.См. Алгоритм Рамера – Дугласа – Пекера .

После выполнения этих изменений ниже приведен результат main ()

int main() 
{
  // create points (each point == square poligon)     
  boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > one, two, three;  
  boost::geometry::read_wkt(         "POLYGON((1 1, 1 0, 0 0, 0 1))", one);  
  boost::geometry::read_wkt(         "POLYGON((2 2, 2 1, 1 1, 1 2))", two); 
  boost::geometry::read_wkt(         "POLYGON((1 1, 1 2, 0 2, 0 1))", three);  
  boost::geometry::correct(one);
  boost::geometry::correct(two);
  boost::geometry::correct(three);

  // create a container for joined points structure  
  boost::geometry::model::multi_polygon< boost::geometry::model::polygon<boost::geometry::model::d2::point_xy<double> > > outputTmp, output, simpl;      
  // join points one by one (because one day we would have many=))    
  boost::geometry::union_(one, two, outputTmp);    
  boost::geometry::union_( outputTmp, three, output);    
  // simplify joined structure  
  boost::geometry::simplify(output, simpl, 0.5);   
  // create an svg image   
  create_svg("make_envelope.svg", simpl, output ); 
}

Boost :: Упрощение объединения геометрии - как это работает?

Обновление

Пожалуйста, войдите или зарегистрируйтесь чтобы ответить на этот вопрос.

1 Ответ

Пожалуйста, войдите или зарегистрируйтесь что бы добавить комментарий.

Boost :: Упрощение объединения геометрии - как это работает?

Обновление

Пожалуйста, войдите или зарегистрируйтесь чтобы ответить на этот вопрос.

1 Ответ

Пожалуйста, войдите или зарегистрируйтесь что бы добавить комментарий.

Похожие темы