| Science observation UT date(s) | 2010 Nov 01 |
| Data filename prefix(es) | S20101101S |
| Science | S20101101S | 87-90 | binning 2x2 |
| Processed Bias | for 20101101 | S20101101S0169_bias.fits | |
| Processed Twilight Flat | for 20101101 | S20101107S0201_flat_r.fits | |
| Processed Fringe Frame | N/A | N/A |
N/A
| Processed: S20101107S0201_flat_r.fits | |||
| Twilight Frames | S20101107S | 201-212 | binning 2x2 |
| Processed Bias | for 20101107 | S20101107S0228_bias.fits | |
| Processed: S20101101S0169_bias.fits | |||
| Biases | S20101101S | 169-172 (173 is bad) | binning 2x2 |
| S20101102S | 131-135 | binning 2x2 | |
| Processed: S20101107S0228_bias.fits | |||
| Biases | S20101107S | 228-231 (232 is bad) | binning 2x2 |
| S20101108S | 26-30 | binning 2x2 | |
The directory paths used in this notebook are assumed to match this structure:
rootdir/
calib/
bias/
flat/
fringe/
raw/
redux/The data processing will happen in the Notebook. It is also convenient to have a PyRAF interactive session ready to use if needed, eg. to imexam. In this step, we will configure the Notebook and the PyRAF session.
This step needs to be repeated everytime the Notebook and/or the PyRAF session is restarted.
In [1]:
import os
import os.path
import shutil
# WARNING: Make sure that the directory structure given in Step 0
# is already in place.
# EDIT HERE - Directories
rootdir = '/data/workspace/DRWorkshopAAO/NGC6872'
# EDIT HERE - File numbers and file names
biases = {
'bias_20101101' : {
'procbias' : 'S20101101S0169_bias.fits',
'frames' : [('S20101101S', '169-172'),
('S20101102S', '131-135'),
]
},
'bias_20101107' : {
'procbias' : 'S20101107S0228_bias.fits',
'frames' : [('S20101107S', '228-231'),
('S20101108S', '26-30')],
}
}
flats = {
'flat_20101107' : {
'procflat' : 'S20101107S0201_flat_r.fits',
'frames' : [('S20101107S', '201-212')],
'flatbias' : biases['bias_20101107']['procbias'],
}
}
science = {
'procsci' : 'ngc6872_r.fits',
'sciroot' : 'S20101101S',
'sciframes' : '87-90',
'scibias' : biases['bias_20101101']['procbias'],
'sciflat' : flats['flat_20101107']['procflat'],
#'scifringe' : fringes['fringe_20101111']['procfringe'],
}
logfile = 'NGC6872_r.log'
# Automatically set directories to match structure from Step 0
raw_path = os.path.join(rootdir, 'raw') + '/'
redux_path = os.path.join(rootdir, 'redux') + '/'
procbias_path = os.path.join(rootdir, 'calib', 'bias') + '/'
procflat_path = os.path.join(rootdir, 'calib', 'flat') + '/'
procfringe_path = os.path.join(rootdir, 'calib', 'fringe') + '/'
products_path = os.path.join(rootdir, 'products') + '/'
# Change to work directory
os.chdir(redux_path)
---> LAUNCH DS9 before proceeding further
Load the package required for the Notebook session and reset tasks to the default parameters.
In [2]:
from pyraf import iraf
iraf.gemini()
iraf.gmos()
iraf.unlearn(iraf.gemini, iraf.gmos, iraf.gemtools)
iraf.gmos.logfile = logfile
iraf.set(stdimage='imtgmos')
This is the initial release of the IRAF FITSUTIL package
to include support for FITS tile compression via 'fpack'.
Please send comments and questions to seaman@noao.edu.
+------------------- Gemini IRAF Package -------------------+
| Version 1.13.1, December 7, 2015 |
| Requires IRAF v2.14.1 or greater |
| Tested with Ureka IRAF v2.16 |
| Gemini Observatory, Hilo, Hawaii |
| Please use the help desk for submission of questions |
| http://www.gemini.edu/sciops/helpdesk/helpdeskIndex.html |
+-----------------------------------------------------------+
Warning setting imtype=fits
Warning setting use_new_imt=no
gemini/:
f2/ gemtools/ gsaoi/ niri/
flamingos/ gmos/ midir/ oscir/
gemlocal/ gnirs/ nifs/ quirc/
gmos/:
gbias gfquick gmosaic gsappwave gssdist
gbpm gfreduce gmosexamples gscalibrate gsskysub
gdark gfresponse gmosinfo gscrmask gsslitfunction
gdisplay gfscatsub gmosinfoifu gscrrej gsstandard
gfapsum gfskysub gmosinfoimag gscut gstransform
gfcombine gftransform gmosinfospec gsdrawslits gswavelength
gfcopy giflat gnscombine gsextract mostools/
gfcube gifringe gnsdark gsflat
gfdisplay giillum gnsskysub gsfquick
gfextract gireduce gprepare gsreduce
gffindblocks girmfringe gqecorr gsscatsub
The commands shown here need to be run from the shell.
Open an xterm or a Terminal, then:
cd /data/workspace/DRWorkshopAAO/NGC6872/redux
pyraf
Then in the PyRAF session:
gemini
gmos
unlearn gemini
unlearn gmos
unlearn gemtools
Get your PyRAF configured. In the PyRAF session:
iraf.gmos.logfile = "NGC6872_r.log"
set stdimage=imtgmos
In [11]:
os.chdir(redux_path)
# delete any pre-existing lists
iraf.delete('*.lis', verify='no')
# biases
for (root, frames) in biases['bias_20101101']['frames']:
iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'bias_20101101.lis')
iraf.delete('tmp_*.lis', verify='no')
for (root, frames) in biases['bias_20101107']['frames']:
iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'bias_20101107.lis')
iraf.delete('tmp_*.lis', verify='no')
# flats
for (root, frames) in flats['flat_20101107']['frames']:
iraf.gemlist(root, frames, Stdout='tmp_'+root+'.lis')
iraf.concat('tmp_*.lis', 'flat_20101107.lis')
iraf.delete('tmp_*.lis', verify='no')
# science
iraf.gemlist(science['sciroot'], science['sciframes'], \
Stdout='sci.lis')
# Get a full list
iraf.concat('*.lis', 'all.lis')
In [5]:
all_files = open('all.lis', 'r')
for line in all_files:
image = line.strip() + '[1]'
print image
iraf.gdisplay(raw_path + image, 1, fl_paste='no')
iraf.sleep(5)
all_files.close()
S20101101S0169[1]
GDISPLAY - Started: 2016-03-23T08:27:14
GDISPLAY: Temporary image name tmpout92482_394S20101101S0169.fits
GMULTIAMP - Started: 2016-03-23T08:27:15
GMULTIAMP - Finished: 2016-03-23T08:27:15
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:16
S20101101S0170[1]
GDISPLAY - Started: 2016-03-23T08:27:21
GDISPLAY: Temporary image name tmpout92482_415S20101101S0170.fits
GMULTIAMP - Started: 2016-03-23T08:27:21
GMULTIAMP - Finished: 2016-03-23T08:27:22
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:23
S20101101S0171[1]
GDISPLAY - Started: 2016-03-23T08:27:28
GDISPLAY: Temporary image name tmpout92482_436S20101101S0171.fits
GMULTIAMP - Started: 2016-03-23T08:27:28
GMULTIAMP - Finished: 2016-03-23T08:27:29
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:30
S20101101S0172[1]
GDISPLAY - Started: 2016-03-23T08:27:35
GDISPLAY: Temporary image name tmpout92482_457S20101101S0172.fits
GMULTIAMP - Started: 2016-03-23T08:27:35
GMULTIAMP - Finished: 2016-03-23T08:27:36
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:37
S20101102S0131[1]
GDISPLAY - Started: 2016-03-23T08:27:42
GDISPLAY: Temporary image name tmpout92482_478S20101102S0131.fits
GMULTIAMP - Started: 2016-03-23T08:27:42
GMULTIAMP - Finished: 2016-03-23T08:27:43
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:44
S20101102S0132[1]
GDISPLAY - Started: 2016-03-23T08:27:49
GDISPLAY: Temporary image name tmpout92482_499S20101102S0132.fits
GMULTIAMP - Started: 2016-03-23T08:27:49
GMULTIAMP - Finished: 2016-03-23T08:27:50
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:51
S20101102S0133[1]
GDISPLAY - Started: 2016-03-23T08:27:56
GDISPLAY: Temporary image name tmpout92482_520S20101102S0133.fits
GMULTIAMP - Started: 2016-03-23T08:27:56
GMULTIAMP - Finished: 2016-03-23T08:27:57
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:27:58
S20101102S0134[1]
GDISPLAY - Started: 2016-03-23T08:28:03
GDISPLAY: Temporary image name tmpout92482_541S20101102S0134.fits
GMULTIAMP - Started: 2016-03-23T08:28:03
GMULTIAMP - Finished: 2016-03-23T08:28:04
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:05
S20101102S0135[1]
GDISPLAY - Started: 2016-03-23T08:28:10
GDISPLAY: Temporary image name tmpout92482_562S20101102S0135.fits
GMULTIAMP - Started: 2016-03-23T08:28:10
GMULTIAMP - Finished: 2016-03-23T08:28:11
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:12
S20101107S0228[1]
GDISPLAY - Started: 2016-03-23T08:28:17
GDISPLAY: Temporary image name tmpout92482_583S20101107S0228.fits
GMULTIAMP - Started: 2016-03-23T08:28:17
GMULTIAMP - Finished: 2016-03-23T08:28:18
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:19
S20101107S0229[1]
GDISPLAY - Started: 2016-03-23T08:28:24
GDISPLAY: Temporary image name tmpout92482_604S20101107S0229.fits
GMULTIAMP - Started: 2016-03-23T08:28:24
GMULTIAMP - Finished: 2016-03-23T08:28:25
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:26
S20101107S0230[1]
GDISPLAY - Started: 2016-03-23T08:28:31
GDISPLAY: Temporary image name tmpout92482_625S20101107S0230.fits
GMULTIAMP - Started: 2016-03-23T08:28:31
GMULTIAMP - Finished: 2016-03-23T08:28:32
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:33
S20101107S0231[1]
GDISPLAY - Started: 2016-03-23T08:28:38
GDISPLAY: Temporary image name tmpout92482_646S20101107S0231.fits
GMULTIAMP - Started: 2016-03-23T08:28:38
GMULTIAMP - Finished: 2016-03-23T08:28:38
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:40
S20101107S0232[1]
GDISPLAY - Started: 2016-03-23T08:28:45
GDISPLAY: Temporary image name tmpout92482_667S20101107S0232.fits
GMULTIAMP - Started: 2016-03-23T08:28:45
GMULTIAMP - Finished: 2016-03-23T08:28:45
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:46
S20101107S0201[1]
GDISPLAY - Started: 2016-03-23T08:28:51
GDISPLAY: Temporary image name tmpout92482_688S20101107S0201.fits
GMULTIAMP - Started: 2016-03-23T08:28:52
GMULTIAMP - Finished: 2016-03-23T08:28:52
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:28:53
S20101107S0202[1]
GDISPLAY - Started: 2016-03-23T08:28:58
GDISPLAY: Temporary image name tmpout92482_709S20101107S0202.fits
GMULTIAMP - Started: 2016-03-23T08:28:59
GMULTIAMP - Finished: 2016-03-23T08:28:59
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:00
S20101107S0203[1]
GDISPLAY - Started: 2016-03-23T08:29:05
GDISPLAY: Temporary image name tmpout92482_730S20101107S0203.fits
GMULTIAMP - Started: 2016-03-23T08:29:05
GMULTIAMP - Finished: 2016-03-23T08:29:06
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:07
S20101107S0204[1]
GDISPLAY - Started: 2016-03-23T08:29:12
GDISPLAY: Temporary image name tmpout92482_751S20101107S0204.fits
GMULTIAMP - Started: 2016-03-23T08:29:12
GMULTIAMP - Finished: 2016-03-23T08:29:13
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:14
S20101107S0205[1]
GDISPLAY - Started: 2016-03-23T08:29:19
GDISPLAY: Temporary image name tmpout92482_772S20101107S0205.fits
GMULTIAMP - Started: 2016-03-23T08:29:19
GMULTIAMP - Finished: 2016-03-23T08:29:20
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:21
S20101107S0206[1]
GDISPLAY - Started: 2016-03-23T08:29:26
GDISPLAY: Temporary image name tmpout92482_793S20101107S0206.fits
GMULTIAMP - Started: 2016-03-23T08:29:26
GMULTIAMP - Finished: 2016-03-23T08:29:27
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:28
S20101107S0207[1]
GDISPLAY - Started: 2016-03-23T08:29:33
GDISPLAY: Temporary image name tmpout92482_814S20101107S0207.fits
GMULTIAMP - Started: 2016-03-23T08:29:33
GMULTIAMP - Finished: 2016-03-23T08:29:34
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:35
S20101107S0208[1]
GDISPLAY - Started: 2016-03-23T08:29:40
GDISPLAY: Temporary image name tmpout92482_835S20101107S0208.fits
GMULTIAMP - Started: 2016-03-23T08:29:40
GMULTIAMP - Finished: 2016-03-23T08:29:41
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:42
S20101107S0209[1]
GDISPLAY - Started: 2016-03-23T08:29:47
GDISPLAY: Temporary image name tmpout92482_856S20101107S0209.fits
GMULTIAMP - Started: 2016-03-23T08:29:47
GMULTIAMP - Finished: 2016-03-23T08:29:48
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:49
S20101107S0210[1]
GDISPLAY - Started: 2016-03-23T08:29:54
GDISPLAY: Temporary image name tmpout92482_877S20101107S0210.fits
GMULTIAMP - Started: 2016-03-23T08:29:54
GMULTIAMP - Finished: 2016-03-23T08:29:55
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:29:56
S20101107S0211[1]
GDISPLAY - Started: 2016-03-23T08:30:01
GDISPLAY: Temporary image name tmpout92482_898S20101107S0211.fits
GMULTIAMP - Started: 2016-03-23T08:30:01
GMULTIAMP - Finished: 2016-03-23T08:30:02
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:30:03
S20101107S0212[1]
GDISPLAY - Started: 2016-03-23T08:30:08
GDISPLAY: Temporary image name tmpout92482_919S20101107S0212.fits
GMULTIAMP - Started: 2016-03-23T08:30:08
GMULTIAMP - Finished: 2016-03-23T08:30:09
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:30:10
S20101101S0087[1]
GDISPLAY - Started: 2016-03-23T08:30:15
GDISPLAY: Temporary image name tmpout92482_940S20101101S0087.fits
GMULTIAMP - Started: 2016-03-23T08:30:15
GMULTIAMP - Finished: 2016-03-23T08:30:16
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:30:17
S20101101S0088[1]
GDISPLAY - Started: 2016-03-23T08:30:22
GDISPLAY: Temporary image name tmpout92482_961S20101101S0088.fits
GMULTIAMP - Started: 2016-03-23T08:30:22
GMULTIAMP - Finished: 2016-03-23T08:30:23
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:30:24
S20101101S0089[1]
GDISPLAY - Started: 2016-03-23T08:30:29
GDISPLAY: Temporary image name tmpout92482_982S20101101S0089.fits
GMULTIAMP - Started: 2016-03-23T08:30:29
GMULTIAMP - Finished: 2016-03-23T08:30:30
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:30:31
S20101101S0090[1]
GDISPLAY - Started: 2016-03-23T08:30:36
GDISPLAY: Temporary image name tmpout92482_1003S20101101S0090.fits
GMULTIAMP - Started: 2016-03-23T08:30:36
GMULTIAMP - Finished: 2016-03-23T08:30:36
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY - Finished: 2016-03-23T08:30:37
This will produce the two biases defined above. Those two biases will be used later to process the science frames, the twilight flats, and the fringe frames.
In [8]:
os.chdir(redux_path)
for biaskey in biases.keys():
iraf.imdelete(biases[biaskey]['procbias'], verify='no')
iraf.imdelete('g//@'+biaskey+'.lis', verify='no')
print 'Creating bias ', biases[biaskey]['procbias']
iraf.gbias('@'+biaskey+'.lis', outbias=biases[biaskey]['procbias'], \
rawpath=raw_path, fl_over='yes', fl_trim='yes', \
fl_vardq='yes')
shutil.copy(biases[biaskey]['procbias'], procbias_path)
# The processed biases are:
os.listdir(procbias_path)
Creating bias S20101101S0169_bias.fits
----------------------------------------------------------------------------
GBIAS -- Wed Mar 23 08:33:17 HST 2016
Input images:
S20101101S0169
S20101101S0170
S20101101S0171
S20101101S0172
S20101102S0131
S20101102S0132
S20101102S0133
S20101102S0134
S20101102S0135
Output bias image: S20101101S0169_bias.fits
----------------------------------------------------------------------------
GIREDUCE -- Wed Mar 23 08:33:17 HST 2016
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:17 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101101S0169
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0169.fits Output gS20101101S0169.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:18 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101101S0170
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0170.fits Output gS20101101S0170.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:19 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101101S0171
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0171.fits Output gS20101101S0171.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:20 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101101S0172
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0172.fits Output gS20101101S0172.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:21 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101102S0131
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101102S0131.fits Output gS20101102S0131.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:22 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101102S0132
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101102S0132.fits Output gS20101102S0132.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:23 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101102S0133
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101102S0133.fits Output gS20101102S0133.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:24 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101102S0134
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101102S0134.fits Output gS20101102S0134.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:33:25 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101102S0135
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101102S0135.fits Output gS20101102S0135.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
Input files:
gS20101101S0169
gS20101101S0170
gS20101101S0171
gS20101101S0172
gS20101102S0131
gS20101102S0132
gS20101102S0133
gS20101102S0134
gS20101102S0135
Output files:
tmpfile92482_1139S20101101S0169
tmpfile92482_1139S20101101S0170
tmpfile92482_1139S20101101S0171
tmpfile92482_1139S20101101S0172
tmpfile92482_1139S20101102S0131
tmpfile92482_1139S20101102S0132
tmpfile92482_1139S20101102S0133
tmpfile92482_1139S20101102S0134
tmpfile92482_1139S20101102S0135
GIREDUCE: Image tmpfile92482_1139S20101101S0169 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101101S0170 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101101S0171 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101101S0172 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101102S0131 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101102S0132 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101102S0133 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101102S0134 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101102S0135 overscan subtracted
GIREDUCE: Image tmpfile92482_1139S20101101S0169 trimmed
GIREDUCE: Image tmpfile92482_1139S20101101S0170 trimmed
GIREDUCE: Image tmpfile92482_1139S20101101S0171 trimmed
GIREDUCE: Image tmpfile92482_1139S20101101S0172 trimmed
GIREDUCE: Image tmpfile92482_1139S20101102S0131 trimmed
GIREDUCE: Image tmpfile92482_1139S20101102S0132 trimmed
GIREDUCE: Image tmpfile92482_1139S20101102S0133 trimmed
GIREDUCE: Image tmpfile92482_1139S20101102S0134 trimmed
GIREDUCE: Image tmpfile92482_1139S20101102S0135 trimmed
GIREDUCE - Cleaning up -- Wed Mar 23 08:34:16 HST 2016
GIREDUCE -- Wed Mar 23 08:34:16 HST 2016
GIREDUCE exit status: good.
----------------------------------------------------------------------------
--------------------------------------------------------------------------------
GEMCOMBINE -- Wed Mar 23 08:34:16 HST 2016
GEMCOMBINE: input = @tmpcombin92482_1136
GEMCOMBINE: output = S20101101S0169_bias.fits
GEMCOMBINE: title = Bias
GEMCOMBINE: bpmfile =
GEMCOMBINE: sci_ext = SCI
GEMCOMBINE: var_ext = VAR
GEMCOMBINE: dq_ext = DQ
Mar 23 8:34: IMCOMBINE
combine = average, scale = none, zero = none, weight = none
reject = avsigclip, mclip = yes, nkeep = 1
lsigma = 3., hsigma = 3.
blank = 0.
statsec = Mar 23 8:34
Images Exp
tmpfile92482_1139S20101101S0169.fits[SCI,1] 0.0
tmpfile92482_1139S20101101S0170.fits[SCI,1] 0.0
tmpfile92482_1139S20101101S0171.fits[SCI,1] 0.0
tmpfile92482_1139S20101101S0172.fits[SCI,1] 0.0
tmpfile92482_1139S20101102S0131.fits[SCI,1] 0.0
tmpfile92482_1139S20101102S0132.fits[SCI,1] 0.0
tmpfile92482_1139S20101102S0133.fits[SCI,1] 0.0
tmpfile92482_1139S20101102S0134.fits[SCI,1] 0.0
tmpfile92482_1139S20101102S0135.fits[SCI,1] 0.0
Output image = tmpcombout92482_1213, ncombine = 9
Number rejected mask = tmpcombdq92482_1214.pl
Sigma image = tmpcombsig92482_1215
Mar 23 8:34: IMCOMBINE
combine = average, scale = none, zero = none, weight = none
reject = avsigclip, mclip = yes, nkeep = 1
lsigma = 3., hsigma = 3.
blank = 0.
statsec = Mar 23 8:34
Images Exp
tmpfile92482_1139S20101101S0169.fits[SCI,2] 0.0
tmpfile92482_1139S20101101S0170.fits[SCI,2] 0.0
tmpfile92482_1139S20101101S0171.fits[SCI,2] 0.0
tmpfile92482_1139S20101101S0172.fits[SCI,2] 0.0
tmpfile92482_1139S20101102S0131.fits[SCI,2] 0.0
tmpfile92482_1139S20101102S0132.fits[SCI,2] 0.0
tmpfile92482_1139S20101102S0133.fits[SCI,2] 0.0
tmpfile92482_1139S20101102S0134.fits[SCI,2] 0.0
tmpfile92482_1139S20101102S0135.fits[SCI,2] 0.0
Output image = tmpcombout92482_1220, ncombine = 9
Number rejected mask = tmpcombdq92482_1221.pl
Sigma image = tmpcombsig92482_1222
Mar 23 8:34: IMCOMBINE
combine = average, scale = none, zero = none, weight = none
reject = avsigclip, mclip = yes, nkeep = 1
lsigma = 3., hsigma = 3.
blank = 0.
statsec = Mar 23 8:34
Images Exp
tmpfile92482_1139S20101101S0169.fits[SCI,3] 0.0
tmpfile92482_1139S20101101S0170.fits[SCI,3] 0.0
tmpfile92482_1139S20101101S0171.fits[SCI,3] 0.0
tmpfile92482_1139S20101101S0172.fits[SCI,3] 0.0
tmpfile92482_1139S20101102S0131.fits[SCI,3] 0.0
tmpfile92482_1139S20101102S0132.fits[SCI,3] 0.0
tmpfile92482_1139S20101102S0133.fits[SCI,3] 0.0
tmpfile92482_1139S20101102S0134.fits[SCI,3] 0.0
tmpfile92482_1139S20101102S0135.fits[SCI,3] 0.0
Output image = tmpcombout92482_1224, ncombine = 9
Number rejected mask = tmpcombdq92482_1225.pl
Sigma image = tmpcombsig92482_1226
GEMCOMBINE -- Finished: Wed Mar 23 08:34:28 HST 2016
GEMCOMBINE -- Exit status: GOOD
--------------------------------------------------------------------------------
GBIAS: Exit status: GOOD.
----------------------------------------------------------------------------
Creating bias S20101107S0228_bias.fits
----------------------------------------------------------------------------
GBIAS -- Wed Mar 23 08:34:29 HST 2016
Input images:
S20101107S0228
S20101107S0229
S20101107S0230
S20101107S0231
S20101107S0232
Output bias image: S20101107S0228_bias.fits
----------------------------------------------------------------------------
GIREDUCE -- Wed Mar 23 08:34:29 HST 2016
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:34:29 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101107S0228
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0228.fits Output gS20101107S0228.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:34:30 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101107S0229
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0229.fits Output gS20101107S0229.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:34:31 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101107S0230
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0230.fits Output gS20101107S0230.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:34:32 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101107S0231
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0231.fits Output gS20101107S0231.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:34:33 HST 2016
Input list = /data/workspace/DRWorkshopAAO/raw/S20101107S0232
Output list =
Output prefix = g
Raw path =
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0232.fits Output gS20101107S0232.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
Input files:
gS20101107S0228
gS20101107S0229
gS20101107S0230
gS20101107S0231
gS20101107S0232
Output files:
tmpfile92482_1234S20101107S0228
tmpfile92482_1234S20101107S0229
tmpfile92482_1234S20101107S0230
tmpfile92482_1234S20101107S0231
tmpfile92482_1234S20101107S0232
GIREDUCE: Image tmpfile92482_1234S20101107S0228 overscan subtracted
GIREDUCE: Image tmpfile92482_1234S20101107S0229 overscan subtracted
GIREDUCE: Image tmpfile92482_1234S20101107S0230 overscan subtracted
GIREDUCE: Image tmpfile92482_1234S20101107S0231 overscan subtracted
GIREDUCE: Image tmpfile92482_1234S20101107S0232 overscan subtracted
GIREDUCE: Image tmpfile92482_1234S20101107S0228 trimmed
GIREDUCE: Image tmpfile92482_1234S20101107S0229 trimmed
GIREDUCE: Image tmpfile92482_1234S20101107S0230 trimmed
GIREDUCE: Image tmpfile92482_1234S20101107S0231 trimmed
GIREDUCE: Image tmpfile92482_1234S20101107S0232 trimmed
GIREDUCE - Cleaning up -- Wed Mar 23 08:35:02 HST 2016
GIREDUCE -- Wed Mar 23 08:35:02 HST 2016
GIREDUCE exit status: good.
----------------------------------------------------------------------------
--------------------------------------------------------------------------------
GEMCOMBINE -- Wed Mar 23 08:35:02 HST 2016
GEMCOMBINE: input = @tmpcombin92482_1231
GEMCOMBINE: output = S20101107S0228_bias.fits
GEMCOMBINE: title = Bias
GEMCOMBINE: bpmfile =
GEMCOMBINE: sci_ext = SCI
GEMCOMBINE: var_ext = VAR
GEMCOMBINE: dq_ext = DQ
WARNING - GEMCOMBINE: five or less images to combine.
Mar 23 8:35: IMCOMBINE
combine = average, scale = none, zero = none, weight = none
reject = avsigclip, mclip = yes, nkeep = 1
lsigma = 3., hsigma = 3.
blank = 0.
statsec = Mar 23 8:35
Images Exp
tmpfile92482_1234S20101107S0228.fits[SCI,1] 0.0
tmpfile92482_1234S20101107S0229.fits[SCI,1] 0.0
tmpfile92482_1234S20101107S0230.fits[SCI,1] 0.0
tmpfile92482_1234S20101107S0231.fits[SCI,1] 0.0
tmpfile92482_1234S20101107S0232.fits[SCI,1] 5.2
Output image = tmpcombout92482_1284, ncombine = 5
Number rejected mask = tmpcombdq92482_1285.pl
Sigma image = tmpcombsig92482_1286
WARNING - GEMCOMBINE: five or less images to combine.
Mar 23 8:35: IMCOMBINE
combine = average, scale = none, zero = none, weight = none
reject = avsigclip, mclip = yes, nkeep = 1
lsigma = 3., hsigma = 3.
blank = 0.
statsec = Mar 23 8:35
Images Exp
tmpfile92482_1234S20101107S0228.fits[SCI,2] 0.0
tmpfile92482_1234S20101107S0229.fits[SCI,2] 0.0
tmpfile92482_1234S20101107S0230.fits[SCI,2] 0.0
tmpfile92482_1234S20101107S0231.fits[SCI,2] 0.0
tmpfile92482_1234S20101107S0232.fits[SCI,2] 5.2
Output image = tmpcombout92482_1291, ncombine = 5
Number rejected mask = tmpcombdq92482_1292.pl
Sigma image = tmpcombsig92482_1293
WARNING - GEMCOMBINE: five or less images to combine.
Mar 23 8:35: IMCOMBINE
combine = average, scale = none, zero = none, weight = none
reject = avsigclip, mclip = yes, nkeep = 1
lsigma = 3., hsigma = 3.
blank = 0.
statsec = Mar 23 8:35
Images Exp
tmpfile92482_1234S20101107S0228.fits[SCI,3] 0.0
tmpfile92482_1234S20101107S0229.fits[SCI,3] 0.0
tmpfile92482_1234S20101107S0230.fits[SCI,3] 0.0
tmpfile92482_1234S20101107S0231.fits[SCI,3] 0.0
tmpfile92482_1234S20101107S0232.fits[SCI,3] 5.2
Output image = tmpcombout92482_1295, ncombine = 5
Number rejected mask = tmpcombdq92482_1296.pl
Sigma image = tmpcombsig92482_1297
GEMCOMBINE -- Finished: Wed Mar 23 08:35:11 HST 2016
GEMCOMBINE -- Exit status: GOOD
--------------------------------------------------------------------------------
GBIAS: Exit status: GOOD.
----------------------------------------------------------------------------
Out[8]:
['S20101101S0169_bias.fits',
'S20101107S0228_bias.fits',
'S20101111S0061_bias.fits']
In [9]:
# WARNING: interactive. It launches imexam. Type 'q' to quit.
for biaskey in biases.keys():
iraf.gdisplay(biases[biaskey]['procbias'], 1, fl_paste='yes')
GDISPLAY - Started: 2016-03-23T08:35:32
GDISPLAY: Temporary image name tmpout92482_1300S20101101S0169_bias.fits
GMULTIAMP - Started: 2016-03-23T08:35:32
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T08:35:34
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T08:35:39
GDISPLAY - Started: 2016-03-23T08:35:39
GDISPLAY: Temporary image name tmpout92482_1334S20101107S0228_bias.fits
GMULTIAMP - Started: 2016-03-23T08:35:40
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T08:35:41
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T08:35:43
This will produce the one twilight flat defined above. That flat will be used later to process the science frames and the fringe frames.
In [12]:
os.chdir(redux_path)
for flatkey in flats.keys():
iraf.imdelete('g//@'+flatkey+'.lis', verify='no')
iraf.imdelete('rg//@'+flatkey+'.lis', verify='no')
flatlist = flatkey+'.lis'
procbias = os.path.join(procbias_path, flats[flatkey]['flatbias'])
print 'Processing raw twilight flats'
# gprepare and bias correct the raw flats
iraf.gireduce('@'+flatlist, rawpath=raw_path, \
fl_over='yes', fl_trim='yes', \
fl_bias='yes', bias=procbias, fl_flat='no', \
fl_vardq='yes')
Processing raw twilight flats
----------------------------------------------------------------------------
GIREDUCE -- Wed Mar 23 08:36:32 HST 2016
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:32 HST 2016
Input list = S20101107S0201
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0201.fits Output gS20101107S0201.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:33 HST 2016
Input list = S20101107S0202
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0202.fits Output gS20101107S0202.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:34 HST 2016
Input list = S20101107S0203
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0203.fits Output gS20101107S0203.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:35 HST 2016
Input list = S20101107S0204
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0204.fits Output gS20101107S0204.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:36 HST 2016
Input list = S20101107S0205
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0205.fits Output gS20101107S0205.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:37 HST 2016
Input list = S20101107S0206
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0206.fits Output gS20101107S0206.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:38 HST 2016
Input list = S20101107S0207
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0207.fits Output gS20101107S0207.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:39 HST 2016
Input list = S20101107S0208
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0208.fits Output gS20101107S0208.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:40 HST 2016
Input list = S20101107S0209
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0209.fits Output gS20101107S0209.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:41 HST 2016
Input list = S20101107S0210
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0210.fits Output gS20101107S0210.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:42 HST 2016
Input list = S20101107S0211
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0211.fits Output gS20101107S0211.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:36:43 HST 2016
Input list = S20101107S0212
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101107S0212.fits Output gS20101107S0212.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
Input files:
gS20101107S0201
gS20101107S0202
gS20101107S0203
gS20101107S0204
gS20101107S0205
gS20101107S0206
gS20101107S0207
gS20101107S0208
gS20101107S0209
gS20101107S0210
gS20101107S0211
gS20101107S0212
Output files:
rgS20101107S0201
rgS20101107S0202
rgS20101107S0203
rgS20101107S0204
rgS20101107S0205
rgS20101107S0206
rgS20101107S0207
rgS20101107S0208
rgS20101107S0209
rgS20101107S0210
rgS20101107S0211
rgS20101107S0212
WARNING - GIREDUCE: Bad Pixel Mask filename is an empty string
Only saturated pixels will be flagged
GIREDUCE: Image rgS20101107S0201 overscan subtracted
GIREDUCE: Image rgS20101107S0202 overscan subtracted
GIREDUCE: Image rgS20101107S0203 overscan subtracted
GIREDUCE: Image rgS20101107S0204 overscan subtracted
GIREDUCE: Image rgS20101107S0205 overscan subtracted
GIREDUCE: Image rgS20101107S0206 overscan subtracted
GIREDUCE: Image rgS20101107S0207 overscan subtracted
GIREDUCE: Image rgS20101107S0208 overscan subtracted
GIREDUCE: Image rgS20101107S0209 overscan subtracted
GIREDUCE: Image rgS20101107S0210 overscan subtracted
GIREDUCE: Image rgS20101107S0211 overscan subtracted
GIREDUCE: Image rgS20101107S0212 overscan subtracted
GIREDUCE: Image rgS20101107S0201 trimmed
GIREDUCE: Image rgS20101107S0202 trimmed
GIREDUCE: Image rgS20101107S0203 trimmed
GIREDUCE: Image rgS20101107S0204 trimmed
GIREDUCE: Image rgS20101107S0205 trimmed
GIREDUCE: Image rgS20101107S0206 trimmed
GIREDUCE: Image rgS20101107S0207 trimmed
GIREDUCE: Image rgS20101107S0208 trimmed
GIREDUCE: Image rgS20101107S0209 trimmed
GIREDUCE: Image rgS20101107S0210 trimmed
GIREDUCE: Image rgS20101107S0211 trimmed
GIREDUCE: Image rgS20101107S0212 trimmed
Output image Bias Flat Dark Scale
rgS20101107S0201 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0201
GIREDUCE: multiplying image rgS20101107S0201 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0202 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0202
GIREDUCE: multiplying image rgS20101107S0202 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0203 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0203
GIREDUCE: multiplying image rgS20101107S0203 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0204 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0204
GIREDUCE: multiplying image rgS20101107S0204 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0205 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0205
GIREDUCE: multiplying image rgS20101107S0205 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0206 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0206
GIREDUCE: multiplying image rgS20101107S0206 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0207 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0207
GIREDUCE: multiplying image rgS20101107S0207 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0208 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0208
GIREDUCE: multiplying image rgS20101107S0208 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0209 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0209
GIREDUCE: multiplying image rgS20101107S0209 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0210 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0210
GIREDUCE: multiplying image rgS20101107S0210 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0211 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0211
GIREDUCE: multiplying image rgS20101107S0211 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101107S0212 /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias INDEF INDEF 0.00
GIREDUCE: Subtracted bias /data/workspace/DRWorkshopAAO/calib/bias/S20101107S0228_bias from rgS20101107S0212
GIREDUCE: multiplying image rgS20101107S0212 by gain
GIREDUCE: output counts in electrons
GIREDUCE - Cleaning up -- Wed Mar 23 08:39:24 HST 2016
GIREDUCE -- Wed Mar 23 08:39:24 HST 2016
GIREDUCE exit status: good.
----------------------------------------------------------------------------
In [13]:
os.chdir(redux_path)
for flatkey in flats.keys():
iraf.imdelete(flats[flatkey]['procflat'], verify='no')
flatlist = flatkey+'.lis'
procflat = flats[flatkey]['procflat']
print 'Creating twilight flat ', procflat
# Create the processed flat
iraf.giflat('rg//@'+flatlist, outflat=procflat, fl_vardq='yes')
# Copy the processed flat to the calibration directory
shutil.copy(procflat, procflat_path)
# The processed flats are:
os.listdir(procflat_path)
Creating twilight flat S20101107S0201_flat_r.fits
----------------------------------------------------------------------------
GIFLAT -- Wed Mar 23 08:56:45 HST 2016
Flat field images = rg//@flat_20101107.lis
Output flat field image = S20101107S0201_flat_r.fits
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0201
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0202
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0203
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0204
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0205
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0206
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0207
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0208
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0209
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0210
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0211
WARNING - GIFLAT: using previously gireduce'd image: rgS20101107S0212
Using input files:
rgS20101107S0201
rgS20101107S0202
rgS20101107S0203
rgS20101107S0204
rgS20101107S0205
rgS20101107S0206
rgS20101107S0207
rgS20101107S0208
rgS20101107S0209
rgS20101107S0210
rgS20101107S0211
rgS20101107S0212
Using STATSEC: [50:900,50:2250] from CCD2
Using NORMSEC: [50:900,50:2250] from CCD2
GIFLAT: Normalization factor 60768.7
GIFLAT exit status: good.
----------------------------------------------------------------------------
Out[13]:
['S20101102S0100_flat_i.fits', 'S20101107S0201_flat_r.fits']
In [14]:
for flatkey in flats.keys():
iraf.gdisplay(flats[flatkey]['procflat'], 1, fl_paste='yes')
GDISPLAY - Started: 2016-03-23T08:57:37
GDISPLAY: Temporary image name tmpout92482_1783S20101107S0201_flat_r.fits
GMULTIAMP - Started: 2016-03-23T08:57:37
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T08:57:39
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T08:57:48
There is no significant fringing in the r-band. No extra fringing correction will be necessary. If you do want to learn about how to correct for fringing, see the i-band notebook.
All we need to do here is the straighforward processing, i.e. the bias and flat correction. (And overscan and trim.)
In [15]:
os.chdir(redux_path)
iraf.imdelete('g//@sci.lis', verify='no')
iraf.imdelete('rg//@sci.lis', verify='no')
iraf.gireduce('@sci.lis', rawpath=raw_path, fl_over='yes', fl_trim='yes', \
fl_bias='yes', bias=science['scibias'], \
fl_flat='yes', flat1=science['sciflat'], \
fl_vardq='yes')
----------------------------------------------------------------------------
GIREDUCE -- Wed Mar 23 08:58:01 HST 2016
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:58:01 HST 2016
Input list = S20101101S0087
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0087.fits Output gS20101101S0087.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:58:02 HST 2016
Input list = S20101101S0088
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0088.fits Output gS20101101S0088.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:58:03 HST 2016
Input list = S20101101S0089
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0089.fits Output gS20101101S0089.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
----------------------------------------------------------------------------
GPREPARE -- Wed Mar 23 08:58:04 HST 2016
Input list = S20101101S0090
Output list =
Output prefix = g
Raw path = /data/workspace/DRWorkshopAAO/raw/
MDF dir = gmos$data/
Add MDF = no
Input /data/workspace/DRWorkshopAAO/raw/S20101101S0090.fits Output gS20101101S0090.fits
GPREPARE exit status: good.
----------------------------------------------------------------------------
Input files:
gS20101101S0087
gS20101101S0088
gS20101101S0089
gS20101101S0090
Output files:
rgS20101101S0087
rgS20101101S0088
rgS20101101S0089
rgS20101101S0090
WARNING - GIREDUCE: Bad Pixel Mask filename is an empty string
Only saturated pixels will be flagged
GIREDUCE: Image rgS20101101S0087 overscan subtracted
GIREDUCE: Image rgS20101101S0088 overscan subtracted
GIREDUCE: Image rgS20101101S0089 overscan subtracted
GIREDUCE: Image rgS20101101S0090 overscan subtracted
GIREDUCE: Image rgS20101101S0087 trimmed
GIREDUCE: Image rgS20101101S0088 trimmed
GIREDUCE: Image rgS20101101S0089 trimmed
GIREDUCE: Image rgS20101101S0090 trimmed
Output image Bias Flat Dark Scale
rgS20101101S0087 S20101101S0169_bias S20101107S0201_flat_r INDEF 0.00
GIREDUCE: Subtracted bias S20101101S0169_bias from rgS20101101S0087
GIREDUCE: Divided image rgS20101101S0087 by flatfield S20101107S0201_flat_r
GIREDUCE: multiplying image rgS20101101S0087 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101101S0088 S20101101S0169_bias S20101107S0201_flat_r INDEF 0.00
GIREDUCE: Subtracted bias S20101101S0169_bias from rgS20101101S0088
GIREDUCE: Divided image rgS20101101S0088 by flatfield S20101107S0201_flat_r
GIREDUCE: multiplying image rgS20101101S0088 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101101S0089 S20101101S0169_bias S20101107S0201_flat_r INDEF 0.00
GIREDUCE: Subtracted bias S20101101S0169_bias from rgS20101101S0089
GIREDUCE: Divided image rgS20101101S0089 by flatfield S20101107S0201_flat_r
GIREDUCE: multiplying image rgS20101101S0089 by gain
GIREDUCE: output counts in electrons
Output image Bias Flat Dark Scale
rgS20101101S0090 S20101101S0169_bias S20101107S0201_flat_r INDEF 0.00
GIREDUCE: Subtracted bias S20101101S0169_bias from rgS20101101S0090
GIREDUCE: Divided image rgS20101101S0090 by flatfield S20101107S0201_flat_r
GIREDUCE: multiplying image rgS20101101S0090 by gain
GIREDUCE: output counts in electrons
GIREDUCE - Cleaning up -- Wed Mar 23 08:59:05 HST 2016
GIREDUCE -- Wed Mar 23 08:59:05 HST 2016
GIREDUCE exit status: good.
----------------------------------------------------------------------------
In [16]:
sciframes = open('sci.lis', 'r')
for line in sciframes:
frame = line.strip()
iraf.gdisplay('rg'+frame, 1, fl_paste='yes')
sciframes.close()
GDISPLAY - Started: 2016-03-23T09:00:38
GDISPLAY: Temporary image name tmpout92482_1875rgS20101101S0087.fits
GMULTIAMP - Started: 2016-03-23T09:00:39
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T09:00:41
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:00:50
GDISPLAY - Started: 2016-03-23T09:00:50
GDISPLAY: Temporary image name tmpout92482_1909rgS20101101S0088.fits
GMULTIAMP - Started: 2016-03-23T09:00:50
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T09:00:52
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:00:53
GDISPLAY - Started: 2016-03-23T09:00:53
GDISPLAY: Temporary image name tmpout92482_1943rgS20101101S0089.fits
GMULTIAMP - Started: 2016-03-23T09:00:54
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T09:00:56
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:00:58
GDISPLAY - Started: 2016-03-23T09:00:58
GDISPLAY: Temporary image name tmpout92482_1977rgS20101101S0090.fits
GMULTIAMP - Started: 2016-03-23T09:00:58
GMULTIAMP: Tiling image...
Image tiled.
GMULTIAMP - Finished: 2016-03-23T09:01:00
GMULTIAMP - Exit status: GOOD
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:01:03
In [17]:
iraf.imdelete('mrg//@sci.lis', verify='no')
iraf.gmosaic('rg//@sci.lis', fl_vardq='yes')
GMOSAIC
GMOSAIC -------------------------------------------------------------------
GMOSAIC Log opened at [Wed 09:01:09 23-Mar-2016]
GMOSAIC
GMOSAIC --------------------
GMOSAIC Input Parameters:
GMOSAIC inimages = rg//@sci.lis
GMOSAIC outimages =
GMOSAIC outpref = m
GMOSAIC fl_paste = no
GMOSAIC fl_vardq = yes
GMOSAIC fl_fixpix = no
GMOSAIC fl_clean = yes
GMOSAIC fl_fulldq = no
GMOSAIC bitflags = all
GMOSAIC dqthresh = 0.1
GMOSAIC geointer = linear
GMOSAIC gap = default
GMOSAIC bpmfile = gmos$data/chipgaps.dat
GMOSAIC statsec = default
GMOSAIC obsmode = IMAGE
GMOSAIC sci_ext = SCI
GMOSAIC var_ext = VAR
GMOSAIC dq_ext = DQ
GMOSAIC key_detsec = DETSEC
GMOSAIC key_ccdsec = CCDSEC
GMOSAIC key_datsec = DATASEC
GMOSAIC key_ccdsum = CCDSUM
GMOSAIC key_obsmode = OBSMODE
GMOSAIC logfile =
GMOSAIC verbose = yes
GMOSAIC fl_real = no
GMOSAIC --------------------
GMOSAIC Input: rgS20101101S0087 Output: mrgS20101101S0087.fits
GMOSAIC Setting rotation Xrot[1]= 0.011000, Yrot[1]= 0.011000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]= 0.012000, Yrot[3]= 0.012000
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[SCI,1][1:1024,1:2304] to image tmpchipsci92482_2029
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[VAR,1][1:1024,1:2304] to image tmpchipvar92482_2030
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[DQ,1][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[SCI,2][1:1024,1:2304] to image tmpchipsci92482_2032
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[VAR,2][1:1024,1:2304] to image tmpchipvar92482_2033
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[SCI,3][1:1024,1:2304] to image tmpchipsci92482_2035
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[VAR,3][1:1024,1:2304] to image tmpchipvar92482_2036
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2022[DQ,3][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC imcopy tmpchipsci92482_2029 mrgS20101101S0087.fits[1,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2032 mrgS20101101S0087.fits[1,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2035 mrgS20101101S0087.fits[1,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2030 mrgS20101101S0087.fits[2,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2033 mrgS20101101S0087.fits[2,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2036 mrgS20101101S0087.fits[2,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2031 mrgS20101101S0087.fits[3,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2034 mrgS20101101S0087.fits[3,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2037 mrgS20101101S0087.fits[3,append][2085:3108,1:2304] INDEF 0.
gmos$data/gmos-s_bpm_EEV_22_3amp_v2_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3828.553
GMOSAIC Input: rgS20101101S0088 Output: mrgS20101101S0088.fits
GMOSAIC Setting rotation Xrot[1]= 0.011000, Yrot[1]= 0.011000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]= 0.012000, Yrot[3]= 0.012000
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[SCI,1][1:1024,1:2304] to image tmpchipsci92482_2045
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[VAR,1][1:1024,1:2304] to image tmpchipvar92482_2046
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[DQ,1][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[SCI,2][1:1024,1:2304] to image tmpchipsci92482_2048
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[VAR,2][1:1024,1:2304] to image tmpchipvar92482_2049
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[SCI,3][1:1024,1:2304] to image tmpchipsci92482_2051
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[VAR,3][1:1024,1:2304] to image tmpchipvar92482_2052
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2038[DQ,3][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC imcopy tmpchipsci92482_2045 mrgS20101101S0088.fits[1,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2048 mrgS20101101S0088.fits[1,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2051 mrgS20101101S0088.fits[1,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2046 mrgS20101101S0088.fits[2,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2049 mrgS20101101S0088.fits[2,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2052 mrgS20101101S0088.fits[2,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2047 mrgS20101101S0088.fits[3,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2050 mrgS20101101S0088.fits[3,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2053 mrgS20101101S0088.fits[3,append][2085:3108,1:2304] INDEF 0.
gmos$data/gmos-s_bpm_EEV_22_3amp_v2_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3784.769
GMOSAIC Input: rgS20101101S0089 Output: mrgS20101101S0089.fits
GMOSAIC Setting rotation Xrot[1]= 0.011000, Yrot[1]= 0.011000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]= 0.012000, Yrot[3]= 0.012000
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[SCI,1][1:1024,1:2304] to image tmpchipsci92482_2061
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[VAR,1][1:1024,1:2304] to image tmpchipvar92482_2062
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[DQ,1][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[SCI,2][1:1024,1:2304] to image tmpchipsci92482_2064
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[VAR,2][1:1024,1:2304] to image tmpchipvar92482_2065
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[SCI,3][1:1024,1:2304] to image tmpchipsci92482_2067
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[VAR,3][1:1024,1:2304] to image tmpchipvar92482_2068
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2054[DQ,3][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC imcopy tmpchipsci92482_2061 mrgS20101101S0089.fits[1,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2064 mrgS20101101S0089.fits[1,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2067 mrgS20101101S0089.fits[1,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2062 mrgS20101101S0089.fits[2,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2065 mrgS20101101S0089.fits[2,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2068 mrgS20101101S0089.fits[2,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2063 mrgS20101101S0089.fits[3,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2066 mrgS20101101S0089.fits[3,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2069 mrgS20101101S0089.fits[3,append][2085:3108,1:2304] INDEF 0.
gmos$data/gmos-s_bpm_EEV_22_3amp_v2_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3802.435
GMOSAIC Input: rgS20101101S0090 Output: mrgS20101101S0090.fits
GMOSAIC Setting rotation Xrot[1]= 0.011000, Yrot[1]= 0.011000
GMOSAIC Setting rotation Xrot[2]= 0.000000, Yrot[2]= 0.000000
GMOSAIC Setting rotation Xrot[3]= 0.012000, Yrot[3]= 0.012000
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[SCI,1][1:1024,1:2304] to image tmpchipsci92482_2077
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[VAR,1][1:1024,1:2304] to image tmpchipvar92482_2078
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[DQ,1][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: 1.02 yshift: 0.21 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[SCI,2][1:1024,1:2304] to image tmpchipsci92482_2080
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[VAR,2][1:1024,1:2304] to image tmpchipvar92482_2081
GMOSAIC xshift: 0.00 yshift: 0.00 xmag: 1.00 ymag: 1.00 xrot: 0.00 yrot: 0.00
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[SCI,3][1:1024,1:2304] to image tmpchipsci92482_2083
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[VAR,3][1:1024,1:2304] to image tmpchipvar92482_2084
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC
GMOSAIC Transforming image tmpimg92482_2070[DQ,3][1:1024,1:2304] to image tmpdq_transimg92482_2018.fits
GMOSAIC xshift: -2.90 yshift: -0.91 xmag: 1.00 ymag: 1.00 xrot: 0.01 yrot: 0.01
GMOSAIC imcopy tmpchipsci92482_2077 mrgS20101101S0090.fits[1,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2080 mrgS20101101S0090.fits[1,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipsci92482_2083 mrgS20101101S0090.fits[1,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2078 mrgS20101101S0090.fits[2,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2081 mrgS20101101S0090.fits[2,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipvar92482_2084 mrgS20101101S0090.fits[2,append][2085:3108,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2079 mrgS20101101S0090.fits[3,append][1:1024,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2082 mrgS20101101S0090.fits[3,append][1043:2066,1:2304] INDEF 0.
GMOSAIC imcopy tmpchipdq92482_2085 mrgS20101101S0090.fits[3,append][2085:3108,1:2304] INDEF 0.
gmos$data/gmos-s_bpm_EEV_22_3amp_v2_mosaic.pl
GMOSAIC Cleaning areas outside imaging field, cleaning value = 3737.172
GMOSAIC All 4 images successfully mosaiced.
GMOSAIC --------------------
GMOSAIC Exit status: SUCCESS
GMOSAIC Log closed at [Wed 09:01:49 23-Mar-2016]
GMOSAIC -------------------------------------------------------------------
GMOSAIC
In [18]:
sciframes = open('sci.lis', 'r')
for line in sciframes:
frame = line.strip()
iraf.gdisplay('mrg'+frame, 1, z2=20000, fl_paste='yes')
sciframes.close()
GDISPLAY - Started: 2016-03-23T09:02:03
WARNING - GDISPLAY: only one science extension found, turning paste off
WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
Switching off the following flags:
ret_roi=no
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:02:11
GDISPLAY - Started: 2016-03-23T09:02:11
WARNING - GDISPLAY: only one science extension found, turning paste off
WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
Switching off the following flags:
ret_roi=no
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:02:13
GDISPLAY - Started: 2016-03-23T09:02:13
WARNING - GDISPLAY: only one science extension found, turning paste off
WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
Switching off the following flags:
ret_roi=no
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:02:14
GDISPLAY - Started: 2016-03-23T09:02:14
WARNING - GDISPLAY: only one science extension found, turning paste off
WARNING - GDISPLAY: Input image has been processed by GMOSAIC.
Switching off the following flags:
ret_roi=no
GDISPLAY: Displaying image...
GDISPLAY: Starting imexamine - quit with q
GDISPLAY - Finished: 2016-03-23T09:02:16
In [19]:
iraf.imdelete(science['procsci'], verify='no')
iraf.imcoadd('mrg//@sci.lis', outimage=science['procsci'], logfile=logfile)
-----------------------------------------------------------------------------
IMCOADD -- Wed Mar 23 09:02:21 HST 2016
Images (and masks) in list
mrgS20101101S0087
mrgS20101101S0088
mrgS20101101S0089
mrgS20101101S0090
Alignment method: wcs
Making individual masks from DQ
Fixpix input images using DQ
Instrument + camera : GMOS-S
Telescope ISS port number: 3
Bad pixel file gmos$data/gmos-s_bpm_EEV_22_3amp_v2_mosaic.pl
Read noise: 3.85 e- Gain: 1.00 e-/ADU
Co-added cleaned output image: ngc6872_r
Pixel scale: 0.1461
Statistics section: [100:3008,100:2204]
Finding objects in mrgS20101101S0087
Sampling for geotran : 20 20
Block size for geotran: 2048 2048
Centering objects in mrgS20101101S0087
Number of objects successfully centered in reference image 485
Alignment method: wcs
Coordinate list: imxymatch.1 Transform: mrgS20101101S0088_trn
Results file:
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.005628878 0.006281923
Coordinate mapping parameters
Mean Xref and Yref: 1554.5 1152.5
Mean Xin and Yin: 1520.273 1152.509
X and Y shift: -34.21567 4.210953E-4 (xin yin)
X and Y scale: 0.9999979 0.9999979 (xin / xref yin / yref)
X and Y axis rotation: 359.99960 359.99960 (degrees degrees)
Centering objects in mrgS20101101S0088
Xmax = 3108 Ymax = 2304
Entering geomap
Fitting geometry rscale
Iterating a maximum of 5 times
Coordinate list: mrgS20101101S0088_trn Transform: mrgS20101101S0088_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 12.51114 11.03705
Coordinate mapping parameters
Mean Xref and Yref: 1598.206 1147.867
Mean Xin and Yin: 1560.464 1148.434
X and Y shift: -34.61303 2.331055 (xin yin)
X and Y scale: 0.999429 0.999429 (xin / xref yin / yref)
X and Y axis rotation: 0.03191 0.03191 (degrees degrees)
Coordinate list: mrgS20101101S0088_trn Transform: mrgS20101101S0088_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 1.950519 1.863659
Coordinate mapping parameters
Mean Xref and Yref: 1606.994 1129.748
Mean Xin and Yin: 1572.842 1130.485
X and Y shift: -34.46964 0.02893066 (xin yin)
X and Y scale: 1.000062 1.000062 (xin / xref yin / yref)
X and Y axis rotation: 359.99524 359.99524 (degrees degrees)
Coordinate list: mrgS20101101S0088_trn Transform: mrgS20101101S0088_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.2339694 0.1931138
Coordinate mapping parameters
Mean Xref and Yref: 1634.163 1125.485
Mean Xin and Yin: 1599.887 1125.839
X and Y shift: -34.18585 0.3146973 (xin yin)
X and Y scale: 0.9999579 0.9999579 (xin / xref yin / yref)
X and Y axis rotation: 359.99829 359.99829 (degrees degrees)
Coordinate list: mrgS20101101S0088_trn Transform: mrgS20101101S0088_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.1612063 0.1291615
Coordinate mapping parameters
Mean Xref and Yref: 1636.073 1098.907
Mean Xin and Yin: 1601.792 1099.232
X and Y shift: -34.16214 0.3497314 (xin yin)
X and Y scale: 0.999943 0.999943 (xin / xref yin / yref)
X and Y axis rotation: 359.99863 359.99863 (degrees degrees)
Coordinate list: imxymatch.1 Transform: mrgS20101101S0089_trn
Results file:
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.005875335 0.006372866
Coordinate mapping parameters
Mean Xref and Yref: 1554.5 1152.5
Mean Xin and Yin: 1554.475 1118.263
X and Y shift: -0.02503858 -34.23117 (xin yin)
X and Y scale: 0.9999979 0.9999979 (xin / xref yin / yref)
X and Y axis rotation: 0.00014 0.00014 (degrees degrees)
Centering objects in mrgS20101101S0089
Xmax = 3108 Ymax = 2304
Entering geomap
Fitting geometry rscale
Iterating a maximum of 5 times
Coordinate list: mrgS20101101S0089_trn Transform: mrgS20101101S0089_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 13.30814 11.71999
Coordinate mapping parameters
Mean Xref and Yref: 1593.81 1151.595
Mean Xin and Yin: 1554.707 1118.522
X and Y shift: -35.34849 -32.28162 (xin yin)
X and Y scale: 0.9997742 0.9997742 (xin / xref yin / yref)
X and Y axis rotation: 0.01718 0.01718 (degrees degrees)
Coordinate list: mrgS20101101S0089_trn Transform: mrgS20101101S0089_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 2.767917 2.361728
Coordinate mapping parameters
Mean Xref and Yref: 1598.868 1147.175
Mean Xin and Yin: 1564.513 1113.621
X and Y shift: -34.93499 -34.30554 (xin yin)
X and Y scale: 1.000112 1.000112 (xin / xref yin / yref)
X and Y axis rotation: 0.00377 0.00377 (degrees degrees)
Coordinate list: mrgS20101101S0089_trn Transform: mrgS20101101S0089_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.2402036 0.1860939
Coordinate mapping parameters
Mean Xref and Yref: 1625.268 1142.746
Mean Xin and Yin: 1590.811 1108.582
X and Y shift: -34.34774 -34.17566 (xin yin)
X and Y scale: 0.9999845 0.9999845 (xin / xref yin / yref)
X and Y axis rotation: 359.99796 359.99796 (degrees degrees)
Coordinate list: mrgS20101101S0089_trn Transform: mrgS20101101S0089_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.146128 0.1150479
Coordinate mapping parameters
Mean Xref and Yref: 1632.61 1149.365
Mean Xin and Yin: 1598.217 1115.231
X and Y shift: -34.34995 -34.16431 (xin yin)
X and Y scale: 0.9999904 0.9999905 (xin / xref yin / yref)
X and Y axis rotation: 359.99838 359.99838 (degrees degrees)
Coordinate list: imxymatch.1 Transform: mrgS20101101S0090_trn
Results file:
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.005748001 0.006625012
Coordinate mapping parameters
Mean Xref and Yref: 1554.5 1152.5
Mean Xin and Yin: 1588.731 1152.488
X and Y shift: 34.24295 -0.02277369 (xin yin)
X and Y scale: 0.999998 0.999998 (xin / xref yin / yref)
X and Y axis rotation: 359.99954 359.99954 (degrees degrees)
Centering objects in mrgS20101101S0090
Xmax = 3108 Ymax = 2304
Entering geomap
Fitting geometry rscale
Iterating a maximum of 5 times
Coordinate list: mrgS20101101S0090_trn Transform: mrgS20101101S0090_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 12.83797 10.97163
Coordinate mapping parameters
Mean Xref and Yref: 1597.235 1147.954
Mean Xin and Yin: 1593.364 1114.441
X and Y shift: -0.9705968 -33.3031 (xin yin)
X and Y scale: 0.9999892 0.9999892 (xin / xref yin / yref)
X and Y axis rotation: 0.00946 0.00946 (degrees degrees)
Coordinate list: mrgS20101101S0090_trn Transform: mrgS20101101S0090_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 2.47633 1.576959
Coordinate mapping parameters
Mean Xref and Yref: 1599.015 1131.845
Mean Xin and Yin: 1599.307 1097.966
X and Y shift: -0.4821177 -34.57532 (xin yin)
X and Y scale: 1.000154 1.000154 (xin / xref yin / yref)
X and Y axis rotation: 359.99533 359.99533 (degrees degrees)
Coordinate list: mrgS20101101S0090_trn Transform: mrgS20101101S0090_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.1917848 0.157337
Coordinate mapping parameters
Mean Xref and Yref: 1622.835 1124.102
Mean Xin and Yin: 1622.594 1089.803
X and Y shift: -0.2125957 -34.31775 (xin yin)
X and Y scale: 1.000015 1.000015 (xin / xref yin / yref)
X and Y axis rotation: 359.99908 359.99908 (degrees degrees)
Coordinate list: mrgS20101101S0090_trn Transform: mrgS20101101S0090_trn
Results file: tmpres92482_2107
Coordinate mapping status
X fit ok. Y fit ok.
Xin and Yin fit rms: 0.1321171 0.1152363
Coordinate mapping parameters
Mean Xref and Yref: 1635.988 1120.243
Mean Xin and Yin: 1635.788 1085.969
X and Y shift: -0.2193163 -34.31055 (xin yin)
X and Y scale: 1.000019 1.000019 (xin / xref yin / yref)
X and Y axis rotation: 359.99954 359.99954 (degrees degrees)
Median sky level for mrgS20101101S0087: 3725.7
Transforming mrgS20101101S0088 to mrgS20101101S0088_trn
Transforming image mrgS20101101S0088[SCI] to image mrgS20101101S0088_trn
Using transform mrgS20101101S0088_trn in database imcoadd.dat
Median sky level for mrgS20101101S0088: 3673.0
Transforming mrgS20101101S0089 to mrgS20101101S0089_trn
Transforming image mrgS20101101S0089[SCI] to image mrgS20101101S0089_trn
Using transform mrgS20101101S0089_trn in database imcoadd.dat
Median sky level for mrgS20101101S0089: 3722.6
Transforming mrgS20101101S0090 to mrgS20101101S0090_trn
Transforming image mrgS20101101S0090[SCI] to image mrgS20101101S0090_trn
Using transform mrgS20101101S0090_trn in database imcoadd.dat
Median sky level for mrgS20101101S0090: 3675.0
Median sky level for mrgS20101101S0087: 3725.7
Median sky level for mrgS20101101S0088: 3673.0
Transforming mrgS20101101S0087badpix +mask to mrgS20101101S0088badpix
Transforming image tmpmsk92482_2101 to image mrgS20101101S0088badpix
Using transform mrgS20101101S0088_trn in database imcoadd.dat
Median sky level for mrgS20101101S0089: 3722.6
Transforming mrgS20101101S0087badpix +mask to mrgS20101101S0089badpix
Transforming image tmpmsk92482_2138 to image mrgS20101101S0089badpix
Using transform mrgS20101101S0089_trn in database imcoadd.dat
Median sky level for mrgS20101101S0090: 3675.0
Transforming mrgS20101101S0087badpix +mask to mrgS20101101S0090badpix
Transforming image tmpmsk92482_2142 to image mrgS20101101S0090badpix
Using transform mrgS20101101S0090_trn in database imcoadd.dat
Mean intensity for mrgS20101101S0087_trn: (N,absolute,relative)= 172 1950178 1.000
Mean intensity for mrgS20101101S0088_trn: (N,absolute,relative)= 172 1949907 1.000
Mean intensity for mrgS20101101S0089_trn: (N,absolute,relative)= 172 1948036 0.999
Mean intensity for mrgS20101101S0090_trn: (N,absolute,relative)= 172 1957529 1.004
Masking cosmic ray events in mrgS20101101S0087_trn
Applying growth radius 1
Masking cosmic ray events in mrgS20101101S0088_trn
Applying growth radius 1
Masking cosmic ray events in mrgS20101101S0089_trn
Applying growth radius 1
Masking cosmic ray events in mrgS20101101S0090_trn
Applying growth radius 1
IMCOADD exit status: good
-----------------------------------------------------------------------------
In [22]:
iraf.display(science['procsci']+'[sci,1]', 1, \
zrange='no', zscale='no', z1=1000, z2=20000)
shutil.copy(science['procsci'], products_path)
# Final reduced image
os.listdir(products_path)
z1=1000. z2=20000.
Out[22]:
['ngc6872_i.fits', 'ngc6872_r.fits']
In [ ]:
Content source: KathleenLabrie/reduxGemini
Similar notebooks: