Repository navigation
PSF photometry on crowded CCD frames #433
Description
Activity
Hi @Gabriel-p!
2.:
I think you may find equivalent python functionalities in astropy/ccdproc.3.:
Photutils has an aperture photometry module, you may check the docs hereCurrently, I'm working on refactoring 5. in photutils. You may check here.
I'm not sure about 1., 4., 6., and 7..
I would say that you should open a PR if you have implemented/improved any of these functionalities you mentioned :)
No need to be sorry.
Reacted by Gabriel Perren@Gabriel-p I would say that
photutilsdoes not have a direct way to estimate the FWHM of PSFs. It would be a great addition tophotutilsto have apsfmeasureequivalent :D
@bsipocz @eteq @hamogu @larrybradley (pinging the jedi masters)Reacted by Gabriel PerrenYes, this would be a nice feature to add to
photutils. I'm pinging @SaraOgaz because I think she was working on a python version ofpsfmeasure.Reacted by Gabriel PerrenHi @Gabriel-p , and all... we have just started an effort at STScI to get down to the nitty gritty of replacing IRAF tasks still used by the community. So first of all, this input is MUCH appreciated. We would like to either have a notebook with python examples for IRAF tasks if they already exist in the python libraries, or if they are still missing, write them. But we're still in the very early stages of this. Long term, we would like to write a psfmeasure replacement if we can't find an existing one.
I will definitely be adding the stellar photometry users guide to our list of examples. If we end up having the time and man-power, we'd like to replace a lot of these guides with the python equivalent walkthrough.
Sorry I don't have a quick solution for you, but just know that we're working on it!
Reacted by Gabriel PerrenThank you all. I'll see if I can come up with some (dirty) script to serve as the base for a proper
psfmeasuretask.I guess this issue can be closed now?
Cheers
@Gabriel-p - a bit more on these steps:
2: can be done with
ccdprocor just directly on image arrays, or for dithered images, you might want to have a look at reproject (http://reproject.readthedocs.io/), which has various ways of putting them all on the same frame without combing.One extra bit here, though: in fact in many cases you don't want to do this at all - the mode we're hoping to implement is instead providing the different exposures separately, but doing a simultaneous fit, rather than coadding and then fitting on the coadds. But with the above tools, essentially both are supported.
3: I think in @mirca's last implementation this is done simultaneously with 5, although you can do it separately if you want.
4, 6, 7: The intent is to provide utilities to do all this. The basic idea would be functions that take the tables that the PSF photometry classes from 5 output and add extra columns with instrumental mags, calibration mags, etc. This is an area where contributions would be especially useful because it's mostly just table operations. (Will try to make an issue for this).
I'll close this now because I think the questions are answered, but will try to make another issue to address the 4/6/7 question above.
Reacted by Gabriel Perren
Hi there. I occasionally need to perform PSF photometry on crowded fields, and I'm looking to get rid of IRAF as soon as possible. This package seems to be ideal replacement, but for what I can gather
photutilsit is still in a heavy development phase.Below are the steps needed to complete a full PSF photometry process, and the packages in IRAF used in each step (as described in A User's Guide to Stellar CCD Photometry with IRAF). I was wondering if these features are in the roadmap to be implemented, or if they are not planned at all.
imalign)imcombine)daophot --> daofind, phot)photcal --> mkcatalog, mknobsfile, mkconfig, fitparams)daophot --> daofind, phot, psf, nstar, substar, allstar)daophot --> mkobsfile, invertfit)Sorry for the long post.