Блог

Я немного изучил эту тему, но есть много мнений, которые точно не дают четкого изображения. Моя проблема заключается в следующем: я разрабатываю приложение на основе GPS для Android, в котором я хочу знать расстояние между моим текущим местоположением, указанным Androids LocationManager, и другим местоположением в режиме реального времени. Я попробовал Haversine formula, формулу Закона косинусов, затем я обнаружил, что Android SDK дает мне простую функцию Location.distanceTo (Местоположение) — я не уверен, на каком методе работает эта функция.

Итак, вопрос в том, какой из них будет полезен для меня в ситуациях, когда реальное расстояние между этими местоположениями в большинстве случаев не будет превышать примерно 100-200 метров? Может быть, мне следует проверить формулы Винсенти? Это действительно так медленно? Может кто-нибудь, пожалуйста, объясните мне, что я должен выбрать?

Не используйте distanceTo. Используйте метод distanceBetween, поскольку, похоже, у вас уже есть координаты, и это все, что вам нужно для этого метода: Location.distanceBetween() Javadoc

Заглядывая в Android-источник для определения distanceTo (местоположения), вы можете видеть, что результат основан на «обратной формуле» геодезии:

Который основан на использовании «Обратной формулы» (раздел 4)

Кроме того, два метода distanceTo и distanceBetween используют один и тот же базовый метод. У них просто есть альтернативные формы ввода / вывода.

Для полноты картины полный источник этого вычисления приведен ниже, но я рекомендую вам самостоятельно проверить класс Location в android.location. (PS Я не проверял правильность вычисления Android. Это было бы хорошим упражнением!)

     private static void computeDistanceAndBearing(double lat1, double lon1,
    double lat2, double lon2, float[] results) {
    // Based on http://www.ngs.noaa.gov/PUBS_LIB/inverse.pdf
    // using the "Inverse Formula" (section 4)

    int MAXITERS = 20;
    // Convert lat/long to radians
    lat1 *= Math.PI / 180.0;
    lat2 *= Math.PI / 180.0;
    lon1 *= Math.PI / 180.0;
    lon2 *= Math.PI / 180.0;

    double a = 6378137.0; // WGS84 major axis
    double b = 6356752.3142; // WGS84 semi-major axis
    double f = (a - b) / a;
    double aSqMinusBSqOverBSq = (a * a - b * b) / (b * b);

    double L = lon2 - lon1;
    double A = 0.0;
    double U1 = Math.atan((1.0 - f) * Math.tan(lat1));
    double U2 = Math.atan((1.0 - f) * Math.tan(lat2));

    double cosU1 = Math.cos(U1);
    double cosU2 = Math.cos(U2);
    double sinU1 = Math.sin(U1);
    double sinU2 = Math.sin(U2);
    double cosU1cosU2 = cosU1 * cosU2;
    double sinU1sinU2 = sinU1 * sinU2;

    double sigma = 0.0;
    double deltaSigma = 0.0;
    double cosSqAlpha = 0.0;
    double cos2SM = 0.0;
    double cosSigma = 0.0;
    double sinSigma = 0.0;
    double cosLambda = 0.0;
    double sinLambda = 0.0;

    double lambda = L; // initial guess
    for (int iter = 0; iter < MAXITERS; iter  ) {
        double lambdaOrig = lambda;
        cosLambda = Math.cos(lambda);
        sinLambda = Math.sin(lambda);
        double t1 = cosU2 * sinLambda;
        double t2 = cosU1 * sinU2 - sinU1 * cosU2 * cosLambda;
        double sinSqSigma = t1 * t1   t2 * t2; // (14)
        sinSigma = Math.sqrt(sinSqSigma);
        cosSigma = sinU1sinU2   cosU1cosU2 * cosLambda; // (15)
        sigma = Math.atan2(sinSigma, cosSigma); // (16)
        double sinAlpha = (sinSigma == 0) ? 0.0 :
            cosU1cosU2 * sinLambda / sinSigma; // (17)
        cosSqAlpha = 1.0 - sinAlpha * sinAlpha;
        cos2SM = (cosSqAlpha == 0) ? 0.0 :
            cosSigma - 2.0 * sinU1sinU2 / cosSqAlpha; // (18)

        double uSquared = cosSqAlpha * aSqMinusBSqOverBSq; // defn
        A = 1   (uSquared / 16384.0) * // (3)
            (4096.0   uSquared *
             (-768   uSquared * (320.0 - 175.0 * uSquared)));
        double B = (uSquared / 1024.0) * // (4)
            (256.0   uSquared *
             (-128.0   uSquared * (74.0 - 47.0 * uSquared)));
        double C = (f / 16.0) *
            cosSqAlpha *
            (4.0   f * (4.0 - 3.0 * cosSqAlpha)); // (10)
        double cos2SMSq = cos2SM * cos2SM;
        deltaSigma = B * sinSigma * // (6)
            (cos2SM   (B / 4.0) *
             (cosSigma * (-1.0   2.0 * cos2SMSq) -
              (B / 6.0) * cos2SM *
              (-3.0   4.0 * sinSigma * sinSigma) *
              (-3.0   4.0 * cos2SMSq)));

        lambda = L  
            (1.0 - C) * f * sinAlpha *
            (sigma   C * sinSigma *
             (cos2SM   C * cosSigma *
              (-1.0   2.0 * cos2SM * cos2SM))); // (11)

        double delta = (lambda - lambdaOrig) / lambda;
        if (Math.abs(delta) < 1.0e-12) {
            break;
        }
    }

    float distance = (float) (b * A * (sigma - deltaSigma));
    results[0] = distance;
    if (results.length > 1) {
        float initialBearing = (float) Math.atan2(cosU2 * sinLambda,
            cosU1 * sinU2 - sinU1 * cosU2 * cosLambda);
        initialBearing *= 180.0 / Math.PI;
        results[1] = initialBearing;
        if (results.length > 2) {
            float finalBearing = (float) Math.atan2(cosU1 * sinLambda,
                -sinU1 * cosU2   cosU1 * sinU2 * cosLambda);
            finalBearing *= 180.0 / Math.PI;
            results[2] = finalBearing;
        }
    }
}