Performance evaluation of MACACO: a multilayer Compton camera

Gamma-ray sources imaging and test-beam results with MACACO III Compton camera

Hadron therapy is a radiotherapy modality which offers a precise energy deposition to the tumors and a dose reduction to healthy tissue as compared to conventional methods. However, methods for real-time monitoring are required to ensure that the radiation dose is deposited on the target. The IRIS group of IFIC-Valencia developed a Compton camera prototype for this purpose, intending to image the Prompt Gammas emitted by the tissue during irradiation. The system detectors are composed of Lanthanum (III) bromide scintillator crystals coupled to silicon photomultipliers. After an initial characterization in the laboratory, in order to assess the system capabilities for future experiments in proton therapy centers, different tests were carried out in two facilities: PARTREC (Groningen, The Netherlands) and the CNA cyclotron (Sevilla, Spain). Characterization studies performed at PARTREC indicated that the detectors linearity was improved with respect to the previous version and an energy resolution of 5.2 % FWHM at 511 keV was achieved. Moreover, the imaging capabilities of the system were evaluated with a line source of 68Ge and a point-like source of 241Am-9Be. Images at 4.439 MeV were obtained from irradiation of a graphite target with an 18 MeV proton beam at CNA, to perform a study of the system potential to detect shifts at different intensities. In this sense, the system was able to distinguish 1 mm variations in the target position at different beam current intensities for measurement times of 1800 and 600 s.

A Portable Three-Layer Compton Camera for Wide-Energy-Range Gamma-ray Imaging: Design, Simulation and Preliminary Testing

(1) Background: The imaging energy range of a typical Compton camera is limited due to the fact that scattered gamma photons are seldom fully absorbed when the incident energies are above 3 MeV. Further improving the upper energy limit of gamma-ray imaging has important application significance in the active interrogation of special nuclear materials and chemical warfare agents, as well as range verification of proton therapy. (2) Methods: To realize gamma-ray imaging in a wide energy range of 0.3~7 MeV, a principle prototype, named a portable three-layer Compton camera, is developed using the scintillation detector that consists of an silicon photomultiplier array coupled with a Gd3Al2Ga3O12:Ce pixelated scintillator array. Implemented in a list-mode maximum likelihood expectation maximization algorithm, a far-field energy-domain imaging method based on the two interaction events is applied to estimate the initial energy and spatial distribution of gamma-ray sources. The simulation model of the detectors is established based on the Monte Carlo simulation toolkit Geant4. The reconstructed images of a 133Ba, a 137Cs and a 60Co point-like sources have been successfully obtained with our prototype in laboratory tests and compared with simulation studies. (3) Results: The proportion of effective imaging events accounts for about 2%, which allows our prototype to realize the reconstruction of the distribution of a 0.05 μSv/h 137Cs source in 10 s. The angular resolution for resolving two 137Cs point-like sources is 15°. Additional simulated imaging of the 6.13 MeV gamma-rays from 14.1 MeV neutron scattering with water preliminarily demonstrates the imaging capability for high incident energy. (4) Conclusions: We conclude that the prototype has a good imaging performance in a wide energy range (0.3~7 MeV), which shows potential in several MeV gamma-ray imaging applications.

Study of Alternative Imaging Methods for In Vivo Boron Neutron Capture Therapy

Boron Neutron Capture Therapy (BNCT) is an innovative and highly selective treatment against cancer. Nowadays, in vivo boron dosimetry is an important method to carry out such therapy in clinical environments. In this work, different imaging methods were tested for dosimetry and tumor monitoring in BNCT based on a Compton camera detector. A dedicated dataset was generated through Monte Carlo tools to study the imaging capabilities. We first applied the Maximum Likelihood Expectation Maximization (MLEM) iterative method to study dosimetry tomography. As well, two methods based on morphological filtering and deep learning techniques with Convolutional Neural Networks (CNN), respectively, were studied for tumor monitoring. Furthermore, clinical aspects such as the dependence on the boron concentration ratio in image reconstruction and the stretching effect along the detector position axis were analyzed. A simulated spherical gamma source was studied in several conditions (different detector distances and boron concentration ratios) using MLEM. This approach proved the possibility of monitoring the boron dose. Tumor monitoring using the CNN method shows promising results that could be enhanced by increasing the training dataset.

