A Prospective Study Measuring Resident and Faculty Contour Concordance: A Potential Tool for Quantitative Assessment of Residents' Performance in Contouring and Target Delineation in Radiation Oncology Residency.

PURPOSE/OBJECTIVE(S): Accurate target delineation (ie, contouring) is essential for radiation treatment planning and radiotherapy efficacy. As a result, improving the quality of target delineation is an important goal in the education of radiation oncology residents. The purpose of this study was to track the concordance of radiation oncology residents' contours with those of faculty physicians over the course of 1 year to assess for patterns. MATERIALS/METHODS: Residents in postgraduate year (PGY) levels 2 to 4 were asked to contour target volumes that were then compared to the finalized, faculty physician-approved contours. Concordance between resident and faculty physician contours was determined by calculating the Jaccard concordance index (JCI), ranging from 0, meaning no agreement, to 1, meaning complete agreement. Multivariate mixed-effect models were used to assess the association of JCI to the fixed effect of PGY level and its interactions with cancer type and other baseline characteristics. Post hoc means of JCI were compared between PGY levels after accounting for multiple comparisons using Tukey's method. RESULTS: In total, 958 structures from 314 patients collected during the 2020-2021 academic year were studied. The mean JCI was 0.77, 0.75, and 0.61 for the PGY-4, PGY-3, and PGY-2 levels, respectively. The JCI score for PGY-2 was found to be lower than those for PGY-3 and PGY-4, respectively (all P < .001). No statistically significant difference of JCI score was found between the PGY-3 and PGY-4 levels. The average JCI score was lowest (0.51) for primary head and/or neck cancers, and it was highest (0.80) for gynecologic cancers. CONCLUSIONS: Tracking and comparing the concordance of resident contours with faculty physician contours is an intriguing method of assessing resident performance in contouring and target delineation and could potentially serve as a quantitative metric, which is lacking currently, in radiation oncology resident evaluation. However, additional study is necessary before this technique can be incorporated into residency assessments.

Cleaning the dose falloff with low modulation in SBRT lung plans.

Purpose.This dosimetric study is intended to lower the modulation factor in lung SBRT plans generated in the Eclipse TPS that could replace highly modulated plans that are prone to the interplay effect.Materials and methods.Twenty clinical lung SBRT plans with high modulation factors (≥4) were replanned in Varian Eclipse TPS version 15.5 utilizing 2 mm craniocaudal and 1 mm axial block margins followed by light optimization in order to reduce modulation. A unique plan optimization methodology, which utilizes a novel shell structure (OptiForR50) for R50%optimization in addition to five consecutive concentric 5 mm shells, was utilized to control dose falloff according to RTOG 0813 and 0915 recommendations. The prescription varied from 34-54 Gy in 1-4 fractions, and the dose objectives were PTV D95%= Rx, PTV Dmax< 140% of Rx, and minimizing the modulation factor. Plan evaluation metrics included modulation factor, CIRTOG, homogeneity index (HI), R50%, D2cm, V105%, and lung V8-12.8Gy(Timmerman Constraint). A random-intercept linear mixed effects model was used with a p ≤ 0.05 threshold to test for statistical significance.Results.The retrospectively generated plans had significantly lower modulation factors (3.65 ± 0.35 versus 4.59 ± 0.54; p < 0.001), lower CIRTOG(0.97 ± 0.02 versus 1.02 ± 0.06; p = 0.001), higher HI (1.35 ± 0.06 versus 1.14 ± 0.04; p < 0.001), lower R50%(4.09 ± 0.45 versus 4.56 ± 0.56; p < 0.001), and lower lungs V8-12.8Gy(Timmerman) (4.61% ± 3.18% versus 4.92% ± 3.37%; p < 0.001). The high dose spillage V105%was borderline significantly lower (0.44% ± 0.49% versus 1.10% ± 1.64%; p = 0.051). The D2cmwas not statistically different (46.06% ± 4.01% versus 46.19% ± 2.80%; p = 0.835).Conclusion.Lung SBRT plans with significantly lower modulation factors can be generated that meet the RTOG constraints, using our planning strategy.

BEAM MATCHING EVALUATION OF TWO SIMILAR LINEAR ACCELERATORS.

