source: Sophya/trunk/SophyaPI/ProgPI/skymapmodule.pic@ 2175

#----------------------
shell rm -f h4*.fits t4*.ppf h8*.fits t8*.ppf clh*.ppf clt*.ppf almh*.ppf almt*.ppf
delobjs *
c++import SkyMap
c++exec \
  double theta,phi; \
  SphereHEALPix<r_4>  h4(256); \
  for(int_4 i=0;i<h4.NbPixels();i++) { \
    h4.PixThetaPhi(i,theta,phi); \
    h4(i) = 0.; \
    double v = (M_PI_2-theta)*(M_PI_2-theta)/(0.50*0.50); \
    v += (phi-M_PI)*(phi-M_PI)/(M_PI*M_PI/2.); \
    v = 100.*exp(-v/2.); \
    v += 10.*NorRand(); \
    if(phi<20.*M_PI/180.||phi>340.*M_PI/180.) continue; \
    if(phi<190.*M_PI/180.&&phi>170.*M_PI/180.) continue; \
    h4(i) = v; \
  } \
  KeepObj(h4); \
  SphereThetaPhi<r_4> t4(400); \
  for(int_4 i=0;i<t4.NbPixels();i++) { \
    t4.PixThetaPhi(i,theta,phi); \
    t4(i) = 0.; \
    double v = (M_PI_2-theta)*(M_PI_2-theta)/(0.50*0.50); \
    v += (phi-M_PI)*(phi-M_PI)/(M_PI*M_PI/2.); \
    v = 100.*exp(-v/2.); \
    v += 10.*NorRand(); \
    if(phi<20.*M_PI/180.||phi>340.*M_PI/180.) continue; \
    if(phi<190.*M_PI/180.&&phi>170.*M_PI/180.) continue; \
    t4(i) = v; \
  } \
  KeepObj(t4); \
  cout<<"End of Job"<<endl;

savefits h4 !h4.fits
saveppf t4 t4.ppf

zone 1 2
disp h4 zoom/4
disp t4 zoom/4
zone

#----------------------
delobjs *
openfits h4.fits

map2double h4
rename h4 h8
savefits h8 !h8.fits

delobjs *
openppf t4.ppf
map2double t4
rename t4 t8
saveppf t8 t8.ppf

#----------------------
delobjs *
openfits h4.fits
openfits h8.fits
openppf t4.ppf
openppf t8.ppf
typemap h4
typemap h8
typemap t4
typemap t8

#----------------------
delobjs *
openfits h8.fits
openppf t8.ppf
map2map h8 t
map2map t8 h

#----------------------
delobjs *
openfits h8.fits
openppf t8.ppf
mapmult h8 1000
mapmult t8 1000

#----------------------
delobjs *
openfits h8.fits
maproj h8 h8p
mapop h8p - h8
disp h8p zoom/4

maproj h8 h8p64 64
disp h8p64 zoom/4
savefits h8p64 !h8p64.fits

delobjs *
openppf t8.ppf
maproj t8 t8p
mapop t8p - t8
disp t8p zoom/4

maproj t8 t8p100 100
disp t8p100 zoom/4
saveppf t8p100 t8p100.ppf

#----------------------
delobjs *
openfits h8p64.fits
map2cl h8p64 clh64
saveppf clh64 clh64.ppf
n/plot clh64.n*(n+1)*val%n ! ! crossmarker5

map2cl h8p64 clh64 50
n/plot clh64.n*(n+1)*val%n ! ! crossmarker5

delobjs *
openppf t8p100.ppf
map2cl t8p100 clt100 191
saveppf clt100 clt100.ppf
n/plot clt100.n*(n+1)*val%n ! ! crossmarker5

#----------------------
delobjs *
openfits h8p64.fits
map2alm h8p64 almh64
saveppf almh64 almh64.ppf
disp almh64

delobjs *
openppf t8p100.ppf
map2alm t8p100 almt100 191
saveppf almt100 almt100.ppf
disp almt100

#----------------------
delobjs *
openppf clh64.ppf
settypemap h
cl2map clh64 h8p64frcl 64
savefits h8p64frcl !h8p64frcl.fits
disp h8p64frcl

delobjs *
openppf clt100.ppf
settypemap t
cl2map clt100 t8p100frcl 100
saveppf t8p100frcl t8p100frcl.ppf
disp t8p100frcl

#----------------------
delobjs *
openppf almh64.ppf
settypemap h
alm2map almh64 h8p64fralm 64
savefits h8p64fralm !h8p64fralm.fits
disp h8p64fralm

delobjs *
openppf almt100.ppf
settypemap t
alm2map almt100 t8p100fralm 100
saveppf t8p100fralm t8p100fralm.ppf
disp t8p100fralm

#----------------------
delobjs *
settypemap h
crmapmask h8m 256 -20,20 100,130 0,1
savefits h8m !h8m.fits
disp h8m zoom/4

delobjs *
settypemap t
crmapmask t8m 400 -20,20 100,130 0,1
saveppf t8m t8m.ppf
disp t8m zoom/4

#----------------------
delobjs *
openfits h8.fits
crmaskfrmap h8fm 256 h8 -1.e-30,1.e-30 0,1
savefits h8fm !h8fm.fits
disp h8fm zoom/4

delobjs *
openppf t8.ppf
crmaskfrmap t8fm 400 t8 -1.e-30,1.e-30 0,1
saveppf t8fm t8fm.ppf
disp t8fm zoom/4

#----------------------
delobjs *
openfits h8.fits
openfits h8m.fits
maskmap h8 h8m
disp h8 zoom/4

delobjs *
openppf t8.ppf
openppf t8m.ppf
maskmap t8 t8m
disp t8 zoom/4

delobjs *
openfits h8.fits
openppf t8m.ppf
maskmap h8 t8m
disp h8 zoom/4

delobjs *
openppf t8.ppf
openfits h8m.fits
maskmap t8 h8m
disp t8 zoom/4

#----------------------
delobjs *
openfits h8m.fits
mapcover h8m 0.9,1 couvh8m
echo $couvh8m

delobjs *
openppf t8m.ppf
mapcover t8m 0.9,1 couvt8m
echo $couvt8m

#----------------------
delobjs *
openfits h8.fits
map2local h8 h8loc 200,300 20,30 0,90 !
disp h8loc
map2local h8 h8loc 200,300 20,30 0,90 ! 45
disp h8loc

delobjs *
openppf t8.ppf
map2local t8 t8loc 200,300 20,30 0,90 !
disp t8loc
map2local t8 t8loc 200,300 20,30 0,90 ! 45
disp t8loc

#----------------------
delobjs *
openfits h8.fits
rename h8 h8save
openfits h8.fits
mapop h8 + h8save - h8save * h8save / h8save
disp h8 zoom/4

openppf t8.ppf
rename t8 t8save
openppf t8.ppf
mapop t8 + t8save - t8save * t8save / t8save
disp t8 zoom/4

mapop h8 - t8save
disp h8 zoom/4
mapop t8 - h8save
disp t8 zoom/4

#----------------------
delobjs *
openfits h8.fits
openfits h8fm.fits
mapstat h8 !    mean sig
echo mean=$mean sig=$sig
mapstat h8 h8fm mean sig
echo mean=$mean sig=$sig

openppf t8.ppf
openppf t8fm.ppf
mapstat t8 !    mean sig
echo mean=$mean sig=$sig
mapstat t8 t8fm mean sig
echo mean=$mean sig=$sig