** This is part of the CTSim program
** Copyright (C) 1983-2000 Kevin Rosenberg
**
-** $Id: reconstruct.cpp,v 1.2 2000/12/06 01:46:43 kevin Exp $
+** $Id: reconstruct.cpp,v 1.5 2000/12/17 22:30:34 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 HAVE_SGP
int nVecFilter = m_pProcessSignal->getNFilterPoints();
- double adPlotXAxis [nVecFilter];
+ double* adPlotXAxis = new double [nVecFilter];
if (nVecFilter > 0 && pSGP) {
double f = m_pProcessSignal->getFilterMin();
ezplot.addCurve (adPlotXAxis, m_pProcessSignal->getFilter(), nVecFilter);
ezplot.plot();
}
- }
+ }\r
+ delete adPlotXAxis;
#endif
}
void
-Reconstructor::reconstructView (int iStartView, int iViewCount, SGP* pSGP)
+Reconstructor::reconstructView (int iStartView, int iViewCount, SGP* pSGP, bool bBackprojectView)
{
double* adFilteredProj = new double [m_nFilteredProjections]; // filtered projections
for (int iView = iStartView; iView < (iStartView + iViewCount); iView++) {
if (m_iTrace == Trace::TRACE_CONSOLE)
- cout <<"Reconstructing view " << iView << " (last = " << m_rProj.nView() - 1 << ")\n";
+ std::cout <<"Reconstructing view " << iView << " (last = " << m_rProj.nView() - 1 << ")\n";
const DetectorArray& rDetArray = m_rProj.getDetectorArray (iView);
const DetectorValue* detval = rDetArray.detValues();
}
#endif
#endif
-
- m_pBackprojector->BackprojectView (adFilteredProj, rDetArray.viewAngle());
+\r
+ if (bBackprojectView)
+ m_pBackprojector->BackprojectView (adFilteredProj, rDetArray.viewAngle());
#ifdef HAVE_SGP
if (m_iTrace >= Trace::TRACE_PLOT && pSGP) {