To evaluate the beam-matching of two Siemens Primus medical linear accelerators (Linacs), the output factor (Sc,p), wedge factor, quality index (TPR20/10), percentage depth dose (PDD) and beam profiles were compared for 6 and 15 MV photon beams. The output factor, the PDD and the beam profile for electron beam compared for 5, 7, 8, 10 and 12 MeV electron beams. The gamma (γ) analysis of 2 mm/2% and 3 mm/3% was performed. According to the measurements, it can be said that 6 MV photon beams in all field sizes (except 4 × 4 cm2) are beam matched. For 15 MV, although the PDDs were matched in all field sizes (except 4 × 4 cm2) for both 2 mm/2% and 3 mm/3% γ criteria, beam profiles in field sizes larger than 10 × 10 cm2 for 3 mm/3% and in field sizes larger than 8 × 8 cm2 for 2 mm/2% were not matched. The electron beams in all applicator sizes (except 5 × 5 cm2 applicator) pass the acceptance γ criteria of 3 mm/3% (γ < 1). Electron beams do not fulfill beam matched in case of the acceptance γ criteria of 2 mm/2%.

Design and performance evaluation of a slit-slat camera for 2D prompt gamma imaging in proton therapy monitoring: A Monte Carlo simulation study.

PURPOSE: We investigated the design of a prompt gamma camera for real-time dose delivery verification and the partial mitigation of range uncertainties. METHODS: A slit slat (SS) camera was optimized using the trade-off between the signal-to-noise ratio and spatial resolution. Then, using the GATE Monte Carlo package, the camera performances were estimated by means of target shifts, beam position quantification, changing the camera distance from the beam, and air cavity inserting. A homogeneous PMMA phantom and the air gaps induced PMMA phantom were used. The air gaps ranged from 5 mm to 30 mm by 5 mm increments were positioned in the middle of the beam range. To reduce the simulation time, phase space scoring was used. The batch method with five realizations was used for stochastic error calculations. RESULTS: The system's detection efficiency was 1.1 × 10 - 4 PGs Emitted PGs ( 1.8 × 10 - 5 $1.1 \times {10}^{-4}\frac{{\rm PGs}}{{\rm Emitted}\ {\rm PGs}}\ (1.8 \times {10}^{-5}$ PGs/proton) for a 10 × 20 cm2 detector (source-to-collimator distance = 15.0 cm). Axial and transaxial resolutions were 23 mm and 18 mm, respectively. The SS camera estimated the range as 69.0 ± 3.4 (relative stochastic error 1-sigma is 5%) and 67.6 ± 1.8 mm (2.6%) for the real range of 67.0 mm for 107 and 108 protons of 100 MeV, respectively. Considering 160 MeV, these values are 155.5 ± 3.1 (2%) and 152.2 ± 2.0 mm (1.3%) for the real range of 152.0 mm for 107 and 108 protons, respectively. Considering phantom shift, for a 100 MeV beam, the precision of the quantification (1-sigma) in the axial and lateral phantom shift estimation is 2.6 mm and 1 mm, respectively. Accordingly, the axial and lateral quantification precisions were 1.3 mm and 1 mm for a 160 MeV beam, respectively. Furthermore, the quantification of an air gap formulated as gap d e t = 0.98 × gap real ${{\rm gap}}_{det}=0.98 \times {{\rm gap}}_{{\rm real}}$ , where gap d e t ${{\rm gap}}_{det}$ and gapreal are the estimated and real air gap, respectively. The precision of the air gap quantification is 1.6 mm (1 sigma). Moreover, 2D PG images show the trajectory of the proton beam through the phantom. CONCLUSION: The proposed slit-slat imaging systems can potentially provide a real-time, in-vivo, and non-invasive treatment monitoring method for proton therapy.

A phantom-based study and clinical implementation of brainlab's treatment planning system for radiosurgical treatments of arteriovenous malformations.

