Averaging of absorbed doses: How matter matters.

BACKGROUND: Dosimetry in radionuclide therapy often requires the calculation of average absorbed doses within and between spatial regions, for example, for voxel-based dosimetry methods, for paired organs, or across multiple tumors. Formation of such averages can be made in different ways, starting from different definitions. PURPOSE: The aim of this study is to formally specify different averaging strategies for absorbed doses, and to compare their results when applied to absorbed dose distributions that are non-uniform within and between regions. METHODS: For averaging within regions, two definitions of the average absorbed dose are considered: the simple average over the region (the region average) and the average when weighting by the mass density (density-weighted region average). The latter is shown to follow from the definition of mean absorbed dose according to the ICRU, and to be consistent with the MIRD formalism. For averaging between different spatial regions, three definitions follow: the volume-weighted, the mass-weighted, and the unweighted average. With respect to characterizing non-uniformity, the different average definitions lead to the use of dose-volume histograms (DVHs) (region average), dose-mass histograms (DMHs) (density-weighted region average), and unweighted histograms (unweighted average). Average absorbed doses are calculated for three worked examples, starting from the different definitions. The first, schematic, example concerns the calculation of the average absorbed dose between two regions with different volumes or mass densities. The second, stylized, example concerns voxel-based dosimetry, for which the average absorbed-dose rate within a region is calculated. The geometries studied include three 177 Lu-filled voxelized spheres, where the sphere masses are held constant while the material compositions, densities, and volumes are varied. For comparison, the mean absorbed-dose rates obtained using unit-density sphere S-values are also included. The third example concerns SPECT/CT-based tumor dosimetry for five patients undergoing therapy with 177 Lu-PSMA and six patients undergoing therapy with 177 Lu-DOTA-TATE, for which the average absorbed-dose rates across multiple tumors are calculated. For the second and third examples, analyses also include representations by histograms. RESULTS: Example 1 shows that the average absorbed doses, calculated using different definitions, can differ considerably if the masses and absorbed doses for two regions are markedly different. From example 2 it is seen that the density-weighted region average is stable under different activity and density distributions and is also in line with results using S-values. In contrast, the region average varies as function of the activity distribution. In example 3, the absorbed dose rates for individual tumors differ by (1.1 ± 4.3)% and (-0.1 ± 0.4)% with maximum deviations of +34.4% and -1.4% for 177 Lu-PSMA and 177 Lu-DOTA-TATE, respectively, when calculated as region averages or density-weighted region averages, with largest deviations obtained when the density is non-uniform. The average absorbed doses calculated across all tumors are similar when comparing mass-weighted and volume-weighted averages but these differ substantially from unweighted averages. CONCLUSION: Different strategies for averaging of absorbed doses within and between regions can lead to substantially different absorbed-dose estimates. At reporting of radionuclide therapy dosimetry, it is important to specify the averaging strategy applied.

ACCURACY OF CT NUMBERS OBTAINED BY DIRA AND MONOENERGETIC PLUS ALGORITHMS IN DUAL-ENERGY COMPUTED TOMOGRAPHY.

Dual-energy computed tomography (CT) can be used in radiotherapy treatment planning for the calculation of absorbed dose distributions. The aim of this work is to evaluate whether there is room for improvement in the accuracy of the Monoenergetic Plus algorithm by Siemens Healthineers. A Siemens SOMATOM Force scanner was used to scan a cylindrical polymethyl methacrylate phantom with four rod-inserts made of different materials. Images were reconstructed using ADMIRE and processed with Monoenergetic Plus. The resulting CT numbers were compared with tabulated values and values simulated by the proof-of-a-concept algorithm DIRA developed by the authors. Both the Monoenergetic Plus and DIRA algorithms performed well; the accuracy of attenuation coefficients was better than about ±1% at the energy of 70 keV. Compared with DIRA, the worse performance of Monoenergetic Plus was caused by its (i) two-material decomposition to iodine and water and (ii) imperfect suppression of the beam hardening artifact in ADMIRE.

OPTIMAL SELECTION OF BASE MATERIALS FOR ACCURATE DUAL-ENERGY COMPUTED TOMOGRAPHY: COMPARISON BETWEEN THE ALVAREZ-MACOVSKI METHOD AND DIRA.

