Optimizing dual energy cone beam CT protocols for preclinical imaging and radiation research

Objective: The aim of this work was to investigate whether quantitative dual-energy CT (DECT) imaging is feasible for small animal irradiators with an integrated cone-beam CT (CBCT) system. Methods: The optimal imaging protocols were determined by analyzing different energy combinations and dose levels. The influence of beam hardening effects and the performance of a beam hardening correction (BHC) were investigated. In addition, two systems from different manufacturers were compared in terms of errors in the extracted effective atomic numbers (Z(eff)) and relative electron densities (rho(e)) for phantom inserts with known elemental compositions and relative electron densities. Results: The optimal energy combination was determined to be 50 and 90kVp. For this combination, Z(eff) and r rho(e) can be extracted with a mean error of 0.11 and 0.010, respectively, at a dose level of 60cGy. Conclusion: Quantitative DECT imaging is feasible for small animal irradiators with an integrated CBCT system. To obtain the best results, optimizing the imaging protocols is required. Well-separated X-ray spectra and a sufficient dose level should be used to minimize the error and noise for Z(eff) and rho(e). When no BHC is applied in the image reconstruction, the size of the calibration phantom should match the size of the imaged object to limit the influence of beam hardening effects. No significant differences in Z(eff) and rho(e) errors are observed between the two systems from different manufacturers. Advances in knowledge: This is the first study that investigates quantitative DECT imaging for small animal irradiators with an integrated CBCT system

Segmentation-Free Statistical Image Reconstruction for Polyenergetic X-Ray Computed Tomography with Experimental Validation

This paper describes a statistical image reconstruction method for x-ray CT that is based on a physical model that accounts for the polyenergetic x-ray source spectrum and the measurement nonlinearities caused by energy-dependent attenuation. Unlike our earlier work, the proposed algorithm does not require pre-segmentation of the object into the various tissue classes (e.g., bone and soft tissue) and allows mixed pixels. The attenuation coefficient of each voxel is modelled as the product of its unknown density and a weighted sum of energy-dependent mass attenuation coefficients. We formulate a penalized-likelihood function for this polyenergetic model and develop an iterative algorithm for estimating the unknown density of each voxel. Applying this method to simulated x-ray CT measurements of objects containing both bone and soft tissue yields images with significantly reduced beam hardening artefacts relative to conventional beam hardening correction methods. We also apply the method to real data acquired from a phantom containing various concentrations of potassium phosphate solution. The algorithm reconstructs an image with accurate density values for the different concentrations, demonstrating its potential for quantitative CT applications.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85911/1/Fessler66.pd

Radiation Metrology of Small Animal Molecular Imaging and Molecular Radiotherapy Using mirco-PET/CT

Genetically engineered animal models of diseases are increasingly recapitulating human diseases. With this, in vivo preclinical imaging of small laboratory animals has emerged as a critical component of biomedical research because of its noninvasive nature allowing serial assay of animal models and monitoring its safety and effectiveness over the history of the disease. The concept of quantitative molecular imaging is to go beyond displaying images in digital form and to consider the image and extract quantitative information that allows for a better understanding of disease progression and treatment. The aim of this work is to demonstrate the need for the metrology of molecular imaging of animal models using micro-PET/CT devices. System characteristics are determined within each subsystem, micro-PET and micro-CT, independent of each other, and as integrated systems. The characterization of tissues, composition and density, by micro-CT was determined along with the noise level of the unit. Moreover, the nominal superficial and deep absorbed doses were estimated to assess the confounding effect of multiple scans in animal studies. The Q value, used to convert counts per milliliter to activity per milliliter, was estimated to assess the observed activity present in the animal. The resolution of the micro-PET subsystem was also estimated using a modified Derenzo phantom to assess the uncertainty of the activity distribution within tissues. Once both modalities were characterized separately the coordinate system of each individual system was checked for spatial accuracy using a cross capillary method. The offset values were then used to establish the same coordinate system for co-registration. Once both micro-PET and micro-CT image data sets had been verified, they were used to generate a voxel image of the subject for use in the Monte Carlo program, MCNP6, where an absorbed dose map was generated for the radiolabeled compound. Two basic examples are given to demonstrate the use of the voxelized absorbed dose maps for calculating the absorbed dose to any segmented organ of interest, across longitudinal studies. In this way, it was shown that an animal-specific model can be used to accurately calculate the absorbed dose for each time point during a study

Developments in PET-MRI for Radiotherapy Planning Applications

