Photutils DAOPhot 不适合星星?
Photutils DAOPhot Not Fitting stars well?
我最近 运行 研究了 PhotUtils 包,并尝试使用它对我拥有的一些图像执行 PSF 光度测量。但是,当我尝试 运行 代码时,得到的结果非常难运行。当我绘制由 get_residual_image() 生成的图像时,星星没有很好地去除。下面显示了一些示例图像。
第一个图像的 sigma 设置为 2.05,因为它在 PhotUtils 文档中的一个示例程序中:
但是,星星似乎只在它们的中心被移除。
第二张图片的 sigma 设置为 5.0。这个特别st运行ge。有些星星被过度移除,有些被移除,一些黑色方块被添加到图像中,等等。
这是我的代码:
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 中绘制了原始图像和最终图像。我使用灰度色图。左边的图像看起来稍微暗一些的原因是它们使用了不同的颜色比例。
也许我有一个参数设置不正确?
有人可以帮我解决这个问题吗?谢谢!
看残像瞬间告诉我背景减法可能是错误的。我可以重现结果并想知道 MMMBackground
是否没有正确完成工作。
仔细看了文档,Getting startet with Photutils终于给出了本质提示:
image -= np.median(image)
我最近 运行 研究了 PhotUtils 包,并尝试使用它对我拥有的一些图像执行 PSF 光度测量。但是,当我尝试 运行 代码时,得到的结果非常难运行。当我绘制由 get_residual_image() 生成的图像时,星星没有很好地去除。下面显示了一些示例图像。
第一个图像的 sigma 设置为 2.05,因为它在 PhotUtils 文档中的一个示例程序中:
但是,星星似乎只在它们的中心被移除。
第二张图片的 sigma 设置为 5.0。这个特别st运行ge。有些星星被过度移除,有些被移除,一些黑色方块被添加到图像中,等等。
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 中绘制了原始图像和最终图像。我使用灰度色图。左边的图像看起来稍微暗一些的原因是它们使用了不同的颜色比例。
也许我有一个参数设置不正确?
有人可以帮我解决这个问题吗?谢谢!
看残像瞬间告诉我背景减法可能是错误的。我可以重现结果并想知道 MMMBackground
是否没有正确完成工作。
仔细看了文档,Getting startet with Photutils终于给出了本质提示:
image -= np.median(image)