/* vim: set et ts=8 sw=8: */
/* gclue-location.c
 *
 * Copyright 2012 Bastien Nocera
 * Copyright 2015 Ankit (Verma)
 *
 * Geoclue 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.
 *
 * Geoclue 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 Geoclue; if not, write to the Free Software Foundation, Inc.,
 * 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
 *
 *    Authors: Bastien Nocera <hadess@hadess.net>
 *             Zeeshan Ali (Khattak) <zeeshanak@gnome.org>
 *             Ankit (Verma) <ankitstarski@gmail.com>
 */

#include "gclue-location.h"
#include "gclue-nmea-utils.h"
#include <math.h>
#include <string.h>
#include <stdlib.h>

#define TIME_DIFF_THRESHOLD (60 * G_USEC_PER_SEC) /* 60 seconds */
#define EARTH_RADIUS_KM 6372.795
#define KNOTS_IN_METERS_PER_SECOND 0.51444
#define RMC_TIME_DIFF_THRESHOLD 5 /* 5 seconds */
#define RMC_DEFAULT_ACCURACY 5    /* 5 meters */

struct _GClueLocationPrivate {
        char   *description;

        gdouble longitude;
        gdouble latitude;
        gdouble altitude;
        gdouble accuracy;
        guint64 timestamp;
        gdouble speed;
        gdouble heading;
};

enum {
        PROP_0,

        PROP_LATITUDE,
        PROP_LONGITUDE,
        PROP_ACCURACY,
        PROP_DESCRIPTION,
        PROP_TIMESTAMP,
        PROP_ALTITUDE,
        PROP_SPEED,
        PROP_HEADING,
};

G_DEFINE_TYPE_WITH_CODE (GClueLocation,
                         gclue_location,
                         G_TYPE_OBJECT,
                         G_ADD_PRIVATE (GClueLocation));

static void
gclue_location_get_property (GObject    *object,
                             guint       property_id,
                             GValue     *value,
                             GParamSpec *pspec)
{
        GClueLocation *location = GCLUE_LOCATION (object);

        switch (property_id) {
        case PROP_DESCRIPTION:
                g_value_set_string (value,
                                    gclue_location_get_description (location));
                break;

        case PROP_LATITUDE:
                g_value_set_double (value,
                                    gclue_location_get_latitude (location));
                break;

        case PROP_LONGITUDE:
                g_value_set_double (value,
                                    gclue_location_get_longitude (location));
                break;

        case PROP_ALTITUDE:
                g_value_set_double (value,
                                    gclue_location_get_altitude (location));
                break;

        case PROP_ACCURACY:
                g_value_set_double (value,
                                    gclue_location_get_accuracy (location));
                break;

        case PROP_TIMESTAMP:
                g_value_set_uint64 (value,
                                    gclue_location_get_timestamp (location));
                break;
        case PROP_SPEED:
                g_value_set_double (value,
                                    gclue_location_get_speed (location));
                break;

        case PROP_HEADING:
                g_value_set_double (value,
                                    gclue_location_get_heading (location));
                break;

        default:
                /* We don't have any other property... */
                G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
                break;
        }
}

static void
gclue_location_set_latitude (GClueLocation *loc,
                             gdouble        latitude)
{
        g_return_if_fail (latitude >= -90.0 && latitude <= 90.0);

        loc->priv->latitude = latitude;
}

static void
gclue_location_set_longitude (GClueLocation *loc,
                              gdouble        longitude)
{
        g_return_if_fail (longitude >= -180.0 && longitude <= 180.0);

        loc->priv->longitude = longitude;
}

static void
gclue_location_set_altitude (GClueLocation *loc,
                             gdouble        altitude)
{
        loc->priv->altitude = altitude;
}

void
gclue_location_set_accuracy (GClueLocation *loc,
                             gdouble        accuracy)
{
        g_return_if_fail (accuracy >= GCLUE_LOCATION_ACCURACY_UNKNOWN);

        loc->priv->accuracy = accuracy;
}

static void
gclue_location_set_timestamp (GClueLocation *loc,
                              guint64        timestamp)
{
        g_return_if_fail (GCLUE_IS_LOCATION (loc));

        loc->priv->timestamp = timestamp;
}

