Dose calculation accuracy of different image value to density tables for cone-beam CT planning in head & neck and pelvic localizations

PurposeWe aimed to identify the most accurate combination of phantom and protocol for image value to density table (IVDT) on volume-modulated arc therapy (VMAT) dose calculation based on kV-Cone-beam CT imaging, for head and neck (H&N) and pelvic localizations.MethodsThree phantoms (Catphan^®600, CIRS^®062M (inner phantom for head and outer phantom for body), and TomoTherapy^® “Cheese” phantom) were used to create IVDT curves of CBCT systems with two different CBCT protocols (Standard-dose Head and Standard Pelvis). Hounsfield Unit (HU) time stability and repeatability for a single On-Board-Imager (OBI) and compatibility of two distinct devices were assessed with Catphan^®600. Images from the anthropomorphic phantom CIRS ATOM^® for both CT and CBCT modalities were used for VMAT dose calculation from different IVDT curves. Dosimetric indices from CT and CBCT imaging were compared.ResultsIVDT curves from CBCT images were highly different depending on phantom used (up to 1000 HU for high densities) and protocol applied (up to 200 HU for high densities). HU time stability was verified over seven weeks. A maximum difference of 3% on the dose calculation indices studied was found between CT and CBCT VMAT dose calculation across the two localizations using appropriate IVDT curves. One IVDT curve per localization can be established with a bi-monthly verification of IVDT-CBCT.ConclusionsThe IVDT-CBCT_CIRS-Head phantom with the Standard-dose Head protocol was the most accurate combination for dose calculation on H&N CBCT images. For pelvic localizations, the IVDT-CBCT_Cheese established with the Standard Pelvis protocol provided the best accuracy.     

Accurate proton treatment planning for pencil beam crossing titanium fixation implants

PurposeTo present a planning strategy for proton pencil-beam scanning when titanium implants need to be crossed by the beam.MethodsWe addressed three issues: the implementation of a CT calibration curve to assign to titanium the correct stopping power; the effect of artefacts on CT images and their reduction by a dedicated algorithm; the differences in dose computation depending on the dose engine, pencil-beam vs Monte-Carlo algorithms. We performed measurement tests on a simple cylinder phantom and on a real implant. These phantoms were irradiated with three geometries (single spots, uniform mono-energetic layer and uniform box), measuring the exit dose either by radio-chromic film or multi-layer ionization chamber. The procedure was then applied on two patients treated for chordoma.ResultsWe had to set in the calibration curve a mass density equal to 4.37 g/cm³ to saturated Hounsfield Units, in order to have the correct stopping power assigned to titanium in TPS. CT artefact reduction algorithm allowed a better reconstruction of the shape and size of the implant. Monte-Carlo resulted accurate in computing the dose distribution whereas the pencil-beam algorithm failed due to sharp density interfaces between titanium and the surrounding material. Finally, the treatment plans obtained on two patients showed the impact of the dose engine algorithm, with 10–20% differences between pencil-beam and Monte-Carlo in small regions distally to the titanium screws.ConclusionThe described combination of CT calibration, artefacts reduction and Monte-Carlo computation provides a reliable methodology to compute dose in patients with titanium implants.     

Novel Monte Carlo dose calculation algorithm for robotic radiosurgery with multi leaf collimator: Dosimetric evaluation

