36 #include "detmon_utils.h"
37 #include "irplib_utils.h"
56 const char * detmon_get_license(
void)
58 const char * detmon_license =
59 "This file is part of the DETMON Instrument Pipeline\n"
60 "Copyright (C) 2002,2003 European Southern Observatory\n"
62 "This program is free software; you can redistribute it and/or modify\n"
63 "it under the terms of the GNU General Public License as published by\n"
64 "the Free Software Foundation; either version 2 of the License, or\n"
65 "(at your option) any later version.\n"
67 "This program is distributed in the hope that it will be useful,\n"
68 "but WITHOUT ANY WARRANTY; without even the implied warranty of\n"
69 "MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the\n"
70 "GNU General Public License for more details.\n"
72 "You should have received a copy of the GNU General Public License\n"
73 "along with this program; if not, write to the Free Software\n"
74 "Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, \n"
76 return detmon_license ;
92 static void expand_region(
const cpl_image * ref,
103 *urx = cpl_image_get_size_x(ref);
105 *ury = cpl_image_get_size_y(ref);
131 detmon_load_frameset_window(
const cpl_frameset * fset, cpl_type type,
132 cpl_size pnum, cpl_size xtnum,
133 cpl_size llx, cpl_size lly,
134 cpl_size urx, cpl_size ury,
135 cpl_size nx, cpl_size ny)
137 cpl_imagelist * list = cpl_imagelist_new();
138 const cpl_size n = cpl_frameset_get_size(fset);
140 if (nx >= 0 && ny >= 0) {
141 error_if(urx - llx + 1 > nx || ury - lly + 1 > ny ,
142 CPL_ERROR_ILLEGAL_INPUT,
143 "window size [%d:%d,%d:%d] larger than output size [%d, %d]",
144 (
int)llx, (
int)urx, (
int)lly, (
int)ury, (
int)nx, (
int)ny);
147 for (cpl_size i = 0; i < n; i++) {
148 const cpl_frame * frm = cpl_frameset_get_position_const(fset, i);
149 const char * fn = cpl_frame_get_filename(frm);
151 const cpl_size next = cpl_frame_get_nextensions(frm);
152 for (cpl_size e = 0; e < next; e++) {
153 cpl_image * image = cpl_image_load_window(fn, type, pnum, e,
155 skip_if(image == NULL);
156 if(nx < 0 || ny < 0 ) {
157 cpl_imagelist_set(list, image, cpl_imagelist_get_size(list));
159 cpl_image * full = cpl_image_new(nx, ny, type);
160 cpl_image_copy(full, image, llx, lly);
161 cpl_image_delete(image);
162 cpl_imagelist_set(list, full, cpl_imagelist_get_size(list));
168 cpl_image * image = cpl_image_load_window(fn, type, pnum, xtnum,
170 skip_if(image == NULL);
171 if(nx < 0 || ny < 0 ) {
172 cpl_imagelist_set(list, image, cpl_imagelist_get_size(list));
174 cpl_image * full = cpl_image_new(nx, ny, type);
175 cpl_image_copy(full, image, llx, lly);
176 cpl_image_delete(image);
177 cpl_imagelist_set(list, full, cpl_imagelist_get_size(list));
184 if (cpl_error_get_code()) {
185 cpl_imagelist_delete(list);
207 cpl_image * detmon_subtract_create_window(
const cpl_image * a,
214 cpl_image * na = cpl_image_extract(a, llx, lly, urx, ury);
215 cpl_image * nb = cpl_image_extract(b, llx, lly, urx, ury);
216 cpl_image_subtract(na, nb);
217 cpl_image_delete(nb);
240 cpl_image * detmon_subtracted_avg(
const cpl_image * on1,
241 const cpl_image * off1,
242 const cpl_image * on2,
243 const cpl_image * off2,
244 cpl_size llx, cpl_size lly,
245 cpl_size urx, cpl_size ury)
247 expand_region(on1, &llx, &lly, &urx, &ury);
248 cpl_image * dif1 = cpl_image_extract(on1, llx, lly, urx, ury);
249 cpl_image * dif2 = cpl_image_extract(on2, llx, lly, urx, ury);
250 cpl_image * b = cpl_image_extract(off1, llx, lly, urx, ury);
253 cpl_image_subtract(dif1, b);
256 b = cpl_image_extract(off2, llx, lly, urx, ury);
257 cpl_image_subtract(dif1, b);
260 cpl_image_subtract(dif2, b);
264 dif_avg = cpl_image_average_create(dif1, dif2);
265 cpl_image_abs(dif_avg);
267 cpl_image_delete(dif1);
268 cpl_image_delete(dif2);