Me hizo buscar un poco en google para "Esférica Índice Espacial". Hay un montón de métodos posibles, utilizando triangular de la descomposición de la esfera, o mosaicos de voronoi. Uno de los métodos que se ve fácilmente realizable sin embargo, es de considerar sus datos en 3d, como en el "3D del Cuadro Delimitador de la sección" aquí:
http://lin-ear-th-inking.blogspot.co.uk/2007/09/geodetic-data-in-postgis-spherical.html
Entonces usted necesita un 3d índice espacial de algún tipo, entonces usted puede encontrar rápidamente todos los puntos dentro de su 1km. Este sería un 1km 3D radio de búsqueda, de modo ligeramente diferente a un 1km de radio a lo largo de la superficie de la tierra, pero para los pequeños de búsqueda, los radios, sería idéntico (hacer la matemáticas para trabajar la corrección).
Si desea absoluta precisión, utilice esto como un primer paso y, a continuación, calcular las distancias a través de un gran círculo de eliminar aquellos que más lejos (distancia a lo largo de una esfera es siempre mayor que la distancia a través de una esfera).