[postgis-users] st_intersection error

[Martin Davis]

Same result in JTS.

As Nicklas points out, this is a problem of precision.  The software is 
working as designed, but it's not designed to use it to perform 
intersections with a tolerance.

Some more sophisticated approach is required, such as vertex snapping 
the lines to a given tolerance, or perhaps using a narrow buffer and 
some sort of segment-based containment tests.

strk wrote:
> On Tue, Sep 15, 2009 at 11:32:48AM +0200, Düster Horst wrote:
>   
>> strk
>>
>> I attached the geometries in my initial mail.
>>
>> Take a look at aline1.sql and aline2.sql
>>     
>
> I confirm a POINT and a LINESTRING in a COLLECTION
> is returned by ST_Intersection() with latest version:
>
> POSTGIS="1.5.0SVN" GEOS="3.2.0-CAPI-1.6.0" PROJ="Rel. 4.6.0, 21 Dec 2007" USE_STATS
>
> Martin, how about JTS ?
>
> --strk; 
>
>  Free GIS & Flash consultant/developer      ()  ASCII Ribbon Campaign
>  http://foo.keybit.net/~strk/services.html  /\  Keep it simple! 
>
>   

-- 
Martin Davis
Senior Technical Architect
Refractions Research, Inc.
(250) 383-3022