[postgis-users] Problem with intersection from postgresql 8.0.7 to 8.1.3 win32

René F. Viancos S. rviancos at gmail.com
Sat Feb 18 07:30:04 PST 2006


i have a problem with the following query.

SELECT DISTINCT(intersection((SELECT collect(the_geom) FROM r13_ejes_32719
WHERE nombre = 'LOS RECUERDOS'),(SELECT collect(the_geom) FROM
r13_ejes_32719 WHERE nombre = 'LOS NOGALES'))) FROM r13_ejes_32719 WHERE
(nombre = 'LOS RECUERDOS' OR nombre = 'LOS NOGALES');

where 'r13_ejes_32719' is the table with the street data, 'nombre' the flied
with the street name.

In postgresql 8.0.3 works fine, but doesn't in postgresql 8.1.7 and i have
executed the postgis_full_version() in both versions;

postgresql 8.0.3, win32 binary package, has the folowing:
POSTGIS="0.9.1" GEOS="2.1.1" PROJ="Rel. 4.4.9, 29 Oct 2004" USE_STATS
DBPROC=" 0.0.1" RELPROC="0.0.1"

postgresql 8.1.7, win32 binary package, has the folowing;
POSTGIS="1.0.4" GEOS="2.1.4" PROJ="Rel. 4.4.9, 29 Oct 2004" USE_STATS
DBPROC=" 0.3.0" RELPROC="0.3.0"



In postgresql 8.0.3 the shape dumper creates this DDL

CREATE TABLE "public"."r13_ejes_32719" (
  "gid" SERIAL,
  "fnode_" BIGINT,
  "tnode_" BIGINT,
  "lpoly_" BIGINT,
  "rpoly_" BIGINT,
  "length" NUMERIC,
  "svial05_" BIGINT,
  "svial05_id" BIGINT,
  "iniizq" NUMERIC(20,0),
  "terizq" NUMERIC(20,0),
  "inider" NUMERIC(20,0),
  "terder" NUMERIC(20,0),
  "nombre" VARCHAR,
  "clase" VARCHAR,
  "prefijo" VARCHAR,
  "observ" VARCHAR,
  "transito" NUMERIC(20,0),
  "id_saf" NUMERIC(20,0),
  "the_geom" "public"."geometry",
  CONSTRAINT "r13_ejes_32719_pkey" PRIMARY KEY("gid"),
  CONSTRAINT "enforce_geotype_the_geom" CHECK ((geometrytype(the_geom) =
'MULTILINESTRING'::text) OR (the_geom IS NULL)),
  CONSTRAINT "enforce_srid_the_geom" CHECK (srid(the_geom) = 32719)
) WITH OIDS;


In postgresql 8.1.7 the shape dumper creates this DDL

CREATE TABLE "public"."r13_ejes_32719" (
  "gid" SERIAL,
  "fnode_" BIGINT,
  "tnode_" BIGINT,
  "lpoly_" BIGINT,
  "rpoly_" BIGINT,
  "length" NUMERIC,
  "svial05_" BIGINT,
  "svial05_id" BIGINT,
  "iniizq" NUMERIC(20,0),
  "terizq" NUMERIC(20,0),
  "inider" NUMERIC(20,0),
  "terder" NUMERIC(20,0),
  "nombre" VARCHAR,
  "clase" VARCHAR,
  "prefijo" VARCHAR,
  "observ" VARCHAR,
  "transito" NUMERIC(20,0),
  "id_saf" NUMERIC(20,0),
  "the_geom" "public"."geometry",
  CONSTRAINT "r13_ejes_32719_pkey" PRIMARY KEY("gid"),
  CONSTRAINT "enforce_dims_the_geom" CHECK (ndims(the_geom) = 2),
  CONSTRAINT "enforce_geotype_the_geom" CHECK ((geometrytype(the_geom) =
'MULTILINESTRING'::text) OR (the_geom IS NULL)),
  CONSTRAINT "enforce_srid_
the_geom" CHECK (srid(the_geom) = 32719)
) WITHOUT OIDS;


Finally, the output for the query in postgresql 8.0.3 is

SRID=-1;POINT( 355514.59375 6290622)     (the intersection point between
"LOS RECUERDOS" and "LOS NOGALES" streets)

and the output in postgresql 8.1.7 is

ERROR:  GEOS Intersection() threw an error! (i don't know why....)


2006/2/18, strk at refractions.net <strk at refractions.net>:
>
> On Sat, Feb 18, 2006 at 12:05:34PM -0300, René F. Viancos S. wrote:
> > Yes, the postgis_full_version() output is for the two cases,
>
> I mean, what is it ? Can you run postgis_full_version() on
> both setups and show it here ?
>
> --strk;
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>



--
René F. Viáncos S.
Director de Geomática y TIC
Vicerrectoría de Investigación y Desarrollo
Universidad de Chile
Tel (56-2) 632 62 09
Cel (56 9) 933 72 66
rviancos at uchile.cl
rviancos at gmail.com
www.investigacion.uchile.cl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20060218/64a434ef/attachment.html>


More information about the postgis-users mailing list