ERIS Pipeline Reference Manual 1.8.10
eris_ifu_extract_spec_static.h
1/* $Id$
2 *
3 * This file is part of the ERIS Pipeline
4 * Copyright (C) 2002,2003 European Southern Observatory
5 *
6 * This program is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation; either version 2 of the License, or
9 * (at your option) any later version.
10 *
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
15 *
16 * You should have received a copy of the GNU General Public License
17 * along with this program; if not, write to the Free Software
18 * Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
19 */
20
21#ifndef ERISP_ERIS_ERIS_IFU_EXTRACT_SPEC_STATIC_H_
22#define ERISP_ERIS_ERIS_IFU_EXTRACT_SPEC_STATIC_H_
23
24#ifdef HAVE_CONFIG_H
25#include <config.h>
26#endif
27
28/*-----------------------------------------------------------------------------
29 Includes
30 -----------------------------------------------------------------------------*/
31
32#include <time.h>
33#include <math.h>
34
35#include <cpl.h>
36#include <hdrl.h>
37#include "eris_ifu_functions.h"
38#include "eris_ifu_vector.h"
39
40/*-----------------------------------------------------------------------------
41 Define
42 -----------------------------------------------------------------------------*/
43
44typedef enum {MASK, POSITION, MAX, FIT, OPTIMAL} mask_method_types;
45struct esParamStruct {
46 mask_method_types mask_method;
47 int center_x;
48 int center_y;
49 double radius;
50 productDepthType productDepth;
51 bool save_mask;
52};
53
54
55struct esSofStruct {
56 hdrl_imagelist *cube;
57 cpl_image *mask;
58// cpl_image *collapsedCube;
59 cpl_propertylist *header;
60 cpl_imagelist *qualImagelist;
61 deqQualityType qualityType;
62 cpl_size nx;
63 cpl_size ny;
64 cpl_size nz;
65};
66
67
68/*-----------------------------------------------------------------------------
69 Functions prototypes
70 -----------------------------------------------------------------------------*/
71
72/*----------------------------------------------------------------------------*/
73void eris_ifu_extract_free_esSofStruct(struct esSofStruct* self);
74
75hdrl_image* eris_ifu_extract_spec_collapse(hdrl_imagelist *cube,
76 cpl_image **contribute);
77
78cpl_image* eris_ifu_extract_spec_create_fit_mask(const hdrl_image *img);
79
80cpl_image* eris_ifu_extract_spec_create_circle_mask(
81 cpl_size center_x,
82 cpl_size center_y,
83 double radius,
84 cpl_size nx,
85 cpl_size ny);
86
87cpl_image* eris_ifu_extract_spec_create_mask(struct esParamStruct params,
88 struct esSofStruct sof,
89 const hdrl_image *collapsedCube,
90 int productDepth);
91
92cpl_bivector* eris_ifu_extract_spectrum(
93 hdrl_imagelist *cube,
94 cpl_image *mask,
95 double startLambda,
96 double deltaLambda,
97 cpl_vector **error,
98 cpl_vector **totalFlux);
99
100//cpl_bivector* eris_ifu_extract_spectrum2(
101// hdrl_imagelist *cube,
102// cpl_image *mask,
103// double startLambda,
104// double deltaLambda,
105// cpl_vector **error,
106// cpl_vector **totalFlux);
107
108//cpl_bivector* eris_ifu_extract_spectrum3(
109// hdrl_imagelist *cube,
110// cpl_image *mask,
111// double startLambda,
112// double deltaLambda,
113// cpl_vector **error,
114// cpl_vector **totalFlux);
115
116cpl_bivector* eris_ifu_optimal_extraction(const hdrl_imagelist *cube,
117 const cpl_imagelist *cube_dqi,
118 const cpl_image *img_mask,
119 double startLambda,
120 double deltaLambda,
121 int productDepth,
122 cpl_vector **error_out);
123
124cpl_error_code eris_ifu_opt_extr_get_center_fwhm(const hdrl_image *hdrl_img,
125 int edge_trim,
126 cpl_size *xcen,
127 cpl_size *ycen,
128 double *fwhm);
129
130cpl_mask* eris_ifu_opt_extr_create_mask(int nx,
131 int ny,
132 int xcen,
133 int ycen,
134 double fwhm);
135
136cpl_error_code eris_ifu_opt_extr_simple_extraction(const hdrl_imagelist *cube,
137 const cpl_mask *mask,
138 cpl_vector **spec,
139 cpl_vector **spec_var);
140
141cpl_bivector * eris_ifu_opt_extr_helper_usepix(const cpl_mask *mask, int *n_usepix);
142
143cpl_error_code eris_ifu_opt_extr_helper_fill_vertical(cpl_image *img, const eris_ifu_vector *vec);
144
145void eris_ifu_opt_extr_helper_fill_horizontal(cpl_image *img, const cpl_image *slice, const cpl_bivector *usepix, int row, int power);
146
147void eris_ifu_opt_extr_helper_set_positive(cpl_image *img);
148
149cpl_vector* eris_ifu_opt_extr_get_col(const cpl_image *img, int colnr);
150
151cpl_vector *eris_ifu_opt_extr_get_row(const cpl_image *img, int rownr);
152
153cpl_error_code eris_ifu_opt_extr_set_row(cpl_image *img,
154 int rownr,
155 const cpl_vector *vec);
156
157cpl_error_code eris_ifu_opt_extr_set_col(cpl_image *img,
158 int colnr,
159 const cpl_vector *vec);
160
161cpl_vector* eris_ifu_opt_extr_create_lambda(int size,
162 double startLambda,
163 double deltaLambda);
164
165cpl_error_code eris_ifu_opt_extr_vector_sqrt(cpl_vector *vec);
166
167cpl_bivector* eris_ifu_opt_extr_doit(const hdrl_imagelist *cube,
168 const cpl_imagelist *cube_dqi,
169 const cpl_mask *mask,
170 const eris_ifu_vector *spec,
171 const eris_ifu_vector *spec_var,
172 double startLambda,
173 double deltaLambda,
174 int productDepth,
175 cpl_vector **error_out);
176
177int eris_ifu_opt_extr_get_not_finite(cpl_vector *vec);
178
179void eris_ifu_opt_extr_convert_0_to_NaN_vec(cpl_vector *vec);
180
181double eris_ifu_opt_extr_estimate_radius(const cpl_image *img,
182 cpl_size center_x,
183 cpl_size center_y,
184 double initial_radius,
185 int productDepth);
186
187//double eris_ifu_opt_extr_median_without_NaN(cpl_vector *vec);
188
189cpl_error_code eris_ifu_cpl_vector_sqrt(cpl_vector *vec);
190
191#endif
192
193
Generated on Fri Sep 19 2025 20:45:25 for ERIS Pipeline Reference Manual by doxygen 1.9.6