Итак, долго ли сказка сказывается, да не скоро дело делается. Но таки разрешил я проблему с определением центра вращения, которая длилась на продолжении раз, два, трех публикаций. Решение оказалось простым донельзя.

А секретов-то и нет! )

Итак, в предыдущих темах я привел пример алгоритма для определения положения центра вращения изображений в предположении, что само изображение не смещается; далее я привел соображения по поводу вычислений с использованием матричной алгебры.

Здесь я опишу реализацию алгоритмического метода и сравню его с методом определения центра вращения по пересечениям радиусов окружности.

Поехали! )

В прошлой записи я пытался при помощи Octave найти центр вращения поля, исходя из координат объектов на снимках, полученных последовательно с некоторым интервалом времени. Тот метод был непригоден в случае, если поле не только вращалось, но и смещалось. Здесь я изложу некоторые соображения, которые позволят (возможно) найти нужные параметры.

Немного алгебры )

Ну, а в следующей заметке я попытаюсь применить изложенное для автоматического определения центра вращения поля и его смещения от кадра к кадру.

Итак, применим Octave для обработки набора снимков звездного поля, сделанного для определения центра вращения поля (при остановленном компенсаторе вращения). Нужно это, чтобы определить, насколько расходится центр вращения поля с центром вращения компенсатора.

Для того, чтобы отождествить одни и те же звезды на соседних кадрах, я разработал простенький алгоритм. Основную массу кода будет занимать как раз он. Ну, а определение центра - задача простая. Для этого достаточно найти уравнения прямых, проходящих через центр изображений одной и той же звезды на соседних кадрах, перпендикулярно отрезку, соединяющему эти изображения; а затем найти точки пересечения всех этих прямых.

Поехали! )

P.S. Я не учитывал здесь то, что в течение наблюдений изображение также может испытывать линейные смещения. Об этом - в следующий раз.

September 2017

S M T W T F S
     1 2
3456789
1011 12 13141516
17181920 212223
24252627282930

Syndicate

RSS Atom

Style Credit

Expand Cut Tags

No cut tags
Page generated Sep. 25th, 2017 04:20 am
Powered by Dreamwidth Studios