A prototype fan-beam optical CT scanner for 3D dosimetry

A cone-beam optical CT based on a convergent light source - Characterization and optimization

PURPOSE

Employing a Fresnel lens and a point-like light source to create a convergent light beam for the camera effectively minimizes stray light and enhances image quality in optical computed tomography (OCT), benefiting 3D dosimetry applications. This study outlines the development of an economical cone-beam optical computed scanner for 3D dosimetry.

METHODS

Optical performance was assessed by calculating modulation transfer function (MTF) with pattern charts. Stray light was evaluated by imaging a cylinder flask and a square grid with 5 mm diameter holes to determine the stray-to-primary ratio. Reconstruction quality was determined using SIRT-TV and compared with spectrophotometry attenuation coefficients, with the best regularization parameter (λ = 0.01) chosen based on contrast-to-noise ratio (CNR). Dosimetry performance was assessed by determining percentage dose depth (PDD) for a 6MV beam with a 5 × 5 cm2 field using FXO-f gel dosimeter, compared with ionization chamber data.

RESULTS

MTF evaluation yielded ≥ 50 % agreement with pattern charts. Stray-to-primary ratio was less than 0.1 or 10 % of the total signal. Reconstruction showed low noise and artifacts, with optimal CNR at λ = 0.01. Attenuation coefficients from optical CT aligned with spectrometer measurements within 1.2 %. PDD calculated with FXO-f gel dosimeter closely matched ionization chamber data (<1.2 % difference), achieving a dose resolution of 0.1 Gy.

CONCLUSION

The built and optimization the de optical-CT based on a convergent beam is read to perform the 3D quality assurance in clinical applications.

Iterative image reconstruction algorithm analysis for optical CT radiochromic gel dosimetry

Background.Modern radiation therapy technologies aim to enhance radiation dose precision to the tumor and utilize hypofractionated treatment regimens. Verifying the dose distributions associated with these advanced radiation therapy treatments remains an active research area due to the complexity of delivery systems and the lack of suitable three-dimensional dosimetry tools. Gel dosimeters are a potential tool for measuring these complex dose distributions. A prototype tabletop solid-tank fan-beam optical CT scanner for readout of gel dosimeters was recently developed. This scanner does not have a straight raypath from source to detector, thus images cannot be reconstructed using filtered backprojection (FBP) and iterative techniques are required.Purpose.To compare a subset of the top performing algorithms in terms of image quality and quantitatively determine the optimal algorithm while accounting for refraction within the optical CT system. The following algorithms were compared: Landweber, superiorized Landweber with the fast gradient projection perturbation routine (S-LAND-FGP), the fast iterative shrinkage/thresholding algorithm with total variation penalty term (FISTA-TV), a monotone version of FISTA-TV (MFISTA-TV), superiorized conjugate gradient with the nonascending perturbation routine (S-CG-NA), superiorized conjugate gradient with the fast gradient projection perturbation routine (S-CG-FGP), superiorized conjugate gradient with with two iterations of CG performed on the current iterate and the nonascending perturbation routine (S-CG-2-NA).Methods.A ray tracing simulator was developed to track the path of light rays as they traverse the different mediums of the optical CT scanner. Two clinical phantoms and several synthetic phantoms were produced and used to evaluate the reconstruction techniques under known conditions. Reconstructed images were analyzed in terms of spatial resolution, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), signal non-uniformity (SNU), mean relative difference (MRD) and reconstruction time. We developed an image quality based method to find the optimal stopping iteration window for each algorithm. Imaging data from the prototype optical CT scanner was reconstructed and analysed to determine the optimal algorithm for this application.Results.The optimal algorithms found through the quantitative scoring metric were FISTA-TV and S-CG-2-NA. MFISTA-TV was found to behave almost identically to FISTA-TV however MFISTA-TV was unable to resolve some of the synthetic phantoms. S-CG-NA showed extreme fluctuations in the SNR and CNR values. S-CG-FGP had large fluctuations in the SNR and CNR values and the algorithm has less noise reduction than FISTA-TV and worse spatial resolution than S-CG-2-NA. S-LAND-FGP had many of the same characteristics as FISTA-TV; high noise reduction and stability from over iterating. However, S-LAND-FGP has worse SNR, CNR and SNU values as well as longer reconstruction time. S-CG-2-NA has superior spatial resolution to all algorithms while still maintaining good noise reduction and is uniquely stable from over iterating.Conclusions.Both optimal algorithms (FISTA-TV and S-CG-2-NA) are stable from over iterating and have excellent edge detection with ESF MTF 50% values of 1.266 mm-1and 0.992 mm-1. FISTA-TV had the greatest noise reduction with SNR, CNR and SNU values of 424, 434 and 0.91 × 10-4, respectively. However, low spatial resolution makes FISTA-TV only viable for large field dosimetry. S-CG-2-NA has better spatial resolution than FISTA-TV with PSF and LSF MTF 50% values of 1.581 mm-1and 0.738 mm-1, but less noise reduction. S-CG-2-NA still maintains good SNR, CNR, and SNU values of 168, 158 and 1.13 × 10-4, respectively. Thus, S-CG-2-NA is a well rounded reconstruction algorithm that would be the preferable choice for small field dosimetry.

