The second polygon is invalid - this is expected behaviour.
GEOS wrote:
#230: polygon intersection exception.
------------------------+---------------------------------------------------
Reporter: anigold | Owner: geos-devel@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@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/geos-devel
--
Martin Davis
Senior Technical Architect
Refractions Research, Inc.
(250) 383-3022
_______________________________________________
geos-devel mailing list
geos-devel@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/geos-devel