Skip to content

nested Patlak and generalized Patlak OSMAPOSL (POSMAPOSL & GPOSMAPOSL) and nested Richardson-Lucy Motion Deconvolution Image Reconstruction (RLMCIR)#1685

Open
nkarakatsanis wants to merge 3445 commits intoUCL:masterfrom
nkarakatsanis:nGPOSMAPOSL_RLMCIR

Conversation

@nkarakatsanis
Copy link

@nkarakatsanis nkarakatsanis commented Mar 2, 2026

Changes in this pull request

This PR introduces nested Patlak and generalized Patlak OSMAPOSL (POSMAPOSL & GPOSMAPOSL) reconstruction.

In addition, nested Richardson-Lucy Motion Deconvolution Image Reconstruction (RLMCIR) is featured, too.

Testing performed

Related issues

Checklist before requesting a review

  • [] I have performed a self-review of my code
  • [] I have added docstrings/doxygen in line with the guidance in the developer guide
  • [] I have implemented unit tests that cover any new or modified functionality (if applicable)
  • [] The code builds and runs on my machine
  • [] documentation/release_XXX.md has been updated with any functionality change (if applicable)

Contribution Notes

Please tick the following:

  • The content of this Pull Request (the Contribution) is intentionally submitted for inclusion in STIR (the Work) under the terms and conditions of the Apache-2.0 License.
  • I (or my institution) have signed the STIR Contribution License Agreement (not required for small changes).

Kris Thielemans added 30 commits February 13, 2012 17:08
- many updates and clarifications, especially related to using CMake
… if all options were given (i.e. if

you attempted to use the defaults, the program crashed). Also, confusingly it printed text with ACF_threshold once it was going, while the command line option is ACF-threshold.
…ite (otherwise it doesn't generate the documentation)
…of using a const member as size for an array (needed for VC2010)
…ing about WIN32 on cygwin is still there.

- add check on system getopt here (not only in buildblock) as stir/getopt.h need to know
…re we called error())

- slightly better error message when reading the matrix failed
…sed to private typedef (as this confuses SWIG). This change makes sense anyway.
…pedef (as this confuses SWIG). This change makes sense anyway.
Kris Thielemans and others added 30 commits January 13, 2014 18:09
- adapt ITKImageInputFileFormat::read_from_file to check if it's a DICOM file, and if so, attempt to read the whole volume (differentiating on series/frame/gate)
The script checks scatter against a previous result. However, there's a random number
generator is used to assign the exact location of the scatter points. With the current
settings, the difference can be larger than the previous threshold, so I increased it to 4% of the max.
- ECAT7 projection data was read incorrectly (wrong order of arguments in a function)
- create_projdata_template sets modality to PET
- handle filesnames with multiple dots
- have sensible voxel-sizes for tangential and axial dimensions
It is highly recommended to use these as opposed to fwdtest and bcktest in your scripts.
Note the different order of the arguments!
- removed a few hard-wired filenames in simulate_data.sh (and use its arguments properly)
- fixed creation of randoms such that the same command would be usable for zero randoms
- moved scripts and *txt files one level up
- moved samples, examples and samples/PET_simulation to new sub-dir examples in the top-level directory
- created CMakeLists.txt in top-level directory (and moved PROJECT and other main statements there from src/CMakeLists.txt).
- moved CMake helper files to src/cmake, adjust many CMakeLists.txt to rely on module path for these files (by no longer specifying location using paths)
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