void
gclue_location_set_description (GClueLocation *loc,
                                const char    *description)
{
        g_return_if_fail (GCLUE_IS_LOCATION (loc));

        g_free (loc->priv->description);
        loc->priv->description = g_strdup (description);
}

static void
gclue_location_set_property (GObject      *object,
                             guint         property_id,
                             const GValue *value,
                             GParamSpec   *pspec)
{
        GClueLocation *location = GCLUE_LOCATION (object);

        switch (property_id) {
        case PROP_DESCRIPTION:
                gclue_location_set_description (location,
                                                g_value_get_string (value));
                break;

        case PROP_LATITUDE:
                gclue_location_set_latitude (location,
                                             g_value_get_double (value));
                break;

        case PROP_LONGITUDE:
                gclue_location_set_longitude (location,
                                              g_value_get_double (value));
                break;

        case PROP_ALTITUDE:
                gclue_location_set_altitude (location,
                                             g_value_get_double (value));
                break;

        case PROP_ACCURACY:
                gclue_location_set_accuracy (location,
                                             g_value_get_double (value));
                break;

        case PROP_TIMESTAMP:
                gclue_location_set_timestamp (location,
                                              g_value_get_uint64 (value));
                break;
        case PROP_SPEED:
                gclue_location_set_speed (location,
                                          g_value_get_double (value));
                break;

        case PROP_HEADING:
                gclue_location_set_heading (location,
                                            g_value_get_double (value));
                break;

        default:
                /* We don't have any other property... */
                G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
                break;
        }
}

static void
gclue_location_constructed (GObject *object)
{
        GClueLocation *location = GCLUE_LOCATION (object);
        gint64 timestamp;

        if (location->priv->timestamp != 0)
                return;

        timestamp = g_get_real_time () / G_USEC_PER_SEC;
        gclue_location_set_timestamp (location, timestamp);
}

static void
gclue_location_finalize (GObject *glocation)
{
        g_clear_pointer (&GCLUE_LOCATION (glocation)->priv->description,
                         g_free);

        G_OBJECT_CLASS (gclue_location_parent_class)->finalize (glocation);
}

static void
gclue_location_class_init (GClueLocationClass *klass)
{
        GObjectClass *glocation_class = G_OBJECT_CLASS (klass);
        GParamSpec *pspec;

        glocation_class->constructed = gclue_location_constructed;
        glocation_class->finalize = gclue_location_finalize;
        glocation_class->get_property = gclue_location_get_property;
        glocation_class->set_property = gclue_location_set_property;

        /**
         * GClueLocation:description:
         *
         * The description of this location.
         */
        pspec = g_param_spec_string ("description",
                                     "Description",
                                     "Description of this location",
                                     NULL,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_DESCRIPTION, pspec);

        /**
         * GClueLocation:latitude:
         *
         * The latitude of this location in degrees.
         */
        pspec = g_param_spec_double ("latitude",
                                     "Latitude",
                                     "The latitude of this location in degrees",
                                     -90.0,
                                     90.0,
                                     0.0,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_LATITUDE, pspec);

        /**
         * GClueLocation:longitude:
         *
         * The longitude of this location in degrees.
         */
        pspec = g_param_spec_double ("longitude",
                                     "Longitude",
                                     "The longitude of this location in degrees",
                                     -180.0,
                                     180.0,
                                     0.0,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_LONGITUDE, pspec);

        /**
         * GClueLocation:altitude:
         *
         * The altitude of this location in meters.
         */
        pspec = g_param_spec_double ("altitude",
                                     "Altitude",
                                     "The altitude of this location in meters",
                                     GCLUE_LOCATION_ALTITUDE_UNKNOWN,
                                     G_MAXDOUBLE,
                                     GCLUE_LOCATION_ALTITUDE_UNKNOWN,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_ALTITUDE, pspec);

        /**
         * GClueLocation:accuracy:
         *
         * The accuracy of this location in meters.
         */
        pspec = g_param_spec_double ("accuracy",
                                     "Accuracy",
                                     "The accuracy of this location in meters",
                                     GCLUE_LOCATION_ACCURACY_UNKNOWN,
                                     G_MAXDOUBLE,
                                     GCLUE_LOCATION_ACCURACY_UNKNOWN,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_ACCURACY, pspec);


        /**
         * GClueLocation:timestamp:
         *
         * A timestamp in seconds since
         * <ulink url="http://en.wikipedia.org/wiki/Unix_epoch">Epoch</ulink>,
         * giving when the location was resolved from an address.
         *
         * A value of 0 (zero) will be interpreted as the current time.
         */
        pspec = g_param_spec_uint64 ("timestamp",
                                     "Timestamp",
                                     "The timestamp of this location "
                                     "in seconds since Epoch",
                                     0,
                                     G_MAXINT64,
                                     0,
                                     G_PARAM_READWRITE |
                                     G_PARAM_CONSTRUCT_ONLY |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_TIMESTAMP, pspec);

        /**
         * GClueLocation:speed
         *
         * The speed in meters per second.
         */
        pspec = g_param_spec_double ("speed",
                                     "Speed",
                                     "Speed in meters per second",
                                     GCLUE_LOCATION_SPEED_UNKNOWN,
                                     G_MAXDOUBLE,
                                     GCLUE_LOCATION_SPEED_UNKNOWN,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_SPEED, pspec);

        /**
         * GClueLocation:heading
         *
         * The positive angle between the direction of movement and the North
         * direction, in clockwise direction. The angle is measured in degrees.
         */
        pspec = g_param_spec_double ("heading",
                                     "Heading",
                                     "The positive Angle between the direction"
                                     " of movement and the North direction, in"
                                     " clockwise direction. The angle is "
                                     "measured in degrees.",
                                     GCLUE_LOCATION_HEADING_UNKNOWN,
                                     G_MAXDOUBLE,
                                     GCLUE_LOCATION_HEADING_UNKNOWN,
                                     G_PARAM_READWRITE |
                                     G_PARAM_STATIC_STRINGS);
        g_object_class_install_property (glocation_class, PROP_HEADING, pspec);
}

