Boost :: Semplificazione dell'unione geometrica: come funziona?

Question

C'è una grande libreria per la geometria in Boost. Permette anche di disegnare immagini SVG. Voglio usarlo in qualche mio progetto ma funziona davvero strano per me (vedi immagine sotto).

Quindi abbiamo 3 punti pixel rappresentati come poligoni quadrati nello spazio 2d

1 1 
0 1 

pic 1

vogliamo ottenere da loro un sindacato e semplificare in modo che quando si scala essa si' d ottenere un triangolo come

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

ma otteniamo questo:

dove gialla linea tratteggiata è l'unione e il verde è la semplificazione.

Sourcecode:

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

richiede almeno spinta 1.47.0 e 3 file da boost/geometry/extensions/io/svg/

Così come farlo semplificare come voglio che significa ottenere forma come pic 2?

Aggiornamento

nuovo codice creato, funziona correttamente, abbastanza collaudato:

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

E questo codice crea tale immagine:

E per 3 punti che restituisce immagini simili @J. Calleja risposta:

Answer 1

Penso che ci siano diversi problemi con il codice:

• I poligoni si sta definendo sono:

Quello is:

three two 
one - 

Così il risultato atteso è diverso da pic2.

• I poligoni devono essere chiusi e orientati in senso orario.

Manca il punto di chiusura e il terzo poligono non è orientato in senso orario. Dai uno sguardo al metodo correct. In questo esempio, dovresti chiamarlo per ogni poligono che definisci.

• Non è possibile usare lo stesso argomento per l'input e l'output quando si utilizza _union.

si dovrebbe usare una variabile temporanea:

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

• Il risultato atteso non può essere il risultato algoritmo.

Dopo aver eseguito il codice corretto, il risultato è:

Questa può essere una valida simplifcation del poligono. Vedi lo Ramer–Douglas–Peucker algorithm.

Dopo l'esecuzione di tali modifiche, di seguito è il principale risultante()

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