SOFA - тестирование, продолжение
Начало -
SOFA.iauAtoc13 - тестирование программы IAU SOFA - The International Astronomical Union's,
Standards Of Fundamental Astronomy,
Международный астрономический союз,
Стандарты Фундаментальной Астрономии
http://www.iausofa.org/index.htmlВсем привет!
Условия тестирования
Исследуются точки наблюдения на поверхности Земли.
Сеть точек наблюдения -
Longitude - от 0 до 360 с дискретностью 10,
Latitude - от -90 до +90 с дискретностью 10 -
37 х 19 = 704 точки
В каждой точке наблюдения
исследуются точки полусферы,
видимой из этой точки наблюдения -
R - от 0 до 360 с дискретностью 10,
D - от 0 до 180 с дискретностью 10.
37 х 19 = 704 точки
Итого: 704 х 704 = 494 209 точек
Для каждой исследуемой точки наблюдения
координаты каждой исследуемой точки полусферы
принимаются исходными тестируемыми
ICRS - координатами.
Тестирование состоит в следующем.
1. Программа SOFA.iauAtco13 -
по координатам точки наблюдения и
ICRS - координатам очередной точки полусферы,
{ исходные ICRS - координаты }
определяются горизонтальные координаты
этой точки полусферы.
2. Программа SOFA.iauAtoc13 -
полученные горизонтальные координаты
преобразовываются обратно в
ICRS - координаты
( результативные )
3. Результативные ICRS - координаты сравниваются
с исходными ICRS - координатами,
разница между ними записывается в
соответствующий протокол.
В идеале разница
между результативными ICRS - координатами
и исходными ICRS - координатами
должна быть минимальной,
обусловленной только недоучетом атмосферных условий.
Однако результаты эксперимента обескураживают -
по полусферам всех точек наблюдения
разница по прямому восхождению колеблется от -360° до 360°,
разница по склонению от -270° до 90° (?),
Ниже представлены три СОКРАЩЕННЫЕ ТАБЛИЦЫ,
с округленными до целых результатами расчетов
в трех точках наблюдения -
на Северном полюсе, на Южном полюсе и в районе Лондона.
На Северном полюсе
разница по ВСЕЙ ПОЛУСФЕРЕ ВИДИМОСТИ составляет
по прямому восхождению
от -360° при
Longitude = 70, Latitude = 90, R_ICRS_Src = 0, D_ICRS_Src = 90,
до 360° при
Longitude = 240, Latitude = 90, R_ICRS_Src = 360, D_ICRS_Src = 0,
по склонению
от -90° до 90°.
NORTH POLE, Longitude = 0°, Latitude = 90°
1 2 3 4 5 6 7 8 9 10 11
1 Ra_ICRS_Src 0 0 0 0 0 0 0 0 0 0 0
2 Dec_ICRS_Src 0 20 40 60 80 90 100 120 140 160 180
3 Ra_Horizont 360 360 0 0 0 171 180 180 180 180 180
4 Dec_Horizont 0 20 40 60 80 90 80 60 40 20 0
5 Ra_ICRS_Res 330 100 100 100 100 289 280 280 280 280 338
6 Dec_ICRS_Res 90 70 50 30 10 0 10 30 50 70 90
7 Ra_ICRS_Delta 330 100 100 100 100 289 280 280 280 280 338
8 Dec_ICRS_Delta 90 50 10 -30 -70 -90 -90 -90 -90 -90 -90
На Южном полюсе
разница по ВСЕЙ ПОЛУСФЕРЕ ВИДИМОСТИ составляет
по прямому восхождению
от -360° при
Longitude = 270, Latitude = -90, R_ICRS_Src = 360, D_ICRS_Src = 90,
до 360° при
Longitude = 270, Latitude = -90, R_ICRS_Src = 0, D_ICRS_Src = 90,
по склонению - от -270° до -90°.
SOUTH POLE, Longitude = 0°, Latitude = -90°
1 2 3 4 5 6 7 8 9 10 11
1 Ra_ICRS_Src 0 0 0 0 0 0 0 0 0 0 0
2 Dec_ICRS_Src 0 20 40 60 80 90 100 120 140 160 180
3 R_Horizont 360 360 0 0 0 171 180 180 180 180 180
4 D_Horizont 0 20 40 60 80 90 80 60 40 20 0
5 R_ICRS_Res 33 280 280 280 280 92 100 100 100 100 28
6 D_ICRS_Res -90 -70 -50 -30 -10 0 -10 -30 -50 -70 -90
7 R_ICRS_Delta 33 280 280 280 280 92 100 100 100 100 28
8 D_ICRS_Delta -90 -90 -90 -90 -90 -90 -110 -150 -190 -230 -270
В районе Лондона
разница по ВСЕЙ ПОЛУСФЕРЕ ВИДИМОСТИ составляет
по прямому восхождению
от -344° при
Longitude = 0, Latitude = 50, R_ICRS_Src = 360, D_ICRS_Src = 40,
до 328° при
Longitude = 0, Latitude = 50, R_ICRS_Src = 30, D_ICRS_Src = 40,
по склонению от -130° до 50°.
LONDON, Longitude = 0°, Latitude = 50°
1 2 3 4 5 6 7 8 9 10 11
1 Ra_ICRS_Src 0 0 0 0 0 0 0 0 0 0 0
2 Dec_ICRS_Src 0 20 40 60 80 90 100 120 140 160 180
3 Ra_Horizont 360 360 0 0 0 171 180 180 180 180 180
4 Dec_Horizont 0 20 40 60 80 90 80 60 40 20 0
5 Ra_ICRS_Res 280 280 16 100 100 292 280 280 280 280 280
6 Dec_ICRS_Res 50 70 90 70 50 -39 -30 -10 10 30 50
7 Ra_ICRS_Delta 280 280 16 100 100 292 280 280 280 280 280
8 Dec_ICRS_Delta 50 50 50 10 -30 -129 -130 -130 -130 -130 -130
Вывод.
Программа SOFA.iauAtco13 работает правильно,
при некоторых незначительных замечаниях,
программа SOFA.iauAtoc13 РАБОТАЕТ НЕПРАВИЛЬНО.
Хотелось бы ошибаться...
У меня задача -
горизонтальные координаты космических объектов
перевести в ICRS - координаты,
именно та, которая не работает.
Можно, конечно, как обычно советуют,
"почитать книжки" и слепить свое.
Но в данном случае кустарщина не желательна...
Попробую SOFA побеспокоить, интересно, что скажут...