// Переменные, используемые для вычисления смещения и расстояния 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)) ;