ВНИМАНИЕ! На форуме началось голосование в конкурсе - астрофотография месяца МАРТ!
0 Пользователей и 1 Гость просматривают эту тему.
Но вопрос в алгоритме, как определить?
void JnowtoJ2000(double appRA, double appDecl, double *r, double *dc){ double ra=0., dec=0., utc1, utc2, tai1, tai2, tt1, tt2, fd, eo, ri; int y, m, d, H, M; DBG("appRa: %g'', appDecl'': %g", appRA, appDecl); appRA *= DS2R; appDecl *= DAS2R;#define SOFA(f, ...) do{if(f(__VA_ARGS__)){WARNX("Error in " #f); goto rtn;}}while(0) // 1. convert system JDate to UTC SOFA(iauJd2cal, JDate, 0., &y, &m, &d, &fd); fd *= 24.; H = (int)fd; fd = (fd - H)*60.; M = (int)fd; fd = (fd - M)*60.; SOFA(iauDtf2d, "UTC", y, m, d, H, M, fd, &utc1, &utc2); SOFA(iauUtctai, utc1, utc2, &tai1, &tai2); SOFA(iauTaitt, tai1, tai2, &tt1, &tt2); iauAtic13(appRA, appDecl, tt1, tt2, &ri, &dec, &eo); ra = iauAnp(ri + eo); ra *= DR2S; dec *= DR2AS; DBG("SOFA: r=%g'', d=%g''", ra, dec);#undef SOFArtn: if(r) *r = ra; if(dc) *dc = dec;}
Всем привет!Дано -Точка наблюдения - F - склонение, широта, A - прямое восхождение, долготаНебесный объект -f - склонение, широта,a - прямое восхождение, долготаДата / Время2021.01.01 00:00:00Определить экваториальные координатынебесного объектасогласно эпохе J2000.0