为什么NaN从欧几里德方法返回

时间:2011-11-08 14:40:40

标签: php euclidean-distance

我对此功能有疑问:

function encludean_method($lat1,$lon1,$lat2,$lon2)
{
    $x1 = $lat1;
    $y1 = $lon1;
    $x2 = $lat2;
    $y2 = $lon2;
    $x = ( pow($x2,2) - pow($x1,2));
    $y = ( pow($y2,2) - pow($y1,2));
    $distance = (sqrt($x + $y));
    return $distance;
}

使用某些值(例如

)调用上述函数时出现的问题
(1.57454123333,103.6200516333,1.57483825,103.619484475)

它返回NaN。谁能告诉我为什么返回NaN以及如何解决它?

4 个答案:

答案 0 :(得分:7)

$ x和$ y可以假设为负值(实际上在您的示例中确实发生了这种情况),总计产生$ x + $ y的负表达式。我不确定你所做的是否正确,不正确的版本更像是:

$x1 = $lat1;
$y1 = $lon1;
$x2 = $lat2;
$y2 = $lon2;
$x = $x2 - $x1;
$y = $y2 - $y1;
$distance = (sqrt(pow($x, 2) + pow($y, 2)));
return $distance;

通过在计算平方根之前进行平方,您可以保证仅对正值进行操作。

答案 1 :(得分:1)

试试这个

function euclidean_method($lat1,$lon1,$lat2,$lon2) {
    $x = $lat1-$lat2;
    $y = $lon1-$lon2;
    return sqrt($x*$x+$y*$y)
}

答案 2 :(得分:0)

问题在于你的x值小于y的负值,因此它会成为你试图占据根的负值。

答案 3 :(得分:0)

由于算法错误。

function encludean_method($lat1,$lon1,$lat2,$lon2)
{
    $x1 = doubleval($lat1);
    $y1 = doubleval($lon1);
    $x2 = doubleval($lat2);
    $y2 = doubleval($lon2);

    $x = doubleval( pow($x2 - $x1, 2.0));
    $y = doubleval( pow($y2 - $y1, 2.0));

    $distance = doubleval(sqrt($x + $y));

    return ($distance);
}

注意:(x 1 2 - x 2 2 )&lt ; (x 1 - x 2 2

这就是您的代码提供x = float(0.00093541820670495)y = float(-0.11753762299304)的原因。因此,sqrt(x + y)太接近于零。

但是如果你试图获得该点在地球表面之间的距离,你应该使用另一个公式 - 你正在使用的公式仅适用于笛卡尔坐标系。

请注意,纬度和经度以度为单位。所以你得到的数字有点神秘=)

所以,谷歌搜索后,我注意到一些有趣的服务和公式。这是您可能尝试调用的代码:

function encludean_method($lat1, $lon1, $lat2, $lon2)
{
    $R = 6372.8; // km

    $lat1 = deg2rad($lat1);
    $lat2 = deg2rad($lat2);
    $lon1 = deg2rad($lon1);
    $lon2 = deg2rad($lon2);

    $a = (cos($lat2) * cos($lon2)) - (cos($lat1) * cos($lon1));
    $b = (cos($lat2) * sin($lon2)) - (cos($lat1) * sin($lon1));
    $c = sin($lat2) - sin($lat1);
    $ch = sqrt(pow($a, 2.0) + pow($b, 2.0) + pow($c, 2.0));
    $phi = 2.0 * asin($ch / 2.0);

    return $R * $phi;
}

现在,encludean_method(1.5745, 103.6200, 1.5748, 103.6195)返回0.0648375378577值,这可能是您要搜索的距离;以公里为单位。

以下是exampleone more(基于this answer)。以下是我用来验证该代码的service

希望这会对你有所帮助! =)