mirror of
https://git.osgeo.org/gitea/postgis/postgis
synced 2024-10-26 01:52:37 +00:00
768a8d838e
git-svn-id: http://svn.osgeo.org/postgis/trunk@477 b70326c6-7e19-0410-871a-916f4a2858ee
3959 lines
108 KiB
MySQL
3959 lines
108 KiB
MySQL
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
--
|
|
-- $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 hte GNU General Public Licence. See the COPYING file.
|
|
--
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
|
|
BEGIN TRANSACTION;
|
|
|
|
-- You might have to define the PL/PgSQL language usually done with the
|
|
-- changelang script.
|
|
|
|
-- Here's some hokey code to test to see if PL/PgSQL is installed
|
|
-- if it is, you get a message "PL/PgSQL is installed"
|
|
-- otherwise it will give a big error message.
|
|
|
|
(select 'PL/PgSQL is installed.' as message from pg_language where lanname='plpgsql') union (select 'You must install PL/PgSQL before running this SQL file,\nor you will get an error. To install PL/PgSQL run:\n\tcreatelang plpgsql <dbname>'::text as message) order by message limit 1;
|
|
|
|
|
|
-----------------------------------------
|
|
-- THESE WERE postgis_sql_VER_start.sql.in
|
|
-----------------------------------------
|
|
|
|
#if USE_VERSION == 71
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- HISTOGRAM2D
|
|
|
|
CREATE FUNCTION histogram2d_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION histogram2d_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE histogram2d (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = histogram2d_in,
|
|
output = histogram2d_out,
|
|
storage = main
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BOX3D
|
|
|
|
CREATE FUNCTION box3d_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION box3d_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE box3d (
|
|
alignment = double,
|
|
internallength = 48,
|
|
input = box3d_in,
|
|
output = box3d_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- SPHEROID
|
|
|
|
CREATE FUNCTION spheroid_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'ellipsoid_in'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION spheroid_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'ellipsoid_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE spheroid (
|
|
alignment = double,
|
|
internallength = 65,
|
|
input = spheroid_in,
|
|
output = spheroid_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
CREATE FUNCTION wkb_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'WKB_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'WKB_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE wkb (
|
|
internallength = variable,
|
|
input = wkb_in,
|
|
output = wkb_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CHIP
|
|
|
|
CREATE FUNCTION chip_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'CHIP_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION chip_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'CHIP_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE chip (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = chip_in,
|
|
output = chip_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- GEOMETRY
|
|
|
|
CREATE FUNCTION geometry_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE geometry (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = geometry_in,
|
|
output = geometry_out,
|
|
storage = main
|
|
);
|
|
|
|
--
|
|
-- GiST Selectivity Function
|
|
--
|
|
|
|
CREATE FUNCTION postgis_gist_sel(oid, oid, int2, opaque, int4)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
#elif USE_VERSION == 72
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- HISTOGRAM2D
|
|
|
|
CREATE FUNCTION histogram2d_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION histogram2d_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE histogram2d (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = histogram2d_in,
|
|
output = histogram2d_out,
|
|
storage = main
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BOX3D
|
|
|
|
CREATE FUNCTION box3d_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION box3d_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE box3d (
|
|
alignment = double,
|
|
internallength = 48,
|
|
input = box3d_in,
|
|
output = box3d_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- SPHEROID
|
|
|
|
CREATE FUNCTION spheroid_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'ellipsoid_in'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION spheroid_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'ellipsoid_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE spheroid (
|
|
alignment = double,
|
|
internallength = 65,
|
|
input = spheroid_in,
|
|
output = spheroid_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
CREATE FUNCTION wkb_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'WKB_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'WKB_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE wkb (
|
|
internallength = variable,
|
|
input = wkb_in,
|
|
output = wkb_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CHIP
|
|
|
|
CREATE FUNCTION chip_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'CHIP_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION chip_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'CHIP_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE chip (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = chip_in,
|
|
output = chip_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- GEOMETRY
|
|
|
|
CREATE FUNCTION geometry_in(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_out(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE geometry (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = geometry_in,
|
|
output = geometry_out,
|
|
storage = main
|
|
);
|
|
|
|
--
|
|
-- GiST Selectivity Function
|
|
--
|
|
|
|
CREATE FUNCTION postgis_gist_sel(opaque, oid, opaque, int4)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
#elif USE_VERSION == 73
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- HISTOGRAM2D
|
|
|
|
CREATE FUNCTION histogram2d_in(cstring)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION histogram2d_out(histogram2d)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE histogram2d (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = histogram2d_in,
|
|
output = histogram2d_out,
|
|
storage = main
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BOX3D
|
|
|
|
CREATE FUNCTION box3d_in(cstring)
|
|
RETURNS box3d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION box3d_out(box3d)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE box3d (
|
|
alignment = double,
|
|
internallength = 48,
|
|
input = box3d_in,
|
|
output = box3d_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- SPHEROID
|
|
|
|
CREATE FUNCTION spheroid_in(cstring)
|
|
RETURNS spheroid
|
|
AS MODULE_FILENAME,'ellipsoid_in'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION spheroid_out(spheroid)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'ellipsoid_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE spheroid (
|
|
alignment = double,
|
|
internallength = 65,
|
|
input = spheroid_in,
|
|
output = spheroid_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
CREATE FUNCTION wkb_in(cstring)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'WKB_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_out(wkb)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'WKB_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE wkb (
|
|
internallength = variable,
|
|
input = wkb_in,
|
|
output = wkb_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CHIP
|
|
|
|
CREATE FUNCTION chip_in(cstring)
|
|
RETURNS chip
|
|
AS MODULE_FILENAME,'CHIP_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION chip_out(chip)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'CHIP_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE chip (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = chip_in,
|
|
output = chip_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- GEOMETRY
|
|
|
|
CREATE FUNCTION geometry_in(cstring)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_out(geometry)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE geometry (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = geometry_in,
|
|
output = geometry_out,
|
|
storage = main
|
|
);
|
|
|
|
--
|
|
-- GiST selectivity function
|
|
--
|
|
|
|
CREATE FUNCTION postgis_gist_sel (internal, oid, internal, int4)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
#elif USE_VERSION == 74
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- HISTOGRAM2D
|
|
|
|
CREATE FUNCTION histogram2d_in(cstring)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION histogram2d_out(histogram2d)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE histogram2d (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = histogram2d_in,
|
|
output = histogram2d_out,
|
|
storage = main
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BOX3D
|
|
|
|
CREATE FUNCTION box3d_in(cstring)
|
|
RETURNS box3d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION box3d_out(box3d)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE box3d (
|
|
alignment = double,
|
|
internallength = 48,
|
|
input = box3d_in,
|
|
output = box3d_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- SPHEROID
|
|
|
|
CREATE FUNCTION spheroid_in(cstring)
|
|
RETURNS spheroid
|
|
AS MODULE_FILENAME,'ellipsoid_in'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION spheroid_out(spheroid)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'ellipsoid_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE spheroid (
|
|
alignment = double,
|
|
internallength = 65,
|
|
input = spheroid_in,
|
|
output = spheroid_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
CREATE FUNCTION wkb_in(cstring)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'WKB_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_out(wkb)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'WKB_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_recv(internal)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'WKB_recv'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CHIP
|
|
|
|
CREATE FUNCTION chip_in(cstring)
|
|
RETURNS chip
|
|
AS MODULE_FILENAME,'CHIP_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION chip_out(chip)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'CHIP_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE chip (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = chip_in,
|
|
output = chip_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- GEOMETRY
|
|
|
|
CREATE FUNCTION geometry_in(cstring)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_out(geometry)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE geometry (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = geometry_in,
|
|
output = geometry_out,
|
|
storage = main
|
|
);
|
|
|
|
--
|
|
-- GiST selectivity function
|
|
--
|
|
|
|
CREATE FUNCTION postgis_gist_sel (internal, oid, internal, int4)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
|
|
#elif USE_VERSION == 75
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- HISTOGRAM2D
|
|
|
|
CREATE FUNCTION histogram2d_in(cstring)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION histogram2d_out(histogram2d)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE histogram2d (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = histogram2d_in,
|
|
output = histogram2d_out,
|
|
storage = main
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BOX3D
|
|
|
|
CREATE FUNCTION box3d_in(cstring)
|
|
RETURNS box3d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION box3d_out(box3d)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE box3d (
|
|
alignment = double,
|
|
internallength = 48,
|
|
input = box3d_in,
|
|
output = box3d_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- SPHEROID
|
|
|
|
CREATE FUNCTION spheroid_in(cstring)
|
|
RETURNS spheroid
|
|
AS MODULE_FILENAME,'ellipsoid_in'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION spheroid_out(spheroid)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'ellipsoid_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE spheroid (
|
|
alignment = double,
|
|
internallength = 65,
|
|
input = spheroid_in,
|
|
output = spheroid_out
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
CREATE FUNCTION wkb_in(cstring)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'WKB_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_out(wkb)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'WKB_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION wkb_recv(internal)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'WKB_recv'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CHIP
|
|
|
|
CREATE FUNCTION chip_in(cstring)
|
|
RETURNS chip
|
|
AS MODULE_FILENAME,'CHIP_in'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION chip_out(chip)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME,'CHIP_out'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE chip (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = chip_in,
|
|
output = chip_out,
|
|
storage = extended
|
|
);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- GEOMETRY
|
|
|
|
CREATE FUNCTION geometry_in(cstring)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_out(geometry)
|
|
RETURNS cstring
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_analyze(internal)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE TYPE geometry (
|
|
alignment = double,
|
|
internallength = variable,
|
|
input = geometry_in,
|
|
output = geometry_out,
|
|
storage = main,
|
|
analyze = geometry_analyze
|
|
);
|
|
|
|
--
|
|
-- GiST selectivity function
|
|
--
|
|
|
|
CREATE FUNCTION postgis_gist_sel (internal, oid, internal, int4)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
|
|
|
|
#endif // USE_VERSION == ...
|
|
|
|
---------------------------------------
|
|
-- THIS WAS postgis_sql_common.sql.in
|
|
---------------------------------------
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- 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 < 75
|
|
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 ) );
|
|
|
|
--
|
|
-- Workaround for old user defined variable length datatype
|
|
-- default value bug. Should not be necessary > 7.2
|
|
--
|
|
|
|
UPDATE pg_type SET typdefault = NULL WHERE typname = 'wkb';
|
|
UPDATE pg_type SET typdefault = NULL WHERE typname = 'geometry';
|
|
UPDATE pg_type SET typdefault = NULL WHERE typname = 'histogram2d';
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- POSTGIS_VERSION()
|
|
|
|
CREATE FUNCTION postgis_version() RETURNS text
|
|
AS 'SELECT 'POSTGIS_VERSION'::text AS version'
|
|
LANGUAGE 'sql';
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIND_SRID( <schema>, <table>, <geom col> )
|
|
|
|
CREATE FUNCTION 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);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- GET_PROJ4_FROM_SRID( <srid> )
|
|
|
|
CREATE FUNCTION get_proj4_from_srid(integer) RETURNS text AS
|
|
'SELECT proj4text::text FROM spatial_ref_sys WHERE srid= $1'
|
|
LANGUAGE 'sql' WITH (iscachable,isstrict);
|
|
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIND_EXTENT( <table name>, <column name> )
|
|
|
|
CREATE FUNCTION find_extent(text,text) RETURNS box3d AS
|
|
'
|
|
DECLARE
|
|
tablename alias for $1;
|
|
columnname alias for $2;
|
|
okay boolean;
|
|
myrec RECORD;
|
|
|
|
BEGIN
|
|
FOR myrec IN EXECUTE ''SELECT extent("''||columnname||''") FROM "''||tablename||''"'' LOOP
|
|
return myrec.extent;
|
|
END LOOP;
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- TRANSFORM ( <geometry>, <srid> )
|
|
--
|
|
-- Test:
|
|
--
|
|
-- trans=# select * from spatial_ref_sys ;
|
|
--
|
|
-- srid | auth_name | auth_srid | srtext | proj4text
|
|
-- ------+---------------+-----------+--------+--------------------------------------------------------------------------
|
|
-- 1 | latlong WGS84 | 1 | | +proj=longlat +datum=WGS84
|
|
-- 2 | BC albers | 2 | | proj=aea ellps=GRS80 lon_0=-126 lat_0=45 lat_1=50 lat_2=58.5 x_0=1000000
|
|
--
|
|
-- select transform( 'SRID=1;POINT(-120.8 50.3)', 2);
|
|
-- -> 'SRID=2;POINT(1370033.37046971 600755.810968684)'
|
|
--
|
|
|
|
CREATE FUNCTION transform_geometry(geometry,text,text,int)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'transform_geom'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION 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);
|
|
|
|
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- COMMON FUNCTIONS
|
|
|
|
CREATE FUNCTION srid(chip)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME,'srid_chip'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION height(chip)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME,'height_chip'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION factor(chip)
|
|
RETURNS FLOAT4
|
|
AS MODULE_FILENAME,'factor_chip'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION width(chip)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME,'width_chip'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION datatype(chip)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME,'datatype_chip'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION compression(chip)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME,'compression_chip'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION setSRID(chip,int4)
|
|
RETURNS chip
|
|
AS MODULE_FILENAME,'setsrid_chip'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION setfactor(chip,float4)
|
|
RETURNS chip
|
|
AS MODULE_FILENAME,'setfactor_chip'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION geometry(CHIP)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'CHIP_to_geom'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION box3d(geometry)
|
|
RETURNS box3d
|
|
AS MODULE_FILENAME,'get_bbox_of_geometry'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION box(geometry)
|
|
RETURNS BOX
|
|
AS MODULE_FILENAME,'geometry2box'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION geometry(box3d)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'get_geometry_of_bbox'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION geometry(text)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_text'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION expand(box3d,float8)
|
|
RETURNS box3d
|
|
AS MODULE_FILENAME,'expand_bbox'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
--
|
|
-- Functions for converting to WKB
|
|
--
|
|
|
|
CREATE FUNCTION asbinary(geometry)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'asbinary_simple'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION asbinary(geometry,TEXT)
|
|
RETURNS wkb
|
|
AS MODULE_FILENAME,'asbinary_specify'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION bytea(wkb)
|
|
RETURNS bytea
|
|
AS MODULE_FILENAME,'WKBtoBYTEA'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION geometry(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'geometryfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION GeomFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'geometryfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION GeomFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'geometryfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION PointFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'PointfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION PointFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'PointfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION LineFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'LinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION LineFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'LinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
|
|
CREATE FUNCTION LinestringFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'LinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION LinestringFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'LinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION PolyFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'PolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION PolyFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'PolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION PolygonFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'PolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION PolygonFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'PolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
|
|
CREATE FUNCTION MPointFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPointfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MPointFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPointfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
|
|
CREATE FUNCTION MultiPointFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPointfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MultiPointFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPointfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MultiLineFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MLinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MultiLineFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MLinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
|
|
CREATE FUNCTION MLineFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MLinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MLineFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MLinefromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MPolyFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MPolyFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MultiPolyFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION MultiPolyFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'MPolyfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
|
|
|
|
CREATE FUNCTION GeomCollFromWKB(wkb,int)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'GCfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
CREATE FUNCTION GeomCollFromWKB(wkb)
|
|
RETURNS GEOMETRY
|
|
AS MODULE_FILENAME,'GCfromWKB_SRID'
|
|
LANGUAGE 'C' WITH (iscachable,isstrict);
|
|
|
|
|
|
-- CREATE FUNCTION index_thing(geometry)
|
|
-- RETURNS BOOL
|
|
-- AS MODULE_FILENAME
|
|
-- LANGUAGE 'C' WITH (isstrict);
|
|
|
|
--
|
|
-- Debugging functions
|
|
--
|
|
|
|
CREATE FUNCTION npoints(geometry)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION nrings(geometry)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict) ;
|
|
|
|
CREATE FUNCTION mem_size(geometry)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
CREATE FUNCTION summary(geometry)
|
|
RETURNS text
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION translate(geometry,float8,float8,float8)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict) ;
|
|
|
|
CREATE FUNCTION dimension(geometry)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict) ;
|
|
|
|
CREATE FUNCTION geometrytype(geometry)
|
|
RETURNS text
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION envelope(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION x(geometry)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'x_point'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION y(geometry)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'y_point'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION z(geometry)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'z_point'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION numpoints(geometry)
|
|
RETURNS integer
|
|
AS MODULE_FILENAME,'numpoints_linestring'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION pointn(geometry,integer)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'pointn_linestring'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION exteriorring(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'exteriorring_polygon'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION numinteriorrings(geometry)
|
|
RETURNS integer
|
|
AS MODULE_FILENAME,'numinteriorrings_polygon'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION interiorringn(geometry,integer)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'interiorringn_polygon'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION numgeometries(geometry)
|
|
RETURNS integer
|
|
AS MODULE_FILENAME,'numgeometries_collection'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometryn(geometry,integer)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometryn_collection'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION max_distance(geometry,geometry)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION optimistic_overlap(geometry,geometry,FLOAT8)
|
|
RETURNS BOOL
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION segmentize(geometry,FLOAT8)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION distance(geometry,geometry)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION astext(geometry)
|
|
RETURNS TEXT
|
|
AS MODULE_FILENAME,'astext_geometry'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION srid(geometry)
|
|
RETURNS int4
|
|
AS MODULE_FILENAME,'srid_geom'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometryfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
|
|
CREATE FUNCTION geometryfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION geomfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION geomfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION polyfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_poly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION polygonfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_poly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION polygonfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_poly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION mpolyfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION linefromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_line'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
|
|
CREATE FUNCTION mlinefromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mline'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION multilinestringfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mline'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION multilinestringfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mline'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION pointfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_point'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION mpointfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoint'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION multipointfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoint'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION multipointfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoint'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION geomcollfromtext(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_gc'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION setSRID(geometry,int4)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION polyfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_poly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
|
|
CREATE FUNCTION mpolyfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION multipolygonfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION multipolygonfromtext(geometry,int)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoly'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION linefromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_line'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION linestringfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_line'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION linestringfromtext(geometry,int)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_line'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION mlinefromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mline'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION pointfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_point'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION mpointfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_mpoint'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION geomcollfromtext(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geometry_from_text_gc'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
|
|
CREATE FUNCTION isempty(geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME,'isempty'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION issimple(geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME,'issimple'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
|
|
CREATE FUNCTION equals(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME,'geomequals'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
|
|
--
|
|
-- Special spheroid functions
|
|
--
|
|
|
|
CREATE FUNCTION length_spheroid(geometry,spheroid)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'length_ellipsoid'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION length3d_spheroid(geometry,spheroid)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'length3d_ellipsoid'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION distance_spheroid(geometry,geometry,spheroid)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'distance_ellipsoid'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
--
|
|
-- Generic operations
|
|
--
|
|
|
|
CREATE FUNCTION length3d(geometry)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION length(geometry)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'length2d'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION area2d(geometry)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION area(geometry)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'area2d'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION perimeter3d(geometry)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION perimeter(geometry)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'perimeter2d'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
---CREATE FUNCTION truly_inside(geometry,geometry)
|
|
--- RETURNS bool
|
|
--- AS MODULE_FILENAME
|
|
--- LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION point_inside_circle(geometry,float8,float8,float8)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION startpoint(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION endpoint(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION isclosed(geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION centroid(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION isring(geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION pointonsurface(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
--
|
|
-- BBox operations
|
|
--
|
|
|
|
CREATE FUNCTION xmin(box3d)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'box3d_xmin'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION ymin(box3d)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'box3d_ymin'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION zmin(box3d)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'box3d_zmin'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION xmax(box3d)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'box3d_xmax'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION ymax(box3d)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'box3d_ymax'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION zmax(box3d)
|
|
RETURNS FLOAT8
|
|
AS MODULE_FILENAME,'box3d_zmax'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
CREATE FUNCTION box3dtobox(box3d)
|
|
RETURNS BOX
|
|
AS MODULE_FILENAME,'box3dtobox'
|
|
LANGUAGE 'C' WITH (isstrict,iscachable);
|
|
|
|
--
|
|
-- Aggregate functions
|
|
--
|
|
|
|
CREATE FUNCTION geom_accum (geometry[],geometry)
|
|
RETURNS geometry[]
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION combine_bbox(box3d,geometry)
|
|
RETURNS box3d
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE AGGREGATE extent(
|
|
sfunc = combine_bbox,
|
|
basetype = geometry,
|
|
stype = box3d
|
|
);
|
|
|
|
CREATE FUNCTION collector(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE AGGREGATE memcollect(
|
|
sfunc = collector,
|
|
basetype = geometry,
|
|
stype = geometry
|
|
);
|
|
|
|
CREATE FUNCTION collect_garray (geometry[])
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE AGGREGATE collect (
|
|
sfunc = geom_accum,
|
|
basetype = geometry,
|
|
stype = geometry[],
|
|
finalfunc = collect_garray
|
|
);
|
|
|
|
|
|
--
|
|
-- Operator definitions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_overleft(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_overright(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_left(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_right(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_contain(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_contained(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_overlap(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_same(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
--
|
|
-- Sorting functions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_lt(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_le(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_gt(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_ge(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_eq(geometry, geometry)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geometry_cmp(geometry, geometry)
|
|
RETURNS integer
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
--
|
|
-- Two dimensional to three dimensional forces
|
|
--
|
|
|
|
CREATE FUNCTION force_2d(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION force_3d(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
--
|
|
-- Force collection
|
|
--
|
|
|
|
CREATE FUNCTION force_collection(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
--
|
|
-- Operator definitions
|
|
--
|
|
|
|
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 = positionjoinsel
|
|
);
|
|
|
|
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
|
|
);
|
|
|
|
--
|
|
-- 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
|
|
);
|
|
|
|
--
|
|
-- GEOS Functions
|
|
--
|
|
|
|
|
|
CREATE FUNCTION intersection(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'intersection'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION buffer(geometry,float8)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'buffer'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION convexhull(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'convexhull'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
CREATE FUNCTION difference(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'difference'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION boundary(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'boundary'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION symdifference(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'symdifference'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
CREATE FUNCTION symmetricdifference(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'symdifference'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
CREATE FUNCTION GeomUnion(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME,'geomunion'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE AGGREGATE MemGeomUnion (
|
|
basetype = geometry,
|
|
sfunc = geomunion,
|
|
stype = geometry
|
|
);
|
|
|
|
CREATE FUNCTION unite_garray (geometry[])
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE AGGREGATE GeomUnion (
|
|
sfunc = geom_accum,
|
|
basetype = geometry,
|
|
stype = geometry[],
|
|
finalfunc = unite_garray
|
|
);
|
|
|
|
|
|
CREATE FUNCTION relate(geometry,geometry)
|
|
RETURNS text
|
|
AS MODULE_FILENAME,'relate_full'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION relate(geometry,geometry,text)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME,'relate_pattern'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION disjoint(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION touches(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION intersects(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION crosses(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION within(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION contains(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION overlaps(geometry,geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION isvalid(geometry)
|
|
RETURNS boolean
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION geosnoop(geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME, 'GEOSnoop'
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
|
|
--
|
|
-- Algorithms
|
|
--
|
|
|
|
CREATE FUNCTION simplify(geometry, float8)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
CREATE FUNCTION line_interpolate_point(geometry, float8)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C' WITH (isstrict);
|
|
|
|
-----------------------------------------
|
|
-- THESE WERE postgis_sql_VER_end.sql.in
|
|
-----------------------------------------
|
|
|
|
#if USE_VERSION == 71
|
|
|
|
--
|
|
-- GiST support functions
|
|
--
|
|
|
|
CREATE FUNCTION ggeometry_consistent(opaque,geometry,int4)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_compress(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_penalty(opaque,opaque,opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_picksplit(opaque, opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_union(bytea, opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_same(opaque, opaque, opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION rtree_decompress(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- R-Tree support functions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_union(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_inter(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_size(geometry,opaque)
|
|
RETURNS float4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
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 = 'ggeometry_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 = 'ggeometry_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 = 'ggeometry_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 = 'rtree_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 = 'ggeometry_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 = 'ggeometry_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 = 'ggeometry_same';
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIX_GEOMETRY_COLUMNS()
|
|
|
|
CREATE FUNCTION fix_geometry_columns() RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
UPDATE geometry_columns SET attrelid = ( SELECT c.oid AS attrelid FROM pg_class c, pg_attribute a WHERE c.relname = geometry_columns.f_table_name::name AND a.attrelid = c.oid AND a.attname = geometry_columns.f_geometry_column::name
|
|
);
|
|
UPDATE geometry_columns SET varattnum = ( SELECT a.attnum FROM pg_class c, pg_attribute a WHERE c.relname = geometry_columns.f_table_name::name AND a.attrelid = c.oid AND a.attname = geometry_columns.f_geometry_column::name
|
|
);
|
|
|
|
RETURN ''geometry_columns table is now linked to the system tables'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN (
|
|
-- <catalog>, <schema>, <table name>, <column name>,
|
|
-- <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.
|
|
--
|
|
CREATE FUNCTION 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;
|
|
|
|
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'')) )
|
|
THEN
|
|
RAISE EXCEPTION ''Invalid type name - valid ones are:
|
|
GEOMETRY, GEOMETRYCOLLECTION, POINT,
|
|
MULTIPOINT, POLYGON, MULTIPOLYGON,
|
|
LINESTRING, or MULTILINESTRING '';
|
|
return ''fail'';
|
|
END IF;
|
|
|
|
IF ( (new_dim >3) or (new_dim <0) ) THEN
|
|
RAISE EXCEPTION ''invalid dimension'';
|
|
return ''fail'';
|
|
END IF;
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
|
|
'' ADD COLUMN '' || quote_ident(column_name) ||
|
|
'' geometry '';
|
|
|
|
EXECUTE ''INSERT INTO geometry_columns VALUES ('' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal(table_name) || '','' ||
|
|
quote_literal(column_name) || '','' ||
|
|
new_dim || '','' || new_srid || '','' ||
|
|
quote_literal(new_type) || '')'';
|
|
|
|
EXECUTE ''select fix_geometry_columns()'';
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
|
|
'' ADD CHECK (SRID('' || quote_ident(column_name) ||
|
|
'') = '' || new_srid || '')'' ;
|
|
|
|
IF (not(new_type = ''GEOMETRY'')) THEN
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
|
|
'' ADD CHECK (geometrytype('' ||
|
|
quote_ident(column_name) || '')='' ||
|
|
quote_literal(new_type) || '' OR ('' ||
|
|
quote_ident(column_name) || '') is null)'';
|
|
END IF;
|
|
|
|
return ''Geometry column '' || column_name || '' added to table ''
|
|
|| table_name || '' WITH a SRID of '' || new_srid ||
|
|
'' and type '' || new_type;
|
|
END;
|
|
' LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN ( <schema>, <table>, <column>, <srid>, <type>, <dim> )
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
CREATE FUNCTION 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);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- DROPGEOMETRYCOLUMN( <schema name>, <table name>, <column name> )
|
|
--
|
|
-- There is no ALTER TABLE DROP COLUMN command in postgresql
|
|
-- There is no ALTER TABLE DROP CONSTRAINT command in postgresql
|
|
-- So, we:
|
|
-- 1. remove the unwanted geom column reference from the
|
|
-- geometry_columns table
|
|
-- 2. update the table so that the geometry column is all NULLS
|
|
-- This is okay since the CHECK srid(geometry) = <srid> is not
|
|
-- checked if geometry is NULL (the isstrict attribute on srid())
|
|
-- 3. add another constraint that the geometry column must be NULL
|
|
-- This, effectively kills the geometry column
|
|
-- (a) its not in the geometry_column table
|
|
-- (b) it only has nulls in it
|
|
-- (c) you cannot add anything to the geom column because it must be NULL
|
|
--
|
|
-- This will screw up if you put a NOT NULL constraint on the geometry
|
|
-- column, so the first thing we must do is remove this constraint (its a
|
|
-- modification of the pg_attribute system table)
|
|
--
|
|
-- We also check to see if the table/column exists in the geometry_columns
|
|
-- table
|
|
|
|
CREATE FUNCTION DropGeometryColumn(varchar,varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
schema_name alias for $1;
|
|
table_name alias for $2;
|
|
column_name alias for $3;
|
|
myrec RECORD;
|
|
okay boolean;
|
|
BEGIN
|
|
-- first we find out if the column is in the geometry_columns table
|
|
okay = ''f'';
|
|
FOR myrec IN SELECT * from geometry_columns where f_table_schema = schema_name 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;
|
|
|
|
-- ensure the geometry column does not have a NOT NULL attribute
|
|
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);
|
|
|
|
-- remove ref from geometry_columns table
|
|
EXECUTE ''delete from geometry_columns where f_table_schema = '' || quote_literal(database_name) ||
|
|
'' and f_table_name = '' || quote_literal(table_name) ||
|
|
'' and f_geometry_column = '' || 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)'';
|
|
|
|
RETURN table_name || ''.'' || column_name ||'' effectively removed.'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
#elif USE_VERSION == 72
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS()
|
|
|
|
CREATE FUNCTION update_geometry_stats() RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
EXECUTE ''update geometry_columns set attrelid = (select pg_class.oid AS attrelid from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_GEOMETRY_column::name), varattnum = (select pg_attribute.attnum from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_GEOMETRY_column::name)'';
|
|
execute ''update geometry_columns set stats = (build_histogram2d( create_histogram2d(find_extent(f_table_name,f_GEOMETRY_column),40 ),f_table_name::text, f_GEOMETRY_column::text)) '';
|
|
return ''done'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS( <table>, <column> )
|
|
|
|
CREATE FUNCTION update_geometry_stats(varchar,varchar) RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
tablename aliAS for $1;
|
|
columnname aliAS for $2;
|
|
|
|
BEGIN
|
|
EXECUTE ''update geometry_columns set attrelid = (select pg_class.oid AS attrelid from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_GEOMETRY_column::name), varattnum = (select pg_attribute.attnum from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_GEOMETRY_column::name)'';
|
|
execute ''update geometry_columns set stats = (build_histogram2d( create_histogram2d(find_extent(''|| quote_literal(tablename) || '',''||quote_literal(columnname) ||''),40 ),''|| quote_literal(tablename) || ''::text,''||quote_literal(columnname) ||''::text )) WHERE f_table_name=''|| quote_literal(tablename) || ''and f_GEOMETRY_column=''||quote_literal(columnname) ;
|
|
return ''done'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CREATE_HISTOGRAM2D( <box>, <size> )
|
|
--
|
|
-- Returns a histgram with 0s in all the boxes.
|
|
|
|
CREATE FUNCTION create_histogram2d(box3d,int)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'create_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BUILD_HISTOGRAM2D( <histogram2d>, <tablename>, <columnname> )
|
|
--
|
|
|
|
CREATE FUNCTION build_histogram2d (HISTOGRAM2D,text,text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'build_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- EXPLODE_HISTOGRAM2D( <histogram2d>, <tablename> )
|
|
--
|
|
|
|
CREATE FUNCTION explode_histogram2d (HISTOGRAM2D,text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'explode_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ESTIMATE_HISTOGRAM2D( <histogram2d>, <box> )
|
|
--
|
|
|
|
CREATE FUNCTION estimate_histogram2d(HISTOGRAM2D,box)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'estimate_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- POSTGISCOSTESTIMATE()
|
|
--
|
|
|
|
CREATE FUNCTION postgisgistcostestimate(opaque,opaque,opaque,opaque,opaque,opaque,opaque,opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'postgisgistcostestimate'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
--
|
|
-- 7.2 GiST support functions
|
|
--
|
|
|
|
CREATE FUNCTION ggeometry_consistent(opaque,geometry,int4)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION ggeometry_compress(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION gbox_penalty(opaque,opaque,opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION gbox_picksplit(opaque, opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION gbox_union(bytea, opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION gbox_same(box, box, opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION rtree_decompress(opaque)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
--
|
|
-- 7.2 RTREE support functions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_union(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION geometry_inter(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
CREATE FUNCTION geometry_size(geometry,opaque)
|
|
RETURNS float4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE'C';
|
|
|
|
--
|
|
-- Create opclass index binding entries.
|
|
--
|
|
|
|
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 = 'box'));
|
|
|
|
-- 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 = 'ggeometry_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 = 'gbox_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 = 'ggeometry_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 = 'rtree_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 = 'gbox_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 = 'gbox_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 = 'gbox_same';
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIX_GEOMETRY_COLUMNS()
|
|
|
|
CREATE FUNCTION fix_geometry_columns() RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
|
|
UPDATE geometry_columns SET attrelid = ( SELECT c.oid AS attrelid FROM pg_class c, pg_attribute a WHERE c.relname = geometry_columns.f_table_name::name AND a.attrelid = c.oid AND a.attname = geometry_columns.f_geometry_column::name
|
|
);
|
|
UPDATE geometry_columns SET varattnum = ( SELECT a.attnum FROM pg_class c, pg_attribute a WHERE c.relname = geometry_columns.f_table_name::name AND a.attrelid = c.oid AND a.attname = geometry_columns.f_geometry_column::name
|
|
);
|
|
|
|
RETURN ''geometry_columns table is now linked to the system tables'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN (
|
|
-- <catalog>, <schema>, <table name>, <column name>,
|
|
-- <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.
|
|
--
|
|
CREATE FUNCTION 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;
|
|
|
|
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'')) )
|
|
THEN
|
|
RAISE EXCEPTION ''Invalid type name - valid ones are:
|
|
GEOMETRY, GEOMETRYCOLLECTION, POINT,
|
|
MULTIPOINT, POLYGON, MULTIPOLYGON,
|
|
LINESTRING, or MULTILINESTRING '';
|
|
return ''fail'';
|
|
END IF;
|
|
|
|
IF ( (new_dim >3) or (new_dim <0) ) THEN
|
|
RAISE EXCEPTION ''invalid dimension'';
|
|
return ''fail'';
|
|
END IF;
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
|
|
'' ADD COLUMN '' || quote_ident(column_name) ||
|
|
'' geometry '';
|
|
|
|
EXECUTE ''INSERT INTO geometry_columns VALUES ('' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal(table_name) || '','' ||
|
|
quote_literal(column_name) || '','' ||
|
|
new_dim || '','' || new_srid || '','' ||
|
|
quote_literal(new_type) || '')'';
|
|
|
|
EXECUTE ''select fix_geometry_columns()'';
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
|
|
'' ADD CHECK (SRID('' || quote_ident(column_name) ||
|
|
'') = '' || new_srid || '')'' ;
|
|
|
|
IF (not(new_type = ''GEOMETRY'')) THEN
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(table_name) ||
|
|
'' ADD CHECK (geometrytype('' ||
|
|
quote_ident(column_name) || '')='' ||
|
|
quote_literal(new_type) || '' OR ('' ||
|
|
quote_ident(column_name) || '') is null)'';
|
|
END IF;
|
|
|
|
return ''Geometry column '' || column_name || '' added to table ''
|
|
|| table_name || '' WITH a SRID of '' || new_srid ||
|
|
'' and type '' || new_type;
|
|
END;
|
|
' LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN ( <schema>, <table>, <column>, <srid>, <type>, <dim> )
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
CREATE FUNCTION 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);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- DROPGEOMETRYCOLUMN( <schema name>, <table name>, <column name> )
|
|
--
|
|
-- There is no ALTER TABLE DROP COLUMN command in postgresql
|
|
-- There is no ALTER TABLE DROP CONSTRAINT command in postgresql
|
|
-- So, we:
|
|
-- 1. remove the unwanted geom column reference from the
|
|
-- geometry_columns table
|
|
-- 2. update the table so that the geometry column is all NULLS
|
|
-- This is okay since the CHECK srid(geometry) = <srid> is not
|
|
-- checked if geometry is NULL (the isstrict attribute on srid())
|
|
-- 3. add another constraint that the geometry column must be NULL
|
|
-- This, effectively kills the geometry column
|
|
-- (a) its not in the geometry_column table
|
|
-- (b) it only has nulls in it
|
|
-- (c) you cannot add anything to the geom column because it must be NULL
|
|
--
|
|
-- This will screw up if you put a NOT NULL constraint on the geometry
|
|
-- column, so the first thing we must do is remove this constraint (its a
|
|
-- modification of the pg_attribute system table)
|
|
--
|
|
-- We also check to see if the table/column exists in the geometry_columns
|
|
-- table
|
|
|
|
CREATE FUNCTION DropGeometryColumn(varchar,varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
schema_name alias for $1;
|
|
table_name alias for $2;
|
|
column_name alias for $3;
|
|
myrec RECORD;
|
|
okay boolean;
|
|
BEGIN
|
|
-- first we find out if the column is in the geometry_columns table
|
|
okay = ''f'';
|
|
FOR myrec IN SELECT * from geometry_columns where f_table_schema = schema_name 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;
|
|
|
|
-- ensure the geometry column does not have a NOT NULL attribute
|
|
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);
|
|
|
|
-- remove ref from geometry_columns table
|
|
EXECUTE ''delete from geometry_columns where f_table_schema = '' || quote_literal(database_name) ||
|
|
'' and f_table_name = '' || quote_literal(table_name) ||
|
|
'' and f_geometry_column = '' || 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)'';
|
|
|
|
RETURN table_name || ''.'' || column_name ||'' effectively removed.'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
#elif USE_VERSION == 73
|
|
|
|
--
|
|
-- 7.3 explicit casting definitions
|
|
--
|
|
|
|
CREATE CAST ( chip AS geometry ) WITH FUNCTION geometry(chip) 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 geometry ) WITH FUNCTION geometry(box3d) AS IMPLICIT;
|
|
CREATE CAST ( text AS geometry) WITH FUNCTION geometry(text) AS IMPLICIT;
|
|
CREATE CAST ( wkb AS bytea ) WITH FUNCTION bytea(wkb) AS IMPLICIT;
|
|
CREATE CAST ( box3d AS box ) WITH FUNCTION box3dtobox(box3d);
|
|
CREATE CAST ( geometry AS text ) WITH FUNCTION astext(geometry);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS()
|
|
|
|
CREATE FUNCTION update_geometry_stats()
|
|
RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
EXECUTE ''update geometry_columns set attrelid = (select pg_class.oid AS attrelid from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name), varattnum = (select pg_attribute.attnum from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name)'';
|
|
execute ''update geometry_columns set stats = (build_histogram2d( create_histogram2d(find_extent(f_table_name,f_geometry_column),40 ),f_table_name::text, f_geometry_column::text)) '';
|
|
return ''done'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS( <table>, <column> )
|
|
|
|
CREATE FUNCTION update_geometry_stats(varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
tablename aliAS for $1;
|
|
columnname aliAS for $2;
|
|
|
|
BEGIN
|
|
EXECUTE ''update geometry_columns set attrelid = (select pg_class.oid AS attrelid from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name), varattnum = (select pg_attribute.attnum from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name)'';
|
|
execute ''update geometry_columns set stats = (build_histogram2d( create_histogram2d(find_extent(''|| quote_literal(tablename) || '',''||quote_literal(columnname) ||''),40 ),''|| quote_literal(tablename) || ''::text,''||quote_literal(columnname) ||''::text )) WHERE f_table_name=''|| quote_literal(tablename) || ''and f_geometry_column=''||quote_literal(columnname) ;
|
|
return ''done'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CREATE_HISTOGRAM2D( <box>, <size> )
|
|
--
|
|
-- Returns a histgram with 0s in all the boxes.
|
|
|
|
CREATE FUNCTION create_histogram2d(box3d,int)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'create_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BUILD_HISTOGRAM2D( <histogram2d>, <tablename>, <columnname> )
|
|
--
|
|
|
|
CREATE FUNCTION build_histogram2d (histogram2d, text, text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'build_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- EXPLODE_HISTOGRAM2D( <histogram2d>, <tablename> )
|
|
--
|
|
|
|
CREATE FUNCTION explode_histogram2d (histogram2d, text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'explode_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ESTIMATE_HISTOGRAM2D( <histogram2d>, <box> )
|
|
--
|
|
|
|
CREATE FUNCTION estimate_histogram2d(histogram2d,box)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'estimate_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
CREATE FUNCTION postgisgistcostestimate(internal,internal,internal,internal,internal,internal,internal,internal)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'postgisgistcostestimate'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
--
|
|
-- 7.2 GiST support functions
|
|
--
|
|
|
|
CREATE FUNCTION ggeometry_consistent(internal,geometry,int4)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_compress(internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_penalty(internal,internal,internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_picksplit(internal, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_union(bytea, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_same(box, box, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION rtree_decompress(internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- 7.2 RTREE support functions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_union(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_inter(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_size(geometry,internal)
|
|
RETURNS float4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- Create opclass index bindings (PG73)
|
|
--
|
|
|
|
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 ggeometry_consistent (internal, geometry, int4),
|
|
FUNCTION 2 gbox_union (bytea, internal),
|
|
FUNCTION 3 ggeometry_compress (internal),
|
|
FUNCTION 4 rtree_decompress (internal),
|
|
FUNCTION 5 gbox_penalty (internal, internal, internal),
|
|
FUNCTION 6 gbox_picksplit (internal, internal),
|
|
FUNCTION 7 gbox_same (box, box, internal);
|
|
|
|
UPDATE pg_opclass
|
|
SET opckeytype = (select oid from pg_type where typname = 'box')
|
|
WHERE opcname = 'gist_geometry_ops';
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIX_GEOMETRY_COLUMNS()
|
|
--
|
|
-- 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 set to the
|
|
-- empty string any schema value not existing in pg_namespace
|
|
--
|
|
|
|
CREATE FUNCTION fix_geometry_columns() RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
-- it would be better to find the correct schema name
|
|
UPDATE geometry_columns SET f_table_schema = ''public''
|
|
WHERE f_table_schema is NULL or f_table_schema NOT IN
|
|
( SELECT nspname::varchar FROM pg_namespace );
|
|
|
|
UPDATE geometry_columns SET attrelid = (
|
|
SELECT c.oid
|
|
FROM pg_class c, pg_attribute a, pg_namespace n
|
|
WHERE c.relname = geometry_columns.f_table_name::name
|
|
AND a.attrelid = c.oid AND c.relnamespace = n.oid
|
|
AND a.attname = geometry_columns.f_geometry_column::name
|
|
AND n.nspname = geometry_columns.f_table_schema::name
|
|
);
|
|
|
|
UPDATE geometry_columns SET varattnum = (
|
|
SELECT a.attnum
|
|
FROM pg_class c, pg_attribute a, pg_namespace n
|
|
WHERE n.nspname = geometry_columns.f_table_schema::name
|
|
AND c.relname = geometry_columns.f_table_name::name
|
|
AND a.attname = geometry_columns.f_geometry_column::name
|
|
AND a.attrelid = c.oid AND c.relnamespace = n.oid
|
|
);
|
|
|
|
RETURN ''geometry_columns table is now linked to the system tables'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN (
|
|
-- <catalog>, <schema>, <table name>, <column name>,
|
|
-- <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.
|
|
--
|
|
CREATE FUNCTION 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;
|
|
|
|
rec RECORD;
|
|
schema_ok bool;
|
|
real_schema name;
|
|
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'')) )
|
|
THEN
|
|
RAISE EXCEPTION ''Invalid type name - valid ones are:
|
|
GEOMETRY, GEOMETRYCOLLECTION, POINT,
|
|
MULTIPOINT, POLYGON, MULTIPOLYGON,
|
|
LINESTRING, or MULTILINESTRING '';
|
|
RETURN ''fail'';
|
|
END IF;
|
|
|
|
IF ( (new_dim >3) or (new_dim <0) ) THEN
|
|
RAISE EXCEPTION ''invalid dimension'';
|
|
RETURN ''fail'';
|
|
END IF;
|
|
|
|
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;
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD COLUMN '' || quote_ident(column_name) ||
|
|
'' geometry '';
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD CHECK (SRID('' || quote_ident(column_name) ||
|
|
'') = '' || new_srid || '')'' ;
|
|
|
|
IF (not(new_type = ''GEOMETRY'')) THEN
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD CHECK (geometrytype('' ||
|
|
quote_ident(column_name) || '')='' ||
|
|
quote_literal(new_type) || '' OR ('' ||
|
|
quote_ident(column_name) || '') is null)'';
|
|
END IF;
|
|
|
|
EXECUTE ''INSERT INTO geometry_columns VALUES ('' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal(real_schema) || '','' ||
|
|
quote_literal(table_name) || '','' ||
|
|
quote_literal(column_name) || '','' ||
|
|
new_dim || '','' || new_srid || '','' ||
|
|
quote_literal(new_type) || '')'';
|
|
|
|
EXECUTE ''select fix_geometry_columns()'';
|
|
--SELECT fix_geometry_columns();
|
|
|
|
return ''Geometry column '' || column_name || '' added to table ''
|
|
|| real_schema || ''.'' || table_name || '' WITH a SRID of '' || new_srid ||
|
|
'' and type '' || new_type;
|
|
END;
|
|
' LANGUAGE 'plpgsql' WITH (isstrict);
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN ( <schema>, <table>, <column>, <srid>, <type>, <dim> )
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
CREATE FUNCTION 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);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- DROPGEOMETRYCOLUMN( <schema name>, <table name>, <column name> )
|
|
--
|
|
-- Removes geometry column reference from geometry_columns table,
|
|
-- and actually drops the column.
|
|
--
|
|
CREATE FUNCTION DropGeometryColumn(varchar,varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
schema_name alias for $1;
|
|
table_name alias for $2;
|
|
column_name alias for $3;
|
|
|
|
myrec RECORD;
|
|
real_schema name;
|
|
okay boolean;
|
|
query text;
|
|
|
|
BEGIN
|
|
|
|
|
|
-- 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;
|
|
|
|
-- first we 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);
|
|
|
|
-- Remove table column
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) || ''.'' ||
|
|
quote_ident(table_name) || '' DROP COLUMN '' ||
|
|
quote_ident(column_name);
|
|
|
|
|
|
RETURN real_schema || ''.'' || table_name || ''.'' || column_name ||'' effectively removed.'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
#elif USE_VERSION == 74
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
-- this has been moved here at _end from _start
|
|
-- because we need the definition of function bytea
|
|
-- found in _common
|
|
CREATE TYPE wkb (
|
|
internallength = variable,
|
|
input = wkb_in,
|
|
output = wkb_out,
|
|
storage = extended,
|
|
send = bytea,
|
|
receive = wkb_recv
|
|
);
|
|
|
|
--
|
|
-- 7.3 explicit casting definitions
|
|
--
|
|
|
|
CREATE CAST ( chip AS geometry ) WITH FUNCTION geometry(chip) 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 geometry ) WITH FUNCTION geometry(box3d) AS IMPLICIT;
|
|
CREATE CAST ( text AS geometry) WITH FUNCTION geometry(text) AS IMPLICIT;
|
|
CREATE CAST ( wkb AS bytea ) WITH FUNCTION bytea(wkb) AS IMPLICIT;
|
|
CREATE CAST ( box3d AS box ) WITH FUNCTION box3dtobox(box3d);
|
|
CREATE CAST ( geometry AS text ) WITH FUNCTION astext(geometry);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS()
|
|
|
|
CREATE FUNCTION update_geometry_stats()
|
|
RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
EXECUTE ''update geometry_columns set attrelid = (select pg_class.oid AS attrelid from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name), varattnum = (select pg_attribute.attnum from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name)'';
|
|
execute ''update geometry_columns set stats = (build_histogram2d( create_histogram2d(find_extent(f_table_name,f_geometry_column),40 ),f_table_name::text, f_geometry_column::text)) '';
|
|
return ''done'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS( <table>, <column> )
|
|
|
|
CREATE FUNCTION update_geometry_stats(varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
tablename aliAS for $1;
|
|
columnname aliAS for $2;
|
|
|
|
BEGIN
|
|
EXECUTE ''update geometry_columns set attrelid = (select pg_class.oid AS attrelid from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name), varattnum = (select pg_attribute.attnum from pg_class,pg_attribute where relname =geometry_columns.f_table_name::name and pg_attribute.attrelid = pg_class.oid and pg_attribute.attname = geometry_columns.f_geometry_column::name)'';
|
|
execute ''update geometry_columns set stats = (build_histogram2d( create_histogram2d(find_extent(''|| quote_literal(tablename) || '',''||quote_literal(columnname) ||''),40 ),''|| quote_literal(tablename) || ''::text,''||quote_literal(columnname) ||''::text )) WHERE f_table_name=''|| quote_literal(tablename) || ''and f_geometry_column=''||quote_literal(columnname) ;
|
|
return ''done'';
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CREATE_HISTOGRAM2D( <box>, <size> )
|
|
--
|
|
-- Returns a histgram with 0s in all the boxes.
|
|
|
|
CREATE FUNCTION create_histogram2d(box3d,int)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'create_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BUILD_HISTOGRAM2D( <histogram2d>, <tablename>, <columnname> )
|
|
--
|
|
|
|
CREATE FUNCTION build_histogram2d (histogram2d, text, text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'build_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- EXPLODE_HISTOGRAM2D( <histogram2d>, <tablename> )
|
|
--
|
|
|
|
CREATE FUNCTION explode_histogram2d (histogram2d, text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'explode_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ESTIMATE_HISTOGRAM2D( <histogram2d>, <box> )
|
|
--
|
|
|
|
CREATE FUNCTION estimate_histogram2d(histogram2d,box)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'estimate_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
CREATE FUNCTION postgisgistcostestimate(internal,internal,internal,internal,internal,internal,internal,internal)
|
|
RETURNS opaque
|
|
AS MODULE_FILENAME,'postgisgistcostestimate'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
--
|
|
-- 7.2 GiST support functions
|
|
--
|
|
|
|
CREATE FUNCTION ggeometry_consistent(internal,geometry,int4)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_compress(internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_penalty(internal,internal,internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_picksplit(internal, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_union(bytea, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_same(box, box, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION rtree_decompress(internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- 7.2 RTREE support functions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_union(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_inter(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_size(geometry,internal)
|
|
RETURNS float4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- Create opclass index bindings (PG74)
|
|
--
|
|
|
|
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 ggeometry_consistent (internal, geometry, int4),
|
|
FUNCTION 2 gbox_union (bytea, internal),
|
|
FUNCTION 3 ggeometry_compress (internal),
|
|
FUNCTION 4 rtree_decompress (internal),
|
|
FUNCTION 5 gbox_penalty (internal, internal, internal),
|
|
FUNCTION 6 gbox_picksplit (internal, internal),
|
|
FUNCTION 7 gbox_same (box, box, internal);
|
|
|
|
UPDATE pg_opclass
|
|
SET opckeytype = (select oid from pg_type where typname = 'box')
|
|
WHERE opcname = 'gist_geometry_ops';
|
|
|
|
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);
|
|
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIX_GEOMETRY_COLUMNS()
|
|
--
|
|
-- 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 set to the
|
|
-- empty string any schema value not existing in pg_namespace
|
|
--
|
|
|
|
CREATE FUNCTION fix_geometry_columns() RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
-- it would be better to find the correct schema name
|
|
UPDATE geometry_columns SET f_table_schema = ''public''
|
|
WHERE f_table_schema is NULL or f_table_schema NOT IN
|
|
( SELECT nspname::varchar FROM pg_namespace );
|
|
|
|
UPDATE geometry_columns SET attrelid = (
|
|
SELECT c.oid
|
|
FROM pg_class c, pg_attribute a, pg_namespace n
|
|
WHERE c.relname = geometry_columns.f_table_name::name
|
|
AND a.attrelid = c.oid AND c.relnamespace = n.oid
|
|
AND a.attname = geometry_columns.f_geometry_column::name
|
|
AND n.nspname = geometry_columns.f_table_schema::name
|
|
);
|
|
|
|
UPDATE geometry_columns SET varattnum = (
|
|
SELECT a.attnum
|
|
FROM pg_class c, pg_attribute a, pg_namespace n
|
|
WHERE n.nspname = geometry_columns.f_table_schema::name
|
|
AND c.relname = geometry_columns.f_table_name::name
|
|
AND a.attname = geometry_columns.f_geometry_column::name
|
|
AND a.attrelid = c.oid AND c.relnamespace = n.oid
|
|
);
|
|
|
|
RETURN ''geometry_columns table is now linked to the system tables'';
|
|
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN (
|
|
-- <catalog>, <schema>, <table name>, <column name>,
|
|
-- <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.
|
|
--
|
|
CREATE FUNCTION 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;
|
|
|
|
rec RECORD;
|
|
schema_ok bool;
|
|
real_schema name;
|
|
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'')) )
|
|
THEN
|
|
RAISE EXCEPTION ''Invalid type name - valid ones are:
|
|
GEOMETRY, GEOMETRYCOLLECTION, POINT,
|
|
MULTIPOINT, POLYGON, MULTIPOLYGON,
|
|
LINESTRING, or MULTILINESTRING '';
|
|
RETURN ''fail'';
|
|
END IF;
|
|
|
|
IF ( (new_dim >3) or (new_dim <0) ) THEN
|
|
RAISE EXCEPTION ''invalid dimension'';
|
|
RETURN ''fail'';
|
|
END IF;
|
|
|
|
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;
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD COLUMN '' || quote_ident(column_name) ||
|
|
'' geometry '';
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD CHECK (SRID('' || quote_ident(column_name) ||
|
|
'') = '' || new_srid || '')'' ;
|
|
|
|
IF (not(new_type = ''GEOMETRY'')) THEN
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD CHECK (geometrytype('' ||
|
|
quote_ident(column_name) || '')='' ||
|
|
quote_literal(new_type) || '' OR ('' ||
|
|
quote_ident(column_name) || '') is null)'';
|
|
END IF;
|
|
|
|
EXECUTE ''INSERT INTO geometry_columns VALUES ('' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal(real_schema) || '','' ||
|
|
quote_literal(table_name) || '','' ||
|
|
quote_literal(column_name) || '','' ||
|
|
new_dim || '','' || new_srid || '','' ||
|
|
quote_literal(new_type) || '')'';
|
|
|
|
EXECUTE ''select fix_geometry_columns()'';
|
|
--SELECT fix_geometry_columns();
|
|
|
|
return ''Geometry column '' || column_name || '' added to table ''
|
|
|| real_schema || ''.'' || table_name || '' WITH a SRID of '' || new_srid ||
|
|
'' and type '' || new_type;
|
|
END;
|
|
' LANGUAGE 'plpgsql' WITH (isstrict);
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN ( <schema>, <table>, <column>, <srid>, <type>, <dim> )
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
CREATE FUNCTION 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);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- DROPGEOMETRYCOLUMN( <schema name>, <table name>, <column name> )
|
|
--
|
|
-- Removes geometry column reference from geometry_columns table,
|
|
-- and actually drops the column.
|
|
--
|
|
CREATE FUNCTION DropGeometryColumn(varchar,varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
schema_name alias for $1;
|
|
table_name alias for $2;
|
|
column_name alias for $3;
|
|
|
|
myrec RECORD;
|
|
real_schema name;
|
|
okay boolean;
|
|
query text;
|
|
|
|
BEGIN
|
|
|
|
|
|
-- 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;
|
|
|
|
-- first we 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);
|
|
|
|
-- Remove table column
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) || ''.'' ||
|
|
quote_ident(table_name) || '' DROP COLUMN '' ||
|
|
quote_ident(column_name);
|
|
|
|
|
|
RETURN real_schema || ''.'' || table_name || ''.'' || column_name ||'' effectively removed.'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
#elif USE_VERSION == 75
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- WKB
|
|
|
|
-- this has been moved here at _end from _start
|
|
-- because we need the definition of function bytea
|
|
-- found in _common
|
|
CREATE TYPE wkb (
|
|
internallength = variable,
|
|
input = wkb_in,
|
|
output = wkb_out,
|
|
storage = extended,
|
|
send = bytea,
|
|
receive = wkb_recv
|
|
);
|
|
|
|
--
|
|
-- 7.3 explicit casting definitions
|
|
--
|
|
|
|
CREATE CAST ( chip AS geometry ) WITH FUNCTION geometry(chip) 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 geometry ) WITH FUNCTION geometry(box3d) AS IMPLICIT;
|
|
CREATE CAST ( text AS geometry) WITH FUNCTION geometry(text) AS IMPLICIT;
|
|
CREATE CAST ( wkb AS bytea ) WITH FUNCTION bytea(wkb) AS IMPLICIT;
|
|
CREATE CAST ( box3d AS box ) WITH FUNCTION box3dtobox(box3d);
|
|
CREATE CAST ( geometry AS text ) WITH FUNCTION astext(geometry);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- UPDATE_GEOMETRY_STATS()
|
|
--
|
|
-- This is no more needed, running ANALYZE will store
|
|
-- geometry attributes statistics in system tables.
|
|
-- It is kept for back compatibility, will raise a message.
|
|
--
|
|
|
|
CREATE FUNCTION update_geometry_stats()
|
|
RETURNS text
|
|
AS
|
|
'
|
|
SELECT ''update_geometry_stats() has been obsoleted. Statistics are automatically built running the ANALYZE command''::text
|
|
'
|
|
LANGUAGE 'sql' ;
|
|
|
|
CREATE FUNCTION update_geometry_stats(varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
SELECT update_geometry_stats();
|
|
'
|
|
LANGUAGE 'sql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- CREATE_HISTOGRAM2D( <box>, <size> )
|
|
--
|
|
-- Returns a histgram with 0s in all the boxes.
|
|
|
|
CREATE FUNCTION create_histogram2d(box3d,int)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'create_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- BUILD_HISTOGRAM2D( <histogram2d>, <tablename>, <columnname> )
|
|
--
|
|
|
|
CREATE FUNCTION build_histogram2d (histogram2d, text, text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'build_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- EXPLODE_HISTOGRAM2D( <histogram2d>, <tablename> )
|
|
--
|
|
|
|
CREATE FUNCTION explode_histogram2d (histogram2d, text)
|
|
RETURNS histogram2d
|
|
AS MODULE_FILENAME,'explode_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ESTIMATE_HISTOGRAM2D( <histogram2d>, <box> )
|
|
--
|
|
|
|
CREATE FUNCTION estimate_histogram2d(histogram2d,box)
|
|
RETURNS float8
|
|
AS MODULE_FILENAME,'estimate_histogram2d'
|
|
LANGUAGE 'C' with (isstrict);
|
|
|
|
--
|
|
-- 7.2 GiST support functions
|
|
--
|
|
|
|
CREATE FUNCTION ggeometry_consistent(internal,geometry,int4)
|
|
RETURNS bool
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION ggeometry_compress(internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_penalty(internal,internal,internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_picksplit(internal, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_union(bytea, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION gbox_same(box, box, internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION rtree_decompress(internal)
|
|
RETURNS internal
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- 7.2 RTREE support functions
|
|
--
|
|
|
|
CREATE FUNCTION geometry_union(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_inter(geometry,geometry)
|
|
RETURNS geometry
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
CREATE FUNCTION geometry_size(geometry,internal)
|
|
RETURNS float4
|
|
AS MODULE_FILENAME
|
|
LANGUAGE 'C';
|
|
|
|
--
|
|
-- Create opclass index bindings (PG75)
|
|
--
|
|
|
|
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 ggeometry_consistent (internal, geometry, int4),
|
|
FUNCTION 2 gbox_union (bytea, internal),
|
|
FUNCTION 3 ggeometry_compress (internal),
|
|
FUNCTION 4 rtree_decompress (internal),
|
|
FUNCTION 5 gbox_penalty (internal, internal, internal),
|
|
FUNCTION 6 gbox_picksplit (internal, internal),
|
|
FUNCTION 7 gbox_same (box, box, internal);
|
|
|
|
UPDATE pg_opclass
|
|
SET opckeytype = (select oid from pg_type where typname = 'box')
|
|
WHERE opcname = 'gist_geometry_ops';
|
|
|
|
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);
|
|
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- FIX_GEOMETRY_COLUMNS()
|
|
--
|
|
-- 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 set to the
|
|
-- empty string any schema value not existing in pg_namespace
|
|
--
|
|
|
|
CREATE FUNCTION fix_geometry_columns() RETURNS text
|
|
AS
|
|
'
|
|
BEGIN
|
|
-- it would be better to find the correct schema name
|
|
UPDATE geometry_columns SET f_table_schema = ''public''
|
|
WHERE f_table_schema is NULL or f_table_schema NOT IN
|
|
( SELECT nspname::varchar FROM pg_namespace );
|
|
|
|
RETURN ''geometry_columns table is now linked to the system tables'';
|
|
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' ;
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN (
|
|
-- <catalog>, <schema>, <table name>, <column name>,
|
|
-- <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.
|
|
--
|
|
CREATE FUNCTION 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;
|
|
|
|
rec RECORD;
|
|
schema_ok bool;
|
|
real_schema name;
|
|
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'')) )
|
|
THEN
|
|
RAISE EXCEPTION ''Invalid type name - valid ones are:
|
|
GEOMETRY, GEOMETRYCOLLECTION, POINT,
|
|
MULTIPOINT, POLYGON, MULTIPOLYGON,
|
|
LINESTRING, or MULTILINESTRING '';
|
|
RETURN ''fail'';
|
|
END IF;
|
|
|
|
IF ( (new_dim >3) or (new_dim <0) ) THEN
|
|
RAISE EXCEPTION ''invalid dimension'';
|
|
RETURN ''fail'';
|
|
END IF;
|
|
|
|
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;
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD COLUMN '' || quote_ident(column_name) ||
|
|
'' geometry '';
|
|
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD CHECK (SRID('' || quote_ident(column_name) ||
|
|
'') = '' || new_srid || '')'' ;
|
|
|
|
IF (not(new_type = ''GEOMETRY'')) THEN
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) ||
|
|
''.'' || quote_ident(table_name) ||
|
|
'' ADD CHECK (geometrytype('' ||
|
|
quote_ident(column_name) || '')='' ||
|
|
quote_literal(new_type) || '' OR ('' ||
|
|
quote_ident(column_name) || '') is null)'';
|
|
END IF;
|
|
|
|
EXECUTE ''INSERT INTO geometry_columns VALUES ('' ||
|
|
quote_literal('''') || '','' ||
|
|
quote_literal(real_schema) || '','' ||
|
|
quote_literal(table_name) || '','' ||
|
|
quote_literal(column_name) || '','' ||
|
|
new_dim || '','' || new_srid || '','' ||
|
|
quote_literal(new_type) || '')'';
|
|
|
|
EXECUTE ''select fix_geometry_columns()'';
|
|
--SELECT fix_geometry_columns();
|
|
|
|
return ''Geometry column '' || column_name || '' added to table ''
|
|
|| real_schema || ''.'' || table_name || '' WITH a SRID of '' || new_srid ||
|
|
'' and type '' || new_type;
|
|
END;
|
|
' LANGUAGE 'plpgsql' WITH (isstrict);
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- ADDGEOMETRYCOLUMN ( <schema>, <table>, <column>, <srid>, <type>, <dim> )
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
CREATE FUNCTION 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);
|
|
|
|
-- - - - - - - - - - - - - - - - - - - - - - - - - - - -
|
|
-- DROPGEOMETRYCOLUMN( <schema name>, <table name>, <column name> )
|
|
--
|
|
-- Removes geometry column reference from geometry_columns table,
|
|
-- and actually drops the column.
|
|
--
|
|
CREATE FUNCTION DropGeometryColumn(varchar,varchar,varchar)
|
|
RETURNS text
|
|
AS
|
|
'
|
|
DECLARE
|
|
schema_name alias for $1;
|
|
table_name alias for $2;
|
|
column_name alias for $3;
|
|
|
|
myrec RECORD;
|
|
real_schema name;
|
|
okay boolean;
|
|
query text;
|
|
|
|
BEGIN
|
|
|
|
|
|
-- 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;
|
|
|
|
-- first we 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);
|
|
|
|
-- Remove table column
|
|
EXECUTE ''ALTER TABLE '' || quote_ident(real_schema) || ''.'' ||
|
|
quote_ident(table_name) || '' DROP COLUMN '' ||
|
|
quote_ident(column_name);
|
|
|
|
|
|
RETURN real_schema || ''.'' || table_name || ''.'' || column_name ||'' effectively removed.'';
|
|
|
|
END;
|
|
'
|
|
LANGUAGE 'plpgsql' WITH (isstrict);
|
|
|
|
|
|
#endif // USE_VERSION == ...
|
|
|
|
END TRANSACTION;
|