File: quantize_coordinates.sql

package info (click to toggle)
postgis 3.5.2%2Bdfsg-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 70,052 kB
  • sloc: ansic: 162,204; sql: 93,950; xml: 53,121; cpp: 12,646; perl: 5,658; sh: 5,369; makefile: 3,434; python: 1,205; yacc: 447; lex: 151; pascal: 58
file content (24 lines) | stat: -rw-r--r-- 1,652 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
-- Returns NULL for NULL geometry
SELECT 't1', ST_QuantizeCoordinates(NULL::geometry, 6) IS NULL;
-- Throws error if no precision specified
SELECT 't2', ST_QuantizeCoordinates('POINT (3 7)', NULL);
-- Preserves SRID
SELECT 't3', ST_SRID(ST_QuantizeCoordinates('SRID=32611;POINT (3 7)', 5)) = 32611;
-- Carries X precision through to other dimensions if not specified
SELECT 't4', ST_X(ST_QuantizeCoordinates('POINT (1.2345678 1.2345678)', 3)) = ST_Y(ST_QuantizeCoordinates('POINT (1.2345678 1.2345678)', 3));
-- Or they can be specified independently
SELECT 't5', ST_X(ST_QuantizeCoordinates('POINT (1.2345678 1.2345678)', 7, 4)) != ST_Y(ST_QuantizeCoordinates('POINT (1.2345678 1.2345678)', 7, 4));
-- Significant digits are preserved
WITH input AS (SELECT ST_MakePoint(1.23456789012345, 0) AS geom)
SELECT 't6', bool_and(abs(ST_X(geom)-ST_X(ST_QuantizeCoordinates(geom, i))) < pow(10, -i))
FROM input, generate_series(1,15) AS i;
-- Test also negative precision
WITH input AS (SELECT ST_MakePoint(-12345.6789012345, 0) AS geom)
SELECT 't7', bool_and(abs(ST_X(geom)-ST_X(ST_QuantizeCoordinates(geom, i))) < pow(10, -i) AND abs(ST_X(geom)-ST_X(ST_QuantizeCoordinates(geom, i))) > pow(10, -i - 3))
FROM input, generate_series(-4,10) AS i;
-- Test NaN is preserved
SELECT 't8', ST_X(ST_QuantizeCoordinates('POINT EMPTY', 7, 4)) IS NULL;
-- Test that precision >= 18 fully preserves the number
SELECT 't9', ST_X(ST_QuantizeCoordinates('POINT (1.234567890123456 0)', 18)) = ST_X('POINT (1.234567890123456 0)');
-- Test very low precision
SELECT 't10', abs(ST_X(ST_QuantizeCoordinates('POINT (1234567890123456 0)', -18)) - 1234567890123456) <= pow(10, 18);