Напишите о результатах.
В методе есть несколько ограничений. Универсальным его назвать вряд ли можно. Но для небольших полей, до 10-15 минут, и не слишком плохой дискретизации, FWHM>2.5-3 пикселей, точность хорошая.
Для относительно ярких звезд (но без перекопа) у нас получается ошибка ~0.005mag. Впрочем тут многое зависит от фокуса/пикселя/камеры.
Можно расширить на большие поля, но тут лучше использовать не Python, а более быстрые решения, вроде IRAF+AstroKit.