PURPOSE: Development of a simple, phantom-based methodology allowing for pilot applications for the Elements TPS cranio-vascular module and clinical implementation prior to AVM treatments. METHODS: A customized phantom was developed to be visible in MRI and CT images. High resolution digital subtraction angiograms (DSAs) and CT images of the phantom were acquired and imported into the Brainlab Elements treatment planning system. A clinical treatment plan with 5 arcs was generated in cranial vascular planning module and delivered to the phantom using a Varian TrueBeam STx Linac equipped with HD-MLCs and Brainlab ExacTrac imaging system for non-coplanar setup verification. The delivered dose was verified using a calibrated ionization chamber placed in the phantom. Upon verification of the TPS workflow, three patients with AVM who have been treated to date at our center using the Brainlab's cranial vascular module for AVM are presented here for retrospective review. RESULTS: The difference between the planed and measured dose by the ionization chamber was found to be less than 1%. Following a successful dose verification study, a clinical workflow was created. Currently, three AVM patients have been treated successfully. Clinical aspects of imaging and treatment planning consideration are presented in retrospective setting. CONCLUSIONS: Dose verification of the Brainlab Elements cranial vascular planning module for intracranial SRS treatments of AVM on Varian TrueBeam was successfully implemented using a custom-made phantom with <1% discrepancy. The Brainlab Elements' cranial vascular module was successfully implemented in clinical workflow to treat patients with AVM. This manuscript provides a guideline for clinical implementation of frameless Linac-based AVM treatment using the Brainlab Elements TPS.

A dosimetric comparative analysis of Brainlab elements and Eclipse RapidArc for spine SBRT treatment planning.

Purpose. This is a dosimetric study comparing stereotactic body radiotherapy (SBRT) plans of spine tumors using Brainlab Elements Spine planning module against Eclipse RapidArc plans. Dose conformity, dose gradient, dose fall-off, and patient-specific quality assurance (QA) metrics were evaluated. Methods:Twenty patients were immobilized in supine position using half Vac-Lok. A prescription dose of 16 Gy in a single fraction was planned for Varian TrueBeam. Conformal arc plans were generated with Pencil beam (PB), MonteCarlo (MC) in Elements, and RapidArc with Acuros XB algorithm in Eclipse using identical treatment geometry.Results. Eclipse, Elements PB, and Elements MC generated dosimetrically conformal plans having Inverse Paddick Conformity Index (IPCI) <1.3. All plans satisfied the dose constraints to target and OARs. Elements PB had a sharper gradient than Elements MC with average GI of 3.67(95% CI: 3.52-3.82) and 4.06 (95% CI: 3.93-4.20) respectively. Eclipse plans were more homogeneous with mean HI = 1.22 (95% CI: 1.20-1.23) that is lower than others. Average maximum clinical target volume (CTV) doses were higher in Elements MC with 22.31 Gy (95% CI: 21.87-22.74), while PB plans have 21.15 Gy (95% CI: 20.36-21.96), respectively. Elements MC and PB plans had lower average dose to 0.35 c.c. of spinal cord (D0.35cc) of 7.60 Gy (95% CI: 7.18-8.02) and 8.42 Gy (95% CI: 7.83-9.01). All plans had >95% points passing the gamma QA criteria at 3%/2 mm.Conclusion. All treatment plans achieved clinically acceptable target coverage >95% and meet spinal cord dose limits. Smart optimization in Brainlab Elements spine module produced dosimetrically superior plans by better spinal cord sparing.

Monte Carlo simulation and analytical calculation methods to investigate the potential of nanoparticles for INTRABEAM® IORT machine.

In the present study, Monte Carlo (MC) simulation and analytical calculation methods were used to investigate the potential of cancer treatment for the combination of IORT with nanoparticles (NPs). The Geant4 MC toolkit was used to simulate ZEISS INTRABEAM® IORT machine and its smallest applicator with 1.5 cm diameter. The dose enhancement effects (DEFs) were obtained for silver (Ag), gold (Au), bismuth (Bi), copper (Cu) and iron (Fe) spherical NPs considering different concentrations. In addition, analytical calculations were performed based on attenuation coefficient formula for sample NPs. Our MC results showed that the use of different NPs led to an increase in DEF up to 40%. Among different NPs, Au had the maximum DEF. In addition, analytical calculations revealed a significant increase, using NPs as well. Our study has suggested that the use of NPs in combination with IORT has the potential to enhance treatment outcomes.

Simultaneous respiratory motion correction and image reconstruction in 4D-multi pinhole small animal SPECT.

