Spatial Databases by Open Standards and Software 5

Spatial Databases by Open Standards and Software 5.

Basics of PostGIS

Created by

XMLmind XSL-FO Converter

. Gábor Nagy

Created by

. Spatial Databases by Open Standards and Software 5.: Basics of PostGIS

Gábor Nagy

Created by

. Lector: Zoltán Siki

Created by

. This module was created within TÁMOP - 4.1.2-08/1/A-2009-0027 "Tananyagfejlesztéssel a GEO-ért" ("Educational material development for GEO") project. The project was funded by the European Union and the Hungarian Government to the amount of HUF 44,706,488.

Publication date 2010

Created by

. Copyright

Created by

. 2010 University of West Hungary Faculty of Geoinformatics

Abstract

PostGIS is a geospatial addition for the PostgreSQL, which is based in OGC 06-103r4 and OGC 04-104r4.

Created by

. The right to this intellectual property is protected by the 1999/LXXVI copyright law. Any unauthorized use of this material is prohibited. No part of this product may be reproduced or transmitted in any form or by any means, electronic or mechanical, including photocopying, recording, or by any information storage and retrieval system without express written permission from the author/publisher.

Created by

. Created by

. Table of Contents

Created by

. 5. Basics of PostGIS

......

Created by

. 1. 5.1 Introduction

......

Created by

2. 5.2 Installing PostGIS

Created by

......

Created by

. 2.1. 5.2.1 Installing from compiled binaries

Created by

......

Created by

. 2.2. 5.2.2 Installing from source

......

Created by

. 2.3. 5.2.3 Installing to a database

......

Created by

3. 5.3 Creating spatial tables

Created by

......

Created by

. 4. 5.4 Input and output functions

......

Created by

. 4.1. 5.4.1 GeomFromText

......

Created by

. 4.2. 5.4.2 GeomFromEWKT

......

Created by

4.3. 5.4.3 AsText

Created by

......

Created by

. 4.4. 5.4.4 Binary input and output

......

Created by

. 4.5. 5.4.5 GML and KML input and output

......

Created by

5. 5.5 Calculating numerical attributes of geometry

Created by

......

Created by

. 5.1. 5.5.1 Area

Created by

......

Created by

. 5.2. 5.5.2 Perimeter

......

Created by

. 5.3. 5.5.3 Length

......

Created by

6. 5.6 Calculating geometry values

Created by

......

Created by

. 6.1. 5.6.1 Buffer

......

Created by

. 6.2. 5.6.2 ConvexHull

......

Created by

6.3. 5.6.3 Centroid

Created by

......

Created by

. 6.4. 5.6.4 Simplify

Created by

......

Created by

. 6.5. 5.6.5 Segmentize

......

Created by

. 7. 5.7 Relation of two geometries

......

Created by

7.1. 5.7.1 Distance

Created by

......

Created by

. 7.2. 5.7.2 Azimuth

Created by

......

Created by

. 7.3. 5.7.3 Equals

......

Created by

7.4. 5.7.4 Contains

Created by

......

Created by

7.5. 5.7.5 Overlaps

Created by

......

Created by

. 7.6. 5.7.6 Intersects

......

Created by

. 7.7. 5.7.7 Touches

......

Created by

7.8. 5.7.8 Relate

Created by

......

Created by

. 8. 5.8 Calculate derived geometry

......

Created by

. 8.1. 5.8.1 Transform

......

Created by

8.2. 5.8.2 Affine

Created by

......

Created by

. 8.3. 5.8.3 SnapToGrid

Created by

......

Created by

. 9. 5.9 Set functions

......

Created by

. 9.1. 5.9.1 ST_Union

......

Created by

9.2. 5.9.2 Intersection

Created by

......

Created by

. 9.3. 5.9.3 Difference

Created by

......

Created by

. 9.4. 5.9.4 SymDifference

......

Created by

10. 5.10 Aggregate functions

Created by

......

Created by

10.1. 5.10.1 ST_Union

Created by

......

Created by

. 10.2. 5.10.2 Collect

......

Created by

. 10.3. 5.10.3 Poligonize

......

Created by

11. 5.11 Operators

Created by

......

Created by

. List of Tables

Created by

. 1. Operators of PostGIS

......

Created by

. Chapter 5. Basics of PostGIS

Created by

. 1. 5.1 Introduction

Created by

. The PostGIS is an additional package to PostgreSQL, which provides geographic objects and functions based on OGC 06-103r4 and OGC 06-104r4. The PostGIS (like PostgreSQL) is an open source software, developed by the Refractions Research.

PostGIS uses some other open source softwares. PROJ.4 for the reprojection, and the GEOS (Geometry Engine, Open Source) for many geospatial functions.

Created by

. PostGIS provides geometry and geography type and functions for these types. Most of these functions are based on OGC 06-103r4 and OGC 06-104r4 standards, but PostGIS has some other functions.

PostGIS 1.5.2 contains 10 types, 785 functions and 17 aggregate functions. The functions and aggregate

Created by

