在健康的Mollview上获取特定值的坐标

时间:2018-08-12 14:28:18

标签: python coordinates masking healpy

是否可以从健康的软体动物视图图中找到匹配的特定值的坐标?

例如: 如果我有一个平滑的匹配图,并且想使用蒙版:

mask = (hitmap > 2.054) & (hitmap < 2.056)

查找匹配图上的所需值。

命中图用于创建像这样的Mollview图:

hp.mollview(hitmap, xsize=2000)

然后可以通过使用遮罩在 mollview图中找到与满足遮罩的像素相对应的坐标(lon,lat)吗?

非常感谢!

1 个答案:

答案 0 :(得分:0)

最好的方法是通过hp.pix2ang()

import healpy as hp
import numpy as hp

# Get the nside and number of pixels in your map
nside = hp.get_nside(mask)
npix = hp.nside2npix(nside)

# Use pix2ang to get the (l, b) coordinates for each pixel
glons, glats = hp.pix2ang(nside, np.arange(npix), lonlat=True)

这将为您提供地图中每个像素的所有坐标(l,b)的列表。要只得到面具中的那些,只需使用:

glons_eff = glons[mask]
glats_eff = glats[mask]