Question

尝试从.gpx文件计算每个gps点之间的距离。我尝试了两种不同的配方。这个应该更准确： dist = 6378.388 * acos（sin（lat1）* sin（lat2）+ cos（lat1）* cos（lat2）* cos（lon2 - lon1））

这个也适用于短距离：距离= sqrt（dx * dx + dy * dy）

每个点之间的距离是最大值。 100米我用我的程序计算了距离，得到了26.8公里。然后我将.gpx文件导入了goole earth，距离我有21公里。我不知道为什么这些价值观差异如此之大。也许你可以帮助我。

这里是代码（只是相关片段）来计算与gpx文件的距离：

...
$xml = file_get_contents($file);
$doc = new DOMDocument();
$doc->loadXML($xml);
$track = array();

/* get the elevation for each gps point */
...

/* get the time for each gps point */
...

/* get the position for each gps point*/                                
$trkpts = $doc->getElementsByTagName('trkpt');
$i=0;
foreach ($trkpts as $trkpt) {
   $track[$i]['position']   = array();
   $track[$i]['position']['lon']    = $trkpt->getAttribute('lat');
   $track[$i]['position']['lat']    = $trkpt->getAttribute('lon');
   $i++;
}

/* calculate the distance for each gps point */
ini_set('precision', '50');
for($i2=0;$i2<($i);$i2++){
    if($i2==0){
      $track[$i2]['distance'] = 0;
      /* get first point */
      $lastlat  = deg2rad($track[$i2]['position']['lat']);
      $lastlon  = deg2rad($track[$i2]['position']['lon']);
      $lastele  = $track[$i2]['ele'];
    }
    else{

       $aktlat = deg2rad($track[$i2]['position']['lat']);
       $aktlon  = deg2rad($track[$i2]['position']['lon']);
       $aktele  = $track[$i2]['ele'];

       /* if coordinates differ calculate distance */                   
       if(($aktlat-$lastlat)!=0 || ($aktlon-$lastlon)!=0){
            $track[$i2]['distance'] = acos(sin($lastlat)*sin($aktlat)+cos($lastlat)*cos($aktlat)*cos($aktlon-$lastlon));

            $track[$i2]['distance'] = $track[$i2]['distance'] * 6378137;

            /* calculate the elevation difference */
            $track[$i2]['distance'] = sqrt(pow($track[$i2]['distance'],2)+pow(($aktele-$lastele),2));

       }
       else{
           $track[$i2]['distance'] = 0;
       }
       /* save act point as last point */
       $lastlat                 = deg2rad($aktlat);
       $lastlon                 = deg2rad($aktlon);
       $lastele             = $aktele;
    }
}

以下是每个想要查看的人的gpx文件：http://mein-sporttagebuch.de/userfiles/tempfiles/runtastic_20130825_1054_Radfahren.gpx

感谢您的任何建议。

修改这里是获取高程的代码（从高程我计算出更多的测量值，如最高点和最深点）：

$eles = $doc->getElementsByTagName('ele');
$i=0;
$minheight  = 0;
$maxheight  = 0;
$asc        = 0;
$desc       = 0;
$lastheight = 0;
foreach ($eles as $ele) {
   if($i>0){
      if($lastheight<$ele->nodeValue){
         $asc +=  $ele->nodeValue-$lastheight;
      }
      else{
         $desc +=  $lastheight-$ele->nodeValue;
      }
   }
   $track[$i] = array();
   $track[$i]['ele'] = $ele->nodeValue;
   if($minheight>$ele->nodeValue || $i==0){
      $minheight = $ele->nodeValue;
   }
   if($maxheight<$ele->nodeValue){
      $maxheight = $ele->nodeValue;
   }

   $lastheight = $ele->nodeValue;
   $i++;
}

Answer 1

循环结束时代码中存在错误。

   /* save act point as last point */
   $lastlat                 = deg2rad($aktlat);
   $lastlon                 = deg2rad($aktlon);

但是aktlat / lon被定义为：

   $aktlat = deg2rad($track[$i2]['position']['lat']);
   $aktlon  = deg2rad($track[$i2]['position']['lon']);

...你已经完成了两次deg2rad转换。它应该是：

   /* save act point as last point */
   $lastlat                 = $aktlat;
   $lastlon                 = $aktlon;

不确定这是不是导致错误的原因，但它无能为力！您还可以使用顶部公式在代码中使用不同的地球半径，但这几乎没有什么区别。

Answer 2

正如我所见，你在计算中取得了提升。可能，谷歌地图不计算其公式中的高程。您是否尝试计算“平坦”方块的距离？

计算gps点之间的距离

2 个答案: