IMAGE QUALITY AND POTENTIAL DOSE REDUCTION USING ADVANCED MODELED ITERATIVE RECONSTRUCTION (ADMIRE) IN ABDOMINAL CT - A REVIEW.

Traditional filtered back projection (FBP) reconstruction methods have served the computed tomography (CT) community well for over 40 years. With the increased use of CT during the last decades, efforts to minimise patient exposure, while maintaining sufficient or improved image quality, have led to the development of model-based iterative reconstruction (MBIR) algorithms from several vendors. The usefulness of the advanced modeled iterative reconstruction (ADMIRE) (Siemens Healthineers) MBIR in abdominal CT is reviewed and its noise suppression and/or dose reduction possibilities explored. Quantitative and qualitative methods with phantom and human subjects were used. Assessment of the quality of phantom images will not always correlate positively with those of patient images, particularly at the higher strength of the ADMIRE algorithm. With few exceptions, ADMIRE Strength 3 typically allows for substantial noise reduction compared to FBP and hence to significant (≈30%) patient dose reductions. The size of the dose reductions depends on the diagnostic task.

Establishing the European diagnostic reference levels for interventional cardiology.

Interventional cardiac procedures may be associated with high patient doses and therefore require special attention to protect the patients from radiation injuries such as skin erythema, cardiovascular tissue reactions or radiation-induced cancer. In this study, patient exposure data is collected from 13 countries (37 clinics and nearly 50 interventional rooms) and for 10 different procedures. Dose data was collected from a total of 14,922 interventional cardiology procedures. Based on these data European diagnostic reference levels (DRL) for air kerma-area product are suggested for coronary angiography (CA, DRL = 35 Gy cm2), percutaneous coronary intervention (PCI, 85 Gy cm2), transcatheter aortic valve implantation (TAVI, 130 Gy cm2), electrophysiological procedures (12 Gy cm2) and pacemaker implantations. Pacemaker implantations were further divided into single-chamber (2.5 Gy cm2) and dual chamber (3.5 Gy cm2) procedures and implantations of cardiac resynchronization therapy pacemaker (18 Gy cm2). Results show that relatively new techniques such as TAVI and treatment of chronic total occlusion (CTO) often produce relatively high doses, and thus emphasises the need for use of an optimization tool such as DRL to assist in reducing patient exposure. The generic DRL presented here facilitate comparison of patient exposure in interventional cardiology.

ACCURATE KAP METER CALIBRATION AS A PREREQUISITE FOR OPTIMISATION IN PROJECTION RADIOGRAPHY.

Modern X-ray units register the air kerma-area product, PKA, with a built-in KAP meter. Some KAP meters show an energy-dependent bias comparable with the maximum uncertainty articulated by the IEC (25 %), adversely affecting dose-optimisation processes. To correct for the bias, a reference KAP meter calibrated at a standards laboratory and two calibration methods described here can be used to achieve an uncertainty of <7 % as recommended by IAEA. A computational model of the reference KAP meter is used to calculate beam quality correction factors for transfer of the calibration coefficient at the standards laboratory, Q0, to any beam quality, Q, in the clinic. Alternatively, beam quality corrections are measured with an energy-independent dosemeter via a reference beam quality in the clinic, Q1, to beam quality, Q Biases up to 35 % of built-in KAP meter readings were noted. Energy-dependent calibration factors are needed for unbiased PKA Accurate KAP meter calibration as a prerequisite for optimisation in projection radiography.

VERIFICATION OF INDICATED SKIN ENTRANCE AIR KERMA FOR CARDIAC X-RAY-GUIDED INTERVENTION USING GAFCHROMIC FILM.

The aim of this work was to verify the indicated maximum entrance surface air kerma (ESAK) using a GE Innova IGS 520 imaging system during cardiac interventional procedures. Gafchromic XR RV3 films were used for the patient measurements to monitor the maximum ESAK. The films were scanned and calibrated to measure maximum ESAK. Thermoluminescent dosemeters were used to measure the backscatter factor from an anthropomorphic thorax phantom. The measured backscatter factor, 1.53, was in good agreement with Monte Carlo simulations but higher than the one used by the imaging system, 1.20. The median of the ratio between indicated maximum ESAK and measured maximum ESAK was 0.68. In this work, the indicated maximum ESAK by the imaging system's dose map model underestimates the measured maximum ESAK by 32 %. The threshold ESAK for follow-up procedures for patient with skin dose in excess of 2 Gy will be reduced to 1.4 Gy.

AUTOMATIC SEGMENTATION OF PELVIS FOR BRACHYTHERAPY OF PROSTATE.

Advanced model-based iterative reconstruction algorithms in quantitative computed tomography (CT) perform automatic segmentation of tissues to estimate material properties of the imaged object. Compared with conventional methods, these algorithms may improve quality of reconstructed images and accuracy of radiation treatment planning. Automatic segmentation of tissues is, however, a difficult task. The aim of this work was to develop and evaluate an algorithm that automatically segments tissues in CT images of the male pelvis. The newly developed algorithm (MK2014) combines histogram matching, thresholding, region growing, deformable model and atlas-based registration techniques for the segmentation of bones, adipose tissue, prostate and muscles in CT images. Visual inspection of segmented images showed that the algorithm performed well for the five analysed images. The tissues were identified and outlined with accuracy sufficient for the dual-energy iterative reconstruction algorithm whose aim is to improve the accuracy of radiation treatment planning in brachytherapy of the prostate.

PARALLELISATION OF THE MODEL-BASED ITERATIVE RECONSTRUCTION ALGORITHM DIRA.

New paradigms for parallel programming have been devised to simplify software development on multi-core processors and many-core graphical processing units (GPU). Despite their obvious benefits, the parallelisation of existing computer programs is not an easy task. In this work, the use of the Open Multiprocessing (OpenMP) and Open Computing Language (OpenCL) frameworks is considered for the parallelisation of the model-based iterative reconstruction algorithm DIRA with the aim to significantly shorten the code's execution time. Selected routines were parallelised using OpenMP and OpenCL libraries; some routines were converted from MATLAB to C and optimised. Parallelisation of the code with the OpenMP was easy and resulted in an overall speedup of 15 on a 16-core computer. Parallelisation with OpenCL was more difficult owing to differences between the central processing unit and GPU architectures. The resulting speedup was substantially lower than the theoretical peak performance of the GPU; the cause was explained.

ASSESSING THE USEFULNESS OF THE QUASI-IDEAL OBSERVER FOR QUALITY CONTROL IN FLUOROSCOPY.

The aim of this work was to evaluate the reliability of the square of the signal-to-noise ratio rate, [Formula: see text], as a precise measurement for quality control test in a digital fluoroscopy system. The quasi-ideal model observer was used to measure [Formula: see text] The dose rate, pulse rate and field of view were varied, and their effect on dose efficiency, defined as [Formula: see text], was evaluated (where [Formula: see text] is the air kerma-area product rate). Measurements were repeated to assess reproducibility. The relative standard deviation in [Formula: see text] over seven consecutive measurements was 5 %. No significant variation in [Formula: see text] was observed across different pulse rates (10-30 pulses s(-1)). The low-dose-rate setting had a superior dose efficiency compared with the medium- and high-dose-rate settings. A smaller field of view resulted in higher dose efficiency.The results show that [Formula: see text] measurements offer the high precision required in quality control constancy tests.

In-situ calibration of clinical built-in KAP meters with traceability to a primary standard using a reference KAP meter.

