Model-based iterative CT reconstruction suffers from its high computational cost. To increase the reconstruction speed, algorithms using a momentum acceleration can be combined with the ordered subset approach. As shown in previous works, the incorporation of ordered subsets strongly improves the reconstruction time but makes the outcome of the reconstruction difficult to predict due to its empirical nature. Specifically, when the number of subsets is too high, divergence or convergence to an unsatisfactory result can easily appear. In this work, we propose a new combination of ordered subsets with momentum to achieve fast reconstruction in a more robust manner. Our approach, referred to as EOS-MFISTA, is an Efficient Ordered Subset strategy based on the steps of MFISTA, a monotonic version of FISTA. In short, EOS-MFISTA uses objective function values to decide on the next update step and does this with little increase in computational effort. The performance of EOS-MFISTA is evaluated on real CT data of a head phantom and a human chest. Starting from a zero image, 200 iterations of EOS-MFISTA are calculated and compared with other algorithms. Quantitative analyses based on the RMSE show that EOS-MFISTA is much more robust than OS-FISTA, especially when the number of subsets is increased to accelerate reconstruction. Visual inspection of results obtained with 50 iterations further supports these findings.
|