Remove stripes Titarenko

Remove stripes Titarenko#

Description

This stripe removal algorithm is developed by Dr. Sofya Titarenko and is a part of the module httomolibgpu.prep.stripe for data correction. This filter works in the sinogram space to minimise stripes leading to Ring Artifacts in the reconstruction. See more about the filter in the original publication [9]. The algorithm uses a priori information about the sinogram, such as its smoothness and symmetry.

Where and how to use it:

Stripe removal method is used for data pre-processing, i.e., before the reconstruction to remove the stripes present in the data. If the stripes are not removed, this will lead to ring artefacts in the reconstruction and further problems with quantification and segmentation.

What are the adjustable parameters:

The important parameter here is beta which is responsible for strength of the filter. Note that lower values increase the filter strength.

Practical example:

In this example we demonstrate how to apply Titarenko’s stripe removal algorithm to the normalised projection data.

Warning

As with any pre-processing tool, one needs to be careful in applying stripe removal algorithm to data. Sub-optimal parameters can potentially lead to other artifacts introduced in the reconstruction.

../../_images/data_stripes_added_sino.png

Fig. 31 Input (sinogram view) with the stripes added#

../../_images/data_stripes_added_proj.png

Fig. 32 Input (projection view) with the stripes added#

../../_images/remove_stripe_ti_sino.png

Fig. 33 After applying remove stripes by Titarenko’s stripe removal (sinogram)#

../../_images/remove_stripe_ti_proj.png

Fig. 34 After applying remove stripes by Titarenko’s stripe removal (projection)#

../../_images/remove_stripe_ti_res_sino.png

Fig. 35 Sinogram view of absolute residual between input and output#

../../_images/remove_stripe_ti_res_proj.png

Fig. 36 Projection view of absolute residual between input and output#