Как преобразовать координаты?

[ORSA]

Как преобразовать экваториальные координаты в азимутальные - поделитесь алгоритмом, а лучше программной библиотекой если есть. Янашел ряд формул в инете но они плохо работают - у тангенса период  PI и нужно корректировать результат вычисления функции ATAN().

Пишу сейчас для себя программку для автоматизации подготовки плана наблюдения:
при отборе из каталога будет учитываться:
- время наблюдения
- геометрия участка (углы на которые закрывается небо строениями и деревьями)
- яркость объектов

[Ivan Mhitarov]

На чем пишем?
Именно такая программка уже есть - AstroPlanner.
Если интересуют формулы - вот кусок из моей программы

Код: [Выделить]

{Эта функция из Монтенбрука}
Function LMST(MJD,Lambda:Extended):Extended;
Var MJD0,T,UT,GMST:Extended;

 Function Frac(X:Extended):Extended;
   Begin  X:=X-Trunc(X); If X<0 Then X:=X+1; Frac:=X; End;

Begin
  MJD0:=Int(MJD);
  UT:=(MJD-MJD0)*24.0;
  T:=(MJD0-51544.5)/36525.0;
  GMST:=6.697374558+1.0027379093*UT
        +(8640184.812866+(0.093104-6.2E-6*T)*T)*T/3600.0;
  LMST:=24.0*Frac((GMST-Lambda/15.0)/24.0);
End;

Function Atn2(Y,X:Extended):Extended;
  Var AX,AY,FI:Extended;
Begin
 If ( X = 0.0 ) AND ( Y = 0.0 )
   Then Atn2:=0.0
   Else
     Begin
       AX:=Abs(X);
       AY:=Abs(Y);
       If ( AX > AY )
         Then FI:=ArcTan(AY/AX)/Rad
         Else FI:=90.0-ArcTan(AX/AY)/Rad;
       If X < 0.0 Then FI:=180.0-FI;
       If Y < 0.0 Then FI:=-FI;
       Atn2:=FI;
     End;
End;



    // Кусок моей программы
    // Sn, Cs - Синус и Косинус от угла в градусах
    // ceLat.Value - широта, ceLong.Value - долгота
    S:=LMST(LocalDateTimeToDateTime(ATime)+15018, -ceLong.Value);
    T:=15*(S-Ra);
    El:=90-Acs(Cs(ceLat.Value)*Cs(Dec)*Cs(T)+Sn(ceLat.Value)*Sn(Dec));
    Az:=Atn2(Cs(Dec)*Sn(T), Cs(Dec)*Sn(ceLat.Value)*Cs(T)-Cs(ceLat.Value)*Sn(Dec));

    Az:=Az+180;   // Отсчитывается от точки юга
    If Az >= 360 Then
      Az:=Az-360;


[recarrion]

Паскаль! До сих пор приятно на листинги смотреть.

[ORSA]

О коллеги - спасибо за листинг. 
Пишу на самом доступном инструменте всех времен и народов ЁКСЕЛЕ-МОКСЕЛЕ :-)

[ORSA]

Фу преобразовал в VBA - работает

Const rad = (2 * 3.14159265358979 / 360)

Function Get_Az(S As Double, Shirota As Double, RA As Double, dec As Double) As Double
Rem Shirota  широта в градусах
Rem S звездное время в часах
Rem RA в в часах
Rem dec = склонение в градусах
Dim t As Double
   t = 15 * (S - RA)
   Az = Atn2(Cos(to_rad(dec)) * Sin(to_rad(t)), Cos(to_rad(dec)) * Sin(to_rad(Shirota)) * Cos(to_rad(t)) -    Cos(to_rad(Shirota)) * Sin(to_rad(dec)))
   Az = Az + 180
   If Az >= 360 Then Az = Az - 360
   Get_Az = Az
End Function

Function Get_Alt(S As Double, Shirota As Double, RA As Double, dec As Double) As Double
Rem Shirota  широта в градусах
Rem S звездное время в часах
Rem RA в в часах
Rem dec = склонение в градусах
Dim t As Double
   t = 15 * (S - RA)
   Get_Alt = 90 - to_deg(acos(Cos(to_rad(Shirota)) * Cos(to_rad(dec)) * Cos(to_rad(t)) + Sin(to_rad(Shirota)) *    Sin(to_rad(dec))))
End Function

Function frac(x As Double) As Double
   x = x - Round(x)
   If x < 0 Then x = x + 1
   frac = x
End Function

Function LMST(MJD As Double, Lambda As Double) As Double
Dim MJD0 As Double, t As Double, UT As Double, GMST As Double

  MJD0 = Int(MJD)
  UT = (MJD - MJD0) * 24#
  t = (MJD0 - 51544.5) / 36525#
  GMST = 6.697374558 + 1.0027379093 * UT + (8640184.812866 + (0.093104 - 0.0000062 * t) * t) * t / 3600#
  LMST = 24# * frac((GMST - Lambda / 15#) / 24#)
 
 
End Function


Function Atn2(Y, x As Double) As Double
Dim AX As Double, AY As Double, FI As Double
If (x = 0#) And (Y = 0#) Then
 Atn2 = 0#
Else
       
    AX = Abs(x)
    AY = Abs(Y)
    If (AX > AY) Then
        FI = Atn(AY / AX) / rad
    Else
        FI = 90# - Atn(AX / AY) / rad
    End If
       
    If x < 0# Then FI = 180# - FI
    If Y < 0# Then FI = -FI
    Atn2 = FI
     
End If

End Function

Function to_rad(x As Double) As Double
   to_rad = x * rad
End Function

Function to_deg(x As Double) As Double
   to_deg = x / rad
End Function

Function acos(x As Double) As Double
   Sheets("s_acos").Cells(5, 1).Value = x
   acos = Sheets("s_acos").Cells(5, 2).Value
End Function

[ORSA]

Ну вот наконец то все и получилось! Чем приятен EXCEL, что теперь программку можно легко доводить и модифицировать под любые желания J