GRAVI Pipeline Reference Manual 1.8.0
Loading...
Searching...
No Matches
gravi_cpl.h
Go to the documentation of this file.
1/* $Id: gravi_cpl.h,v 1.12 2014/11/12 06:10:40 nazouaoui Exp $
2 *
3 * This file is part of the GRAVI 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., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
19 */
20
21#ifndef GRAVI_CPL_H_
22#define GRAVI_CPL_H_
23
24/*-----------------------------------------------------------------------------
25 Includes
26 -----------------------------------------------------------------------------*/
27
28#include <cpl.h>
29
30/*-----------------------------------------------------------------------------
31 Prototypes
32 -----------------------------------------------------------------------------*/
33
34cpl_array ** gravi_array_new_list (int n, cpl_type type, int size);
35
36double ** gravi_table_get_data_array_double (cpl_table * table, const char * name);
37float ** gravi_table_get_data_array_float (cpl_table * table, const char * name);
38float complex ** gravi_table_get_data_array_float_complex (cpl_table * table, const char * name);
39double complex ** gravi_table_get_data_array_double_complex (cpl_table * table, const char * name);
40int ** gravi_table_get_data_array_int (cpl_table * table, const char * name);
41
42double gravi_table_get_column_mean (cpl_table * table, const char * name, int base, int nbase);
43double gravi_table_get_column_flagged_mean (cpl_table * table, const char * name);
44double gravi_table_get_column_std (cpl_table * table, const char * name, int base, int nbase);
45cpl_array * gravi_table_get_column_sum_array (cpl_table * table, const char * name, int base, int nbase);
46cpl_array * gravi_table_get_column_mean_array (cpl_table * table, const char * name, int base, int nbase);
47
48#define gravi_table_get_value(table,name,row,value) cpl_array_get (cpl_table_get_array (table,name,row), value, NULL)
49#define gravi_table_set_value(table,name,row,value,val) cpl_array_set (cpl_table_get_data_array (table,name)[row], value, val);
50
51int gravi_array_threshold_min (cpl_array * array, double lo_cut);
52
53cpl_array * gravi_array_init_double (long n , double value);
54cpl_array * gravi_array_init_int (long n, int value);
55cpl_array * gravi_array_init_double_complex (long n, double complex value);
56cpl_array * gravi_array_init_float_complex (long n, float complex value);
57
58cpl_array * gravi_array_wrap_complex (cpl_array * input_re, cpl_array * input_im);
59cpl_array * gravi_array_wrap_float_complex (cpl_array * input_re, cpl_array * input_im);
60cpl_array * gravi_array_compute_norm2 (cpl_array * input_re, cpl_array * input_im);
61cpl_array * gravi_array_smooth (cpl_array * input, int nsmooth);
62cpl_array * gravi_array_rebin (const cpl_array * input, const cpl_array * errs,
63 cpl_table * oi_wave_sc, cpl_table * oi_wave_ft);
64cpl_array * gravi_array_create_inverse (cpl_array *input);
65cpl_array * gravi_array_cexp (double complex factor, const cpl_array * input);
66
67cpl_error_code gravi_array_normalize_complex (cpl_array * input);
68cpl_error_code gravi_array_phase_unwrap (cpl_array * input);
69cpl_error_code gravi_array_phase_wrap (cpl_array * input);
70cpl_error_code gravi_array_multiply_phasor (cpl_array * input, double complex factor, cpl_array * phase);
71cpl_error_code gravi_array_add_phasor (cpl_array * input, double complex factor, cpl_array * phase);
72cpl_error_code gravi_array_add_phase (cpl_array * input, double factor, cpl_array * phase);
73cpl_error_code gravi_array_multiply_conj (cpl_array * input1, cpl_array * input2);
74cpl_error_code gravi_array_add_phasors (cpl_array * input, cpl_array * add, cpl_array * sub);
75
76cpl_error_code gravi_table_set_array_double_complex (cpl_table * table, const char * name,
77 cpl_size row, cpl_array * visR, cpl_array * visI);
78cpl_error_code gravi_table_set_array_phase (cpl_table * table, const char * name,
79 cpl_size row, cpl_array * phase);
80cpl_error_code gravi_table_set_string_fixlen (cpl_table *table, const char *name,
81 int row, const char *value, int len);
82int gravi_table_are_equal (cpl_table * first, cpl_table * second);
83
84cpl_error_code gravi_table_new_column (cpl_table * table, const char * name, const char * unit, cpl_type type);
85cpl_error_code gravi_table_new_column_array (cpl_table * table, const char * name, const char * unit, cpl_type type, cpl_size size);
86cpl_error_code gravi_table_init_column_array (cpl_table * table, const char * name, const char * unit, cpl_type type, cpl_size size);
87
88cpl_error_code gravi_table_interpolate_column (cpl_table * to_table,
89 const char * to_x,
90 const char * to_y,
91 const cpl_table * from_table,
92 const char * from_x,
93 const char * from_y);
94
95cpl_error_code gravi_table_add_columns (cpl_table * oi_vis1, const char *name1,
96 cpl_table * oi_vis2, const char *name2);
97cpl_error_code gravi_table_runint_column (cpl_table * oi_vis, const char *input_name,
98 const char *output_name, int nsmooth, int nbase);
99cpl_error_code gravi_table_smooth_column (cpl_table * oi_vis, const char *input_name,
100 const char *output_name, int nsmooth, int nbase);
101
102cpl_error_code gravi_table_multiply_scalar (cpl_table * table, const char * name, int base, int nbase, double value);
103cpl_error_code gravi_table_add_scalar (cpl_table * table, const char * name, int base, int nbase, double value);
104
105cpl_array * gravi_table_create_sigma_array (cpl_table * oi_wave);
106cpl_array * gravi_table_create_wave_array (cpl_table * oi_wave);
107
108cpl_error_code gravi_array_get_group_delay_loop (cpl_array ** input,
109 cpl_array ** flag,
110 cpl_array * sigma,
111 double * gd, cpl_size nrow,
112 double max_width, int verbose);
113
114cpl_error_code gravi_table_compute_group_delay (cpl_table * table,
115 const char *input,
116 const char *flag,
117 const char *output, cpl_table * oi_wave);
118
119cpl_matrix * get_matrix_from_vector (cpl_vector * , cpl_vector * );
120
121cpl_imagelist * gravi_imagelist_from_column (cpl_table * table_data, const char * data_x);
122cpl_imagelist * gravi_imagelist_wrap_column (cpl_table * table_data, const char * data_x);
123cpl_error_code gravi_imagelist_unwrap_images (cpl_imagelist * imglist);
124double gravi_image_get_noise_window (cpl_image *img,
125 cpl_size llx, cpl_size lly,
126 cpl_size urx, cpl_size ury);
127
128cpl_vector * gravi_table_get_vector (cpl_table * , cpl_size , const char * );
129cpl_vector * gravi_table_get_vector_scalar (cpl_table * table,
130 const char * name,
131 cpl_size base,
132 cpl_size nbase);
133
134cpl_vector * gravi_table_get_vector_diff (cpl_table * spectrum_data, int index,
135 const char * regname1,
136 const char * regname2);
137
138cpl_array * gravi_array_wrap_image (cpl_image *);
139
140cpl_image * gravi_image_collapse_median_x (cpl_image * img,
141 cpl_size drop_from,
142 cpl_size drop_to);
143
144cpl_error_code gravi_image_subtract_collapse (cpl_image * img,
145 const cpl_image * collapse,
146 int direction);
147
148cpl_error_code gravi_image_subtract_window (cpl_image * img1, const cpl_image * img2,
149 cpl_size llx, cpl_size lly,
150 cpl_size urx, cpl_size ury,
151 cpl_size llx2, cpl_size lly2);
152
153cpl_error_code gravi_image_replace_window (cpl_image * img1, const cpl_image * img2,
154 cpl_size llx, cpl_size lly,
155 cpl_size urx, cpl_size ury,
156 cpl_size llx2, cpl_size lly2);
157
158cpl_error_code gravi_image_fill (cpl_image * img, double value);
159cpl_image * gravi_image_wrap_matrix (cpl_matrix * matrix);
160cpl_image * gravi_image_from_matrix (cpl_matrix * matrix);
161cpl_image * gravi_image_wrap_vector (cpl_vector * vector);
162cpl_image * gravi_image_from_vector (cpl_vector * vector);
163cpl_image * gravi_image_from_column (cpl_table * table_data, const char * data_x, cpl_size row);
164
165double gravi_image_get_quantile (const cpl_image * img, double value);
166double gravi_array_get_quantile (cpl_array * arr, double value);
167
168cpl_vector * gravi_vector_median (const cpl_vector * vector, cpl_size hw);
169cpl_error_code gravi_vector_abs (cpl_vector * vector);
170cpl_size gravi_vector_get_maxpos (cpl_vector * vector);
171cpl_vector * gravi_vector_extract (const cpl_vector * vector, int start, int step);
172cpl_error_code gravi_vector_unwrap_with_guess (cpl_vector * vector, cpl_vector * ref, double ref_to_phase);
173
174double gravi_vector_get_mean_clip (cpl_vector * vector_in,
175 double percent, double nsigma);
176
177cpl_array * gravi_table_get_column_dimension (const cpl_table * table, const char * name);
178
179cpl_matrix * gravi_matrix_interpolate_col (cpl_matrix * matrix,
180 cpl_vector * xref,
181 cpl_vector * xout);
182
183cpl_matrix * gravi_matrix_invertSV_create (cpl_matrix * a_in);
184
185cpl_table * gravi_table_extract_time_interval (cpl_table *table, double start, double end);
186
187
188#endif /* GRAVI_CPL_H_ */
189
190
cpl_array * gravi_array_init_int(long n, int value)
Definition: gravi_cpl.c:540
cpl_error_code gravi_table_set_array_double_complex(cpl_table *table, const char *name, cpl_size row, cpl_array *visR, cpl_array *visI)
Definition: gravi_cpl.c:646
cpl_array * gravi_array_smooth(cpl_array *input, int nsmooth)
Definition: gravi_cpl.c:1217
cpl_table * gravi_table_extract_time_interval(cpl_table *table, double start, double end)
Extract rows from table based on the TIME column.
Definition: gravi_cpl.c:3193
cpl_error_code gravi_table_smooth_column(cpl_table *oi_vis, const char *input_name, const char *output_name, int nsmooth, int nbase)
Definition: gravi_cpl.c:936
cpl_error_code gravi_table_add_columns(cpl_table *oi_vis1, const char *name1, cpl_table *oi_vis2, const char *name2)
Definition: gravi_cpl.c:772
float complex ** gravi_table_get_data_array_float_complex(cpl_table *table, const char *name)
Definition: gravi_cpl.c:453
cpl_error_code gravi_table_runint_column(cpl_table *oi_vis, const char *input_name, const char *output_name, int nsmooth, int nbase)
Definition: gravi_cpl.c:843
cpl_vector * gravi_vector_extract(const cpl_vector *vector, int start, int step)
Extract part of a vector.
Definition: gravi_cpl.c:2628
double gravi_table_get_column_flagged_mean(cpl_table *table, const char *name)
Function to compute the mean of a column table with arrays.
Definition: gravi_cpl.c:217
double gravi_image_get_noise_window(cpl_image *img, cpl_size llx, cpl_size lly, cpl_size urx, cpl_size ury)
Estimate the median noise in a window. This noise is estimated as the median{img**2}**0....
Definition: gravi_cpl.c:1771
cpl_array * gravi_array_create_inverse(cpl_array *input)
Definition: gravi_cpl.c:1003
cpl_array * gravi_table_get_column_dimension(const cpl_table *table, const char *name)
Return an array ready for cpl_table_set_column_dimension.
Definition: gravi_cpl.c:2046
double gravi_table_get_column_mean(cpl_table *table, const char *name, int base, int nbase)
Definition: gravi_cpl.c:275
cpl_image * gravi_image_wrap_matrix(cpl_matrix *matrix)
Wrap matrix into an image (data not duplicated)
Definition: gravi_cpl.c:2347
cpl_error_code gravi_table_new_column(cpl_table *table, const char *name, const char *unit, cpl_type type)
Definition: gravi_cpl.c:1588
double complex ** gravi_table_get_data_array_double_complex(cpl_table *table, const char *name)
Definition: gravi_cpl.c:478
cpl_image * gravi_image_from_vector(cpl_vector *vector)
Definition: gravi_cpl.c:2383
cpl_array * gravi_array_cexp(double complex factor, const cpl_array *input)
Compute the complex exponention of an array: cexp (factor * input)
Definition: gravi_cpl.c:1069
cpl_vector * gravi_table_get_vector_scalar(cpl_table *table, const char *name, cpl_size base, cpl_size nbase)
Definition: gravi_cpl.c:2223
cpl_vector * gravi_table_get_vector_diff(cpl_table *spectrum_data, int index, const char *regname1, const char *regname2)
Create a vector from the row index of the column regname.
Definition: gravi_cpl.c:2288
cpl_error_code gravi_table_compute_group_delay(cpl_table *table, const char *input, const char *flag, const char *output, cpl_table *oi_wave)
Definition: gravi_cpl.c:1467
cpl_array * gravi_array_init_float_complex(long n, float complex value)
Definition: gravi_cpl.c:562
cpl_image * gravi_image_from_column(cpl_table *table_data, const char *data_x, cpl_size row)
Create an image from a column array in table.
Definition: gravi_cpl.c:1984
cpl_error_code gravi_image_fill(cpl_image *img, double value)
Fill entire image with value.
Definition: gravi_cpl.c:2330
cpl_matrix * gravi_matrix_invertSV_create(cpl_matrix *a_in)
Invers a matrix with singular value decomposition.
Definition: gravi_cpl.c:2853
cpl_imagelist * gravi_imagelist_from_column(cpl_table *table_data, const char *data_x)
Create an imagelist from a column array in table.
Definition: gravi_cpl.c:2016
cpl_error_code gravi_array_phase_unwrap(cpl_array *input)
Definition: gravi_cpl.c:1020
cpl_error_code gravi_table_init_column_array(cpl_table *table, const char *name, const char *unit, cpl_type type, cpl_size size)
Definition: gravi_cpl.c:1624
cpl_error_code gravi_array_add_phase(cpl_array *input, double factor, cpl_array *phase)
Add a REAL phase to a REAL phase array, in-place: input = input + factor * phase.
Definition: gravi_cpl.c:1173
cpl_error_code gravi_vector_abs(cpl_vector *vector)
Return the running median of a vector, with special care for the boundaray, that are filled with the ...
Definition: gravi_cpl.c:2520
double gravi_image_get_quantile(const cpl_image *img, double value)
Compute the quantile of an image.
Definition: gravi_cpl.c:2409
cpl_error_code gravi_vector_unwrap_with_guess(cpl_vector *vector, cpl_vector *ref, double ref_to_phase)
Unwrap a phase vector following a guess vector. The difference is actually unwrap and shall thus be s...
Definition: gravi_cpl.c:2651
cpl_array * gravi_table_create_wave_array(cpl_table *oi_wave)
Definition: gravi_cpl.c:969
cpl_array * gravi_array_wrap_image(cpl_image *)
Wrap the data of na image into an array.
Definition: gravi_cpl.c:2073
cpl_error_code gravi_table_set_array_phase(cpl_table *table, const char *name, cpl_size row, cpl_array *phase)
Definition: gravi_cpl.c:670
cpl_array * gravi_array_init_double(long n, double value)
Definition: gravi_cpl.c:529
double ** gravi_table_get_data_array_double(cpl_table *table, const char *name)
Definition: gravi_cpl.c:405
cpl_error_code gravi_array_normalize_complex(cpl_array *input)
Definition: gravi_cpl.c:985
cpl_imagelist * gravi_imagelist_wrap_column(cpl_table *table_data, const char *data_x)
Wrap a column array of a table into an imagelist.
Definition: gravi_cpl.c:1689
cpl_vector * gravi_vector_median(const cpl_vector *vector, cpl_size hw)
Return the running median of a vector, with special care for the boundaray, that are filled with the ...
Definition: gravi_cpl.c:2479
cpl_error_code gravi_image_subtract_window(cpl_image *img1, const cpl_image *img2, cpl_size llx, cpl_size lly, cpl_size urx, cpl_size ury, cpl_size llx2, cpl_size lly2)
Definition: gravi_cpl.c:1903
double gravi_array_get_quantile(cpl_array *arr, double value)
Compute the value of the vector corresponding to the quantile 'thr' (0 < thr < 1)
Definition: gravi_cpl.c:2443
cpl_error_code gravi_imagelist_unwrap_images(cpl_imagelist *imglist)
Unwrap an imagelist an all its images.
Definition: gravi_cpl.c:1659
cpl_error_code gravi_array_get_group_delay_loop(cpl_array **input, cpl_array **flag, cpl_array *sigma, double *gd, cpl_size nrow, double max_width, int verbose)
Optimized computation of GDELAY for a list of arrays.
Definition: gravi_cpl.c:1335
cpl_array * gravi_table_get_column_sum_array(cpl_table *table, const char *name, int base, int nbase)
Definition: gravi_cpl.c:359
cpl_image * gravi_image_from_matrix(cpl_matrix *matrix)
Definition: gravi_cpl.c:2355
cpl_error_code gravi_array_add_phasor(cpl_array *input, double complex factor, cpl_array *phase)
Add a REAL phase to a COMPLEX array, in-place: input = input + cexp (factor * phase)
Definition: gravi_cpl.c:1117
cpl_matrix * gravi_matrix_interpolate_col(cpl_matrix *matrix, cpl_vector *xref, cpl_vector *xout)
Linear interpolation of matrix column.
Definition: gravi_cpl.c:2799
cpl_vector * gravi_table_get_vector(cpl_table *, cpl_size, const char *)
Create a vector from the row index of the column regname.
Definition: gravi_cpl.c:2167
cpl_error_code gravi_table_set_string_fixlen(cpl_table *table, const char *name, int row, const char *value, int len)
Set string in table, ensuring fixed length (right space padding). see cpl_table_set_string.
Definition: gravi_cpl.c:701
cpl_array * gravi_array_wrap_float_complex(cpl_array *input_re, cpl_array *input_im)
Definition: gravi_cpl.c:596
cpl_array * gravi_array_compute_norm2(cpl_array *input_re, cpl_array *input_im)
Definition: gravi_cpl.c:620
int gravi_array_threshold_min(cpl_array *array, double lo_cut)
Definition: gravi_cpl.c:71
cpl_error_code gravi_image_replace_window(cpl_image *img1, const cpl_image *img2, cpl_size llx, cpl_size lly, cpl_size urx, cpl_size ury, cpl_size llx2, cpl_size lly2)
Definition: gravi_cpl.c:1936
double gravi_table_get_column_std(cpl_table *table, const char *name, int base, int nbase)
Definition: gravi_cpl.c:320
cpl_array * gravi_table_get_column_mean_array(cpl_table *table, const char *name, int base, int nbase)
Definition: gravi_cpl.c:381
cpl_error_code gravi_image_subtract_collapse(cpl_image *img, const cpl_image *collapse, int direction)
Definition: gravi_cpl.c:1853
cpl_image * gravi_image_wrap_vector(cpl_vector *vector)
Wrap vector into an image (data not duplicated)
Definition: gravi_cpl.c:2375
cpl_array * gravi_array_rebin(const cpl_array *input, const cpl_array *errs, cpl_table *oi_wave_sc, cpl_table *oi_wave_ft)
Definition: gravi_cpl.c:724
cpl_array * gravi_table_create_sigma_array(cpl_table *oi_wave)
Definition: gravi_cpl.c:956
double gravi_vector_get_mean_clip(cpl_vector *vector_in, double percent, double nsigma)
Return the mean of a vector after extrema and RMS clipping.
Definition: gravi_cpl.c:2574
cpl_error_code gravi_table_multiply_scalar(cpl_table *table, const char *name, int base, int nbase, double value)
Multiply scalar or array column by scalar.
Definition: gravi_cpl.c:2698
cpl_matrix * get_matrix_from_vector(cpl_vector *, cpl_vector *)
Copy the content of two vector into a newly allocated 2D matrix.
Definition: gravi_cpl.c:2117
cpl_error_code gravi_array_phase_wrap(cpl_array *input)
Definition: gravi_cpl.c:1051
cpl_array * gravi_array_init_double_complex(long n, double complex value)
Definition: gravi_cpl.c:551
int ** gravi_table_get_data_array_int(cpl_table *table, const char *name)
Definition: gravi_cpl.c:502
cpl_array ** gravi_array_new_list(int n, cpl_type type, int size)
Allocate a list of arrays, pre-filled with 0.0.
Definition: gravi_cpl.c:93
cpl_error_code gravi_array_multiply_phasor(cpl_array *input, double complex factor, cpl_array *phase)
Multiply a REAL phase to a COMPLEX array, in-place: input = input * cexp (factor * phase)
Definition: gravi_cpl.c:1091
cpl_size gravi_vector_get_maxpos(cpl_vector *vector)
Return the index of maximum in a vector. If several indexes exists with the maximum value,...
Definition: gravi_cpl.c:2546
cpl_error_code gravi_table_add_scalar(cpl_table *table, const char *name, int base, int nbase, double value)
Multply scalar or array column by scalar.
Definition: gravi_cpl.c:2748
cpl_error_code gravi_array_add_phasors(cpl_array *input, cpl_array *add, cpl_array *sub)
Add a pair of COMPLEX arrays to a COMPLEX array, in-place: input = input + add*conj(sub)
Definition: gravi_cpl.c:1143
float ** gravi_table_get_data_array_float(cpl_table *table, const char *name)
Definition: gravi_cpl.c:429
cpl_array * gravi_array_wrap_complex(cpl_array *input_re, cpl_array *input_im)
Definition: gravi_cpl.c:573
cpl_error_code gravi_table_interpolate_column(cpl_table *to_table, const char *to_x, const char *to_y, const cpl_table *from_table, const char *from_x, const char *from_y)
Definition: gravi_cpl.c:121
cpl_error_code gravi_table_new_column_array(cpl_table *table, const char *name, const char *unit, cpl_type type, cpl_size size)
Definition: gravi_cpl.c:1610
int gravi_table_are_equal(cpl_table *first, cpl_table *second)
Check if two tables have the same content.
Definition: gravi_cpl.c:1517
cpl_error_code gravi_array_multiply_conj(cpl_array *input1, cpl_array *input2)
Definition: gravi_cpl.c:1194
cpl_image * gravi_image_collapse_median_x(cpl_image *img, cpl_size drop_from, cpl_size drop_to)
collapse image along the x direction, droping a part of the image
Definition: gravi_cpl.c:1816