PURPOSE: Respiratory motion in the chest region during single photon emission computed tomography (SPECT) is a major degrading factor that reduces the accuracy of image quantification. This effect is more notable when the tumor is very small, or the spatial resolution of the imaging system is less than the respiratory motion amplitude. Small animals imaging systems with sub-millimeter spatial resolution need more attention to the respiratory motion for quantitative studies. We developed a motion-embedded four-dimensional (4D)-multi pinhole SPECT (MPS) reconstruction algorithm for respiratory motion correction. This algorithm makes full use of projection statistics for reconstruction of every individual frame. METHODS: The ROBY phantom with small tumors in liver was generated in eight different phases during one respiratory cycle. The MPS projections were modeled using a fast ray tracing method simulating an MPS acquisition. Individual frames were reconstructed and used for motion estimation. The Demons non-rigid registration algorithm was used to calculate deformation vector fields (DVFs) for simultaneous motion correction and image reconstruction. A motion-embedded 4D-MPS method was used to reconstruct images using all the projections and corresponding DVFs, simultaneously. The 4D-MPS reconstructed images were compared to the low-count single frame (LCSF) reconstructed image, the three-dimensional (3D)-MPS images reconstructed using individual frames, and post reconstruction registration (PRR) that aligns all individual phases to a reference frame using Demons-derived DVFs. The tumor volume relative error (TVE), tumor contrast relative error (TCE), and dice index (DI) for 2, 3, and 4 mm liver were calculated and compared for different reconstruction methods. RESULTS: For the 4D-MPS reconstruction method, TVE was reduced and DI was higher compared to PRR, 3D-MPS, and LCSF. The extent of the improvement was higher for the small tumor size (i.e. 2 mm). For the biggest tumor in contrast 3 (i.e. 4 mm) TVE for 4D-MPS, PRR, 3D-MPS and, LCSF were 1.33%, 8%, 8%, and 14.67%, respectively. CONCLUSIONS: The results suggest that motion-embedded 4D-MPS method is an effective and practical way for respiratory motion correction. It reconstructs high quality gated frames while using all projection data to reconstruct each frame.

Evaluation of normal lung tissue complication probability in gated and conventional radiotherapy using the 4D XCAT digital phantom.

PURPOSE: The present study was conducted to investigate normal lung tissue complication probability in gated and conventional radiotherapy (RT) as a function of diaphragm motion, lesion size, and its location using 4D-XCAT digital phantom in a simulation study. MATERIALS AND METHODS: Different time series of 3D-CT images were generated using the 4D-XCAT digital phantom. The binary data obtained from this phantom were then converted to the digital imaging and communication in medicine (DICOM) format using an in-house MATLAB-based program to be compatible with our treatment planning system (TPS). The 3D-TPS with superposition computational algorithm was used to generate conventional and gated plans. Treatment plans were generated for 36 different XCAT phantom configurations. These included four diaphragm motions of 20, 25, 30 and 35 mm, three lesion sizes of 3, 4, and 5 cm in diameter and each tumor was placed in four different lung locations (right lower lobe, right upper lobe, left lower lobe and left upper lobe). The complication of normal lung tissue was assessed in terms of mean lung dose (MLD), the lung volume receiving ≥20 Gy (V20), and normal tissue complication probability (NTCP). RESULTS: The results showed that the gated RT yields superior outcomes in terms of normal tissue complication compared to the conventional RT. For all cases, the gated radiation therapy technique reduced the mean dose, V20, and NTCP of lung tissue by up to 5.53 Gy, 13.38%, and 23.89%, respectively. CONCLUSIONS: The results of this study showed that the gated RT provides significant advantages in terms of the normal lung tissue complication, compared to the conventional RT, especially for the lesions near the diaphragm.

Iterative reconstruction with boundary detection for carbon ion computed tomography.

In heavy ion radiation therapy, improving the accuracy in range prediction of the ions inside the patient's body has become essential. Accurate localization of the Bragg peak provides greater conformity of the tumor while sparing healthy tissues. We investigated the use of carbon ions directly for computed tomography (carbon CT) to create the relative stopping power map of a patient's body. The Geant4 toolkit was used to perform a Monte Carlo simulation of the carbon ion trajectories, to study their lateral and angular deflections and the most likely paths, using a water phantom. Geant4 was used to create carbonCT projections of a contrast and spatial resolution phantom, with a cone beam of 430 MeV/u carbon ions. The contrast phantom consisted of cranial bone, lung material, and PMMA inserts while the spatial resolution phantom contained bone and lung material inserts with line pair (lp) densities ranging from 1.67 lp cm-1 through 5 lp cm-1. First, the positions of each carbon ion on the rear and front trackers were used for an approximate reconstruction of the phantom. The phantom boundary was extracted from this approximate reconstruction, by using the position as well as angle information from the four tracking detectors, resulting in the entry and exit locations of the individual ions on the phantom surface. Subsequent reconstruction was performed by the iterative algebraic reconstruction technique coupled with total variation minimization (ART-TV) assuming straight line trajectories for the ions inside the phantom. The influence of number of projections was studied with reconstruction from five different sets of projections: 15, 30, 45, 60 and 90. Additionally, the effect of number of ions on the image quality was investigated by reducing the number of ions/projection while keeping the total number of projections at 60. An estimation of carbon ion range using the carbonCT image resulted in improved range prediction compared to the range calculated using a calibration curve.

