Skip to content

Conversation

@saarlemo
Copy link
Contributor

@saarlemo saarlemo commented Mar 27, 2025

Orthogonal distance ray tracing projector for SPECT

  • Assumes circular holes for now
  • A voxel is weighed by the parallel and orthogonal distance from the detector element center normal
  • The weight is calculated from a Gaussian evaluated at the orthogonal distance, with variance dependent on the parallel distance
  • CPU, OpenCL and CUDA

SPECT reconstructions with matrix operations ($Ax$ forward projection and $A^Ty$ backprojection)

  • MATLAB and Python
  • The examples include a single pixel backprojection and single voxel forward projection

Rotation-based projector FOV resolution is no longer constrained by sinogram resolution

  • Sinograms are resized and resampled to match FOV XZ-plane
  • Cubical voxels still necessary

The function volume3Dviewer now has support for PET/CT and SPECT/CT multimodal overlays (MATLAB only as Octave has no support for the alpha function).

Miscellaneous improvements

  • Fixed typos
  • Macros for OpenCL / CUDA status checks
  • Reformatted the wall of bullet points in readme

saarlemo added 30 commits March 4, 2025 17:21
@villekf villekf merged commit 09c9254 into villekf:master May 13, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants