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; // potential 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