Come ottenere rapidamente 13 righe di record vicino a un punto (lon, lat) postgis

È possibile applicare la funzione distanza utilizzando l'operatore <-> direttamente nel ORDER BY clausola. In questo modo il pianificatore utilizzerà l'indice gist:

EXPLAIN (ANALYSE,COSTS OFF)
SELECT * FROM data
ORDER BY geog <-> ST_Point(6.5, 48.7) 
LIMIT 13;

                            QUERY PLAN
----------------------------------------------------------------------
Limit (actual time=15.019..15.213 rows=13 loops=1)
  ->  Index Scan using idx_data_point on data (actual time=15.017..15.210 rows=13 loops=1)
        Order By: (geog <-> '0101000020E61000000000000000001A409A99999999594840'::geography)
Planning Time: 0.372 ms
Execution Time: 15.309 ms

Demo:db<>fiddle