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