Du berechnest die Koordinaten des Vierecks(!) in dem die Orte liegen dürfen und gibst aus der Datenbank alle Orte zurück, welche diese Bedingung erfüllen.
Im zweiten Schritt berechnest Du die Entfernung dieser Orte zum Ausgangspunkt und gibst die Orte zurück, die innerhalb des Umkreises liegen.
Das muss man ausprobieren, ob eine näherung mit einer Bounding-Box schneller ist als gleich direkt zu berechnen.