[postgis-users] ERROR: Operation on mixed SRID geometries after upgrading to recent postgres/postgis

Sebastian Reitenbach sebastia at l00-bugdead-prods.de
Fri Apr 4 06:32:21 PDT 2008


Hi all,

I'm in the phase of evaluating an upgrade of postgres/postgis from 
8.1.4/1.1.2 to 8.3.X/1.3.2.
But while testing, I ran into trouble with one of the queries that worked 
before

Now I have postgresql-8.3.1.tar.bz2 or postgresql-8.3.0.tar.bz2, and 
proj-4.6.0.tar.gz, postgis-1.3.2.tar.gz, and geos-3.0.0.tar.bz2 installed. 
With the old database/postgis version constellation (geos-2.2.2.tar.bz2, 
postgis-1.1.2.tar.gz, proj-4.4.9.tar.gz, and postgresql-8.1.4.tar.bz2) the 
query below worked:

select "TILE_ID", innerst_shape
from
(select original_tile."UNIQUE_ID", 
original_tile."TILE_ID", 
original_tile."TILE_INTERIOR_SHAPE", 
original_tile."TILE_POLY_SHAPE", 
difference(original_tile."TILE_INTERIOR_SHAPE", 
geomunion(intersect_tile."TILE_POLY_SHAPE")) as 
innerst_shape 
FROM "tile" original_tile 
JOIN "tile" intersect_tile 
ON (intersect_tile."TILE_ID"!= original_tile."TILE_ID" AND 
intersect_tile."TILE_POLY_SHAPE" && original_tile."TILE_INTERIOR_SHAPE" and 
Intersects(intersect_tile."TILE_POLY_SHAPE", 
original_tile."TILE_INTERIOR_SHAPE") and 
round(intersect_tile."TILE_ID"/100000)<=round(original_tile."TILE_ID"/100000)) 
WHERE original_tile."TILE_ID" in 
(select "TILE_ID" 
from "tile" 
where "TILE_POLY_SHAPE"&& GeomFromText('POLYGON((75.7491342 
31.0239757,75.7473798 30.7984119,76.0086717 30.7966414,76.0110403 
31.02219,75.7491342 31.0239757))', 4326) and 
Intersects ("TILE_POLY_SHAPE", GeomFromText('POLYGON((75.7491342 
31.0239757,75.7473798 30.7984119,76.0086717 30.7966414,76.0110403 
31.02219,75.7491342 31.0239757))', 4326)) ) 
GROUP BY original_tile."UNIQUE_ID", original_tile."TILE_ID", 
original_tile."TILE_INTERIOR_SHAPE", 
original_tile."TILE_POLY_SHAPE") as subquery 
where 
not IsEmpty(innerst_shape) and 
Intersects(innerst_shape, GeomFromText('POLYGON((75.7491342 
31.0239757,75.7473798 30.7984119,76.0086717 30.7966414,76.0110403 
31.02219,75.7491342 31.0239757))', 4326))
and not Touches(innerst_shape, GeomFromText('POLYGON((75.7491342 
31.0239757,75.7473798 30.7984119,76.0086717 30.7966414,76.0110403 
31.02219,75.7491342 31.0239757))', 4326))


ERROR:  Operation on mixed SRID geometries

********** Error **********

ERROR: Operation on mixed SRID geometries
SQL state: XX000

this worked with the old versin, without glitch. 
So anybody has an idea what this error message mean, and more importantly 
how to overcome it?

kind regards
Sebastian