Opened 14 years ago

Closed 14 years ago

#386 closed defect (invalid)

Within Function gives wrong result, Intersection threw an error.

Reported by: svenanders Owned by: pramsey
Priority: medium Milestone: PostGIS 1.4.2
Component: postgis Version:
Keywords: Cc:

Description

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.

Attachments (4)

schema.sql (307 bytes ) - added by svenanders 14 years ago.
Schema File for example Table
62771 (526.9 KB ) - added by svenanders 14 years ago.
OSM Relation 62771 Lower Saxony in SQL
62676 (63.1 KB ) - added by svenanders 14 years ago.
OSM Relation 62676 Cuxhaven in SQL Syntax
bug.sql (200 bytes ) - added by svenanders 14 years ago.
SQL Commands that fails

Download all attachments as: .zip

Change History (8)

by svenanders, 14 years ago

Attachment: schema.sql added

Schema File for example Table

by svenanders, 14 years ago

Attachment: 62771 added

OSM Relation 62771 Lower Saxony in SQL

by svenanders, 14 years ago

Attachment: 62676 added

OSM Relation 62676 Cuxhaven in SQL Syntax

by svenanders, 14 years ago

Attachment: bug.sql added

SQL Commands that fails

comment:1 by pramsey, 14 years ago

Milestone: PostGIS 1.4.1PostGIS 1.4.2

comment:2 by pramsey, 14 years ago

Resolution: invalid
Status: newclosed

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)

comment:3 by svenanders, 14 years ago

Resolution: invalid
Status: closedreopened

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?

comment:4 by svenanders, 14 years ago

Resolution: invalid
Status: reopenedclosed

Okay i found the error its my fault….

Note: See TracTickets for help on using tickets.