'Returns all Locations whose geometry intersects the given feature.
Relies on Shapely.
@ToDo: provide an option to use PostGIS/Spatialite'

@staticmethod
def get_features_by_feature(feature):

    from shapely.wkt import loads as wkt_loads
    shape = wkt_loads(feature.wkt)
    return current.gis.get_features_by_shape(shape)