. functions are started with “ST_” prefix, but most of them may be used without this prefix. The descriptions of the functions use the name with prefix, where we can not use the function without “ST_” (for example

ST_UNION

Geography type is similar to geometry type. The objects of geography type are determined on a sphere or

Created by

. spheroid. It is different from a geometry type with spherical reference system, because it uses spherical topology instead of plane.

2. 5.2 Installing PostGIS

Created by

. 2.1. 5.2.1 Installing from compiled binaries

Created by

. Under Linux operating system, the simplest solution is to use the package manager to install the PostGIS package of the distribution. For example use the

apt-get install postgis

command in a Debian based system. (The name of the required package may be different.)

Created by

. When we would like to install PostgreSQL to a Windows operating system, we can use more different installer programs. Some installers contain the PostGIS package, and it is selectable during the installation process.

The PostgreSQL and PostGIS could be installed from a ZIP file too. Both of the downloaded files must be unpacked to a directory. The PostgreSQL programs are available from this directory.

Created by

. 2.2. 5.2.2 Installing from source

The PostgreSQL and the PostGIS are open source programs. The source code of these programs can be downloaded from the following home pages (http://www.postgresql.org and http://www.postgis.org).

Created by

. Both of these programs have regular build system, and could be installed by

./configure

Created by

. make

make install

commands. This work needs some other software. We can install this software by the package manager (need the development packages, with -dev ended name) or build from source also.

Created by

. 2.3. 5.2.3 Installing to a database

The installation of PostGIS to the database server is not sufficient, because PostGIS contains stored procedures. These stored procedures (and types) must be available in each database, where we should use the geospatial function of PostGIS.

Created by

. All of the necessary SQL commands, which create the required database objects, are in the

postgis.sql

file. This file is in the

Created by

. share/contrib/postgis-version

directory, where version is the version of PostGIS.

This file can be run from the psql console. For example:

Created by

. \i /usr/local/share/contrib/postgis-1.5/postgis.sql

We can use the -f option of the psql command from the shell:

Created by

. psql -f /usr/local/share/contrib/postgis-1.5/postgis.sql gisdatabase gisuser

After the execution of these commands, PostGIS is available in the database.

Created by

spatial_ref_sys

table was created by the

postgis.sql

Created by

. , but the table is empty. We can fill the descriptions of the spatial reference systems to this table by the

spatial_ref_sys.sql

file, which locates in the directory of the

postgis.sql

Created by

. 3. 5.3 Creating spatial tables

Created by

. PostGIS provides geometry and geography types. Geometry and geography columns can be created by the

CREATE TABLE

SQL command, together with the other columns of the table:

Created by

. gisdb=# CREATE TABLE gistable

(id serial PRIMARY KEY,

name char(20),

Created by

. geom geometry);

CREATE TABLE

Created by

. This is not a perfect solution, because many applications need one record in the geometry

_columns

table for each geometry type columns of the tables of the database. Use the

AddGeomemertyColumns

Created by

. function instead to create the geometry type columns of a table, after the table has been created. This function can be called by a

SELECT

command. For example:

Created by

. gisdb=# CREATE TABLE gistable

(id serial PRIMARY KEY,

name char(20));

Created by

. CREATE TABLE

gisdb=# SELECT AddGeometryColumn('gistable', 'geom', 23700, 'POLYGON', 2);

Created by

. addgeometrycolumn

------

public.gistable.geom SRID:23700 TYPE:POLYGON DIMS:2

Created by

. (1 row)

The first parameter of the function is the name of the table, the second parameter is the name of the new geometry column. The third parameter is the SRID (the identification number of the reference system), the fourth prameter is the type of the geometry, fifth parameter is the dimension of the geometry.

Created by

. AddGeomemertyColumn

function inserts the necessary record into the geometry_columns table, and set constrains to check the geometry type, dimension and SRID (identification number of the spatial reference system).

Created by

. The spatial index is very important for effective spatial queries. GIST type indexes can be built for the geometry type columns:

gisdb=# CREATE INDEX gistable_geom ON gistable USING GIST (geom);

Created by

. CREATE INDEX

The GIST (Generalized Search Tree) index provides an R-tree like data structure for fast queries on the basis of the bounding rectangles. The creation of this index is recommended for each geometry column, because it is very useful in many geospatial operations.

Created by

DropGeometryColumn

function can be used to remove a geometry column from a table. For example:

Created by

. SELECT DropGeometryColumn('gistable','geom');

4. 5.4 Input and output functions

Created by

. These functions provide the creation of geometry objects from text and binary data, and dump the geometry to various text and binary formats.

Created by

. 4.1. 5.4.1 GeomFromText

GeomFromText

Created by

GeometryFromText

function returns a geometry value.

GeomFromText('POLYGON((4 4,4 7,7 7,7 4,4 4))', 23700)

Created by

. The first parameter is the WKT representation of the geometry, the second parameter is the identification number of the spatial reference system (SRID).

Created by

ST_GeogFromText

is similar to the

GeomFromText

, but creates geography values. This function needs only one parameter, the WKT representation of the geography; the SRID will be 4326 (WGS84).

Created by

. 4.2. 5.4.2 GeomFromEWKT

Created by

GeomFromEWKT

is similar to the

GeomFromText

function, but uses EWKT (Extended WKT) representation.

Created by

. GeomFromEWKT('SRID=23700;POLYGON((4 4,4 7,7 7,7 4,4 4))')

This function has only one parameter: The EWKT representation text of the geometry, which contains the SRID too.

Created by

. 4.3. 5.4.3 AsText

Created by

. AsTex

t function returns the WKT representation of the geometry. The

AsEWKT

is very similar, but the result contains the SRID of the geometry too. Both of these functions have one parameter: one geometry value.

Created by

. 4.4. 5.4.4 Binary input and output

PostGIS has functions for the binary input and output. These functions are very similar to the text functions, but

Created by

. use WKB or EWKB format for the input (

GeomFromWKB

GeomFromEWBK

) or the output (

Created by

. AsBinary

AsEWKB

Created by

. 4.5. 5.4.5 GML and KML input and output

PostGIS can handle GML and KML data. The

Created by

. and the

functions return the representation of the geometry in these XML based formats. The result is a text value. This text value contains the XML description of the geometry in the required format:

Created by

. gisdb=# SELECT AsGML(GeomFromText('POLYGON((4 4,4 7,7 7,7 4,4 4))', 4326));

Created by

. ------

Created by

. 4,4 4,7 7,7 7,4 4,4

Created by

(1 row)

gisdb=# SELECT AsKML(GeomFromText('POLYGON((4 4,4 7,7 7,7 4,4 4))', 4326));

Created by

. askml

------

Created by

4,4 4,7 7,7 7,4 4,4

Created by

. (1 row)

Created by

. GeomFromGML

GeomFromKML

functions provide the data input from these formats.

Created by

. 5. 5.5 Calculating numerical attributes of geometry

Created by

. 5.1. 5.5.1 Area

function returns the area of a Polygon or MultiPolygon geometry or the total area of the polygon elements of a GeometryCollection. If the parameter of the function has a different geometry type the result will be zero.

Created by

. gisdb=# SELECT ST_Area('POLYGON((1 1,1 4,4 4,4 1,1 1))');

st_area

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT ST_Area('POLYGON((1 1,1 4,4 4,4 1,1 1),

Created by

. (2 2,2 3,3 3,3 2,2 2))');

st_area

------

Created by

(1 row)

Created by

. gisdb=# SELECT ST_Area('LINESTRING(4 5,7 9,11 2)');

st_area

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT ST_Area('GEOMETRYCOLLECTION(

Created by

. POLYGON((1 1,1 4,4 4,4 1,1 1)),POINT(7 8))');

st_area

------

Created by

(1 row)

Created by

. The parameter of the

ST_Area

function (as in the examples above) may be a text value with the WKT representation of the geometry too. This solution will be applied in some of the following examples of other functions.

Created by

. 5.2. 5.5.2 Perimeter

Perimeter

Created by

. function returns the perimeter (length of the boundary) of a Polygon or MultiPolygon geometry or the the total perimeter of the polygon elements of a GeometryCollection. If the parameter of the function has a different geometry type the result will be zero.

gisdb=# SELECT Perimeter('POLYGON((1 1,1 4,4 4,4 1,1 1))');

Created by

. perimeter

------

Created by

. (1 row)

gisdb=# SELECT Perimeter('POLYGON((1 1,1 4,4 4,4 1,1 1),

Created by

. (2 2,2 3,3 3,3 2,2 2))');

perimeter

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Perimeter('LINESTRING(4 5,7 9,11 2)');

Created by

. perimeter

------

Created by

. (1 row)

gisdb=# SELECT Perimeter('GEOMETRYCOLLECTION(

Created by

. POLYGON((1 1,1 4,4 4,4 1,1 1)),POINT(7 8))');

perimeter

Created by

. ------

(1 row)

Created by

Perimeter3D

function works similarly to the Perimeter function, but calculates the 3D perimeter, if the geometry has Z coordinates. The

Created by

. Perimeter2D

function always returns the 2D perimeter.

gisdb=# SELECT Perimeter2D('POLYGON((1 1 0,1 4 0,4 4 4,

Created by

. 4 1 0,1 1 0))');

