]> git.openstreetmap.org Git - nominatim.git/blob - lib-sql/functions/utils.sql
update Ubuntu22 install script for pip install
[nominatim.git] / lib-sql / functions / utils.sql
1 -- SPDX-License-Identifier: GPL-2.0-only
2 --
3 -- This file is part of Nominatim. (https://nominatim.org)
4 --
5 -- Copyright (C) 2022 by the Nominatim developer community.
6 -- For a full list of authors see the git log.
7
8 -- Assorted helper functions for the triggers.
9
10 CREATE OR REPLACE FUNCTION get_center_point(place GEOMETRY)
11   RETURNS GEOMETRY
12   AS $$
13 DECLARE
14   geom_type TEXT;
15 BEGIN
16   geom_type := ST_GeometryType(place);
17   IF geom_type = ' ST_Point' THEN
18     RETURN place;
19   END IF;
20   IF geom_type = 'ST_LineString' THEN
21     RETURN ST_LineInterpolatePoint(place, 0.5);
22   END IF;
23
24   RETURN ST_PointOnSurface(place);
25 END;
26 $$
27 LANGUAGE plpgsql IMMUTABLE;
28
29
30 CREATE OR REPLACE FUNCTION geometry_sector(partition INTEGER, place GEOMETRY)
31   RETURNS INTEGER
32   AS $$
33 BEGIN
34   RETURN (partition*1000000) + (500-ST_X(place)::INTEGER)*1000 + (500-ST_Y(place)::INTEGER);
35 END;
36 $$
37 LANGUAGE plpgsql IMMUTABLE;
38
39
40
41 CREATE OR REPLACE FUNCTION array_merge(a INTEGER[], b INTEGER[])
42   RETURNS INTEGER[]
43   AS $$
44 DECLARE
45   i INTEGER;
46   r INTEGER[];
47 BEGIN
48   IF array_upper(a, 1) IS NULL THEN
49     RETURN b;
50   END IF;
51   IF array_upper(b, 1) IS NULL THEN
52     RETURN a;
53   END IF;
54   r := a;
55   FOR i IN 1..array_upper(b, 1) LOOP  
56     IF NOT (ARRAY[b[i]] <@ r) THEN
57       r := r || b[i];
58     END IF;
59   END LOOP;
60   RETURN r;
61 END;
62 $$
63 LANGUAGE plpgsql IMMUTABLE;
64
65 -- Return the node members with a given label from a relation member list
66 -- as a set.
67 --
68 -- \param members      Member list in osm2pgsql middle format.
69 -- \param memberLabels Array of labels to accept.
70 --
71 -- \returns Set of OSM ids of nodes that are found.
72 --
73 CREATE OR REPLACE FUNCTION get_rel_node_members(members TEXT[],
74                                                 memberLabels TEXT[])
75   RETURNS SETOF BIGINT
76   AS $$
77 DECLARE
78   i INTEGER;
79 BEGIN
80   FOR i IN 1..ARRAY_UPPER(members,1) BY 2 LOOP
81     IF members[i+1] = ANY(memberLabels)
82        AND upper(substring(members[i], 1, 1))::char(1) = 'N'
83     THEN
84       RETURN NEXT substring(members[i], 2)::bigint;
85     END IF;
86   END LOOP;
87
88   RETURN;
89 END;
90 $$
91 LANGUAGE plpgsql IMMUTABLE;
92
93
94 CREATE OR REPLACE FUNCTION get_rel_node_members(members JSONB, memberLabels TEXT[])
95   RETURNS SETOF BIGINT
96   AS $$
97 DECLARE
98   member JSONB;
99 BEGIN
100   FOR member IN SELECT * FROM jsonb_array_elements(members)
101   LOOP
102     IF member->>'type' = 'N' and member->>'role' = ANY(memberLabels) THEN
103         RETURN NEXT (member->>'ref')::bigint;
104     END IF;
105   END LOOP;
106
107   RETURN;
108 END;
109 $$
110 LANGUAGE plpgsql IMMUTABLE;
111
112
113 -- Copy 'name' to or from the default language.
114 --
115 -- \param country_code     Country code of the object being named.
116 -- \param[inout] name      List of names of the object.
117 --
118 -- If the country named by country_code has a single default language,
119 -- then a `name` tag is copied to `name:<country_code>` if this tag does
120 -- not yet exist and vice versa.
121 CREATE OR REPLACE FUNCTION add_default_place_name(country_code VARCHAR(2),
122                                                   INOUT name HSTORE)
123   AS $$
124 DECLARE
125   default_language VARCHAR(10);
126 BEGIN
127   IF name is not null AND array_upper(akeys(name),1) > 1 THEN
128     default_language := get_country_language_code(country_code);
129     IF default_language IS NOT NULL THEN
130       IF name ? 'name' AND NOT name ? ('name:'||default_language) THEN
131         name := name || hstore(('name:'||default_language), (name -> 'name'));
132       ELSEIF name ? ('name:'||default_language) AND NOT name ? 'name' THEN
133         name := name || hstore('name', (name -> ('name:'||default_language)));
134       END IF;
135     END IF;
136   END IF;
137 END;
138 $$
139 LANGUAGE plpgsql IMMUTABLE;
140
141
142 -- Find the nearest artificial postcode for the given geometry.
143 -- TODO For areas there should not be more than two inside the geometry.
144 CREATE OR REPLACE FUNCTION get_nearest_postcode(country VARCHAR(2), geom GEOMETRY)
145   RETURNS TEXT
146   AS $$
147 DECLARE
148   outcode TEXT;
149   cnt INTEGER;
150 BEGIN
151     -- If the geometry is an area then only one postcode must be within
152     -- that area, otherwise consider the area as not having a postcode.
153     IF ST_GeometryType(geom) in ('ST_Polygon','ST_MultiPolygon') THEN
154         SELECT min(postcode), count(*) FROM
155               (SELECT postcode FROM location_postcode
156                 WHERE ST_Contains(geom, location_postcode.geometry) LIMIT 2) sub
157           INTO outcode, cnt;
158
159         IF cnt = 1 THEN
160             RETURN outcode;
161         ELSE
162             RETURN null;
163         END IF;
164     END IF;
165
166     SELECT postcode FROM location_postcode
167      WHERE ST_DWithin(geom, location_postcode.geometry, 0.05)
168           AND location_postcode.country_code = country
169      ORDER BY ST_Distance(geom, location_postcode.geometry) LIMIT 1
170     INTO outcode;
171
172     RETURN outcode;
173 END;
174 $$
175 LANGUAGE plpgsql STABLE;
176
177
178 CREATE OR REPLACE FUNCTION get_country_code(place geometry)
179   RETURNS TEXT
180   AS $$
181 DECLARE
182   nearcountry RECORD;
183   countries TEXT[];
184 BEGIN
185 -- RAISE WARNING 'get_country_code, start: %', ST_AsText(place);
186
187   -- Try for a OSM polygon
188   SELECT array_agg(country_code) FROM location_area_country
189     WHERE country_code is not null and st_covers(geometry, place)
190     INTO countries;
191
192   IF array_length(countries, 1) = 1 THEN
193     RETURN countries[1];
194   END IF;
195
196   IF array_length(countries, 1) > 1 THEN
197     -- more than one country found, confirm against the fallback data what to choose
198     FOR nearcountry IN
199         SELECT country_code FROM country_osm_grid
200           WHERE ST_Covers(geometry, place) AND country_code = ANY(countries)
201           ORDER BY area ASC
202     LOOP
203         RETURN nearcountry.country_code;
204     END LOOP;
205     -- Still nothing? Choose the country code with the smallest partition number.
206     -- And failing that, just go by the alphabet.
207     FOR nearcountry IN
208         SELECT cc,
209                (SELECT partition FROM country_name WHERE country_code = cc) as partition
210         FROM unnest(countries) cc
211         ORDER BY partition, cc
212     LOOP
213         RETURN nearcountry.cc;
214     END LOOP;
215
216     -- Should never be reached.
217     RETURN countries[1];
218   END IF;
219
220 -- RAISE WARNING 'osm fallback: %', ST_AsText(place);
221
222   -- Try for OSM fallback data
223   -- The order is to deal with places like HongKong that are 'states' within another polygon
224   FOR nearcountry IN
225     SELECT country_code from country_osm_grid
226     WHERE st_covers(geometry, place) order by area asc limit 1
227   LOOP
228     RETURN nearcountry.country_code;
229   END LOOP;
230
231 -- RAISE WARNING 'near osm fallback: %', ST_AsText(place);
232
233   RETURN NULL;
234 END;
235 $$
236 LANGUAGE plpgsql STABLE;
237
238
239 CREATE OR REPLACE FUNCTION get_country_language_code(search_country_code VARCHAR(2))
240   RETURNS TEXT
241   AS $$
242 DECLARE
243   nearcountry RECORD;
244 BEGIN
245   FOR nearcountry IN
246     SELECT distinct country_default_language_code from country_name
247     WHERE country_code = search_country_code limit 1
248   LOOP
249     RETURN lower(nearcountry.country_default_language_code);
250   END LOOP;
251   RETURN NULL;
252 END;
253 $$
254 LANGUAGE plpgsql STABLE;
255
256
257 CREATE OR REPLACE FUNCTION get_partition(in_country_code VARCHAR(10))
258   RETURNS INTEGER
259   AS $$
260 DECLARE
261   nearcountry RECORD;
262 BEGIN
263   FOR nearcountry IN
264     SELECT partition from country_name where country_code = in_country_code
265   LOOP
266     RETURN nearcountry.partition;
267   END LOOP;
268   RETURN 0;
269 END;
270 $$
271 LANGUAGE plpgsql STABLE;
272
273
274 -- Find the parent of an address with addr:street/addr:place tag.
275 --
276 -- \param token_info Naming info with the address information.
277 -- \param partition  Partition where to search the parent.
278 -- \param centroid   Location of the address.
279 --
280 -- \return Place ID of the parent if one was found, NULL otherwise.
281 CREATE OR REPLACE FUNCTION find_parent_for_address(token_info JSONB,
282                                                    partition SMALLINT,
283                                                    centroid GEOMETRY)
284   RETURNS BIGINT
285   AS $$
286 DECLARE
287   parent_place_id BIGINT;
288 BEGIN
289   -- Check for addr:street attributes
290   parent_place_id := getNearestNamedRoadPlaceId(partition, centroid, token_info);
291   IF parent_place_id is not null THEN
292     {% if debug %}RAISE WARNING 'Get parent from addr:street: %', parent_place_id;{% endif %}
293     RETURN parent_place_id;
294   END IF;
295
296   -- Check for addr:place attributes.
297   parent_place_id := getNearestNamedPlacePlaceId(partition, centroid, token_info);
298   {% if debug %}RAISE WARNING 'Get parent from addr:place: %', parent_place_id;{% endif %}
299   RETURN parent_place_id;
300 END;
301 $$
302 LANGUAGE plpgsql STABLE;
303
304
305 CREATE OR REPLACE FUNCTION delete_location(OLD_place_id BIGINT)
306   RETURNS BOOLEAN
307   AS $$
308 DECLARE
309 BEGIN
310   DELETE FROM location_area where place_id = OLD_place_id;
311 -- TODO:location_area
312   RETURN true;
313 END;
314 $$
315 LANGUAGE plpgsql;
316
317 -- Create a bounding box with an extent computed from the radius (in meters)
318 -- which in turn is derived from the given search rank.
319 CREATE OR REPLACE FUNCTION place_node_fuzzy_area(geom GEOMETRY, rank_search INTEGER)
320   RETURNS GEOMETRY
321   AS $$
322 DECLARE
323   radius FLOAT := 500;
324 BEGIN
325   IF rank_search <= 16 THEN -- city
326     radius := 15000;
327   ELSIF rank_search <= 18 THEN -- town
328     radius := 4000;
329   ELSIF rank_search <= 19 THEN -- village
330     radius := 2000;
331   ELSIF rank_search  <= 20 THEN -- hamlet
332     radius := 1000;
333   END IF;
334
335   RETURN ST_Envelope(ST_Collect(
336                      ST_Project(geom::geography, radius, 0.785398)::geometry,
337                      ST_Project(geom::geography, radius, 3.9269908)::geometry));
338 END;
339 $$
340 LANGUAGE plpgsql IMMUTABLE;
341
342
343 CREATE OR REPLACE FUNCTION add_location(place_id BIGINT, country_code varchar(2),
344                                         partition INTEGER, keywords INTEGER[],
345                                         rank_search INTEGER, rank_address INTEGER,
346                                         in_postcode TEXT, geometry GEOMETRY,
347                                         centroid GEOMETRY)
348   RETURNS BOOLEAN
349   AS $$
350 DECLARE
351   locationid INTEGER;
352   secgeo GEOMETRY;
353   postcode TEXT;
354 BEGIN
355   PERFORM deleteLocationArea(partition, place_id, rank_search);
356
357   -- add postcode only if it contains a single entry, i.e. ignore postcode lists
358   postcode := NULL;
359   IF in_postcode is not null AND in_postcode not similar to '%(,|;)%' THEN
360       postcode := upper(trim (in_postcode));
361   END IF;
362
363   IF ST_GeometryType(geometry) in ('ST_Polygon','ST_MultiPolygon') THEN
364     FOR secgeo IN select split_geometry(geometry) AS geom LOOP
365       PERFORM insertLocationAreaLarge(partition, place_id, country_code, keywords, rank_search, rank_address, false, postcode, centroid, secgeo);
366     END LOOP;
367
368   ELSEIF ST_GeometryType(geometry) = 'ST_Point' THEN
369     secgeo := place_node_fuzzy_area(geometry, rank_search);
370     PERFORM insertLocationAreaLarge(partition, place_id, country_code, keywords, rank_search, rank_address, true, postcode, centroid, secgeo);
371
372   END IF;
373
374   RETURN true;
375 END;
376 $$
377 LANGUAGE plpgsql;
378
379
380 CREATE OR REPLACE FUNCTION quad_split_geometry(geometry GEOMETRY, maxarea FLOAT,
381                                                maxdepth INTEGER)
382   RETURNS SETOF GEOMETRY
383   AS $$
384 DECLARE
385   xmin FLOAT;
386   ymin FLOAT;
387   xmax FLOAT;
388   ymax FLOAT;
389   xmid FLOAT;
390   ymid FLOAT;
391   secgeo GEOMETRY;
392   secbox GEOMETRY;
393   seg INTEGER;
394   geo RECORD;
395   area FLOAT;
396   remainingdepth INTEGER;
397   added INTEGER;
398 BEGIN
399
400 --  RAISE WARNING 'quad_split_geometry: maxarea=%, depth=%',maxarea,maxdepth;
401
402   IF (ST_GeometryType(geometry) not in ('ST_Polygon','ST_MultiPolygon') OR NOT ST_IsValid(geometry)) THEN
403     RETURN NEXT geometry;
404     RETURN;
405   END IF;
406
407   remainingdepth := maxdepth - 1;
408   area := ST_AREA(geometry);
409   IF remainingdepth < 1 OR area < maxarea THEN
410     RETURN NEXT geometry;
411     RETURN;
412   END IF;
413
414   xmin := st_xmin(geometry);
415   xmax := st_xmax(geometry);
416   ymin := st_ymin(geometry);
417   ymax := st_ymax(geometry);
418   secbox := ST_SetSRID(ST_MakeBox2D(ST_Point(ymin,xmin),ST_Point(ymax,xmax)),4326);
419
420   -- if the geometry completely covers the box don't bother to slice any more
421   IF ST_AREA(secbox) = area THEN
422     RETURN NEXT geometry;
423     RETURN;
424   END IF;
425
426   xmid := (xmin+xmax)/2;
427   ymid := (ymin+ymax)/2;
428
429   added := 0;
430   FOR seg IN 1..4 LOOP
431
432     IF seg = 1 THEN
433       secbox := ST_SetSRID(ST_MakeBox2D(ST_Point(xmin,ymin),ST_Point(xmid,ymid)),4326);
434     END IF;
435     IF seg = 2 THEN
436       secbox := ST_SetSRID(ST_MakeBox2D(ST_Point(xmin,ymid),ST_Point(xmid,ymax)),4326);
437     END IF;
438     IF seg = 3 THEN
439       secbox := ST_SetSRID(ST_MakeBox2D(ST_Point(xmid,ymin),ST_Point(xmax,ymid)),4326);
440     END IF;
441     IF seg = 4 THEN
442       secbox := ST_SetSRID(ST_MakeBox2D(ST_Point(xmid,ymid),ST_Point(xmax,ymax)),4326);
443     END IF;
444
445     IF st_intersects(geometry, secbox) THEN
446       secgeo := st_intersection(geometry, secbox);
447       IF NOT ST_IsEmpty(secgeo) AND ST_GeometryType(secgeo) in ('ST_Polygon','ST_MultiPolygon') THEN
448         FOR geo IN select quad_split_geometry(secgeo, maxarea, remainingdepth) as geom LOOP
449           IF NOT ST_IsEmpty(geo.geom) AND ST_GeometryType(geo.geom) in ('ST_Polygon','ST_MultiPolygon') THEN
450             added := added + 1;
451             RETURN NEXT geo.geom;
452           END IF;
453         END LOOP;
454       END IF;
455     END IF;
456   END LOOP;
457
458   RETURN;
459 END;
460 $$
461 LANGUAGE plpgsql IMMUTABLE;
462
463
464 CREATE OR REPLACE FUNCTION split_geometry(geometry GEOMETRY)
465   RETURNS SETOF GEOMETRY
466   AS $$
467 DECLARE
468   geo RECORD;
469 BEGIN
470   -- 10000000000 is ~~ 1x1 degree
471   FOR geo IN select quad_split_geometry(geometry, 0.25, 20) as geom LOOP
472     RETURN NEXT geo.geom;
473   END LOOP;
474   RETURN;
475 END;
476 $$
477 LANGUAGE plpgsql IMMUTABLE;
478
479 CREATE OR REPLACE FUNCTION simplify_large_polygons(geometry GEOMETRY)
480   RETURNS GEOMETRY
481   AS $$
482 BEGIN
483   IF ST_GeometryType(geometry) in ('ST_Polygon','ST_MultiPolygon')
484      and ST_MemSize(geometry) > 3000000
485   THEN
486     geometry := ST_SimplifyPreserveTopology(geometry, 0.0001);
487   END IF;
488   RETURN geometry;
489 END;
490 $$
491 LANGUAGE plpgsql IMMUTABLE;
492
493
494 CREATE OR REPLACE FUNCTION place_force_delete(placeid BIGINT)
495   RETURNS BOOLEAN
496   AS $$
497 DECLARE
498     osmid BIGINT;
499     osmtype character(1);
500     pclass text;
501     ptype text;
502 BEGIN
503   SELECT osm_type, osm_id, class, type FROM placex WHERE place_id = placeid INTO osmtype, osmid, pclass, ptype;
504   DELETE FROM import_polygon_delete where osm_type = osmtype and osm_id = osmid and class = pclass and type = ptype;
505   DELETE FROM import_polygon_error where osm_type = osmtype and osm_id = osmid and class = pclass and type = ptype;
506   -- force delete by directly entering it into the to-be-deleted table
507   INSERT INTO place_to_be_deleted (osm_type, osm_id, class, type, deferred)
508          VALUES(osmtype, osmid, pclass, ptype, false);
509   PERFORM flush_deleted_places();
510
511   RETURN TRUE;
512 END;
513 $$
514 LANGUAGE plpgsql;
515
516
517 CREATE OR REPLACE FUNCTION place_force_update(placeid BIGINT)
518   RETURNS BOOLEAN
519   AS $$
520 DECLARE
521   placegeom GEOMETRY;
522   geom GEOMETRY;
523   diameter FLOAT;
524   rank SMALLINT;
525 BEGIN
526   UPDATE placex SET indexed_status = 2 WHERE place_id = placeid;
527
528   SELECT geometry, rank_address INTO placegeom, rank
529     FROM placex WHERE place_id = placeid;
530
531   IF placegeom IS NOT NULL AND ST_IsValid(placegeom) THEN
532     IF ST_GeometryType(placegeom) in ('ST_Polygon','ST_MultiPolygon')
533        AND rank > 0
534     THEN
535       FOR geom IN SELECT split_geometry(placegeom) LOOP
536         UPDATE placex SET indexed_status = 2
537          WHERE ST_Intersects(geom, placex.geometry)
538                and indexed_status = 0
539                and ((rank_address = 0 and rank_search > rank) or rank_address > rank)
540                and (rank_search < 28 or name is not null or (rank >= 16 and address ? 'place'));
541       END LOOP;
542     ELSE
543         diameter := update_place_diameter(rank);
544         IF diameter > 0 THEN
545           IF rank >= 26 THEN
546             -- roads may cause reparenting for >27 rank places
547             update placex set indexed_status = 2 where indexed_status = 0 and rank_search > rank and ST_DWithin(placex.geometry, placegeom, diameter);
548           ELSEIF rank >= 16 THEN
549             -- up to rank 16, street-less addresses may need reparenting
550             update placex set indexed_status = 2 where indexed_status = 0 and rank_search > rank and ST_DWithin(placex.geometry, placegeom, diameter) and (rank_search < 28 or name is not null or address ? 'place');
551           ELSE
552             -- for all other places the search terms may change as well
553             update placex set indexed_status = 2 where indexed_status = 0 and rank_search > rank and ST_DWithin(placex.geometry, placegeom, diameter) and (rank_search < 28 or name is not null);
554           END IF;
555         END IF;
556     END IF;
557     RETURN TRUE;
558   END IF;
559
560   RETURN FALSE;
561 END;
562 $$
563 LANGUAGE plpgsql;
564
565 CREATE OR REPLACE FUNCTION flush_deleted_places()
566   RETURNS INTEGER
567   AS $$
568 BEGIN
569   -- deleting large polygons can have a massive effect on the system - require manual intervention to let them through
570   INSERT INTO import_polygon_delete (osm_type, osm_id, class, type)
571     SELECT osm_type, osm_id, class, type FROM place_to_be_deleted WHERE deferred;
572
573   -- delete from place table
574   ALTER TABLE place DISABLE TRIGGER place_before_delete;
575   DELETE FROM place USING place_to_be_deleted
576     WHERE place.osm_type = place_to_be_deleted.osm_type
577           and place.osm_id = place_to_be_deleted.osm_id
578           and place.class = place_to_be_deleted.class
579           and place.type = place_to_be_deleted.type
580           and not deferred;
581   ALTER TABLE place ENABLE TRIGGER place_before_delete;
582
583   -- Mark for delete in the placex table
584   UPDATE placex SET indexed_status = 100 FROM place_to_be_deleted
585     WHERE placex.osm_type = 'N' and place_to_be_deleted.osm_type = 'N'
586           and placex.osm_id = place_to_be_deleted.osm_id
587           and placex.class = place_to_be_deleted.class
588           and placex.type = place_to_be_deleted.type
589           and not deferred;
590   UPDATE placex SET indexed_status = 100 FROM place_to_be_deleted
591     WHERE placex.osm_type = 'W' and place_to_be_deleted.osm_type = 'W'
592           and placex.osm_id = place_to_be_deleted.osm_id
593           and placex.class = place_to_be_deleted.class
594           and placex.type = place_to_be_deleted.type
595           and not deferred;
596   UPDATE placex SET indexed_status = 100 FROM place_to_be_deleted
597     WHERE placex.osm_type = 'R' and place_to_be_deleted.osm_type = 'R'
598           and placex.osm_id = place_to_be_deleted.osm_id
599           and placex.class = place_to_be_deleted.class
600           and placex.type = place_to_be_deleted.type
601           and not deferred;
602
603    -- Mark for delete in interpolations
604    UPDATE location_property_osmline SET indexed_status = 100 FROM place_to_be_deleted
605     WHERE place_to_be_deleted.osm_type = 'W'
606           and place_to_be_deleted.class = 'place'
607           and place_to_be_deleted.type = 'houses'
608           and location_property_osmline.osm_id = place_to_be_deleted.osm_id
609           and not deferred;
610
611    -- Clear todo list.
612    TRUNCATE TABLE place_to_be_deleted;
613
614    RETURN NULL;
615 END;
616 $$ LANGUAGE plpgsql;