Multimodality assessment of the coronary microvasculature with TIMI frame count versus perfusion PET highlights coronary changes characteristic of coronary microvascular disease.

Background: The diagnosis of coronary microvascular disease (CMVD) remains challenging. Perfusion PET-derived myocardial blood flow (MBF) reserve (MBFR) can quantify CMVD but is not widely available. Thrombolysis in Myocardial Infarction (TIMI) frame count (TFC) is an angiography-based method that has been proposed as a measure of CMVD. Here, we compare TFC and PET-derived MBF measurements to establish the role of TFC in assessing for CMVD. We use coronary modeling to elucidate the relationship between MBFR and TFC and propose TFC thresholds for identifying CMVD. Methods: In a cohort of 123 individuals (age 58 ± 12.1, 63% women, 41% Caucasian) without obstructive coronary artery disease who had undergone perfusion PET and coronary angiography for clinical indications, we compared TFC and perfusion PET parameters using Pearson correlation (PCC) and linear regression modeling. We used mathematical modeling of the coronary circulation to understand the relationship between these parameters and performed Receiver Operating Curve (ROC) analysis. Results: We found a significant negative correlation between TFC and MBFR. Sex, race and ethnicity, and nitroglycerin administration impact this relationship. Coronary modeling showed an uncoupling between TFC and flow in epicardial vessels. In ROC analysis, TFC performed well in women (AUC 0.84-0.89) and a moderately in men (AUC 0.68-0.78). Conclusions: We established an inverse relationship between TFC and PET-derived MBFR, which is affected by patient selection and procedural factors. TFC represents a measure of the volume of the epicardial coronary compartment, which is increased in patients with CMVD, and performs well in identifying women with CMVD.

Improving Generalizability of PET DL Algorithms: List-Mode Reconstructions Improve DOTATATE PET Hepatic Lesion Detection Performance.

Deep learning (DL) algorithms used for DOTATATE PET lesion detection typically require large, well-annotated training datasets. These are difficult to obtain due to low incidence of gastroenteropancreatic neuroendocrine tumors (GEP-NETs) and the high cost of manual annotation. Furthermore, networks trained and tested with data acquired from site specific PET/CT instrumentation, acquisition and processing protocols have reduced performance when tested with offsite data. This lack of generalizability requires even larger, more diverse training datasets. The objective of this study is to investigate the feasibility of improving DL algorithm performance by better matching the background noise in training datasets to higher noise, out-of-domain testing datasets. 68Ga-DOTATATE PET/CT datasets were obtained from two scanners: Scanner1, a state-of-the-art digital PET/CT (GE DMI PET/CT; n = 83 subjects), and Scanner2, an older-generation analog PET/CT (GE STE; n = 123 subjects). Set1, the data set from Scanner1, was reconstructed with standard clinical parameters (5 min; Q.Clear) and list-mode reconstructions (VPFXS 2, 3, 4, and 5-min). Set2, data from Scanner2 representing out-of-domain clinical scans, used standard iterative reconstruction (5 min; OSEM). A deep neural network was trained with each dataset: Network1 for Scanner1 and Network2 for Scanner2. DL performance (Network1) was tested with out-of-domain test data (Set2). To evaluate the effect of training sample size, we tested DL model performance using a fraction (25%, 50% and 75%) of Set1 for training. Scanner1, list-mode 2-min reconstructed data demonstrated the most similar noise level compared that of Set2, resulting in the best performance (F1 = 0.713). This was not significantly different compared to the highest performance, upper-bound limit using in-domain training for Network2 (F1 = 0.755; p-value = 0.103). Regarding sample size, the F1 score significantly increased from 25% training data (F1 = 0.478) to 100% training data (F1 = 0.713; p < 0.001). List-mode data from modern PET scanners can be reconstructed to better match the noise properties of older scanners. Using existing data and their associated annotations dramatically reduces the cost and effort in generating these datasets and significantly improves the performance of existing DL algorithms. List-mode reconstructions can provide an efficient, low-cost method to improve DL algorithm generalizability.

