** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: backprojectors.cpp,v 1.20 2000/12/18 09:31:26 kevin Exp $
+** $Id: backprojectors.cpp,v 1.22 2001/01/02 16:02:13 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
{
{"Nearest"},
{"Linear"},
- {"Frequency Preinterpolationj"},
+ {"Frequency Preinterpolation"},
#if HAVE_BSPLINE_INTERP
{"B-Spline"},
{"B-Spline 1st Order"},
iDetCenter = (nDet - 1) / 2; // index refering to L=0 projection
rotScale = proj.rotInc();
- if (proj.geometry() == Scanner::GEOMETRY_PARALLEL)\r
- rotScale /= (proj.nView() * proj.rotInc() / PI); // scale by number of PI rotations\r
- else if (proj.geometry() == Scanner::GEOMETRY_EQUIANGULAR || proj.geometry() == Scanner::GEOMETRY_EQUILINEAR)\r
- rotScale /= (proj.nView() * proj.rotInc() / (2 * PI)); // scale by number of 2PI rotations\r
- else\r
- sys_error (ERR_SEVERE, "Invalid geometry type %d [Backproject::Backproject]", proj.geometry());\r
+ if (proj.geometry() == Scanner::GEOMETRY_PARALLEL)
+ rotScale /= (proj.nView() * proj.rotInc() / PI); // scale by number of PI rotations
+ else if (proj.geometry() == Scanner::GEOMETRY_EQUIANGULAR || proj.geometry() == Scanner::GEOMETRY_EQUILINEAR)
+ rotScale /= (proj.nView() * proj.rotInc() / (2 * PI)); // scale by number of 2PI rotations
+ else
+ sys_error (ERR_SEVERE, "Invalid geometry type %d [Backproject::Backproject]", proj.geometry());
v = im.getArray();
nx = im.nx();
}
void Backproject::errorIndexOutsideDetector (int ix, int iy, double theta, double L, int iDetPos)
-{\r
+{
#if 1
std::ostringstream os;
os << "ix=" << ix << ", iy=" << iy << ", theta=" << theta << ", L=" << L << ", detinc=" << detInc << "\n";
os << "yMin=" << yMin << ", yMax=" << yMax << ", yInc=" << yInc << "\n";
os << "iDetPos index outside bounds: " << iDetPos << " [backprojector]";;
- sys_error (ERR_WARNING, os.str().c_str());\r
+ sys_error (ERR_WARNING, os.str().c_str());
#endif
}
*pImCol++ += filteredProj[iDetPos] + (detRemainder * deltaFilteredProj[iDetPos]);
} // end for iy
} //end linear
- } // end for ix\r
-\r
+ } // end for ix
+
delete deltaFilteredProj;
}