static void
gclue_location_init (GClueLocation *location)
{
        location->priv = gclue_location_get_instance_private (location);

        location->priv->altitude = GCLUE_LOCATION_ALTITUDE_UNKNOWN;
        location->priv->accuracy = GCLUE_LOCATION_ACCURACY_UNKNOWN;
        location->priv->speed = GCLUE_LOCATION_SPEED_UNKNOWN;
        location->priv->heading = GCLUE_LOCATION_HEADING_UNKNOWN;
}

static gdouble
get_accuracy_from_hdop (gdouble hdop)
{
        /* FIXME: These are really just rough estimates based on:
         *        http://en.wikipedia.org/wiki/Dilution_of_precision_%28GPS%29#Meaning_of_DOP_Values
         */
        if (hdop <= 1)
                return 0;
        else if (hdop <= 2)
                return 1;
        else if (hdop <= 5)
                return 3;
        else if (hdop <= 10)
                return 50;
        else if (hdop <= 20)
                return 100;
        else
                return 300;
}

static gdouble
parse_coordinate_string (const char *coordinate,
                         const char *direction)
{
        gdouble minutes, degrees, out;
        g_autofree gchar *degrees_str = NULL;
        const char *c, *dot_str, *coordinate_end;
        char *conversion_end;
        gint dot_offset;

        if (!coordinate || !direction || coordinate[0] == '\0' || direction[0] == '\0')
                return INVALID_COORDINATE;

        if (direction[0] != 'N' &&
            direction[0] != 'S' &&
            direction[0] != 'E' &&
            direction[0] != 'W') {
                g_warning ("Unknown direction '%s' for coordinates, ignoring..",
                           direction);
                return INVALID_COORDINATE;
        }

        for (c = coordinate; *c != '\0'; c++) {
                if (!(g_ascii_isdigit (*c) || *c == '.')) {
                        g_warning ("Invalid coordinate string '%s'", coordinate);
                        return INVALID_COORDINATE;
                }
        }
        coordinate_end = c;

        dot_str = g_strstr_len (coordinate, 6, ".");
        dot_offset = dot_str - coordinate;
        if (!dot_str || dot_offset < 3) {
                g_warning ("Invalid coordinate string '%s'", coordinate);
                return INVALID_COORDINATE;
        }

        degrees_str = g_strndup (coordinate, dot_offset - 2);
        degrees = g_ascii_strtod (degrees_str, NULL);

        minutes = g_ascii_strtod (dot_str - 2, &conversion_end);
        if (conversion_end != coordinate_end) {
                g_warning ("Invalid coordinate string '%s'", coordinate);
                return INVALID_COORDINATE;
        }

        /* Include the minutes as part of the degrees */
        out = degrees + (minutes / 60.0);

        if (direction[0] == 'S' || direction[0] == 'W')
                out = 0 - out;

        return out;
}