The choice of the material base to which the material decomposition is performed in dual-energy computed tomography may affect the quality of reconstructed images. The aim of this work is to investigate how the commonly used bases (water, bone), (water, iodine) and (photoelectric effect, Compton scattering) affect the reconstructed linear attenuation coefficient in the case of the Alvarez-Macovski method. The performance of this method is also compared with the performance of the Dual-energy Iterative Reconstruction Algorithm (DIRA). In both cases, the study is performed using simulations. The results show that the Alvarez-Macovski method produced artefacts when iodine was present in the phantom together with human tissues since this method can only work with one doublet. It was shown that these artefacts could be avoided with DIRA using the (water, bone) doublet for tissues and the (water, iodine) doublet for the iodine solution.

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%.

Investigation of a synthetic diamond detector response in kilovoltage photon beams.

PURPOSE: An important characteristic of radiation dosimetry detectors is their energy response which consists of absorbed-dose and intrinsic energy responses. The former can be characterized using Monte Carlo (MC) simulations, whereas the latter (i.e., detector signal per absorbed dose to detector) is extracted from experimental data. Such a characterization is especially relevant when detectors are used in nonrelative measurements at a beam quality that differs from the calibration beam quality. Having in mind the possible application of synthetic diamond detectors (microDiamond PTW 60019, Freiburg, Germany) for nonrelative dosimetry of low-energy brachytherapy (BT) beams, we determined their intrinsic and absorbed-dose energy responses in 25-250 kV beams relative to a 60 Co beam, which is usually the reference beam quality for detector calibration in radiotherapy. MATERIAL AND METHODS: Three microDiamond detectors and, for comparison, two silicon diodes (PTW 60017) were calibrated in terms of air-kerma free in air in six x-ray beam qualities (from 25 to 250 kV) and in terms of absorbed dose to water in a 60 Co beam at the national metrology laboratory in Sweden. The PENELOPE/penEasy MC radiation transport code was used to calculate the absorbed-dose energy response of the detectors (modeled based on blueprints) relative to air and water depending on calibration conditions. The MC results were used to extract the relative intrinsic energy response of the detectors from the overall energy response. Measurements using an independent setup with a single ophthalmic BEBIG I25.S16 125 I BT seed (effective photon energy of 28 keV) were used as a qualitative check of the extracted intrinsic energy response correction factors. Additionally, the impact of the thickness of the active volume as well as the presence of extra-cameral components on the absorbed-dose energy response of a microDiamond detector was studied using MC simulations. RESULTS: The relative intrinsic energy response of the microDiamond detectors was higher by a factor of 2 in 25 and 50 kV beams compared to the 60 Co beam. The variation in the relative intrinsic energy response of silicon diodes was within 10% over the investigated photon energy range. The use of relative intrinsic energy response correction factors improved the agreement among the absorbed dose to water values determined using microDiamond detectors and silicon diodes, as well as with the TG-43 formalism-based calculations for the 125 I seed. MC study of microDiamond detector design features provided a possible explanation for inter-detector response variation at low-energy photon beams by differences in the effective thickness of the active volume. CONCLUSIONS: MicroDiamond detectors had a non-negligible variation in the relative intrinsic energy response (factor of 2) which was comparable to that in the absorbed-dose energy response relative to water at low-energy photon beams. Silicon diodes, in contrast, had an absorbed-dose energy dependence on photon energy that varied by a factor of 6, whereas the intrinsic energy dependence on beam quality was within 10%. It is important to decouple these two responses for a full characterization of detector energy response especially when the user and reference beam qualities differ significantly, and MC alone is not enough.

Suitability of microDiamond detectors for the determination of absorbed dose to water around high-dose-rate 192 Ir brachytherapy sources.

