Вычисление расстояния между точкой и виртуальной линией двух ширин / широт

Пожалуйста, обратитесь к примеру и соответствующему изображению.

Я хотел бы добиться следующего: обеспечить два места (лат / LNG), которые показаны ниже как A и B . Из этого будет проведена виртуальная линия, а затем будет вычислено расстояние между этой линией и C (при любом измерении).

Я достиг этого в настоящее время в Google Maps API v3, но хотел бы также иметь возможность выполнять это за кулисами на моем языке по своему выбору. Любые советы / идеи будут с благодарностью!

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

Самая короткая длина - это расстояние, которое вам нужно, если я не ошибаюсь?

Может быть, я делаю это слишком сложным, но вам нужно расстояние от точки до прямой. Это расстояние от точки вдоль AB, которая связывает AB с C линией, ортогональной к AB. Этот вектор, перпендикулярный AB, определяется как

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

(Я использовал квадратные скобки, чтобы определить вектор или массив из двух элементов). Расстояние между C [xp, yp] и точкой A

u=[x1-xp, y1-xp]

Расстояние между линией и C - это просто проекция u на v. Если предположить, что mod (v) = 1 (просто нормализовать его), то

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

Единственное осложнение заключается в том, что вы, вероятно, хотите убедиться, что ваши координаты не являются парами WGS84 lat / log, а спроецированы (или используются геодезические координаты). Вы можете использовать OGR или Proj4 для этого.

Будучи немного отвращенным ко всей этой математике, я бы подошел к этому с другой стороны. Я бы сделал это «реальной» строкой, а не виртуальной строкой, а затем использовал бы существующие инструменты.

Если A и B имеют общий атрибут, вы можете соединить их, нарисовав линию (в Kosmo GIS есть инструмент, который будет создавать линии из точек, и я считаю, что для этого есть также плагин QGIS). Если у вас есть линии, функция «ближний» на точечном слое «С» даст вам расстояние до линии. Пусть программное обеспечение справится с математикой для вас!

Если вы использовали Java на Android, это только одна строка с библиотечной функцией

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

distanceToLine:

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

Вычисляет расстояние на сфере между точкой p и отрезком от начала до конца.

Параметры: p - точка для измерения

начало - начало отрезка

конец - конец отрезка

Возвращает: расстояние в метрах (при условии сферической земли)

Просто добавьте библиотеку в свой

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