问题标签 [healpy]

我正在尝试在我创建的 venv 中安装 healpy。我成功安装了 numpy 和 astropy 但是当我运行以下命令时：

或者

安装失败并显示一条很长的消息，该消息以一些开头

并以红色文本结尾：

有人可以帮我在我的venv中正确安装healpy吗？

我想知道在 healpy 地图上绘制矢量或定向线段的最简单方法，例如本文的图 22：https ://arxiv.org/abs/1807.06212 ，截图如下。matplotlib 的箭袋之类的东西是否可用于在 healpy 中绘制例如正交或 mollweide 投影？

我healpy.query_polygon用来获取多边形内的 healpix 索引列表。根据文档：

vertices：顶点数组，包含多边形的顶点，形状为 (N, 3)。

但是当我尝试从以下多边形中获取所有索引时，会RuntimeError: Unknown exception出现一个：

在[1]：

输出[1]：

在[2]：

输出[2]：

在这里，您可以看到这些点位于球体上的可视化效果。

只是为了好玩，我尝试转置输入数组，所以它的形状变成了 (3, 4)。Unknown exception问题消失了。但是这样的输入与文档相矛盾，所以我不相信。

在[1]：

输出[1]：

在[2]：

输出[2]：

我将不胜感激任何建议。

我有一个高分辨率的 healpix 贴图（nside = 4096），我想在给定半径的磁盘中平滑，比如说 10 arcmin。

作为一个非常新的healpy并阅读了文档，我发现一个 - 不太好的 - 方法是执行“锥形搜索”，即在每个像素周围找到磁盘内的像素，平均它们并给出这个中心像素的新值。然而，这是非常耗时的。

我认为 healpy.sphtfunc.smoothing 函数可能会有所帮助，因为它声明您可以输入任何自定义光束窗口函数，但我根本不明白它是如何工作的......

非常感谢你的帮助 ！

假设我有一个包含 WCS 标头的 FITS 文件，所以我可以这样做：

listofhdus[0].data 将包含一个 2D numpy 数组（纽约乘以 NX），它对应于银河坐标系中完整天空的一小部分。

如果我想将该天空图叠加到以下带有 Mollweide 投影的散点图中，那么将该 2d numpy 转换为 healpix 格式（称为 WCS）的最佳方法是什么？

其中经度和纬度是两个二维数组（使用 numpy.meshgrid 生成），其中包含我感兴趣的一些天体物理源的银河坐标？

我想我可以尝试将我的 healpix 像素转换为坐标，以某种方式将它们与我的天空图中可用的像素匹配并从那里插入值，但必须有一些更简单、更优雅和精确的东西，对吧？

我有一张包含大约 250 万个元素的 healpy 地图。我想将地图顺时针旋转约 45 度。我找不到一个简单的方法来执行此操作。我想我可以将所有像素转换为 theta/phi 然后旋转，但是每个像素都有一个相应的强度值，需要在旋转过程中保持完整。

有没有更直接的方法来做到这一点？

在 IDL 版本的 HEALPix (pix2vec_nest) 中有一个很好的函数，它可以返回一个带有像素顶点坐标的向量 (n,3,4)。C 版本中的相同功能不会做同样的事情，看起来你可以只知道像素中心的坐标。有没有办法计算顶点的坐标？

这个东西的几何和数学对我来说有点难写一个 HEALPix 像素化的生成器。我试图理解 IDL 代码以便将功能添加到 C 库中，但它看起来非常混乱，带有难以辨认的变量名称。

这段代码在某种程度上看起来很简单，但文件“Healpix_3.50\src\idl\toolkitpix2vec_nest.pro”中还有其他代码可以使用所有这些参数调用此函数..如果您有想法..提前致谢。

我有一张健康地图，其中每个像素值都是一些有意义的物理测量值。我想将地图分成边长为 5 度的“正方形”。基本上我想平均每个 5 deg^2 bin 中的像素值。有没有直接的方法可以做到这一点？

我有一系列 theta 和 phi 使用 healpy pix2ang 命令，然后转换为 RA, Decl.::

我只是想将这些投影到例如 Aitoff 类型的投影上，但是对于我的生活来说，我无法通过以下方式弄清楚如何做到这一点： https ://healpy.readthedocs.io/en/latest/generated/healpy.visufunc .projplot.html

并没有真正做任何事情。

基本上我想要的只是 hp.cartview 的功能，但我不希望我的机器在每次调用 cartview 函数时都浪费内存来绘制实际地图。如何以二维数组的形式获得 healpy 的笛卡尔投影，而不必每次都绘制投影？