PurposeAt introduction in 2014, dose calculation for the first MLC on a robotic SRS/SBRT platform was limited to a correction-based Finite-Size Pencil Beam (FSPB) algorithm. We report on the dosimetric accuracy of a novel Monte Carlo (MC) dose calculation algorithm for this MLC, included in the Precision™ treatment planning system.MethodsA phantom was built of one slab (5.0 cm) of lung-equivalent material (0.09…0.29 g/cc) enclosed by 3.5 cm (above) and 5 cm (below) slabs of solid water (1.045 g/cc). This was irradiated using rectangular (15.4 × 15.4 mm² to 53.8 × 53.7 mm²) and two irregular MLC-fields. Radiochromic film (EBT3) was positioned perpendicular to the slabs and parallel to the beam. Calculated dose distributions were compared to film measurements using line scans and 2D gamma analysis.ResultsMeasured and MC calculated percent depth dose curves showed a characteristic dose drop within the low-density region, which was not correctly reproduced by FSPB. Superior average gamma pass rates (2%/1 mm) were found for MC (91.2 ± 1.5%) compared to FSPB (55.4 ± 2.7%). However, MC calculations exhibited localized anomalies at mass density transitions around 0.15 g/cc, which were traced to a simplification in electron transport. Absence of these anomalies was confirmed in a modified build of the MC engine, which increased gamma pass rates to 96.6 ± 1.2%.ConclusionsThe novel MC algorithm greatly improves dosimetric accuracy in heterogeneous tissue, potentially expanding the clinical use of robotic radiosurgery with MLC. In-depth, independent validation is paramount to identify and reduce the residual uncertainties in any software solution.     

Development of a deformable phantom for experimental verification of 4D Monte Carlo simulations in a deforming anatomy

PurposeTo verify the accuracy of 4D Monte Carlo (MC) simulations, using the 4DdefDOSXYZnrc user code, in a deforming anatomy. We developed a tissue-equivalent and reproducible deformable lung phantom and evaluated 4D simulations of delivered dose to the phantom by comparing calculations against measurements.MethodsA novel deformable phantom consisting of flexible foam, emulating lung tissue, inside a Lucite external body was constructed. A removable plug, containing an elastic tumor that can hold film and other dosimeters, was inserted in the phantom. Point dose and position measurements were performed inside and outside the tumor using RADPOS 4D dosimetry system. The phantom was irradiated on an Elekta Infinity linac in both stationary and moving states. The dose delivery was simulated using delivery log files and the phantom motion recorded with RADPOS.ResultsReproducibility of the phantom motion was determined to be within 1 mm. The phantom motion presented realistic features like hysteresis. MC calculations and measurements agreed within 2% at the center of tumor. Outside the tumor agreements were better than 5% which were within the positional/dose reading uncertainties at the measurement points. More than 94% of dose points from MC simulations agreed within 2%/2 mm compared to film measurements.ConclusionThe deformable lung phantom presented realistic and reproducible motion characteristics and its use for verification of 4D dose calculations was demonstrated. Our 4DMC method is capable of accurate calculations of the realistic dose delivered to a moving and deforming anatomy during static and dynamic beam delivery techniques.     

Measuring the dose in bone for spine stereotactic body radiotherapy

PurposeCurrent quality assurance of radiotherapy involving bony regions generally utilises homogeneous phantoms and dose calculations, ignoring the challenges of heterogeneities with dosimetry problems likely occurring around bone. Anthropomorphic phantoms with synthetic bony materials enable realistic end-to-end testing in clinical scenarios. This work reports on measurements and calculated corrections required to directly report dose in bony materials in the context of comprehensive end-to-end dosimetry audit measurements (63 plans, 6 planning systems).Materials and methodsRadiochromic film and microDiamond measurements were performed in an anthropomorphic spine phantom containing bone equivalent materials. Medium dependent correction factors, k_med, were established using 6 MV and 10 MV Linear Accelerator Monte Carlo simulations to account for the detectors being calibrated in water, but measuring in regions of bony material. Both cortical and trabecular bony material were investigated for verification of dose calculations in dose-to-medium (D_m,m) and dose-to-water (D_w,w) scenarios.ResultsFor D_m,m calculations, modelled correction factors for cortical and trabecular bone in film measurements, and for trabecular bone in microDiamond measurements were 0.875(±0.1%), 0.953(±0.3%) and 0.962(±0.4%), respectively. For D_w,w calculations, the corrections were 0.920(±0.1%), 0.982(±0.3%) and 0.993(±0.4%), respectively. In the audit, application of the correction factors improves the mean agreement between treatment plans and measured microDiamond dose from −2.4%(±3.9%) to 0.4%(±3.7%).ConclusionMonte Carlo simulations provide a method for correcting the dose measured in bony materials allowing more accurate comparison with treatment planning system doses. In verification measurements, algorithm specific correction factors should be applied to account for variations in bony material for calculations based on D_m,m and D_w,w.     

