A validation framework to assess performance of commercial deformable image registration in lung radiotherapy

IntroductionDeformable image registration (DIR) can play an important role in the context of adaptive radiotherapy. The AAPM Task Group 132 (TG-132) has described several quantitative measures for DIR error assessment but they can only be accurately defined when there is a ground-truth present in high-contrast regions. This work aims to set out a framework to obtain optimal results for CT-CT lung DIR in clinical setting for a commercially available system by quantifying the DIR performance in both low- and high-contrast regions.MethodsFive publicly available thorax datasets were used to assess the DIR quality. A “Ghost fiducial” method was implemented by windowing the contrast in a new feature provided by Varian Velocity v4.1. Target registration error (TRE) of the landmarks and Dice-similarity coefficient of the tumour were calculated at three different contrast settings to assess the algorithm in high- and low-contrast scenarios.ResultsFor the original unedited dataset, higher resolution DIR methods showed best performance acceptable within the recommended limit according to TG-132, when actual displacements were less than 10 mm. The relation of the actual displacement of the landmarks and TRE shows the limited capacity of the algorithm to deal with movements larger than 10 mm.ConclusionThis work found the performance of DIR methods and settings available in Varian Velocity v4.1 to be a function of contrast level as well as extent of motion. This highlights the need for multiple metrics to assess different aspects of DIR performance for various applications related to low-contrast and/or high-contrast regions.     

Evaluation of the Block Matching deformable registration algorithm in the field of head-and-neck adaptive radiotherapy

Background and purposeTo compare the accuracy of the Block Matching deformable registration (DIR) against rigid image registration (RIR) for head-and-neck multi-modal images CT to cone-beam CT (CBCT) registration.Material and methodsPlanning-CT and weekly CBCT of 10 patients were used for this study. Several volumes, including medullary canal (MC), thyroid cartilage (TC), hyoid bone (HB) and submandibular gland (SMG) were transposed from CT to CBCT images using either DIR or RIR. Transposed volumes were compared with the manual delineation of these volumes on every CBCT. The parameters of similarity used for analysis were: Dice Similarity Index (DSI), 95%-Hausdorff Distance (95%-HD) and difference of volumes (cc).ResultsWith DIR, the major mean difference of volumes was −1.4 cc for MC, revealing limited under-segmentation. DIR limited variability of DSI and 95%-HD. It significantly improved DSI for TC and HB and 95%-HD for all structures but SMG. With DIR, mean 95%-HD (mm) was 3.01 ± 0.80, 5.33 ± 2.51, 4.99 ± 1.69, 3.07 ± 1.31 for MC, TC, HB and SMG, respectively. With RIR, it was 3.92 ± 1.86, 6.94 ± 3.98, 6.44 ± 3.37 and 3.41 ± 2.25, respectively.ConclusionBlock Matching is a valid algorithm for deformable multi-modal CT to CBCT registration. Values of 95%-HD are useful for ongoing development of its application to the cumulative dose calculation.     

Commissioning and validation of commercial deformable image registration software for adaptive contouring

PurposeTo report the commissioning and validation of deformable image registration(DIR) software for adaptive contouring.MethodsDIR (SmartAdapt®v13.6) was validated using two methods namely contour propagation accuracy and landmark tracking, using physical phantoms and clinical images of various disease sites. Five in-house made phantoms with various known deformations and a set of 10 virtual phantoms were used. Displacement in lateral, anterio-posterior (AP) and superior-inferior (SI) direction were evaluated for various organs and compared with the ground truth. Four clinical sites namely, brain (n = 5), HN (n = 9), cervix (n = 18) and prostate (n = 23) were used. Organs were manually delineated by a radiation oncologist, compared with the deformable image registration (DIR) generated contours. 3D slicer v4.5.0.1 was used to analyze Dice Similarity Co-efficient (DSC), shift in centre of mass (COM) and Hausdorff distances Hf_95%/avg.ResultsMean (SD) DSC, Hf_95% (mm), Hf_avg (mm) and COM of all the phantoms 1–5 were 0.84 (0.2) mm, 5.1 (7.4) mm, 1.6 (2.2) mm, and 1.6 (0.2) mm respectively. Phantom-5 had the largest deformation as compared to phantoms 1–4, and hence had suboptimal indices. The virtual phantom resulted in consistent results for all the ROIs investigated. Contours propagated for brain patients were better with a high DSC score (0.91 (0.04)) as compared to other sites (HN: 0.84, prostate: 0.81 and cervix 0.77). A similar trend was seen in other indices too. The accuracy of propagated contours is limited for complex deformations that include large volume and shape change of bladder and rectum respectively. Visual validation of the propagated contours is recommended for clinical implementation.ConclusionThe DIR algorithm was commissioned and validated for adaptive contouring.     

