การวัดระยะทางระหว่างสองพิกัดใน PHP


146

สวัสดีฉันต้องการคำนวณระยะทางระหว่างจุดสองจุดที่มีละติจูดและลองจิจูด

ฉันต้องการหลีกเลี่ยงการเรียก API ภายนอกใด ๆ

ฉันพยายามใช้สูตร Haversine ใน PHP:

นี่คือรหัส:

class CoordDistance
 {
    public $lat_a = 0;
    public $lon_a = 0;
    public $lat_b = 0;
    public $lon_b = 0;

    public $measure_unit = 'kilometers';

    public $measure_state = false;

    public $measure = 0;

    public $error = '';



    public function DistAB()

      {
          $delta_lat = $this->lat_b - $this->lat_a ;
          $delta_lon = $this->lon_b - $this->lon_a ;

          $earth_radius = 6372.795477598;

          $alpha    = $delta_lat/2;
          $beta     = $delta_lon/2;
          $a        = sin(deg2rad($alpha)) * sin(deg2rad($alpha)) + cos(deg2rad($this->lat_a)) * cos(deg2rad($this->lat_b)) * sin(deg2rad($beta)) * sin(deg2rad($beta)) ;
          $c        = asin(min(1, sqrt($a)));
          $distance = 2*$earth_radius * $c;
          $distance = round($distance, 4);

          $this->measure = $distance;

      }
    }

ทดสอบด้วยคะแนนที่กำหนดซึ่งมีระยะทางสาธารณะฉันไม่ได้ผลลัพธ์ที่น่าเชื่อถือ

ฉันไม่เข้าใจว่ามีข้อผิดพลาดในสูตรดั้งเดิมหรือในการนำไปใช้ของฉัน


4
ฉันพบรหัสการทำงานที่นี่ในหลายภาษารวมถึง php geodatasource.com/developers/php
krishna

คำตอบ:


273

ไม่นานที่ผ่านมาฉันเขียนตัวอย่างของสูตรแฮเวอรีนและเผยแพร่บนเว็บไซต์ของฉัน:

/**
 * Calculates the great-circle distance between two points, with
 * the Haversine formula.
 * @param float $latitudeFrom Latitude of start point in [deg decimal]
 * @param float $longitudeFrom Longitude of start point in [deg decimal]
 * @param float $latitudeTo Latitude of target point in [deg decimal]
 * @param float $longitudeTo Longitude of target point in [deg decimal]
 * @param float $earthRadius Mean earth radius in [m]
 * @return float Distance between points in [m] (same as earthRadius)
 */
function haversineGreatCircleDistance(
  $latitudeFrom, $longitudeFrom, $latitudeTo, $longitudeTo, $earthRadius = 6371000)
{
  // convert from degrees to radians
  $latFrom = deg2rad($latitudeFrom);
  $lonFrom = deg2rad($longitudeFrom);
  $latTo = deg2rad($latitudeTo);
  $lonTo = deg2rad($longitudeTo);

  $latDelta = $latTo - $latFrom;
  $lonDelta = $lonTo - $lonFrom;

  $angle = 2 * asin(sqrt(pow(sin($latDelta / 2), 2) +
    cos($latFrom) * cos($latTo) * pow(sin($lonDelta / 2), 2)));
  return $angle * $earthRadius;
}

$earthRadius➽หมายเหตุว่าคุณได้รับกลับระยะทางในหน่วยเดียวกับที่คุณผ่านในพารามิเตอร์ ค่าเริ่มต้นคือ 6371000 เมตรดังนั้นผลลัพธ์จะเป็น [m] ด้วย ในการรับผลลัพธ์เป็นไมล์คุณสามารถผ่าน 3959 ไมล์ได้$earthRadiusและผลลัพธ์จะเป็น [mi] ในความคิดของฉันมันเป็นนิสัยที่ดีที่จะยึดติดกับหน่วย SI หากไม่มีเหตุผลใดที่จะทำอย่างอื่น

แก้ไข:

