@@ -66,7 +66,7 @@ from astropy.stats import sigma_clipped_stats
6666# Photometry tools
6767import sep
6868from photutils.detection import DAOStarFinder
69- from photutils.psf import PSFPhotometry, IterativePSFPhotometry, IntegratedGaussianPRF , make_psf_model_image
69+ from photutils.psf import PSFPhotometry, IterativePSFPhotometry, CircularGaussianSigmaPRF , make_psf_model_image
7070from photutils.background import LocalBackground, MMMBackground
7171
7272# Firefly
@@ -337,7 +337,7 @@ psf_fwhm = 0.16 # PSF FWHM in arcsec
337337pixscale = 0.1 # VIS pixel scale in arcsec/px
338338
339339init_params = Table([objects[" x" ],objects[" y" ]] , names = [" x" ," y" ]) # initial positions
340- psf_model = IntegratedGaussianPRF (flux = 1 , sigma = psf_fwhm/ pixscale / 2.35 )
340+ psf_model = CircularGaussianSigmaPRF (flux = 1 , sigma = psf_fwhm/ pixscale / 2.35 )
341341psf_model.x_0.fixed = True
342342psf_model.y_0.fixed = True
343343psf_model.sigma.fixed = False
@@ -457,7 +457,7 @@ psf_fwhm = 0.3 # arcsec
457457pixscale = 0.3 # arcsec/px
458458
459459init_params = Table([xy[0 ],xy[1 ]] , names = [" x" ," y" ]) # initial positions
460- psf_model = IntegratedGaussianPRF (flux = 1 , sigma = psf_fwhm/ pixscale / 2.35 )
460+ psf_model = CircularGaussianSigmaPRF (flux = 1 , sigma = psf_fwhm/ pixscale / 2.35 )
461461psf_model.x_0.fixed = True
462462psf_model.y_0.fixed = True
463463psf_model.sigma.fixed = False
0 commit comments