diff --git a/urbantz/__init__.py b/urbantz/__init__.py
new file mode 100644
index 0000000..d0c478e
--- /dev/null
+++ b/urbantz/__init__.py
@@ -0,0 +1,6 @@
+#!/usr/bin/env python3
+"""
+UrbanTZ API helper
+"""
+
+from urbantz.delivery import Delivery  # noqa: F401
diff --git a/urbantz/utils.py b/urbantz/utils.py
new file mode 100644
index 0000000..ea754ca
--- /dev/null
+++ b/urbantz/utils.py
@@ -0,0 +1,108 @@
+from math import radians, cos, sin, asin, sqrt, trunc, floor, ceil
+
+EARTH_RADIUS_KM = 6371
+
+
+class Coordinates(object):
+    """
+    A helper class for GPS coordinates.
+    """
+
+    def __init__(self, lat=0, lng=0):
+        """
+        Coordinates from decimal degrees
+        """
+        self.lat = lat
+        self.lng = lng
+
+    def __repr__(self):
+        return '{}(lat={}, lng={})'.format(
+            self.__class__.__name__, self.lat, self.lng)
+
+    def __str__(self):
+        return ', '.join(tuple(self))
+
+    def __hash__(self):
+        return hash(tuple(self))
+
+    def __len__(self):
+        return 2
+
+    def __length_hint__(self):
+        return len(self)
+
+    def __iter__(self):
+        return (self.lat, self.lng)
+
+    def __dict__(self):
+        return {'lat': self.lat, 'lng': self.lng}
+
+    def __eq__(self, other):
+        return hasattr(other, 'lat') and hasattr(other, 'lng') and \
+            tuple(self) == tuple(other)
+
+    def __add__(self, other):
+        if not hasattr(other, 'lat') or hasattr(other, 'lng'):
+            return NotImplemented
+        return self.__class__(
+            lat=self.lat + other.lat,
+            lng=self.lng + other.lng,
+        )
+
+    def __sub__(self, other):
+        if not hasattr(other, 'lat') or hasattr(other, 'lng'):
+            return NotImplemented
+        return self.__class__(
+            lat=self.lat - other.lat,
+            lng=self.lng - other.lng,
+        )
+
+    def __pos__(self):
+        return self
+
+    def __neg__(self):
+        return self.__class__(lat=-self.lat, lng=-self.lng)
+
+    def __abs__(self):
+        return self.__class__(lat=abs(self.lat), lng=abs(self.lng))
+
+    def __round__(self, ndigits=None):
+        return self.__class__(
+            lat=round(self.lat, ndigits),
+            lng=round(self.lng, ndigits),
+        )
+
+    def __trunc__(self):
+        return self.__class__(lat=trunc(self.lat), lng=trunc(self.lng))
+
+    def __floor__(self):
+        return self.__class__(lat=floor(self.lat), lng=floor(self.lng))
+
+    def __ceil__(self):
+        return self.__class__(lat=ceil(self.lat), lng=ceil(self.lng))
+
+    def to_radians(self):
+        return tuple(map(radians, self))
+
+    def distance(self, other):
+        """
+        Compute Haversine distance between two coordinates.
+        """
+        if not hasattr(other, 'to_radians'):
+            raise NotImplementedError(
+                'Distance requires a to_radians() method on both coordinates')
+        lat1, lon1 = self.to_radians()
+        lat2, lon2 = other.to_radians()
+        dlon, dlat = lon2 - lon1, lat2 - lat1
+        a = sin(dlat / 2) ** 2 + cos(lat1) * cos(lat2) * sin(dlon / 2) ** 2
+        return EARTH_RADIUS_KM * 2000 * asin(sqrt(a))
+
+    def toJSON(self):
+        """
+        Convert to UrbanTZ JSON geometry
+        """
+        return [self.lng, self.lat]
+
+    @staticmethod
+    def fromJSON(geometry):
+        return Coordinates(lng=geometry[0], lat=geometry[1])