00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029 #ifdef HAVE_CONFIG_H
00030 # include <config.h>
00031 #endif
00032
00033
00034
00035
00036
00037 #include "sinfo_badsky_ini_by_cpl.h"
00038 #include "sinfo_raw_types.h"
00039 #include "sinfo_pro_types.h"
00040 #include "sinfo_hidden.h"
00041 #include "sinfo_globals.h"
00042 #include "sinfo_functions.h"
00043 #include "sinfo_file_handling.h"
00044
00045
00046
00047
00048 static void
00049 parse_section_frames(badsky_config * cfg, cpl_parameterlist* cpl_cfg,
00050 cpl_frameset* sof, cpl_frameset** raw, int* status);
00051 static void
00052 parse_section_badpix(badsky_config * cfg, cpl_parameterlist* cpl_cfg);
00053 static void
00054 parse_section_thresh(badsky_config * cfg, cpl_parameterlist* cpl_cfg);
00060
00061
00062
00074 badsky_config *
00075 sinfo_parse_cpl_input_badsky(cpl_parameterlist * cpl_cfg,
00076 cpl_frameset* sof,
00077 cpl_frameset** raw)
00078 {
00079 badsky_config * cfg ;
00080 int status = 0;
00081
00082 cfg = sinfo_badsky_cfg_create();
00083 parse_section_badpix (cfg, cpl_cfg);
00084 parse_section_thresh (cfg, cpl_cfg);
00085 parse_section_frames (cfg, cpl_cfg, sof, raw,&status);
00086 if(status>0) {
00087 sinfo_msg_error("parsing cpl input");
00088 sinfo_badsky_cfg_destroy(cfg);
00089 cfg = NULL ;
00090 return NULL ;
00091 }
00092 return cfg ;
00093
00094 }
00095
00107 static void
00108 parse_section_frames(badsky_config * cfg,
00109 cpl_parameterlist * cpl_cfg,
00110 cpl_frameset * sof,
00111 cpl_frameset ** raw,
00112 int* status)
00113 {
00114 int i=0;
00115 int nraw = 0;
00116 char * tag=NULL;
00117
00118 int nraw_good = 0;
00119 cpl_frame* frame=NULL;
00120 cpl_parameter *p;
00121 char spat_res[FILE_NAME_SZ];
00122 char lamp_status[FILE_NAME_SZ];
00123 char band[FILE_NAME_SZ];
00124 int ins_set=0;
00125
00126 *raw=cpl_frameset_new();
00127
00128
00129 sinfo_contains_frames_type(sof,raw,RAW_SKY);
00130 nraw=cpl_frameset_get_size(*raw);
00131 if (nraw < 1) {
00132 sinfo_msg_error("Too few (%d) raw frames %s present in"
00133 "frameset!Aborting...",nraw,RAW_SKY);
00134 (*status)++;
00135 return;
00136 }
00137
00138
00139
00140 cfg->framelist = cpl_malloc(nraw * sizeof(char*));
00141
00142
00143 for (i=0 ; i<nraw ; i++) {
00144 frame = cpl_frameset_get_frame(*raw,i);
00145 tag = (char*)cpl_frame_get_tag(frame) ;
00146 if(sinfo_file_exists((char*) cpl_frame_get_filename(frame))==1)
00147 {
00148
00149 cfg->framelist[i]=cpl_strdup(cpl_frame_get_filename(frame));
00150 nraw_good++;
00151 }
00152 }
00153
00154 if(nraw_good<1) {
00155 sinfo_msg_error("Error: no good raw frame in input, something wrong!");
00156 (*status)++;
00157 return;
00158 }
00159
00160
00161 if(NULL != cpl_frameset_find(sof,PRO_MASTER_DARK)) {
00162 frame = cpl_frameset_find(sof,PRO_MASTER_DARK);
00163 strcpy(cfg -> sinfo_dark, cpl_strdup(cpl_frame_get_filename(frame)));
00164 } else {
00165 sinfo_msg_warning("Frame %s not found!", PRO_MASTER_DARK);
00166 }
00167
00168
00169 cfg->nframes = nraw_good ;
00170
00171
00172 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.out_filename");
00173 strcpy(cfg -> outName, cpl_parameter_get_string(p));
00174
00175 frame = cpl_frameset_get_frame(*raw,0);
00176 sinfo_get_spatial_res(frame,spat_res);
00177
00178
00179 switch(sinfo_frame_is_on(frame))
00180 {
00181 case 0:
00182 strcpy(lamp_status,"on");
00183 break;
00184 case 1:
00185 strcpy(lamp_status,"off");
00186 break;
00187 case -1:
00188 strcpy(lamp_status,"undefined");
00189 break;
00190 default:
00191 strcpy(lamp_status,"undefined");
00192 break;
00193 }
00194
00195 sinfo_get_band(frame,band);
00196 sinfo_msg("spatial resolution: %s lamp status: %s band: %s \n",
00197 spat_res, lamp_status, band);
00198
00199
00200 sinfo_get_ins_set(band,&ins_set);
00201 return;
00202
00203 }
00204
00213 static void
00214 parse_section_badpix(badsky_config * cfg, cpl_parameterlist * cpl_cfg)
00215 {
00216 cpl_parameter *p;
00217
00218 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.sigma_factor");
00219 cfg -> sigmaFactor = cpl_parameter_get_double(p);
00220
00221 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.method_index");
00222 cfg -> methodInd = cpl_parameter_get_int(p);
00223
00224 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.factor");
00225 cfg -> factor = cpl_parameter_get_double(p);
00226
00227 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.iterations");
00228 cfg -> iterations = cpl_parameter_get_int(p);
00229
00230 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.low_rejection");
00231 cfg -> loReject = cpl_parameter_get_double(p);
00232
00233 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.high_rejection");
00234 cfg -> hiReject = cpl_parameter_get_double(p);
00235
00236 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.llx");
00237 cfg -> llx = cpl_parameter_get_int(p);
00238
00239 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.lly");
00240 cfg -> lly = cpl_parameter_get_int(p);
00241
00242 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.urx");
00243 cfg -> urx = cpl_parameter_get_int(p);
00244
00245 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.ury");
00246 cfg -> ury = cpl_parameter_get_int(p);
00247
00248 }
00249
00258 static void
00259 parse_section_thresh(badsky_config * cfg, cpl_parameterlist * cpl_cfg)
00260 {
00261 cpl_parameter *p;
00262
00263 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.threshold_index");
00264 cfg -> threshInd = cpl_parameter_get_bool(p);
00265
00266 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.mean_factor");
00267 cfg -> meanfactor = cpl_parameter_get_double(p);
00268
00269
00270 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.min_cut");
00271 cfg -> mincut = cpl_parameter_get_double(p);
00272
00273 p = cpl_parameterlist_find(cpl_cfg, "sinfoni.bp_sky.max_cut");
00274 cfg -> maxcut = cpl_parameter_get_double(p);
00275
00276 }
00277
00285 void
00286 sinfo_badsky_free(badsky_config * cfg)
00287 {
00288 cpl_free(cfg->framelist);
00289 sinfo_badsky_cfg_destroy(cfg);
00290
00291 return;
00292
00293 }