---

RWindowPcaOp.hpp

#ifndef _RWindowPcaOp_H
#define _RWindowPcaOp_H

#include <stdio.h>

#include "AnImage.hpp"
#include "AbsRImageFilter.hpp"

#include "RImage.hpp"           
#include "Real3CImage.hpp"

class RWindowPcaOp : public AbsRImageFilter
{
private:

int                             width, height;

int                             pcaChannel;                     // wanted PCA channel

int                             nWindows;
int                             widthW, heightW;                // Window size
int                             nwW, nhW, dimensionW;

double*                 meanSample;
double*                 covariance;


double*                 eigenValue;
double*                 triDiagonalCov;
double*                 eigenVectors;

double*                 e;                              // temporary eigenvalue computation

protected:

void computeWindowCovariance(RImage* g);

void triDiagonalisation();

// Public householder tridiagonalisation in C
void tri_diag(double* a,double* d,double* e,double* z,int n,double tol);
int calc_eigenstructure(double* d,double* e,double* z,int n,double macheps);

public:

RWindowPcaOp(int widthWindow=8, int heightWindow=8, int pcaChannel=0);
virtual ~RWindowPcaOp();

// Default is principal component for each channel
virtual RImage* filter(RImage* source, RImage* destination=0);

// All channels recovered
// virtual Real3CImage* singleFilter(AnImage* source, Real3CImage* destination);

virtual void report(FILE* file){}
virtual void report(){}
};

#endif

Restoreinpaint sourceforge project `C++/Java Image Processing, Restoration, Inpainting Project'.

Bernard De Cuyper: Open Project Leader: Concept, design and development. Bernard De Cuyper & Eddy Fraiha 2002, 2003. Bernard De Cuyper 2004. Open and free, for friendly usage only. Modifications on Belgium ground of this piece of artistic work, by governement institutions or companies, must be notified to Bernard De Cuyper.