Measuring the distance between two coordinates in PHP

Question

Hi I have the need to calculate the distance between two points having the lat and long   I would like to avoid any call to external API   I tried to implement the Haversine Formula in PHP   Here is the code   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- gt lat b -  this- gt lat a              delta lon    this- gt lon b -  this- gt lon a               earth radius   6372 795477598              alpha       delta lat 2             beta        delta lon 2             a          sin deg2rad  alpha     sin deg2rad  alpha     cos deg2rad  this- gt lat a     cos deg2rad  this- gt lat b     sin deg2rad  beta     sin deg2rad  beta                c          asin min 1  sqrt  a                distance   2  earth radius    c             distance   round  distance  4               this- gt measure    distance                   Testing it with some given points which have public distances I don t get a reliable result   I don t understand if there is an error in the original formula or in my implementation

User · Answer

For exact values do it like that   public function DistAB            delta lat    this- gt lat b -  this- gt lat a          delta lon    this- gt lon b -  this- gt lon a           a   pow sin  delta lat 2   2          a    cos deg2rad  this- gt lat a9     cos deg2rad  this- gt 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- gt measure    distance      Hmm I think that should do it     Edit   For formulars and at least JS-implementations try  http   www movable-type co uk scripts latlong html  Dare me    I forgot to deg2rad all the values in the circle-functions

User · Answer

I found this code which is giving me reliable results   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          results    echo distance 32 9697  -96 80322  29 46786  -98 53506   M       Miles lt br gt    echo distance 32 9697  -96 80322  29 46786  -98 53506   K       Kilometers lt br gt    echo distance 32 9697  -96 80322  29 46786  -98 53506   N       Nautical Miles lt br gt

User · Answer

For the ones who like shorter and faster not calling 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

User · Answer

Try this function out to calculate distance between to points of latitude and longitude  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                            Then use the fucntion as  echo calculateDistanceBetweenTwoPoints  11 657740   77 766270   11 074820   77 002160   ML  true 5     Hope it helps

User · Answer

The multiplier is changed at every coordinate because of the great circle distance theory as written here     http   en wikipedia org wiki Great-circle distance  and you can calculate the nearest value using this formula described here   http   en wikipedia org wiki Great-circle distance Worked example  the key is converting each degree - minute - second value to all degree value   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

User · Answer

It s just addition to  martinstoeckli and  Janith Chinthana answers  For those who curious about which algorithm is fastest i wrote the performance test  Best performance result shows optimized function from 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      Here is test results   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

User · Answer

Try this gives awesome results      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

User · Answer

Not long ago I wrote an example of the haversine formula  and published it on my website          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        Note that you get the distance back in the same unit as you pass in with the parameter  earthRadius  The default value is 6371000 meters so the result will be in  m  too  To get the result in miles  you could e g  pass 3959 miles as  earthRadius and the result would be in  mi   In my opinion it is a good habit to stick with the SI units  if there is no particular reason to do otherwise   Edit   As TreyA correctly pointed out  the Haversine formula has weaknesses with antipodal points because of rounding errors  though it is stable for small distances   To get around them  you could use the Vincenty formula instead          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

User · Answer

Here the simple and perfect code for calculating the distance between two latitude and longitude  The following code have been found from here - 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

User · Answer

Hello here Code For Get Distance and Time Using Two Different Lat and Long     url   https   maps googleapis com maps api distancematrix json units imperial amp origins 16 538048 80 613266 amp 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     You can check Example Below Link get time between two different locations using latitude and longitude in php

User · Answer

One of the easiest ways is    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    It will round off up to 2 decimal points

User · Answer

Quite old question  but for those interested in a PHP code that returns the same results as Google Maps  the following does the job         Computes the distance between two coordinates        Implementation based on reverse engineering of     lt code gt google maps geometry spherical computeDistanceBetween   lt  code gt          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     I ve tested with various coordinates and it works perfectly  I think it should be faster then some alternatives too  But didn t test that  Hint  Google Maps uses 6378137 as Earth radius  So using it with other algorithms might work as well

[php] Measuring the distance between two coordinates in PHP

Examples related to php

Examples related to coordinates

Examples related to geocoding

Examples related to haversine