Simplified patient-specific renal dosimetry in 177Lu therapy: a proof of concept

PurposeThe aim of this proof-of-concept study is to propose a simplified personalized kidney dosimetry procedure in ¹⁷⁷Lu peptide receptor radionuclide therapy (PRRT) for neuroendocrine tumors and metastatic prostate cancer. It relies on a single quantitative SPECT/CT acquisition and multiple radiometric measurements executed with a collimated external probe, properly directed on kidneys.MethodsWe conducted a phantom study involving external count-rate measurements in an abdominal phantom setup filled with activity concentrations of ^99mTc, reproducing patient-relevant organ effective half-lives occurring in ¹⁷⁷Lu PRRT. GATE Monte Carlo (MC) simulations of the experiment, using ^99mTc and ¹⁷⁷Lu as sources, were performed. Furthermore, we tested this method via MC on a clinical case of ¹⁷⁷Lu-DOTATATE PRRT with SPECT/CT images at three time points (2, 20 and 70 hrs), comparing a simplified kidney dosimetry, employing a single SPECT/CT and probe measurements at three time points, with the complete MC dosimetry.ResultsThe experimentally estimated kidney half-life with background subtraction applied was compatible within 3% with the expected value. The MC simulations of the phantom study, both with ^99mTc and ¹⁷⁷Lu, confirmed a similar level of accuracy. Concerning the clinical case, the simplified dosimetric method led to a kidney dose estimation compatible with the complete MC dosimetry within 6%, 12% and 2%, using respectively the SPECT/CT at 2, 20 and 70 hrs.ConclusionsThe proposed simplified procedure provided a satisfactory accuracy and would reduce the imaging required to derive the kidney absorbed dose to a unique quantitative SPECT/CT, with consequent benefits in terms of clinic workflows and patient comfort.     

Combination of dual-energy computed tomography and iterative metal artefact reduction to increase general quality of imaging for radiotherapy patients with high dense materials. Phantom study

PurposeTo evaluate the use of pseudo-monoenergetic reconstructions (PMR) from dual-energy computed tomography, combined with the iterative metal artefact reduction (iMAR) method.MethodsPseudo-monoenergetic CT images were obtained using the dual-energy mode on the Siemens Somatom Definition AS scanner. A range of PMR combinations (70–130 keV) were used with and without iMAR. A Virtual Water™ phantom was used for quantitative assessment of error in the presence of high density materials: titanium, alloys 330 and 600. The absolute values of CT number differences (AD) and normalised standard deviations (NSD) were calculated for different phantom positions. Image quality was assessed using an anthropomorphic pelvic phantom with an embedded hip prosthesis. Image quality was scored blindly by five observers.ResultsAD and NSD values revealed differences in CT number errors between tested sets. AD and NSD were reduced in the vicinity of metal for images with iMAR (p < 0.001 for AD/NSD). For ROIs away from metal, with and without iMAR, 70 keV PMR and pCT AD values were lower than for the other reconstructions (p = 0.039). Similarly, iMAR NSD values measured away from metal were lower for 130 keV and 70 keV PMR (p = 0.002). Image quality scores were higher for 70 keV and 130 keV PMR with iMAR (p = 0.034).ConclusionThe use of 70 keV PMR with iMAR allows for significant metal artefact reduction and low CT number errors observed in the vicinity of dense materials. It is therefore an attractive alternative to high keV imaging when imaging patients with metallic implants, especially in the context of radiotherapy planning.     

