You cannot select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
241 lines
7.9 KiB
PHP
241 lines
7.9 KiB
PHP
<?php
|
|
namespace AnthonyMartin\GeoLocation;
|
|
|
|
/**
|
|
* <p>Represents a point on the surface of a sphere. (The Earth is almost
|
|
* spherical.)</p>
|
|
*
|
|
* <p>To create an instance, call one of the static methods fromDegrees() or
|
|
* fromRadians().</p>
|
|
*
|
|
* <p>This is a PHP port of Java code that was originally published at
|
|
* <a href="http://JanMatuschek.de/LatitudeLongitudeBoundingCoordinates#Java">
|
|
* http://JanMatuschek.de/LatitudeLongitudeBoundingCoordinates#Java</a>.</p>
|
|
*
|
|
* Many thanks to the original author: Jan Philip Matuschek
|
|
*
|
|
* @author Anthony Martin
|
|
* @version November 21 2012
|
|
*/
|
|
class GeoLocation {
|
|
|
|
private $radLat; // latitude in radians
|
|
private $radLon; // longitude in radians
|
|
|
|
private $degLat; // latitude in degrees
|
|
private $degLon; // longitude in degrees
|
|
|
|
private $angular; // angular radius
|
|
|
|
const EARTHS_RADIUS_KM = 6371.01;
|
|
const EARTHS_RADIUS_MI = 3958.762079;
|
|
|
|
protected static $MIN_LAT; // -PI/2
|
|
protected static $MAX_LAT; // PI/2
|
|
protected static $MIN_LON; // -PI
|
|
protected static $MAX_LON; // PI
|
|
|
|
public function __construct() {
|
|
self::$MIN_LAT = deg2rad(-90); // -PI/2
|
|
self::$MAX_LAT = deg2rad(90); // PI/2
|
|
self::$MIN_LON = deg2rad(-180); // -PI
|
|
self::$MAX_LON = deg2rad(180); // PI
|
|
}
|
|
|
|
/**
|
|
* @param double $latitude the latitude, in degrees.
|
|
* @param double $longitude the longitude, in degrees.
|
|
* @return GeoLocation
|
|
*/
|
|
public static function fromDegrees($latitude, $longitude) {
|
|
$location = new GeoLocation();
|
|
$location->radLat = deg2rad($latitude);
|
|
$location->radLon = deg2rad($longitude);
|
|
$location->degLat = $latitude;
|
|
$location->degLon = $longitude;
|
|
$location->checkBounds();
|
|
return $location;
|
|
}
|
|
|
|
/**
|
|
* @param double $latitude the latitude, in radians.
|
|
* @param double $longitude the longitude, in radians.
|
|
* @return GeoLocation
|
|
*/
|
|
public static function fromRadians($latitude, $longitude) {
|
|
$location = new GeoLocation();
|
|
$location->radLat = $latitude;
|
|
$location->radLon = $longitude;
|
|
$location->degLat = rad2deg($latitude);
|
|
$location->degLon = rad2deg($longitude);
|
|
$location->checkBounds();
|
|
return $location;
|
|
}
|
|
|
|
protected function checkBounds() {
|
|
if ($this->radLat < self::$MIN_LAT || $this->radLat > self::$MAX_LAT ||
|
|
$this->radLon < self::$MIN_LON || $this->radLon > self::$MAX_LON)
|
|
throw new \Exception("Invalid Argument");
|
|
}
|
|
|
|
/**
|
|
* Computes the great circle distance between this GeoLocation instance
|
|
* and the location argument.
|
|
* @param GeoLocation $location
|
|
* @param string $unit_of_measurement
|
|
* @internal param float $radius the radius of the sphere, e.g. the average radius for a
|
|
* spherical approximation of the figure of the Earth is approximately
|
|
* 6371.01 kilometers.
|
|
* @return double the distance, measured in the same unit as the radius
|
|
* argument.
|
|
*/
|
|
public function distanceTo(GeoLocation $location, $unit_of_measurement) {
|
|
$radius = $this->getEarthsRadius($unit_of_measurement);
|
|
|
|
return acos(sin($this->radLat) * sin($location->radLat) +
|
|
cos($this->radLat) * cos($location->radLat) *
|
|
cos($this->radLon - $location->radLon)) * $radius;
|
|
}
|
|
|
|
/**
|
|
* @return double the latitude, in degrees.
|
|
*/
|
|
public function getLatitudeInDegrees() {
|
|
return $this->degLat;
|
|
}
|
|
|
|
/**
|
|
* @return double the longitude, in degrees.
|
|
*/
|
|
public function getLongitudeInDegrees() {
|
|
return $this->degLon;
|
|
}
|
|
|
|
/**
|
|
* @return double the latitude, in radians.
|
|
*/
|
|
public function getLatitudeInRadians() {
|
|
return $this->radLat;
|
|
}
|
|
|
|
/**
|
|
* @return double the longitude, in radians.
|
|
*/
|
|
public function getLongitudeInRadians() {
|
|
return $this->radLon;
|
|
}
|
|
|
|
/**
|
|
* @return double angular radius.
|
|
*/
|
|
public function getAngular() {
|
|
return $this->angular;
|
|
}
|
|
|
|
public function __toString() {
|
|
return "(" . $this->degLat . ", " . $this->degLon . ") = (" .
|
|
$this->radLat . " rad, " . $this->radLon . " rad";
|
|
}
|
|
|
|
|
|
/**
|
|
* <p>Computes the bounding coordinates of all points on the surface
|
|
* of a sphere that have a great circle distance to the point represented
|
|
* by this GeoLocation instance that is less or equal to the distance
|
|
* argument.</p>
|
|
* <p>For more information about the formulae used in this method visit
|
|
* <a href="http://JanMatuschek.de/LatitudeLongitudeBoundingCoordinates">
|
|
* http://JanMatuschek.de/LatitudeLongitudeBoundingCoordinates</a>.</p>
|
|
*
|
|
* @param double $distance the distance from the point represented by this
|
|
* GeoLocation instance. Must me measured in the same unit as the radius
|
|
* argument.
|
|
* @param string $unit_of_measurement
|
|
* @throws \Exception
|
|
* @internal param radius the radius of the sphere, e.g. the average radius for a
|
|
* spherical approximation of the figure of the Earth is approximately
|
|
* 6371.01 kilometers.
|
|
* @return GeoLocation[] an array of two GeoLocation objects such that:<ul>
|
|
* <li>The latitude of any point within the specified distance is greater
|
|
* or equal to the latitude of the first array element and smaller or
|
|
* equal to the latitude of the second array element.</li>
|
|
* <li>If the longitude of the first array element is smaller or equal to
|
|
* the longitude of the second element, then
|
|
* the longitude of any point within the specified distance is greater
|
|
* or equal to the longitude of the first array element and smaller or
|
|
* equal to the longitude of the second array element.</li>
|
|
* <li>If the longitude of the first array element is greater than the
|
|
* longitude of the second element (this is the case if the 180th
|
|
* meridian is within the distance), then
|
|
* the longitude of any point within the specified distance is greater
|
|
* or equal to the longitude of the first array element
|
|
* <strong>or</strong> smaller or equal to the longitude of the second
|
|
* array element.</li>
|
|
* </ul>
|
|
*/
|
|
public function boundingCoordinates($distance, $unit_of_measurement) {
|
|
$radius = $this->getEarthsRadius($unit_of_measurement);
|
|
if ($radius < 0 || $distance < 0) throw new \Exception('Arguments must be greater than 0.');
|
|
|
|
// angular distance in radians on a great circle
|
|
$this->angular = $distance / $radius;
|
|
|
|
$minLat = $this->radLat - $this->angular;
|
|
$maxLat = $this->radLat + $this->angular;
|
|
|
|
$minLon = 0;
|
|
$maxLon = 0;
|
|
if ($minLat > self::$MIN_LAT && $maxLat < self::$MAX_LAT) {
|
|
$deltaLon = asin(sin($this->angular) /
|
|
cos($this->radLat));
|
|
$minLon = $this->radLon - $deltaLon;
|
|
if ($minLon < self::$MIN_LON) $minLon += 2 * pi();
|
|
$maxLon = $this->radLon + $deltaLon;
|
|
if ($maxLon > self::$MAX_LON) $maxLon -= 2 * pi();
|
|
} else {
|
|
// a pole is within the distance
|
|
$minLat = max($minLat, self::$MIN_LAT);
|
|
$maxLat = min($maxLat, self::$MAX_LAT);
|
|
$minLon = self::$MIN_LON;
|
|
$maxLon = self::$MAX_LON;
|
|
}
|
|
|
|
return array(
|
|
GeoLocation::fromRadians($minLat, $minLon),
|
|
GeoLocation::fromRadians($maxLat, $maxLon)
|
|
);
|
|
}
|
|
|
|
protected function getEarthsRadius($unit_of_measurement) {
|
|
$u = $unit_of_measurement;
|
|
if($u == 'miles' || $u == 'mi')
|
|
return $radius = self::EARTHS_RADIUS_MI;
|
|
elseif($u == 'kilometers' || $u == 'km')
|
|
return $radius = self::EARTHS_RADIUS_KM;
|
|
|
|
else throw new \Exception('You must supply a valid unit of measurement');
|
|
}
|
|
/**
|
|
* Retrieves Geocoding information from Google
|
|
* eg. $response = GeoLocation::getGeocodeFromGoogle($location);
|
|
* $latitude = $response->results[0]->geometry->location->lng;
|
|
* $longitude = $response->results[0]->geometry->location->lng;
|
|
* @param string $location address, city, state, etc.
|
|
* @return \stdClass
|
|
*/
|
|
public static function getGeocodeFromGoogle($location) {
|
|
$url = 'http://maps.googleapis.com/maps/api/geocode/json?address='.urlencode($location).'&sensor=false';
|
|
$ch = curl_init();
|
|
curl_setopt($ch, CURLOPT_URL,$url);
|
|
curl_setopt($ch, CURLOPT_RETURNTRANSFER, true);
|
|
return json_decode(curl_exec($ch));
|
|
}
|
|
public static function MilesToKilometers($miles) {
|
|
return $miles * 1.6093439999999999;
|
|
}
|
|
public static function KilometersToMiles($km) {
|
|
return $km * 0.621371192237334;
|
|
}
|
|
}
|
|
|