Потратил неделю на анализ данных от RainDog.
ВОт определенные мной по его фото координаты TYC 4629-37-1, на которую он указал:
1 045 274
2 061 281
3 078 286
4 094 291
5 111 294
6 128 297
7 145 298
8 162 298
9 179 297
10 196 295
11 213 292
12 230 287
13 246 282
14 262 275
15 277 268
16 292 259
17 306 250
18 320 239
Начал с определения масштаба.
ПОследнее положение TYC 4629-37-1 (320 239) .
Последнее положение ПОлярной (631 120).
Расстояние между ними в пикселях = sqrt((631-320)^2+(239-120)^2)=333.
Координаты на 2000г
TYC 4629-37-1 4ч42м46,2с 89*37'48,3"
Полярная 2ч31м48,7с 89*15'51"
Вычислил угол =1690".
Масштаб получился довольно крупный = 5,075" на пиксель (прошлый раз у RainDog было 1,215, почти в 4 раза мельче).
Далее - составление программ и "автоматизированный" анализ массива данных.
Разбиваю все 18 точек на 3 блока по 6, в каждом блоке составляю 20 комбинаций по 3 точки, определяющие радиус трека и координаты центра.
Затем вычисляю средний радиус трека и координаты центра.
Например, для среднего блока с 7 по 12 получил (на первую колонку не смотрите, это просто код комбинации в блоке, X Y в пикселях, радиус в секундах дуги, все числа - с множителем 1000):
1.0e+003 * X Y R"
0.1230 0.1530 0.0080 1.4723
0.1240 0.1530 0.0070 1.4774
0.1250 0.1530 0.0060 1.4825
0.1260 0.1530 0.0290 1.3658
0.1340 0.1530 0.0060 1.4825
0.1350 0.1530 0.0040 1.4926
0.1360 0.1540 0.0360 1.3304
0.1450 0.1530 0.0030 1.4977
0.1460 0.1550 0.0480 1.2698
0.1560 0.1600 0.0800 1.1090
0.2340 0.1530 0.0050 1.4877
0.2350 0.1530 0.0030 1.4978
0.2360 0.1550 0.0410 1.3048
0.2450 0.1520 0.0010 1.5081
0.2460 0.1570 0.0560 1.2284
0.2560 0.1630 0.0930 1.0404
0.3450 0.1520 0.0000 1.5135
0.3460 0.1600 0.0690 1.1611
0.3560 0.1690 0.1110 0.9453
0.4560 0.1770 0.1390 0.7976
srednee =
1.0e+003 * X Y R"
0.2153 0.1566 0.0372 1.3232
Видно, что в триадах получился большой разброс центра по Y (от 0 до 139 пикселей) и по радиусу (от 797" до 1508") .
Это из-за того, что точки на траектории получены довольно близко.
Поэтому взял "черезстрочную развертку" - из всех четных точек сделал два блока: 2-12 и 8-18.
ВОт что получилось:
точки 2-4-6-8-10-12
1.0e+003 * X Y R"
0.1230 0.1610 0.0100 1.4660
0.1240 0.1550 0.0290 1.3650
0.1250 0.1550 0.0290 1.3650
0.1260 0.1540 0.0320 1.3489
0.1340 0.1520 0.0460 1.2789
0.1350 0.1540 0.0380 1.3205
0.1360 0.1540 0.0390 1.3157
0.1450 0.1550 0.0280 1.3697
0.1460 0.1540 0.0340 1.3394
0.1560 0.1540 0.0400 1.3110
0.2340 0.1510 0.0620 1.1976
0.2350 0.1540 0.0460 1.2801
0.2360 0.1540 0.0450 1.2850
0.2450 0.1550 0.0280 1.3702
0.2460 0.1540 0.0360 1.3295
0.2560 0.1540 0.0430 1.2949
0.3450 0.1530 0.0070 1.4772
0.3460 0.1520 0.0260 1.3807
0.3560 0.1540 0.0420 1.3008
0.4560 0.1570 0.0560 1.2284
srednee =
1.0e+003 * X Y R"
0.2153 0.1543 0.0358 1.3312
-----------------
точки 8-10-12-14-16-18
1.0e+003 * X Y R"
0.1230 0.1570 0.0550 1.2336
0.1240 0.1560 0.0460 1.2794
0.1250 0.1560 0.0430 1.2946
0.1260 0.1560 0.0410 1.3048
0.1340 0.1540 0.0370 1.3254
0.1350 0.1540 0.0370 1.3254
0.1360 0.1540 0.0370 1.3254
0.1450 0.1540 0.0370 1.3254
0.1460 0.1540 0.0370 1.3254
0.1560 0.1540 0.0360 1.3304
0.2340 0.1510 0.0270 1.3791
0.2350 0.1510 0.0310 1.3591
0.2360 0.1520 0.0320 1.3533
0.2450 0.1530 0.0340 1.3424
0.2460 0.1530 0.0350 1.3374
0.2560 0.1530 0.0350 1.3374
0.3450 0.1540 0.0380 1.3212
0.3460 0.1540 0.0370 1.3261
0.3560 0.1540 0.0360 1.3309
0.4560 0.1530 0.0350 1.3377
srednee =
1.0e+003 * X Y R"
0.2153 0.1538 0.0373 1.3247
РАзброс меду отдельными вычислениями заметно уменьшился, а средние величины почти не изменились, что говорит о практической эффективности усреднения всех 20 сочетаний даже для "полной ахинеи" вначале при малом расстоянии между точками.
Итого, получился радиус трека от 1323" (22'03") до 1331" (22'11").
Это хорошо в сравнении с координатами 2000г (полярное расстояние 22'12"), но вряд ли хорошо для текущих координат с учетом прецессии и нутации Земной Оси, которые пока найти не могу, потому что в Carte de Sel и Selesta TYC4629-37-1 нет, а в СТеллариуме она есть, но текущие координаты только средние
5ч38м1с 89*38'56" (с учетом прецессии но без нутации, а ведь это +-17"), этому соответствует радиус трека 21'04" что на минуту меньше рассчитанного по фото.
Привожу расчитанный мной график изменеия радиуса трека (полярного расстояния) с учетом аберрации. По горизонтали время года в месяцах, по вертикали аберрационное отклонение от среднего расстояния до Полюса.
Поскольку звезда почти на 6 часов, график по фазе противоположен известному графику для Гаммы Дракона, которая на 18 часов.
Для конца сентября должно быть максимальное отклонение ОТ Полюса (максимальный радиус треков).
Качественно вроде бы правильно, а количество пока дает повод для разночтений.