X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=include%2Fbackprojectors.h;h=d8f9bb56d5715b3d26054c03530e94af42638c8e;hb=5a6caa64e880f613b82e516031028d02fd127257;hp=bd30b9cf47c0a4739da2f301e023505506d0f26f;hpb=9776c9a12ba53419d34563a5ec57c90e3d6798f4;p=ctsim.git

diff --git a/include/backprojectors.h b/include/backprojectors.h
index bd30b9c..d8f9bb5 100644
--- 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.25 2001/03/11 15:27:30 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
@@ -37,6 +37,7 @@
 class Backproject;
 class ImageFile;
 class Projections;
+struct ReconstructionROI;
 
 class Backprojector
 {
@@ -59,11 +60,13 @@ class Backprojector
   static const int INTERP_3BSPLINE;
 #endif
 
-  Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName, const int interpFactor);
+  Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, 
+    const char* const interpName, const int interpFactor, const ReconstructionROI* pROI);
 
   ~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;}
@@ -100,24 +103,26 @@ class Backprojector
   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);
+  bool initBackprojector (const Projections& proj, ImageFile& im, const char* const backprojName,
+    const char* const interpName, const int interpFactor, const ReconstructionROI* pROI);
 };
 
 
 class Backproject
 {
  public:
-    Backproject (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
+    Backproject (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, 
+      const ReconstructionROI* pROI);
 
     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 +137,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);
@@ -142,8 +149,8 @@ class Backproject
 class BackprojectTrig : public Backproject
 {
  public:
-  BackprojectTrig (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
-      : Backproject (proj, im, interpID, interpFactor)
+  BackprojectTrig (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, const ReconstructionROI* pROI)
+      : Backproject (proj, im, interpID, interpFactor, pROI)
       {}
 
   void BackprojectView (const double* const t, const double view_angle);
@@ -153,10 +160,11 @@ class BackprojectTrig : public Backproject
 class BackprojectTable : public Backproject
 {
  public:
-  BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
+  BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, const ReconstructionROI* pROI);
   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;
@@ -169,10 +177,11 @@ class BackprojectTable : public Backproject
 class BackprojectDiff : public Backproject
 {
  public:
-  BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
+  BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, const ReconstructionROI* pROI);
   ~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;
@@ -184,8 +193,8 @@ class BackprojectDiff : public Backproject
 class BackprojectIntDiff : public BackprojectDiff
 {
  public:
-  BackprojectIntDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
-    :  BackprojectDiff (proj, im, interpID, interpFactor)
+  BackprojectIntDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, const ReconstructionROI* pROI)
+    :  BackprojectDiff (proj, im, interpID, interpFactor, pROI)
     {}
   
   void BackprojectView (const double* const t, const double view_angle);
@@ -194,8 +203,8 @@ class BackprojectIntDiff : public BackprojectDiff
 class BackprojectEquilinear : public BackprojectTable
 {
  public:
-  BackprojectEquilinear (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
-      : BackprojectTable (proj, im, interpID, interpFactor)
+  BackprojectEquilinear (const Projections& proj, ImageFile& im, int interpID, const int interpFactor,  const ReconstructionROI* pROI)
+      : BackprojectTable (proj, im, interpID, interpFactor, pROI)
       {}
 
   void BackprojectView (const double* const t, const double view_angle);
@@ -207,8 +216,8 @@ class BackprojectEquilinear : public BackprojectTable
 class BackprojectEquiangular : public BackprojectTable
 {
  public:
-  BackprojectEquiangular (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
-      : BackprojectTable (proj, im, interpID, interpFactor)
+  BackprojectEquiangular (const Projections& proj, ImageFile& im, int interpID, const int interpFactor,  const ReconstructionROI* pROI)
+      : BackprojectTable (proj, im, interpID, interpFactor, pROI)
       {}
 
   void BackprojectView (const double* const t, const double view_angle);