The hybridization of magnetic resonance imaging (MRI) and positron emission tomography (PET) provides the benefit of soft-tissue contrast and specific molecular information in a simultaneous acquisition. The applications of PET-MRI in radiotherapy are only starting to be realised. However, quantitative accuracy of PET relies on accurate attenuation correction (AC) of, not only the patient anatomy but also MRI hardware and current methods, which are prone to artefacts caused by dense materials. Quantitative accuracy of PET also relies on full characterization of patient motion during the scan. The simultaneity of PET-MRI makes it especially suited for motion correction. However, quality assurance (QA) procedures for such corrections are lacking. Therefore, a dynamic phantom that is PET and MR compatible is required. Additionally, respiratory motion characterization is needed for conformal radiotherapy of lung. 4D-CT can provide 3D motion characterization but suffers from poor soft-tissue contrast. In this thesis, I examine these problems, and present solutions in the form of improved MR-hardware AC techniques, a PET/MRI/CT-compatible tumour respiratory motion phantom for QA measurements, and a retrospective 4D-PET-MRI technique to characterise respiratory motion. Chapter 2 presents two techniques to improve upon current AC methods that use a standard helical CT scan for MRI hardware in PET-MRI. One technique uses a dual-energy computed tomography (DECT) scan to construct virtual monoenergetic image volumes and the other uses a tomotherapy linear accelerator to create CT images at megavoltage energies (1.0 MV) of the RF coil. The DECT-based technique reduced artefacts in the images translating to improved μ-maps. The MVCT-based technique provided further improvements in artefact reduction, resulting in artefact free μ-maps. This led to more AC of the breast coil. In chapter 3, I present a PET-MR-CT motion phantom for QA of motion-correction protocols. This phantom is used to evaluate a clinically available real-time dynamic MR images and a respiratory-triggered PET-MRI protocol. The results show the protocol to perform well under motion conditions. Additionally, the phantom provided a good model for performing QA of respiratory-triggered PET-MRI. Chapter 4 presents a 4D-PET/MRI technique, using MR sequences and PET acquisition methods currently available on hybrid PET/MRI systems. This technique is validated using the motion phantom presented in chapter 3 with three motion profiles. I conclude that our 4D-PET-MRI technique provides information to characterise tumour respiratory motion while using a clinically available pulse sequence and PET acquisition method

Reducing Radiation Dose to the Female Breast during CT Coronary Angiography: A Simulation Study Comparing Breast Shielding, Angular Tube Current Modulation, Reduced kV, and Partial Angle Protocols Using an Unknown-location Signal-detectability Metric

Purpose: The authors compared the performance of five protocols intended to reduce dose to the breast during computed tomography (CT) coronary angiography scans using a model observer unknown-location signal-detectability metric. Methods: The authors simulated CT images of an anthropomorphic female thorax phantom for a 120 kV reference protocol and five “dose reduction” protocols intended to reduce dose to the breast: 120 kV partial angle (posteriorly centered), 120 kV tube-current modulated (TCM), 120 kV with shielded breasts, 80 kV, and 80 kV partial angle (posteriorly centered). Two image quality tasks were investigated: the detection and localization of 4-mm, 3.25 mg/ml and 1-mm, 6.0 mg/ml iodine contrast signals randomly located in the heart region. For each protocol, the authors plotted the signal detectability, as quantified by the area under the exponentially transformed free response characteristic curve estimator (AˆFE), as well as noise and contrast-to-noise ratio (CNR) versus breast and lung dose. In addition, the authors quantified each protocol\u27s dose performance as the percent difference in dose relative to the reference protocol achieved while maintaining equivalentAˆFE. Results: For the 4-mm signal-size task, the 80 kV full scan and 80 kV partial angle protocols decreased dose to the breast (80.5% and 85.3%, respectively) and lung (80.5% and 76.7%, respectively) withAˆFE= 0.96, but also resulted in an approximate three-fold increase in image noise. The 120 kV partial protocol reduced dose to the breast (17.6%) at the expense of increased lung dose (25.3%). The TCM algorithm decreased dose to the breast (6.0%) and lung (10.4%). Breast shielding increased breast dose (67.8%) and lung dose (103.4%). The 80 kV and 80 kV partial protocols demonstrated greater dose reductions for the 4-mm task than for the 1-mm task, and the shielded protocol showed a larger increase in dose for the 4-mm task than for the 1-mm task. In general, the CNR curves indicate a similar relative ranking of protocol performance as the correspondingAˆFEcurves, however, the CNR metric overestimated the performance of the shielded protocol for both tasks, leading to corresponding underestimates in the relative dose increases compared to those obtained when using theAˆFEmetric. Conclusions: The 80 kV and 80 kV partial angle protocols demonstrated the greatest reduction to breast and lung dose, however, the subsequent increase in image noise may be deemed clinically unacceptable. Tube output for these protocols can be adjusted to achieve a more desirable noise level with lesser breast dose savings. Breast shielding increased breast and lung dose when maintaining equivalentAˆFE. The results demonstrated that comparisons of dose performance depend on both the image quality metric and the specific task, and that CNR may not be a reliable metric of signal detectability

Multi-Energy Blended CBCT Spectral Imaging Using a Spectral Modulator with Flying Focal Spot (SMFFS)

