Zwiększenie :: Uproszczenie geometrii - jak to działa?

Question

W bibliotece Boost znajduje się świetna biblioteka do geometrii. Pozwala także rysować obrazy SVG. Chcę go użyć w moim projekcie, ale działa to dla mnie naprawdę dziwnie (patrz zdjęcie poniżej).

Więc mamy 3 punkty pixel reprezentowane kwadratowych poligons w 2D przestrzeni

1 1 
0 1 

pic 1

chcemy uzyskać od nich unii i uprościć je tak, że kiedy skala to my” d dostać trójkąt jak

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 

pic 2

ale mamy to:

gdzie żółty Doted linia jest unia i zielony to uproszczenie.

źródłowego:

#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); 
} 

wymaga co najmniej impuls 1.47.0 i 3 pliki z boost/geometry/extensions/io/svg/

Więc jak zrobić to jak ja chcą uprościć czyli aby uzyskać kształt jak pic 2?

Aktualizacja

Utworzono nowy kod, działa poprawnie, całkiem testowane:

#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); 
} 

I ten kod tworzy taki obraz:

I za 3 punkty zwraca obrazów alike @J. Calleja odpowiedź:

Answer 1

Myślę, że istnieje kilka problemów z kodem:

• Wielokąty są definiowane są:

Że jest:

three two 
one - 

Więc oczekiwany wynik jest różny od pic2.

• Wielokąty powinny być zamknięte i skierowane zgodnie z ruchem wskazówek zegara.

Brakuje punktu zamknięcia, a trzeci wielokąt nie jest skierowany zgodnie z ruchem wskazówek zegara. Spójrz na metodę correct. W tym przykładzie powinieneś go nazwać dla każdego zdefiniowanego wielokąta.

• Nie można przy użyciu _union użyć tego samego argumentu dla wejścia i wyjścia.

Należy użyć zmiennej tymczasowej:

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

• Twój oczekiwany wynik nie może być wynikiem algorytmu.

Po wykonaniu poprawiony kod, wynik jest:

To może być ważny simplifcation swojego wielokąta. Zobacz Ramer–Douglas–Peucker algorithm.

Po dokonaniu tych zmian, poniżej jest otrzymany główny()

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); 
}