Skip to content

Latest commit

 

History

History
47 lines (30 loc) · 2.32 KB

README.md

File metadata and controls

47 lines (30 loc) · 2.32 KB
Raw

pymosaics

PyPI version Build Status codecov Total alerts Language grade: Python

CAUTION: this is a simple functional project, although it should work as generally expected, do make sanity checks and make sure you understand what the mosaic does - as usual. In particular, the re-pixelization approach (mapping/interpolation) is suitable for maps with flux density, as long as the pixels are smaller than PSF

Merge FITS file mosaic

  • maps all images in the same pixels (taken from one of the images, or healpix)

  • sums with a flexible algorithm

  • output mosaics format is compatible with input - allowing map-reduce operations

    note that pixelization used for mosaic and for the output is not necessarily the same, i.e. for healpix it's useful to have "regular" projection output.

    clearly, healpix pixels are useful primarily for all-sky.

Installation

$ pip install pymosaic-fits

Example:

$ mosaic tests/data/isgri_sky_ima_{1,2}.fits.gz out.fits # just a normal mosaic, pixels/output from first image
$ mosaic tests/data/isgri_sky_ima_{1,2}.fits.gz out.fits --mock # mock, to show assembly
$ mosaic tests/data/isgri_sky_ima_{1,2}.fits.gz out.fits --pixels healpix # healpix pixels, projected output
import mosaic                                                            
mosaic.mosaic_fn_list(["tests/data/isgri_sky_ima_1.fits.gz", "tests/data/isgri_sky_ima_2.fits.gz"], "out.fits")

Additional mosaic analysis

Using https://github.com/astromatic/sextractor/ we provide some analysis, suitable for example for INTEGRAL ISGRI and JEMX mosaics:

image