# vypsani spektra ze snimku
# spoustet v terminalu jako:
#   $ python slice.py > tabulka
# a vysledna tabulka se ulozi

import pyfits

# tato funkce si bere jako argument nazev snimku
def slice(filename):

    # otevreme jej
    f = pyfits.open(filename)

    # vytahneme si z nej obrazek
    obr = f[0].data
 
    # zjistme jeho rozmery
    m,n = obr.shape

    # zvolime si souradnice (ja mam 2100)
    for i in range(0,m):
       s = sum(obr[i,2100:2200])/100
       print("{0} {1} {2}".format(i,obr[i,2100],s))

    # a vypsana tabulka obsahuje sloupecky:
    # pixel     hodnota v obr[i,2100]   prumerna hodnota obr[i,2100-2200]


slice('IMG_3480.fits')