The air kerma-area product (KAP) is used for settings of diagnostic reference levels. The International Atomic Energy Agency (IAEA) recommends that doses in diagnostic radiology (including the KAP values) be estimated with an accuracy of at least ± 7% (k = 2). Industry standards defined by the International Electrotechnical Commission (IEC) specify that the uncertainty of KAP meter measurements should be less than ± 25% (k = 2). Medical physicists willing to comply with the IAEA's recommendation need to apply correction factors to KAP values reported by x-ray units. The aim of this work is to present and evaluate a calibration method for built-in KAP meters on clinical x-ray units. The method is based on (i) a tandem calibration method, which uses a reference KAP meter calibrated to measure the incident radiation, (ii) measurements using an energy-independent ionization chamber to correct for the energy dependence of the reference KAP meter, and (iii) Monte Carlo simulations of the beam quality correction factors that correct for differences between beam qualities at a standard laboratory and the clinic. The method was applied to the KAP meter in a Siemens Aristos FX plus unit. It was found that values reported by the built-in KAP meter differed from the more accurate values measured by the reference KAP meter by more than 25% for high tube voltages (more than 140 kV) and heavily filtered beams (0.3 mm Cu). Associated uncertainties were too high to claim that the IEC's limit of 25% was exceeded. Nevertheless the differences were high enough to justify the need for a more accurate calibration of built-in KAP meters.

Quantifying the potential for dose reduction with visual grading regression.

OBJECTIVES: To propose a method to study the effect of exposure settings on image quality and to estimate the potential for dose reduction when introducing dose-reducing measures. METHODS: Using the framework of visual grading regression (VGR), a log(mAs) term is included in the ordinal logistic regression equation, so that the effect of reducing the dose can be quantitatively related to the effect of adding post-processing. In the ordinal logistic regression, patient and observer identity are treated as random effects using generalised linear latent and mixed models. The potential dose reduction is then estimated from the regression coefficients. The method was applied in a single-image study of coronary CT angiography (CTA) to evaluate two-dimensional (2D) adaptive filters, and in an image-pair study of abdominal CT to evaluate 2D and three-dimensional (3D) adaptive filters. RESULTS: For five image quality criteria in coronary CTA, dose reductions of 16-26% were predicted when adding 2D filtering. Using five image quality criteria for abdominal CT, it was estimated that 2D filtering permits doses were reduced by 32-41%, and 3D filtering by 42-51%. CONCLUSIONS: VGR including a log(mAs) term can be used for predictions of potential dose reduction that may be useful for guiding researchers in designing subsequent studies evaluating diagnostic value. With appropriate statistical analysis, it is possible to obtain direct numerical estimates of the dose-reducing potential of novel acquisition, reconstruction or post-processing techniques.

Quantifying the potential for dose reduction with visual grading regression.

Objectives To propose a method to study the effect of exposure settings on image quality and to estimate the potential for dose reduction when introducing dose-reducing measures. Methods Using the framework of visual grading regression (VGR), a log(mAs) term is included in the ordinal logistic regression equation, so that the effect of reducing the dose can be quantitatively related to the effect of adding post-processing. In the ordinal logistic regression, patient and observer identity are treated as random effects using generalised linear latent and mixed models. The potential dose reduction is then estimated from the regression coefficients. The method was applied in a single-image study of coronary CT angiography (CTA) to evaluate two-dimensional (2D) adaptive filters, and in an image-pair study of abdominal CT to evaluate 2D and three-dimensional (3D) adaptive filters. Results For five image quality criteria in coronary CTA, dose reductions of 16-26% were predicted when adding 2D filtering. Using five image quality criteria for abdominal CT, it was estimated that 2D filtering permits doses were reduced by 32-41%, and 3D filtering by 42-51%. Conclusions VGR including a log(mAs) term can be used for predictions of potential dose reduction that may be useful for guiding researchers in designing subsequent studies evaluating diagnostic value. With appropriate statistical analysis, it is possible to obtain direct numerical estimates of the dose-reducing potential of novel acquisition, reconstruction or post-processing techniques.

Patient organ radiation doses during treatment for aneurysmal subarachnoid hemorrhage.