Mixed-beam approach for high-risk prostate cancer: Carbon-ion boost followed by photon intensity-modulated radiotherapy. Dosimetric and geometric evaluations (AIRC IG-14300)

Background and purposeThe aim was to evaluate dosimetric uncertainties of a mixed beam approach for patients with high-risk prostate cancer (PCa). The treatment consists of a carbon ion radiotherapy (CIRT) boost followed by whole-pelvis intensity-modulated RT (IMRT).Materials and methodsPatients were treated with a CIRT boost of 16.6 Gy/4 fractions followed by whole-pelvis IMRT of 50 Gy/25 fractions, with consequent long term androgen deprivation therapy. Deformable computed tomography image registration (DIR) was performed and corresponding doses were used for plan sum. A comparative IMRT photon plan was obtained as whole-pelvis IMRT of 50 Gy/25 fractions followed by a boost of 28 Gy/14 fractions. DIR performances were evaluated through structure-related and image characteristics parameters.ResultsUntil now, five patients out of ten total enrolled ended the treatment. Dosimetric parameters were lower in CIRT + IMRT than IMRT-only plans for all organs at risk (OARs) except femoral heads.Regarding DIR evaluation, femoral heads were the less deformed OAR. Penile bulb, bladder and anal canal showed intermediate deformation. Rectum was the most deformed. DIR algorithms were patient (P)-dependent, as performances were the highest for P3 and P4, intermediate for P2 and P5, and the lowest for P1.ConclusionsCIRT allows better OARs sparing while increasing the efficacy due to the higher radio-biological effect of carbon ions. However, a mixed beam approach could introduce DIR problems in multi-centric treatments with different operative protocols. The development of this prospective trial will lead to more mature data concerning the clinical impact of implementing DIR procedures in dose accumulation applications for high-risk PCa treatments.     

A quantitative evaluation of deformable image registration based on MV cone beam CT images: Impact of deformation magnitudes and image modalities

Background and purposeTo evaluate the impact of deformation magnitude and image modality on deformable-image-registration (DIR) accuracy using Halcyon megavoltage cone beam CT images (MVCBCT).Materials and methodsPlanning CT images of an anthropomorphic Head phantom were aligned rigidly with MVCBCT and re-sampled to achieve the same resolution, denoted as pCT. MVCBCT was warped with twenty simulated pre-known virtual deformation fields (T_i, i = 1–20) with increasing deformation magnitudes, yielding warped CBCT (wCBCT). The pCT and MVCBCT were registered to wCBCT respectively (Multi-modality and Uni-modality DIR), generating deformation vector fields V_i and V_i′ (i = 1–20). V_i and V_i′ were compared with T_i respectively to assess the DIR accuracy geometrically. In addition, V_i, T_i, and V_i′ were applied to pCT, generating deformed CT (dCT_i), ground-truth CT (G_i) and deformed CT′ (dCT_i′) respectively. The Hounsfield Unit (HU) on these virtual CT images were also compared.ResultsThe mean errors of vector displacement increased with the deformation magnitude. For deformation magnitudes between 2.82 mm and 7.71 mm, the errors of uni-modality DIR were 1.16 mm ~ 1.73 mm smaller than that of multi-modality (p = 0.0001, Wilcoxon signed rank test). DIR could reduce the maximum signed and absolute HU deviations from 70.8 HU to 11.4 HU and 208 HU to 46.2 HU respectively.ConclusionsAs deformation magnitude increases, DIR accuracy continues to deteriorate and uni-modality DIR consistently outperformed multi-modality DIR. DIR-based adaptive radiotherapy utilizing the noisy MVCBCT images is only conditionally applicable with caution.     

Dose warping performance in deformable image registration in lung

