Healpy如何在每个像素上写标号

如何在Healpy天图上标上像素的坐标，示例如下：

import matplotlib
import healpy as hp
matplotlib.rcParams.update({'font.size': 8})
nside = 4
mmap = np.arange(hp.nside2npix(nside))
hp.mollview(mmap, title='', cbar=None)
for i in range(len(testmap)):
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi, '%d'%i)
plt.show()

效果图：

orthview视图：

import matplotlib
import healpy as hp
matplotlib.rcParams.update({'font.size': 8})
nside = 4
mmap = np.arange(hp.nside2npix(nside))
hp.orthview(mmap, title='', cbar=None)
for i in range(len(testmap)):
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi, '%d'%i)
plt.show()

效果图：

orthview半视图

import matplotlib
import healpy as hp
matplotlib.rcParams.update({'font.size': 8})
nside = 4
mmap = np.arange(hp.nside2npix(nside))
hp.orthview(mmap, title='', cbar=None, half_sky=True)
for i in range(len(testmap)):
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi, '%d'%i)
plt.show()

效果图：

选择某个几个像素值：

import matplotlib
matplotlib.rcParams.update({'font.size': 8})

nside_temp = 1
mmap_temp = np.zeros(hp.nside2npix(nside_temp))
mmap_temp[4] = 1
mmap_temp4 = hp.ud_grade(mmap_temp, 4)
nozeros = np.array(np.where(mmap_temp4 != 0))
nside = 4

mmap = np.arange(hp.nside2npix(nside))
hp.mollview(mmap, title='', cbar=None)

sss = 0
for i in range(len(mmap)):
    if i in nozeros:
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi + 0.12, '%d(%d)'%(i, sss), color='r')
        sss = sss + 1
    else:
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi, '%d'%i)
plt.show()

效果图：

rothview视图：

import matplotlib
matplotlib.rcParams.update({'font.size': 8})

nside_temp = 1
mmap_temp = np.zeros(hp.nside2npix(nside_temp))
mmap_temp[4] = 1
mmap_temp4 = hp.ud_grade(mmap_temp, 4)
nozeros = np.array(np.where(mmap_temp4 != 0))
nside = 4

mmap = np.arange(hp.nside2npix(nside))
hp.orthview(mmap, title='', cbar=None, half_sky=True)

sss = 0
for i in range(len(mmap)):
    if i in nozeros:
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi + 0.12, '%d(%d)'%(i, sss), color='r')
        sss = sss + 1
    else:
        theta, phi = hp.pix2ang(nside, i)
        hp.projtext(theta, phi, '%d'%i)
plt.show()

效果图：

在每个像素上画点：

import healpy as hp
import numpy as np
import matplotlib.pyplot as plt
import warnings
warnings.filterwarnings('ignore')

import matplotlib as mpl

# cmap = mpl.cm.summer

m = np.zeros(hp.nside2npix(1))
mrang = np.arange(hp.nside2npix(4))
m[3] = 1
m4 = hp.ud_grade(m, 4)
mmmm = m4*np.NAN
pixs4 = np.array(np.nonzero(m4))

mmmm[pixs4] = mrang[pixs4]
hp.orthview(mmmm, title="", half_sky=True, rot=(-20,30), cbar=None, hold=True, norm='hist')
hp.graticule()
loc = hp.pix2ang(4, np.arange(hp.nside2npix(4)), lonlat=True)
locpoint4 = hp.pix2ang(4, pixs4, lonlat=True)
hp.projscatter(loc, lonlat=True, c='r', s=5)
for i in np.arange(16):
    hp.projtext(locpoint4[0][0][i],  locpoint4[1][0][i], "%d"%i, lonlat=True, fontsize = 16)
plt.savefig('gaussian.png', dpi = 100, transparent = True, bbox_inches='tight')
plt.show()

效果图：