Calcular la distancia entre un punto y una línea virtual de dos lat / lngs

Consulte el ejemplo y la imagen correspondiente.

Me gustaría lograr lo siguiente: proporcionar dos ubicaciones (lat / lng), que se muestran a continuación como A y B . A partir de esto, se trazaría una línea virtual y luego se calcularía la distancia entre esta línea y C (en cualquier medida).

Lo he logrado actualmente en la API v3 de Google Maps, pero me gustaría poder realizar esto detrás de escena en el idioma que elija. Cualquier sugerencia / idea sería muy apreciada!

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 longitud más corta es la distancia que necesita, a menos que me equivoque?

Tal vez lo estoy haciendo demasiado complicado, pero lo que quieres es la distancia de un punto a una línea. Esa es la distancia desde un punto a lo largo de AB que une AB con C con una línea ortogonal a AB. Este vector perpendicular a AB viene dado por

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

(He utilizado los corchetes para definir un vector o una matriz de dos elementos). La distancia entre C [xp, yp] y el punto A es

u=[x1-xp, y1-xp]

La distancia entre la línea y C es solo la proyección de u sobre v. Si suponemos que mod (v) = 1 (solo lo normalizamos), entonces

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

La única complicación es que probablemente desee asegurarse de que sus coordenadas no sean pares lat / log WGS84, sino proyectadas (o utilice coordenadas geodésicas). Puede usar OGR o Proj4 para esto.

Siendo un poco reacio a todas estas matemáticas también, lo abordaría desde un ángulo diferente. Lo convertiría en una línea 'real', en lugar de una línea virtual, y luego usaría las herramientas existentes.

Si A y B comparten un atributo, puede conectarlos dibujando una línea (Kosmo GIS tiene una herramienta que creará líneas a partir de puntos, y creo que también hay un complemento QGIS para esto). Una vez que tenga las líneas, una función 'cerca' en la capa de puntos 'C' le dará la distancia a la línea. ¡Deje que el software maneje las matemáticas por usted!

Si estabas usando Java en Android, es solo una línea con la función de biblioteca

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

distanceToLine:

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

Calcula la distancia en la esfera entre el punto p y el segmento de línea de principio a fin.

Parámetros: p - el punto a medir

inicio: el comienzo del segmento de línea

end - el final del segmento de línea

Retornos: la distancia en metros (suponiendo tierra esférica)

Simplemente agregue la biblioteca a su

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