Performance evaluation of resolution and sensitivity of C-SPECT's variable slat-stack collimator.

BACKGROUND: Myocardial perfusion imaging is commonly performed using SPECT, where both general-purpose and dedicated scanners are available. A limitation with general-purpose systems has been the inability to image dynamically since different projections are obtained far apart in time due to scanner rotation. Dedicated systems can have this capability since they acquire completely sampled projections (i.e., those with enough angular views for reconstruction) with short time frames. C-SPECT, does not need any scanner or patient motion to obtain complete projections, allowing fast dynamics. When imaging fast dynamics, the optimal collimator settings are not necessarily the same as for static imaging, where longer acquisitions can be utilized. Thus, C-SPECT offers adaptive collimation in the transverse and axial directions. PURPOSE: The performance of adaptation in the axial direction was characterized herein. METHODS: The ratio of the resolution metric in high-sensitivity mode to that in the high-resolution mode, termed resolution boost factor, was determined. Analogously, the sensitivity boost factor was also determined. Comparisons were made with theory and simulations. RESULTS: The boost factors for resolution and sensitivity, averaged over the 14 modules of the system, were determined to be 1.72 and 1.75, respectively. CONCLUSIONS: The boost factors, which ideally would be two, were between 10% and 15% below optimal values and tracked each other, suggesting mechanical challenges in the apparatus, such as incomplete closure of adjacent slats, but show reasonably successful adaptation between modes.

Modeling contrast-to-noise ratio from list mode reconstructions of 68Ga DOTATATE PET/CT: predicting detectability of hepatic metastases in shorter acquisition PET reconstructions.

BACKGROUND: Deep learning (DL) algorithms have shown promise in identifying and quantifying lesions in PET/CT. However, the accuracy and generalizability of these algorithms relies on large, diverse datasets which are time and labor intensive to curate. Modern PET/CT scanners may acquire data in list mode, allowing for multiple reconstructions of the same datasets with different parameters and imaging times. These reconstructions may provide a wide range of image characteristics to increase the size and diversity of datasets. Training algorithms with shorter imaging times and higher noise properties requires that lesions remain detectable. The purpose of this study is to model and predict the contrast-to-noise ratio (CNR) for shorter imaging times based on CNR from longer duration, lower noise images for 68Ga DOTATATE PET hepatic lesions and identify a threshold above which lesions remain detectable. METHODS: 68Ga DOTATATE subjects (n=20) with hepatic lesions were divided into two subgroups. The "Model" group (n=4 subjects; n=9 lesions; n=36 datapoints) was used to identify the relationship between CNR and imaging time. The "Test" group (n=16 subjects; n=44 lesions; n=176 datapoints) was used to evaluate the prediction provided by the model. RESULTS: CNR plotted as a function of imaging time for a subset of identified subjects was very well fit with a quadratic model. For the remaining subjects, the measured CNR showed a very high linear correlation with the predicted CNR for these lesions (R2 > 0.97) for all imaging durations. From the model, a threshold CNR=6.9 at 5-minutes predicted CNR > 5 at 2-minutes. Visual inspection of lesions in 2-minute images with CNR above the threshold in 5-minute images were assessed and rated as a 4 or 5 (probably positive or definitely positive) confirming 100% lesion detectability on the shorter 2-minute PET images. CONCLUSIONS: CNR for shorter DOTATATE PET imaging times may be accurately predicted using list mode reconstructions of longer acquisitions. A threshold CNR may be applied to longer duration images to ensure lesion detectability of shorter duration reconstructions. This method can aid in the selection of lesions to include in novel data augmentation techniques for deep learning.

A Novel Mouse Model of Radiation-Induced Cardiac Injury Reveals Biological and Radiological Biomarkers of Cardiac Dysfunction with Potential Clinical Relevance.

