You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardexpand all lines: docs/source/algorithms.rst
+9-9
Original file line number
Diff line number
Diff line change
@@ -122,24 +122,24 @@ MLEM/OSEM
122
122
While only OSEM is selectable, MLEM is enabled if no subsets are used. This method can be used for PET, SPECT or CT data, or any other Poisson-based data. Note that CT uses its own transmission tomography based formula, while
123
123
PET and SPECT use the linear model. Useful algorithm for PET and SPECT, but not particularly recommended for CT. Use OSL_OSEM for regularized version (see below).
124
124
125
-
|EM MLEM based on: https://doi.org/10.1111/j.2517-6161.1977.tb01600.x
126
-
|EM OSEM: https://doi.org/10.1109/42.363108
125
+
|Emission tomography (ET) MLEM based on: https://doi.org/10.1111/j.2517-6161.1977.tb01600.x
126
+
|ET OSEM: https://doi.org/10.1109/42.363108
127
127
128
128
RAMLA
129
129
^^^^^
130
130
131
131
Similar to OSEM, but has guaranteed convergence and is dependent on the relaxation parameter lambda (or lambdaN in Python), see RELAXATION PARAMETER. Slower to converge than OSEM. Can be used with or without subsets.
132
132
Note that the default lambda values might not work with RAMLA. Not recommended for CT but has transmission tomography based version implemented. See BSREM for regularized version.
133
133
134
-
EM version based on: https://doi.org/10.1109/42.538946
134
+
ET version based on: https://doi.org/10.1109/42.538946
135
135
136
136
MRAMLA
137
137
^^^^^^
138
138
139
-
Unregularized version of the MBSREM. Almost identical to RAMLA, i.e. requires lambda, but supports preconditioners. EM preconditioner is also highly recommended! Has some additional steps to guarantee convergence.
139
+
Unregularized version of the MBSREM. Almost identical to RAMLA, i.e. requires lambda, but supports preconditioners. ET preconditioner is also highly recommended! Has some additional steps to guarantee convergence.
140
140
Also has dedicated transmission tomography version. Useful for any Poisson-based data, if regularization is not used.
141
141
142
-
EM version based on: https://doi.org/10.1109/TMI.2003.812251
142
+
ET version based on: https://doi.org/10.1109/TMI.2003.812251
143
143
144
144
ROSEM
145
145
^^^^^
@@ -229,18 +229,18 @@ OSL based on: https://doi.org/10.1109/42.52985
229
229
MBSREM
230
230
^^^^^^
231
231
232
-
Regularized version of MRAMLA. Requires relaxation parameter lambda, and supports preconditioners. EM preconditioner is also highly recommended! Has some additional steps to guarantee convergence.
232
+
Regularized version of MRAMLA. Requires relaxation parameter lambda, and supports preconditioners. ET preconditioner is also highly recommended! Has some additional steps to guarantee convergence.
233
233
Also has dedicated transmission tomography version. Useful for any Poisson-based data, if regularization is used.
234
234
235
-
EM version based on: https://doi.org/10.1109/TMI.2003.812251
235
+
ET version based on: https://doi.org/10.1109/TMI.2003.812251
236
236
237
237
BSREM
238
238
^^^^^
239
239
240
240
Regularized version of RAMLA. However, unlike MBSREM, BSREM handles the regularization differently. While MBSREM computes the regularization after every subset, BSREM does it only after one full iteration (epoch). This can
241
241
sometimes be useful as less regularization steps might be used. Requires relaxation parameter lambda. Also has dedicated transmission tomography version.
242
242
243
-
EM version based on: https://doi.org/10.1109/42.921477
243
+
ET version based on: https://doi.org/10.1109/42.921477
244
244
245
245
ROSEM-MAP
246
246
^^^^^^^^^
@@ -268,7 +268,7 @@ Note that for PET and SPECT data the relaxation parameter can safely begin at 1,
268
268
until it is of the right magnitude. Too high values will cause quick divergence while too low values will cause slow convergence. For CT, the default value of 1 is divided by 10000 (when you use default values). This should work for
269
269
most CT applications, but it might not be optimal. To fix this, a proper normalization would be required for the backprojection (suggestions are welcome!).
270
270
271
-
EM version based on: https://doi.org/10.1109/TMI.2019.2898271
271
+
ET version based on: https://doi.org/10.1109/TMI.2019.2898271
0 commit comments