Я ищу функцию, которая возвращает точку (широта, долгота) между двумя точками (где я также указываю их широту и долготу), и эта точка основана на проценте расстояния.
Итак, я указываю Lat1, Lon1, Lat2, Lon2 и% в функции, и она возвращает точку, например, удаленную на 20% от первой точки до второй.
Предполагая, что координата является десятичным числом. Вы можете использовать это уравнение.
function midpoint(lat1, long1, lat2, long2, per) {
return [lat1 + (lat2 - lat1) * per, long1 + (long2 - long1) * per];
}
Верните новую желаемую координату [широта, долгота] на основе процента (например, per = 0,2 для 20%).
Вот ссылка, которая очень поможет (проверьте внизу)
http://www.movable-type.co.uk/scripts/latlong.html
Промежуточная точка
Также может быть вычислена промежуточная точка в любой части пути большого круга между двумя точками.
Формула:
a = sin((1−f)⋅δ) / sin δ b = sin(f⋅δ) / sin δ x = a ⋅ cos φ1 ⋅ cos λ1 + b ⋅ cos φ2 ⋅ cos λ2 y = a ⋅ cos φ1 ⋅ sin λ1 + b ⋅ cos φ2 ⋅ sin λ2 z = a ⋅ sin φ1 + b ⋅ sin φ2 φi = atan2(z, √x² + y²) λi = atan2(y, x)
где f - дробь вдоль маршрута большого круга (f = 0 - точка 1, f = 1 - точка 2), δ < / em> - угловое расстояние d / R между двумя точками.
Оба ответа могут быть полезны для определенных пользователей, но я хотел бы указать на некоторые проблемы.
Решение lguiel подходит для игр или небольших расстояний, но не работает для геометрических вычислений на Земле как глобусе.
Ответ спирографа правильный, но, поскольку он носит теоретический характер, его может быть слишком сложно запрограммировать.
Я перевел второй ответ на практический язык программирования, чтобы вы могли использовать его в своих проектах.
Полный код | Запускать в Интернете
// Original calculation from https://www.movable-type.co.uk/scripts/latlong.html
LatLng calculateIntermediatePoint(LatLng point1, LatLng point2, double perc) {
//const φ1 = this.lat.toRadians(), λ1 = this.lon.toRadians();
//const φ2 = point.lat.toRadians(), λ2 = point.lon.toRadians();
double lat1 = degreesToRadians(point1.latitude);
double lng1 = degreesToRadians(point1.longitude);
double lat2 = degreesToRadians(point2.latitude);
double lng2 = degreesToRadians(point2.longitude);
//const Δφ = φ2 - φ1;
//const Δλ = λ2 - λ1;
double deltaLat = lat2 - lat1;
double deltaLng = lng2 - lng1;
//const a = Math.sin(Δφ/2) * Math.sin(Δφ/2) + Math.cos(φ1) * Math.cos(φ2) * Math.sin(Δλ/2) * Math.sin(Δλ/2);
//const δ = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
double calcA = sin(deltaLat / 2) * sin(deltaLat / 2) +
cos(lat1) * cos(lat2) * sin(deltaLng / 2) * sin(deltaLng / 2);
double calcB = 2 * atan2(sqrt(calcA), sqrt(1 - calcA));
//const A = Math.sin((1-fraction)*δ) / Math.sin(δ);
//const B = Math.sin(fraction*δ) / Math.sin(δ);
double A = sin((1 - perc) * calcB) / sin(calcB);
double B = sin(perc * calcB) / sin(calcB);
//const x = A * Math.cos(φ1) * Math.cos(λ1) + B * Math.cos(φ2) * Math.cos(λ2);
//const y = A * Math.cos(φ1) * Math.sin(λ1) + B * Math.cos(φ2) * Math.sin(λ2);
//const z = A * Math.sin(φ1) + B * Math.sin(φ2);
double x = A * cos(lat1) * cos(lng1) + B * cos(lat2) * cos(lng2);
double y = A * cos(lat1) * sin(lng1) + B * cos(lat2) * sin(lng2);
double z = A * sin(lat1) + B * sin(lat2);
//const φ3 = Math.atan2(z, Math.sqrt(x*x + y*y));
//const λ3 = Math.atan2(y, x);
double lat3 = atan2(z, sqrt(x * x + y * y));
double lng3 = atan2(y, x);
//const lat = φ3.toDegrees();
//const lon = λ3.toDegrees();
return LatLng(radiansToDegrees(lat3), radiansToDegrees(lng3));
}
