X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fbackprojectors.h;h=b36dfb83cefbb7ef776ef5d4e2da3871e378a9ff;hp=08b72d450fbfb121a78532d1a1b37e1bbf9bb948;hb=a05f3cb550877e94aa118cc04b361c0c8fdb3dc3;hpb=2d39e823ba389fc68e5317c422b55be006094252 diff --git a/include/backprojectors.h b/include/backprojectors.h index 08b72d4..b36dfb8 100644 --- a/include/backprojectors.h +++ b/include/backprojectors.h @@ -9,7 +9,7 @@ ** This is part of the CTSim program ** Copyright (C) 1983-2000 Kevin Rosenberg ** -** $Id: backprojectors.h,v 1.4 2000/06/22 10:17:28 kevin Exp $ +** $Id: backprojectors.h,v 1.14 2000/08/31 08:38:58 kevin Exp $ ** ** This program is free software; you can redistribute it and/or modify ** it under the terms of the GNU General Public License (version 2) as @@ -26,95 +26,113 @@ ******************************************************************************/ +#ifndef __BACKPROJECTORS_H +#define __BACKPROJECTORS_H + #undef HAVE_BSPLINE_INTERP +#include "imagefile.h" + + class Backproject; +class ImageFile; +class Projections; class Backprojector { public: - typedef enum { - BPROJ_INVALID, - BPROJ_TRIG, - BPROJ_TABLE, - BPROJ_DIFF, - BPROJ_DIFF2, - BPROJ_IDIFF2 -} BackprojectID; - - typedef enum { - INTERP_INVALID, - INTERP_NEAREST, // Nearest neighbor + static const int BPROJ_INVALID; + static const int BPROJ_TRIG; + static const int BPROJ_TABLE; + static const int BPROJ_DIFF; + static const int BPROJ_DIFF2; + static const int BPROJ_IDIFF2; + static const int BPROJ_IDIFF3; + + static const int INTERP_INVALID; + static const int INTERP_NEAREST; + static const int INTERP_LINEAR; + static const int INTERP_FREQ_PREINTERPOLATION; #if HAVE_BSPLINE_INTERP - I_BSPLINE, - I_1BSPLINE, // 1st order B-Spline - I_2BSPLINE, - I_3BSPLINE, + static const int INTERP_BSPLINE; + static const int INTERP_1BSPLINE; + static const int INTERP_2BSPLINE; + static const int INTERP_3BSPLINE; #endif - INTERP_LINEAR // Linear interpolation - } InterpolationID; - Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName); + Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName, const int interpFactor); - ~Backprojector (void); + ~Backprojector (); void BackprojectView (const double* const viewData, const double viewAngle); - bool fail(void) const {return m_fail;} + bool fail() const {return m_fail;} + const string& failMessage() const {return m_failMessage;} + + static const int getBackprojectCount() {return s_iBackprojectCount;} + static const char** getBackprojectNameArray() {return s_aszBackprojectName;} + static const char** getBackprojectTitleArray() {return s_aszBackprojectTitle;} + static int convertBackprojectNameToID (const char* const bprojName); + static const char* convertBackprojectIDToName (const int bprojID); + static const char* convertBackprojectIDToTitle (const int bprojID); + + static const int getInterpCount() {return s_iInterpCount;} + static const char** getInterpNameArray() {return s_aszInterpName;} + static const char** getInterpTitleArray() {return s_aszInterpTitle;} + static int convertInterpNameToID (const char* const interpName); + static const char* convertInterpIDToName (const int interpID); + static const char* convertInterpIDToTitle (const int interpID); + private: string m_nameBackproject; string m_nameInterpolation; - BackprojectID m_idBackproject; - InterpolationID m_idInterpolation; - bool m_fail; + int m_idBackproject; + int m_idInterpolation; Backproject* m_pBackprojectImplem; + bool m_fail; + string m_failMessage; + + static const char* s_aszBackprojectName[]; + static const char* s_aszBackprojectTitle[]; + static const int s_iBackprojectCount; - static const char BPROJ_TRIG_STR[]= "trig"; - static const char BPROJ_TABLE_STR[]= "table"; - static const char BPROJ_DIFF_STR[]= "diff"; - static const char BPROJ_DIFF2_STR[]= "diff2"; - static const char BPROJ_IDIFF2_STR[]= "idiff2"; - - static const char INTERP_NEAREST_STR[]= "nearest"; - static const char INTERP_LINEAR_STR[]= "linear"; - static const char INTERP_BSPLINE_STR[]= "bspline"; - - static const InterpolationID convertInterpolationNameToID (const char* const interpName); - static const char* convertInterpolationIDToName (const InterpolationID interpID); - static const BackprojectID convertBackprojectNameToID (const char* const bprojName); - static const char* convertBackprojectIDToName (const BackprojectID bprojID); - - bool initBackprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName); + static const char* s_aszInterpName[]; + static const char* s_aszInterpTitle[]; + static const int s_iInterpCount; + + bool initBackprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName, const int interpFactor); }; class Backproject { public: - Backproject (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType); + Backproject (const Projections& proj, ImageFile& im, int interpID, const int interpFactor); virtual ~Backproject (); virtual void BackprojectView (const double* const viewData, const double viewAngle) {}; protected: - void ScaleImageByRotIncrement (void); + void ScaleImageByRotIncrement (); void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double r, double phi, double L, int ni); void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double L, int ni); - Backprojector::InterpolationID interpType; const Projections& proj; ImageFile& im; + int interpType; ImageFileArray v; - kuint32 nx; - kuint32 ny; + kint32 nx; + kint32 ny; double detInc; double rotInc; int iDetCenter; // index refering to L=0 projection int nDet; double xMin, xMax, yMin, yMax; // Retangular coords of phantom double xInc, yInc; // size of cells + int m_interpFactor; + double m_dFocalLength; private: Backproject (const Backproject& rhs); @@ -125,8 +143,8 @@ class Backproject class BackprojectTrig : public Backproject { public: - BackprojectTrig (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType) - : Backproject::Backproject (proj, im, interpType) + BackprojectTrig (const Projections& proj, ImageFile& im, int interpID, const int interpFactor) + : Backproject::Backproject (proj, im, interpID, interpFactor) {} void BackprojectView (const double* const t, double view_angle); @@ -136,12 +154,12 @@ class BackprojectTrig : public Backproject class BackprojectTable : public Backproject { public: - BackprojectTable (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType); + BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor); ~BackprojectTable (); void BackprojectView (const double* const t, double view_angle); - private: + protected: Array2d arrayR; Array2d arrayPhi; kfloat64** r; @@ -152,7 +170,7 @@ class BackprojectTable : public Backproject class BackprojectDiff : public Backproject { public: - BackprojectDiff (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType); + BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor); ~BackprojectDiff (); void BackprojectView (const double* const t, double view_angle); @@ -166,8 +184,8 @@ class BackprojectDiff : public Backproject class BackprojectDiff2 : public BackprojectDiff { public: - BackprojectDiff2 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType) - : BackprojectDiff::BackprojectDiff (proj, im, interpType) + BackprojectDiff2 (const Projections& proj, ImageFile& im, int interpID, const int interpFactor) + : BackprojectDiff::BackprojectDiff (proj, im, interpID, interpFactor) {} void BackprojectView (const double* const t, double view_angle); @@ -176,11 +194,43 @@ class BackprojectDiff2 : public BackprojectDiff class BackprojectIntDiff2 : public BackprojectDiff { public: - BackprojectIntDiff2 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType) - : BackprojectDiff::BackprojectDiff (proj, im, interpType) + BackprojectIntDiff2 (const Projections& proj, ImageFile& im, int interpID, const int interpFactor) + : BackprojectDiff::BackprojectDiff (proj, im, interpID, interpFactor) {} void BackprojectView (const double* const t, double view_angle); }; +class BackprojectIntDiff3 : public BackprojectDiff +{ + public: + BackprojectIntDiff3 (const Projections& proj, ImageFile& im, int interpID, const int interpFactor) + : BackprojectDiff::BackprojectDiff (proj, im, interpID, interpFactor) + {} + + void BackprojectView (const double* const t, double view_angle); +}; + +class BackprojectEquilinear : public BackprojectTable +{ + public: + BackprojectEquilinear (const Projections& proj, ImageFile& im, int interpID, const int interpFactor) + : BackprojectTable::BackprojectTable (proj, im, interpID, interpFactor) + {} + + void BackprojectView (const double* const t, double view_angle); +}; + +class BackprojectEquiangular : public BackprojectTable +{ + public: + BackprojectEquiangular (const Projections& proj, ImageFile& im, int interpID, const int interpFactor) + : BackprojectTable::BackprojectTable (proj, im, interpID, interpFactor) + {} + + void BackprojectView (const double* const t, double view_angle); +}; + + +#endif