In [1]:

    

import subprocess
from astropy.io import fits
import healpy as hp
from astropy.utils.data import download_file
from astropy.utils.data import Conf 
Conf.remote_timeout = 10000 ## In case there's a slow connection
%pylab
%matplotlib inline


    

    




Using matplotlib backend: MacOSX
Populating the interactive namespace from numpy and matplotlib

In [2]:

    

##Planck first

#mapfiles = [
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/LFI_SkyMap_030_1024_R2.01_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/LFI_SkyMap_044_1024_R2.01_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/LFI_SkyMap_070_1024_R2.01_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/LFI_SkyMap_070_2048_R2.01_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/HFI_SkyMap_100_2048_R2.02_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/HFI_SkyMap_143_2048_R2.02_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/HFI_SkyMap_217_2048_R2.02_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/HFI_SkyMap_353_2048_R2.02_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/HFI_SkyMap_545_2048_R2.02_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/HFI_SkyMap_857_2048_R2.02_full.fits",
#"http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/maps/component-maps/cmb/COM_CMB_IQU-commander-field-Int_2048_R2.01_full.fits"
#]
#maps = [download_file(file, cache = True) for file in mapfiles]
#healmaps = [hp.read_map(pixelmap, memmap = True, h = False) for pixelmap in maps]


    

In [3]:

    

### wmap
mapfiles = [
    'http://lambda.gsfc.nasa.gov/data/map/dr5/skymaps/9yr/raw/wmap_band_imap_r9_9yr_K_v5.fits',
    #'http://lambda.gsfc.nasa.gov/data/map/dr5/skymaps/9yr/raw/wmap_band_imap_r9_9yr_Ka_v5.fits',
    #'http://lambda.gsfc.nasa.gov/data/map/dr5/skymaps/9yr/raw/wmap_band_imap_r9_9yr_Q_v5.fits',
    #'http://lambda.gsfc.nasa.gov/data/map/dr5/skymaps/9yr/raw/wmap_band_imap_r9_9yr_V_v5.fits',
    #'http://lambda.gsfc.nasa.gov/data/map/dr5/skymaps/9yr/raw/wmap_band_imap_r9_9yr_W_v5.fits'
]
maps = [download_file(file, cache = True) for file in mapfiles]
healmaps = [hp.read_map(pixelmap, field = (0,1), memmap = True, h = False) for pixelmap in maps]


    

    




NSIDE = 512
ORDERING = NESTED in fits file
INDXSCHM = IMPLICIT






    




/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/healpy/fitsfunc.py:335: UserWarning: No INDXSCHM keyword in header file : assume IMPLICIT
  "assume {}".format(schm))






    




Ordering converted to RING
Ordering converted to RING

In [9]:

    

hp.mollview(healmaps[0][0], xsize = 800, title = 'temperature');

hp.mollview(healmaps[0][1], xsize = 800, title = 'Nobs');


    

In [ ]: