/home/kopmann/git-mirror/fdhwlib/fdhwlib/akutil/analysematrix.h

Go to the documentation of this file.

/***************************************************************************
    analyseMatrix.h  -  description

    begin                : Tue Apr 18 2006
    copyright            : (C) 2006 by Andreas Kopmann
    email                : kopmann@hpe.fzk.de
 ***************************************************************************/


#ifndef ANALYSEMATRIX_H
#define ANALYSEMATRIX_H

#include <cstdio>


class analyseMatrix {
public: 
  analyseMatrix();
  ~analyseMatrix();

  void filterBy2x2(int filter[2][2], int col, int row, 
           double matrix[20][22], double res[20][22], int mask[20][22]);

  void displayMatrix(FILE *fout, double matrix[20][22], const char *title);

  void displayMatrix(FILE *fout, int matrix[20][22], const char *title, int norm = 0, 
                      int value = 0, const char *tag = 0);

  int displayPixelList(FILE *fout, int matrix[20][22]);

  void classifyByThreshold(double thresh, double matrix[20][22], int res[20][22], int up = 1);

  void setMatrix(double value, double matrix[20][22]);
  void setMatrix(int value, int matrix[20][22]);

  void multMatrix(double value, double matrix[20][22]);  

  void multMask(int value, int mask[20][22]);
  
  void divMatrix(double res[20][22], double matrix[20][22]);
  
  void divMatrix(double res[20][22], int matrix[20][22]);
  
  void meanMatrix(double matrix[20][22], 
              double *mean, double *var, int *n, int mask[20][22]);

  void meanMatrix(double matrix[20][22], 
              double *mean, double *var, int *n, int mask[20][22], int err[20][22]);

  void addMask(int res[20][22], int mask[20][22]);

  void addMatrix(double res[20][22], double matrix[20][22]);
  
  void copyMatrix(double res[20][22], double matrix[20][22]);
  
    
  void orMask(int res[20][22], int mask[20][22]);

  void invertMask(int mask[20][22]); 

  int nMask(int mask[20][22]);

  void copyMask(int res[20][22], int mask[20][22]);

  void maxMask(int max, int mask[20][22]);

  void minMask(int min, int mask[20][22]);

  
  void findHighNoise1(FILE *fout, double matrix[20][22], int res[20][22], int err[20][22]);

  void findHighNoise2(FILE *fout, double matrix[20][22], int res[20][22], int err[20][22]);
  
  // ======= Function working with the local arrays ============
        
  void clearNoisyPixel(); 
  
  void clearFaultyPixel();

  void  addNoisyPixel(int mask[20][22]);
  
  void findHighNoise(FILE *fout, double matrix[20][22]);

  void meanOfNormalPixel(double matrix[20][22], 
              double *mean, double *var, int *n);
      
  void meanOfNoisyPixel(double matrix[20][22], 
              double *mean, double *var, int *n);

  int nNoisyPixel();      
  
  void displayNoisyPixelList(FILE *fout);       

                  
  // ========== Local arrays =======
      
  int noisyPixel[20][22];

  int noisyPixelLast[20][22];
  
  int faultyPixel[20][22];

  int noisyRuns;


  double noisyNVariances;

  int noisyMinVariance;

  double noisyNeighborDiff;
 
  double noisyAbsThresh; 

  double noisyLowThresh;  

  
  int debug;
  
};

#endif