Influence of the background in Compton camera images for proton therapy treatment monitoring

Objective. Background events are one of the most relevant contributions to image degradation in Compton camera imaging for hadron therapy treatment monitoring. A study of the background and its contribution to image degradation is important to define future strategies to reduce the background in the system.Approach. In this simulation study, the percentage of different kinds of events and their contribution to the reconstructed image in a two-layer Compton camera have been evaluated. To this end, GATE v8.2 simulations of a proton beam impinging on a PMMA phantom have been carried out, for different proton beam energies and at different beam intensities.Main results. For a simulated Compton camera made of Lanthanum (III) Bromide monolithic crystals, coincidences caused by neutrons arriving from the phantom are the most common type of background produced by secondary radiations in the Compton camera, causing between 13% and 33% of the detected coincidences, depending on the beam energy. Results also show that random coincidences are a significant cause of image degradation at high beam intensities, and their influence in the reconstructed images is studied for values of the time coincidence windows from 500 ps to 100 ns.Significance. Results indicate the timing capabilities required to retrieve the fall-off position with good precision. Still, the noise observed in the image when no randoms are considered make us consider further background rejection methods.

Hybrid PET/Compton-camera imaging: an imager for the next generation

Compton cameras can offer advantages over gamma cameras for some applications, since they are well suited for multitracer imaging and for imaging high-energy radiotracers, such as those employed in radionuclide therapy. While in conventional clinical settings state-of-the-art Compton cameras cannot compete with well-established methods such as PET and SPECT, there are specific scenarios in which they can constitute an advantageous alternative. The combination of PET and Compton imaging can benefit from the improved resolution and sensitivity of current PET technology and, at the same time, overcome PET limitations in the use of multiple radiotracers. Such a system can provide simultaneous assessment of different radiotracers under identical conditions and reduce errors associated with physical factors that can change between acquisitions. Advances are being made both in instrumentation developments combining PET and Compton cameras for multimodal or three-gamma imaging systems, and in image reconstruction, addressing the challenges imposed by the combination of the two modalities or the new techniques. This review article summarizes the advances made in Compton cameras for medical imaging and their combination with PET.

Development of a two-layer dense-pixel LYSO Compton camera prototype for prompt gamma imaging

Objective.Lutetium-yttrium orthosilicate (LYSO)-based Compton camera (CC) has been proposed for prompt gamma imaging due to its high detection efficiency and position resolution. However, very few LYSO CC prototypes have been built and used for practical evaluation. In this study, we built a lightweight dense-pixel silicon photomultiplier-based two-layer LYSO CC prototype for future prompt gamma imaging.Approach.We attempt the first-ever effort to use the double-encoding with the thick light guide and coding circuit structure for 46 × 46 dense-pixel LYSO detectors construction and use pixel segmentation based on centroid mapping to obtain 4232 spectral calibrations. We also present a framework for list-mode projection data acquisition based on the decoding of the time series data obtained by data acquisition card in this study. Finally, the standard source calibration, ring-like²²Na source with non-uniform intensity, and mixed point-like source with a wide energy spectrum experiments were implemented to evaluate the resolution metrics and imaging performance of the prototype.Main results.The lateral position resolution of the prototype was 1 mm, and the maximum measurement deviation is 2.5 mm and 5 mm in the depth direction for the scatterer and absorber, respectively. In the experiments, the measured energy resolution was 9.63% @ 1.33 MeV for the scatterer and 10.8% @ 1.33 MeV for the absorber. And the detection efficiency of the prototype for a spherical⁶⁰Co source with a diameter of 2.8 mm at 10 cm far was 5.7 × 10^-3@ 1.33 MeV and the full width at half maximum of the reconstruction was 5.5 mm. Besides, the spatial position offset within 2 mm of the radioactive source at 10 cm can be distinguished.Signification.The developed two-layer dense-pixel LYSO CC contributes to incorporating Compton imaging techniques for prompt gamma detection and multiple energy sources into nuclear medical imaging.

