Measure aperture photometry using the mido.sh routine and other assorted scripts. An easy way to run this is to use apphot_interactive.
% apphot --help Usage: apphot 20191022T025524.2_acm_sci.fits aps1.cdfp 3.0 AP1 arg1 - FITS image to be measured arg2 - name of the cdfp file containing aperture position (Xcent,Ycent) arg3 - aperture radius in units of arcseconds arg4 - basename for the output table file Additional options: -v = print verbose comments --help = view a help document
Note that the cdfp file specifiying the aapertures can have any name (i.e. it is not tied to the basename of the input image). This was done so that one cdfp file can be used to compute magnitude sets using different aperture sizes. The output from this routine is a table file with the basename specified by arg4. The regions file for these apertures is stored using the table basename.
This routine was assembled from a lot of pre-existing routines in April 2020. I collected some working notes for this period.
Below I show a more verbose example of an apphot run.
% ls aps1.cdfp % apphot ../local_red/WCS/20191022T025524.2_acm_sci.fits aps1.cdfp 3.0 APER1 **** Please wait, mido.sh is computing images parameters. % ls 20191022T025524.2_acm_sci.fits APER1.parlab APER1.table APER1.cdfp APER1.reg aps1.cdfp The APER1.t table file contains: % cat APER1.parlab Xpix X pixels (aperture X center) Ypix Y pixels (aperture Y center) Xcen X center, pixels (intensity weighted centroid) Ycen Y center, pixels (intensity weighted centroid) Mag30 magnitude (assuming zp=30) MeanAp Average value per pixel in aperture MeanSky Average value per pixel in sky annulus Npix Number of pixels in aperture Fpeak Peak pixel value (no bkg-sub, image ADU) Mtype marker code (1=circle,2=box,3=ellipse) The Ra,Dec values are pasted into the table file to get a cdfp file for use in psoitional cross-matching in the, for instance, the zeropoint determination later: % head -5 APER1.cdfp Ra_hrs,Dec_deg,Xpix,Ypix,Xcen,Ycen,Mag30,MeanAp,MeanSky,Npix,Fpeak,Mtype, 19.740600586 26.221525192 737.84 209.03 737.74 209.10 12.322 31865.816 1263.165 385.0 64219.45 1 19.739345551 26.207641602 472.52 90.99 472.54 91.06 15.031 2632.667 115.620 386.0 13957.64 1 19.736711502 26.249803543 168.15 756.53 168.20 756.59 15.885 1236.209 81.096 383.0 5352.90 1 19.739686966 26.216396332 562.44 186.81 562.46 186.85 16.129 1002.018 78.761 383.0 5290.61 1