The wcsf_ series of scripts drive the processes I use to derive and install WCS header information in FITS images. Here I briefly review the 4 basic scripts and record a lot ofthe steps I use to process 3 images.
# I set up the sample reduction
Where I did this work: /home/sco/A_wcsf/run3
% cat ./S/list.0
/home/sco/A_wcsf/images/20180206T020222.8_acm_sci_proc.fits
/home/sco/A_wcsf/images/20180206T022042.6_acm_sci_proc.fits
/home/sco/A_wcsf/images/20180206T032156.5_acm_sci_proc.fits
% bigds9 ./S/list.0 4 4
**** Things look reasonable.
# I install rough WCS
% wcsf_rough ./S/list.0 N
Usage: wcsf_rough image.list Y
arg1 - Name of file with list of FITS images (can be full path)
arg2 - run in debug mode (Y/N)
**** For 3 images this took about 5min (from a cold start, not having run this in a few days).
% pwd
/home/sco/A_wcsf/run3
% ls
20180206T020222.8_acm_sci_proc.fits big.NOMARK fits_review.pass0 m2.crpix Store_DSS/
20180206T022042.6_acm_sci_proc.fits ds9_open.Size fits_review.pass0_Explain m3.crpix
20180206T032156.5_acm_sci_proc.fits fits_review.params_used m1.crpix S/
# I test the WCS (and make USNO catalog files)
% wcsf_viscat ./S/list.0 N
Would you also like to run in verbose mode? (Y/N): N
Parameter values set with params_4_fits_review calls:
Vmag,radius,circrad (for wcs_refine) = 18.0 3.00 4.0
magname,magmax (for usno_look_wcs) = Red 20.0
Skip the use of wcs_refine? (Y/N): Y
% ls
20180206T020222.8_acm_sci_proc.fits 20180206T022042.6_acm_sci_proc-usno.cdfp ds9_open.Size m2.crpix
20180206T020222.8_acm_sci_proc.png 20180206T032156.5_acm_sci_proc.fits fits_review.params_used m3.crpix
20180206T020222.8_acm_sci_proc-usno.cdfp 20180206T032156.5_acm_sci_proc.png fits_review.pass0 S/
20180206T022042.6_acm_sci_proc.fits 20180206T032156.5_acm_sci_proc-usno.cdfp fits_review.pass0_Explain Store_DSS/
20180206T022042.6_acm_sci_proc.png big.NOMARK m1.crpix
**** This took about 20 sec! The png files are not so important, but the *-usno.cdfp will be
used later.
# I compile image catalogs
**** I change the adress of the nput images - I want to use the
images that have rough WCS installed.
% cat S/list.1
/home/sco/A_wcsf/run3/20180206T020222.8_acm_sci_proc.fits
/home/sco/A_wcsf/run3/20180206T022042.6_acm_sci_proc.fits
/home/sco/A_wcsf/run3/20180206T032156.5_acm_sci_proc.fits
% wcsf_imgcat ./S/list.1 N
Would you also like to run in verbose mode? (Y/N): N
% ls
big.NOMARK fits_review.defaults fits_review.pass0 local_red/ m2.crpix S/
ds9_open.Size fits_review.params_used fits_review.pass0_Explain m1.crpix m3.crpix Store_DSS/
# Notice how thinsg cleaned up! Where did the files go?
% ls local_red/
IMGCAT0/
% ls local_red/IMGCAT0/
20180206T020222.8_acm_sci_proc_bkg1.fits 20180206T022042.6_acm_sci_proc_bkg1.fits 20180206T032156.5_acm_sci_proc_bkg1.fits
20180206T020222.8_acm_sci_proc.cat0 20180206T022042.6_acm_sci_proc.cat0 20180206T032156.5_acm_sci_proc.cat0
20180206T020222.8_acm_sci_proc.cdfp 20180206T022042.6_acm_sci_proc.cdfp 20180206T032156.5_acm_sci_proc.cdfp
20180206T020222.8_acm_sci_proc_detsig.fits 20180206T022042.6_acm_sci_proc_detsig.fits 20180206T032156.5_acm_sci_proc_detsig.fits
20180206T020222.8_acm_sci_proc.fits 20180206T022042.6_acm_sci_proc.fits 20180206T032156.5_acm_sci_proc.fits
20180206T020222.8_acm_sci_proc_label0.fits 20180206T022042.6_acm_sci_proc_label0.fits 20180206T032156.5_acm_sci_proc_label0.fits
20180206T020222.8_acm_sci_proc.png 20180206T022042.6_acm_sci_proc.png 20180206T032156.5_acm_sci_proc.png
20180206T020222.8_acm_sci_proc_sigma.fits 20180206T022042.6_acm_sci_proc_sigma.fits 20180206T032156.5_acm_sci_proc_sigma.fits
20180206T020222.8_acm_sci_proc-usno.cdfp 20180206T022042.6_acm_sci_proc-usno.cdfp 20180206T032156.5_acm_sci_proc-usno.cdfp
*** This took about 20 seconds!
# Compute the final WCS headers
% wcsf_final ./S/list.1 N
Would you also like to run in verbose mode? (Y/N): N
*** This took about a minute, and mostly just answered "Y" and "2lin".
% ls
big.NOMARK fits_review.defaults fits_review.pass0 local_red/ m2.crpix S/
ds9_open.Size fits_review.params_used fits_review.pass0_Explain m1.crpix m3.crpix Store_DSS/
% ls local_red/
IMGCAT0/ WCS/
% ls local_red/WCS
20180206T020222.8_acm_sci_proc.fits 20180206T022042.6_acm_sci_proc.fits 20180206T032156.5_acm_sci_proc.fits
I could compile a nice summary list:
% gethead ./local_red/WCS/*fits RMSRA RMSDEC NUMWCS FILTER EXPTIME
20180206T020222.8_acm_sci_proc.fits 0.7800 0.8660 51 B 5
20180206T022042.6_acm_sci_proc.fits 0.3950 0.4000 22 B 5
20180206T032156.5_acm_sci_proc.fits 0.2650 0.4130 7 B 5
April 17 reduction of 14 images
===============================
% ls -1 /home/sco/A_wcsf/images/20*fits > list.1 # stored list.1 in ./S
% bigds9 ./S/list.1 4 4 # 10:55 am
% wcsf_rough ./S/list.1 N # took 10:55am to 11:06am
% wcsf_viscat ./S/list.1 N # took 11:08am to 11:10am
% ls -1 /home/sco/A_wcsf/Run_apr17/20*fits >list.2
% cp list.2 S
% wcsf_imgcat ./S/list.2 N # Use defaults, 11:14am to 11:15am
*** All images and producst are in ./local_red/IMGCAT0/
*** I can use the original list.1 since I am just using list.1 to get image
basenames. The input images will come from .//local_red/IMGCAT0/
% wcsf_final ./S/list.1 N # start 11:21am to 11:28 am
At this point after about 30min I have 14 WCS-calibrated images in
/home/sco/A_wcsf/Run_apr17/local_red/WCS
I could compile a nice summary list:
% gethead ./local_red/WCS/*fits RMSRA RMSDEC NUMWCS FILTER EXPTIME
20180206T020222.8_acm_sci_proc.fits 0.9050 0.7340 50 B 5
20180206T022042.6_acm_sci_proc.fits 0.5410 0.4040 23 B 5
20180206T032156.5_acm_sci_proc.fits 0.7420 0.4010 8 B 5
20180206T035223.8_acm_sci_proc.fits 0.7380 0.4380 8 B 5
20180206T061218.4_acm_sci_proc.fits 0.3490 0.5040 17 r` 5
20180206T070902.2_acm_sci_proc.fits 0.4880 0.5910 14 r` 5
20180206T072507.6_acm_sci_proc.fits 0.7150 1.3590 6 r` 5
20180206T083937.3_acm_sci_proc.fits 0.2200 0.3490 11 r` 5
20180206T092518.0_acm_sci_proc.fits 0.6330 0.5700 7 r` 5
20180206T100533.3_acm_sci_proc.fits 0.3580 0.2580 9 r` 5
20180206T101258.1_acm_sci_proc.fits 0.2650 0.1770 10 r` 5
20180206T102445.6_acm_sci_proc.fits 0.4900 0.7590 14 r` 5
20180206T104913.4_acm_sci_proc.fits 0.2050 0.2630 10 r` 5
20180206T105144.7_acm_sci_proc.fits 0.3480 0.2780 11 r` 5
So, the first 4 images are B filter, the next 10 are Gunn r