Есть лог GPS-логера с координатами. Нужно рассчитать сколько проехала машина.
Поиски привели к "формуле гаверсинуса" и "расстояние большого круга"
http://gis-lab.info/qa/great-circles.html
http://ru.w3support.net/index.php?db=so&id=27928
http://js-php.ru/web-development/distance-from-dot-to-dot/

Вот только как это на 1С перевести...

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

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

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

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

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

atan2 нет в 1С

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

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

Типа этого? http://megocode3.wordpress.com/2008/02/05/haversine-formula-in-c/
и по ADO вызывать?

(9) Достаточно http://www.1cpp.ru/docum/icpp/html/Math.html

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

(10) Восьмерка

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

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

(13) (14) Примеры можно?

(15) В поиске забанили? http://yandex.ua/yandsearch?text=java+math&lr=10347

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

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

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

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

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

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

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

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

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

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

По Формуле гаверсинусов
http://www.movable-type.co.uk/scripts/latlong.html и пример в Excel из первой ссылки

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

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

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

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

http://clip2net.com/s/vgcX

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

http://www.delphisources.ru/pages/faq/base/calculate_distance.html

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

Пля

Скрипт=Новый COMОбъект("MSScriptControl.ScriptControl");
Скрипт.language="javascript";
Математика=Скрипт.Eval("Math");
Сообщить(Математика.atan2(10,5));


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

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

(0)
http://vasnake.blogspot.com/2006/12/blog-post_15.html

+(38) http://dxdy.ru/topic19762.html

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

Может я че то не понимаю, но

Результат = sqrt((x2-x1)*(x2-x1)) + (у2-у1)(у2-у1))

+(39) В общем думаю самый простой вариант это взять длину дуги градуса широты и долготы для местности, погрешность будет, но не большая.
Тогда длина дуги между двумя точками будет такая:

длина пути = ((lon_1*d_Lat)^2+(lon_1*d_Lon)^2)^(1/2)

где:
lat_1 - длина одного градуса широты для местности
lon_1 - длина одного градуса долготы для местности
d_Lat=lat1-lon2 - разница широт
d_Lon=lon1-lon2 - разница долгот.

(18) ну и какие там проблемы вычисления ?

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)

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

(43) Не сходится

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

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

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

(35) Ругается

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

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

http://maps.google.com/ - проверь линейкой

(51) с (27)

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

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

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

Уряяяяя.


Функция 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))));
   Сообщить(Расстояние);
   
КонецПроцедуры

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

http://api.yandex.ru/maps/jsapi/doc/ref/reference/geocoordsystem.xml
rulerDistance

{Number} rulerDistance(point1, point2)


Возвращает расстояние между двумя заданными точками (в метрах), если двигаться вдоль прямой соединяющей точки на карте. Используется для расчета расстояний в элементе управления "Линейка".

Параметры:Имя    Тип    Описание
point1    YMaps.GeoPoint    Первая точка.
point2    YMaps.GeoPoint    Вторая точка.


Возвращает:Тип    Описание
Number    Расстояние "по линейке" между двумя заданными точками.

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

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

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

Вот из (34) переделал. вроде работает

Процедура РасстояниеПоКоординатам(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);
   
КонецПроцедуры

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

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

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

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

тогда

Скрипт=Новый COMОбъект("MSScriptControl.ScriptControl");
Скрипт.language="javascript";
Сообщить(Скрипт.Eval("Math.atan2(10,5)"));
Сообщить(Скрипт.Eval("Math.PI"));


(65) На продажную версию пароль. Демку качай. В ней нет возможности добавления регионов, а в остальном и математика и хранение и слив с устройств да всё остальное присутствуют. Модально.
http://freefin.ucoz.org/1C81/Demo/1CGPS8demo.zip

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

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