---

MDistanceFMM.hpp

#ifndef __MDistanceFMM__
#define __MDistanceFMM__



#include "AnImage.hpp"
#include "RImage.hpp"

#include "AnImageOp.hpp"

#include "MPsi2DValue.hpp"
#include "MDomainOfAction.hpp"
#include "MCandidateHeapList.hpp"
#include "MNeigbourSet.hpp"


class MDistanceFMM : public AnImageOp
{
public:

const static int        lookupTable[256][3];

protected:

int                                     frontMax;
double*                         speedFactor;

int                                     maskFlag;       // 1= mask area, -1 not mask, 0 = both

int                                     iteration;

double                          selectedDistance;

double                          root2;

double                          slowness;
MNeigbourSet*                   neigbourSet;

MCandidateHeapList*             candidateList;

protected:

int                                     iterMax;
int                                     width;
int                                     height;

MDomainOfAction*                        dmap;                           

int                                     xSelected, ySelected, currentFront;

protected:


virtual double                  getPotential(int i, int j)  
                {return (dmap->getLabel(i,j) != Unknown ? (slowness ) : 1e12); }

// selection
virtual MPsi2DValue*            select(){ return candidateList->getTop(); }
virtual void                    getSelectionNeighbour(int x, int y);



// resolving the quadratic system (See Sethian & Kimmel)
virtual void                    resolve(int x, int y);  

public:

MDistanceFMM(int amaskFlag=1);
virtual ~MDistanceFMM(){ delete dmap; delete candidateList; delete neigbourSet; delete speedFactor; }

virtual void                    setMaskActivity(int amaskFlag){ maskFlag= amaskFlag; }

virtual double                  getDistance(int x, int y){return dmap->getValue(x, y);}
virtual int                             getFront(int x, int y){return dmap->getFront(x, y);}
virtual StatePsi                        getLabel(int x, int y){return dmap->getLabel(x, y);}
virtual double                  getSpeedCoeff(int fronNumber){return speedFactor[fronNumber];}


virtual void                    setData(int w, int h, double v0, int x0, int y0, double v1, int x1, int y1);
virtual void                    setData(int w, int h, int n, int* frontNum, double* v, int* xMask, int* yMask);
virtual void                    setData(int w, int h, int n0, double v0, int* xMask0, int* yMask0, 
                                                                        int n1, double v1, int* xMask1, int* yMask1);
virtual void                    setData(AnImage* im);
        

virtual int                             getMaxIterations(){ return iterMax; } 
virtual int                             numberOfIteration(){ return iteration; }
virtual void                    setIteration(int i){ iteration=i; }

virtual bool                    iterate();


virtual bool                    iterate(int* x, int* y)                         
                                        { 
                                        bool rc= iterate();

                                        *x= xSelected;
                                        *y= ySelected;

                                        return rc;
                                        }


virtual bool    iterate(int* x, int* y, double* rate)                                   
                                { 
                                bool rc= iterate(x,y);

                                *rate= dmap->getValue(xSelected,ySelected);

                                return rc;
                                }

virtual RImage*                 compute(RImage* dest=0);

virtual void                    sortedArea(AnImage* mask, 
                                                                   int* npoints, int** xposition, int** yposition, bool inside=true); 
virtual void                    sortedArea(AnImage* mask, int* npoints, 
                                                                int** xposition, int** yposition, double** rating, bool inside=true);

virtual AnImage*                filter(AnImage* mask, AnImage* dest=0);
virtual AnImage*                selectedFilter(AnImage* mask, AnImage* src, AnImage* dest=0);

virtual RImage*         rfilter(AnImage* mask, RImage* dest=0);

virtual AnImage*                greyLabelFilter(AnImage* src, AnImage* dest=0);
virtual AnImage*                colorLabelFilter(AnImage* src, AnImage* dest=0);

// euclidian Paths
virtual AnImage*                bidirPath(AnImage* src, double v0, int x0, int y0, double v1, int x1, int y1, AnImage* dest=0);

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.