如何绘制HEALPix分箱？基于Healpy实现天球等面积分箱可视化

Healpy提供了直接绘制HEALPix像素边界的内置函数，无需手动计算坐标，下面是实现类似你现有可视化效果的简单方法：

核心实现步骤 #

• 使用healpy.draw_edges()函数直接绘制指定分辨率下的HEALPix像素边界
• 结合matplotlib的Aitoff投影（和你当前代码一致）完成天球可视化
• 调整样式参数匹配你的现有图效果

完整代码示例 #

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

# 设置HEALPix分辨率参数nside（必须是2的幂，值越大分辨率越高）
nside = 8
# 计算对应总像素数
npix = hp.nside2npix(nside)

# 创建Aitoff投影子图
plt.subplot(111, projection="aitoff")
plt.title("HEALPix Equal Binning", fontsize=30)

# 绘制HEALPix像素边界，设置颜色和线宽
hp.draw_edges(nside, color='k', linewidth=1)

# 调整坐标轴样式，和你的可视化对齐
ax = plt.gca()
ax.set_xticklabels(ax.get_xticklabels()[::-1])
plt.xticks(fontsize=20)
plt.yticks(fontsize=20)
plt.xlabel("L", fontsize=25)
plt.ylabel("B", fontsize=25)

plt.show()

代码说明 #

• nside：HEALPix的分辨率参数，必须为2的整数次幂（如1、2、4、8等），nside越大，像素数量越多，分辨率越高
• hp.draw_edges()：Healpy内置函数，自动计算并绘制所有HEALPix像素的边界，无需手动处理天球坐标转换
• 坐标轴调整部分和你原代码保持一致，确保可视化风格匹配

这段代码会生成和你现有等面积分箱图风格一致的HEALPix分箱可视化，HEALPix的像素在赤道附近近似方形，两极附近近似菱形，保证每个像素面积相等。

内容的提问来源于stack exchange，提问作者Kyle