ในฐานะที่เป็น TreyA ชี้ได้อย่างถูกต้องออกสูตร Haversine มีจุดอ่อนกับจุดที่ตรงกันข้ามกับเท้าเนื่องจากข้อผิดพลาดในการปัดเศษ (แม้ว่ามันจะเป็นมีเสถียรภาพสำหรับระยะทางขนาดเล็ก) คุณสามารถใช้สูตร Vincentyแทนได้

/**
 * Calculates the great-circle distance between two points, with
 * the Vincenty formula.
 * @param float $latitudeFrom Latitude of start point in [deg decimal]
 * @param float $longitudeFrom Longitude of start point in [deg decimal]
 * @param float $latitudeTo Latitude of target point in [deg decimal]
 * @param float $longitudeTo Longitude of target point in [deg decimal]
 * @param float $earthRadius Mean earth radius in [m]
 * @return float Distance between points in [m] (same as earthRadius)
 */
public static function vincentyGreatCircleDistance(
  $latitudeFrom, $longitudeFrom, $latitudeTo, $longitudeTo, $earthRadius = 6371000)
{
  // convert from degrees to radians
  $latFrom = deg2rad($latitudeFrom);
  $lonFrom = deg2rad($longitudeFrom);
  $latTo = deg2rad($latitudeTo);
  $lonTo = deg2rad($longitudeTo);

  $lonDelta = $lonTo - $lonFrom;
  $a = pow(cos($latTo) * sin($lonDelta), 2) +
    pow(cos($latFrom) * sin($latTo) - sin($latFrom) * cos($latTo) * cos($lonDelta), 2);
  $b = sin($latFrom) * sin($latTo) + cos($latFrom) * cos($latTo) * cos($lonDelta);

  $angle = atan2(sqrt($a), $b);
  return $angle * $earthRadius;
}

1
@TreyA - มีหลายเวอร์ชั่นที่เป็นไปได้รุ่นนี้ใช้สูตรในWikipediaและผ่านการทดสอบอย่างดี $ angle หมายถึงมุมที่อยู่ตรงกลางของโลกเป็นเรเดียนดังนั้นเราจึงสามารถคูณมันด้วยรัศมีของโลก ฉันยังสามารถแสดงตัวอย่างของสูตร Vincenty ที่ซับซ้อนมากขึ้นถ้ามีคนสนใจ
martinstoeckli

@TreyA - ใช่ฉันรู้ฉันไม่แน่ใจว่าคุณต้องการพูดอะไร คุณได้ทดสอบฟังก์ชั่นแล้วมันคำนวณผลที่ผิดหรือเปล่า? คุณเคยดูสูตรใน Wikipedia ไหม? คุณควรทำแบบทดสอบของคุณเองและให้ตัวอย่างสิ่งที่คุณคิดว่าผิด
martinstoeckli

ขออภัยฉันต้องอธิบายบางสิ่งบางอย่างในขณะนี้ 1) คำถามเกี่ยวกับสูตร Haversine คุณควรบอกเราว่าคุณแนะนำให้ใช้สูตรอื่นหรือไม่ 2) สูตร Haversine มีจุดอ่อนรอบเสา แต่เป็นที่ถูกต้องสำหรับระยะทางขนาดเล็ก (มันเป็นปัญหาของสูตร arccosine ก) 3) คุณระบุว่ามีขั้นตอนที่ขาดหายไปจากการคำนวณ $ angle นั่นเป็นความผิดพลาดเพียงอย่างเดียวมันไม่สามารถปรับปรุงผลลัพธ์ได้โปรดช่วยทดสอบด้วย! 4) ฉันยอมรับว่ามันจะดีกว่าถ้าใช้สูตรมั่นคง Vincenty ฉันเสนอให้เป็นตัวอย่างแล้ว บางทีคุณอาจจะเขียนคำตอบด้วยก็ได้
martinstoeckli

@martinstoekli - คุณถูกต้องคุณไม่มีขั้นตอนที่ขาดหายไปในสูตร Haversine ของคุณ ฉันลบความคิดเห็นของฉันเพื่อไม่ให้ผู้อ่านในอนาคตสับสน
TreyA

