MB-DECTNet: a model-based unrolling network for accurate 3D dual-energy CT reconstruction from clinically acquired helical scans.

Objective.Over the past several decades, dual-energy CT (DECT) imaging has seen significant advancements due to its ability to distinguish between materials. DECT statistical iterative reconstruction (SIR) has exhibited potential for noise reduction and enhanced accuracy. However, its slow convergence and substantial computational demands render the elapsed time for 3D DECT SIR often clinically unacceptable. The objective of this study is to accelerate 3D DECT SIR while maintaining subpercentage or near-subpercentage accuracy.Approach.We incorporate DECT SIR into a deep-learning model-based unrolling network for 3D DECT reconstruction (MB-DECTNet), which can be trained end-to-end. This deep learning-based approach is designed to learn shortcuts between initial conditions and the stationary points of iterative algorithms while preserving the unbiased estimation property of model-based algorithms. MB-DECTNet comprises multiple stacked update blocks, each containing a data consistency layer (DC) and a spatial mixer layer, with the DC layer functioning as a one-step update from any traditional iterative algorithm.Main results.The quantitative results indicate that our proposed MB-DECTNet surpasses both the traditional image-domain technique (MB-DECTNet reduces average bias by a factor of 10) and a pure deep learning method (MB-DECTNet reduces average bias by a factor of 8.8), offering the potential for accurate attenuation coefficient estimation, akin to traditional statistical algorithms, but with considerably reduced computational costs. This approach achieves 0.13% bias and 1.92% mean absolute error and reconstructs a full image of a head in less than 12 min. Additionally, we show that the MB-DECTNet output can serve as an initializer for DECT SIR, leading to further improvements in results.Significance.This study presents a model-based deep unrolling network for accurate 3D DECT reconstruction, achieving subpercentage error in estimating virtual monoenergetic images for a full head at 60 and 150 keV in 30 min, representing a 40-fold speedup compared to traditional approaches. These findings have significant implications for accelerating DECT SIR and making it more clinically feasible.

Motion-compensated scheme for sequential scanned statistical iterative dual-energy CT reconstruction.

Objective.Dual-energy computed tomography (DECT) has been widely used to reconstruct numerous types of images due its ability to better discriminate tissue properties. Sequential scanning is a popular dual-energy data acquisition method as it requires no specialized hardware. However, patient motion between two sequential scans may lead to severe motion artifacts in DECT statistical iterative reconstructions (SIR) images. The objective is to reduce the motion artifacts in such reconstructions.Approach.We propose a motion-compensation scheme that incorporates a deformation vector field into any DECT SIR. The deformation vector field is estimated via the multi-modality symmetric deformable registration method. The precalculated registration mapping and its inverse or adjoint are then embedded into each iteration of the iterative DECT algorithm.Main results.Results from a simulated and clinical case show that the proposed framework is capable of reducing motion artifacts in DECT SIRs. Percentage mean square errors in regions of interest in the simulated and clinical cases were reduced from 4.6% to 0.5% and 6.8% to 0.8%, respectively. A perturbation analysis was then performed to determine errors in approximating the continuous deformation by using the deformation field and interpolation. Our findings show that errors in our method are mostly propagated through the target image and amplified by the inverse matrix of the combination of the Fisher information and Hessian of the penalty term.Significance.We have proposed a novel motion-compensation scheme to incorporate a 3D registration method into the joint statistical iterative DECT algorithm in order to reduce motion artifacts caused by inter-scan motion, and successfully demonstrate that interscan motion corrections can be integrated into the DECT SIR process, enabling accurate imaging of radiological quantities on conventional SECT scanners, without significant loss of either computational efficiency or accuracy.

Towards subpercentage uncertainty proton stopping-power mapping via dual-energy CT: Direct experimental validation and uncertainty analysis of a statistical iterative image reconstruction method.

