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.16 2000/12/03 15: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 #ifndef __BACKPROJECTORS_H
30 #define __BACKPROJECTORS_H
32 #undef HAVE_BSPLINE_INTERP
34 #include "imagefile.h"
44 static const int BPROJ_INVALID;
45 static const int BPROJ_TRIG;
46 static const int BPROJ_TABLE;
47 static const int BPROJ_DIFF;
48 static const int BPROJ_DIFF2;
49 static const int BPROJ_IDIFF2;
50 static const int BPROJ_IDIFF3;
52 static const int INTERP_INVALID;
53 static const int INTERP_NEAREST;
54 static const int INTERP_LINEAR;
55 static const int INTERP_FREQ_PREINTERPOLATION;
56 #if HAVE_BSPLINE_INTERP
57 static const int INTERP_BSPLINE;
58 static const int INTERP_1BSPLINE;
59 static const int INTERP_2BSPLINE;
60 static const int INTERP_3BSPLINE;
63 Backprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName, const int interpFactor);
67 void BackprojectView (const double* const viewData, const double viewAngle);
69 bool fail() const {return m_fail;}
70 const string& failMessage() const {return m_failMessage;}
72 static const int getBackprojectCount() {return s_iBackprojectCount;}
73 static const char** getBackprojectNameArray() {return s_aszBackprojectName;}
74 static const char** getBackprojectTitleArray() {return s_aszBackprojectTitle;}
75 static int convertBackprojectNameToID (const char* const bprojName);
76 static const char* convertBackprojectIDToName (const int bprojID);
77 static const char* convertBackprojectIDToTitle (const int bprojID);
79 static const int getInterpCount() {return s_iInterpCount;}
80 static const char** getInterpNameArray() {return s_aszInterpName;}
81 static const char** getInterpTitleArray() {return s_aszInterpTitle;}
82 static int convertInterpNameToID (const char* const interpName);
83 static const char* convertInterpIDToName (const int interpID);
84 static const char* convertInterpIDToTitle (const int interpID);
88 string m_nameBackproject;
89 string m_nameInterpolation;
91 int m_idInterpolation;
92 Backproject* m_pBackprojectImplem;
96 static const char* s_aszBackprojectName[];
97 static const char* s_aszBackprojectTitle[];
98 static const int s_iBackprojectCount;
100 static const char* s_aszInterpName[];
101 static const char* s_aszInterpTitle[];
102 static const int s_iInterpCount;
104 bool initBackprojector (const Projections& proj, ImageFile& im, const char* const backprojName, const char* const interpName, const int interpFactor);
111 Backproject (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
113 virtual ~Backproject ();
115 virtual void BackprojectView (const double* const viewData, const double viewAngle) {};
118 void ScaleImageByRotIncrement ();
119 void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double r, double phi, double L, int ni);
120 void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double L, int ni);
122 const Projections& proj;
130 int iDetCenter; // index refering to L=0 projection
132 double xMin, xMax, yMin, yMax; // Retangular coords of phantom
133 double xInc, yInc; // size of cells
135 double m_dFocalLength;
138 Backproject (const Backproject& rhs);
139 Backproject& operator= (const Backproject& rhs);
143 class BackprojectTrig : public Backproject
146 BackprojectTrig (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
147 : Backproject::Backproject (proj, im, interpID, interpFactor)
150 void BackprojectView (const double* const t, double view_angle);
154 class BackprojectTable : public Backproject
157 BackprojectTable (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
158 virtual ~BackprojectTable ();
160 void BackprojectView (const double* const t, double view_angle);
163 Array2d<kfloat64> arrayR;
164 Array2d<kfloat64> arrayPhi;
170 class BackprojectDiff : public Backproject
173 BackprojectDiff (const Projections& proj, ImageFile& im, int interpID, const int interpFactor);
176 void BackprojectView (const double* const t, double view_angle);
181 double im_xinc, im_yinc;
184 class BackprojectDiff2 : public BackprojectDiff
187 BackprojectDiff2 (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
188 : BackprojectDiff::BackprojectDiff (proj, im, interpID, interpFactor)
191 void BackprojectView (const double* const t, double view_angle);
194 class BackprojectIntDiff2 : public BackprojectDiff
197 BackprojectIntDiff2 (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
198 : BackprojectDiff::BackprojectDiff (proj, im, interpID, interpFactor)
201 void BackprojectView (const double* const t, double view_angle);
205 class BackprojectIntDiff3 : public BackprojectDiff
208 BackprojectIntDiff3 (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
209 : BackprojectDiff::BackprojectDiff (proj, im, interpID, interpFactor)
212 void BackprojectView (const double* const t, double view_angle);
215 class BackprojectEquilinear : public BackprojectTable
218 BackprojectEquilinear (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
219 : BackprojectTable::BackprojectTable (proj, im, interpID, interpFactor)
222 void BackprojectView (const double* const t, double view_angle);
224 virtual ~BackprojectEquilinear()
228 class BackprojectEquiangular : public BackprojectTable
231 BackprojectEquiangular (const Projections& proj, ImageFile& im, int interpID, const int interpFactor)
232 : BackprojectTable::BackprojectTable (proj, im, interpID, interpFactor)
235 void BackprojectView (const double* const t, double view_angle);
237 virtual ~BackprojectEquiangular()