Offset a LatLon object by a heading (in degrees) and distance (in km)
        to return a new LatLon object

def offset(self, heading_initial, distance, ellipse = 'WGS84'):
        '''
        Offset a LatLon object by a heading (in degrees) and distance (in km)
        to return a new LatLon object
        '''
        lat1, lon1 = self.lat.decimal_degree, self.lon.decimal_degree
        g = pyproj.Geod(ellps = ellipse)
        distance = distance * 1000 # Convert km to meters
        lon2, lat2, back_bearing = g.fwd(lon1, lat1, heading_initial, distance, radians = False)
        return LatLon(Latitude(lat2), Longitude(lon2))