Remove stripes combo

Remove stripes combo#

Description

This stripe removal algorithm is developed by Dr. Nghia Vo 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 [12]. Also more information about the method can be found on the author’s software Sarepy.

Remove stripes combination uses four algorithms applied sequentially to the data in order to remove different types of stripe artefacts as presented in the paper [12].

Where and how to use it:

Remove stripes combination should be applied when different types of stripe artefacts present in the data (e.g. small or large or partial stripes). It is useful to have a look at this classification of different ring/stripe artefacts beforehand. The method can be used blindly, however, applying multiple pre-processing algorithms can distort the data, resulting in artefacts in the reconstruction. Meaning that if there are not to many stripes in the data and the method is still applied, could result in suboptimal reconstructions.

What are the adjustable parameters:

As this method based on multiple algorithms, including Remove stripes by sorting, it has a number of free parameters. The snr parameter is used to detect and remove large stripes. The associated to that method parameter la_size is used to remove large stripes, also by means of sorting. Parameter sm_size defines the size of the median filter when sorting is performed to remove thinner stripes.

Practical example:

In this example we demonstrate how to apply remove stripes combo 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. 25 Input (sinogram view) with the stripes added#

../../_images/data_stripes_added_proj.png

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

../../_images/remove_all_stripe_sino.png

Fig. 27 After applying remove stripes combo (sinogram)#

../../_images/remove_all_stripe_proj.png

Fig. 28 After applying remove stripes combo (projection)#

../../_images/remove_all_stripe_res_sino.png

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

../../_images/remove_all_stripe_res_proj.png

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