perimeter2d

Created by

. ------

(1 row)

Created by

. gisdb=# SELECT Perimeter3D('POLYGON((1 1 0,1 4 0,4 4 4,

4 1 0,1 1 0))');

Created by

. perimeter3d

------

Created by

(1 row)

Created by

. 5.3. 5.5.3 Length

Length

function returns the length of a LineString or MultiLineString geometry or the the total length of the linestring

Created by

. elements of a GeometryCollection. If the parameter of the function has a different geometry type the result will be zero.

gisdb=# SELECT ST_Length('LINESTRING(1 2,5 2,9 5)');

Created by

. st_length

------

Created by

. (1 row)

gisdb=# SELECT ST_Length('MULTILINESTRING((1 2,5 2,9 5),

Created by

. (1 4,3 4))');

st_length

Created by

. ------

(1 row)

Created by

. gisdb=# SELECT ST_Length('POLYGON((1 1,1 4,4 4,4 1,1 1))');

st_length

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT ST_Length('GEOMETRYCOLLECTION(

Created by

. LINESTRING(1 4,1 10),POINT(7 8))');

st_length

Created by

. ------

(1 row)

Created by

. The “length” is also the name of the function, which returns the length of a text. PostgreSQL makes a distinction between the two length functions by the type of the parameter. We have to use

