From: Kevin M. Rosenberg Date: Wed, 7 Mar 2001 19:02:44 +0000 (+0000) Subject: r616: *** empty log message *** X-Git-Tag: debian-4.5.3-3~401 X-Git-Url: http://git.kpe.io/?p=ctsim.git;a=commitdiff_plain;h=b7cbad1230585741079fad4881a8d67be2592100 r616: *** empty log message *** --- diff --git a/libctsim/backprojectors.cpp b/libctsim/backprojectors.cpp index bab1dc7..bfd594e 100644 --- a/libctsim/backprojectors.cpp +++ b/libctsim/backprojectors.cpp @@ -8,7 +8,7 @@ ** This is part of the CTSim program ** Copyright (c) 1983-2001 Kevin Rosenberg ** -** $Id: backprojectors.cpp,v 1.29 2001/03/01 07:30:49 kevin Exp $ +** $Id: backprojectors.cpp,v 1.30 2001/03/07 19:02:44 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 @@ -365,7 +365,7 @@ BackprojectTrig::BackprojectView (const double* const filteredProj, const double double frac = p - pFloor; // fraction distance from det if (iDetPos >= 0 && iDetPos < nDet - 1) v[ix][iy] += rotScale * ((1-frac) * filteredProj[iDetPos] + frac * filteredProj[iDetPos+1]); - } else if (interpType = Backprojector::INTERP_CUBIC) { + } else if (interpType == Backprojector::INTERP_CUBIC) { double p = iDetCenter + (L / detInc); // position along detector if (p >= 0 && p < nDet) v[ix][iy] += rotScale * pCubicInterp->interpolate (p); @@ -441,7 +441,7 @@ BackprojectTable::BackprojectView (const double* const filteredProj, const doubl double frac = dPos - dPosFloor; // fraction distance from det if (iDetPos >= 0 && iDetPos < nDet - 1) pImCol[iy] += ((1-frac) * filteredProj[iDetPos] + frac * filteredProj[iDetPos+1]); - } else if (interpType = Backprojector::INTERP_CUBIC) { + } else if (interpType == Backprojector::INTERP_CUBIC) { double p = iDetCenter + (L / detInc); // position along detector if (p >= 0 && p < nDet) pImCol[iy] += pCubicInterp->interpolate (p); @@ -518,7 +518,7 @@ BackprojectDiff::BackprojectView (const double* const filteredProj, const double double frac = curDetPos - detPosFloor; // fraction distance from det if (iDetPos > 0 && iDetPos < nDet - 1) *pImCol++ += filteredProj[iDetPos] + (frac * (filteredProj[iDetPos+1] - filteredProj[iDetPos])); - } else if (interpType = Backprojector::INTERP_CUBIC) { + } else if (interpType == Backprojector::INTERP_CUBIC) { double p = iDetCenter + curDetPos; // position along detector if (p >= 0 && p < nDet) *pImCol++ += pCubicInterp->interpolate (p); @@ -589,7 +589,7 @@ BackprojectIntDiff::BackprojectView (const double* const filteredProj, const dou if (iDetPos >= 0 && iDetPos <= iLastDet) *pImCol++ += filteredProj[iDetPos] + (detRemainder * deltaFilteredProj[iDetPos]); } // end for iy - } else if (interpType = Backprojector::INTERP_CUBIC) { + } else if (interpType == Backprojector::INTERP_CUBIC) { for (int iy = 0; iy < ny; iy++, curDetPos += det_dy) { *pImCol++ += pCubicInterp->interpolate (static_cast(curDetPos) / 65536); }