PURPOSE: To assess the potential of a joint dual-energy computerized tomography (CT) reconstruction process (statistical image reconstruction method built on a basis vector model (JSIR-BVM)) implemented on a 16-slice commercial CT scanner to measure high spatial resolution stopping-power ratio (SPR) maps with uncertainties of less than 1%. METHODS: JSIR-BVM was used to reconstruct images of effective electron density and mean excitation energy from dual-energy CT (DECT) sinograms for 10 high-purity samples of known density and atomic composition inserted into head and body phantoms. The measured DECT data consisted of 90 and 140 kVp axial sinograms serially acquired on a Philips Brilliance Big Bore CT scanner without beam-hardening corrections. The corresponding SPRs were subsequently measured directly via ion chamber measurements on a MEVION S250 superconducting synchrocyclotron and evaluated theoretically from the known sample compositions and densities. Deviations of JSIR-BVM SPR values from their theoretically calculated and directly measured ground-truth values were evaluated for our JSIR-BVM method and our implementation of the Hünemohr-Saito (H-S) DECT image-domain decomposition technique for SPR imaging. A thorough uncertainty analysis was then performed for five different scenarios (comparison of JSIR-BVM stopping-power ratio/stopping power (SPR/SP) to International Commission on Radiation Measurements and Units benchmarks; comparison of JSIR-BVM SPR to measured benchmarks; and uncertainties in JSIR-BVM SPR/SP maps for patients of unknown composition) per the Joint Committee for Guides in Metrology and the Guide to Expression of Uncertainty in Measurement, including the impact of uncertainties in measured photon spectra, sample composition and density, photon cross section and I-value models, and random measurement uncertainty. Estimated SPR uncertainty for three main tissue groups in patients of unknown composition and the weighted proportion of each tissue type for three proton treatment sites were then used to derive a composite range uncertainty for our method. RESULTS: Mean JSIR-BVM SPR estimates deviated by less than 1% from their theoretical and directly measured ground-truth values for most inserts and phantom geometries except for high-density Delrin and Teflon samples with SPR error relative to proton measurements of 1.1% and -1.0% (head phantom) and 1.1% and -1.1% (body phantom). The overall root-mean-square (RMS) deviations over all samples were 0.39% and 0.52% (head phantom) and 0.43% and 0.57% (body phantom) relative to theoretical and directly measured ground-truth SPRs, respectively. The corresponding RMS (maximum) errors for the image-domain decomposition method were 2.68% and 2.73% (4.68% and 4.99%) for the head phantom and 0.71% and 0.87% (1.37% and 1.66%) for the body phantom. Compared to H-S SPR maps, JSIR-BVM yielded 30% sharper and twofold sharper images for soft tissues and bone-like surrogates, respectively, while reducing noise by factors of 6 and 3, respectively. The uncertainty (coverage factor k = 1) of the DECT-to-benchmark values comparison ranged from 0.5% to 1.5% and is dominated by scanning-beam photon-spectra uncertainties. An analysis of the SPR uncertainty for patients of unknown composition showed a JSIR-BVM uncertainty of 0.65%, 1.21%, and 0.77% for soft-, lung-, and bony-tissue groups which led to a composite range uncertainty of 0.6-0.9%. CONCLUSIONS: Observed JSIR-BVM SPR estimation errors were all less than 50% of the estimated k = 1 total uncertainty of our benchmarking experiment, demonstrating that JSIR-BVM high spatial resolution, low-noise SPR mapping is feasible and is robust to variations in the geometry of the scanned object. In contrast, the much larger H-S SPR estimation errors are dominated by imaging noise and residual beam-hardening artifacts. While the uncertainties characteristic of our current JSIR-BVM implementation can be as large as 1.5%, achieving < 1% total uncertainty is feasible by improving the accuracy of scanner-specific scatter-profile and photon-spectrum estimates. With its robustness to beam-hardening artifact, image noise, and variations in phantom size and geometry, JSIR-BVM has the potential to achieve high spatial-resolution SPR mapping with subpercentage accuracy and estimated uncertainty in the clinical setting.

Dosimetric impact of range uncertainty in passive scattering proton therapy.

PURPOSE: The objective of this study was to investigate the dosimetric impact of range uncertainty in a large cohort of patients receiving passive scatter proton therapy. METHODS: A cohort of 120 patients were reviewed in this study retrospectively, of which 61 were brain, 39 lung, and 20 prostate patients. Range uncertainties of ±3.5% (overshooting and undershooting by 3.5%, respectively) were added and recalculated on the original plans, which had been planned according to our clinical planning protocol while keeping beamlines, apertures, compensators, and dose grids intact. Changes in the coverage on CTV and DVH for critical organs were compared and analyzed. Correlation between dose change and minimal distance between CTV and critical organs were also investigated. RESULTS: Although CTV coverages and maximum dose to critical organs were largely maintained for most brain patients, large variations over 5% were still observed sporadically. Critical organs, such as brainstem and chiasm, could still be affected by range uncertainty at 4 cm away from CTV. Coverage and OARs in lung and prostate patients were less likely to be affected by range uncertainty with very few exceptions. CONCLUSION: The margin recipe in modern TPS leads to clinically acceptable OAR doses in the presence of range uncertainties. However, range uncertainties still pose a noticeable challenge for small but critical serial organs near tumors, and occasionally for large parallel organs that are located distal to incident proton beams.

Performance comparison of a dedicated total breast PET system with a clinical whole-body PET system: a simulation study.

This paper presents a novel PET geometry for breast cancer imaging. The scanner consists of a 'stadium' (a rectangle with two semi-circles on opposite sides) shaped ring, along with anterior and posterior panels to provide high sensitivity and high spatial resolution for an imaging field-of-view (FOV) that include both breasts, mediastinum and axilla. We simulated this total-breast PET system using GATE and reconstructed the coincidence events using a GPU-based list-mode image reconstruction implementing maximum likelihood expectation-maximization (ML-EM) algorithm. The rear-panel is made up of a single layer of LSO crystals (3.2 × 3.2 × 20 mm3each), while the 'stadium'-shaped elongated ring and the anterior panel are made with dual-layered LSO crystals (1.6 × 1.6 × 6 mm3each). The energy resolution and coincidence resolving time of all detectors are assumed to be 12% and 250 ps full-width-at-half-maximum, respectively. Various sized simulated lesions (4, 5, 6 mm) having 4:1, 5:1, and 6:1 lesion-to-background radioactivity concentration ratios, mimicking different biological uptakes, were strategically located throughout a volumetric torso phantom. We compared system sensitivity and lesion detectability of the dedicated total-breast PET system to a state-of-the-art clinical whole-body PET scanner. The mean sensitivity of the total-breast PET system is 3.21 times greater than that of a whole-body PET scanner in the breast regions. The total-breast PET system also provides better contrast-recovery coefficients for lesions of all sizes and lesion-to-background ratios in the breast when compared to a reference clinical whole-body PET scanner. Receiver operating characteristics (ROC) study shows the area under the ROC curve is 0.948 and 0.924 for the total-breast system and the whole-body PET scanner, respectively, in the detection of 4 mm diameter lesions with 4:1 lesion-to-background ratio. This study demonstrates our novel geometry can provide an imaging FOV larger than conventional PEM systems to simultaneously image both breasts, chest wall and axillae with significantly improved lesion detectability in the breasts when compared to a whole-body PET scanner.

