File: rt_convexhull.sql

package info (click to toggle)
postgis 2.5.1%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 75,792 kB
  • sloc: ansic: 139,314; sql: 136,281; xml: 48,954; sh: 4,906; perl: 4,509; makefile: 2,897; python: 1,198; yacc: 441; cpp: 305; lex: 132
file content (136 lines) | stat: -rw-r--r-- 3,397 bytes parent folder | download | duplicates (11)
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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
DROP TABLE IF EXISTS test_raster_convexhull;

CREATE TABLE test_raster_convexhull AS
	SELECT ST_AddBand(ST_AddBand(ST_MakeEmptyRaster(9, 9, 0, 0, 1, -1, 0, 0, 0), 1, '8BUI', 0, 0), 2, '8BUI', 1, 0) AS rast;

SELECT
	ST_AsText(ST_ConvexHull(rast)),
	ST_AsText(ST_MinConvexHull(rast)),
	ST_AsText(ST_MinConvexHull(rast, 1)),
	ST_AsText(ST_MinConvexHull(rast, 2))
FROM test_raster_convexhull;

UPDATE test_raster_convexhull SET
	rast = ST_SetValues(
		rast, 1, 1, 1,
		ARRAY[
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 1, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0]
		]::double precision[][]
	);

SELECT
	ST_AsText(ST_MinConvexHull(rast, 1))
FROM test_raster_convexhull;

UPDATE test_raster_convexhull SET
	rast = ST_SetValues(
		rast, 1, 1, 1,
		ARRAY[
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 1, 0, 0, 0, 0, 0],
			[0, 0, 0, 1, 1, 0, 0, 0, 0],
			[0, 0, 0, 1, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0]
		]::double precision[][]
	);

SELECT
	ST_AsText(ST_MinConvexHull(rast, 1))
FROM test_raster_convexhull;

UPDATE test_raster_convexhull SET
	rast = ST_SetValues(
		rast, 1, 1, 1,
		ARRAY[
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 1, 0, 1, 0, 0, 0],
			[0, 0, 0, 1, 1, 0, 0, 0, 0],
			[0, 0, 0, 1, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0]
		]::double precision[][]
	);

SELECT
	ST_AsText(ST_MinConvexHull(rast, 1))
FROM test_raster_convexhull;

UPDATE test_raster_convexhull SET
	rast = ST_SetValues(
		rast, 1, 1, 1,
		ARRAY[
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 1, 0, 0, 0, 0, 1],
			[0, 0, 0, 1, 1, 0, 0, 0, 0],
			[0, 0, 0, 1, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0]
		]::double precision[][]
	);

SELECT
	ST_AsText(ST_MinConvexHull(rast, 1))
FROM test_raster_convexhull;

UPDATE test_raster_convexhull SET
	rast = ST_SetValues(
		rast, 2, 1, 1,
		ARRAY[
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[1, 0, 0, 0, 0, 1, 0, 0, 0],
			[0, 0, 0, 0, 1, 1, 0, 0, 0],
			[0, 0, 0, 0, 0, 1, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0]
		]::double precision[][]
	);

SELECT
	ST_AsText(ST_MinConvexHull(rast, 2)),
	ST_AsText(ST_MinConvexHull(rast))
FROM test_raster_convexhull;

UPDATE test_raster_convexhull SET
	rast = ST_SetValues(
		rast, 2, 1, 1,
		ARRAY[
			[0, 0, 0, 0, 0, 0, 0, 0, 1],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[1, 0, 0, 0, 0, 1, 0, 0, 0],
			[0, 0, 0, 0, 1, 1, 0, 0, 0],
			[0, 0, 0, 0, 0, 1, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 0, 0, 0, 0, 0, 0, 0],
			[0, 0, 1, 0, 0, 0, 0, 0, 0]
		]::double precision[][]
	);

SELECT
	ST_AsText(ST_MinConvexHull(rast, 1)),
	ST_AsText(ST_MinConvexHull(rast, 2)),
	ST_AsText(ST_MinConvexHull(rast))
FROM test_raster_convexhull;

DROP TABLE IF EXISTS test_raster_convexhull;