问题标签 [healpy]

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

示例图

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

非常感谢您的帮助！

我想使用 healpy 在一个图中绘制许多子图。如何：

1. 设置颜色条的位置？
2. 设置颜色条的刻度和刻度标签？
3. 设置子图的位置和大小？

我想根据一般坐标生成如图 1 所示的图，该图是在 MATLAB 中绘制的

现在，我只使用healpy将其绘制如下：

生成图 3（类似于图 2）的类似代码如下：

代码生成图 3

我正在尝试使用healpy中的rotate_skymap函数在python 3中旋转天空图，如下所示：

当我运行代码时，虽然我收到了这个错误：

这是什么意思，我该如何解决？

是否可以从healpy mollview plot中找到命中图特定值的坐标？

例如：如果我有一个平滑的命中图并想使用遮罩：

在命中图上找到所需的值。

命中图用于创建如下所示的 mollview 图：

那么是否可以通过使用掩码找到与mollview 图中满足掩码的像素对应的坐标 (lon,lat) ？

提前谢谢了！

我正在尝试使用healpy以下脚本以 HEALPix 格式编写 FITS 文件：

这会引发以下错误：

任何帮助表示赞赏！

我浏览了 astropy 文档并得出结论，如果给定时间，没有可以转换为赤道坐标 RA、dec 的具有 lon、lat 坐标的本地 ECEF（地心、地球固定）框架。这是真的？

最终，我想使用以下方法创建地图：

map = HEALPix(nside=NSIDE, order='nested', frame=MY_REF_FRAME())

MY_REF_FRAME = ITRS显然不是一种选择。

如果有人可以帮助我找到是否可以使用最新版本的 astropy 来定义这种情况，我将不胜感激。

谢谢！

埃里克

当我沿经度或纬度旋转 Healpix 地图时，我得到了错误的行为。我可能在这里遗漏了一些明显的东西，但到目前为止，我找不到什么。

见演示：

和结果：

• 原图
• 旋转经度
• 旋转纬度

地图中沿经度旋转的点是垂直平移的，而沿纬度旋转的则是水平平移的。我期待相反的情况。

关于这里有什么问题的任何提示？

E.

我正在从输入角功率谱 Cl 生成随机 healpix 图。如果我使用healpy.synalm，然后使用healpy.alm2map，最后通过在生成的map上运行healpy.anafast来测试map，输出和输入功率谱不一致，尤其是在高l时，输出功率谱高于输入（见下面代码产生的情节）。如果我直接使用healpy.synfast，我会得到一个与输入一致的输出功率谱。如果我使用healpy.synfast 中的alms 并使用healpy.alm2map 从synfast alms 生成地图，这同样适用。当我查看synfast的源代码时，似乎只是调用了synalm和alm2map，所以我不明白，为什么他们的结果不一致。我的测试代码如下所示：

NSIDE 越低，偏移越大。

非常感谢您的帮助。

我在笛卡尔投影中有一张带有球坐标系的地图。头文件读取 CTYPE1 = 'GLON-CAR' 和 CTYPE2 = 'GLAT-CAR'。我想将此地图转换为 healpix 投影。我知道我需要使用工具“ang2pix”将球坐标中的地图转换为 healpix 像素。但我无法理解如何用笛卡尔投影中可用的值填充像素。

我正在尝试使用healpy 中的 anafast从蒙版地图中计算角功率谱。与F90 版本相比，anafast 的 python 版本如何处理掩码的效果，后者明确将掩码文件作为可选输入？

谢谢！