Ticket #8310: geom.cpp

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

using namespace boost::geometry;

typedef model::d2::point_xy<float> point;
typedef model::polygon<point, true, false> polygon;


void print( const polygon &poly )
{
        bool first = true;
        std::cout << "POLYGON(( ";
        for( const point &p : poly.outer() ) {
                if( first ) {
                        first = false;
                }
                else {
                        std::cout << ", ";
                }
                std::cout << p.x() << " " << p.y();
        }
        std::cout << " ))" << std::endl;
}
void print( const std::vector< polygon > &out )
{
        for( const polygon &poly : out ) {
                std::cout << "\t";
                print( poly );
        }
        std::cout << std::endl;
}

void testDiff( const polygon &a, const polygon &b )
{
        std::vector< polygon > out;
        boost::geometry::difference( a, b, out );

        print( out );
}
void testSymDiff( const polygon &a, const polygon &b )
{
        std::vector< polygon > out;
        boost::geometry::sym_difference( a, b, out );
        print( out );
}
void testIntersection( const polygon &a, const polygon &b )
{
        std::vector< polygon > out;
        boost::geometry::intersection( a, b, out );
        print( out );
}
void testUnion( const polygon &a, const polygon &b )
{
        std::vector< polygon > out;
        boost::geometry::union_( a, b, out );
        print( out );
}
void test( const polygon &a, const polygon &b )
{
        std::cout << "a: ";
        print( a );
        std::cout << "b: ";
        print( b );
        std::cout << std::endl;

        std::cout << "a - b:" << std::endl;
        testDiff( a, b );
        std::cout << "b - a:" << std::endl;
        testDiff( b, a );
        std::cout << "b ^ a:" << std::endl;
        testSymDiff( a, b );
        std::cout << "b & a:" << std::endl;
        testIntersection( b, a );
        std::cout << "b | a:" << std::endl;
        testUnion( b, a );
        std::cout << "disjoint: " << boost::geometry::disjoint( a, b ) << std::endl;
        std::cout << "intersects: " << boost::geometry::intersects( a, b ) << std::endl;
}
void test( const std::string &aStr, const std::string &bStr )
{
        polygon a, b;

        boost::geometry::read_wkt( aStr, a);
        boost::geometry::correct( a );
        boost::geometry::read_wkt( bStr, b );
        boost::geometry::correct( b );

        test( a, b );
}

int main( )
{
        std::cout.precision( 15 );
        std::cout << std::fixed;

        test( "POLYGON(( -2.559375047683716 -0.615625500679016, -2.559375047683716 0.384374797344208, 7.940625190734863 0.384374588727951, 7.940625190734863 -0.615625441074371 ))",
                  "POLYGON(( 1.000000000000000 0.384374707937241, 1.000000000000000 0.000000000000000, 0.000000000000000 0.000000000000000, 0.000000000000000 0.384374737739563 ))" );
        std::cout << " = = = = = " << std::endl << std::endl;

        test( "POLYGON(( 0.000000000000000 0.373437941074371, 1.000000000000000 0.373437792062759, 1.000000000000000 0.000000000000000, 0.000000000000000 0.000000000000000 ))",
                  "POLYGON(( -2.592187881469727 -0.626561701297760, -2.592187643051147 0.373438000679016, 7.907812595367432 0.373437851667404, 7.907812595367432 -0.626561224460602 ))" );
        std::cout << " = = = = = " << std::endl << std::endl;

        test( "POLYGON(( 1.000000000000000 2.531045913696289, 1.000000000000000 3.000000000000000, 2.000000000000000 3.000000000000000, 2.000000000000000 2.531045913696289 ))",
                  "POLYGON(( 5.204249382019043 3.531043529510498, 5.204247951507568 2.531045675277710, -5.295750617980957 2.531046152114868, -5.295751094818115 3.531045913696289 ))" );
        std::cout << " = = = = = " << std::endl << std::endl;

        std::cout << "This one is working as expected...." << std::endl;
        test( "POLYGON(( 1. 2., 1. 3., 2. 3., 2. 2. ))",
                  "POLYGON(( 5. 3., 5. 2., -5. 2., -5. 3. ))" );

        /*

        test( "POLYGON((  ))",
                  "POLYGON((  ))" );
        std::cout << " = = = = = " << std::endl << std::endl;
        */
        return 0;
}