地理坐标周围的边界

时间:2014-02-22 17:10:26

标签: php google-maps google-maps-api-3 geocoding

我在Google地图上显示了一组地理坐标。如何确定最外面的点来构建包含所有点的多边形边界?

我知道如何创建一个多边形。但我需要确定所有外部点。

所有地理坐标都在mysql数据库中。我想使用PHP来分隔创建多边形所需的坐标。

我不知道如何做一个边界框。

1 个答案:

答案 0 :(得分:0)

你可以使用delaunay三角剖分并连接凸包上的所有顶点,即。找不到邻居的顶点。你可以下载我的php class convex hull @ phpclasses.org。您还可以查看我的网站http://www.phpdevpad.de/geofence作为delaunay三角测量和基于lat-lng对的凹壳的示例。你可以使用我的例子号。 5找到lat-lng对的凸壳:

require_once("convex-hull.php"); 
//example5 
$mapPadding  = 100; 
$mapWidth    = 500; 
$mapHeight   = 500; 
$mapLonLeft  =1000; 
$mapLatBottom=1000; 
$mapLonRight =   0; 
$mapLatTop   =   0; 
$set=array(); 
$geocoord = array ("8.6544487,50.1005233", 
                   "8.7839489,50.0907496", 
                   "8.1004734,50.2002273", 
                   "8.4117234,50.0951493", 
                   "8.3508367,49.4765982", 
                   "9.1828630,48.7827027", 
                   "9.1686483,48.7686426", 
                   "9.2118466,48.7829101", 
                   "8.9670738,48.9456327"); 

foreach ($geocoord as $key => $arr) 
{ 
    list($lon,$lat) = explode(",",$arr); 
    $mapLonLeft = min($mapLonLeft,$lon); 
    $mapLonRight = max($mapLonRight,$lon); 
    $mapLatBottom = min($mapLatBottom,$lat); 
    $mapLatTop = max($mapLatTop,$lat); 
    $set[]=array($lon,$lat); 
} 

$mapLonDelta = $mapLonRight-$mapLonLeft; 
$mapLatDelta = $mapLatTop-$mapLatBottom; 
$mapLatTopY=$mapLatTop*(M_PI/180); 
$worldMapWidth=(($mapWidth/$mapLonDelta)*360)/(2*M_PI); 
$LatBottomSin=min(max(sin($mapLatBottom*(M_PI/180)),-0.9999),0.9999); 
$mapOffsetY=$worldMapWidth/2 * log((1+$LatBottomSin)/(1-$LatBottomSin)); 
$LatTopSin=min(max(sin($mapLatTop*(M_PI/180)),-0.9999),0.9999); 
$mapOffsetTopY=$worldMapWidth/2 * log((1+$LatTopSin)/(1-$LatTopSin)); 
$mapHeightD=$mapOffsetTopY-$mapOffsetY; 
$mapRatioH=$mapHeight/$mapHeightD; 
$newWidth=$mapWidth*($mapHeightD/$mapHeight); 
$mapRatioW=$mapWidth/$newWidth; 

foreach ($set as $key => $arr) 
{ 
    list($lon,$lat) = $arr; 
    $tx = ($lon - $mapLonLeft) * ($newWidth/$mapLonDelta)*$mapRatioW; 
    $f = sin($lat*M_PI/180); 
    $ty = ($mapHeightD-(($worldMapWidth/2 * log((1+$f)/(1-$f)))-$mapOffsetY)); 
} 

$chull=new convexhull(); 
$chull->main($set,$mapWidth,$mapHeightD); 

然后凸包在数组中,您需要从超级角上移除顶点并检查顶点是否在lat-lng对中:

foreach ($chull->convexhull as $key => $arr)
      {
     foreach ($arr as $ikey => $iarr)
     {
        list($x1,$y1,$x2,$y2) = $iarr;
        if (abs($x1) != SUPER_TRIANGLE && abs($y1) != SUPER_TRIANGLE && abs($x2) != SUPER_TRIANGLE && abs($y2) != SUPER_TRIANGLE)
        {
           $ok=0;
           foreach ($chull->pointset as $iikey => $iiarr)
           {
          if ($iiarr==array($x1,$y1))
          {
             $ok=1;
          }
           }
           if ($ok)
           {
          solution[]=set[$iikey];  
           }
        }
     }
      }