Browse
Recent Submissions
Item Open Access Assessment of ventricular function using gated blood pool planar and - SPECT imaging: a phantom study(University of the Free State, 2023) Pieters, Hané; van Staden, J. A.In the field of Nuclear Medicine, Gated Blood Pool (GBP) investigations are essential in offering vital insights into cardiac function, specifically the Left Ventricular Ejection Fraction (LVEF). The evaluation of ventricle volume changes during the end-diastolic (ED) and end-systolic (ES) phases plays a critical role in detecting, diagnosing, and managing various cardiac diseases. While the longstanding preference for Gated Blood Pool Planar (GBP-P) methods lies in their validation, non-invasiveness, and straightforward application, challenges such as anatomical overlap and the need for background correction persist. The theoretical superiority of three-dimensional (3D) analogues, specifically Gated Blood Pool SPECT (GBP-S) studies, promises to overcome GBP-P challenges by offering true volumetric representation without the need for background correction. However, this transition introduces complexity in algorithms and processing software needed for GBP-S studies. Accuracy and precision in determining LVEF is paramount in both GBP-P and GBP-S methods to avoid misdiagnosis, improper treatment, or negligence. Rigorous testing and comparison to known or true values are imperative for validating GBP processing software programs to meet set standards. A key advancement in validating these software programs is the use of digital hybrid phantoms, notably the advanced 4D-XCAT model. The model, blending voxelised and mathematical elements, mimics human anatomy and physiology. Paired with the Monte Carlo (MC) code, SIMIND, these 4D-XCAT models can be used to simulate clinically realistic GBP images, generating a database for testing, and validating various GBP software packages. Importantly, this approach avoids radiation exposure to patients and researchers and enhances the reliability of outcomes by providing benchmark input parameters for software evaluation. The primary aim of this investigation was to assess ventricular function using MC-simulated GBP-P and GBP-S images of digital patient phantom studies based on 4D-XCAT models with varying cardiac volumes and functions. The aim was achieved by considering three objectives. Firstly, a Monte Carlo simulated cardiac phantom for planar and SPECT studies was validated. The modelled gamma camera was verified using routine quality control procedures outlined by the National Electrical Manufacturers Association (NEMA). Furthermore, 3D cardiac phantoms were printed and imaged according to GBP-P and GBP-S imaging guidelines. Simulated images of these phantoms were generated using the MC code SIMIND and compared to the gamma camera-acquired images. The successful verification of these simulated images led to the next step, namely verifying the use of the 4D-XCAT model in image simulations. By simulating GBP-P and GBP-S studies of a single 4D-XCAT model, the study demonstrates excellent agreement between known and calculated ventricular parameters, confirming the reliability of the 4D-XCAT model in simulating cardiac imaging. Building on the successful simulation of the 4D-XCAT model, the second objective involved the creation of a comprehensive database comprising 64 clinically realistic GBP patient models. GBP-P simulated images from these models were utilised to evaluate four commercially available GBP-P processing software programs. The study yielded a strong correlation between calculated LVEF values and known values, thereby affirming the reliability and accuracy of the GBP-P processing software. Lastly, the research was extended to include clinically realistic GBP-S studies as part of the third objective. The phase focused on assessing the commercially available GBP-S processing software, Quantitative Blood Pool SPECT (QBS), from Cedars-Sinai. GBP-S images of the 4D-XCAT GBP database were simulated according to established imaging guidelines. The study assessed the accuracy of QBS in terms of LV volumes and EF values for two different reconstruction techniques and found a strong correlation between the calculated and known LV volumes and EF values. This paved the way for future multi-centre studies to validate other GBP-S processing software packages.Item Open Access Quantification accuracy for I-123 SPECT/CT studies using LEHR and me collimators: a Monte Carlo study(University of the Free State, 2021) Richards, Anneray; Van Staden, J. A.; Du Raan, H.The accurate quantification of Nuclear Medicine single photon emission tomography (SPECT) plays an important part in radiopharmaceutical therapy. Accurately quantifying SPECT images of a diagnostic radionuclide such as I-123 is desirable, although not a straightforward process as it is hindered by the complex decay scheme. I-123 has low-energy primary emissions of 159 keV, and performing acquisitions with a low-energy resolution (LEHR) collimator result in images with high resolution. However, I-123 also has high-energy photon emissions which degrade image contrast and quantification accuracy. This degradation can be reduced by using medium-energy collimators (ME); however, at the expense of spatial resolution. Most clinical facilities have access to LEHR collimators, but not necessarily ME collimators. The aim of this study was to evaluate the quantification accuracy of I-123 LEHR and ME collimated SPECT images when an optimised OSEM reconstruction protocol is applied. To accomplish the aim three objectives were identified: 1) validation of a SIMIND modelled gamma camera fitted with LEHR and ME collimators, 2) optimisation of the iterative reconstruction algorithm in terms of equivalent iterations and SPECT corrections, and based on these results, 3) evaluation of the quantification accuracy of I-123 LEHR and ME SPECT images. The first objective of this study, to validate the SIMIND modelled gamma camera fitted with LEHR and ME collimators for I-123, involved comparing measured and simulated I-123 data. Results of measured and simulated planar performance tests (system energy resolution, system spatial resolution, and system sensitivity) were compared for both collimators. The validation included a visual comparison of reconstructed SPECT images of a quality control phantom in terms of uniformity, cold contrast, resolution, and linearity. The measured and simulated planar results for system energy resolution, system spatial resolution and system sensitivity differed by 3.4%, 6.4% and 5.3%, respectively. The visual comparison performed on the reconstructed SPECT images showed good agreement between the measured and simulated data. The second objective was to optimise the OSEM iterative reconstruction algorithm concerning the number of iterations and SPECT corrections. SPECT images of voxel-based phantoms of spherical objects and image quality phantoms were simulated and reconstructed with different numbers of effective iterations. The count density recovery, image noise, contrast and resolution were evaluated. The image quality phantom was also reconstructed with different corrections (attenuation, scatter and collimator-detector response (CDR)) and compared. The optimal number of equivalent iterations was selected as 64 and the contribution of the different corrections was appreciated. When septal penetration and scatter was compensated for as part of the CDR correction, the LEHR collimator results were comparable to that obtained with the ME collimator. This led to the aim of the final objective: to determine the quantification accuracy of I-123 SPECT studies in patient phantoms acquired with LEHR and ME collimators. Using voxel based patient phantoms, the quantification accuracy was assessed for LEHR and ME SPECT images of spherical objects. Quantification errors smaller than 3.8% were recorded for both the LEHR and ME collimators when attenuation, scatter and CDR (including septal penetration and scatter) corrections were applied. Therefore, to conclude, when appropriate SPECT corrections were applied during the reconstruction of I 123 LEHR and ME SPECT images, the image quality between the collimators were comparable and quantification accuracy of up to 3.8% was achievable.Item Open Access Computed tomography radiomic texture features dependence upon imaging parameters(University of the Free State, 2019-11) Makosa, Frank; Rae, W. I. D.; Court, L. E.; Acho, S. N. N.Introduction and Aim: Few studies have been carried out to determine the influence of Computed Tomography (CT) acquisition parameters (slice thickness, tube potential difference (kVp), and tube current time product (mAs)) on the quantitative image features in radiomics studies. There is little evidence in the published literature, of studies that use mathematics to establish radiomic texture features that are independent of the CT scan technique parameters. The stability of radiomic texture features may have a great impact on the diagnosis and treatment of cancers. Robust texture features can be used to track radiotherapy treatment response. In this study radiomic texture features were investigated to identify features that did not depend on the CT technique parameters. Methodology: The credence cartridge radiomics (CCR) phantom was imaged at four CT units at the Universitas Academic and the National District hospitals. The tube current-time product (mAs) was varied from 75 to 400 mAs in steps of 25mAs while the kilovoltage peak and slice thickness were kept set at 120kVp and 5 mm respectively. The CT tube potential was investigated at 80, 100, 120 and 135 kVp whilst mAs and slice thickness was kept set at 300 mAs and 5 mm respectively. The slice thickness was varied from 1 mm to 5 mm whilst the mAs and kVp was kept constant at 300 mAs and 120kVp respectively. The acquisition field of view (FOV) and pitch were kept constant. The images obtained were processed using PyRadiomics software platform of 3D Slicer and the Matlab 2017a package. PyRadiomics was used to segment and extract a total of 105 radiomics texture features for each region of interest (ROI) delineated on an image. The 105 radiomic features included 13 shape features, 18 first order statistics features, 23 grey-level co-occurrence matrix, 14 grey level difference matrix, 16 grey-level run length matrix, 16 grey level size zone matrix and 5 neighbourhood grey tone difference matrix features. For each 10 CCR phantom inserts, 16 ROI of 2cm diameter was segmented by aligning the centre of the ROI at the centre of the insert. The Matlab package was used to segment and extract image matrices that were used to perform hand GLCM calculations. A kV Cone Beam Computed Tomography (kV CBCT) acquired cervical cancer data-set was used to establish the robust radiomic texture features response to radiotherapy treatment. The kV CBCT images were acquired first day and weekly during the 25 treatment fractions. Results: Five first order statistic radiomic features and six grey level co-occurrence matrix features were identified in the experimental test and mathematical manual calculations tests to vary with coefficients of variance of less than or equal to 10 % when the slice thickness was varied. Most of the radiomic texture features were weak and unstable (coefficients of variance above 10%) at very small slice thickness (≥2.5 mm) and robust at medium (≥2.5 mm) to large slice thickness (3.75 mm and 5 mm) (coefficients of variance ≤ 10 %). The above was attributed to an averaging effect (image smoothening) on the images when the slice thickness of image acquisition is increased. The image noise was observed to be less in large slice thickness when compared to noise at small slice thickness. Radiomics features were independent and stable to the tube potential at greater than 100 kV. At high tube potential the radiation attenuated signal detected at the CT detector was higher cancelling the noise effects. The robustness of these radiomic features depended on the material comprising the insert analysed. The extent of mAs dependence observed for the dense cork and plaster resin materials inserts was low compared to the dependence on the solid acrylic material insert. All the other phantom inserts (rubber particles, natural cork and the 3 acrylonitrile butadiene styrene plastic) data plots showed smaller variations around the central axis (zero feature value) of the skewness, uniformity, entropy and kurtosis features graphs. Irrespective of the mAs changes, the radiomic texture feature values obtained from all of the ABS materials inserts, rubber particles and natural cork inserts were consistently smaller, closer to zero. A general decrease in image noise as the mAs of image acquisition was increased in images of uniform or relatively uniform material was also observed. The patient tumour analysis showed some radiomic texture features response to radiotherapy treatment. This was shown by the changes observed on the inverse difference, inverse difference moment, entropy and difference variance texture features. The texture features had their values decrease from start of treatment (first fraction) to the last treatment fraction. The decrease was not smooth along the treatment period, there were some anomalies on the trends. This decrease was ascribed to the change in the heterogeneity of the tissues within the treatment region of interest evaluated. Conclusion: Overall, using theoretical analysis and a practical approach, robust radiomic features that were independent of the CT scan parameters were observed. The experimental approach showed that the phantom insert materials had influence on radiomic texture feature values obtained in investigations. Radiomic texture features demonstrated that tumours had a variation of heterogeneity between them. The observation agrees with other clinical studies that showed that tumours exhibit some extensive genetic and phenotypic variations. Radiomic texture features can be utilised to depict tumour texture changes along the treatment timeline as shown in this study. A great challenge would be to associate the radiomic texture feature changes to the clinical biological changes. For future robust radiomic feature studies, the use of phantoms with tissue like materials was proposed.Item Open Access Development of a brachytherapy treatment planning module for cervix cancer utilising biological dose metrics(University of the Free State, 2020-12) Van der Walt, Hester Catharina; Shaw, WilliamThe contouring uncertainties associated with the use of computed tomography imaging for brachytherapy planning creates the need to investigate an alternative planning method for Image-Guided Adaptative Brachytherapy. This alternative method needs to be more robust against imaging and contouring uncertainties compared to the original GEC-ESTRO prescription regarding dose-volume histogram criteria. This study evaluates the utilisation of biological dose metrics (equivalent uniform dose (EUD)) during Image-Guided Adaptive Brachytherapy (IGABT) treatment planning and optimisation. A retrospective planning study was conducted. The eighteen patients that were included in the planning study received CT-based Image-Guided Brachytherapy (IGBT) in combination with external beam radiotherapy (EBRT) between 2014 and 2015. A novel biological optimisation model was developed and used to efficiently and automatically optimise brachytherapy (BT) plans by utilising either dose-volume, biological metrics/indexes or both for fast treatment plan generation. The module was refined to allow forward and inverse planning and optimisation of combined interstitial and intracavitary brachytherapy. Additionally, the utilisation of OAR total dose constraints during the planning process was applied and evaluated. The results of the study showed that the inverse optimisation tool could produce better target volume coverage compared to the manual optimisation tool. In some cases, the inverse optimisation tool led to higher OAR doses; however, the values recorded were still within set constraints. When comparing the conventional and biological planning methods, the biological planning produced superior CTV-T doses and dose distributions within the CTV-T. The inverse biological approach reported significantly higher average CTV-THR D98 % , D100 % values and CTV-TIR EUD, D90 % D98 % and D100 % values compared to inverse conventional planning approach. With this, the inverse biological approach also had the ability to record significant lower average bladder wall EUD and D0.1cm3D0.1cm3 values. Even though the inverse conventional planning approach reported significant lower average small bowel D2cm3 values compared to the IBG approach, both approaches were still well below the D2cm3 hard constraint of 75 Gy. Dose escalation was achieved in the CTV-T with a reduction in OAR dose with the combination of interstitial/intracavitary brachytherapy. It was concluded from the study that the incorporation and utilisation of biological metrics, which incorporates the entire dose distribution in the organ of interest, is the preferred approach when compared to the conventional physical dose-volume approach.Item Open Access Comparison between measured and simulated activity using Gafchromic™ film with radionuclides(University of the Free State, 2020-08) Joubert, Maria Magdalena; Du Plessis, F. C. P.; Van Staden, J. A.In this study, Gafchromic™ film XR-QA2 and RT-QA2 were used to characterise the film energy response against various radionuclides. The film response was investigated with respect to different backscatter materials. The sensitivity of the two types of films was compared, and a film stack method was tested to allow the user to obtain sequential, cumulative doses at different time points. Monte Carlo (MC) simulations were used to link optical density (OD) values from measurements to the absorbed dose in the film. This was achieved by using conversion factors obtained by BEAMDP, BEAMnrc and DOSXYZnrc simulations to get the absorbed dose in the film. A neutron depletion theoretical model was introduced that can describe film response as a function of cumulated activity and absorbed dose. Background: Gafchromic™ film has been used for quality assurance in various studies but not in nuclear medicine applications. Once the OD has been determined after film exposure to a radionuclide, it can be linked to the absorbed dose using the air kerma rate constant at distances that approximates point sources and the dose in water can be linked to the dose in film using MC simulations to get conversion factors. MC simulations are known as a gold standard to get the absorbed dose in materials. Materials and Methods: XR-QA2 and RT-QA2 Gafchromic™ film were irradiated with the following radionuclides: Am-241, Cs-137, Tc-99m and I-131. The OD was calculated, and a function describing the relationship between the OD and the time-activity was derived based on the neutron depletion model. Different backscatter materials such as Corrugated fibreboard carton (CFC) or air equivalent material, polystyrene, Polymethyl Methacrylate (PMMA or perspex) and lead were used to investigate the effect it has on film response. The sensitivity of each film was investigated and compared. BEAMDP, BEAMnrc and DOSXYZnrc simulations were used to link the film response, OD, to the absorbed dose. The MC simulations were done replicating the exact geometry as with the physical measurements to get the absorbed dose in the film. Results: The new neutron depletion model fitted the OD vs cumulative activity accurately as well as the OD vs absorbed dose. The XR-QA2 Gafchromic™ film has shown to be the most sensitive film when using air equivalent material with radionuclides, especially with low energy radionuclides such as Am-241. When using more than one layer, the OD sensitivity of the film can be increased as well. The film stack method investigated also showed to be less time consuming when relating stacked film data to single film data. The fluence obtained from BEAMDP confirmed that the radionuclide containers have an effect on the radionuclide spectra’s. Lead was also the backscatter material which showed higher OD change but lower absorbed dose values. Conclusions: The neutron depletion theoretical model is more accurate than higher-order polynomial fits because it contains less free parameters. The XR-QA2 Gafchromic™ is better to use in nuclear medicine because of its sensitivity at low energies and because the sensitivity can be increased by using multiple layers of film. Film stack methods can be used to decrease experiment times. BEAMnrc can be used to accurately model radionuclides within their containers to evaluate the container effects. Lead showed a higher induced OD with lower absorbed dose, and the air equivalent material showed the lower OD change but higher absorbed dose.Item Open Access Estimation of the eye lens doses in a catheterization laboratory from available image parameters(University of the Free State, 2020-07) Phutheho, Mokete; Acho, Sussan; Rae, William; Rose, AndreBackground and objective: New data on eye lens dosimetry supports the theory that the threshold of radiation-induced cataracts could be substantially lower than previously believed with some investigators arguing that cataracts could be classified as a stochastic rather than a deterministic effect. Based on these new data, the International Commission on Radiological Protection (ICRP) has reduced the occupational eye lens dose limit from 150 mSv to 20 mSv averaged over a defined period of 5 years, with no single year exceeding 50 mSv. The new reduction in the annual dose limit will have considerable implications particularly in high exposure environments such as interventional cardiology and radiology. It is therefore imperative that strategies for effective dose reduction, radiation protection, eye dose monitoring, and dosimetry be implemented in countries that have already adopted the new eye dose limit. The main aim of this study was to develop methods that can be applied to estimate eye dose equivalent from the available imaging parameters and whole-body equivalent measured over the lead apron at the chest level. The study also aimed to establish a method to estimate eye lens dose based on the workload of interventionalists. Material and methods: The study included four interventional cardiologists. A total of 127 procedures were performed in a period of three months. The procedures were categorised into diagnostic (CA) and therapeutic (CA+PCI) procedures. During these procedures, two different active dosimeters were used to measure scatter dose (one attached on the canthus of the protective eyewear to measure eye lens dose (ELD) and the other at the chest level to measure whole-body dose) to the cardiologists. The dose area product (DAP), air kerma (Ka,r), fluoroscopic time, total cine images were recorded after every procedure. The efficacy of the protective eyewear used at Universitas Hospital was evaluated in a separate study. Results: Average eye dose per CA and CA+PCI procedures were 195.1±112 and 391.8±202.9 µSv, respectively. The average dose per procedure obtained by combining all the monitored procedures was 250.9±168.3 µSv. The minimum workload necessary to exceed the annual eye lens dose limit calculated using an equation established in this study was 80 procedures. The dose reduction factor of the protective eyewear was ~2. Applying this factor increased the minimum procedures necessary for a doctor to exceed the limit to 160 procedures per year. Excellent correlation was found between ELD and DAP (R2 = 0.78). Excellent correlation was also found between ELD and Ka,r (R2 = 0.72). A poor but significant correlation was found between ELD and chest dose (R2 = 0.45). Three methods based on the ratios of ELD to DAP, Ka,r and chest dose were established. The calculation error using the methods based on DAP and Ka,r was ±20%. The respective calculation error was ±37% using the method based on chest dose. Conclusion: The accumulated eye dose of interventional cardiologists working at the Universitas Hospital can easily surpass the newly set annual eye lens dose limit after performing relatively low numbers of interventional procedures. The high average dose per procedure reported in this study highlights immediate need for implementation of radiation optimization strategies to mitigate the risk of radiation-induced cataracts. This is the first study in South Africa to establish methods that can be used to estimate eye lens doses at any time. More research is needed in the South African context to further investigate eye lens dose in interventional suits. This will allow for comparison of results obtained at different institutions and improvement in accuracy of estimation methods.Item Open Access Characterization of small megavoltage photon beams for radiography(University of the Free State, 2017-07) Setilo, Itumeleng; Du Plessis, F. C. P.English: Introduction The landscape of radiation treatment techniques is ever evolving in pursuit of improved target coverage. The latest techniques such as IMRT, SBRT, SRS and VMAT, provide improved target coverage by controlling the intensity of the given dose through the use of multiple small fields in contrast to large fields in conventional treatments. The advantage of using these large fields is that, their characteristics are fully understood. The introduction of small fields leads to improved coverage, but the physics of these fields are not fully understood. So, when used in patient treatment, it resulted in unaccounted radiation exposure due to inaccurate commissioning and inaccurate absolute dose calibration at these field sizes. The errors were due to incorrect detectors used for data collection, and incorrect application of factors when performing absolute dose calibration. This report investigated the characteristics of these small fields using different detectors whilst varying the SSD and the incident photon beam energy. The measurements included beam profiles, percentage depth dose (PDD) curves as well as the relative output factors (ROF). Materials and Methods The photon energies, 6 MV, 10 MV and 15 MV were delivered using the Synergy LINAC, which is equipped with Agility multileaf collimators (MLCs). The detectors that were investigated were the CC01 ion chamber, EFD-3G diode, PTW60019 microdiamond, EBT2 radiochromic film and the EDR2 radiographic film. Measurements were carried out using water as a medium for the CC01 ion chamber, EFD-3G diode and the PTW60019. Films were placed in between water equivalent RW3 phantom slabs. These measurements were carried out at 90 cm, 95 cm, 100 cm and 110 cm source to surface distances (SSD). The field sizes that were investigated were 1×1 cm², 2×2 cm², 3×3 cm², 4×4 cm², 5×5 cm² and 10×10 cm², these fields sizes were set using Jaws and MLCs. The 10×10 cm² field size was included as a reference field. Results and Discussion The results showed that the beam profiles were insignificantly different at the various SSDs for the detectors. The EBT2 film showed the sharpest penumbra, with the EDR2 and the CC01 showing broad penumbrae, but the difference was negligible. The PDD measurements showed that the difference between the detectors after Depth of maximum dose (Dmax) were insignificant. The films differed significantly at shallower depths, and this can be attributed to setup, as well as the artefacts that showed up when the films were being analyzed. The PDD measurements indicated that the setup used for the films was not adequate for measuring the 1 cm square field sizes and below. Dmax was used to compare the detectors, though it did not vary greatly for the detectors, it was shown that there is a change in the manner in which this factor changes with field size. Below a certain field size, 2 cm for the 6 MV and 10 MV and 3 cm for the 15 MV, the Dmax would start shifting back to the surface instead of moving deeper as expected. The relative output factor (ROF) increased with energy, and this is true for all the fields which had lateral electronic equilibrium (LEE). This relation broke down as the field sizes decreased due to the onset of lateral electronic disequilibrium (LED). The high-density detector, PTW60019 gave the highest ROF for the different energies, with the less dense CC01 giving the lowest ROFs. This showed that the density of the detector had an effect on the output factor measured. Conclusion The fields were characterized with the different detectors, barring the artefacts experienced with film measurements in some instances, these detectors can be used safely for the small fields. The ROFs can be measured at longer SSDs as they showed little variation due to increased SSDs.Item Open Access Commissioning and optimization of a total skin electron therapy technique using a high rate electron facility(University of the Free State, 2007-06) Yousif, Yousif Abd Alla Mohammed; Willemse, C. A.English: Total skin electron therapy (TSET) is the treatment of choice for several malignant diseases of the skin (Kaposi sarcoma, mycosis fungoides). Several different techniques have been developed in various centers, in order to achieve homogeneous dose distribution over a large irradiation field (200 x 80 ern"). However, to implement a TSET technique one has to account for a variety of parameters, from geometric (room design, space constrains) to physical (number, angle and energy of the beams). To obtain the most acceptable dose distribution an extensive set of measurements and a large number of calculations have to be performed. Therefore Monte Carlo simulation of TSET can facilitate optimization of this technique. In this study we implemented and optimized a TSET technique using 4 and 6 MeVelectron beams. The dosimetrie procedure intended to obtain adequate dose uniformity over the entire surface of the patient, and to reduce the patient treatment time using a high dose rate facility on the Elekta Precise accelerator. The EGS4/BEAM code package running on a Windows based platform was used for the MC simulation. Percentage depth-dose curves and beam profiles were calculated and measured experimentally for the 40x40 cm2 nominal field at both 100 cm SSD and at the patient surface at the treatment plane (SSD 350 cm) for a single beam. The accuracy of the simulated beam was validated by the good correspondence (within less than 2%) between measured beam characteristic parameters (Rso, dmax, Rp) and Monte Carlo calculated results. To obtain a uniform profile vertically, two vertical angles of incidence were used. The angle between the two beams that gave best uniformity was considered the optimum angle. The patient is to be placed on a rotating platform perpendicular to the beam and rotated through 60 degree increments to obtain six horizontal directions of beam incidence. The doses expected in the patient were measured with Kodak EDR2 films positioned at different levels between slices of a Rando phantom. TLDs were placed on the surface to relate the film measurements to dose. The delivered doses in the treatment plane were compared to simulated data that was obtained from the MC simulation. The penetration depth of the dose distribution varied over various scanning directions between 2-3 mm and 3- 4 mm for 4 and 6 MeV respectively. This information is useful when treatment of lesions of different thickness are being considered. The composite percentage depth dose of all six dual fields for both 4 and 6 MeV yielded an 80 % dose at - 7 mm and - 9 mm depth, respectively. Good dose uniformity was achieved for both energies and it was about ± 5% for 4 MeVand about ± 3% for 6 MeVover a range of - 100 to +100 cm. The bremsstrahlung contamination was 0.9 and 1.3 %. Generally there was good agreement between the dose distribution calculated with MC and measured with films, thus validating our MC calculations. The dose distributions in phantom were found to comply with the guidelines described in the AAPM TG-23 protocol, showing the suitability of this technique for treatments of the skin diseases. The HDRE is a useful operational mode providing reasonable output, field size, and Xray contamination. Use of a dual field technique produces reasonable beam uniformity over an area large enough to allow total skin electron therapy in a conventional treatment room. Monte Carlo techniques provided a guiding principle to assist the verification of the beam characterization of a TSET technique. The absolute calibration of dose to the patient required the measurement of the ratio "skin dose to calibration point dose"; thiswas achieved by measurements with a parallel plate ionization chamber and TLDs.Item Open Access Kwantifisering van radionukliedverspreiding deur planare beelding met 'n sintillasiekamera(University of the Free State, 2000-06) Marais, Johan; Van Aswegen, A.English: In this study attenuation and scatter correction techniques in planar imaging were investigated as well as the practical implementation of these techniques in clinical studies. The attenuation correction technique that was investigated was the geometrical mean method where an iterative method was used to determine the attenuation coefficient. This adaptation made the attenuation correction method independent of depth. In this study the depth dependence of the geometrical mean method was solved which is an important contribution. The three-energy window scatter correction technique was investigated with Monte Carlo methods for planar imaging. These simulations showed that the scatter correction method under corrects by approximately 10%. The scatter correction technique cannot be performed on it's own but must be used in conjunction with the attenuation correction method. The overall effect of the quantitation (attenuation and scatter correction) was determined in two ways. First a known amount of activity was imaged. Then the correction was applied to the images. Subsequently the corrected activity obtained from the latter images was compared to the known value. The calculated activity was 98,6 % of the known activity value. Secondly clinical studies were performed where the bladder activity was quantitated in vivo and correlated to the in vitro measured value. In this case the correlation coefficient was 0,996. The quantitation technique was applied to the clinical determination of the glomerular filtration rate. The quantitation technique was compared to the routine method that requires handling of urine. The correlation coefficient was 0.96 and the linear regression line had a slope ofO.91 and a y-axis intercept of4.8l. The quantitation of the bladder activity leads to a simplification of the clinical study procedure since the physical handling of urine can be eliminated. The accuracy of this method depends on the conversion factor that compensates for the difference between the sensitivity of the camera and the well counter. This factor must be checked on a regular basis to ensure the accuracy of this method. The bladder as source organ was used in this study. The bladder is a relative isolated organ and therefor the influence of radioactive uptake in other organs is relatively small. If these quantitations are applied to organs that are not isolated, i.e. the heart where the liver can influence it, this problem has to be address. In this project a quantitative planar imaging technique was investigated that corrects for attenuation and scatter and thereby facilitates depth independent corrections.Item Open Access An evaluation of the effect of scatter and attenuation correction of gamma photons on the reconstructed radionuclide distribution in the myocardial wall during spect imaging(University of the Free State, 2000-11) Mdletshe, Nhlakanipho; Van Aswegen, A.; Du Raan, H.The purpose of this study was firstly to evaluate the selection of reconstruction parameters (i.e. the number of subsets and the number of iterations) based on phantom studies. The second aim was to evaluate the effect of the non-uniform attenuation and scatter correction on myocardial perfusion studies performed on healthy volunteers as well as patients with proven inferior wall perfusion defects. The quality of the images from the phantom studies showed that 16 subsets with 2 iterations gave the best results if considering image noise and image resolution. These number of subsets and iterations were therefore used as reconstruction parameters in the patient studies The application of an attenuation correction to the emission data required that attenuation coefficient maps of the subjects were obtained from transmission images. I39Ce was chosen as the transmission source and used in conjunction with 99mTcas the emission source. The emission data were corrected for scatter according to the triple energy window method. In the healthy male and female volunteers, the attenuation and scatter corrected myocardial SPECT images showed an improvement in the homogeneity of the counts distribution compared to the uncorrected images. The counts distribution in the inferior region improved after the attenuation correction was applied, however it exceeded the counts in the anterior region. After applying a scatter as well as an attenuation correction to the emission data, the counts in the inferior region of the myocardium were slightly reduced. This was a result of the scatter correction eliminating scattered counts in the inferior region originating mainly from the liver. The apparent lower counts in the anterior region could be a result of too little compensation for scatter in the inferior wall, and needs to be investigated further. The defects in the three unhealthy patients, were not obscured after applying the scatter and attenuation correction to the emission data. The correction technique did not introduce false negative results in these patients. The application of scatter and attenuation correction techniques shows promising results for the interpretation of myocardial perfusion studies. These correction algorithms however need to be investigated thoroughly before being used in the routine clinical practice to avoid the introduction of artefacts.Item Open Access Using Monte Carlo techniques to evaluate the dose distributions from a radiotherapy treatment planning system(University of the Free State, 2000-06) Awusi, Kavuma; Lotter, M. G.; Willemse, C. A.English: In this study we used Monte Carlo techniques to simulate the SL25 linear accelerator treatment head using the BEAM Code. The main purpose of study was to evaluate the dose distributions obtained by the CADPLAN treatment planning system (TPS) for 8 MV photon beams of a SL25 linear accelerator in realistic patient models. Simulation of the treatment head involves modeling of the main components of the treatment head that have influence on the absorption and scattering of radiation. Simulation of the accelerator was done in two parts to minimize the simulation time. Analysis of the data generated by the BEAM code was carried out using BEAMDP, another subsidiary of the BEAM code. We calculated the beam characteristics which are difficult to measure experimentally, such as angular distributions, spectral distributions, planar fluence and planar energy fluence at a plane located just above the jaws of the treatment head. The phase space files at the isocenter were used as source input for DOSXYZ, a MC code to calculate 3D dose distributions in water or CT based phantoms. The DOSXYZ code was used to calculate depth dose and cross plane profiles in a water phantom. The data obtained with Monte Carlo methods were compared with that obtained by ionization chamber measurements. Depth dose and cross plane profiles obtained by Monte Carlo methods and ionization chamber measurements generally agreed within 2%. We created patient models from CT data of real patients using the CTCREATE option of the DOSXYZ program. Dose distributions for a number of field sizes and different anatomical sites were calculated with the DOSXYZ code and compared with corresponding dose distributions calculated by the TPS. The modified BATHO and ETAR inhomogeneity correction methods used in the TPS were evaluated. Results show that Monte Carlo methods can accurately reproduce ion chamber measurements in a water phantom. Monte Carlo techniques are very useful for evaluating the accuracy of dose distributions generated by treatment planning systems in patient based models where measurements are impossible. The BATHO and ETAR methods showed comparable results to the Monte Carlo results. This could be due to the inefficiency of the method (visualization of the dose distributions) that we used for the comparison of the results. A more quantitative method like the use of the dose difference volume histogram could give a more comprehensive evaluation.Item Open Access Development of a Monte Carlo simulation method for the evaluation of dose distribution calculations of radiotherapy treatment planning systems(University of the Free State, 1999-11) Du Plessis, Frederik Carl Philippus; Willemse, C. A.; Lotter, M. G.In this study a method is described whereby the dose distributions calculated by any treatment planning system (TPS) could be evaluated using dose distributions calculated with Monte Carlo simulations. The Monte Carlo dose simulations can be regarded as the golden standard. The method developed in this study involved the Monte Carlo simulation of a Philips SL75/14 based generic accelerator with the BEAM code. This was done to obtain beam information stored in phase space files that were characteristic of the generic accelerator. This beam data were then used for the simulation of dose distributions in a mathematical water phantom using the Monte Carlo code, DOSXYZ. The same beam data were used to generate the data base for the TPS that uses it for dose calculations in CT based patient models. The BATHO and ETAR inhomogeneity correction algorithms implemented on a CADPLAN TPS were evaluated. The CT slices that make up the patient model, on the TPS, were transformed to material data. Each of these materials (57 in total) covered a discrete CT interval in a total CT number range of 3000 CT numbers. Each of the 57 materials was represented in the preprocessor code (pEGS4) to allow dose simulations in realistic patient models with the DOSXYZ code. Dose distributions were calculated in a maxillary sinus (head), lung and prostate patient for photon beams with size 2x2, 5x5 and 10xlO cm2 . These dose distributions were calculated on the TPS using the BATHO and ET AR methods. The DOSXYZ dose distributions were scaled to the TPS calculated dose distributions by normalization to the dose in water at 2 cm depth on the beam central axis. Dose difference volume histograms, percentage depth dose curves and 2D dose distributions were obtained to evaluate these dose distributions. The BATHO and ET AR methods cannot model lateral and longitudinal electron transport through complex media. These effects were apparent in large inhomogeneities such as in the lung model where the Monte Carlo dose simulation gave a wider beam penumbra for the large field, and in the deviation of the TPS dose distributions in these regions for the small field size. The method developed in this study could also be applied to any IPS that uses more sophisticated models. Manufacturers of IPS's in particular could use the methods described in this study to evaluate their dose calculation algorithms. Key words: Monte Carlo, CT based patient model, DOSXYZ, BEAM, Treatment planning, dose distributions, lung, maxillary sinus, water phantom, lateral electron transport, TPS, inhomogeneity.Item Open Access The assessment of potential radiation hazards from gold mines in the Free State Goldfields to members of the public(University of the Free State, 1998-11) Ellis, Jozua Francois; Botha, J. C.; Van Aswegen, A.The gold mines In the Free State Goldfields extract and process ore from underground, which contains naturally radioactive uranium and its associated decay products. This assessment aimed to cost effectively determine the major potential radiation hazards to the public from the gold mines in the area. The potential exposure sources from the mines are radon gas, radioactive dust, contaminated water and external gamma radiation. The assessment focussed mainly on the public's potential exposure to radon gas emanating from tailings dams, waste rock dumps and upcast shafts from underground workings. The rate of radon emanation from the dams was measured using several different techniques, and the potential dispersion of the radon was modelled using internationally accepted modelling codes and local weather data for the Free State Goldfields. A maximum potential contribution to the natural background radon levels of 6 Bq m-3 was calculated. This is a small increment to the background levels in the order of 25 to 35 Bq m-3". Environmental measurements of outdoor radon concentrations confirmed the modelling results to the extent that no significantly high radon concentration could be detected in the environment. Background radon levels in towns outside the Free State Goldfields are in the same order as those measured around the mines. Similar environmental measurements of airborne dust and water sources around the mines indicated relatively low levels of radiation. A conservative estimate of the total potential exposure of the public in the Free State Goldfields is in the order of 130 to 250 µSv/a. This can be interpreted as well within the internationally accepted public dose limit of 1000 µSv/a.Item Open Access Development and validation of a molecular assay and evaluation of the GeneXpert® MTB/RIF assay for the rapid detection of genital tuberculosis(University of the Free State, 2016-02) Sokhela, Maxwell Cebolenkosi; Goedhals, D.; Hoosen, A. A.English: Tuberculosis (TB) is a communicable disease which is caused by the bacterium Mycobacterium tuberculosis (MTB). According to the World Health Organization, globally in 2015 there were 10.4 million new cases and 1.4 million deaths due to TB. TB is one of the leading causes of death in South Africa resulting in approximately 8.4% of deaths in 2015. The most common manifestation of TB involves the lungs, defined as pulmonary TB (PTB), while TB affecting other organs is defined as extrapulmonary TB (EPTB). EPTB accounts for only 20% of all TB cases in human immunodeficiency virus negative individuals. Approximately 1.8% of all TB cases have a genitourinary site, with the prevalence of genital TB (GTB) in South Africa reported to range from 6.2-21.0%. One of the leading symptoms of GTB in females is infertility, usually resulting from the involvement of the fallopian tubes and endometrium. Approximately 40-80% of women with GTB will become infertile. The detection of microorganisms through microscopy is the oldest technique for laboratory diagnosis. While microscopy is rapid and inexpensive, it requires a high bacterial load which is not present in paucibacillary EPTB samples. Culture of MTB is widely regarded as the gold standard for TB diagnosis. While culture has a long turnaround time, culture remains important since it is more sensitive than microscopy. In addition, growth is required for species identification, drug susceptibility testing and genotyping of cultured organisms may be useful for epidemiological studies. Little is known regarding which technique is best for the detection of GTB from clinical samples apart from culture. Molecular based techniques hold the promise of a more rapid and accurate diagnosis of EPTB. The aim of this project was the development and validation of an in-house nested PCR assay and the validation of the GeneXpert® MTB/RIF (GeneXpert) assay for the laboratory diagnosis of GTB. In total 54 samples were submitted for GTB screening from women being investigated for infertility at the Unit for Human Reproduction, Universitas Academic Hospital, Bloemfontein. This included 44 endometrial tissue samples and 10 menstrual fluid samples. All samples underwent testing with the GeneXpert, the in-house nested PCR and culture. The nested PCR was designed targeting the insertion sequence element 6110 (IS6110) found in members of the MTB complex. The analytical sensitivity/limit of detection (LOD) for the GeneXpert was determined to be 250pg while the LOD for the nested polymerase chain reaction (PCR) was 62.5fg. Both assays displayed excellent analytical specificity by discriminating TB deoxyribose nucleic acid (DNA) from other bacterial and nontuberculous mycobacterial DNA. The diagnostic sensitivity and specificity was determined using culture as the reference method. Culture was able to detect GTB in 2 of the 54 samples including one menstrual fluid and one endometrial tissue sample, thus indicating a GTB prevalence of 3.7%. The GeneXpert detected 1 of the 54 samples as positive indicating a sensitivity of 50% and a specificity of 100%. The nested PCR detected both positive samples resulting in a sensitivity and specificity of 100%. The GeneXpert obtained a positive predictive value (PPV) of 100% and a negative predictive value (NPV) of 98.1%, while the nested PCR obtained a PPV and NPV of 100%. The two GTB isolates underwent genotyping using spoligotyping and mycobacterial interspersed repetitive unit – variable number of tandem repeats (MIRU-VNTR). The menstrual fluid isolate was characterised as a Beijing strain and the endometrial tissue isolate as an X3 strain. The nested PCR showed a greater sensitivity than the GeneXpert as a result of the better LOD. Despite this, both techniques could be implemented for GTB screening in combination with culture. Screening of menstrual fluid samples using the GeneXpert assay would be well suited for GTB screening in resource limited areas.Item Open Access Monte Carlo study on megavolt x-ray therapy for development of suitable targets for the evaluation of nano particle dose enhancement(University of the Free State, 2016-09) Mutsakanyi, Stalyn; Du Plessis, F. C. P.English: INTRODUCTION: Radiation dose enhancement with nanoparticles is a treatment technique involving the irradiation of tumour seeded with high atomic number (high Z) material. This work describes the generation of x-ray beams using a 6 MeV Elekta Precise linac head using low-Z Bremsstrahlung target materials, water and carbon combined with tungsten. The aim of the study was to simulate photon energy spectra appropriate for high-Z nanoparticles dose enhancement in tumour using EGSnrc MC codes. MATERIALS AND METHOD: BEAMnrc Monte Carlo (MC) code successfully modelled the treatment head components of a flattening filter free 6 MV Elekta Precise linear accelerator. Simulations were run using suitable histories to generate high energy x-ray beams of differing quality from electron spectra obtained using 6 MeV electron beam. Water and carbon layers were the primary target which were inserted in the path of the 6 MeV electron pencil beam before it hits the tungsten Bremsstrahlung target to act as moderators that slow down electron before they hit a tungsten layer. The electron spectra obtained just after the primary target was used as the incident beam to the tungsten target which acts as the secondary target to generate x-ray photon beams. Therefore the x-ray beam source target was either water/tungsten or carbon/tungsten combination. Different photon spectra were obtained for investigation in nanoparticles (NPs) based photon therapy. An original linac using a normal tungsten target of 0.3 cm thickness was also simulated to benchmark the results. The photon spectra obtained below X,Y jaws were used as input sources in the DOSXYZnrc MC code to simulate dose distribution in water and a patient CT phantom. The simulations were carried out using source 2 in DOSXYZnrc. A 40 x 40 x 40 cm3 water phantom was simulated at 100 cm SSD using a range of field sizes to characterize the beams. The phantom had voxel size of 0.2 × 0.2 × 0.2 cm3. The photon beams were characterised in terms of percentage depth doses and beam profiles. These x-ray beams were then used to quantify the variation of tumour dose enhancement in a constructed patient CT phantom. The prostate tumour was used as the planning target volume (PTV). The PTV composition was either a tumour only or a tumour volume seeded with atoms of gold nanoparticles with concentration of 7mg/g of tumour. These tumour/NPs model was manually drawn on to the CT dataset from actual CT images of the patient using MCSHOW graphical user interface (GUI). The tumour composition was made part of the patient CT data set using a locally-developed Interactive Data Language (IDL) code that converts the density of the drawn volume into the desired tumour density. The 3DCRT was used as the treatment strategy and 4, 5 and 6 field beams were investigated. With this model, we were able to estimate more accurately the effect of altered beams on NPs radiation dose enhancement. For both simulations using BEAMnrc and DOSXYZnrc the electron cut-off energy (ECUT) and photon cut-off energy (PCUT) was 0.521 MeV and 0.01 MeV respectively. The number of histories was chosen so that the statistical uncertainty along the CAX had an average value of 1% at 0 – 30 cm depth. RESULTS AND CONCLUSION: The results showed that the use of electron moderators in generating x-ray beams for use in NPs seeded tumours can lead to a significant dose enhancement. Photon spectra obtained with water/tungsten or carbon/tungsten Bremsstrahlung targets combinations showed significant changes at various target thickness. There is a significant dependence of dose enhancement factors (DEF) on the mean energy of the x-ray beams as well as the target thickness. DEFs ranging from 0.05% to 7.5% were obtained at various Bremsstrahlung target combinations. Based on the results, carbon is more efficient at moderating the electron beam to generate photon beams for dose enhancement at lower thickness (approximately 1.4 cm) compared to water (approximately 2.5 cm), although water can just be as good at larger thickness. At these thicknesses the mean photon beam energy is approximately 0.4 MV. In summary, the results of this work indicate that the use of photon beams from low-Z Bremsstrahlung targets as electron slowing down medium could enable significant clinical dose enhancement during external beam radiotherapy for NPs seeded tumours. MC techniques showed to be valuable tools for dose calculations in both water and patient CT phantom.Item Open Access Monte Carlo evaluation of the dose perturbation effect of various hip prostheses during pelvic megavoltage photon radiotherapy(University of the Free State, 2016-11) Mahuvava, Courage; Du Plessis, F. C. P.English: Introduction: Hip prostheses (HPs) are routinely used in hip augmentation surgery to replace painful or dysfunctional hip joints, especially in the elderly population. A number of patients with HPs are undergoing pelvic radiotherapy (RT) for localised prostate cancer. However, radiographic discrepancies between high-density and high-atomic-number (Z) inserts and surrounding tissue may cause considerable dose perturbations within the target volume and in regions where tissues interface with the prosthetic device. Furthermore, conventional treatment planning systems (TPSs) do not accurately predict dose effects incurred around metallic implants. Therefore, concerns regarding dose inhomogeneities near the prosthesis always arise, especially in patients with bilateral hip prostheses (bHPs) who require teletherapy of prostate cancer, where the tumour typically lies between the prostheses. The aim of this study was to evaluate the dosimetric effect of various HPs during 3D conformal prostate RT using Monte Carlo (MC) simulations. Materials and methods: The MC radiation transport simulation user-code BEAMnrc was used to simulate an Elekta Precise linear accelerator (linac) head, based on the manufacturer’s specifications. The MC linac model was validated by comparing dosimetric features including depth dose and dose profile data simulated in a cubic water tank (WT) with measured values. DOSXYZnrc was used to calculate 3D absorbed dose distributions in a CT based phantom (patient model) with and without HPs. Simulations were performed for 6, 10, 15 and 20 MV conformal photon beams using different beam arrangements. Three treatment plans were generated by XiO TPS and incorporated into MC simulations: a four–field (4F) box plan, a five–field (5F) plan and a six–field (6F) plan. The planning target volume (PTV) was generated by a 1 cm expansion of the prostate alone. The HP materials used were stainless steel (SS316L), titanium (Ti6Al4V) and ultra-high-molecular-weight-polyethylene (UHMWPE). These prosthetic models were manually drawn into the CT dataset from actual CT images of the patient pelvis using MCSHOW graphical user interface (GUI). The prosthesis was made part of the patient using a locally-developed Interactive Data Language (IDL) code that converts the density of the drawn volume into the desired HP material density. Both unilateral and bilateral models were considered in the simulations and dose perturbation factors (DPFs) were calculated on the proximal and distal interfaces of the implant. The dose reduction in the PTV as well as the dose to critical organs was also evaluated. Results: Results indicated that the central axis depth dose within and beyond the inhomogeneity drops significantly due to beam attenuation. For patients with bHPs, the dose contribution from lateral ports at 6 MV was attenuated by up to 23% and 17% for SS316L and Ti6Al4V, respectively. For a unilateral HP (uHP), the respective dose attenuations were 19% and 12%. The dose perturbation was always < 1% for a patient fitted with UHMWPE. Up to 38% dose increase was found at the proximal bone–HP interface due to backscattered electrons from the metal implant. There was a weak dependence of dose distribution on beam energy at the target isocenter, with the maximum dose reduction ranging only from 22.8 to 16.9% from 6 to 20 MV in a patient with bilateral steel HPs. However, interface effects were more pronounced at higher beam energies. However, increasing the number of treatment beams improved the plan quality. The greatest PTV dose perturbation was observed in a 4F box and lowest in a 6F plan. Production of scatter radiation was found to be larger for backscatter compared to forward scatter in this study. Conclusions: The dose perturbation effect of metallic HPs is significant and must be taken into account during treatment planning. UHMWPE poses no significant dose perturbation in the shadow of the implant and on the interface with tissue or bone. The use of MC–based TPSs is recommended for treatments using beam portals passing through HPs. MCSHOW allows the addition of HP contours in the virtual phantom from CT dataset of a patient without a HP. This allows one to carry out MC calculations for several implant models without metal artefacts. Results also highlight the significant influence of the implant’s composition and the beam position relative to the HP as well as beam energy on the dose distribution. Increasing the beam energy may help overcome the attenuation effects of metallic HPs and to improve target coverage. Therefore this study recommends plans with a larger number of beams that would allow avoiding the hip inhomogeneity in order to effectively compensate for dose attenuated in fields passing through HPs. 1It is also evident from the results that the shadowing effect is density-dependent, and its maximum value is for the SS316L HP. A more sophisticated, non-coplanar beam orientation may be necessary to avoid the HPs whilst sparing organs at risk (OARs) and giving sufficient target coverage.Item Open Access Evaluation 99mTc and 123I quantification using SPECT/CT(University of the Free State, 2015-02) Mongane, Modisenyane Simon; Van Staden, J. A.; Du Raan, H.English: A review of Single Photon Emission Tomography (SPECT) quantification shows that different protocols and phantoms are available to evaluate SPECT quantification accuracy. This study was necessitated by the lack of standardized protocols and the widespread use of a variety of non-standard phantoms. The aim of this work was to evaluate the influence of the geometry of a radionuclide distribution on SPECT quantification accuracy for 99mTc and 123I isotopes in an abdominal phantom. In order to achieve the aim, the following steps were taken: The preparatory phase of the study was to design and construct an abdominal phantom, verify the accuracy of the attenuation coefficients obtained with the Computed Tomography (CT) scanner, determine the accuracy of the source calibrator used in this study and then obtain a calibration factor in order to convert image counts to activity. During the quantification phase SPECT data were acquired, the influence of not applying scatter correction explicitly was evaluated and the final quantification was performed using the proposed standard clinical reconstruction protocols. The influence of different tumour sizes and locations in the abdominal phantom relative to a high uptake organ on the quantification accuracy was evaluated. Finally, parameters in the Ordered Subset Expectation Maximization (OSEM) reconstruction protocol were altered in order to investigate the influence of number of subsets and iterations on the quantified data. The non-standard Density Phantom with five different compounds was used for the verification of the 99mTc and 123I attenuation coefficients. The percentage difference between the measured and theoretical attenuation coefficients values were < 3%, except for Polystyrene (85% and 65% respectively). The SPECT calibration factor was determined for both 99mTc (11.0 ± 1.3 cpm/kBq) and 123I (10.8 ± 0.3 cpm/kBq) using the Cylindrical Phantom. The in-house built Abdominal Phantom was used to evaluate the tumour activity quantification accuracy. The quantification accuracy of 99mTc and 123I was found to change significantly (p < 0.05) as a function of tumour size after corrections for “spill out” counts due to the partial volume effect, scatter and attenuation were applied. On the other hand, there was no significant difference in the quantification accuracy (p > 0.05) for each tumour at different tumour-liver distances when appropriate scatter and attenuation corrections were applied. The influence of OSEM parameters showed no dependence on the tumour-liver distance and no significant difference (p > 0.05) between quantification with background activity as compared to no background activity. In conclusion, the study showed that the quantification accuracy for 99mTc and 123I was comparable to other published studies. It was found that the tumour quantification accuracy is not influenced by proximity of high uptake organs when appropriate correction factors were applied. Tumour size influenced the accuracy of SPECT quantification for both radionuclides. The results of this study also showed that at least 128 Maximum Likelihood Expectation Maximization (MLEM) equivalent iterations were needed during iterative reconstruction to achieve convergence and consistent SPECT quantification accuracy. Finally, it is recommended that the evaluated quantification protocol may be used in our nuclear medicine clinic for 99mTc and 123I quantification.Item Open Access Optimisation of delivery efficiency in prostate intensity modulated radiotherapy planning(University of the Free State, 2016-02) Fourie, Nicola Sieglinde; Muhammed, Omer Abdul-Aziz Ali; Rae, William Ian DuncombeEnglish: Evidence that supports dose escalation for prostate cancer is growing and with Intensity Modulated Radiation Therapy (IMRT) higher conformal target doses can be delivered. With more segments and higher monitor units (MU’s), target conformity can be improved, however this results in longer delivery times, which makes it difficult to ensure accurate dose delivery, as intra-fractional as well as target movement plays an increasing role. Evidence from the literature indicates that secondary radiation-induced cancer risk is proportional to the beam-on time (thus the MU’s). Improvements in IMRT delivery efficiency while maintaining plan quality can be achieved by reducing the complexity of an IMRT plan. This can be done by changing the optimization parameters during the optimization process. Less “complex” prostate IMRT plans will require fewer MU’s by using less segments resulting in shorter delivery times and therefore reduced risk of secondary cancers. The goal of this study was to recommend a set of optimization parameter values that will improve the delivery efficiency of prostate IMRT treatment plan while maintaining plan quality. Fifteen clinical prostate IMRT plans (15 MV), already used for treatment, were re-optimized, using a XiO treatment planning system (TPS). Changes in total MU’s and segments were evaluated for changes in some of the optimization parameter values. Eleven optimization parameters (some of them used more than once with different values) were used to generate 15 new IMRT combination plans (ICP’s) for each patient for both 6 and 15 MV, resulting in 450 plans being assessed. One parameter was changed at a time while all other variables were kept constant. Plan quality was evaluated in terms of four variables: MU, number of segments, homogeneity index and conformity index while the delivery efficiency was evaluated in terms of delivery time. To our knowledge no time delivery model has been proposed for a Siemens® ARTISTETM Linear Accelerator (Linac). Using the principles given in the literature we derived such a time delivery model by adding the radio frequency wave component and Multi Leaf Collimator delay time. K-means clustering was then used to analyse the data in terms of the five variables and the top 10 ICP’s in 3 patients in terms of a faster more conformal, delivered plan were identified. To confirm the delivery efficiency and accuracy, the fluences of these top 10 ICP’s were measured on a Siemens® ARTISTETM Linac with the step and shoot method and compared to the treatment planning system’s fluences. The evaluation criteria chosen were 3% and 3 mm, distance to agreement. A 3 dimensional dose volume histogram program was used to determine the percentage pass rates on the planned target volumes and the organs at risk. The optimization parameters such as the minimum MU’s per segment, intensity level, minimum segment size and minimum segment area; demonstrated the greatest influence on the total number of segments, while the total MU’s was most greatly influenced by the filters and intensity level optimization parameter. Controversy exists regarding which energy should be used, 6 MV or 15 MV, when treating prostate cancer. Both energies were considered here during the optimization process and it was concluded that the optimization parameters are not greatly influenced by the beam energy. However, it was seen that beam arrangement has an influence on optimization parameter behaviour. A limitation of this study is that the beam angle distribution was not investigated. Thus recommendations could be made in terms of which ICP demonstrated the most improved delivery efficiency of a prostate IMRT treatment plan while maintaining plan quality. The optimisation parameter which was introduced to the optimization process was a General High filter. Gaining knowledge about the behaviour of the optimization parameters during optimization makes it easier to advise and assist treatment planners preparing complex IMRT plans.Item Open Access Manual quantification for right and left ventricular function using cardiac magnetic resonance imaging(University of the Free State, 2012-07) Willemse, Nanette; Rae, W. I. D.; Herbst, C. P.Abstract not availableItem Open Access Development of a particle source model for a synergy linear accelerator to be used in Monte Carlo radiation dose calculations for cancer therapy(University of the Free State, 2014-05) Van Eeden, Dete; Du Plessis, F. C. P.English: In oncology patients are treated for cancer with various methods such as surgery, chemo therapy and radiation therapy. Accurate radiation treatment planning and dose delivery to the tumour is necessary for the successful outcome of cancer treatment. In order to achieve this goal accurate radiation dose calculation codes must be utilized. EGSnrc based Monte Carlo (MC) codes such as BEAMnrc and DOSXYZnrc have been developed for just this purpose. The problem that arises in using these MC codes is that they lack suitable x-ray beam source models. These models must be accurate in order to replicate the true clinical x-ray beam emanating from the linear accelerator. One such machine for which radiation source data must be derived is currently being used at the Oncology department in Universitas Hospital Annex. It is desirable to model this linear accelerator in order to perform MC based dose calculations for radiation treatment. The use of MC based dose calculations is certainly not new in the radiation physics environment. Various authors have studied the replication of radiation beam characteristics using source models to simulate the phase-space parameters of particles produced by the linear accelerator. These parameters include the charge, energy, direction, and position of each particle as it crosses a certain reference plane below the linear accelerator. An accurate source model should be able to re-generate particles with the exact set of above-mentioned parameters as would be produced by the real linear accelerator. Sources can be very simple such as a single point from which the particles are radiating with a single invariant energy spectrum. Studies have shown that these beam models can yield accurate beam data over relatively small field sizes and is not general enough to use over a whole range of clinically useful field sizes. A graphical user interface (GUI) was developed that can assist in the construction of the source model. The source model can describe energy and fluence distributions for photons and electrons as separate point sources each with their own SSD. The accuracy of the model was validated by comparing simulated profiles with measured data for an Elekta Synergy linear accelerator. The modified Schiff formula was used to derive the bremsstrahlung spectra emanating from the target. The x-ray fluence Gaussian distribution consisted of the primary fluence from the target, which was modified by the primary collimator, secondary collimators as well as the multileaf collimators. The truncation and beam scatter caused by the face of the collimators were modelled with error functions. Exponential functions were used to model off-axis collimator transmission. Profiles and percentage depth dose curves were obtained with the source for square field sizes of 1 × 1 cm2 up to a 40 × 40 cm2. Offset fields for 10 × 10 cm2, 15 × 15 cm2 and 20 × 20 cm2, rectangular fields as well as wedged fields were included. Irregular field shapes were simulated to evaluate the source model‘s capability of reproducing complex treatment fields. Film dose verification was done in an anthropomorphic Rando® phantom and compared with the MC source model for 6 MV x-ray beams. A criterion of 2% / 2 mm was used to compare MC data and measured data. This study demonstrated that a diversity of field sizes and percentage depth dose curves can be modelled within 2% / 2 mm. The model can replicate irregular field sizes used for complex treatments. Minor discrepancies were found for the relative dose comparisons between the MC and film data for the anthropomorphic phantom.