X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=blobdiff_plain;f=include%2Fbackprojectors.h;h=5e897cd22795ed259d00317ffcf0cc91e0ef4951;hp=c28a11b750c4ea15dc6260428650ecc62430061b;hb=1a050c98763fbbc0662731b0b76953acede6f5d7;hpb=f4a23943110823118f35756dd41fbd6707f04511 diff --git a/include/backprojectors.h b/include/backprojectors.h index c28a11b..5e897cd 100644 --- a/include/backprojectors.h +++ b/include/backprojectors.h @@ -1,25 +1,15 @@ /***************************************************************************** ** FILE IDENTIFICATION ** -** Name: backproject.h +** Name: backproject.h ** Purpose: Backprojection classes -** Programmer: Kevin Rosenberg -** Date Started: June 2000 +** Programmer: Kevin Rosenberg +** Date Started: June 2000 ** ** This is part of the CTSim program -** Copyright (C) 1983-2000 Kevin Rosenberg -** -** $Id: backprojectors.h,v 1.3 2000/06/17 20:12:14 kevin Exp $ -** $Log: backprojectors.h,v $ -** Revision 1.3 2000/06/17 20:12:14 kevin -** Converted Scanner and Projections to C++ -** -** Revision 1.2 2000/06/13 16:20:31 kevin -** finished c++ conversions -** -** Revision 1.1 2000/06/10 23:00:17 kevin -** *** empty log message *** +** Copyright (c) 1983-2001 Kevin Rosenberg ** +** $Id$ ** ** 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,32 +25,120 @@ ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA ******************************************************************************/ + +#ifndef __BACKPROJECTORS_H +#define __BACKPROJECTORS_H + +#undef HAVE_BSPLINE_INTERP + +#include "imagefile.h" + + +class Backproject; +class ImageFile; +class Projections; +struct ReconstructionROI; + +class Backprojector +{ + public: + static const int BPROJ_INVALID; + static const int BPROJ_TRIG; + static const int BPROJ_TABLE; + static const int BPROJ_DIFF; + static const int BPROJ_IDIFF; + + static const int INTERP_INVALID; + static const int INTERP_NEAREST; + static const int INTERP_LINEAR; + static const int INTERP_CUBIC; + static const int INTERP_FREQ_PREINTERPOLATION; +#if HAVE_BSPLINE_INTERP + static const int INTERP_BSPLINE; + static const int INTERP_1BSPLINE; + static const int INTERP_2BSPLINE; + static const int INTERP_3BSPLINE; +#endif + + 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;} + + static const int getBackprojectCount() {return s_iBackprojectCount;} + 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* 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); + + + private: + std::string m_nameBackproject; + std::string m_nameInterpolation; + int m_idBackproject; + int m_idInterpolation; + Backproject* m_pBackprojectImplem; + bool m_fail; + std::string m_failMessage; + + static const char* const s_aszBackprojectName[]; + static const char* const s_aszBackprojectTitle[]; + static const int s_iBackprojectCount; + + 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, const ReconstructionROI* pROI); +}; + + class Backproject { public: - Backproject (const Projections& proj, ImageFile& im, InterpolationType interpType); + Backproject (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, + const ReconstructionROI* pROI); virtual ~Backproject (); - virtual void BackprojectView (const double* const t, const double view_angle) {}; + virtual void BackprojectView (const double* const viewData, const double viewAngle) = 0; + virtual void PostProcessing (); // call after backprojecting all views protected: - void ScaleImageByRotIncrement (void); - 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); - - InterpolationType interpType; + void ScaleImageByRotIncrement (); + void errorIndexOutsideDetector (int ix, int iy, double theta, double r, double phi, double L, int ni); + void errorIndexOutsideDetector (int ix, int iy, double theta, double L, int ni); const Projections& proj; ImageFile& im; + int interpType; ImageFileArray v; - kuint32 nx; - kuint32 ny; + kint32 nx; + kint32 ny; double detInc; - double rotInc; - int iDetCenter; // index refering to L=0 projection + double rotScale; + 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 + double xInc, yInc; // size of cells + int m_interpFactor; + double m_dFocalLength; + double m_dSourceDetectorLength; + bool m_bPostProcessingDone; private: Backproject (const Backproject& rhs); @@ -71,23 +149,24 @@ class Backproject class BackprojectTrig : public Backproject { public: - BackprojectTrig (const Projections& proj, ImageFile& im, InterpolationType interpType) - : Backproject::Backproject (proj, im, interpType) + 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, double view_angle); + void BackprojectView (const double* const t, const double view_angle); }; class BackprojectTable : public Backproject { public: - BackprojectTable (const Projections& proj, ImageFile& im, InterpolationType interpType); - ~BackprojectTable (); + BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, const ReconstructionROI* pROI); + virtual ~BackprojectTable (); - void BackprojectView (const double* const t, double view_angle); + virtual void BackprojectView (const double* const t, const double view_angle); + virtual void PostProcessing (); // call after backprojecting all views - private: + protected: Array2d arrayR; Array2d arrayPhi; kfloat64** r; @@ -98,10 +177,11 @@ class BackprojectTable : public Backproject class BackprojectDiff : public Backproject { public: - BackprojectDiff (const Projections& proj, ImageFile& im, InterpolationType interpType); + BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor, const ReconstructionROI* pROI); ~BackprojectDiff (); - void BackprojectView (const double* const t, 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; @@ -109,25 +189,42 @@ class BackprojectDiff : public Backproject double im_xinc, im_yinc; }; -class BackprojectDiff2 : public BackprojectDiff + +class BackprojectIntDiff : public BackprojectDiff { public: - BackprojectDiff2 (const Projections& proj, ImageFile& im, InterpolationType interpType) - : BackprojectDiff::BackprojectDiff (proj, im, interpType) + 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, double view_angle); + void BackprojectView (const double* const t, const double view_angle); }; -class BackprojectIntDiff2 : public BackprojectDiff +class BackprojectEquilinear : public BackprojectTable { public: - BackprojectIntDiff2 (const Projections& proj, ImageFile& im, InterpolationType interpType) - : BackprojectDiff::BackprojectDiff (proj, im, interpType) - {} - - void BackprojectView (const double* const t, double view_angle); + 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); + + virtual ~BackprojectEquilinear() + {} +}; + +class BackprojectEquiangular : public BackprojectTable +{ + public: + 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); + + virtual ~BackprojectEquiangular() + {} }; -Backproject* selectBackprojector (BackprojType type, const Projections& proj, ImageFile& im, InterpolationType interpType); +#endif