Machine log file-based dose verification using novel iterative CBCT reconstruction algorithm in commercial software during volumetric modulated arc therapy for prostate cancer patients

PurposeTo evaluate the utility of the use of iterative cone-beam computed tomography (CBCT) for machine log file-based dose verification during volumetric modulated arc therapy (VMAT) for prostate cancer patients.MethodsAll CBCT acquisition data were used to reconstruct images with the Feldkamp-Davis-Kress algorithm (FDK-CBCT) and the novel iterative algorithm (iCBCT). The Hounsfield unit (HU)-electron density curves for CBCT images were created using the Advanced Electron Density Phantom. The I’mRT and anthropomorphic phantoms were irradiated with VMAT after CBCT registration. Subsequently, fourteen prostate cancer patients received VMAT after CBCT registration. Machine log files and both CBCT images were exported to the PerFRACTION software, and a 3D patient dose was reconstructed. Mean dose for planning target volume (PTV), the bladder, and rectum and the 3D gamma analysis were evaluated.ResultsFor the phantom studies, the variation of HU values was observed at the central position surrounding the bones in FDK-CBCT. There were almost no changes in the difference of doses at the isocenter between measurement and reconstructed dose for planning CT (pCT), FDK-CBCT, and iCBCT. Mean dose differences of PTV, rectum, and bladder between iCBCT and pCT were approximately 2% lower than those between FDK-CBCT and pCT. For the clinical study, average gamma analysis for 2%/2 mm was 98.22% ± 1.07 and 98.81% ± 1.25% in FDK-CBCT and iCBCT, respectively.ConclusionsA similar machine log file-based dose verification accuracy is obtained for FDK-CBCT and iCBCT during VMAT for prostate cancer patients.     

Experimental benchmarking of RayStation proton dose calculation algorithms inside and outside the target region in heterogeneous phantom geometries

PurposeThe aim of the presented study was to complement existing literature on benchmarking proton dose by comparing dose calculations with experimental measurements in heterogeneous phantom. Points of interest inside and outside the target were considered to quantify the magnitude of calculation uncertainties in current and previous proton therapy practice that might especially have an impact on the dose in organs at risk (OARs).MethodsThe RayStation treatment planning system (RaySearch Laboratories), offering two dose calculation algorithms for pencil beam scanning in proton therapy, i.e., Pencil Beam (PB) and Monte Carlo (MC), was utilized. Treatment plans for a target located behind the interface of the heterogeneous tissues were generated. Dose measurements within and behind the target were performed in a water phantom with embedded slabs of various tissue equivalent materials and 24 PinPoint ionization chambers (PTW). In total 12 test configurations encompassing two different target depths, oblique beam incidence of 30 degrees and range shifter, were considered.ResultsPB and MC calculated doses agreed equally well with the measurements for all test geometries within the target, including the range shifter (mean dose differences ± 3%). Outside the target, the maximum dose difference of 9% (19%) was observed for MC (PB) for the oblique beam incidence and inserted range shifter.ConclusionThe accuracy of MC dose algorithm was superior compared to the PB algorithm, especially outside the target volumes. MC based dose calculation should therefore be preferred in treatment scenarios with heterogeneities, especially to reduce clinically relevant uncertainties for OARs.     

能谱CT优化胃肿瘤扫描辐射剂量对肾上腺嗜铬细胞瘤的诊断价值