PURPOSE: Radiation-induced cardiotoxicity is a significant concern in thoracic oncology patients. However, the basis for this disease pathology is not well characterized. We developed a novel mouse model of radiation-induced cardiotoxicity to investigate pathophysiologic mechanisms and identify clinically targetable biomarkers of cardiac injury. EXPERIMENTAL DESIGN: Single radiation doses of 20, 40, or 60 Gy were delivered to the cardiac apex of female C57BL/6 mice ages 9-11 weeks, with or without adjacent lung tissue, using conformal radiotherapy. Cardiac tissue was harvested up to 24 weeks post-radiotherapy for histologic analysis. Echocardiography and Technetium-99m sestamibi single photon emission computed tomography (SPECT) at 8 and 16 weeks post-radiotherapy were implemented to evaluate myocardial function and perfusion. Mouse cardiac tissue and mouse and human plasma were harvested for biochemical studies. RESULTS: Histopathologically, radiotherapy resulted in perivascular fibrosis 8 and 24 (P < 0.05) weeks post-radiotherapy. Apical perfusion deficits on SPECT and systolic and diastolic dysfunction on echocardiography 8 and 16 weeks post-radiotherapy were also observed (P < 0.05). Irradiated cardiac tissue and plasma showed significant increases in placental growth factor (PlGF), IL6, and TNFα compared with nonradiated matched controls, with greater increases in cardiac cytokine levels when radiotherapy involved lung. Human plasma showed increased PlGF (P = 0.021) and TNFα (P = 0.036) levels after thoracic radiotherapy. PlGF levels demonstrated a strong correlation (r = 0.89, P = 0.0001) with mean heart dose. CONCLUSIONS: We developed and characterized a pathophysiologically relevant mouse model of radiation-induced cardiotoxicity involving in situ irradiation of the cardiac apex. The model can be used to integrate radiomic and biochemical markers of cardiotoxicity to inform early therapeutic intervention and human translational studies.

Diagnostic accuracy of SPECT and PET myocardial perfusion imaging in patients with left bundle branch block or ventricular-paced rhythm.

BACKGROUND: The difference in diagnostic accuracy of coronary artery disease (CAD) between vasodilator SPECT and PET myocardial perfusion imaging (MPI) in patients with left bundle branch block (LBBB) or ventricular-paced rhythm (VPR) is unknown. METHODS: We identified patients with LBBB or VPR who underwent either vasodilator SPECT or PET MPI and subsequent coronary angiography. LBBB/VPR-related septal and anteroseptal defects were defined as perfusion defects involving those regions in the absence of obstructive CAD in the left anterior descending artery or left main coronary artery. RESULTS: Of the 55 patients who underwent coronary angiography, 38 (69%) underwent SPECT and 17 patients (31%) underwent PET. PET compared to SPECT demonstrated higher sensitivity (88% vs 60%), specificity (56% vs 14%), positive predictive value (64% vs 20%), negative predictive value (83% vs 50%), and overall superior diagnostic accuracy (AUC .72 (95% CI .50-.93) vs .37 (95% CI .20-.54), P = .01) to detect obstructive CAD. LBBB/VPR-related septal and anteroseptal defects were more common with SPECT compared to PET (septal: 72% vs 17%, P = .001; anteroseptal: 47% vs 8%, P = .02). CONCLUSIONS: PET has higher diagnostic accuracy when compared to SPECT for the detection of obstructive CAD in patients with LBBB or VPR.

Special Section Guest Editorial: Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine.

Guest editors Scott D. Metzler, Samuel Matej, and J. Webster Stayman provide an introduction to the Special Section on Three-Dimensional Image Reconstruction in Nuclear Medicine, PET, and CT.

Quantification of defect contrast in microSPECT imaging of a myocardial phantom.