A double photon coincidence detection method for medical gamma-ray imaging

Cascade nuclides emit two or more gamma rays successively through an intermediate state. The coincidence detection of cascade gamma rays provides several advantages in gamma-ray imaging. In this review article, three applications of the double photon coincidence method are reviewed. Double-photon emission imaging with mechanical collimators and Compton double-photon emission imaging can identify radioactive source positions with their angular-resolving detectors, and reduce the crosstalk between nuclides. In addition, a novel method of coincidence Compton imaging is proposed by taking coincidence detection between a Compton event and a photopeak events. Although this type of coincidence Compton imaging cannot specify the location, it can be useful in multi-nuclide Compton imaging.

The OpenGATE ecosystem for Monte Carlo simulation in medical physics

This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.

Joint image reconstruction algorithm in Compton cameras

Objective. To demonstrate the benefits of using an joint image reconstruction algorithm based on the List Mode Maximum Likelihood Expectation Maximization that combines events measured in different channels of information of a Compton camera. Approach. Both simulations and experimental data are employed to show the algorithm performance. Main results. The obtained joint images present improved image quality and yield better estimates of displacements of high-energy gamma-ray emitting sources. The algorithm also provides images that are more stable than any individual channel against the noisy convergence that characterizes Maximum Likelihood based algorithms. Significance. The joint reconstruction algorithm can improve the quality and robustness of Compton camera images. It also has high versatility, as it can be easily adapted to any Compton camera geometry. It is thus expected to represent an important step in the optimization of Compton camera imaging.

Machine learning-based event recognition in SiFi Compton camera imaging for proton therapy monitoring

Objective. Online monitoring of dose distribution in proton therapy is currently being investigated with the detection of prompt gamma (PG) radiation emitted from a patient during irradiation. The SiPM and scintillation Fiber based Compton Camera (SiFi-CC) setup is being developed for this aim. Approach. A machine learning approach to recognize Compton events is proposed, reconstructing the PG emission profile during proton therapy. The proposed method was verified on pseudo-data generated by a Geant4 simulation for a single proton beam impinging on a polymethyl methacrylate (PMMA) phantom. Three different models including the boosted decision tree (BDT), multilayer perception (MLP) neural network, and k-nearest neighbour (k-NN) were trained using 10-fold cross-validation and then their performances were assessed using the receiver operating characteristic (ROI) curves. Subsequently, after event selection by the most robust model, a software based on the List-Mode Maximum Likelihood Estimation Maximization (LM-MLEM) algorithm was applied for the reconstruction of the PG emission distribution profile. Main results. It was demonstrated that the BDT model excels in signal/background separation compared to the other two. Furthermore, the reconstructed PG vertex distribution after event selection showed a significant improvement in distal falloff position determination. Significance. A highly satisfactory agreement between the reconstructed distal edge position and that of the simulated Compton events was achieved. It was also shown that a position resolution of 3.5 mm full width at half maximum (FWHM) in distal edge position determination is feasible with the proposed setup.

Compton imaging for medical applications

Compton imaging exploits inelastic scattering, known as Compton scattering, using a Compton camera consisting of a scatterer detector in the front layer and an absorber detector in the back layer. This method was developed for astronomy, and in recent years, research and development for environmental and medical applications has been actively conducted. Compton imaging can discriminate gamma rays over a wide energy range from several hundred keV to several MeV. Therefore, it is expected to be applied to the simultaneous imaging of multiple nuclides in nuclear medicine and prompt gamma ray imaging for range verification in particle therapy. In addition, multiple gamma coincidence imaging is expected to be realized, which allows the source position to be determined from a single coincidence event using nuclides that emit multiple gamma rays simultaneously, such as nuclides that emit a single gamma ray simultaneously with positron decay. This review introduces various efforts toward the practical application of Compton imaging in the medical field, including in vivo studies, and discusses its prospects.

