Рассчитать расстояние по координатам в 1С #501045

Есть лог GPS-логера с координатами. Нужно рассчитать сколько проехала машина. Поиски привели к "формуле гаверсинуса" и "расстояние большого круга" Вот только как это на 1С перевести...

о боже, сложить длины отрезков - погрешность менее 5% будет...

у тебя же координаты на каждую секунду должны быть - так там вообще за 0% погреность принять можно

Не ругайся а пальцем ткни ;)

Школа была давно, институт не окончен ...

А что мешает? В третьей ссылке все написано. sin, cos и arctan в 1с есть.

atan2 нет в 1С

что такое atan2? и чем просто atan не устраивает?

Если скуль, то вывести расчет в хранимку - и быстрее и ограничений в формулах нет.

Типа этого? и по ADO вызывать?

Достаточно

В код невчитывался, но суть в передаче в скуль файла (или таблицы) и получении назад результата. Можно даже с 1срр заморочиться.

Восьмерка

все указанное в так же можно получить из объекта Math джабыскрипт

тогда через жабаскрипт

Примеры можно?

В поиске забанили?

пример коорднат

lat="47.877987" lon="35.071194" lat="47.877289" lon="35.072243" lat="47.876556" lon="35.072796"

рисуешь прямоугольник в голове, потом считаешь диагональ как в 5м классе, и умножаешь на 111 км/градус

А atan2 в какой из 3-х ссылок упоминается? Что-то не увидел. Вся эта "сферичность" нужна только для немаленьких расстояний, когда сказывается сферичность Земли. Для небольших расстояний можно Землю считать плоской+сделать поправку (1 градус - 111 км).

(|lat1-lat2|^2+|lon1-lon2|^2) ^  1/2

+20. Или формулу в ряд разложить. Наверное, будет 2-й член (разница в квадрате).

Получается 1,26 + 0,92 = 2,18

Расстояние между двумя первыми точками 110,23 метра

это по какой формуле? или это в реале?

в первой ссылке классный пример в Экселе, пробовал вводить разные данные 1я формула отлично считает, а там нет никаких ATAN2

По Формуле гаверсинусов и пример в Excel из первой ссылки

Есть еще подъем в гору, или спуск с горы. Или например ты проехал 50 метров прямо, повернул на 90 грудосов и проехал еще 50 метров, т.е. реально ты проехал - 100 метров, а по координатам - 70 метров (длина гипотенузы < длины суммы катетов)

у него тока координаты, датчика уровня над водой нет я подозреваю)

Из testcase? Пытаюсь перевести

У меня в екселе есть эта функция

Есть и высота. Но расчетов с высотой не видел

Здесь на Delphi написано. может поможет

Пля

:))) я думал только я такой умный

Да у меня тоже есть )), там в примере 3 формулы расчета 1я через ACOS,остальные ATAN2 и ATAN. Пробовал вводить разные координаты-результаты совершенно одинаковые, при условии, что первая функция проще. Так что можно без джавы, стандарной 8 функцией.

+

всё интересное пропустил :)

Может я че то не понимаю, но Результат = sqrt((x2-x1)*(x2-x1)) + (у2-у1)(у2-у1))

+ В общем думаю самый простой вариант это взять длину дуги градуса широты и долготы для местности, погрешность будет, но не большая. Тогда длина дуги между двумя точками будет такая: длина пути = ((lon_1*d_Lat)^2+(lon_1*d_Lon)^2)^(1/2) где: lat_1 - длина одного градуса широты для местности lon_1 - длина одного градуса долготы для местности d_Lat=lat1-lon2 - разница широт d_Lon=lon1-lon2 - разница долгот.

ну и какие там проблемы вычисления ? lat="47.877987" lon="35.071194" lat="47.877289" lon="35.072243" lat="47.876556" lon="35.072796" расстояние == sqrt((47.877987 - 47.877289)^2 + (35.071194 - 35.072243)^2) + sqrt((47.877289 - 47.876556)^2 + (35.072243 - 35.072796)^2)

На радиус Земли в этой местности нужно еще умножить.

Не сходится

блин, так это ж долгота и широта :)

через яндекс.api можно...Он по своим картам может и маршрут проложить и расстояние вернуть

всё просто : вводим вектор-функцию и берём интеграл =)

Ругается

Получится елси перевести из градусов в радианы и воспользоваться

с чем не сходится???

- проверь линейкой

с

+ 1) Длина дуги окружности = радиус * угол в радианах. 2) На малых маштабах сферичностью можно принебречь

IronDemon, посмотр решение от систеха, там расстояние считается.

Можно и так, только брать интеграл не аналитически, а численно.

Уряяяяя. Функция ABS(Число)    Если Число<0 Тогда        Возврат -1*Число;    Иначе        Возврат Число;    КонецЕсли; КонецФункции Процедура КнопкаВыполнитьНажатие(Кнопка)    Пи=3.14;    РадиусЗ = 6372795;    Ширина1=47.877987;    Долгота1=35.071194;    Ширина2=47.877289;    Долгота2=35.072243;        Расстояние = РадиусЗ*ATAN(Sqrt(Pow(COS(Пи*Ширина2/180)*SIN(ABS(Пи*Долгота2/180-Пи*Долгота1/180)),2)+Pow(COS(Пи*Ширина1/180)*SIN(Пи*Ширина2/180)-SIN(Пи*Ширина1/180)*COS(Пи*Ширина2/180)*COS(ABS(Пи*Долгота2/180-Пи*Долгота1/180)),2))/(SIN(Пи*Ширина1/180)*SIN(Пи*Ширина2/180)+COS(Пи*Ширина1/180)*COS(Пи*Ширина2/180)*COS(ABS(Пи*Долгота2/180-Пи*Долгота1/180))));    Сообщить(Расстояние);     КонецПроцедуры

