Towards Accurate Partial Volume Correction - Perturbation for SPECT Resolution Estimation

The accuracy of quantitative SPECT imaging is limited by the Partial Volume Effect as a result of the relatively poor spatial resolution. There is currently no consensus on the optimal Partial Volume Correction (PVC) algorithm in the application of SPECT oncology imaging. Several promising candidates require information on the reconstructed resolution - usually in the form of the Point Spread Function (PSF). A particular challenge that SPECT poses for PVC is that the resolution is known to vary with position in the field-of-view, as well as with activity distribution and reconstruction method. In this work, we assessed the potential benefit of using perturbation to measure case-specific resolution for PVC. A small point source was used to measure the resolution in phantoms designed to replicate the issues encountered in oncology imaging, including anthropomorphic phantoms which had not previously been examined in perturbation applications. Results demonstrate that, provided that a sufficient number of iterations is used for image reconstruction, perturbation can be used to measure a case-specific PSF. When PVC is applied with this case-specific PSF, quantitative accuracy is improved compared with no correction or applying PVC with an inappropriate PSF

Towards accurate partial volume correction in (99m}^Tc oncology SPECT: perturbation for case-specific resolution estimation

BACKGROUND: Currently, there is no consensus on the optimal partial volume correction (PVC) algorithm for oncology imaging. Several existing PVC methods require knowledge of the reconstructed resolution, usually as the point spread function (PSF)-often assumed to be spatially invariant. However, this is not the case for SPECT imaging. This work aimed to assess the accuracy of SPECT quantification when PVC is applied using a case-specific PSF. METHODS: Simulations of SPECT [Formula: see text]Tc imaging were performed for a range of activity distributions, including those replicating typical clinical oncology studies. Gaussian PSFs in reconstructed images were estimated using perturbation with a small point source. Estimates of the PSF were made in situations which could be encountered in a patient study, including; different positions in the field of view, different lesion shapes, sizes and contrasts, noise-free and noisy data. Ground truth images were convolved with the perturbation-estimated PSF, and with a PSF reflecting the resolution at the centre of the field of view. Both were compared with reconstructed images and the root-mean-square error calculated to assess the accuracy of the estimated PSF. PVC was applied using Single Target Correction, incorporating the perturbation-estimated PSF. Corrected regional mean values were assessed for quantitative accuracy. RESULTS: Perturbation-estimated PSF values demonstrated dependence on the position in the Field of View and the number of OSEM iterations. A lower root mean squared error was observed when convolution of the ground truth image was performed with the perturbation-estimated PSF, compared with convolution using a different PSF. Regional mean values following PVC using the perturbation-estimated PSF were more accurate than uncorrected data, or data corrected with PVC using an unsuitable PSF. This was the case for both simple and anthropomorphic phantoms. For the simple phantom, regional mean values were within 0.7% of the ground truth values. Accuracy improved after 5 or more OSEM iterations (10 subsets). For the anthropomorphic phantoms, post-correction regional mean values were within 1.6% of the ground truth values for noise-free uniform lesions. CONCLUSION: Perturbation using a simulated point source could potentially improve quantitative SPECT accuracy via the application of PVC, provided that sufficient reconstruction iterations are used

An update on computational anthropomorphic anatomical models

The prevalent availability of high-performance computing coupled with validated computerized simulation platforms as open-source packages have motivated progress in the development of realistic anthropomorphic computational models of the human anatomy. The main application of these advanced tools focused on imaging physics and computational internal/external radiation dosimetry research. This paper provides an updated review of state-of-the-art developments and recent advances in the design of sophisticated computational models of the human anatomy with a particular focus on their use in radiation dosimetry calculations. The consolidation of flexible and realistic computational models with biological data and accurate radiation transport modeling tools enables the capability to produce dosimetric data reflecting actual setup in clinical setting. These simulation methodologies and results are helpful resources for the medical physics and medical imaging communities and are expected to impact the fields of medical imaging and dosimetry calculations profoundly.</p

Consistent and invertible deformation vector fields for a breathing anthropomorphic phantom: a post-processing framework for the XCAT phantom.

Breathing motion is challenging for radiotherapy planning and delivery. This requires advanced four-dimensional (4D) imaging and motion mitigation strategies and associated validation tools with known deformations. Numerical phantoms such as the XCAT provide reproducible and realistic data for simulation-based validation. However, the XCAT generates partially inconsistent and non-invertible deformations where tumours remain rigid and structures can move through each other. We address these limitations by post-processing the XCAT deformation vector fields (DVF) to generate a breathing phantom with realistic motion and quantifiable deformation. An open-source post-processing framework was developed that corrects and inverts the XCAT-DVFs while preserving sliding motion between organs. Those post-processed DVFs are used to warp the first XCAT-generated image to consecutive time points providing a 4D phantom with a tumour that moves consistently with the anatomy, the ability to scale lung density as well as consistent and invertible DVFs. For a regularly breathing case, the inverse consistency of the DVFs was verified and the tumour motion was compared to the original XCAT. The generated phantom and DVFs were used to validate a motion-including dose reconstruction (MIDR) method using isocenter shifts to emulate rigid motion. Differences between the reconstructed doses with and without lung density scaling were evaluated. The post-processing framework produced DVFs with a maximum [Formula: see text]-percentile inverse-consistency error of 0.02 mm. The generated phantom preserved the dominant sliding motion between the chest wall and inner organs. The tumour of the original XCAT phantom preserved its trajectory while deforming consistently with the underlying tissue. The MIDR was compared to the ground truth dose reconstruction illustrating its limitations. MIDR with and without lung density scaling resulted in small dose differences up to 1 Gy (prescription 54 Gy). The proposed open-source post-processing framework overcomes important limitations of the original XCAT phantom and makes it applicable to a wider range of validation applications within radiotherapy