Quantitative 4D-PET reconstruction for small animal using SMEIR-reconstructed 4D-CBCT.

Respiratory motions in small animals PET cause image degradation during reconstruction. This work aims to develop a motion compensated 4D-PET reconstruction method using accurate motion corrections and attenuation corrections from 4D-CBCT images reconstructed using a simultaneous motion estimation and image reconstruction (SMEIR) method. Projections of 4D-CBCT were calculated using a ray-tracing method on a digital 4D rat phantom, and list-mode data of 4D-PET with matched respiratory phases were simulated using the GATE Monte Carlo package. The respiratory rate was set at 1.0 second per cycle with 10 phases of 30 projection images each. 4D-CBCT images were reconstructed using the SMEIR method and motion information and linear attenuation from 4D-CBCT were subsequently used for motion compensated 4D-PET reconstruction and attenuation corrections. We quantitatively evaluate the reconstructed 4D-PET using the errors of tumor volume and standard uptake values of tumors with different sizes. The tumor motion was successfully reconstructed and showed good agreement with the original phantom. The proposed method reduced tumor volume errors and standard uptake value errors. For tumor diameters of 3.0, 4.5, and 6.0 mm, the tumor volume errors are 32.5%, 29.2% and 19.4% respectively with motion compensation and 45.1%, 37.5% and 20.2% respectively without compensation.

Internal dosimetry for radioembolization therapy with Yttrium-90 microspheres.

The absorbed doses in the liver and adjacent viscera in Yttrium-90 radioembolization therapy for metastatic liver lesions are not well-documented. We sought for a clinically practical way to determine the dosimetry of this advent treatment. Six different female XCAT BMIs and seven different male XCAT BMIs were generated. Using Monte Carlo GATE code simulation, the total of 100MBq 90 Y was deposited uniformly in the source organ, liver. Self-irradiation and absorbed doses in lung, kidney and bone marrow were calculated. The mean energy of Yittrium-90 (i.e., 0.937 MeV) was used. The S-values and equivalent doses in target organs were estimated. The dose absorbed in the liver was between 84 and 53 Gy and below the target of 80 to 150 Gy. The absorbed dose in the bone marrow, lungs, and kidneys are very low and below 0.1 , 0.4, and 0.5 Gy respectively. Our study indicates that larger activities than the conventional dose of 3 GBq may be both required and safe. Further confirmations in clinical settings are needed.

The importance of BMI in dosimetry of 153Sm-EDTMP bone pain palliation therapy: A Monte Carlo study.

Using digital phantoms as an atlas compared to acquiring CT data for internal radionuclide dosimetry decreases patient overall radiation dose and reduces the required analysis effort and time for organ segmentation. The drawback is that the phantom may not match exactly with the patient. We assessed the effect of varying BMIs on dosimetry results for a bone pain palliation agent, 153Sm-EDTMP. The simulation was done using the GATE Monte Carlo code. Female XCAT phantoms with the following different BMIs were employed: 18.6, 20.8, 22.1, 26.8, 30.3 and 34.7kg/m2. S-factors (mGy/MBq.s) and SAFs (kg-1) were calculated for the dosimetry of the radiation from major source organs including spine, ribs, kidney and bladder into different target organs as well as whole body dosimetry from spine. The differences in dose estimates from different phantoms compared to those from the phantom with BMI of 26.8kg/m2 as the reference, were calculated for both gamma and beta radiations. The relative differences (RD) of the S-factors or SAFs from the values of reference phantom were calculated. RDs greater than 10% and 100% were frequent in radiations to organs for photon and beta particles, respectively. The relative differences in whole body SAFs from the reference phantom were 15.4%, 7%, 4.2%, -9.8% and -1.4% for BMIs of 18.6, 20.8, 22.1, 30.3 and 34.7kg/m2, respectively. The differences in whole body S-factors for the phantoms with BMIs of 18.6, 20.8, 22.1, 30.3 and 34.7kg/m2 were 39.5%, 19.4%, 8.8%, -7.9% and -4.3%, respectively. The dosimetry of the gamma photons and beta particles changes substantially with the use of phantoms with different BMIs. The change in S-factors is important for dose calculation and can change the prescribed therapeutic dose of 153Sm-EDTMP. Thus a phantom with BMI better matched to the patient is suggested for therapeutic purposes where dose estimates closer to those in the actual patient are required.