Validation of complex radiotherapy techniques using polymer gel dosimetry

Modern radiotherapy delivers highly conformal dose distributions to irregularly shaped target volumes while sparing the surrounding normal tissue. Due to the complex planning and delivery techniques, dose verification and validation of the whole treatment workflow by end-to-end tests became much more important and polymer gel dosimeters are one of the few possibilities to capture the delivered dose distribution in 3D. The basic principles and formulations of gel dosimetry and its evaluation methods are described and the available studies validating device-specific geometrical parameters as well as the dose delivery by advanced radiotherapy techniques, such as 3D-CRT/IMRT and stereotactic radiosurgery treatments, the treatment of moving targets, online-adaptive magnetic resonance-guided radiotherapy as well as proton and ion beam treatments, are reviewed. The present status and limitations as well as future challenges of polymer gel dosimetry for the validation of complex radiotherapy techniques are discussed.

Iterative image reconstruction with polar coordinate discretized system matrix for optical CT radiochromic gel dosimetry

BACKGROUND

Gel dosimeters are a potential tool for measuring the complex dose distributions that characterize modern radiotherapy. A prototype tabletop solid-tank fan-beam optical CT scanner for readout of gel dosimeters was recently developed. This scanner does not have a straight raypath from source to detector, thus images cannot be reconstructed using filtered backprojection (FBP) and iterative techniques are required. Iterative image reconstruction requires a system matrix that describes the geometry of the imaging system. Stored system matrices can become immensely large, making them impractical for storage on a typical desktop computer.

PURPOSE

Here we develop a method to reduce the storage size of optical CT system matrices through use of polar coordinate discretization while accounting for the refraction in optical CT systems.

METHODS

A ray tracing simulator was developed to track the path of light rays as they traverse the different mediums of the optical CT scanner. Cartesian coordinate discretized system matrices (CCDSMs) and polar coordinate discretized system matrices (PCDSMs) were generated by discretizing the reconstruction area of the optical CT scanner into a Cartesian pixel grid and a polar coordinate pixel grid, respectively. The length of each ray through each pixel was calculated and used to populate the system matrices. To ensure equal weighting during iterative reconstruction, the radial rings of PCDSMs were asymmetrically spaced such that the area of each polar pixel was constant. Two clinical phantoms and several synthetic phantoms were produced and used to evaluate the reconstruction techniques under known conditions. Reconstructed images were analyzed in terms of spatial resolution, signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), signal nonuniformity (SNU), and Gamma map pass percentage.

RESULTS

A storage size reduction of 99.72% was found when comparing a PCDSM to a CCDSM with the same total number of pixels. Images reconstructed with a PCDSM were found to have superior SNR, CNR, SNU, and Gamma (1 mm, 1%) pass percentage compared to those reconstructed with a CCDSM. Increasing spatial resolution in the radial direction with increasing radial distance was found in both PCDSM and CCDSM reconstructions due to the outer regions refracting light more severely. Images reconstructed with a PCDSM showed a decrease in spatial resolution in the azimuthal directions as radial distance increases, due to the widening of the polar pixels. However, this can be mitigated with only a slight increase in storage size by increasing the number of projections. A loss of spatial resolution in the radial direction within 5 mm radially from center was found when reconstructing with a PCDSM, due to the large innermost pixels. However, this was remedied by increasing the number of radial rings within the PCDSM, yielding radial spatial resolution on par with images reconstructed with a CCDSM and a storage size reduction of 99.26%.

