-- Use case example an address at the intersection of 2 streets:
-- SELECT pprint_addy(addy), st_astext(geomout),rating FROM geocode_intersection('School St', 'Washington St', 'MA', 'Boston','02117');
--DROP FUNCTION tiger.geocode_intersection(text,text,text,text,text,integer);
-CREATE OR REPLACE FUNCTION tiger.geocode_intersection(IN roadway1 text, IN roadway2 text, IN in_state text, IN in_city text DEFAULT '', IN in_zip text DEFAULT '',
+CREATE OR REPLACE FUNCTION geocode_intersection(IN roadway1 text, IN roadway2 text, IN in_state text, IN in_city text DEFAULT '', IN in_zip text DEFAULT '',
IN num_results integer DEFAULT 10, OUT ADDY NORM_ADDY,
OUT GEOMOUT GEOMETRY,
- OUT RATING INTEGER) RETURNS SETOF record AS
+ OUT RATING INTEGER) RETURNS SETOF record AS
$$
DECLARE
var_na_road norm_addy;
var_zip varchar(5)[];
in_statefp varchar(2) ;
var_debug boolean := false;
+ results record;
BEGIN
IF COALESCE(roadway1,'') = '' OR COALESCE(roadway2,'') = '' THEN
-- not enough to give a result just return
RETURN ;
ELSE
- var_na_road := normalize_address(roadway1 || ', ' || COALESCE(in_city,'') || ', ' || in_state || ' ' || in_zip);
- var_na_inter1 := normalize_address(roadway2 || ', ' || COALESCE(in_city,'') || ', ' || in_state || ' ' || in_zip);
+ var_na_road := normalize_address('0 ' || roadway1 || ', ' || COALESCE(in_city,'') || ', ' || in_state || ' ' || in_zip);
+ var_na_inter1 := normalize_address('0 ' || roadway2 || ', ' || COALESCE(in_city,'') || ', ' || in_state || ' ' || in_zip);
END IF;
in_statefp := statefp FROM state_lookup As s WHERE s.abbrev = in_state;
IF COALESCE(in_zip,'') > '' THEN -- limit search to 2 plus or minus the input zip
END IF;
IF var_zip IS NULL AND in_city > '' THEN
- var_zip := array_agg(zip ORDER BY zip) FROM zip_lookup_base WHERE statefp = in_statefp AND lower(city) = lower(in_city);
+ var_zip := array_agg(zip) FROM zip_lookup_base WHERE statefp = in_statefp AND lower(city) = lower(in_city);
END IF;
IF var_zip IS NULL THEN
- var_zip := array_agg(zip ORDER BY zip) FROM zip_lookup_base WHERE statefp = in_statefp AND lower(city) LIKE lower(in_city) || '%' ;
+ var_zip := array_agg(zip) FROM zip_lookup_base WHERE statefp = in_statefp AND lower(city) LIKE lower(in_city) || '%' ;
END IF;
IF var_debug THEN
RAISE NOTICE 'var_zip: %, city: %', quote_nullable(var_zip), quote_nullable(in_city);
END IF;
var_sql := '
- WITH a AS (SELECT * FROM addr WHERE statefp = $1 AND zip = ANY($5::text[]) ),
- e1 AS (SELECT e.the_geom, f.*, a.fromhn, a.tohn,
- CASE WHEN a.side = ''L'' THEN e.tfidl ELSE e.tfidr END AS tfid, a.zip, a.side
- FROM a
+ WITH
+ f1 AS (SELECT * FROM featnames
+ WHERE statefp = $1 AND ( lower(name) = $2 ' ||
+ CASE WHEN length(var_na_road.streetName) > 5 THEN ' or lower(fullname) LIKE $6 || ''%'' ' ELSE '' END || ')'
+ || ') ,
+ f2 AS (SELECT * FROM featnames
+ WHERE statefp = $1 AND ( lower(name) = $4 ' ||
+ CASE WHEN length(var_na_inter1.streetName) > 5 THEN ' or lower(fullname) LIKE $7 || ''%'' ' ELSE '' END || ')'
+ || ' ) ,
+ a1 AS (SELECT f.*, addr.fromhn, addr.tohn, addr.side , addr.zip
+ FROM addr INNER JOIN f1 AS f ON addr.tlid = f.tlid
+ WHERE addr.statefp = $1 AND addr.zip = ANY($5::text[])
+ ),
+ a2 AS (SELECT f.*, addr.fromhn, addr.tohn, addr.side , addr.zip
+ FROM addr INNER JOIN f2 AS f ON addr.tlid = f.tlid
+ WHERE addr.statefp = $1 AND addr.zip = ANY($5::text[])
+ ),
+ e1 AS (SELECT e.the_geom, e.tnidf, e.tnidt, a.*,
+ CASE WHEN a.side = ''L'' THEN e.tfidl ELSE e.tfidr END AS tfid
+ FROM a1 As a
INNER JOIN edges AS e ON a.tlid = e.tlid
- INNER JOIN
- featnames AS f ON (e.tlid = f.tlid AND f.statefp = $1 AND e.statefp = f.statefp)
- WHERE e.statefp = $1 AND (lower(e.fullname) = $6 or lower(f.name) LIKE $2 || ''%'')
- ),
- e2 AS (SELECT e.the_geom, f.* , a.fromhn, a.tohn,
- CASE WHEN a.side = ''L'' THEN e.tfidl ELSE e.tfidr END AS tfid, a.zip, a.side
- FROM e1 INNER JOIN edges AS e ON (e.statefp = $1 AND ST_Intersects(e.the_geom, e1.the_geom) )
- INNER JOIN a ON a.tlid = e.tlid
- INNER JOIN
- (SELECT * FROM featnames WHERE statefp = $1) AS f ON (e.tlid = f.tlid)
- WHERE e.statefp = $1 AND (lower(e.fullname) = $7 or lower(f.name) LIKE $4 || ''%'')
- )
- SELECT ROW(address, predirabrv , fename, fetype,fedirs,NULL,place, $8,zip, true)::norm_addy, geomout, rating
- FROM (SELECT DISTINCT ON(e1.tlid, e1.side)
- CASE WHEN ST_Intersects(ST_StartPoint(e1.the_geom), e2.the_geom) THEN
+ WHERE e.statefp = $1 ) ,
+ e2 AS (SELECT e.the_geom, e.tnidf, e.tnidt, a.*,
+ CASE WHEN a.side = ''L'' THEN e.tfidl ELSE e.tfidr END AS tfid
+ FROM e1 INNER JOIN edges AS e ON (e.statefp = $1 AND ST_Intersects(e.the_geom, e1.the_geom)
+ AND ARRAY[e.tnidf, e.tnidt] && ARRAY[e1.tnidf, e1.tnidt] )
+ INNER JOIN a2 AS a ON a.tlid = e.tlid
+ WHERE e.statefp = $1 AND (lower(e.fullname) = $7 or lower(a.name) LIKE $4 || ''%'')
+ ),
+ segs AS (SELECT DISTINCT ON(e1.tlid, e1.side)
+ CASE WHEN e1.tnidf = e2.tnidf OR e1.tnidf = e2.tnidt THEN
e1.fromhn
- WHEN ST_Intersects(ST_EndPoint(e1.the_geom), e2.the_geom) THEN
- e1.tohn END As address, e1.predirabrv , COALESCE(e1.prequalabr || '' '','''' ) || e1.name As fename,
+ ELSE
+ e1.tohn END As address, e1.predirabrv As fedirp, COALESCE(e1.prequalabr || '' '','''' ) || e1.name As fename,
COALESCE(e1.suftypabrv,e1.pretypabrv) As fetype, e1.sufdirabrv AS fedirs,
p.name As place, e1.zip,
CASE WHEN ST_Intersects(ST_StartPoint(e1.the_geom), e2.the_geom) THEN
ST_StartPoint(e1.the_geom)
WHEN ST_Intersects(ST_EndPoint(e1.the_geom), e2.the_geom) THEN
ST_EndPoint(e1.the_geom)
- ELSE ST_EndPoint(e1.the_geom) END AS geomout ,
+ ELSE ST_EndPoint(e1.the_geom) END AS geom ,
CASE WHEN lower(p.name) = $3 THEN 0 ELSE 1 END
+ levenshtein_ignore_case(e1.name || COALESCE('' '' || e1.sufqualabr, ''''),$2) +
CASE WHEN e1.fullname = $6 THEN 0 ELSE levenshtein_ignore_case(e1.fullname, $6) END +
+ levenshtein_ignore_case(e2.name || COALESCE('' '' || e2.sufqualabr, ''''),$4)
- AS rating
+ + CASE WHEN e2.tfid = e1.tfid THEN 0 ELSE 3 END
+ AS a_rating
FROM e1
- INNER JOIN (SELECT * FROM faces WHERE statefp = $1) As fa1 ON (e1.tfid = fa1.tfid)
INNER JOIN e2 ON (
ST_Intersects(e1.the_geom, e2.the_geom) )
+ INNER JOIN (SELECT * FROM faces WHERE statefp = $1) As fa1 ON (e1.tfid = fa1.tfid )
LEFT JOIN place AS p ON (fa1.placefp = p.placefp AND p.statefp = $1 )
- ORDER BY e1.tlid, e1.side, rating LIMIT 20) As foo ORDER BY rating LIMIT
-
- ' || num_results::text;
-
- /**var_sql := replace(replace(
- replace(replace(replace(
- replace(
- replace(
- replace(var_sql, '$1', quote_nullable(in_statefp)),
- '$2', quote_nullable(lower(var_na_road.streetName) ) ),
- '$3', quote_nullable(lower(in_city)) ),
- '$4', quote_nullable(lower(var_na_inter1.streetName) ) ),
- '$5', quote_nullable(var_zip) ),
- '$6', quote_nullable(lower(var_na_road.streetName || ' ' || COALESCE(var_na_road.streetTypeAbbrev,'') )) ) ,
- '$7', quote_nullable(lower(var_na_inter1.streetName || ' ' || COALESCE(var_na_inter1.streetTypeAbbrev,'') ) ,
- '$8', quote_nullable(in_state ) ) ); **/
+ ORDER BY e1.tlid, e1.side, a_rating LIMIT $9*4 )
+ SELECT address, fedirp , fename, fetype,fedirs,place, zip , geom, a_rating
+ FROM segs ORDER BY a_rating LIMIT $9';
IF var_debug THEN
- RAISE NOTICE 'sql: %', replace(replace(
- replace(replace(replace(
+ RAISE NOTICE 'sql: %', replace(replace(replace(
+ replace(replace(replace(
replace(
replace(
replace(var_sql, '$1', quote_nullable(in_statefp)),
'$5', quote_nullable(var_zip) ),
'$6', quote_nullable(lower(var_na_road.streetName || ' ' || COALESCE(var_na_road.streetTypeAbbrev,'') )) ) ,
'$7', quote_nullable(lower(var_na_inter1.streetName || ' ' || COALESCE(var_na_inter1.streetTypeAbbrev,'') ) ) ) ,
- '$8', quote_nullable(in_state ) );
+ '$8', quote_nullable(in_state ) ), '$9', num_results::text );
END IF;
- RETURN QUERY EXECUTE var_sql USING in_statefp, lower(var_na_road.streetName), lower(in_city), lower(var_na_inter1.streetName), var_zip,
+ FOR results IN EXECUTE var_sql USING in_statefp, lower(var_na_road.streetName), lower(in_city), lower(var_na_inter1.streetName), var_zip,
lower(var_na_road.streetName || ' ' || COALESCE(var_na_road.streetTypeAbbrev,'')),
- lower(var_na_inter1.streetName || ' ' || COALESCE(var_na_inter1.streetTypeAbbrev,'')), in_state ;
+ lower(var_na_inter1.streetName || ' ' || COALESCE(var_na_inter1.streetTypeAbbrev,'')), in_state, num_results LOOP
+ ADDY.preDirAbbrev := results.fedirp;
+ ADDY.streetName := results.fename;
+ ADDY.streetTypeAbbrev := results.fetype;
+ ADDY.postDirAbbrev := results.fedirs;
+ ADDY.location := results.place;
+ ADDY.stateAbbrev := in_state;
+ ADDY.zip := results.zip;
+ ADDY.parsed := TRUE;
+ ADDY.address := results.address;
+
+ GEOMOUT := results.geom;
+ RATING := results.a_rating;
+ RETURN NEXT;
+ END LOOP;
+ RETURN;
END;
$$
LANGUAGE plpgsql IMMUTABLE
COST 1000
ROWS 10;
-
-
projects / postgis / commitdiff