Converti una query SQL complessa in SQLAlchemy

Il tuo HAVING è gestito correttamente, ma gli stai passando l'espressione sbagliata. Sembra che tu stia usando Python 2, dal momento che il confronto relazionale tra una stringa e un intero

'distance' < 25

non genera un'eccezione, ma restituisce False invece. In altre parole la tua query è uguale a

locations = db.session.query(...).having(False).all()

il che spiega perché ottieni zero risultati:tutte le righe sono esplicitamente filtrate dalla clausola HAVING, come si vede nella versione stampata:

...
HAVING false = 1  -- remove all rows

Una soluzione è usare un costrutto adatto, come column() , per produrre l'espressione:

locations = db.session.query(...).having(column('distance') < 25).all()

Non dovresti racchiudere l'espressione complessa dell'elemento dell'elenco di selezione in un select() , che rappresenta un'istruzione SELECT. O etichetta il text() così com'è:

text('( 6371 * acos( cos( radians("53.6209798282177") ) * '
     'cos( radians( lat ) ) * cos( radians( lng ) - radians("13.96948162900808") ) + '
     'sin( radians("53.6209798282177") ) * sin( radians( lat ) ) ) ) '
     'AS distance')

oppure costruisci l'espressione usando il modello:

(6371 *
 func.acos(func.cos(func.radians(53.6209798282177)) *
           func.cos(func.radians(Location.lat)) *
           func.cos(func.radians(Location.lng) - func.radians(13.96948162900808)) +
           func.sin(func.radians(53.6209798282177)) *
           func.sin(func.radians(Location.lat)))).label('distance')

Potresti migliorare la leggibilità della costruzione della tua query creando una funzione per la distanza del cerchio grande , e con un po' di lavoro potresti implementare un metodo ibrido su Location :

import math

def gc_distance(lat1, lng1, lat2, lng2, math=math):
    ang = math.acos(math.cos(math.radians(lat1)) *
                    math.cos(math.radians(lat2)) *
                    math.cos(math.radians(lng2) -
                             math.radians(lng1)) +
                    math.sin(math.radians(lat1)) *
                    math.sin(math.radians(lat2)))

    return 6371 * ang

class Location(db.Model):
    ...
    @hybrid_method
    def distance(self, lat, lng):
        return gc_distance(lat, lng, self.lat, self.lng)

    @distance.expression
    def distance(cls, lat, lng):
        return gc_distance(lat, lng, cls.lat, cls.lng, math=func)

locations = db.session.query(
        Location,
        Location.distance(53.6209798282177,
                          13.96948162900808).label('distance')).\
    having(column('distance') < 25).\
    order_by('distance').\
    all()

Nota che il modo in cui usi HAVING per eliminare le righe non di gruppo non è portabile. Ad esempio in Postgresql la presenza della clausola HAVING trasforma una query in una query raggruppata, anche senza una clausola GROUP BY. Potresti invece usare una sottoquery:

stmt = db.session.query(
        Location,
        Location.distance(53.6209798282177,
                          13.96948162900808).label('distance')).\
    subquery()

location_alias = db.aliased(Location, stmt)

locations = db.session.query(location_alias).\
    filter(stmt.c.distance < 25).\
    order_by(stmt.c.distance).\
    all()