r591: Added Center-Detector length to scanning and reconstruction

[ctsim.git] / include / backprojectors.h

diff --git a/include/backprojectors.h b/include/backprojectors.h
index 8fbb06b86382652f2af15fbc3cefd328ef646120..0a0cb1275dfdd8034db2b30df0d8c4a73d562c9c 100644 (file)
--- a/include/backprojectors.h
+++ b/include/backprojectors.h
@@ -9,7 +9,7 @@
 **  This is part of the CTSim program
 **  Copyright (c) 1983-2001 Kevin Rosenberg
 **
-**  $Id: backprojectors.h,v 1.21 2001/02/11 04:56:37 kevin Exp $
+**  $Id: backprojectors.h,v 1.24 2001/03/01 07:30:49 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
@@ -64,20 +64,21 @@ class Backprojector
   ~Backprojector ();
                 
   void BackprojectView (const double* const viewData, const double viewAngle);
+  void PostProcessing();
 
   bool fail() const {return m_fail;}
   const std::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 const char* const* getBackprojectNameArray() {return s_aszBackprojectName;}
+  static const char* const* 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 const char* const * getInterpNameArray() {return s_aszInterpName;}
+  static const char* const * 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);
@@ -92,12 +93,12 @@ class Backprojector
   bool m_fail;
   std::string m_failMessage;
 
-  static const char* s_aszBackprojectName[];
-  static const char* s_aszBackprojectTitle[];
+  static const char* const s_aszBackprojectName[];
+  static const char* const s_aszBackprojectTitle[];
   static const int s_iBackprojectCount;
 
-  static const char* s_aszInterpName[];
-  static const char* s_aszInterpTitle[];
+  static const char* const s_aszInterpName[];
+  static const char* const 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);
@@ -111,13 +112,13 @@ class Backproject
 
     virtual ~Backproject ();
 
-    virtual void BackprojectView (const double* const viewData, const double viewAngle) {};
+    virtual void BackprojectView (const double* const viewData, const double viewAngle) = 0;
+    virtual void PostProcessing (); // call after backprojecting all views
 
  protected:
     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;
     int interpType;
@@ -132,6 +133,8 @@ class Backproject
     double xInc, yInc; // size of cells
     int m_interpFactor;
     double m_dFocalLength;
+    double m_dSourceDetectorLength;
+    bool m_bPostProcessingDone;
 
  private:
     Backproject (const Backproject& rhs);
@@ -156,7 +159,8 @@ class BackprojectTable : public Backproject
   BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
   virtual ~BackprojectTable ();
 
-  void BackprojectView (const double* const t, const double view_angle);
+  virtual void BackprojectView (const double* const t, const double view_angle);
+  virtual void PostProcessing (); // call after backprojecting all views
 
  protected:
   Array2d<kfloat64> arrayR;
@@ -172,7 +176,8 @@ class BackprojectDiff : public Backproject
   BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
   ~BackprojectDiff ();
 
-  void BackprojectView (const double* const t, const double view_angle);
+  virtual void BackprojectView (const double* const t, const double view_angle);
+  virtual void PostProcessing (); // call after backprojecting all views
 
  protected:
   double start_r;