Mercurial > gemma
view schema/isrs_tests.sql @ 2535:73c8762cee60
Stretch area generation: Throw errors instead of returning NULL
It is now an error if not both distance marks for the given stretch
can be found or no contiguous axis can be constructed between them.
author | Tom Gottfried <tom@intevation.de> |
---|---|
date | Thu, 07 Mar 2019 12:20:00 +0100 |
parents | 00cac7890574 |
children | dc4fae4bdb8f |
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 throws_ok($$ SELECT ISRSrange_points(isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 99999)::isrs)) $$, 'P0002', NULL, 'ISRSrange_points fails if not both distance marks can be found'); SELECT throws_ok($$ SELECT ISRSrange_axis(isrsrange( ('AT', 'XXX', '00001', '00000', 0)::isrs, ('AT', 'XXX', '00001', '00000', 2)::isrs), 0) $$, 'P0002', NULL, 'ISRSrange_axis fails if no contiguous axis can be constructed'); 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');