Utilities for galfitm (mostly python)

##galfitm-visualisations

Work in progress...

import galfitm_visualisations as gv

A couple of examples...

Figures are output in a folder named plots in the folder from which the code is run, so you may need to create this first.

Create a pdf figure of the input, model and residual from a single-band galfitm output fits file. Automatic scaling should be reasonable, but you can also specify various options manually.

gv.plotimg('imgblock.fits')

To include nonparametric output,

gv.plotimg('imgblock.fits', nonparam=1)

If you have several different fits you want to compare, say single-sersic, bulge+disk, bulge+disk+bar:

gv.plotimg(('ss.fits', 'bd.fits', 'bdb.fits'))

If two of these have nonparametric images enabled:

gv.plotimg(('ss.fits', 'bd.fits', 'bdb.fits'), nonparam=2)

If you have performed multi-band fits, then just specify the band labels you want to plot:

gv.plotimg(('ss.fits', 'bd.fits', 'bdb.fits'), bands=('g', 'r', 'i'))

To make an equaivalent figure based on single-band fits, you can specify the location of the band name in the filename with {}:

gv.plotimg(('ss_{}.fits', 'bd_{}.fits', 'bdb_{}.fits'), bands=('g', 'r', 'i'))

This functionality can be repurposed to make a single figure for several different galaxies. Say you have a list of object IDs in ids, corresponding to a set of normal (named fits/1C_sersic/1C_{}.fits) and nonparametric (fits/1C_sersic/1C_{}.fits) galfitm output FITS files for single-band fits to SDSS galaxies. Then you can could create a nice figure with:

gv.plotimg(('fits/1C_sersic/1C_{}.fits', 'fits/1C_sersic/1C_{}_nonparam.fits'),
                  bands=ids, zoom=2, uservmin=-10, uservmax=1000, usernorm=20,
                  uservrange=1000, nonparam=1, label=True, ids=('normal', 'nonparam'), 
                  name='1C_sersic')