PurposeIt is unclear that spatial accuracy can reflect the impact of deformed dose distribution. In this study, we used dosimetric parameters to compare an in-house deformable image registration (DIR) system using NiftyReg, with two commercially available systems, MIM Maestro (MIM) and Velocity AI (Velocity).MethodsFor 19 non-small-cell lung cancer patients, the peak inspiration (0%)-4DCT images were deformed to the peak expiration (50%)-4DCT images using each of the three DIR systems, which included computation of the deformation vector fields (DVF). The 0%-gross tumor volume (GTV) and the 0%-dose distribution were also then deformed using the DVFs. The agreement in the dose distributions for the GTVs was evaluated using generalized equivalent uniform dose (gEUD), mean dose (D_mean), and three-dimensional (3D) gamma index (criteria: 3 mm/3%). Additionally, a Dice similarity coefficient (DSC) was used to measure the similarity of the GTV volumes.ResultsD_mean and gEUD demonstrated good agreement between the original and deformed dose distributions (differences were generally less than 3%) in 17 of the patients. In two other patients, the Velocity system resulted in differences in gEUD of 50.1% and 29.7% and in D_mean of 11.8% and 4.78%. The gamma index comparison showed statistically significant differences for the in-house DIR vs. MIM, and MIM vs. Velocity.ConclusionsThe finely tuned in-house DIR system could achieve similar spatial and dose accuracy to the commercial systems. Care must be taken, as we found errors of more than 5% for D_mean and 30% for gEUD, even with a commercially available DIR tool.     

Accuracy of target delineation by positron emission tomography-based auto-segmentation methods after deformable image registration: A phantom study

PurposeDiagnostic positron emission tomography and computed tomography (PET/CT) images can be fused to the planning CT images by a deformable image registration (DIR). The aim of this study was to evaluate the standardized uptake value (SUV) and target delineation on deformed PET images.MethodsWe used a cylindrical phantom and removable inserts of four spheres (16–38 mm in diameter) and three ellipsoids with a volume equal to the 38-mm-diameter sphere (S38) in each. S38 was filled with 18F-fluorodeoxyglucose activity, and then PET/CT images were acquired. The contours of S38 were generated using original PET images by PET auto-segmentation (PET-AS) methods of (1) SUV2.5, (2) 40% of maximum SUV (SUV40%max), and (3) gradient-based (GB), and were deformed to the other inserts by DIR. We compared the volumes and the SUV_max with the generated contours using the deformed PET images.ResultsThe SUV_max was slightly decreased by DIR; the mean absolute difference was −0.10 ± 0.04. For SUV2.5 and SUV40%max, the differences in S38 volumes between the original and deformed PET images were less than 5%, regardless of deformation type. For the GB, the contoured volumes obtained from deformed PET images were larger than those of the original PET images for the deformation type of ellipsoids. When the S38 was deformed to the 16-mm-diameter sphere, the maximum volume difference was −22.8%.ConclusionsAlthough SUV fluctuations by DIR were negligible, the target delineation on deformed PET images by the GB should be carefully considered owing to the distortion of intensity profiles.     

Monte Carlo dose calculation in presence of low-density media: Application to lung SBRT treated during DIBH

PurposeCommercial algorithms used in Radiotherapy include approximations that are generally acceptable. However their limits can be seen when confronted with small fields and low-density media. These conditions exist during the treatment of lung cancers with Stereotactic Body Radiation Therapy (SBRT) achieved with the “Deep Inspiration Breath Hold” (DIBH) technique. A Monte Carlo (MC) model of a linear accelerator was used to assess the performance of two algorithms (Varian Acuros and AAA) in these conditions. This model is validated using phantoms with different densities. Lastly, results for SBRT cases are compared to both Acuros and AAA.MethodsA Varian TrueBeam linac was modeled using GATE/Geant4 and validated by comparing dose distributions for simple fields to measurements in water and in heterogeneous phantoms composed of PMMA and two types of cork (corresponding to lung densities during free-breathing and DIBH). Experimental measurements are also compared to AAA and Acuros. Finally, results of Acuros/AAA are compared to MC for a clinical case (SBRT during DIBH).ResultsBased on 1D gamma index comparisons with measurements in water, the TrueBeam model was validated (>97% of points passed this test). In heterogeneous phantoms, and in particular for small field sizes, very low density (0.12 g.cm⁻³) and at the edge of the field, MC model was still in good agreement with measurements whilst AAA and Acuros showed discrepancies. With the patient CT, similar differences between MC and AAA/Acuros were observed for static fields but disappeared using an SBRT arc field.ConclusionsOur MC model is validated and limits of commercial algorithms are shown in very low densities.     

