C:/Users/albertlr/projects/TarzaNN/TarzaNN/Filter.h

Go to the documentation of this file.

#ifndef FILTER_H_INCLUDED
#define FILTER_H_INCLUDED

#include "common.h"
#include "Matrix.h"
#include <math.h>
#include "nowarn.h"
#include <QMutex>
#include <QLinkedList>
#include <QTextStream>

#define EVEN_FILTERS

class GatingUnits;

/*how to extend images for convolution*/
enum
{
    convolutionMirror           =       2,      /* Mirror border pixels, width determined by filter size*/
    convolutionFill                     =       4,      /* Use the background color for missing pixels (?)*/
    convolutionEdgeExtend       =       8,      /* Repeat edge pixels*/
};

class FilterOperation {
public:
        FilterOperation(QString theName,QString theValue) {
                name = theName;
                value = theValue;
        };
        FilterOperation() {};
        QString name;
        QString value;
};

class Filter: public Matrix {
        friend class UI_FilterProp;
        friend class UI_DesignerPropDisplay;
public:
    Filter(int w, int h);
        Filter(QString fileName, int fileType); 
        float* getKernel();
    void _convolve( Matrix* input, Matrix* output);
    void convolve( Matrix* input, Matrix* output, bool scale=true);
    void STMconvolve( Matrix* input, GatingUnits* gu);
    void setFromConvolveWindow(int _x1,int _y1, int _w, int _h);

    static void _convolveRetinal( Matrix* input, Matrix* output, int padding, int distPixels);
    static void destroyRetinalImageFilters();
    static Filter* retinalImageFilters[1024];
    static int retinalImageFilterNo;
    static int retinalImageFilterSize;
    static QMutex retinalMutex;
    
    void setPaddingType(int _type);
    int getPaddingType();
    void separateKernel();
    static void svd(double **a, int m, int n, double w[], double **v);
    static double pythag(double a,double b) ;
    ~Filter();

        void addOperation(const QString &opName, const QString &opValue);

        void setParameters(paramMap params);
        void XMLSerialize(int tabs, QTextStream *buf);
        static int getType(const QString & name);
        static QString getType(int value);
        
#ifdef ST_CONVOLVE
        static double compute_skip( int input, int output, int filter){
            return (double)(input-filter)/(double)(output-1);
        }
#else
        static int computeSkip(int input, int output){
                return (int)(input/output);
        }
#endif//ST_CONVOLVE

    // should be protected, Gfilter friend, etc
    float l1_norm, l2_norm;
    int env_x, env_y;
    int pos_x, pos_y;
    void adjust_l(float l1, float l2);
        QImage* getPreviewImage();

    // filter type
    int fType;
        
protected:
    int separable_status;
    enum {SEPARABLE_NOT_CHECKED, SEPARABLE_NO, SEPARABLE_YES};
    Filter **separated_filters; //for now only2

    void find_limits( int w, int h, float sigma, float aspect, float orientation, float& h_m_x, float& h_m_y);
    Filter* nextFilter;
        
        paramMap initParams;
        QLinkedList<FilterOperation> operations;

private:
    Filter();
        float* kernel;
    bool useConvolveWindow;
    int convolve_x1, convolve_y1;
    int convolve_w, convolve_h;

    int padding_type;

#ifdef VIMAGE_H
    void* tempConvBuffer;
        void* tempScaleBuffer;
#endif  
};

#endif //FILTER_H_INCLUDED