X-Git-Url: https://codewiz.org/gitweb?a=blobdiff_plain;f=examples%2Fdevelgps%2Fcompass.c;fp=examples%2Fdevelgps%2Fcompass.c;h=e1edcfba4eec50e7c6b52d5588e48c70302e62c9;hb=ab03b5b37726fb8b39a6f6c915f8668bddd779ee;hp=0000000000000000000000000000000000000000;hpb=ec4804af8993aaf53193c71abdf890f7283fe5fb;p=bertos.git

diff --git a/examples/develgps/compass.c b/examples/develgps/compass.c
new file mode 100644
index 00000000..e1edcfba
--- /dev/null
+++ b/examples/develgps/compass.c
@@ -0,0 +1,99 @@
+/**
+ * \file
+ * <!--
+ * This file is part of BeRTOS.
+ *
+ * Bertos is free software; you can redistribute it and/or modify
+ * it under the terms of the GNU General Public License as published by
+ * the Free Software Foundation; either version 2 of the License, or
+ * (at your option) any later version.
+ *
+ * This program is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
+ *
+ * As a special exception, you may use this file as part of a free software
+ * library without restriction.  Specifically, if other files instantiate
+ * templates or use macros or inline functions from this file, or you compile
+ * this file and link it with other files to produce an executable, this
+ * file does not by itself cause the resulting executable to be covered by
+ * the GNU General Public License.  This exception does not however
+ * invalidate any other reasons why the executable file might be covered by
+ * the GNU General Public License.
+ *
+ * Copyright 2010 Develer S.r.l. (http://www.develer.com/)
+ *
+ * -->
+ *
+ * \brief DevelGPS: compass routines.
+ *
+ * \author Andrea Righi <arighi@develer.com>
+ */
+
+#include <cfg/compiler.h>
+#include <cfg/debug.h>
+#include <math.h>
+#include "compass.h"
+
+static const char *_compass_heading[] =
+{
+	"N", "NNE", "NE", "ENE",
+	"E", "ESE", "SE", "SSE",
+	"S", "SSW", "SW", "WSW",
+	"W", "WNW", "NW", "NNW",
+};
+
+/**
+ * Use the Haversine formula to calculate great-circle distances between the
+ * two points.
+ *
+ * The Haversine formula remains particularly well-conditioned for numerical
+ * computation even at small distances, unlike calculations based on the
+ * spherical law of cosines.
+ */
+float distance(float lat1, float lon1, float lat2, float lon2)
+{
+	const float PLANET_RADIUS = 6371000;
+	float d_lat = deg2rad(lat2 - lat1);
+	float d_lon = deg2rad(lon2 - lon1);
+
+	float a = sin(d_lat / 2) * sin(d_lat / 2) +
+			cos(deg2rad(lat1)) * cos(deg2rad(lat2)) *
+			sin(d_lon / 2) * sin(d_lon / 2);
+	float c = 2 * atan2(sqrt(a), sqrt(1 - a));
+
+	return PLANET_RADIUS * c;
+}
+
+/**
+ * Evaluate the bearing (also known as forward azimuth) using spherical law
+ * coordinates.
+ *
+ * The bearing is a straight line along a great-circle arc from the start point
+ * to the destination point.
+ */
+int bearing(float lat1, float lon1, float lat2, float lon2)
+{
+	float res;
+
+	res = rad2deg(atan2(sin(deg2rad(lon2 - lon1)) *
+		cos(deg2rad(lat2)), cos(deg2rad(lat1)) *
+		sin(deg2rad(lat2)) - sin(deg2rad(lat1)) *
+		cos(deg2rad(lat2)) * cos(deg2rad(lon2) -
+		deg2rad(lon1))));
+	return ((int)res + 360) % 360;
+}
+
+const char *compass_heading(int bearing)
+{
+	ASSERT(bearing >= 0 && bearing < 360);
+	/*
+	 * bearing / 22.5
+	 */
+	return _compass_heading[(bearing << 4) / 360];
+}