CONCLUSIONS

Discretizing the image pixel elements in polar coordinates achieved a system matrix storage size reduction of 99.26% with only minimal reduction in the image quality.

Commissioning of a solid tank design for fan-beam optical CT based 3D radiation dosimetry

Objective.Optical computed tomography (CT) is one of the leading modalities for imaging gel dosimeters used in the verification of complex radiotherapy treatments. In previous work, a novel fan-beam optical CT scanner design was proposed that could significantly reduce the volume of the refractive index baths that are commonly found in optical CT systems. Here, the proposed scanner has been manufactured and commissioned.Approach.Image reconstruction is performed through algebraic reconstruction technique and iterated using the fast iterative shrinkage-thresholding algorithm (FISTA) algorithm. Ray tracing for algebraic reconstruction was performed using an in-house developed ray tracing simulator. A set of Sylgard® 184 phantoms were created to commission spatial resolution, geometric deformity, contrast-to-noise ratio (CNR), and scan settings.Main Results.The scanner is capable of a 0.929 mm-1spatial resolution, observed at 200 iterations, although the spatial resolution is highly dependent on the number of iterations. The geometric distortion, measured by scanning a needle phantom with the prototype scanner as well as a conventional x-ray CT was found to be within <0.25 mm. The CNR was found to peak between 65 and 190 occurring between 50 and 100 iterations and was highly dependent on the region chosen for background noise calculation. The proposed scanner is capable of scanning and reading out slices in less than 1 min per slice.Significance.This work displays the viability of a fan-beam optical CT scanner with minimal index matching using ray-traced algebraic reconstruction.

Radiation Dosimetry by Use of Radiosensitive Hydrogels and Polymers: Mechanisms, State-of-the-Art and Perspective from 3D to 4D

Gel dosimetry was developed in the 1990s in response to a growing need for methods to validate the radiation dose distribution delivered to cancer patients receiving high-precision radiotherapy. Three different classes of gel dosimeters were developed and extensively studied. The first class of gel dosimeters is the Fricke gel dosimeters, which consist of a hydrogel with dissolved ferrous ions that oxidize upon exposure to ionizing radiation. The oxidation results in a change in the nuclear magnetic resonance (NMR) relaxation, which makes it possible to read out Fricke gel dosimeters by use of quantitative magnetic resonance imaging (MRI). The radiation-induced oxidation in Fricke gel dosimeters can also be visualized by adding an indicator such as xylenol orange. The second class of gel dosimeters is the radiochromic gel dosimeters, which also exhibit a color change upon irradiation but do not use a metal ion. These radiochromic gel dosimeters do not demonstrate a significant radiation-induced change in NMR properties. The third class is the polymer gel dosimeters, which contain vinyl monomers that polymerize upon irradiation. Polymer gel dosimeters are predominantly read out by quantitative MRI or X-ray CT. The accuracy of the dosimeters depends on both the physico-chemical properties of the gel dosimeters and on the readout technique. Many different gel formulations have been proposed and discussed in the scientific literature in the last three decades, and scanning methods have been optimized to achieve an acceptable accuracy for clinical dosimetry. More recently, with the introduction of the MR-Linac, which combines an MRI-scanner and a clinical linear accelerator in one, it was shown possible to acquire dose maps during radiation, but new challenges arise.

Optimization of solid tank design for fan-beam optical CT based 3D radiation dosimetry

