Calcolo della distanza tra un punto e una linea virtuale di due lat / lng

Fare riferimento all'esempio e all'immagine corrispondente.

Desidero di raggiungere i seguenti: fornire due riprese (lat / lng), che sono di seguito riportate come A e B . Da questo, verrebbe tracciata una linea virtuale e quindi la distanza tra questa linea e C verrebbe calcolata (in qualsiasi misura).

L'ho raggiunto attualmente nell'API v3 di Google Maps, ma vorrei poterlo eseguire anche dietro le quinte nella mia lingua prescelta. Eventuali suggerimenti / idee sarebbero molto apprezzati!

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 lunghezza più breve è la distanza richiesta, a meno che non mi sbagli?

Forse lo sto rendendo troppo complicato, ma quello che vuoi è la distanza da un punto a una linea. Questa è la distanza da un punto lungo AB che collega AB con C con una linea ortogonale ad AB. Questo vettore perpendicolare ad AB è dato da

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

(Ho usato le parentesi quadre per definire un vettore o una matrice a due elementi). La distanza tra C [xp, yp] e il punto A è

u=[x1-xp, y1-xp]

La distanza tra la linea e C è solo la proiezione di u su v. Se assumiamo mod (v) = 1 (semplicemente normalizzalo), allora

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

L'unica complicazione è che probabilmente vuoi assicurarti che le tue coordinate non siano coppie lat / log WGS84, ma proiettate (o usa coordinate geodetiche). È possibile utilizzare OGR o Proj4 per questo.

Essendo un po 'contrario anche a tutta questa matematica, ci arriverei da una prospettiva diversa. Vorrei renderlo una linea "reale", piuttosto che una linea virtuale, e quindi utilizzare gli strumenti esistenti.

Se A e B condividono un attributo, potresti collegarli tracciando una linea (Kosmo GIS ha uno strumento che creerà linee da punti, e credo che ci sia anche un plugin QGIS per questo). Una volta che hai le linee, una funzione 'vicino' sul livello del punto 'C' ti darà la distanza dalla linea. Lascia che il software gestisca la matematica per te!

Se stavi usando Java su Android, è solo una riga con la funzione di libreria

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

distanceToLine:

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

Calcola la distanza sulla sfera tra il punto p e il segmento di linea dall'inizio alla fine.

Parametri: p - il punto da misurare

inizio - l'inizio del segmento di linea

end - la fine del segmento di linea

Restituisce: la distanza in metri (assumendo la terra sferica)

Aggiungi libreria al tuo

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