docs: 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))