spirit posteets tagged distance  [ Profile ]

Sort by: Date / Title /

  1. 10 years ago
    Great circle distance
    Great circle distance
    Find the distance in kilometres between two points on the surface of the earth. This is just the sort of problem stored functions were made for. For a first order approximation, ignore deviations of the earth's surface from the perfectly spherical. Then the distance in radians is given by a number of trigonometric formulas. ACOS and COS behave reasonably:
    
                 COS(lat1-lat2)*(1+COS(lon1-lon2)) - COS(lat1+lat2)*(1-COS(lon1-lon2))
    rads = ACOS( --------------------------------------------------------------------- )
                                                  2
    
    We need to convert degrees latitude and longitude to radians, and we need to know the length in km of one radian on the earth's surface, which is 6378.388. The function:
    
    set log_bin_trust_function_creators=TRUE;
    
    DROP FUNCTION IF EXISTS GeoDistKM;
    DELIMITER |
    CREATE FUNCTION GeoDistKM( lat1 FLOAT, lon1 FLOAT, lat2 FLOAT, lon2 FLOAT ) RETURNS float
    BEGIN
      DECLARE pi, q1, q2, q3 FLOAT;
      DECLARE rads FLOAT DEFAULT 0;
      SET pi = PI();
      SET lat1 = lat1 * pi / 180;
      SET lon1 = lon1 * pi / 180;
      SET lat2 = lat2 * pi / 180;
      SET lon2 = lon2 * pi / 180;
      SET q1 = COS(lon1-lon2);
      SET q2 = COS(lat1-lat2);
      SET q3 = COS(lat1+lat2);
      SET rads = ACOS( 0.5*((1.0+q1)*q2 - (1.0-q1)*q3) ); 
      RETURN 6378.388 * rads;
    END;
    |
    DELIMITER ;
    
    -- toronto to montreal (505km):
    select geodistkm(43.6667,-79.4167,45.5000,-73.5833);
    +----------------------------------------------+
    | geodistkm(43.6667,-79.4167,45.5000,-73.5833) |
    +----------------------------------------------+
    |                           505.38836669921875 |
    +----------------------------------------------+
    Paste this in your website: <script type="text/javascript" src="http://www.posteet.com/embed/1555"></script>
  2. 10 years ago
    1. function distance($lat1, $lon1, $lat2, $lon2, $unit) {
    2.  
    3.   $theta = $lon1 - $lon2;
    4.   $dist = sin(deg2rad($lat1)) * sin(deg2rad($lat2))cos(deg2rad($lat1)) * cos(deg2rad($lat2)) * cos(deg2rad($theta));
    5.   $dist = acos($dist);
    6.   $dist = rad2deg($dist);
    7.   $miles = $dist * 60 * 1.1515;
    8.   $unit = strtoupper($unit);
    9.  
    10.   if ($unit == "K") {
    11.     return ($miles * 1.609344);
    12.   } else if ($unit == "N") {
    13.       return ($miles * 0.8684);
    14.     } else {
    15.         return $miles;
    16.       }
    17. }
    18.  
    19. echo distance(32.9697, -96.80322, 29.46786, -98.53506, "m") . " miles<br>";
    20. echo distance(32.9697, -96.80322, 29.46786, -98.53506, "k") . " kilometers<br>";
    21. echo distance(32.9697, -96.80322, 29.46786, -98.53506, "n") . " nautical miles<br>";
    Paste this in your website: <script type="text/javascript" src="http://www.posteet.com/embed/1554"></script>

First / Previous / Next / Last / Page 1 of 1 (2 posteets)