[postgis-users] Intersects() and large datasets: error?

[strk at refractions.net]

Hubert,
GEOS 1.0.0 was not robust and that brought to those kind of errors.
Try GEOS 2.0.0.
--strk;

On Wed, Sep 22, 2004 at 04:32:29PM +0200, Hubert Fröhlich wrote:
> Hi list,
> 
> I am working with PostgreSQL 7.4.1, Postgis 0.8.2 and GEOS 1.0.0 .
> 
> I have a table with some columns and a geometry column for MULTIPOLYGONs
> 
> I try to test the scenario "give me the geometry containing a point with 
> given coordinates (borders included)" with some SQL statement like
> 
> select col1, col2  from table1 where pgpoly && 
> geometryfromtext('POINT(4427268.73 5477865.96))', 31494) and 
> intersects(geometryfromtext('POINT(4427268.73 5477865.96))', 31494), 
> pgpoly);
> 
> When the table is filled with a sub-dataset (ca. 80 000 multipolygons) 
> of my final data, everything works fine.
> 
> When I use my dataset (ca. 10 mill. multipolygons) , I get errors like
> 
> NOTICE:  TopologyException: side location conflict (4.42661e+06,5.47702e+06)
> ERROR:  GEOS intersects() threw an error!
> 
> The box overlap alone
> 
> select col1, col2  from table1 where pgpoly && 
> geometryfromtext('POINT(4427268.73 5477865.96))', 31494)
> 
> works properly.
> What has gone wrong?
> 
> Regards,
> Hubert
> 
> -- 
> -------------------------------------------------------------------------------
> Dr.-Ing. Hubert Fröhlich			
> Bezirksfinanzdirektion München 			
> Alexandrastr. 3, D-80538 München, GERMANY
> Tel. :+49 (0)89 / 2190 - 2980
> Fax  :+49 (0)89 / 2190 - 2997
> hubert dot froehlich at bvv dot bayern dot de
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users