[postgis-devel] strange problem with operators

[Carl Anderson]

Working with postgresql-7.4.3 and  GEOS (cvs) and PostGIS (cvs)

I have found a condition where

select difference (geometry(shp1),geometry(shp2)) from zoning_test.bad;
NOTICE:  TopologyException: no outgoing dirEdge found 
(2.23561e+06,1.42152e+06)
ERROR:  GEOS difference() threw an error!

BUT

select difference (geometry(astext(shp1)),geometry(astext(shp2))) from 
zoning_test.bad;

succeeds

Any hints as how to extract the geometries for analysis



I got there by

shp1 = difference (shp1,shp0)
then
shp2 = difference( shp1,shp2)    -- fails                             




BTW
PostGIS (cvs) depends on Proj 4.4.8  (pj_get_release)
Can a note be added to the Makefile or CHANGES

C.


-- 
Carl Anderson
GIS Manager Fulton County, Georgia
carl.anderson at co.fulton.ga.us
404.730.8026