Note
A new set of map and cube classes is being developed in gammapy.maps
and long-term will replace the existing gammapy.image.SkyImage
and
gammapy.cube.SkyCube
classes. Please consider trying out gammapy.maps
and changing your scripts to use those new classes. See Data Structures for Images and Cubes (gammapy.maps).
Cube Style Analysis (gammapy.cube
)¶
Introduction¶
The cube
module bundles functionality for combined spatial and
spectral analysis (cube style analysis) of gamma-ray sources.
Some information on cube style analysis in gamma-ray astronomy can be found here:
Getting Started¶
Use SkyCube
to read a Fermi-LAT diffuse model cube:
>>> from gammapy.cube import SkyCube
>>> filename = '$GAMMAPY_EXTRA/test_datasets/unbundled/fermi/gll_iem_v02_cutout.fits'
>>> cube = SkyCube.read(filename, format='fermi-background')
>>> print(cube)
Sky cube flux with shape=(30, 21, 61) and unit=1 / (cm2 MeV s sr):
n_lon: 61 type_lon: GLON-CAR unit_lon: deg
n_lat: 21 type_lat: GLAT-CAR unit_lat: deg
n_energy: 30 unit_energy: MeV
Use the cube methods to do computations:
import astropy.units as u
emin, emax = [1, 10] * u.GeV
image = cube.sky_image_integral(emin=emin, emax=emax)
image.show('ds9')
Using gammapy.cube
¶
Gammapy tutorial notebooks that show examples using gammapy.cube
:
Reference/API¶
gammapy.cube Package¶
Sky cubes (3-dimensional: energy, lon, lat).
Functions¶
compute_npred_cube (flux_cube, exposure_cube, …) |
Compute predicted counts cube. |
compute_npred_cube_simple (flux_cube, …) |
Compute predicted counts cube (using a simple method). |
make_background_cube (pointing, obstime, bkg, …) |
Calculate background predicted counts cube. |
make_exposure_cube (pointing, livetime, aeff, …) |
Calculate exposure cube. |
Classes¶
CombinedModel3D (spatial_model, spectral_model) |
Combine spatial and spectral model into a 3D model. |
SingleObsCubeMaker (obs, empty_cube_images, …) |
Compute SkyCube images for one observation. |
SkyCube ([name, data, wcs, energy_axis, meta]) |
Sky cube with dimensions lon, lat and energy. |
SkyCubeImages ([name, images, wcs, energy, meta]) |
Class to represent connection between SkyImage and SkyCube . |
StackedObsCubeMaker (empty_cube_images[, …]) |
Compute stacked cubes for many observations. |