Результат какой?

rulerDistance {Number} rulerDistance(point1, point2) Возвращает расстояние между двумя заданными точками (в метрах), если двигаться вдоль прямой соединяющей точки на карте. Используется для расчета расстояний в элементе управления "Линейка". Параметры:Имя    Тип    Описание point1    YMaps.GeoPoint    Первая точка. point2    YMaps.GeoPoint    Вторая точка. Возвращает:Тип    Описание Number    Расстояние "по линейке" между двумя заданными точками.

Расстояние = 110,20373939206339422

в закладочку =)

интересно, а северном полюсе взлетит функция?

Вот из переделал. вроде работает Процедура РасстояниеПоКоординатам(StartLat,StartLong,EndLat,EndLong,Distance,Bearing)  // Передаваемые широта/долгота в градусах и сотых долях  //StartLat: double; // Начальная широта  //StartLong: double; // Начальная долгота  //EndLat: double; // Конечная широта  //EndLong: double; // Конечная долгота  //// Переменные, используемые для вычисления смещения и расстояния  //fPhimean: Double; // Средняя широта  //fdLambda: Double; // Разница между двумя значениями долготы  //fdPhi: Double; // Разница между двумя значениями широты  //fAlpha: Double; // Смещение  //fRho: Double; // Меридианский радиус кривизны  //fNu: Double; // Поперечный радиус кривизны  //fR: Double; // Радиус сферы Земли  //fz: Double; // Угловое расстояние от центра сфероида  //fTemp: Double; // Временная переменная, использующаяся в вычислениях  //Distance: Double; // Вычисленное расстояния в метрах  //Bearing: Double; // Вычисленное от и до смещение  // Константы, используемые для вычисления смещения и расстояния  D2R = 0.017453; // Константа для преобразования градусов в радианы  R2D = 57.295781; // Константа для преобразования радиан в градусы  a = 6378137.0; // Основные полуоси  b = 6356752.314245; // Неосновные полуоси  e2 = 0.006739496742337; // Квадрат эксцентричности эллипсоида  f = 0.003352810664747; // Выравнивание эллипсоида  fdLambda = (StartLong - EndLong) * D2R;  fdPhi = (StartLat - EndLat) * D2R;  fPhimean = ((StartLat + EndLat) / 2.0) * D2R;  // Вычисляем меридианные и поперечные радиусы кривизны средней широты  fTemp = 1 - e2 * (Pow(Sin(fPhimean), 2));  fRho = (a * (1 - e2)) / Pow(fTemp, 1.5);  fNu = a / (Sqrt(1 - e2 * (Sin(fPhimean) * Sin(fPhimean))));  // Вычисляем угловое расстояние  fz = Sqrt(Pow(Sin(fdPhi / 2.0), 2) + Cos(EndLat * D2R) * Cos(StartLat * D2R) * Pow(Sin(fdLambda / 2.0), 2));  fz = 2 * ASin(fz);  // Вычисляем смещение  fAlpha = Cos(EndLat * D2R) * Sin(fdLambda) * 1 / Sin(fz);  fAlpha = ASin(fAlpha);  // Вычисляем радиус Земли  fR = (fRho * fNu) / ((fRho * Pow(Sin(fAlpha), 2)) + (fNu * Pow(Cos(fAlpha), 2)));  // Получаем смещение и расстояние  Distance = (fz * fR);  AbsAD = fAlpha * R2D;  if AbsAD < 0 then      AbsAD = AbsAD * (-1);  endif;    if ((StartLat < EndLat) and (StartLong < EndLong)) then      Bearing =  AbsAD;  elsif ((StartLat < EndLat) and (StartLong > EndLong)) then      Bearing = 360 - AbsAD;  elsif ((StartLat > EndLat) and (StartLong > EndLong)) then      Bearing = 180 + AbsAD;  elsif ((StartLat > EndLat) and (StartLong < EndLong)) then      Bearing = 180 - AbsAD;  endif;   КонецПроцедуры Процедура КнопкаВыполнитьНажатие(Кнопка)    Distance = 0;    Bearing = 0;        //StartLat = 55.791444; // Начальная широта    //StartLong = 37.665028; // Начальная долгота    //EndLat = 55.791589; // Конечная широта    //EndLong = 37.669148; // Конечная долгота        StartLat = 47.877987; // Начальная широта    StartLong = 35.071194; // Начальная долгота    EndLat = 47.877289; // Конечная широта    EndLong = 35.072243; // Конечная долгота            РасстояниеПоКоординатам(StartLat,StartLong,EndLat,EndLong,Distance,Bearing);        Сообщить("Distance (расстояние) = Расстояние в метрах " + Distance);    Сообщить("Bearing (смещение) = Смещение в градусах " + Bearing);     КонецПроцедуры

можешь сдуть с проекта: freefin.ucoz.org в восьмерошной демке модули форм не закрыты.

Пароля нет на 1CGPS8distr.rar  :(

вообщето в все сказано и за атаны и даже за ПИ

чего ругается, типо нету метода ?

тогда

На продажную версию пароль. Демку качай. В ней нет возможности добавления регионов, а в остальном и математика и хранение и слив с устройств да всё остальное присутствуют. Модально.

.........  А чего это вы тут делаете?!.........

Плюшки считаем