Не очень понял. Там все сильно навороченно, сотни функций с мало чего говорящими названиями. Существует какая-нибудь интерактивная, общепонятная программа без наворотов?
Ведь насколько я понимаю, задача решается просто:
1. Определить fwhm снимка.
2. Выделить на нем все звездные объекты.
3. Определить среднюю яркость пиксела в областях, где никаких объектов не выделилось.
4. Опознать несколько звезд или взять уже опознанные.
5. Определить сумму интенсивности пикселов в пределах fwhm от центроида звезды, умножить на 2. Получится интенсивность, соответствующая ее звездной величине. Усреднить по всем опознанным звездам. То есть по сути, выведем "отдачу" снимка - какой звездной величине соответствует 1 единица интенсивности на снимке.
6. Предельная звездная величина вычисляется как равная звездной величине звезды, которая дает суммарную интенсивность в пределах fwhm, равную 3 (т.к. три сигмы) / 2 (так как fwhm это по определению радиус в который попадает половина всего света звезды) * площадь попадающую в радиус fwhm * средняя интенсивность пиксела без звезд, определенная в (3).
Я верно понимаю процесс? Есть какая-то прога его реализующая, в которой можно разобраться легко?