00001 #ifndef BP_SKY_H 00002 #define BP_SKY_H 00003 /******************************************************************************* 00004 * E.S.O. - VLT project 00005 * 00006 * "@(#) $Id: bp_sky.h,v 1.2 2005/04/29 12:18:58 amodigli Exp $" 00007 * 00008 * who when what 00009 * -------- -------- ---------------------------------------------- 00010 * schreib 06/05/03 created 00011 */ 00012 00013 /************************************************************************ 00014 * bp_sky.h 00015 * routine to search for bad pixels 00016 *---------------------------------------------------------------------- 00017 */ 00018 00019 /* 00020 * header files 00021 */ 00022 #include <cpl.h> 00023 00024 /*---------------------------------------------------------------------------- 00025 * Function ANSI C prototypes 00026 *--------------------------------------------------------------------------*/ 00027 00028 /*---------------------------------------------------------------------------- 00029 Function : bp_sky() 00030 In : ini_file: file name of according .ini file 00031 Out : integer (0 if it worked, -1 if it doesn't) 00032 Job : 00033 00034 this function searches for static bad pixels in stacks of flatfield frames, 00035 that means it stacks the flatfield in a data cube and takes the median 00036 along the z-axis to be sure to have no cosmics, then the intensity tilt of 00037 each column is removed. 00038 00039 The standard deviation and mean of the pixel values within a defined 00040 rectangular zone is determined in a way that the extreme low and high 00041 values are cut off. 00042 00043 The next step is to indicate each pixel as bad that has a deviation from 00044 the mean greater than a user defined factor times the standard deviation. 00045 This step can be leaped over if wished. 00046 00047 The next step is to calculate the median of the 8 nearest neighbors for 00048 each pixel and to determine the deviation of each pixel value from this 00049 median. 00050 00051 If this deviation is greater than a defined factor times the standard 00052 deviation the pixel value is replaced by the median. The last step is 00053 repeated to be able to consider also clusters of bad pixels. Then the 00054 resulting image is compared with the input image with the column intensity 00055 tilt removed and each changed pixel is indicated as bad. Finally, a bad 00056 pixel mask is produced that means each good pixel is marked with 1 and 00057 each bad pixel with 0. 00058 00059 ---------------------------------------------------------------------------*/ 00060 int bp_sky (cpl_parameterlist* config, cpl_frameset* set, char* procatg); 00061 #endif 00062 00063 /*--------------------------------------------------------------------------*/
1.2.13.1 written by Dimitri van Heesch,
© 1997-2001