[postgis-users] ST_contains() throws an error

[James DeMichele]

Here't he full version:

 

                                postgis_full_version

------------------------------------------------------------------------
-------------

 POSTGIS="1.3.1" GEOS="3.0.0rc4-CAPI-1.3.3" PROJ="Rel. 4.5.0, 22 Oct
2006" USE_STATS

 

-James

 

From: postgis-users-bounces at postgis.refractions.net
[mailto:postgis-users-bounces at postgis.refractions.net] On Behalf Of Obe,
Regina
Sent: Monday, December 31, 2007 3:46 AM
To: PostGIS Users Discussion
Subject: RE: [postgis-users] ST_contains() throws an error

 

Not quite enough information there to guess.

 

select  postgis_full_version();

 

instead of postgis_version()

 

That will give the exact Geos version you are running.

 

Hope that helps,

Regina

 

________________________________

From: postgis-users-bounces at postgis.refractions.net
[mailto:postgis-users-bounces at postgis.refractions.net] On Behalf Of
James DeMichele
Sent: Friday, December 28, 2007 2:25 PM
To: postgis-users at postgis.refractions.net
Subject: [postgis-users] ST_contains() throws an error

Hi,

I have a table with 27000 records and each record has 2 geometry
columns. The first column has an original MULTIPOLYGON (named
original_poly), and the second column has some tweaks (buffer and
simplification) applied to the original data, and is also stored as a
MULTIPOLYGON (named new_poly).

 

I want to get the list of these rows where the new_poly does NOT contain
the original_poly. The first time I ran it, I got this: "ERROR:  GEOS
contains() threw an error!"

 

So, then I thought that maybe some of the polygons were not valid. There
were 6 invalid polygons, so then I tried running this query:

 

select count(*) 

from temp_places_buffer 

where(st_isvalid(new_poly) = true and st_isvalid(redfin_poly) = true) 

and not st_contains(t1.new_poly, t1.redfin_poly);

 

Unfortunately, I still get this error:

 

ERROR:  GEOS contains() threw an error!

 

What are possible reasons that this would still be throwing an error,
other than invalid polygons, since the query should now only be
calculating st_contains on valid polygons?

 

I am running this version:

 

            postgis_version

---------------------------------------

 1.3 USE_GEOS=1 USE_PROJ=1 USE_STATS=1

 

Thanks.

 

-James

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20071231/d8c79d82/attachment.html>