An accurate probabilistic model with detector resolution and Doppler broadening correction in list-mode MLEM reconstruction for Compton camera

Objective. The Compton cameras have been researched for medical applications and radioactive material detection. It is challenging for the Compton camera to realize high-resolution reconstruction when the incident photon energy is below 200 keV. However, multiple kinds of nuclear medical radionuclides are in this energy range, such as 201Tl, 67Ga, 99mTc, and 123I. In this work, we propose an improved probabilistic model with correction of detector energy resolution, detector spatial resolution, and Doppler broadening effect. The proposed model is used for numerical calculation of the system matrix in the list-mode maximum likelihood expectation maximization (LM-MLEM) algorithm. Approach. The model can improve the imaging resolution of LM-MLEM reconstruction by taking Doppler broadening effect into account. It performs well, especially in the following situations: low-energy photon incidence below 200 keV or (and) small distance between scattering and absorbing positions. Main results. Firstly, three main factors that affect the angular resolution of the Compton camera are theoretically analyzed and quantitatively calculated. The results of the analysis indicate the necessity of including the Doppler broadening effect in the model. Secondly, the details and derivation of the proposed probabilistic model are described. Thirdly, both Monte Carlo (MC) simulations and experiments are carried out to verify the performance of the proposed algorithm. The simulations focus on the low-energy reconstruction in which 201Tl (70 keV) and 99mTc (141 keV) are simulated. And the experiments are based on a single-layer Compton camera composed of a Timepix3 detector. Significance. The results of the simulations and the Timepix3-based experiments are presented to verify the effectiveness of the proposed algorithm. The model improves the Compton imaging resolution when the photon energy is below 200 keV.

Self-Collimating SPECT With Multi-Layer Interspaced Mosaic Detectors

Conventional single photon emission computed tomography (SPECT) relies on mechanical collimation whose resolution and sensitivity are interdependent, the best performance a SPECT system can attain is only a compromise of these two equally desired properties. To simultaneously achieve high resolution and sensitivity, we propose to use sensitive detectors constructed in a multi-layer in ter spaced mosaicdetectors (MATRICES) architecture to accomplish part of the collimation needed. We name this new approach self-collimation. We evaluate three self-collimating SPECT systems and report their imaging performance: 1) A simulated human brain SPECT achieves 3.88% sensitivity, it clearly resolves 0.5-mm and 1.0-mm hot-rod patterns at noise-free and realistic count-levels, respectively; 2) a simulated mouse SPECT achieves 1.25% sensitivity, it clearly resolves 50- [Formula: see text] and 100- [Formula: see text] hot-rod patterns at noise-free and realistic count-levels, respectively; 3) a SPECT prototype achieves 0.14% sensitivity and clearly separates 0.3-mm-diameter point sources of which the center-to-center neighbor distance is also 0.3 mm. Simulated contrast phantom studies show excellent resolution and signal-to-noise performance. The unprecedented system performance demonstrated by these 3 SPECT scanners is a clear manifestation of the superiority of the self-collimating approach over conventional mechanical collimation. It represents a potential paradigm shift in SPECT technology development.

A time-of-flight-based reconstruction for real-time prompt-gamma imaging in proton therapy

We propose a novel prompt-gamma (PG) imaging modality for real-time monitoring in proton therapy: PG time imaging (PGTI). By measuring the time-of-flight (TOF) between a beam monitor and a PG detector, our goal is to reconstruct the PG vertex distribution in 3D. In this paper, a dedicated, non-iterative reconstruction strategy is proposed (PGTI reconstruction). Here, it was resolved under a 1D approximation to measure a proton range shift along the beam direction. In order to show the potential of PGTI in the transverse plane, a second method, based on the calculation of the centre of gravity (COG) of the TIARA pixel detectors' counts was also explored. The feasibility of PGTI was evaluated in two different scenarios. Under the assumption of a 100 ps (rms) time resolution (achievable in single proton regime), MC simulations showed that a millimetric proton range shift is detectable at 2σwith 10⁸incident protons in simplified simulation settings. With the same proton statistics, a potential 2 mm sensitivity (at 2σwith 10⁸incident protons) to beam displacements in the transverse plane was found using the COG method. This level of precision would allow to act in real-time if the treatment does not conform to the 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σwith 10⁹incident protons. By showing the feasibility of a time-based algorithm for the reconstruction of the PG vertex distribution for a simplified anatomy, this work poses a theoretical basis for the future development of a PG imaging detector based on the measurement of particle TOF.

