Copyright © Michael Richmond.
This work is licensed under a Creative Commons License.
Your job today is to create master dark frames and master flatfield frames for the images we collected March 30, 2002. You will use these master frames to reduce the images of 1107 Lictoria.
First, you must combine the 1-second dark frames to create a master dark suitable for subtracting from the 1-second flatfield frames. The raw images are called dark1.001, dark1.002, etc.
chmod 644 dark1.*
box 1 sr=100 sc=100 nr=50 nc=50
box 1 show
abx dark1.001 box=1
Now, combine all the 1-second dark frames via the median command to create a master dark frame, and compare it to one of the individual dark frames.
median dark1.001 dark1.002 ... dark1.010 outfile=dark1.fts nomean verbose
tv dark1.001 box=1 zoom=5 z=150 l=100
tv dark1.fts box=1 zoom=5 z=150 l=100
What do you see?
Write down the differences and similarities of the
individual and master dark frames.
Okay, you have a master dark frame for the 1-second exposures. Make sure it is called dark1.fts.
Now, repeat the procedure to create a master dark frame suitable for the 15-second target exposures we took later in the night. You need to combine the dark frames called dark15.001, dark15.002, and so on, to create a master dark frame called dark15.fts. As before, calculate the mean and stdev of pixel values in box 1 (the same box) of one individual 15-second dark frame, and the master 15-second dark frame. Compare the individual and master darks like so:
tv dark15.001 box=1 zoom=5 z=150 l=200
tv dark15.fts box=1 zoom=5 z=150 l=200
Describe in your own words the difference between these two frames.
Let's try to calculate how many "counts" are generated by the dark current each second. Find an area in the 15-second master dark frame which is free from hot pixels, and exactly 10 x 10 pixels in size. Use the cursor to pick the coordinates of the upper left-hand corner -- say you choose row=116, col=121. Then
box 2 sr=116 sc=121 nr=10 nc=10
We have a series of 1-second, V-band images of the twilight sky, called flatv.001, flatv.002, etc. You must subtract the dark current from each frame, and the combine them via the median command to create a master flatfield frame.
/bin/rm filename
sub flatv.001 dark1.fts
sub flatv.002 dark1.fts
etc.
median flatv.001 flatv.002 ... outfile=flatv.fts verbose
Note that the nomean option does NOT appear this time.
The median command will print out some stuff to the screen.
Explain what it is doing, and why.
Copyright © Michael Richmond.
This work is licensed under a Creative Commons License.