7

我有一个用 HEALPY 制作的 HealPix 图,如Healpy:从数据到 Healpix 地图(像素较少,例如采用 nside=2,请参见下面的代码)。

import healpy as hp
import numpy as np
import matplotlib.pyplot as plt

# Set the number of sources and the coordinates for the input
nsources = int(1.e4)
nside = 2
npix = hp.nside2npix(nside)

# Coordinates and the density field f
thetas = np.random.random(nsources) * np.pi
phis = np.random.random(nsources) * np.pi * 2.
fs = np.random.randn(nsources)

# Go from HEALPix coordinates to indices
indices = hp.ang2pix(nside, thetas, phis)

# Initate the map and fill it with the values
hpxmap = np.zeros(npix, dtype=np.float)
hpxmap[indices] += fs[indices]

# Inspect the map
hp.mollview(hpxmap)

示例图

如何在绘图上的每个 HEALPIx 的中心写一个带有值的文本?例如,如何使用类似的数组为每个“像素”编写标识符 range(len(hpxmap))

非常感谢您的帮助!

4

0 回答 0