Influence of sub-nanosecond time of flight resolution for online range verification in proton therapy using the line-cone reconstruction in Compton imaging

Online ion range monitoring in hadron therapy can be performed via detection of secondary radiation, such as promptγ-rays, emitted during treatment. The promptγemission profile is correlated with the ion depth-dose profile and can be reconstructed via Compton imaging. The line-cone reconstruction, using the intersection between the primary beam trajectory and the cone reconstructed via a Compton camera, requires negligible computation time compared to iterative algorithms. A recent report hypothesised that time of flight (TOF) based discrimination could improve the precision of theγfall-off position (FOP) measured via line-cone reconstruction, where TOF comprises both the proton transit time from the phantom entrance untilγemission, and the flight time of theγ-ray to the detector. The aim of this study was to implement such a method and investigate the influence of temporal resolution on the precision of the FOP. Monte Carlo simulations of a 160 MeV proton beam incident on a homogeneous PMMA phantom were performed using GATE. The Compton camera consisted of a silicon-based scatterer and CeBr₃scintillator absorber. The temporal resolution of the detection system (absorber + beam trigger) was varied between 0.1 and 1.3 ns rms and a TOF-based discrimination method applied to eliminate unlikely solution(s) from the line-cone reconstruction. The FOP was obtained for varying temporal resolutions and its precision obtained from its shift across 100 independentγemission profiles compared to a high statistics reference profile. The optimal temporal resolution for the given camera geometry and 10⁸primary protons was 0.2 ns where a precision of 2.30 ± 0.15 mm (1σ) on the FOP was found. This precision is comparable to current state-of-the-art Compton imaging using iterative reconstruction methods or 1D imaging with mechanically collimated devices, and satisfies the requirement of being smaller than the clinical safety margins.

Advanced Monte Carlo simulations of emission tomography imaging systems with GATE

Built on top of the Geant4 toolkit, GATE is collaboratively developed for more than 15 years to design Monte Carlo simulations of nuclear-based imaging systems. It is, in particular, used by researchers and industrials to design, optimize, understand and create innovative emission tomography systems. In this paper, we reviewed the recent developments that have been proposed to simulate modern detectors and provide a comprehensive report on imaging systems that have been simulated and evaluated in GATE. Additionally, some methodological developments that are not specific for imaging but that can improve detector modeling and provide computation time gains, such as Variance Reduction Techniques and Artificial Intelligence integration, are described and discussed.

Proton range verification with MACACO II Compton camera enhanced by a neural network for event selection

The applicability extent of hadron therapy for tumor treatment is currently limited by the lack of reliable online monitoring techniques. An active topic of investigation is the research of monitoring systems based on the detection of secondary radiation produced during treatment. MACACO, a multi-layer Compton camera based on LaBr₃ scintillator crystals and SiPMs, is being developed at IFIC-Valencia for this purpose. This work reports the results obtained from measurements of a 150 MeV proton beam impinging on a PMMA target. A neural network trained on Monte Carlo simulations is used for event selection, increasing the signal to background ratio before image reconstruction. Images of the measured prompt gamma distributions are reconstructed by means of a spectral reconstruction code, through which the 4.439 MeV spectral line is resolved. Images of the emission distribution at this energy are reconstructed, allowing calculation of the distal fall-off and identification of target displacements of 3 mm.

Design of 4π Directional Radiation Detector based on Compton Scattering Effect