Development of a dynamic deformable thorax phantom for the quality management of deformable image registration

The purpose of this study was to develop a novel dynamic deformable thorax phantom for deformable image registration (DIR) quality assurance (QA) and to verify as a tool for commissioning and DIR QA.The phantom consists of a base phantom, an inner phantom, and a motor-derived piston. The base phantom is an acrylic cylinder phantom with a diameter of 180 mm. The inner phantom consists of deformable, 20 mm thick disk-shaped sponges. To evaluate the physical characteristics of the phantom, we evaluated its image quality and deformation. DIR accuracies were evaluated using the three types of commercially DIR software (MIM, RayStation, and Velocity AI) to test the feasibility of this phantom. We used different DIR parameters to test the impact of parameters on DIR accuracy in various phantom settings. To evaluate DIR accuracy, a target registration error (TRE) was calculated using the anatomical landmark points.The three locations (i.e., distal, middle, and proximal positions) had different displacement amounts. This result indicated that the inner phantom was not moved but deformed. In cases with different phantom settings and marker settings, the ranges of the average TRE were 0.63–15.60 mm (MIM). In cases with different DIR parameters settings, the ranges of the average TRE were as follows: 0.73–7.10 mm (MIM), 8.25–8.66 mm (RayStation), and 8.26–8.43 mm (Velocity). These results suggest that our phantom could evaluate the detailed DIR behaviors with TRE. Therefore, this is indicative of the potential usefulness of our phantom in DIR commissioning and QA.     

Review of skin dose calculation software in interventional cardiology

PurposeIn interventional cardiology, patients may be exposed to high doses to the skin resulting in skin burns following single or multiple procedures. Reviewing and analysing available software (online or offline) may help medical physicists assessing the maximum skin dose to the patient together with the dose distribution during (or after) these procedures.Method and resultsCapabilities and accuracy of available software were analysed through an extensive bibliography search and contacts with both vendor and authors. Their markedly differed among developers.In total, 22 software were identified and reviewed according to their algorithms and their capabilities. Special attention was dedicated to their main features and limitations of interest for the intended clinical use.While the accuracy of the 12 software products validated with measurements on phantoms was acceptable (within ± 25%), the agreement was poor for the two products validated on patients (within ± 43% and ± 76%, respectively). In addition, no software has been validated on angiographic units from all manufacturers, though several software developers claimed vendor-independent transportability. Only one software allows for multiple procedures dose calculation.ConclusionLarge differences among vendors made it clear that work remains to be done before an accurate and reliable skin dose mapping is available for all patients.     

Validation of a method for “dose of the day” calculation in head-neck tomotherapy by using planning ct-to-MVCT deformable image registration

PurposeThe aim of this study was to test the feasibility and dosimetric accuracy of a method that employs planning CT-to-MVCT deformable image registration (DIR) for calculation of the daily dose for head and neck (HN) patients treated with Helical Tomotherapy (HT).MethodsFor each patient, the planning kVCT (CTplan) was deformably registered to the MVCT acquired at the 15th therapy session (MV15) with a B-Spline Free Form algorithm using Mattes mutual information (open-source software 3D Slicer), resulting in a deformed CT (CTdef). On the same day as MVCT15, a kVCT was acquired with the patient in the same treatment position (CT15). The original HT plans were recalculated both on CTdef and CT15, and the corresponding dose distributions were compared; local dose differences <2% of the prescribed dose (DD2%) and 2D/3D gamma-index values (2%-2 mm) were assessed respectively with Mapcheck SNC Patient software (Sun Nuclear) and with 3D-Slicer.ResultsOn average, 87.9% ± 1.2% of voxels were found for DD2% (on average 27 slices available for each patient) and 94.6% ± 0.8% of points passed the 2D gamma analysis test while the 3D gamma test was satisfied in 94.8% ± 0.8% of body’s voxels.ConclusionsThis study represents the first demonstration of the dosimetric accuracy of kVCT-to-MVCT DIR for dose of the day computations. The suggested method is sufficiently fast and reliable to be used for daily delivered dose evaluations in clinical strategies for adaptive Tomotherapy of HN cancer.     

