TemplateSpatialModel#

class gammapy.modeling.models.TemplateSpatialModel(map, meta=None, normalize=True, interp_kwargs=None, filename=None, copy_data=True)[source]#

Bases: gammapy.modeling.models.spatial.SpatialModel

Spatial sky map template model.

For more information see Template spatial model.

Parameters
mapMap

Map template.

metadict, optional

Meta information, meta[‘filename’] will be used for serialization

normalizebool

Normalize the input map so that it integrates to unity.

interp_kwargsdict

Interpolation keyword arguments passed to gammapy.maps.Map.interp_by_coord. Default arguments are {‘method’: ‘linear’, ‘fill_value’: 0}.

Filenamestr

Name of the map file

copy_databool

Create a deepcopy of the map data or directly use the original. True by default, can be turned to False to save memory in case of large maps.

Attributes Summary

covariance

default_parameters

evaluation_bin_size_min

evaluation_radius

Evaluation radius (Angle).

evaluation_region

Evaluation region

frame

frozen

Frozen status of a model, True if all parameters are frozen

is_energy_dependent

map

Template map (Map)

parameters

Parameters (Parameters)

phi_0

position

SkyCoord

position_error

Get 95% containment position error as (EllipseSkyRegion)

position_lonlat

Spatial model center position (lon, lat) in rad and frame of the model

tag

type

Methods Summary

__call__(lon, lat[, energy])

Call evaluate method

copy(**kwargs)

evaluate(lon, lat[, energy])

Evaluate the model at given coordinates.

evaluate_geom(geom)

Evaluate model on Geom

freeze()

Freeze all parameters

from_dict(data)

from_parameters(parameters, **kwargs)

Create model from parameter list

from_position(position, **kwargs)

Define the position of the model using a sky coord

integrate_geom(geom[, oversampling_factor])

Integrate model on Geom or RegionGeom.

plot([ax, geom])

Plot spatial model.

plot_error([ax])

Plot position error

plot_grid([geom])

Plot spatial model energy slices in a grid.

plot_interative([ax, geom])

Plot spatial model.

read(filename[, normalize])

Read spatial template model from FITS image.

reassign(datasets_names, new_datasets_names)

Reassign a model from one dataset to another

to_dict([full_output])

Create dict for YAML serilisation

to_region(**kwargs)

Model outline from template map boundary (RectangleSkyRegion).

unfreeze()

Restore parameters frozen status to default

write([overwrite])

Attributes Documentation

covariance#
default_parameters = <gammapy.modeling.parameter.Parameters object>#
evaluation_bin_size_min#
evaluation_radius#

Evaluation radius (Angle).

Set to half of the maximal dimension of the map.

evaluation_region#

Evaluation region

frame#
frozen#

Frozen status of a model, True if all parameters are frozen

is_energy_dependent#
map#

Template map (Map)

parameters#

Parameters (Parameters)

phi_0#
position#

SkyCoord

position_error#

Get 95% containment position error as (EllipseSkyRegion)

position_lonlat#

Spatial model center position (lon, lat) in rad and frame of the model

tag = ['TemplateSpatialModel', 'template']#
type#

Methods Documentation

__call__(lon, lat, energy=None)#

Call evaluate method

copy(**kwargs)#
evaluate(lon, lat, energy=None)[source]#

Evaluate the model at given coordinates. Note that, if the map data assume negative values, these are clipped to zero.

evaluate_geom(geom)#

Evaluate model on Geom

Parameters
geomWcsGeom
Returns
Map
freeze()#

Freeze all parameters

classmethod from_dict(data)[source]#
classmethod from_parameters(parameters, **kwargs)#

Create model from parameter list

Parameters
parametersParameters

Parameters for init

Returns
modelModel

Model instance

classmethod from_position(position, **kwargs)#

Define the position of the model using a sky coord

Parameters
positionSkyCoord

Position

Returns
modelSpatialModel

Spatial model

integrate_geom(geom, oversampling_factor=None)#

Integrate model on Geom or RegionGeom.

Integration is performed by simple rectangle approximation, the pixel center model value is multiplied by the pixel solid angle. An oversampling factor can be used for precision. By default, this parameter is set to None and an oversampling factor is automatically estimated based on the model estimation maximal bin width.

For a RegionGeom, the model is integrated on a tangent WCS projection in the region.

Parameters
geomWcsGeom or RegionGeom

The geom on which the integration is performed

oversampling_factorint or None

The oversampling factor to use for integration. Default is None: the factor is estimated from the model minimimal bin size

Returns
Map or gammapy.maps.RegionNDMap, containing

the integral value in each spatial bin.

plot(ax=None, geom=None, **kwargs)[source]#

Plot spatial model.

Parameters
axAxes, optional

Axis

geomWcsGeom, optional

Geom to use for plotting.

**kwargsdict

Keyword arguments passed to plot()

Returns
axAxes, optional

Axis

plot_error(ax=None, **kwargs)#

Plot position error

Parameters
axAxes, optional

Axis

**kwargsdict

Keyword arguments passed to plot()

Returns
axAxes, optional

Axis

plot_grid(geom=None, **kwargs)#

Plot spatial model energy slices in a grid.

Parameters
geomWcsGeom, optional

Geom to use for plotting.

**kwargsdict

Keyword arguments passed to plot()

Returns
axAxes, optional

Axis

plot_interative(ax=None, geom=None, **kwargs)[source]#

Plot spatial model.

Parameters
axAxes, optional

Axis

geomWcsGeom, optional

Geom to use for plotting.

**kwargsdict

Keyword arguments passed to plot()

Returns
axAxes, optional

Axis

classmethod read(filename, normalize=True, **kwargs)[source]#

Read spatial template model from FITS image. If unit is not given in the FITS header the default is sr-1.

Parameters
filenamestr

FITS image filename.

normalizebool

Normalize the input map so that it integrates to unity.

kwargsdict

Keyword arguments passed to Map.read().

reassign(datasets_names, new_datasets_names)#

Reassign a model from one dataset to another

Parameters
datasets_namesstr or list

Name of the datasets where the model is currently defined

new_datasets_namesstr or list

Name of the datasets where the model should be defined instead. If multiple names are given the two list must have the save length, as the reassignment is element-wise.

Returns
modelModel

Reassigned model.

to_dict(full_output=False)[source]#

Create dict for YAML serilisation

to_region(**kwargs)[source]#

Model outline from template map boundary (RectangleSkyRegion).

unfreeze()#

Restore parameters frozen status to default

write(overwrite=False)[source]#