]> git.openstreetmap.org Git - nominatim.git/blob - lib-sql/functions/placex_triggers.sql
move transliteration of housenumbers into indexing
[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         i := getorcreate_housenumber_id(make_standard_name(NEW.address->'conscriptionnumber'));
670         IF NEW.address ? 'streetnumber' THEN
671             i := getorcreate_housenumber_id(make_standard_name(NEW.address->'streetnumber'));
672             NEW.housenumber := (NEW.address->'conscriptionnumber') || '/' || (NEW.address->'streetnumber');
673         ELSE
674             NEW.housenumber := NEW.address->'conscriptionnumber';
675         END IF;
676       ELSEIF NEW.address ? 'streetnumber' THEN
677         NEW.housenumber := NEW.address->'streetnumber';
678         i := getorcreate_housenumber_id(make_standard_name(NEW.address->'streetnumber'));
679       ELSEIF NEW.address ? 'housenumber' THEN
680         NEW.housenumber := NEW.address->'housenumber';
681         i := getorcreate_housenumber_id(make_standard_name(NEW.housenumber));
682       END IF;
683       NEW.housenumber := transliteration(NEW.housenumber);
684
685       addr_street := NEW.address->'street';
686       addr_place := NEW.address->'place';
687
688       IF NEW.address ? 'postcode' and NEW.address->'postcode' not similar to '%(:|,|;)%' THEN
689         i := getorcreate_postcode_id(NEW.address->'postcode');
690       END IF;
691   END IF;
692
693   NEW.postcode := null;
694
695   -- recalculate country and partition
696   IF NEW.rank_search = 4 AND NEW.address is not NULL AND NEW.address ? 'country' THEN
697     -- for countries, believe the mapped country code,
698     -- so that we remain in the right partition if the boundaries
699     -- suddenly expand.
700     NEW.country_code := lower(NEW.address->'country');
701     NEW.partition := get_partition(lower(NEW.country_code));
702     IF NEW.partition = 0 THEN
703       NEW.country_code := lower(get_country_code(NEW.centroid));
704       NEW.partition := get_partition(NEW.country_code);
705     END IF;
706   ELSE
707     IF NEW.rank_search >= 4 THEN
708       NEW.country_code := lower(get_country_code(NEW.centroid));
709     ELSE
710       NEW.country_code := NULL;
711     END IF;
712     NEW.partition := get_partition(NEW.country_code);
713   END IF;
714   {% if debug %}RAISE WARNING 'Country updated: "%"', NEW.country_code;{% endif %}
715
716   -- waterway ways are linked when they are part of a relation and have the same class/type
717   IF NEW.osm_type = 'R' and NEW.class = 'waterway' THEN
718       FOR relation_members IN select members from planet_osm_rels r where r.id = NEW.osm_id and r.parts != array[]::bigint[]
719       LOOP
720           FOR i IN 1..array_upper(relation_members, 1) BY 2 LOOP
721               IF relation_members[i+1] in ('', 'main_stream', 'side_stream') AND substring(relation_members[i],1,1) = 'w' THEN
722                 {% if debug %}RAISE WARNING 'waterway parent %, child %/%', NEW.osm_id, i, relation_members[i];{% endif %}
723                 FOR linked_node_id IN SELECT place_id FROM placex
724                   WHERE osm_type = 'W' and osm_id = substring(relation_members[i],2,200)::bigint
725                   and class = NEW.class and type in ('river', 'stream', 'canal', 'drain', 'ditch')
726                   and ( relation_members[i+1] != 'side_stream' or NEW.name->'name' = name->'name')
727                 LOOP
728                   UPDATE placex SET linked_place_id = NEW.place_id WHERE place_id = linked_node_id;
729                   {% if 'search_name' in db.tables %}
730                     DELETE FROM search_name WHERE place_id = linked_node_id;
731                   {% endif %}
732                 END LOOP;
733               END IF;
734           END LOOP;
735       END LOOP;
736       {% if debug %}RAISE WARNING 'Waterway processed';{% endif %}
737   END IF;
738
739   NEW.importance := null;
740   SELECT wikipedia, importance
741     FROM compute_importance(NEW.extratags, NEW.country_code, NEW.osm_type, NEW.osm_id)
742     INTO NEW.wikipedia,NEW.importance;
743
744 {% if debug %}RAISE WARNING 'Importance computed from wikipedia: %', NEW.importance;{% endif %}
745
746   -- ---------------------------------------------------------------------------
747   -- For low level elements we inherit from our parent road
748   IF NEW.rank_search > 27 THEN
749
750     {% if debug %}RAISE WARNING 'finding street for % %', NEW.osm_type, NEW.osm_id;{% endif %}
751     NEW.parent_place_id := null;
752
753     -- if we have a POI and there is no address information,
754     -- see if we can get it from a surrounding building
755     inherited_address := ''::HSTORE;
756     IF NEW.osm_type = 'N' AND addr_street IS NULL AND addr_place IS NULL
757        AND NEW.housenumber IS NULL THEN
758       FOR location IN
759         -- The additional && condition works around the misguided query
760         -- planner of postgis 3.0.
761         SELECT address from placex where ST_Covers(geometry, NEW.centroid)
762             and geometry && NEW.centroid
763             and (address ? 'housenumber' or address ? 'street' or address ? 'place')
764             and rank_search > 28 AND ST_GeometryType(geometry) in ('ST_Polygon','ST_MultiPolygon')
765             limit 1
766       LOOP
767         NEW.housenumber := location.address->'housenumber';
768         addr_street := location.address->'street';
769         addr_place := location.address->'place';
770         inherited_address := location.address;
771       END LOOP;
772     END IF;
773
774     -- We have to find our parent road.
775     NEW.parent_place_id := find_parent_for_poi(NEW.osm_type, NEW.osm_id,
776                                                NEW.partition,
777                                                ST_Envelope(NEW.geometry),
778                                                addr_street, addr_place);
779
780     -- If we found the road take a shortcut here.
781     -- Otherwise fall back to the full address getting method below.
782     IF NEW.parent_place_id is not null THEN
783
784       -- Get the details of the parent road
785       SELECT p.country_code, p.postcode, p.name FROM placex p
786        WHERE p.place_id = NEW.parent_place_id INTO location;
787
788       IF addr_street is null and addr_place is not null THEN
789         -- Check if the addr:place tag is part of the parent name
790         SELECT count(*) INTO i
791           FROM svals(location.name) AS pname WHERE pname = addr_place;
792         IF i = 0 THEN
793           NEW.address = NEW.address || hstore('_unlisted_place', addr_place);
794         END IF;
795       END IF;
796
797       NEW.country_code := location.country_code;
798       {% if debug %}RAISE WARNING 'Got parent details from search name';{% endif %}
799
800       -- determine postcode
801       IF NEW.address is not null AND NEW.address ? 'postcode' THEN
802           NEW.postcode = upper(trim(NEW.address->'postcode'));
803       ELSE
804          NEW.postcode := location.postcode;
805       END IF;
806       IF NEW.postcode is null THEN
807         NEW.postcode := get_nearest_postcode(NEW.country_code, NEW.geometry);
808       END IF;
809
810       IF NEW.name is not NULL THEN
811           NEW.name := add_default_place_name(NEW.country_code, NEW.name);
812           name_vector := make_keywords(NEW.name);
813
814           IF NEW.rank_search <= 25 and NEW.rank_address > 0 THEN
815             result := add_location(NEW.place_id, NEW.country_code, NEW.partition,
816                                    name_vector, NEW.rank_search, NEW.rank_address,
817                                    upper(trim(NEW.address->'postcode')), NEW.geometry,
818                                    NEW.centroid);
819             {% if debug %}RAISE WARNING 'Place added to location table';{% endif %}
820           END IF;
821
822       END IF;
823
824       {% if not db.reverse_only %}
825       IF array_length(name_vector, 1) is not NULL
826          OR inherited_address is not NULL OR NEW.address is not NULL
827       THEN
828         SELECT * INTO name_vector, nameaddress_vector
829           FROM create_poi_search_terms(NEW.place_id,
830                                        NEW.partition, NEW.parent_place_id,
831                                        inherited_address || NEW.address,
832                                        NEW.country_code, NEW.housenumber,
833                                        name_vector, NEW.centroid);
834
835         IF array_length(name_vector, 1) is not NULL THEN
836           INSERT INTO search_name (place_id, search_rank, address_rank,
837                                    importance, country_code, name_vector,
838                                    nameaddress_vector, centroid)
839                  VALUES (NEW.place_id, NEW.rank_search, NEW.rank_address,
840                          NEW.importance, NEW.country_code, name_vector,
841                          nameaddress_vector, NEW.centroid);
842           {% if debug %}RAISE WARNING 'Place added to search table';{% endif %}
843         END IF;
844       END IF;
845       {% endif %}
846
847       RETURN NEW;
848     END IF;
849
850   END IF;
851
852   -- ---------------------------------------------------------------------------
853   -- Full indexing
854   {% if debug %}RAISE WARNING 'Using full index mode for % %', NEW.osm_type, NEW.osm_id;{% endif %}
855   SELECT * INTO location FROM find_linked_place(NEW);
856   IF location.place_id is not null THEN
857     {% if debug %}RAISE WARNING 'Linked %', location;{% endif %}
858
859     -- Use the linked point as the centre point of the geometry,
860     -- but only if it is within the area of the boundary.
861     geom := coalesce(location.centroid, ST_Centroid(location.geometry));
862     IF geom is not NULL AND ST_Within(geom, NEW.geometry) THEN
863         NEW.centroid := geom;
864     END IF;
865
866     {% if debug %}RAISE WARNING 'parent address: % rank address: %', parent_address_level, location.rank_address;{% endif %}
867     IF location.rank_address > parent_address_level
868        and location.rank_address < 26
869     THEN
870       NEW.rank_address := location.rank_address;
871     END IF;
872
873     -- merge in the label name
874     IF NOT location.name IS NULL THEN
875       NEW.name := location.name || NEW.name;
876     END IF;
877
878     -- merge in extra tags
879     NEW.extratags := hstore('linked_' || location.class, location.type)
880                      || coalesce(location.extratags, ''::hstore)
881                      || coalesce(NEW.extratags, ''::hstore);
882
883     -- mark the linked place (excludes from search results)
884     UPDATE placex set linked_place_id = NEW.place_id
885       WHERE place_id = location.place_id;
886     -- ensure that those places are not found anymore
887     {% if 'search_name' in db.tables %}
888       DELETE FROM search_name WHERE place_id = location.place_id;
889     {% endif %}
890     PERFORM deleteLocationArea(NEW.partition, location.place_id, NEW.rank_search);
891
892     SELECT wikipedia, importance
893       FROM compute_importance(location.extratags, NEW.country_code,
894                               'N', location.osm_id)
895       INTO linked_wikipedia,linked_importance;
896
897     -- Use the maximum importance if one could be computed from the linked object.
898     IF linked_importance is not null AND
899        (NEW.importance is null or NEW.importance < linked_importance)
900     THEN
901       NEW.importance = linked_importance;
902     END IF;
903   ELSE
904     -- No linked place? As a last resort check if the boundary is tagged with
905     -- a place type and adapt the rank address.
906     IF NEW.rank_address > 0 and NEW.extratags ? 'place' THEN
907       SELECT address_rank INTO place_address_level
908         FROM compute_place_rank(NEW.country_code, 'A', 'place',
909                                 NEW.extratags->'place', 0::SMALLINT, False, null);
910       IF place_address_level > parent_address_level and
911          place_address_level < 26 THEN
912         NEW.rank_address := place_address_level;
913       END IF;
914     END IF;
915   END IF;
916
917   -- Initialise the name vector using our name
918   NEW.name := add_default_place_name(NEW.country_code, NEW.name);
919   name_vector := make_keywords(NEW.name);
920
921   -- make sure all names are in the word table
922   IF NEW.admin_level = 2
923      AND NEW.class = 'boundary' AND NEW.type = 'administrative'
924      AND NEW.country_code IS NOT NULL AND NEW.osm_type = 'R'
925   THEN
926     PERFORM create_country(NEW.name, lower(NEW.country_code));
927     {% if debug %}RAISE WARNING 'Country names updated';{% endif %}
928
929     -- Also update the list of country names. Adding an additional sanity
930     -- check here: make sure the country does overlap with the area where
931     -- we expect it to be as per static country grid.
932     FOR location IN
933       SELECT country_code FROM country_osm_grid
934        WHERE ST_Covers(geometry, NEW.centroid) and country_code = NEW.country_code
935        LIMIT 1
936     LOOP
937       {% if debug %}RAISE WARNING 'Updating names for country '%' with: %', NEW.country_code, NEW.name;{% endif %}
938       UPDATE country_name SET name = name || NEW.name WHERE country_code = NEW.country_code;
939     END LOOP;
940   END IF;
941
942   -- For linear features we need the full geometry for determining the address
943   -- because they may go through several administrative entities. Otherwise use
944   -- the centroid for performance reasons.
945   IF ST_GeometryType(NEW.geometry) in ('ST_LineString', 'ST_MultiLineString') THEN
946     geom := NEW.geometry;
947   ELSE
948     geom := NEW.centroid;
949   END IF;
950
951   IF NEW.rank_address = 0 THEN
952     max_rank := geometry_to_rank(NEW.rank_search, NEW.geometry, NEW.country_code);
953     -- Rank 0 features may also span multiple administrative areas (e.g. lakes)
954     -- so use the geometry here too. Just make sure the areas don't become too
955     -- large.
956     IF NEW.class = 'natural' or max_rank > 10 THEN
957       geom := NEW.geometry;
958     END IF;
959   ELSEIF NEW.rank_address > 25 THEN
960     max_rank := 25;
961   ELSE
962     max_rank = NEW.rank_address;
963   END IF;
964
965   SELECT * FROM insert_addresslines(NEW.place_id, NEW.partition, max_rank,
966                                     NEW.address, geom, NEW.country_code)
967     INTO NEW.parent_place_id, NEW.postcode, nameaddress_vector;
968
969   {% if debug %}RAISE WARNING 'RETURN insert_addresslines: %, %, %', NEW.parent_place_id, NEW.postcode, nameaddress_vector;{% endif %}
970
971   IF NEW.address is not null AND NEW.address ? 'postcode' 
972      AND NEW.address->'postcode' not similar to '%(,|;)%' THEN
973     NEW.postcode := upper(trim(NEW.address->'postcode'));
974   END IF;
975
976   IF NEW.postcode is null AND NEW.rank_search > 8 THEN
977     NEW.postcode := get_nearest_postcode(NEW.country_code, NEW.geometry);
978   END IF;
979
980   -- if we have a name add this to the name search table
981   IF NEW.name IS NOT NULL THEN
982
983     IF NEW.rank_search <= 25 and NEW.rank_address > 0 THEN
984       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);
985       {% if debug %}RAISE WARNING 'added to location (full)';{% endif %}
986     END IF;
987
988     IF NEW.rank_search between 26 and 27 and NEW.class = 'highway' THEN
989       result := insertLocationRoad(NEW.partition, NEW.place_id, NEW.country_code, NEW.geometry);
990       {% if debug %}RAISE WARNING 'insert into road location table (full)';{% endif %}
991     END IF;
992
993     result := insertSearchName(NEW.partition, NEW.place_id, name_vector,
994                                NEW.rank_search, NEW.rank_address, NEW.geometry);
995     {% if debug %}RAISE WARNING 'added to search name (full)';{% endif %}
996
997     {% if not db.reverse_only %}
998         INSERT INTO search_name (place_id, search_rank, address_rank,
999                                  importance, country_code, name_vector,
1000                                  nameaddress_vector, centroid)
1001                VALUES (NEW.place_id, NEW.rank_search, NEW.rank_address,
1002                        NEW.importance, NEW.country_code, name_vector,
1003                        nameaddress_vector, NEW.centroid);
1004     {% endif %}
1005
1006   END IF;
1007
1008   {% if debug %}RAISE WARNING 'place update % % finsihed.', NEW.osm_type, NEW.osm_id;{% endif %}
1009
1010   RETURN NEW;
1011 END;
1012 $$
1013 LANGUAGE plpgsql;
1014
1015
1016 CREATE OR REPLACE FUNCTION placex_delete()
1017   RETURNS TRIGGER
1018   AS $$
1019 DECLARE
1020   b BOOLEAN;
1021   classtable TEXT;
1022 BEGIN
1023   -- RAISE WARNING 'placex_delete % %',OLD.osm_type,OLD.osm_id;
1024
1025   IF OLD.linked_place_id is null THEN
1026     update placex set linked_place_id = null, indexed_status = 2 where linked_place_id = OLD.place_id and indexed_status = 0;
1027     {% if debug %}RAISE WARNING 'placex_delete:01 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1028     update placex set linked_place_id = null where linked_place_id = OLD.place_id;
1029     {% if debug %}RAISE WARNING 'placex_delete:02 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1030   ELSE
1031     update placex set indexed_status = 2 where place_id = OLD.linked_place_id and indexed_status = 0;
1032   END IF;
1033
1034   IF OLD.rank_address < 30 THEN
1035
1036     -- mark everything linked to this place for re-indexing
1037     {% if debug %}RAISE WARNING 'placex_delete:03 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1038     UPDATE placex set indexed_status = 2 from place_addressline where address_place_id = OLD.place_id 
1039       and placex.place_id = place_addressline.place_id and indexed_status = 0 and place_addressline.isaddress;
1040
1041     {% if debug %}RAISE WARNING 'placex_delete:04 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1042     DELETE FROM place_addressline where address_place_id = OLD.place_id;
1043
1044     {% if debug %}RAISE WARNING 'placex_delete:05 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1045     b := deleteRoad(OLD.partition, OLD.place_id);
1046
1047     {% if debug %}RAISE WARNING 'placex_delete:06 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1048     update placex set indexed_status = 2 where parent_place_id = OLD.place_id and indexed_status = 0;
1049     {% if debug %}RAISE WARNING 'placex_delete:07 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1050     -- reparenting also for OSM Interpolation Lines (and for Tiger?)
1051     update location_property_osmline set indexed_status = 2 where indexed_status = 0 and parent_place_id = OLD.place_id;
1052
1053   END IF;
1054
1055   {% if debug %}RAISE WARNING 'placex_delete:08 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1056
1057   IF OLD.rank_address < 26 THEN
1058     b := deleteLocationArea(OLD.partition, OLD.place_id, OLD.rank_search);
1059   END IF;
1060
1061   {% if debug %}RAISE WARNING 'placex_delete:09 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1062
1063   IF OLD.name is not null THEN
1064     {% if 'search_name' in db.tables %}
1065       DELETE from search_name WHERE place_id = OLD.place_id;
1066     {% endif %}
1067     b := deleteSearchName(OLD.partition, OLD.place_id);
1068   END IF;
1069
1070   {% if debug %}RAISE WARNING 'placex_delete:10 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1071
1072   DELETE FROM place_addressline where place_id = OLD.place_id;
1073
1074   {% if debug %}RAISE WARNING 'placex_delete:11 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1075
1076   -- remove from tables for special search
1077   classtable := 'place_classtype_' || OLD.class || '_' || OLD.type;
1078   SELECT count(*)>0 FROM pg_tables WHERE tablename = classtable and schemaname = current_schema() INTO b;
1079   IF b THEN
1080     EXECUTE 'DELETE FROM ' || classtable::regclass || ' WHERE place_id = $1' USING OLD.place_id;
1081   END IF;
1082
1083   {% if debug %}RAISE WARNING 'placex_delete:12 % %',OLD.osm_type,OLD.osm_id;{% endif %}
1084
1085   RETURN OLD;
1086
1087 END;
1088 $$
1089 LANGUAGE plpgsql;