static gboolean
coordinates_ok (double latitude,
                double longitude)
{
        if (latitude > 90.0 || latitude < -90 ||
            longitude > 180.0 || longitude < -180.0)
                return FALSE;
        else
                return TRUE;
}

static gdouble
parse_altitude_string (const char *altitude,
                       const char *unit)
{
        if (altitude[0] == '\0' || unit[0] == '\0')
                return GCLUE_LOCATION_ALTITUDE_UNKNOWN;

        if (unit[0] != 'M') {
                g_warning ("Unknown unit '%s' for altitude, ignoring..",
                           unit);

                return GCLUE_LOCATION_ALTITUDE_UNKNOWN;
        }

        return g_ascii_strtod (altitude, NULL);
}

/* Return a timestamp derived from the NMEA timestamp and system date as
 * seconds since epoch.
 * If system time cannot be retrieved, return 0.
 * If timestamp parsing fails, return system time.
 * If the parsed time is in the future when compared to the system time,
 * return the parsed time yesterday.
 */
static gint64
parse_nmea_timestamp (const char *nmea_ts)
{
        g_autoptr(GDateTime) now = NULL;
        g_autoptr(GDateTime) midnight = NULL;
        g_autoptr(GDateTime) ts = NULL;
        GTimeSpan timespan;

        now = g_date_time_new_now_utc ();
        if (now == NULL) {
                g_warning ("Failed to get the current UTC time");
                return 0;
        }

        if (!nmea_ts || !*nmea_ts) {  /* Empty timestamp, no warning */
                return g_date_time_to_unix (now);
        }

        timespan = gclue_nmea_timestamp_to_timespan (nmea_ts);
        if (timespan < 0) {
                g_warning ("Failed to parse NMEA timestamp '%s'", nmea_ts);
                return g_date_time_to_unix (now);
        }

        midnight = g_date_time_new_utc (g_date_time_get_year (now),
                                        g_date_time_get_month (now),
                                        g_date_time_get_day_of_month (now),
                                        0,
                                        0,
                                        0);
        ts = g_date_time_add (midnight, timespan);

        if (g_date_time_difference (ts, now) > TIME_DIFF_THRESHOLD) {
                g_debug ("NMEA timestamp '%s' in future. Assuming yesterday's.",
                         nmea_ts);
                g_clear_pointer (&ts, g_date_time_unref);
                ts = g_date_time_add (midnight, timespan - G_TIME_SPAN_DAY);
        }

        return g_date_time_to_unix (ts);
}

/**
 * gclue_location_new:
 * @latitude: a valid latitude
 * @longitude: a valid longitude
 * @accuracy: accuracy of location in meters
 * @description: a description for the location
 *
 * Creates a new #GClueLocation object.
 *
 * Returns: a new #GClueLocation object. Use g_object_unref() when done.
 **/
GClueLocation *
gclue_location_new (gdouble latitude,
                    gdouble longitude,
                    gdouble accuracy,
                    const char *description)
{
        return g_object_new (GCLUE_TYPE_LOCATION,
                             "latitude", latitude,
                             "longitude", longitude,
                             "accuracy", accuracy,
                             "description", description,
                             NULL);
}

/**
 * gclue_location_new_full:
 * @latitude: a valid latitude
 * @longitude: a valid longitude
 * @accuracy: accuracy of location in meters
 * @speed: speed in meters per second
 * @heading: heading in degrees
 * @altitude: altitude of location in meters
 * @timestamp: timestamp in seconds since the Epoch
 * @description: a description for the location
 *
 * Creates a new #GClueLocation object.
 *
 * Returns: a new #GClueLocation object. Use g_object_unref() when done.
 **/