目的:探讨能谱CT优化胃肿瘤扫描辐射剂量对肾上腺嗜铬细胞瘤的诊断价值。方法:采用回顾性、抽样、随机研究方法选择2012年9月到2017年2月在我院诊治的肾上腺嗜铬细胞瘤患者59例作为研究对象,所有患者都给予常规CT扫描与能谱CT优化胃肿瘤扫描,记录和比较辐射剂量与图像质量。结果:所有病例包膜均完整,边缘清楚,肿瘤内见单发或多发低密度区,肿瘤实质区呈不均匀显著强化。常规CT与能谱CT的图像质量主观评分分别为3.89±0.45分和4.54±0.34分;常规CT与能谱CT图像的胃肿瘤CT值分别为31.94±6.39HU和35.29±5.19HU,对比都有显著差异(P0.05)。能谱CT图像的膀胱和皮下脂肪图像噪声值都显著低于常规CT图像,对比差异都有统计学意义(P0.05);能谱CT扫描的CTDIvol和DLP分别为12.39±3.48mGy和624.10±39.19mGy.cm,都显著低于常规CT扫描的14.09±4.13mGy和653.92±56.29mGy.cm(P0.05)。结论:能谱CT优化胃肿瘤扫描在肾上腺嗜铬细胞瘤诊断中的应用能有效减少辐射剂量与图像噪声,提高图像CT值与主观质量,临床应用价值更高。     

3D-printed patient-specific pelvis phantom for dosimetry measurements for prostate stereotactic radiotherapy with dominant intraprostatic lesion boost

AimDeveloping and assessing the feasibility of using a three-dimensional (3D) printed patient-specific anthropomorphic pelvis phantom for dose calculation and verification for stereotactic ablative radiation therapy (SABR) with dose escalation to the dominant intraprostatic lesions.Material and methodsA 3D-printed pelvis phantom, including bone-mimicking material, was fabricated based on the computed tomography (CT) images of a prostate cancer patient. To compare the extent to which patient and phantom body and bones overlapped, the similarity Dice coefficient was calculated. Modular cylindrical inserts were created to encapsulate radiochromic films and ionization chamber for absolute dosimetry measurements at the location of prostate and at the boost region. Gamma analysis evaluation with 2%/2mm criteria was performed to compare treatment planning system calculations and measured dose when delivering a 10 flattening filter free (FFF) SABR plan and a 10FFF boost SABR plan.ResultsDice coefficients of 0.98 and 0.91 were measured for body and bones, respectively, demonstrating agreement between patient and phantom outlines. For the boost plans the gamma analysis yielded 97.0% of pixels passing 2%/2mm criteria and these results were supported by the chamber average dose difference of 0.47 ± 0.03%. These results were further improved when overriding the bone relative electron density: 97.3% for the 2%/2mm gamma analysis, and 0.05 ± 0.03% for the ionization chamber average dose difference.ConclusionsThe modular patient-specific 3D-printed pelvis phantom has proven to be a highly attractive and versatile tool to validate prostate SABR boost plans using multiple detectors.     

Measurement of the influence of titanium hip prosthesis on therapeutic electron beam dose distributions in a novel pelvic phantom

PurposeTo investigate the degree of 18 and 22 MeV electron beam dose perturbations caused by unilateral hip titanium (Ti) prosthesis.MethodsMeasurements were acquired using Gafchromic EBT2 film in a novel pelvic phantom made out of Nylon-12 slices in which a Ti-prosthesis is embedded. Dose perturbations were measured and compared using depth doses for 8 × 8, 10 × 10 and 11 × 11 cm² applicator-defined field sizes at 95 cm source-surface-distance (SSD). Comparisons were also made between film data at 100 cm SSD for a 10 × 10 cm² field and dose calculations made on CMS XiO treatment planning system utilizing the pencil beam algorithm. The extent of dose deviations caused by the Ti prosthesis based on film data was quantified through the dose enhancement factor (DEF), defined as the ratio of the dose influenced by the prosthesis and the unchanged beam.ResultsAt the interface between Nylon-12 and the Ti implant on the prosthesis entrance side, the dose increased to values of 21 ± 1% and 23 ± 1% for 18 and 22 MeV electron beams, respectively. DEFs increased with increasing electron energy and field size, and were found to fall off quickly with distance from the nylon-prosthesis interface. A comparison of film and XiO depth dose data for 18 and 22 MeV gave relative errors of 20% and 25%, respectively.ConclusionThis study outlines the lack of accuracy of the XiO TPS for electron planning in highly heterogeneous media. So a dosimetric error of 20–25% could influence clinical outcome.     

