lib-php/PlaceLookup.php

   1 <?php
   2 /**
   3  * SPDX-License-Identifier: GPL-2.0-only
   4  *
   5  * This file is part of Nominatim. (https://nominatim.org)
   6  *
   7  * Copyright (C) 2022 by the Nominatim developer community.
   8  * For a full list of authors see the git log.
   9  */
  10
  11 namespace Nominatim;
  12
  13 require_once(CONST_LibDir.'/AddressDetails.php');
  14 require_once(CONST_LibDir.'/Result.php');
  15
  16 class PlaceLookup
  17 {
  18     protected $oDB;
  19
  20     protected $aLangPrefOrderSql = "''";
  21
  22     protected $bAddressDetails = false;
  23     protected $bExtraTags = false;
  24     protected $bNameDetails = false;
  25
  26     protected $bIncludePolygonAsText = false;
  27     protected $bIncludePolygonAsGeoJSON = false;
  28     protected $bIncludePolygonAsKML = false;
  29     protected $bIncludePolygonAsSVG = false;
  30     protected $fPolygonSimplificationThreshold = 0.0;
  31
  32     protected $sAnchorSql = null;
  33     protected $sAddressRankListSql = null;
  34     protected $sAllowedTypesSQLList = null;
  35     protected $bDeDupe = true;
  36
  37
  38     public function __construct(&$oDB)
  39     {
  40         $this->oDB =& $oDB;
  41     }
  42
  43     public function doDeDupe()
  44     {
  45         return $this->bDeDupe;
  46     }
  47
  48     public function setIncludeAddressDetails($b)
  49     {
  50         $this->bAddressDetails = $b;
  51     }
  52
  53     public function loadParamArray($oParams, $sGeomType = null)
  54     {
  55         $aLangs = $oParams->getPreferredLanguages();
  56         $this->aLangPrefOrderSql =
  57             'ARRAY['.join(',', $this->oDB->getDBQuotedList($aLangs)).']';
  58
  59         $this->bExtraTags = $oParams->getBool('extratags', false);
  60         $this->bNameDetails = $oParams->getBool('namedetails', false);
  61
  62         $this->bDeDupe = $oParams->getBool('dedupe', $this->bDeDupe);
  63
  64         if ($sGeomType === null || $sGeomType == 'geojson') {
  65             $this->bIncludePolygonAsGeoJSON = $oParams->getBool('polygon_geojson');
  66         }
  67
  68         if ($oParams->getString('format', '') !== 'geojson') {
  69             if ($sGeomType === null || $sGeomType == 'text') {
  70                 $this->bIncludePolygonAsText = $oParams->getBool('polygon_text');
  71             }
  72             if ($sGeomType === null || $sGeomType == 'kml') {
  73                 $this->bIncludePolygonAsKML = $oParams->getBool('polygon_kml');
  74             }
  75             if ($sGeomType === null || $sGeomType == 'svg') {
  76                 $this->bIncludePolygonAsSVG = $oParams->getBool('polygon_svg');
  77             }
  78         }
  79         $this->fPolygonSimplificationThreshold
  80             = $oParams->getFloat('polygon_threshold', 0.0);
  81
  82         $iWantedTypes =
  83             ($this->bIncludePolygonAsText ? 1 : 0) +
  84             ($this->bIncludePolygonAsGeoJSON ? 1 : 0) +
  85             ($this->bIncludePolygonAsKML ? 1 : 0) +
  86             ($this->bIncludePolygonAsSVG ? 1 : 0);
  87         if ($iWantedTypes > CONST_PolygonOutput_MaximumTypes) {
  88             if (CONST_PolygonOutput_MaximumTypes) {
  89                 userError('Select only '.CONST_PolygonOutput_MaximumTypes.' polgyon output option');
  90             } else {
  91                 userError('Polygon output is disabled');
  92             }
  93         }
  94     }
  95
  96     public function getMoreUrlParams()
  97     {
  98         $aParams = array();
  99
 100         if ($this->bAddressDetails) {
 101             $aParams['addressdetails'] = '1';
 102         }
 103         if ($this->bExtraTags) {
 104             $aParams['extratags'] = '1';
 105         }
 106         if ($this->bNameDetails) {
 107             $aParams['namedetails'] = '1';
 108         }
 109
 110         if ($this->bIncludePolygonAsText) {
 111             $aParams['polygon_text'] = '1';
 112         }
 113         if ($this->bIncludePolygonAsGeoJSON) {
 114             $aParams['polygon_geojson'] = '1';
 115         }
 116         if ($this->bIncludePolygonAsKML) {
 117             $aParams['polygon_kml'] = '1';
 118         }
 119         if ($this->bIncludePolygonAsSVG) {
 120             $aParams['polygon_svg'] = '1';
 121         }
 122
 123         if ($this->fPolygonSimplificationThreshold > 0.0) {
 124             $aParams['polygon_threshold'] = $this->fPolygonSimplificationThreshold;
 125         }
 126
 127         if (!$this->bDeDupe) {
 128             $aParams['dedupe'] = '0';
 129         }
 130
 131         return $aParams;
 132     }
 133
 134     public function setAnchorSql($sPoint)
 135     {
 136         $this->sAnchorSql = $sPoint;
 137     }
 138
 139     public function setAddressRankList($aList)
 140     {
 141         $this->sAddressRankListSql = '('.join(',', $aList).')';
 142     }
 143
 144     public function setAllowedTypesSQLList($sSql)
 145     {
 146         $this->sAllowedTypesSQLList = $sSql;
 147     }
 148
 149     public function setLanguagePreference($aLangPrefOrder)
 150     {
 151         $this->aLangPrefOrderSql = $this->oDB->getArraySQL(
 152             $this->oDB->getDBQuotedList($aLangPrefOrder)
 153         );
 154     }
 155
 156     private function addressImportanceSql($sGeometry, $sPlaceId)
 157     {
 158         if ($this->sAnchorSql) {
 159             $sSQL = 'ST_Distance('.$this->sAnchorSql.','.$sGeometry.')';
 160         } else {
 161             $sSQL = '(SELECT max(ai_p.importance * (ai_p.rank_address + 2))';
 162             $sSQL .= '   FROM place_addressline ai_s, placex ai_p';
 163             $sSQL .= '   WHERE ai_s.place_id = '.$sPlaceId;
 164             $sSQL .= '     AND ai_p.place_id = ai_s.address_place_id ';
 165             $sSQL .= '     AND ai_s.isaddress ';
 166             $sSQL .= '     AND ai_p.importance is not null)';
 167         }
 168
 169         return $sSQL.' AS addressimportance,';
 170     }
 171
 172     private function langAddressSql($sHousenumber)
 173     {
 174         if ($this->bAddressDetails) {
 175             return ''; // langaddress will be computed from address details
 176         }
 177
 178         return 'get_address_by_language(place_id,'.$sHousenumber.','.$this->aLangPrefOrderSql.') AS langaddress,';
 179     }
 180
 181     public function lookupOSMID($sType, $iID)
 182     {
 183         $sSQL = 'select place_id from placex where osm_type = :type and osm_id = :id';
 184         $iPlaceID = $this->oDB->getOne($sSQL, array(':type' => $sType, ':id' => $iID));
 185
 186         if (!$iPlaceID) {
 187             return null;
 188         }
 189
 190         $aResults = $this->lookup(array($iPlaceID => new Result($iPlaceID)), 0, 30, true);
 191
 192         return empty($aResults) ? null : reset($aResults);
 193     }
 194
 195     public function lookup($aResults, $iMinRank = 0, $iMaxRank = 30, $bAllowLinked = false)
 196     {
 197         Debug::newFunction('Place lookup');
 198
 199         if (empty($aResults)) {
 200             return array();
 201         }
 202         $aSubSelects = array();
 203
 204         $sPlaceIDs = Result::joinIdsByTable($aResults, Result::TABLE_PLACEX);
 205         if ($sPlaceIDs) {
 206             Debug::printVar('Ids from placex', $sPlaceIDs);
 207             $sSQL  = 'SELECT ';
 208             $sSQL .= '    osm_type,';
 209             $sSQL .= '    osm_id,';
 210             $sSQL .= '    class,';
 211             $sSQL .= '    type,';
 212             $sSQL .= '    admin_level,';
 213             $sSQL .= '    rank_search,';
 214             $sSQL .= '    rank_address,';
 215             $sSQL .= '    min(place_id) AS place_id,';
 216             $sSQL .= '    min(parent_place_id) AS parent_place_id,';
 217             $sSQL .= '    -1 as housenumber,';
 218             $sSQL .= '    country_code,';
 219             $sSQL .= $this->langAddressSql('-1');
 220             $sSQL .= '    get_name_by_language(name,'.$this->aLangPrefOrderSql.') AS placename,';
 221             $sSQL .= "    get_name_by_language(name, ARRAY['ref']) AS ref,";
 222             if ($this->bExtraTags) {
 223                 $sSQL .= 'hstore_to_json(extratags)::text AS extra,';
 224             }
 225             if ($this->bNameDetails) {
 226                 $sSQL .= 'hstore_to_json(name)::text AS names,';
 227             }
 228             $sSQL .= '    avg(ST_X(centroid)) AS lon, ';
 229             $sSQL .= '    avg(ST_Y(centroid)) AS lat, ';
 230             $sSQL .= '    COALESCE(importance,0.75-(rank_search::float/40)) AS importance, ';
 231             $sSQL .= $this->addressImportanceSql(
 232                 'ST_Collect(centroid)',
 233                 'min(CASE WHEN placex.rank_search < 28 THEN placex.place_id ELSE placex.parent_place_id END)'
 234             );
 235             $sSQL .= "    COALESCE(extratags->'place', extratags->'linked_place') AS extra_place ";
 236             $sSQL .= ' FROM placex';
 237             $sSQL .= " WHERE place_id in ($sPlaceIDs) ";
 238             $sSQL .= '   AND (';
 239             $sSQL .= "        placex.rank_address between $iMinRank and $iMaxRank ";
 240             if (14 >= $iMinRank && 14 <= $iMaxRank) {
 241                 $sSQL .= "    OR (extratags->'place') = 'city'";
 242             }
 243             if ($this->sAddressRankListSql) {
 244                 $sSQL .= '    OR placex.rank_address in '.$this->sAddressRankListSql;
 245             }
 246             $sSQL .= '       ) ';
 247             if ($this->sAllowedTypesSQLList) {
 248                 $sSQL .= 'AND placex.class in '.$this->sAllowedTypesSQLList;
 249             }
 250             if (!$bAllowLinked) {
 251                 $sSQL .= '    AND linked_place_id is null ';
 252             }
 253             $sSQL .= ' GROUP BY ';
 254             $sSQL .= '     osm_type, ';
 255             $sSQL .= '     osm_id, ';
 256             $sSQL .= '     class, ';
 257             $sSQL .= '     type, ';
 258             $sSQL .= '     admin_level, ';
 259             $sSQL .= '     rank_search, ';
 260             $sSQL .= '     rank_address, ';
 261             $sSQL .= '     housenumber,';
 262             $sSQL .= '     country_code, ';
 263             $sSQL .= '     importance, ';
 264             if (!$this->bDeDupe) {
 265                 $sSQL .= 'place_id,';
 266             }
 267             if (!$this->bAddressDetails) {
 268                 $sSQL .= 'langaddress, ';
 269             }
 270             $sSQL .= '     placename, ';
 271             $sSQL .= '     ref, ';
 272             if ($this->bExtraTags) {
 273                 $sSQL .= 'extratags, ';
 274             }
 275             if ($this->bNameDetails) {
 276                 $sSQL .= 'name, ';
 277             }
 278             $sSQL .= '     extra_place ';
 279
 280             $aSubSelects[] = $sSQL;
 281         }
 282
 283         // postcode table
 284         $sPlaceIDs = Result::joinIdsByTable($aResults, Result::TABLE_POSTCODE);
 285         if ($sPlaceIDs) {
 286             Debug::printVar('Ids from location_postcode', $sPlaceIDs);
 287             $sSQL = 'SELECT';
 288             $sSQL .= "  'P' as osm_type,";
 289             $sSQL .= '  (SELECT osm_id from placex p WHERE p.place_id = lp.parent_place_id) as osm_id,';
 290             $sSQL .= "  'place' as class, 'postcode' as type,";
 291             $sSQL .= '  null::smallint as admin_level, rank_search, rank_address,';
 292             $sSQL .= '  place_id, parent_place_id,';
 293             $sSQL .= '  -1 as housenumber,';
 294             $sSQL .= '  country_code,';
 295             $sSQL .= $this->langAddressSql('-1');
 296             $sSQL .= '  postcode as placename,';
 297             $sSQL .= '  postcode as ref,';
 298             if ($this->bExtraTags) {
 299                 $sSQL .= 'null::text AS extra,';
 300             }
 301             if ($this->bNameDetails) {
 302                 $sSQL .= 'null::text AS names,';
 303             }
 304             $sSQL .= '  ST_x(geometry) AS lon, ST_y(geometry) AS lat,';
 305             $sSQL .= '  (0.75-(rank_search::float/40)) AS importance, ';
 306             $sSQL .= $this->addressImportanceSql('geometry', 'lp.parent_place_id');
 307             $sSQL .= '  null::text AS extra_place ';
 308             $sSQL .= 'FROM location_postcode lp';
 309             $sSQL .= " WHERE place_id in ($sPlaceIDs) ";
 310             $sSQL .= "   AND lp.rank_address between $iMinRank and $iMaxRank";
 311
 312             $aSubSelects[] = $sSQL;
 313         }
 314
 315         // All other tables are rank 30 only.
 316         if ($iMaxRank == 30) {
 317             // TIGER table
 318             if (CONST_Use_US_Tiger_Data) {
 319                 $sPlaceIDs = Result::joinIdsByTable($aResults, Result::TABLE_TIGER);
 320                 if ($sPlaceIDs) {
 321                     Debug::printVar('Ids from Tiger table', $sPlaceIDs);
 322                     $sHousenumbers = Result::sqlHouseNumberTable($aResults, Result::TABLE_TIGER);
 323                     // Tiger search only if a housenumber was searched and if it was found
 324                     // (realized through a join)
 325                     $sSQL = ' SELECT ';
 326                     $sSQL .= "     'T' AS osm_type, ";
 327                     $sSQL .= '     (SELECT osm_id from placex p WHERE p.place_id=blub.parent_place_id) as osm_id, ';
 328                     $sSQL .= "     'place' AS class, ";
 329                     $sSQL .= "     'house' AS type, ";
 330                     $sSQL .= '     null::smallint AS admin_level, ';
 331                     $sSQL .= '     30 AS rank_search, ';
 332                     $sSQL .= '     30 AS rank_address, ';
 333                     $sSQL .= '     place_id, ';
 334                     $sSQL .= '     parent_place_id, ';
 335                     $sSQL .= '     housenumber_for_place as housenumber,';
 336                     $sSQL .= "     'us' AS country_code, ";
 337                     $sSQL .= $this->langAddressSql('housenumber_for_place');
 338                     $sSQL .= '     null::text AS placename, ';
 339                     $sSQL .= '     null::text AS ref, ';
 340                     if ($this->bExtraTags) {
 341                         $sSQL .= 'null::text AS extra,';
 342                     }
 343                     if ($this->bNameDetails) {
 344                         $sSQL .= 'null::text AS names,';
 345                     }
 346                     $sSQL .= '     st_x(centroid) AS lon, ';
 347                     $sSQL .= '     st_y(centroid) AS lat,';
 348                     $sSQL .= '     -1.15 AS importance, ';
 349                     $sSQL .= $this->addressImportanceSql('centroid', 'blub.parent_place_id');
 350                     $sSQL .= '     null::text AS extra_place ';
 351                     $sSQL .= ' FROM (';
 352                     $sSQL .= '     SELECT place_id, ';    // interpolate the Tiger housenumbers here
 353                     $sSQL .= '         CASE WHEN startnumber != endnumber';
 354                     $sSQL .= '              THEN ST_LineInterpolatePoint(linegeo, (housenumber_for_place-startnumber::float)/(endnumber-startnumber)::float)';
 355                     $sSQL .= '              ELSE ST_LineInterpolatePoint(linegeo, 0.5) END AS centroid, ';
 356                     $sSQL .= '         parent_place_id, ';
 357                     $sSQL .= '         housenumber_for_place';
 358                     $sSQL .= '     FROM (';
 359                     $sSQL .= '            location_property_tiger ';
 360                     $sSQL .= '            JOIN (values '.$sHousenumbers.') AS housenumbers(place_id, housenumber_for_place) USING(place_id)) ';
 361                     $sSQL .= '     WHERE ';
 362                     $sSQL .= '         housenumber_for_place >= startnumber';
 363                     $sSQL .= '         AND housenumber_for_place <= endnumber';
 364                     $sSQL .= ' ) AS blub'; //postgres wants an alias here
 365
 366                     $aSubSelects[] = $sSQL;
 367                 }
 368             }
 369
 370             // osmline - interpolated housenumbers
 371             $sPlaceIDs = Result::joinIdsByTable($aResults, Result::TABLE_OSMLINE);
 372             if ($sPlaceIDs) {
 373                 Debug::printVar('Ids from interpolation', $sPlaceIDs);
 374                 $sHousenumbers = Result::sqlHouseNumberTable($aResults, Result::TABLE_OSMLINE);
 375                 // interpolation line search only if a housenumber was searched
 376                 // (realized through a join)
 377                 $sSQL = 'SELECT ';
 378                 $sSQL .= "  'W' AS osm_type, ";
 379                 $sSQL .= '  osm_id, ';
 380                 $sSQL .= "  'place' AS class, ";
 381                 $sSQL .= "  'house' AS type, ";
 382                 $sSQL .= '  null::smallint AS admin_level, ';
 383                 $sSQL .= '  30 AS rank_search, ';
 384                 $sSQL .= '  30 AS rank_address, ';
 385                 $sSQL .= '  place_id, ';
 386                 $sSQL .= '  parent_place_id, ';
 387                 $sSQL .= '  housenumber_for_place as housenumber,';
 388                 $sSQL .= '  country_code, ';
 389                 $sSQL .= $this->langAddressSql('housenumber_for_place');
 390                 $sSQL .= '  null::text AS placename, ';
 391                 $sSQL .= '  null::text AS ref, ';
 392                 if ($this->bExtraTags) {
 393                     $sSQL .= 'null::text AS extra, ';
 394                 }
 395                 if ($this->bNameDetails) {
 396                     $sSQL .= 'null::text AS names, ';
 397                 }
 398                 $sSQL .= '  st_x(centroid) AS lon, ';
 399                 $sSQL .= '  st_y(centroid) AS lat, ';
 400                 // slightly smaller than the importance for normal houses
 401                 $sSQL .= '  -0.1 AS importance, ';
 402                 $sSQL .= $this->addressImportanceSql('centroid', 'blub.parent_place_id');
 403                 $sSQL .= '  null::text AS extra_place ';
 404                 $sSQL .= '  FROM (';
 405                 $sSQL .= '     SELECT ';
 406                 $sSQL .= '         osm_id, ';
 407                 $sSQL .= '         place_id, ';
 408                 $sSQL .= '         country_code, ';
 409                 $sSQL .= '         CASE ';             // interpolate the housenumbers here
 410                 $sSQL .= '           WHEN startnumber != endnumber ';
 411                 $sSQL .= '           THEN ST_LineInterpolatePoint(linegeo, (housenumber_for_place-startnumber::float)/(endnumber-startnumber)::float) ';
 412                 $sSQL .= '           ELSE linegeo ';
 413                 $sSQL .= '         END as centroid, ';
 414                 $sSQL .= '         parent_place_id, ';
 415                 $sSQL .= '         housenumber_for_place ';
 416                 $sSQL .= '     FROM (';
 417                 $sSQL .= '            location_property_osmline ';
 418                 $sSQL .= '            JOIN (values '.$sHousenumbers.') AS housenumbers(place_id, housenumber_for_place) USING(place_id)';
 419                 $sSQL .= '          ) ';
 420                 $sSQL .= '     WHERE housenumber_for_place >= 0 ';
 421                 $sSQL .= '  ) as blub'; //postgres wants an alias here
 422
 423                 $aSubSelects[] = $sSQL;
 424             }
 425         }
 426
 427         if (empty($aSubSelects)) {
 428             return array();
 429         }
 430
 431         $sSQL = join(' UNION ', $aSubSelects);
 432         Debug::printSQL($sSQL);
 433         $aPlaces = $this->oDB->getAll($sSQL, null, 'Could not lookup place');
 434
 435         foreach ($aPlaces as &$aPlace) {
 436             $aPlace['importance'] = (float) $aPlace['importance'];
 437             if ($this->bAddressDetails) {
 438                 // to get addressdetails for tiger data, the housenumber is needed
 439                 $aPlace['address'] = new AddressDetails(
 440                     $this->oDB,
 441                     $aPlace['place_id'],
 442                     $aPlace['housenumber'],
 443                     $this->aLangPrefOrderSql
 444                 );
 445                 $aPlace['langaddress'] = $aPlace['address']->getLocaleAddress();
 446             }
 447
 448             if ($this->bExtraTags) {
 449                 if ($aPlace['extra']) {
 450                     $aPlace['sExtraTags'] = json_decode($aPlace['extra'], true);
 451                 } else {
 452                     $aPlace['sExtraTags'] = (object) array();
 453                 }
 454             }
 455
 456             if ($this->bNameDetails) {
 457                 $aPlace['sNameDetails'] = $this->extractNames($aPlace['names']);
 458             }
 459
 460             $aPlace['addresstype'] = ClassTypes\getLabelTag(
 461                 $aPlace,
 462                 $aPlace['country_code']
 463             );
 464
 465             $aResults[$aPlace['place_id']] = $aPlace;
 466         }
 467
 468         $aResults = array_filter(
 469             $aResults,
 470             function ($v) {
 471                 return !($v instanceof Result);
 472             }
 473         );
 474
 475         Debug::printVar('Places', $aResults);
 476
 477         return $aResults;
 478     }
 479
 480
 481     private function extractNames($sNames)
 482     {
 483         if (!$sNames) {
 484             return (object) array();
 485         }
 486
 487         $aFullNames = json_decode($sNames, true);
 488         $aNames = array();
 489
 490         foreach ($aFullNames as $sKey => $sValue) {
 491             if (strpos($sKey, '_place_') === 0) {
 492                 $sSubKey = substr($sKey, 7);
 493                 if (array_key_exists($sSubKey, $aFullNames)) {
 494                     $aNames[$sKey] = $sValue;
 495                 } else {
 496                     $aNames[$sSubKey] = $sValue;
 497                 }
 498             } else {
 499                 $aNames[$sKey] = $sValue;
 500             }
 501         }
 502
 503         return $aNames;
 504     }
 505
 506
 507     /* returns an array which will contain the keys
 508      *   aBoundingBox
 509      * and may also contain one or more of the keys
 510      *   asgeojson
 511      *   askml
 512      *   assvg
 513      *   astext
 514      *   lat
 515      *   lon
 516      */
 517     public function getOutlines($iPlaceID, $fLon = null, $fLat = null, $fRadius = null, $fLonReverse = null, $fLatReverse = null)
 518     {
 519
 520         $aOutlineResult = array();
 521         if (!$iPlaceID) {
 522             return $aOutlineResult;
 523         }
 524
 525         // Get the bounding box and outline polygon
 526         $sSQL = 'select place_id,0 as numfeatures,st_area(geometry) as area,';
 527         if ($fLonReverse != null && $fLatReverse != null) {
 528             $sSQL .= ' ST_Y(closest_point) as centrelat,';
 529             $sSQL .= ' ST_X(closest_point) as centrelon,';
 530         } else {
 531             $sSQL .= ' ST_Y(centroid) as centrelat, ST_X(centroid) as centrelon,';
 532         }
 533         $sSQL .= ' ST_YMin(geometry) as minlat,ST_YMax(geometry) as maxlat,';
 534         $sSQL .= ' ST_XMin(geometry) as minlon,ST_XMax(geometry) as maxlon';
 535         if ($this->bIncludePolygonAsGeoJSON) {
 536             $sSQL .= ',ST_AsGeoJSON(geometry) as asgeojson';
 537         }
 538         if ($this->bIncludePolygonAsKML) {
 539             $sSQL .= ',ST_AsKML(geometry) as askml';
 540         }
 541         if ($this->bIncludePolygonAsSVG) {
 542             $sSQL .= ',ST_AsSVG(geometry) as assvg';
 543         }
 544         if ($this->bIncludePolygonAsText) {
 545             $sSQL .= ',ST_AsText(geometry) as astext';
 546         }
 547         if ($fLonReverse != null && $fLatReverse != null) {
 548             $sFrom = ' from (SELECT * , CASE WHEN (class = \'highway\') AND (ST_GeometryType(geometry) = \'ST_LineString\') THEN ';
 549             $sFrom .=' ST_ClosestPoint(geometry, ST_SetSRID(ST_Point('.$fLatReverse.','.$fLonReverse.'),4326))';
 550             $sFrom .=' ELSE centroid END AS closest_point';
 551             $sFrom .= ' from placex where place_id = '.$iPlaceID.') as plx';
 552         } else {
 553             $sFrom = ' from placex where place_id = '.$iPlaceID;
 554         }
 555         if ($this->fPolygonSimplificationThreshold > 0) {
 556             $sSQL .= ' from (select place_id,centroid,ST_SimplifyPreserveTopology(geometry,'.$this->fPolygonSimplificationThreshold.') as geometry'.$sFrom.') as plx';
 557         } else {
 558             $sSQL .= $sFrom;
 559         }
 560
 561         $aPointPolygon = $this->oDB->getRow($sSQL, null, 'Could not get outline');
 562
 563         if ($aPointPolygon && $aPointPolygon['place_id']) {
 564             if ($aPointPolygon['centrelon'] !== null && $aPointPolygon['centrelat'] !== null) {
 565                 $aOutlineResult['lat'] = $aPointPolygon['centrelat'];
 566                 $aOutlineResult['lon'] = $aPointPolygon['centrelon'];
 567             }
 568
 569             if ($this->bIncludePolygonAsGeoJSON) {
 570                 $aOutlineResult['asgeojson'] = $aPointPolygon['asgeojson'];
 571             }
 572             if ($this->bIncludePolygonAsKML) {
 573                 $aOutlineResult['askml'] = $aPointPolygon['askml'];
 574             }
 575             if ($this->bIncludePolygonAsSVG) {
 576                 $aOutlineResult['assvg'] = $aPointPolygon['assvg'];
 577             }
 578             if ($this->bIncludePolygonAsText) {
 579                 $aOutlineResult['astext'] = $aPointPolygon['astext'];
 580             }
 581
 582             if (abs($aPointPolygon['minlat'] - $aPointPolygon['maxlat']) < 0.0000001) {
 583                 $aPointPolygon['minlat'] = $aPointPolygon['minlat'] - $fRadius;
 584                 $aPointPolygon['maxlat'] = $aPointPolygon['maxlat'] + $fRadius;
 585             }
 586
 587             if (abs($aPointPolygon['minlon'] - $aPointPolygon['maxlon']) < 0.0000001) {
 588                 $aPointPolygon['minlon'] = $aPointPolygon['minlon'] - $fRadius;
 589                 $aPointPolygon['maxlon'] = $aPointPolygon['maxlon'] + $fRadius;
 590             }
 591
 592             $aOutlineResult['aBoundingBox'] = array(
 593                                                (string)$aPointPolygon['minlat'],
 594                                                (string)$aPointPolygon['maxlat'],
 595                                                (string)$aPointPolygon['minlon'],
 596                                                (string)$aPointPolygon['maxlon']
 597                                               );
 598         }
 599
 600         // as a fallback we generate a bounding box without knowing the size of the geometry
 601         if ((!isset($aOutlineResult['aBoundingBox'])) && isset($fLon)) {
 602             $aBounds = array(
 603                         'minlat' => $fLat - $fRadius,
 604                         'maxlat' => $fLat + $fRadius,
 605                         'minlon' => $fLon - $fRadius,
 606                         'maxlon' => $fLon + $fRadius
 607                        );
 608
 609             $aOutlineResult['aBoundingBox'] = array(
 610                                                (string)$aBounds['minlat'],
 611                                                (string)$aBounds['maxlat'],
 612                                                (string)$aBounds['minlon'],
 613                                                (string)$aBounds['maxlon']
 614                                               );
 615         }
 616         return $aOutlineResult;
 617     }
 618 }