#!/usr/bin/env python
import sys,pyfits
from equalize import histeq
from kapteyn import maputils
from matplotlib import pyplot as plt
from scipy.ndimage.filters import gaussian_filter
import warnings
warnings.simplefilter("ignore",DeprecationWarning)
if __name__ == "__main__":
HST_img = pyfits.open(sys.argv[1])
HST_hdr = HST_img[1].header
filtered_data,_ = histeq(HST_img[1].data + gaussian_filter(HST_img[1].data,2))
f = maputils.FITSimage(externalheader=HST_hdr,externaldata=filtered_data)
mplim = f.Annotatedimage(cmap="bone")
mplim.Image()
mplim.Graticule()
mplim.plot()
plt.savefig(sys.argv[2])
HST_img.close()
#!/usr/bin/env python
import sys, pyfits
from equalize import histeq
from kapteyn import maputils
from matplotlib import pyplot as plt
# from scipy.ndimage.filters import gaussian_filter
import warnings
warnings.simplefilter("ignore", DeprecationWarning)
if __name__ == "__main__":
HST_img = pyfits.open(sys.argv[1])
HST_hdr = HST_img["PRIMARY"].header
data, _ = histeq(HST_img["PRIMARY"].data)
f = maputils.FITSimage(externalheader=HST_hdr, externaldata=data)
mplim = f.Annotatedimage(cmap="bone")
mplim.Image()
mplim.Graticule()
mplim.plot()
plt.savefig(sys.argv[2], dpi=300)
HST_img.close()