Impact of bowtie filter and detector collimation on multislice CT scatter profiles: A simulation study.

PURPOSE: To investigate via Monte Carlo simulations, the impact of scan subject size, antiscatter grid (ASG), collimator size, and bowtie filter on the distribution of scatter radiation in a typical realistically modeled third generation 16 slice diagnostic computed tomography (CT) scanner. METHODS: Full radiation transport was simulated with Geant4 in a realistic CT scanner geometric model, including the imaging phantom, bowtie filter (BTF), collimators and detector assembly, except for the ASGs. An analytical method was employed to quantify the probable transmission through the ASG of each photon intersecting the detector array. Normalized scatter profiles (NSP) and scatter-to-primary-ratio (SPR) profiles were simulated for 90 and 140 kVp beams for different size phantoms and slice thicknesses. The impact of CT scatter on the reconstructed attenuation coefficient factor was also studied as were the modulating effects of phantom- and patient-tissue heterogeneities on scatter profiles. A method to characterize the relative spatial frequency content of sinogram signals was developed to assess the latter. RESULTS: For the 21.4-cm diameter phantom, NSP and SPR increase linearly with collimator opening for both tube potentials, with the 90 kVp scan exhibiting slightly larger NSP and SPR. The BTF modestly modulates scatter under the phantom center, reducing the prominent off-axis lobes by factors of 1.1-1.3. The ASG reduces scatter on the central axis NSP threefold, and reduces scatter at the detectors outside the phantom shadow by factors of 25 to 500. For the phantoms with diameters of 27 and 32 cm, the scatter increases roughly three- and fourfold, respectively, demonstrating that scatter monotonically increases with phantom size, despite deployment of the ASG and BTF. In the absence of a scan subject, the ASG reduces the signal profile arising photons scattered by the BTF. Without ASG, the in-air scatter profile is relatively flat compared to the scatter profile when the ASG is present. For both 90 and 140 kVp photon spectra, the calculated attenuation coefficient decreases linearly with increasing collimation size. For both homogeneous and heterogeneous objects, NSPs are dominated by low spatial frequency content compared to the primary signal. However, the SPR, which quantifies the local magnitude of nonlinear detector response and is dominated by the high frequency content of the primary profile, can contribute strongly to high-spatial frequency streaking artifacts near high-density structures in reconstructed image artifacts. CONCLUSION: Public-domain Monte Carlo codes, Geant-4 in particular, is a feasible method for characterizing CT detector response to scattered- and off-focal radiation. Our study demonstrates that the ASG substantially reduces the scatter radiation and reshapes scatter-radiation profiles and affects the accuracy with which the detector array can measure narrow-beam attenuation due its inability to distinguish between true uncollided primary and narrow-angle coherently scattered photons. Hence, incorporating the impact of detector array collimation into the forward-projection signal formation models used by iterative reconstruction algorithms is necessary to use CT for accurately characterizing material properties. While tissue heterogeneities exercise a modest influence on local NPS shape and magnitude, they do not add significant high spatial frequency content.

A second-generation virtual-pinhole PET device for enhancing contrast recovery and improving lesion detectability of a whole-body PET/CT scanner.

PURPOSE: We have developed a second-generation virtual-pinhole (VP) positron emission tomography (PET) device that can position a flat-panel PET detector around a patient's body using a robotic arm to enhance the contrast recovery coefficient (CRC) and detectability of lesions in any region-of-interest using a whole-body PET/computed tomography (CT) scanner. METHODS: We constructed a flat-panel VP-PET device using 32 high-resolution detectors, each containing a 4  ×  4 MPPC array and 16  ×  16 LYSO crystals of 1.0  ×  1.0  ×  3.0 mm3 each. The flat-panel detectors can be positioned around a patient's body anywhere in the imaging field-of-view (FOV) of a Siemens Biograph 40 PET/CT scanner by a robotic arm. New hardware, firmware and software have been developed to support the additional detector signals without compromising a scanner's native functions. We stepped a 22 Na point source across the axial FOV of the scanner to measure the sensitivity profile of the VP-PET device. We also recorded the coincidence events measured by the scanner detectors and by the VP-PET detectors when imaging phantoms of different sizes. To assess the improvement in the CRC of small lesions, we imaged an elliptical torso phantom measuring 316  ×  228  ×  162 mm3 that contains spherical tumors with diameters ranging from 3.3 to 11.4 mm with and without the VP-PET device. Images were reconstructed using a list mode Maximum-Likelihood Estimation-Maximization algorithm implemented on multiple graphics processing units (GPUs) to support the unconventional geometries enabled by a VP-PET system. The mean and standard deviation of the CRC were calculated for tumors of different sizes. Monte Carlo simulation was also conducted to image clusters of lesions in a torso phantom using a PET/CT scanner alone or the same scanner equipped with VP-PET devices. Receiver operating characteristic (ROC) curves were analyzed for three system configurations to evaluate the improvement in lesion detectability by the VP-PET device over the native PET/CT scanner. RESULTS: The repeatability in positioning the flat-panel detectors using a robotic arm is better than 0.15 mm in all three directions. Experimental results show that the average CRC of 3.3, 4.3, and 6.0 mm diameter tumors was 0.82%, 2.90%, and 5.25%, respectively, when measured by the native scanner. The corresponding CRC was 2.73%, 6.21% and 10.13% when imaged by the VP-PET insert device with the flat-panel detector under the torso phantom. These values may be further improved to 4.31%, 9.65% and 18.01% by a future dual-panel VP-PET insert device if DOI detectors are employed to triple its detector efficiency. Monte Carlo simulation results show that the tumor detectability can be improved by a VP-PET device that has a single flat-panel detector. The improvement is greater if the VP-PET device employs a dual-panel design. CONCLUSIONS: We have developed a prototype flat-panel VP-PET device and integrated it with a clinical PET/CT scanner. It significantly enhances the contrast of lesions, especially for those that are borderline detectable by the native scanner, within regions-of-interest specified by users. Simulation demonstrated the enhancement in lesion detectability with the VP-PET device. This technology may become a cost-effective solution for organ-specific imaging tasks.