Cone-beam CT (CBCT) spectral imaging has great potential in medical and industrial applications, but it is very challenging as scatter and spectral effects are seriously twisted. In this work, we present the first attempt to develop a stationary spectral modulator with flying focal spot (SMFFS) technology as a promising, low-cost approach to accurately solving the X-ray scattering problem and physically enabling spectral imaging in a unified framework, and with no significant misalignment in data sampling of spectral projections. Based on an in-depth analysis of optimal energy separation from different combinations of modulator materials and thicknesses, we present a practical design of a mixed two-dimensional spectral modulator that can generate multi-energy blended CBCT spectral projections. To deal with the twisted scatter-spectral challenge, we propose a novel scatter-decoupled material decomposition (SDMD) method by taking advantage of a scatter similarity in SMFFS. A Monte Carlo simulation is conducted to validate the strong similarity of X-ray scatter distributions across the flying focal spot positions. Both numerical simulations using a clinical abdominal CT dataset, and physics experiments on a tabletop CBCT system using a GAMMEX multi-energy CT phantom, are carried out to demonstrate the feasibility of our proposed SDMD method for CBCT spectral imaging with SMFFS. In the physics experiments, the mean relative errors in selected ROI for virtual monochromatic image (VMI) are 0.9\% for SMFFS, and 5.3\% and 16.9\% for 80/120 kV dual-energy cone-beam scan with and without scatter correction, respectively. Our preliminary results show that SMFFS can effectively improve the quantitative imaging performance of CBCT.Comment: 10 pages, 13 figure

Advances in dual-energy computed tomography imaging of radiological properties

Dual-energy computed tomography (DECT) has shown great potential in the reduction of uncertainties of proton ranges and low energy photon cross section estimation used in radiation therapy planning. The work presented herein investigated three contributions for advancing DECT applications. 1) A linear and separable two-parameter DECT, the basis vector model (BVM) was used to estimate proton stopping power. Compared to other nonlinear two-parameter models in the literature, the BVM model shows a comparable accuracy achieved for typical human tissues. This model outperforms other nonlinear models in estimations of linear attenuation coefficients. This is the first study to clearly illustrate the advantages of linear model not only in accurately mapping radiological quantities for radiation therapy, but also in providing a unique model for accurate linear forward projection modelling, which is needed by the statistical iterative reconstruction (SIR) and other advanced DECT reconstruction algorithms. 2) Accurate DECT requires knowledge of x-ray beam properties. Using the Birch-Marshall1 model and beam hardening correction coefficients encoded in a CT scanner’s sinogram header files, an efficient and accurate way to estimate the x-ray spectrum is proposed. The merits of the proposed technique lie in requiring no physical transmission measurement after a one-time calibration against an independently measured spectrum. This technique can also be used in monitoring the aging of x-ray CT tubes. 3) An iterative filtered back projection with anatomical constraint (iFBP-AC) algorithm was also implemented on a digital phantom to evaluate its ability in mitigating beam hardening effects and supporting accurate material decomposition for in vivo imaging of photon cross section and proton stopping power. Compared to iFBP without constraints, both algorithms demonstrate high efficiency of convergence. For an idealized digital phantom, similar accuracy was observed under a noiseless situation. With clinically achievable noise level added to the sinograms, iFBP-AC greatly outperforms iFBP in prediction of photon linear attenuation at low energy, i.e., 28 keV. The estimated mean errors of iFBP and iFBP-AC for cortical bone are 1% and 0.7%, respectively; the standard deviations are 0.6% and 5%, respectively. The achieved accuracy of iFBP-AC shows robustness versus contrast level. Similar mean errors are maintained for muscle tissue. The standard deviation achieved by iFBP-AC is 1.2%. In contrast, the standard deviation yielded by iFBP is about 20.2%. The algorithm of iFBP-AC shows potential application of quantitative measurement of DECT. The contributions in this thesis aim to improve the clinical performance of DECT

Mammographic density. Measurement of mammographic density

Mammographic density has been strongly associated with increased risk of breast cancer. Furthermore, density is inversely correlated with the accuracy of mammography and, therefore, a measurement of density conveys information about the difficulty of detecting cancer in a mammogram. Initial methods for assessing mammographic density were entirely subjective and qualitative; however, in the past few years methods have been developed to provide more objective and quantitative density measurements. Research is now underway to create and validate techniques for volumetric measurement of density. It is also possible to measure breast density with other imaging modalities, such as ultrasound and MRI, which do not require the use of ionizing radiation and may, therefore, be more suitable for use in young women or where it is desirable to perform measurements more frequently. In this article, the techniques for measurement of density are reviewed and some consideration is given to their strengths and limitations

Optimization of the energy for Breast monochromatic absorption X-ray Computed Tomography

The limits of mammography have led to an increasing interest on possible alternatives such as the breast Computed Tomography (bCT). The common goal of all X-ray imaging techniques is to achieve the optimal contrast resolution, measured through the Contrast to Noise Ratio (CNR), while minimizing the radiological risks, quantified by the dose. Both dose and CNR depend on the energy and the intensity of the X-rays employed for the specific imaging technique. Some attempts to determine an optimal energy for bCT have suggested the range 22keV\u201334keV, some others instead suggested the range 50keV\u201360keV depending on the parameters considered in the study. Recent experimental works, based on the use of monochromatic radiation and breast specimens, show that energies around 32keV give better image quality respect to setups based on higher energies. In this paper we report a systematic study aiming at defining the range of energies that maximizes the CNR at fixed dose in bCT. The study evaluates several compositions and diameters of the breast and includes various reconstruction algorithms as well as different dose levels. The results show that a good compromise between CNR and dose is obtained using energies around 28keV