Description
Computes the average of a set of processed and aligned images.
Usage
stack.py input-file-list-file-name output-file-name
The input file list is a text file of the paths of images to stack. Output file is another FITS file.
Examples
The file 'brgmlist1.txt' contains the following entries:
072407/processed/im0101.a.fits 072407/processed/im0102.a.fits 072407/processed/im0103.a.fits 072407/processed/im0104.a.fits 072407/processed/im0105.a.fits 072407/processed/im0106.a.fits 072407/processed/im0107.a.fits 072407/processed/im0108.a.fits 072407/processed/im0109.a.fits 072407/processed/im0110.a.fits 072407/processed/im0111.a.fits 072407/processed/im0112.a.fits 072407/processed/im0113.a.fits 072407/processed/im0114.a.fits 072407/processed/im0115.a.fits 072407/processed/im0116.a.fits 072407/processed/im0117.a.fits 072407/processed/im0118.a.fits 072407/processed/im0119.a.fits 072407/processed/im0120.a.fits
Calling 'stack.py' creates the file 'brgm1.fits':
2: ~/Projects/Python/AstroPy > stack.py 072407/intermediates/brgmlist1.txt 072407/intermediates/brgm1.fits reading 072407/processed/im0101.a.fits reading 072407/processed/im0102.a.fits reading 072407/processed/im0103.a.fits reading 072407/processed/im0104.a.fits reading 072407/processed/im0105.a.fits reading 072407/processed/im0106.a.fits reading 072407/processed/im0107.a.fits reading 072407/processed/im0108.a.fits reading 072407/processed/im0109.a.fits reading 072407/processed/im0110.a.fits reading 072407/processed/im0111.a.fits reading 072407/processed/im0112.a.fits reading 072407/processed/im0113.a.fits reading 072407/processed/im0114.a.fits reading 072407/processed/im0115.a.fits reading 072407/processed/im0116.a.fits reading 072407/processed/im0117.a.fits reading 072407/processed/im0118.a.fits reading 072407/processed/im0119.a.fits reading 072407/processed/im0120.a.fits stack is (20, 512, 512) av is (512, 512), mean = 968.938119635, std = 1412.00217144 writing 072407/intermediates/brgm1.fits
The file 'brgmlist2.txt' contains the following entries:
072407/processed/im0277.a.fits 072407/processed/im0279.a.fits 072407/processed/im0280.a.fits 072407/processed/im0281.a.fits 072407/processed/im0282.a.fits 072407/processed/im0283.a.fits 072407/processed/im0284.a.fits 072407/processed/im0285.a.fits 072407/processed/im0286.a.fits 072407/processed/im0287.a.fits 072407/processed/im0288.a.fits 072407/processed/im0289.a.fits 072407/processed/im0290.a.fits 072407/processed/im0291.a.fits 072407/processed/im0292.a.fits 072407/processed/im0293.a.fits 072407/processed/im0294.a.fits 072407/processed/im0295.a.fits 072407/processed/im0296.a.fits
Image 278 was not included, as it is a multiple image taken while the planet is moving in the frame. Calling 'stack.py' creates the file 'brgm2.fits':
5: ~/Projects/Python/AstroPy > stack.py 072407/intermediates/brgmlist2.txt 072407/intermediates/brgm2.fits reading 072407/processed/im0277.a.fits reading 072407/processed/im0279.a.fits reading 072407/processed/im0280.a.fits reading 072407/processed/im0281.a.fits reading 072407/processed/im0282.a.fits reading 072407/processed/im0283.a.fits reading 072407/processed/im0284.a.fits reading 072407/processed/im0285.a.fits reading 072407/processed/im0286.a.fits reading 072407/processed/im0287.a.fits reading 072407/processed/im0288.a.fits reading 072407/processed/im0289.a.fits reading 072407/processed/im0290.a.fits reading 072407/processed/im0291.a.fits reading 072407/processed/im0292.a.fits reading 072407/processed/im0293.a.fits reading 072407/processed/im0294.a.fits reading 072407/processed/im0295.a.fits reading 072407/processed/im0296.a.fits stack is (19, 512, 512) av is (512, 512), mean = 954.750451992, std = 1127.93224344 writing 072407/intermediates/brgm2.fits
Here are the resulting stacked images ('brgm1.fits' = left, 'brgm2.fits' = right) displayed at scales of 0-2000 and 0-10000:
İSky Coyote 2007