GClueLocation *
gclue_location_new_full (gdouble     latitude,
                         gdouble     longitude,
                         gdouble     accuracy,
                         gdouble     speed,
                         gdouble     heading,
                         gdouble     altitude,
                         guint64     timestamp,
                         const char *description)
{
        return g_object_new (GCLUE_TYPE_LOCATION,
                             "latitude", latitude,
                             "longitude", longitude,
                             "accuracy", accuracy,
                             "speed", speed,
                             "heading", heading,
                             "altitude", altitude,
                             "timestamp", timestamp,
                             "description", description,
                             NULL);
}

static GClueLocation *
gclue_location_create_from_gga (const char *gga)
{
        GClueLocation *location;
        gdouble latitude, longitude, accuracy, altitude;
        gdouble hdop; /* Horizontal Dilution Of Precision */
        guint64 timestamp;
        g_auto(GStrv) parts = NULL;

        parts = g_strsplit (gga, ",", -1);
        if (g_strv_length (parts) < 11) {
                g_warning ("Received short NMEA GGA sentence, discarding.");
                return NULL;
        }

        if (g_ascii_strtoll (parts[6], NULL, 10) == 0) {
                /* No fix, ignore. */
                return NULL;
        }

        /* For syntax of GGA sentences:
         * https://gpsd.gitlab.io/gpsd/NMEA.html#_gga_global_positioning_system_fix_data
         */
        timestamp = parse_nmea_timestamp (parts[1]);
        latitude = parse_coordinate_string (parts[2], parts[3]);
        longitude = parse_coordinate_string (parts[4], parts[5]);
        if (!coordinates_ok (latitude, longitude)) {
                g_warning ("Invalid coordinates (%f, %f) on NMEA GGA sentence.",
                           latitude, longitude);
                return NULL;
        }

        altitude = parse_altitude_string (parts[9], parts[10]);

        hdop = g_ascii_strtod (parts[8], NULL);
        accuracy = get_accuracy_from_hdop (hdop);

        location = g_object_new (GCLUE_TYPE_LOCATION,
                                 "latitude", latitude,
                                 "longitude", longitude,
                                 "accuracy", accuracy,
                                 "timestamp", timestamp,
                                 "description", "GPS GGA",
                                 NULL);
        if (altitude != GCLUE_LOCATION_ALTITUDE_UNKNOWN)
                g_object_set (location, "altitude", altitude, NULL);

        return location;
}

static GClueLocation *
gclue_location_create_from_rmc (const char     *rmc,
                                GClueLocation  *prev_location)
{
        GClueLocation *location;
        g_auto(GStrv) parts = NULL;
        gdouble accuracy;
        gdouble altitude;

        parts = g_strsplit (rmc, ",", -1);
        if (g_strv_length (parts) < 12) {
                g_warning ("Invalid NMEA RMC sentence.");
                return NULL;
        }

        /* RMC sentence is invalid */
        if (g_strcmp0 (parts[2], "A") != 0) {
                return NULL;
        }

        guint64 timestamp = parse_nmea_timestamp (parts[1]);
        gdouble lat = parse_coordinate_string (parts[3], parts[4]);
        gdouble lon = parse_coordinate_string (parts[5], parts[6]);
        if (!coordinates_ok (lat, lon)) {
                g_warning ("Invalid coordinates (%f, %f) on NMEA RMC sentence.",
                           lat, lon);
                return NULL;
        }

        gdouble speed = GCLUE_LOCATION_SPEED_UNKNOWN;
        if (parts[7][0] != '\0')
                speed = g_ascii_strtod (parts[7], NULL) * KNOTS_IN_METERS_PER_SECOND;

        gdouble heading = GCLUE_LOCATION_HEADING_UNKNOWN;
        if (parts[8][0] != '\0')
                heading = g_ascii_strtod (parts[8], NULL);

        /* Some receivers use '0.0,0.0' as invalid speed and heading */
        if (speed == 0.0 && heading == 0.0) {
                speed = GCLUE_LOCATION_SPEED_UNKNOWN;
                heading = GCLUE_LOCATION_HEADING_UNKNOWN;
        }

        accuracy = RMC_DEFAULT_ACCURACY;
        altitude = GCLUE_LOCATION_ALTITUDE_UNKNOWN;
        if (prev_location != NULL) {
                guint64 prev_loc_timestamp;

                prev_loc_timestamp = gclue_location_get_timestamp (prev_location);

                /* Sentence is older then previous location, reject */
                if (timestamp < prev_loc_timestamp)
                        return NULL;

                if (timestamp - prev_loc_timestamp < RMC_TIME_DIFF_THRESHOLD) {
                        accuracy = gclue_location_get_accuracy (prev_location);
                        altitude = gclue_location_get_altitude (prev_location);
                }
        }

        location = g_object_new (GCLUE_TYPE_LOCATION,
                                 "latitude", lat,
                                 "longitude", lon,
                                 "timestamp", timestamp,
                                 "speed", speed,
                                 "heading", heading,
                                 "description", "GPS RMC",
                                 "accuracy", accuracy,
                                 "altitude", altitude,
                                 NULL);

        return location;
}