PURPOSE: Experimental dosimetry of high-dose-rate (HDR) 192 Ir brachytherapy (BT) sources is complicated due to high dose and dose-rate gradients, and softening of photon energy spectrum with depth. A single crystal synthetic diamond detector microDiamond (PTW 60019, Freiburg, Germany) has a small active volume, high sensitivity, direct readout, and nearly water-equivalent active volume. The purpose of this study was to evaluate the suitability of microDiamond detectors for the determination of absorbed dose to water around HDR 192 Ir BT sources. Three microDiamond detectors were used, allowing for the comparison of their properties. METHODS: In-phantom measurements were performed using microSelectron and VariSource iX HDR 192 Ir BT treatment units. Their treatment planning systems (TPSs), Oncentra (v. 4.3) and BrachyVision (v. 13.6), respectively, were used to create irradiation plans for a cubic PMMA phantom with the microDiamond positioned at one of three source-to-detector distances (SDDs) (1.5, 2.5, and 5.5 cm) at a time. The source was stepped in increments of 0.5 cm over a total length of 6 cm to yield absorbed dose of 2 Gy at the nominal reference-point of the detector. Detectors were calibrated in 60 Co beam in terms of absorbed dose to water, and Monte Carlo (MC) calculated beam quality correction factors were applied to account for absorbed-dose energy dependence. Phantom correction factors were applied to account for differences in dimensions between the measurement phantom and a water phantom used for absorbed dose calculations made with a TPS. The same measurements were made with all three of the detectors. Additionally, dose-rate dependence and stability of the detectors were evaluated in 60 Co beam. RESULTS: The percentage differences between experimentally determined and TPS-calculated absorbed doses to water were from -1.3% to +2.9%. The values agreed to within experimental uncertainties, which were from 1.9% to 4.3% (k = 2) depending on the detector, SDD and treatment delivery unit. No dose-rate or intrinsic energy dependence corrections were applied. All microDiamonds were comparable in terms of preirradiation dose, stability of the readings and energy response, and showed a good agreement. CONCLUSIONS: The results indicate that the microDiamond is potentially suitable for the determination of absorbed dose to water around HDR 192 Ir BT sources and may be used for independent verification of TPS's calculations, as well as for QA measurements of HDR 192 Ir BT treatment delivery units at clinical sites.

A model-based iterative reconstruction algorithm DIRA using patient-specific tissue classification via DECT for improved quantitative CT in dose planning.

PURPOSE: To develop and evaluate-in a proof-of-concept configuration-a novel iterative reconstruction algorithm (DIRA) for quantitative determination of elemental composition of patient tissues for application to brachytherapy with low energy (< 50 keV) photons and proton therapy. METHODS: DIRA was designed as a model-based iterative reconstruction algorithm, which uses filtered backprojection, automatic segmentation and multimaterial tissue decomposition. The evaluation was done for a phantom derived from the voxelized ICRP 110 male phantom. Soft tissues were decomposed to the lipid, protein and water triplet, bones were decomposed to the compact bone and bone marrow doublet. Projections were derived using the Drasim simulation code for an axial scanning configuration resembling a typical DECT (dual-energy CT) scanner with 80 kV and Sn140 kV x-ray spectra. The iterative loop produced mono-energetic images at 50 and 88 keV without beam hardening artifacts. Different noise levels were considered: no noise, a typical noise level in diagnostic imaging and reduced noise level corresponding to tenfold higher doses. An uncertainty analysis of the results was performed using type A and B evaluations. The two approaches were compared. RESULTS: Linear attenuation coefficients averaged over a region were obtained with relative errors less than 0.5% for all evaluated regions. Errors in average mass fractions of the three-material decomposition were less than 0.04 for no noise and reduced noise levels and less than 0.11 for the typical noise level. Mass fractions of individual pixels were strongly affected by noise, which slightly increased after the first iteration but subsequently stabilized. Estimates of uncertainties in mass fractions provided by the type B evaluation differed from the type A estimates by less than 1.5% for most cases. The algorithm was fast, the results converged after 5 iterations. The algorithmic complexity of forward polyenergetic projection calculation was much reduced by using material doublets and triplets. CONCLUSIONS: The simulations indicated that DIRA is capable of determining elemental composition of tissues, which are needed in brachytherapy with low energy (< 50 keV) photons and proton therapy. The algorithm provided quantitative monoenergetic images with beam hardening artifacts removed. Its convergence was fast, image sharpness expressed via the modulation transfer function was maintained, and image noise did not increase with the number of iterations.

