Complute and collect WCS data from ACAM headers processed with wcs_doall.

 

% acam_wcs_stats.sh
Usage: acam_wcs_stats.sh 20160302T034710.0_acm_sci.fits 
arg1 - WCS-calibrated FIT image name

Suppose you have WCS-calibrated a set of ACAM images and you want to see if the image properties are consistant. Usually, the first thing I do (see below) is run the gethead routine to check the RMSRA and RMSDEC values in the header to insurethat I have processed all of the images with wcs_doall (i.e. the values are NOT 1.0,1.0 exactly). Then, because acam_wcs_stats presently only runs on single images, I build a RUN script (see RUN_STATS below). This creates the desired output table I need for plotting (see acam_plots) and other analysis.

* Be sure you have solutions:
% gethead RMSRA RMSDEC 20*fits
20160301T023620.4_acm_sci.fits 0.2410 0.1170
20160301T030347.7_acm_sci.fits 0.0020 0.0060
20160301T031055.4_acm_sci.fits 0.2770 0.1960
20160301T031747.2_acm_sci.fits 0.4370 0.2350
20160301T032232.1_acm_sci.fits 0.0670 0.1990
20160301T032733.5_acm_sci.fits 0.1970 0.2940
20160301T033146.9_acm_sci.fits 0.2280 0.1820
20160301T034459.0_acm_sci.fits 0.1460 0.3100
20160301T035356.8_acm_sci.fits 0.1980 0.4610
20160301T040708.6_acm_sci.fits 0.4860 0.1840
20160301T041917.6_acm_sci.fits 0.2530 0.0770
20160301T043419.8_acm_sci.fits 0.0020 0.0060
20160301T044319.0_acm_sci.fits 0.2800 0.2570
20160301T045456.8_acm_sci.fits 0.4170 0.3530

* Build an input file 
% ls -1 20*fits > RUN_STATS 

**** I edit the RUN_STATS file to make is a script

% chmod 777 RUN_STATS 
% cat RUN_STATS 
#
acam_wcs_stats.sh 20160301T023620.4_acm_sci.fits
acam_wcs_stats.sh 20160301T030347.7_acm_sci.fits
acam_wcs_stats.sh 20160301T031055.4_acm_sci.fits
acam_wcs_stats.sh 20160301T031747.2_acm_sci.fits
acam_wcs_stats.sh 20160301T032232.1_acm_sci.fits
acam_wcs_stats.sh 20160301T032733.5_acm_sci.fits
acam_wcs_stats.sh 20160301T033146.9_acm_sci.fits
acam_wcs_stats.sh 20160301T034459.0_acm_sci.fits
acam_wcs_stats.sh 20160301T035356.8_acm_sci.fits
acam_wcs_stats.sh 20160301T040708.6_acm_sci.fits
acam_wcs_stats.sh 20160301T041917.6_acm_sci.fits
acam_wcs_stats.sh 20160301T043419.8_acm_sci.fits
acam_wcs_stats.sh 20160301T044319.0_acm_sci.fits
acam_wcs_stats.sh 20160301T045456.8_acm_sci.fits


% RUN_STATS
     0.268    88.619    0.26947     180.000    20160301T023620.4_acm_sci.fits          
     0.006    88.142    0.27072      21.452    20160301T030347.7_acm_sci.fits          
     0.339    87.917    0.27156      21.452    20160301T031055.4_acm_sci.fits          
     0.496    88.028    0.27109      21.452    20160301T031747.2_acm_sci.fits          
     0.210    88.073    0.27105      21.452    20160301T032232.1_acm_sci.fits          
     0.354    88.155    0.27009      21.452    20160301T032733.5_acm_sci.fits          
     0.292    88.279    0.27213      21.452    20160301T033146.9_acm_sci.fits          
     0.343    88.186    0.27051      66.856    20160301T034459.0_acm_sci.fits          
     0.502    88.274    0.27211      94.089    20160301T035356.8_acm_sci.fits          
     0.520    88.609    0.27252     115.276    20160301T040708.6_acm_sci.fits          
     0.264    88.392    0.27085     137.787    20160301T041917.6_acm_sci.fits          
     0.006    89.052    0.26887     162.610    20160301T043419.8_acm_sci.fits          
     0.380    86.651    0.27047     189.374    20160301T044319.0_acm_sci.fits          
     0.546    88.101    0.27107     214.103    20160301T045456.8_acm_sci.fits      

% cat acam_wcs_stats.explain

col01 = RMS, combined Ra,Dec residual rms in arcsec
col02 = tz, tz = Acam_PosAngle - HETQ  (degrees)
col03 = plate scale (arcsec/pixel_) 
col04 = HETQ, HET parallactic angle 
col05 = image name  
      RMS       tz        ps          HETQ     Image