/*ZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZZ
	Функция определяет расстояние между точками 1 и 2 по их координатам.
	X - Широта, Y - Долгота, 6371000 - радиус земли в метрах
	31.03.2015
*/
ALTER FUNCTION [dbo].[fnAddressGeoDist] ( @Pos1X DECIMAL(9,6), 
					  @Pos1Y DECIMAL(9,6), 
					  @Pos2X DECIMAL(9,6), 
					  @Pos2Y DECIMAL(9,6)
) RETURNS int
AS
BEGIN
	DECLARE @Dist int
/*
	cos(d) = sin(φА)·sin(φB) + cos(φА)·cos(φB)·cos(λА − λB),
		где φА и φB — широты, λА, λB — долготы данных пунктов, d — расстояние между пунктами, измеряемое в радианах длиной дуги большого круга земного шара. 
	Расстояние между пунктами, измеряемое в километрах, определяется по формуле:
		L = d·R, где R = 6371 км — средний радиус земного шара.
*/
	SET @Dist = 6371000 * 2 * ASIN(SQRT(
      POWER(SIN((@Pos1X - ABS(@Pos2X)) * PI()/180 / 2), 2) +
      COS(@Pos1X * PI()/180) *
      COS(ABS(@Pos2X) * PI()/180) *
      POWER(SIN((@Pos1Y - @Pos2Y) * PI()/180 / 2), 2)
    ))
	RETURN @Dist
END