Тематические рабочие группы > Вопросы астрометрии
Программа для проверки точности астрометрии на короткой дуге
ctac:
Есть и интернет вариант: http://www.fitsblink.net/residuals/index.html
ctac:
Загружаю сюда программу astrostat
ctac:
--- Цитата: Vit от 05 июля 2007, 23:47:18 ---А насколько можно доверять этой программе?
Я тут загнал в нее наблюдения кометы 2006VZ13, вот что она выдала:
2006vz13.txt
Number of observations 10
UTC,hms RA,hms dRA,s DEC,g`" dDEC,", delta,"
30.96212 20 54 8.490 -0.022 66 25 0.00 -0.90 0.91
30.96321 20 54 7.320 0.074 66 25 5.30 0.22 0.50
30.96411 20 54 6.160 -0.039 66 25 8.90 0.37 0.43
30.96486 20 54 5.250 -0.075 66 25 12.40 0.98 1.08
30.96567 20 54 4.460 0.075 66 25 14.90 0.37 0.58
30.96647 20 54 3.500 0.047 66 25 16.70 -0.91 0.95
30.96730 20 54 2.500 0.011 66 25 20.60 -0.21 0.22
30.96802 20 54 1.530 -0.123 66 25 23.80 0.23 0.77
30.96890 20 54 0.620 -0.008 66 25 27.10 0.15 0.15
30.96975 20 53 59.700 0.059 66 25 29.90 -0.30 0.47
------ ------
Sigma, arcsec 0.40 0.59
Честно говоря, Астрометрика дает для этих наблюдений максимальную погрешность, к примеру, по Decl. 0.35”. Сравнивая свои наблюдательные данные по астероидам с Насавским калькулятором эфемерид, я не отмечал даже на самых плохих снимках расхождений больше чем в 0.5”, а для выборочных это значение укладывалось в 0.2”. Но чтобы ошибка достигала 1”, такого не было. Мне кажется, что выдаваемые данные в этой программе сильно зависят от количества наблюдений и (или) длины дуги, неплохо бы оценить минимальные значения этих параметров.
--- Конец цитаты ---
Может тут ошибка в том, что не учтен косинус угла склонения? Надо будет проверить.
ctac:
Для проверки решил попробовать чужую астрометрию:
--- Код: --- K09SH2S C2009 09 29.27435 00 50 16.87 +20 27 17.3 20.1 V ES141G96
K09SH2S C2009 09 29.28127 00 50 15.03 +20 27 50.6 20.4 V ES141G96
K09SH2S C2009 09 29.28819 00 50 13.23 +20 28 24.0 20.6 V ES141G96
K09SH2S C2009 09 29.29513 00 50 11.36 +20 28 57.2 20.5 V ES141G96
K09SH2S C2009 09 29.35716 00 49 55.05 +20 33 52.4 19.8 V ES141G96
K09SH2S C2009 09 29.35780 00 49 54.87 +20 33 55.6 20.1 V ES141G96
K09SH2S C2009 09 29.35843 00 49 54.73 +20 33 58.5 20.3 V ES141G96
K09SH2S C2009 09 29.35905 00 49 54.55 +20 34 01.7 20.4 V ES141G96
--- Конец кода ---
http://www.cfa.harvard.edu/mpec/K09/K09SE1.html
Вот что получил с нашей программой:
--- Код: --- Number of observations 8
UTC,hms RA,hms dRA,s DEC,g`" dDEC,", delta,"
29.27435 0 50 16.870 0.011 20 27 17.30 -0.43 0.46
29.28127 0 50 15.030 -0.006 20 27 50.60 -0.12 0.14
29.28819 0 50 13.230 0.017 20 28 24.00 0.28 0.37
29.29513 0 50 11.360 -0.026 20 28 57.20 0.40 0.54
29.35716 0 49 55.050 0.001 20 33 52.40 -0.16 0.16
29.35780 0 49 54.870 -0.011 20 33 55.60 0.00 0.15
29.35843 0 49 54.730 0.015 20 33 58.50 -0.11 0.24
29.35905 0 49 54.550 -0.001 20 34 1.70 0.13 0.14
------ ------
Sigma, arcsec 0.20 0.26
--- Конец кода ---
--- Код: ---
Residuals in seconds of arc
090929 G96 0.0 0.0
090929 G96 0.2- 0.0
090929 G96 0.1+ 0.1+
090929 G96 0.4- 0.1-
090929 G96 0.1+ 0.3-
090929 G96 0.0 0.1-
090929 G96 0.3+ 0.2-
090929 G96 0.1+ 0.1+
--- Конец кода ---
Плохо коррелирует. Надо еще в других программах проверить.
Теперь из программы FIND_ORB:
.00+ .00+
.30- .02+
.02- .18+
.61- .09+
.04+ .35-
.13- .18-
.24+ .26-
.00+ .00+
Разница тут вот в чем:
- МРС считает всю орбиту и отклонения от ее дуги
- aststat считает полином и отклонения от него
- FIND_ORB считает это короткой дугой, а не прямой линией
Vit:
--- Цитировать ---Теперь из программы FIND_ORB:
--- Конец цитаты ---
Стас, сравнение не совсем корректно, т.к. ты взял из МРС резедуалы со всей дуги, а в Findorb загнал только часть.
Вот что получается если использовать всю дугу:
--- Код: ---09 09 29.27435 G96 00 50 16.87 +20 27 17.3 .00 .00 0.1312 1.1264
09 09 29.28127 G96 00 50 15.03 +20 27 50.6 .23- .05- 0.1312 1.1264
09 09 29.28819 G96 00 50 13.23 +20 28 24.0 .11+ .05+ 0.1312 1.1264
09 09 29.29513 G96 00 50 11.36 +20 28 57.2 .44- .09- 0.1313 1.1265
09 09 29.35716 G96 00 49 55.05 +20 33 52.4 .12+ .31- 0.1315 1.1267
09 09 29.35780 G96 00 49 54.87 +20 33 55.6 .05- .14- 0.1316 1.1267
09 09 29.35843 G96 00 49 54.73 +20 33 58.5 .31+ .22- 0.1316 1.1267
09 09 29.35905 G96 00 49 54.55 +20 34 01.7 .07+ .06+ 0.1316 1.1267
09 09 29.36020 H36 00 49 53.45 +20 33 59.3 .46+ 1.1+ 0.1316 1.1267
09 09 29.36439 H36 00 49 52.36 +20 34 18.0 .34+ .15+ 0.1316 1.1268
09 09 29.36857 H36 00 49 51.22 +20 34 36.8 .54- .65- 0.1316 1.1268
09 09 29.51466 G96 00 49 14.85 +20 46 03.5 .29- .41- 0.1323 1.1274
09 09 29.51530 G96 00 49 14.73 +20 46 06.9 .22+ .06+ 0.1323 1.1274
09 09 29.51593 G96 00 49 14.54 +20 46 09.6 .28- .12- 0.1323 1.1274
09 09 29.54976 F65 00 49 07.66 +20 49 10.6 .06+ .12+ 0.1325 1.1275
09 09 29.55071 F65 00 49 07.42 +20 49 14.9 .13+ .06+ 0.1325 1.1275
09 09 29.55160 F65 00 49 07.16 +20 49 19.3 .30- .37+ 0.1325 1.1275
09 09 29.55246 F65 00 49 06.98 +20 49 22.9 .28+ .03+ 0.1325 1.1275
09 09 29.82471 J95 00 48 08.74 +21 09 43.9 .11+ .46+ 0.1338 1.1286
09 09 29.83048 J95 00 48 07.38 +21 10 10.1 .43+ .21- 0.1338 1.1287
09 09 29.84083 J95 00 48 04.82 +21 10 58.1 .56- .34- 0.1339 1.1287
--- Конец кода ---
Как видно дастаточно корректно с МРС.
Навигация
Перейти к полной версии