Тематические рабочие группы > Вопросы астрометрии
Программа для проверки точности астрометрии на короткой дуге
ctac:
Программа aststat (360Кб)
Автор: Ирина Степановна Гусева (ГАО РАН, Пулково)
октябрь 2006г
Программа работает на основе стандартного формата астрометрических данных MPC
Пример1:
--- Код: --- K07L00F C2007 06 07.89636 14 16 43.33 +17 23 50.7 B16
K07L00F C2007 06 07.89830 14 16 45.40 +17 23 47.3 B16
K07L00F C2007 06 07.89852 14 16 45.64 +17 23 46.7 B16
K07L00F C2007 06 07.89896 14 16 46.12 +17 23 46.2 B16
K07L00F C2007 06 07.90026 14 16 47.54 +17 23 43.9 B16
K07L00F C2007 06 07.90396 14 16 51.45 +17 23 37.5 B16
K07L00F C2007 06 07.90504 14 16 52.57 +17 23 35.8 B16
K07L00F C2007 06 07.90570 14 16 53.30 +17 23 34.6 B16
K07L00F C2007 06 07.90635 14 16 53.97 +17 23 33.3 B16
K07L00F C2007 06 07.90700 14 16 54.67 +17 23 32.9 B16
K07L00F C2007 06 07.90983 14 16 57.69 +17 23 28.3 B16
K07L00F C2007 06 07.91352 14 17 01.65 +17 23 21.5 B16
K07L00F C2007 06 07.91439 14 17 02.57 +17 23 19.5 B16
K07L00F C2007 06 07.91482 14 17 03.04 +17 23 19.1 B16
K07L00F C2007 06 07.91678 14 17 05.11 +17 23 15.8 B16
K07L00F C2007 06 07.91938 14 17 07.85 +17 23 11.2 B16
K07L00F C2007 06 07.91982 14 17 08.35 +17 23 10.9 B16
K07L00F C2007 06 07.92025 14 17 08.81 +17 23 09.6 B16
K07L00F C2007 06 07.92310 14 17 11.88 +17 23 04.6 B16
K07L00F C2007 06 07.92332 14 17 12.16 +17 23 04.6 B16
K07L00F C2007 06 07.92352 14 17 12.18 +17 23 03.5 B16
K07L00F C2007 06 07.92396 14 17 12.78 +17 23 03.5 B16
K07L00F C2007 06 07.92462 14 17 13.44 +17 23 02.6 B16
K07L00F C2007 06 07.92505 14 17 14.03 +17 23 01.8 B16
K07L00F C2007 06 07.92744 14 17 16.46 +17 22 56.8 B16
K07L00F C2007 06 07.92940 14 17 18.68 +17 22 55.9 B16
K07L00F C2007 06 07.92962 14 17 18.89 +17 22 52.7 B16
K07L00F C2007 06 07.93005 14 17 19.25 +17 22 53.3 B16
K07L00F C2007 06 07.93093 14 17 20.15 +17 22 51.1 B16
K07L00F C2007 06 07.93180 14 17 21.16 +17 22 49.8 B16
K07L00F C2007 06 07.93202 14 17 21.36 +17 22 49.0 B16
K07L00F C2007 06 07.93333 14 17 22.77 +17 22 47.5 B16
K07L00F C2007 06 07.93442 14 17 23.93 +17 22 45.2 B16
K07L00F C2007 06 07.93550 14 17 25.10 +17 22 43.1 B16
K07L00F C2007 06 07.93707 14 17 26.69 +17 22 40.7 B16
K07L00F C2007 06 07.94207 14 17 32.09 +17 22 32.3 B16
K07L00F C2007 06 07.94425 14 17 34.46 +17 22 28.1 B16
--- Конец кода ---
Расшифровка по номеру столбца:
1-5 - пусто
6-12 - обозначение объекта съемки (программой не учитывается)
13-14 - пусто
15-19 - год
20 - пусто
21-22 - месяц
23 - пусто
24-25 - день
26-31-32 - доли суток(иногда применяют расширеную запись)
32 - пусто
33-34 - RA h
35 - пусто
36-37 - RA m
38 - пусто
39-43-44 - RA sec (иногда применяют расширеную запись)
44 - пусто
45 - знак DEC
46-47 - DEC d
48 - пусто
49-50 - DEC '
51 - пусто
52-55-56 - DEC " (иногда применяют расширеную запись)
56-77 - пусто (хотя там может встречаться и оценки блеска, но я их опускаю за ненадобностью) (программой не учитывается)
78-80 - код обсерватории по списку MPC (программой не учитывается)
А вот что дает в результате программа aststat.
Пример2:
--- Код: --- 2007LF.txt
Number of observations 37
UTC,hms RA,hms dRA,s DEC,g`" dDEC,", delta,"
7.89636 14 16 43.330 -0.004 17 23 50.70 -0.05 0.07
7.89830 14 16 45.400 -0.003 17 23 47.30 -0.12 0.13
7.89852 14 16 45.640 0.002 17 23 46.70 -0.34 0.34
7.89896 14 16 46.120 0.013 17 23 46.20 -0.08 0.20
7.90026 14 16 47.540 0.046 17 23 43.90 -0.15 0.68
7.90396 14 16 51.450 0.009 17 23 37.50 -0.19 0.23
7.90504 14 16 52.570 -0.023 17 23 35.80 -0.04 0.33
7.90570 14 16 53.300 0.002 17 23 34.60 -0.10 0.11
7.90635 14 16 53.970 -0.021 17 23 33.30 -0.29 0.41
7.90700 14 16 54.670 -0.014 17 23 32.90 0.43 0.47
7.90983 14 16 57.690 -0.013 17 23 28.30 0.69 0.72
7.91352 14 17 1.650 0.011 17 23 21.50 0.23 0.28
7.91439 14 17 2.570 0.002 17 23 19.50 -0.27 0.27
7.91482 14 17 3.040 0.014 17 23 19.10 0.07 0.21
7.91678 14 17 5.110 -0.007 17 23 15.80 0.13 0.17
7.91938 14 17 7.850 -0.041 17 23 11.20 0.00 0.59
7.91982 14 17 8.350 -0.010 17 23 10.90 0.46 0.48
7.92025 14 17 8.810 -0.009 17 23 9.60 -0.10 0.16
7.92310 14 17 11.880 0.021 17 23 4.60 -0.21 0.36
7.92332 14 17 12.160 0.066 17 23 4.60 0.17 0.97
7.92352 14 17 12.180 -0.127 17 23 3.50 -0.59 1.91
7.92396 14 17 12.780 0.003 17 23 3.50 0.17 0.18
7.92462 14 17 13.440 -0.041 17 23 2.60 0.41 0.71
7.92505 14 17 14.030 0.091 17 23 1.80 0.34 1.35
7.92744 14 17 16.460 -0.029 17 22 56.80 -0.55 0.69
7.92940 14 17 18.680 0.100 17 22 55.90 1.92 2.39
7.92962 14 17 18.890 0.076 17 22 52.70 -0.91 1.41
7.93005 14 17 19.250 -0.023 17 22 53.30 0.43 0.55
7.93093 14 17 20.150 -0.062 17 22 51.10 -0.25 0.92
7.93180 14 17 21.160 0.020 17 22 49.80 -0.06 0.29
7.93202 14 17 21.360 -0.015 17 22 49.00 -0.48 0.53
7.93333 14 17 22.770 -0.002 17 22 47.50 0.27 0.27
7.93442 14 17 23.930 -0.005 17 22 45.20 -0.16 0.17
7.93550 14 17 25.100 0.013 17 22 43.10 -0.40 0.44
7.93707 14 17 26.690 -0.072 17 22 40.70 -0.11 1.03
7.94207 14 17 32.090 -0.006 17 22 32.30 0.09 0.12
7.94425 14 17 34.460 0.039 17 22 28.10 -0.37 0.66
------ ------
Sigma, arcsec 0.61 0.46
--- Конец кода ---
Так что делает это ПО ?
Она берет все те точки, астрометрии, что вы ей даете и строит среднюю линию по ним (полином, если мне не изменяет память ::) ) и затем оценивает разброс точек относительно этой средней точки. Убирая выбросы (смотрите колонку delta," ), что привышают определенную планку (для этих измерений я взял планку в 0.5") можно улучшить общий результат.
Например вот до такого вида:
Пример3:
--- Код: --- K07L00F C2007 06 07.89636 14 16 43.33 +17 23 50.7 B16
K07L00F C2007 06 07.89830 14 16 45.40 +17 23 47.3 B16
K07L00F C2007 06 07.89852 14 16 45.64 +17 23 46.7 B16
K07L00F C2007 06 07.89896 14 16 46.12 +17 23 46.2 B16
K07L00F C2007 06 07.90396 14 16 51.45 +17 23 37.5 B16
K07L00F C2007 06 07.90504 14 16 52.57 +17 23 35.8 B16
K07L00F C2007 06 07.90570 14 16 53.30 +17 23 34.6 B16
K07L00F C2007 06 07.90635 14 16 53.97 +17 23 33.3 B16
K07L00F C2007 06 07.91352 14 17 01.65 +17 23 21.5 B16
K07L00F C2007 06 07.91439 14 17 02.57 +17 23 19.5 B16
K07L00F C2007 06 07.91482 14 17 03.04 +17 23 19.1 B16
K07L00F C2007 06 07.91678 14 17 05.11 +17 23 15.8 B16
K07L00F C2007 06 07.92025 14 17 08.81 +17 23 09.6 B16
K07L00F C2007 06 07.92310 14 17 11.88 +17 23 04.6 B16
K07L00F C2007 06 07.92396 14 17 12.78 +17 23 03.5 B16
K07L00F C2007 06 07.93180 14 17 21.16 +17 22 49.8 B16
K07L00F C2007 06 07.93333 14 17 22.77 +17 22 47.5 B16
K07L00F C2007 06 07.93442 14 17 23.93 +17 22 45.2 B16
K07L00F C2007 06 07.93550 14 17 25.10 +17 22 43.1 B16
K07L00F C2007 06 07.94207 14 17 32.09 +17 22 32.3 B16
--- Конец кода ---
И!
Пример4:
--- Код: --- 2007LF-3.txt
Number of observations 20
UTC,hms RA,hms dRA,s DEC,g`" dDEC,", delta,"
7.89636 14 16 43.330 -0.002 17 23 50.70 0.08 0.08
7.89830 14 16 45.400 -0.002 17 23 47.30 0.01 0.03
7.89852 14 16 45.640 0.003 17 23 46.70 -0.22 0.22
7.89896 14 16 46.120 0.014 17 23 46.20 0.04 0.20
7.90396 14 16 51.450 0.009 17 23 37.50 -0.08 0.15
7.90504 14 16 52.570 -0.023 17 23 35.80 0.07 0.33
7.90570 14 16 53.300 0.003 17 23 34.60 0.00 0.04
7.90635 14 16 53.970 -0.021 17 23 33.30 -0.19 0.35
7.91352 14 17 1.650 0.010 17 23 21.50 0.31 0.34
7.91439 14 17 2.570 0.001 17 23 19.50 -0.20 0.20
7.91482 14 17 3.040 0.012 17 23 19.10 0.14 0.23
7.91678 14 17 5.110 -0.009 17 23 15.80 0.20 0.24
7.92025 14 17 8.810 -0.011 17 23 9.60 -0.04 0.16
7.92310 14 17 11.880 0.018 17 23 4.60 -0.15 0.30
7.92396 14 17 12.780 0.001 17 23 3.50 0.22 0.22
7.93180 14 17 21.160 0.016 17 22 49.80 -0.03 0.23
7.93333 14 17 22.770 -0.006 17 22 47.50 0.29 0.30
7.93442 14 17 23.930 -0.009 17 22 45.20 -0.14 0.19
7.93550 14 17 25.100 0.008 17 22 43.10 -0.39 0.40
7.94207 14 17 32.090 -0.011 17 22 32.30 0.08 0.18
------ ------
Sigma, arcsec 0.17 0.18
--- Конец кода ---
Точность измерений в среднем увеличилась в 3.6 раз по RA и в 2.5 по DEC! Не плохо. :)
Как работать с программой aststat? Просто и легко.
1) Нужно в ту же папку, где лежит у вас aststat.exe закачать условный файл скажем asteroid.txt (важно, что бы не было лишних пробелов, отступов и сторочек! по сравнению со стандартным форматом);
2) Запустить ПО aststat и написать имя файла, что выжелаете обработать: asteroid.txt , нажать Enter и дать имя новому файлу, в который будут сохранены результаты оценки точности. Скажем rez-asteroid.txt. Нажать Enter.
3) Программа закроется и через пару секунд создаст файл rez-asteroid.txt.
4) Далее удаляем все выбросы, из файла asteroid.txt и сохраняем его под новым именем, например asteroid-2.txt .
5) Далее заново прогоняем шаги 1-3 (изменяя имена естественно) и смотрим насколько нас устраивает результат. Улучшать результат можно долго. Только смотрите, что бы не осталось 3-4 точки в результате ::) Думаю, что для нормальной оценки нужно не менее 10 точек.
Vit:
А насколько можно доверять этой программе?
Я тут загнал в нее наблюдения кометы 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:
--- Цитата: Vit от 05 июля 2007, 23:47:18 ---А насколько можно доверять этой программе?
Честно говоря, Астрометрика дает для этих наблюдений максимальную погрешность, к примеру, по Decl. 0.35”. Сравнивая свои наблюдательные данные по астероидам с Насавским калькулятором эфемерид, я не отмечал даже на самых плохих снимках расхождений больше чем в 0.5”, а для выборочных это значение укладывалось в 0.2”. Но чтобы ошибка достигала 1”, такого не было. Мне кажется, что выдаваемые данные в этой программе сильно зависят от количества наблюдений и (или) длины дуги, неплохо бы оценить минимальные значения этих параметров.
--- Конец цитаты ---
Я ей доверяю на 100%.
Виталий, вышли сюда, пожалуйста, твои наблюдения астероидов, что ты отправил на регистрацию в MPC - для них точность пересчитают в Пулково (в MPC в таком формате не дадут показатели точности). А потом можешь пронаблюдать любой NEO и увидеть на сайте http://unicorn.eis.uva.es/cgi-bin/neodys/neoibo реальные погрешности твоих наблюдений с учетом всех наблюдений.
andyp:
А исходник у этой программы есть ?
Хотелось бы использовать его в OrbSearch.
ctac:
--- Цитата: andyp от 23 июля 2007, 14:03:55 ---А исходник у этой программы есть ?
Хотелось бы использовать его в OrbSearch.
--- Конец цитаты ---
Все что у меня есть - я выложил. Думаю, что там ни чего сложного нет - просто полином строится на коротком отрезке дуги, где ее можно принять за прямую, а объект движется без ускорений.
Навигация
Перейти к полной версии