import matplotlib.pyplot as plt 
from astropy.io import fits 
from astropy.wcs import WCS 
from astropy import units as u 
from astropy.visualization.wcsaxes.frame import EllipticalFrame
from reproject import reproject_interp
from spectral_cube import SpectralCube

fn = '/srv/data/dragons/dragons-canfar-Nov2025/frequency_cubes/dragons_freqK_gal_conv_i.car.fits'
freqIdat = fits.getdata(fn)
freqIhdr = fits.getheader(fn)

indat = fn

freqIhdr

SIMPLE  =                    T / conforms to FITS standard                      
BITPIX  =                  -64 / array data type                                
NAXIS   =                    3 / number of array dimensions                     
NAXIS1  =                  720                                                  
NAXIS2  =                  360                                                  
NAXIS3  =                  680                                                  
EXTEND  =                    T                                                  
DATE    = '2025-11-14'                                                          
INSTRUME= 'DRAO-15 '                                                            
EPOCH   = '2000    '                                                            
CDELT1  =                 -0.5                                                  
CRPIX1  =                360.5                                                  
CRVAL1  =                    0                                                  
CTYPE1  = 'GLON-CAR'                                                            
CUNIT1  = 'deg     '                                                            
CDELT2  =                  0.5                                                  
CRPIX2  =                180.5                                                  
CRVAL2  =                    0                                                  
CTYPE2  = 'GLAT-CAR'                                                            
CUNIT2  = 'deg     '                                                            
CDELT3  =                  1.0                                                  
CRPIX3  =                    1                                                  
CRVAL3  =           350.520833                                                  
CTYPE3  = 'FREQ    '                                                            
CUNIT3  = 'MHz     '                                                            
BUNIT   = 'K       '                                                            
BMEXTENS=                    T                                                  
CASAMBM =                    T                                                  

#need a target header,