问题标签 [healpy]

我最近正在使用healpy。我使用“projplot”制作了以下散点图，并使用了标记大小来定义大小。这是情节： 这里的圆圈大小与值相近

现在我想根据值与颜色条制作相同的图形。我想要银河平面中的以下内容：这里的颜色条显示了值。

我想我需要使用'markercolor'，但我不确定。你能告诉我如何用颜色条生成散点图的代码（如第二张图：每个圆圈都应该根据值获得颜色）吗？先感谢您。

我正在尝试用 healpy 绘制 3 个 numpy 数组。其中两个对应于角度位置theta和phi，另一个是温度。我已经使用了上一个答案中的信息：Plotting a numpy array in healpy来做到这一点，我做了如下：

Where theta, phiandTemp是我的数组（每个都有近 100 万个元素，但看起来与这条线没有不兼容m_sm[indx] = Temp）。然后，我使用以下方式进行投影mollview：

但是我得到了以下地图，其中有一条奇怪的白线图案穿过所有地图。

银河坐标中的温度图，带有奇怪的白线图案

即使NSIDE = 256我获得了相同的模式。我不知道这是否可能是因为我正在使用立方体网格，并且对于网格的每个单元格，我从笛卡尔位置计算相对于盒子中心的角度位置。我能做些什么？或者我的代码有什么问题？提前致谢。

我正在尝试使用healpy HEALPiX 的python 实现。我在 Ubuntu 上，使用 anaconda2。我 pip install healpy==1.9.1 并使用 conda list 检查它是否已安装，并且它使用的是 conda pip，而不是系统。当我“导入healpy”时，我得到：

从我发现的问题似乎与 gcc 编译器有关？我以前在这台计算机上安装过这个版本，所以这似乎是 anaconda 找到正确库的问题。

我已经检查并安装了 gcc 版本 4.8.5。我对连接anaconda和库以取得进展的了解还不够，也找不到太多信息！

使用python healpy，我想显示（并导出）用于天空镶嵌（在不同级别）的 healpix 网格。我该怎么做呢？理想情况下，我需要中点、网格的顶点以及连接它们的线。

任何帮助深表感谢。

我正在使用 hp.mollview() 绘制整个投影未完全填充的地图。这使得地图的边界在哪里变得不清楚，因为白色背景只是融入了图形区域的白色。是否可以在地图区域周围绘制边框？

我有一组大于 10M 的大型对象，带有 RA 和偏角的文件。我想使用 healpix/healpy 制作这些的对数密度全天地图。我当前的代码如下所示：

我得到了错误：

我究竟做错了什么？？

谢谢，尼克

我有一个带有银河坐标的 healpix 地图。我想将该地图转换为赤道坐标系。我知道我可以使用 mollview() 函数在赤道坐标中绘制地图。你有什么办法可以保存这种转换后的地图吗？

谢谢维努

我有一个星系列表要绘制到 healpix 地图上（我正在使用 healpy 来做）每个星系都有一个固定的通量，我需要以这样的方式绘制它们，以便每个星系的通量在地图上保持不变.

这是我的代码：

我得到的结果是 140 而不是 30，这是通量的真实总和。

我知道发生了什么，并且相同的通量分布在多个像素上（第一个星系为 6 个像素，第二个星系为 4 个像素），我知道我可以这样做：

这将保留总光子数，但不一定会保留我需要的每个星系的光子数。即这种方法最终导致第一个星系贡献了 12.86 的通量，而第二个星系贡献了 17.14 的通量。

有没有办法在每个坐标占用多少像素之前计算出，然后根据这个改变倾倒的通量？

提前致谢！

我正在使用 healpy.sphtfunc.smoothing 对我的二进制蒙版进行变迹，我遇到的问题是，如果我的温度降低了 100K 并且我制作了一个与切割相对应的二进制蒙版，那么在使用上述例程对蒙版进行变迹后，当我应用它时在我的地图上，我得到 120K 或超过 100K 的数字。所以我很困惑是否对二进制掩码或 (map*Binary_mask) 进行切趾

这是没有二值掩码变迹的结果图

这是二进制掩码 5 度变迹后的结果图

我使用的代码需要使用 map2alm 和 alm2map（主要是后者）在 healpy 中进行频繁的 spin-2 转换。但是，当我越来越高时，这种变换不可避免地会变慢。

我正在使用带有openmp的healpy，并且由于必须进行大量的转换，这仍然会导致完成一项工作需要很长时间。有没有进一步的方法来尝试加速可能适用于这种情况的转换？例如，由于我只对部分天空感兴趣，有没有办法在使用 alm2map 时只重建指定像素中的值而不是整个天空，这样会更快吗？

谢谢你的帮助