No entendí muy bien tu código, pero creo haber entendido la pregunta.
A veces puede parecer "complicado" trabajar con coordenadas, pues la tierra no es plana, por lo que no se debe aplicar geometría euclideana, sino la geometría esférica.
Análisis matemático del problema
Por definición, para saber si un punto se encuentra dentro del área circular, basta con saber si la distancia del punto al centro del círculo es menor o igual al radio.
Ahora, el problema es entonces saber cuál es la distancia del punto al centro del círculo.
Aplicar la fórmula de distancia entre puntos que siempre hemos conocido no sirve, pues es sólo aplicable para un espacio euclideano, y evidentemente no estamos trabajando en un espacio euclideano.
Y de hecho, la ruta más corta para llegar a un determinado punto, no es una recta como sucede en el espacio euclideano.
En la Figura 1 se puede observar eso.
Ahora entonces, supongamos que de la Figura 1 queremos calcular la distancia c entre los puntos w y v, para ello utilizaremos la fórmula de Haversine que nos permite calcular el ángulo central entre dos puntos.
¿Pero para qué queremos el ángulo, si buscamos la distancia? Porque si recordamos un poco, la longitud de arco (lo que estamos buscando) divido entre el radio del círculo (en este caso esfera, pero lo reducimos a dos dimensiones en donde se encuentran v, w y el centro de la esfera pues no necesitamos nada más) es igual al ángulo central.
Donde big-theta es el ángulo central, S es la longitud de arco y r es el radio
Entonces por lo tanto, conociendo el radio y el ángulo podemos sacar la distancia entre los puntos.
El radio es el radio de la tierra porque pues, estamos trabajando con coordenadas GPS.
Entonces, en conclusión debemos hacer los siguientes sencillos pasos para calcular la distancia esférica entre dos puntos en una esfera:
- Obtener la función haversine del ángulo central entre los dos puntos con la fórmula de haversine. Revisa el link de Wikipedia para saber los detalles de la fórmula.
- Aplicar la función inversa a lo obtenido en el paso anterior para obtener el ángulo central entre los dos puntos.
- Calcular la distancia entre los puntos con la fórmula de longitud de arco, pues ya conocemos el ángulo central y el radio.
Implementación en JS
Conociendo lo anterior es sencillo programarlo, sin embargo debemos de tomar en cuenta que las coordenadas que proporcionaste son coordenadas decimales y se requieren coordenadas en radianes, entonces basta simplemente con convertir de grados decimales a radianes. Para ello está la siguiente función
const deg2Rad = degrees => degrees * Math.PI / 180;
Luego la fórmula 4 de haversine queda así
const hav = theta => (1 - Math.cos(theta)) / 2;
Y con esas dos funciones ya podemos programar otra función que realice el paso 1. Dicha función es la siguiente.
const havPoints = (point1, point2) => {
// latitudes in radians
const lat1 = deg2Rad(point1.latitude);
const lat2 = deg2Rad(point2.latitude);
// longitudes in radians
const long1 = deg2Rad(point1.longitude);
const long2 = deg2Rad(point2.longitude);
// haversine formula
return hav(lat2 - lat1) + Math.cos(lat1) * Math.cos(lat2) * hav(long2 - long1);
};
Por último, para hacer el paso 2 y 3, utilizamos la siguiente función
const calculateDistance = (point1, point2, radius=6_371) => {
// hav(central angle)
const hav_central_angle = havPoints(point1, point2);
// now that we have the haversine of the central angle its time to get the central angle by applying the inverse
// central angle = archav(hav(central angle)) = 2arcsin(hav(central angle))
const central_angle = 2 * Math.asin(Math.sqrt(hav_central_angle));
// now just apply the well known formula d = r * theta
return radius * central_angle;
};
Espero que eso te haya servido de ayuda, disculpa los comentarios en inglés en el código, usualmente siempre programo en inglés, pero a final de cuentas creo que el cómo sacar la distancia entre dos puntos quedó claro.
Teniendo estas funciones, es fácil programar la función que filtré las coordenadas que están dentro del círculo, simplemente hay que revisar si la distancia del centro del círculo al punto es menor o igual que el radio del círculo, si es así, entonces, obviamente está dentro.
Snippet completo
Te dejo un snippet para que lo veas en acción
/**
* Converts the given degrees to radians
* @param degrees degrees :)
* @return the equivalent in radians
*/
const deg2Rad = degrees => degrees * Math.PI / 180;
/**
* Calculates de haversine of the central angle between two points
* @param point1 an object with properties latitude and longitude in DEGREES
* @param point2 an object with properties latitude and longitude in DEGREES
* @return the middle angle (in radians) between the two points
*/
const havPoints = (point1, point2) => {
// latitudes in radians
const lat1 = deg2Rad(point1.latitude);
const lat2 = deg2Rad(point2.latitude);
// longitudes in radians
const long1 = deg2Rad(point1.longitude);
const long2 = deg2Rad(point2.longitude);
// haversine formula
return hav(lat2 - lat1) + Math.cos(lat1) * Math.cos(lat2) * hav(long2 - long1);
};
/**
* Calculates the haversine of the given angle
* @param theta the ange in RADIANS
* @return the haversine in radians
*/
const hav = theta => (1 - Math.cos(theta)) / 2;
/**
* Uses the haversine function to calculate the SPHERICAL distance between two points
* @param point1 an object with properties latitude and longitude in DEGREES
* @param point2 an object with properties latitude and longitude in DEGREES
* @param radius the radius of the sphere in km. Default is the radius of the earth -> 6,371km
* @return the distance in km between the two points
*/
const calculateDistance = (point1, point2, radius=6_371) => {
// hav(central angle)
const hav_central_angle = havPoints(point1, point2);
// now that we have the haversine of the central angle its time to get the central angle by applying the inverse
// central angle = archav(hav(central angle)) = 2arcsin(hav(central angle))
const central_angle = 2 * Math.asin(Math.sqrt(hav_central_angle));
// now just apply the well known formula d = r * theta
return radius * central_angle;
};
const filterCoordsByDistance = (center, points, minAcceptableDistance) => points.filter(
point => calculateDistance(center, point) <= minAcceptableDistance
);
const transformData = () => {
let coords = {
"coordOne": [4.3005119999999994, -75.69080319999999],
"coordTwo": [7.8005119999999994,-5.69080319999999],
"coordThree": [4.8005119999999994,-75.66080319999999],
}
return Object.values(coords).map(coords => {
return {
latitude: coords[0],
longitude: coords[1]
}
});
}
const center = {
latitude: 4.8005119999999994,
longitude: -75.69080319999999
};
const data = transformData();
console.log(`Distancias de los puntos al centro, en orden: ${data.map(point => calculateDistance(center, point)).sort()}`);
const circleRadius = 100;
const pointsInRange = filterCoordsByDistance(center, transformData(), circleRadius);
console.log(`Puntos que están en un radio de ${circleRadius}km del centro: ${JSON.stringify(pointsInRange)}`);
Si lo ejecutaste, habrás visto que sólo un punto no está dentro del círculo con radio de 100km, si quieres comprobarlo ve a Google Maps o a alguno otro similar y calcula la distancia entre los tres puntos con el centro del círculo y verás que en efecto, las distancias son correctas y sólo deben existir dos puntos en el círculo de radio 100km.
Últimas notas
- Nuevamente perdón por el código en inglés, espero la explicación que fue dada sea suficiente.
- El código no está 100% testeado y listo para producción obviamente.
- Mencionas que se debe filtrar de acuerdo al área circular, supongo que te refieres al radio, en cualquier caso puedes obtener el radio fácilmente.
userCoords
) representa el centro de la circunferencia. El radio de la misma estaría dado en km, por el valor dearea
. Y lo que necesitas determinar es si los puntos representados por las coordenadas encoordOne
,coordTwo
ycoordThree
se encuentran dentro o fuera de la circunferencia definida anteriormente. ¿Es esto correcto?