lib-sql/functions/place_triggers.sql

   1 CREATE OR REPLACE FUNCTION place_insert()
   2   RETURNS TRIGGER
   3   AS $$
   4 DECLARE
   5   i INTEGER;
   6   existing RECORD;
   7   existingplacex RECORD;
   8   existingline RECORD;
   9   existinggeometry GEOMETRY;
  10   existingplace_id BIGINT;
  11   result BOOLEAN;
  12   partition INTEGER;
  13 BEGIN
  14
  15   {% if debug %}
  16     RAISE WARNING '-----------------------------------------------------------------------------------';
  17     RAISE WARNING 'place_insert: % % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type,st_area(NEW.geometry);
  18   {% endif %}
  19   -- filter wrong tupels
  20   IF ST_IsEmpty(NEW.geometry) OR NOT ST_IsValid(NEW.geometry) OR ST_X(ST_Centroid(NEW.geometry))::text in ('NaN','Infinity','-Infinity') OR ST_Y(ST_Centroid(NEW.geometry))::text in ('NaN','Infinity','-Infinity') THEN
  21     INSERT INTO import_polygon_error (osm_type, osm_id, class, type, name, country_code, updated, errormessage, prevgeometry, newgeometry)
  22       VALUES (NEW.osm_type, NEW.osm_id, NEW.class, NEW.type, NEW.name, NEW.address->'country', now(), ST_IsValidReason(NEW.geometry), null, NEW.geometry);
  23 --    RAISE WARNING 'Invalid Geometry: % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type;
  24     RETURN null;
  25   END IF;
  26
  27   -- decide, whether it is an osm interpolation line => insert intoosmline, or else just placex
  28   IF NEW.class='place' and NEW.type='houses' and NEW.osm_type='W' and ST_GeometryType(NEW.geometry) = 'ST_LineString' THEN
  29     -- Have we already done this place?
  30     select * from place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type INTO existing;
  31
  32     -- Get the existing place_id
  33     select * from location_property_osmline where osm_id = NEW.osm_id INTO existingline;
  34
  35     -- Handle a place changing type by removing the old data (this trigger is executed BEFORE INSERT of the NEW tupel)
  36     IF existing.osm_type IS NULL THEN
  37       DELETE FROM place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class;
  38     END IF;
  39
  40     DELETE from import_polygon_error where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
  41     DELETE from import_polygon_delete where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
  42
  43     -- update method for interpolation lines: delete all old interpolation lines with same osm_id (update on place) and insert the new one(s) (they can be split up, if they have > 2 nodes)
  44     IF existingline.osm_id IS NOT NULL THEN
  45       delete from location_property_osmline where osm_id = NEW.osm_id;
  46     END IF;
  47
  48     -- for interpolations invalidate all nodes on the line
  49     update placex p set indexed_status = 2
  50       from planet_osm_ways w
  51       where w.id = NEW.osm_id and p.osm_type = 'N' and p.osm_id = any(w.nodes);
  52
  53
  54     INSERT INTO location_property_osmline (osm_id, address, linegeo)
  55       VALUES (NEW.osm_id, NEW.address, NEW.geometry);
  56
  57
  58     IF existing.osm_type IS NULL THEN
  59       return NEW;
  60     END IF;
  61
  62     IF coalesce(existing.address, ''::hstore) != coalesce(NEW.address, ''::hstore)
  63        OR (coalesce(existing.extratags, ''::hstore) != coalesce(NEW.extratags, ''::hstore))
  64        OR existing.geometry::text != NEW.geometry::text
  65        THEN
  66
  67       update place set
  68         name = NEW.name,
  69         address = NEW.address,
  70         extratags = NEW.extratags,
  71         admin_level = NEW.admin_level,
  72         geometry = NEW.geometry
  73         where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type;
  74     END IF;
  75
  76     RETURN NULL;
  77
  78   ELSE -- insert to placex
  79
  80     -- Patch in additional country names
  81     IF NEW.admin_level = 2 AND NEW.type = 'administrative'
  82           AND NEW.address is not NULL AND NEW.address ? 'country' THEN
  83         SELECT name FROM country_name WHERE country_code = lower(NEW.address->'country') INTO existing;
  84         IF existing.name IS NOT NULL THEN
  85             NEW.name = existing.name || NEW.name;
  86         END IF;
  87     END IF;
  88
  89     -- Have we already done this place?
  90     select * from place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type INTO existing;
  91
  92     -- Get the existing place_id
  93     select * from placex where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type INTO existingplacex;
  94
  95     -- Pure postcodes are never queried from placex so we don't add them.
  96     -- location_postcodes is filled from the place table directly.
  97     IF NEW.class = 'place' AND NEW.type = 'postcode' THEN
  98       -- Remove old placex entry if the type changed to postcode.
  99       IF existingplacex.type IS NOT NULL AND existingplacex.type != 'postcode' THEN
 100         DELETE FROM placex where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
 101       END IF;
 102       RETURN NEW;
 103     END IF;
 104
 105     -- Handle a place changing type by removing the old data
 106     -- My generated 'place' types are causing havok because they overlap with real keys
 107     -- TODO: move them to their own special purpose key/class to avoid collisions
 108     IF existing.osm_type IS NULL THEN
 109       DELETE FROM place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class;
 110     END IF;
 111
 112     {% if debug %}RAISE WARNING 'Existing: %',existing.osm_id;{% endif %}
 113     {% if debug %}RAISE WARNING 'Existing PlaceX: %',existingplacex.place_id;{% endif %}
 114
 115     -- Log and discard
 116     IF existing.geometry is not null AND st_isvalid(existing.geometry)
 117       AND st_area(existing.geometry) > 0.02
 118       AND ST_GeometryType(NEW.geometry) in ('ST_Polygon','ST_MultiPolygon')
 119       AND st_area(NEW.geometry) < st_area(existing.geometry)*0.5
 120       THEN
 121       INSERT INTO import_polygon_error (osm_type, osm_id, class, type, name, country_code, updated, errormessage, prevgeometry, newgeometry)
 122         VALUES (NEW.osm_type, NEW.osm_id, NEW.class, NEW.type, NEW.name, NEW.address->'country', now(),
 123         'Area reduced from '||st_area(existing.geometry)||' to '||st_area(NEW.geometry), existing.geometry, NEW.geometry);
 124       RETURN null;
 125     END IF;
 126
 127     DELETE from import_polygon_error where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
 128     DELETE from import_polygon_delete where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
 129
 130     -- To paraphrase, if there isn't an existing item, OR if the admin level has changed
 131     IF existingplacex.osm_type IS NULL OR
 132         (existingplacex.class = 'boundary' AND
 133           ((coalesce(existingplacex.admin_level, 15) != coalesce(NEW.admin_level, 15) AND existingplacex.type = 'administrative') OR
 134           (existingplacex.type != NEW.type)))
 135     THEN
 136
 137       {% if config.get_bool('LIMIT_REINDEXING') %}
 138       IF existingplacex.osm_type IS NOT NULL THEN
 139         -- sanity check: ignore admin_level changes on places with too many active children
 140         -- or we end up reindexing entire countries because somebody accidentally deleted admin_level
 141         SELECT count(*) INTO i FROM
 142           (SELECT 'a' FROM placex, place_addressline
 143             WHERE address_place_id = existingplacex.place_id
 144                   and placex.place_id = place_addressline.place_id
 145                   and indexed_status = 0 and place_addressline.isaddress LIMIT 100001) sub;
 146         IF i > 100000 THEN
 147           RETURN null;
 148         END IF;
 149       END IF;
 150       {% endif %}
 151
 152       IF existing.osm_type IS NOT NULL THEN
 153         -- pathological case caused by the triggerless copy into place during initial import
 154         -- force delete even for large areas, it will be reinserted later
 155         UPDATE place set geometry = ST_SetSRID(ST_Point(0,0), 4326) where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type;
 156         DELETE from place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type;
 157       END IF;
 158
 159       -- No - process it as a new insertion (hopefully of low rank or it will be slow)
 160       insert into placex (osm_type, osm_id, class, type, name,
 161                           admin_level, address, extratags, geometry)
 162         values (NEW.osm_type, NEW.osm_id, NEW.class, NEW.type, NEW.name,
 163                 NEW.admin_level, NEW.address, NEW.extratags, NEW.geometry);
 164
 165       {% if debug %}RAISE WARNING 'insert done % % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type,NEW.name;{% endif %}
 166
 167       RETURN NEW;
 168     END IF;
 169
 170     -- Special case for polygon shape changes because they tend to be large and we can be a bit clever about how we handle them
 171     IF existing.geometry::text != NEW.geometry::text
 172        AND ST_GeometryType(existing.geometry) in ('ST_Polygon','ST_MultiPolygon')
 173        AND ST_GeometryType(NEW.geometry) in ('ST_Polygon','ST_MultiPolygon')
 174        THEN
 175
 176       -- Get the version of the geometry actually used (in placex table)
 177       select geometry from placex where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type into existinggeometry;
 178
 179       -- Performance limit
 180       IF st_area(NEW.geometry) < 0.000000001 AND st_area(existinggeometry) < 1 THEN
 181
 182         -- re-index points that have moved in / out of the polygon, could be done as a single query but postgres gets the index usage wrong
 183         update placex set indexed_status = 2 where indexed_status = 0
 184             AND ST_Intersects(NEW.geometry, placex.geometry)
 185             AND NOT ST_Intersects(existinggeometry, placex.geometry)
 186             AND rank_search > existingplacex.rank_search AND (rank_search < 28 or name is not null);
 187
 188         update placex set indexed_status = 2 where indexed_status = 0
 189             AND ST_Intersects(existinggeometry, placex.geometry)
 190             AND NOT ST_Intersects(NEW.geometry, placex.geometry)
 191             AND rank_search > existingplacex.rank_search AND (rank_search < 28 or name is not null);
 192
 193       END IF;
 194
 195     END IF;
 196
 197
 198     IF coalesce(existing.name::text, '') != coalesce(NEW.name::text, '')
 199        OR coalesce(existing.extratags::text, '') != coalesce(NEW.extratags::text, '')
 200        OR coalesce(existing.address, ''::hstore) != coalesce(NEW.address, ''::hstore)
 201        OR coalesce(existing.admin_level, 15) != coalesce(NEW.admin_level, 15)
 202        OR existing.geometry::text != NEW.geometry::text
 203        THEN
 204
 205       update place set
 206         name = NEW.name,
 207         address = NEW.address,
 208         extratags = NEW.extratags,
 209         admin_level = NEW.admin_level,
 210         geometry = NEW.geometry
 211         where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type;
 212
 213
 214       IF NEW.class = 'boundary' AND NEW.type = 'postal_code' THEN
 215           IF NEW.address is NULL OR NOT NEW.address ? 'postcode' THEN
 216               -- postcode was deleted, no longer retain in placex
 217               DELETE FROM placex where place_id = existingplacex.place_id;
 218               RETURN NULL;
 219           END IF;
 220
 221           NEW.name := hstore('ref', NEW.address->'postcode');
 222       END IF;
 223
 224       IF NEW.class in ('boundary')
 225          AND ST_GeometryType(NEW.geometry) not in ('ST_Polygon','ST_MultiPolygon') THEN
 226           DELETE FROM placex where place_id = existingplacex.place_id;
 227           RETURN NULL;
 228       END IF;
 229
 230       update placex set
 231         name = NEW.name,
 232         address = NEW.address,
 233         parent_place_id = null,
 234         extratags = NEW.extratags,
 235         admin_level = NEW.admin_level,
 236         indexed_status = 2,
 237         geometry = NEW.geometry
 238         where place_id = existingplacex.place_id;
 239       -- if a node(=>house), which is part of a interpolation line, changes (e.g. the street attribute) => mark this line for reparenting
 240       -- (already here, because interpolation lines are reindexed before nodes, so in the second call it would be too late)
 241       IF NEW.osm_type='N'
 242          and (coalesce(existing.address, ''::hstore) != coalesce(NEW.address, ''::hstore)
 243              or existing.geometry::text != NEW.geometry::text)
 244       THEN
 245           result:= osmline_reinsert(NEW.osm_id, NEW.geometry);
 246       END IF;
 247
 248       -- linked places should get potential new naming and addresses
 249       IF existingplacex.linked_place_id is not NULL THEN
 250         update placex x set
 251           name = p.name,
 252           extratags = p.extratags,
 253           indexed_status = 2
 254         from place p
 255         where x.place_id = existingplacex.linked_place_id
 256               and x.indexed_status = 0
 257               and x.osm_type = p.osm_type
 258               and x.osm_id = p.osm_id
 259               and x.class = p.class;
 260       END IF;
 261
 262     END IF;
 263
 264     -- Abort the add (we modified the existing place instead)
 265     RETURN NULL;
 266   END IF;
 267
 268 END;
 269 $$ LANGUAGE plpgsql;
 270
 271
 272 CREATE OR REPLACE FUNCTION place_delete()
 273   RETURNS TRIGGER
 274   AS $$
 275 DECLARE
 276   has_rank BOOLEAN;
 277 BEGIN
 278
 279   {% if debug %}RAISE WARNING 'delete: % % % %',OLD.osm_type,OLD.osm_id,OLD.class,OLD.type;{% endif %}
 280
 281   -- deleting large polygons can have a massive effect on the system - require manual intervention to let them through
 282   IF st_area(OLD.geometry) > 2 and st_isvalid(OLD.geometry) THEN
 283     SELECT bool_or(not (rank_address = 0 or rank_address > 25)) as ranked FROM placex WHERE osm_type = OLD.osm_type and osm_id = OLD.osm_id and class = OLD.class and type = OLD.type INTO has_rank;
 284     IF has_rank THEN
 285       insert into import_polygon_delete (osm_type, osm_id, class, type) values (OLD.osm_type,OLD.osm_id,OLD.class,OLD.type);
 286       RETURN NULL;
 287     END IF;
 288   END IF;
 289
 290   -- mark for delete
 291   UPDATE placex set indexed_status = 100 where osm_type = OLD.osm_type and osm_id = OLD.osm_id and class = OLD.class and type = OLD.type;
 292
 293   -- interpolations are special
 294   IF OLD.osm_type='W' and OLD.class = 'place' and OLD.type = 'houses' THEN
 295     UPDATE location_property_osmline set indexed_status = 100 where osm_id = OLD.osm_id; -- osm_id = wayid (=old.osm_id)
 296   END IF;
 297
 298   RETURN OLD;
 299 END;
 300 $$
 301 LANGUAGE plpgsql;
 302