Optical computed tomography (CT) is one of the leading modalities for imaging gel dosimeters for 3D radiation dosimetry. There exist multiple scanner designs that have showcased excellent 3D dose verification capabilities of optical CT gel dosimetry. However, due to multiple experimental and reconstruction based factors there is currently no single scanner that has become a preferred standard. A significant challenge with setup and maintenance can be attributed to maintaining a large refractive index bath (1-15 l). In this work, a prototype solid 'tank' optical CT scanner is proposed that minimizes the volume of refractive index bath to between 10 and 35 ml. A ray-path simulator was created to optimize the design such that the solid tank geometry maximizes light collection across the detector array, maximizes the volume of the dosimeter scanned, and maximizes the collected signal dynamic range. An objective function was created to score possible geometries, and was optimized to find a local maximum geometry score from a set of possible design parameters. The design parameters optimized include the block length x _bl , bore position x _bc , fan-laser position x _lp , lens block face semi-major axis length x _ma , and the lens block face eccentricity x _be . For the proposed design it was found that each of these parameters can have a significant effect on the signal collection efficacy within the scanner. Simulations scores are specific to the attenuation characteristics and refractive index of a simulated dosimeter. It was found that for a FlexyDos3D dosimeter, the ideal values for each of the five variables were: x _bl = 314 mm, x _bc = 6.5 mm, x _lp = 50 mm, x _ma = 66 mm, and x _be = 0. In addition, a ClearView^™ dosimeter was found to have ideal values at: x _bl = 204 mm, x _bc = 13 mm, x _lp = 58 mm, x _ma = 69 mm, and x _be = 0. The ray simulator can also be used for further design and testing of new, unique and purpose-built optical CT geometries.

Determination of paramagnetic ferrous gel sensitivity in low energy x-ray beam produced by a miniature accelerator

The INTRABEAM Carl Zeiss Surgical system (Oberkochen, Germany) is a miniature accelerator producing low energy photons (50 keV maximum). The published dosimetric characterization of the INTRABEAM was based on detectors (radiochromic films or ionization chambers) not allowing measuring the absorbed dose in the first millimeters of the irradiated medium, where the dose is actually prescribed. This study aims at determining with Magnetic Resonance Imaging (MRI) the sensitivity of a paramagnetic gel in order to measure the dose deposit produced with the INTRABEAM from 0 to 20 mm. Although spherical applicators are mostly used with the INTRABEAM system for breast applications, this study focuses on surface applicators that are of interest for cutaneous carcinomas. The irradiations at 12 different dose levels (between 2 Gy and 50 Gy at the gel surface) were performed with the INTRABEAM and a 4 cm surface applicator. The gel used in this study is a new « sensitive » material. In order to compare gel sensitivity at low energy with high energy, gels were irradiated by an 18 MV photon beam produced by a Varian Clinac 2100 CD. T₂ weighted multi echo MRI sequences were performed with 16 echo times. The T₂ signal versus echo times was fitted with a mono-exponential function with 95% confidence interval. The calibration curve determined at low energy is a linear function (r² = 0.9893) with a sensitivity of 0.0381 s^-1.Gy^-1, a similar linear function was obtained at high energy (0.0372 s^-1.Gy^-1 with r² = 0.9662). The calibration curve at low energy was used to draw isodose maps from the MR images. The PDD (Percent Depth Dose) determined in the gel is within 5%-1mm of the ionization chamber PDD except for one point. The dosimetric sensitivity of this new paramagnetic ferrous gel was determined with MRI measurements. It allowed measuring the dose distribution specifically in the first millimeters for an irradiation with the INTRABEAM miniature accelerator equipped with a surface applicator.

Fixed, object-specific intensity compensation for cone beam optical CT radiation dosimetry

Optical cone beam computed tomography (CT) scanning of radiochromic gel dosimeters, using a CCD camera and a low stray light convergent source, provides fast, truly 3D radiation dosimetry with high accuracy. However, a key limiting factor in radiochromic gel dosimetry at large (⩾10 cm diameter) volumes is the initial attenuation of the dosimeters. It is not unusual to observe a 5-10×  difference in signal intensity through the dosimeter center versus through the surrounding medium in pre-irradiation images. Thus, all dosimetric information in a typical experiment is measured within the lower 10%-20% of the camera sensor's range, and re-use of gels is often not possible due to a lack of transmission. To counteract this, in this note we describe a simple method to create source compensators by printing on transparent films. This technique, which is easily implemented and inexpensive, is an optical analogue to the bowtie filter in x-ray CT. We present transmission images and solution phantom reconstructions to demonstrate that (1) placing compensators beyond the focal zone of the imaging lens prevents high spatial frequency features of the printed films from generating reconstruction artifacts, and (2) object-specific compensation considerably reduces the range of intensities measured in projection images. This will improve the measurable dose range in optical CT dosimetry, and will enable imaging of larger gel volumes (∼15 cm diameter). Additionally, it should enable re-use of dosimeters by printing a new compensator for a second experiment.