Image quality characteristics for virtual monoenergetic images using dual-layer spectral detector CT: Comparison with conventional tube-voltage images

PurposeTo investigate the image quality characteristics for virtual monoenergetic images compared with conventional tube-voltage image with dual-layer spectral CT (DLCT).MethodsHelical scans were performed using a first-generation DLCT scanner, two different sizes of acrylic cylindrical phantoms, and a Catphan phantom. Three different iodine concentrations were inserted into the phantom center. The single-tube voltage for obtaining virtual monoenergetic images was set to 120 or 140 kVp. Conventional 120- and 140-kVp images and virtual monoenergetic images (40–200-keV images) were reconstructed from slice thicknesses of 1.0 mm. The CT number and image noise were measured for each iodine concentration and water on the 120-kVp images and virtual monoenergetic images. The noise power spectrum (NPS) was also calculated.ResultsThe iodine CT numbers for the iodinated enhancing materials were similar regardless of phantom size and acquisition method. Compared with the iodine CT numbers of the conventional 120-kVp images, those for the monoenergetic 40-, 50-, and 60-keV images increased by approximately 3.0-, 1.9-, and 1.3-fold, respectively. The image noise values for each virtual monoenergetic image were similar (for example, 24.6 HU at 40 keV and 23.3 HU at 200 keV obtained at 120 kVp and 30-cm phantom size). The NPS curves of the 70-keV and 120-kVp images for a 1.0-mm slice thickness over the entire frequency range were similar.ConclusionVirtual monoenergetic images represent stable image noise over the entire energy spectrum and improved the contrast-to-noise ratio than conventional tube voltage using the dual-layer spectral detector CT.     

Comparison and Combination of Dual-Energy- and Iterative-Based Metal Artefact Reduction on Hip Prosthesis and Dental Implants

Purpose

To compare and combine dual-energy based and iterative metal artefact reduction on hip prosthesis and dental implants in CT.

Material and Methods

A total of 46 patients (women:50%,mean age:63±15years) with dental implants or hip prostheses (n = 30/20) were included and examined with a second-generation Dual Source Scanner. 120kV equivalent mixed-images were derived from reconstructions of the 100/Sn140kV source images using no metal artefact reduction (NOMAR) and iterative metal artefact reduction (IMAR). We then generated monoenergetic extrapolations at 130keV from source images without IMAR (DEMAR) or from source images with IMAR, (IMAR+DEMAR). The degree of metal artefact was quantified for NOMAR, IMAR, DEMAR and IMAR+DEMAR using a Fourier-based method and subjectively rated on a five point Likert scale by two independent readers.

Results

In subjects with hip prosthesis, DEMAR and IMAR resulted in significantly reduced artefacts compared to standard reconstructions (33% vs. 56%; for DEMAR and IMAR; respectively, p<0.005), but the degree of artefact reduction was significantly higher for IMAR (all p<0.005). In contrast, in subjects with dental implants only IMAR showed a significant reduction of artefacts whereas DEMAR did not (71%, vs. 8% p<0.01 and p = 0.1; respectively). Furthermore, the combination of IMAR with DEMAR resulted in additionally reduced artefacts (Hip prosthesis: 47%, dental implants 18%; both p<0.0001).

Conclusion

IMAR allows for significantly higher reduction of metal artefacts caused by hip prostheses and dental implants, compared to a dual energy based method. The combination of DE-source images with IMAR and subsequent monoenergetic extrapolation provides an incremental benefit compared to both single methods.     

Automatic exposure control at single- and dual-heartbeat CTCA on a 320-MDCT volume scanner: Effect of heart rate,exposure phase window setting,and reconstruction algorithm

