[Boost-bugs] [Boost C++ Libraries] #8366: "Overlay invalid input exception" (3)

Subject: [Boost-bugs] [Boost C++ Libraries] #8366: "Overlay invalid input exception" (3)
From: Boost C++ Libraries (noreply_at_[hidden])
Date: 2013-04-02 08:19:25

#8366: "Overlay invalid input exception" (3)
----------------------------------------------------------------+-----------
 Reporter:  Volker Schöch <vschoech@…>                          |       Owner:  barendgehrels
     Type:  Bugs                                                |      Status:  new          
Milestone:  To Be Determined                                    |   Component:  geometry     
  Version:  Boost 1.52.0                                        |    Severity:  Problem      
 Keywords:  polygon, multi polygon, difference, sym_difference  |  
----------------------------------------------------------------+-----------
 Please find below some code that triggers "Overlay invalid input
 exception" (in the last statement of the example). I have a couple of
 different reproductions, and since I'm not sure if they all share the same
 root cause, I'm filing them in separate tickets.

 As always, my polygon type is oriented '''counter-clockwise''' and '''not
 closed''', my point type is based on '''int'''.

 This is the data that triggers the exception:

 {{{
 _intPolygon polygon( "MULTIPOLYGON(((529 998,5337 998,5337 3475,529
 3475)))" );
 {
        _intPolygon const polygonB = _intPolygon("MULTIPOLYGON(((528
 3314,1734 2054,2934 1670,4140 1754,5340 2072,5340 32767,528 32767)))") ^
 _intPolygon("MULTIPOLYGON(((528 3218,1734 1784,2934 1400,4140 2582,5340
 1832,5340 32767,528 32767)))");
        polygon -= polygonB;
 }
 {
        _intPolygon const polygonB = _intPolygon("MULTIPOLYGON(((528
 3218,1734 1784,2934 1400,4140 2582,5340 1832,5340 32767,528 32767)))") ^
 _intPolygon("MULTIPOLYGON(((528 2498,1734 1406,2934 1574,4140 3002,5340
 1178,5340 32767,528 32767)))");
        polygon -= polygonB;
 }
 {
        _intPolygon const polygonB = _intPolygon("MULTIPOLYGON(((528
 2498,1734 1406,2934 1574,4140 3002,5340 1178,5340 32767,528 32767)))") ^
 _intPolygon("MULTIPOLYGON(((528 2420,1734 2186,2934 2378,4140 2750,5340
 1250,5340 32767,528 32767)))");
        polygon -= polygonB;
 }
 {
        _intPolygon const polygonB = _intPolygon("MULTIPOLYGON(((528
 2420,1734 2186,2934 2378,4140 2750,5340 1250,5340 32767,528 32767)))") ^
 _intPolygon("MULTIPOLYGON(((528 1724,1734 2552,2934 1640,4140 2396,5340
 1460,5340 32767,528 32767)))");
        polygon -= polygonB;
 }
 }}}

 This is my code that wraps boost::geometry to implement the operators used
 above:

 {{{
 template<typename T>
 template<typename Geometry>
 _TPolygon< T > _TPolygon< T >::operator-(Geometry const& geometry) const
 {
        // should not be necessary
        //if( boost::geometry::area(geometry)==0 ) return *this;

        _TPolygon< T > polygonOut;
        boost::geometry::difference(*this, geometry, polygonOut);

        // should not be necessary
        //boost::geometry::correct( polygonOut );

        return polygonOut;
 }

 template<typename T>
 template<typename Geometry>
 _TPolygon<T>& _TPolygon< T >::operator-=(Geometry const& geometry)
 {
        // boost::geometry::difference cannot operate in-place
        // http://lists.boost.org/geometry/2012/02/1796.php

        *this = *this - geometry;
        return *this;
 }

 template<typename T>
 template<typename Geometry>
 _TPolygon< T > _TPolygon< T >::operator^(Geometry const& geometry) const
 {
        _TPolygon< T > polygonOut;
        boost::geometry::sym_difference(*this, geometry, polygonOut);

        // should not be necessary
        //boost::geometry::correct( polygonOut );

        return polygonOut;
 }
 }}}

-- 
Ticket URL: <https://svn.boost.org/trac/boost/ticket/8366>
Boost C++ Libraries <http://www.boost.org/>
Boost provides free peer-reviewed portable C++ source libraries.