Ischemic heart disease remains a significant public health concern, accentuating the importance of basic research and therapeutic studies of small animals in which myocardial changes can be reproducibly detected and quantified. Few or no studies have investigated the performance of microSPECT in quantifying myocardial lesions. We utilized three versions of a multi-compartment phantom containing two left ventricular myocardial compartments (one uniform and one with a transmural 'cold' defect), a ventricular blood pool, and a background compartment, where each version had a different myocardial wall thickness (0.75, 1.0 and 1.25 mm). Each compartment was imaged separately while acquiring list-mode data. The separate compartment data were manipulated into a single data set with a known defect contrast, blood-pool and background activity. Data were processed with background-free defect-contrast values of 0 (no defect), -0.25, -0.5, -0.75, and -1.0 (all defect), three ratios of blood-pool to myocardial activity, 0 (no blood pool activity), 0.1, and 0.2 (20% of the activity in the healthy myocardial compartment), and three ratios of uniform background 0 (no background activity), 0.1 and 0.2, relative to the healthy myocardial compartment. For each wall thickness, defect contrast, blood-pool, and background activity combination, 25 list-mode noise realizations were generated and reconstructed. Volumes of interest were drawn and used to determine mean contrast recovery coefficients (CRCs) over the noise ensembles. We developed a slope-analysis procedure to estimate a single CRC over all contrast levels, with resulting CRC values (for no blood-pool and no background) of 0.848, 0.946, and 0.834 for the 0.75, 1.0, and 1.25 mm wall thicknesses, respectively. We also determined and validated a reprocessing method to calculate an ideal CRC. This work demonstrates the quantitative abilities of microSPECT for myocardial-defect imaging utilizing CRC and establishes a framework for evaluating defect-imaging capabilities in other systems.

Analytic Determination of Rectangular-Pinhole Sensitivity With Penetration.

Modern small-animal SPECT systems use multiple pinhole collimators per detector to increase sensitivity while still maintaining high resolution. This resolution is a combination of aperture resolution combined with detector resolution, which is mitigated by magnification. Higher magnification results in better resolution, but fewer apertures per detector. When multiple pinhole collimators project onto the same detector, those with a rectangular field of view (FOV) can be packed more tightly than those with a circular FOV. In addition, a rectangular aperture can be used to obtain different resolution-sensitivity tradeoffs in the two orthogonal directions. Thus, these rectangular-pinhole collimators can have independent FOVs and independent resolution values in the two directions of the rectangular aperture. Previous work has determined the amount of penetration for circular pinholes (i.e., circular apertures with circular FOVs), where the pinhole walls were modeled as cones. In this work, a formula for the penetrative sensitivity for rectangular apertures with a rectangular FOV is determined. The formula was validated using numerical calculations for various combinations of acceptance angles, aperture sizes, linear attenuation coefficients, and incidence angles.

Quantification of myocardial uptake rate constants in dynamic small-animal SPECT using a cardiac phantom.

Myocardial blood flow and myocardial blood flow reserve (MBFR) measurements are often used clinically to quantify coronary microvascular function. Developing imaging-based methods to measure MBFR for research in mice would be advantageous for evaluating new treatment methods for coronary microvascular disease (CMVD), yet this is more challenging in mice than in humans. This work investigates microSPECT's quantitative capabilities of cardiac imaging by utilizing a multi-part cardiac phantom and applying a known kinetic model to synthesize kinetic data from static data, allowing for assessment of kinetic modeling accuracy. The phantom was designed with four main components: two left-ventricular (LV) myocardial sections and two LV blood-pool sections, sized for end-systole (ES) and end-diastole (ED). Each section of the phantom was imaged separately while acquiring list-mode data. These static, separate-compartment data were manipulated into synthetic dynamic data using a kinetic model representing the myocardium and blood-pool activity concentrations over time and then combined into a set of dynamic image frames and reconstructed. Regions of interest were drawn on the resulting images, and kinetic parameters were estimated. This process was performed for three tracer uptake values (K 1), three myocardial wall thicknesses, ten filter parameters, and 20 iterations for 25 noise ensembles. The degree of filtering and iteration number were optimized to minimize the root mean-squared error (RMSE) of K 1 values, with the largest number of iterations and minimal filtering yielding the lowest error. Using the optimized parameters, K 1 was determined with reasonable error (~3% RMSE) over all wall thicknesses and K 1 input values. This work demonstrates that accurate and precise measurements of K 1 are possible for the U-SPECT+ system used in this study, for several different uptake rates and LV dimensions. Additionally, it allows for future investigation utilizing other imaging systems, including PET studies with any radiotracer, as well as with additional phantom parts containing lesions.

