Point spread function photometry with astropy/photutils

Lightning talk at PyAstro 2017

Transcript

1. Point Spread Function Photometry with astropy/photutils: a Google Summer of Code project Z´ e Vin´ ıcius mirca.github.io Lorentz Center, Leiden May 10, 2017
2. • Google Summer of Code 1
3. • Google Summer of Code • astropy/photutils 1
4. Google Summer of Code 2
5. Google Summer of Code • Project Idea: PSF photometry for fitting overlapping sources simultaneously 2
6. Google Summer of Code • Project Idea: PSF photometry for fitting overlapping sources simultaneously • Astropy Mentors 2
7. Google Summer of Code • Project Idea: PSF photometry for fitting overlapping sources simultaneously • Astropy Mentors • Student 2
8. DAOPHOT-like PSF Photometry of a Single Star • Data: Z = znm, n = 1, 2, ..., N, m = 1, 2, ..., M • Parameters: θ = (x0, y0, F) • Model: g(x, y, θ) • Objective: θ = arg min θ∈Θ n,m (znm − g(n, m, θ))2 3
9. DAOPHOT-like PSF Photometry of Multiple Overlapping Stars • Data: Z = znm, n = 1, 2, ..., N, m = 1, 2, ..., M • Parameters: θ(1) = x(1) 0 , y(1) 0 , F(1) θ(2) = x(2) 0 , y(2) 0 , F(2) ... θ(K) = x(K) 0 , y(K) 0 , F(K) • Model: g x, y, θ(k) • Objective: θ = arg min θ∈Θ n,m znm − k g n, m, θ(k) 2 4
10. DAOPHOT-like PSF Photometry stars = find(image) while len(stars) > 0: groups = group(stars) fitted_groups = fit(groups, model) residual = subtract(image, fitted_groups) stars = find(residual) basic_phot = BasicPSFPhotometry(...) itr_phot = IterativelySubtractedPSFPhotometry(...) daophot = DAOPhotPSFPhotometry(...) results = daophot(image) 5
11. DAOPHOT-like PSF Photometry 0 10 20 30 40 50 60 0 10 20 30 40 50 60 6
12. DAOPHOT-like PSF Photometry 0 10 20 30 40 50 60 0 10 20 30 40 50 60 7
13. astropy/photutils 8