问题标签 [healpy]

我正在应用卷积技术来卷积 2 个数据集，一个 nside = 256 的 healpix 图和一个形状的主光束 (256, 256)，以便测量卷积后的 healpix 图的总强度。我的问题是，在将我的地图与主光束卷积后，我在卷积后的地图中得到了环。我尝试使用 lanczos 或 Gaussian 内核对其进行归一化以处理环，但所有这些方法都失败了。

在下面的代码中，我使用 scipy 中的查询函数在给定半径内的 healpix 地图中搜索最近的像素，并使用地图坐标获取主光束中相应像素的乘积之和。我得到的最终图像中有环。请问谁能帮我解决这个问题？提前致谢。

我有一个healpy图，它是在CMB的特定透镜模拟期间获得的偏转（透镜势的梯度）。如果可能的话，我想使用 healpy 获得透镜潜力的地图。我注意到有一个 healpy 函数alm2map_der1()会给我一个 healpy 地图及其给定地图的施舍的一阶导数。我假设这个一阶导数是地图的梯度 - 如果我错了，请纠正我。本质上我想知道我是否可以使用healpy来做这个的倒退过程。我想去除渐变，我只想要透镜潜力。

到目前为止，我一直在尝试使用偏转和透镜势能谱之间的关系。偏转的 Cls = l(l+1) * 透镜势的 Cls，并将其重新排列为：透镜势的 Cls = 偏转的 Cls / l(l+1)，然后使用 synfast 将其转换回映射。我似乎没有得到正确的地图。

有没有更好的方法来做我想做的事情？也许甚至不使用healpy？

I have a healpy full-sky map, and need to chose some sky patches centered around given pixels. hp.query_disc finds the pixels within a given disc, which is nice. But I would like to be able to add/subtract on different sky patches, and so the fact that different discs cannot be the same size causes a problem.

I have also used

hp.gnomview(return_projected_map=True)

to look at my sky patches, and found that each returned 2D array has the same size. However since I would probably do stuff on many submaps, I do not want the maps to be drawn every time. Is it possible to obtain only the 2D array? And if so, is it fine to use plt.imshow after adding/subtracting those arrays?

I feel like there should be more elegant ways to achieve what I want -- Are there any existing routines that I overlooked? Thanks!

我有一张来自 AKARI Far Infrared Surveyor 数据库（公开发布）的 HEALPIx 全天地图。我曾尝试使用“平滑”地图healpy，但结果看起来很奇怪。有没有更好的办法？然而，我的问题与任何全天 HEALPIx 地图（即 IRAS、Planck、WISE、WMAP）有关。

我的目标是将此 AKARI 地图的有效点扩散函数“平滑”到 1 度的角分辨率（原始数据的 PSF 约为 1 弧分）。这样我就可以将远红外 AKARI 地图与分辨率较低的微波地图（特别是异常微波前景的那些）进行比较。

在下面的示例中，我使用的是降级版本的地图，因此它足够小，可以上传到 Github。这意味着像素约为 3.42 角分。在 PSF 平滑之前，我不会过多地降低像素比例，但这只是一个例子：

我已经尝试过 healpy.sphtfunct.smoothing 例程（https://healpy.readthedocs.org/en/latest/generated/healpy.sphtfunc.smoothing.html#healpy.sphtfunc.smoothing）。据我所知，smoothing转换映射成球谐函数，然后与高斯卷积，然后将其转换回空间映射。

我已经在 github 存储库中保存了 ipython 笔记本以及低分辨率 .FITS HEALpix 地图，这里：

https://github.com/aaroncnb/healpy_smoothing_test

（您需要healpy安装软件包）

通过在笔记本中运行代码，您可以轻松地可视化我遇到的问题 - 在对地图进行平滑之后，有一些奇怪的“伪影”，就好像像素已经被迭代箱平均，而不是使用圆形高斯平滑轮廓。我期望看到的只是输入图的模糊版本。

在平滑完成之前，我认为我缺少关于转换为球谐函数的一些基本知识。

有没有人尝试过在 HEALPIx 地图上进行这种全天空平滑？