1
@PratikCJoshi - ในที่สุดก็พบเวลาที่จะเพิ่มบันทึกเกี่ยวกับการใช้หน่วยที่แตกต่างกัน
martinstoeckli

64

ฉันพบรหัสนี้ซึ่งให้ผลลัพธ์ที่เชื่อถือได้กับฉัน

function distance($lat1, $lon1, $lat2, $lon2, $unit) {

  $theta = $lon1 - $lon2;
  $dist = sin(deg2rad($lat1)) * sin(deg2rad($lat2)) +  cos(deg2rad($lat1)) * cos(deg2rad($lat2)) * cos(deg2rad($theta));
  $dist = acos($dist);
  $dist = rad2deg($dist);
  $miles = $dist * 60 * 1.1515;
  $unit = strtoupper($unit);

  if ($unit == "K") {
      return ($miles * 1.609344);
  } else if ($unit == "N") {
      return ($miles * 0.8684);
  } else {
      return $miles;
  }
}

ผล :

echo distance(32.9697, -96.80322, 29.46786, -98.53506, "M") . " Miles<br>";
echo distance(32.9697, -96.80322, 29.46786, -98.53506, "K") . " Kilometers<br>";
echo distance(32.9697, -96.80322, 29.46786, -98.53506, "N") . " Nautical Miles<br>";

2
สิ่งที่ดีผมพยายามนี้และยัง google แผนที่แสดงระยะทางเดียวกันการเปลี่ยนแปลงเพียงทศนิยมที่นี่และมี ..
Zohair

ถ้าคุณต้องการคำนวณระยะทางระหว่างจุดสามจุด
kexxcream

3
เรียกใช้ฟังก์ชั่นนี้สองครั้งและเพิ่มพวกเขาอีกวิธีหนึ่งคุณเปลี่ยนฟังก์ชั่นตามนั้น
Janith Chinthana

ส่งคืน NaN ในบางเงื่อนไขstackoverflow.com/questions/37184259/…
Zahur Sh

23

มันเป็นเพียงส่วนเติมจาก@martinstoeckliและ@Janith Chinthanaตอบ สำหรับผู้ที่อยากรู้อยากเห็นเกี่ยวกับการที่อัลกอริทึมที่เร็วที่สุดที่ฉันเขียนการทดสอบประสิทธิภาพ ผลลัพธ์ที่ดีที่สุดแสดงฟังก์ชั่นที่ได้รับการปรับปรุงจากcodexworld.com :

/**
 * Optimized algorithm from http://www.codexworld.com
 *
 * @param float $latitudeFrom
 * @param float $longitudeFrom
 * @param float $latitudeTo
 * @param float $longitudeTo
 *
 * @return float [km]
 */
function codexworldGetDistanceOpt($latitudeFrom, $longitudeFrom, $latitudeTo, $longitudeTo)
{
    $rad = M_PI / 180;
    //Calculate distance from latitude and longitude
    $theta = $longitudeFrom - $longitudeTo;
    $dist = sin($latitudeFrom * $rad) 
        * sin($latitudeTo * $rad) +  cos($latitudeFrom * $rad)
        * cos($latitudeTo * $rad) * cos($theta * $rad);

    return acos($dist) / $rad * 60 *  1.853;
}

นี่คือผลการทดสอบ:

Test name       Repeats         Result          Performance     
codexworld-opt  10000           0.084952 sec    +0.00%
codexworld      10000           0.104127 sec    -22.57%
custom          10000           0.107419 sec    -26.45%
custom2         10000           0.111576 sec    -31.34%
custom1         10000           0.136691 sec    -60.90%
vincenty        10000           0.165881 sec    -95.26%

ในรหัสของคุณตัวคูณสำหรับอัลกอริทึม codexworld คือ 1.852 ในขณะที่ต้นฉบับจริงคือ 1.1515 ทำไมนี้ ทำไมถึงแตกต่าง
GotBatteries

@GotBatteries mulitplier ดั้งเดิมสำหรับไมล์ ฟังก์ชันที่ปรับให้เหมาะสมจะส่งคืนผลลัพธ์เป็นกิโลเมตร 1.1515 * 1.609344 = 1.853. ขอขอบคุณแก้ไขที่ 1.853
Alexander Yancharuk

