Search

X-ray phase-contrast imaging enhances soft tissue visualization by leveraging the phase shift of X-rays passing through materials. It permits to minimize radiation exposure due to high contrast, as well as high resolution imaging limited by the wavelength of the X-rays. Phase retrieval extracts the phase shift computationally, but simulated images fail to recreate low-frequency noise observed in experimental images. To this end, we propose a new method to simulate phase contrast images using the Wigner Distribution Function. This permits the simulation of wave and particle effects simultaneously and simulates images photon by photon. Here, we give a first demonstration of the method by simulating the Gaussian double-slit experiment. It has the potential for realistic simulation of low-dose imaging., Comment: 4 pages, 5 figures

Proton computed tomography (CT) is similar to x-ray CT but relies on protons rather than photons to form an image. In its most common operation mode, the measured quantity is the amount of energy that a proton has lost while traversing the imaged object from which a relative stopping power map can be obtained via tomographic reconstruction. To this end, a calorimeter which measures the energy deposited by protons downstream of the scanned object has been studied or implemented as energy detector in several proton CT prototypes. An alternative method is to measure the proton's residual velocity and thus its kinetic energy via the time of flight (TOF) between at least two sensor planes. In this work, we study the precision, i.e. image noise, which can be expected from TOF proton CT systems. We rely on physics models on the one hand and statistical models of the relevant uncertainties on the other to derive closed form expressions for the noise in projection images. The TOF measurement error scales with the distance between the TOF sensor planes and is reported as velocity error in ps/m. We use variance reconstruction to obtain noise maps of a water cylinder phantom given the scanner characteristics and additionally reconstruct noise maps for a calorimeter-based proton CT system as reference. We find that TOF proton CT with 30 ps/m velocity error reaches similar image noise as a calorimeter-based proton CT system with 1% energy error (1 sigma error). A TOF proton CT system with a 50 ps/m velocity error produces slightly less noise than a 2% calorimeter system. Noise in a reconstructed TOF proton CT image is spatially inhomogeneous with a marked increase towards the object periphery. This systematic study of image noise in TOF proton CT can serve as a guide for future developments of this alternative solution for estimating the residual energy of protons after the scanned object.

We are currently conceiving, through (MC) simulation, a multi-channel gamma detector array (TIARA for Time-of-flight Imaging ARrAy) for the online monitoring of protontherapy treatments. By measuring the Time-Of-Flight (TOF) between a beam monitor placed upstream and the Prompt-Gamma (PG) detector, our goal is to reconstruct the PG vertex distribution to detect a possible deviation of proton beam delivery. In this paper, two non-iterative reconstruction strategies are proposed. The first is based on the resolution of an analytical formula describing the PG vertex distribution in 3D. Here, it was resolved under a one-dimensional approximation in order to measure a potential proton range shift along the beam direction. The second is based on the calculation of the Centre-Of-Gravity (COG) of the TIARA pixel detectors counts and also provides 3D information on a possible beam displacement. The PG vertex reconstruction was evaluated in two different scenarios. A coincidence time resolution of 100 ps (rms) can be attained in single proton regime (operating a reduction of the beam current) and using an external beam monitor to provide a start trigger for the TOF measurement. Under these conditions, MC simulations have shown that a millimetric proton range shift sensitivity can be achieved at 2$\sigma$ with 10$^{8}$ incident protons. This level of accuracy would allow to act in real-time if the treatment does not conform to treatment plan. A worst case scenario of a 1 ns (rms) TOF resolution was also considered to demonstrate that a degraded timing information can be compensated by increasing the acquisition statistics: in this case, a 2 mm range shift would be detectable at 2$\sigma$ with 10$^{9}$ incident protons. The COG method has shown excellent capabilities of detecting lateral beam displacements: a 2 mm sensitivity was found at 2$\sigma$ with 10$^{8}$ incident protons.

The accuracy of hadrontherapy treatment is currently limited by ion-range uncertainties. In order to fully exploit the potential of this technique, we propose the development of a novel system for online control of particle therapy, based on TOF-resolved (time-of-flight) Prompt Gamma (PG) imaging with 100 ps time resolution. Our aim is to detect a possible deviation of the proton range with respect to treatment planning within the first few irradiation spots at the beginning of the session. The system consists of a diamond-based beam hodoscope for single proton tagging, operated in time coincidence with one or more gamma detectors placed downstream of the patient. The TOF between the proton time of arrival in the hodoscope and the PG detection time provides an indirect measurement of the proton range in the patient with a precision strictly related to the system time resolution. With a single ~38 cm$^{3}$ BaF2 detector placed at 15 cm from a heterogeneous PMMA target, we obtained a coincidence time resolution of 101 ps (rms). This system allowed us to measure the thickness and position of an air cavity within a PMMA target, and the associated proton range shift: a 3 mm shift can be detected at 2$\sigma$ confidence level within a single large irradiation spot (~10$^{8}$ protons). We are currently conceiving a multi-channel PG timing detector with 3D target coverage. Each pixel will provide the PG detection time and its hit position, that can be used to reconstruct the 3D distribution of PG vertices in the patient. Our approach does not require collimation and allows to dramatically increase the detection efficiency. Since both signal detection and background rejection are based on TOF, the constraints on energy resolution can be relaxed to further improve time resolution., Comment: 4 pages, 6 figures

A method is proposed and evaluated to model large and inconvenient phase space files used in Monte Carlo simulations by a compact Generative Adversarial Network (GAN). The GAN is trained based on a phase space dataset to create a neural network, called Generator (G), allowing G to mimic the multidimensional data distribution of the phase space. At the end of the training process, G is stored with about 0.5 million weights, around 10 MB, instead of few GB of the initial file. Particles are then generated with G to replace the phase space dataset. This concept is applied to beam models from linear accelerators (linacs) and from brachytherapy seed models. Simulations using particles from the reference phase space on one hand and those generated by the GAN on the other hand were compared. 3D distributions of deposited energy obtained from source distributions generated by the GAN were close to the reference ones, with less than 1% of voxel-by-voxel relative difference. Sharp parts such as the brachytherapy emission lines in the energy spectra were not perfectly modeled by the GAN. Detailed statistical properties and limitations of the GAN-generated particles still require further investigation, but the proposed exploratory approach is already promising and paves the way for a wide range of applications., Comment: Phys Med Biol (Accepted Manuscript online 30 August 2019)