Optimisation of quantitative lung SPECT applied to mild COPD: a software phantom simulation study.

BACKGROUND: The amount of inhomogeneities in a (99m)Tc Technegas single-photon emission computed tomography (SPECT) lung image, caused by reduced ventilation in lung regions affected by chronic obstructive pulmonary disease (COPD), is correlated to disease advancement. A quantitative analysis method, the CVT method, measuring these inhomogeneities was proposed in earlier work. To detect mild COPD, which is a difficult task, optimised parameter values are needed. METHODS: In this work, the CVT method was optimised with respect to the parameter values of acquisition, reconstruction and analysis. The ordered subset expectation maximisation (OSEM) algorithm was used for reconstructing the lung SPECT images. As a first step towards clinical application of the CVT method in detecting mild COPD, this study was based on simulated SPECT images of an advanced anthropomorphic lung software phantom including respiratory and cardiac motion, where the mild COPD lung had an overall ventilation reduction of 5%. RESULTS: The best separation between healthy and mild COPD lung images as determined using the CVT measure of ventilation inhomogeneity and 125 MBq (99m)Tc was obtained using a low-energy high-resolution collimator (LEHR) and a power 6 Butterworth post-filter with a cutoff frequency of 0.6 to 0.7 cm(-1). Sixty-four reconstruction updates and a small kernel size should be used when the whole lung is analysed, and for the reduced lung a greater number of updates and a larger kernel size are needed. CONCLUSIONS: A LEHR collimator and 125 (99m)Tc MBq together with an optimal combination of cutoff frequency, number of updates and kernel size, gave the best result. Suboptimal selections of either cutoff frequency, number of updates and kernel size will reduce the imaging system's ability to detect mild COPD in the lung phantom.

A system for remote dosimetry audit of 3D-CRT, IMRT and VMAT based on lithium formate dosimetry.

The aim of this work was to develop and test a remote end-to-end audit system using lithium formate EPR dosimeters. Four clinics were included in a pilot study, absorbed doses determined in the PTV agreed with TPS calculated doses within ±5% for 3D-CRT and ±7% (k=1) for IMRT/VMAT dose plans.

Investigation of signal fading in lithium formate EPR dosimeters using a new sensitive method.

The aim of this study was to investigate signal fading in lithium formate electron paramagnetic resonance (EPR) dosimeters used for clinical applications in radiotherapy. A new experimental method for determination of signal fading, designed to resolve small changes in signal from slowly decaying unstable radicals, was used. Possible signal fading in lithium formate due to different storage temperatures was also tested. Air humidity was kept at a constant level of 33% throughout the experiments. The conclusion drawn from the investigations was that the EPR signal from lithium formate is stable during at least 1 month after irradiation and is not sensitive to variations in storage temperature <40 °C when kept at a relative air humidity of 33%. This makes lithium formate a suitable dosimeter for transfer dosimetry in clinical audits.

Determination of absorbed dose to water around a clinical HDR (192)Ir source using LiF:Mg,Ti TLDs demonstrates an LET dependence of detector response.

