Jae H Lee, Yushu Yao, Uttam Shrestha, Grant T Gullberg, Youngho Seo
{"title":"Handling Big Data in Medical Imaging: Iterative Reconstruction with Large-Scale Automated Parallel Computation.","authors":"Jae H Lee, Yushu Yao, Uttam Shrestha, Grant T Gullberg, Youngho Seo","doi":"10.1109/NSSMIC.2014.7430758","DOIUrl":"https://doi.org/10.1109/NSSMIC.2014.7430758","url":null,"abstract":"<p><p>The primary goal of this project is to implement the iterative statistical image reconstruction algorithm, in this case maximum likelihood expectation maximum (MLEM) used for dynamic cardiac single photon emission computed tomography, on Spark/GraphX. This involves porting the algorithm to run on large-scale parallel computing systems. Spark is an easy-to- program software platform that can handle large amounts of data in parallel. GraphX is a graph analytic system running on top of Spark to handle graph and sparse linear algebra operations in parallel. The main advantage of implementing MLEM algorithm in Spark/GraphX is that it allows users to parallelize such computation without any expertise in parallel computing or prior knowledge in computer science. In this paper we demonstrate a successful implementation of MLEM in Spark/GraphX and present the performance gains with the goal to eventually make it useable in clinical setting.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2014 ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2014-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2014.7430758","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34463545","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Uttam Shrestha, Elias H Botvinick, Yerem Yeghiazarians, Youngho Seo, Grant T Gullberg
{"title":"Quantitative Signature of <i>Coronary Steal</i> in a Patient with Occluded Coronary Arteries Supported by Collateral Circulation Using Dynamic SPECT.","authors":"Uttam Shrestha, Elias H Botvinick, Yerem Yeghiazarians, Youngho Seo, Grant T Gullberg","doi":"10.1109/NSSMIC.2014.7430923","DOIUrl":"https://doi.org/10.1109/NSSMIC.2014.7430923","url":null,"abstract":"<p><p><i>Coronary steal</i> (CS) is a physiological process that induces absolute decrease in blood flow in collateralized myocardium compared to resting flow during coronary vasodilation due to redistribution of blood away from collateral-dependent myocardium. Although, CS has been well known for decades, there are very few noninvasive perfusion studies in humans that quantitatively predict the existence of CS. In this study, we show that the quantitative measurement of absolute value of regional myocardial blood flow (MBF) and coronary flow reserve (CFR) using dynamic single photon emitted computed tomography (SPECT) can help estimate the presence of CS in myocardium with obstructed coronary artery and collateral circulation.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2014 ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2014-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2014.7430923","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34463547","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Shih-Ying Huang, Dragana Savic, Jaewon Yang, Uttam Shrestha, Youngho Seo
{"title":"The Effect of Magnetic Field on Positron Range and Spatial Resolution in an Integrated Whole-Body Time-Of-Flight PET/MRI System.","authors":"Shih-Ying Huang, Dragana Savic, Jaewon Yang, Uttam Shrestha, Youngho Seo","doi":"10.1109/NSSMIC.2014.7431006","DOIUrl":"https://doi.org/10.1109/NSSMIC.2014.7431006","url":null,"abstract":"<p><p>Simultaneous imaging systems combining positron emission tomography (PET) and magnetic resonance imaging (MRI) have been actively investigated. A PET/MR imaging system (GE Healthcare) comprised of a time-of-flight (TOF) PET system utilizing silicon photomultipliers (SiPMs) and 3-tesla (3T) MRI was recently installed at our institution. The small-ring (60 cm diameter) TOF PET subsystem of this PET/MRI system can generate images with higher spatial resolution compared with conventional PET systems. We have examined theoretically and experimentally the effect of uniform magnetic fields on the spatial resolution for high-energy positron emitters. Positron emitters including <sup>18</sup>F, <sup>124</sup>I, and <sup>68</sup>Ga were simulated in water using the Geant4 Monte Carlo toolkit in the presence of a uniform magnetic field (0, 3, and 7 Tesla). The positron annihilation position was tracked to determine the 3D spatial distribution of the 511-keV gammy ray emission. The full-width at tenth maximum (FWTM) of the positron point spread function (PSF) was determined. Experimentally, <sup>18</sup>F and <sup>68</sup>Ga line source phantoms in air and water were imaged with an investigational PET/MRI system and a PET/CT system to investigate the effect of magnetic field on the spatial resolution of PET. The full-width half maximum (FWHM) of the line spread function (LSF) from the line source was determined as the system spatial resolution. Simulations and experimental results show that the in-plane spatial resolution was slightly improved at field strength as low as 3 Tesla, especially when resolving signal from high-energy positron emitters in the air-tissue boundary.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2014 ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2014-11-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2014.7431006","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34460225","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Luca Caucci, Abhinav K Jha, Lars R Furenlid, Eric W Clarkson, Matthew A Kupinski, Harrison H Barrett
{"title":"Image Science with Photon-Processing Detectors.","authors":"Luca Caucci, Abhinav K Jha, Lars R Furenlid, Eric W Clarkson, Matthew A Kupinski, Harrison H Barrett","doi":"10.1109/NSSMIC.2013.6829331","DOIUrl":"https://doi.org/10.1109/NSSMIC.2013.6829331","url":null,"abstract":"<p><p>We introduce and discuss photon-processing detectors and we compare them with photon-counting detectors. By estimating a relatively small number of attributes for each collected photon, photon-processing detectors may help understand and solve a fundamental theoretical problem of any imaging system based on photon-counting detectors, namely null functions. We argue that photon-processing detectors can improve task performance by estimating position, energy, and time of arrival for each collected photon. We consider a continuous-to-continuous linear operator to relate the object being imaged to the collected data, and discuss how this operator can be analyzed to derive properties of the imaging system. Finally, we derive an expression for the characteristic functional of an imaging system that produces list-mode data.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2013 ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2013-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2013.6829331","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"9871818","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Rahul Khare, Guillaume Sala, Paul Kinahan, Giuseppe Esposito, Filip Banovac, Kevin Cleary, Andinet Enquobahrie
{"title":"Experimental Evaluation of a Deformable Registration Algorithm for Motion Correction in PET-CT Guided Biopsy.","authors":"Rahul Khare, Guillaume Sala, Paul Kinahan, Giuseppe Esposito, Filip Banovac, Kevin Cleary, Andinet Enquobahrie","doi":"10.1109/NSSMIC.2013.6829037","DOIUrl":"https://doi.org/10.1109/NSSMIC.2013.6829037","url":null,"abstract":"<p><p>Positron emission tomography computed tomography (PET-CT) images are increasingly being used for guidance during percutaneous biopsy. However, due to the physics of image acquisition, PET-CT images are susceptible to problems due to respiratory and cardiac motion, leading to inaccurate tumor localization, shape distortion, and attenuation correction. To address these problems, we present a method for motion correction that relies on respiratory gated CT images aligned using a deformable registration algorithm. In this work, we use two deformable registration algorithms and two optimization approaches for registering the CT images obtained over the respiratory cycle. The two algorithms are the BSpline and the symmetric forces Demons registration. In the first optmization approach, CT images at each time point are registered to a single reference time point. In the second approach, deformation maps are obtained to align each CT time point with its adjacent time point. These deformations are then composed to find the deformation with respect to a reference time point. We evaluate these two algorithms and optimization approaches using respiratory gated CT images obtained from 7 patients. Our results show that overall the BSpline registration algorithm with the reference optimization approach gives the best results.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2013-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2013.6829037","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"33085117","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
William C J Hunter, Robert S Miyaoka, Lawrence MacDonald, Wendy McDougald, Thomas K Lewellen
{"title":"Light-Sharing Interface for dMiCE Detectors using Sub-Surface Laser Engraving.","authors":"William C J Hunter, Robert S Miyaoka, Lawrence MacDonald, Wendy McDougald, Thomas K Lewellen","doi":"10.1109/NSSMIC.2013.6829407","DOIUrl":"https://doi.org/10.1109/NSSMIC.2013.6829407","url":null,"abstract":"<p><p>We have previously reported on dMiCE, a method of resolving depth or interaction (DOI) in a pair of discrete crystals by encoding light sharing properties as a function of depth in the interface of this crystal-element pair. A challenge for this method is the cost and repeatability of interface treatment for a crystal pair. In this work, we report our preliminary results on using sub-surface laser engraving (SSLE) as a means of forming this depth-dependent interface in a dMiCE detector. A surplus first-generation SSLE system was used to create a partially reflective layer 100-microns thick at the boundary between two halves of a 1.4-by-2.9-by-20 mmˆ3 LYSO crystal. The boundary of these paired crystal elements was positioned between two 3-mm wide Geiger-Müller avalanche photodiodes from Hamamatsu. The responses of these two photodetectors were acquired for an ensemble of 511-keV photons collimated to interact at a fixed depth in just one crystal element. Interaction position was then varied to measure detector response as a function of depth, which was then used to maximum-likelihood positions events. Despite use of sub-optimal SSLE processing we found an average DOI resolution of 3.4 mm for front-sided readout and 3.9 mm for back-sided readout. We expect DOI resolution can be improved significantly by optimizing the SSLE process and pattern.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2013 ","pages":"1-7"},"PeriodicalIF":0.0,"publicationDate":"2013-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2013.6829407","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"32908906","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Fenghua Weng, Srijeeta Bagchi, Qiu Huang, Youngho Seo
{"title":"Design Studies of a CZT-based Detector Combined with a Pixel-Geometry-Matching Collimator for SPECT Imaging.","authors":"Fenghua Weng, Srijeeta Bagchi, Qiu Huang, Youngho Seo","doi":"10.1109/NSSMIC.2013.6829458","DOIUrl":"https://doi.org/10.1109/NSSMIC.2013.6829458","url":null,"abstract":"<p><p>Single Photon Emission Computed Tomography (SPECT) suffers limited efficiency due to the need for collimators. Collimator properties largely decide the data statistics and image quality. Various materials and configurations of collimators have been investigated in many years. The main thrust of our study is to evaluate the design of pixel-geometry-matching collimators to investigate their potential performances using Geant4 Monte Carlo simulations. Here, a pixel-geometry-matching collimator is defined as a collimator which is divided into the same number of pixels as the detector's and the center of each pixel in the collimator is a one-to-one correspondence to that in the detector. The detector is made of Cadmium Zinc Telluride (CZT), which is one of the most promising materials for applications to detect hard X-rays and <i>γ</i>-rays due to its ability to obtain good energy resolution and high light output at room temperature. For our current project, we have designed a large-area, CZT-based gamma camera (20.192 cm×20.192 cm) with a small pixel pitch (1.60 mm). The detector is pixelated and hence the intrinsic resolution can be as small as the size of the pixel. Materials of collimator, collimator hole geometry, detection efficiency, and spatial resolution of the CZT detector combined with the pixel-matching collimator were calculated and analyzed under different conditions. From the simulation studies, we found that such a camera using rectangular holes has promising imaging characteristics in terms of spatial resolution, detection efficiency, and energy resolution.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2013 ","pages":"1-4"},"PeriodicalIF":0.0,"publicationDate":"2013-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2013.6829458","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"32800457","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Soo Mee Kim, Adam M Alessio, Bruno De Man, Evren Asma, Paul E Kinahan
{"title":"Direct Reconstruction of CT-based Attenuation Correction Images for PET with Cluster-Based Penalties.","authors":"Soo Mee Kim, Adam M Alessio, Bruno De Man, Evren Asma, Paul E Kinahan","doi":"10.1109/NSSMIC.2013.6829245","DOIUrl":"https://doi.org/10.1109/NSSMIC.2013.6829245","url":null,"abstract":"<p><p>Extremely low-dose CT acquisitions for the purpose of PET attenuation correction will have a high level of noise and biasing artifacts due to factors such as photon starvation. This work explores a priori knowledge appropriate for CT iterative image reconstruction for PET attenuation correction. We investigate the maximum a posteriori (MAP) framework with cluster-based, multinomial priors for the direct reconstruction of the PET attenuation map. The objective function for direct iterative attenuation map reconstruction was modeled as a Poisson log-likelihood with prior terms consisting of quadratic (Q) and mixture (M) distributions. The attenuation map is assumed to have values in 4 clusters: air+background, lung, soft tissue, and bone. Under this assumption, the MP was a mixture probability density function consisting of one exponential and three Gaussian distributions. The relative proportion of each cluster was jointly estimated during each voxel update of direct iterative coordinate decent (dICD) method. Noise-free data were generated from NCAT phantom and Poisson noise was added. Reconstruction with FBP (ramp filter) was performed on the noise-free (ground truth) and noisy data. For the noisy data, dICD reconstruction was performed with the combination of different prior strength parameters (β and γ) of Q- and M-penalties. The combined quadratic and mixture penalties reduces the RMSE by 18.7% compared to post-smoothed iterative reconstruction and only 0.7% compared to quadratic alone. For direct PET attenuation map reconstruction from ultra-low dose CT acquisitions, the combination of quadratic and mixture priors offers regularization of both variance and bias and is a potential method to derive attenuation maps with negligible patient dose. However, the small improvement in quantitative accuracy relative to the substantial increase in algorithm complexity does not currently justify the use of mixture-based PET attenuation priors for reconstruction of CT images for PET attenuation correction.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":" ","pages":""},"PeriodicalIF":0.0,"publicationDate":"2013-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2013.6829245","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"34017309","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ryeojin Park, Brian W Miller, Abhinav K Jha, Lars R Furenlid, William C J Hunter, Harrison H Barrett
{"title":"A Prototype Detector for a Novel High-Resolution PET System: BazookaPET.","authors":"Ryeojin Park, Brian W Miller, Abhinav K Jha, Lars R Furenlid, William C J Hunter, Harrison H Barrett","doi":"10.1109/NSSMIC.2012.6551486","DOIUrl":"https://doi.org/10.1109/NSSMIC.2012.6551486","url":null,"abstract":"<p><p>We have designed and are developing a novel proof-of-concept PET system called BazookaPET. In order to complete the PET configuration, at least two detector elements are required to detect positron-electron annihilation events. Each detector element of the BazookaPET has two independent data acquisition channels. One side of the scintillation crystal is optically coupled to a 4×4 silicon photomultiplier (SiPM) array and the other side is a CCD-based gamma camera. Using these two separate channels, we can obtain data with high energy, temporal and spatial resolution data by associating the data outputs via several maximum-likelihood estimation (MLE) steps. In this work, we present the concept of the system and the prototype detector element. We focus on characterizing individual detector channels, and initial experimental calibration results are shown along with preliminary performance-evaluation results. We measured energy resolution and the integrated traces of the slit-beam images from both detector channel outputs. A photo-peak energy resolution of ~5.3% FWHM was obtained from the SiPM and ~48% FWHM from the CCD at 662 keV. We assumed SiPM signals follow Gaussian statistics and estimated the 2D interaction position using MLE. Based on our the calibration experiments, we computed the Cramér-Rao bound (CRB) for the SiPM detector channel and found that the CRB resolution is better than 1 mm in the center of the crystal.</p>","PeriodicalId":73298,"journal":{"name":"IEEE Nuclear Science Symposium conference record. Nuclear Science Symposium","volume":"2012 ","pages":"2123-2127"},"PeriodicalIF":0.0,"publicationDate":"2012-10-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://sci-hub-pdf.com/10.1109/NSSMIC.2012.6551486","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"33959564","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}