Within Function gives wrong result, Intersection threw an error.

[svenanders]

I have two boundarys , first one is Lower Saxony

the second one is Cuxhaven.

There are both from OpenStreetMap.

Cuxhaven is in Lower Saxony but the Within function tells is always false.

When I try to get the intersection i got the following error.

TopologyException: no outgoing dirEdge found at 8.49226 53.4721 GEOS Intersection() threw an error!

POSTGIS="1.4.1" GEOS="3.2.0-CAPI-1.6.0" PROJ="Rel. 4.6.0, 21 Dec 2007" USE_STATS

Postgresql 8.3

Debian lenny.

[svenanders]

Schema File for example Table

[svenanders]

OSM Relation 62771 Lower Saxony in SQL

[svenanders]

OSM Relation 62676 Cuxhaven in SQL Syntax

[svenanders]

SQL Commands that fails

[pramsey]

These polygons are invalid

geog=# select isvalid(g) from boundary_test;
NOTICE:  Holes are nested at or near point 8.55575 53.5249
NOTICE:  Holes are nested at or near point 8.65388 53.6025
 isvalid 
---------
 f
 f
(2 rows)

[svenanders]

Sorry, i beleve that at least the Polygon for cuxhaven is valid. I have looked at our map at:

​http://www.openstreetmap.org/?mlon=8.65388&mlat=53.6025

Also in the source polygon there is no point: 8.65388 53.6025

Is postigs rounding coords?

[svenanders]

Okay i found the error its my fault….