/**
 * gclue_location_create_from_nmeas:
 * @nmea: A NULL terminated array NMEA sentence strings
 * @prev_location: Previous location provided from the location source
 *
 * Creates a new #GClueLocation object by combining data from multiple NMEA
 * sentences.
 *
 * Returns: a new #GClueLocation object if GGA or RMC sentences are found,
 * a %NULL on all other cases and errors. Unref using #g_object_unref() when
 * done with it.
 **/
GClueLocation *
gclue_location_create_from_nmeas (const char     *nmeas[],
                                  GClueLocation  *prev_location)
{
        GClueLocation *gga_loc = NULL;
        GClueLocation *rmc_loc = NULL;
        const char **iter;

        for (iter = nmeas; *iter != NULL; iter++) {
                if (!gga_loc && gclue_nmea_type_is (*iter, "GGA"))
                        gga_loc = gclue_location_create_from_gga (*iter);
                if (!rmc_loc && gclue_nmea_type_is (*iter, "RMC"))
                        rmc_loc = gclue_location_create_from_rmc
                                (*iter, prev_location);
                if (gga_loc && rmc_loc)
                    break;
        }

        if (gga_loc && rmc_loc) {
                gclue_location_set_speed
                        (gga_loc, gclue_location_get_speed(rmc_loc));
                gclue_location_set_heading
                        (gga_loc, gclue_location_get_heading(rmc_loc));
                g_object_set (gga_loc, "description", "GPS GGA+RMC", NULL);
                g_object_unref (rmc_loc);

                return gga_loc;
        }
        if (gga_loc)
                return gga_loc;
        if (rmc_loc)
                return rmc_loc;

        g_debug ("Valid NMEA GGA or RMC sentence not found");
        return NULL;
}

/**
 * gclue_location_duplicate:
 * @location: the #GClueLocation instance to duplicate.
 *
 * Creates a new copy of @location object (with the same timestamp).
 *
 * Returns: a new #GClueLocation object. Use g_object_unref() when done.
 **/
GClueLocation *
gclue_location_duplicate (GClueLocation *location)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (location), NULL);

        return g_object_new
                (GCLUE_TYPE_LOCATION,
                 "latitude", location->priv->latitude,
                 "longitude", location->priv->longitude,
                 "accuracy", location->priv->accuracy,
                 "altitude", location->priv->altitude,
                 "timestamp", location->priv->timestamp,
                 "speed", location->priv->speed,
                 "heading", location->priv->heading,
                 "description", location->priv->description,
                 NULL);
}

/**
 * gclue_location_duplicate_fresh:
 * @location: the #GClueLocation instance to duplicate.
 *
 * Creates a new copy of @location object with a refreshed timestamp.
 *
 * Returns: a new #GClueLocation object. Use g_object_unref() when done.
 **/
GClueLocation *
gclue_location_duplicate_fresh (GClueLocation *location)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (location), NULL);

        return g_object_new
                (GCLUE_TYPE_LOCATION,
                 "latitude", location->priv->latitude,
                 "longitude", location->priv->longitude,
                 "accuracy", location->priv->accuracy,
                 "altitude", location->priv->altitude,
                 "speed", location->priv->speed,
                 "heading", location->priv->heading,
                 "description", location->priv->description,
                 NULL);
}

const char *
gclue_location_get_description (GClueLocation *loc)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (loc), NULL);

        return loc->priv->description;
}

/**
 * gclue_location_get_latitude:
 * @loc: a #GClueLocation
 *
 * Gets the latitude of location @loc.
 *
 * Returns: The latitude of location @loc.
 **/
