1 /*****************************************************************************
5 ** Purpose: Backprojection classes
6 ** Programmer: Kevin Rosenberg
7 ** Date Started: June 2000
9 ** This is part of the CTSim program
10 ** Copyright (C) 1983-2000 Kevin Rosenberg
12 ** $Id: backprojectors.h,v 1.7 2000/07/09 08:16:17 kevin Exp $
14 ** This program is free software; you can redistribute it and/or modify
15 ** it under the terms of the GNU General Public License (version 2) as
16 ** published by the Free Software Foundation.
18 ** This program is distributed in the hope that it will be useful,
19 ** but WITHOUT ANY WARRANTY; without even the implied warranty of
20 ** MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
21 ** GNU General Public License for more details.
23 ** You should have received a copy of the GNU General Public License
24 ** along with this program; if not, write to the Free Software
25 ** Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
26 ******************************************************************************/
29 #undef HAVE_BSPLINE_INTERP
48 INTERP_NEAREST, // Nearest neighbor
49 #if HAVE_BSPLINE_INTERP
51 I_1BSPLINE, // 1st order B-Spline
55 INTERP_LINEAR, // Linear interpolation
56 INTERP_FREQ_PREINTERPOLATE,
59 static const char BPROJ_TRIG_STR[]= "trig";
60 static const char BPROJ_TABLE_STR[]= "table";
61 static const char BPROJ_DIFF_STR[]= "diff";
62 static const char BPROJ_DIFF2_STR[]= "diff2";
63 static const char BPROJ_IDIFF2_STR[]= "idiff2";
64 static const char BPROJ_IDIFF3_STR[]= "idiff3";
66 static const char INTERP_NEAREST_STR[]= "nearest";
67 static const char INTERP_LINEAR_STR[]= "linear";
68 static const char INTERP_BSPLINE_STR[]= "bspline";
69 static const char INTERP_FREQ_PREINTERPOLATE_STR[]= "freq_preinterpolate";
71 Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName);
73 ~Backprojector (void);
75 void BackprojectView (const double* const viewData, const double viewAngle);
77 bool fail(void) const {return m_fail;}
78 const string& failMessage(void) const {return m_failMessage;}
81 string m_nameBackproject;
82 string m_nameInterpolation;
83 BackprojectID m_idBackproject;
84 InterpolationID m_idInterpolation;
85 Backproject* m_pBackprojectImplem;
89 static const InterpolationID convertInterpolationNameToID (const char* const interpName);
90 static const char* convertInterpolationIDToName (const InterpolationID interpID);
91 static const BackprojectID convertBackprojectNameToID (const char* const bprojName);
92 static const char* convertBackprojectIDToName (const BackprojectID bprojID);
94 bool initBackprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName);
101 Backproject (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType);
103 virtual ~Backproject ();
105 virtual void BackprojectView (const double* const viewData, const double viewAngle) {};
108 void ScaleImageByRotIncrement (void);
109 void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double r, double phi, double L, int ni);
110 void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double L, int ni);
112 Backprojector::InterpolationID interpType;
113 const Projections& proj;
120 int iDetCenter; // index refering to L=0 projection
122 double xMin, xMax, yMin, yMax; // Retangular coords of phantom
123 double xInc, yInc; // size of cells
126 Backproject (const Backproject& rhs);
127 Backproject& operator= (const Backproject& rhs);
131 class BackprojectTrig : public Backproject
134 BackprojectTrig (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType)
135 : Backproject::Backproject (proj, im, interpType)
138 void BackprojectView (const double* const t, double view_angle);
142 class BackprojectTable : public Backproject
145 BackprojectTable (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType);
146 ~BackprojectTable ();
148 void BackprojectView (const double* const t, double view_angle);
151 Array2d<kfloat64> arrayR;
152 Array2d<kfloat64> arrayPhi;
158 class BackprojectDiff : public Backproject
161 BackprojectDiff (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType);
164 void BackprojectView (const double* const t, double view_angle);
169 double im_xinc, im_yinc;
172 class BackprojectDiff2 : public BackprojectDiff
175 BackprojectDiff2 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType)
176 : BackprojectDiff::BackprojectDiff (proj, im, interpType)
179 void BackprojectView (const double* const t, double view_angle);
182 class BackprojectIntDiff2 : public BackprojectDiff
185 BackprojectIntDiff2 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType)
186 : BackprojectDiff::BackprojectDiff (proj, im, interpType)
189 void BackprojectView (const double* const t, double view_angle);
193 class BackprojectIntDiff3 : public BackprojectDiff
196 BackprojectIntDiff3 (const Projections& proj, ImageFile& im, Backprojector::InterpolationID interpType)
197 : BackprojectDiff::BackprojectDiff (proj, im, interpType)
200 void BackprojectView (const double* const t, double view_angle);