From: Sarah Hoffmann Date: Thu, 30 Apr 2015 19:48:25 +0000 (+0200) Subject: Merge branch 'reverse-lookup-including-tiger-lookup' of https://github.com/mtmail... X-Git-Tag: v.2.5.0~60 X-Git-Url: https://git.openstreetmap.org/nominatim.git/commitdiff_plain/cf36e8d2fba81351889208e4499e8d28e7d9593c?hp=cc4070618851aa4194f6263cdf3de63fcdf8074d Merge branch 'reverse-lookup-including-tiger-lookup' of https://github.com/mtmail/Nominatim --- diff --git a/AUTHORS b/AUTHORS index ea48a49c..b90e7057 100644 --- a/AUTHORS +++ b/AUTHORS @@ -15,3 +15,5 @@ Nominatim was written by: Kurt Roeckx Rodolphe Quiédeville Marc Tobias Metten + +Reverse geocoding using Tiger data feature made possible with support from Guru Labs diff --git a/ChangeLog b/ChangeLog index f8c4bb69..d6002d43 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,3 +1,6 @@ + * reverse geocoding looking includes looking up Tiger data + + 2.4 * drop support for postgres 8.4 diff --git a/lib/PlaceLookup.php b/lib/PlaceLookup.php index 6da82fd8..84c68a99 100644 --- a/lib/PlaceLookup.php +++ b/lib/PlaceLookup.php @@ -5,6 +5,8 @@ protected $iPlaceID; + protected $bIsTiger = false; + protected $aLangPrefOrder = array(); protected $bAddressDetails = false; @@ -35,6 +37,11 @@ $this->iPlaceID = $this->oDB->getOne($sSQL); } + function setIsTiger($b = false) + { + $this->bIsTiger = $b; + } + function lookup() { if (!$this->iPlaceID) return null; @@ -49,12 +56,29 @@ $sSQL .= " (case when centroid is null then st_y(st_centroid(geometry)) else st_y(centroid) end) as lat,"; $sSQL .= " (case when centroid is null then st_x(st_centroid(geometry)) else st_x(centroid) end) as lon"; $sSQL .= " from placex where place_id = ".(int)$this->iPlaceID; + + + if ($this->bIsTiger) + { + $sSQL = "select place_id,partition, 'T' as osm_type, place_id as osm_id, 'place' as class, 'house' as type, null as admin_level, housenumber, null as street, null as isin, postcode,"; + $sSQL .= " 'us' as country_code, null as extratags, parent_place_id, null as linked_place_id, 30 as rank_address, 30 as rank_search,"; + $sSQL .= " coalesce(null,0.75-(30::float/40)) as importance, null as indexed_status, null as indexed_date, null as wikipedia, 'us' as calculated_country_code, "; + $sSQL .= " get_address_by_language(place_id, $sLanguagePrefArraySQL) as langaddress,"; + $sSQL .= " null as placename,"; + $sSQL .= " null as ref,"; + $sSQL .= " st_y(centroid) as lat,"; + $sSQL .= " st_x(centroid) as lon"; + $sSQL .= " from location_property_tiger where place_id = ".(int)$this->iPlaceID; + } + $aPlace = $this->oDB->getRow($sSQL); + if (PEAR::IsError($aPlace)) { failInternalError("Could not lookup place.", $sSQL, $aPlace); } + if (!$aPlace['place_id']) return null; if ($this->bAddressDetails) @@ -63,6 +87,7 @@ $aPlace['aAddress'] = $aAddress; } + $aClassType = getClassTypes(); $sAddressType = ''; $sClassType = $aPlace['class'].':'.$aPlace['type'].':'.$aPlace['admin_level']; diff --git a/lib/ReverseGeocode.php b/lib/ReverseGeocode.php index cb8c0c6b..ae83af1e 100644 --- a/lib/ReverseGeocode.php +++ b/lib/ReverseGeocode.php @@ -69,12 +69,15 @@ { $sPointSQL = 'ST_SetSRID(ST_Point('.$this->fLon.','.$this->fLat.'),4326)'; $iMaxRank = $this->iMaxRank; + $iMaxRank_orig = $this->iMaxRank; // Find the nearest point $fSearchDiam = 0.0004; $iPlaceID = null; $aArea = false; $fMaxAreaDistance = 1; + $bIsInUnitedStates = false; + $bPlaceIsTiger = false; while(!$iPlaceID && $fSearchDiam < $fMaxAreaDistance) { $fSearchDiam = $fSearchDiam * 2; @@ -90,7 +93,7 @@ if ($fSearchDiam > 0.008 && $iMaxRank > 22) $iMaxRank = 22; if ($fSearchDiam > 0.001 && $iMaxRank > 26) $iMaxRank = 26; - $sSQL = 'select place_id,parent_place_id,rank_search from placex'; + $sSQL = 'select place_id,parent_place_id,rank_search,calculated_country_code from placex'; $sSQL .= ' WHERE ST_DWithin('.$sPointSQL.', geometry, '.$fSearchDiam.')'; $sSQL .= ' and rank_search != 28 and rank_search >= '.$iMaxRank; $sSQL .= ' and (name is not null or housenumber is not null)'; @@ -107,12 +110,52 @@ } $iPlaceID = $aPlace['place_id']; $iParentPlaceID = $aPlace['parent_place_id']; + $bIsInUnitedStates = ($aPlace['calculated_country_code'] == 'us'); + } + + // Only street found? If it's in the US we can check TIGER data for nearest housenumber + if ($bIsInUnitedStates && $iMaxRank_orig >= 28 && $iPlaceID && ($aPlace['rank_search'] == 26 || $aPlace['rank_search'] == 27 )) + { + $fSearchDiam = 0.001; + $sSQL = 'SELECT place_id,parent_place_id,30 as rank_search '; + if (CONST_Debug) { $sSQL .= ', housenumber, ST_distance('.$sPointSQL.', centroid) as distance, st_y(centroid) as lat, st_x(centroid) as lon'; } + $sSQL .= ' FROM location_property_tiger WHERE parent_place_id = '.$iPlaceID; + $sSQL .= ' AND ST_DWithin('.$sPointSQL.', centroid, '.$fSearchDiam.')'; + $sSQL .= ' ORDER BY ST_distance('.$sPointSQL.', centroid) ASC limit 1'; + + + // print all house numbers in the parent (street) + if (CONST_Debug) + { + $sSQL = preg_replace('/limit 1/', 'limit 100', $sSQL); + var_dump($sSQL); + + $aAllHouses = $this->oDB->getAll($sSQL); + foreach($aAllHouses as $i) + { + echo $i['housenumber'] . ' | ' . $i['distance'] * 1000 . ' | ' . $i['lat'] . ' | ' . $i['lon']. ' | '. "
\n"; + } + } + + $aPlaceTiger = $this->oDB->getRow($sSQL); + if (PEAR::IsError($aPlace)) + { + failInternalError("Could not determine closest Tiger place.", $sSQL, $aPlaceTiger); + } + if ($aPlaceTiger) + { + if (CONST_Debug) var_dump('found Tiger place', $aPlaceTiger); + $bPlaceIsTiger = true; + $aPlace = $aPlaceTiger; + $iPlaceID = $aPlaceTiger['place_id']; + $iParentPlaceID = $aPlaceTiger['parent_place_id']; // the street + } } // The point we found might be too small - use the address to find what it is a child of if ($iPlaceID && $iMaxRank < 28) { - if ($aPlace['rank_search'] > 28 && $iParentPlaceID) + if ($aPlace['rank_search'] > 28 && $iParentPlaceID && !$bPlaceIsTiger) { $iPlaceID = $iParentPlaceID; } @@ -132,6 +175,7 @@ $oPlaceLookup->setLanguagePreference($this->aLangPrefOrder); $oPlaceLookup->setIncludeAddressDetails($this->bAddressDetails); $oPlaceLookup->setPlaceId($iPlaceID); + $oPlaceLookup->setIsTiger($bPlaceIsTiger); return $oPlaceLookup->lookup(); } diff --git a/tests/features/api/reverse.feature b/tests/features/api/reverse.feature index b46ab2cd..2486b013 100644 --- a/tests/features/api/reverse.feature +++ b/tests/features/api/reverse.feature @@ -11,3 +11,15 @@ Feature: Reverse geocoding | ID | country | 0 | Deutschland + @Tiger + Scenario: TIGER house number + Given the request parameters + | addressdetails + | 1 + When looking up coordinates 40.6863624710666,-112.060005720023 + And exactly 1 result is returned + And result addresses contain + | ID | house_number | road | postcode | country_code + | 0 | 7094 | Kings Estate Drive | 84128 | us + And result 0 has not attributes osm_id,osm_type +