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 bd30b9cf47c0a4739da2f301e023505506d0f26f..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.22 2001/02/22 18:22:40 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,6 +64,7 @@ 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;}
@@ -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;