Comment mesurer la distance et créer un cadre de sélection basé sur deux points de latitude et de longitude en Java?

StackOverflow https://stackoverflow.com/questions/120283

Question

Je souhaite trouver la distance entre deux points différents. Je sais que cela peut être accompli avec la distance du grand cercle. http://www.meridianworlddata.com/Distance-calculation.asp

Une fois terminé, avec un point et une distance, j'aimerais trouver le point qui se trouve au nord et à l'est afin de créer une boîte autour du point.

Était-ce utile?

La solution

Nous avons eu quelques succès en utilisant OpenMap pour tracer de nombreuses données de position. Il existe une classe LatLonPoint qui possède des bases fonctionnalité, y compris la distance.

Autres conseils

Voici une implémentation Java de la formule Haversine . J'utilise ceci dans un projet pour calculer la distance en miles entre lat / longs.

public static double distFrom(double lat1, double lng1, double lat2, double lng2) {
    double earthRadius = 3958.75; // miles (or 6371.0 kilometers)
    double dLat = Math.toRadians(lat2-lat1);
    double dLng = Math.toRadians(lng2-lng1);
    double sindLat = Math.sin(dLat / 2);
    double sindLng = Math.sin(dLng / 2);
    double a = Math.pow(sindLat, 2) + Math.pow(sindLng, 2)
            * Math.cos(Math.toRadians(lat1)) * Math.cos(Math.toRadians(lat2));
    double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));
    double dist = earthRadius * c;

    return dist;
    }

Vous pouvez également utiliser SimpleLatLng . Apache 2.0 est sous licence et utilisé dans un système de production que je connais: le mien.

Petite histoire:

Je recherchais une géo-bibliothèque simple et ne parvenais pas à en trouver une pour répondre à mes besoins. Et qui veut écrire, tester et déboguer ces petits outils géographiques encore et encore dans chaque application? Il doit y avoir un meilleur moyen!

So SimpleLatLng est né comme moyen de stocker des données de latitude-longitude, d’effectuer des calculs de distance et de créer des limites profilées.

Je sais que je suis deux ans trop tard pour aider l’affiche originale, mais mon but est d’aider les personnes comme moi qui trouvent cette question dans une recherche. J'aimerais que certaines personnes l'utilisent et contribuent aux tests et à la vision de ce petit utilitaire léger.

Pour une distance plus précise (0,5 mm), vous pouvez également utiliser l'approximation de Vincenty:

/**
 * Calculates geodetic distance between two points specified by latitude/longitude using Vincenty inverse formula
 * for ellipsoids
 * 
 * @param lat1
 *            first point latitude in decimal degrees
 * @param lon1
 *            first point longitude in decimal degrees
 * @param lat2
 *            second point latitude in decimal degrees
 * @param lon2
 *            second point longitude in decimal degrees
 * @returns distance in meters between points with 5.10<sup>-4</sup> precision
 * @see <a href="http://www.movable-type.co.uk/scripts/latlong-vincenty.html">Originally posted here</a>
 */
public static double distVincenty(double lat1, double lon1, double lat2, double lon2) {
    double a = 6378137, b = 6356752.314245, f = 1 / 298.257223563; // WGS-84 ellipsoid params
    double L = Math.toRadians(lon2 - lon1);
    double U1 = Math.atan((1 - f) * Math.tan(Math.toRadians(lat1)));
    double U2 = Math.atan((1 - f) * Math.tan(Math.toRadians(lat2)));
    double sinU1 = Math.sin(U1), cosU1 = Math.cos(U1);
    double sinU2 = Math.sin(U2), cosU2 = Math.cos(U2);

    double sinLambda, cosLambda, sinSigma, cosSigma, sigma, sinAlpha, cosSqAlpha, cos2SigmaM;
    double lambda = L, lambdaP, iterLimit = 100;
    do {
        sinLambda = Math.sin(lambda);
        cosLambda = Math.cos(lambda);
        sinSigma = Math.sqrt((cosU2 * sinLambda) * (cosU2 * sinLambda)
                + (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda) * (cosU1 * sinU2 - sinU1 * cosU2 * cosLambda));
        if (sinSigma == 0)
            return 0; // co-incident points
        cosSigma = sinU1 * sinU2 + cosU1 * cosU2 * cosLambda;
        sigma = Math.atan2(sinSigma, cosSigma);
        sinAlpha = cosU1 * cosU2 * sinLambda / sinSigma;
        cosSqAlpha = 1 - sinAlpha * sinAlpha;
        cos2SigmaM = cosSigma - 2 * sinU1 * sinU2 / cosSqAlpha;
        if (Double.isNaN(cos2SigmaM))
            cos2SigmaM = 0; // equatorial line: cosSqAlpha=0 (§6)
        double C = f / 16 * cosSqAlpha * (4 + f * (4 - 3 * cosSqAlpha));
        lambdaP = lambda;
        lambda = L + (1 - C) * f * sinAlpha
                * (sigma + C * sinSigma * (cos2SigmaM + C * cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM)));
    } while (Math.abs(lambda - lambdaP) > 1e-12 && --iterLimit > 0);

    if (iterLimit == 0)
        return Double.NaN; // formula failed to converge

    double uSq = cosSqAlpha * (a * a - b * b) / (b * b);
    double A = 1 + uSq / 16384 * (4096 + uSq * (-768 + uSq * (320 - 175 * uSq)));
    double B = uSq / 1024 * (256 + uSq * (-128 + uSq * (74 - 47 * uSq)));
    double deltaSigma = B
            * sinSigma
            * (cos2SigmaM + B
                    / 4
                    * (cosSigma * (-1 + 2 * cos2SigmaM * cos2SigmaM) - B / 6 * cos2SigmaM
                            * (-3 + 4 * sinSigma * sinSigma) * (-3 + 4 * cos2SigmaM * cos2SigmaM)));
    double dist = b * A * (sigma - deltaSigma);

    return dist;
}

