source: corel/pre.py @ 2

#!/bin/env python
# coding: utf-8

VERSION='1.0'

import pyfits
from optparse import OptionParser

import numpy


parser = OptionParser()

(options, args) = parser.parse_args ()

if (len(args) !=2) :
    print "Usage : pre.py input.fits output.fits"

# Ouverture de l'image fits en argument
hdulist = pyfits.open(args[0])
scidata = hdulist[0].data

# suppression prescan/overscan/bord de plage

gauche = scidata[:2048,50:1074]
droite = scidata[:2048,1087:2111]

scidata = numpy.hstack((gauche,droite))
dim = scidata.shape

# Remplacement bete des mauvaises colonnes du E0
# PremiÚre colonne indiquée, largeur 4 colonnes

BadColumns = [92,119,431,715,842,1144,1160,1511]
for c in BadColumns:
    scidata[:,c:c+4]=scidata[:,c-5:c-1]

# Suppresion du fond (version brutale)

scidata -= 7000
scidata = numpy.where(scidata<0,0,scidata)


prihdr = hdulist[0].header
prihdr.add_history('Suppression os/ps')
prihdr.add_history('Interpolation mauvaises colonnes')
hdulist[0].data = scidata
hdulist.writeto(args[1],clobber=True)




hdulist.close()