Attenuation correction in 4D-PET using a single-phase attenuation map and rigidity-adaptive deformable registration.

PURPOSE: Four-dimensional positron emission tomography (4D-PET) imaging is a potential solution to the respiratory motion effect in the thoracic region. Computed tomography (CT)-based attenuation correction (AC) is an essential step toward quantitative imaging for PET. However, due to the temporal difference between 4D-PET and a single attenuation map from CT, typically available in routine clinical scanning, motion artifacts are observed in the attenuation-corrected PET images, leading to errors in tumor shape and uptake. We introduced a practical method to align single-phase CT with all other 4D-PET phases for AC. METHODS: A penalized non-rigid Demons registration between individual 4D-PET frames without AC provides the motion vectors to be used for warping single-phase attenuation map. The non-rigid Demons registration was used to derive deformation vector fields (DVFs) between PET matched with the CT phase and other 4D-PET images. While attenuated PET images provide useful data for organ borders such as those of the lung and the liver, tumors cannot be distinguished from the background due to loss of contrast. To preserve the tumor shape in different phases, an ROI-covering tumor was excluded from nonrigid transformation. Instead the mean DVF of the central region of the tumor was assigned to all voxels in the ROI. This process mimics a rigid transformation of the tumor along with a nonrigid transformation of other organs. A 4D-XCAT phantom with spherical lung tumors, with diameters ranging from 10 to 40 mm, was used to evaluate the algorithm. The performance of the proposed hybrid method for attenuation map estimation was compared to (a) the Demons nonrigid registration only and (b) a single attenuation map based on quantitative parameters in individual PET frames. RESULTS: Motion-related artifacts were significantly reduced in the attenuation-corrected 4D-PET images. When a single attenuation map was used for all individual PET frames, the normalized root-mean-square error (NRMSE) values in tumor region were 49.3% (STD: 8.3%), 50.5% (STD: 9.3%), 51.8% (STD: 10.8%) and 51.5% (STD: 12.1%) for 10-mm, 20-mm, 30-mm, and 40-mm tumors, respectively. These errors were reduced to 11.9% (STD: 2.9%), 13.6% (STD: 3.9%), 13.8% (STD: 4.8%), and 16.7% (STD: 9.3%) by our proposed method for deforming the attenuation map. The relative errors in total lesion glycolysis (TLG) values were -0.25% (STD: 2.87%) and 3.19% (STD: 2.35%) for 30-mm and 40-mm tumors, respectively, in proposed method. The corresponding values for Demons method were 25.22% (STD: 14.79%) and 18.42% (STD: 7.06%). Our proposed hybrid method outperforms the Demons method especially for larger tumors. For tumors smaller than 20 mm, nonrigid transformation could also provide quantitative results. CONCLUSION: Although non-AC 4D-PET frames include insignificant anatomical information, they are still useful to estimate the DVFs to align the attenuation map for accurate AC. The proposed hybrid method can recover the AC-related artifacts and provide quantitative AC-PET images.

Practical Nuclear Medicine and Utility of Phantoms for Internal Dosimetry: XCAT Compared with Zubal.

PURPOSE: The absorbed doses for two radioisotopes, 99mTc and 131I, between previously validated Zubal phantom and the recently developed XCAT phantom were compared. MATERIALS AND METHODS: GATE Monte Carlo code was used to simulate the statistical process of radiation. A XCAT phantom with voxel and matrix sizes similar to a standard Zubal phantom was generated. According to Medical International Radiation Dose formalism, specific absorbed fraction (SAF) values for photons and S-factors for beta particles were tabulated. The amounts of absorbed doses were calculated and compared in different organs. RESULTS: The differences of gamma radiation doses, SAFs, between Zubal and XCAT are >50% in adrenal from adrenal, pancreas from pancreas and thyroid from thyroid, in lung from kidney, kidneys from lungs and in kidneys from thyroid and thyroid from kidneys. The beta radiation doses differences between Zubal and XCAT are >50% in thyroid from thyroid, bladder from bladder, kidney from kidney, liver from bladder, thyroid from bladder and kidney from thyroid. The size and distances of the organs were different between XCAT and Zubal phantoms. Denoted differences of SAF and S-factors correspond to the different organ geometries in phantoms. CONCLUSION: The results of absorbed doses in Zubal and XCAT phantoms are different. The variations prohibit easy comparison or interchangeability of dosimetry between these phantoms.