ทำไมคุณไม่ใช้ M_PI / 180 และ $ rad * 60 * 1.853 เป็นค่าคงที่เพื่อประสิทธิภาพที่ดีกว่า
Evren Yurtesen

@EvrenYurtesen เป็นความคิดที่ดีหากความสำคัญของคุณคือประสิทธิภาพ แต่การบำรุงรักษาและการอ่านจะกลายเป็นความซับซ้อนมากขึ้นฉันคิดว่า
Alexander Yancharuk

เพียงใส่ความคิดเห็นในบรรทัดก่อนหน้าและพูดว่า // M_PI / 180 ... ฯลฯ ฉันไม่รู้ว่าทำไมมันคงยากที่จะรักษา ไม่ใช่สิ่งที่คุณจะเปลี่ยน
Evren Yurtesen

10

นี่คือรหัสที่ง่ายและสมบูรณ์แบบสำหรับการคำนวณระยะทางระหว่างละติจูดและลองจิจูดสองอัน พบรหัสต่อไปนี้จากที่นี่ - http://www.codexworld.com/distance-between-two-addresses-google-maps-api-php/

$latitudeFrom = '22.574864';
$longitudeFrom = '88.437915';

$latitudeTo = '22.568662';
$longitudeTo = '88.431918';

//Calculate distance from latitude and longitude
$theta = $longitudeFrom - $longitudeTo;
$dist = sin(deg2rad($latitudeFrom)) * sin(deg2rad($latitudeTo)) +  cos(deg2rad($latitudeFrom)) * cos(deg2rad($latitudeTo)) * cos(deg2rad($theta));
$dist = acos($dist);
$dist = rad2deg($dist);
$miles = $dist * 60 * 1.1515;

$distance = ($miles * 1.609344).' km';

5

สำหรับคนที่ชอบสั้นและเร็วกว่า (ไม่ใช่การเรียก deg2rad ())

function circle_distance($lat1, $lon1, $lat2, $lon2) {
  $rad = M_PI / 180;
  return acos(sin($lat2*$rad) * sin($lat1*$rad) + cos($lat2*$rad) * cos($lat1*$rad) * cos($lon2*$rad - $lon1*$rad)) * 6371;// Kilometers
}

2

ลองสิ่งนี้จะให้ผลลัพธ์ที่ยอดเยี่ยม

function getDistance($point1_lat, $point1_long, $point2_lat, $point2_long, $unit = 'km', $decimals = 2) {
        // Calculate the distance in degrees
        $degrees = rad2deg(acos((sin(deg2rad($point1_lat))*sin(deg2rad($point2_lat))) + (cos(deg2rad($point1_lat))*cos(deg2rad($point2_lat))*cos(deg2rad($point1_long-$point2_long)))));

        // Convert the distance in degrees to the chosen unit (kilometres, miles or nautical miles)
        switch($unit) {
            case 'km':
                $distance = $degrees * 111.13384; // 1 degree = 111.13384 km, based on the average diameter of the Earth (12,735 km)
                break;
            case 'mi':
                $distance = $degrees * 69.05482; // 1 degree = 69.05482 miles, based on the average diameter of the Earth (7,913.1 miles)
                break;
            case 'nmi':
                $distance =  $degrees * 59.97662; // 1 degree = 59.97662 nautic miles, based on the average diameter of the Earth (6,876.3 nautical miles)
        }
        return round($distance, $decimals);
    }

2

ค่อนข้างเก่าคำถาม แต่สำหรับผู้ที่สนใจในรหัส PHP ที่ส่งคืนผลลัพธ์เดียวกันกับ Google Maps ต่อไปนี้ทำงาน:

/**
 * Computes the distance between two coordinates.
 *
 * Implementation based on reverse engineering of
 * <code>google.maps.geometry.spherical.computeDistanceBetween()</code>.
 *
 * @param float $lat1 Latitude from the first point.
 * @param float $lng1 Longitude from the first point.
 * @param float $lat2 Latitude from the second point.
 * @param float $lng2 Longitude from the second point.
 * @param float $radius (optional) Radius in meters.
 *
 * @return float Distance in meters.
 */
