Add OpenLayers/Spherical.js

lib/OpenLayers/Spherical.js

@@ -0,0 +1,59 @@
+/* Copyright (c) 2006-2012 by OpenLayers Contributors (see authors.txt for
+ * full list of contributors). Published under the Clear BSD license.
+ * See http://svn.openlayers.org/trunk/openlayers/license.txt for the
+ * full text of the license. */
+
+/**
+ * @see http://www.movable-type.co.uk/scripts/latlong.html
+ * @see http://code.google.com/apis/maps/documentation/javascript/reference.html#spherical
+ */
+
+
+/**
+ * Namespace: Spherical
+ */
+OpenLayers.Spherical = OpenLayers.Spherical || {};
+
+
+OpenLayers.Spherical.DEFAULT_RADIUS = 6378137;
+
+
+/**
+ * APIFunction: computeDistanceBetween
+ * Computes the distance between two LonLats.
+ *
+ * Parameters:
+ * from   - {<OpenLayers.LonLat>}
+ * to     - {<OpenLayers.LonLat>}
+ * radius - {Float}
+ *
+ * Returns:
+ * {Float} The distance in meters.
+ */
+OpenLayers.Spherical.computeDistanceBetween = function(from, to, radius) {
+  var R = radius || OpenLayers.Spherical.DEFAULT_RADIUS;
+  var sinHalfDeltaLon = Math.sin(Math.PI * (to.lon - from.lon) / 360);
+  var sinHalfDeltaLat = Math.sin(Math.PI * (to.lat - from.lat) / 360);
+  var a = sinHalfDeltaLat * sinHalfDeltaLat +
+      sinHalfDeltaLon * sinHalfDeltaLon * Math.cos(Math.PI * from.lat / 180) * Math.cos(Math.PI * to.lat / 180); 
+  return 2 * R * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a)); 
+};
+
+
+/**
+ * APIFunction: computeHeading
+ * Computes the heading from one LonLat to another LonLat.
+ *
+ * Parameters:
+ * from - {<OpenLayers.LonLat>}
+ * to   - {<OpenLayers.LonLat>}
+ *
+ * Returns:
+ * {Float} The heading in degrees.
+ */
+OpenLayers.Spherical.computeHeading = function(from, to) {
+    var y = Math.sin(Math.PI * (from.lon - to.lon) / 180) * Math.cos(Math.PI * to.lat / 180);
+    var x = Math.cos(Math.PI * from.lat / 180) * Math.sin(Math.PI * to.lat / 180) -
+        Math.sin(Math.PI * from.lat / 180) * Math.cos(Math.PI * to.lat / 180) * Math.cos(Math.PI * (from.lon - to.lon) / 180);
+    return 180 * Math.atan2(y, x) / Math.PI;
+};