Geometric inaccuracy and co-registration errors for CT,DynaCT and MRI images used in robotic stereotactic radiosurgery treatment planning

PurposeTo measure the combined errors due to geometric inaccuracy and image co-registration on secondary images (dynamic CT angiography (dCTA), 3D DynaCT angiography (DynaCTA), and magnetic resonance images (MRI)) that are routinely used to aid in target delineation and planning for stereotactic radiosurgery (SRS).MethodsThree phantoms (one commercial and two in-house built) and two different analysis approaches (commercial and MATLAB based) were used to quantify the magnitude of geometric image distortion and co-registration errors for different imaging modalities within CyberKnife’s MultiPlan treatment planning software. For each phantom, the combined errors were reported as a mean target registration error (TRE). The mean TRE’s for different intramodality imaging parameters (e.g., mAs, kVp, and phantom set-ups) and for dCTA, DynaCTA, and MRI systems were measured.ResultsOnly X-ray based imaging can be performed with the commercial phantom, and the mean TRE ± standard deviation values were large compared to the in-house analysis using MATLAB. With the 3D printed phantom, even drastic changes in treatment planning CT imaging protocols did not greatly influence the mean TRE (<0.5 mm for a 1 mm slice thickness CT). For all imaging modalities, the largest mean TRE was found on DynaCT, followed by T2-weighted MR images (albeit all <1 mm).ConclusionsThe user may overestimate the mean TRE if the commercial phantom and MultiPlan were used solely. The 3D printed phantom design is a sensitive and suitable quality assurance tool for measuring 3D geometric inaccuracy and co-registration errors across all imaging modalities.     

Accuracy of dose calculation algorithms for static and rotational IMRT of lung cancer: A phantom study

PurposeTo evaluate the dosimetric accuracy of Pencil beam (PB), Anisotropic Analytical Algorithm (AAA) and Collapsed Cone Convolution Superposition (CCCS) in thoracic tumours for various IMRT techniques.MethodsStep-and-shoot Linac IMRT (IMRT), arc volumetric RapidArc (RA) and Helical Tomotherapy (HT) lung treatments for different clinical situations (mediastinum tumour, single metastasis and multiple metastases) were simulated and calculated with PB/AAA, AAA, CCCS, respectively. Delivery quality assurance plans were first verified in homogeneous media (Cheese phantom and ArcCHECK); then several low-density inhomogeneous phantoms were used: the Multiplug ArcCHECK, the commercial ArcCHECK slightly modified with a low density lung–shape insert and a custom-made slab heterogeneous phantom simulating the thorax region. Absolute doses and planar dose maps were checked to assess the agreement between measured and calculated dose distributions.ResultsIn total, data referred to 195 point dose measurements and 189 planar measurements were considered. Average point absolute deviations <3% were found for all the delivery techniques/dose algorithms. In small targets completely embedded in very low density media, deviations up to 7–10% and 4–5% were found for PB and AAA/CCCS respectively. Excellent results were found for planar measurements in ArcCHECK configurations, where ≥95% of points satisfy the 3%/3 mm acceptance criteria for all the algorithms.ConclusionsA satisfactory agreement (<2%) between planned and measured doses was generally found for CCCS and AAA, excepting the very critical situation of a small tumour completely embedded in air. A significant dose overestimation (from few to 5–7%) was confirmed for PB in complex inhomogeneous arrangements.     

Accuracy of 2-hydroxyglutarate quantification by short-echo proton-MRS at 3 T: A phantom study