function computeDistance($lat1, $lng1, $lat2, $lng2, $radius = 6378137)
{
    static $x = M_PI / 180;
    $lat1 *= $x; $lng1 *= $x;
    $lat2 *= $x; $lng2 *= $x;
    $distance = 2 * asin(sqrt(pow(sin(($lat1 - $lat2) / 2), 2) + cos($lat1) * cos($lat2) * pow(sin(($lng1 - $lng2) / 2), 2)));

    return $distance * $radius;
}

ฉันทดสอบด้วยพิกัดต่าง ๆ และมันทำงานได้อย่างสมบูรณ์

ฉันคิดว่ามันน่าจะเร็วกว่านี้แล้วมีทางเลือกอื่นด้วย แต่ไม่ได้ทำการทดสอบ

คำแนะนำ: Google Maps ใช้ 6378137 เป็นรัศมีของโลก ดังนั้นการใช้กับอัลกอริทึมอื่นอาจใช้งานได้เช่นกัน


1

สำหรับค่าที่แน่นอนทำเช่นนั้น:

public function DistAB()
{
      $delta_lat = $this->lat_b - $this->lat_a ;
      $delta_lon = $this->lon_b - $this->lon_a ;

      $a = pow(sin($delta_lat/2), 2);
      $a += cos(deg2rad($this->lat_a9)) * cos(deg2rad($this->lat_b9)) * pow(sin(deg2rad($delta_lon/29)), 2);
      $c = 2 * atan2(sqrt($a), sqrt(1-$a));

      $distance = 2 * $earth_radius * $c;
      $distance = round($distance, 4);

      $this->measure = $distance;
}

อืมฉันคิดว่าควรจะทำ ...

แก้ไข:

สำหรับ formulars และอย่างน้อย JS-implementations ลอง: http://www.movable-type.co.uk/scripts/latlong.html

กล้าฉัน ... ฉันลืม deg2rad ค่าทั้งหมดในฟังก์ชันวงกลม ...


ขอบคุณสำหรับคำตอบ. ฉันได้ตรวจสอบการใช้งานนี้ด้วยการคำนวณอย่างง่ายระหว่างpointA (42,12) และ pointB (43,12)โดยใช้ $ earth_radius = 6372.795477598 ฉันได้รับผลลัพธ์ 12745.591 เมื่อมันควรจะอยู่ที่ประมาณ 110,94
maxdangelo

1

สวัสดีรหัสที่นี่เพื่อรับระยะทางและเวลาโดยใช้ Lat และ Long สองแบบที่แตกต่างกัน

$url ="https://maps.googleapis.com/maps/api/distancematrix/json?units=imperial&origins=16.538048,80.613266&destinations=23.0225,72.5714";



    $ch = curl_init();
    // Disable SSL verification

    curl_setopt($ch, CURLOPT_SSL_VERIFYPEER, false);
    // Will return the response, if false it print the response
    curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
    // Set the url
    curl_setopt($ch, CURLOPT_URL,$url);
    // Execute
    $result=curl_exec($ch);
    // Closing
    curl_close($ch);

    $result_array=json_decode($result);
print_r($result_array);

คุณสามารถตรวจสอบตัวอย่างด้านล่างลิงก์รับเวลาระหว่างสถานที่สองแห่งโดยใช้ละติจูดและลองจิจูดใน php


6
อาจไม่จำเป็นต้องเรียก api สำหรับบางสิ่งที่สามารถพบได้โดยใช้คณิตศาสตร์
Ivotje50

1

ลองใช้ฟังก์ชันนี้เพื่อคำนวณระยะทางระหว่างถึงจุดละติจูดและลองจิจูด

