trigger-example.sql

   1 -- PHOROS -- Photogrammetric Road Survey
   2 -- Copyright (C) 2011 Bert Burgemeister
   3 --
   4 -- This program is free software; you can redistribute it and/or modify
   5 -- it under the terms of the GNU General Public License as published by
   6 -- the Free Software Foundation; either version 2 of the License, or
   7 -- (at your option) any later version.
   8 --
   9 -- This program is distributed in the hope that it will be useful,
  10 -- but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 -- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12 -- GNU General Public License for more details.
  13 --
  14 -- You should have received a copy of the GNU General Public License along
  15 -- with this program; if not, write to the Free Software Foundation, Inc.,
  16 -- 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
  17 ----------------------------------------------------------------------
  18
  19
  20 -- This is a working example for a PL/pgSQL trigger function body.
  21 -- Once deployed, it is fired on changes to the user point table.  It
  22 -- can be deployed with your presentation project like this:
  23 --
  24 -- ./phoros --redefine-trigger-function=<your-presentation-project> --plpgsql-body=<name-of-this-file> ...
  25 --
  26 -- The code threads points from the user point table, grouped by
  27 -- column description, to linestrings which are stored into the user
  28 -- line table.  The strings ~0@*~A and ~1@*~A are placeholders for the
  29 -- names of user point table and user line table, respectively.
  30
  31 DECLARE
  32   polyline GEOMETRY;
  33   new_point point_bag%ROWTYPE;
  34   tried_point point_bag%ROWTYPE;
  35   previous_point point_bag%ROWTYPE;
  36   starting_point GEOMETRY;
  37   point_group RECORD;
  38   old_description text;
  39   new_description text;
  40   -- Set max_bend = 180 in order not to discard any points, and to get
  41   -- crazy zigzag lines.
  42   max_bend DOUBLE PRECISION DEFAULT 91; -- degrees
  43
  44 BEGIN
  45
  46   -- Muffle warnings about implicitly created stuff:
  47   SET client_min_messages TO ERROR;
  48
  49   IF TG_OP = 'UPDATE'
  50   THEN
  51     old_description := OLD.description;
  52     new_description := NEW.description;
  53   END IF;
  54
  55   IF TG_OP = 'INSERT'
  56   THEN
  57     new_description := NEW.description;
  58   END IF;
  59
  60   IF TG_OP = 'DELETE'
  61   THEN
  62     old_description := OLD.description;
  63   END IF;
  64
  65   DROP TABLE IF EXISTS point_bag_table;
  66
  67   CREATE TEMPORARY TABLE point_bag_table
  68     (id SERIAL primary key, coordinates GEOMETRY, description TEXT)
  69     ON COMMIT DROP;
  70
  71   <<thread_one_line>>
  72   FOR point_group
  73   IN SELECT description FROM ~0@*~A
  74   WHERE description = old_description OR description = new_description
  75   GROUP BY description
  76   LOOP
  77     INSERT INTO point_bag_table (coordinates, description)
  78       SELECT coordinates, description
  79         FROM ~0@*~A
  80         WHERE attribute = 'polyline'
  81               AND
  82               description = point_group.description;
  83
  84     starting_point :=
  85       (SELECT st_centroid(st_collect(coordinates))
  86          FROM point_bag_table);
  87
  88     previous_point :=
  89       (SELECT ROW(id, coordinates)
  90          FROM point_bag_table
  91          ORDER BY st_distance(point_bag_table.coordinates, starting_point)
  92          LIMIT 1);
  93
  94     DELETE FROM point_bag_table WHERE id = previous_point.id;
  95
  96     new_point :=
  97       (SELECT ROW(id, coordinates)
  98          FROM point_bag_table
  99          ORDER BY st_distance(point_bag_table.coordinates,
 100                               previous_point.coordinates)
 101          LIMIT 1);
 102
 103     polyline := st_makeline(previous_point.coordinates,
 104                         new_point.coordinates);
 105
 106     DELETE FROM point_bag_table WHERE id = new_point.id;
 107
 108     <<add_or_discard_point>>
 109     LOOP
 110       previous_point.coordinates := st_pointn(polyline,1);
 111
 112       new_point :=
 113         (SELECT ROW(id, coordinates)
 114            FROM point_bag_table
 115            ORDER BY st_distance(coordinates, previous_point.coordinates)
 116            LIMIT 1);
 117
 118       EXIT WHEN new_point IS NULL;
 119
 120       IF bendedness(st_pointn(polyline, 2),
 121                     st_pointn(polyline, 1),
 122                     new_point.coordinates)
 123          < bendedness(st_pointn(polyline, st_npoints(polyline) - 1),
 124                       st_pointn(polyline, st_npoints(polyline)),
 125                       new_point.coordinates)
 126          AND
 127          bendedness(st_pointn(polyline, 2), st_pointn(polyline, 1),
 128                      new_point.coordinates)
 129          < radians(max_bend)
 130       THEN
 131           polyline := st_addpoint(polyline, new_point.coordinates, 0);
 132           DELETE FROM point_bag_table WHERE id = new_point.id;
 133       END IF;
 134
 135       polyline := st_reverse(polyline);
 136
 137       DELETE FROM point_bag_table WHERE id = tried_point.id;
 138
 139       tried_point := new_point;
 140     END LOOP add_or_discard_point;
 141
 142     -- RAISE NOTICE '%', st_astext(polyline);
 143
 144     DELETE FROM ~1@*~A
 145            WHERE description = point_group.description;
 146
 147     IF polyline IS NOT NULL
 148     THEN
 149       INSERT INTO ~1@*~A
 150              (description, line)
 151              VALUES (point_group.description, polyline);
 152     END IF;
 153
 154   END LOOP thread_one_line;
 155 END;