Тематические рабочие группы > Вопросы астрометрии

Программа для проверки точности астрометрии на короткой дуге

(1/3) > >>

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.

--- Конец цитаты ---

 Все что у меня есть - я выложил. Думаю, что там ни чего сложного нет - просто полином строится на коротком отрезке дуги, где ее можно принять за прямую, а объект движется без ускорений.

Навигация

[0] Главная страница сообщений

[#] Следующая страница

Перейти к полной версии