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