PURPOSE: The aim of this retrospective study was to estimate risk organ doses and to estimate radiation risks during the imaging work-up and treatment for aneurysmal subarachnoid hemorrhage (SAH). METHODS: The imaging procedures comprised computed tomography and digital subtraction angiography studies for diagnosis or endovascular interventional procedures in 50 consecutive patients. Equivalent organ doses (H(T)) to skin, brain, eye lens, salivary glands, thyroid and oral mucosa were measured using thermoluminescence dosimeters in an anthropomorphic head phantom. Picture archiving and communication system (PACS) and radiological information system (RIS) records were analyzed and the frequency of each imaging procedure was recorded as well as the registered individual kerma-length product (P(KL)) and the kerma-area product (P(KA)). The doses were computed by multiplying the recorded P(KL) and P(KA) values by the conversion coefficients H(T)/P(KL) and H(T)/P(KA) from the head phantom. RESULTS: The mean fluoroscopy time, P(KL) and P(KA) were 38 min, 7269 mGy cm and 286 Gy cm(2), respectively. The estimated mean equivalent doses were as follows: skin 2.51 Sv, brain 0.92 Sv, eye lens 0.43 Sv and salivary glands 0.23 Sv. Maximum organ doses were 2.3-3.5 times higher than the mean. Interventional procedures contributed 66 % to skin dose, 55 % to brain dose and 25 % to eye lens dose. Of the patients with an estimated skin dose exceeding 6 Sv, only 1 developed temporary epilation. CONCLUSION: The risk for radiation-induced cancer for SAH patients is low (2-3 cases per 1,000 patients, of which 90 % are expected to be benign types) compared with the risk of tissue reactions on the head such as skin erythema and epilation (1 temporary epilation per 50 patients).

Breast dosimetry using high-resolution voxel phantoms.

A computer model of X-ray mammography has been developed, which uses quasi-realistic high-resolution voxel phantoms to simulate the breast. The phantoms have 400 microm voxels and simulate the three-dimensional distributions of adipose and fibro-glandular tissues, Cooper's ligaments, ducts and skin and allow the estimation of dose to individual tissues. Calculations of the incident air kerma to mean glandular dose conversion factor, g, were made using a Mo/Mo spectrum at 28 kV for eight phantoms in the thickness range 40-80 mm and of varying glandularity. The values differed from standard tabulations used for breast dosimetry by up to 43%, because of the different spatial distribution of glandular tissue within the breast. To study this further, additional voxel phantoms were constructed, which gave variations of between 9 and 59% compared with standard values. For accurate breast dosimetry, it is therefore very important to take the distribution of glandular tissues into account.

Calculation of the properties of digital mammograms using a computer simulation.

A Monte Carlo computer model of mammography has been developed to study and optimise the performance of digital mammographic systems. The program uses high-resolution voxel phantoms to model the breast, which simulate the adipose and fibroglandular tissues, Cooper's ligaments, ducts and skin in three dimensions. The model calculates the dose to each tissue, and also the quantities such as energy imparted to image pixels, noise per image pixel and scatter-to-primary (S/P) ratios. It allows studies of the dependence of image properties on breast structure and on position within the image. The program has been calibrated by calculating and measuring the pixel values and noise for a digital mammographic system. The thicknesses of two components of this system were unknown, and were adjusted to obtain a good agreement between measurement and calculation. The utility of the program is demonstrated with the calculations of the variation of the S/P ratio with and without a grid, and of the image contrast across the image of a 50-mm-thick breast phantom.

Monte Carlo simulation of a mammographic test phantom.

A test phantom, including a wide range of mammographic tissue equivalent materials and test details, was imaged on a digital mammographic system. In order to quantify the effect of scatter on the contrast obtained for the test details, calculations of the scatter-to-primary ratio (S/P) have been made using a Monte Carlo simulation of the digital mammographic imaging chain, grid and test phantom. The results show that the S/P values corresponding to the imaging conditions used were in the range 0.084-0.126. Calculated and measured pixel values in different regions of the image were compared as a validation of the model and showed excellent agreement. The results indicate the potential of Monte Carlo methods in the image quality-patient dose process optimisation, especially in the assessment of imaging conditions not available on standard mammographic units.

