I have a FITS file named 'my_cube.fits' with WCS. The file has spatial information on axis 1 and 2 (X and Y) and spectral information on axis 3 (Z). When I load it using astropy.io.fits, the spectr...

I write script for calibration of image (dark frame and flat field)...Here is part of code for n in range(len(img)): with pyfits.open(img[n], mode='update', memmap=True) as im: imgg = im[0].data...

I've looked and looked for a solution to this problem and am turning up nothing. I'm generating rectangular FITS images through matplotlib and subsequently applying WCS coordinates to them using ...