Model-Based Normalization of a Fractional-Crystal Collimator for Small-Animal PET Imaging.

Previously, we proposed to use a coincidence collimator to achieve fractional-crystal resolution in PET imaging. We have designed and fabricated a collimator prototype for a small-animal PET scanner, A-PET. To compensate for imperfections in the fabricated collimator prototype, collimator normalization, as well as scanner normalization, is required to reconstruct quantitative and artifact-free images. In this study, we develop a normalization method for the collimator prototype based on the A-PET normalization using a uniform cylinder phantom. We performed data acquisition without the collimator for scanner normalization first, and then with the collimator from eight different rotation views for collimator normalization. After a reconstruction without correction, we extracted the cylinder parameters from which we generated expected emission sinograms. Single scatter simulation was used to generate the scattered sinograms. We used the least-squares method to generate the normalization coefficient for each LOR based on measured, expected and scattered sinograms. The scanner and collimator normalization coefficients were factorized by performing two normalizations separately. The normalization methods were also verified using experimental data acquired from A-PET with and without the collimator. In summary, we developed a model-base collimator normalization that can significantly reduce variance and produce collimator normalization with adequate statistical quality within feasible scan time.

Design of a dual-resolution collimator for preclinical cardiac SPECT with a stationary triple-detector system.

PURPOSE: One approach to preclinical single-photon emission computed tomography (SPECT) imaging that provides both high resolution and high sensitivity is based on imaging a mouse inside a collimating tube; many magnified pinhole projection images from a small target region, e.g., the heart, can be recorded simultaneously on multiple detectors with little multiplexing since each pinhole aperture's opening angle is restricted to view mostly the target organ. However, to obtain complete data for reconstruction, it may be necessary to scan the mouse through the target region of the tube. The authors are developing a different approach based on acquisition and reconstruction of both low-resolution and high-resolution projection data acquired sequentially through many pinholes embedded in two tungsten tube sections of different diameters, a "scout" section and a high-resolution section, placed end-to-end along the axis of a triple-head clinical SPECT scanner. This paper describes the design procedures used to determine the geometric parameters of two new collimator-tube sections, as well as one approach for joint reconstruction of data acquired from both sections. METHODS: The high-resolution section was designed by projecting as many pinhole views of a simulated mouse heart as possible over each detector's camera, with no overlapping of heart projections and minimal overlapping between adjacent "hot" organ and cardiac projections. The authors then jointly optimized the geometric design of the scout section for a triple-detector camera system, as well as the number of maximum-likelihood expectation maximization (MLEM) iterations required to provide minimum mean-squared error of reconstructed voxel counts throughout a 7-cm axial range, with the constraints of fixed, 2.4-mm scout system resolution at the tube center for all apertures, limited multiplexing, and no detector motion. Simulated mouse projection data from both tube sections were then reconstructed to illustrate a simple approach for using high-resolution data to improve the whole-body scout images within a cylindrical region surrounding the heart. RESULTS: The 2-cm-inner-radius high-resolution tube section accommodated 87 platinum-iridium pinhole inserts, each with a 0.3-mm square aperture; their radial distances from the centerline of the system ranged from 2.2 to 3.0 cm. The optimal radial distance to the closest scout pinhole and optimal number of MLEM iterations were 4.4 cm and 35 iterations, respectively, and the radial distances of the 39 scout pinholes ranged from 4.4 to 4.8 cm; aperture sizes ranged from 1.1 to 1.7 mm transaxially and 0.9-1.5 mm axially. After including data from the high-resolution section viewing the heart region into whole-body mouse reconstructions from scout data, the authors obtained high-resolution images of the heart, embedded within lower resolution images of the body, with minimal artifacts. CONCLUSIONS: The authors have optimized a dual-resolution collimator tube that provides both whole-body projections of a mouse and more targeted projections centered on the heart that can be jointly reconstructed to obtain high-resolution images of the heart embedded within lower-resolution whole-body images.

