if not 'fe' in locals(): import pyfits #fe_NS_mosaic.fits H2_NS_mosaic.fits H2_mosaic.fits KS_mosaic.fits Ks_NS_mosaic.fits # fe_montage_georgs_southeast.fits # altair_big_h2.fits # fe_montage_georgs_southeast.fits fef = pyfits.open('altair_big_fe.fits') fe = fef[0].data if os.path.exists('altair_big_fe_unsharp.fits'): feunsharp = pyfits.getdata('altair_big_fe_unsharp.fits') fenormed = pyfits.getdata('altair_big_fe_normed.fits') else: #kernel = astropy.nddata.convolution.make_kernel.make_kernel([151,151],kernelwidth=50) #fesmooth = astropy.nddata.convolve(fe,kernel) fesmooth = agpy.smooth(mask_image(fe,downsample=4),250, fft_pad=False, interpolate_nan=True, psf_pad=False, ignore_edge_zeros=True, normalize_kernel=True, use_numpy_fft=True, nthreads=1, use_rfft=True, complextype=np.float32, silent=False, boundary='fill') feunsharp = fe for ii in xrange(4): for jj in xrange(4): shape = feunsharp[ii::4,jj::4].shape feunsharp[ii::4,jj::4] -= fesmooth[:shape[0],:shape[1]] fef[0].data = feunsharp fef.writeto("altair_big_fe_unsharp.fits",clobber=True) fenormed = fe for ii in xrange(4): for jj in xrange(4): shape = fenormed[ii::4,jj::4].shape
import agpy import numpy as np import PIL, ImageEnhance prefix = "" if not 'fe' in locals(): import pyfits fef = pyfits.open('big_mosaic_feii.fits') fe = fef[0].data if os.path.exists('big_mosaic_feii_unsharp.fits'): feunsharp = pyfits.getdata('big_mosaic_feii_unsharp.fits') fenormed = pyfits.getdata('big_mosaic_feii_normed.fits') else: #kernel = astropy.nddata.convolution.make_kernel.make_kernel([151,151],kernelwidth=50) #fesmooth = astropy.nddata.convolve(fe,kernel) mimg = mask_image(fe, downsample=4) fesmooth = agpy.smooth(mimg, 250, fft_pad=False, interpolate_nan=True, psf_pad=False, ignore_edge_zeros=True, normalize_kernel=True, use_numpy_fft=True, nthreads=1, use_rfft=True, complextype=np.float32, silent=False, boundary='fill') feunsharp = fe for ii in xrange(4):
import FITS_tools prefix="" if not 'fe' in locals(): fef = pyfits.open('big_mosaic_feii.fits') header = fef[0].header fe = fef[0].data avm = pyavm.AVM.from_header(header) mywcs = avm.to_wcs() if os.path.exists('big_mosaic_feii_unsharp.fits'): feunsharp = pyfits.getdata('big_mosaic_feii_unsharp.fits') fenormed = pyfits.getdata('big_mosaic_feii_normed.fits') else: #kernel = astropy.nddata.convolution.make_kernel.make_kernel([151,151],kernelwidth=50) #fesmooth = astropy.nddata.convolve(fe,kernel) mimg = mask_image(fe,downsample=4) fesmooth = convolve_fft(mimg, Gaussian2DKernel(250), fft_pad=False, interpolate_nan=True, psf_pad=False, ignore_edge_zeros=True, normalize_kernel=True, use_numpy_fft=True, nthreads=1, use_rfft=True, complextype=np.float32, silent=False, boundary='fill') feunsharp = fe for ii in range(4): for jj in range(4): shape = feunsharp[ii::4,jj::4].shape feunsharp[ii::4,jj::4] -= fesmooth[:shape[0],:shape[1]] fef[0].data = feunsharp fef.writeto("big_mosaic_feii_unsharp.fits",clobber=True) fenormed = fe for ii in range(4):