PURPOSE: Experimental radiation dosimetry with thermoluminescent dosimeters (TLDs), calibrated in a (60)Co or megavoltage (MV) photon beam, is recommended by AAPM TG-43U1for verification of Monte Carlo calculated absorbed doses around brachytherapy sources. However, it has been shown by Carlsson Tedgren et al. [Med. Phys. 38, 5539-5550 (2011)] that for TLDs of LiF:Mg,Ti, detector response was 4% higher in a (137)Cs beam than in a (60)Co one. The aim of this work was to investigate if similar over-response exists when measuring absorbed dose to water around (192)Ir sources, using LiF:Mg,Ti dosimeters calibrated in a 6 MV photon beam. METHODS: LiF dosimeters were calibrated to measure absorbed dose to water in a 6 MV photon beam and used to measure absorbed dose to water at distances of 3, 5, and 7 cm from a clinical high dose rate (HDR) (192)Ir source in a polymethylmethacrylate (PMMA) phantom. Measured values were compared to values of absorbed dose to water calculated using a treatment planning system (TPS) including corrections for the difference in energy absorption properties between calibration quality and the quality in the users' (192)Ir beam and for the use of a PMMA phantom instead of the water phantom underlying dose calculations in the TPS. RESULTS: Measured absorbed doses to water around the (192)Ir source were overestimated by 5% compared to those calculated by the TPS. Corresponding absorbed doses to water measured in a previous work with lithium formate electron paramagnetic resonance (EPR) dosimeters by Antonovic et al. [Med. Phys. 36, 2236-2247 (2009)], using the same irradiation setup and calibration procedure as in this work, were 2% lower than those calculated by the TPS. The results obtained in the measurements in this work and those obtained using the EPR lithium formate dosimeters were, within the expanded (k = 2) uncertainty, in agreement with the values derived by the TPS. The discrepancy between the results using LiF:Mg,Ti TLDs and the EPR lithium formate dosimeters was, however, statistically significant and in agreement with the difference in relative detector responses found for the two detector systems by Carlsson Tedgren et al. [Med. Phys. 38, 5539-5550 (2011)] and by Adolfsson et al. [Med. Phys. 37, 4946-4959 (2010)]. CONCLUSIONS: When calibrated in (60)Co or MV photon beams, correction for the linear energy transfer (LET) dependence of LiF:Mg,Ti detector response will be needed as to measure absorbed doses to water in a (192)Ir beam with highest accuracy. Such corrections will depend on the manufacturing process (MTS-N Poland or Harshaw TLD-100) and details of the annealing and read-out schemes used.

Estimating statistical uncertainty of Monte Carlo efficiency-gain in the context of a correlated sampling Monte Carlo code for brachytherapy treatment planning with non-normal dose distribution.

Correlated sampling Monte Carlo methods can shorten computing times in brachytherapy treatment planning. Monte Carlo efficiency is typically estimated via efficiency gain, defined as the reduction in computing time by correlated sampling relative to conventional Monte Carlo methods when equal statistical uncertainties have been achieved. The determination of the efficiency gain uncertainty arising from random effects, however, is not a straightforward task specially when the error distribution is non-normal. The purpose of this study is to evaluate the applicability of the F distribution and standardized uncertainty propagation methods (widely used in metrology to estimate uncertainty of physical measurements) for predicting confidence intervals about efficiency gain estimates derived from single Monte Carlo runs using fixed-collision correlated sampling in a simplified brachytherapy geometry. A bootstrap based algorithm was used to simulate the probability distribution of the efficiency gain estimates and the shortest 95% confidence interval was estimated from this distribution. It was found that the corresponding relative uncertainty was as large as 37% for this particular problem. The uncertainty propagation framework predicted confidence intervals reasonably well; however its main disadvantage was that uncertainties of input quantities had to be calculated in a separate run via a Monte Carlo method. The F distribution noticeably underestimated the confidence interval. These discrepancies were influenced by several photons with large statistical weights which made extremely large contributions to the scored absorbed dose difference. The mechanism of acquiring high statistical weights in the fixed-collision correlated sampling method was explained and a mitigation strategy was proposed.

Comparison of clinical and physical measures of image quality in chest and pelvis computed radiography at different tube voltages.

The aim of this work was to study the dependence of image quality in digital chest and pelvis radiography on tube voltage, and to explore correlations between clinical and physical measures of image quality. The effect on image quality of tube voltage in these two examinations was assessed using two methods. The first method relies on radiologists' observations of images of an anthropomorphic phantom, and the second method was based on computer modeling of the imaging system using an anthropomorphic voxel phantom. The tube voltage was varied within a broad range (50-150 kV), including those values typically used with screen-film radiography. The tube charge was altered so that the same effective dose was achieved for each projection. Two x-ray units were employed using a computed radiography (CR) image detector with standard tube filtration and antiscatter device. Clinical image quality was assessed by a group of radiologists using a visual grading analysis (VGA) technique based on the revised CEC image criteria. Physical image quality was derived from a Monte Carlo computer model in terms of the signal-to-noise ratio, SNR, of anatomical structures corresponding to the image criteria. Both the VGAS (visual grading analysis score) and SNR decrease with increasing tube voltage in both chest PA and pelvis AP examinations, indicating superior performance if lower tube voltages are employed. Hence, a positive correlation between clinical and physical measures of image quality was found. The pros and cons of using lower tube voltages with CR digital radiography than typically used in analog screen-film radiography are discussed, as well as the relevance of using VGAS and quantum-noise SNR as measures of image quality in pelvis and chest radiography.