Reconstruction of multiple-pinhole micro-SPECT data using origin ensembles.

PURPOSE: The authors are currently developing a dual-resolution multiple-pinhole microSPECT imaging system based on three large NaI(Tl) gamma cameras. Two multiple-pinhole tungsten collimator tubes will be used sequentially for whole-body "scout" imaging of a mouse, followed by high-resolution (hi-res) imaging of an organ of interest, such as the heart or brain. Ideally, the whole-body image will be reconstructed in real time such that data need only be acquired until the area of interest can be visualized well-enough to determine positioning for the hi-res scan. The authors investigated the utility of the origin ensemble (OE) algorithm for online and offline reconstructions of the scout data. This algorithm operates directly in image space, and can provide estimates of image uncertainty, along with reconstructed images. Techniques for accelerating the OE reconstruction were also introduced and evaluated. METHODS: System matrices were calculated for our 39-pinhole scout collimator design. SPECT projections were simulated for a range of count levels using the MOBY digital mouse phantom. Simulated data were used for a comparison of OE and maximum-likelihood expectation maximization (MLEM) reconstructions. The OE algorithm convergence was evaluated by calculating the total-image entropy and by measuring the counts in a volume-of-interest (VOI) containing the heart. Total-image entropy was also calculated for simulated MOBY data reconstructed using OE with various levels of parallelization. RESULTS: For VOI measurements in the heart, liver, bladder, and soft-tissue, MLEM and OE reconstructed images agreed within 6%. Image entropy converged after ∼2000 iterations of OE, while the counts in the heart converged earlier at ∼200 iterations of OE. An accelerated version of OE completed 1000 iterations in <9 min for a 6.8M count data set, with some loss of image entropy performance, whereas the same dataset required ∼79 min to complete 1000 iterations of conventional OE. A combination of the two methods showed decreased reconstruction time and no loss of performance when compared to conventional OE alone. CONCLUSIONS: OE-reconstructed images were found to be quantitatively and qualitatively similar to MLEM, yet OE also provided estimates of image uncertainty. Some acceleration of the reconstruction can be gained through the use of parallel computing. The OE algorithm is useful for reconstructing multiple-pinhole SPECT data and can be easily modified for real-time reconstruction.

A unified Fourier theory for time-of-flight PET data.

Fully 3D time-of-flight (TOF) PET scanners offer the potential of previously unachievable image quality in clinical PET imaging. TOF measurements add another degree of redundancy for cylindrical PET scanners and make photon-limited TOF-PET imaging more robust than non-TOF PET imaging. The data space for 3D TOF-PET data is five-dimensional with two degrees of redundancy. Previously, consistency equations were used to characterize the redundancy of TOF-PET data. In this paper, we first derive two Fourier consistency equations and Fourier-John equation for 3D TOF PET based on the generalized projection-slice theorem; the three partial differential equations (PDEs) are the dual of the sinogram consistency equations and John's equation. We then solve the three PDEs using the method of characteristics. The two degrees of entangled redundancy of the TOF-PET data can be explicitly elicited and exploited by the solutions of the PDEs along the characteristic curves, which gives a complete understanding of the rich structure of the 3D x-ray transform with TOF measurement. Fourier rebinning equations and other mapping equations among different types of PET data are special cases of the general solutions. We also obtain new Fourier rebinning and consistency equations (FORCEs) from other special cases of the general solutions, and thus we obtain a complete scheme to convert among different types of PET data: 3D TOF, 3D non-TOF, 2D TOF and 2D non-TOF data. The new FORCEs can be used as new Fourier-based rebinning algorithms for TOF-PET data reduction, inverse rebinnings for designing fast projectors, or consistency conditions for estimating missing data. Further, we give a geometric interpretation of the general solutions--the two families of characteristic curves can be obtained by respectively changing the azimuthal and co-polar angles of the biorthogonal coordinates in Fourier space. We conclude the unified Fourier theory by showing that the Fourier consistency equations are necessary and sufficient for 3D x-ray transform with TOF measurement. Finally, we give numerical examples of inverse rebinning for a 3D TOF PET and Fourier-based rebinning for a 2D TOF PET using the FORCEs to show the efficacy of the unified Fourier solutions.