A fast dual wavelength laser beam fluid-less optical CT scanner for radiotherapy 3D gel dosimetry I: design and development

Three dimensional dosimetry by optical CT readout of radiosensitive gels or solids has previously been indicated as a solution for measurement of radiotherapy 3D dose distributions. The clinical uptake of these dosimetry methods has been limited, partly due to impracticalities of the optical readout such as the expertise and labour required for refractive index fluid matching. In this work a fast laser beam optical CT scanner is described, featuring fluid-less and dual wavelength operation. A second laser with a different wavelength is used to provide an alternative reference scan to the commonly used pre-irradiation scan. Transmission data for both wavelengths is effectively acquired simultaneously, giving a single scan process. Together with the elimination of refractive index fluid matching issues, scanning practicality is substantially improved. Image quality and quantitative accuracy were assessed for both dual and single wavelength methods. The dual wavelength scan technique gave improvements in uniformity of reconstructed optical attenuation coefficients in the sample 3D volume. This was due to a reduction of artefacts caused by scan to scan changes. Optical attenuation measurement accuracy was similar for both dual and single wavelength modes of operation. These results established the basis for further work on dosimetric performance.

Three-dimensional radiation dosimetry using polymer gel and solid radiochromic polymer: From basics to clinical applications

Accurate dose measurement tools are needed to evaluate the radiation dose delivered to patients by using modern and sophisticated radiation therapy techniques. However, the adequate tools which enable us to directly measure the dose distributions in three-dimensional (3D) space are not commonly available. One such 3D dose measurement device is the polymer-based dosimeter, which changes the material property in response to radiation. These are available in the gel form as polymer gel dosimeter (PGD) and ferrous gel dosimeter (FGD) and in the solid form as solid plastic dosimeter (SPD). Those are made of a continuous uniform medium which polymerizes upon irradiation. Hence, the intrinsic spatial resolution of those dosimeters is very high, and it is only limited by the method by which one converts the dose information recorded by the medium to the absorbed dose. The current standard methods of the dose quantification are magnetic resonance imaging, optical computed tomography, and X-ray computed tomography. In particular, magnetic resonance imaging is well established as a method for obtaining clinically relevant dosimetric data by PGD and FGD. Despite the likely possibility of doing 3D dosimetry by PGD, FGD or SPD, the tools are still lacking wider usages for clinical applications. In this review article, we summarize the current status of PGD, FGD, and SPD and discuss the issue faced by these for wider acceptance in radiation oncology clinic and propose some directions for future development.

Role of gel dosimeters in boron neutron capture therapy

Gel dosimeters have acquired a unique status in radiotherapy, especially with the advent of the new techniques in which there is a need for three-dimensional dose measurement with high spatial resolution. One of the techniques in which the use of gel dosimeters has drawn the attention of the researchers is the boron neutron capture therapy. Exploring the history of gel dosimeters, this paper sets out to study their role in the boron neutron capture therapy dosimetric process.

Investigations into the feasibility of optical-CT 3D dosimetry with minimal use of refractively matched fluids

PURPOSE

In optical-CT, the use of a refractively matched polyurethane solid-tank in place of a fluid bath has the potential to greatly increase practical convenience, reduce cost, and possibly improve the efficacy of flood corrections. This work investigates the feasibility of solid-tank optical-CT imaging for 3D dosimetry through computer simulation.

METHODS

A matlab ray-tracing simulation platform, ScanSim, was used to model a parallel-source telecentric optical-CT imaging system through a polyurethane solid-tank containing a central cylindrical hollow into which PRESAGE radiochromic dosimeters can be placed. A small amount of fluid fills the 1-5 mm gap between the dosimeter and the walls of the tank. The use of the solid-tank reduces the required amount of fluid by approximately 97%. To characterize the efficacy of solid-tank, optical-CT scanning simulations investigated sensitivity to refractive index (RI) mismatches between dosimeter, solid-tank, and fluid, for a variety of dosimeter (RI = 1.5-1.47) and fluid (RI = 1.55-1.0) combinations. Efficacy was evaluated through the usable radius (ru) metric, defined as the fraction of the radius of the dosimeter where measured dose is predicted to be within 2% of the ground truth entered into the simulation. Additional simulations examined the effect of increasing gap size (1-5 mm) between the dosimeter and solid-tank well. The effects of changing the lens tolerance (0.5°-5.0°) were also investigated.

