Technical note: SpekPy Web-online x-ray spectrum calculations using an interface to the SpekPy toolkit.

Knowledge of the photon spectrum emitted from an x-ray tube is frequently needed in imaging and dosimetry contexts. As the spectrum characteristics are influenced by several parameters and routine measurement of a spectrum is often impractical, a variety of software programs have been developed over the decades for convenient calculations. SpekPy is a state-of-the-art software package containing several spectrum models, and was created to estimate photon spectra originating from x-ray tubes using a small set of input parameters (e.g., anode material, anode angle, tube potential, filtration, etc.). SpekPy is distributed as a Python toolkit and is available free of charge. The toolkit does, however, lack a graphical user interface and a user is required to write a Python script to make use of it. In this work this limitation is addressed by introducing a web application called SpekPy Web: a graphical user interface together with an application programmable interface (API). These developments both make the SpekPy spectrum models accessible to a broader set of users and increases the ease of use for existing users. SpekPy Web is hosted at: https://spekpy.smile.ki.se. The functionality of the software is demonstrated, using its API, by estimating first half-value layers (HVLs) for 15 standard beam qualities from the International Bureau of Weights and Measures (BIPM). The estimated HVLs were found to all be within 3.5% agreement when compared to experimental values, with an average calculation time of 2.5 s per spectrum. half-value-layer, software, x-ray spectrum.

Accuracy of the dose-shift approximation in estimating the delivered dose in SBRT of lung tumors considering setup errors and breathing motions.

INTRODUCTION: Geometrical uncertainties can result in a delivered dose to the tumor different from that estimated in the static treatment plan. The purpose of this project was to investigate the accuracy of the dose calculated to the clinical target volume (CTV) with the dose-shift approximation, in stereotactic body radiation therapy (SBRT) of lung tumors considering setup errors and breathing motion. The dose-shift method was compared with a beam-shift method with dose recalculation. MATERIAL AND METHODS: Included were 10 patients (10 tumors) selected to represent a variety of SBRT-treated lung tumors in terms of tumor location, CTV volume, and tumor density. An in-house developed toolkit within a treatment planning system allowed the shift of either the dose matrix or a shift of the beam isocenter with dose recalculation, to simulate setup errors and breathing motion. Setup shifts of different magnitudes (up to 10 mm) and directions as well as breathing with different peak-to-peak amplitudes (up to 10:5:5 mm) were modeled. The resulting dose-volume histograms (DVHs) were recorded and dose statistics were extracted. RESULTS: Generally, both the dose-shift and beam-shift methods resulted in calculated doses lower than the static planned dose, although the minimum (D98%) dose exceeded the prescribed dose in all cases, for setup shifts up to 5 mm. The dose-shift method also generally underestimated the dose compared with the beam-shift method. For clinically realistic systematic displacements of less than 5 mm, the results demonstrated that in the minimum dose region within the CTV, the dose-shift method was accurate to 2% (root-mean-square error). Breathing motion only marginally degraded the dose distributions. CONCLUSIONS: Averaged over the patients and shift directions, the dose-shift approximation was determined to be accurate to approximately 2% (RMS) within the CTV, for clinically relevant geometrical uncertainties for SBRT of lung tumors.

Comparison of CT number calibration techniques for CBCT-based dose calculation.

PURPOSE: The aim of this work was to compare and validate various computed tomography (CT) number calibration techniques with respect to cone beam CT (CBCT) dose calculation accuracy. METHODS: CBCT dose calculation accuracy was assessed for pelvic, lung, and head and neck (H&N) treatment sites for two approaches: (1) physics-based scatter correction methods (CBCTr); (2) density override approaches including assigning water density to the entire CBCT (W), assignment of either water or bone density (WB), and assignment of either water or lung density (WL). Methods for CBCT density assignment within a commercially available treatment planning system (RSauto), where CBCT voxels are binned into six density levels, were assessed and validated. Dose-difference maps and dose-volume statistics were used to compare the CBCT dose distributions with the ground truth of a planning CT acquired the same day as the CBCT. RESULTS: For pelvic cases, all CTN calibration methods resulted in average dose-volume deviations below 1.5 %. RSauto provided larger than average errors for pelvic treatments for patients with large amounts of adipose tissue. For H&N cases, all CTN calibration methods resulted in average dose-volume differences below 1.0 % with CBCTr (0.5 %) and RSauto (0.6 %) performing best. For lung cases, WL and RSauto methods generated dose distributions most similar to the ground truth. CONCLUSION: The RSauto density override approach is an attractive option for CTN adjustments for a variety of anatomical sites. RSauto methods were validated, resulting in dose calculations that were consistent with those calculated on diagnostic-quality CT images, for CBCT images acquired of the lung, for patients receiving pelvic RT in cases without excess adipose tissue, and for H&N cases.