GeomFromText

function in the

Length

Created by

. function, or use

ST_Length

function.

Created by

. Length3D

function works similarly to the Length function, but calculates the 3D perimeter, if the geometry has Z coordinates.

gisdb=# SELECT ST_Length('LINESTRING(1 3 2,5 15 5)');

Created by

. st_length

------

Created by

. 12.6491106406735

(1 row)

Created by

. gisdb=# SELECT ST_Length2D('LINESTRING(1 3 2,5 15 5)');

st_length2d

------

Created by

. 12.6491106406735

(1 row)

Created by

. gisdb=# SELECT ST_Length3D('LINESTRING(1 3 2,5 15 5)');

st_length3d

Created by

. ------

Created by

. (1 row)

Length_Spheroid

Lebgth3D_Spheroid

Created by

. ) function can calculate the length of the geometry on a spheroid (sphere or ellipsoid).

gisdb=# SELECT ST_Length_Spheroid('LINESTRING(18.41281 47.18458,

Created by

. 18.41883 47.18878,18.42706 47.19097)',

'SPHEROID["WGS84",6378137,298.257223563]');

st_length_spheroid

Created by

. ------

1322.39548875418

Created by

. (1 row)

6. 5.6 Calculating geometry values

Created by

. These functions return a new geometry, which are built from the original geometries and sometimes other parameters.

Created by

. 6.1. 5.6.1 Buffer

Buffer

function returns a Polygon or MultiPolygon geometry that represents the set of the points whose distance from a geometry specified by the first parameter are less than or equal to the distance specified by the second parameter.

Created by

. gisdb=# SELECT AsText(Buffer('POINT(3 4)',1));

astext

Created by

. ------

POLYGON((4 4,3.981 3.805,3.924 3.617,3.831 3.444,3.707 3.293,

Created by

. 3.556 3.169,3.383 3.076,3.195 3.019,3 3,

2.805 3.019,2.617 3.076,2.444 3.169,2.293 3.293,

2.169 3.444,2.076 3.617,2.019 3.805,2 4,

Created by

. 2.019 4.195,2.076 4.383,2.169 4.556,2.293 4.707,

2.444 4.831,2.617 4.924,2.805 4.981,3 5,

Created by

. 3.195 4.981,3.383 4.924,3.556 4.831,3.707 4.707,

3.831 4.556,3.924 4.383,3.981 4.195,4 4))

(1 row)

Created by

. By definition, the boundary of the shape of the buffer contains arc sections, but the Polygon and MultiPolygon geometries are limited to have linear rings. The arc sections will be replaced by more linear sections. By default, a quarter of a circle will be replaced by 8 line segments. The optionally third parameter of the function specifies the number of line segments or other parameters of the buffer generation.

Created by

. gisdb=# SELECT AsText(ST_Buffer('POINT(3 4)',1,'quad_segs=2'));

astext

Created by

. ------

POLYGON((4 4,3.707 3.293,3 3,2.293 3.293,2 4,

2.293 4.707,3 5,3.707 4.707,4 4))

Created by

. (1 row)

6.2. 5.6.2 ConvexHull

Created by

. ConvexHull

function returns a geometry that represents the least convex point set, which contains all the points of a geometry, which were specified by the parameter of the function.

Created by

. gisdb=# SELECT AsText(ConvexHull('POINT(3 5)'));

astext

Created by

. ------

POINT(3 5)

(1 row)

Created by

. gisdb=# SELECT AsText(ConvexHull('LINESTRING(2 3,5 8)'));

astext

Created by

. ------

LINESTRING(2 3,5 8)

Created by

. (1 row)

gisdb=# SELECT AsText(ConvexHull('MULTIPOINT(2 2,4 5)'));

Created by

. astext

------

LINESTRING(2 2,4 5)

Created by

. (1 row)

Created by

. gisdb=# SELECT AsText(ConvexHull('MULTIPOINT(2 2,4 5,3 4)'));

astext

------

Created by

. POLYGON((2 2,3 4,4 5,2 2))

(1 row)

Created by

. gisdb=# SELECT AsText(ConvexHull('POLYGON((1 1,1 4,4 4,4 1,1 1),

(2 2,2 3,3 3,3 2,2 2))'));

Created by

. astext

------

POLYGON((1 1,1 4,4 4,4 1,1 1))

Created by

. (1 row)

gisdb=# SELECT AsText(ConvexHull('POLYGON((1 5,6 9,4 6,8 2,1 5))'));

Created by

. astext

------

Created by

. POLYGON((8 2,1 5,6 9,8 2))

(1 row)

Created by

. 6.3. 5.6.3 Centroid

Centroid

Created by

. function returns the geometric centre of the geometry as a Point geometry.

gisdb=# SELECT AsText(Centroid('POLYGON((1 1,1 4,4 4,4 1,1 1))'));

astext

Created by

. ------

POINT(2.5 2.5)

Created by

. (1 row)

gisdb=# SELECT AsText(Centroid('POLYGON((1 1,1 4,4 4,4 1,1 1),

Created by

. (2 2,2 3,3 3,3 2,2 2))'));

astext

Created by

. ------

