Calcul de la distance entre un point et une ligne virtuelle de deux lat / lngs

Veuillez vous référer à l'exemple et à l'image correspondante.

Je voudrais atteindre les objectifs suivants: fournir deux emplacements (lat / lng), qui sont présentés ci - dessous en A et B . À partir de cela, une ligne virtuelle serait tracée, puis la distance entre cette ligne et C serait calculée (dans n'importe quelle mesure).

J'ai atteint cet objectif actuellement dans Google Maps API v3, mais je voudrais également pouvoir effectuer cela en coulisses dans la langue de mon choix. Tous les conseils / idées seraient grandement appréciés!

def get_perp( X1, Y1, X2, Y2, X3, Y3):
    """************************************************************************************************ 
    Purpose - X1,Y1,X2,Y2 = Two points representing the ends of the line segment
              X3,Y3 = The offset point 
    'Returns - X4,Y4 = Returns the Point on the line perpendicular to the offset or None if no such
                        point exists
    '************************************************************************************************ """
    XX = X2 - X1 
    YY = Y2 - Y1 
    ShortestLength = ((XX * (X3 - X1)) + (YY * (Y3 - Y1))) / ((XX * XX) + (YY * YY)) 
    X4 = X1 + XX * ShortestLength 
    Y4 = Y1 + YY * ShortestLength
    if X4 < X2 and X4 > X1 and Y4 < Y2 and Y4 > Y1:
        return X4,Y4
    return None

La longueur la plus courte est la distance dont vous avez besoin, sauf erreur de ma part?

Peut-être que je rends les choses trop compliquées, mais ce que vous voulez, c'est la distance d'un point à une ligne. C'est la distance d'un point le long de AB qui relie AB à C avec une ligne orthogonale à AB. Ce vecteur perpendiculaire à AB est donné par

v=[x2-x1, -(y2-y1)] # Point A is [x1,y1] Point B is [x2,y2]

(J'ai utilisé les crochets pour définir un vecteur ou un tableau à deux éléments). La distance entre C [xp, yp] et le point A est

u=[x1-xp, y1-xp]

La distance entre la ligne et C n'est que la projection de u sur v. Si nous supposons que mod (v) = 1 (il suffit de le normaliser), alors

distance = u*v = abs( (x2-x1)*(y1-yp) - (x1-xp)*(y2-y1) )

La seule complication est que vous voulez probablement vous assurer que vos coordonnées ne sont pas des paires lat / log WGS84, mais projetées (ou utilisez des coordonnées géodésiques). Vous pouvez utiliser OGR ou Proj4 pour cela.

Étant un peu opposé à toutes ces mathématiques aussi, j'y viendrais sous un angle différent. Je voudrais en faire une ligne «réelle» plutôt qu'une ligne virtuelle, puis utiliser les outils existants.

Si A et B partagent un attribut, vous pouvez les connecter en traçant une ligne (Kosmo GIS dispose d'un outil qui créera des lignes à partir de points, et je pense qu'il existe également un plugin QGIS pour cela). Une fois que vous avez les lignes, une fonction «proche» sur la couche de points «C» vous donnera la distance à la ligne. Laissez le logiciel gérer les mathématiques pour vous!

Si vous utilisiez Java sur Android, sa seule ligne avec la fonction de bibliothèque

import static com.google.maps.android.PolyUtil.distanceToLine;

distanceToLine:

public static double distanceToLine(LatLng p, LatLng start,LatLng end)

Calcule la distance sur la sphère entre le point p et le segment de ligne du début à la fin.

Paramètres: p - le point à mesurer

start - le début du segment de ligne

end - la fin du segment de ligne

Renvoie: la distance en mètres (en supposant la terre sphérique)

Ajoutez simplement une bibliothèque à votre

dependencies {
    compile 'com.google.maps.android:android-maps-utils:0.5+'
}