gdouble
gclue_location_get_latitude (GClueLocation *loc)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (loc), 0.0);

        return loc->priv->latitude;
}

/**
 * gclue_location_get_longitude:
 * @loc: a #GClueLocation
 *
 * Gets the longitude of location @loc.
 *
 * Returns: The longitude of location @loc.
 **/
gdouble
gclue_location_get_longitude (GClueLocation *loc)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (loc), 0.0);

        return loc->priv->longitude;
}

/**
 * gclue_location_get_altitude:
 * @loc: a #GClueLocation
 *
 * Gets the altitude of location @loc.
 *
 * Returns: The altitude of location @loc.
 **/
gdouble
gclue_location_get_altitude (GClueLocation *loc)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (loc),
                              GCLUE_LOCATION_ALTITUDE_UNKNOWN);

        return loc->priv->altitude;
}

/**
 * gclue_location_get_accuracy:
 * @loc: a #GClueLocation
 *
 * Gets the accuracy (in meters) of location @loc.
 *
 * Returns: The accuracy of location @loc.
 **/
gdouble
gclue_location_get_accuracy (GClueLocation *loc)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (loc),
                              GCLUE_LOCATION_ACCURACY_UNKNOWN);

        return loc->priv->accuracy;
}

/**
 * gclue_location_get_timestamp:
 * @loc: a #GClueLocation
 *
 * Gets the timestamp (in seconds since the Epoch) of location @loc. See
 * #GClueLocation:timestamp.
 *
 * Returns: The timestamp of location @loc.
 **/
guint64
gclue_location_get_timestamp (GClueLocation *loc)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (loc), 0);

        return loc->priv->timestamp;
}

/**
 * gclue_location_get_speed:
 * @location: a #GClueLocation
 *
 * Gets the speed in meters per second.
 *
 * Returns: The speed, or %GCLUE_LOCATION_SPEED_UNKNOWN if speed in unknown.
 **/
gdouble
gclue_location_get_speed (GClueLocation *location)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (location),
                              GCLUE_LOCATION_SPEED_UNKNOWN);

        return location->priv->speed;
}

/**
 * gclue_location_set_speed:
 * @location: a #GClueLocation
 * @speed: speed in meters per second
 *
 * Sets the speed.
 **/
void
gclue_location_set_speed (GClueLocation *location,
                          gdouble        speed)
{
        location->priv->speed = speed;

        g_object_notify (G_OBJECT (location), "speed");
}

/**
 * gclue_location_set_speed_from_prev_location:
 * @location: a #GClueLocation
 * @prev_location: a #GClueLocation
 *
 * Calculates the speed based on provided previous location @prev_location
 * and sets it on @location.
 **/
void
gclue_location_set_speed_from_prev_location (GClueLocation *location,
                                             GClueLocation *prev_location)
{
        gdouble speed;
        guint64 timestamp, prev_timestamp;

        g_return_if_fail (GCLUE_IS_LOCATION (location));
        g_return_if_fail (prev_location == NULL ||
                          GCLUE_IS_LOCATION (prev_location));

        if (prev_location == NULL) {
               speed = GCLUE_LOCATION_SPEED_UNKNOWN;

               goto out;
        }

        timestamp = gclue_location_get_timestamp (location);
        prev_timestamp = gclue_location_get_timestamp (prev_location);

        if (timestamp <= prev_timestamp) {
               speed = GCLUE_LOCATION_SPEED_UNKNOWN;

               goto out;
        }

        speed = gclue_location_get_distance_from (location, prev_location) /
                (timestamp - prev_timestamp);

out:
        location->priv->speed = speed;

        g_object_notify (G_OBJECT (location), "speed");
}

/**
 * gclue_location_get_heading:
 * @location: a #GClueLocation
 *
 * Gets the positive angle between direction of movement and North direction.
 * The angle is measured in degrees.
 *
 * Returns: The heading, or %GCLUE_LOCATION_HEADING_UNKNOWN if heading is
 *          unknown.
 **/
gdouble
gclue_location_get_heading (GClueLocation *location)
{
        g_return_val_if_fail (GCLUE_IS_LOCATION (location),
                              GCLUE_LOCATION_HEADING_UNKNOWN);

        return location->priv->heading;
}

