** Date Started: Aug 84
**
** This is part of the CTSim program
-** Copyright (C) 1983-2000 Kevin Rosenberg
+** Copyright (c) 1983-2001 Kevin Rosenberg
**
-** $Id: reconstruct.cpp,v 1.3 2000/12/16 02:44:26 kevin Exp $
+** $Id: reconstruct.cpp,v 1.9 2001/01/28 19:10:18 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
adPlotXAxis[i] = f;
if (m_pProcessSignal->getFilter()) {
- EZPlot ezplot (*pSGP);
+ EZPlot ezplot;
ezplot.ezset ("title Filter Response");
ezplot.addCurve (adPlotXAxis, m_pProcessSignal->getFilter(), nVecFilter);
- ezplot.plot();
+ ezplot.plot (pSGP);
}
- }\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());
+ if (bBackprojectView)
+ m_pBackprojector->BackprojectView (adFilteredProj, rDetArray.viewAngle());
#ifdef HAVE_SGP
if (m_iTrace >= Trace::TRACE_PLOT && pSGP) {
- EZPlot ezplotProj (*pSGP);
+ EZPlot ezplotProj;
ezplotProj.ezset ("clear");
ezplotProj.ezset ("title Raw Projection");
ezplotProj.ezset ("xticks major 5");
ezplotProj.ezset ("xlabel ");
ezplotProj.ezset ("ylabel ");
- ezplotProj.ezset ("yporigin .5");
- ezplotProj.ezset ("ylength .5");
+ ezplotProj.ezset ("yporigin 0.1");
+ ezplotProj.ezset ("ylength 0.45");
+ ezplotProj.ezset ("xlength 0.66");
ezplotProj.ezset ("box.");
ezplotProj.ezset ("grid.");
ezplotProj.addCurve (m_adPlotXAxis, detval, m_rProj.nDet());
- ezplotProj.plot();
+ pSGP->setTextPointSize(12);
+ ezplotProj.plot (pSGP);
ezplotProj.ezset ("clear");
ezplotProj.ezset ("title Filtered Projection");
ezplotProj.ezset ("xticks major 5");
ezplotProj.ezset ("xlabel ");
ezplotProj.ezset ("ylabel ");
- ezplotProj.ezset ("ylength .5");
+ ezplotProj.ezset ("yporigin 0.55");
+ ezplotProj.ezset ("ylength 0.45");
+ ezplotProj.ezset ("xlength 0.66");
ezplotProj.ezset ("box");
ezplotProj.ezset ("grid");
ezplotProj.addCurve (m_adPlotXAxis, adFilteredProj, m_nFilteredProjections);
- ezplotProj.plot();
+ pSGP->setTextPointSize(12);
+ ezplotProj.plot (pSGP);
}
#endif //HAVE_SGP
}
-\r
- delete adFilteredProj;\r
+
+ delete adFilteredProj;
}