Quality control of CT systems by automated monitoring of key performance indicators: a two-year study.

The purpose of this study was to develop a method of performing routine periodical quality controls (QC) of CT systems by automatically analyzing key performance indicators (KPIs), obtainable from images of manufacturers' quality assurance (QA) phantoms. A KPI pertains to a measurable or determinable QC parameter that is influenced by other underlying fundamental QC parameters. The established KPIs are based on relationships between existing QC parameters used in the annual testing program of CT scanners at the Karolinska University Hospital in Stockholm, Sweden. The KPIs include positioning, image noise, uniformity, homogeneity, the CT number of water, and the CT number of air. An application (MonitorCT) was developed to automatically evaluate phantom images in terms of the established KPIs. The developed methodology has been used for two years in clinical routine, where CT technologists perform daily scans of the manufacturer's QA phantom and automatically send the images to MonitorCT for KPI evaluation. In the cases where results were out of tolerance, actions could be initiated in less than 10 min. 900 QC scans from two CT scanners have been collected and analyzed over the two-year period that MonitorCT has been active. Two types of errors have been registered in this period: a ring artifact was discovered with the image noise test, and a calibration error was detected multiple times with the CT number test. In both cases, results were outside the tolerances defined for MonitorCT, as well as by the vendor. Automated monitoring of KPIs is a powerful tool that can be used to supplement established QC methodologies. Medical physicists and other professionals concerned with the performance of a CT system will, using such methods, have access to comprehensive data on the current and historical (trend) status of the system such that swift actions can be taken in order to ensure the quality of the CT examinations, patient safety, and minimal disruption of service.

Stopping-power ratio estimation for proton radiotherapy using dual-energy computed tomography and prior-image constrained denoising.

BACKGROUND: Dual-energy computed tomography (DECT) is a promising technique for estimating stopping-power ratio (SPR) for proton therapy planning. It is known, however, that deriving electron density (ED) and effective atomic number (EAN) from DECT data can cause noise amplification in the resulting SPR images. This can negate the benefits of DECT. PURPOSE: This work introduces a new algorithm for estimating SPR from DECT with noise suppression, using a pair of CT scans with spectral separation. The method is demonstrated using phantom measurements. MATERIALS AND METHODS: An iterative algorithm is presented, reconstructing ED and EAN with noise suppression, based on Prior Image Constrained Denoising (PIC-D). The algorithm is tested using a Siemens Definition AS+ CT scanner (Siemens Healthcare, Forchheim, Germany). Three phantoms are investigated: a calibration phantom (CIRS 062M), a QA phantom (CATPHAN 700), and an anthropomorphic head phantom (CIRS 731-HN). A task-transfer function (TTF) and the noise power spectrum are derived from SPR images of the QA phantom for the evaluation of image quality. Comparisons of accuracy and noise for ED, EAN, and SPR are made for various versions of the algorithm in comparison to a solution based on Siemens syngo.via Rho/Z software and the current clinical standard of a single-energy CT stoichiometric calibration. A gamma analysis is also applied to the SPR images of the head phantom and water-equivalent distance (WED) is evaluated in a treatment planning system for a proton treatment field. RESULTS: The algorithm is effective at suppressing noise in both ED and EAN and hence also SPR. The noise is tunable to a level equivalent to or lower than that of the syngo.via Rho/Z software. The spatial resolution (10% and 50% frequencies in the TTF) does not degrade even for the highest noise suppression investigated, although the average spatial frequency of noise does decrease. The PIC-D algorithm showed better accuracy than syngo.via Rho/Z for low density materials. In the calibration phantom, it was superior even when excluding lung substitutes, with root-mean-square deviations for ED and EAN less than 0.3% and 2%, respectively, compared to 0.5% and 3%. In the head phantom, however, the SPR accuracy of the PIC-D algorithm was comparable (excluding sinus tissue) to that derived from syngo.via Rho/Z: less than 1% error for soft tissue, brain, and trabecular bone substitutes and 5-7% for cortical bone, with the larger error for the latter likely related to the phantom geometry. Gamma evaluation showed that PIC-D can suppress noise in a patient-like geometry without introducing substantial errors in SPR. The absolute pass rates were almost identical for PIC-D and syngo.via Rho/Z. In the WED evaluations no general differences were shown. CONCLUSIONS: The PIC-D DECT algorithm provides scanner-specific calibration and tunable noise suppression. It is vendor agnostic and applicable to any pair of CT scans with spectral separation. Improved accuracy to current methods was not clearly demonstrated for the complex geometry of a head phantom, but the suppression of noise without spatial resolution degradation and the possibility of incorporating constraints on image properties, suggests the usefulness of the approach.

