transformer la latitude de longitude en mètres

Il y a des bibliothèques sur jstott.me.uk pour PHP , Java et JavaScript lequel faisiez-le, par exemple

var lld1 = new LatLng(40.718119, -73.995667); // New York
document.write("New York Lat/Long: " + lld1.toString() + "<br />");
var lld2 = new LatLng(51.499981, -0.125313);  // London
document.write("London Lat/Long: " + lld2.toString() + "<br />");
var d = lld1.distance(lld2);
document.write("Surface Distance between New York and London: " + d + "km");

Pour commencer, je pense que vous avez votre latitude et longitude inversée. Mesures de longitude x et mesures de latitude Y.

La latitude est facile à transformer en une distance nord-sud. Nous savons que 360 ​​degrés est un cercle complet autour de la terre à travers les poteaux et Cette distance est 40008000 mètres . Tant que vous n'avez pas besoin de prendre en compte les erreurs dues à la Terre étant PAS parfaitement sphérique, la formule est Deltatalatitude * 40008000/360 code>. P>

La partie délicate est la conversion Longitude to x, comme vous le soupçonnez. Comme cela dépend de la latitude dont vous avez besoin pour choisir la latitude que vous allez utiliser - vous pouvez choisir la latitude de votre origine, la latitude de votre destination ou un point arbitraire entre les deux. La circonférence de l'équateur (latitude 0) est de 40075160 mètres. La circonférence d'un cercle à une latitude donnée sera proportionnelle au cosinus, la formule sera donc Deltalongitude * 40075160 * COS (latitude) / 360 Code>. P>

Modifier : strong> Votre commentaire indique que vous avez eu des problèmes avec la formule de longitude; Vous avez peut-être utilisé des diplômes au lieu de radians dans l'appel à cos code>, c'est une erreur commune des recrues. Pour vous assurer qu'il n'y a pas d'ambiguïté, voici le code de travail de Python. P>

def asRadians(degrees):
    return degrees * pi / 180

def getXYpos(relativeNullPoint, p):
    """ Calculates X and Y distances in meters.
    """
    deltaLatitude = p.latitude - relativeNullPoint.latitude
    deltaLongitude = p.longitude - relativeNullPoint.longitude
    latitudeCircumference = 40075160 * cos(asRadians(relativeNullPoint.latitude))
    resultX = deltaLongitude * latitudeCircumference / 360
    resultY = deltaLatitude * 40008000 / 360
    return resultX, resultY