Task Equivalence for Model and Human-Observer Comparisons in SPECT Localization Studies.

While mathematical model observers are intended for efficient assessment of medical imaging systems, their findings should be relevant for human observers as the primary clinical end users. We have investigated whether pursuing equivalence between the model and human-observer tasks can help ensure this goal. A localization ROC (LROC) study tested prostate lesion detection in simulated In-111 SPECT imaging with anthropomorphic phantoms. The test images were 2D slices extracted from reconstructed volumes. The iterative OSEM reconstruction method was used with Gaussian postsmoothing. Variations in the number of iterations and the level of postfiltering defined the test strategies in the study. Human-observer performance was compared with that of a visual-search (VS) observer, a scanning channelized Hotelling observer, and a scanning nonprewhitening (CNPW) observer. These model observers were applied with precise information about the target regions of interest (ROIs). ROI knowledge was a study variable for the human observers. In one study format, the humans read the SPECT image alone. With a dual-modality format, the SPECT image was presented alongside an anatomical image slice extracted from the density map of the phantom. Performance was scored by area under the LROC curve. The human observers performed significantly better with the dual-modality format, and correlation with the model observers was also improved. Given the human-observer data from the SPECT study format, the Pearson correlation coefficients for the model observers were 0.58 (VS), -0.12 (CH), and -0.23 (CNPW). The respective coefficients based on the human-observer data from the dual-modality study were 0.72, 0.27, and -0.11. These results point towards the continued development of the VS observer for enhancing task equivalence in model-observer studies.

Respiratory motion correction in 4D-PET by simultaneous motion estimation and image reconstruction (SMEIR).

In conventional 4D positron emission tomography (4D-PET), images from different frames are reconstructed individually and aligned by registration methods. Two issues that arise with this approach are as follows: (1) the reconstruction algorithms do not make full use of projection statistics; and (2) the registration between noisy images can result in poor alignment. In this study, we investigated the use of simultaneous motion estimation and image reconstruction (SMEIR) methods for motion estimation/correction in 4D-PET. A modified ordered-subset expectation maximization algorithm coupled with total variation minimization (OSEM-TV) was used to obtain a primary motion-compensated PET (pmc-PET) from all projection data, using Demons derived deformation vector fields (DVFs) as initial motion vectors. A motion model update was performed to obtain an optimal set of DVFs in the pmc-PET and other phases, by matching the forward projection of the deformed pmc-PET with measured projections from other phases. The OSEM-TV image reconstruction was repeated using updated DVFs, and new DVFs were estimated based on updated images. A 4D-XCAT phantom with typical FDG biodistribution was generated to evaluate the performance of the SMEIR algorithm in lung and liver tumors with different contrasts and different diameters (10-40 mm). The image quality of the 4D-PET was greatly improved by the SMEIR algorithm. When all projections were used to reconstruct 3D-PET without motion compensation, motion blurring artifacts were present, leading up to 150% tumor size overestimation and significant quantitative errors, including 50% underestimation of tumor contrast and 59% underestimation of tumor uptake. Errors were reduced to less than 10% in most images by using the SMEIR algorithm, showing its potential in motion estimation/correction in 4D-PET.

A New Approach for Scatter Removal and Attenuation Compensation from SPECT/CT Images.

OBJECTIVE(S): In SPECT, the sinogram contains scatter and lack of attenuated counts that degrade the reconstructed image quality and quantity. Many techniques for attenuation and scatter correction have been proposed. An acceptable method of correction is to incorporate effects into an iterative statistical reconstruction. Here, we propose new Maximum Likelihood Expectation Maximization (MLEM) formula to correct scattering and attenuating photons during reconstruction. MATERIALS AND METHODS: In this work, scatters are estimated through Klein-Nishina formula in all iterations and CT images are used for accurate attenuation correction. Reconstructed images resulted from different MLEM reconstruction formula have been compared considering profile agreement, contrast, mean square error, signal-to-noise ratio, contrast-to-noise ratio and computational time. RESULTS: The proposed formula has a good profile agreement, increased contrast, signal-to-noise (SNR) & contrast-to-noise ratio (CNR), computational time and decreased mean square error (MSE) compared with uncorrected images and/or images from conventional formula. CONCLUSION: In conclusion, by applying the proposed formula we were able to correct attenuation and scatter via MLEM and improve the image quality, which is a necessary step for both qualitative and quantitative SPECT images.

