Healpix mollview outlines

Use of Cobaya. camb, CLASS, cosmomc, compilers, etc.
Post Reply
Samuel Leach
Posts: 16
Joined: October 15 2004
Affiliation: SISSA, Trieste
Contact:

Healpix mollview outlines

Post by Samuel Leach » January 24 2007

I want to make a Healpix mollview plot overlaying several "outline" structures, something similar to this:

http://www.astro.caltech.edu/~lgg/bicep/fig3.jpg

I know how to insert a single outline structure using the outline= option, but I haven't figured out to insert several outlines at the same time.

Can someobody help me do that ?

Thanks in advance.

Ben Gold
Posts: 81
Joined: September 25 2004
Affiliation: University of Minnesota
Contact:

Healpix mollview outlines

Post by Ben Gold » January 25 2007

Mollview should accept multiple outlines if you make them as separate structures and give them to mollview using the "outline=" option. But if you have trouble with that, I think you can just directly call "outline_coord2uv", which is the routine that actually draws the outlines. You might have to look at the source code to get all the arguments right.

Samuel Leach
Posts: 16
Joined: October 15 2004
Affiliation: SISSA, Trieste
Contact:

Healpix mollview outlines

Post by Samuel Leach » January 25 2007

Thanks Ben, you are right that outline_coord2uv can be used. Nicolas Ponthieu provided me with this example code to do that:

Code: Select all

boomerang_ra=[72,72,92,92,72]
boomerang_dec=[-40,-50,-50,-40,-40]
boomerang_area= {coord:'C',ra:boomerang_ra,dec:boomerang_dec,line:0}

rot_ang=[0,-90,0]
eul_mat = euler_matrix_new(rot_ang(0), -rot_ang(1), rot_ang(2), /Deg, /ZYX)
mollview,'dust.fits',graticule=[45,20],rot=rot_ang,min=0,max=30,units='!4l!3K',$
         title='Dust at 150GHz',factor=1e6
outline_coord2uv,boomerang_area,coord_out,eul_mat,projection='moll',/show,thick=3,col=0

Post Reply