POINT(2.5 2.5)

(1 row)

Created by

. 6.4. 5.6.4 Simplify

Created by

. Simplify

function returns a simplified variant of the geometry, which is the first parameter of the function. This function uses the Douglas-Peuker algorithm, with a tolerance, which is the second parameter.

gisdb=# SELECT AsText(Simplify('LINESTRING(2 4,4 7,7 9)', 1));

Created by

. astext

------

Created by

. LINESTRING(2 4,7 9)

(1 row)

Created by

. gisdb=# SELECT AsText(Simplify('LINESTRING(2 4,4 7,7 9)', 0.5));

astext

------

Created by

. LINESTRING(2 4,4 7,7 9)

(1 row)

Created by

. ST_SiplifyPreserveTopology

function is the same as the

Simplify

, but preserves the topology relations of parts of the object (linear rings, linestring segments).

Created by

. SELECT AsText(Simplify(ST_GeomFromText(

'MULTILINESTRING((1 5,5 7,9 5),(1 1,5 6,9 1))'),2));

Created by

. astext

------

MULTILINESTRING((1 5,9 5),(1 1,5 6,9 1))

Created by

. (1 row)

SELECT AsText(ST_SimplifyPreserveTopology(ST_GeomFromText(

Created by

. 'MULTILINESTRING((1 5,5 7,9 5),(1 1,5 6,9 1))'),2));

astext

------

Created by

. MULTILINESTRING((1 5,5 7,9 5),(1 1,5 6,9 1))

(1 row)

Created by

. 6.5. 5.6.5 Segmentize

Segmentize

Created by

. function returns a modified variant of the geometry (the first parameter), which has no longer segments than a given distance (the second parameter).

gisdb=# SELECT AsText(Segmentize('LINESTRING(1 3,1 6)',3));

Created by

. astext

------

LINESTRING(1 3,1 6)

Created by

. (1 row)

Created by

. gisdb=# SELECT AsText(Segmentize('LINESTRING(1 3,1 9)',3));

astext

------

Created by

. LINESTRING(1 3,1 6,1 9)

(1 row)

Created by

. gisdb=# SELECT AsText(Segmentize('LINESTRING(1 3,1 9)',2));

astext

Created by

. ------

LINESTRING(1 3,1 5,1 7,1 9)

(1 row)

Created by

. gisdb=# SELECT AsText(Segmentize(

'POLYGON((3 2,3 8,5 8,5 2,3 2))',3));

Created by

. astext

------

Created by

. POLYGON((3 2,3 5,3 8,5 8,5 5,5 2,3 2))

(1 row)

Created by

. 7. 5.7 Relation of two geometries

Created by

. 7.1. 5.7.1 Distance

Distance

function returns the shortest distance between two geometries, which are the parameters of the function. The result is zero, if the geometries have a common point.

Created by

. gisdb=# SELECT Distance('POINT(1 2)','POINT(4 6)');

distance

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Distance('LINESTRING(2 2,8 2)','POINT(4 6)');

Created by

. distance

------

Created by

. (1 row)

Created by

. gisdb=# SELECT Distance('POLYGON((2 2,2 3,8 3,8 2,2 2))',

'POINT(4 6)');

distance

Created by

. ------

Created by

. (1 row)

Distance

function can be used between geography objects too.

Created by

. gisdb=# SELECT ST_Distance(

ST_GeogFromText('SRID=4326;POINT(18.41883 47.18878)'),

Created by

. ST_GeogFromText('SRID=4326;POINT(18.41281 47.18458)'));

st_distance

------

Created by

. 652.839475007109

(1 row)

Created by

. ST_ShortestLine

function returns the shortest line segment between two geometries as a LineString geometry. The length of this LineString equals to the result of the

Distance

function.

Created by

. gisdb=# SELECT AsText(ST_ShortestLine('POINT(1 2)','POINT(4 6)'));

astext

Created by

. ------

LINESTRING(1 2,4 6)

Created by

. (1 row)

gisdb=# SELECT AsText(ST_ShortestLine('LINESTRING(2 2,8 2)',

Created by

. 'POINT(4 6)'));

astext

------

Created by

(1 row)

Created by

. gisdb=# SELECT AsText(ST_ShortestLine(

'POLYGON((2 2,2 3,8 3,8 2,2 2))','POINT(4 6)'));

Created by

. astext

------

Created by

(1 row)

Created by

. 7.2. 5.7.2 Azimuth

Azimuth

function returns the azimuth in radians from the point specified in the first parameter to the point specified in the second parameter. If any parameter has no Point geometry, the function returns NULL.

Created by

. gisdb=# SELECT Azimuth('POINT(1 1)','POINT(1 2)');

azimuth

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Azimuth('POINT(1 1)','POINT(2 2)');

Created by

. azimuth

------

Created by

. 0.785398163397448

(1 row)

Created by

. 7.3. 5.7.3 Equals

Equals

function returns true value, if the point set of the geometries (specified by the parameters of the function) are spatially equal.

Created by

. gisdb=# SELECT Equals('POINT(2 4)','POINT(3 5)');

equals

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Equals('POINT(2 4)','POINT(2 4)');

Created by

. equals

------

Created by

. (1 row)

