]> git.openstreetmap.org Git - nominatim.git/blob - lib-sql/functions/placex_triggers.sql
b47758df3b3bfd3f1c066e75b59859961ea3cae8
[nominatim.git] / lib-sql / functions / placex_triggers.sql
1 -- Trigger functions for the placex table.
2
3 -- Find the parent road of a POI.
4 --
5 -- \returns Place ID of parent object or NULL if none
6 --
7 -- Copy data from linked items (POIs on ways, addr:street links, relations).
8 --
9 CREATE OR REPLACE FUNCTION find_parent_for_poi(poi_osm_type CHAR(1),
10                                                poi_osm_id BIGINT,
11                                                poi_partition SMALLINT,
12                                                bbox GEOMETRY,
13                                                addr_street TEXT,
14                                                addr_place TEXT,
15                                                fallback BOOL = true)
16   RETURNS BIGINT
17   AS $$
18 DECLARE
19   parent_place_id BIGINT DEFAULT NULL;
20   location RECORD;
21   parent RECORD;
22 BEGIN
23     {% if debug %}RAISE WARNING 'finding street for % %', poi_osm_type, poi_osm_id;{% endif %}
24
25     -- Is this object part of an associatedStreet relation?
26     FOR location IN
27       SELECT members FROM planet_osm_rels
28       WHERE parts @> ARRAY[poi_osm_id]
29         and members @> ARRAY[lower(poi_osm_type) || poi_osm_id]
30         and tags @> ARRAY['associatedStreet']
31     LOOP
32       FOR i IN 1..array_upper(location.members, 1) BY 2 LOOP
33         IF location.members[i+1] = 'street' THEN
34           FOR parent IN
35             SELECT place_id from placex
36              WHERE osm_type = 'W' and osm_id = substring(location.members[i],2)::bigint
37                and name is not null
38                and rank_search between 26 and 27
39           LOOP
40             RETURN parent.place_id;
41           END LOOP;
42         END IF;
43       END LOOP;
44     END LOOP;
45
46     parent_place_id := find_parent_for_address(addr_street, addr_place,
47                                                poi_partition, bbox);
48     IF parent_place_id is not null THEN
49       RETURN parent_place_id;
50     END IF;
51
52     IF poi_osm_type = 'N' THEN
53       -- Is this node part of an interpolation?
54       FOR parent IN
55         SELECT q.parent_place_id
56           FROM location_property_osmline q, planet_osm_ways x
57          WHERE q.linegeo && bbox and x.id = q.osm_id
58                and poi_osm_id = any(x.nodes)
59          LIMIT 1
60       LOOP
61         {% if debug %}RAISE WARNING 'Get parent from interpolation: %', parent.parent_place_id;{% endif %}
62         RETURN parent.parent_place_id;
63       END LOOP;
64
65       -- Is this node part of any other way?
66       FOR location IN
67         SELECT p.place_id, p.osm_id, p.rank_search, p.address,
68                coalesce(p.centroid, ST_Centroid(p.geometry)) as centroid
69           FROM placex p, planet_osm_ways w
70          WHERE p.osm_type = 'W' and p.rank_search >= 26
71                and p.geometry && bbox
72                and w.id = p.osm_id and poi_osm_id = any(w.nodes)
73       LOOP
74         {% if debug %}RAISE WARNING 'Node is part of way % ', location.osm_id;{% endif %}
75
76         -- Way IS a road then we are on it - that must be our road
77         IF location.rank_search < 28 THEN
78           {% if debug %}RAISE WARNING 'node in way that is a street %',location;{% endif %}
79           return location.place_id;
80         END IF;
81
82         SELECT find_parent_for_poi('W', location.osm_id, poi_partition,
83                                    location.centroid,
84                                    location.address->'street',
85                                    location.address->'place',
86                                    false)
87           INTO parent_place_id;
88         IF parent_place_id is not null THEN
89           RETURN parent_place_id;
90         END IF;
91       END LOOP;
92     END IF;
93
94     IF fallback THEN
95       IF addr_street is null and addr_place is not null THEN
96         -- The address is attached to a place we don't know.
97         -- Instead simply use the containing area with the largest rank.
98         FOR location IN
99           SELECT place_id FROM placex
100             WHERE bbox && geometry AND _ST_Covers(geometry, ST_Centroid(bbox))
101                   AND rank_address between 5 and 25
102             ORDER BY rank_address desc
103         LOOP
104             RETURN location.place_id;
105         END LOOP;
106       ELSEIF ST_Area(bbox) < 0.005 THEN
107         -- for smaller features get the nearest road
108         SELECT getNearestRoadPlaceId(poi_partition, bbox) INTO parent_place_id;
109         {% if debug %}RAISE WARNING 'Checked for nearest way (%)', parent_place_id;{% endif %}
110       ELSE
111         -- for larger features simply find the area with the largest rank that
112         -- contains the bbox, only use addressable features
113         FOR location IN
114           SELECT place_id FROM placex
115             WHERE bbox && geometry AND _ST_Covers(geometry, ST_Centroid(bbox))
116                   AND rank_address between 5 and 25
117             ORDER BY rank_address desc
118         LOOP
119             RETURN location.place_id;
120         END LOOP;
121       END IF;
122     END IF;
123
124     RETURN parent_place_id;
125 END;
126 $$
127 LANGUAGE plpgsql STABLE;
128
129 -- Try to find a linked place for the given object.
130 CREATE OR REPLACE FUNCTION find_linked_place(bnd placex)
131   RETURNS placex
132   AS $$
133 DECLARE
134   relation_members TEXT[];
135   rel_member RECORD;
136   linked_placex placex%ROWTYPE;
137   bnd_name TEXT;
138 BEGIN
139   IF bnd.rank_search >= 26 or bnd.rank_address = 0
140      or ST_GeometryType(bnd.geometry) NOT IN ('ST_Polygon','ST_MultiPolygon')
141      or bnd.type IN ('postcode', 'postal_code')
142   THEN
143     RETURN NULL;
144   END IF;
145
146   IF bnd.osm_type = 'R' THEN
147     -- see if we have any special relation members
148     SELECT members FROM planet_osm_rels WHERE id = bnd.osm_id INTO relation_members;
149     {% if debug %}RAISE WARNING 'Got relation members';{% endif %}
150
151     -- Search for relation members with role 'lable'.
152     IF relation_members IS NOT NULL THEN
153       FOR rel_member IN
154         SELECT get_rel_node_members(relation_members, ARRAY['label']) as member
155       LOOP
156         {% if debug %}RAISE WARNING 'Found label member %', rel_member.member;{% endif %}
157
158         FOR linked_placex IN
159           SELECT * from placex
160           WHERE osm_type = 'N' and osm_id = rel_member.member
161             and class = 'place'
162         LOOP
163           {% if debug %}RAISE WARNING 'Linked label member';{% endif %}
164           RETURN linked_placex;
165         END LOOP;
166
167       END LOOP;
168     END IF;
169   END IF;
170
171   IF bnd.name ? 'name' THEN
172     bnd_name := make_standard_name(bnd.name->'name');
173     IF bnd_name = '' THEN
174       bnd_name := NULL;
175     END IF;
176   END IF;
177
178   -- If extratags has a place tag, look for linked nodes by their place type.
179   -- Area and node still have to have the same name.
180   IF bnd.extratags ? 'place' and bnd_name is not null THEN
181     FOR linked_placex IN
182       SELECT * FROM placex
183       WHERE make_standard_name(name->'name') = bnd_name
184         AND placex.class = 'place' AND placex.type = bnd.extratags->'place'
185         AND placex.osm_type = 'N'
186         AND placex.linked_place_id is null
187         AND placex.rank_search < 26 -- needed to select the right index
188         AND _st_covers(bnd.geometry, placex.geometry)
189     LOOP
190       {% if debug %}RAISE WARNING 'Found type-matching place node %', linked_placex.osm_id;{% endif %}
191       RETURN linked_placex;
192     END LOOP;
193   END IF;
194
195   IF bnd.extratags ? 'wikidata' THEN
196     FOR linked_placex IN
197       SELECT * FROM placex
198       WHERE placex.class = 'place' AND placex.osm_type = 'N'
199         AND placex.extratags ? 'wikidata' -- needed to select right index
200         AND placex.extratags->'wikidata' = bnd.extratags->'wikidata'
201         AND placex.linked_place_id is null
202         AND placex.rank_search < 26
203         AND _st_covers(bnd.geometry, placex.geometry)
204       ORDER BY make_standard_name(name->'name') = bnd_name desc
205     LOOP
206       {% if debug %}RAISE WARNING 'Found wikidata-matching place node %', linked_placex.osm_id;{% endif %}
207       RETURN linked_placex;
208     END LOOP;
209   END IF;
210
211   -- Name searches can be done for ways as well as relations
212   IF bnd_name is not null THEN
213     {% if debug %}RAISE WARNING 'Looking for nodes with matching names';{% endif %}
214     FOR linked_placex IN
215       SELECT placex.* from placex
216       WHERE make_standard_name(name->'name') = bnd_name
217         AND ((bnd.rank_address > 0
218               and bnd.rank_address = (compute_place_rank(placex.country_code,
219                                                          'N', placex.class,
220                                                          placex.type, 15::SMALLINT,
221                                                          false, placex.postcode)).address_rank)
222              OR (bnd.rank_address = 0 and placex.rank_search = bnd.rank_search))
223         AND placex.osm_type = 'N'
224         AND placex.linked_place_id is null
225         AND placex.rank_search < 26 -- needed to select the right index
226         AND _st_covers(bnd.geometry, placex.geometry)
227     LOOP
228       {% if debug %}RAISE WARNING 'Found matching place node %', linked_placex.osm_id;{% endif %}
229       RETURN linked_placex;
230     END LOOP;
231   END IF;
232
233   RETURN NULL;
234 END;
235 $$
236 LANGUAGE plpgsql STABLE;
237
238
239 -- Insert address of a place into the place_addressline table.
240 --
241 -- \param obj_place_id  Place_id of the place to compute the address for.
242 -- \param partition     Partition number where the place is in.
243 -- \param maxrank       Rank of the place. All address features must have
244 --                      a search rank lower than the given rank.
245 -- \param address       Address terms for the place.
246 -- \param geometry      Geometry to which the address objects should be close.
247 --
248 -- \retval parent_place_id  Place_id of the address object that is the direct
249 --                          ancestor.
250 -- \retval postcode         Postcode computed from the address. This is the
251 --                          addr:postcode of one of the address objects. If
252 --                          more than one of has a postcode, the highest ranking
253 --                          one is used. May be NULL.
254 -- \retval nameaddress_vector  Search terms for the address. This is the sum
255 --                             of name terms of all address objects.
256 CREATE OR REPLACE FUNCTION insert_addresslines(obj_place_id BIGINT,
257                                                partition SMALLINT,
258                                                maxrank SMALLINT,
259                                                address HSTORE,
260                                                geometry GEOMETRY,
261                                                country TEXT,
262                                                OUT parent_place_id BIGINT,
263                                                OUT postcode TEXT,
264                                                OUT nameaddress_vector INT[])
265   AS $$
266 DECLARE
267   address_havelevel BOOLEAN[];
268
269   location_isaddress BOOLEAN;
270   current_boundary GEOMETRY := NULL;
271   current_node_area GEOMETRY := NULL;
272
273   parent_place_rank INT := 0;
274   addr_place_ids BIGINT[];
275
276   location RECORD;
277 BEGIN
278   parent_place_id := 0;
279   nameaddress_vector := '{}'::int[];
280
281   address_havelevel := array_fill(false, ARRAY[maxrank]);
282
283   FOR location IN
284     SELECT * FROM get_places_for_addr_tags(partition, geometry,
285                                                    address, country)
286     ORDER BY rank_address, distance, isguess desc
287   LOOP
288     {% if not db.reverse_only %}
289       nameaddress_vector := array_merge(nameaddress_vector,
290                                         location.keywords::int[]);
291     {% endif %}
292
293     IF location.place_id is not null THEN
294       location_isaddress := not address_havelevel[location.rank_address];
295       IF not address_havelevel[location.rank_address] THEN
296         address_havelevel[location.rank_address] := true;
297         IF parent_place_rank < location.rank_address THEN
298           parent_place_id := location.place_id;
299           parent_place_rank := location.rank_address;
300         END IF;
301       END IF;
302
303       INSERT INTO place_addressline (place_id, address_place_id, fromarea,
304                                      isaddress, distance, cached_rank_address)
305         VALUES (obj_place_id, location.place_id, not location.isguess,
306                 true, location.distance, location.rank_address);
307
308       addr_place_ids := array_append(addr_place_ids, location.place_id);
309     END IF;
310   END LOOP;
311
312   FOR location IN
313     SELECT * FROM getNearFeatures(partition, geometry, maxrank)
314     WHERE addr_place_ids is null or not addr_place_ids @> ARRAY[place_id]
315     ORDER BY rank_address, isguess asc,
316              distance *
317                CASE WHEN rank_address = 16 AND rank_search = 15 THEN 0.2
318                     WHEN rank_address = 16 AND rank_search = 16 THEN 0.25
319                     WHEN rank_address = 16 AND rank_search = 18 THEN 0.5
320                     ELSE 1 END ASC
321   LOOP
322     -- Ignore all place nodes that do not fit in a lower level boundary.
323     CONTINUE WHEN location.isguess
324                   and current_boundary is not NULL
325                   and not ST_Contains(current_boundary, location.centroid);
326
327     -- If this is the first item in the rank, then assume it is the address.
328     location_isaddress := not address_havelevel[location.rank_address];
329
330     -- Further sanity checks to ensure that the address forms a sane hierarchy.
331     IF location_isaddress THEN
332       IF location.isguess and current_node_area is not NULL THEN
333         location_isaddress := ST_Contains(current_node_area, location.centroid);
334       END IF;
335       IF not location.isguess and current_boundary is not NULL
336          and location.rank_address != 11 AND location.rank_address != 5 THEN
337         location_isaddress := ST_Contains(current_boundary, location.centroid);
338       END IF;
339     END IF;
340
341     IF location_isaddress THEN
342       address_havelevel[location.rank_address] := true;
343       parent_place_id := location.place_id;
344
345       -- Set postcode if we have one.
346       -- (Returned will be the highest ranking one.)
347       IF location.postcode is not NULL THEN
348         postcode = location.postcode;
349       END IF;
350
351       -- Recompute the areas we need for hierarchy sanity checks.
352       IF location.rank_address != 11 AND location.rank_address != 5 THEN
353         IF location.isguess THEN
354           current_node_area := place_node_fuzzy_area(location.centroid,
355                                                      location.rank_search);
356         ELSE
357           current_node_area := NULL;
358           SELECT p.geometry FROM placex p
359               WHERE p.place_id = location.place_id INTO current_boundary;
360         END IF;
361       END IF;
362     END IF;
363
364     -- Add it to the list of search terms
365     {% if not db.reverse_only %}
366       nameaddress_vector := array_merge(nameaddress_vector,
367                                         location.keywords::integer[]);
368     {% endif %}
369
370     INSERT INTO place_addressline (place_id, address_place_id, fromarea,
371                                      isaddress, distance, cached_rank_address)
372         VALUES (obj_place_id, location.place_id, not location.isguess,
373                 location_isaddress, location.distance, location.rank_address);
374   END LOOP;
375 END;
376 $$
377 LANGUAGE plpgsql;
378
379
380 CREATE OR REPLACE FUNCTION placex_insert()
381   RETURNS TRIGGER
382   AS $$
383 DECLARE
384   postcode TEXT;
385   result BOOLEAN;
386   is_area BOOLEAN;
387   country_code VARCHAR(2);
388   diameter FLOAT;
389   classtable TEXT;
390 BEGIN
391   {% if debug %}RAISE WARNING '% % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type;{% endif %}
392
393   NEW.place_id := nextval('seq_place');
394   NEW.indexed_status := 1; --STATUS_NEW
395
396   NEW.country_code := lower(get_country_code(NEW.geometry));
397
398   NEW.partition := get_partition(NEW.country_code);
399   NEW.geometry_sector := geometry_sector(NEW.partition, NEW.geometry);
400
401   IF NEW.osm_type = 'X' THEN
402     -- E'X'ternal records should already be in the right format so do nothing
403   ELSE
404     is_area := ST_GeometryType(NEW.geometry) IN ('ST_Polygon','ST_MultiPolygon');
405
406     IF NEW.class in ('place','boundary')
407        AND NEW.type in ('postcode','postal_code')
408     THEN
409       IF NEW.address IS NULL OR NOT NEW.address ? 'postcode' THEN
410           -- most likely just a part of a multipolygon postcode boundary, throw it away
411           RETURN NULL;
412       END IF;
413
414       NEW.name := hstore('ref', NEW.address->'postcode');
415
416     ELSEIF NEW.class = 'highway' AND is_area AND NEW.name is null
417            AND NEW.extratags ? 'area' AND NEW.extratags->'area' = 'yes'
418     THEN
419         RETURN NULL;
420     ELSEIF NEW.class = 'boundary' AND NOT is_area
421     THEN
422         RETURN NULL;
423     ELSEIF NEW.class = 'boundary' AND NEW.type = 'administrative'
424            AND NEW.admin_level <= 4 AND NEW.osm_type = 'W'
425     THEN
426         RETURN NULL;
427     END IF;
428
429     SELECT * INTO NEW.rank_search, NEW.rank_address
430       FROM compute_place_rank(NEW.country_code,
431                               CASE WHEN is_area THEN 'A' ELSE NEW.osm_type END,
432                               NEW.class, NEW.type, NEW.admin_level,
433                               (NEW.extratags->'capital') = 'yes',
434                               NEW.address->'postcode');
435
436     -- a country code make no sense below rank 4 (country)
437     IF NEW.rank_search < 4 THEN
438       NEW.country_code := NULL;
439     END IF;
440
441   END IF;
442
443   {% if debug %}RAISE WARNING 'placex_insert:END: % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type;{% endif %}
444
445 {% if not disable_diff_updates %}
446   -- The following is not needed until doing diff updates, and slows the main index process down
447
448   IF NEW.osm_type = 'N' and NEW.rank_search > 28 THEN
449       -- might be part of an interpolation
450       result := osmline_reinsert(NEW.osm_id, NEW.geometry);
451   ELSEIF NEW.rank_address > 0 THEN
452     IF (ST_GeometryType(NEW.geometry) in ('ST_Polygon','ST_MultiPolygon') AND ST_IsValid(NEW.geometry)) THEN
453       -- Performance: We just can't handle re-indexing for country level changes
454       IF st_area(NEW.geometry) < 1 THEN
455         -- mark items within the geometry for re-indexing
456   --    RAISE WARNING 'placex poly insert: % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type;
457
458         UPDATE placex SET indexed_status = 2
459          WHERE ST_Intersects(NEW.geometry, placex.geometry)
460                and indexed_status = 0
461                and ((rank_address = 0 and rank_search > NEW.rank_address)
462                     or rank_address > NEW.rank_address
463                     or (class = 'place' and osm_type = 'N')
464                    )
465                and (rank_search < 28
466                     or name is not null
467                     or (NEW.rank_address >= 16 and address ? 'place'));
468       END IF;
469     ELSE
470       -- mark nearby items for re-indexing, where 'nearby' depends on the features rank_search and is a complete guess :(
471       diameter := update_place_diameter(NEW.rank_search);
472       IF diameter > 0 THEN
473   --      RAISE WARNING 'placex point insert: % % % % %',NEW.osm_type,NEW.osm_id,NEW.class,NEW.type,diameter;
474         IF NEW.rank_search >= 26 THEN
475           -- roads may cause reparenting for >27 rank places
476           update placex set indexed_status = 2 where indexed_status = 0 and rank_search > NEW.rank_search and ST_DWithin(placex.geometry, NEW.geometry, diameter);
477           -- reparenting also for OSM Interpolation Lines (and for Tiger?)
478           update location_property_osmline set indexed_status = 2 where indexed_status = 0 and ST_DWithin(location_property_osmline.linegeo, NEW.geometry, diameter);
479         ELSEIF NEW.rank_search >= 16 THEN
480           -- up to rank 16, street-less addresses may need reparenting
481           update placex set indexed_status = 2 where indexed_status = 0 and rank_search > NEW.rank_search and ST_DWithin(placex.geometry, NEW.geometry, diameter) and (rank_search < 28 or name is not null or address ? 'place');
482         ELSE
483           -- for all other places the search terms may change as well
484           update placex set indexed_status = 2 where indexed_status = 0 and rank_search > NEW.rank_search and ST_DWithin(placex.geometry, NEW.geometry, diameter) and (rank_search < 28 or name is not null);
485         END IF;
486       END IF;
487     END IF;
488   END IF;
489
490
491    -- add to tables for special search
492    -- Note: won't work on initial import because the classtype tables
493    -- do not yet exist. It won't hurt either.
494   classtable := 'place_classtype_' || NEW.class || '_' || NEW.type;
495   SELECT count(*)>0 FROM pg_tables WHERE tablename = classtable and schemaname = current_schema() INTO result;
496   IF result THEN
497     EXECUTE 'INSERT INTO ' || classtable::regclass || ' (place_id, centroid) VALUES ($1,$2)' 
498     USING NEW.place_id, ST_Centroid(NEW.geometry);
499   END IF;
500
501 {% endif %} -- not disable_diff_updates
502
503   RETURN NEW;
504
505 END;
506 $$
507 LANGUAGE plpgsql;
508
509 CREATE OR REPLACE FUNCTION placex_update()
510   RETURNS TRIGGER
511   AS $$
512 DECLARE
513   i INTEGER;
514   location RECORD;
515   relation_members TEXT[];
516
517   geom GEOMETRY;
518   parent_address_level SMALLINT;
519   place_address_level SMALLINT;
520
521   addr_street TEXT;
522   addr_place TEXT;
523
524   max_rank SMALLINT;
525
526   name_vector INTEGER[];
527   nameaddress_vector INTEGER[];
528   addr_nameaddress_vector INTEGER[];
529
530   inherited_address HSTORE;
531
532   linked_node_id BIGINT;
533   linked_importance FLOAT;
534   linked_wikipedia TEXT;
535
536   result BOOLEAN;
537 BEGIN
538   -- deferred delete
539   IF OLD.indexed_status = 100 THEN
540     {% if debug %}RAISE WARNING 'placex_update delete % %',NEW.osm_type,NEW.osm_id;{% endif %}
541     delete from placex where place_id = OLD.place_id;
542     RETURN NULL;
543   END IF;
544
545   IF NEW.indexed_status != 0 OR OLD.indexed_status = 0 THEN
546     RETURN NEW;
547   END IF;
548
549   {% if debug %}RAISE WARNING 'placex_update % % (%)',NEW.osm_type,NEW.osm_id,NEW.place_id;{% endif %}
550
551   NEW.indexed_date = now();
552
553   {% if 'search_name' in db.tables %}
554     DELETE from search_name WHERE place_id = NEW.place_id;
555   {% endif %}
556   result := deleteSearchName(NEW.partition, NEW.place_id);
557   DELETE FROM place_addressline WHERE place_id = NEW.place_id;
558   result := deleteRoad(NEW.partition, NEW.place_id);
559   result := deleteLocationArea(NEW.partition, NEW.place_id, NEW.rank_search);
560   UPDATE placex set linked_place_id = null, indexed_status = 2
561          where linked_place_id = NEW.place_id;
562   -- update not necessary for osmline, cause linked_place_id does not exist
563
564   NEW.extratags := NEW.extratags - 'linked_place'::TEXT;
565   NEW.address := NEW.address - '_unlisted_place'::TEXT;
566
567   IF NEW.linked_place_id is not null THEN
568     {% if debug %}RAISE WARNING 'place already linked to %', NEW.linked_place_id;{% endif %}
569     RETURN NEW;
570   END IF;
571
572   -- Postcodes are just here to compute the centroids. They are not searchable
573   -- unless they are a boundary=postal_code.
574   -- There was an error in the style so that boundary=postal_code used to be
575   -- imported as place=postcode. That's why relations are allowed to pass here.
576   -- This can go away in a couple of versions.
577   IF NEW.class = 'place'  and NEW.type = 'postcode' and NEW.osm_type != 'R' THEN
578     RETURN NEW;
579   END IF;
580
581   -- Speed up searches - just use the centroid of the feature
582   -- cheaper but less acurate
583   NEW.centroid := ST_PointOnSurface(NEW.geometry);
584   {% if debug %}RAISE WARNING 'Computing preliminary centroid at %',ST_AsText(NEW.centroid);{% endif %}
585
586   -- recompute the ranks, they might change when linking changes
587   SELECT * INTO NEW.rank_search, NEW.rank_address
588     FROM compute_place_rank(NEW.country_code,
589                             CASE WHEN ST_GeometryType(NEW.geometry)
590                                         IN ('ST_Polygon','ST_MultiPolygon')
591                             THEN 'A' ELSE NEW.osm_type END,
592                             NEW.class, NEW.type, NEW.admin_level,
593                             (NEW.extratags->'capital') = 'yes',
594                             NEW.address->'postcode');
595   -- We must always increase the address level relative to the admin boundary.
596   IF NEW.class = 'boundary' and NEW.type = 'administrative'
597      and NEW.osm_type = 'R' and NEW.rank_address > 0
598   THEN
599     -- First, check that admin boundaries do not overtake each other rank-wise.
600     parent_address_level := 3;
601     FOR location IN
602       SELECT rank_address,
603              (CASE WHEN extratags ? 'wikidata' and NEW.extratags ? 'wikidata'
604                         and extratags->'wikidata' = NEW.extratags->'wikidata'
605                    THEN ST_Equals(geometry, NEW.geometry)
606                    ELSE false END) as is_same
607       FROM placex
608       WHERE osm_type = 'R' and class = 'boundary' and type = 'administrative'
609             and admin_level < NEW.admin_level and admin_level > 3
610             and rank_address > 0
611             and geometry && NEW.centroid and _ST_Covers(geometry, NEW.centroid)
612       ORDER BY admin_level desc LIMIT 1
613     LOOP
614       IF location.is_same THEN
615         -- Looks like the same boundary is replicated on multiple admin_levels.
616         -- Usual tagging in Poland. Remove our boundary from addresses.
617         NEW.rank_address := 0;
618       ELSE
619         parent_address_level := location.rank_address;
620         IF location.rank_address >= NEW.rank_address THEN
621           IF location.rank_address >= 24 THEN
622             NEW.rank_address := 25;
623           ELSE
624             NEW.rank_address := location.rank_address + 2;
625           END IF;
626         END IF;
627       END IF;
628     END LOOP;
629
630     IF NEW.rank_address > 9 THEN
631         -- Second check that the boundary is not completely contained in a
632         -- place area with a higher address rank
633         FOR location IN
634           SELECT rank_address FROM placex
635           WHERE class = 'place' and rank_address < 24
636                 and rank_address > NEW.rank_address
637                 and geometry && NEW.geometry
638                 and geometry ~ NEW.geometry -- needed because ST_Relate does not do bbox cover test
639                 and ST_Relate(geometry, NEW.geometry, 'T*T***FF*') -- contains but not equal
640           ORDER BY rank_address desc LIMIT 1
641         LOOP
642           NEW.rank_address := location.rank_address + 2;
643         END LOOP;
644     END IF;
645   ELSEIF NEW.class = 'place' and NEW.osm_type = 'N'
646      and NEW.rank_address between 16 and 23
647   THEN
648     -- If a place node is contained in a admin boundary with the same address level
649     -- and has not been linked, then make the node a subpart by increasing the
650     -- address rank (city level and above).
651     FOR location IN
652         SELECT rank_address FROM placex
653         WHERE osm_type = 'R' and class = 'boundary' and type = 'administrative'
654               and rank_address = NEW.rank_address
655               and geometry && NEW.centroid and _ST_Covers(geometry, NEW.centroid)
656         LIMIT 1
657     LOOP
658       NEW.rank_address = NEW.rank_address + 2;
659     END LOOP;
660   ELSE
661     parent_address_level := 3;
662   END IF;
663
664   {% if debug %}RAISE WARNING 'Copy over address tags';{% endif %}
665   -- housenumber is a computed field, so start with an empty value
666   NEW.housenumber := NULL;
667   IF NEW.address is not NULL THEN
668       IF NEW.address ? 'conscriptionnumber' THEN
669         IF NEW.address ? 'streetnumber' THEN
670             NEW.housenumber := (NEW.address->'conscriptionnumber') || ';' || (NEW.address->'streetnumber');
671         ELSE
672             NEW.housenumber := NEW.address->'conscriptionnumber';
673         END IF;
674       ELSEIF NEW.address ? 'streetnumber' THEN
675         NEW.housenumber := NEW.address->'streetnumber';
676       ELSEIF NEW.address ? 'housenumber' THEN
677         NEW.housenumber := NEW.address->'housenumber';
678       END IF;
679       NEW.housenumber := create_housenumber_id(NEW.housenumber);
680
681       addr_street := NEW.address->'street';
682       addr_place := NEW.address->'place';
683
684       IF NEW.address ? 'postcode' and NEW.address->'postcode' not similar to '%(:|,|;)%' THEN
685         i := getorcreate_postcode_id(NEW.address->'postcode');
686       END IF;
687   END IF;
688
689   NEW.postcode := null;
690
691   -- recalculate country and partition
692   IF NEW.rank_search = 4 AND NEW.address is not NULL AND NEW.address ? 'country' THEN
693     -- for countries, believe the mapped country code,
694     -- so that we remain in the right partition if the boundaries
695     -- suddenly expand.
696     NEW.country_code := lower(NEW.address->'country');
697     NEW.partition := get_partition(lower(NEW.country_code));
698     IF NEW.partition = 0 THEN
699       NEW.country_code := lower(get_country_code(NEW.centroid));
700       NEW.partition := get_partition(NEW.country_code);
701     END IF;
702   ELSE
703     IF NEW.rank_search >= 4 THEN
704       NEW.country_code := lower(get_country_code(NEW.centroid));
705     ELSE
706       NEW.country_code := NULL;
707     END IF;
708     NEW.partition := get_partition(NEW.country_code);
709   END IF;
710   {% if debug %}RAISE WARNING 'Country updated: "%"', NEW.country_code;{% endif %}
711
712   -- waterway ways are linked when they are part of a relation and have the same class/type
713   IF NEW.osm_type = 'R' and NEW.class = 'waterway' THEN
714       FOR relation_members IN select members from planet_osm_rels r where r.id = NEW.osm_id and r.parts != array[]::bigint[]
715       LOOP
716           FOR i IN 1..array_upper(relation_members, 1) BY 2 LOOP
717               IF relation_members[i+1] in ('', 'main_stream', 'side_stream') AND substring(relation_members[i],1,1) = 'w' THEN
718                 {% if debug %}RAISE WARNING 'waterway parent %, child %/%', NEW.osm_id, i, relation_members[i];{% endif %}
719                 FOR linked_node_id IN SELECT place_id FROM placex
720                   WHERE osm_type = 'W' and osm_id = substring(relation_members[i],2,200)::bigint
721                   and class = NEW.class and type in ('river', 'stream', 'canal', 'drain', 'ditch')
722                   and ( relation_members[i+1] != 'side_stream' or NEW.name->'name' = name->'name')
723                 LOOP
724                   UPDATE placex SET linked_place_id = NEW.place_id WHERE place_id = linked_node_id;
725                   {% if 'search_name' in db.tables %}
726                     DELETE FROM search_name WHERE place_id = linked_node_id;
727                   {% endif %}
728                 END LOOP;
729               END IF;
730           END LOOP;
731       END LOOP;
732       {% if debug %}RAISE WARNING 'Waterway processed';{% endif %}
733   END IF;
734
735   NEW.importance := null;
736   SELECT wikipedia, importance
737     FROM compute_importance(NEW.extratags, NEW.country_code, NEW.osm_type, NEW.osm_id)
738     INTO NEW.wikipedia,NEW.importance;
739
740 {% if debug %}RAISE WARNING 'Importance computed from wikipedia: %', NEW.importance;{% endif %}
741
742   -- ---------------------------------------------------------------------------
743   -- For low level elements we inherit from our parent road
744   IF NEW.rank_search > 27 THEN
745
746     {% if debug %}RAISE WARNING 'finding street for % %', NEW.osm_type, NEW.osm_id;{% endif %}
747     NEW.parent_place_id := null;
748
749     -- if we have a POI and there is no address information,
750     -- see if we can get it from a surrounding building
751     inherited_address := ''::HSTORE;
752     IF NEW.osm_type = 'N' AND addr_street IS NULL AND addr_place IS NULL
753        AND NEW.housenumber IS NULL THEN
754       FOR location IN
755         -- The additional && condition works around the misguided query
756         -- planner of postgis 3.0.
757         SELECT address from placex where ST_Covers(geometry, NEW.centroid)
758             and geometry && NEW.centroid
759             and (address ? 'housenumber' or address ? 'street' or address ? 'place')
760             and rank_search > 28 AND ST_GeometryType(geometry) in ('ST_Polygon','ST_MultiPolygon')
761             limit 1
762       LOOP
763         NEW.housenumber := location.address->'housenumber';
764         addr_street := location.address->'street';
765         addr_place := location.address->'place';
766         inherited_address := location.address;
767       END LOOP;
768     END IF;
769
770     -- We have to find our parent road.
771     NEW.parent_place_id := find_parent_for_poi(NEW.osm_type, NEW.osm_id,
772                                                NEW.partition,
773                                                ST_Envelope(NEW.geometry),
774                                                addr_street, addr_place);
775
776     -- If we found the road take a shortcut here.
777     -- Otherwise fall back to the full address getting method below.
778     IF NEW.parent_place_id is not null THEN
779
780       -- Get the details of the parent road
781       SELECT p.country_code, p.postcode, p.name FROM placex p
782        WHERE p.place_id = NEW.parent_place_id INTO location;
783
784       IF addr_street is null and addr_place is not null THEN
785         -- Check if the addr:place tag is part of the parent name
786         SELECT count(*) INTO i
787           FROM svals(location.name) AS pname WHERE pname = addr_place;
788         IF i = 0 THEN
789           NEW.address = NEW.address || hstore('_unlisted_place', addr_place);
790         END IF;
791       END IF;
792
793       NEW.country_code := location.country_code;
794       {% if debug %}RAISE WARNING 'Got parent details from search name';{% endif %}
795
796       -- determine postcode
797       IF NEW.address is not null AND NEW.address ? 'postcode' THEN
798           NEW.postcode = upper(trim(NEW.address->'postcode'));
799       ELSE
800          NEW.postcode := location.postcode;
801       END IF;
802       IF NEW.postcode is null THEN
803         NEW.postcode := get_nearest_postcode(NEW.country_code, NEW.geometry);
804       END IF;
805
806       IF NEW.name is not NULL THEN
807           NEW.name := add_default_place_name(NEW.country_code, NEW.name);
808           name_vector := make_keywords(NEW.name);
809
810           IF NEW.rank_search <= 25 and NEW.rank_address > 0 THEN
811             result := add_location(NEW.place_id, NEW.country_code, NEW.partition,
812                                    name_vector, NEW.rank_search, NEW.rank_address,
813                                    upper(trim(NEW.address->'postcode')), NEW.geometry,
814                                    NEW.centroid);
815             {% if debug %}RAISE WARNING 'Place added to location table';{% endif %}
816           END IF;
817
818       END IF;
819
820       {% if not db.reverse_only %}
821       IF array_length(name_vector, 1) is not NULL
822          OR inherited_address is not NULL OR NEW.address is not NULL
823       THEN
824         SELECT * INTO name_vector, nameaddress_vector
825           FROM create_poi_search_terms(NEW.place_id,
826                                        NEW.partition, NEW.parent_place_id,
827                                        inherited_address || NEW.address,
828                                        NEW.country_code, NEW.housenumber,
829                                        name_vector, NEW.centroid);
830
831         IF array_length(name_vector, 1) is not NULL THEN
832           INSERT INTO search_name (place_id, search_rank, address_rank,
833                                    importance, country_code, name_vector,
834                                    nameaddress_vector, centroid)
835                  VALUES (NEW.place_id, NEW.rank_search, NEW.rank_address,
836                          NEW.importance, NEW.country_code, name_vector,
837                          nameaddress_vector, NEW.centroid);
838           {% if debug %}RAISE WARNING 'Place added to search table';{% endif %}
839         END IF;
840       END IF;
841       {% endif %}
842
843       RETURN NEW;
844     END IF;
845
846   END IF;
847
848   -- ---------------------------------------------------------------------------
849   -- Full indexing
850   {% if debug %}RAISE WARNING 'Using full index mode for % %', NEW.osm_type, NEW.osm_id;{% endif %}
851   SELECT * INTO location FROM find_linked_place(NEW);
852   IF location.place_id is not null THEN
853     {% if debug %}RAISE WARNING 'Linked %', location;{% endif %}
854
855     -- Use the linked point as the centre point of the geometry,
856     -- but only if it is within the area of the boundary.
857     geom := coalesce(location.centroid, ST_Centroid(location.geometry));
858     IF geom is not NULL AND ST_Within(geom, NEW.geometry) THEN
859         NEW.centroid := geom;
860     END IF;
861
862     {% if debug %}RAISE WARNING 'parent address: % rank address: %', parent_address_level, location.rank_address;{% endif %}
863     IF location.rank_address > parent_address_level
864        and location.rank_address < 26
865     THEN
866       NEW.rank_address := location.rank_address;
867     END IF;
868
869     -- merge in the label name
870     IF NOT location.name IS NULL THEN
871       NEW.name := location.name || NEW.name;
872     END IF;
873
874     -- merge in extra tags
875     NEW.extratags := hstore('linked_' || location.class, location.type)
876                      || coalesce(location.extratags, ''::hstore)
877                      || coalesce(NEW.extratags, ''::hstore);
878
879     -- mark the linked place (excludes from search results)
880     UPDATE placex set linked_place_id = NEW.place_id
881       WHERE place_id = location.place_id;
882     -- ensure that those places are not found anymore
883     {% if 'search_name' in db.tables %}
884       DELETE FROM search_name WHERE place_id = location.place_id;
885     {% endif %}
886     PERFORM deleteLocationArea(NEW.partition, location.place_id, NEW.rank_search);
887
888     SELECT wikipedia, importance
889       FROM compute_importance(location.extratags, NEW.country_code,
890                               'N', location.osm_id)
891       INTO linked_wikipedia,linked_importance;
892
893     -- Use the maximum importance if one could be computed from the linked object.
894     IF linked_importance is not null AND
895        (NEW.importance is null or NEW.importance < linked_importance)
896     THEN
897       NEW.importance = linked_importance;
898     END IF;
899   ELSE
900     -- No linked place? As a last resort check if the boundary is tagged with
901     -- a place type and adapt the rank address.
902     IF NEW.rank_address > 0 and NEW.extratags ? 'place' THEN
903       SELECT address_rank INTO place_address_level
904         FROM compute_place_rank(NEW.country_code, 'A', 'place',
905                                 NEW.extratags->'place', 0::SMALLINT, False, null);
906       IF place_address_level > parent_address_level and
907          place_address_level < 26 THEN
908         NEW.rank_address := place_address_level;
909       END IF;
910     END IF;
911   END IF;
912
913   -- Initialise the name vector using our name
914   NEW.name := add_default_place_name(NEW.country_code, NEW.name);
915   name_vector := make_keywords(NEW.name);
916
917   -- make sure all names are in the word table
918   IF NEW.admin_level = 2
919      AND NEW.class = 'boundary' AND NEW.type = 'administrative'
920      AND NEW.country_code IS NOT NULL AND NEW.osm_type = 'R'
921   THEN
922     PERFORM create_country(NEW.name, lower(NEW.country_code));
923     {% if debug %}RAISE WARNING 'Country names updated';{% endif %}
924
925     -- Also update the list of country names. Adding an additional sanity
926     -- check here: make sure the country does overlap with the area where
927     -- we expect it to be as per static country grid.
928     FOR location IN
929       SELECT country_code FROM country_osm_grid
930        WHERE ST_Covers(geometry, NEW.centroid) and country_code = NEW.country_code
931        LIMIT 1
932     LOOP
933       {% if debug %}RAISE WARNING 'Updating names for country '%' with: %', NEW.country_code, NEW.name;{% endif %}
934       UPDATE country_name SET name = name || NEW.name WHERE country_code = NEW.country_code;
935     END LOOP;
936   END IF;
937
938   -- For linear features we need the full geometry for determining the address
939   -- because they may go through several administrative entities. Otherwise use
940   -- the centroid for performance reasons.
941   IF ST_GeometryType(NEW.geometry) in ('ST_LineString', 'ST_MultiLineString') THEN
942     geom := NEW.geometry;
943   ELSE
944     geom := NEW.centroid;
945   END IF;
946
947   IF NEW.rank_address = 0 THEN
948     max_rank := geometry_to_rank(NEW.rank_search, NEW.geometry, NEW.country_code);
949     -- Rank 0 features may also span multiple administrative areas (e.g. lakes)
950     -- so use the geometry here too. Just make sure the areas don't become too
951     -- large.
952     IF NEW.class = 'natural' or max_rank > 10 THEN
953       geom := NEW.geometry;
954     END IF;
955   ELSEIF NEW.rank_address > 25 THEN
956     max_rank := 25;
957   ELSE
958     max_rank = NEW.rank_address;
959   END IF;
960
961   SELECT * FROM insert_addresslines(NEW.place_id, NEW.partition, max_rank,
962                                     NEW.address, geom, NEW.country_code)
963     INTO NEW.parent_place_id, NEW.postcode, nameaddress_vector;
964
965   {% if debug %}RAISE WARNING 'RETURN insert_addresslines: %, %, %', NEW.parent_place_id, NEW.postcode, nameaddress_vector;{% endif %}
966
967   IF NEW.address is not null AND NEW.address ? 'postcode' 
968      AND NEW.address->'postcode' not similar to '%(,|;)%' THEN
969     NEW.postcode := upper(trim(NEW.address->'postcode'));
970   END IF;
971
972   IF NEW.postcode is null AND NEW.rank_search > 8 THEN
973     NEW.postcode := get_nearest_postcode(NEW.country_code, NEW.geometry);
974   END IF;
975
976   -- if we have a name add this to the name search table
977   IF NEW.name IS NOT NULL THEN
978
979     IF NEW.rank_search <= 25 and NEW.rank_address > 0 THEN
980       result := add_location(NEW.place_id, NEW.country_code, NEW.partition, name_vector, NEW.rank_search, NEW.rank_address, upper(trim(NEW.address->'postcode')), NEW.geometry, NEW.centroid);
981       {% if debug %}RAISE WARNING 'added to location (full)';{% endif %}
982     END IF;
983
984     IF NEW.rank_search between 26 and 27 and NEW.class = 'highway' THEN
985       result := insertLocationRoad(NEW.partition, NEW.place_id, NEW.country_code, NEW.geometry);
986       {% if debug %}RAISE WARNING 'insert into road location table (full)';{% endif %}
987     END IF;
988
989     result := insertSearchName(NEW.partition, NEW.place_id, name_vector,
990                                NEW.rank_search, NEW.rank_address, NEW.geometry);
991     {% if debug %}RAISE WARNING 'added to search name (full)';{% endif %}
992
993     {% if not db.reverse_only %}
994         INSERT INTO search_name (place_id, search_rank, address_rank,
995                                  importance, country_code, name_vector,
996                                  nameaddress_vector, centroid)
997                VALUES (NEW.place_id, NEW.rank_search, NEW.rank_address,
998                        NEW.importance, NEW.country_code, name_vector,
999                        nameaddress_vector, NEW.centroid);
1000     {% endif %}
1001
1002   END IF;
1003
1004   {% if debug %}RAISE WARNING 'place update % % finsihed.', NEW.osm_type, NEW.osm_id;{% endif %}
1005
1006   RETURN NEW;
1007 END;
1008 $$
1009 LANGUAGE plpgsql;
1010
1011
1012 CREATE OR REPLACE FUNCTION placex_delete()
1013   RETURNS TRIGGER
1014   AS $$
1015 DECLARE
1016   b BOOLEAN;
1017   classtable TEXT;
1018 BEGIN
1019   -- RAISE WARNING 'placex_delete % %',OLD.osm_type,OLD.osm_id;
1020
1021   IF OLD.linked_place_id is null THEN
1022     update placex set linked_place_id = null, indexed_status = 2 where linked_place_id = OLD.place_id and indexed_status = 0;
1023     {% if debug %}RAISE WARNING 'placex_delete:01 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1024     update placex set linked_place_id = null where linked_place_id = OLD.place_id;
1025     {% if debug %}RAISE WARNING 'placex_delete:02 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1026   ELSE
1027     update placex set indexed_status = 2 where place_id = OLD.linked_place_id and indexed_status = 0;
1028   END IF;
1029
1030   IF OLD.rank_address < 30 THEN
1031
1032     -- mark everything linked to this place for re-indexing
1033     {% if debug %}RAISE WARNING 'placex_delete:03 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1034     UPDATE placex set indexed_status = 2 from place_addressline where address_place_id = OLD.place_id 
1035       and placex.place_id = place_addressline.place_id and indexed_status = 0 and place_addressline.isaddress;
1036
1037     {% if debug %}RAISE WARNING 'placex_delete:04 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1038     DELETE FROM place_addressline where address_place_id = OLD.place_id;
1039
1040     {% if debug %}RAISE WARNING 'placex_delete:05 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1041     b := deleteRoad(OLD.partition, OLD.place_id);
1042
1043     {% if debug %}RAISE WARNING 'placex_delete:06 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1044     update placex set indexed_status = 2 where parent_place_id = OLD.place_id and indexed_status = 0;
1045     {% if debug %}RAISE WARNING 'placex_delete:07 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1046     -- reparenting also for OSM Interpolation Lines (and for Tiger?)
1047     update location_property_osmline set indexed_status = 2 where indexed_status = 0 and parent_place_id = OLD.place_id;
1048
1049   END IF;
1050
1051   {% if debug %}RAISE WARNING 'placex_delete:08 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1052
1053   IF OLD.rank_address < 26 THEN
1054     b := deleteLocationArea(OLD.partition, OLD.place_id, OLD.rank_search);
1055   END IF;
1056
1057   {% if debug %}RAISE WARNING 'placex_delete:09 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1058
1059   IF OLD.name is not null THEN
1060     {% if 'search_name' in db.tables %}
1061       DELETE from search_name WHERE place_id = OLD.place_id;
1062     {% endif %}
1063     b := deleteSearchName(OLD.partition, OLD.place_id);
1064   END IF;
1065
1066   {% if debug %}RAISE WARNING 'placex_delete:10 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1067
1068   DELETE FROM place_addressline where place_id = OLD.place_id;
1069
1070   {% if debug %}RAISE WARNING 'placex_delete:11 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1071
1072   -- remove from tables for special search
1073   classtable := 'place_classtype_' || OLD.class || '_' || OLD.type;
1074   SELECT count(*)>0 FROM pg_tables WHERE tablename = classtable and schemaname = current_schema() INTO b;
1075   IF b THEN
1076     EXECUTE 'DELETE FROM ' || classtable::regclass || ' WHERE place_id = $1' USING OLD.place_id;
1077   END IF;
1078
1079   {% if debug %}RAISE WARNING 'placex_delete:12 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1080
1081   RETURN OLD;
1082
1083 END;
1084 $$
1085 LANGUAGE plpgsql;