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

[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]

А насколько можно доверять этой программе?

Честно говоря, Астрометрика дает для этих наблюдений максимальную погрешность, к примеру, по Decl. 0.35”. Сравнивая свои наблюдательные данные по астероидам с Насавским калькулятором эфемерид, я не отмечал даже на самых плохих снимках расхождений больше чем в 0.5”, а для выборочных это значение укладывалось в 0.2”. Но чтобы ошибка достигала 1”, такого не было. Мне кажется, что выдаваемые данные в этой программе сильно зависят от количества наблюдений и (или) длины дуги, неплохо бы оценить минимальные значения этих параметров.

 Я ей доверяю на 100%.

 Виталий, вышли сюда, пожалуйста, твои наблюдения астероидов, что ты отправил на регистрацию в MPC - для них точность пересчитают в Пулково (в MPC в таком формате не дадут показатели точности). А потом можешь пронаблюдать любой NEO и увидеть на сайте http://unicorn.eis.uva.es/cgi-bin/neodys/neoibo реальные погрешности твоих наблюдений с учетом всех наблюдений.

[andyp]

А исходник у этой программы есть ?
Хотелось бы использовать его в  OrbSearch.

[ctac]

А исходник у этой программы есть ?
Хотелось бы использовать его в  OrbSearch.

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

[ctac]

Есть и интернет вариант: http://www.fitsblink.net/residuals/index.html

[ctac]

Загружаю сюда программу astrostat

[ctac]

А насколько можно доверять этой программе?

Я тут загнал в нее наблюдения кометы 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
Как видно дастаточно корректно с МРС.

[ched]

Я считаю, что подобная программа не очень-то нужна, так как за один "приём" незачем получать много позиций.


МРС рекомендует получать 2-3 позиции на объект за наблюдательную сессию, а не 10 и не 20. Если высылать в МРС помногу астрометрии, полученной для одного объекта за одну наблюдательную сессию, то вы ставите наблюдения других обсерваторий по данному объекту в неравное положение. У вас же могут быть и систематические ошибки!

[ctac]

Я считаю, что подобная программа не очень-то нужна, так как за один "приём" незачем получать много позиций.


МРС рекомендует получать 2-3 позиции на объект за наблюдательную сессию, а не 10 и не 20. Если высылать в МРС помногу астрометрии, полученной для одного объекта за одну наблюдательную сессию, то вы ставите наблюдения других обсерваторий по данному объекту в неравное положение. У вас же могут быть и систематические ошибки!

 Вот ситуация, когда нужно выбирать точки: снимаю астероид для построения кривой блеска: 300 точек за 6 часов наблюдений.
 Какие высылать в MPC? Крайние и среднюю - в слепую?
 Такая программа нужна

[ctac]

Теперь из программы FIND_ORB:

Стас, сравнение не совсем корректно, т.к. ты взял из МРС резедуалы со всей дуги, а в Findorb загнал только часть.
Вот что получается если использовать всю дугу:

Как видно дастаточно корректно с МРС.

 Да, Виталий - согласен. Но тут я сейчас проверял именно по короткой дуге

[ctac]

Вот для примера, что прислал Дима Честнов только что:
"Евпатория послала 518 позиций 2009 AP10 с 15 по 23 сентября
пятьсот восемнадцать!!!!!!! *WALL*"

[ched]

Программа может быть полезна для проверки астрометрии нового объекта перед отправкой её в MPC. Если есть хоть одна позиция, отклоняющаяся от общего курса более чем на 1,5", то робот MPC вообще не примет астрометрию данного объекта.