ImageProfileEstimator

class gammapy.maps.ImageProfileEstimator(x_edges=None, method='sum', axis='lon', center=None)

Bases: object

Estimate profile from image.

Parameters

x_edgesAngle

Coordinate edges to define a custom measument grid (optional).

method[‘sum’, ‘mean’]

Compute sum or mean within profile bins.

axis[‘lon’, ‘lat’, ‘radial’]

Along which axis to estimate the profile.

centerSkyCoord

Center coordinate for the radial profile option.

Examples

This example shows how to compute a counts profile for the Fermi galactic center region:

import matplotlib.pyplot as plt
from gammapy.maps import ImageProfileEstimator
from gammapy.maps import Map
from astropy import units as u

# load example data
filename = '$GAMMAPY_DATA/fermi-3fhl-gc/fermi-3fhl-gc-counts.fits.gz'
fermi_cts = Map.read(filename)

# set up profile estimator and run
p = ImageProfileEstimator(axis='lon', method='sum')
profile = p.run(fermi_cts)

# smooth profile and plot
smoothed = profile.smooth(kernel='gauss')
smoothed.peek()
plt.show()

Methods Summary

run(self, image[, image_err, mask])

Run image profile estimator.

Methods Documentation

run(self, image, image_err=None, mask=None)

Run image profile estimator.

Parameters

imageMap

Input image to run profile estimator on.

image_errMap

Input error image to run profile estimator on.

maskMap

Optional mask to exclude regions from the measurement.

Returns

profileImageProfile

Result image profile object.