Improved iterative reconstruction method for Compton imaging using median filter.

A Compton camera is a device for imaging a radio-source distribution without using a mechanical collimator. Ordered-subset expectation-maximization (OS-EM) is widely used to reconstruct Compton images. However, the OS-EM algorithm tends to over-concentrate and amplify noise in the reconstructed image. It is, thus, necessary to optimize the number of iterations to develop high-quality images, but this has not yet been achieved. In this paper, we apply a median filter to an OS-EM algorithm and introduce a median root prior expectation-maximization (MRP-EM) algorithm to overcome this problem. In MRP-EM, the median filter is used to update the image in each iteration. We evaluated the quality of images reconstructed by our proposed method and compared them with those reconstructed by conventional algorithms using mathematical phantoms. The spatial resolution was estimated using the images of two point sources. Reproducibility was evaluated on an ellipsoidal phantom by calculating the residual sum of squares, zero-mean normalized cross-correlation, and mutual information. In addition, we evaluated the semi-quantitative performance and uniformity on the ellipsoidal phantom. MRP-EM reduces the generated noise and is robust with respect to the number of iterations. An evaluation of the reconstructed image quality using some statistical indices shows that our proposed method delivers better results than conventional techniques.

Annihilation gamma imaging for carbon ion beam range monitoring using Si/CdTe Compton camera.

In this study, we performed on-beam monitoring of 511 keV annihilation gamma emissions using a Compton camera. Beam monitoring experiments were conducted using carbon ion beams of 290 MeV/u irradiated on a polymethyl methacrylate (PMMA) phantom. The intensity of the beams was 3 × 109 particles per pulse, with 20 pulses per minute. A Compton camera based on a silicon/cadmium telluride (Si/CdTe) detector was used to monitor the annihilation gamma rays emitted from the phantom. We successfully reconstructed the energy events of 511 keV annihilation gamma rays and developed Compton images using a simple back-projection method. The distribution of the annihilation gamma ray generation traced the beam trajectory and the peak intensity position was a few millimeters shorter than the Bragg peak position. Moreover, the effect of the beam range shifter with 30, 60, and 90 mm water equivalent thickness (WET) was clearly visualized in the reconstructed Compton images. The experimentally measured values of the corresponding range shifts in the PMMA phantom (28.70 mm, 52.49 mm, and 76.77 mm, respectively) were consistent with the shifts of the Bragg peak position (25.50 mm, 51.30 mm and 76.70 mm, respectively) evaluated by Monte Carlo simulation. The results show that the Si/CdTe Compton camera has strong potential for on-beam monitoring of annihilation gamma rays in particle therapy in clinical situations.

Effect of number of views on cross-sectional Compton imaging: A fundamental study with backprojection.

INTRODUCTION: We have been developing a medical imaging technique using a Compton camera, which is expected to reconstruct three-dimensional images. If the number of views is not sufficient, star-shaped artifacts (streak artifacts) could arise in cross-sectional images. Therefore, we estimated the point spread function (PSF) of cross-sectional Compton images and the effect of the number of views by Monte Carlo simulations and experimental studies. MATERIALS AND METHODS: A cross-sectional Compton image was reconstructed using a dataset comprising 719 view directions and PSF was analyzed using a radial distribution. The peak height, full width at half maximum (FWHM), background (BG), and residual sum of squares (RSS) were calculated from the obtained PSF. In addition, RSSs were plotted against the number of views to estimate the required number to suppress star-shaped artifacts. RESULTS: There was no correlation found between the number of views and both FWHM (16 mm) and peak/BG ratio (∼1 × 104). RSSs were reduced with the number of views and approached the minimum asymptotically. Correlation was observed between the required number of views and the number of Compton events used for image reconstruction. CONCLUSION: We determined the PSF of cross-sectional Compton images and the effect of the number of views on the images. The required number of views to suppress the star-shaped artifact is related to the square root of the number of Compton events used to reconstruct the image. From this study, we concluded that 21 or more views are required for clinical purposes.

Imaging of monochromatic beams by measuring secondary electron bremsstrahlung for carbon-ion therapy using a pinhole x-ray camera.

