Dr Lucia Florescu

Broken ray transforms (BRTs) are typically considered to be reciprocal, meaning that the transform is independent of the direction in which a photon travels along a given broken ray. However, if the photon can change its energy (or be absorbed and re-radiated at a different frequency) at the vertex of the ray, then reciprocity is lost. In optics, non-reciprocal BRTs are applicable to imaging problems with fluorescent contrast agents. In the case of x-ray imaging, problems with single Compton scattering also give rise to non-reciprocal BRTs. In this paper, we focus on tomographic optical fluorescence imaging and show that, by reversing the path of a photon and using the non-reciprocity of the data function, we can reconstruct simultaneously and independently all optical properties of the medium (the intrinsic attenuation coefficients at the excitation and the fluorescence frequency and the concentration of the contrast agent). Our results are also applicable to inverting BRTs that arise due to single Compton scattering.

Numerical experiments were performed to analyse the effect of data loss at the edges of the sample on the accuracy of optical-CT reconstruction, in the context of applications to radiation dosimetry.

We address the interior problem of computed tomography that occurs when projection data is only available for a region in the interior of the sample. In this case, it is not possible to accurately reconstruct the attenuation function even in the interior domain. We consider an algorithm for correcting the interior tomography reconstruction which is based on prior knowledge in the interior domain. This correction algorithm is evaluated by performing numerical experiments with the Shepp-Logan phantom for various amounts of data loss, noise in the available projection data, various values of the attenuation function known a priori, and various positions within the sample where the prior information is available. Good performance of the algorithm based on prior knowledge at one point is demonstrated in the case of noiseless data. In the presence of noise in the projection data, improvements in the reconstructed attenuation function are obtained based on prior knowledge at a number of points in the interior domain. The robustness of the correction algorithm to errors in the values of the attenuation function used as prior knowledge was also investigated.

We present a tomographic imaging technique based on angularly-selective measurements of fluorescent light that enables for the first time simultaneous reconstruction of the attenuation coefficient at two energies and of the contrast-agent concentration.

Optical methods of biomedical tomographic imaging are of considerable interest due to their non-invasive nature and sensitivity to physiologically important markers. Similarly to other imaging modalities, optical methods can be enhanced by utilizing extrinsic contrast agents. Typically, these are fluorescent molecules, which can aggregate in regions of interest due to various mechanisms. In the current approaches to imaging, the intrinsic (related to the tissue) and extrinsic (related to the contrast agent) optical parameters are determined separately. This can result in errors, in particular, due to using simplified heuristic models for the spectral dependence of the optical parameters. Recently, we have developed the theory of non-reciprocal broken-ray tomography (NRBRT) for fluorescence imaging of weakly scattering systems. NRBRT enables simultaneous reconstruction of the fluorophore concentration as well as of the intrinsic optical attenuation coefficient at both the excitation and the emission wavelengths. Importantly, no assumption about the spectral dependence of the tissue optical properties is made in NRBRT. In this study, we perform numerical validation of NRBRT under realistic conditions using the Monte Carlo method to generate forward data. We demonstrate that NRBRT can be used for tomographic imaging of samples of up to four scattering lengths in size. The effects of physical characteristics of the detectors such as the area and the acceptance angle are also investigated.

We perform numerical experiments based on Monte Carlo simulations and clinical CT data to investigate Cherenkov light emission in molecular radiation therapy of hyperthyroidism, and demonstrate that Cherenkov light-based dosimetry could be feasible.

Numerical experiments based on Monte Carlo simulations and clinical CT data are performed to investigate the spatial and spectral characteristics of Cherenkov light emission and the relationship between Cherenkov light intensity and deposited dose in molecular radiotherapy of hyperthyroidism and papillary thyroid carcinoma. It is found that Cherenkov light is emitted mostly in the treatment volume, the spatial distribution of Cherenkov light at the surface of the patient presents high-value regions at locations that depend on the symmetry and location of the treatment volume, and the surface light in the near-infrared spectral region originates from the treatment site. The effect of inter-patient variability in the tissue optical parameters and radioisotope uptake on the linear relationship between the dose absorbed by the treatment volume and Cherenkov light intensity at the surface of the patient is investigated, and measurements of surface light intensity for which this effect is minimal are identified. The use of Cherenkov light measurements at the patient surface for molecular radiation therapy dosimetry is also addressed.