'Returns a generator of locations whose shape intersects the given LatLon.
Relies on Shapely.
@todo: provide an option to use PostGIS/Spatialite'

@staticmethod
def get_features_by_latlon(lat, lon):

    from shapely.geometry import point
    return current.gis.get_features_by_shape(point.Point(lon, lat))
