35 #include "muse_ampl_z.h"
69 static const char *muse_ampl_help =
70 "This recipe combines several separate amplifier images (flat-fields with special FITS headers containing pico amplifier measurements) into one master image file and computes the instrumental throughput per IFU (and slice). Processing trims the raw data and records the overscan statistics, subtracts the bias (taking account of the overscan, if overscan is not &none&), and optionally, the dark from each raw input image, converts them from adu to count, scales them according to their exposure time, and combines them using input parameters. To compute the throughput the image is converted into a pixel table, the flux is then integrated over the given filter band, and the ratio of the expected flux (FITS header INS.AMPL2.CURR) to measured flux is taken, in the same units. If a geometry table was given as input, the relative area of the IFUs is taken into account when computing the flux per unit area. The resulting ratio is the instrument efficiency (throughput) and saved as QC parameters for the whole input image and per slice in the output pixel table.";
72 static const char *muse_ampl_help_esorex =
73 "\n\nInput frames for raw frame tag \"AMPL\":\n"
74 "\n Frame tag Type Req #Fr Description"
75 "\n -------------------- ---- --- --- ------------"
76 "\n AMPL raw Y >=3 Special raw flat-field taken with pico-amplifier readings in the FITS header"
77 "\n MASTER_BIAS calib Y 1 Master bias"
78 "\n MASTER_DARK calib . 1 Master dark"
79 "\n BADPIX_TABLE calib . 1 Known bad pixels"
80 "\n GEOMETRY_TABLE calib . 1 Relative positions of the slices in the field of view"
81 "\n FILTER_LIST calib Y 1 Filter definitions; here, it has to contain the filter curve for the filter given in INS.AMPL2.FILTER"
82 "\n TRACE_TABLE calib Y 1 Trace table"
83 "\n WAVECAL_TABLE calib Y 1 Wavelength calibration table"
84 "\n\nProduct frames for raw frame tag \"AMPL\":\n"
85 "\n Frame tag Level Description"
86 "\n -------------------- -------- ------------"
87 "\n MASTER_AMPL final Combined master AMPL image, written if --savemaster=true"
88 "\n TABLE_AMPL final Output table with computations for each CCD pixel, written if --savetable=true."
89 "\n AMPL_CONVOLVED final Combined and convolved master AMPL image";
100 static cpl_recipeconfig *
101 muse_ampl_new_recipeconfig(
void)
103 cpl_recipeconfig *recipeconfig = cpl_recipeconfig_new();
107 cpl_recipeconfig_set_tag(recipeconfig, tag, 3, -1);
108 cpl_recipeconfig_set_input(recipeconfig, tag,
"MASTER_BIAS", 1, 1);
109 cpl_recipeconfig_set_input(recipeconfig, tag,
"MASTER_DARK", -1, 1);
110 cpl_recipeconfig_set_input(recipeconfig, tag,
"BADPIX_TABLE", -1, 1);
111 cpl_recipeconfig_set_input(recipeconfig, tag,
"GEOMETRY_TABLE", -1, 1);
112 cpl_recipeconfig_set_input(recipeconfig, tag,
"FILTER_LIST", 1, 1);
113 cpl_recipeconfig_set_input(recipeconfig, tag,
"TRACE_TABLE", 1, 1);
114 cpl_recipeconfig_set_input(recipeconfig, tag,
"WAVECAL_TABLE", 1, 1);
115 cpl_recipeconfig_set_output(recipeconfig, tag,
"MASTER_AMPL");
116 cpl_recipeconfig_set_output(recipeconfig, tag,
"TABLE_AMPL");
117 cpl_recipeconfig_set_output(recipeconfig, tag,
"AMPL_CONVOLVED");
133 static cpl_error_code
134 muse_ampl_prepare_header(
const char *aFrametag, cpl_propertylist *aHeader)
136 cpl_ensure_code(aFrametag, CPL_ERROR_NULL_INPUT);
137 cpl_ensure_code(aHeader, CPL_ERROR_NULL_INPUT);
138 if (!strcmp(aFrametag,
"MASTER_AMPL")) {
141 "Number of saturated pixels in raw image i in input list");
144 "Number of saturated pixels in output master image");
145 }
else if (!strcmp(aFrametag,
"TABLE_AMPL")) {
148 "[ph] Integrated number of photons detected on the CCD");
151 "[W] Integrated power detected on the CCD");
154 "[%] Channel throughput compared to photodiode a");
157 "[%] Estimated error of throughput compared to photodiode 2");
160 "[%] Slice j throughput compared to photodiode 2");
161 }
else if (!strcmp(aFrametag,
"AMPL_CONVOLVED")) {
164 "Number of saturated pixels in raw image i in input list");
167 "Number of saturated pixels in output master image");
170 "[ph] Integrated number of photons detected on the CCD");
173 "[W] Integrated power detected on the CCD");
176 "[%] Channel throughput compared to photodiode a");
179 "[%] Estimated error of throughput compared to photodiode 2");
182 "[%] Slice j throughput compared to photodiode 2");
184 cpl_msg_warning(__func__,
"Frame tag %s is not defined", aFrametag);
185 return CPL_ERROR_ILLEGAL_INPUT;
187 return CPL_ERROR_NONE;
200 static cpl_frame_level
201 muse_ampl_get_frame_level(
const char *aFrametag)
204 return CPL_FRAME_LEVEL_NONE;
206 if (!strcmp(aFrametag,
"MASTER_AMPL")) {
207 return CPL_FRAME_LEVEL_FINAL;
209 if (!strcmp(aFrametag,
"TABLE_AMPL")) {
210 return CPL_FRAME_LEVEL_FINAL;
212 if (!strcmp(aFrametag,
"AMPL_CONVOLVED")) {
213 return CPL_FRAME_LEVEL_FINAL;
215 return CPL_FRAME_LEVEL_NONE;
229 muse_ampl_get_frame_mode(
const char *aFrametag)
234 if (!strcmp(aFrametag,
"MASTER_AMPL")) {
237 if (!strcmp(aFrametag,
"TABLE_AMPL")) {
240 if (!strcmp(aFrametag,
"AMPL_CONVOLVED")) {
258 muse_ampl_create(cpl_plugin *aPlugin)
262 if (cpl_plugin_get_type(aPlugin) == CPL_PLUGIN_TYPE_RECIPE) {
263 recipe = (cpl_recipe *)aPlugin;
271 muse_ampl_new_recipeconfig(),
272 muse_ampl_prepare_header,
273 muse_ampl_get_frame_level,
274 muse_ampl_get_frame_mode);
279 cpl_msg_set_time_on();
283 recipe->parameters = cpl_parameterlist_new();
288 p = cpl_parameter_new_range(
"muse.muse_ampl.nifu",
290 "IFU to handle. If set to 0, all IFUs are processed serially. If set to -1, all IFUs are processed in parallel.",
295 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"nifu");
296 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"nifu");
298 cpl_parameterlist_append(recipe->parameters, p);
301 p = cpl_parameter_new_value(
"muse.muse_ampl.overscan",
303 "If this is \"none\", stop when detecting discrepant overscan levels (see ovscsigma), for \"offset\" it assumes that the mean overscan level represents the real offset in the bias levels of the exposures involved, and adjusts the data accordingly; for \"vpoly\", a polynomial is fit to the vertical overscan and subtracted from the whole quadrant.",
305 (
const char *)
"vpoly");
306 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"overscan");
307 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"overscan");
309 cpl_parameterlist_append(recipe->parameters, p);
312 p = cpl_parameter_new_value(
"muse.muse_ampl.ovscreject",
314 "This influences how values are rejected when computing overscan statistics. Either no rejection at all (\"none\"), rejection using the DCR algorithm (\"dcr\"), or rejection using an iterative constant fit (\"fit\").",
316 (
const char *)
"dcr");
317 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"ovscreject");
318 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"ovscreject");
320 cpl_parameterlist_append(recipe->parameters, p);
323 p = cpl_parameter_new_value(
"muse.muse_ampl.ovscsigma",
325 "If the deviation of mean overscan levels between a raw input image and the reference image is higher than |ovscsigma x stdev|, stop the processing. If overscan=\"vpoly\", this is used as sigma rejection level for the iterative polynomial fit (the level comparison is then done afterwards with |100 x stdev| to guard against incompatible settings). Has no effect for overscan=\"offset\".",
328 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"ovscsigma");
329 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"ovscsigma");
331 cpl_parameterlist_append(recipe->parameters, p);
334 p = cpl_parameter_new_value(
"muse.muse_ampl.ovscignore",
336 "The number of pixels of the overscan adjacent to the data region of the CCD that are ignored when computing statistics or fits.",
339 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"ovscignore");
340 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"ovscignore");
342 cpl_parameterlist_append(recipe->parameters, p);
345 p = cpl_parameter_new_enum(
"muse.muse_ampl.combine",
347 "Type of combination to use",
349 (
const char *)
"sigclip",
351 (
const char *)
"average",
352 (
const char *)
"median",
353 (
const char *)
"minmax",
354 (
const char *)
"sigclip");
355 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"combine");
356 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"combine");
358 cpl_parameterlist_append(recipe->parameters, p);
361 p = cpl_parameter_new_value(
"muse.muse_ampl.nlow",
363 "Number of minimum pixels to reject with minmax",
366 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"nlow");
367 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"nlow");
369 cpl_parameterlist_append(recipe->parameters, p);
372 p = cpl_parameter_new_value(
"muse.muse_ampl.nhigh",
374 "Number of maximum pixels to reject with minmax",
377 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"nhigh");
378 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"nhigh");
380 cpl_parameterlist_append(recipe->parameters, p);
383 p = cpl_parameter_new_value(
"muse.muse_ampl.nkeep",
385 "Number of pixels to keep with minmax",
388 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"nkeep");
389 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"nkeep");
391 cpl_parameterlist_append(recipe->parameters, p);
394 p = cpl_parameter_new_value(
"muse.muse_ampl.lsigma",
396 "Low sigma for pixel rejection with sigclip",
399 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"lsigma");
400 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"lsigma");
402 cpl_parameterlist_append(recipe->parameters, p);
405 p = cpl_parameter_new_value(
"muse.muse_ampl.hsigma",
407 "High sigma for pixel rejection with sigclip",
410 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"hsigma");
411 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"hsigma");
413 cpl_parameterlist_append(recipe->parameters, p);
416 p = cpl_parameter_new_value(
"muse.muse_ampl.fbeam",
418 "Factor to describe the widening of the beam from the focal plane to photo diode 2.",
421 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"fbeam");
422 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"fbeam");
424 cpl_parameterlist_append(recipe->parameters, p);
427 p = cpl_parameter_new_value(
"muse.muse_ampl.temp",
429 "Lamp temperature [K] used to create the black body function.",
432 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"temp");
433 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"temp");
435 cpl_parameterlist_append(recipe->parameters, p);
438 p = cpl_parameter_new_value(
"muse.muse_ampl.savemaster",
440 "Save the processed and combined master image before any concolution is done.",
443 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"savemaster");
444 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"savemaster");
446 cpl_parameterlist_append(recipe->parameters, p);
449 p = cpl_parameter_new_value(
"muse.muse_ampl.savetable",
451 "Save the table with all the processed pixel values.",
454 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"savetable");
455 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"savetable");
457 cpl_parameterlist_append(recipe->parameters, p);
460 p = cpl_parameter_new_value(
"muse.muse_ampl.merge",
462 "Merge output products from different IFUs into a common file.",
465 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CFG,
"merge");
466 cpl_parameter_set_alias(p, CPL_PARAMETER_MODE_CLI,
"merge");
468 cpl_parameterlist_append(recipe->parameters, p);
488 cpl_ensure_code(aParams, CPL_ERROR_NULL_INPUT);
489 cpl_ensure_code(aParameters, CPL_ERROR_NULL_INPUT);
492 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.nifu");
493 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
494 aParams->
nifu = cpl_parameter_get_int(p);
496 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.overscan");
497 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
498 aParams->
overscan = cpl_parameter_get_string(p);
500 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.ovscreject");
501 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
502 aParams->
ovscreject = cpl_parameter_get_string(p);
504 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.ovscsigma");
505 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
506 aParams->
ovscsigma = cpl_parameter_get_double(p);
508 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.ovscignore");
509 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
510 aParams->
ovscignore = cpl_parameter_get_int(p);
512 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.combine");
513 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
514 aParams->
combine_s = cpl_parameter_get_string(p);
516 (!strcasecmp(aParams->
combine_s,
"average")) ? MUSE_AMPL_PARAM_COMBINE_AVERAGE :
517 (!strcasecmp(aParams->
combine_s,
"median")) ? MUSE_AMPL_PARAM_COMBINE_MEDIAN :
518 (!strcasecmp(aParams->
combine_s,
"minmax")) ? MUSE_AMPL_PARAM_COMBINE_MINMAX :
519 (!strcasecmp(aParams->
combine_s,
"sigclip")) ? MUSE_AMPL_PARAM_COMBINE_SIGCLIP :
520 MUSE_AMPL_PARAM_COMBINE_INVALID_VALUE;
521 cpl_ensure_code(aParams->
combine != MUSE_AMPL_PARAM_COMBINE_INVALID_VALUE,
522 CPL_ERROR_ILLEGAL_INPUT);
524 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.nlow");
525 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
526 aParams->
nlow = cpl_parameter_get_int(p);
528 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.nhigh");
529 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
530 aParams->
nhigh = cpl_parameter_get_int(p);
532 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.nkeep");
533 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
534 aParams->
nkeep = cpl_parameter_get_int(p);
536 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.lsigma");
537 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
538 aParams->
lsigma = cpl_parameter_get_double(p);
540 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.hsigma");
541 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
542 aParams->
hsigma = cpl_parameter_get_double(p);
544 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.fbeam");
545 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
546 aParams->
fbeam = cpl_parameter_get_double(p);
548 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.temp");
549 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
550 aParams->
temp = cpl_parameter_get_double(p);
552 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.savemaster");
553 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
554 aParams->
savemaster = cpl_parameter_get_bool(p);
556 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.savetable");
557 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
558 aParams->
savetable = cpl_parameter_get_bool(p);
560 p = cpl_parameterlist_find(aParameters,
"muse.muse_ampl.merge");
561 cpl_ensure_code(p, CPL_ERROR_DATA_NOT_FOUND);
562 aParams->
merge = cpl_parameter_get_bool(p);
576 muse_ampl_exec(cpl_plugin *aPlugin)
578 if (cpl_plugin_get_type(aPlugin) != CPL_PLUGIN_TYPE_RECIPE) {
581 cpl_recipe *recipe = (cpl_recipe *)aPlugin;
582 cpl_msg_set_threadid_on();
584 cpl_frameset *usedframes = cpl_frameset_new(),
585 *outframes = cpl_frameset_new();
587 muse_ampl_params_fill(¶ms, recipe->parameters);
589 cpl_errorstate prestate = cpl_errorstate_get();
591 if (params.
nifu < -1 || params.
nifu > kMuseNumIFUs) {
592 cpl_msg_error(__func__,
"Please specify a valid IFU number (between 1 and "
593 "%d), 0 (to process all IFUs consecutively), or -1 (to "
594 "process all IFUs in parallel) using --nifu.", kMuseNumIFUs);
598 cpl_boolean donotmerge = CPL_FALSE;
600 if (params.
nifu > 0) {
603 rc = muse_ampl_compute(proc, ¶ms);
604 cpl_frameset_join(usedframes, proc->
usedframes);
605 cpl_frameset_join(outframes, proc->
outframes);
607 donotmerge = CPL_TRUE;
608 }
else if (params.
nifu < 0) {
609 int *rcs = cpl_calloc(kMuseNumIFUs,
sizeof(
int));
611 #pragma omp parallel for default(none) \
612 shared(outframes, params, rcs, recipe, usedframes)
613 for (nifu = 1; nifu <= kMuseNumIFUs; nifu++) {
619 int *rci = rcs + (nifu - 1);
620 *rci = muse_ampl_compute(proc, pars);
621 if (rci && cpl_error_get_code() == MUSE_ERROR_CHIP_NOT_LIVE) {
625 #pragma omp critical(muse_processing_used_frames)
626 cpl_frameset_join(usedframes, proc->
usedframes);
627 #pragma omp critical(muse_processing_output_frames)
628 cpl_frameset_join(outframes, proc->
outframes);
633 for (nifu = 1; nifu <= kMuseNumIFUs; nifu++) {
634 if (rcs[nifu-1] != 0) {
640 for (params.
nifu = 1; params.
nifu <= kMuseNumIFUs && !rc; params.
nifu++) {
643 rc = muse_ampl_compute(proc, ¶ms);
644 if (rc && cpl_error_get_code() == MUSE_ERROR_CHIP_NOT_LIVE) {
647 cpl_frameset_join(usedframes, proc->
usedframes);
648 cpl_frameset_join(outframes, proc->
outframes);
652 UNUSED_ARGUMENT(donotmerge);
654 if (!cpl_errorstate_is_equal(prestate)) {
658 cpl_msg_set_level(CPL_MSG_INFO);
665 if (params.
merge && !donotmerge) {
666 muse_utils_frameset_merge_frames(outframes);
674 cpl_frameset_join(recipe->frames, usedframes);
675 cpl_frameset_join(recipe->frames, outframes);
676 cpl_frameset_delete(usedframes);
677 cpl_frameset_delete(outframes);
690 muse_ampl_destroy(cpl_plugin *aPlugin)
694 if (cpl_plugin_get_type(aPlugin) == CPL_PLUGIN_TYPE_RECIPE) {
695 recipe = (cpl_recipe *)aPlugin;
701 cpl_parameterlist_delete(recipe->parameters);
718 cpl_plugin_get_info(cpl_pluginlist *aList)
720 cpl_recipe *recipe = cpl_calloc(1,
sizeof *recipe);
721 cpl_plugin *plugin = &recipe->interface;
725 helptext = cpl_sprintf(
"%s%s", muse_ampl_help,
726 muse_ampl_help_esorex);
728 helptext = cpl_sprintf(
"%s", muse_ampl_help);
732 cpl_plugin_init(plugin, CPL_PLUGIN_API, MUSE_BINARY_VERSION,
733 CPL_PLUGIN_TYPE_RECIPE,
735 "Determine the instrumental throughput from exposures taken with the pico-amplifier / photo diode readings.",
743 cpl_pluginlist_append(aList, plugin);
void muse_processing_delete(muse_processing *aProcessing)
Free the muse_processing structure.
int nifu
IFU to handle. If set to 0, all IFUs are processed serially. If set to -1, all IFUs are processed in ...
Structure to hold the parameters of the muse_ampl recipe.
muse_cplframework_type muse_cplframework(void)
Return the CPL framework the recipe is run under.
int nhigh
Number of maximum pixels to reject with minmax.
double ovscsigma
If the deviation of mean overscan levels between a raw input image and the reference image is higher ...
double temp
Lamp temperature [K] used to create the black body function.
cpl_frameset * usedframes
muse_processing * muse_processing_new(const char *aName, cpl_recipe *aRecipe)
Create a new processing structure.
const char * muse_get_license(void)
Get the pipeline copyright and license.
int nlow
Number of minimum pixels to reject with minmax.
int combine
Type of combination to use.
int ovscignore
The number of pixels of the overscan adjacent to the data region of the CCD that are ignored when com...
const char * overscan
If this is "none", stop when detecting discrepant overscan levels (see ovscsigma), for "offset" it assumes that the mean overscan level represents the real offset in the bias levels of the exposures involved, and adjusts the data accordingly; for "vpoly", a polynomial is fit to the vertical overscan and subtracted from the whole quadrant.
int merge
Merge output products from different IFUs into a common file.
int savetable
Save the table with all the processed pixel values.
int nkeep
Number of pixels to keep with minmax.
void muse_cplerrorstate_dump_some(unsigned aCurrent, unsigned aFirst, unsigned aLast)
Dump some CPL errors.
void muse_processinginfo_delete(cpl_recipe *)
Clear all information from the processing info and from the recipe config.
double hsigma
High sigma for pixel rejection with sigclip.
const char * ovscreject
This influences how values are rejected when computing overscan statistics. Either no rejection at al...
double fbeam
Factor to describe the widening of the beam from the focal plane to photo diode 2.
const char * combine_s
Type of combination to use (as string)
cpl_error_code muse_cplframeset_erase_duplicate(cpl_frameset *aFrames)
Erase all duplicate frames from a frameset.
cpl_error_code muse_cplframeset_erase_all(cpl_frameset *aFrames)
Erase all frames in a frameset.
int savemaster
Save the processed and combined master image before any concolution is done.
void muse_processinginfo_register(cpl_recipe *, cpl_recipeconfig *, muse_processing_prepare_header_func *, muse_processing_get_frame_level_func *, muse_processing_get_frame_mode_func *)
Register extended functionalities for MUSE recipes.
double lsigma
Low sigma for pixel rejection with sigclip.
cpl_error_code muse_processing_prepare_property(cpl_propertylist *, const char *, cpl_type, const char *)
Prepare and check the specified property.