]> git.openstreetmap.org Git - nominatim.git/blob - lib-sql/tiger_import_start.sql
Merge pull request #2562 from lonvia/copyright-headers
[nominatim.git] / lib-sql / tiger_import_start.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 DROP TABLE IF EXISTS location_property_tiger_import;
8 CREATE TABLE location_property_tiger_import (linegeo GEOMETRY, place_id BIGINT, partition INTEGER, parent_place_id BIGINT, startnumber INTEGER, endnumber INTEGER, interpolationtype TEXT, postcode TEXT);
9
10 CREATE OR REPLACE FUNCTION tiger_line_import(linegeo GEOMETRY, in_startnumber INTEGER,
11                                              in_endnumber INTEGER, interpolationtype TEXT,
12                                              token_info JSONB, in_postcode TEXT) RETURNS INTEGER
13   AS $$
14 DECLARE
15   startnumber INTEGER;
16   endnumber INTEGER;
17   stepsize INTEGER;
18   numberrange INTEGER;
19   place_centroid GEOMETRY;
20   out_partition INTEGER;
21   out_parent_place_id BIGINT;
22   location RECORD;
23
24 BEGIN
25
26   IF in_endnumber > in_startnumber THEN
27     startnumber = in_startnumber;
28     endnumber = in_endnumber;
29   ELSE
30     startnumber = in_endnumber;
31     endnumber = in_startnumber;
32   END IF;
33
34   IF startnumber < 0 THEN
35     RAISE WARNING 'Negative house number range (% to %)', startnumber, endnumber;
36     RETURN 0;
37   END IF;
38
39   numberrange := endnumber - startnumber;
40
41   IF (interpolationtype = 'odd' AND startnumber % 2 = 0) OR (interpolationtype = 'even' AND startnumber % 2 = 1) THEN
42     startnumber := startnumber + 1;
43     stepsize := 2;
44   ELSE
45     IF (interpolationtype = 'odd' OR interpolationtype = 'even') THEN
46       stepsize := 2;
47     ELSE -- everything else assumed to be 'all'
48       stepsize := 1;
49     END IF;
50   END IF;
51
52   -- Filter out really broken tiger data
53   IF numberrange > 0 AND (numberrange::float/stepsize::float > 500)
54                     AND ST_length(linegeo)/(numberrange::float/stepsize::float) < 0.000001 THEN
55     RAISE WARNING 'Road too short for number range % to % (%)',startnumber,endnumber,
56                   ST_length(linegeo)/(numberrange::float/stepsize::float);
57     RETURN 0;
58   END IF;
59
60   place_centroid := ST_Centroid(linegeo);
61   out_partition := get_partition('us');
62
63   out_parent_place_id := getNearestNamedRoadPlaceId(out_partition, place_centroid,
64                                                     token_info);
65
66   IF out_parent_place_id IS NULL THEN
67     SELECT getNearestParallelRoadFeature(out_partition, linegeo)
68       INTO out_parent_place_id;
69   END IF;
70
71   IF out_parent_place_id IS NULL THEN
72     SELECT getNearestRoadPlaceId(out_partition, place_centroid)
73       INTO out_parent_place_id;
74   END IF;
75
76 --insert street(line) into import table
77 insert into location_property_tiger_import (linegeo, place_id, partition, parent_place_id, startnumber, endnumber, interpolationtype, postcode)
78 values (linegeo, nextval('seq_place'), out_partition, out_parent_place_id, startnumber, endnumber, interpolationtype, in_postcode);
79
80   RETURN 1;
81 END;
82 $$
83 LANGUAGE plpgsql;