/**
 * gclue_location_set_heading:
 * @location: a #GClueLocation
 * @heading: heading in degrees
 *
 * Sets the heading.
 **/
void
gclue_location_set_heading (GClueLocation *location,
                            gdouble        heading)
{
        location->priv->heading = heading;

        g_object_notify (G_OBJECT (location), "heading");
}

/**
 * gclue_location_set_heading_from_prev_location:
 * @location: a #GClueLocation
 * @prev_location: a #GClueLocation
 *
 * Calculates the heading direction in degrees with respect to North direction
 * based on provided @prev_location and sets it on @location.
 **/
void
gclue_location_set_heading_from_prev_location (GClueLocation *location,
                                               GClueLocation *prev_location)
{
        gdouble dlat, dlon, x, y, angle, lat, lon, prev_lat, prev_lon;

        g_return_if_fail (GCLUE_IS_LOCATION (location));
        g_return_if_fail (prev_location == NULL ||
                          GCLUE_IS_LOCATION (prev_location));

        if (prev_location == NULL) {
               location->priv->heading = GCLUE_LOCATION_HEADING_UNKNOWN;

               return;
        }

        lat = gclue_location_get_latitude (location);
        lon = gclue_location_get_longitude (location);
        prev_lat = gclue_location_get_latitude (prev_location);
        prev_lon = gclue_location_get_longitude (prev_location);

        if (lat == prev_lat && lon == prev_lon) {
               location->priv->heading = GCLUE_LOCATION_HEADING_UNKNOWN;

               return;
        }

        /* atan2(y, x) is a function which takes in coordinate values of
         * a 2D point and returns the angle of line from origin to that
         * coordinate makes with the positive X-axis, in the range (-PI,+PI].
         *
         * We want to match the inclusive end (+PI) of atan2's range to our
         * desired zero direction (north), so we take the vector pointing
         * south as our X-axis and the vector pointing east as our Y-axis
         * and measure our movement vector coordinates */
        dlon = (lon - prev_lon);
        dlat = (lat - prev_lat);
        /* in that basis.
         *
         * Latitudes decrease towards the south and longitudes
         * increase towards the east, so: */
        x = -dlat;
        y = dlon;
        angle = atan2(y, x) * 180.0 / M_PI;
        /* We now have the angle of our movement vector measured from the
         * vector pointing south in degrees, with the positive direction
         * being counterclockwise.
         *
         * We want to convert it to heading: Degrees from from north with
         * positive values only and positive direction being clockwise.
         *
         * If angle is positive (the movement is more towards east than
         * west), we need to subtract it from the heading of our reference
         * vector (south == 180 deg). If the angle is negative, we need to
         * add its negative to the heading of the reference vector. Both
         * cases result in */
        location->priv->heading = 180.0 - angle;

        g_object_notify (G_OBJECT (location), "heading");
}

/**
 * gclue_location_get_distance_from:
 * @loca: a #GClueLocation
 * @locb: a #GClueLocation
 *
 * Calculates the distance in meters, along the curvature of the Earth,
 * between 2 locations. Note that altitude changes are not
 * taken into account.
 *
 * Returns: a distance in meters.
 **/
double
gclue_location_get_distance_from (GClueLocation *loca,
                                  GClueLocation *locb)
{
        gdouble dlat, dlon, lat1, lat2;
        gdouble a, c;

        g_return_val_if_fail (GCLUE_IS_LOCATION (loca), 0.0);
        g_return_val_if_fail (GCLUE_IS_LOCATION (locb), 0.0);

        /* Algorithm from:
         * http://www.movable-type.co.uk/scripts/latlong.html */

        dlat = (locb->priv->latitude - loca->priv->latitude) * M_PI / 180.0;
        dlon = (locb->priv->longitude - loca->priv->longitude) * M_PI / 180.0;
        lat1 = loca->priv->latitude * M_PI / 180.0;
        lat2 = locb->priv->latitude * M_PI / 180.0;

        a = sin (dlat / 2) * sin (dlat / 2) +
            sin (dlon / 2) * sin (dlon / 2) * cos (lat1) * cos (lat2);
        c = 2 * atan2 (sqrt (a), sqrt (1-a));
        return 1000.0 * EARTH_RADIUS_KM * c;
}
