Real-time CBCT reconstructions using Krylov solvers in repeated scanning procedures

This work introduces a new efficient iterative solver for the reconstruction of real-time cone-beam computed tomography, which is based on the prior image constrained compressed sensing (PICCS) regularization and leverages the efficiency of Krylov subspace methods. In particular, we focus on the setting where a sequence of under-sampled computed tomography (CT) scans are taken on the same object with only local changes (e.g. changes in a tumour size or the introduction of a surgical tool). This is very common, for example, in image-guided surgery, where the amount of measurements is limited to ensure the safety of the patient. In this case, we can also typically assume that a (good) initial reconstruction for the solution exists, coming from a previously over-sampled scan, so we can use this information to aid the subsequent reconstructions. The effectiveness of this method is demonstrated in both a synthetic scan and using real CT data, where it can be observed that the PICCS framework is very effective for the reduction of artefacts, and that the new method is faster than other common alternatives used in the same setting.