Fourier rebinning and consistency equations for time-of-flight PET planograms.

Due to the unique geometry, dual-panel PET scanners have many advantages in dedicated breast imaging and on-board imaging applications since the compact scanners can be combined with other imaging and treatment modalities. The major challenges of dual-panel PET imaging are the limited-angle problem and data truncation, which can cause artifacts due to incomplete data sampling. The time-of-flight (TOF) information can be a promising solution to reduce these artifacts. The TOF planogram is the native data format for dual-panel TOF PET scanners, and the non-TOF planogram is the 3D extension of linogram. The TOF planograms is five-dimensional while the objects are three-dimensional, and there are two degrees of redundancy. In this paper, we derive consistency equations and Fourier-based rebinning algorithms to provide a complete understanding of the rich structure of the fully 3D TOF planograms. We first derive two consistency equations and John's equation for 3D TOF planograms. By taking the Fourier transforms, we obtain two Fourier consistency equations and the Fourier-John equation, which are the duals of the consistency equations and John's equation, respectively. We then solve the Fourier consistency equations and Fourier-John equation using the method of characteristics. The two degrees of entangled redundancy of the 3D TOF data can be explicitly elicited and exploited by the solutions along the characteristic curves. As the special cases of the general solutions, we obtain Fourier rebinning and consistency equations (FORCEs), and thus we obtain a complete scheme to convert among different types of PET planograms: 3D TOF, 3D non-TOF, 2D TOF and 2D non-TOF planograms. The FORCEs can be used as Fourier-based rebinning algorithms for TOF-PET data reduction, inverse rebinnings for designing fast projectors, or consistency conditions for estimating missing data. As a byproduct, we show the two consistency equations are necessary and sufficient for 3D TOF planograms. Finally, we give numerical examples of implementation of a fast 2D TOF planogram projector and Fourier-based rebinning for a 2D TOF planograms using the FORCEs to show the efficacy of the Fourier-based solutions.

Investigation of imaging properties for submillimeter rectangular pinholes.

PURPOSE: Recently, a multipinhole collimator with inserts that have both rectangular apertures and rectangular fields of view (FOVs) has been proposed for SPECT imaging since it can tile the projection onto the detector efficiently and the FOVs in transverse and axial directions become separable. The purpose of this study is to investigate the image properties of rectangular-aperture pinholes with submillimeter apertures sizes. METHODS: In this work, the authors have conducted sensitivity and FOV experiments for 18 replicates of a prototype insert fabricated in platinum/iridium (Pt/Ir) alloy with submillimeter square-apertures. A sin(q)θ fit to the experimental sensitivity has been performed for these inserts. For the FOV measurement, the authors have proposed a new formula to calculate the projection intensity of a flood image on the detector, taking into account the penumbra effect. By fitting this formula to the measured projection data, the authors obtained the acceptance angles. RESULTS: The mean (standard deviation) of fitted sensitivity exponents q and effective edge lengths we were, respectively, 10.8 (1.8) and 0.38 mm (0.02 mm), which were close to the values, 7.84 and 0.396 mm, obtained from Monte Carlo calculations using the parameters of the designed inserts. For the FOV measurement, the mean (standard deviation) of the transverse and axial acceptances were 35.0° (1.2°) and 30.5° (1.6°), which are in good agreement with the designed values (34.3° and 29.9°). CONCLUSIONS: These results showed that the physical properties of the fabricated inserts with submillimeter aperture size matched our design well.

Transmission-less attenuation estimation from time-of-flight PET histo-images using consistency equations.