RESULTS

As the RI mismatch between the dosimeter and solid-tank increased from 0 to 0.02, the usable radius decreased from 97.6% to 50.2%. The optimal fluid RI decreased nonlinearly from 1.5 to 1.34 as the mismatch increased and was up to 9% lower than the tank. Media mismatches between the dosimeter and solid-tank also exacerbate the effects of changing the gap size, with no easily quantifiable relationship with usable radius. Generally, the optimal fluid RI value increases as gap size increases and is closely matched to the dosimeter at large gap sizes (> 3 mm). Increasing the telecentric lens tolerance increases the usable radius for all refractive media combinations and improves the maximum usable radius of mismatched media to that of perfectly matched media for tolerances > 5.0°. The maximum usable radius can be improved up to a factor of 2 when lens tolerances are small (< 1.0°).

CONCLUSIONS

Dry solid-tank optical-CT imaging in a telecentric system is feasible if the dosimeter RI is a close match with the solid-tank (< 0.01 difference), providing accurate dose measurements within ± 2% of true dose to over 80% of the dosimeter volume. In order to achieve accurate measurements over 96% of the dosimeter volume (representing out to 2 mm from the dosimeter edge), the dosimeter-tank RI mismatch must be less than 0.005. Optimal results occur when the RI of the dosimeter and tank is the same, in which case the fluid will have the same RI. If mismatches between the tank and dosimeter RI occur, the RI of the matching fluid needs to be fine tuned to achieve the highest usable radius.

Incorporating multislice imaging into x-ray CT polymer gel dosimetry

PURPOSE

To evaluate multislice computed tomography (CT) scanning for fast and reliable readout of radiation therapy (RT) dose distributions using CT polymer gel dosimetry (PGD) and to establish a baseline assessment of image noise and uniformity in an unirradiated gel dosimeter.

METHODS

A 16-slice CT scanner was used to acquire images through a 1 L cylinder filled with water. Additional images were collected using a single slice machine. The variability in CT number (NCT) associated with the anode heel effect was evaluated and used to define a new slice-by-slice background subtraction artifact removal technique for CT PGD. Image quality was assessed for the multislice system by evaluating image noise and uniformity. The agreement in NCT for slices acquired simultaneously using the multislice detector array was also examined. Further study was performed to assess the effects of increasing x-ray tube load on the constancy of measured NCT and overall scan time. In all cases, results were compared to the single slice machine. Finally, images were collected throughout the volume of an unirradiated gel dosimeter to quantify image noise and uniformity before radiation is delivered.

RESULTS

Slice-by-slice background subtraction effectively removes the variability in NCT observed across images acquired simultaneously using the multislice scanner and is the recommended background subtraction method when using a multislice CT system. Image noise was higher for the multislice system compared to the single slice scanner, but overall image quality was comparable between the two systems. Further study showed NCT was consistent across image slices acquired simultaneously using the multislice detector array for each detector configuration of the slice thicknesses examined. In addition, the multislice system was found to eliminate variations in NCT due to increasing x-ray tube load and reduce scanning time by a factor of 4 when compared to imaging a large volume using a single slice scanner. Images acquired through an unirradiated, active gel revealed NCT varies between the top and bottom of the 1 L cylinder as well as across the diameter of the cylinder by up to 7 HU.

CONCLUSIONS

Multislice CT imaging has been evaluated for CT PGD and found to be the superior technique compared to single slice imaging in terms of the time required to complete a scan and the tube load characteristics associated with each scanning method. The implementation of multislice scanning is straightforward and expected to facilitate routine gel dosimetry measurements for complex dose distributions in modern RT centers.

Effects of refractive index mismatch in optical CT imaging of polymer gel dosimeters

PURPOSE

Proposing an image reconstruction technique, algebraic reconstruction technique-refraction correction (ART-rc). The proposed method takes care of refractive index mismatches present in gel dosimeter scanner at the boundary, and also corrects for the interior ray refraction. Polymer gel dosimeters with high dose regions have higher refractive index and optical density compared to the background medium, these changes in refractive index at high dose results in interior ray bending.

METHODS

