In [2]:
%matplotlib inline 
import healpy as hp
import numpy as np

This crashes the kernel:

Bad pixels as NaN


In [ ]:
# Import data from FITS file
irc9 = hp.read_map('/work1/users/aaronb/Databrary/HEALPix/AKARI_HEALPix_orig/4096_nside/akari_9_4096_unseen.fits', 
                   nest=False)

# Replace UNSEEN with NaNs:

irc9[irc9==hp.UNSEEN] = np.nan

# Do the smoothing
irc9_smooth = hp.sphtfunc.smoothing(irc9, 
                                    fwhm = 0.00285, 
                                    iter=1)
hp.mollview(irc9_smooth,
           norm='hist',
           nest=False)

This doesn't:

Bad pixels as UNSEEN (-1.6375e+30)


In [1]:
# Import data from FITS file
irc9 = hp.read_map('/work1/users/aaronb/Databrary/HEALPix/AKARI_HEALPix_orig/4096_nside/akari_9_4096_unseen.fits', 
                   nest=False)

# Replace UNSEEN with NaNs:

irc9[irc9==hp.UNSEEN] = np.nan

# Do the smoothing
irc9_smooth = hp.sphtfunc.smoothing(irc9, 
                                    fwhm = 0.00285, 
                                    iter=1)
hp.mollview(irc9_smooth,
           norm='hist',
           nest=False)


NSIDE = 4096
ORDERING = RING in fits file
INDXSCHM = IMPLICIT
Sigma is 4.160648 arcmin (0.001210 rad) 
-> fwhm is 9.797578 arcmin
Sigma is 0.000000 arcmin (0.000000 rad) 
-> fwhm is 0.000000 arcmin
/work1/users/aaronb/Softbrary/Anaconda/lib/python2.7/site-packages/healpy/projaxes.py:859: MaskedArrayFutureWarning: setting an item on a masked array which has a shared mask will not copy the mask and also change the original mask array in the future.
Check the NumPy 1.11 release notes for more information.
  result[np.isinf(val.data)] = -np.inf