UCL / PETPVC

Partial Volume Correction in PET
Apache License 2.0
50 stars 14 forks source link

Richardson-Lucy can create very high voxel values at alternating iterations #81

Closed KrisThielemans closed 3 years ago

KrisThielemans commented 3 years ago

We sometimes see that the max in the image at even (or odd) iterations becomes very large (1E20 or more), while in the next iteration it is fine. This seems due to 1 or 2 voxels. The rest of the image is not affected.

Presumably this is because of numerical problems in the division of RL.