The inclusion of the effects of refraction is an important step in reconstruction of optical density in gel dosimeters. The proposed ray tracing algorithm models the interior multiple refraction at the inhomogeneities. Jacob's ray tracing algorithm has been modified to calculate the pathlengths of the ray that traverses through the higher dose regions. The algorithm computes the length of the ray in each pixel along its path and is used as the weight matrix. Algebraic reconstruction technique and pixel based reconstruction algorithms are used for solving the reconstruction problem. The proposed method is tested with numerical phantoms for various noise levels. The experimental dosimetric results are also presented.

RESULTS

The results show that the proposed scheme ART-rc is able to reconstruct optical density inside the dosimeter better than the results obtained using filtered backprojection and conventional algebraic reconstruction approaches. The quantitative improvement using ART-rc is evaluated using gamma-index. The refraction errors due to regions of different refractive indices are discussed. The effects of modeling of interior refraction in the dose region are presented.

CONCLUSIONS

The errors propagated due to multiple refraction effects have been modeled and the improvements in reconstruction using proposed model is presented. The refractive index of the dosimeter has a mismatch with the surrounding medium (for dry air or water scanning). The algorithm reconstructs the dose profiles by estimating refractive indices of multiple inhomogeneities having different refractive indices and optical densities embedded in the dosimeter. This is achieved by tracking the path of the ray that traverses through the dosimeter. Extensive simulation studies have been carried out and results are found to be matching that of experimental results.

Radiation-induced refraction artifacts in the optical CT readout of polymer gel dosimeters

PURPOSE

The objective of this work is to demonstrate imaging artifacts that can occur during the optical computed tomography (CT) scanning of polymer gel dosimeters due to radiation-induced refractive index (RI) changes in polyacrylamide gels.

METHODS

A 1 L cylindrical polyacrylamide gel dosimeter was irradiated with 3 × 3 cm(2) square beams of 6 MV photons. A prototype fan-beam optical CT scanner was used to image the dosimeter. Investigative optical CT scans were performed to examine two types of rayline bending: (i) bending within the plane of the fan-beam and (ii) bending out the plane of the fan-beam. To address structured errors, an iterative Savitzky-Golay (ISG) filtering routine was designed to filter 2D projections in sinogram space. For comparison, 2D projections were alternatively filtered using an adaptive-mean (AM) filter.

RESULTS

In-plane rayline bending was most notably observed in optical CT projections where rays of the fan-beam confronted a sustained dose gradient that was perpendicular to their trajectory but within the fan-beam plane. These errors caused distinct streaking artifacts in image reconstructions due to the refraction of higher intensity rays toward more opaque regions of the dosimeter. Out-of-plane rayline bending was observed in slices of the dosimeter that featured dose gradients perpendicular to the plane of the fan-beam. These errors caused widespread, severe overestimations of dose in image reconstructions due to the higher-than-actual opacity that is perceived by the scanner when light is bent off of the detector array. The ISG filtering routine outperformed AM filtering for both in-plane and out-of-plane rayline errors caused by radiation-induced RI changes. For in-plane rayline errors, streaks in an irradiated region (>7 Gy) were as high as 49% for unfiltered data, 14% for AM, and 6% for ISG. For out-of-plane rayline errors, overestimations of dose in a low-dose region (∼50 cGy) were as high as 13 Gy for unfiltered data, 10 Gy for AM, and 3.1 Gy for ISG. The ISG routine also addressed unrelated artifacts that previously needed to be manually removed in sinogram space. However, the ISG routine blurred reconstructions, causing losses in spatial resolution of ∼5 mm in the plane of the fan-beam and ∼8 mm perpendicular to the fan-beam.

CONCLUSIONS

This paper reveals a new category of imaging artifacts that can affect the optical CT readout of polyacrylamide gel dosimeters. Investigative scans show that radiation-induced RI changes can cause significant rayline errors when rays confront a prolonged dose gradient that runs perpendicular to their trajectory. In fan-beam optical CT, these errors manifested in two ways: (1) distinct streaking artifacts caused by in-plane rayline bending and (2) severe overestimations of opacity caused by rays bending out of the fan-beam plane and missing the detector array. Although the ISG filtering routine mitigated these errors better than an adaptive-mean filtering routine, it caused unacceptable losses in spatial resolution.