Feasibility study of a point-of-care positron emission tomography system with interactive imaging capability.

PURPOSE: We investigated the feasibility of a novel positron emission tomography (PET) system that provides near real-time feedback to an operator who can interactively scan a patient to optimize image quality. The system should be compact and mobile to support point-of-care (POC) molecular imaging applications. In this study, we present the key technologies required and discuss the potential benefits of such new capability. METHODS: The core of this novel PET technology includes trackable PET detectors and a fully three-dimensional, fast image reconstruction engine implemented on multiple graphics processing units (GPUs) to support dynamically changing geometry by calculating the system matrix on-the-fly using a tube-of-response approach. With near real-time image reconstruction capability, a POC-PET system may comprise a maneuverable front PET detector and a second detector panel which can be stationary or moved synchronously with the front detector such that both panels face the region-of-interest (ROI) with the detector trajectory contoured around a patient's body. We built a proof-of-concept prototype using two planar detectors each consisting of a photomultiplier tube (PMT) optically coupled to an array of 48 × 48 lutetium-yttrium oxyorthosilicate (LYSO) crystals (1.0 × 1.0 × 10.0 mm3 each). Only 38 × 38 crystals in each arrays can be clearly re-solved and used for coincidence detection. One detector was mounted to a robotic arm which can position it at arbitrary locations, and the other detector was mounted on a rotational stage. A cylindrical phantom (102 mm in diameter, 150 mm long) with nine spherical lesions (8:1 tumor-to-background activity concentration ratio) was imaged from 27 sampling angles. List-mode events were reconstructed to form images without or with time-of-flight (TOF) information. We conducted two Monte Carlo simulations using two POC-PET systems. The first one uses the same phantom and detector setup as our experiment, with the detector coincidence re-solving time (CRT) ranging from 100 to 700 ps full-width-at-half-maximum (FWHM). The second study simulates a body-size phantom (316 × 228 × 160 mm3 ) imaged by a larger POC-PET system that has 4 × 6 modules (32 × 32 LYSO crystals/module, four in axial and six in transaxial directions) in the front panel and 3 × 8 modules (16 × 16 LYSO crystals/module, three in axial and eight in transaxial directions) in the back panel. We also evaluated an interactive scanning strategy by progressively increasing the number of data sets used for image reconstruction. The updated images were analyzed based on the number of data sets and the detector CRT. RESULTS: The proof-of-concept prototype re-solves most of the spherical lesions despite a limited number of coincidence events and incomplete sampling. TOF information reduces artifacts in the reconstructed images. Systems with better timing resolution exhibit improved image quality and reduced artifacts. We observed a reconstruction speed of 0.96 × 106 events/s/iteration for 600 × 600 × 224 voxel rectilinear space using four GPUs. A POC-PET system with significantly higher sensitivity can interactively image a body-size object from four angles in less than 7 min. CONCLUSIONS: We have developed GPU-based fast image reconstruction capability to support a PET system with arbitrary and dynamically changing geometry. Using TOF PET detectors, we demonstrated the feasibility of a PET system that can provide timely visual feedback to an operator who can scan a patient interactively to support POC imaging applications.

Low-dose photon counting CT reconstruction bias reduction with multi-energy alternating minimization algorithm.

Photon counting CT (PCCT) is an x-ray imaging technique that has undergone great development in the past decade. PCCT has the potential to improve dose efficiency and low-dose performance. In this paper, we propose a statistics-based iterative algorithm to perform a direct reconstruction of material-decomposed images. Compared with the conventional sinogram-based decomposition method which has degraded performance in low-dose scenarios, the multi-energy alternating minimization algorithm for photon counting CT (MEAM-PCCT) can generate accurate material-decomposed image with much smaller biases.

Experimental implementation of a joint statistical image reconstruction method for proton stopping power mapping from dual-energy CT data.

