projects / ctsim.git / blobdiff
? search:
summary | shortlog | log | commit | commitdiff | tree
raw | inline | side by side
r186: *** empty log message ***
[ctsim.git] / include / backprojectors.h
diff --git a/include/backprojectors.h b/include/backprojectors.h
index 598c2449ff26ec56edd94fff430e7a7d3bbd5b4c..b36dfb83cefbb7ef776ef5d4e2da3871e378a9ff 100644 (file)
--- 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.11 2000/07/20 11:17:31 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
@@ -35,79 +35,71 @@
 
 
 class Backproject;
+class ImageFile;
+class Projections;
 
 class Backprojector
 {
  public:
-  typedef  enum {
-    BPROJ_INVALID,
-    BPROJ_TRIG,
-    BPROJ_TABLE,
-    BPROJ_DIFF,
-    BPROJ_DIFF2,
-    BPROJ_IDIFF2,
-    BPROJ_IDIFF3
-} 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 
-    INTERP_FREQ_PREINTERPOLATION,
-  } InterpolationID;
-
-  static const char BPROJ_TRIG_STR[];
-  static const char BPROJ_TABLE_STR[];
-  static const char BPROJ_DIFF_STR[];
-  static const char BPROJ_DIFF2_STR[];
-  static const char BPROJ_IDIFF2_STR[];
-  static const char BPROJ_IDIFF3_STR[];
-
-  static const char BPROJ_TRIG_TITLE_STR[];
-  static const char BPROJ_TABLE_TITLE_STR[];
-  static const char BPROJ_DIFF_TITLE_STR[];
-  static const char BPROJ_DIFF2_TITLE_STR[];
-  static const char BPROJ_IDIFF2_TITLE_STR[];
-  static const char BPROJ_IDIFF3_TITLE_STR[];
- 
-  static const char INTERP_NEAREST_STR[];
-  static const char INTERP_LINEAR_STR[];
-  static const char INTERP_BSPLINE_STR[];
-  static const char INTERP_FREQ_PREINTERPOLATION_STR[];
-
-  static const char INTERP_NEAREST_TITLE_STR[];
-  static const char INTERP_LINEAR_TITLE_STR[];
-  static const char INTERP_BSPLINE_TITLE_STR[];
-  static const char INTERP_FREQ_PREINTERPOLATION_TITLE_STR[];
 
   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;}
-  const string& failMessage(void) const {return m_failMessage;}
+  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;
+  int m_idBackproject;
+  int m_idInterpolation;
   Backproject* m_pBackprojectImplem;
   bool m_fail;
   string m_failMessage;
 
-  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);
+  static const char* s_aszBackprojectName[];
+  static const char* s_aszBackprojectTitle[];
+  static const int s_iBackprojectCount;
+
+  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);
 };
@@ -116,20 +108,20 @@ class Backprojector
 class Backproject
 {
  public:
-    Backproject (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor);
+    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);
 
     const Projections& proj;
     ImageFile& im;
-    Backprojector::InterpolationID interpType;
+    int interpType;
     ImageFileArray v;
     kint32 nx;
     kint32 ny;
@@ -140,6 +132,7 @@ class Backproject
     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);
@@ -150,8 +143,8 @@ class Backproject
 class BackprojectTrig : public Backproject
 {
  public:
-  BackprojectTrig (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor)
-      : Backproject::Backproject (proj, im, interpType, interpFactor)
+  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);
@@ -161,12 +154,12 @@ class BackprojectTrig : public Backproject
 class BackprojectTable : public Backproject
 {
  public:
-  BackprojectTable (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor);
+  BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
   ~BackprojectTable ();
 
   void BackprojectView (const double* const t, double view_angle);
 
- private:
+ protected:
   Array2d<kfloat64> arrayR;
   Array2d<kfloat64> arrayPhi;
   kfloat64** r;
@@ -177,7 +170,7 @@ class BackprojectTable : public Backproject
 class BackprojectDiff : public Backproject
 {
  public:
-  BackprojectDiff (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor);
+  BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
   ~BackprojectDiff ();
 
   void BackprojectView (const double* const t, double view_angle);
@@ -191,8 +184,8 @@ class BackprojectDiff : public Backproject
 class BackprojectDiff2 : public BackprojectDiff
 {
  public:
-  BackprojectDiff2 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor)
-    :  BackprojectDiff::BackprojectDiff (proj, im, interpType, interpFactor)
+  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);
@@ -201,8 +194,8 @@ class BackprojectDiff2 : public BackprojectDiff
 class BackprojectIntDiff2 : public BackprojectDiff
 {
  public:
-  BackprojectIntDiff2 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor)
-    :  BackprojectDiff::BackprojectDiff (proj, im, interpType, interpFactor)
+  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);
@@ -212,12 +205,32 @@ class BackprojectIntDiff2 : public BackprojectDiff
 class BackprojectIntDiff3 : public BackprojectDiff
 {
  public:
-  BackprojectIntDiff3 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor)
-    :  BackprojectDiff::BackprojectDiff (proj, im, interpType, interpFactor)
+  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
Computed Tomography simulator and reconstructor