r156: *** empty log message ***

[ctsim.git] / include / backprojectors.h

diff --git a/include/backprojectors.h b/include/backprojectors.h
index 0a2a8e12c68b75ed0c91ffc17b42a97cb1e3a659..598c2449ff26ec56edd94fff430e7a7d3bbd5b4c 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.6 2000/07/07 15:30:59 kevin Exp $
+**  $Id: backprojectors.h,v 1.11 2000/07/20 11:17:31 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,8 +26,14 @@
 ******************************************************************************/
 
 
+#ifndef __BACKPROJECTORS_H
+#define __BACKPROJECTORS_H
+
 #undef HAVE_BSPLINE_INTERP
 
+#include "imagefile.h"
+
+
 class Backproject;
 
 class Backprojector
@@ -52,22 +58,35 @@ class Backprojector
     I_2BSPLINE,
     I_3BSPLINE,
 #endif
-    INTERP_LINEAR        // Linear interpolation 
+    INTERP_LINEAR,        // Linear interpolation 
+    INTERP_FREQ_PREINTERPOLATION,
   } InterpolationID;
 
-  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 BPROJ_IDIFF3_STR[]=   "idiff3";
+  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[]=  "nearest";
-  static const char INTERP_LINEAR_STR[]=   "linear";
-  static const char INTERP_BSPLINE_STR[]=  "bspline";
+  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);
+  Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName, const int interpFactor);
 
   ~Backprojector (void);
                 
@@ -90,14 +109,14 @@ class Backprojector
   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);
+  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, Backprojector::InterpolationID interpType, const int interpFactor);
 
     virtual ~Backproject ();
 
@@ -108,18 +127,19 @@ class Backproject
     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;
+    Backprojector::InterpolationID 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;
 
  private:
     Backproject (const Backproject& rhs);
@@ -130,8 +150,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, Backprojector::InterpolationID interpType, const int interpFactor)
+      : Backproject::Backproject (proj, im, interpType, interpFactor)
       {}
 
   void BackprojectView (const double* const t, double view_angle);
@@ -141,7 +161,7 @@ class BackprojectTrig : public Backproject
 class BackprojectTable : public Backproject
 {
  public:
-  BackprojectTable (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType);
+  BackprojectTable (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor);
   ~BackprojectTable ();
 
   void BackprojectView (const double* const t, double view_angle);
@@ -157,7 +177,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, Backprojector::InterpolationID interpType, const int interpFactor);
   ~BackprojectDiff ();
 
   void BackprojectView (const double* const t, double view_angle);
@@ -171,8 +191,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, Backprojector::InterpolationID interpType, const int interpFactor)
+    :  BackprojectDiff::BackprojectDiff (proj, im, interpType, interpFactor)
     {}
 
   void BackprojectView (const double* const t, double view_angle);
@@ -181,8 +201,8 @@ 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, Backprojector::InterpolationID interpType, const int interpFactor)
+    :  BackprojectDiff::BackprojectDiff (proj, im, interpType, interpFactor)
     {}
   
   void BackprojectView (const double* const t, double view_angle);
@@ -192,11 +212,12 @@ class BackprojectIntDiff2 : public BackprojectDiff
 class BackprojectIntDiff3 : public BackprojectDiff
 {
  public:
-  BackprojectIntDiff3 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType)
-    :  BackprojectDiff::BackprojectDiff (proj, im, interpType)
+  BackprojectIntDiff3 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType, const int interpFactor)
+    :  BackprojectDiff::BackprojectDiff (proj, im, interpType, interpFactor)
     {}
   
   void BackprojectView (const double* const t, double view_angle);
 };
 
 
+#endif