Created by

. gisdb=# SELECT Equals('MULTIPOINT(2 4,5 8)','MULTIPOINT(5 8,2 4)');

equals

------

Created by

(1 row)

Created by

. gisdb=# SELECT Equals('LINESTRING(2 4,5 8)','LINESTRING(5 8,2 4)');

equals

Created by

. ------

(1 row)

Created by

. gisdb=# SELECT Equals('POLYGON((1 1,1 5,5 1,1 1))',

'POLYGON((1 5,1 1,5 1,1 5))');

Created by

. equals

------

Created by

(1 row)

Created by

. gisdb=# SELECT Equals('POLYGON((2 2,2 4,4 4,4 2,2 2))',

'POLYGON((2 2,2 3,2 4,4 4,4 2,2 2))');

Created by

. equals

------

Created by

. (1 row)

7.4. 5.7.4 Contains

Created by

. Contains

function returns true value, if the geometry specified by the first parameter contains all the points of the geometry specified by the second parameter.

Created by

. gisdb=# SELECT Contains('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POINT(3 3)');

Created by

. contains

------

Created by

. (1 row)

gisdb=# SELECT Contains('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'POINT(8 4)');

contains

Created by

. ------

(1 row)

Created by

'LINESTRING(2 2,3 4)');

Created by

. contains

------

Created by

(1 row)

Created by

. contains

------

Created by

. (1 row)

gisdb=# SELECT Contains('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'POLYGON((2 2,2 3,3 3,3 2,2 2))');

contains

Created by

. ------

(1 row)

Created by

. 'POLYGON((4 4,4 7,7 7,7 4,4 4))');

contains

------

Created by

(1 row)

Created by

. gisdb=# SELECT Contains('LINESTRING(2 3,4 3)',

'POINT(3 3)');

Created by

. contains

------

Created by

. (1 row)

gisdb=# SELECT Contains('LINESTRING(2 3,7 3)',

Created by

. 'LINESTRING(3 3,4 3)');

contains

Created by

. ------

Created by

. (1 row)

7.5. 5.7.5 Overlaps

Created by

. Overlaps

function returns true value, if the common set of the points of the two geometries has the same dimension as the objects. This function returns true value, if the two geometries have common points, but both of the geometries have points, which are not contained by the other object. This function does not work with GeometryCollection objects.

Created by

. gisdb=# SELECT Overlaps('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((2 2,2 3,3 3,3 2,2 2))');

overlaps

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Overlaps('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'POLYGON((4 4,4 7,7 7,7 4,4 4))');

overlaps

Created by

. ------

(1 row)

Created by

. gisdb=# SELECT Overlaps('LINESTRING(2 3,7 3)',

Created by

. overlaps

------

Created by

(1 row)

Created by

. gisdb=# SELECT Overlaps('LINESTRING(2 3,7 3)',

overlaps

Created by

. ------

Created by

. (1 row)

7.6. 5.7.6 Intersects

Created by

. Intersects

function returns true value, if the two geometries have at least one common point.

Created by

. gisdb=# SELECT Intersects('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((4 4,4 7,7 7,7 4,4 4))');

intersects

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Intersects('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'LINESTRING(2 4,6 4)'); intersects

------

Created by

(1 row)

Created by

. gisdb=# SELECT Intersects('POLYGON((1 1,1 5,5 5,5 1,1 1))',

intersects

Created by

. ------

Created by

. (1 row)

7.7. 5.7.7 Touches

Created by

. Touches

function returns true value, if the two geometries have common points, but these points are located only on the boundaries of the objects.

Created by

. gisdb=# SELECT Touches('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((5 2,5 4,7 4,7 2,5 2))');

touches

Created by

. ------

Created by

. (1 row)

gisdb=# SELECT Touches('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'POLYGON((5 2,4 4,7 4,7 2,5 2))');

touches

Created by

. ------

(1 row)

Created by

. 7.8. 5.7.8 Relate

Created by

. Relate

function returns a text, which contains the element of the DE-9IM.

gisdb=# SELECT Relate('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'POLYGON((5 2,5 4,7 4,7 2,5 2))');

relate

------

Created by

. FF2F11212

(1 row)

Created by

. If the

Relate

function has three parameters, the third parameter is a DE-9IM pattern. The function returns true, if the pattern matches the relation of the two geometries, given in the first and second parameters.

Created by

. gisdb=# SELECT Relate('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((5 2,5 4,7 4,7 2,5 2))','****1****');

Created by

. relate

------

Created by

. (1 row)

8. 5.8 Calculate derived geometry

Created by

. These functions return a modified geometry preserving the topology of the original object.

Created by

. 8.1. 5.8.1 Transform

PostGIS can transform geometries between spatial reference systems. The geometries contain the identification number of spatial reference system, further on SRID. The spatial_ref_sys table contains the descriptions of spatial reference systems identified by the SRID.

Created by

. Transform

function returns a transformed geometry from the original reference system (specified as the first parameter, SRID is involved in the geometry) to a new reference system, which SRID is specified by the second parameter.

Created by

. gisdb=# SELECT AsText(Transform(

'SRID=4326;POINT(18.4183 47.1888)',23700));

Created by

. astext

------