In positron emission tomography (PET) imaging, attenuation correction with accurate attenuation estimation is crucial for quantitative patient studies. Recent research showed that the attenuation sinogram can be determined up to a scaling constant utilizing the time-of-flight information. The TOF-PET data can be naturally and efficiently stored in a histo-image without information loss, and the radioactive tracer distribution can be efficiently reconstructed using the DIRECT approaches. In this paper, we explore transmission-less attenuation estimation from TOF-PET histo-images. We first present the TOF-PET histo-image formation and the consistency equations in the histo-image parameterization, then we derive a least-squares solution for estimating the directional derivatives of the attenuation factors from the measured emission histo-images. Finally, we present a fast solver to estimate the attenuation factors from their directional derivatives using the discrete sine transform and fast Fourier transform while considering the boundary conditions. We find that the attenuation histo-images can be uniquely determined from the TOF-PET histo-images by considering boundary conditions. Since the estimate of the attenuation directional derivatives can be inaccurate for LORs tangent to the patient boundary, external sources, e.g. a ring or annulus source, might be needed to give an accurate estimate of the attenuation gradient for such LORs. The attenuation estimation from TOF-PET emission histo-images is demonstrated using simulated 2D TOF-PET data.

Review of SPECT collimator selection, optimization, and fabrication for clinical and preclinical imaging.

In single photon emission computed tomography, the choice of the collimator has a major impact on the sensitivity and resolution of the system. Traditional parallel-hole and fan-beam collimators used in clinical practice, for example, have a relatively poor sensitivity and subcentimeter spatial resolution, while in small-animal imaging, pinhole collimators are used to obtain submillimeter resolution and multiple pinholes are often combined to increase sensitivity. This paper reviews methods for production, sensitivity maximization, and task-based optimization of collimation for both clinical and preclinical imaging applications. New opportunities for improved collimation are now arising primarily because of (i) new collimator-production techniques and (ii) detectors with improved intrinsic spatial resolution that have recently become available. These new technologies are expected to impact the design of collimators in the future. The authors also discuss concepts like septal penetration, high-resolution applications, multiplexing, sampling completeness, and adaptive systems, and the authors conclude with an example of an optimization study for a parallel-hole, fan-beam, cone-beam, and multiple-pinhole collimator for different applications.

LOR-interleaving image reconstruction for PET imaging with fractional-crystal collimation.

Positron emission tomography (PET) has become an important modality in medical and molecular imaging. However, in most PET applications, the resolution is still mainly limited by the physical crystal sizes or the detector's intrinsic spatial resolution. To achieve images with better spatial resolution in a central region of interest (ROI), we have previously proposed using collimation in PET scanners. The collimator is designed to partially mask detector crystals to detect lines of response (LORs) within fractional crystals. A sequence of collimator-encoded LORs is measured with different collimation configurations. This novel collimated scanner geometry makes the reconstruction problem challenging, as both detector and collimator effects need to be modeled to reconstruct high-resolution images from collimated LORs. In this paper, we present a LOR-interleaving (LORI) algorithm, which incorporates these effects and has the advantage of reusing existing reconstruction software, to reconstruct high-resolution images for PET with fractional-crystal collimation. We also develop a 3D ray-tracing model incorporating both the collimator and crystal penetration for simulations and reconstructions of the collimated PET. By registering the collimator-encoded LORs with the collimator configurations, high-resolution LORs are restored based on the modeled transfer matrices using the non-negative least-squares method and EM algorithm. The resolution-enhanced images are then reconstructed from the high-resolution LORs using the MLEM or OSEM algorithm. For validation, we applied the LORI method to a small-animal PET scanner, A-PET, with a specially designed collimator. We demonstrate through simulated reconstructions with a hot-rod phantom and MOBY phantom that the LORI reconstructions can substantially improve spatial resolution and quantification compared to the uncollimated reconstructions. The LORI algorithm is crucial to improve overall image quality of collimated PET, which can have significant implications in preclinical and clinical ROI imaging applications.