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

Question

使用 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")

输出

有没有办法只显示地图的某些区域？例如，只有上半球，还是只有左侧？

我的想法是只查看一小块天空以查看小点源，或者类似于来自 LSST"half-sky" 的投影

Answer 1

你可以使用mask，它是一个相同大小的boolean map，其中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)

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

How can I partially plot a Healpix map using Healpy?

python

healpy