Bedside monitoring of CBF with xenon-CT and a mobile scanner: a novel method in neurointensive care.

Combining previously independently established techniques our objective was to develop and evaluate a method for bedside qualitative assessment of cerebral blood flow in neurointensive care (NICU) patients. The CT-protocol was optimized using phantoms and comparing a mobile CT-scanner (Tomoscan-M, Philips) with two stationary CT scanners. Thirty-two per cent xenon was delivered with standard equipment (Enhancer 3000). Mean cortical flow in volunteers was 48 ml/min/100 g, with the mean vascular territorial flow varying between 45 and 66 ml/min/100 g. The potential clinical usefulness was illustrated in three patients with vasospasm following subarachnoid haemorrhage. Our conclusion is that quantitative bedside measurements of CBF can be repeatedly performed in an easy and safe way in a standard NICU-setting, using xenon-inhalation and a mobile CT-scanner. The method is useful for the decision-making, and is a good example of how the quality of multi-modality monitoring in the NICU can be developed and further diversified.

Influence of the characteristic curve on the clinical image quality of lumbar spine and chest radiographs.

The "European Guidelines on Quality Criteria for Diagnostic Radiographic Images" do not address the choice of the film characteristic (H&D) curve, which is an important parameter for the description of a radiographic screen-film system. The image contrast of clinical lumbar spine and chest radiographs was altered by digital image processing techniques, simulating images with different H&D curves, both steeper and flatter than the original. The manipulated images were printed on film for evaluation. Seven experienced radiologists evaluated the clinical image quality by analysing the fulfilment of the European Image Criteria (ICS) and by visual grading analysis (VGA) of in total 224 lumbar spine and 360 chest images. A parallel study of the effect of the H&D curve has also been made using a theoretical model. The contrast (DeltaOD) of relevant anatomical details was calculated, using a Monte Carlo simulation-model of the complete imaging system including a 3D voxel phantom of a patient. Correlations between the calculated contrast and the radiologists' assessment by VGA were sought. The results of the radiologists' assessment show that the quality in selected regions of lumbar spine and chest images can be significantly improved by the use of films with a steeper H&D curve compared with the standard latitude film. Significant (p<0.05) correlations were found between the VGA results and the calculations of the contrast of transverse processes and trabecular details in the lumbar spine vertebrae, and with the contrast of blood vessels in the retrocardiac area of the chest.

A study and optimization of lumbar spine X-ray imaging systems.

A Monte Carlo program has been developed that incorporates a voxel phantom of an adult patient in a model of the complete X-ray imaging system, including the anti-scatter grid and screen-film receptor. This allows the realistic estimation of patient dose and the corresponding image (optical density map) for a wide range of equipment configurations. This paper focuses on the application of the program to lumbar spine anteroposterior and lateral screen-film examinations. The program has been applied to study the variation of physical image quality measures and effective dose for changing system parameters such as tube voltage, grid design and screen-film system speed. These variations form the basis for optimization of these system parameters. In our approach to optimization, the best systems are those that can match (or come close to) the calculated image quality measure of systems preferred in a recent European clinical trial, but with lower patient dose. The largest dose savings found were 21% for a 400 speed class system with a grid having a strip density of 40 cm(-1) and a grid ratio of 16. A further dose saving of 13% was possible when a 600 speed class system was employed. The best systems found from the optimization correspond to those recommended by the European Commission guidelines on image quality criteria for diagnostic radiographic images.

Schemes for the optimization of chest radiography using a computer model of the patient and x-ray imaging system.