Spatial resolution, noise properties, and detectability index of a deep learning reconstruction algorithm for dual-energy CT of the abdomen.

BACKGROUND: Iterative reconstruction (IR) has increasingly replaced traditional reconstruction methods in computed tomography (CT). The next paradigm shift in image reconstruction is likely to come from artificial intelligence, with deep learning reconstruction (DLR) solutions already entering the clinic. An enduring disadvantage to IR has been a change in noise texture, which can affect diagnostic confidence. DLR has demonstrated the potential to overcome this issue and has recently become available for dual-energy CT. PURPOSE: To evaluate the spatial resolution, noise properties, and detectability index of a commercially available DLR algorithm for dual-energy CT of the abdomen and compare it to single-energy (SE) CT. METHODS: An oval 25 cm x 35 cm custom-made phantom was scanned on a GE Revolution CT scanner (GE Healthcare, Waukesha, WI) at two dose levels (13 and 5 mGy) and two iodine concentrations (8 and 2 mg/mL), using three typical abdominal scan protocols: dual-energy (DE), SE 80 kV (SE-80 kV) and SE 120 kV (SE-120 kV). Reconstructions were performed with three strengths of IR (ASiR-V: AR0%, AR50%, AR100%) and three strengths of DLR (TrueFidelity: low, medium, high). The DE acquisitions were reconstructed as mono-energetic images between 40 and 80 keV. The noise power spectrum (NPS), task transfer function (TTF), and detectability index (d') were determined for the reconstructions following the recommendations of AAPM Task Group 233. RESULTS: Noise magnitude reductions (relative to AR0%) for the SE protocols were on average (-29%, -21%) for (AR50%, TF-M), while for DE-70 keV were (-28%, -43%). There was less reduction in mean frequency (fav ) for DLR than for IR, with similar results for SE and DE imaging. There was, however, a substantial change in the NPS shape when using DE with DLR, quantifiable by a marked reduction in the peak frequency (fpeak ) that was absent in SE mode. All protocols and reconstructions (including AR0%) exhibited slight to moderate shifts towards lower spatial frequencies at the lower dose (<12% in fav ). Spatial resolution was consistently superior for DLR compared to IR for SE but not for DE. All protocols and reconstructions (including AR0%) showed decreased resolution with reduced dose and iodine concentration, with less decrease for DLR compared to IR. DLR displayed a higher d' than IR. The effect of energy was large: d' increased with lower keV, and SE-80 kV had higher d' than SE-120 kV. Using DE with DLR could provide higher d' than SE-80 kV at the higher dose but not at lower dose. CONCLUSIONS: DE imaging with DLR maintained spatial resolution and reduced noise magnitude while displaying less change in noise texture than IR. The d' was also higher with DLR than IR, suggesting superiority in detectability of iodinated contrast. Despite these trends being consistent with those previously established for SE imaging, there were some noteworthy differences. For DE imaging there was no improvement in resolution compared to IR and a change in noise texture. DE imaging with low keV and DLR had superior detectability to SE DLR at the high dose but was not better than SE-80 kV at low dose.

Ventriculoperitoneal shunt survey in pediatric patients in an emergency setting: a comparison between ultra-low-dose CT using a tin filter and plain radiography.

OBJECTIVE: The aim of this study was to investigate the diagnostic performance, effective radiation dose, and examination time of ventriculoperitoneal shunt evaluation using full-body ultra-low-dose CT (ULD CT) with a tin filter compared with digital plain radiography in a pediatric population. METHODS: A retrospective cross-sectional study was conducted in an emergency setting. Data from 143 children were collected. Sixty were examined with ULD CT with a tin filter and 83 with digital plain radiography methods. Effective doses and times were compared between the two methods. Two observers in pediatric radiology evaluated the patient images. Clinical findings and results from shunt revision, if it was performed, were used to evaluate the diagnostic performance between modalities. An examination-room simulation was performed of the two methods to estimate representative examination times. RESULTS: The mean effective radiation dose for ULD CT with the tin filter was estimated to be 0.29 ± 0.16 mSv compared with 0.16 ± 0.19 mSv for digital plain radiography, with both examinations associated with a very low lifetime attributable risk (< 0.01%). The shunt tip could be more reliably located with ULD CT. ULD CT also allowed assessment of additional findings to explain patient symptoms, such as a cyst at the tip of the shunt catheter and the presence of an obstructing rubber nipple in the duodenum that could not have been observed on a plain radiograph. The examination time with ULD CT of the shunt was estimated to be 20 minutes. The examination time of the shunt with digital plain radiography, including the examination itself time and transfer of the patient between rooms, was estimated to be 60 minutes. CONCLUSIONS: ULD CT using a tin filter allows good visualization of the position or disruption of the shunt catheter that is comparable or superior to plain radiography, at a higher dose, while providing additional findings and reduced patient discomfort.

Technical note: suppression of artifacts arising from simultaneous cone-beam imaging and RF transponder tracking in prostate radiotherapy.

PURPOSE: Artifacts in treatment-room cone-beam reconstructions have been observed at the authors' center when cone-beam acquisition is simultaneous with radio frequency (RF) transponder tracking using the Calypso 4D system (Calypso Medical, Seattle, WA). These artifacts manifest as CT-number modulations and increased CT-noise. The authors present a method for the suppression of the artifacts. METHODS: The authors propose a three-stage postprocessing technique that can be applied to image volumes previously reconstructed by a cone-beam system. The stages are (1) segmentation of voxels into air, soft-tissue, and bone; (2) application of a 2D spatial-filter in the axial plane to the soft-tissue voxels; and (3) normalization to remove streaking along the axial-direction. The algorithm was tested on patient data acquired with Synergy XVI cone-beam CT systems (Elekta, Crawley, United Kingdom). RESULTS: The computational demands of the suggested correction are small, taking less than 15 s per cone-beam reconstruction on a desktop PC. For a moderate loss of spatial-resolution, the artifacts are strongly suppressed and low-contrast visibility is improved. CONCLUSIONS: The correction technique proposed is fast and effective in removing the artifacts caused by simultaneous cone-beam imaging and RF-transponder tracking.

Technical Note: SpekPy v2.0-a software toolkit for modeling x-ray tube spectra.

PURPOSE: SpekPy is a free toolkit for modeling x-ray tube spectra with the Python programming language. In this article, the advances in version 2.0 (v2) of the software are described, including additional target materials and more accurate modeling of the heel effect. Use of the toolkit is also demonstrated. METHODS: The predictions of SpekPy are illustrated in comparison to experimentally determined spectra: three radiation quality reference (RQR) series tungsten spectra and one mammography spectrum with a molybdenum target. The capability of the software to correctly model changes in tube output with tube potential is also assessed, using the example of a GE RevolutionTM CT scanner (GE Healthcare, Waukesha, WI, USA) and specifications in the system's Technical Reference Manual. Furthermore, we note that there are several physics models available in SpekPy. These are compared on and off the central axis, to illustrate the differences. RESULTS: SpekPy agrees closely with the experimental spectra over a wide range of tube potentials, both visually and in terms of first and second half-value layers (HVLs) (within 2% here). The CT scanner spectrum output (normalized to 120 kV tube potential) agreed within 4% over the range of 70 to 140 kV. The default physics model (casim) is adequate in most situations. The advanced option (kqp) should be used if high accuracy is desired for modeling the anode heel effect, as it fully includes the effects of bremsstrahlung anisotropy. CONCLUSIONS: SpekPy v2 can reliably predict on- and off-axis spectra for tungsten and molybdenum targets. SpekPy's open-source MIT license allows users the freedom to incorporate this powerful toolkit into their own projects.

Estimation of delivered dose to lung tumours considering setup uncertainties and breathing motion in a cohort of patients treated with stereotactic body radiation therapy.

INTRODUCTION: Dose-response relationships for local control of lung tumours treated with stereotactic body radiotherapy (SBRT) have proved ambiguous, however, these have been based on the prescribed or planned dose. Delivered dose to the target may be a better predictor for local control. In this study, the probability of the delivered minimum dose to the clinical target volume (CTV) in relation to the prescribed dose was estimated for a cohort of patients, considering geometrical uncertainties. MATERIALS AND METHODS: Delivered doses were retrospectively simulated for 50 patients treated with SBRT for lung tumours, comparing two image-guidance techniques: pre-treatment verification computed tomography (IG1) and online cone-beam computed tomography (IG2). The prescribed dose was typically to the 67% isodose line of the treatment plan. Simulations used in-house software that shifted the static planned dose according to a breathing motion and sampled setup/matching errors. Each treatment was repeatedly simulated, generating a multiplicity of dose-volume histograms (DVH). From these, tumour-specific and population-averaged statistics were derived. RESULTS: For IG1, the probability that the minimum CTV dose (D98%) exceeded 100% of the prescribed dose was 90%. With IG2, this probability increased to 99%. CONCLUSIONS: Doses below the prescribed dose were delivered to a considerably larger part of the population prior to the introduction of online soft-tissue image-guidance. However, there is no clear evidence that this impacts local control, when compared to previous published data.

Performance characterization of a prototype dual-layer cone-beam computed tomography system.

PURPOSE: Conventional cone-beam computed tomography CT (CBCT) provides limited discrimination between low-contrast tissues. Furthermore, it is limited to full-spectrum energy integration. A dual-energy CBCT system could be used to separate photon energy spectra with the potential to increase the visibility of clinically relevant features and acquire additional information relevant in a multitude of clinical imaging applications. In this work, the performance of a novel dual-layer dual-energy CBCT (DL-DE-CBCT) C-arm system is characterized for the first time. METHODS: A prototype dual-layer detector was fitted into a commercial interventional C-arm CBCT system to enable DL-DE-CBCT acquisitions. DL-DE reconstructions were derived from material-decomposed Compton scatter and photoelectric base functions. The modulation transfer function (MTF) of the prototype DL-DE-CBCT was compared to that of a commercial CBCT. Noise and uniformity characteristics were evaluated using a cylindrical water phantom. Effective atomic numbers and electron densities were estimated in clinically relevant tissue substitutes. Iodine quantification was performed (for 0.5-15 mg/ml concentrations) and virtual noncontrast (VNC) images were evaluated. Finally, contrast-to-noise ratios (CNR) and CT number accuracies were estimated. RESULTS: The prototype and commercial CBCT showed similar spatial resolution, with a mean 10% MTF of 5.98 cycles/cm and 6.28 cycles/cm, respectively, using a commercial standard reconstruction. The lowest noise was seen in the 80 keV virtual monoenergetic images (VMI) (7.40 HU) and the most uniform images were seen at VMI 60 keV (4.74 HU) or VMI 80 keV (1.98 HU), depending on the uniformity measure used. For all the tissue substitutes measured, the mean accuracy in effective atomic number was 98.2% (SD 1.2%) and the mean accuracy in electron density was 100.3% (SD 0.9%). Iodine quantification images showed a mean difference of -0.1 (SD 0.5) mg/ml compared to the true iodine concentration for all blood and iodine-containing objects. For VNC images, all blood substitutes containing iodine averaged a CT number of 43.2 HU, whereas a blood-only substitute measured 44.8 HU. All water-containing iodine substitutes measured a mean CT number of 2.6 in the VNC images. A noise-suppressed dataset showed a CNR peak at VMI 40 keV and low at VMI 120 keV. In the same dataset without noise suppression applied, a peak in CNR was obtained at VMI 70 keV and a low at VMI 120 keV. The estimated CT numbers of various clinically relevant objects were generally very close to the calculated CT number. CONCLUSIONS: The performance of a prototype dual-layer dual-energy C-arm CBCT system was characterized. Spatial resolution and noise were comparable with a commercially available C-arm CBCT system, while offering dual-energy capability. Iodine quantifications, effective atomic numbers, and electron densities were in good agreement with expected values, indicating that the system can be used to reliably evaluate the material composition of clinically relevant tissues. The VNC and monoenergetic images indicate a consistent ability to separate clinically relevant tissues. The results presented indicate that the system could find utility in diagnostic, interventional, and radiotherapy planning settings.

Age-specific and gender-specific radiation risks in paediatric angiography and interventional cardiology: conversion coefficients and risk reference values.

OBJECTIVES: To estimate risk for exposure-induced cancer death (REID), organ-specific risks of exposure-induced cancer death (REIDHT) and associated conversion coefficients (CCREID:KAP=REID/kerma-area product (KAP), CCREIDHT:KAP=REIDHT/KAP) in paediatric cardiac catheterizations using data from radiation dose structured reports (RDSR). A novel risk surveillance tool consisting of age-specific and gender-specific risk reference values (RRVs) related to population cancer risk is suggested. METHODS: The PCXMC v.2.0 code is used together with exposure-related information from RDSR from a cohort of 238 children to assess cancer risks and related conversion coefficients. The KAP corresponding to 1 in 1000 of increased REID is used to define age-specific and gender-specific KAP values to monitor risk in such patient cohorts, here denoted as RRVs. RESULTS: The REID estimates ranged from below 1 up to 300 in 100,000, and the RRVs for the different age groups and gender ranged from 0.77 Gycm2 and 2.1 Gycm2 for neonates (female, male) to 11 Gycm2 and 25 Gycm2 for 15-year-olds (female, male). The CCREID:KAP and CCREIDHT:KAP decreased biexponentially with increased age, being notably higher for female patients. CONCLUSIONS: Prominent risk contributing organs were the lungs and the (female) breast. The concept of age-specific and gender-specific RRVs related to population cancer risk is introduced and is intended to be used as a supporting tool for physicians performing such interventions. ADVANCES IN KNOWLEDGE: Age-related and gender-related conversion coefficients for radiation risk, CCREID:KAP and CCREIDHT:KAP, are introduced and a novel risk surveillance concept, the RRV, is suggested for paediatric cardiac catheterizations.

A model for the energy and angular distribution of x rays emitted from an x-ray tube. Part I. Bremsstrahlung production.

PURPOSE: To develop an analytical model for bremsstrahlung production in a thick x-ray target (i.e., the x-ray tube anode) that takes into account the intrinsic bremsstrahlung angular distribution. METHODS: X-ray spectrum models developed from theoretical principles have traditionally treated the angular distribution of the bremsstrahlung production as spherically uniform. This assumption stems from the rationale that electrons promptly attain a diffuse directional distribution in an x-ray target due to multiple scattering, thereby effectively masking the intrinsic bremsstrahlung angular distribution. In this work, a model that explicitly accounts for the angular distribution of the bremsstrahlung production is presented. The model combines Monte Carlo-calculated depth, energy, and angular distributions of electrons penetrating the x-ray target, and incorporates theoretical results for the differential bremsstrahlung cross section. The effects of using different simplified model assumptions for the electron penetration and the intrinsic bremsstrahlung angular distribution are analyzed for tungsten and molybdenum targets in the energy range 20-300 keV. RESULTS: Typical assumptions of previous models are shown to introduce errors in calculated spectra. Particularly, it is shown that predictions of fluence and air kerma free-in-air can be overestimated by 15-30% (2-3% in aluminum half-value layer thickness) for clinically relevant beam qualities. The present model is able to reproduce comprehensive Monte Carlo calculations of the bremsstrahlung production generally to within 1%. CONCLUSIONS: The bremsstrahlung model developed in this work is an improvement over previous models in that the main features of the electron penetration and the resulting bremsstrahlung are considered in detail. The model can be used for more accurate predictions of the energy and angular distribution of x rays emitted from an x-ray tube.

A validation of SpekPy: A software toolkit for modelling X-ray tube spectra.

PURPOSE: To validate the SpekPy software toolkit that has been developed to estimate the spectra emitted from tungsten anode X-ray tubes. The model underlying the toolkit introduces improvements upon a well-known semi-empirical model of X-ray emission. MATERIALS AND METHODS: Using the same theoretical framework as the widely-used SpekCalc software, new electron penetration data was simulated using the Monte Carlo (MC) method, alternative bremsstrahlung cross-sections were applied, L-line characteristic emissions were included, and improvements to numerical methods implemented. The SpekPy toolkit was developed with the Python programming language. The toolkit was validated against other popular X-ray spectrum models (50 to 120 kVp), X-ray spectra estimated with MC (30 to 150 kVp) as well as reference half value layers (HVL) associated with numerous radiation qualities from standard laboratories (20 to 300 kVp). RESULTS: The toolkit can be used to estimate X-ray spectra that agree with other popular X-ray spectrum models for typical configurations in diagnostic radiology as well as with MC spectra over a wider range of conditions. The improvements over SpekCalc are most evident at lower incident electron energies for lightly and moderately filtered radiation qualities. Using the toolkit, estimations of the HVL over a large range of standard radiation qualities closely match reference values. CONCLUSIONS: A toolkit to estimate X-ray spectra has been developed and extensively validated for central-axis spectra. This toolkit can provide those working in Medical Physics and beyond with a powerful and user-friendly way of estimating spectra from X-ray tubes.

A model for the energy and angular distribution of x rays emitted from an x-ray tube. Part II. Validation of x-ray spectra from 20 to 300 kV.

PURPOSE: To present and validate a complete x-ray emission model (bremsstrahlung and characteristic x-ray emission) for the energy range 20-300 kV. METHODS: An analytical x-ray spectrum model that combines the bremsstrahlung emission model developed in Part I with a previously developed characteristic x-ray emission model is validated by comparison with Monte Carlo calculations, published measured spectra, and models developed by other authors. Furthermore, the assumptions and limitations of previous spectrum models are summarized, and their predictions are compared with results obtained by Monte Carlo simulations of x rays emitted from tungsten and molybdenum targets. RESULTS: The model is able to reproduce narrow-beam Monte Carlo calculations to within 0.5% in terms of the first and second aluminum half-value layer thickness (HVL). Compared with measured spectra, the difference in HVL is < 2% for typical diagnostic and therapeutic beam qualities available at primary standard laboratories. Compared with previous spectrum models, the present model performs especially well for low kilovoltage x-ray beams (below 50 kV), and is reliable for a wider range of take-off angles, that is, the angle between the target surface and the direction of emission. The difference in model and Monte Carlo predictions of the energy-fluence weighted air kerma (i.e., the photon energy absorption in air) is < 0.5% using the present model, while previous spectrum models can differ by more than 10%. CONCLUSIONS: The x-ray emission model developed in this work has been validated against Monte Carlo calculations and measured results. The model provides an efficient alternative to comprehensive Monte Carlo simulations and is an improvement over previous models. The model can be used to predict both central- and off-axis spectra, as well as off-axis effects such as the (anode) heel effect.

Experimental assessment of a phosphor model for estimating the relative extrinsic efficiency in radioluminescent detectors.

Optimising phosphor screens in dose detectors or imaging sensor designs is a cumbersome and time- consuming work normally involving specialised measuring equipment and advanced modelling. It is known that crucial optical parameters of the same phosphor may vary within a wide range of values. The aim of this work was to experimentally assess a simple previously published model where the case specific optical parameters (scattering and absorption) are instead represented by a fixed, single parameter, the light extinction factor, ξ. The term extrinsic efficiency, N, of a phosphor is also introduced, differing from the common denotation "absolute efficiency", after noting that unknown factors (such as temperature dependence) can have an influence during efficiency estimations and hence difficult to claim absoluteness. N is expressed as the ratio of light energy emitted per unit area at the phosphor surface to incident x-ray energy fluence. By focusing on ratios and relative changes in this study, readily available instruments in a Medical Physics Department (i.e. a photometer) could be used. The varying relative extrinsic efficiency for an extended range of particle sizes (7.5 and 25 µm) and layer thicknesses (220 to 830 µm) were calculated in the model from the input parameters: the mean particle size of the phosphor, the layer thickness, the light extinction factor and the calculated energy imparted to the layer. In-house manufactured screens (Gd2O2S:Tb) were used for better control of design parameters. The model provided good qualitative agreement to experiment with quantitative deviations in relative extrinsic efficiency within approximately 2%.

Calculating organ and effective doses in paediatric interventional cardiac radiology based on DICOM structured reports - Is detailed examination data critical to dose estimates?

PURPOSE: To estimate effective dose (E), equivalent organ doses (HT) and associated conversion coefficients (CCE:KAP = E/KAP, CCHT:KAP = HT/KAP; KAP = Kerma-area product) in paediatric cardiac interventions, using detailed exposure data from radiation dose structured reports (RDSR). These "RDSR dose estimations" have been compared with estimations performed using the approach currently implemented in the clinic that is based on a simplified assumptions method (SAM). METHODS: The Monte Carlo system PCXMC, incorporated into a previously developed framework, was used to calculate E and HT for 202 children. The calculations were performed with input values from RDSR, and also using simplified assumptions, including fixed nominal values for the focus-skin distance, collimated beam size, irradiation geometry and patient size (age, weight and height). RESULTS: Mean HT to critical organs were: 5-25 mSv (lungs), 5-8 mSv (breasts) and 5-22 mSv (heart), with the lower and upper end of the doses associated with the neonatal and 15 years group, respectively. The associated mean CCHT:KAP for the different age groups were: 9.4-1.6 mSv/Gycm2 (lungs), 8.9-0.54 mSv/Gycm2 (breasts) and 9.3-1.4 mSv/Gycm2 (heart). CONCLUSIONS: The extension of the concept of a conversion coefficient for HT is introduced and CCHT:KAP values for paediatric cardiac interventions divided in age groups are presented. This method of linking the KAP to HT is intended for use in epidemiological/cohort studies or in clinics that do not have access to RDSR. Further, the population-averaged conversion coefficients for the critical organs estimated from RDSR, displayed no statistically significant difference compared with the SAM approach.

The synthetic localizer radiograph - A new CT scan planning method.

OBJECTIVE: To investigate if the conventional localizer radiograph (LR) can be replaced by a synthetic LR (SLR), generated from a low-dose spiral CT scan, for CT scan planning with minimal changes to current clinical workflows. METHODS: A dosimetric comparison of SLRs and LRs was made using Monte Carlo methods. Water equivalent diameters (WEDs) of a centered and mis-centered phantom were estimated from low-dose spiral CT scans and LRs acquired at different angles. Body sizes, in the form of two lengths and two diameters obtained from SLRs and LRs, were compared for 10 patients (4 men and 6 women with a mean age of 74.8 and 76.2 years respectively) undergoing CT of thorax and abdomen. The image quality of SLRs for CT scan planning relative to LRs was rated using a 5-grade scale by four radiologists and two CT radiographers. RESULTS: An SLR can be obtained at a comparable effective dose to that of traditionally acquired LRs: 0.14 mSv. WEDs from LRs were more affected by mis-centering than WEDs calculated from low-dose spiral scans. One significant discrepancy of estimated body sizes was observed, the broadest part of the patient that on lateral localizers showed a mean deviation of 17.7 mm (range: 7.3-28.7 mm, p < 0.001). The anteroposterior/posteroanterior SLR image quality was assessed as better compared to an LR while the same could not be shown for lateral localizers. CONCLUSIONS: SLRs based on low-dose spiral scans can replace LRs for CT planning.

Comparison of acquisition protocols for ventilation/perfusion SPECT-a Monte Carlo study.

One of the most commonly used imaging techniques for diagnosing pulmonary embolism (PE) is ventilation/perfusion (V/P) scintigraphy. The aim of this study was to evaluate the performance of the currently used imaging protocols for V/P single photon emission computed tomography (V/P SPECT) at two nuclear medicine department sites and to investigate the effect of altering important protocol parameters. The Monte Carlo technique was used to simulate 4D digital phantoms with perfusion defects. Six imaging protocols were included in the study and a total of 72 digital patients were simulated. Six dually trained radiologists/nuclear medicine physicians reviewed the images and reported all perfusion mismatch findings. The radiologists also visually graded the image quality. No statistically significant differences in diagnostic performance were found between the studied protocols, but visual grading analysis pointed out one protocol as significantly superior to four of the other protocols. Considering the study results, we have decided to harmonize our clinical protocols for imaging patients with suspected PE. The administered Technegas and macro aggregated albumin activities have been altered, a low energy all purpose collimator is used instead of a low energy high resolution collimator and the acquisition times have been lowered.