Ce code a été librement adapté de http: //www.movable- type.co.uk/scripts/latlong-vincenty.html

Formule Haversine Distance corrigée ....

public static double HaverSineDistance(double lat1, double lng1, double lat2, double lng2) 
{
    // mHager 08-12-2012
    // http://en.wikipedia.org/wiki/Haversine_formula
    // Implementation

    // convert to radians
    lat1 = Math.toRadians(lat1);
    lng1 = Math.toRadians(lng1);
    lat2 = Math.toRadians(lat2);
    lng2 = Math.toRadians(lng2);

    double dlon = lng2 - lng1;
    double dlat = lat2 - lat1;

    double a = Math.pow((Math.sin(dlat/2)),2) + Math.cos(lat1) * Math.cos(lat2) * Math.pow(Math.sin(dlon/2),2);

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

    return EARTH_RADIUS * c;
}   

http://www.movable-type.co.uk/scripts/ latlong.html

public static Double distanceBetweenTwoLocationsInKm(Double latitudeOne, Double longitudeOne, Double latitudeTwo, Double longitudeTwo) {
        if (latitudeOne == null || latitudeTwo == null || longitudeOne == null || longitudeTwo == null) {
            return null;
        }

        Double earthRadius = 6371.0;
        Double diffBetweenLatitudeRadians = Math.toRadians(latitudeTwo - latitudeOne);
        Double diffBetweenLongitudeRadians = Math.toRadians(longitudeTwo - longitudeOne);
        Double latitudeOneInRadians = Math.toRadians(latitudeOne);
        Double latitudeTwoInRadians = Math.toRadians(latitudeTwo);
        Double a = Math.sin(diffBetweenLatitudeRadians / 2) * Math.sin(diffBetweenLatitudeRadians / 2) + Math.cos(latitudeOneInRadians) * Math.cos(latitudeTwoInRadians) * Math.sin(diffBetweenLongitudeRadians / 2)
                * Math.sin(diffBetweenLongitudeRadians / 2);
        Double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
        return (earthRadius * c);
    }

Vous pouvez utiliser le bibliothèque de géodésie Java pour le GPS , il utilise les les formules de Vincenty qui tiennent compte de la surface de la Terre courbure.

La mise en œuvre se déroule comme suit:

import org.gavaghan.geodesy.*;
...
GeodeticCalculator geoCalc = new GeodeticCalculator();
Ellipsoid reference = Ellipsoid.WGS84;
GlobalPosition pointA = new GlobalPosition(latitude, longitude, 0.0);
GlobalPosition userPos = new GlobalPosition(userLat, userLon, 0.0);
double distance = geoCalc.calculateGeodeticCurve(reference, userPos, pointA).getEllipsoidalDistance();

La distance résultante est en mètres.

Je sais qu’il existe de nombreuses réponses, mais en effectuant des recherches sur ce sujet, j’ai constaté que la plupart des réponses utilisées ici utilisaient la formule Haversine, mais la formule de Vincenty est en réalité plus précise. Un article a adapté le calcul à partir d'une version Javascript, mais il est très difficile à manier. J'ai trouvé une version supérieure car:

  1. Il possède également une licence ouverte.
  2. Il utilise les principes de la POO.
  3. Il est plus facile de choisir l'ellipsoïde que vous souhaitez utiliser.
  4. Plusieurs méthodes sont disponibles pour permettre différents calculs à l'avenir.
  5. C'est bien documenté.

VincentyDistanceCalculator

Cette méthode vous aiderait à trouver la distance entre un emplacement géographique en km.

private double getDist(double lat1, double lon1, double lat2, double lon2)
{
    int R = 6373; // radius of the earth in kilometres
    double lat1rad = Math.toRadians(lat1);
    double lat2rad = Math.toRadians(lat2);
    double deltaLat = Math.toRadians(lat2-lat1);
    double deltaLon = Math.toRadians(lon2-lon1);

    double a = Math.sin(deltaLat/2) * Math.sin(deltaLat/2) +
            Math.cos(lat1rad) * Math.cos(lat2rad) *
            Math.sin(deltaLon/2) * Math.sin(deltaLon/2);
    double c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a));

    double d = R * c;
    return d;
}

J'utilise généralement MATLAB avec la Mapping Toolbox , puis j'utilise le code dans Java. Utilisation de la JA MATLAB Builder. Cela me simplifie grandement la vie. Etant donné que la plupart des écoles disposent d’un accès gratuit pour les étudiants, vous pouvez l’essayer (ou obtenir la version d’essai pour en finir avec votre travail).

Pour Android, il existe une approche simple.

 public static float getDistanceInMeter(LatLng start, LatLng end) { 
    float[] results = new float[1];
    Location.distanceBetween(start.latitude, start.longitude, end.latitude, end.longitude, results);
    return results[0];

}

;

https://developer.android.com/reference/android/location/Location#distanceBetween (lat1, lng1, lat2, lng2, sortie [])

Licencié sous: CC-BY-SA avec attribution
Non affilié à StackOverflow
scroll top