diff options
author | Willem Jan Palenstijn <Willem.Jan.Palenstijn@cwi.nl> | 2019-01-24 13:48:49 +0100 |
---|---|---|
committer | Willem Jan Palenstijn <Willem.Jan.Palenstijn@cwi.nl> | 2019-01-24 14:18:03 +0100 |
commit | e2d5375ecc5c8fc45b796e0bd9d948cd729abcc9 (patch) | |
tree | 774f65447d591ce27abbaccc205940a010775c7b /src/SparseMatrixProjector2D.cpp | |
parent | b9563df34df98480d35a53caa2b81d998440f9c5 (diff) | |
download | astra-e2d5375ecc5c8fc45b796e0bd9d948cd729abcc9.tar.gz astra-e2d5375ecc5c8fc45b796e0bd9d948cd729abcc9.tar.bz2 astra-e2d5375ecc5c8fc45b796e0bd9d948cd729abcc9.tar.xz astra-e2d5375ecc5c8fc45b796e0bd9d948cd729abcc9.zip |
Remove largely unimplemented CProjector2D::projectPoint method
This includes the astra_mex_projector('splat') matlab function.
Diffstat (limited to 'src/SparseMatrixProjector2D.cpp')
-rw-r--r-- | src/SparseMatrixProjector2D.cpp | 41 |
1 files changed, 0 insertions, 41 deletions
diff --git a/src/SparseMatrixProjector2D.cpp b/src/SparseMatrixProjector2D.cpp index 6fded6a..46302f3 100644 --- a/src/SparseMatrixProjector2D.cpp +++ b/src/SparseMatrixProjector2D.cpp @@ -174,44 +174,3 @@ void CSparseMatrixProjector2D::computeSingleRayWeights(int _iProjectionIndex, _iStoredPixelCount = p.getStoredPixelCount(); } -//---------------------------------------------------------------------------------------- -// Splat a single point -std::vector<SDetector2D> CSparseMatrixProjector2D::projectPoint(int _iRow, int _iCol) -{ - unsigned int iVolumeIndex = _iCol * m_pVolumeGeometry->getGridRowCount() + _iRow; - - // NOTE: This is very slow currently because we don't have the - // sparse matrix stored in an appropriate form for this function. - std::vector<SDetector2D> ret; - - const CSparseMatrix* pMatrix = dynamic_cast<CSparseMatrixProjectionGeometry2D*>(m_pProjectionGeometry)->getMatrix(); - - for (int iAngle = 0; iAngle < m_pProjectionGeometry->getProjectionAngleCount(); ++iAngle) - { - for (int iDetector = 0; iDetector < m_pProjectionGeometry->getDetectorCount(); ++iDetector) - { - int iRayIndex = iAngle * m_pProjectionGeometry->getDetectorCount() + iDetector; - const unsigned int* piColIndices; - const float32* pfValues; - unsigned int iSize; - - pMatrix->getRowData(iRayIndex, iSize, pfValues, piColIndices); - - for (unsigned int i = 0; i < iSize; ++i) { - if (piColIndices[i] == iVolumeIndex) { - SDetector2D s; - s.m_iIndex = iRayIndex; - s.m_iAngleIndex = iAngle; - s.m_iDetectorIndex = iDetector; - ret.push_back(s); - break; - } else if (piColIndices[i] > iVolumeIndex) { - break; - } - } - } - } - return ret; -} - -//---------------------------------------------------------------------------------------- |