在 healpy mollview 上获取特定值的坐标
Get coordinates for specific value on healpy mollview
是否可以从 healpy mollview plot 中找到 hitmap 特定值的坐标?
例如:
如果我有一个平滑的命中图并想使用遮罩:
mask = (hitmap > 2.054) & (hitmap < 2.056)
在命中图上找到所需的值。
命中图用于创建这样的 mollview 图:
hp.mollview(hitmap, xsize=2000)
然后是否可以通过mask找到mollview plot中满足mask的像素点对应的坐标(lon,lat)?
非常感谢!
最好的方法是通过 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]
是否可以从 healpy mollview plot 中找到 hitmap 特定值的坐标?
例如: 如果我有一个平滑的命中图并想使用遮罩:
mask = (hitmap > 2.054) & (hitmap < 2.056)
在命中图上找到所需的值。
命中图用于创建这样的 mollview 图:
hp.mollview(hitmap, xsize=2000)
然后是否可以通过mask找到mollview plot中满足mask的像素点对应的坐标(lon,lat)?
非常感谢!
最好的方法是通过 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]