Отправляет email-рассылки с помощью сервиса Sendsay

Delphi для профессионалов

Рассылка закрыта

При закрытии подписчики были переданы в рассылку "Delphi - проблемы и решения" на которую и рекомендуем вам подписаться.

Вы можете найти рассылки сходной тематики в Каталоге рассылок.

  Октябрь 2002   →
 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
24.10.2002
15:05:04
23:21:18
25
26
27
28
29
30
31
Автор
Ведущий

Статистика

8.943 подписчиков
-22 за неделю
←   Все выпуски   →

Delphi для профессионалов


Информационный Канал Subscribe.Ru
eServer.ru :: Хостинг без ограничений!
Рассылка: Delphi для профессионалов

Количество подписчиков: 12
Рассылка 'Delphi для профессионалов'

Вычисление расстояния между двуми точками, имея широту и долготу

 

Входные данные:
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;

Лирическое отступление автора:
в качестве входных параметров используются ШИРОТЫ (в множественном числе, ударение на втором слоге), ведь их две.
Но хмммм.... долгота(ы???) тоже две, а как будет звучать множественное число? Загадка.
Наверное не существует такой формы. (P.S. зато я знаю как будет множественное число от слова ДНО! Слабо?)

Все вопросы присылайте по адресу: delphi-ask@igp.org.ua
Все ответы присылайте по адресу: delphi-reply@igp.org.ua
Все подсказки присылайте по адресу: delphi-hint@igp.org.ua
Со всем списком советов можно познакомиться по адресу: http://igp.org.ua/delphine/.
Наш проект Delphine можно найти по следующей ссылке: http://igp.org.ua/products/delphine/.

Заходите в наш форум: Форум Delphi-профессионалов. Адрес: http://igp.org.ua/cgi-bin/yabb/yabb.pl.
Здесь Вы найдете ответы на все Ваши вопросы!
http://subscribe.ru/
E-mail: ask@subscribe.ru 		Отписаться
Убрать рекламу
В избранное