Answered: How to create a HEALpix map with Python?

Fri, 09 Jun 2023 10:23:16 +0000

healpy has a function for that. Here is an example:

import healpy as hp
from healpy.newvisufunc import projview, newprojplot

import matplotlib.pyplot as plt
import matplotlib as mpl
import numpy as np

def coords_to_hpx(lon, lat, nside, radec=False):
    npix = hp.nside2npix(nside)
    if radec:
        eq = SkyCoord(lon, lat, 'icrs', unit='deg')
        l, b = eq.galactic.l.value, eq.galactic.b.value
    else:
        l, b = lon, lat

    theta = np.radians(90. - b)
    phi = np.radians(l)
    indices = hp.ang2pix(nside, theta, phi)
    idx, counts = np.unique(indices, return_counts=True)
    hpx_map = np.zeros(npix, dtype=int)
    hpx_map[idx] = counts
    return hpx_map

data = np.loadtxt('photon_lists/p2209fb.txt')
l = data[:,3]
b = data[:,4]
l[l >= 180] -= 360

hpx = coords_to_hpx(l, b, nside=8, radec=False)
hpx[hpx > 50] = 50
projview(
        hpx,
        coord=["G"],
        graticule=True,
        graticule_labels=True,
        unit="counts",
        xlabel="longitude",
        ylabel="latitude",
        cb_orientation="horizontal",
        cbar_ticks=np.round(np.linspace(min(hpx), max(hpx), 10, endpoint=True)),
        projection_type="aitoff",
        )

plt.show()

CMS/Wiki for a colloboration

Fri, 09 Jun 2023 10:17:33 +0000

Is there a good resource on how to choose and set up a content management system or wiki for a collaboration? We need to migrate from Confluence.

How to download Meerkat data?

Mon, 22 Nov 2021 17:31:13 +0000

hallo

How to analyze HESS data with gammapy?

Fri, 19 Nov 2021 12:58:38 +0000

Given are the event lists of a selection of HESS runs. How to analyze the data with gammapy?

PUNCH4NFDI Marketplace - Recent questions and answers in Tutorials

Answered: How to create a HEALpix map with Python?

CMS/Wiki for a colloboration

How to download Meerkat data?

How to analyze HESS data with gammapy?