Добрый день. Извиняюсь, что уже с набившей оскомину темой. Есть входящие координаты в градусах на WGS 84, их пересчитала в Гаусса-Крюгера также в градусах.
Конечный вариант нужен в метрах и тут возникли сложности. Пользовалась формулами ГОСТА Р 51794-2001 (25-26 формулы).
Конечный вид в T-SQL следующий
set @Lr= Pi()*(@Ld-@dL/3600)/180 'перевод градусов в радианы
set @Br= Pi()*(@Bd-@dB/3600)/180 'перевод градусов в радианы
set @a= (6+@Lr)/6
set @n= convert(int,@a)
set @l=(@Lr-(3+6*(@n-1)))/57.29577951
set @Bm= 6367558.4968*@Br - sin(2*@Br)*(16002.89+66.9607*power(sin(@Br),2) + 0.3515*power(sin(@Br),4) - (power(@l,2)*(1594561.25 + 5336.535*power(sin(@Br),2) + 26.79*power(sin(@Br),4) + 0.149*power(sin(@Br),6) + power ((@l),2)*(672483.4 - 811219.9*power(sin(@Br),2) + 5420*power(sin(@Br),4)- 10.6*power(sin(@Br),6) + power(@l,2)*(278194 -(830174*power(sin(@Br),2)) + 572434*power(sin(@Br),4)-16010*power(sin(@Br),6) + power((@l),2)*(109500 - 574700*power(sin(@Br),2) + 863700*power(sin(@Br),4) - 398600*power(sin(@Br),6)))))))
set @Lm= (5+10*@n)*1000000 + @l*cos(@Br)*(6378245 +
21346.1415*power(sin(@Br),2) +107.159*power(sin(@Br),4) + 0.5977*power(sin(@Br),6) + power((@l),2)*(1070204.16 -
2136826.66*power(sin(@Br),2) +17.98*power(sin(@Br),4) -
11.99*power(sin(@Br),6) + power((@l),2)*(270806 - 1523417*power(sin(@Br),2) + 1327645*power(sin(@Br),4) -
21701*power((@Br),6) + power((@l),2)*(79690 -
866190*power(sin(@Br),2) + 1730360*power((@Br),4) -
945460*power((@Br),6)))))
return @Bm
return @Lm
В результате получается несоответствие. В чем ошибка, может быть изначально не те формулы брала?
Спасибо за помощь!!!