PURPOSE: To experimentally commission a dual-energy CT (DECT) joint statistical image reconstruction (JSIR) method, which is built on a linear basis vector model (BVM) of material characterization, for proton stopping power ratio (SPR) estimation. METHODS: The JSIR-BVM method builds on the relationship between the energy-dependent photon attenuation coefficients and the proton stopping power via a pair of BVM component weights. The two BVM component images are simultaneously reconstructed from the acquired DECT sinograms and then used to predict the electron density and mean excitation energy (I-value), which are required by the Bethe equation for SPR computation. A post-reconstruction image-based DECT method, which utilizes the two separate CT images reconstructed via the scanner's software, was implemented for comparison. The DECT measurement data were acquired on a Philips Brilliance scanner at 90 and 140 kVp for two phantoms of different sizes. Each phantom contains 12 different soft and bony tissue surrogates with known compositions. The SPR estimation results were compared to the reference values computed from the known compositions. The difference of the computed water equivalent path lengths (WEPL) across the phantoms between the two methods was also compared. RESULTS: The overall root-mean-square (RMS) of SPR estimation error of the JSIR-BVM method are 0.33% and 0.37% for the head- and body-sized phantoms, respectively, and all SPR estimates of the test samples are within 0.7% of the reference ground truth. The image-based method achieves overall RMS errors of 2.35% and 2.50% for the head- and body-sized phantoms, respectively. The JSIR-BVM method also reduces the pixel-wise random variation by 4-fold to 6-fold within homogeneous regions compared to the image-based method. The average differences between the JSIR-BVM method and the image-based method are 0.54% and 1.02% for the head- and body-sized phantoms, respectively. CONCLUSION: By taking advantage of an accurate polychromatic CT data model and a model-based DECT statistical reconstruction algorithm, the JSIR-BVM method accounts for both systematic bias and random noise in the acquired DECT measurement data. Therefore, the JSIR-BVM method achieves good accuracy and precision on proton SPR estimation for various tissue surrogates and object sizes. In contrast, the experimentally achievable accuracy of the image-based method may be limited by the uncertainties in the image formation process. The result suggests that the JSIR-BVM method has the potential for more accurate SPR prediction compared to post-reconstruction image-based methods in clinical settings.

Impact of joint statistical dual-energy CT reconstruction of proton stopping power images: Comparison to image- and sinogram-domain material decomposition approaches.

PURPOSE: The purpose of this study was to assess the performance of a novel dual-energy CT (DECT) approach for proton stopping power ratio (SPR) mapping that integrates image reconstruction and material characterization using a joint statistical image reconstruction (JSIR) method based on a linear basis vector model (BVM). A systematic comparison between the JSIR-BVM method and previously described DECT image- and sinogram-domain decomposition approaches is also carried out on synthetic data. METHODS: The JSIR-BVM method was implemented to estimate the electron densities and mean excitation energies (I-values) required by the Bethe equation for SPR mapping. In addition, image- and sinogram-domain DECT methods based on three available SPR models including BVM were implemented for comparison. The intrinsic SPR modeling accuracy of the three models was first validated. Synthetic DECT transmission sinograms of two 330 mm diameter phantoms each containing 17 soft and bony tissues (for a total of 34) of known composition were then generated with spectra of 90 and 140 kVp. The estimation accuracy of the reconstructed SPR images were evaluated for the seven investigated methods. The impact of phantom size and insert location on SPR estimation accuracy was also investigated. RESULTS: All three selected DECT-SPR models predict the SPR of all tissue types with less than 0.2% RMS errors under idealized conditions with no reconstruction uncertainties. When applied to synthetic sinograms, the JSIR-BVM method achieves the best performance with mean and RMS-average errors of less than 0.05% and 0.3%, respectively, for all noise levels, while the image- and sinogram-domain decomposition methods show increasing mean and RMS-average errors with increasing noise level. The JSIR-BVM method also reduces statistical SPR variation by sixfold compared to other methods. A 25% phantom diameter change causes up to 4% SPR differences for the image-domain decomposition approach, while the JSIR-BVM method and sinogram-domain decomposition methods are insensitive to size change. CONCLUSION: Among all the investigated methods, the JSIR-BVM method achieves the best performance for SPR estimation in our simulation phantom study. This novel method is robust with respect to sinogram noise and residual beam-hardening effects, yielding SPR estimation errors comparable to intrinsic BVM modeling error. In contrast, the achievable SPR estimation accuracy of the image- and sinogram-domain decomposition methods is dominated by the CT image intensity uncertainties introduced by the reconstruction and decomposition processes.

Multi-GPU Acceleration of Branchless Distance Driven Projection and Backprojection for Clinical Helical CT.

