XProfan

/**
* Converts degrees/minutes/seconds to degrees
*
* Converts a string which represents a latitude/longitude as degree/minutes/seconds
* to a float degree value. If no valid string is passed it will return 0.
*
* @access  public
* @param   string  $dms  latitude/longitude as degree/minutes/seconds
* @return  float   degree
*/

function dms2deg($dms, $language = GEO_LANGUAGE_DEFAULT) {

    global $cfgStrings;
     $negativeSigns = array($cfgStrings[GEO_ORIENTATION_SHORT][$language][4], $cfgStrings[GEO_ORIENTATION_SHORT][$language][6], -);
     $negativeSignsString = $cfgStrings[GEO_ORIENTATION_SHORT][$language][4].$cfgStrings[GEO_ORIENTATION_SHORT][$language][6];

    if (strlen($dms) == 6) {

         $dms = 0.$dms;

    } elseif (strlen($dms) == 5) {

         $dms = 00.$dms;

    }

     $searchPattern = |s*([$negativeSignsString-+]?)s*(d{1,3})[°s]*(d{1,2})[\s]*(d{1,2})([,.]*)(d*)[\s]*([$negativeSignsString-+]?)|i;

    if(preg_match($searchPattern, $dms, $result)) {

        if (in_array(strtoupper($result[1]), $negativeSigns) || in_array(strtoupper($result[7]), $negativeSigns)) {

             $algSign = -1.;

        } else {

             $algSign = 1.;

        }

        if (((1. * $result[2]) > 360) || ($result[3] >= 60) || ($result[4] >= 60)) {

            return PEAR::raiseError(Values out of range);

        }

        return $algSign * ($result[2] +(($result[3] + (($result[4]...$result[6]) * 10/6)/100)*10/6)/100);

    } else {

        return PEAR::raiseError(No DMS-Format (Like 51° 24 32.123\ W));

    }

}

/**
* Converts a float value to degrees/minutes/seconds
*
* Converts a float value to degrees/minutes/second (e.g. 50.1833300 to 50° 10 60)
* The seconds could contain the number of decimal places one passes to the optional
* parameter $decPlaces. The direction (N, S, W, E) must be added manually
* (e.g. $output = E .deg2dms(7.441944); )
*
* @access  public
* @param   float   $degFloat
* @param   int     $decPlaces
* @return  string  degrees minutes seconds
*/

function deg2dms($degFloat, $decPlaces = 0) {

     $degree = abs(floor($degFloat));
     $minSec = 60 * ($degFloat - $degree);
     $minutes = floor($minSec);
     $seconds = round(($minSec - $minutes) * 60, $decPlaces);
    return $degree.° $minutes $seconds;

}

/**
* Returns the radius of the earth
*
* Returns the radius of the earth in the given unit.
* GEO_EARTH_RADIUS is set to the mean value: 6371. km
* equatorial radius as of WGS84: 6378.137 km
*
* @access  public
* @param   int     $unit  use the GEO_UNIT_* constants
*/

function getEarthRadius($unit = GEO_UNIT_DEFAULT) {

    switch ($unit) {

        case GEO_UNIT_KM:// Kilometer
        return GEO_EARTH_RADIUS;
        case GEO_UNIT_MI:// Meilen
        return GEO_EARTH_RADIUS * 0.621371;
        case GEO_UNIT_IN:// Zoll/Inch
        return GEO_EARTH_RADIUS * 39370.08;
        case GEO_UNIT_SM:// See-Meilen
        return GEO_EARTH_RADIUS * 0.5399568;
        default:
        return GEO_EARTH_RADIUS;

    }

}