Obliczanie odległości między punktem a wirtualną linią dwóch lat / lngs

Proszę odnieść się do przykładu i odpowiedniego obrazu.

Ja zamierzamy uzyskać następujące: dostarczenie dwóch miejscach (łac / ng), które są pokazane poniżej jako A i B . Na tej podstawie zostanie narysowana linia wirtualna, a następnie obliczona zostanie odległość między tą linią a C (w dowolnym pomiarze).

Osiągnąłem to obecnie w Google Maps API v3, ale chciałbym również móc to zrobić za kulisami w moim wybranym języku. Wszelkie wskazówki / pomysły będą mile widziane!

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

Najkrótsza długość to wymagany dystans, chyba że się mylę?

Może robię to zbyt skomplikowane, ale to, czego chcesz, to odległość od punktu do linii. Jest to odległość od punktu wzdłuż AB, który łączy AB z C linią prostopadłą do AB. Ten wektor prostopadły do ​​AB podano przez

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

(Użyłem nawiasów kwadratowych do zdefiniowania wektora lub tablicy dwuelementowej). Odległość między C [xp, yp] a punktem A wynosi

u=[x1-xp, y1-xp]

Odległość między linią a C jest tylko rzutem u na v. Jeśli założymy mod (v) = 1 (po prostu normalizujemy), to

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

Jedyną komplikacją jest to, że prawdopodobnie chcesz się upewnić, że twoje współrzędne nie są parami lat / log WGS84, ale są rzutowane (lub używasz współrzędnych geodezyjnych). Możesz do tego użyć OGR lub Proj4 .

Będąc trochę niechętnym całej tej matematyce, podchodziłbym do tego z innej perspektywy. Uczyniłbym to „rzeczywistą” linią, a nie linią wirtualną, a następnie używałbym istniejących narzędzi.

Jeśli A i B mają wspólny atrybut, możesz połączyć je poprzez narysowanie linii (Kosmo GIS ma narzędzie, które utworzy linie z punktów, i sądzę, że jest do tego również wtyczka QGIS). Po uzyskaniu linii funkcja „near” na warstwie punktowej „C” da ci odległość do linii. Niech oprogramowanie zajmie się matematyką!

Jeśli korzystasz z java na Androidzie, jest to tylko jedna linia z funkcją biblioteki

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

distanceToLine:

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

Oblicza odległość na kuli między punktem p a segmentem linii zaczynającym się kończyć.

Parametry: p - punkt, który należy zmierzyć

start - początek odcinka linii

end - koniec segmentu linii

Zwraca: odległość w metrach (przy założeniu kulistej ziemi)

Po prostu dodaj bibliotekę do swojego

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