void c_tbegtd (const IRAFPointer tp, const IRAFPointer cp, int row, double *buffer) {
/* Read a value of type double from a table column.
arguments:
IRAFPointer tp i: table descriptor
IRAFPointer cp i: column descriptor
int row i: row number (one indexed)
double *buffer o: value read from table
*/
TableDescr *tbl_descr;
ColumnDescr *col_descr;
int anynul=0;
long firstelem=1, nelem=1;
double nulval=IRAF_INDEFD;
int status = 0;
tbl_descr = (TableDescr *)tp;
col_descr = (ColumnDescr *)cp;
if (col_descr->datatype < 0) {
char *value;
int maxch;
maxch = col_descr->width + 5;
value = (char *)calloc (maxch, sizeof(char));
c_tbegtt (tp, cp, row, value, maxch);
if (strcmp (value, "INDEF") == 0)
*buffer = IRAF_INDEFD;
else
*buffer = atof (value);
free (value);
} else if (col_descr->datatype == IRAF_INT) {
int i_value;
c_tbegti (tp, cp, row, &i_value);
if (i_value == IRAF_INDEFI)
*buffer = IRAF_INDEFD;
else
*buffer = i_value;
} else if (col_descr->datatype == IRAF_SHORT) {
short si_value;
c_tbegts (tp, cp, row, &si_value);
if (si_value == IRAF_INDEFS)
*buffer = IRAF_INDEFD;
else
*buffer = si_value;
} else {
/* fits_read_col_dbl = ffgcvd */
fits_read_col_dbl (tbl_descr->fptr, col_descr->colnum,
(long)row, firstelem, nelem, nulval,
buffer, &anynul, &status);
if (status != 0)
setError (status, "c_tbegtd: error reading element");
}
}
void c_tbegti (const IRAFPointer tp, const IRAFPointer cp, int row, int *buffer) {
/* Read an integer value from a table column.
arguments:
IRAFPointer tp i: table descriptor
IRAFPointer cp i: column descriptor
int row i: row number (one indexed)
int *buffer o: value read from table
*/
int anynul=0;
long firstelem=1, nelem=1;
int nulval=IRAF_INDEFI;
int status = 0;
TableDescr *tbl_descr = (TableDescr *)tp;
ColumnDescr *col_descr = (ColumnDescr *)cp;
if (col_descr->datatype < 0) {
char *value;
int maxch;
maxch = col_descr->width + 5;
value = (char *)calloc (maxch, sizeof(char));
c_tbegtt (tp, cp, row, value, maxch);
if (strcmp (value, "INDEF") == 0)
*buffer = IRAF_INDEFI;
else
*buffer = atoi (value);
free (value);
} else if (col_descr->datatype == IRAF_BOOL) {
Bool b_value;
c_tbegtb (tp, cp, row, &b_value);
if (b_value)
*buffer = 1;
else
*buffer = 0;
} else if (col_descr->datatype == IRAF_SHORT) {
short si_value;
c_tbegts (tp, cp, row, &si_value);
if (si_value == IRAF_INDEFS)
*buffer = IRAF_INDEFI;
else
*buffer = si_value;
} else {
/* fits_read_col_int = ffgcvk */
fits_read_col_int (tbl_descr->fptr, col_descr->colnum,
(long)row, firstelem, nelem, nulval,
buffer, &anynul, &status);
if (status != 0)
setError (status, "c_tbegti: error reading element");
}
}
static int ReadBpixTab (TblInfo *tabinfo, int row, TblRow *tabrow) {
/* OPT_ELEM is an optional column. */
if (tabinfo->cp_opt_elem == 0) {
strcpy (tabrow->opt_elem, "ANY");
} else {
c_tbegtt (tabinfo->tp, tabinfo->cp_opt_elem, row,
tabrow->opt_elem, STIS_CBUF);
if (c_iraferr())
return (TABLE_ERROR);
}
c_tbegti (tabinfo->tp, tabinfo->cp_xstart, row, &tabrow->xstart);
if (c_iraferr())
return (TABLE_ERROR);
tabrow->xstart--;
c_tbegti (tabinfo->tp, tabinfo->cp_ystart, row, &tabrow->ystart);
if (c_iraferr())
return (TABLE_ERROR);
tabrow->ystart--;
c_tbegti (tabinfo->tp, tabinfo->cp_length, row, &tabrow->length);
if (c_iraferr())
return (TABLE_ERROR);
c_tbegti (tabinfo->tp, tabinfo->cp_axis, row, &tabrow->axis);
if (c_iraferr())
return (TABLE_ERROR);
c_tbegts (tabinfo->tp, tabinfo->cp_flag, row, &tabrow->flag);
if (c_iraferr())
return (TABLE_ERROR);
if (tabrow->axis !=1 && tabrow->axis != 2) {
printf ("ERROR Axis = %d in BPIXTAB, but it must be 1 or 2.\n",
tabrow->axis);
c_tbtclo (tabinfo->tp);
return (TABLE_ERROR);
}
if (tabrow->length <= 0) {
printf (
"ERROR Length = %d in BPIXTAB, but it must be positive.\n",
tabrow->length);
c_tbtclo (tabinfo->tp);
return (TABLE_ERROR);
}
return (0);
}
/* rejpar_in -- Read parameters either from user input or table.
Description:
------------
Reads CL parameters and does necessary checkings
Input parameters from crrej reference table:
-------------------------------------------
Col. Name Parameter
"skysub" sky Sky levels subtraction scheme
"crsigmas" sigmas Rejection thresholds
"crradius" radius Radius (in pixels) to propagate the cosmic ray
"crthresh" thresh Propagation factor
"initgues" initgues Scheme of computing initial-guess image
"scalense" scalense multiplicative noise in percents
"badinpdq" badinpdq Data quality pset
"crmask" mask flag CR-rejected pixels in input files?
Input parameters from input image primary header:
------------------------------------------------
Date Author Description
---- ------ -----------
24-Sep-1998 W.J. Hack Initial ACS Version
*/
int rejpar_in (clpar *par, int newpar[], int nimgs, float exptot, int *niter, float sigma[])
{
extern int status;
IRAFPointer tp;
IRAFPointer colptr, colptr1;
int i, nrows, nmatch, row;
int crsplit_in, crsplit, maxcrsplit;
float exp_in, meanexp, mindiff, diff;
char maskstr[ACS_CBUF+1];
void PrRefInfo (char *, char *, char *, char *, char *);
void WhichError (int);
/* -------------------------------- begin ---------------------------------- */
crsplit_in = nimgs;
exp_in = exptot / (float) crsplit_in;
par->meanexp = exp_in;
/* if all parameters are specified by the user, no need to open the
reference CRREJ table */
if (newpar[0] < MAX_PAR) {
tp = c_tbtopn (par->tbname, IRAF_READ_ONLY, 0);
if (c_iraferr() != 0) {
sprintf (MsgText,"CRREJTAB table '%s' does not exist", par->tbname);
trlerror (MsgText);
return (status = TABLE_ERROR);
}
nrows = c_tbpsta (tp, TBL_NROWS);
/* read the columns CRSPLIT and MEANEXP */
c_tbcfnd1 (tp, "crsplit", &colptr);
if (colptr == 0) {
trlerror ("column CRSPLIT does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbcfnd1 (tp, "meanexp", &colptr1);
if (colptr1 == 0) {
trlerror ("column MEANEXP does not exist in CRREJTAB\n");
return (status = COLUMN_NOT_FOUND);
}
nmatch = 0;
/* find the largest value in the CRSPLIT column */
for (i = 1; i <= nrows; i++) {
c_tbegti (tp, colptr, i, &crsplit);
if (i == 1) maxcrsplit = crsplit;
if (crsplit > maxcrsplit)
maxcrsplit = crsplit;
}
if (crsplit_in >= maxcrsplit) crsplit_in = maxcrsplit;
/* find the matching row in CRREJTAB to use */
for (i = 1; i <= nrows; i++) {
c_tbegti (tp, colptr, i, &crsplit);
c_tbegtr (tp, colptr1, i, &meanexp);
diff = meanexp - exp_in;
if (crsplit_in == crsplit && diff >= 0.) {
nmatch++;
if (nmatch == 1) mindiff = diff;
if (diff <= mindiff) {
row = i;
mindiff = diff;
}
}
}
if (nmatch == 0) {
trlerror (" No matching CRSPLIT and MEANEXP in CRREJTAB");
return (status = ROW_NOT_FOUND);
}
/* read the sigmas parameter */
if (newpar[CRSIGMAS] == 0) {
c_tbcfnd1 (tp, "crsigmas", &colptr);
if (colptr == 0) {
trlerror ("column CRSIGMAS does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegtt (tp, colptr, row, par->sigmas, ACS_LINE);
}
/* read other parameters */
if (newpar[SKYSUB] == 0) {
c_tbcfnd1 (tp, "skysub", &colptr);
if (colptr == 0) {
trlerror ("column SKYSUB does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegtt (tp, colptr, row, par->sky, ACS_FITS_REC);
}
/* CR propagation parameters */
if (newpar[CRRADIUS] == 0) {
c_tbcfnd1 (tp, "crradius", &colptr);
if (colptr == 0) {
trlerror ("column CRRADIUS does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegtr (tp, colptr, row, &par->radius);
}
if (newpar[CRTHRESH] == 0) {
c_tbcfnd1 (tp, "crthresh", &colptr);
if (colptr == 0) {
trlerror ("column CRTHRESH does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegtr (tp, colptr, row, &par->thresh);
}
/* figure out how to do initial comparison image */
if (newpar[INITGUES] == 0) {
c_tbcfnd1 (tp, "initgues", &colptr);
if (colptr == 0) {
trlerror ("column INITGUES does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegtt (tp, colptr, row, par->initgues, ACS_FITS_REC);
}
/* read the noise model */
if (newpar[SCALENSE] == 0) {
c_tbcfnd1 (tp, "scalense", &colptr);
if (colptr == 0) {
trlerror ("column SCALENSE does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegtr (tp, colptr, row, &par->scalense);
}
if (newpar[BADINPDQ] == 0) {
c_tbcfnd1 (tp, "badinpdq", &colptr);
if (colptr == 0) {
trlerror ("column BADINPDQ does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegts (tp, colptr, row, &par->badinpdq);
}
if (newpar[CRMASK] == 0) {
c_tbcfnd1 (tp, "crmask", &colptr);
if (colptr == 0) {
trlerror ("column CRMASK does not exist in CRREJTAB");
return (status = COLUMN_NOT_FOUND);
}
c_tbegti (tp, colptr, row, &par->mask);
}
c_tbtclo (tp);
}
PrRefInfo ("crrejtab", par->tbname, "", "", "");
/* parse the sigmas string into numbers */
*niter = strtor (par->sigmas, sigma);
if (status != ACS_OK) {
WhichError (status);
return (status);
}
if (*niter > MAX_ITER) {
sprintf (MsgText,"No more than %d iterations permitted.", MAX_ITER);
trlerror (MsgText);
return (status = ERROR_RETURN);
}
if (*niter <= 0) {
trlerror ("Number of iterations is ZERO.");
return (status = ERROR_RETURN);
}
/* other fixed (for now) parameters */
par->crval = (short) DATAREJECT;
par->fillval = 0.;
/* print out which parameter are used */
if (par->verbose) {
sprintf (MsgText,"\n number of images = %d", nimgs);
trlmessage (MsgText);
sprintf (MsgText," CRREJ ref table used: %s", par->tbname);
trlmessage (MsgText);
sprintf (MsgText," initial guess method: %s", par->initgues);
trlmessage (MsgText);
sprintf (MsgText," total exposure time = %0.1f", exptot);
trlmessage (MsgText);
sprintf (MsgText," sigmas used: %s", par->sigmas);
trlmessage (MsgText);
sprintf (MsgText," sky subtraction used: %s", par->sky);
trlmessage (MsgText);
sprintf (MsgText," rejection radius = %0.1f", par->radius);
trlmessage (MsgText);
sprintf (MsgText," propagation threshold = %0.1f", par->thresh);
trlmessage (MsgText);
sprintf (MsgText," scale noise = %0.1f%%", par->scalense);
trlmessage (MsgText);
sprintf (MsgText," input bad bits value = %d", par->badinpdq);
trlmessage (MsgText);
if (par->mask == 1) {
strcpy (maskstr,"YES");
} else {
strcpy (maskstr, "NO");
}
sprintf (MsgText," reset crmask = %s\n", maskstr);
trlmessage (MsgText);
}
return (status);
}
void c_tbegtt (const IRAFPointer tp, const IRAFPointer cp, int row, char *buffer,
int maxch) {
/* Read a string value from a table column.
arguments:
IRAFPointer tp i: table descriptor
IRAFPointer cp i: column descriptor
int row i: row number (one indexed)
char *buffer o: value read from table
int maxch i: maximum length of the string (not incl NULL)
*/
TableDescr *tbl_descr;
ColumnDescr *col_descr;
int anynul=0;
long firstelem=1, nelem=1;
char *value;
int len;
int status = 0;
tbl_descr = (TableDescr *)tp;
col_descr = (ColumnDescr *)cp;
if (col_descr->datatype < 0) {
if (col_descr->width >= maxch)
len = col_descr->width;
else
len = maxch;
} else {
len = SZ_FITS_STR;
}
value = (char *)calloc (len+1, sizeof(char));
if (col_descr->datatype == IRAF_BOOL) {
Bool b_value;
c_tbegtb (tp, cp, row, &b_value);
if (b_value == True)
strcpy (value, "yes");
else
strcpy (value, "no");
} else if (col_descr->datatype == IRAF_DOUBLE) {
double d_value;
c_tbegtd (tp, cp, row, &d_value);
if (d_value == IRAF_INDEFD)
strcpy (value, "INDEF");
else
sprintf (value, "%.16g", d_value);
} else if (col_descr->datatype == IRAF_REAL) {
float r_value;
c_tbegtr (tp, cp, row, &r_value);
if (r_value >= 0.99999 * IRAF_INDEFR &&
r_value <= 1.00001 * IRAF_INDEFR) {
strcpy (value, "INDEF");
} else {
sprintf (value, "%.7g", r_value);
}
} else if (col_descr->datatype == IRAF_INT) {
int i_value;
c_tbegti (tp, cp, row, &i_value);
if (i_value == IRAF_INDEFI)
strcpy (value, "INDEF");
else
sprintf (value, "%d", i_value);
} else if (col_descr->datatype == IRAF_SHORT) {
short si_value;
c_tbegts (tp, cp, row, &si_value);
if (si_value == IRAF_INDEFS)
strcpy (value, "INDEF");
else
sprintf (value, "%hd", si_value);
} else {
/* fits_read_col_str = ffgcvs */
fits_read_col_str (tbl_descr->fptr, col_descr->colnum,
(long)row, firstelem, nelem, "INDEF",
&value, &anynul, &status);
if (status != 0)
setError (status, "c_tbegtt: error reading element");
}
copyString (buffer, value, maxch);
free (value);
}
void c_tbegtb (const IRAFPointer tp, const IRAFPointer cp, int row, Bool *buffer) {
/* Read a boolean value from a table column.
arguments:
IRAFPointer tp i: table descriptor
IRAFPointer cp i: column descriptor
int row i: row number (one indexed)
Bool *buffer o: value (True or False) read from table
*/
TableDescr *tbl_descr;
ColumnDescr *col_descr;
int anynul=0;
long firstelem=1, nelem=1;
int nulval=0;
char s_value[11]={'\0'};
int status = 0;
tbl_descr = (TableDescr *)tp;
col_descr = (ColumnDescr *)cp;
if (col_descr->datatype < 0) {
char *value;
int maxch;
maxch = col_descr->width + 5;
value = (char *)calloc (maxch, sizeof(char));
c_tbegtt (tp, cp, row, value, maxch);
str_lower (value, value);
if (strcmp (value, "1") == 0) {
*buffer = True;
} else if (strcmp (value, "0") == 0) {
*buffer = False;
} else if (strncmp (value, "yes", 3) == 0 ||
strncmp (value, "true", 4) == 0) {
*buffer = True;
} else {
*buffer = False;
}
free (value);
} else if (col_descr->datatype == IRAF_INT) {
int i_value;
c_tbegti (tp, cp, row, &i_value);
if (i_value == IRAF_INDEFI)
*buffer = False;
else if (i_value)
*buffer = True;
else
*buffer = False;
} else if (col_descr->datatype == IRAF_SHORT) {
short si_value;
c_tbegts (tp, cp, row, &si_value);
if (si_value == IRAF_INDEFS)
*buffer = False;
else if (si_value)
*buffer = True;
else
*buffer = False;
} else {
/* fits_read_col_log = ffgcvl */
fits_read_col_log (tbl_descr->fptr, col_descr->colnum,
(long)row, firstelem, nelem, nulval,
s_value, &anynul, &status);
if (status != 0)
setError (status, "c_tbegtb: error reading element");
if (s_value[0])
*buffer = True;
else
*buffer = False;
}
}
