from astropy.io import fits as _fits
import matplotlib.pyplot as _plt
#from matplotlib.pyplot import xticks as _xticks
import numpy as _np
from numpy import linalg as _LA
from scipy.optimize import curve_fit as _curve_fit
from scipy.special import gamma as _gamma
from scipy.stats import chi2 as _chi2

#============================================================================================================================
#plt.rc('text', usetex=True)
#plt.rc('font', family='serif')         Adds LaTeX Font, must have LaTeX installed
#plt.rcParams.update({'font.size':20})
#============================================================================================================================

def OpenFits(file,contrast=2**16-1,BoW=1):
    """Open fit/fits file. Example OpenFits('Atlas.fits')
        d is an array of the size of the CCD eg. 1020x1530"""
    f = _fits.open(file)
    d = f[0].data
    if BoW == 0:
        _plt.imshow(d, vmax=contrast, cmap='Greys', origin='lowerleft')
    if BoW == 1:
        _plt.imshow(d, vmax=contrast, cmap='gray', origin='lowerleft')
    _plt.show()
    return d