PurposeWe set out to investigate the potential confounding effect of variable concentration of N-acetyl-l-aspartate (NAA) and Glutamate (Glu) on measurement of the brain oncometabolite 2-hydroxyglutarate (2HG) using a standard MRS protocol. This issue may arise due to spectral overlap at clinical magnetic field strengths and thus complicate the usage of 2HG as a putative biomarker of gliomas bearing mutations of the isocitrate dehydrogenase (IDH) 1 and 2 genes.MethodsSpectra from 25 phantoms (50 mL falcon test tubes) containing a range of known concentrations of 2HG, NAA and Glu were acquired using a clinical 3 T scanner with a quadrature head coil, single-voxel point-resolved spectroscopy sequence with TE = 30 ms. Metabolite concentrations were estimated by linear combination analysis and a simulated basis set.ResultsNAA and Glu concentrations can have a significant confounding effect on 2HG measurements, whereby the negative changes in concentration of these metabolites typically observed in (peri)lesional areas can lead to under-estimation of 2HG concentration with respect to spectra acquired in presence of physiological levels of NAA and Glu.ConclusionThe confounding effect of NAA and Glu concentration changes needs to be considered: in patients, it may mask the presence of 2HG at low concentrations, however it is not expected to lead to false positives. 2HG data acquired using standard short echo-time MRS protocols should be considered with caution.     

A novel framework for spatial normalization of dose distributions in voxel-based analyses of brain irradiation outcomes

PurposeTo devise a novel Spatial Normalization framework for Voxel-based analysis (VBA) in brain radiotherapy. VBAs rely on accurate spatial normalization of different patients’ planning CTs on a common coordinate system (CCS). The cerebral anatomy, well characterized by MRI, shows instead poor contrast in CT, resulting in potential inaccuracies in VBAs based on CT alone.MethodsWe analyzed 50 meningioma patients treated with proton-therapy, undergoing planning CT and T1-weighted (T1w) MRI. The spatial normalization pipeline based on MR and CT images consisted in: intra-patient registration of CT to T1w, inter-patient registration of T1w to MNI space chosen as CCS, doses propagation to MNI.The registration quality was compared with that obtained by Statistical Parametric Mapping software (SPM), used as benchmark. To evaluate the accuracy of dose normalization, the dose organ overlap (DOO) score was computed on gray matter, white matter and cerebrospinal fluid before and after normalization. In addition, the trends in the DOOs distribution were investigated by means of cluster analysis.ResultsThe registration quality was higher for the proposed method compared to SPM (p < 0.001). The DOO scores showed a significant improvement after normalization (p < 0.001). The cluster analysis highlighted 2 clusters, with one of them including the majority of data and exhibiting acceptable DOOs.ConclusionsOur study presents a robust tool for spatial normalization, specifically tailored for brain dose VBAs. Furthermore, the cluster analysis provides a formal criterion for patient exclusion in case of non-acceptable normalization results. The implemented framework lays the groundwork for future reliable VBAs in brain irradiation studies.     

Rapid pedobarographic image registration based on contour curvature and optimization

Image registration, the process of optimally aligning homologous structures in multiple images, has recently been demonstrated to support automated pixel-level analysis of pedobarographic images and, subsequently, to extract unique and biomechanically relevant information from plantar pressure data. Recent registration methods have focused on robustness, with slow but globally powerful algorithms. In this paper, we present an alternative registration approach that affords both speed and accuracy, with the goal of making pedobarographic image registration more practical for near-real-time laboratory and clinical applications. The current algorithm first extracts centroid-based curvature trajectories from pressure image contours, and then optimally matches these curvature profiles using optimization based on dynamic programming. Special cases of disconnected images (that occur in high-arched subjects, for example) are dealt with by introducing an artificial spatially linear bridge between adjacent image clusters. Two registration algorithms were developed: a ‘geometric’ algorithm, which exclusively matched geometry, and a ‘hybrid’ algorithm, which performed subsequent pseudo-optimization. After testing the two algorithms on 30 control image pairs considered in a previous study, we found that, when compared with previously published results, the hybrid algorithm improved overlap ratio (p=0.010), but both current algorithms had slightly higher mean-squared error, assumedly because they did not consider pixel intensity. Nonetheless, both algorithms greatly improved the computational efficiency (25±8 and 53±9 ms per image pair for geometric and hybrid registrations, respectively). These results imply that registration-based pixel-level pressure image analyses can, eventually, be implemented for practical clinical purposes.     

Accurate dosimetric measurement of large extended SSD fields for comparison to TPS models

