File: union.sql

package info (click to toggle)
postgis 3.5.3%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 69,528 kB
  • sloc: ansic: 162,229; sql: 93,970; xml: 53,139; cpp: 12,646; perl: 5,658; sh: 5,369; makefile: 3,435; python: 1,205; yacc: 447; lex: 151; pascal: 58
file content (90 lines) | stat: -rw-r--r-- 2,654 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
CREATE TABLE geoms (geom geometry);

-- Create 20x20 square polygons
WITH coords AS (SELECT * FROM generate_series(0, 19) AS x, generate_series(0, 19) AS y)
INSERT INTO geoms
SELECT ST_Square(1.0, x, y, ST_MakePoint(0, 0))
FROM coords;

WITH u AS (SELECT ST_Union(geom) AS g FROM geoms)
SELECT ST_Area(g), ST_XMin((g)), ST_YMin(g), ST_XMax(g), ST_YMax(g) from u;

-- NOTE: using gridSize > 0 requires GEOS-3.9 or higher
WITH u AS (SELECT ST_Union(geom, -1.0) AS g FROM geoms)
SELECT ST_Area(g), ST_XMin((g)), sT_YMin(g), ST_XMax(g), ST_YMax(g) from u;

TRUNCATE TABLE geoms;

-- Empty table
SELECT ST_Union(geom) IS NULL FROM geoms;

-- NULLs
INSERT INTO geoms (geom) VALUES (NULL), (NULL), (NULL), (NULL), (NULL);
SELECT ST_Union(geom) IS NULL FROM geoms;

-- Add empty point
INSERT INTO geoms SELECT ST_GeomFromText('POINT EMPTY') AS geom;
SELECT ST_AsText(ST_Union(geom)) FROM geoms;

-- Add empty polygon
INSERT INTO geoms SELECT ST_GeomFromText('POLYGON EMPTY') AS geom;
SELECT ST_AsText(ST_Union(geom)) FROM geoms;

/** PG16 force_parallel_mode was renamed to debug_parallel_query, thus the need for this conditional **/
CREATE OR REPLACE PROCEDURE p_force_parellel_mode(param_state text) language plpgsql AS
$$
BEGIN
	IF (_postgis_pgsql_version())::integer < 160 THEN
		IF param_state = 'on' THEN
			SET force_parallel_mode=on;
		ELSE
			SET force_parallel_mode=off;
		END IF;
	ELSE
		IF param_state = 'on' THEN
			SET debug_parallel_query=on;
		ELSE
			SET debug_parallel_query=off;
		END IF;
	END IF;
END;
$$;

-- Test in parallel mode

TRUNCATE TABLE geoms;

CALL p_force_parellel_mode('on'::text);
SET parallel_setup_cost = 0;
SET parallel_tuple_cost = 0;
SET min_parallel_table_scan_size = 0;
SET max_parallel_workers_per_gather = 4;

WITH coords AS (SELECT * FROM generate_series(0, 19) AS x, generate_series(0, 19) AS y)
INSERT INTO geoms
SELECT ST_Square(1.0, x, y, ST_MakePoint(0, 0))
FROM coords;

WITH u AS (SELECT ST_Union(geom) AS g FROM geoms)
SELECT ST_Area(g), ST_XMin((g)), ST_YMin(g), ST_XMax(g), ST_YMax(g) from u;

WITH u AS (SELECT ST_Union(geom, -1.0) AS g FROM geoms)
SELECT ST_Area(g), ST_XMin((g)), sT_YMin(g), ST_XMax(g), ST_YMax(g) from u;

TRUNCATE TABLE geoms;

SELECT ST_Union(geom) IS NULL FROM geoms;

INSERT INTO geoms (geom) VALUES (NULL), (NULL), (NULL), (NULL), (NULL);
SELECT ST_Union(geom) IS NULL FROM geoms;

INSERT INTO geoms SELECT ST_GeomFromText('POINT EMPTY') AS geom;
SELECT ST_AsText(ST_Union(geom)) FROM geoms;

INSERT INTO geoms SELECT ST_GeomFromText('POLYGON EMPTY') AS geom;
SELECT ST_AsText(ST_Union(geom)) FROM geoms;


DROP TABLE geoms;

DROP PROCEDURE IF EXISTS p_force_parellel_mode(text);