Domenico Tiziani / gammalib

Added GWcsSIN class.

# Conflicts:
#	ChangeLog
#	NEWS

In the introduction of the GSkyMap::shape() method, a change was introduced that led to writing single maps as 3D images. The code was modified so that the original behaviour of writing 2D images was recovered.

This replaces two arguments:
- naxis: the dimension of the image
- naxes: a pointer to an array containing the number of pixels for each dimension

…id pointer being free'd (not clear why this is the case)

The shape() and ndim() methods allow the handling of images with more than three dimensions in GSkyMap. Now higher dimensional images are properly read and written by the class, and the shape() method allows setting of higher dimensional map parameters.

…udes a modification to GSkyRegionCircle class

Add a warning that the method should not be used during setup.

Before, always the first map has been accessed.

The check done in GSkymap::contains() was not correct, and pixel indices that were beyond the edge of the map by one pixel were still accepted. This lead in some cases to artefacts when using the GSkymap::operator() sky direction access operator. This is now corrected.

The flux() method now computes an approximated flux value for HealPix projections. This flux estimate is accurate to 1-3% level, dependent on pixel and Nside parameter.

The following methods have been added:
- set_z_phi
- loc2dir
- xyf2loc
- pix2xyf
- max_pixrad
- boundaries
- neighbours

The wrap around in the X axis is now correctly dealt with in the interpolation operator.

Upon reading of the WCS map the size of the map is checked. Exceptions are thrown in case that the map is larger then 360 in X and larger than 180 in Y.

The flux() methods return the correct flux in a pixel considering that a bilinear interpolation is performed in the map. Flux is defined as intensity times solid angle.

This works so far only properly for WCS projections. For HealPix maps, only the pixel intensity times the solid angle is returned.

This exception may occur if the skymap projection cannot represent the specified sky coordinate.