我相信另一种选择是将地图转换为标准的矩形数组，然后进行平滑处理。但是，我仍然对在不离开 HEALPIx 格式的情况下解决问题感到好奇。

有许多 Healpix IDL 例程旨在查找属于由其顶点定义的某个几何区域（例如球形三角形、球形多边形）的像素索引。有query_*例程（例如query_triangle）。请参阅此处的文档：

http://healpix.jpl.nasa.gov/html/idlnode45.htm

我希望在我的 healpy 程序中使用这些像素索引。任何一个

query_*(A) 我可以将 IDL例程中的像素索引输出列表保存为一种data_file.save格式。然后，您可以.sav使用各种模块将此像素索引文件导入 Python，例如http://www.astropython.org/packages/idlsave94/

(B) 以某种方式根本不使用 IDL 会方便得多！Healpy 有几个与像素相关的功能，但似乎没有办法query_*单独使用 healpy 来“转换”IDL 例程。

有没有办法query_polygon使用healpy？是否有可能做到这一点？

I'm trying to use others colormaps on healpy.mollview I succeded with this code

but I get an unexpected blue background and there is no way I can set it to white

我一直在尝试使用 CAMB 和 healpy 制作 E 和 B 模式地图，但是到目前为止我制作的东西看起来并不多。

我正在使用 CAMB 在以下链接生成模拟地图：http: //lambda.gsfc.nasa.gov/toolbox/tb_camb_form.cfm，除了选中张量 cl 的框和不同的 r 之外，我保留了所有默认参数.

有了这个，我一直在使用以下代码来生成地图：

快速的谷歌搜索显示与我得到的结果截然不同的结果，您可以在以下链接中看到： http://tinypic.com/view.php ?pic=ifvrpt &s=9# .VmY7Jx8So8o http://tinypic.com /view.php?pic=a5ngx&s=9#.VmY7KR8So8o

我期待在 E 和 B 中看到某些图案，因为它们是渐变和卷曲成分，这在我在网上找到的图像中很明显，但在我的产品中却没有。

所以我有三个问题：

1. 我用来生成地图的逻辑是否有错误？

2. 是否需要额外的过滤才能看到预期的 E 和 B 模式？

   • 类似于以下链接： http ://cosmologist.info/polar/EBsupport.html

3. 在 healpy 以外的 healpix 版本中，mollview 中有一个名为 POLARIZATION 的参数（http://healpix.jpl.nasa.gov/html/idlnode33.htm#idl:mollview:polarization），其中设置 POLARIZATION = 3 允许您过度绘制极化信息作为无头矢量。这就是我认为他们在这里所做的事情：http: //bicepkeck.org/media/b_over_b_rect_BICEP2.png 和这里的第一张图片： https ://inspirehep.net/record/787860/plots ，有没有办法做到这一点在愈合？由于 POLARIZATION 参数似乎不包含在此版本中。

使用 healpy 时，我可以使用在 Mollview 中绘制 Healpix 地图

或在教程中

哪个输出

有没有办法只显示地图的某些区域？例如，只有上半球，还是只有左侧？

我想到的是只查看一小块天空以查看小点源，或者类似于 LSST 的“半天空”投影

I have tried debugging my code and I've realised that ultimately it breaks down when I try to save my AltAz coordinates into a .csv file because its not a numpy array, its a SkyCoord object. Could someone suggest a simple of way of converting a large table of Equatorial coordinates to AltAz or how I can get my code to save to file.

我有一张银河坐标地图，我需要将它保存在另一个文件的赤道坐标中。我知道我可以使用：

这应该返回一个numpy存储在map_rot. 但是当我阅读时map_rot，我发现它是一个仅填充了-inf值的 masked_array 和mask=False, fill_value=-1.6735e+30 （因此，显然-inf不是掩码）。此外，元素的总数与我期望的地图（ ）map_rot的像素数不匹配。npix=12*nside**2例如，如果nside=256我希望获得npix=786432, whilemap_rot有400*800=320000元素。这是怎么回事？

（我已经看过这篇文章，但是我有一个极化地图，所以我需要旋转斯托克斯的参数。既然mollview知道怎么做，我试图直接从获取新地图mollview。）