Mercurial > gemma
view schema/isrs_tests.sql @ 2519:2fa12229c946 octree-diff
Merged default into octree-diff branch.
author | Sascha L. Teichmann <sascha.teichmann@intevation.de> |
---|---|
date | Tue, 05 Mar 2019 17:43:43 +0100 |
parents | 00cac7890574 |
children | 73c8762cee60 |
line wrap: on
line source
-- This is Free Software under GNU Affero General Public License v >= 3.0 -- without warranty, see README.md and license for details. -- SPDX-License-Identifier: AGPL-3.0-or-later -- License-Filename: LICENSES/AGPL-3.0.txt -- Copyright (C) 2018 by via donau -- – Österreichische Wasserstraßen-Gesellschaft mbH -- Software engineering by Intevation GmbH -- Author(s): -- * Tom Gottfried <tom@intevation.de> -- -- pgTAP test script for ISRS location code types and functions -- SELECT results_eq($$ SELECT isrs_fromText('DEBON03901G007906548') $$, $$ SELECT CAST(('DE', 'BON', '03901', 'G0079', 6548) AS isrs) $$, 'Correct ISRS text input gives respective ISRS location code'); SELECT throws_ok($$ SELECT isrs_fromText('DEUXXX039000000005023') $$, 22023, NULL, 'ISRS text input needs to have correct length'); SELECT ok( 'DEBON03901G007906548' = isrs_asText(isrs_fromText('DEBON03901G007906548')) , 'isrs_asText() is the inverse of isrs_fromText()'); SELECT ok( ISRSrange_area(isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 1)::isrs), ST_SetSRID('POLYGON((0 1, 0 2, 1 2, 1 1, 0 1))'::geometry, 4326) ) IS NULL, 'ISRSrange_area returns NULL, if given area does not intersect with axis'); SELECT results_eq($$ SELECT every(ST_DWithin( ST_Boundary(ISRSrange_area( isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 1)::isrs), (SELECT ST_Collect(CAST(area AS geometry)) FROM waterway.waterway_area))), geom, 1)) FROM waterway.distance_marks_virtual WHERE location_code IN( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 1)::isrs) $$, $$ SELECT true $$, 'Resulting polygon almost ST_Touches points corresponding to stretch'); \set test_area 'POLYGON((-1 1, 2 1, 2 -1, -1 -1, -1 1))' SELECT ok( 2 = ST_NumGeometries( ISRSrange_area( isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 1)::isrs), ST_SetSRID(ST_Collect( ST_Translate(:'test_area', (ST_XMax(:'test_area'::geometry) - ST_XMin(:'test_area'::geometry))/2 + 0.1, 0), ST_Translate(:'test_area', -((ST_XMax(:'test_area'::geometry) - ST_XMin(:'test_area'::geometry))/2 + 0.1), 0)), 4326))), 'Two polygons intersecting the axis lead to two polygons in result'); SELECT ok( 1 = ST_NumGeometries( ISRSrange_area( isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 1)::isrs), ST_SetSRID(ST_Collect( ST_Translate(:'test_area', (ST_XMax(:'test_area'::geometry) - ST_XMin(:'test_area'::geometry))/3, 0), ST_Translate(:'test_area', -((ST_XMax(:'test_area'::geometry) - ST_XMin(:'test_area'::geometry))/3), 0)), 4326))), 'Self-intersecting multipolygon leads to one polygon in result'); SELECT results_eq($$ SELECT every(ST_DWithin( ST_Boundary(ISRSrange_area( isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 2)::isrs), (SELECT ST_Collect(CAST(area AS geometry)) FROM waterway.waterway_area))), geom, 1)) FROM waterway.distance_marks_virtual WHERE location_code IN( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 2)::isrs) $$, $$ SELECT true $$, 'Area generated from non-matching distance mark and non-contiguous axis');