Obtaining directional information is required in many applications such as nuclear homeland security, contamination mapping after a nuclear incident and radiological events, or during the decontamination work. However, many directional radiation detectors are based on directional shielding, made of lead or tungsten collimators, introducing two main drawbacks. The first is the size and weight, making those detectors too heavy and irrelevant for utilization in handheld devices, drone mapping, or space applications. The second drawback is the limited field of view, which requires multiple detectors to cover the whole required field of view or machinery to rotate the narrow field of view detector. We propose a novel 4π directional detector based on a segmented hollow cubic detector, which uses the Compton effect interactions with no heavy collimators. The symmetrical cubical design provides both higher efficiency and 4π detection ability. Instead of traditional two types of detectors (scatterer and absorber) structure, we use the same type of detector, based on GAGG(Ce) scintillator coupled to silicon photomultiplier. Additional advantage of the proposed detector obtained by locating the photon sensors inside the detector, behind the scintillators, which improves the radiation hardness required for space applications. Furthermore, such arrangement flattens the temperature variation across the detector, providing better gain stability. The main advantage of the proposed detector is the ability of 4pi radiation detection for high energy gamma-rays without the use of heavy collimators.

Roadmap toward the 10 ps time-of-flight PET challenge

Since the seventies, positron emission tomography (PET) has become an invaluable medical molecular imaging modality with an unprecedented sensitivity at the picomolar level, especially for cancer diagnosis and the monitoring of its response to therapy. More recently, its combination with x-ray computed tomography (CT) or magnetic resonance (MR) has added high precision anatomic information in fused PET/CT and PET/MR images, thus compensating for the modest intrinsic spatial resolution of PET. Nevertheless, a number of medical challenges call for further improvements in PET sensitivity. These concern in particular new treatment opportunities in the context personalized (also called precision) medicine, such as the need to dynamically track a small number of cells in cancer immunotherapy or stem cells for tissue repair procedures. A better signal-to-noise ratio (SNR) in the image would allow detecting smaller size tumours together with a better staging of the patients, thus increasing the chances of putting cancer in complete remission. Moreover, there is an increasing demand for reducing the radioactive doses injected to the patients without impairing image quality. There are three ways to improve PET scanner sensitivity: improving detector efficiency, increasing geometrical acceptance of the imaging device and pushing the timing performance of the detectors. Currently, some pre-localization of the electron-positron annihilation along a line-of-response (LOR) given by the detection of a pair of annihilation photons is provided by the detection of the time difference between the two photons, also known as the time-of-flight (TOF) difference of the photons, whose accuracy is given by the coincidence time resolution (CTR). A CTR of about 10 picoseconds FWHM will ultimately allow to obtain a direct 3D volume representation of the activity distribution of a positron emitting radiopharmaceutical, at the millimetre level, thus introducing a quantum leap in PET imaging and quantification and fostering more frequent use of 11C radiopharmaceuticals. The present roadmap article toward the advent of 10 ps TOF-PET addresses the status and current/future challenges along the development of TOF-PET with the objective to reach this mythic 10 ps frontier that will open the door to real-time volume imaging virtually without tomographic inversion. The medical impact and prospects to achieve this technological revolution from the detection and image reconstruction point-of-views, together with a few perspectives beyond the TOF-PET application are discussed.

Prompt gamma imaging for the identification of regional proton range deviations due to anatomic change in a heterogeneous region

OBJECTIVES

Prompt gamma (PG) imaging has previously been demonstrated for use in proton range verification of a brain treatment with a homogeneous target region. In this study, the feasibility of PG imaging to detect anatomic change within a heterogeneous region is presented.

METHODS

A prompt gamma camera recorded several fractions of a patient treatment to the base of skull. An evaluation CT revealed a decrease in sinus cavity filling during the treatment course. Comparison of PG profiles between measurement and simulation was performed to investigate range variations between planned and measured pencil beam spot positions.

RESULTS

For one field, an average over range of 3 mm due to the anatomic change could be detected for a subset of spots traversing the sinus cavity region. The two other fields appeared less impacted by the change but predicted range variations could not be detected. These results were partially consistent with the simulations of the evaluation CT.

