[geos-devel] troubles with gPolygonize

[Masucci, A]

Hi,
gPolygonize doesn't work with some simple geometries.
Here is an example:

> LS = list(
+ readWKT("LINESTRING (425963 576719, 425980 576703)"),
+ readWKT("LINESTRING (425963 576719, 425882 577073)"),
+ readWKT("LINESTRING (425980 576703, 426082 577072)"),
+ readWKT("LINESTRING (425882 577073, 426082 577072)"),
+ readWKT("LINESTRING (426138 577068, 426082 577072)"),
+ readWKT("LINESTRING (426138 577068, 426420 577039)"),
+ readWKT("LINESTRING (426420 577039, 426554 576990)"),
+ readWKT("LINESTRING (426751 576924, 426776 576823)"),
+ readWKT("LINESTRING (426751 576924, 426783 576919)"),
+ readWKT("LINESTRING (426751 576924, 426714 576953)"),
+ readWKT("LINESTRING (426776 576823, 426783 576919)"),
+ readWKT("LINESTRING (426658 576966, 426554 576990)"),
+ readWKT("LINESTRING (426658 576966, 426667 577031)"),
+ readWKT("LINESTRING (426658 576966, 426714 576953)"),
+ readWKT("LINESTRING (426667 577031, 426714 576953)")
+ )
> gPolygonize(LS)

Traceback:
 1: .Call("rgeos_polygonize", .RGEOS_HANDLE, splist, id, p4s, getCutEdges,     PACKAGE = "rgeos")
 2: gPolygonize(LS)
Error in gPolygonize(LS) : caught access violation - continue with care

I hope you can help me
Kindest regards

Paolo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/geos-devel/attachments/20130205/b666fe52/attachment.html>