CODEWITHSUNDEEP
plothealpy
bluebird343
bluebird343

Reputation: 143

Healpy Plotting

I'm very new to this, but I have been given a .fits file that I've plotted with healpy using this code.

map = hp.read_map('.fits')
hp.mollview(map, title=fileName,min=-1,max=1)
hp.graticule()
plt.show()

I'd like to zoom into the coordinates RA = [-50,50] and Dec = [-70,-45].

How can I do this? I only see the options for rotating the map, not zooming in on an area. Am I misunderstanding something?

Thank you for your time!

Upvotes: 1

Views: 2248

Answers (1)

Andrea Zonca
Andrea Zonca

Reputation: 8773

You cannot zoom in Mollweide projection, you should use a Gnomonic projection instead:

hp.gnomview(map, rot=[colatitude_deg, longitude_deg])

you can change xsize and reso to change the width of the projected area.

Upvotes: 1

Related Questions