function calculateDistanceBetweenTwoPoints($latitudeOne='', $longitudeOne='', $latitudeTwo='', $longitudeTwo='',$distanceUnit ='',$round=false,$decimalPoints='')
    {
        if (empty($decimalPoints)) 
        {
            $decimalPoints = '3';
        }
        if (empty($distanceUnit)) {
            $distanceUnit = 'KM';
        }
        $distanceUnit = strtolower($distanceUnit);
        $pointDifference = $longitudeOne - $longitudeTwo;
        $toSin = (sin(deg2rad($latitudeOne)) * sin(deg2rad($latitudeTwo))) + (cos(deg2rad($latitudeOne)) * cos(deg2rad($latitudeTwo)) * cos(deg2rad($pointDifference)));
        $toAcos = acos($toSin);
        $toRad2Deg = rad2deg($toAcos);

        $toMiles  =  $toRad2Deg * 60 * 1.1515;
        $toKilometers = $toMiles * 1.609344;
        $toNauticalMiles = $toMiles * 0.8684;
        $toMeters = $toKilometers * 1000;
        $toFeets = $toMiles * 5280;
        $toYards = $toFeets / 3;


              switch (strtoupper($distanceUnit)) 
              {
                  case 'ML'://miles
                         $toMiles  = ($round == true ? round($toMiles) : round($toMiles, $decimalPoints));
                         return $toMiles;
                      break;
                  case 'KM'://Kilometers
                        $toKilometers  = ($round == true ? round($toKilometers) : round($toKilometers, $decimalPoints));
                        return $toKilometers;
                      break;
                  case 'MT'://Meters
                        $toMeters  = ($round == true ? round($toMeters) : round($toMeters, $decimalPoints));
                        return $toMeters;
                      break;
                  case 'FT'://feets
                        $toFeets  = ($round == true ? round($toFeets) : round($toFeets, $decimalPoints));
                        return $toFeets;
                      break;
                  case 'YD'://yards
                        $toYards  = ($round == true ? round($toYards) : round($toYards, $decimalPoints));
                        return $toYards;
                      break;
                  case 'NM'://Nautical miles
                        $toNauticalMiles  = ($round == true ? round($toNauticalMiles) : round($toNauticalMiles, $decimalPoints));
                        return $toNauticalMiles;
                      break;
              }


    }

จากนั้นใช้ fucntion เป็น

echo calculateDistanceBetweenTwoPoints('11.657740','77.766270','11.074820','77.002160','ML',true,5);

หวังว่ามันจะช่วย


ตรวจสอบกับสถานการณ์จริงทำงานที่สมบูรณ์แบบในกรณีของฉัน
Daxesh Vekariya

1
ใช้เวลาเกือบ 5 ชั่วโมงในการเขียนและตรวจสอบในสถานการณ์จริง
Manojkiran.A

0

ตัวคูณจะถูกเปลี่ยนในทุกพิกัดเนื่องจากทฤษฎีระยะทางวงกลมใหญ่ดังที่เขียนไว้ที่นี่:

http://en.wikipedia.org/wiki/Great-circle_distance

และคุณสามารถคำนวณค่าที่ใกล้ที่สุดโดยใช้สูตรนี้อธิบายที่นี่:

http://en.wikipedia.org/wiki/Great-circle_distance#Worked_example

กุญแจสำคัญคือการแปลงค่าแต่ละองศา - นาที - วินาทีให้เป็นค่าทั้งหมด:

N 36°7.2', W 86°40.2'  N = (+) , W = (-), S = (-), E = (+) 
referencing the Greenwich meridian and Equator parallel

(phi)     36.12° = 36° + 7.2'/60' 

(lambda)  -86.67° = 86° + 40.2'/60'

0

หนึ่งในวิธีที่ง่ายที่สุดคือ:

$my_latitude = "";
$my_longitude = "";
$her_latitude = "";
$her_longitude = "";

$distance = round((((acos(sin(($my_latitude*pi()/180)) * sin(($her_latitude*pi()/180))+cos(($my_latitude*pi()/180)) * cos(($her_latitude*pi()/180)) * cos((($my_longitude- $her_longitude)*pi()/180))))*180/pi())*60*1.1515*1.609344), 2);
echo $distance;

มันจะปัดเศษทศนิยมไม่เกิน 2

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.