问题标签 [healpy]

根据我对量子物理学和大学数学的记忆，对于每种模式l，都存在m = l-1, l, l+1。为什么 healpix （在我的情况下，特别是healpy）球谐函数例程，例如healpy.sphtfunc.alm2map，返回相同长度的l和数组？m

I'm just starting to work with Healpy and have noticed that if I use a map to get alm's and then use those alm's to generate a new map, I do not get the map I started with. Here's what I'm looking at:

As you can see, new_map doesn't match the input map, m. I imagine there's some subtlety to these functions that I'm missing. Any idea?

我对healpy和python很陌生，但我想在银河坐标系中用银河平面和赤道平面绘制天空图。我遇到的第一个问题是我无法划清界限：

我试图遵循这个文档： https ://healpy.readthedocs.org/en/1.5.0/generated/healpy.visufunc.projplot.html

另外我不确定在哪里可以获得银河平面和赤道平面的坐标。

我正在尝试构建一个 Python 库 Healpy。我的第一次尝试是使用 pip。输入后：

我得到了结果：

之后我尝试自己构建库 - 我下载了源代码并使用命令构建了 cfitsio：

之后，我得到了一个libcfitsio.a文件，pkg-config包含 .pc 文件的文件夹和包含文件夹中的一些 .h 文件。但输入后：

我得到了同样的错误。在 Windows 中构建 lib 后，我得到了同样的错误。有没有人有类似的问题？对于任何建议，我都会非常感激。

这个问题有两个部分：

第 1 部分：我遇到的主要问题是尝试使用 healpy mollview 绘制半边天图。我正在使用的当前解决方法是：

这感觉有点笨拙，所以我正在寻找一种更好的方法来解决这个问题。

第 2 部分：第二个问题是标题。有没有办法制作完全没有标题的地图？

提前致谢

我有宇宙微波背景的健康地图。每张地图都有三个字段。我在代码中编写了以下函数来计算一般地图数量的第一个字段的功率谱：

'files' 输入是映射文件名的数组作为字符串。上面的代码为 ipython 中的 10 个映射返回以下输出：

这在一个列表中有 10 个数组。每个阵列具有特定地图的第一场功率谱数据。我想要做的是在我的代码中调用此列表中的各个数组，以便对每个单独的数据数组执行一个函数。

我看了个数组里面的pages Array？和访问python中数组列表中的元素，但没有发现任何一个像我希望的那样有用。

我正在尝试创建单极子 (l=0, m=0)、偶极子 (l=1, m=-1,0,1) 和四极子 (l=2, m=-2,-1, 0,1,2) 的球谐函数使用 Healpy。到目前为止，我在 healpy 教程的帮助下构建了这个非常简单的代码：

但这只是返回单极子 (x) 的值以及单极子和偶极子 (y) 的值：

当我尝试 hp.mollview(x) 时，出现以下错误：

如何创建单极子和偶极子的地图。另外，我如何创建四极杆的地图？

我正在使用以下脚本来计算银河坐标（以度为单位）到天体坐标的银河中心 (GC) 位置：

这些是我得到的结果：

首先，天体坐标和银河坐标中的数字都是错误的。有可能我使用了错误的函数（我希望如此），或者函数本身有问题。有人知道出了什么问题吗？

第二：看来，我得到了以弧度表示的数字，对吗？

我有一系列具有两个不同索引 i 和 j 的地图。让它像 map_series[i][j] 一样被索引。

编辑 1/21：一个最小的工作示例是这样的 map_series=np.array([np.array([np.arange(12) + 0.1*(i+1) + 0.01*(j+1) for j in range(3)]) for i in range(5)])

我想对每个人应用相同的面具；如果map_series是一维的，这些都起作用。

我可以想象应用这些地图的几种不同方式：

(A) 将掩码应用于整个阵列：

(B1) 将掩码应用于数组的每个元素：

(B2) 将掩码应用于数组的每个元素：

(C) 以 Python 方式枚举列表：

然而，所有这些都无法提供我想要的输出。

在尝试 (A) 或 (C) 时，我在第一步后得到一个错误，告诉我 TypeError: bad number of pixel。

在尝试 (B1) 时，我没有收到错误，但我也没有 maps_series_ma 的元素有掩码；事实上，它们甚至看起来都不是 hp.ma 对象。不过，奇怪的是：当我返回 temp 时，它确实有适当的掩码。

在尝试（B2）时，我收到错误 AttributeError: 'numpy.ndarray' object has no attribute 'mask' （在查看我的语法后，我完全明白了！）

我有点困惑如何去做。(A) 和 (B1) 对我来说似乎都可以接受......

非常感谢任何帮助，谢谢，山姆

我想知道 healpy 中的 alm2fits 是什么。

假设我在 healpy 中有 3 个用于 T、E 和 B 的施舍，我如何导出它们以便通过 fit2alm 的 fortran 代码读取它们？

hp.fitsfunc.write_alm('alms.fits',[alm_T,alm_E,alm_B])

似乎不喜欢列表和

hp.fitsfunc.mwrfits(filename= 'alm.fits',data=[alm_T,alm_E,alm_B],colnames=['T','E','B'],keys=None)

似乎以错误的顺序保存它们，因此 fit2alm 引发 'Inconsistent l^2+l+m+1 -> l,m mapping' 错误

谢谢

一个。