POINT(602233.50006042 205129.440451882)

Created by

. (1 row)

4326 is the SRID of WGS84, 23700 is the SRID of Hungarian Projection System (EOV).

Created by

. The SRID of an object can be modified without transformation by the

SetSRID

function.

Created by

. gisdb=# SELECT AsText(ST_SetSRID

('SRID=4326;POINT(18.4183 47.1888)',23700));

Created by

. astext

------

POINT(18.4183 47.1888)

Created by

. (1 row)

8.2. 5.8.2 Affine

Created by

. Affine

function returns an affine transformed geometry from the input geometry specified as the first parameter. The other 6 or 12 parameters of the function contain the parameters of a 2D or 3D transformation.

Created by

. 8.3. 5.8.3 SnapToGrid

SnapToGrid

Created by

. function returns modified geometry snapping points to a grid. The resolution of the grid is specified by the second parameter. The coordinates of a geometry can be rounded up by this function.

gisdb=# SELECT AsText(SnapToGrid(

Created by

. 'POINT(602233.50006042 205129.440451882)',0.01));

astext

------

Created by

. POINT(602233.5 205129.44)

(1 row)

Created by

. 9. 5.9 Set functions

These functions return a geometry which represents a point set which is the result of a set operation between two geometries.

Created by

. 9.1. 5.9.1 ST_Union

Created by

function returns a geometry that represents the union of the two geometries, which are specified as the parameters of the function.

Created by

. gisdb=# SELECT AsText(ST_Union('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((3 3,3 7,7 7,7 3,3 3))'));

Created by

. astext

------

POLYGON((1 1,1 5,3 5,3 7,7 7,7 3,5 3,5 1,1 1))

Created by

. (1 row)

Created by

. Figure 1. Union of two POLYGON geometries

gisdb=# SELECT AsText(ST_Union('LINESTRING(2 3,5 3)',

'LINESTRING(4 3,7 3)')); astext

Created by

. ------

MULTILINESTRING((2 3,4 3),(4 3,5 3),(5 3,7 3))

Created by

. (1 row)

gisdb=# SELECT AsText(ST_Union('POINT(2 4)','POINT(7 3)'));

Created by

. astext

------

Created by

. MULTIPOINT(2 4,7 3)

(1 row)

Created by

. 9.2. 5.9.2 Intersection

Intersection

function returns a geometry that represents the common parts of the two geometries, which are specified as the parameters of the function.

Created by

. gisdb=# SELECT AsText(Intersection('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((3 3,3 7,7 7,73,3 3))'));

Created by

. astext

------

Created by

. POLYGON((3 5,5 5,5 3,3 3,3 5))

(1 row)

Created by

. Figure 2. Intersection of two POLYGON geometries

gisdb=# SELECT AsText(Intersection('LINESTRING(2 3,5 3)',

Created by

. 'LINESTRING(4 3,7 3)')); astext

------

Created by

(1 row)

Created by

. gisdb=# SELECT AsText(Intersection('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'LINESTRING(3 4,8 9)'));

astext

Created by

. ------

LINESTRING(3 4,4 5)

Created by

. (1 row)

gisdb=# SELECT AsText(Intersection('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'MULTIPOINT(2 4,6 2)'));

astext

Created by

. ------

POINT(2 4)

(1 row)

Created by

. 9.3. 5.9.3 Difference

Created by

. Difference

function returns a geometry that represents that part of the first geometry, which is not in the second geometry.

gisdb=# SELECT AsText(Difference('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'POLYGON((3 3,3 7,7 7,7 3,3 3))'));

astext

------

Created by

. POLYGON((1 1,1 5,3 5,3 3,5 3,5 1,1 1))

(1 row)

Created by

. Figure 3. Difference of two POLYGON geometries

gisdb=# SELECT AsText(Difference('LINESTRING(2 3,5 3)',

'LINESTRING(4 3,7 3)')); astext

Created by

. ------

LINESTRING(2 3,4 3)

Created by

. (1 row)

gisdb=# SELECT AsText(Difference('POLYGON((1 1,1 5,5 5,5 1,1 1))',

Created by

. 'LINESTRING(3 4,8 9)')); astext

------

Created by

. POLYGON((1 1,1 5,4 5,5 5,5 1,1 1))

(1 row)

Created by

. gisdb=# SELECT AsText(Difference('POLYGON((1 1,1 5,5 5,5 1,1 1))',

'MULTIPOINT(2 4,6 2)')); astext

------

Created by

. POLYGON((1 1,1 5,5 5,5 1,1 1))

(1 row)

Created by

. 9.4. 5.9.4 SymDifference

SymDifference

Created by

. function returns a geometry that represents that part of the two input geometry, which belongs to only one of them.

gisdb=# SELECT AsText(SymDifference(

Created by

. 'POLYGON((1 1,1 5,5 5,5 1,1 1))',

'POLYGON((3 3,3 7,7 7,7 3,3 3))'));

astext

Created by

. ------

MULTIPOLYGON(((1 1,1 5,3 5,3 3,5 3,5 1,1 1)),

Created by

. ((5 3,5 5,3 5,3 7,7 7,7 3,5 3)))

(1 row)

Created by

. Figure 4. Symmetric difference of two POLYGON geometries

