如何使用Healpy部分绘制Healpix地图?

时间:2016-01-06 18:08:09

标签: python healpy

使用healpy时,我可以使用

在Mollview中绘制Healpix地图
import healpy
map = 'filename.fits'
healpy.visufunc.mollview(map)

或在教程中

>>> import numpy as np
>>> import healpy as hp
>>> NSIDE = 32
>>> m = np.arange(hp.nside2npix(NSIDE))
>>> hp.mollview(m, title="Mollview image RING")

输出

enter image description here

有没有办法只显示地图的某些区域?例如,只有上半球,或只有左侧?

我的想法是只观看天空中的小片以观察小点源,或者像半天空"来自LSST的投影

enter image description here

1 个答案:

答案 0 :(得分:2)

你可以使用一个掩码,它是一个相同大小的布尔映射,其中1被掩盖,0被掩盖:

http://healpy.readthedocs.org/en/latest/tutorial.html#masked-map-partial-maps

示例:

import numpy as np
import healpy as hp
NSIDE = 32
m = hp.ma(np.arange(hp.nside2npix(NSIDE), dtype=np.double))
mask = np.zeros(hp.nside2npix(NSIDE), dtype=np.bool)
pixel_theta, pixel_phi = hp.pix2ang(NSIDE, np.arange(hp.nside2npix(NSIDE)))
mask[pixel_theta > np.pi/2] = 1
m.mask = mask
hp.mollview(m)

enter image description here