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.
  99       DELETE FROM placex where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
 100       RETURN NEW;
 101     END IF;
 102
 103     -- Handle a place changing type by removing the old data
 104     -- My generated 'place' types are causing havok because they overlap with real keys
 105     -- TODO: move them to their own special purpose key/class to avoid collisions
 106     IF existing.osm_type IS NULL THEN
 107       DELETE FROM place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class;
 108     END IF;
 109
 110     {% if debug %}RAISE WARNING 'Existing: %',existing.osm_id;{% endif %}
 111     {% if debug %}RAISE WARNING 'Existing PlaceX: %',existingplacex.place_id;{% endif %}
 112
 113     -- Log and discard
 114     IF existing.geometry is not null AND st_isvalid(existing.geometry)
 115       AND st_area(existing.geometry) > 0.02
 116       AND ST_GeometryType(NEW.geometry) in ('ST_Polygon','ST_MultiPolygon')
 117       AND st_area(NEW.geometry) < st_area(existing.geometry)*0.5
 118       THEN
 119       INSERT INTO import_polygon_error (osm_type, osm_id, class, type, name, country_code, updated, errormessage, prevgeometry, newgeometry)
 120         VALUES (NEW.osm_type, NEW.osm_id, NEW.class, NEW.type, NEW.name, NEW.address->'country', now(),
 121         'Area reduced from '||st_area(existing.geometry)||' to '||st_area(NEW.geometry), existing.geometry, NEW.geometry);
 122       RETURN null;
 123     END IF;
 124
 125     DELETE from import_polygon_error where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
 126     DELETE from import_polygon_delete where osm_type = NEW.osm_type and osm_id = NEW.osm_id;
 127
 128     -- To paraphrase, if there isn't an existing item, OR if the admin level has changed
 129     IF existingplacex.osm_type IS NULL OR
 130         (existingplacex.class = 'boundary' AND
 131           ((coalesce(existingplacex.admin_level, 15) != coalesce(NEW.admin_level, 15) AND existingplacex.type = 'administrative') OR
 132           (existingplacex.type != NEW.type)))
 133     THEN
 134
 135       {% if config.get_bool('LIMIT_REINDEXING') %}
 136       IF existingplacex.osm_type IS NOT NULL THEN
 137         -- sanity check: ignore admin_level changes on places with too many active children
 138         -- or we end up reindexing entire countries because somebody accidentally deleted admin_level
 139         SELECT count(*) INTO i FROM
 140           (SELECT 'a' FROM placex, place_addressline
 141             WHERE address_place_id = existingplacex.place_id
 142                   and placex.place_id = place_addressline.place_id
 143                   and indexed_status = 0 and place_addressline.isaddress LIMIT 100001) sub;
 144         IF i > 100000 THEN
 145           RETURN null;
 146         END IF;
 147       END IF;
 148       {% endif %}
 149
 150       IF existing.osm_type IS NOT NULL THEN
 151         -- pathological case caused by the triggerless copy into place during initial import
 152         -- force delete even for large areas, it will be reinserted later
 153         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;
 154         DELETE from place where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type;
 155       END IF;
 156
 157       -- No - process it as a new insertion (hopefully of low rank or it will be slow)
 158       insert into placex (osm_type, osm_id, class, type, name,
 159                           admin_level, address, extratags, geometry)
 160         values (NEW.osm_type, NEW.osm_id, NEW.class, NEW.type, NEW.name,
 161                 NEW.admin_level, NEW.address, NEW.extratags, NEW.geometry);
 162
 163       {% if debug %}RAISE WARNING 'insert done % % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type,NEW.name;{% endif %}
 164
 165       RETURN NEW;
 166     END IF;
 167
 168     -- Special case for polygon shape changes because they tend to be large and we can be a bit clever about how we handle them
 169     IF existing.geometry::text != NEW.geometry::text
 170        AND ST_GeometryType(existing.geometry) in ('ST_Polygon','ST_MultiPolygon')
 171        AND ST_GeometryType(NEW.geometry) in ('ST_Polygon','ST_MultiPolygon')
 172        THEN
 173
 174       -- Get the version of the geometry actually used (in placex table)
 175       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;
 176
 177       -- Performance limit
 178       IF st_area(NEW.geometry) < 0.000000001 AND st_area(existinggeometry) < 1 THEN
 179
 180         -- 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
 181         update placex set indexed_status = 2 where indexed_status = 0
 182             AND ST_Intersects(NEW.geometry, placex.geometry)
 183             AND NOT ST_Intersects(existinggeometry, placex.geometry)
 184             AND rank_search > existingplacex.rank_search AND (rank_search < 28 or name is not null);
 185
 186         update placex set indexed_status = 2 where indexed_status = 0
 187             AND ST_Intersects(existinggeometry, placex.geometry)
 188             AND NOT ST_Intersects(NEW.geometry, placex.geometry)
 189             AND rank_search > existingplacex.rank_search AND (rank_search < 28 or name is not null);
 190
 191       END IF;
 192
 193     END IF;
 194
 195
 196     IF coalesce(existing.name::text, '') != coalesce(NEW.name::text, '')
 197        OR coalesce(existing.extratags::text, '') != coalesce(NEW.extratags::text, '')
 198        OR coalesce(existing.address, ''::hstore) != coalesce(NEW.address, ''::hstore)
 199        OR coalesce(existing.admin_level, 15) != coalesce(NEW.admin_level, 15)
 200        OR existing.geometry::text != NEW.geometry::text
 201        THEN
 202
 203       update place set
 204         name = NEW.name,
 205         address = NEW.address,
 206         extratags = NEW.extratags,
 207         admin_level = NEW.admin_level,
 208         geometry = NEW.geometry
 209         where osm_type = NEW.osm_type and osm_id = NEW.osm_id and class = NEW.class and type = NEW.type;
 210
 211
 212       IF NEW.class = 'boundary' AND NEW.type = 'postal_code' THEN
 213           IF NEW.address is NULL OR NOT NEW.address ? 'postcode' THEN
 214               -- postcode was deleted, no longer retain in placex
 215               DELETE FROM placex where place_id = existingplacex.place_id;
 216               RETURN NULL;
 217           END IF;
 218
 219           NEW.name := hstore('ref', NEW.address->'postcode');
 220       END IF;
 221
 222       IF NEW.class in ('boundary')
 223          AND ST_GeometryType(NEW.geometry) not in ('ST_Polygon','ST_MultiPolygon') THEN
 224           DELETE FROM placex where place_id = existingplacex.place_id;
 225           RETURN NULL;
 226       END IF;
 227
 228       update placex set
 229         name = NEW.name,
 230         address = NEW.address,
 231         parent_place_id = null,
 232         extratags = NEW.extratags,
 233         admin_level = NEW.admin_level,
 234         indexed_status = 2,
 235         geometry = NEW.geometry
 236         where place_id = existingplacex.place_id;
 237       -- if a node(=>house), which is part of a interpolation line, changes (e.g. the street attribute) => mark this line for reparenting
 238       -- (already here, because interpolation lines are reindexed before nodes, so in the second call it would be too late)
 239       IF NEW.osm_type='N'
 240          and (coalesce(existing.address, ''::hstore) != coalesce(NEW.address, ''::hstore)
 241              or existing.geometry::text != NEW.geometry::text)
 242       THEN
 243           result:= osmline_reinsert(NEW.osm_id, NEW.geometry);
 244       END IF;
 245
 246       -- linked places should get potential new naming and addresses
 247       IF existingplacex.linked_place_id is not NULL THEN
 248         update placex x set
 249           name = p.name,
 250           extratags = p.extratags,
 251           indexed_status = 2
 252         from place p
 253         where x.place_id = existingplacex.linked_place_id
 254               and x.indexed_status = 0
 255               and x.osm_type = p.osm_type
 256               and x.osm_id = p.osm_id
 257               and x.class = p.class;
 258       END IF;
 259
 260     END IF;
 261
 262     -- Abort the add (we modified the existing place instead)
 263     RETURN NULL;
 264   END IF;
 265
 266 END;
 267 $$ LANGUAGE plpgsql;
 268
 269
 270 CREATE OR REPLACE FUNCTION place_delete()
 271   RETURNS TRIGGER
 272   AS $$
 273 DECLARE
 274   has_rank BOOLEAN;
 275 BEGIN
 276
 277   {% if debug %}RAISE WARNING 'delete: % % % %',OLD.osm_type,OLD.osm_id,OLD.class,OLD.type;{% endif %}
 278
 279   -- deleting large polygons can have a massive effect on the system - require manual intervention to let them through
 280   IF st_area(OLD.geometry) > 2 and st_isvalid(OLD.geometry) THEN
 281     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;
 282     IF has_rank THEN
 283       insert into import_polygon_delete (osm_type, osm_id, class, type) values (OLD.osm_type,OLD.osm_id,OLD.class,OLD.type);
 284       RETURN NULL;
 285     END IF;
 286   END IF;
 287
 288   -- mark for delete
 289   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;
 290
 291   -- interpolations are special
 292   IF OLD.osm_type='W' and OLD.class = 'place' and OLD.type = 'houses' THEN
 293     UPDATE location_property_osmline set indexed_status = 100 where osm_id = OLD.osm_id; -- osm_id = wayid (=old.osm_id)
 294   END IF;
 295
 296   RETURN OLD;
 297 END;
 298 $$
 299 LANGUAGE plpgsql;
 300