PurposeThere is little evidence in the literature which quantifies the accuracy of Treatment Planning Systems (TPSs) using large fields at extended SSD (eSSD). This paper introduces the approach taken at Christchurch Hospital, New Zealand to validate the use of the Monaco TPS for Total Body Irradiation (TBI) treatments.MethodsA purpose-built device for allowing precise movements of block-like phantoms called a Phantom Mobility Device (PMD) was used for collecting measurements at eSSD. These measurements were used for determining the ability of the Monaco TPS (originally validated for SSDs between 80 and 110 cm) to accurately model dose distributions for TBI treatments at Christchurch Hospital on either treatment machine one (T1) or two (T2) with SSD values of 341 and 432.6 and clinically useful field sizes of 120 and 170 cm, respectively.ResultsWe found that within the limits of measurement uncertainty the PMD contributed no determinable scatter to the measurements and proved a reliable approach for eSSD dose measurements. Additionally, by applying depth and off-axis distance constraints of use for TPS information it is possible to use the existing Monaco CCC model at eSSD for block phantom geometries. Dose Difference (DD) analysis showed a clinically acceptable agreement between the CCC model and measured data over a range of depths and off-axis distances.ConclusionsThe PMD was determined to be a useful tool for accurate measurement of extended SSD treatment fields. Monaco TPS CCC model agreed well for block phantoms so future comparisons to anthropomorphic phantoms or patient data are feasible.     

Calculation of tumour and normal tissue biological effective dose in 90Y liver radioembolization with different dosimetric methods

PurposeRadioembolization with ⁹⁰Y microspheres is an effective treatment for unresectable liver tumours. Two types of microspheres are available: resin (SIR-Spheres®) and glass (Theraspheres®). The aim of this study is to compare biological effective dose (BED) values obtained with three different dosimetric methods.Methods29 HCC patients were included in this study: 15 were treated with resin(mean injected activity 1.5 GBq, range 0.8–2.7 GBq) and 14 with glass microspheres (2.6 GBq, range 1.3–4.1 GBq). Average doses to tumours and normal liver tissues were calculated with AAPM, multi-compartmental MIRD and Voxel-based methods and consequently the BED values were obtained. Planar images were used for the AAPM method: ^99mTc-MAA SPECT-CT attenuation and scatter corrected images (resin) and ⁹⁹m Tc-MAA SPECT attenuation corrected (glass) were employed for the other two methods.ResultsRegardless of type of microspheres, both for tumours and normal liver tissues, no significant statistical differences were found between MIRD and Voxel for both doses and BED values. Conversely AAPM gave discordant results with respect to the other two methods (Mann-Whitney p-values ⩽ 0.01). For resin spheres the calculated tumour-to-normal tissue ratios on planar images were on average 14 times greater than those obtained on SPECT-CT images, while they were 4 times greater on glass. A linear correlation was observed between MIRD and Voxel BEDs.ConclusionsThe AAPM method appears to be less precise for absorbed dose and BED estimation, while MIRD and voxel based dosimetry are more confident each other.     

Polyvinyl alcohol cryogel based vessel mimicking material for modelling the progression of atherosclerosis

PurposeThe stiffness of Polyvinyl-alcohol cryogel can be adjusted through application of consecutive freeze-thaw cycles. This material has potential applications in the production of tissue mimicking phantoms in diagnostic ultrasound. The aim of this study was to use PVA-c to produce a range of geometrically and acoustically identical vessel phantoms modelling stages of atherosclerosis which could be verified through mechanical testing, thus allowing for more precision in quantitative in-vitro flow analysis of atherosclerosis.MethodsA series of anatomically realistic walled renal artery flow phantoms were constructed using PVA-c. In order to ensure precise modelling of atherosclerosis, a modified procedure of ISO27:2017 was used to compare the mechanical properties of PVA-c. Results were compared for the standard “dumbbell” test object and a modified vessel test object. The geometric accuracy and reproducibility of the vessel models were tested before and after implantation in flow phantoms.ResultsNo significant difference was found between the mechanical properties of the dumbbell test samples and the vessels for any number of freeze thaw cycles, with a correlation coefficient of R² = 0.9767 across the dataset, indicating that a direct comparison between the mechanical properties of the dumbbell test samples and the phantom vessels was established. The geometric reproducibility showed that before and after implantation there was no significant difference between individual vessel geometries (p = 0.337 & p = 0.176 respectively).ConclusionsPolyvinyl-alcohol cryogel is a useful material for the production of arterial flow phantoms. Care should be taken when using this material to ensure its mechanical properties have been correctly characterised. The guidelines of ISO37:2017 potentially provide the best procedure to ensure this.