[postgis-users] Intersection doesn't work on some data

[David Mitchell]

I have got two MULTIPOLYGON features, each containing a single POLYGON. 
(attached). When I ask postgis to find the intersection of these, it 
gives the wrong result. It should give:

MULTILINESTRING ((-121.331929 41.18388, -121.34693 41.18388), 
(-121.34693 41.18388, -121.36293 41.18388), (-121.36293 41.18388, 
-121.36544 41.18388), (-121.36544 41.18388, -121.37179 41.18388), 
(-121.37179 41.18388, -121.37608 41.18401), (-121.37608
41.184019, -121.44593 41.18384))

But in fact gives:

MULTILINESTRING((-121.33193 41.18388,-121.34693 41.18388),(-121.34693 
41.18388,-121.36293 41.18388),(-121.36293 41.18388,-121.36545 
41.18388),(-121.36545 41.18388,-121.37179 41.18388),(-121.37179 
41.18388,-121.37608 41.18402),(-121.37608 41.18402,1.7e-308 
-1.00226318765399e-13))

Note the final line string has very tiny numbers for its second 
coordinate. I have written a program to run these geometries directly 
through geos and the output is correct, so it must be a postgis problem.

Any help would be appreciated.

-- 
David Mitchell
Software Engineer
Telogis

NOTICE:
This message (including any attachments) contains CONFIDENTIAL
INFORMATION intended for a specific individual and purpose, and
is protected by law.  If you are not the intended recipient,
you should delete this message and are hereby notified that any
disclosure, copying, or distribution of this message, or the
taking of any action based on it, is strictly prohibited.
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: counties.wkt
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20040324/e4b8397b/attachment.ksh>