** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: projections.cpp,v 1.13 2000/07/11 10:32:44 kevin Exp $
+** $Id: projections.cpp,v 1.15 2000/07/15 08:36: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
if (! fs)
return false;
- off_t testPos;
fs.readInt16 (_hsize);
fs.readInt16 (_signature);
fs.readInt32 (_nView);
darray.setViewAngle (view_angle);
// darray.setNDet ( nDet);
- for (int i = 0; i < nDet; i++) {
+ for (unsigned int i = 0; i < nDet; i++) {
kfloat32 detval;
fs.readFloat32 (detval);
detval_ptr[i] = detval;
fs.writeFloat64 (view_angle);
fs.writeInt32 (nDet);
- for (int i = 0; i < nDet; i++) {
+ for (unsigned int i = 0; i < nDet; i++) {
kfloat32 detval = detval_ptr[i];
fs.writeFloat32 (detval);
}
*/
bool
-Projections::reconstruct (ImageFile& im, const char* const filterName, double filt_param, const char* const filterMethodName, const int zeropad, const char* const interpName, int interpFactor, const char* const backprojectName, const int trace)
+Projections::reconstruct (ImageFile& im, const char* const filterName, double filt_param, const char* const filterMethodName, const int zeropad, const char* const interpName, int interpFactor, const char* const backprojectName, const int trace) const
{
- int nview = m_nView;
double detInc = m_detInc;
int n_filteredProj = m_nDet * interpFactor;
double filteredProj [n_filteredProj]; // filtered projections
if (trace >= TRACE_TEXT)
printf ("Reconstructing view %d (last = %d)\n", iview, m_nView - 1);
- DetectorArray& darray = getDetectorArray (iview);
- DetectorValue* detval = darray.detValues();
+ const DetectorArray& darray = getDetectorArray (iview);
+ const DetectorValue* detval = darray.detValues();
filter.filterSignal (detval, filteredProj);