Суббота, 23 Ноя 2024, 01:35
Uchi.ucoz.ru
Меню сайта
Форма входа

Категории раздела
Высшая математика [11]
Экономическая социология [95]
Основы Менеджмента [64]
Бухгалтерский учёт [157]
Философия [163]
Мировая Экономика [603]
Бизнес планирование [29]
Финансирование и кредитование инвест [105]
Ценообразование [46]
Гражданское право [196]
Права Человека [173]
Основы Маркетинга [207]
Основы энергосбережения [55]
Информатика [0]
Экология и устойчивое развитие [0]
Физика для студентов [0]
Основы права [0]
Политология [0]
Не стандартные примеры на Delphi [169]
Примеры на Delphi7 [108]
Алгоритмы [94]
API [110]
Pascal [152]
Базы Данных [6]
Новости
Чего не хватает сайту?
500
Статистика
Зарегистрировано на сайте:
Всего: 51646


Онлайн всего: 14
Гостей: 14
Пользователей: 0
Яндекс.Метрика
Рейтинг@Mail.ru

Каталог статей


Главная » Статьи » Студентам » Алгоритмы

Как вычислить расстояние, имея широту и долготу.
Попробуйте следующий код. Я им пользуюсь продолжительное
время.
Входные данные:
StartLat (начальная широта) = Градусы и сотые
доли
StartLong (начальная долгота) = Градусы и сотые доли
EndLat (конечная
широта) = Градусы и сотые доли
EndLong (конечная долгота) = Градусы и сотые
доли

Выходные данные:
Distance (расстояние) = Расстояние в метрах
Bearing
(смещение) = Смещение в градусах

Не забудьте включить модуль Math в список используемых (USES) модулей.

var

// Передаваемые широта/долгота в градусах и
сотых долях
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; //
Вычисленное от и до смещениеEnd

const

// Константы, используемые для вычисления
смещения и расстояния
D2R: Double = 0.017453; // Константа для преобразования градусов в радианы
R2D: Double = 57.295781; // Константа для преобразования радиан в градусы
a: Double = 6378137.0; // Основные полуоси
b: Double = 6356752.314245; // Неосновные полуоси
e2: Double = 0.006739496742337; // Квадрат эксцентричности эллипсоида
f: Double = 0.003352810664747; // Выравнивание эллипсоида
begin

// Вычисляем разницу между двумя долготами и
широтами и получаем среднюю широту
fdLambda := (StartLong - EndLong) * D2R;
fdPhi := (StartLat - EndLat) * D2R;
fPhimean := ((StartLat + EndLat) / 2.0) *
D2R;

// Вычисляем меридианные и поперечные радиусы
кривизны средней широты
fTemp := 1 - e2 *
(Power(Sin(fPhimean),2));
fRho := (a * (1 - e2)) / Power(fTemp,
1.5);
fNu := a / (Sqrt(1 - e2 * (Sin(fPhimean) *
Sin(fPhimean))));

// Вычисляем угловое расстояние
fz :=Sqrt(Power(Sin(fdPhi/2.0),2)+Cos(EndLat*D2R)*Cos(StartLat*D2R)*Power(Sin(fdLambda/2.0),2)) ;

fz := 2 * ArcSin(fz);

// Вычисляем смещение
fAlpha := Cos(EndLat * D2R) * Sin(fdLambda) * 1 / Sin(fz);
fAlpha := ArcSin(fAlpha);

// Вычисляем радиус Земли
fR := (fRho * fNu) / ((fRho * Power(Sin(fAlpha),2)) + (fNu * Power(Cos(fAlpha),2)));

// Получаем смещение и расстояние
Distance := (fz * fR);
if((StartLat < EndLat) and (StartLong <
EndLong)) then
Bearing := Abs(fAlpha * R2D)
else if ((StartLat < EndLat) and (StartLong >
EndLong)) then
Bearing := 360 - Abs(fAlpha * R2D)
else if ((StartLat > EndLat) and (StartLong >
EndLong)) then
Bearing := 180 + Abs(fAlpha * R2D)
else if ((StartLat > EndLat) and (StartLong <
EndLong)) then
Bearing := 180 - Abs(fAlpha *
R2D);end;
Категория: Алгоритмы | Добавил: Lerka (21 Ноя 2012)
Просмотров: 1126 | Рейтинг: 1.3/ 7 Оштрафовать | Жаловаться на материал
Похожие материалы
Всего комментариев: 0

Для блога (HTML)


Для форума (BB-Code)


Прямая ссылка

Профиль
Суббота
23 Ноя 2024
01:35


Вы из группы: Гости
Вы уже дней на сайте
У вас: непрочитанных сообщений
Добавить статью
Прочитать сообщения
Регистрация
Вход
Улучшенный поиск
Поиск по сайту Поиск по всему интернету
Наши партнеры
Интересное
Популярное статьи
Портфолио ученика начальной школы
УХОД ЗА ВОЛОСАМИ ОЧЕНЬ ПРОСТ — ХОЧУ Я ЭТИМ ПОДЕЛИТ...
Диктанты 2 класс
Детство Л.Н. Толстого
Библиографический обзор литературы о музыке
Авторская программа элективного курса "Практи...
Контрольная работа по теме «Углеводороды»
Поиск
Главная страница
Используются технологии uCoz