CONCLUSION

We report the first clinical application of PG imaging that detected some of the expected small regional proton range deviations due to anatomic change in a heterogeneous region. However, several limitations exist with the technology that may limit its sensitivity to detect range deviations in heterogeneous regions.

ADVANCES IN KNOWLEDGE

We report on the first detection of range variations due to anatomic change in a heterogeneous region using PGI. The results confirm the feasibility of using PG-based range verification in highly heterogeneous target regions to identify deviations from the treatment plan.

Image reconstruction for a multi-layer Compton telescope: an analytical model for three interaction events

Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enable in vivo range verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.

A spectral reconstruction algorithm for two-plane Compton cameras

One factor limiting the current applicability extent of hadron therapy is the lack of a reliable method for real time treatment monitoring. The use of Compton imaging systems as monitors requires the correct reconstruction of the distribution of prompt gamma productions during patient irradiation. In order to extract the maximum information from all the measurable events, we implemented a spectral reconstruction method that assigns to all events a probability of being either partial or total energy depositions. The method, implemented in a list-mode maximum likelihood expectation maximization algorithm, generates a four dimensional image in the joint spatial-spectral domain, in which the voxels containing the emission positions and energies are obtained. The analytical model used for the system response function is also employed to derive an analytical expression for the sensitivity, which is calculated via Monte Carlo integration. The performance of the method is evaluated through reconstruction of various experimental and simulated sources with different spatial and energy distributions. The results show that the proposed method can recover the spectral and spatial information simultaneously, but only under the assumption of ideal measurements. The analysis of the Monte Carlo simulations has led to the identification of two important degradation sources: the mispositioning of the gamma interaction point and the missing energy recorded in the interaction. Both factors are related to the high energy transferred to the recoil electrons, which can travel far from the interaction point and even escape the detector. These effects prevent the direct application of the current method in more realistic scenarios. Nevertheless, experimental point-like sources have been accurately reconstructed and the spatial distributions and spectral emission of complex simulated phantoms can be identified.

Electron-tracking Compton camera imaging of technetium-95m

Imaging was conducted using an electron tracking-Compton camera (ETCC), which measures γ-rays with energies in the range of 200-900 keV from 95mTc. 95mTc was produced by the 95Mo(p, n)95mTc reaction on a 95Mo-enriched target. A method for recycling 95Mo-enriched molybdenum trioxide was employed, and the recycled yield of 95Mo was 70%-90%. Images were obtained with the gate of three energies. The results showed that the spatial resolution increases with increasing γ-ray energy, and suggested that the ETCC with high-energy γ-ray emitters such as 95mTc is useful for the medical imaging of deep tissue and organs in the human body.

3D prompt gamma imaging for proton beam range verification

We tested the ability of a single Compton camera (CC) to produce 3-dimensional (3D) images of prompt gammas (PGs) emitted during the irradiation of a tissue-equivalent plastic phantom with proton pencil beams for clinical doses delivered at clinical dose rates. PG measurements were made with a small prototype CC placed at three different locations along the proton beam path. We evaluated the ability of the CC to produce images at each location for two clinical scenarios: (1) the delivery of a single 2 Gy pencil beam from a hypo-fractionated treatment (~9  ×  10⁸ protons), and (2) a single pencil beam from a standard treatment (~1  ×  10⁸ protons). Additionally, the data measured at each location were combined to simulate measurements with a larger scale, clinical CC and its ability to image shifts in the Bragg peak (BP) range for both clinical scenarios. With our prototype CC, the location of the distal end of the BP could be seen with the CC placed up to 4 cm proximal or distal to the BP distal falloff. Using the data from the simulated full scale clinical CC, 3D images of the PG emission were produced with the delivery of as few as 1  ×  10⁸ protons, and shifts in the proton beam range as small as 2 mm could be detected for delivery of a 2 Gy spot. From these results we conclude that 3D PG imaging for proton range verification under clinical beam delivery conditions is possible with a single CC.