Недавно я столкнулся с пакетом PhotUtils и пытаюсь использовать его для выполнения PSF-фотометрии на некоторых изображениях, которые у меня есть. Однако, когда я пытаюсь запустить код, я получаю очень странные результаты. Когда я строю изображение, сгенерированное get_residual_image (), звезды удаляются не очень хорошо. Некоторые примеры изображений показаны ниже.
На первом изображении сигма установлена на 2.05, как в одной из примеров программ в документации PhotUtils: 
Однако звезды появляются только должны быть удалены в их центре.
На втором изображении сигма установлена на 5,0. Этот особенно странный. Некоторые звезды слишком удалены, некоторые удалены, некоторые черные квадраты добавлены к изображению и т. Д. c.
Вот мой код:
import photutils
from photutils.psf import DAOPhotPSFPhotometry as DAOP
from photutils.psf import IntegratedGaussianPRF as PRF
from photutils.background import MMMBackground
bkg = MMMBackground()
background = 2.5*bkg(img)
gaussian_prf = PRF(sigma=5.0)
gaussian_prf.sigma.fixed = False
photTester = DAOP(8,background,5,gaussian_prf,31)
photResults = photTester(imgStars)
finalImg = photTester.get_residual_image()
После этого я просто строю исходное и окончательное изображение в MatPlotLib. Я использую цветовую карту в оттенках серого. Причина, по которой левые изображения выглядят немного темнее, заключается в том, что они используют другую цветовую шкалу.
Возможно, я неправильно установил один из параметров?
Может кто-нибудь помочь мне с этим? Спасибо!