Sliding motion is a challenge for deformable image registration because it leads to discontinuities in the sought deformation. In this paper, we present a method to handle sliding motion using ...multiple B-spline transforms. The proposed method decomposes the sought deformation into sliding regions to allow discontinuities at their interfaces, but prevents unrealistic solutions by forcing those interfaces to match. The method was evaluated on 16 lung cancer patients against a single B-spline transform approach and a multi B-spline transforms approach without the sliding constraint at the interface. The target registration error (TRE) was significantly lower with the proposed method (TRE = 1.5 mm) than with the single B-spline approach (TRE = 3.7 mm) and was comparable to the multi B-spline approach without the sliding constraint (TRE = 1.4 mm). The proposed method was also more accurate along region interfaces, with 37% less gaps and overlaps when compared to the multi B-spline transforms without the sliding constraint.
This paper studies the impact of tiny changes in region-of-interest (ROI) tomography system matrices on the variance of the reconstructed ROI. In small-scale and medium-scale examples, the variance ...in the reconstructed ROI was estimated for different system matrices. The results revealed a striking and counterintuitive phenomenon: a tiny change in the system matrix can dramatically affect the variance of the ROI estimate. In one of our examples, a decrease of 0.1% in one element out of hundreds of thousands of the system matrix resulted in a systematic reduction of the variance inside the ROI, and by a factor of 5 to 10 for some pixels. Our results agree with a recently proven theorem about the ability of additional measurements to reduce the variance in ROI tomography.
Proton transmission imaging has been proposed and investigated as imaging modality complementary to x-ray based techniques in proton beam therapy. In particular, it addresses the issue of range ...uncertainties due to the conversion of an x-ray patient computed tomography (CT) image expressed in Hounsfield Units (HU) to relative stopping power (RSP) needed as input to the treatment planning system. One approach to exploit a single proton radiographic projection is to perform a patient-specific calibration of the CT to RSP conversion curve by optimising the match between a measured and a numerically integrated proton radiography. In this work, we develop the mathematical tools needed to perform such an optimisation in an efficient and robust way. Our main focus lies on set-ups which combine pencil beam scanning with a range telescope detector, although most of our methods can be employed in combination with other set-ups as well. Proton radiographies are simulated in Monte Carlo using an idealised detector and applying the same data processing chain used with experimental data. This approach allows us to have a ground truth CT-RSP curve to compare the optimisation results with. Our results show that the parameters of the CT-RSP curve are strongly correlated when using a pencil beam based set-up, which leads to unrealistic variation in the optimised CT-RSP curves. To address this issue, we introduce a regularisation procedure which guarantees a plausible degree of smoothness in the optimised CT-RSP curves. We investigate three different methods to perform the numerical projection operation needed to generate a proton digitally reconstructed radiography. We find that the approximate and computationally faster method performs as well as the more accurate but more demanding method. We perform a Monte Carlo experiment based on a head and neck patient to evaluate the range accuracy achievable with the optimised CT-RSP curves and find an agreement with the ground truth expectation of better than Formula: see text. Our results further indicate that the region in the patient in which the proton radiography is acquired does not necessarily have to correspond to the treatment volume to achieve this accuracy. This is important as the imaged region could be freely chosen, e.g. in order to spare organs at risk.
Proton computed tomography (pCT) has high accuracy and dose efficiency in producing spatial maps of the relative stopping power (RSP) required for treatment planning in proton therapy. With ...fluence-modulated pCT (FMpCT), prescribed noise distributions can be achieved, which allows to decrease imaging dose by employing object-specific dynamically modulated fluence during the acquisition. For FMpCT acquisitions we divide the image into region-of-interest (ROI) and non-ROI volumes. In proton therapy, the ROI volume would encompass all treatment beams. An optimization algorithm then calculates dynamically modulated fluence that achieves low prescribed noise inside the ROI and high prescribed noise elsewhere. It also produces a planned noise distribution, which is the expected noise map for that fluence, as calculated with a Monte Carlo simulation. The optimized fluence can be achieved by acquiring pCT images with grids of intensity modulated pencil beams. In this work, we interfaced the control system of a clinical proton beam line to deliver the optimized fluence. Using three phantoms we acquired images with uniform fluence, with a constant noise prescription, and with an FMpCT task. Image noise distributions as well as fluence maps were compared to the corresponding planned distributions as well as to the prescription. Furthermore, we propose a correction method that removes image artifacts stemming from the acquisition with pencil beams having a spatially varying energy distribution that is not seen in clinical operation. RSP accuracy of FMpCT scans was compared to uniform scans and was found to be comparable to standard pCT scans. While we identified technical improvements for future experimental acquisitions, in particular related to an unexpected pencil beam size reduction and a misalignment of the fluence pattern, agreement with the planned noise was satisfactory and we conclude that FMpCT optimized for specific image noise prescriptions is experimentally feasible.
Protons undergo many small angle deflections when traversing a medium, such as a patient. This effect, known as multiple Coulomb scattering (MCS), leads to degraded image resolution in proton ...radiography and computed tomography (CT) and to lateral spreading of the dose distribution in proton therapy. To optimally account for MCS in proton imaging, the most likely path (MLP) of a proton is estimated based on its position and propagation angle measured in front of and behind the object. In this work, we propose a functional which quantifies the likelihood of a proton trajectory and study how it can be used to model proton trajectories in a homogeneous medium. We focus on two aspects: first, we present an analytical method to quickly generate proton trajectories in a homogeneous medium based on the likelihood functional and validate it through Monte Carlo simulations. It could be used for fast generation of proton CT images without a full Monte Carlo simulation, or potentially to complement the components in a treatment planning Monte Carlo which simulate MCS. Second, by maximising the likelihood functional, we derive an expression for the MLP which is equivalent to the conventional ones reported in the literature yet computationally more convenient. Moreover, we show that the MLP is strictly a polynomial function if the protons' energy loss in the medium is approximated as a polynomial and that the orders of both are linked. We validate our MLP through Monte Carlo simulations and compare proton CT images reconstructed with our expression and with the conventional one. We find that an MLP polynomial of orders larger than five do not lead to increased spatial resolution compared to lower order expressions.
We propose the Reconstruction Toolkit (RTK, http://www.openrtk.org), an open-source toolkit for fast cone-beam CT reconstruction, based on the Insight Toolkit (ITK) and using GPU code extracted from ...Plastimatch. RTK is developed by an open consortium (see affiliations) under the non-contaminating Apache 2.0 license. The quality of the platform is daily checked with regression tests in partnership with Kitware, the company supporting ITK. Several features are already available: Elekta, Varian and IBA inputs, multi-threaded Feldkamp-David-Kress reconstruction on CPU and GPU, Parker short scan weighting, multi-threaded CPU and GPU forward projectors, etc. Each feature is either accessible through command line tools or C++ classes that can be included in independent software. A MIDAS community has been opened to share CatPhan datasets of several vendors (Elekta, Varian and IBA). RTK will be used in the upcoming cone-beam CT scanner developed by IBA for proton therapy rooms. Many features are under development: new input format support, iterative reconstruction, hybrid Monte Carlo / deterministic CBCT simulation, etc. RTK has been built to freely share tomographic reconstruction developments between researchers and is open for new contributions.
Scattering proton CT Krah, N; Quiñones, C T; Létang, J M ...
Physics in medicine & biology,
11/2020, Letnik:
65, Številka:
22
Journal Article
Recenzirano
Odprti dostop
Proton computed tomography (CT) is an imaging modality investigated mainly in the context of proton therapy as a complement to x-ray CT. It uses protons with high enough energy to fully traverse the ...imaged object. Common prototype systems measure each proton's position and direction upstream and downstream of the object as well as the energy loss which can be converted into the water equivalent thickness. A reconstruction algorithm then produces a map of the relative stopping power in the object. As an alternative to energy-loss proton CT, it has been proposed to reconstruct a map of the object's scattering power based on the protons' angular dispersion which can be estimated from the measured directions. As in energy-loss proton CT, reconstruction should best be performed considering the non-linear shape of proton trajectories due to multiple Coulomb scattering (MCS), but no algorithm to achieve this is so far available in the literature. In this work, we propose a filtered backprojection algorithm with distance-driven binning to account for the protons' most likely path. Furthermore, we present a systematic study of scattering proton CT in terms of inherent noise and spatial resolution and study the artefacts which arise from the physics of MCS. Our analysis is partly based on analytical models and partly on Monte Carlo simulations. Our results show that the proposed algorithm performs well in reconstructing relative scattering power maps, i.e. scattering power relative to that of water. Spatial resolution is improved by almost a factor of three compared to straight line projection and is comparable to energy-loss proton CT. Image noise, on the other hand, is inherently much higher. For example, in a water cylinder of 20 cm diameter, representative of a human head, noise in the central image pixel is about 40 times higher in scattering proton CT than in energy-loss proton CT. Relative scattering power in dense regions such as bone inserts is systematically underestimated by a few percent, depending on beam energy and phantom geometry.
Proton computed tomography (CT) has been described as a solution for imaging the proton stopping power of patient tissues, therefore reducing the uncertainty of the conversion of x-ray CT images to ...relative stopping power (RSP) maps and its associated margins. This study aimed to investigate this assertion under the assumption of ideal detection systems. We have developed a Monte Carlo framework to assess proton CT performances for the main steps of a proton therapy treatment planning, i.e. proton or x-ray CT imaging, conversion to RSP maps based on the calibration of a tissue phantom, and proton dose simulations. Irradiations of a computational phantom with pencil beams were simulated on various anatomical sites and the proton range was assessed on the reference, the proton CT-based and the x-ray CT-based material maps. Errors on the tissue's RSP reconstructed from proton CT were found to be significantly smaller and less dependent on the tissue distribution. The imaging dose was also found to be much more uniform and conformal to the primary beam. The mean absolute deviation for range calculations based on x-ray CT varies from 0.18 to 2.01 mm depending on the localization, while it is smaller than 0.1 mm for proton CT. Under the assumption of a perfect detection system, proton range predictions based on proton CT are therefore both more accurate and more uniform than those based on x-ray CT.
The use of ion computed tomography (CT) promises to yield improved relative stopping power (RSP) estimation as input to particle therapy treatment planning. Recently, proton CT (pCT) has been shown ...to yield RSP accuracy on par with state-of-the-art x-ray dual energy CT. There are however concerns that the lower spatial resolution of pCT compared to x-ray CT may limit its potential, which has spurred interest in the use of helium ion CT (HeCT). The goal of this study was to investigate image quality of pCT and HeCT in terms of noise, spatial resolution, RSP accuracy and imaging dose using a detailed Monte Carlo (MC) model of an existing ion CT prototype.
Three phantoms were used in simulated pCT and HeCT scans allowing estimation of noise, spatial resolution and the scoring of dose. An additional phantom was used to evaluate RSP accuracy. The imaging dose required to achieve the same image noise in a water and a head phantom was estimated at both native spatial resolution, and in a scenario where the HeCT spatial resolution was reduced and matched to that of pCT using Hann windowing of the reconstruction filter. A variance reconstruction formalism was adapted to account for Hann windowing.
We confirmed that the scanner prototype would produce higher spatial resolution for HeCT than pCT by a factor 1.8 (0.86 lp mm
versus 0.48 lp mm
at the center of a 20 cm water phantom). At native resolution, HeCT required a factor 2.9 more dose than pCT to achieve the same noise, while at matched resolution, HeCT required only 38% of the pCT dose. Finally, RSP mean absolute percent error (MAPE) was found to be 0.59% for pCT and 0.67% for HeCT.
This work compared the imaging performance of pCT and HeCT when using an existing scanner prototype, with the spatial resolution advantage of HeCT coming at the cost of increased dose. When matching spatial resolution via Hann windowing, HeCT had a substantial dose advantage. Both modalities provided state-of-the-art RSP MAPE. HeCT might therefore help reduce the dose exposure of patients with comparable image noise to pCT, enhanced spatial resolution and acceptable RSP accuracy at the same time.