A computer program has been developed to model chest radiography. It incorporates a voxel phantom of an adult and includes antiscatter grid, radiographic screen, and film. Image quality is quantified by calculating the contrast (deltaOD) and the ideal observer signal-to-noise ratio (SNR(I)) for a number of relevant anatomical details at various positions in the anatomy. Detector noise and system unsharpness are modeled and their influence on image quality is considered. A measure of useful dynamic range is computed and defined as the fraction of the image that is reproduced at an optical density such that the film gradient exceeds a preset value. The effective dose is used as a measure of the radiation risk for the patient. A novel approach to patient dose and image quality optimization has been developed and implemented. It is based on a reference system acknowledged to yield acceptable image quality in a clinical trial. Two optimizations schemes have been studied, the first including the contrast of vessels as measure of image quality and the second scheme using also the signal-to-noise ratio of calcifications. Both schemes make use of our measure of useful dynamic range as a key quantity. A large variety of imaging conditions was simulated by varying the tube voltage, antiscatter device, screen-film system, and maximum optical density in the computed image. It was found that the optical density is crucial in screen-film chest radiography. Significant dose savings (30%-50%) can be accomplished without sacrificing image quality by using low-atomic-number grids with a low grid ratio or an air gap and more sensitive screen-film system. Dose-efficient configurations proposed by the model agree well with the example of good radiographic technique suggested by the European Commission.

Demonstration of correlations between clinical and physical image quality measures in chest and lumbar spine screen-film radiography.

The ability to predict clinical image quality from physical measures is useful for optimization in diagnostic radiology. In this work, clinical and physical assessments of image quality are compared and correlations between the two are derived. Clinical assessment has been made by a group of expert radiologists who evaluated fulfillment of the European image criteria for chest and lumbar spine radiography using two scoring methods: image criteria score (ICS) and visual grading analysis score (VGAS). Physical image quality measures were calculated using a Monte Carlo simulation model of the complete imaging system. This model includes a voxelized male anatomy and was used to calculate contrast and signal-to-noise ratio of various important anatomical details and measures of dynamic range. Correlations between the physical image quality measures on the one hand and the ICS and VGAS on the other were sought. 16 chest and 4 lumbar spine imaging system configurations were compared in frontal projection. A statistically significant correlation with clinical image quality was found in chest posteroanterior radiography for the contrast of blood vessels in the retrocardiac area and a measure of useful dynamic range. In lumbar spine anteroposterior radiography, a similar significant correlation with clinical image quality was found between the contrast and signal-to-noise ratio of the trabecular structures in the L1-L5 vertebrae. The significant correlation shows that clinical image quality can, at least in some cases, be predicted from appropriate measures of physical image quality.

Influence of anode/filter material and tube potential on contrast, signal-to-noise ratio and average absorbed dose in mammography: a Monte Carlo study.

The comparative performance of mammographic X-ray systems that use different anode/filter combinations has been assessed for screen-film and digital imaging. Monte Carlo techniques have been used to calculate average glandular dose as well as contrast and signal-to-noise ratio for imaging two test details. Five anode/filter combinations have been studied to establish the potential for dose saving or image quality improvement. For screen-film mammography, it was found that little benefit is gained by changing from a standard 28 kV molybdenum/molybdenum spectrum for breasts up to 6 cm thick. For thicker breasts, where the tube potential for the standard technique might be increased, 20% improvement in contrast can be achieved without dose penalty using molybdenum/rhodium or rhodium/rhodium spectra, whereas dose savings of more than 50% can be attained whilst maintaining contrast using tungsten/rhodium or rhodium/aluminium spectra. In digital mammography, a molybdenum/molybdenum spectrum delivers the lowest dose for a 2 cm breast, but gives the highest dose for thicker breasts. Tungsten/rhodium or rhodium/aluminium spectra provide the lowest doses at greater thicknesses. It is concluded that for screen-film mammography, molybdenum/molybdenum is the spectrum of choice for all but the thickest or most glandular breasts. In digital mammography, an alternative spectrum is preferable for breasts thicker than 2 cm.