Mercurial > gemma
view schema/isrs_tests.sql @ 2431:48495bd3081d
Construct stretch area between distance marks also from non-contiguous axis
Before, both distance marks were required to be next to the same linestring
of the waterway axis.
author | Tom Gottfried <tom@intevation.de> |
---|---|
date | Thu, 28 Feb 2019 19:44:38 +0100 |
parents | 8481e6266691 |
children | 00cac7890574 |
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 ok( ST_DWithin( (SELECT geom FROM waterway.distance_marks_virtual WHERE location_code = ('AT', 'XXX', '00001', '00000', 0)::isrs), 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))), 1) AND ST_DWithin( (SELECT geom FROM waterway.distance_marks_virtual WHERE location_code = ('AT', 'XXX', '00001', '00000', 1)::isrs), 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))), 1), '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 ok( 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)) IS NOT NULL, 'Area generated from non-matching distance mark and non-contiguous axis');