Использование формулы Haversine в Javascript

Я пытаюсь использовать формулу расстояния Haversine (как найдено здесь:http://www.movable-type.co.uk/scripts/latlong.html) но я не могу заставить его работать, см. следующий код

    function test() { 
    var lat2 = 42.741; 
    var lon2 = -71.3161; 
    var lat1 = 42.806911; 
    var lon1 = -71.290611; 

    var R = 6371; // km 
    //has a problem with the .toRad() method below.
    var dLat = (lat2-lat1).toRad();  
    var dLon = (lon2-lon1).toRad();  
    var a = Math.sin(dLat/2) * Math.sin(dLat/2) + 
                    Math.cos(lat1.toRad()) * Math.cos(lat2.toRad()) * 
                    Math.sin(dLon/2) * Math.sin(dLon/2);  
    var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
    var d = R * c; 

    alert(d); 
}

и ошибка:

Uncaught TypeError: Object -0.06591099999999983 has no method 'toRad'

Я понимаю, что он должен сделать следующее:

Number.prototype.toRad = function() {
return this * Math.PI / 180;
}

но когда я помещаю это ниже функции, он все равно возвращается с тем же сообщением об ошибке. Как заставить его использовать метод helper? Или есть ли альтернативный способ получить этот код, чтобы заставить его работать? Спасибо!

25
haversine javascript

8 ответов

этот код работает:

Number.prototype.toRad = function() {
   return this * Math.PI / 180;
}

var lat2 = 42.741; 
var lon2 = -71.3161; 
var lat1 = 42.806911; 
var lon1 = -71.290611; 

var R = 6371; // km 
//has a problem with the .toRad() method below.
var x1 = lat2-lat1;
var dLat = x1.toRad();  
var x2 = lon2-lon1;
var dLon = x2.toRad();  
var a = Math.sin(dLat/2) * Math.sin(dLat/2) + 
                Math.cos(lat1.toRad()) * Math.cos(lat2.toRad()) * 
                Math.sin(dLon/2) * Math.sin(dLon/2);  
var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
var d = R * c; 

alert(d);

обратите внимание, как я определил x1 и x2. Играть с ним в: https://tinker.io/3f794

38

вот рефакторинг функции на основе 3 другие ответы!

обратите внимание, что аргументами coords являются [долгота, широта].

function haversineDistance(coords1, coords2, isMiles) {
  function toRad(x) {
    return x * Math.PI / 180;
  }

  var lon1 = coords1[0];
  var lat1 = coords1[1];

  var lon2 = coords2[0];
  var lat2 = coords2[1];

  var R = 6371; // km

  var x1 = lat2 - lat1;
  var dLat = toRad(x1);
  var x2 = lon2 - lon1;
  var dLon = toRad(x2)
  var a = Math.sin(dLat / 2) * Math.sin(dLat / 2) +
    Math.cos(toRad(lat1)) * Math.cos(toRad(lat2)) *
    Math.sin(dLon / 2) * Math.sin(dLon / 2);
  var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
  var d = R * c;

  if(isMiles) d /= 1.60934;

  return d;
}
20

почему бы не попробовать прямое решение? Вместо расширения прототипа номера просто определите toRad как регулярную функцию:

function toRad(x) {
   return x * Math.PI / 180;
}

и затем вызов toRad везде:

var dLat = toRad(lat2-lat1);

расширение прототипа номера не всегда работает так, как ожидалось. Например, вызов 123.toRad() не работает. Я думаю, что если вы var x1 = lat2 - lat1; x1.toRad(); работает лучше, чем делать (lat2-lat1).toRad()

13

ЕС6 на JavaScript/NodeJS рефакторингу версия:

/**
 * Calculates the haversine distance between point A, and B.
 * @param {number[]} latlngA [lat, lng] point A
 * @param {number[]} latlngB [lat, lng] point B
 * @param {boolean} isMiles If we are using miles, else km.
 */
const haversineDistance = (latlngA, latlngB, isMiles) => {
  const toRad = x => (x * Math.PI) / 180;
  const R = 6371; // km

  const dLat = toRad(latlngB[1] - latlngA[1]);
  const dLatSin = Math.sin(dLat / 2);
  const dLon = toRad(latlngB[0] - latlngA[0]);
  const dLonSin = Math.sin(dLon / 2);

  const a = (dLatSin * dLatSin) +
            (Math.cos(toRad(latlngA[1])) * Math.cos(toRad(latlngB[1])) * dLonSin * dLonSin);
  const c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
  let distance = R * c;

  if (isMiles) distance /= 1.60934;

  return distance;
}
8

когда я помещаю это ниже функции

вам нужно только поставить его выше точки, где вы называете test(). Где test сама функция объявлена не имеет значения.

1

перед вызовом этих расширений в функции необходимо расширить прототип номера.

поэтому просто убедитесь 

Number.prototype.toRad = function() {
  return this * Math.PI / 180;
}

называется до ваша функция вызывается.

1

другой вариант для уменьшения избыточности, а также совместим с объектами Google LatLng:

  function haversine_distance(coords1, coords2) {

     function toRad(x) {
         return x * Math.PI / 180;
    }

  var dLat = toRad(coords2.latitude - coords1.latitude);
  var dLon = toRad(coords2.longitude - coords1.longitude)

  var a = Math.sin(dLat / 2) * Math.sin(dLat / 2) +
          Math.cos(toRad(coords1.latitude)) * 
          Math.cos(toRad(coords2.latitude)) *
          Math.sin(dLon / 2) * Math.sin(dLon / 2);

  return 12742 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
}
1

это Java-имплеметация решения talkol выше. Его или ее решение очень хорошо для нас. Я не пытаюсь ответить на вопрос, так как исходный вопрос был для javascript. Я просто делюсь нашей Java-реализацией данного решения javascript, если другие найдут его использование.

// this was a pojo class we used internally...
public class GisPostalCode {

    private String country;
    private String postalCode;
    private double latitude;
    private double longitude;

    // getters/setters, etc.
}


public static double distanceBetweenCoordinatesInMiles2(GisPostalCode c1, GisPostalCode c2) {

    double lat2 = c2.getLatitude();
    double lon2 = c2.getLongitude();
    double lat1 = c1.getLatitude();
    double lon1 = c1.getLongitude();

    double R = 6371; // km
    double x1 = lat2 - lat1;
    double dLat = x1 * Math.PI / 180;
    double x2 = lon2 - lon1;
    double dLon = x2 * Math.PI / 180;

    double a = Math.sin(dLat/2) * Math.sin(dLat/2) +
        Math.cos(lat1*Math.PI/180) * Math.cos(lat2*Math.PI/180) *
        Math.sin(dLon/2) * Math.sin(dLon/2);

    double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
    double d = R * c;

    // convert to miles
    return d / 1.60934;
}
0