Model-based image reconstruction (MBIR) techniques have the potential to generate high quality images from noisy measurements and a small number of projections which can reduce the x-ray dose in patients. These MBIR techniques rely on projection and backprojection to refine an image estimate. One of the widely used projectors for these modern MBIR based technique is called branchless distance driven (DD) projection and backprojection. While this method produces superior quality images, the computational cost of iterative updates keeps it from being ubiquitous in clinical applications. In this paper, we provide several new parallelization ideas for concurrent execution of the DD projectors in multi-GPU systems using CUDA programming tools. We have introduced some novel schemes for dividing the projection data and image voxels over multiple GPUs to avoid runtime overhead and inter-device synchronization issues. We have also reduced the complexity of overlap calculation of the algorithm by eliminating the common projection plane and directly projecting the detector boundaries onto image voxel boundaries. To reduce the time required for calculating the overlap between the detector edges and image voxel boundaries, we have proposed a pre-accumulation technique to accumulate image intensities in perpendicular 2D image slabs (from a 3D image) before projection and after backprojection to ensure our DD kernels run faster in parallel GPU threads. For the implementation of our iterative MBIR technique we use a parallel multi-GPU version of the alternating minimization (AM) algorithm with penalized likelihood update. The time performance using our proposed reconstruction method with Siemens Sensation 16 patient scan data shows an average of 24 times speedup using a single TITAN X GPU and 74 times speedup using 3 TITAN X GPUs in parallel for combined projection and backprojection.

Technical Note: On the accuracy of parametric two-parameter photon cross-section models in dual-energy CT applications.

PURPOSE: To evaluate and compare the theoretically achievable accuracy of two families of two-parameter photon cross-section models: basis vector model (BVM) and modified parametric fit model (mPFM). METHOD: The modified PFM assumes that photoelectric absorption and scattering cross-sections can be accurately represented by power functions in effective atomic number and/or energy plus the Klein-Nishina cross-section, along with empirical corrections that enforce exact prediction of elemental cross-sections. Two mPFM variants were investigated: the widely used Torikoshi model (tPFM) and a more complex "VCU" variant (vPFM). For 43 standard soft and bony tissues and phantom materials, all consisting of elements with atomic number less than 20 (except iodine), we evaluated the theoretically achievable accuracy of tPFM and vPFM for predicting linear attenuation, photoelectric absorption, and energy-absorption coefficients, and we compared it to a previously investigated separable, linear two-parameter model, BVM. RESULTS: For an idealized dual-energy computed tomography (DECT) imaging scenario, the cross-section mapping process demonstrates that BVM more accurately predicts photon cross-sections of biological mixtures than either tPFM or vPFM. Maximum linear attenuation coefficient prediction errors were 15% and 5% for tPFM and BVM, respectively. The root-mean-square (RMS) prediction errors of total linear attenuation over the 20 keV to 1000 keV energy range of tPFM and BVM were 0.93% (tPFM) and 0.1% (BVM) for adipose tissue, 0.8% (tPFM) and 0.2% (BVM) for muscle tissue, and 1.6% (tPFM) and 0.2% (BVM) for cortical bone tissue. With exception of the thyroid and Teflon, the RMS error for photoelectric absorption and scattering coefficient was within 4% for the tPFM and 2% for the BVM. Neither model predicts the photon cross-sections of thyroid tissue accurately, exhibiting relative errors as large as 20%. For the energy-absorption coefficients prediction error, RMS errors for the BVM were less than 1.5%, while for the tPFM, the RMS errors were as large as 16%. CONCLUSION: Compared to modified PFMs, BVM shows superior potential to support dual-energy CT cross-section mapping. In addition, the linear, separable BVM can be more efficiently deployed by iterative model-based DECT image-reconstruction algorithms.

A simultaneous beta and coincidence-gamma imaging system for plant leaves.

Positron emitting isotopes, such as (11)C, (13)N, and (18)F, can be used to label molecules. The tracers, such as (11)CO2, are delivered to plants to study their biological processes, particularly metabolism and photosynthesis, which may contribute to the development of plants that have a higher yield of crops and biomass. Measurements and resulting images from PET scanners are not quantitative in young plant structures or in plant leaves due to poor positron annihilation in thin objects. To address this problem we have designed, assembled, modeled, and tested a nuclear imaging system (simultaneous beta-gamma imager). The imager can simultaneously detect positrons ([Formula: see text]) and coincidence-gamma rays (γ). The imaging system employs two planar detectors; one is a regular gamma detector which has a LYSO crystal array, and the other is a phoswich detector which has an additional BC-404 plastic scintillator for beta detection. A forward model for positrons is proposed along with a joint image reconstruction formulation to utilize the beta and coincidence-gamma measurements for estimating radioactivity distribution in plant leaves. The joint reconstruction algorithm first reconstructs beta and gamma images independently to estimate the thickness component of the beta forward model and afterward jointly estimates the radioactivity distribution in the object. We have validated the physics model and reconstruction framework through a phantom imaging study and imaging a tomato leaf that has absorbed (11)CO2. The results demonstrate that the simultaneously acquired beta and coincidence-gamma data, combined with our proposed joint reconstruction algorithm, improved the quantitative accuracy of estimating radioactivity distribution in thin objects such as leaves. We used the structural similarity (SSIM) index for comparing the leaf images from the simultaneous beta-gamma imager with the ground truth image. The jointly reconstructed images yield SSIM indices of 0.69 and 0.63, whereas the separately reconstructed beta alone and gamma alone images had indices of 0.33 and 0.52, respectively.

Cardiorespiratory interactions: Noncontact assessment using laser Doppler vibrometry.

