照片DAOPhot不适合星星吗？

Question

我最近遇到了PhotUtils包，并试图使用它对我的一些图像进行PSF测光。然而，当我试图运行代码时，我会得到非常奇怪的结果。当我绘制由get_residual_image()生成的图像时，星体没有被很好地移除。下面显示了一些示例图像。

第一个图像的西格玛设置为2.05，正如PhotUtils文档中的一个示例程序所示：

然而，恒星似乎只在它们的中心被移除。

第二图像将西格玛设置为5.0。这个特别奇怪。有些恒星被移除，有些被移除，一些黑色方块被添加到图像中，等等。

​

这是我的代码：

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中绘制原始和最后的图像。我用的是灰阶彩色地图。左图像稍微变暗的原因是它们使用了不同的颜色缩放。

也许我把其中一个参数设置错了？

有人能帮我解决这个问题吗？谢谢!

Answer 1

看着残差图像，立即告诉我，背景减法可能是错误的。我可以复制结果，并想知道，如果MMMBackground没有正确地完成这项工作。在仔细研究了文档之后，用照片来获得明星终于给出了基本的提示：

image -= np.median(image)