00001 #ifndef __EikonalFastSweepingOp__
00002 #define __EikonalFastSweepingOp__
00003 
00004 
00005 
00034 #include <stdio.h>
00035 #include <math.h>
00036 
00037 #include "AbsEikonalFastSweepingOp.hpp"
00038 
00039 
00040 class EikonalFastSweepingOp : public  AbsEikonalFastSweepingOp
00041 {
00042 protected:
00043 
00044 RImage*         P;              
00045 
00046 protected:
00047 
00048 virtual double getPotential_h(int i, int j){ return P->get(i,j)*h; }   
00049 
00050 public:
00051         
00052 EikonalFastSweepingOp(int iter=1, RImage* imageMap=0, double hx=1.0):AbsEikonalFastSweepingOp(iter, hx)
00053         { P= imageMap; }
00054 virtual ~EikonalFastSweepingOp(){}
00055 
00056 virtual void setPotential(RImage* imageMap){P= imageMap}                
00057 
00058 virtual void report(FILE* file){ fprintf(file,"Distance Fast Sweeping Operator\n"); }
00059 };
00060 
00061 
00062 #endif