The application of a noncontact physiological recording technique, based on the method of laser Doppler vibrometry (LDV), is described. The effectiveness of the LDV method as a physiological recording modality lies in the ability to detect very small movements of the skin, associated with internal mechanophysiological activities. The method is validated for a range of cardiovascular variables, extracted from the contour of the carotid pulse waveform as a function of phase of the respiration cycle. Data were obtained from 32 young healthy participants, while resting and breathing spontaneously. Individual beats were assigned to four segments, corresponding with inspiration and expiration peaks and transitional periods. Measures relating to cardiac and vascular dynamics are shown to agree with the pattern of effects seen in the substantial body of literature based on human and animal experiments, and with selected signals recorded simultaneously with conventional sensors. These effects include changes in heart rate, systolic time intervals, and stroke volume. There was also some evidence for vascular adjustments over the respiration cycle. The effectiveness of custom algorithmic approaches for extracting the key signal features was confirmed. The advantages of the LDV method are discussed in terms of the metrological properties and utility in psychophysiological research. Although used here within a suite of conventional sensors and electrodes, the LDV method can be used on a stand-alone, noncontact basis, with no requirement for skin preparation, and can be used in harsh environments including the MR scanner.

Line Integral Alternating Minimization Algorithm for Dual-Energy X-Ray CT Image Reconstruction.

We propose a new algorithm, called line integral alternating minimization (LIAM), for dual-energy X-ray CT image reconstruction. Instead of obtaining component images by minimizing the discrepancy between the data and the mean estimates, LIAM allows for a tunable discrepancy between the basis material projections and the basis sinograms. A parameter is introduced that controls the size of this discrepancy, and with this parameter the new algorithm can continuously go from a two-step approach to the joint estimation approach. LIAM alternates between iteratively updating the line integrals of the component images and reconstruction of the component images using an image iterative deblurring algorithm. An edge-preserving penalty function can be incorporated in the iterative deblurring step to decrease the roughness in component images. Images from both simulated and experimentally acquired sinograms from a clinical scanner were reconstructed by LIAM while varying the regularization parameters to identify good choices. The results from the dual-energy alternating minimization algorithm applied to the same data were used for comparison. Using a small fraction of the computation time of dual-energy alternating minimization, LIAM achieves better accuracy of the component images in the presence of Poisson noise for simulated data reconstruction and achieves the same level of accuracy for real data reconstruction.

A generalized reconstruction framework for unconventional PET systems.

PURPOSE: Quantitative estimation of the radionuclide activity concentration in positron emission tomography (PET) requires precise modeling of PET physics. The authors are focused on designing unconventional PET geometries for specific applications. This work reports the creation of a generalized reconstruction framework, capable of reconstructing tomographic PET data for systems that use right cuboidal detector elements positioned at arbitrary geometry using a regular Cartesian grid of image voxels. METHODS: The authors report on a variety of design choices and optimization for the creation of the generalized framework. The image reconstruction algorithm is maximum likelihood-expectation-maximization. System geometry can be specified using a simple script. Given the geometry, a symmetry seeking algorithm finds existing symmetry in the geometry with respect to the image grid to improve the memory usage/speed. Normalization is approached from a geometry independent perspective. The system matrix is computed using the Siddon's algorithm and subcrystal approach. The program is parallelized through open multiprocessing and message passing interface libraries. A wide variety of systems can be modeled using the framework. This is made possible by modeling the underlying physics and data correction, while generalizing the geometry dependent features. RESULTS: Application of the framework for three novel PET systems, each designed for a specific application, is presented to demonstrate the robustness of the framework in modeling PET systems of unconventional geometry. Three PET systems of unconventional geometry are studied. (1) Virtual-pinhole half-ring insert integrated into Biograph-40: although the insert device improves image quality over conventional whole-body scanner, the image quality varies depending on the position of the insert and the object. (2) Virtual-pinhole flat-panel insert integrated into Biograph-40: preliminary results from an investigation into a modular flat-panel insert are presented. (3) Plant PET system: a reconfigurable PET system for imaging plants, with resolution of greater than 3.3 mm, is shown. Using the automated symmetry seeking algorithm, the authors achieved a compression ratio of the storage and memory requirement by a factor of approximately 50 for the half-ring and flat-panel systems. For plant PET system, the compression ratio is approximately five. The ratio depends on the level of symmetry that exists in different geometries. CONCLUSIONS: This work brings the field closer to arbitrary geometry reconstruction. A generalized reconstruction framework can be used to validate multiple hypotheses and the effort required to investigate each system is reduced. Memory usage/speed can be improved with certain optimizations.

A dedicated high-resolution PET imager for plant sciences.

PET provides an in vivo molecular and functional imaging capability that could be valuable for studying the interaction of plants in changing environments at the whole-plant level. We have developed a dedicated plant PET imager housed in a plant growth chamber (PGC), which provides a fully controlled environment. The system currently contains two types of scintillation detector modules from commercial small animal PET scanners: 84 microPET® detectors, which are made with scintillation crystal arrays of 2.2 mm(3) × 2.2 mm(3) × 10 mm(3) crystals to provide a large detection area; and 32 Inveon™ detectors, which are made with scintillation crystal arrays of 1.5 mm(3) × 1.5 mm(3) × 10 mm(3) crystals to provide higher spatial resolution. The detector modules are configured to form two half-rings, which provide a 15 cm-diameter trans-axial field of view (FOV) for dynamic tomographic imaging of small plants. Alternatively, the Inveon detectors can be reconfigured to form quarter-rings, which provide a 25 cm FOV using step-and-shoot motion. The imager contains two linear stages that move detectors vertically at different heights for multisection scanning, and two rotation stages to collect coincidence events from all angles when using the step-and-shoot acquisition. The detector modules and mechanical components of the imager are housed inside a PGC that regulates the environmental parameters. The system has a typical energy resolution of 15% for the Inveon detectors and 24% for the microPET detectors, timing resolution of 1.8 ns, and sensitivity of 1.3%, 1.4% and 3.0% measured at the center of the FOV, 5 cm off to the larger half-ring and 5 cm off to the smaller half-ring, respectively (with a 350-650 keV energy window and 3.1 ns timing window). The system's spatial resolution is capable of resolving rod sources of 1.25 mm diameter spaced 2.5 mm apart (center to center) using the ML-EM reconstruction algorithm. Preliminary imaging experiments using soybean and wild type and mutant maize labeled with (11)CO2 produced high-quality dynamic PET images that reveal the translocation and distribution patterns of photoassimilates. This system can be used to provide an in vivo molecular and functional imaging capability for plant research.

