#include <math.h>
#include <xsh_drl.h>
#include <xsh_drl_check.h>
#include <xsh_badpixelmap.h>
#include <xsh_data_pre.h>
#include <xsh_data_order.h>
#include <xsh_data_wavemap.h>
#include <xsh_data_localization.h>
#include <xsh_data_rec.h>
#include <xsh_dfs.h>
#include <xsh_pfits.h>
#include <xsh_error.h>
#include <xsh_utils_image.h>
#include <xsh_msg.h>
#include <xsh_fit.h>
#include <cpl.h>
#include <gsl/gsl_bspline.h>
#include <gsl/gsl_linalg.h>
#include <f2c.h>
#include <clapack.h>
Go to the source code of this file.
Defines | |
#define | REGDEBUG_MEDIAN_SPLINE 0 |
Functions | |
static void | fit_spline (xsh_wavemap_list *wave_list, int idx, int nbkpts, int order, int niter, float kappa, float ron2, float gain, const char *sarm) |
static xsh_wavemap_list * | xsh_wavemap_list_new (cpl_image *wavemap, cpl_image *slitmap, xsh_localization *list, xsh_order_list *order_table, xsh_pre *pre_sci, int nbkpts, cpl_frame *usr_defined_break_points_frame, xsh_subtract_sky_single_param *sky_par, xsh_instrument *instrument) |
static xsh_rec_list * | xsh_wavemap_list_build_sky (xsh_wavemap_list *wave_list, xsh_pre *pre_mflat, xsh_instrument *instr) |
Build the sky spectrum list. | |
static int | wave_compare (const void *un, const void *deux) |
static cpl_frame * | xsh_wavelist_subtract_sky (xsh_pre *pre_sci, xsh_pre *pre_mflat, xsh_rec_list *sky_list, xsh_wavemap_list *wave_list, xsh_instrument *instrument, const char *rec_prefix) |
subtract sky | |
static int | xsh_bspline_eval_all (const double x, gsl_vector *B, size_t *idx, gsl_bspline_workspace *w, size_t *start_i) |
static size_t | xsh_bspline_find_interval (const double x, int *flag, gsl_bspline_workspace *w, size_t *start_i) |
static cpl_error_code | xsh_mflat_normalize (xsh_order_list *order_table, xsh_pre *pre_mflat) |
Normalize the master flat to the value at the center of the slit on each order and Y position. | |
cpl_frame * | xsh_subtract_sky_single (cpl_frame *sci_frame, cpl_frame *order_table_frame, cpl_frame *slitmap_frame, cpl_frame *wavemap_frame, cpl_frame *loc_table_frame, cpl_frame *mflat_frame, cpl_frame *usr_defined_break_points_frame, xsh_instrument *instrument, int nbkpts, xsh_subtract_sky_single_param *sky_par, cpl_frame **sky_spectrum, cpl_frame **sky_spectrum_eso, const char *rec_prefix, const int clean_tmp) |
Subtract the sky background for single frame. If sky_spectrum is NOT NULL it is saved as a product, of type REC. | |
cpl_frame * | xsh_save_sky_model (cpl_frame *obj_frame, cpl_frame *sub_sky_frame, const char *sky_tag, xsh_instrument *instrument) |
cpl_frame * | xsh_add_sky_model (cpl_frame *subsky_frame, cpl_frame *sky_frame, xsh_instrument *instrument, const char *prefix) |
Variables | |
double | xsh_nbkpts_uvb [XSH_ORDERS_UVB] |
Create the wavemap list. | |
double | xsh_nbkpts_vis [XSH_ORDERS_VIS] |
double | xsh_nbkpts_nir [XSH_ORDERS_NIR] |