Towards optimization in digital chest radiography using Monte Carlo modelling.

A Monte Carlo based computer model of the x-ray imaging system was used to investigate how various image quality parameters of interest in chest PA radiography and the effective dose E vary with tube voltage (90-150 kV), additional copper filtration (0-0.5 mm), anti-scatter method (grid ratios 8-16 and air gap lengths 20-40 cm) and patient thickness (20-28 cm) in a computed radiography (CR) system. Calculated quantities were normalized to a fixed value of air kerma (5.0 microGy) at the automatic exposure control chambers. Soft-tissue nodules were positioned at different locations in the anatomy and calcifications in the apical region. The signal-to-noise ratio, SNR, of the nodules and the nodule contrast relative to the contrast of bone (C/C(B)) as well as relative to the dynamic range in the image (C(rel)) were used as image quality measures. In all anatomical regions, except in the densest regions in the thickest patients, the air gap technique provides higher SNR and contrast ratios than the grid technique and at a lower effective dose E. Choice of tube voltage depends on whether quantum noise (SNR) or the contrast ratios are most relevant for the diagnostic task. SNR increases with decreasing tube voltage while C/C(B) increases with increasing tube voltage.

Evaluation of image quality of lumbar spine images: a comparison between FFE and VGA.

PURPOSE: The aim of the present study is to compare two different methods for evaluation of the quality of clinical X-ray images. METHODS: Based on fifteen lumbar spine radiographs, two new sets of images were created. A hybrid image set was created by adding two distributions of artificial lesions to each original image. The image quality parameters spatial resolution and noise were manipulated and a total of 210 hybrid images were created. A set of 105 disease-free images was created by applying the same combinations of spatial resolution and noise to the original images. The hybrid images were evaluated with the free-response forced error experiment (FFE) and the normal images with visual grading analysis (VGA) by nine experienced radiologists. RESULTS: In the VGA study, images with low noise were preferred over images with higher noise levels. The alteration of the MTF had a limited influence on the VGA score. For the FFE study, the visibility of the lesions was independent of the sharpness and the noise level. No correlation was found between the two image quality measures. CONCLUSIONS: FFE is a precise method for evaluation of image quality, but the results are only valid for the type of lesion used in the study, whereas VGA is a more general method for clinical image quality assessment. The results of the FFE study indicate that there might be a potential to lower the dose levels in lumbar spine radiography without losing important diagnostic information.

Inter-observer variation in masked and unmasked images for quality evaluation of clinical radiographs.

PURPOSE: To investigate the influence of masking on the inter-observer variation in image quality evaluation of clinical radiographs of chest and lumbar spine. BACKGROUND: Inter-observer variation is a big problem in image quality evaluation since this variation is often much bigger than the variation in image quality between, for example, two radiographic systems. In this study, we have evaluated the effect of masking on the inter-observer variation. The idea of the masking was to force every observer to view exactly the same part of the image and to avoid the effect of the overall 'first impression' of the image. A discussion with a group of European expert radiologists before the study indicated that masking might be a good way to reduce the inter-observer variation. METHODS: Five chest and five lumbar spine radiographs were collected together with detailed information regarding exposure conditions. The radiographs were digitised with a high-performance scanner and five different manipulations were performed, simulating five different exposure conditions. The contrast, noise and spatial resolution were manipulated by this method. The images were printed onto the film and the individual masks were produced for each film, showing only the parts of the images that were necessary for the image quality evaluation. The quality of the images was evaluated on ordinary viewing boxes by a large group of experienced radiologists. The images were examined with and without the masks with a set of image criteria (if fulfilled, 1 point; and not fulfilled, 0 point), and the mean score was calculated for each simulated exposure condition. RESULTS: The results of this study indicate that-contrary to what was supposed-the inter-observer variation increased when the images were masked. In some cases, especially for chest, this increase was statistically significant. CONCLUSIONS: Based on the results of this study, image masking in studies of fulfilment of image criteria cannot be recommended.