PurposeTo investigate whether electrocardiogram (ECG)-gated single- and dual-heartbeat computed tomography coronary angiography (CTCA) with automatic exposure control (AEC) yields images with uniform image noise at reduced radiation doses.Materials and methodsUsing an anthropomorphic chest CT phantom we performed prospectively ECG-gated single- and dual-heartbeat CTCA on a second-generation 320-multidetector CT volume scanner. The exposure phase window was set at 75%, 70–80%, 40–80%, and 0–100% and the heart rate at 60 or 80 or corr80 bpm; images were reconstructed with filtered back projection (FBP) or iterative reconstruction (IR, adaptive iterative dose reduction 3D). We applied AEC and set the image noise level to 20 or 25 HU. For each technique we determined the image noise and the radiation dose to the phantom center.ResultsWith half-scan reconstruction at 60 bpm, a 70–80% phase window- and a 20-HU standard deviation (SD) setting, the imagenoise level and -variation along the z axis manifested similar curves with FBP and IR. With half-scan reconstruction, the radiation dose to the phantom center with 70–80% phase window was 18.89 and 12.34 mGy for FBP and 4.61 and 3.10 mGy for IR at an SD setting SD of 20 and 25 HU, respectively. At 80 bpm with two-segment reconstruction the dose was approximately twice that of 60 bpm at both SD settings. However, increasing radiation dose at corr80 bpm was suppressed to 1.39 times compared to 60 bpm.ConclusionAEC at ECG-gated single- and dual-heartbeat CTCA controls the image noise at different radiation dose.     

Dosimetric effect of tissue heterogeneity for 125I prostate implants

AimTo use Monte Carlo (MC) together with voxel phantoms to analyze the tissue heterogeneity effect in the dose distributions and equivalent uniform dose (EUD) for ¹²⁵I prostate implants.BackgroundDose distribution calculations in low dose-rate brachytherapy are based on the dose deposition around a single source in a water phantom. This formalism does not take into account tissue heterogeneities, interseed attenuation, or finite patient dimensions effects. Tissue composition is especially important due to the photoelectric effect.Materials and methodsThe computed tomographies (CT) of two patients with prostate cancer were used to create voxel phantoms for the MC simulations. An elemental composition and density were assigned to each structure. Densities of the prostate, vesicles, rectum and bladder were determined through the CT electronic densities of 100 patients. The same simulations were performed considering the same phantom as pure water. Results were compared via dose–volume histograms and EUD for the prostate and rectum.ResultsThe mean absorbed doses presented deviations of 3.3–4.0% for the prostate and of 2.3–4.9% for the rectum, when comparing calculations in water with calculations in the heterogeneous phantom. In the calculations in water, the prostate D₉₀ was overestimated by 2.8–3.9% and the rectum D_0.1cc resulted in dose differences of 6–8%. The EUD resulted in an overestimation of 3.5–3.7% for the prostate and of 7.7–8.3% for the rectum.ConclusionsThe deposited dose was consistently overestimated for the simulation in water. In order to increase the accuracy in the determination of dose distributions, especially around the rectum, the introduction of the model-based algorithms is recommended.     

Accuracy of dose calculation on iterative CBCT for head and neck radiotherapy

