/*
 				field.h

*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
*
*	Part of:	SWarp
*
*	Author:		E.BERTIN (IAP)
*
*	Contents:	Handling of field structures.
*
*	Last modify:	07/03/2001
*
*%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
*/

#ifndef _FITSCAT_H_
#include "fitscat.h"
#endif

#ifndef _FIELD_H_
#define _FIELD_H_

#define		MAXINFIELD	8192	/* Maximum number of input files */

/*------------------------------ field flags --------------------------------*/
#define		DETECT_FIELD	0x01	/* Detection */
#define		MEASURE_FIELD	0x02	/* Measurement */
#define		FLAG_FIELD	0x04	/* Flagging */
#define		RMS_FIELD	0x08	/* Weighting with std deviations */
#define		VAR_FIELD	0x10	/* Weighting with variances */
#define		WEIGHT_FIELD	0x20	/* Weighting with weights */
#define		BACKRMS_FIELD	0x40	/* Weighting from a backrms matrix */

#define		FIELD_READ	0x100	/* Field available for read access */
#define		FIELD_WRITE	0x200	/* Field available for read access */

/*--------------------------------- typedefs --------------------------------*/
typedef enum	{BACK_RELATIVE, BACK_ABSOLUTE}	backenum;

typedef struct field
  {
  char		filename[MAXCHAR];	/* image filename */
  char		*rfilename;		/* pointer to the reduced image name */
  char		hfilename[MAXCHAR];	/* header filename */
  char		ident[MAXCHAR];		/* field identifier (read from FITS) */
  catstruct	*cat;			/* cat structure */
  tabstruct	*tab;			/* tab structure */
/* ---- main image parameters */
  int		frameno;		/* pos in Multi-extension FITS file */
  int		width, height;		/* x,y size of the field */
  size_t	npix;			/* total number of pixels */
  double	gain;			/* conversion factor e-/ADU */
  double	fscale;			/* relative photometric scale */
  double	ngamma;			/* normalized photo gamma */
/* ---- background parameters */
  float		*back;			/* ptr to the background map in mem */
  float		*dback;			/* ptr to the background deriv. map */
  float		*sigma;			/* ptr to the sigma map */
  float		*dsigma;		/* Ptr to the sigma deriv. map */
  int		backw, backh;		/* x,y size of a bkgnd mesh */
  int		nbackp;			/* total nb of pixels per bkgnd mesh */
  int		nbackx, nbacky;		/* x,y number of bkgnd meshes */
  int		nback;			/* total number of bkgnd meshes */
  int		nbackfx, nbackfy;	/* x,y size of bkgnd filtering mask */
  float		backdefault;		/* default background value */
  float		backmean;		/* median bkgnd value in image */
  float		backsig;		/* median bkgnd rms in image */
  float		sigfac;			/* scaling RMS factor (for WEIGHTs) */
  PIXTYPE	*pix;			/* pixel data */
  PIXTYPE	*backline;		/* current interpolated bkgnd line */
  backenum     	back_type;		/* background type */
/* ---- dynamic compression parameters */
  PIXTYPE	dyncomp_fac;		/* dynamic compression factor */
/* ---- astrometric parameters */
  struct structwcs	*wcs;		/* astrometric data */
  int		flags;			/* flags defining the field type */
  int		cflags;			/* flags defining the conversion */
  double	pixscale;		/* pixel size in arcsec.pix-1 */
  double	epoch;			/* epoch of coordinates */
  double	*rawmin;		/* Starting pixel for coaddition */
  double	*rawmax;		/* Ending pixel for coaddition */
/* ---- image interpolation */
  PIXTYPE	weight_thresh;		/* interpolation threshold */
  struct field	*reffield;	       	/* pointer to a reference field */
  }	fieldstruct;

/*------------------------------- functions ---------------------------------*/

extern fieldstruct	*inherit_field(char *filename, fieldstruct *reffield,
					int fflags),
			*init_field(fieldstruct **infield,
				fieldstruct **inwfield, int ninput,
				char *filename),
			*load_field(char *filename, int frameno, int fieldno);

extern int		read_aschead(char *filename, int frameno,
					tabstruct *tab);

extern void		end_field(fieldstruct *field),
			printinfo_field(fieldstruct *field,
					fieldstruct *wfield),
			scale_field(fieldstruct *field, fieldstruct *reffield),
			writefitsinfo_field(fieldstruct *field);

#endif