Derive the half-light radius, r50, and other quartile parameters of sources marked in the regions file associated with an input FITS image.
% find_r50 --help
usage: find_r50 Rsco2043.fits [-v] [-h] [-r1 Rin] [-r2 Rout]
Additional options:
-v = print verbose comments and run in debug mode
-h = show a help file and then exit
-r1 = inner bkg annulus radius in pixel units (can be "none")
-r2 = outer bkg annulus radius in pixel units (can be "none")
In many cases, I use an image source catalog such as that derived with
image_catalogs to define
apertures for bright sources and then derive the r50. Here is an exampe
where we use a source catalog:
% ls local_red/ S/ % imgcat0_region_files /home/sco/AUTO_match_for_wcs/S/20191018T024545.3_acm_sci.fits -n 8 -far 16 % ls 20191018T024545.3_acm_sci_apertures.cdfp 20191018T024545.3_acm_sci.cdfp 20191018T024545.3_acm_sci.fits S/ 20191018T024545.3_acm_sci_apertures.cdfp.reg 20191018T024545.3_acm_sci.cdfp.reg local_red/ % cp 20191018T024545.3_acm_sci_apertures.cdfp.reg 20191018T024545.3_acm_sci.reg % find_r50 20191018T024545.3_acm_sci.fits R50 = 4.236500 -+ 0.097633 (pixels) N=26 % ls 20191018T024545.3_acm_sci_apertures.cdfp 20191018T024545.3_acm_sci.fits local_red/ profiles.table 20191018T024545.3_acm_sci_apertures.cdfp.reg 20191018T024545.3_acm_sci.info profiles_full.parlab r25.stats 20191018T024545.3_acm_sci_bkgannuli.reg 20191018T024545.3_acm_sci.profpars profiles_full.table r50.stats 20191018T024545.3_acm_sci.cdfp 20191018T024545.3_acm_sci.reg profiles.params r75.stats 20191018T024545.3_acm_sci.cdfp.reg ERROR.find_r50 profiles.parlab S/Notice that a lot of peripheral files are built, but some particularly useful ones are: