[Rtk-users] RTK Scatter Correction

Lotte Schyns lotte.schyns at maastro.nl
Fri Feb 23 11:37:25 EST 2018


Thank you Simon. I see that the implementation of the Boellaard method is indeed not as complete as in the paper. I will have a go with the other methods. Especially the MC one seems interesting.

On 21-02-18 08:50, Simon Rit wrote:
Dear Lotte,
These two options address different issues. The first one deals with the patient scatter and the second one with the detector glare which might come from the detector scatter (but maybe not (only)). See the work of, e.g., Poludniowski to understand the difference
dx.doi.org/10.1088/0031-9155/54/12/016<http://dx.doi.org/10.1088/0031-9155/54/12/016>
dx.doi.org/10.1088/0031-9155/56/6/019<http://dx.doi.org/10.1088/0031-9155/56/6/019>
You'll see that the second paper suggests that both must be corrected.
I haven't used those implementations a lot but I think the detector glare follows closely Poludniowski's paper but Boellaard does not. Boellaard implementation is the simplest option when only a constant is subtracted and I actually believe that it's not working properly (but this would need to be checked, it's been on my todo list for a long while).
There are other options:
- auto-detection of I0 in air does something similar to Boellaard's but differently
http://www.openrtk.org/Doxygen/classrtk_1_1I0EstimationProjectionFilter.html
- empirical cupping correction:
https://doi.org/10.1118/1.2188076
- Monte Carlo correction if you have a CT image of your object. The tools to do this are open source, see
https://doi.org/10.1016/j.phro.2017.09.002
I hope this helps,
Simon

On Tue, Feb 20, 2018 at 3:14 PM, Lotte Schyns <lotte.schyns at maastro.nl<mailto:lotte.schyns at maastro.nl>> wrote:
Hello,

I'm investigating the possibility to perform scatter corrections in RTK.
From what I could find, there seem to be two options  (correct me if I'm
wrong):

1) rtk::BoellaardScatterCorrectionImageFilter (Boellaard paper)
2) rtk::ScatterGlareCorrectionImageFilter (Poludniowski paper)

Since both methods are based on a deconvolution approach using the
edge-spread function, I was wondering what the difference in
implementation is and in which cases one method would be preferred over
the other.

Lotte
_______________________________________________
Rtk-users mailing list
Rtk-users at public.kitware.com<mailto:Rtk-users at public.kitware.com>
https://public.kitware.com/mailman/listinfo/rtk-users


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://public.kitware.com/pipermail/rtk-users/attachments/20180223/17fd31e8/attachment-0001.html>


More information about the Rtk-users mailing list