Measurement of radiation dose in dental radiology.

Patient dose audit is an important tool for quality control and it is important to have a well-defined and easy to use method for dose measurements. In dental radiology, the most commonly used dose parameters for the setting of diagnostic reference levels (DRLs) are the entrance surface air kerma (ESAK) for intraoral examinations and dose width product (DWP) for panoramic examinations. DWP is the air kerma at the front side of the secondary collimator integrated over the collimator width and an exposure cycle. ESAK or DWP is usually measured in the absence of the patient but with the same settings of tube voltage (kV), tube current (mA) and exposure time as with the patient present. Neither of these methods is easy to use, and, in addition, DWP is not a risk related quantity. A better method of monitoring patient dose would be to use a dose area product (DAP) meter for all types of dental examinations. In this study, measurements with a DAP meter are reported for intraoral and panoramic examinations. The DWP is also measured with a pencil ionisation chamber and the product of DWP and the height H (DWP x H) of the secondary collimator (measured using film) was compared to DAP. The results show that it is feasible to measure DAP using a DAP meter for both intraoral and panoramic examinations. The DAP is therefore recommended for the setting of DRLs.

The influence of patient thickness and imaging system on patient dose and physical image quality in digital chest imaging.

The aim of this work was to study the influence of patient thickness, tube voltage and image detector on patient dose, contrast and ideal observer signal-to-noise ratio (SNR(I)), for pathological details positioned at different regions in the image in posterior-anterior (PA) chest radiology. A Monte Carlo computational model was used to compute measures of physical image quality (contrast, SNR(I)) and patient effective dose, E. Two metastasis-like details positioned in the central right lung and right lung near the spine, respectively, were studied. The tube voltage was varied between 100 and 150 kV and the patient thickness between 20 and 28 cm. Both, a computed radiography (CR) system and a direct radiography (DR) system, were investigated. The DR system provides both lower doses and better image quality compared with the CR system. The SNR(I)2/E is approximately 2.9 times higher for the DR system compared with the CR system.

Effect of scatter on reconstructed image quality in cone beam computed tomography: evaluation of a scatter-reduction optimisation function.

The effect of scatter on reconstructed image quality in cone beam computed tomography was investigated and a function which can be used in scatter-reduction optimisation tasks was tested. Projections were calculated using the Monte Carlo method in an axially symmetric cone beam geometry consisting of a point source, water phantom and a single row of detector elements. Image reconstruction was performed using the filtered backprojection method. Image quality was assessed by the L2-norm-based difference relative to a reference image derived from (1) weighted linear attenuation coefficients and (2) projections by primary photons. It was found that the former function was strongly affected by the beam hardening artefact and did not properly reflect the amount of scatter but the latter function increased with increasing beam width, was higher for the larger phantom and exhibited properties which made it a good candidate for scatter-reduction optimisation tasks using polyenergetic beams.

Distributions of scatter-to-primary and signal-to-noise ratios per pixel in digital chest imaging.

The aim of this work was to calculate distributions of scatter-to-primary ratios (epsilon(s)/epsilon(p)) and signal-to-noise ratios per pixel (SNRp) in chest images. Such distributions may provide useful information on how physical image quality (contrast, SNR) is distributed over the posterior/anterior (PA) chest image. A Monte Carlo computer program was used for the calculations, including a model of both the patient (voxel phantom) and the imaging system (X-ray tube, anti-scatter grid and image detector). The calculations were performed for three PA thicknesses 20, 24 and 28 cm. For a 24 cm patient, the epsilon(s)/epsilon(p) varies between 0.5 in the lung to 2.5 behind the spine and heart. The corresponding variation of the SNRp is a factor of 3, with the highest values in the lung. Increasing the patient thickness from 20 to 28 cm increases the epsilon(s)/epsilon(p) by a factor of 2.2 behind the spine and heart.