Ticket #386 (closed defect: invalid)

Opened 3 years ago

Last modified 3 years ago

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

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

Change History

Changed 3 years ago by svenanders

Schema File for example Table

Changed 3 years ago by svenanders

OSM Relation 62771 Lower Saxony in SQL

Changed 3 years ago by svenanders

OSM Relation 62676 Cuxhaven in SQL Syntax

Changed 3 years ago by svenanders

SQL Commands that fails

Changed 3 years ago by pramsey

  • milestone changed from PostGIS 1.4.1 to PostGIS 1.4.2

Changed 3 years ago by pramsey

  • status changed from new to closed
  • resolution set to invalid

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)

Changed 3 years ago by svenanders

  • status changed from closed to reopened
  • resolution invalid deleted

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?

Changed 3 years ago by svenanders

  • status changed from reopened to closed
  • resolution set to invalid

Okay i found the error its my fault....

Note: See TracTickets for help on using tickets.