summaryrefslogtreecommitdiffstats
path: root/src/SparseMatrixProjector2D.cpp
diff options
context:
space:
mode:
authorWillem Jan Palenstijn <Willem.Jan.Palenstijn@cwi.nl>2019-01-24 13:48:49 +0100
committerWillem Jan Palenstijn <Willem.Jan.Palenstijn@cwi.nl>2019-01-24 14:18:03 +0100
commite2d5375ecc5c8fc45b796e0bd9d948cd729abcc9 (patch)
tree774f65447d591ce27abbaccc205940a010775c7b /src/SparseMatrixProjector2D.cpp
parentb9563df34df98480d35a53caa2b81d998440f9c5 (diff)
downloadastra-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.cpp41
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;
-}
-
-//----------------------------------------------------------------------------------------