The influence of resolution recovery by using collimator detector response during 3D OSEM image reconstruction on (99m)Tc-ECD brain SPET images.

Partial volume effect, due to the poor spatial resolution of single photon emission tomography (SPET), significantly restricts the absolute quantification of the regional brain uptake and limits the accuracy of the absolute measurement of blood flow. In this study the importance of compensation for the collimator-detector response (CDR) in the technetium-99m ethyl cysteinate dimer ((99m)Tc-ECD) brain SPET was assessed, by incorporating system response in the ordered-subsets expectation maximization (OSEM) reconstruction algorithm. By placing a point source of (99m)Tc at different distances from the face of the collimator, CDR were found and modeled using Gaussian functions. A fillable slice of the brain phantom was designed and filled by (99m)Tc. Projections acquired from the phantom and also 4 patients who underwent the (99m)Tc-ECD brain SPET were used in this study. To reconstruct the images, 3D OSEM algorithm was used. System blurring functions were modeled, during the reconstruction in both projection and backprojection steps. Our results were compared with the conventional resolution recovery using Metz filter in filtered backprojection (FBP). Visual inspection of the images was performed by six nuclear medicine specialists. Quantitative analysis was also studied by calculating the contrast and the count density of the reconstructed images. For the phantom images, background counts and noise were decreased by 3D OSEM compared to the FBP-Metz method. Quantitatively, the ratio of the counts of the occupied hot region to that of the cold region of the reconstructed by FBP-Metz images was 1.14. This value was decreased from 1.12 to 0.86 for 3D OSEM of 2 and 30 iterations respectively. The reference value was 0.85 for the planar image. For clinical images, hot to cold regions (grey to white matter), the count ratio was increased from 1.44 in FBP-Metz to 3.2 and 4 in 3D OSEM with 10 and 20 iterations respectively. Based on the interpretability of images, the best scores (3.79±0.51) by the physicians were given to the images reconstructed by 3D OSEM and 10 iterations. This value was 0.63±0.77 for FBP-Metz images. In conclusion, by incorporating the distance dependent CDR during 3D OSEM, it was possible to reconstruct the brain images with much higher resolution and contrast as compared to the conventional resolution recovery method, which used FBP-Metz. It was however important to make a trade-off between noise and resolution by determining an optimum iterations number.

Quantification and reduction of the collimator-detector response effect in SPECT by applying a system model during iterative image reconstruction: a simulation study.

BACKGROUND: Detector blurring and non-ideal collimation decrease the spatial resolution of the single-photon emission computed tomography (SPECT) images. Iterative reconstruction algorithms such as ordered subsets expectation maximization (OSEM) can incorporate degrading factors during reconstruction. We investigated the quantitative errors associated with poor SPECT resolution and evaluated the importance of two-dimensional (2D) and three-dimensional (3D) resolution recovery by modelling system response during iterative image reconstruction. METHODS: Different phantoms consisted of the NURBS-based cardiac-torso (NCAT) liver phantom with small tumors, the Zubal brain phantom and the NCAT heart phantom were used in this study. Monte Carlo simulation was used to create SPECT projections. Gaussian functions were used to model collimator detector response (CDR). Modeled CDRs were applied during OSEM. Both noise-free and noisy projections were created. RESULTS: Even with noise-free projections, conventional OSEM algorithm provided limited quantitative accuracy compared to both 2D and 3D resolution recovery. The 3D implementation of resolution recovery, however, yielded superior results compared to its 2D implementation. For the liver phantom, the ability to distinguish small tumors in both transverse and axial planes was improved. For the brain phantom, gray to white matter activity ratio was increased from 3.14 ± 0.04 in simple OSEM to 3.84 ± 0.06 in 3D OSEM. For the NCAT heart phantom, 3D resolution recovery, results in images with thinner wall and higher contrast for different noise levels. CONCLUSION: There are considerable quantitative errors associated with CDR, especially when the size of the target is comparable with the spatial resolution of the system. Between different reconstruction algorithms, 3D OSEM that consider the 3D nature of CDR, improve both the visual quality and the quantitative accuracy of any SPECT studies.