带有 healpy 的 2-dim ndarray 的 Gnomonic 投影答案

【问题标题】：Gnomonic projection into 2-dim ndarray with healpy带有 healpy 的 2-dim ndarray 的 Gnomonic 投影
【发布时间】：2013-05-02 15:16:40
【问题描述】：

Gnomview 允许可视化球面地图的 gnomonic 投影。如何将投影存储在二维 numpy.ndarray 中？ healpy 相当于 IDL 的 gnomview 的 map_out 选项？

【问题讨论】：

标签： python numpy astronomy healpy

【解决方案1】：

healpy 不支持这个，所以我实现了，你可以从github仓库查看healpy的开发版本，然后：

In [14]: hp.gnomview(np.arange(12), return_projected_map=True)
Out[14]: 
masked_array(data =
 [[ 4.  4.  4. ...,  4.  4.  4.]
 [ 4.  4.  4. ...,  4.  4.  4.]
 [ 4.  4.  4. ...,  4.  4.  4.]
 ..., 
 [ 4.  4.  4. ...,  4.  4.  4.]
 [ 4.  4.  4. ...,  4.  4.  4.]
 [ 4.  4.  4. ...,  4.  4.  4.]],
         mask = False,
   fill_value = -1.6375e+30)

返回的投影图是一个 2d 掩码的 numpy 数组。这将在下一个healpy 版本中提供。

【讨论】：

对不起，我没听懂。如何使用此选项？ Helpy 告诉我：gnomview() 得到了一个意外的关键字参数 'return_projected_map'
但是有没有什么办法可以通过 healpy 获得该投影？也许与 GnomonicProj 的projmap 一起使用？我试过了，但我无法使用该功能，并且文档帮助不大
我在healpy的开发版中添加了这个，只需下载它：github.com/healpy/healpy/archive/master.zip，解压，cd进入文件夹并运行python setup.py install（必要时使用sudo）
嗨，我在github上写了关于如何构建开发版本的更详细说明：github.com/healpy/healpy/issues/105#issuecomment-17485534
虽然线程已经有点老了 - 是否还有一个函数可以将投影地图写入 FITS 文件，或者我必须自己编写，例如使用 AstroPy？如果包含，我没有找到它...在此先感谢。