File: edge_cases.pg

package info (click to toggle)
pgrouting 4.0.1-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 17,332 kB
  • sloc: cpp: 21,315; sql: 10,419; ansic: 9,795; perl: 1,142; sh: 919; javascript: 314; xml: 182; makefile: 29
file content (89 lines) | stat: -rw-r--r-- 2,297 bytes parent folder | download
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
/* :file: This file is part of the pgRouting project.
:copyright: Copyright (c) 2025-2026 pgRouting developers
:license: Creative Commons Attribution-Share Alike 3.0 https://creativecommons.org/licenses/by-sa/3.0 */

BEGIN;

SELECT CASE WHEN min_version('3.8.0') THEN plan(6) ELSE plan(1) END;

CREATE TABLE grid(
      id serial NOT NULL primary key,
      dir text,
      the_geom geometry
);

-- creating the grid geometry
INSERT INTO grid(dir, the_geom)
SELECT
  CASE WHEN s1 % 2 = 0 THEN 'FT' ELSE 'TF' END,
  ST_AsTexT(ST_MakeLine(ST_MakePoint(1,s1), ST_MakePoint(6,s1)))
FROM (SELECT generate_series(1,6) AS s1) AS foo

UNION ALL

SELECT
  CASE WHEN s1 % 2 = 0 THEN 'FT' ELSE 'TF' END,
  ST_AsTexT(ST_MakeLine(ST_MakePoint(s1,1), ST_MakePoint(s1,6)))
  FROM (SELECT generate_series(1,6) AS s1) AS foo;


CREATE OR REPLACE FUNCTION issue_test()
RETURNS SETOF TEXT AS
$BODY$
BEGIN

  IF NOT min_version('3.8.0') THEN
    RETURN QUERY SELECT skip(1, 'pgr_separateTouching new on 3.8.0');
    RETURN;
  END IF;

  RETURN QUERY
  SELECT is((SELECT count(*)::INTEGER FROM grid), 12, 'There are 12 original edges');

  CREATE TABLE grid_r AS
  SELECT *, NULL::INTEGER source, NULL::INTEGER target, NULL::TEXT dir
  FROM pgr_separateTouching('SELECT id, the_geom AS geom FROM grid', 0.001);

  RETURN QUERY
  SELECT has_total_edges_vertices('grid_r', 20, 20);

  RETURN QUERY
  SELECT set_eq(
    'SELECT id, count(*) FROM grid_r GROUP BY id ORDER BY id',
    'VALUES (1,5),(6,5),(7,5),(12,5)',
    'For each original edge there are 5 subedges');

  -- set the source information
  UPDATE grid_r AS e
  SET source = rv.id
  FROM tmp_rv rv
  WHERE ST_StartPoint(e.geom) = rv.geom;

  -- set the target information
  UPDATE grid_r AS e
  SET target = rv.id
  FROM tmp_rv rv
  WHERE ST_EndPoint(e.geom) = rv.geom;

  -- se the dir
  UPDATE grid_r AS e
  SET dir = grid.dir
  FROM grid
  WHERE grid.id = e.id;

  RETURN QUERY
  SELECT is(
    (SELECT DISTINCT count(*)::INTEGER FROM grid GROUP BY dir),
    6, '6 is the count of different TF and FT in grid');
  RETURN QUERY
  SELECT is(
    (SELECT DISTINCT count(*)::INTEGER FROM grid_r GROUP BY dir),
    10, '10 is the count of different TF and FT in grid_r (4*5)');

END;
$BODY$ LANGUAGE plpgsql;

SELECT issue_test();

SELECT finish();
ROLLBACK;