self-hosted/postgis
self-hosted/postgis
1
0
Fork 0
mirror of https://git.osgeo.org/gitea/postgis/postgis synced 2024-10-25 17:42:38 +00:00
Code Issues Packages Projects Releases Wiki Activity
postgis/lwgeom/lwpostgis.sql.in
Sandro Santilli b01cca1a22 Enforced schema use in UpdateGeometrySRID 
git-svn-id: http://svn.osgeo.org/postgis/trunk@1151 b70326c6-7e19-0410-871a-916f4a2858ee
2004-12-15 12:59:40 +00:00

3319 lines
94 KiB
MySQL
Raw Blame History

-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
--
-- $Id$
--
-- PostGIS - Spatial Types for PostgreSQL
-- http://postgis.refractions.net
-- Copyright 2001-2003 Refractions Research Inc.
--
-- This is free software; you can redistribute and/or modify it under
-- the terms of the GNU General Public Licence. See the COPYING file.
--
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
#if USE_VERSION > 71
#define CREATEFUNCTION CREATE OR REPLACE FUNCTION
#else
#define CREATEFUNCTION CREATE FUNCTION
#endif
BEGIN;
-------------------------------------------------------------------
-- HISTOGRAM2D TYPE (lwhistogram2d)
-------------------------------------------------------------------
#if USE_VERSION < 73
# define HISTOGRAM_IN_REP opaque
# define HISTOGRAM_OUT_REP opaque
#else
# define HISTOGRAM_IN_REP histogram2d
# define HISTOGRAM_OUT_REP cstring
#endif
CREATEFUNCTION histogram2d_in(HISTOGRAM_OUT_REP)
RETURNS HISTOGRAM_IN_REP
AS '@MODULE_FILENAME@', 'lwhistogram2d_in'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION histogram2d_out(HISTOGRAM_IN_REP)
RETURNS HISTOGRAM_OUT_REP
AS '@MODULE_FILENAME@', 'lwhistogram2d_out'
LANGUAGE 'C' WITH (isstrict);
CREATE TYPE histogram2d (
alignment = double,
internallength = variable,
input = histogram2d_in,
output = histogram2d_out,
storage = main
);
-------------------------------------------------------------------
-- SPHEROID TYPE
-------------------------------------------------------------------
#if USE_VERSION < 73
# define SPHEROID_IN_REP opaque
# define SPHEROID_OUT_REP opaque
#else
# define SPHEROID_IN_REP spheroid
# define SPHEROID_OUT_REP cstring
#endif
CREATEFUNCTION spheroid_in(SPHEROID_OUT_REP)
RETURNS SPHEROID_IN_REP
AS '@MODULE_FILENAME@','ellipsoid_in'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION spheroid_out(SPHEROID_IN_REP)
RETURNS SPHEROID_OUT_REP
AS '@MODULE_FILENAME@','ellipsoid_out'
LANGUAGE 'C' WITH (isstrict);
CREATE TYPE spheroid (
alignment = double,
internallength = 65,
input = spheroid_in,
output = spheroid_out
);
-------------------------------------------------------------------
-- GEOMETRY TYPE (lwgeom)
-------------------------------------------------------------------
#if USE_VERSION < 73
# define GEOMETRY_IN_REP opaque
# define GEOMETRY_OUT_REP opaque
#else
# define GEOMETRY_IN_REP geometry
# define GEOMETRY_OUT_REP cstring
#endif
CREATEFUNCTION geometry_in(GEOMETRY_OUT_REP)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_in'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_out(GEOMETRY_IN_REP)
RETURNS cstring
AS '@MODULE_FILENAME@','LWGEOM_out'
LANGUAGE 'C' WITH (isstrict,iscachable);
#if USE_VERSION >= 80
CREATEFUNCTION geometry_analyze(internal)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_analyze'
LANGUAGE 'C' WITH (isstrict);
#endif
#if USE_VERSION > 73
CREATEFUNCTION geometry_recv(internal)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_recv'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION geometry_send(geometry)
RETURNS bytea
AS '@MODULE_FILENAME@','WKBFromLWGEOM'
LANGUAGE 'C' WITH (isstrict);
#endif
CREATE TYPE geometry (
internallength = variable,
input = geometry_in,
output = geometry_out,
#if USE_VERSION > 73
send = geometry_send,
receive = geometry_recv,
#endif
#if USE_VERSION >= 80
analyze = geometry_analyze,
#endif // USE_VERSION >= 80
storage = main
);
-------------------------------------------------------------------
-- BOX3D TYPE
-------------------------------------------------------------------
#if USE_VERSION < 73
# define BOX3D_IN_REP opaque
# define BOX3D_OUT_REP opaque
#else
# define BOX3D_IN_REP box3d
# define BOX3D_OUT_REP cstring
#endif
CREATEFUNCTION box3d_in(BOX3D_OUT_REP)
RETURNS BOX3D_IN_REP
AS '@MODULE_FILENAME@', 'BOX3D_in'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION box3d_out(BOX3D_IN_REP)
RETURNS BOX3D_OUT_REP
AS '@MODULE_FILENAME@', 'BOX3D_out'
LANGUAGE 'C' WITH (isstrict);
CREATE TYPE box3d (
alignment = double,
internallength = 48,
input = box3d_in,
output = box3d_out
);
CREATEFUNCTION xmin(box3d)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','BOX3D_xmin'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION ymin(box3d)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','BOX3D_ymin'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION zmin(box3d)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','BOX3D_zmin'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION xmax(box3d)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','BOX3D_xmax'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION ymax(box3d)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','BOX3D_ymax'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION zmax(box3d)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','BOX3D_zmax'
LANGUAGE 'C' WITH (isstrict,iscachable);
-------------------------------------------------------------------
-- CHIP TYPE
-------------------------------------------------------------------
#if USE_VERSION < 73
# define CHIP_IN_REP opaque
# define CHIP_OUT_REP opaque
#else
# define CHIP_IN_REP chip
# define CHIP_OUT_REP cstring
#endif
CREATEFUNCTION chip_in(CHIP_OUT_REP)
RETURNS CHIP_IN_REP
AS '@MODULE_FILENAME@','CHIP_in'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION chip_out(CHIP_IN_REP)
RETURNS CHIP_OUT_REP
AS '@MODULE_FILENAME@','CHIP_out'
LANGUAGE 'C' WITH (isstrict);
CREATE TYPE chip (
alignment = double,
internallength = variable,
input = chip_in,
output = chip_out,
storage = extended
);
-----------------------------------------------------------------------
-- BOX2D
-----------------------------------------------------------------------
#if USE_VERSION < 73
# define BOX2D_IN_REP opaque
# define BOX2D_OUT_REP opaque
#else
# define BOX2D_IN_REP box2d
# define BOX2D_OUT_REP cstring
#endif
CREATEFUNCTION box2d_in(BOX2D_OUT_REP)
RETURNS box2d
AS '@MODULE_FILENAME@','BOX2DFLOAT4_in'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_out(BOX2D_IN_REP)
RETURNS cstring
AS '@MODULE_FILENAME@','BOX2DFLOAT4_out'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATE TYPE box2d (
internallength = 16,
input = box2d_in,
output = box2d_out,
storage = plain
);
---- BOX2D support functions
CREATEFUNCTION box2d_overleft(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_overleft'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_overright(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_overright'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_left(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_left'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_right(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_right'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_contain(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_contain'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_contained(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_contained'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_overlap(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_overlap'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_same(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_same'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d_intersects(box2d, box2d)
RETURNS bool
AS '@MODULE_FILENAME@', 'BOX2D_intersects'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION xmin(box2d)
RETURNS FLOAT4
AS '@MODULE_FILENAME@','BOX2DFLOAT4_xmin'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION ymin(box2d)
RETURNS FLOAT4
AS '@MODULE_FILENAME@','BOX2DFLOAT4_ymin'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION xmax(box2d)
RETURNS FLOAT4
AS '@MODULE_FILENAME@','BOX2DFLOAT4_xmax'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION ymax(box2d)
RETURNS FLOAT4
AS '@MODULE_FILENAME@','BOX2DFLOAT4_ymax'
LANGUAGE 'C' WITH (isstrict,iscachable);
-- lwgeom operator support functions
-------------------------------------------------------------------
-- BTREE indexes
-------------------------------------------------------------------
CREATEFUNCTION geometry_lt(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'lwgeom_lt'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION geometry_le(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'lwgeom_le'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION geometry_gt(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'lwgeom_gt'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION geometry_ge(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'lwgeom_ge'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION geometry_eq(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'lwgeom_eq'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION geometry_cmp(geometry, geometry)
RETURNS integer
AS '@MODULE_FILENAME@', 'lwgeom_cmp'
LANGUAGE 'C' WITH (isstrict);
--
-- Sorting operators for Btree
--
CREATE OPERATOR < (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_lt,
COMMUTATOR = '>', NEGATOR = '>=',
RESTRICT = contsel, JOIN = contjoinsel
);
CREATE OPERATOR <= (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_le,
COMMUTATOR = '>=', NEGATOR = '>',
RESTRICT = contsel, JOIN = contjoinsel
);
CREATE OPERATOR = (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_eq,
COMMUTATOR = '=', -- we might implement a faster negator here
RESTRICT = contsel, JOIN = contjoinsel
);
CREATE OPERATOR >= (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_ge,
COMMUTATOR = '<=', NEGATOR = '<',
RESTRICT = contsel, JOIN = contjoinsel
);
CREATE OPERATOR > (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_gt,
COMMUTATOR = '<', NEGATOR = '<=',
RESTRICT = contsel, JOIN = contjoinsel
);
#if USE_VERSION >= 74
CREATE OPERATOR CLASS btree_geometry_ops
DEFAULT FOR TYPE geometry USING btree AS
OPERATOR 1 < ,
OPERATOR 2 <= ,
OPERATOR 3 = ,
OPERATOR 4 >= ,
OPERATOR 5 > ,
FUNCTION 1 geometry_cmp (geometry, geometry);
#endif // USE_VERSION >= 74
-------------------------------------------------------------------
-- GiST indexes
-------------------------------------------------------------------
#if USE_VERSION == 71
CREATEFUNCTION postgis_gist_sel(oid, oid, int2, opaque, int4)
#elif USE_VERSION == 72
CREATEFUNCTION postgis_gist_sel(opaque, oid, opaque, int4)
#else
CREATEFUNCTION postgis_gist_sel (internal, oid, internal, int4)
#endif
RETURNS float8
AS '@MODULE_FILENAME@', 'LWGEOM_gist_sel'
LANGUAGE 'C';
#if USE_VERSION == 71
CREATEFUNCTION postgis_gist_joinsel(oid, oid, int2, opaque, smallint)
#elif USE_VERSION == 72
CREATEFUNCTION postgis_gist_joinsel(opaque, oid, opaque, smallint)
#else
CREATEFUNCTION postgis_gist_joinsel(internal, oid, internal, smallint)
#endif
RETURNS float8
AS '@MODULE_FILENAME@', 'LWGEOM_gist_joinsel'
LANGUAGE 'C';
CREATEFUNCTION geometry_overleft(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_overleft'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_overright(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_overright'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_left(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_left'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_right(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_right'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_contain(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_contain'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_contained(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_contained'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_overlap(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_overlap'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry_same(geometry, geometry)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_same'
LANGUAGE 'C' WITH (isstrict,iscachable);
-- GEOMETRY operators
CREATE OPERATOR << (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_left,
COMMUTATOR = '>>',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR &< (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_overleft,
COMMUTATOR = '&>',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR && (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_overlap,
COMMUTATOR = '&&',
RESTRICT = postgis_gist_sel, JOIN = postgis_gist_joinsel
);
CREATE OPERATOR &> (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_overright,
COMMUTATOR = '&<',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR >> (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_right,
COMMUTATOR = '<<',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR ~= (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_same,
COMMUTATOR = '~=',
RESTRICT = eqsel, JOIN = eqjoinsel
);
CREATE OPERATOR @ (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_contained,
COMMUTATOR = '~',
RESTRICT = contsel, JOIN = contjoinsel
);
CREATE OPERATOR ~ (
LEFTARG = geometry, RIGHTARG = geometry, PROCEDURE = geometry_contain,
COMMUTATOR = '@',
RESTRICT = contsel, JOIN = contjoinsel
);
-- BOX2D operators
CREATE OPERATOR << (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_left,
COMMUTATOR = '>>',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR &< (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_overleft,
COMMUTATOR = '&>',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR && (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_overlap,
COMMUTATOR = '&&',
RESTRICT = postgis_gist_sel, JOIN = postgis_gist_joinsel
);
CREATE OPERATOR &> (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_overright,
COMMUTATOR = '&<',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR >> (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_right,
COMMUTATOR = '<<',
RESTRICT = positionsel, JOIN = positionjoinsel
);
CREATE OPERATOR ~= (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_same,
COMMUTATOR = '~=',
RESTRICT = eqsel, JOIN = eqjoinsel
);
CREATE OPERATOR @ (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_contained,
COMMUTATOR = '~',
RESTRICT = contsel, JOIN = contjoinsel
);
CREATE OPERATOR ~ (
LEFTARG = box2d, RIGHTARG = box2d, PROCEDURE = box2d_contain,
COMMUTATOR = '@',
RESTRICT = contsel, JOIN = contjoinsel
);
-- gist support functions
#if USE_VERSION < 73
#define OPAQUE_TYPE opaque
#else
#define OPAQUE_TYPE internal
#endif
CREATEFUNCTION LWGEOM_gist_consistent(OPAQUE_TYPE,geometry,int4)
RETURNS bool
AS '@MODULE_FILENAME@' ,'LWGEOM_gist_consistent'
LANGUAGE 'C';
CREATEFUNCTION LWGEOM_gist_compress(OPAQUE_TYPE)
RETURNS OPAQUE_TYPE
AS '@MODULE_FILENAME@','LWGEOM_gist_compress'
LANGUAGE 'C';
CREATEFUNCTION LWGEOM_gist_penalty(OPAQUE_TYPE,OPAQUE_TYPE,OPAQUE_TYPE)
RETURNS OPAQUE_TYPE
AS '@MODULE_FILENAME@' ,'LWGEOM_gist_penalty'
LANGUAGE 'C';
CREATEFUNCTION LWGEOM_gist_picksplit(OPAQUE_TYPE, OPAQUE_TYPE)
RETURNS OPAQUE_TYPE
AS '@MODULE_FILENAME@' ,'LWGEOM_gist_picksplit'
LANGUAGE 'C';
CREATEFUNCTION LWGEOM_gist_union(bytea, OPAQUE_TYPE)
RETURNS OPAQUE_TYPE
AS '@MODULE_FILENAME@' ,'LWGEOM_gist_union'
LANGUAGE 'C';
CREATEFUNCTION LWGEOM_gist_same(box2d, box2d, OPAQUE_TYPE)
RETURNS OPAQUE_TYPE
AS '@MODULE_FILENAME@' ,'LWGEOM_gist_same'
LANGUAGE 'C';
CREATEFUNCTION LWGEOM_gist_decompress(OPAQUE_TYPE)
RETURNS OPAQUE_TYPE
AS '@MODULE_FILENAME@' ,'LWGEOM_gist_decompress'
LANGUAGE 'C';
-------------------------------------------
-- GIST opclass index binding entries.
-------------------------------------------
#if USE_VERSION == 71
INSERT INTO pg_opclass (opcname, opcdeftype)
SELECT 'gist_geometry_ops', oid
FROM pg_type
WHERE typname = 'geometry';
SELECT o.oid AS opoid, o.oprname
INTO TABLE rt_ops_tmp
FROM pg_operator o, pg_type t
WHERE o.oprleft = t.oid AND t.typname = 'geometry';
-- box_left
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 1
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '<<';
-- box_overleft
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 2
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '&<';
-- box_overlap
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 3
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '&&';
-- box_overright
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 4
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '&>';
-- box_right
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 5
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '>>';
-- box_same
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 6
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '~=';
-- box_contains
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 7
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '~';
-- box_contained
INSERT INTO pg_amop (amopid, amopclaid, amopopr, amopstrategy)
SELECT am.oid, opcl.oid, c.opoid, 8
FROM pg_am am, pg_opclass opcl, rt_ops_tmp c
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND c.oprname = '@';
DROP TABLE rt_ops_tmp;
--
-- Add the entries to amproc for the support methods.
-- Note the amprocnum numbers associated with each are specific!
--
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 1
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_consistent';
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 2
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_union';
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 3
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_compress';
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 4
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_decompress';
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 5
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_penalty';
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 6
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_picksplit';
INSERT INTO pg_amproc (amid, amopclaid, amproc, amprocnum)
SELECT am.oid, opcl.oid, pro.oid, 7
FROM pg_am am, pg_opclass opcl, pg_proc pro
WHERE amname = 'gist' AND opcname = 'gist_geometry_ops'
AND proname = 'LWGEOM_gist_same';
#else
#if USE_VERSION == 72
--
-- Create opclass index binding entries for PG72
--
INSERT INTO pg_opclass (opcamid, opcname, opcintype, opcdefault, opckeytype)
VALUES (
(SELECT oid FROM pg_am WHERE amname = 'gist'),
'gist_geometry_ops',
(SELECT oid FROM pg_type WHERE typname = 'geometry'),
true,
(SELECT oid FROM pg_type WHERE typname = 'box2d'));
-- drop table rt_ops_tmp;
SELECT o.oid AS opoid, o.oprname
INTO TABLE rt_ops_tmp
FROM pg_operator o, pg_type t
WHERE o.oprleft = t.oid
AND t.typname = 'geometry';
-- poly_left
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 1, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '<<';
-- poly_overleft
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 2, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '&<';
-- poly_overlap
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 3, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '&&';
-- poly_overright
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 4, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '&>';
-- poly_right
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 5, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '>>';
-- poly_same
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 6, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '~=';
-- poly_contains
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 7, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '~';
-- poly_contained
INSERT INTO pg_amop (amopclaid, amopstrategy, amopreqcheck, amopopr)
SELECT opcl.oid, 8, true, c.opoid
FROM pg_opclass opcl, rt_ops_tmp c
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and c.oprname = '@';
DROP TABLE rt_ops_tmp;
-- add the entries to amproc for the support methods
-- note the amprocnum numbers associated with each are specific!
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 1, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_consistent';
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 2, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_union';
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 3, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_compress';
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 4, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_decompress';
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 5, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_penalty';
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 6, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_picksplit';
INSERT INTO pg_amproc (amopclaid, amprocnum, amproc)
SELECT opcl.oid, 7, pro.oid
FROM pg_opclass opcl, pg_proc pro
WHERE
opcamid = (SELECT oid FROM pg_am WHERE amname = 'gist')
and opcname = 'gist_geometry_ops'
and proname = 'LWGEOM_gist_same';
#else // USE_VERSION >= 73
--
-- Create opclass index bindings for PG>=73
--
CREATE OPERATOR CLASS gist_geometry_ops
DEFAULT FOR TYPE geometry USING gist AS
OPERATOR 1 << RECHECK,
OPERATOR 2 &< RECHECK,
OPERATOR 3 && RECHECK,
OPERATOR 4 &> RECHECK,
OPERATOR 5 >> RECHECK,
OPERATOR 6 ~= RECHECK,
OPERATOR 7 ~ RECHECK,
OPERATOR 8 @ RECHECK,
FUNCTION 1 LWGEOM_gist_consistent (internal, geometry, int4),
FUNCTION 2 LWGEOM_gist_union (bytea, internal),
FUNCTION 3 LWGEOM_gist_compress (internal),
FUNCTION 4 LWGEOM_gist_decompress (internal),
FUNCTION 5 LWGEOM_gist_penalty (internal, internal, internal),
FUNCTION 6 LWGEOM_gist_picksplit (internal, internal),
FUNCTION 7 LWGEOM_gist_same (box2d, box2d, internal);
UPDATE pg_opclass
SET opckeytype = (select oid from pg_type where typname = 'box2d')
WHERE opcname = 'gist_geometry_ops';
-- TODO: add btree binding...
#endif // USE_VERSION >= 73
#endif // USE_VERSION != 71
-- other lwgeom functions
CREATEFUNCTION addBBOX(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_addBBOX'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION getSRID(geometry)
RETURNS int4
AS '@MODULE_FILENAME@','LWGEOM_getSRID'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION getbbox(geometry)
RETURNS box2d
AS '@MODULE_FILENAME@','LWGEOM_to_BOX2DFLOAT4'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION translate(geometry,float8,float8,float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_translate'
LANGUAGE 'C' WITH (isstrict) ;
CREATEFUNCTION translate(geometry,float8,float8)
RETURNS geometry
AS '
SELECT translate($1, $2, $3, 0)
' LANGUAGE 'SQL' WITH (isstrict);
--- CHIP functions
CREATEFUNCTION srid(chip)
RETURNS int4
AS '@MODULE_FILENAME@','CHIP_getSRID'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION height(chip)
RETURNS int4
AS '@MODULE_FILENAME@','CHIP_getHeight'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION factor(chip)
RETURNS FLOAT4
AS '@MODULE_FILENAME@','CHIP_getFactor'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION width(chip)
RETURNS int4
AS '@MODULE_FILENAME@','CHIP_getWidth'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION datatype(chip)
RETURNS int4
AS '@MODULE_FILENAME@','CHIP_getDatatype'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION compression(chip)
RETURNS int4
AS '@MODULE_FILENAME@','CHIP_getCompression'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION setSRID(chip,int4)
RETURNS chip
AS '@MODULE_FILENAME@','CHIP_setSRID'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION setFactor(chip,float4)
RETURNS chip
AS '@MODULE_FILENAME@','CHIP_setFactor'
LANGUAGE 'C' WITH (isstrict,iscachable);
------------------------------------------------------------------------
-- DEBUG
------------------------------------------------------------------------
CREATEFUNCTION mem_size(geometry)
RETURNS int4
AS '@MODULE_FILENAME@', 'LWGEOM_mem_size'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION summary(geometry)
RETURNS text
AS '@MODULE_FILENAME@', 'LWGEOM_summary'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION npoints(geometry)
RETURNS int4
AS '@MODULE_FILENAME@', 'LWGEOM_npoints'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION nrings(geometry)
RETURNS int4
AS '@MODULE_FILENAME@', 'LWGEOM_nrings'
LANGUAGE 'C' WITH (isstrict);
------------------------------------------------------------------------
-- OGC defined
------------------------------------------------------------------------
CREATEFUNCTION NumPoints(geometry)
RETURNS int4
AS '@MODULE_FILENAME@', 'LWGEOM_numpoints_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION NumGeometries(geometry)
RETURNS int4
AS '@MODULE_FILENAME@', 'LWGEOM_numgeometries_collection'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION GeometryN(geometry,integer)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_geometryn_collection'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION Dimension(geometry)
RETURNS int4
AS '@MODULE_FILENAME@', 'LWGEOM_dimension'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION ExterioRring(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_exteriorring_polygon'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION NumInteriorRings(geometry)
RETURNS integer
AS '@MODULE_FILENAME@','LWGEOM_numinteriorrings_polygon'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION InteriorRingN(geometry,integer)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_interiorringn_polygon'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION GeometryType(geometry)
RETURNS text
AS '@MODULE_FILENAME@', 'LWGEOM_getTYPE'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION PointN(geometry,integer)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_pointn_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION X(geometry)
RETURNS float8
AS '@MODULE_FILENAME@','LWGEOM_x_point'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION Y(geometry)
RETURNS float8
AS '@MODULE_FILENAME@','LWGEOM_y_point'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION Z(geometry)
RETURNS float8
AS '@MODULE_FILENAME@','LWGEOM_z_point'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION StartPoint(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_startpoint_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION EndPoint(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_endpoint_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION IsClosed(geometry)
RETURNS boolean
AS '@MODULE_FILENAME@', 'LWGEOM_isclosed_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION IsEmpty(geometry)
RETURNS boolean
AS '@MODULE_FILENAME@', 'LWGEOM_isempty'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION SRID(geometry)
RETURNS int4
AS '@MODULE_FILENAME@','LWGEOM_getSRID'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION SetSRID(geometry,int4)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_setSRID'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION AsBinary(geometry)
RETURNS bytea
AS '@MODULE_FILENAME@','WKBFromLWGEOM'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION AsBinary(geometry,text)
RETURNS bytea
AS '@MODULE_FILENAME@','WKBFromLWGEOM'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION AsText(geometry)
RETURNS TEXT
AS '@MODULE_FILENAME@','LWGEOM_asText'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION GeometryFromText(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_from_text'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION GeometryFromText(geometry, int4)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOM_from_text'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION GeomFromText(geometry)
RETURNS geometry AS '
SELECT geometryfromtext($1)
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION GeomFromText(geometry, int4)
RETURNS geometry AS '
SELECT geometryfromtext($1, $2)
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PointFromText(geometry)
RETURNS geometry
AS ' SELECT CASE WHEN geometrytype($1) = ''POINT'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PointFromText(geometry, int4)
RETURNS geometry
AS ' SELECT CASE WHEN geometrytype($1) = ''POINT'' THEN setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LineFromText(geometry)
RETURNS geometry
AS ' SELECT CASE WHEN geometrytype($1) = ''LINESTRING'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LineFromText(geometry, int4)
RETURNS geometry
AS ' SELECT CASE WHEN geometrytype($1) = ''LINESTRING'' THEN setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LineStringFromText(geometry)
RETURNS geometry
AS 'SELECT LineFromText($1)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LineStringFromText(geometry, int4)
RETURNS geometry
AS 'SELECT LineFromText($1, $2)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolyFromText(geometry)
RETURNS geometry
AS ' SELECT CASE WHEN geometrytype($1) = ''POLYGON'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolyFromText(geometry, int4)
RETURNS geometry
AS ' SELECT CASE WHEN geometrytype($1) = ''POLYGON'' THEN setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolygonFromText(geometry, int4)
RETURNS geometry
AS 'SELECT PolyFromText($1, $2)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolygonFromText(geometry)
RETURNS geometry
AS 'SELECT PolyFromText($1)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MLineFromText(geometry, int4)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''MULTILINESTRING'' THEN setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MLineFromText(geometry)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''MULTILINESTRING'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiLineStringFromText(geometry)
RETURNS geometry
AS 'SELECT MLineFromText($1)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiLineStringFromText(geometry, int4)
RETURNS geometry
AS 'SELECT MLineFromText($1, $2)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPointFromText(geometry, int4)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''MULTIPOINT'' THEN setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPointFromText(geometry)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''MULTIPOINT'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPointFromText(geometry, int4)
RETURNS geometry
AS 'SELECT MPointFromText($1, $2)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPointFromText(geometry)
RETURNS geometry
AS 'SELECT MPointFromText($1)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPolyFromText(geometry, int4)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''MULTIPOLYGON'' THEN setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPolyFromText(geometry)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''MULTIPOLYGON'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPolygonFromText(geometry, int4)
RETURNS geometry
AS 'SELECT MPolyFromText($1, $2)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPolygonFromText(geometry)
RETURNS geometry
AS 'SELECT MPolyFromText($1)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION GeomCollFromText(geometry, int4)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''GEOMETRYCOLLECTION'' THEN
setSRID($1,$2)
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION GeomCollFromText(geometry)
RETURNS geometry
AS 'SELECT CASE WHEN geometrytype($1) = ''GEOMETRYCOLLECTION'' THEN $1
ELSE NULL END
' LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION GeomFromWKB(bytea)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOMFromWKB'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION GeomFromWKB(bytea, int)
RETURNS geometry
AS 'SELECT setSRID(GeomFromWKB($1), $2)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PointFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''POINT'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PointFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''POINT'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LineFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''LINESTRING'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LineFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''LINESTRING'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LinestringFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''LINESTRING'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION LinestringFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''LINESTRING'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolyFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''POLYGON'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolyFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''POLYGON'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolygonFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''POLYGON'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION PolygonFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''POLYGON'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPointFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOINT'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPointFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOINT'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPointFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOINT'' THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPointFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOINT'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiLineFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTILINESTRING''
THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiLineFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTILINESTRING'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MLineFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTILINESTRING''
THEN setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MLineFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTILINESTRING'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPolyFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOLYGON'' THEN
setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MPolyFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOLYGON'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPolyFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOLYGON'' THEN
setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION MultiPolyFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''MULTIPOLYGON'' THEN $1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION GeomCollFromWKB(geometry, int)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''GEOMETRYCOLLECTION'' THEN
setSRID($1, $2)
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION GeomCollFromWKB(geometry)
RETURNS geometry
AS '
SELECT CASE WHEN geometrytype($1) = ''GEOMETRYCOLLECTION'' THEN
$1
ELSE NULL END
'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
------------------------------------------------------------------------
-- Misures
------------------------------------------------------------------------
-- this is a fake (for back-compatibility)
-- uses 3d if 3d is available, 2d otherwise
CREATEFUNCTION length3d(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_length_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION length2d(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_length2d_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION length(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_length_linestring'
LANGUAGE 'C' WITH (isstrict);
-- this is a fake (for back-compatibility)
-- uses 3d if 3d is available, 2d otherwise
CREATEFUNCTION length3d_spheroid(geometry, spheroid)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','LWGEOM_length_ellipsoid_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION length_spheroid(geometry, spheroid)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','LWGEOM_length_ellipsoid_linestring'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION length2d_spheroid(geometry, spheroid)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','LWGEOM_length2d_ellipsoid_linestring'
LANGUAGE 'C' WITH (isstrict);
-- this is a fake (for back-compatibility)
-- uses 3d if 3d is available, 2d otherwise
CREATEFUNCTION perimeter3d(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_perimeter_poly'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION perimeter2d(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_perimeter2d_poly'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION perimeter(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_perimeter_poly'
LANGUAGE 'C' WITH (isstrict);
-- this is an alias for 'area(geometry)'
-- there is nothing such an 'area3d'...
CREATEFUNCTION area2d(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@', 'LWGEOM_area_polygon'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION area(geometry)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','LWGEOM_area_polygon'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION distance_spheroid(geometry,geometry,spheroid)
RETURNS FLOAT8
AS '@MODULE_FILENAME@','LWGEOM_distance_ellipsoid_point'
LANGUAGE 'C' WITH (isstrict);
-- Minimum distance. 2d only.
CREATEFUNCTION distance(geometry,geometry)
RETURNS float8
AS '@MODULE_FILENAME@', 'LWGEOM_mindistance2d'
LANGUAGE 'C' WITH (isstrict,iscachable);
-- Maximum distance between linestrings. 2d only. Very bogus.
CREATEFUNCTION max_distance(geometry,geometry)
RETURNS float8
AS '@MODULE_FILENAME@', 'LWGEOM_maxdistance2d_linestring'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION point_inside_circle(geometry,float8,float8,float8)
RETURNS bool
AS '@MODULE_FILENAME@', 'LWGEOM_inside_circle_point'
LANGUAGE 'C' WITH (isstrict);
------------------------------------------------------------------------
-- MISC
------------------------------------------------------------------------
CREATEFUNCTION force_2d(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_2d'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION force_3dz(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_3dz'
LANGUAGE 'C' WITH (isstrict);
-- an alias for force_3dz
CREATEFUNCTION force_3d(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_3dz'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION force_3dm(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_3dm'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION force_4d(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_4d'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION force_collection(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_collection'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION multi(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_force_multi'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION collector(geometry, geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_collect'
LANGUAGE 'C';
CREATEFUNCTION collect(geometry, geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_collect'
LANGUAGE 'C';
CREATE AGGREGATE memcollect(
sfunc = collect,
basetype = geometry,
stype = geometry
);
CREATEFUNCTION geom_accum (geometry[],geometry)
RETURNS geometry[]
AS '@MODULE_FILENAME@', 'LWGEOM_accum'
LANGUAGE 'C';
CREATEFUNCTION collect_garray (geometry[])
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_collect_garray'
LANGUAGE 'C';
CREATE AGGREGATE collect (
sfunc = geom_accum,
basetype = geometry,
stype = geometry[],
finalfunc = collect_garray
);
CREATEFUNCTION expand(box3d,float8)
RETURNS box3d
AS '@MODULE_FILENAME@', 'BOX3D_expand'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION expand(box2d,float8)
RETURNS box2d
AS '@MODULE_FILENAME@', 'BOX2DFLOAT4_expand'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION expand(geometry,float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_expand'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION envelope(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_envelope'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION reverse(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_reverse'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION ForceRHR(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_forceRHR_poly'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION noop(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_noop'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION zmflag(geometry)
RETURNS smallint
AS '@MODULE_FILENAME@', 'LWGEOM_zmflag'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION ndims(geometry)
RETURNS smallint
AS '@MODULE_FILENAME@', 'LWGEOM_ndims'
LANGUAGE 'C' WITH (iscachable,isstrict);
------------------------------------------------------------------------
-- CONSTRUCTORS
------------------------------------------------------------------------
CREATEFUNCTION makePoint(float8, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_makepoint'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION makePoint(float8, float8, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_makepoint'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION makePoint(float8, float8, float8, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_makepoint'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION makePointM(float8, float8, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_makepoint3dm'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION makeBox2d(geometry, geometry)
RETURNS box2d
AS '@MODULE_FILENAME@', 'BOX2DFLOAT4_construct'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION makeBox3d(geometry, geometry)
RETURNS box3d
AS '@MODULE_FILENAME@', 'BOX3D_construct'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION makeline_garray (geometry[])
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_makeline_garray'
LANGUAGE 'C';
CREATEFUNCTION LineFromMultiPoint(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_line_from_mpoint'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION MakeLine(geometry, geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_makeline'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION AddPoint(geometry, geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_addpoint'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATEFUNCTION AddPoint(geometry, geometry, integer)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_addpoint'
LANGUAGE 'C' WITH (iscachable,isstrict);
CREATE AGGREGATE makeline (
sfunc = geom_accum,
basetype = geometry,
stype = geometry[],
finalfunc = makeline_garray
);
CREATEFUNCTION polygonize_garray (geometry[])
RETURNS geometry
AS '@MODULE_FILENAME@', 'GEOS_polygonize_garray'
LANGUAGE 'C';
CREATE AGGREGATE polygonize (
sfunc = geom_accum,
basetype = geometry,
stype = geometry[],
finalfunc = polygonize_garray
);
------------------------------------------------------------------------
--
-- Aggregate functions
--
CREATEFUNCTION combine_bbox(box2d,geometry)
RETURNS box2d
AS '@MODULE_FILENAME@', 'BOX2DFLOAT4_combine'
LANGUAGE 'C';
CREATE AGGREGATE extent(
sfunc = combine_bbox,
basetype = geometry,
stype = box2d
);
CREATEFUNCTION combine_bbox(box3d,geometry)
RETURNS box3d
AS '@MODULE_FILENAME@', 'BOX3D_combine'
LANGUAGE 'C';
CREATE AGGREGATE extent3d(
sfunc = combine_bbox,
basetype = geometry,
stype = box3d
);
-----------------------------------------------------------------------
-- CREATE_HISTOGRAM2D( <box2d>, <size> )
-----------------------------------------------------------------------
--
-- Returns a histgram with 0s in all the boxes.
--
-----------------------------------------------------------------------
CREATEFUNCTION create_histogram2d(box2d,int)
RETURNS histogram2d
AS '@MODULE_FILENAME@','create_lwhistogram2d'
LANGUAGE 'C' with (isstrict);
-----------------------------------------------------------------------
-- BUILD_HISTOGRAM2D( <histogram2d>, <tablename>, <columnname> )
-----------------------------------------------------------------------
CREATEFUNCTION build_histogram2d (histogram2d,text,text)
RETURNS histogram2d
AS '@MODULE_FILENAME@','build_lwhistogram2d'
LANGUAGE 'C' with (isstrict);
#if USE_VERSION >= 73
-----------------------------------------------------------------------
-- BUILD_HISTOGRAM2D(<histogram2d>,<schema>,<tablename>,<columnname>)
-----------------------------------------------------------------------
-- This is a wrapper to the omonimous schema unaware function,
-- thanks to Carl Anderson for the idea.
-----------------------------------------------------------------------
CREATEFUNCTION build_histogram2d (histogram2d,text,text,text)
RETURNS histogram2d
AS '
BEGIN
EXECUTE ''SET local search_path = ''||$2||'',public'';
RETURN public.build_histogram2d($1,$3,$4);
END
' LANGUAGE 'plpgsql' with (isstrict);
#endif // USE_VERSION >= 73
-----------------------------------------------------------------------
-- EXPLODE_HISTOGRAM2D( <histogram2d>, <tablename> )
-----------------------------------------------------------------------
CREATEFUNCTION explode_histogram2d (histogram2d,text)
RETURNS histogram2d
AS '@MODULE_FILENAME@','explode_lwhistogram2d'
LANGUAGE 'C' with (isstrict);
-----------------------------------------------------------------------
-- ESTIMATE_HISTOGRAM2D( <histogram2d>, <box> )
-----------------------------------------------------------------------
CREATEFUNCTION estimate_histogram2d(histogram2d,box2d)
RETURNS float8
AS '@MODULE_FILENAME@','estimate_lwhistogram2d'
LANGUAGE 'C' with (isstrict);
#if USE_VERSION >= 80
-----------------------------------------------------------------------
-- ESTIMATED_EXTENT( <schema name>, <table name>, <column name> )
-----------------------------------------------------------------------
CREATEFUNCTION estimated_extent(text,text,text) RETURNS box2d AS
'@MODULE_FILENAME@', 'LWGEOM_estimated_extent'
LANGUAGE 'C' with (isstrict);
-----------------------------------------------------------------------
-- ESTIMATED_EXTENT( <table name>, <column name> )
-----------------------------------------------------------------------
CREATEFUNCTION estimated_extent(text,text) RETURNS box2d AS
'SELECT estimated_extent(current_schema()::text, $1, $2)'
LANGUAGE 'SQL' WITH (isstrict);
#endif // USE_VERSION >= 80
-----------------------------------------------------------------------
-- FIND_EXTENT( <schema name>, <table name>, <column name> )
-----------------------------------------------------------------------
CREATEFUNCTION find_extent(text,text,text) RETURNS box2d AS
'
DECLARE
schemaname alias for $1;
tablename alias for $2;
columnname alias for $3;
myrec RECORD;
BEGIN
FOR myrec IN EXECUTE ''SELECT extent("''||columnname||''") FROM "''||schemaname||''"."''||tablename||''"'' LOOP
return myrec.extent;
END LOOP;
END;
'
LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- FIND_EXTENT( <table name>, <column name> )
-----------------------------------------------------------------------
CREATEFUNCTION find_extent(text,text) RETURNS box2d AS
'
DECLARE
tablename alias for $1;
columnname alias for $2;
myrec RECORD;
BEGIN
FOR myrec IN EXECUTE ''SELECT extent("''||columnname||''") FROM "''||tablename||''"'' LOOP
return myrec.extent;
END LOOP;
END;
'
LANGUAGE 'plpgsql' WITH (isstrict);
-------------------------------------------------------------------
-- SPATIAL_REF_SYS
-------------------------------------------------------------------
CREATE TABLE spatial_ref_sys (
srid integer not null primary key,
auth_name varchar(256),
auth_srid integer,
srtext varchar(2048),
proj4text varchar(2048)
);
-------------------------------------------------------------------
-- GEOMETRY_COLUMNS
-------------------------------------------------------------------
CREATE TABLE geometry_columns (
f_table_catalog varchar(256) not null,
f_table_schema varchar(256) not null,
f_table_name varchar(256) not null,
f_geometry_column varchar(256) not null,
coord_dimension integer not null,
srid integer not null,
type varchar(30) not null,
#if USE_VERSION < 80
attrelid oid,
varattnum int,
stats histogram2d,
#endif
CONSTRAINT geometry_columns_pk primary key (
f_table_catalog,
f_table_schema,
f_table_name,
f_geometry_column ) );
-----------------------------------------------------------------------
-- RENAME_GEOMETRY_TABLE_CONSTRAINTS()
-----------------------------------------------------------------------
-- This function has been obsoleted for the difficulty in
-- finding attribute on which the constraint is applied.
-- AddGeometryColumn will name the constraints in a meaningful
-- way, but nobody can rely on it since old postgis versions did
-- not do that.
-----------------------------------------------------------------------
CREATEFUNCTION rename_geometry_table_constraints() RETURNS text
AS
'
SELECT ''rename_geometry_table_constraint() is obsoleted''::text
' LANGUAGE 'SQL';
-----------------------------------------------------------------------
-- FIX_GEOMETRY_COLUMNS()
-----------------------------------------------------------------------
-- This function will:
--
-- o try to fix the schema of records with an invalid one
-- (for PG>=73)
--
-- o link records to system tables through attrelid and varattnum
-- (for PG<75)
--
-- o delete all records for which no linking was possible
-- (for PG<75)
--
--
-----------------------------------------------------------------------
CREATEFUNCTION fix_geometry_columns() RETURNS text
AS
'
DECLARE
mislinked record;
result text;
linked integer;
deleted integer;
#if USE_VERSION >= 73
foundschema integer;
#endif
BEGIN
#if USE_VERSION >= 73
-- Since 7.3 schema support has been added.
-- Previous postgis versions used to put the database name in
-- the schema column. This needs to be fixed, so we try to
-- set the correct schema for each geometry_colums record
-- looking at table, column, type and srid.
UPDATE geometry_columns SET f_table_schema = n.nspname
FROM pg_namespace n, pg_class c, pg_attribute a,
pg_constraint sridcheck, pg_constraint typecheck
WHERE ( f_table_schema is NULL
OR f_table_schema = ''''
OR f_table_schema NOT IN (
SELECT nspname::varchar
FROM pg_namespace nn, pg_class cc, pg_attribute aa
WHERE cc.relnamespace = nn.oid
AND cc.relname = f_table_name::name
AND aa.attrelid = cc.oid
AND aa.attname = f_geometry_column::name))
AND f_table_name::name = c.relname
AND c.oid = a.attrelid
AND c.relnamespace = n.oid
AND f_geometry_column::name = a.attname
AND sridcheck.conrelid = c.oid
--AND sridcheck.conname = ''$1''
AND sridcheck.consrc LIKE ''(srid(% = %)''
AND typecheck.conrelid = c.oid
--AND typecheck.conname = ''$2''
AND typecheck.consrc LIKE
''((geometrytype(%) = ''''%''''::text) OR (% IS NULL))''
AND sridcheck.consrc ~ textcat('' = '', srid::text)
AND typecheck.consrc ~ textcat('' = '''''', type::text)
AND NOT EXISTS (
SELECT oid FROM geometry_columns gc
WHERE c.relname::varchar = gc.f_table_name
#if USE_VERSION >= 73
AND n.nspname::varchar = gc.f_table_schema
#endif
AND a.attname::varchar = gc.f_geometry_column
);
GET DIAGNOSTICS foundschema = ROW_COUNT;
#endif
#if USE_VERSION >= 80
-- no linkage to system table needed
return ''fixed:''||foundschema::text;
#endif
-- fix linking to system tables
SELECT 0 INTO linked;
FOR mislinked in
SELECT gc.oid as gcrec,
a.attrelid as attrelid, a.attnum as attnum
FROM geometry_columns gc, pg_class c,
#if USE_VERSION >= 73
pg_namespace n, pg_attribute a
#else
pg_attribute a
#endif
WHERE ( gc.attrelid IS NULL OR gc.attrelid != a.attrelid
OR gc.varattnum IS NULL OR gc.varattnum != a.attnum)
#if USE_VERSION >= 73
AND n.nspname = gc.f_table_schema::name
AND c.relnamespace = n.oid
#endif
AND c.relname = gc.f_table_name::name
AND a.attname = f_geometry_column::name
AND a.attrelid = c.oid
LOOP
UPDATE geometry_columns SET
attrelid = mislinked.attrelid,
varattnum = mislinked.attnum,
stats = NULL
WHERE geometry_columns.oid = mislinked.gcrec;
SELECT linked+1 INTO linked;
END LOOP;
-- remove stale records
DELETE FROM geometry_columns WHERE attrelid IS NULL;
GET DIAGNOSTICS deleted = ROW_COUNT;
result =
#if USE_VERSION >= 73
''fixed:'' || foundschema::text ||
#endif
'' linked:'' || linked::text ||
'' deleted:'' || deleted::text;
return result;
END;
'
LANGUAGE 'plpgsql' ;
-----------------------------------------------------------------------
-- PROBE_GEOMETRY_COLUMNS()
-----------------------------------------------------------------------
-- Fill the geometry_columns table with values probed from the system
-- catalogues. 3d flag can not be probed, it defaults to 2
--
-- Note that bogus records already in geometry_columns are not
-- overridden (a check for schema.table.column is performed), so
-- to have a fresh probe backup your geometry_column, delete from
-- it and probe.
-----------------------------------------------------------------------
CREATEFUNCTION probe_geometry_columns() RETURNS text AS
'
DECLARE
inserted integer;
oldcount integer;
probed integer;
stale integer;
BEGIN
SELECT count(*) INTO oldcount FROM geometry_columns;
SELECT count(*) INTO probed
FROM pg_class c, pg_attribute a, pg_type t,
#if USE_VERSION >= 73
pg_namespace n,
#endif
pg_constraint sridcheck, pg_constraint typecheck
WHERE t.typname = ''geometry''
AND a.atttypid = t.oid
AND a.attrelid = c.oid
#if USE_VERSION >= 73
AND c.relnamespace = n.oid
AND sridcheck.connamespace = n.oid
AND typecheck.connamespace = n.oid
#endif
AND sridcheck.conrelid = c.oid
--AND sridcheck.conname = ''$1''
AND sridcheck.consrc LIKE ''(srid(% = %)''
AND typecheck.conrelid = c.oid
--AND typecheck.conname = ''$2'';
AND typecheck.consrc LIKE
''((geometrytype(%) = ''''%''''::text) OR (% IS NULL))''
;
INSERT INTO geometry_columns SELECT
''''::varchar as f_table_catalogue,
#if USE_VERSION >= 73
n.nspname::varchar as f_table_schema,
#else
''''::varchar as f_table_schema,
#endif
c.relname::varchar as f_table_name,
a.attname::varchar as f_geometry_column,
2 as coord_dimension,
trim(both '' =)'' from substr(sridcheck.consrc,
strpos(sridcheck.consrc, ''='')))::integer as srid,
trim(both '' =)'''''' from substr(typecheck.consrc,
strpos(typecheck.consrc, ''=''),
strpos(typecheck.consrc, ''::'')-
strpos(typecheck.consrc, ''='')
))::varchar as type
#if USE_VERSION < 80
, a.attrelid,
a.attnum as varattnum,
null::histogram2d as stats
#endif
FROM pg_class c, pg_attribute a, pg_type t,
#if USE_VERSION >= 73
pg_namespace n,
#endif
pg_constraint sridcheck, pg_constraint typecheck
WHERE t.typname = ''geometry''
AND a.atttypid = t.oid
AND a.attrelid = c.oid
#if USE_VERSION >= 73
AND c.relnamespace = n.oid
AND sridcheck.connamespace = n.oid
AND typecheck.connamespace = n.oid
#endif
AND sridcheck.conrelid = c.oid
--AND sridcheck.conname = ''$1''
AND sridcheck.consrc LIKE ''(srid(% = %)''
AND typecheck.conrelid = c.oid
--AND typecheck.conname = ''$2''
AND typecheck.consrc LIKE
''((geometrytype(%) = ''''%''''::text) OR (% IS NULL))''
AND NOT EXISTS (
SELECT oid FROM geometry_columns gc
WHERE c.relname::varchar = gc.f_table_name
#if USE_VERSION >= 73
AND n.nspname::varchar = gc.f_table_schema
#endif
AND a.attname::varchar = gc.f_geometry_column
);
GET DIAGNOSTICS inserted = ROW_COUNT;
IF oldcount > probed THEN
stale = oldcount-probed;
ELSE
stale = 0;
END IF;
RETURN ''probed:''||probed||
'' inserted:''||inserted||
'' conflicts:''||probed-inserted||
'' stale:''||stale;
END
' LANGUAGE 'plpgsql';
-----------------------------------------------------------------------
-- ADDGEOMETRYCOLUMN
-- <catalogue>, <schema>, <table>, <column>, <srid>, <type>, <dim>
-----------------------------------------------------------------------
--
-- Type can be one of geometry, GEOMETRYCOLLECTION, POINT, MULTIPOINT, POLYGON,
-- MULTIPOLYGON, LINESTRING, or MULTILINESTRING.
--
-- Types (except geometry) are checked for consistency using a CHECK constraint
-- uses SQL ALTER TABLE command to add the geometry column to the table.
-- Addes a row to geometry_columns.
-- Addes a constraint on the table that all the geometries MUST have the same
-- SRID. Checks the coord_dimension to make sure its between 0 and 3.
-- Should also check the precision grid (future expansion).
-- Calls fix_geometry_columns() at the end.
--
-----------------------------------------------------------------------
CREATEFUNCTION AddGeometryColumn(varchar,varchar,varchar,varchar,integer,varchar,integer)
RETURNS text
AS
'
DECLARE
catalog_name alias for $1;
schema_name alias for $2;
table_name alias for $3;
column_name alias for $4;
new_srid alias for $5;
new_type alias for $6;
new_dim alias for $7;
#if USE_VERSION >= 73
rec RECORD;
schema_ok bool;
real_schema name;
#endif
fixgeomres text;
BEGIN
IF ( not ( (new_type =''GEOMETRY'') or
(new_type =''GEOMETRYCOLLECTION'') or
(new_type =''POINT'') or
(new_type =''MULTIPOINT'') or
(new_type =''POLYGON'') or
(new_type =''MULTIPOLYGON'') or
(new_type =''LINESTRING'') or
(new_type =''MULTILINESTRING'') or
(new_type =''GEOMETRYCOLLECTIONM'') or
(new_type =''POINTM'') or
(new_type =''MULTIPOINTM'') or
(new_type =''POLYGONM'') or
(new_type =''MULTIPOLYGONM'') or
(new_type =''LINESTRINGM'') or
(new_type =''MULTILINESTRINGM'')) )
THEN
RAISE EXCEPTION ''Invalid type name - valid ones are:
GEOMETRY, GEOMETRYCOLLECTION, POINT,
MULTIPOINT, POLYGON, MULTIPOLYGON,
LINESTRING, MULTILINESTRING,
GEOMETRYCOLLECTIONM, POINTM,
MULTIPOINTM, POLYGONM, MULTIPOLYGONM,
LINESTRINGM, or MULTILINESTRINGM '';
return ''fail'';
END IF;
IF ( (new_dim >4) or (new_dim <0) ) THEN
RAISE EXCEPTION ''invalid dimension'';
return ''fail'';
END IF;
IF ( (new_type LIKE ''%M'') and (new_dim!=3) ) THEN
RAISE EXCEPTION ''TypeM needs 3 dimensions'';
return ''fail'';
END IF;
#if USE_VERSION >= 73
IF ( schema_name != '''' ) THEN
schema_ok = ''f'';
FOR rec IN SELECT nspname FROM pg_namespace WHERE text(nspname) = schema_name LOOP
schema_ok := ''t'';
END LOOP;
if ( schema_ok <> ''t'' ) THEN
RAISE NOTICE ''Invalid schema name - using current_schema()'';
SELECT current_schema() into real_schema;
ELSE
real_schema = schema_name;
END IF;
ELSE
SELECT current_schema() into real_schema;
END IF;
#endif
-- Add geometry column
EXECUTE ''ALTER TABLE '' ||
#if USE_VERSION >= 73
quote_ident(real_schema) || ''.'' || quote_ident(table_name)
#else
quote_ident(table_name)
#endif
|| '' ADD COLUMN '' || quote_ident(column_name) ||
'' geometry '';
-- Delete stale record in geometry_column (if any)
EXECUTE ''DELETE FROM geometry_columns WHERE
f_table_catalog = '' || quote_literal('''') ||
'' AND f_table_schema = '' ||
#if USE_VERSION >= 73
quote_literal(real_schema) ||
#else
quote_literal('''') ||
#endif
'' AND f_table_name = '' || quote_literal(table_name) ||
'' AND f_geometry_column = '' || quote_literal(column_name);
-- Add record in geometry_column
EXECUTE ''INSERT INTO geometry_columns VALUES ('' ||
quote_literal('''') || '','' ||
#if USE_VERSION >= 73
quote_literal(real_schema) || '','' ||
#else
quote_literal('''') || '','' ||
#endif
quote_literal(table_name) || '','' ||
quote_literal(column_name) || '','' ||
new_dim || '','' || new_srid || '','' ||
quote_literal(new_type) || '')'';
-- Add table checks
EXECUTE ''ALTER TABLE '' ||
#if USE_VERSION >= 73
quote_ident(real_schema) || ''.'' || quote_ident(table_name)
#else
quote_ident(table_name)
#endif
|| '' ADD CONSTRAINT "enforce_srid_'' ||
column_name || ''" CHECK (SRID('' || quote_ident(column_name) ||
'') = '' || new_srid || '')'' ;
EXECUTE ''ALTER TABLE '' ||
#if USE_VERSION >= 73
quote_ident(real_schema) || ''.'' || quote_ident(table_name)
#else
quote_ident(table_name)
#endif
|| '' ADD CONSTRAINT "enforce_dims_'' ||
column_name || ''" CHECK (ndims('' || quote_ident(column_name) ||
'') = '' || new_dim || '')'' ;
IF (not(new_type = ''GEOMETRY'')) THEN
EXECUTE ''ALTER TABLE '' ||
#if USE_VERSION >= 73
quote_ident(real_schema) || ''.'' || quote_ident(table_name)
#else
quote_ident(table_name)
#endif
|| '' ADD CONSTRAINT "enforce_geotype_'' ||
column_name || ''" CHECK (geometrytype('' ||
quote_ident(column_name) || '')='' ||
quote_literal(new_type) || '' OR ('' ||
quote_ident(column_name) || '') is null)'';
END IF;
SELECT fix_geometry_columns() INTO fixgeomres;
return
#if USE_VERSION >= 73
real_schema || ''.'' ||
#endif
table_name || ''.'' || column_name ||
'' SRID:'' || new_srid ||
'' TYPE:'' || new_type ||
'' DIMS:'' || new_dim || ''\n '' ||
''geometry_column '' || fixgeomres;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
----------------------------------------------------------------------------
-- ADDGEOMETRYCOLUMN ( <schema>, <table>, <column>, <srid>, <type>, <dim> )
----------------------------------------------------------------------------
--
-- This is a wrapper to the real AddGeometryColumn, for use
-- when catalogue is undefined
--
----------------------------------------------------------------------------
CREATEFUNCTION AddGeometryColumn(varchar,varchar,varchar,integer,varchar,integer) RETURNS text AS '
DECLARE
ret text;
BEGIN
SELECT AddGeometryColumn('''',$1,$2,$3,$4,$5,$6) into ret;
RETURN ret;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
----------------------------------------------------------------------------
-- ADDGEOMETRYCOLUMN ( <table>, <column>, <srid>, <type>, <dim> )
----------------------------------------------------------------------------
--
-- This is a wrapper to the real AddGeometryColumn, for use
-- when catalogue and schema are undefined
--
----------------------------------------------------------------------------
CREATEFUNCTION AddGeometryColumn(varchar,varchar,integer,varchar,integer) RETURNS text AS '
DECLARE
ret text;
BEGIN
SELECT AddGeometryColumn('''','''',$1,$2,$3,$4,$5) into ret;
RETURN ret;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- DROPGEOMETRYCOLUMN
-- <catalogue>, <schema>, <table>, <column>
-----------------------------------------------------------------------
--
-- Removes geometry column reference from geometry_columns table.
-- Drops the column with pgsql >= 73.
-- Make some silly enforcements on it for pgsql < 73
--
-----------------------------------------------------------------------
CREATEFUNCTION DropGeometryColumn(varchar, varchar,varchar,varchar)
RETURNS text
AS
'
DECLARE
catalog_name alias for $1;
schema_name alias for $2;
table_name alias for $3;
column_name alias for $4;
myrec RECORD;
okay boolean;
real_schema name;
BEGIN
#if USE_VERSION >= 73
-- Find, check or fix schema_name
IF ( schema_name != '''' ) THEN
okay = ''f'';
FOR myrec IN SELECT nspname FROM pg_namespace WHERE text(nspname) = schema_name LOOP
okay := ''t'';
END LOOP;
IF ( okay <> ''t'' ) THEN
RAISE NOTICE ''Invalid schema name - using current_schema()'';
SELECT current_schema() into real_schema;
ELSE
real_schema = schema_name;
END IF;
ELSE
SELECT current_schema() into real_schema;
END IF;
#else
real_schema = schema_name;
#endif // USE_VERSION >= 73
-- Find out if the column is in the geometry_columns table
okay = ''f'';
FOR myrec IN SELECT * from geometry_columns where f_table_schema = text(real_schema) and f_table_name = table_name and f_geometry_column = column_name LOOP
okay := ''t'';
END LOOP;
IF (okay <> ''t'') THEN
RAISE EXCEPTION ''column not found in geometry_columns table'';
RETURN ''f'';
END IF;
-- Remove ref from geometry_columns table
EXECUTE ''delete from geometry_columns where f_table_schema = '' ||
quote_literal(real_schema) || '' and f_table_name = '' ||
quote_literal(table_name) || '' and f_geometry_column = '' ||
quote_literal(column_name);
#if USE_VERSION < 73
-- Remove not-null constraint to table column
EXECUTE ''update pg_attribute set attnotnull = false from pg_class where pg_attribute.attrelid = pg_class.oid and pg_class.relname = '' || quote_literal(table_name) ||'' and pg_attribute.attname = '' || quote_literal(column_name);
-- update the given table/column so that it it all NULLS
EXECUTE ''update "''||table_name||''" set "''||column_name||''"= NULL'';
-- add = NULL constraint to given table/column
EXECUTE ''ALTER TABLE "''||table_name||''" ADD CHECK ("''||column_name||''" IS NULL)'';
#else
-- Remove table column
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) || ''.'' ||
quote_ident(table_name) || '' DROP COLUMN '' ||
quote_ident(column_name);
#endif
RETURN real_schema || ''.'' || table_name || ''.'' || column_name ||'' effectively removed.'';
END;
'
LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- DROPGEOMETRYCOLUMN
-- <schema>, <table>, <column>
-----------------------------------------------------------------------
--
-- This is a wrapper to the real DropGeometryColumn, for use
-- when catalogue is undefined
--
-----------------------------------------------------------------------
CREATEFUNCTION DropGeometryColumn(varchar,varchar,varchar)
RETURNS text
AS
'
DECLARE
ret text;
BEGIN
SELECT DropGeometryColumn('''',$1,$2,$3) into ret;
RETURN ret;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- DROPGEOMETRYCOLUMN
-- <table>, <column>
-----------------------------------------------------------------------
--
-- This is a wrapper to the real DropGeometryColumn, for use
-- when catalogue and schema is undefined.
--
-----------------------------------------------------------------------
CREATEFUNCTION DropGeometryColumn(varchar,varchar)
RETURNS text
AS
'
DECLARE
ret text;
BEGIN
SELECT DropGeometryColumn('''','''',$1,$2) into ret;
RETURN ret;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- DROPGEOMETRYTABLE
-- <catalogue>, <schema>, <table>
-----------------------------------------------------------------------
--
-- Drop a table and all its references in geometry_columns
--
-----------------------------------------------------------------------
CREATEFUNCTION DropGeometryTable(varchar, varchar,varchar)
RETURNS text
AS
'
DECLARE
catalog_name alias for $1;
schema_name alias for $2;
table_name alias for $3;
real_schema name;
BEGIN
#if USE_VERSION >= 73
IF ( schema_name = '''' ) THEN
SELECT current_schema() into real_schema;
ELSE
real_schema = schema_name;
END IF;
#endif // USE_VERSION >= 73
-- Remove refs from geometry_columns table
EXECUTE ''DELETE FROM geometry_columns WHERE '' ||
#if USE_VERSION >= 73
''f_table_schema = '' || quote_literal(real_schema) ||
'' AND '' ||
#endif
'' f_table_name = '' || quote_literal(table_name);
-- Remove table
EXECUTE ''DROP TABLE ''
#if USE_VERSION >= 73
|| quote_ident(real_schema) || ''.'' ||
#endif
quote_ident(table_name);
RETURN
#if USE_VERSION >= 73
real_schema || ''.'' ||
#endif
table_name ||'' dropped.'';
END;
'
LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- DROPGEOMETRYTABLE
-- <schema>, <table>
-----------------------------------------------------------------------
--
-- Drop a table and all its references in geometry_columns
--
-----------------------------------------------------------------------
CREATEFUNCTION DropGeometryTable(varchar,varchar) RETURNS text AS
'SELECT DropGeometryTable('''',$1,$2)'
LANGUAGE 'sql' WITH (isstrict);
-----------------------------------------------------------------------
-- DROPGEOMETRYTABLE
-- <table>
-----------------------------------------------------------------------
--
-- Drop a table and all its references in geometry_columns
-- For PG>=73 use current_schema()
--
-----------------------------------------------------------------------
CREATEFUNCTION DropGeometryTable(varchar) RETURNS text AS
'SELECT DropGeometryTable('''','''',$1)'
LANGUAGE 'sql' WITH (isstrict);
-----------------------------------------------------------------------
-- UPDATEGEOMETRYSRID
-- <catalogue>, <schema>, <table>, <column>, <srid>
-----------------------------------------------------------------------
--
-- Change SRID of all features in a spatially-enabled table
--
-----------------------------------------------------------------------
CREATEFUNCTION UpdateGeometrySRID(varchar,varchar,varchar,varchar,integer)
RETURNS text
AS
'
DECLARE
catalog_name alias for $1;
schema_name alias for $2;
table_name alias for $3;
column_name alias for $4;
new_srid alias for $5;
myrec RECORD;
okay boolean;
cname varchar;
real_schema name;
BEGIN
#if USE_VERSION >= 73
-- Find, check or fix schema_name
IF ( schema_name != '''' ) THEN
okay = ''f'';
FOR myrec IN SELECT nspname FROM pg_namespace WHERE text(nspname) = schema_name LOOP
okay := ''t'';
END LOOP;
IF ( okay <> ''t'' ) THEN
RAISE EXCEPTION ''Invalid schema name'';
ELSE
real_schema = schema_name;
END IF;
ELSE
SELECT INTO real_schema current_schema()::text;
END IF;
#endif // USE_VERSION >= 73
-- Find out if the column is in the geometry_columns table
okay = ''f'';
FOR myrec IN SELECT * from geometry_columns where f_table_schema = text(real_schema) and f_table_name = table_name and f_geometry_column = column_name LOOP
okay := ''t'';
END LOOP;
IF (okay <> ''t'') THEN
RAISE EXCEPTION ''column not found in geometry_columns table'';
RETURN ''f'';
END IF;
-- Update ref from geometry_columns table
EXECUTE ''UPDATE geometry_columns SET SRID = '' || new_srid ||
'' where f_table_schema = '' ||
quote_literal(real_schema) || '' and f_table_name = '' ||
quote_literal(table_name) || '' and f_geometry_column = '' ||
quote_literal(column_name);
-- Make up constraint name
cname = ''enforce_srid_'' || column_name;
-- Drop enforce_srid constraint
#if USE_VERSION < 73
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
#else
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
''.'' || quote_ident(table_name) ||
#endif
'' DROP constraint '' || quote_ident(cname);
-- Update geometries SRID
#if USE_VERSION < 73
EXECUTE ''UPDATE '' || quote_ident(table_name) ||
#else
EXECUTE ''UPDATE '' || quote_ident(real_schema) ||
''.'' || quote_ident(table_name) ||
#endif
'' SET '' || quote_ident(column_name) ||
'' = setSRID('' || quote_ident(column_name) ||
'', '' || new_srid || '')'';
-- Reset enforce_srid constraint
#if USE_VERSION < 73
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
#else
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
''.'' || quote_ident(table_name) ||
#endif
'' ADD constraint '' || quote_ident(cname) ||
'' CHECK (srid('' || quote_ident(column_name) ||
'') = '' || new_srid || '')'';
RETURN real_schema || ''.'' || table_name || ''.'' || column_name ||'' SRID changed to '' || new_srid;
END;
'
LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- UPDATEGEOMETRYSRID
-- <schema>, <table>, <column>, <srid>
-----------------------------------------------------------------------
CREATEFUNCTION UpdateGeometrySRID(varchar,varchar,varchar,integer)
RETURNS text
AS '
DECLARE
ret text;
BEGIN
SELECT UpdateGeometrySRID('''',$1,$2,$3,$4) into ret;
RETURN ret;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- UPDATEGEOMETRYSRID
-- <table>, <column>, <srid>
-----------------------------------------------------------------------
CREATEFUNCTION UpdateGeometrySRID(varchar,varchar,integer)
RETURNS text
AS '
DECLARE
ret text;
BEGIN
SELECT UpdateGeometrySRID('''','''',$1,$2,$3) into ret;
RETURN ret;
END;
' LANGUAGE 'plpgsql' WITH (isstrict);
-----------------------------------------------------------------------
-- UPDATE_GEOMETRY_STATS()
-----------------------------------------------------------------------
--
-- Only meaningful for PG<75.
-- Gather statisticts about geometry columns for use
-- with cost estimator.
--
-- It is defined also for PG>=75 for back-compatibility
--
-----------------------------------------------------------------------
#if USE_VERSION >= 80
CREATEFUNCTION update_geometry_stats() RETURNS text
AS ' SELECT ''update_geometry_stats() has been obsoleted. Statistics are automatically built running the ANALYZE command''::text' LANGUAGE 'sql';
#else // USE_VERSION < 80
CREATEFUNCTION update_geometry_stats()
RETURNS text
AS
'
DECLARE
result text;
stated integer;
fixres text;
BEGIN
SELECT fix_geometry_columns() INTO fixres;
UPDATE geometry_columns SET
#if USE_VERSION >= 73
stats = (build_histogram2d(create_histogram2d(
find_extent(f_table_schema, f_table_name, f_geometry_column), 40), f_table_schema, f_table_name, f_geometry_column))
FROM pg_class c, pg_attribute a, pg_namespace n
WHERE n.nspname = f_table_schema::name
AND c.relname = f_table_name::name
AND c.relnamespace = n.oid
#else // USE_VERSION < 73
stats = (build_histogram2d(create_histogram2d(
find_extent(f_table_name, f_geometry_column),
40), f_table_name, f_geometry_column))
FROM pg_class c, pg_attribute a
WHERE c.relname = f_table_name::name
#endif
AND a.attname = f_geometry_column::name
AND a.attrelid = c.oid
AND geometry_columns.attrelid is not null;
GET DIAGNOSTICS stated = ROW_COUNT;
result = fixres || '' stats:'' || stated::text;
return result;
END;
'
LANGUAGE 'plpgsql' ;
#endif // USE_VERSION < 80
-----------------------------------------------------------------------
-- UPDATE_GEOMETRY_STATS( <table>, <column> )
-----------------------------------------------------------------------
--
-- Only meaningful for PG<75.
-- Gather statisticts about a geometry column for use
-- with cost estimator.
--
-- It is defined also for PG>=75 for back-compatibility
--
-----------------------------------------------------------------------
#if USE_VERSION >= 80
CREATEFUNCTION update_geometry_stats(varchar,varchar) RETURNS text
AS 'SELECT update_geometry_stats();' LANGUAGE 'sql' ;
#else
CREATEFUNCTION update_geometry_stats(varchar,varchar) RETURNS text
AS
'
DECLARE
tablename aliAS for $1;
columnname aliAS for $2;
stated integer;
result text;
fixres text;
BEGIN
SELECT fix_geometry_columns() INTO fixres;
EXECUTE ''UPDATE geometry_columns SET
#if USE_VERSION >= 73
stats = (build_histogram2d(create_histogram2d(
find_extent(f_table_schema,
f_table_name,
f_geometry_column), 40),
f_table_schema, f_table_name,
f_geometry_column))
FROM pg_class c, pg_attribute a, pg_namespace n
WHERE n.nspname = f_table_schema::name
AND c.relname = f_table_name::name
AND a.attname = f_geometry_column::name
AND c.relnamespace = n.oid
AND a.attrelid = c.oid
#else // USE_VERSION < 73
stats = (build_histogram2d(create_histogram2d(
find_extent(f_table_name, f_geometry_column),
40), f_table_name, f_geometry_column))
FROM pg_class c, pg_attribute a
WHERE c.relname = f_table_name::name
AND a.attname = f_geometry_column::name
AND a.attrelid = c.oid
#endif
AND f_table_name = '' || quote_literal(tablename) || ''
AND f_geometry_column = '' || quote_literal(columnname)
|| '' AND geometry_columns.attrelid is not null'';
GET DIAGNOSTICS stated = ROW_COUNT;
result = fixres || '' stats:'' || stated::text;
return result;
END;
'
LANGUAGE 'plpgsql' ;
#endif // USE_VERSION < 80
-----------------------------------------------------------------------
-- FIND_SRID( <schema>, <table>, <geom col> )
-----------------------------------------------------------------------
CREATEFUNCTION find_srid(varchar,varchar,varchar) RETURNS int4 AS
'DECLARE
schem text;
tabl text;
sr int4;
BEGIN
IF $1 IS NULL THEN
RAISE EXCEPTION ''find_srid() - schema is NULL!'';
END IF;
IF $2 IS NULL THEN
RAISE EXCEPTION ''find_srid() - table name is NULL!'';
END IF;
IF $3 IS NULL THEN
RAISE EXCEPTION ''find_srid() - column name is NULL!'';
END IF;
schem = $1;
tabl = $2;
-- if the table contains a . and the schema is empty
-- split the table into a schema and a table
-- otherwise drop through to default behavior
IF ( schem = '''' and tabl LIKE ''%.%'' ) THEN
schem = substr(tabl,1,strpos(tabl,''.'')-1);
tabl = substr(tabl,length(schem)+2);
ELSE
schem = schem || ''%'';
END IF;
select SRID into sr from geometry_columns where f_table_schema like schem and f_table_name = tabl and f_geometry_column = $3;
IF NOT FOUND THEN
RAISE EXCEPTION ''find_srid() - couldnt find the corresponding SRID - is the geometry registered in the GEOMETRY_COLUMNS table? Is there an uppercase/lowercase missmatch?'';
END IF;
return sr;
END;
'
LANGUAGE 'plpgsql' WITH (iscachable);
---------------------------------------------------------------
-- PROJ support
---------------------------------------------------------------
CREATEFUNCTION get_proj4_from_srid(integer) RETURNS text AS
'SELECT proj4text::text FROM spatial_ref_sys WHERE srid= $1'
LANGUAGE 'sql' WITH (iscachable,isstrict);
CREATEFUNCTION transform_geometry(geometry,text,text,int)
RETURNS geometry
AS '@MODULE_FILENAME@','transform_geom'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION transform(geometry,integer) RETURNS geometry AS
'BEGIN
RETURN transform_geometry( $1 , get_proj4_from_srid(SRID( $1 ) ), get_proj4_from_srid( $2 ), $2 );
END;'
LANGUAGE 'plpgsql' WITH (iscachable,isstrict);
-----------------------------------------------------------------------
-- POSTGIS_VERSION()
-----------------------------------------------------------------------
CREATEFUNCTION postgis_version() RETURNS text
AS 'SELECT \'@POSTGIS_VERSION@\'::text AS version'
LANGUAGE 'sql';
CREATEFUNCTION postgis_proj_version() RETURNS text
AS '@MODULE_FILENAME@'
LANGUAGE 'C';
CREATEFUNCTION postgis_scripts_installed() RETURNS text
AS 'SELECT \'@POSTGIS_SCRIPTS_VERSION@\'::text AS version'
LANGUAGE 'sql';
CREATEFUNCTION postgis_lib_version() RETURNS text
AS '@MODULE_FILENAME@'
LANGUAGE 'C';
CREATEFUNCTION postgis_scripts_released() RETURNS text
AS '@MODULE_FILENAME@'
LANGUAGE 'C';
CREATEFUNCTION postgis_uses_stats() RETURNS bool
AS '@MODULE_FILENAME@'
LANGUAGE 'C';
CREATEFUNCTION postgis_geos_version() RETURNS text
AS '@MODULE_FILENAME@'
LANGUAGE 'C';
CREATEFUNCTION postgis_full_version() RETURNS text
AS '
DECLARE
libver text;
projver text;
geosver text;
usestats bool;
dbproc text;
relproc text;
fullver text;
BEGIN
SELECT postgis_lib_version() INTO libver;
SELECT postgis_proj_version() INTO projver;
SELECT postgis_geos_version() INTO geosver;
SELECT postgis_uses_stats() INTO usestats;
SELECT postgis_scripts_installed() INTO dbproc;
SELECT postgis_scripts_released() INTO relproc;
fullver = \'POSTGIS="\' || libver || \'"\';
IF geosver IS NOT NULL THEN
fullver = fullver || \' GEOS="\' || geosver || \'"\';
END IF;
IF projver IS NOT NULL THEN
fullver = fullver || \' PROJ="\' || projver || \'"\';
END IF;
IF usestats THEN
fullver = fullver || \' USE_STATS\';
END IF;
fullver = fullver || \' DBPROC="\' || dbproc || \'"\';
fullver = fullver || \' RELPROC="\' || relproc || \'"\';
IF dbproc != relproc THEN
fullver = fullver || \' (needs proc upgrade)\';
END IF;
RETURN fullver;
END
' LANGUAGE 'plpgsql';
---------------------------------------------------------------
-- CASTS
---------------------------------------------------------------
CREATEFUNCTION box2d(geometry)
RETURNS box2d
AS '@MODULE_FILENAME@','LWGEOM_to_BOX2DFLOAT4'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box3d(geometry)
RETURNS box3d
AS '@MODULE_FILENAME@','LWGEOM_to_BOX3D'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box(geometry)
RETURNS box
AS '@MODULE_FILENAME@','LWGEOM_to_BOX'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box2d(box3d)
RETURNS box2d
AS '@MODULE_FILENAME@','BOX3D_to_BOX2DFLOAT4'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box3d(box2d)
RETURNS box3d
AS '@MODULE_FILENAME@','BOX2DFLOAT4_to_BOX3D'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION box(box3d)
RETURNS box
AS '@MODULE_FILENAME@','BOX3D_to_BOX'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION text(geometry)
RETURNS text
AS '@MODULE_FILENAME@','LWGEOM_to_text'
LANGUAGE 'C' WITH (isstrict,iscachable);
-- this is kept for backward-compatibility
CREATEFUNCTION box3dtobox(box3d)
RETURNS box
AS 'SELECT box($1)'
LANGUAGE 'SQL' WITH (isstrict,iscachable);
CREATEFUNCTION geometry(box2d)
RETURNS geometry
AS '@MODULE_FILENAME@','BOX2DFLOAT4_to_LWGEOM'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry(box3d)
RETURNS geometry
AS '@MODULE_FILENAME@','BOX3D_to_LWGEOM'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry(text)
RETURNS geometry
AS '@MODULE_FILENAME@','parse_WKT_lwgeom'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry(chip)
RETURNS geometry
AS '@MODULE_FILENAME@','CHIP_to_LWGEOM'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION geometry(bytea)
RETURNS geometry
AS '@MODULE_FILENAME@','LWGEOMFromWKB'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION bytea(geometry)
RETURNS bytea
AS '@MODULE_FILENAME@','WKBFromLWGEOM'
LANGUAGE 'C' WITH (isstrict,iscachable);
-- 7.3+ explicit casting definitions
#if USE_VERSION >= 73
CREATE CAST (geometry AS box2d) WITH FUNCTION box2d(geometry) AS IMPLICIT;
CREATE CAST (geometry AS box3d) WITH FUNCTION box3d(geometry) AS IMPLICIT;
CREATE CAST (geometry AS box) WITH FUNCTION box(geometry) AS IMPLICIT;
CREATE CAST (box3d AS box2d) WITH FUNCTION box2d(box3d) AS IMPLICIT;
CREATE CAST (box2d AS box3d) WITH FUNCTION box3d(box2d) AS IMPLICIT;
CREATE CAST (box2d AS geometry) WITH FUNCTION geometry(box2d) AS IMPLICIT;
CREATE CAST (box3d AS box) WITH FUNCTION box(box3d) AS IMPLICIT;
CREATE CAST (box3d AS geometry) WITH FUNCTION geometry(box3d) AS IMPLICIT;
CREATE CAST (text AS geometry) WITH FUNCTION geometry(text) AS IMPLICIT;
CREATE CAST (geometry AS text) WITH FUNCTION text(geometry) AS IMPLICIT;
CREATE CAST (chip AS geometry) WITH FUNCTION geometry(chip) AS IMPLICIT;
CREATE CAST (bytea AS geometry) WITH FUNCTION geometry(bytea) AS IMPLICIT;
CREATE CAST (geometry AS bytea) WITH FUNCTION bytea(geometry) AS IMPLICIT;
#endif
---------------------------------------------------------------
-- Algorithms
---------------------------------------------------------------
CREATEFUNCTION simplify(geometry, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_simplify2d'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION line_interpolate_point(geometry, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_line_interpolate_point'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION segmentize(geometry, float8)
RETURNS geometry
AS '@MODULE_FILENAME@', 'LWGEOM_segmentize2d'
LANGUAGE 'C' WITH (isstrict,iscachable);
---------------------------------------------------------------
-- GEOS
---------------------------------------------------------------
CREATEFUNCTION intersection(geometry,geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','intersection'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION buffer(geometry,float8)
RETURNS geometry
AS '@MODULE_FILENAME@','buffer'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION buffer(geometry,float8,integer)
RETURNS geometry
AS '@MODULE_FILENAME@','buffer'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION convexhull(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','convexhull'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION difference(geometry,geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','difference'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION boundary(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','boundary'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION symdifference(geometry,geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','symdifference'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION symmetricdifference(geometry,geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','symdifference'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION GeomUnion(geometry,geometry)
RETURNS geometry
AS '@MODULE_FILENAME@','geomunion'
LANGUAGE 'C' WITH (isstrict);
CREATE AGGREGATE MemGeomUnion (
basetype = geometry,
sfunc = geomunion,
stype = geometry
);
CREATEFUNCTION unite_garray (geometry[])
RETURNS geometry
AS '@MODULE_FILENAME@'
LANGUAGE 'C';
CREATE AGGREGATE GeomUnion (
sfunc = geom_accum,
basetype = geometry,
stype = geometry[],
finalfunc = unite_garray
);
CREATEFUNCTION relate(geometry,geometry)
RETURNS text
AS '@MODULE_FILENAME@','relate_full'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION relate(geometry,geometry,text)
RETURNS boolean
AS '@MODULE_FILENAME@','relate_pattern'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION disjoint(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION touches(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION intersects(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION crosses(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION within(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION contains(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION overlaps(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION IsValid(geometry)
RETURNS boolean
AS '@MODULE_FILENAME@', 'isvalid'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION GEOSnoop(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@', 'GEOSnoop'
LANGUAGE 'C' WITH (isstrict);
-- This is also available w/out GEOS
CREATEFUNCTION Centroid(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION IsRing(geometry)
RETURNS boolean
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION PointOnSurface(geometry)
RETURNS geometry
AS '@MODULE_FILENAME@'
LANGUAGE 'C' WITH (isstrict);
CREATEFUNCTION IsSimple(geometry)
RETURNS boolean
AS '@MODULE_FILENAME@', 'issimple'
LANGUAGE 'C' WITH (isstrict,iscachable);
CREATEFUNCTION Equals(geometry,geometry)
RETURNS boolean
AS '@MODULE_FILENAME@','geomequals'
LANGUAGE 'C' WITH (isstrict,iscachable);
-----------------------------------------------------------------------
-- SVG OUTPUT
-----------------------------------------------------------------------
CREATEFUNCTION AsSVG(geometry,int4,int4)
RETURNS TEXT
AS '@MODULE_FILENAME@','assvg_geometry'
LANGUAGE 'C';
CREATEFUNCTION AsSVG(geometry,int4)
RETURNS TEXT
AS '@MODULE_FILENAME@','assvg_geometry'
LANGUAGE 'C';
CREATEFUNCTION AsSVG(geometry)
RETURNS TEXT
AS '@MODULE_FILENAME@','assvg_geometry'
LANGUAGE 'C';
-----------------------------------------------------------------------
-- GML OUTPUT
-----------------------------------------------------------------------
-- AsGML(geom, precision, version)
CREATEFUNCTION AsGML(geometry, int4, int4)
RETURNS TEXT
AS '@MODULE_FILENAME@','LWGEOM_asGML'
LANGUAGE 'C';
-- AsGML(geom, precision) / version=2
CREATEFUNCTION AsGML(geometry, int4)
RETURNS TEXT
AS '@MODULE_FILENAME@','LWGEOM_asGML'
LANGUAGE 'C';
-- AsGML(geom) / precision=15 version=2
CREATEFUNCTION AsGML(geometry)
RETURNS TEXT
AS '@MODULE_FILENAME@','LWGEOM_asGML'
LANGUAGE 'C';
---------------------------------------------------------------
-- END
---------------------------------------------------------------
COMMIT;
Reference in a new issue View git blame Copy permalink