A feasibility study on the imaging of monochromatic carbon-ion beams for carbon-ion therapy was performed. The evaluation was based on Monte Carlo simulations and beam-irradiation experiments, using a pinhole x-ray camera, which measured secondary electron bremsstrahlung (SEB). The simulation results indicated that the trajectories of the carbon-ion beams with injection energies of 278, 249 and 218 MeV/u in a water phantom, were clearly imaged by measuring the SEB with energies from 30 to 60 keV, using a pinhole camera. The Bragg-peak positions for these three injection energies were located at the positions where the ratios of the counts of SEB acquisitions to the maximum counts were approximately 0.23, 0.26 and 0.29, respectively. Moreover, we experimentally demonstrated that it was possible to identify the Bragg-peak positons, at the positions where the ratios coincided with the simulation results. However, the estimated Bragg-peak positions for the injection energies of 278 and 249 MeV/u were slightly deeper than the expected positions. In conclusion, for both the simulations and experiments, we found that the 25 mm shifts in the Bragg-peak positions can be observed by this method.

Monte Carlo simulation of photon emission below a few hundred kiloelectronvolts for beam monitoring in carbon ion therapy.

The purpose of this study is to determine whether the main component of the low-energy (63-68 keV) particles emitted perpendicularly to the 12C beam from the 12C-irradiated region in a water phantom is secondary electron bremsstrahlung (SEB). Monte Carlo simulations of a 12C-beam (290 MeV/u) irradiated on a water phantom were performed. A detector was placed beside the water phantom with a lead collimator between the phantom and the detector. To move the Bragg-peak position, a binary filter was placed in an upper stream of the phantom. The energy distributions of the particles incident on the detector and those deposited in the detector were analyzed. The simulation was also performed with suppressed delta-ray and/or bremsstrahlung generation to identify the SEB components. It was found that the particles incident on the detector were predominantly photons and neutrons. The yields of the photons and energy deposition decreased with the suppression of SEB generation. It is concluded that one of the predominant components of the yields in the regions shallower than the Bragg-peak position is due to SEB generation, and these components become significantly smaller in regions deeper than the Bragg-peak position.

Beam range estimation by measuring bremsstrahlung.

We describe a new method for estimating the beam range in heavy-ion radiation therapy by measuring the ion beam bremsstrahlung. We experimentally confirm that the secondary electron bremsstrahlung process provides the dominant bremsstrahlung contribution. A Monte Carlo simulation shows that the number of background photons from annihilation gamma rays is about 1% of the bremsstrahlung strength in the low-energy region used in our estimation (63-68 keV). Agreement between the experimental results and the theoretical prediction for the characteristic shape of the bremsstrahlung spectrum validates the effectiveness of our new method in estimating the ion beam range.

Targeted delivery of chemotherapeutic agents using improved radiosensitive liquid core microcapsules and assessment of their antitumor effect.

PURPOSE: Radiation-sensitive microcapsules composed of alginate and hyaluronic acid are being developed. We report the development of improved microcapsules that were prepared using calcium- and yttrium-induced polymerization. We previously reported on the combined antitumor effect of carboplatin-containing microcapsules and radiotherapy. METHODS AND MATERIALS: We mixed a 0.1% (wt/vol) solution of hyaluronic acid with a 0.2% alginate solution. Carboplatin (l mg) and indocyanine green (12.5 microg) were added to this mixture, and the resultant material was used for capsule preparation. The capsules were prepared by spraying the material into a mixture containing a 4.34% CaCl(2) solution supplemented with 0-0.01% yttrium. These capsules were irradiated with single doses of 0.5, 1.0, 1.5, or 2 Gy (60)Co gamma-rays. Immediately after irradiation, the frequency of microcapsule decomposition was determined using a microparticle-induced X-ray emission camera. The amount of core content released was estimated by particle-induced X-ray emission and colorimetric analysis with 0.25% indocyanine green. The antitumor effect of the combined therapy was determined by monitoring its effects on the diameter of an inoculated Meth A fibrosarcoma. RESULTS: Microcapsules that had been polymerized using a 4.34% CaCl(2) solution supplemented with 5.0 x 10(-3)% (10(-3)% meant or 10%(-3)) yttrium exhibited the maximal decomposition, and the optimal release of core content occurred after 2-Gy irradiation. The microcapsules exhibited a synergistic antitumor effect combined with 2-Gy irradiation and were associated with reduced adverse effects. CONCLUSION: The results of our study have shown that our liquid core microcapsules can be used in radiotherapy for targeted delivery of chemotherapeutic agents.