gisdb=# SELECT AsText(SymDifference('LINESTRING(2 3,5 3)',

Created by

. 'LINESTRING(4 3,7 3)'));

astext

Created by

. ------

MULTILINESTRING((2 3,4 3),(5 3,7 3))

(1 row)

Created by

. gisdb=# SELECT AsText(SymDifference(

'POLYGON((1 1,1 5,5 5,5 1,1 1))','LINESTRING(3 4,8 9)'));

Created by

. astext

------

Created by

. GEOMETRYCOLLECTION(LINESTRING(4 5,8 9),

POLYGON((1 1,1 5,4 5,5 5,5 1,1 1)))

(1 row)

Created by

. gisdb=# SELECT AsText(SymDifference(

Created by

. 'POLYGON((1 1,1 5,5 5,5 1,1 1))','MULTIPOINT(2 4,6 2)'));

astext

------

Created by

. GEOMETRYCOLLECTION(POINT(6 2),

POLYGON((1 1,1 5,5 5,5 1,1 1)))

Created by

. (1 row)

10. 5.10 Aggregate functions

Created by

. PostGIS’ aggregate functions create geometry from more geometries. These are used, where the usual SQL aggregate functions, for example in queries containing GROUP BY statement.

10.1. 5.10.1 ST_Union

Created by

. ST_Union

aggregate returns a geometry that represents the union of all input geometries.

Created by

. gisdb=# SELECT id, AsText(geom) FROM polygons; id | astext

----+------

Created by

. 1 | POLYGON((1 1,1 5,5 5,5 1,1 1))

2 | POLYGON((3 3,3 7,7 7,7 3,3 3))

3 | POLYGON((8 8,8 9,9 9,9 8,8 8))

Created by

. (3 rows)

gisdb=# SELECT AsText(ST_Union(geom)) FROM polygons;

Created by

. astext

------

Created by

. MULTIPOLYGON(((8 8,8 9,9 9,9 8,8 8)),

((1 1,1 5,3 5,3 7,7 7,7 3,5 3,5 1,1 1)))

Created by

. (1 row)

10.2. 5.10.2 Collect

Created by

. Collect

aggregate returns a Multi* geometry or a GeometryCollection, which contains the elements of the original geometries. In contrast to the

ST_Union

, this aggregate does not do any spatial operation between the geometries.

Created by

. gisdb=# SELECT AsText(Collect(geom)) FROM polygons;

astext

Created by

. ------

MULTIPOLYGON(((1 1,1 5,5 5,5 1,1 1)),((3 3,3 7,7 7,7 3,3 3)),

Created by

. ((8 8,8 9,9 9,9 8,8 8)))

(1 row)

Created by

. 10.3. 5.10.3 Poligonize

Poligonize

aggregate returns a GeometryCollection that contains the possible Polygons, which can be created from the input LineString and MultiLineString objects.

Created by

. gisdb=# SELECT id, AsText(geom) from linework;

id | astext

Created by

. ----+------

1 | LINESTRING(1 1,1 2)

Created by

. 2 | LINESTRING(1 2,2 2)

Created by

. (8 rows)

gisdb=# SELECT AsText(Polygonize(geom)) from linework;

Created by

. astext

------

GEOMETRYCOLLECTION(

Created by

. POLYGON((1 1,1 2,2 2,2 1,1 1)),

POLYGON((3 2,2 2,2 3,3 3,3 2)))

Created by

. (1 row)

11. 5.11 Operators

Created by

. PostGIS has some operators that return true or false values depending on the spatial relationship of the bounding boxes of the operands. These operators use spatial indexes of any operands.

Created by

. Table 1. Operators of PostGIS

returns TRUE if A's bounding box overlaps B's A && B

Created by

. A &< B returns TRUE if A's bounding box overlaps or is to the left of B's

A &<| returns TRUE if A's bounding box overlaps or is below B's B

Created by

. A &> B returns TRUE if A' bounding box overlaps or is to the right of B's

A << B returns TRUE if A's bounding box is strictly to the left of B's

Created by

. A <<| returns TRUE if A's bounding box is strictly below B's B

A = B returns TRUE if A's bounding box is the same as B's

Created by

. A >> B returns TRUE if A's bounding box is strictly to the right of B's

A @ B returns TRUE if A's bounding box is contained by B's

Created by

. B returns TRUE if A's bounding box overlaps or is above B's

A |>> returns TRUE if A's bounding box is strictly above B's B

A ~ B returns TRUE if A's bounding box contains B's

Created by

. A ~= B returns TRUE if A's bounding box is the same as B's

These operators are very useful for fast pre-screening in spatial queries.

Created by

. Bibliography

Created by

. PostgreSQL Global Development Group:

Created by

. 1996-2010.

Refractions Research Inc.:

Created by

. PostGIS 1.5.2 manual,

Created by

. Open Geospatial Consortium:

OpenGIS Implementation Standard for Geographic information - Simple feature access - Part 1: Common architecture, OGC 06-103r4,

Created by

. Open Geospatial Consortium:

OpenGIS Implementation Standard for Geographic information - Simple feature access - Part 2: SQL

Created by

. option, OGC 06-104r4,

Created by