PurposeTo evaluate the feasibility of the use of iterative cone-beam computed tomography (CBCT) for dose calculation in the head and neck region.MethodsThis study includes phantom and clinical studies. All acquired CBCT images were reconstructed with Feldkamp–Davis–Kress algorithm-based CBCT (FDK-CBCT) and iterative CBCT (iCBCT) algorithm. The Hounsfield unit (HU) consistency between the head and body phantoms was determined in both reconstruction techniques. Volumetric modulated arc therapy (VMAT) plans were generated for 16 head and neck patients on a planning CT scan, and the doses were recalculated on FDK-CBCT and iCBCT with Anisotropic Analytical Algorithm (AAA) and Acuros XB (AXB). As a comparison of the accuracy of dose calculations, the absolute dosimetric difference and 1%/1 mm gamma passing rate analysis were analyzed.ResultsThe difference in the mean HU values between the head and body phantoms was larger for FDK-CBCT (max value: 449.1 HU) than iCBCT (260.0 HU). The median dosimetric difference from the planning CT were <1.0% for both FDK-CBCT and iCBCT but smaller differences were found with iCBCT (planning target volume D_50%: 0.38% (0.15–0.59%) for FDK-CBCT, 0.28% (0.13–0.49%) for iCBCT, AAA; 0.14% (0.04–0.19%) for FDK-CBCT, 0.07% (0.02–0.20%) for iCBCT). The mean gamma passing rate was significantly better in iCBCT than FDK-CBCT (AAA: 98.7% for FDK-CBCT, 99.4% for iCBCT; AXB: 96.8% for FDK_CBCT, 97.5% for iCBCT).ConclusionThe iCBCT-based dose calculation in VMAT for head and neck cancer was accurate compared to FDK-CBCT.     

Three-dimensional printer-generated patient-specific phantom for artificial in vivo dosimetry in radiotherapy quality assurance

Pretreatment intensity-modulated radiotherapy quality assurance is performed using simple rectangular or cylindrical phantoms; thus, the dosimetric errors caused by complex patient-specific anatomy are absent in the evaluation objects. In this study, we construct a system for generating patient-specific three-dimensional (3D)-printed phantoms for radiotherapy dosimetry. An anthropomorphic head phantom containing the bone and hollow of the paranasal sinus is scanned by computed tomography (CT). Based on surface rendering data, a patient-specific phantom is formed using a fused-deposition-modeling-based 3D printer, with a polylactic acid filament as the printing material. Radiophotoluminescence glass dosimeters can be inserted in the 3D-printed phantom. The phantom shape, CT value, and absorbed doses are compared between the actual and 3D-printed phantoms. The shape difference between the actual and printed phantoms is less than 1 mm except in the bottom surface region. The average CT value of the infill region in the 3D-printed phantom is −6 ± 18 Hounsfield units (HU) and that of the vertical shell region is 126 ± 18 HU. When the same plans were irradiated, the dose differences were generally less than 2%. These results demonstrate the feasibility of the 3D-printed phantom for artificial in vivo dosimetry in radiotherapy quality assurance.     

PETSTEP: Generation of synthetic PET lesions for fast evaluation of segmentation methods

PurposeThis work describes PETSTEP (PET Simulator of Tracers via Emission Projection): a faster and more accessible alternative to Monte Carlo (MC) simulation generating realistic PET images, for studies assessing image features and segmentation techniques.MethodsPETSTEP was implemented within Matlab as open source software. It allows generating three-dimensional PET images from PET/CT data or synthetic CT and PET maps, with user-drawn lesions and user-set acquisition and reconstruction parameters. PETSTEP was used to reproduce images of the NEMA body phantom acquired on a GE Discovery 690 PET/CT scanner, and simulated with MC for the GE Discovery LS scanner, and to generate realistic Head and Neck scans. Finally the sensitivity (S) and Positive Predictive Value (PPV) of three automatic segmentation methods were compared when applied to the scanner-acquired and PETSTEP-simulated NEMA images.ResultsPETSTEP produced 3D phantom and clinical images within 4 and 6 min respectively on a single core 2.7 GHz computer. PETSTEP images of the NEMA phantom had mean intensities within 2% of the scanner-acquired image for both background and largest insert, and 16% larger background Full Width at Half Maximum. Similar results were obtained when comparing PETSTEP images to MC simulated data. The S and PPV obtained with simulated phantom images were statistically significantly lower than for the original images, but led to the same conclusions with respect to the evaluated segmentation methods.ConclusionsPETSTEP allows fast simulation of synthetic images reproducing scanner-acquired PET data and shows great promise for the evaluation of PET segmentation methods.