Virtual clinical trials in medical imaging: a review

The accelerating complexity and variety of medical imaging devices and methods have outpaced the ability to evaluate and optimize their design and clinical use. This is a significant and increasing challenge for both scientific investigations and clinical applications. Evaluations would ideally be done using clinical imaging trials. These experiments, however, are often not practical due to ethical limitations, expense, time requirements, or lack of ground truth. Virtual clinical trials (VCTs) (also known as in silico imaging trials or virtual imaging trials) offer an alternative means to efficiently evaluate medical imaging technologies virtually. They do so by simulating the patients, imaging systems, and interpreters. The field of VCTs has been constantly advanced over the past decades in multiple areas. We summarize the major developments and current status of the field of VCTs in medical imaging. We review the core components of a VCT: computational phantoms, simulators of different imaging modalities, and interpretation models. We also highlight some of the applications of VCTs across various imaging modalities

CT breast dose reduction with the use of breast positioning and organ‐based tube current modulation

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/136262/1/mp12076.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/136262/2/mp12076_am.pd

Inclusion of quasi-vertex views in a brain-dedicated multi-pinhole SPECT system for improved imaging performance

With brain-dedicated multi-detector systems employing pinhole apertures the usage of detectors facing the top of the patient\u27s head (i.e., quasi-vertex views) can provide the advantage of additional viewing from close to the brain for improved detector coverage. In this paper, we report the results of simulation and reconstruction studies to investigate the impact of the quasi-vertex views on the imaging performance of AdaptiSPECT-C, a brain-dedicated stationary SPECT system under development. In this design, both primary and scatter photons from regions located inferior to the brain can contribute to SPECT projections acquired by the quasi-vertex views, and thus degrade AdaptiSPECT-C imaging performance. In this work, we determined the proportion, origin, and nature (i.e., primary, scatter, and multiple-scatter) of counts emitted from structures within the head and throughout the body contributing to projections from the different AdaptiSPECT-C detector rings, as well as from a true vertex view detector. We simulated phantoms used to assess different aspects of image quality (i.e., uniform sphere and Derenzo), as well as anthropomorphic phantoms with multiple count levels emulating clinical(123)I activity distributions (i.e., DaTscan and perfusion). We determined that attenuation and scatter in the patient\u27s body greatly diminish the probability of the photons emitted outside the volume of interest reaching to detectors and being recorded within the 15% photopeak energy window. In addition, we demonstrated that the inclusion of the residual of such counts in the system acquisition does not degrade visual interpretation or quantitative analysis. The addition of the quasi-vertex detectors increases volumetric sensitivity, angular sampling, and spatial resolution leading to significant enhancement in image quality, especially in the striato-thalamic and superior regions of the brain. Besides, the use of quasi-vertex detectors improves the recovery of clinically relevant metrics such as the striatal binding ratio and mean activity in selected cerebral structures

Personalised body counter calibration using anthropometric parameters

This book describes the development of a new method for personalisation of efficiency factors in partial body counting. Its achieved goal is the quantification of uncertainties in those factors due to variation in anatomy of the measured persons, and their reduction by correlation with anthropometric parameters. The method was applied to a detector system at the In Vivo Measurement Laboratory at Karlsruhe Institute of Technology using Monte Carlo simulation and computational phantoms

Joint Activity and Attenuation Reconstruction From Multiple Energy Window Data With Photopeak Scatter Re-Estimation in Non-TOF 3-D PET

Estimation of attenuation from positron emission tomography (PET) data only is of interest for hybrid PET-MR and systems where CT is not available or recommended. However, when using data from a single energy window, emission-based non-time-of-flight (TOF) PET attenuation correction (AC) methods suffer from “cross-talk” artifacts. Based on earlier work, this article explores the hypothesis that cross-talk can be reduced by using more than one energy window. We propose an algorithm for the simultaneous estimation of both activity and attenuation images, as well as, the scatter component of the measured data from a PET acquisition, using multiple energy windows. The model for the measurements is 3-D and accounts for the finite energy resolution of PET detectors; it is restricted to single scatter. The proposed energy-based simultaneous maximum likelihood reconstruction of activity and attenuation with photopeak scatter re-estimation algorithm is compared with simultaneous estimation from a single energy window simultaneous maximum likelihood reconstruction of activity and attenuation with photopeak scatter re-estimation. The evaluation is based on simulations using the characteristics of the Siemens mMR scanner. Phantoms of different complexity were investigated. In particular, a 3-D XCAT torso phantom was used to assess the inpainting of attenuation values within the lung region. Results show that the cross-talk present in non-TOF maximum likelihood reconstruction of activity and attenuation reconstructions is significantly reduced when using multiple energy windows and indicate that the proposed approach warrants further investigation