Compressed sampling strategies for tomography.

We investigate new sampling strategies for projection tomography, enabling one to employ fewer measurements than expected from classical sampling theory without significant loss of information. Inspired by compressed sensing, our approach is based on the understanding that many real objects are compressible in some known representation, implying that the number of degrees of freedom defining an object is often much smaller than the number of pixels/voxels. We propose a new approach based on quasi-random detector subsampling, whereas previous approaches only addressed subsampling with respect to source location (view angle). The performance of different sampling strategies is considered using object-independent figures of merit, and also based on reconstructions for specific objects, with synthetic and real data. The proposed approach can be implemented using a structured illumination of the interrogated object or the detector array by placing a coded aperture/mask at the source or detector side, respectively. Advantages of the proposed approach include (i) for structured illumination of the detector array, it leads to fewer detector pixels and allows one to integrate detectors for scattered radiation in the unused space; (ii) for structured illumination of the object, it leads to a reduced radiation dose for patients in medical scans; (iii) in the latter case, the blocking of rays reduces scattered radiation while keeping the same energy in the transmitted rays, resulting in a higher signal-to-noise ratio than that achieved by lowering exposure times or the energy of the source; (iv) compared to view-angle subsampling, it allows one to use fewer measurements for the same image quality, or leads to better image quality for the same number of measurements. The proposed approach can also be combined with view-angle subsampling.

Prospects for in vivo estimation of photon linear attenuation coefficients using postprocessing dual-energy CT imaging on a commercial scanner: comparison of analytic and polyenergetic statistical reconstruction algorithms.

PURPOSE: Accurate patient-specific photon cross-section information is needed to support more accurate model-based dose calculation for low energy photon-emitting modalities in medicine such as brachytherapy and kilovoltage x-ray imaging procedures. A postprocessing dual-energy CT (pDECT) technique for noninvasive in vivo estimation of photon linear attenuation coefficients has been experimentally implemented on a commercial CT scanner and its accuracy assessed in idealized phantom geometries. METHODS: Eight test materials of known composition and density were used to compare pDECT-estimated linear attenuation coefficients to NIST reference values over an energy range from 10 keV to 1 MeV. As statistical image reconstruction (SIR) has been shown to reconstruct images with less random and systematic error than conventional filtered backprojection (FBP), the pDECT technique was implemented with both an in-house polyenergetic SIR algorithm, alternating minimization (AM), as well as a conventional FBP reconstruction algorithm. Improvement from increased spectral separation was also investigated by filtering the high-energy beam with an additional 0.5 mm of tin. The law of propagated uncertainty was employed to assess the sensitivity of the pDECT process to errors in reconstructed images. RESULTS: Mean pDECT-estimated linear attenuation coefficients for the eight test materials agreed within 1% of NIST reference values for energies from 1 MeV down to 30 keV, with mean errors rising to between 3% and 6% at 10 keV, indicating that the method is unbiased when measurement and calibration phantom geometries are matched. Reconstruction with FBP and AM algorithms conferred similar mean pDECT accuracy. However, single-voxel pDECT estimates reconstructed on a 1 × 1 × 3 mm(3) grid are shown to be highly sensitive to reconstructed image uncertainty; in some cases pDECT attenuation coefficient estimates exhibited standard deviations on the order of 20% around the mean. Reconstruction with the statistical AM algorithm led to standard deviations roughly 40% to 60% less than FBP reconstruction. Additional tin filtration of the high energy beam exhibits similar pDECT estimation accuracy as the unfiltered beam, even when scanning with only 25% of the dose. Using the law of propagated uncertainty, low Z materials are found to be more sensitive to image reconstruction errors than high Z materials. Furthermore, it is estimated that reconstructed CT image uncertainty must be limited to less than 0.25% to achieve a target linear-attenuation coefficient estimation uncertainty of 3% at 28 keV. CONCLUSIONS: That pDECT supports mean linear attenuation coefficient measurement accuracies of 1% of reference values for energies greater than 30 keV is encouraging. However, the sensitivity of the pDECT measurements to noise and systematic errors in reconstructed CT images warrants further investigation in more complex phantom geometries. The investigated statistical reconstruction algorithm, AM, reduced random measurement uncertainty relative to FBP owing to improved noise performance. These early results also support efforts to increase DE spectral separation, which can further reduce the pDECT sensitivity to measurement uncertainty.