Python get_input_mapsの例

コード例 #1

ファイル: smallfield_0016.py プロジェクト: rwogburn/smallfield

#!/bin/env python

import smallfield
import healmap
import healpy as hp
import numpy as np

if __name__ == "__main__":
  use_split = 'D'
  use_nside = 512
  out_nside = 32
  outfname = 'smallfield_0016.fits'
  dl_lcdm,lc = smallfield.calc_lcdm_dl()
  print dl_lcdm

  h1,h2 = healmap.get_input_maps(use_split,use_nside=use_nside,dir='/n/bicepfs3/general/input_maps/planck_pr2/')
  lat,lon,ipix = smallfield.gen_field_centers(nside=out_nside,minlat=-1.0)
  i0 = 0
  if i0==0:
    rb = np.zeros([12*out_nside*out_nside])
    sb = 0.0 * rb
    re = 0.0 * rb
    se = 0.0 * rb
  else:
    rb = hp.read_map(outfname,field=0)
    sb = hp.read_map(outfname,field=1)
    re = hp.read_map(outfname,field=2)
    se = hp.read_map(outfname,field=3)

  for i in xrange(i0,lat.size):
    if (rb[ipix[i]]!=0.) or (sb[ipix[i]]!=0.):

コード例 #2

ファイル: splitcomp.py プロジェクト: rwogburn/smallfield

#!/opt/local/bin/python2.7

import smallfield
import healmap
import healpy as hp
import numpy as np

if __name__ == "__main__":
  use_split = 'D'
  use_nside = 512

  h1,h2 = healmap.get_input_maps(use_split,use_nside=use_nside,dir='..')
  sfac = 0.0395
  squ = 10
  st = 200
  rot=[0,90,0]
  #
  hp.orthview(h1.map[0,:]*sfac,rot=rot,nest=False,title="Split 1 T",min=0,max=st)
  hp.orthview(h1.map[1,:]*sfac,rot=rot,nest=False,title="Split 1 Q",min=-squ,max=squ)
  hp.orthview(h1.map[2,:]*sfac,rot=rot,nest=False,title="Split 1 U",min=-squ,max=squ)
  #
  hp.orthview(h2.map[0,:]*sfac,rot=rot,nest=False,title="Split 2 T",min=0,max=st)
  hp.orthview(h2.map[1,:]*sfac,rot=rot,nest=False,title="Split 2 Q",min=-squ,max=squ)
  hp.orthview(h2.map[2,:]*sfac,rot=rot,nest=False,title="Split 2 U",min=-squ,max=squ)
  #
  plt.show()