[geos-devel] [GEOS] #230: polygon intersection exception.
Martin Davis
mbdavis at refractions.net
Wed Feb 25 11:47:01 EST 2009
The second polygon is invalid - this is expected behaviour.
GEOS wrote:
> #230: polygon intersection exception.
> ------------------------+---------------------------------------------------
> Reporter: anigold | Owner: geos-devel at lists.osgeo.org
> Type: defect | Status: new
> Priority: major | Milestone: 3.1.0
> Component: Default | Version: 3.0.0
> Severity: Unassigned | Keywords:
> ------------------------+---------------------------------------------------
> {{{
> std::string str = "POLYGON ((112.3903710739794900
> 22.3851129617344210, 114.1261324827294900 22.3851129617344210,
> 114.1261324827294900 23.5479374992368590, 112.3903710739794900
> 23.5479374992368590, 112.3903710739794900 22.3851129617344210)) ";
> geos::geom::Geometry* pgeom0 = reader.read( str );
>
> str = "POLYGON ((114.0876590000000000 22.5249449999999990,
> 114.1853440000000000 22.5531790000000000, 114.1704040000000000
> 22.5438380000000010, 114.1423630000000000 22.5351350000000000,
> 114.0874290000000000 22.5245199999999990, 114.0892680000000000
> 22.5255820000000000, 114.0936350000000100 22.5266430000000000,
> 114.0876590000000000 22.5249449999999990))";
> geos::geom::Geometry* pgeom1 = reader.read( str );
>
> geos::geom::Geometry* pint = pgeom0->intersection( pgeom1 );
>
> }}}
>
> when execute the last line code, the util::TopologyException will be
> throwed.
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> geos-devel mailing list
> geos-devel at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/geos-devel
--
Martin Davis
Senior Technical Architect
Refractions Research, Inc.
(250) 383-3022
More information about the geos-devel
mailing list