Reducing the influence of the partial volume effect on SPECT activity quantitation with 3D modelling of spatial resolution in iterative reconstruction

What scans we will read: imaging instrumentation trends in clinical oncology

Oncological diseases account for a significant portion of the burden on public healthcare systems with associated costs driven primarily by complex and long-lasting therapies. Through the visualization of patient-specific morphology and functional-molecular pathways, cancerous tissue can be detected and characterized non-invasively, so as to provide referring oncologists with essential information to support therapy management decisions. Following the onset of stand-alone anatomical and functional imaging, we witness a push towards integrating molecular image information through various methods, including anato-metabolic imaging (e.g., PET/CT), advanced MRI, optical or ultrasound imaging.This perspective paper highlights a number of key technological and methodological advances in imaging instrumentation related to anatomical, functional, molecular medicine and hybrid imaging, that is understood as the hardware-based combination of complementary anatomical and molecular imaging. These include novel detector technologies for ionizing radiation used in CT and nuclear medicine imaging, and novel system developments in MRI and optical as well as opto-acoustic imaging. We will also highlight new data processing methods for improved non-invasive tissue characterization. Following a general introduction to the role of imaging in oncology patient management we introduce imaging methods with well-defined clinical applications and potential for clinical translation. For each modality, we report first on the status quo and, then point to perceived technological and methodological advances in a subsequent status go section. Considering the breadth and dynamics of these developments, this perspective ends with a critical reflection on where the authors, with the majority of them being imaging experts with a background in physics and engineering, believe imaging methods will be in a few years from now.Overall, methodological and technological medical imaging advances are geared towards increased image contrast, the derivation of reproducible quantitative parameters, an increase in volume sensitivity and a reduction in overall examination time. To ensure full translation to the clinic, this progress in technologies and instrumentation is complemented by advances in relevant acquisition and image-processing protocols and improved data analysis. To this end, we should accept diagnostic images as "data", and - through the wider adoption of advanced analysis, including machine learning approaches and a "big data" concept - move to the next stage of non-invasive tumour phenotyping. The scans we will be reading in 10 years from now will likely be composed of highly diverse multi-dimensional data from multiple sources, which mandate the use of advanced and interactive visualization and analysis platforms powered by Artificial Intelligence (AI) for real-time data handling by cross-specialty clinical experts with a domain knowledge that will need to go beyond that of plain imaging.

SPECT/CT: an update on technological developments and clinical applications

Functional nuclear medicine imaging with single-photon emission CT (SPECT) in combination with anatomical CT has been commercially available since the beginning of this century. The combination of the two modalities has improved both the sensitivity and specificity of many clinical applications and CT in conjunction with SPECT that allows for spatial overlay of the SPECT data on good anatomy images. Introduction of diagnostic CT units as part of the SPECT/CT system has also potentially allowed for a more cost-efficient use of the equipment. Most of the SPECT systems available are based on the well-known Anger camera principle with NaI(Tl) as a scintillation material, parallel-hole collimators and multiple photomultiplier tubes, which, from the centroid of the scintillation light, determine the position of an event. Recently, solid-state detectors using cadmium-zinc-telluride became available and clinical SPECT cameras employing multiple pinhole collimators have been developed and introduced in the market. However, even if new systems become available with better hardware, the SPECT reconstruction will still be affected by photon attenuation and scatter and collimator response. Compensation for these effects is needed even for qualitative studies to avoid artefacts leading to false positives. This review highlights the recent progress for both new SPECT cameras systems as well as for various data-processing and compensation methods.

Collimator optimization and collimator-detector response compensation in myocardial perfusion SPECT using the ideal observer with and without model mismatch and an anthropomorphic model observer

The collimator is the primary factor that determines the spatial resolution and noise tradeoff in myocardial perfusion SPECT images. In this paper, the goal was to find the collimator that optimizes the image quality in terms of a perfusion defect detection task. Since the optimal collimator could depend on the level of approximation of the collimator-detector response (CDR) compensation modeled in reconstruction, we performed this optimization for the cases of modeling the full CDR (including geometric, septal penetration and septal scatter responses), the geometric CDR, or no model of the CDR. We evaluated the performance on the detection task using three model observers. Two observers operated on data in the projection domain: the Ideal Observer (IO) and IO with Model-Mismatch (IO-MM). The third observer was an anthropomorphic Channelized Hotelling Observer (CHO), which operated on reconstructed images. The projection-domain observers have the advantage that they are computationally less intensive. The IO has perfect knowledge of the image formation process, i.e. it has a perfect model of the CDR. The IO-MM takes into account the mismatch between the true (complete and accurate) model and an approximate model, e.g. one that might be used in reconstruction. We evaluated the utility of these projection domain observers in optimizing instrumentation parameters. We investigated a family of 8 parallel-hole collimators, spanning a wide range of resolution and sensitivity tradeoffs, using a population of simulated projection (for the IO and IO-MM) and reconstructed (for the CHO) images that included background variability. We simulated anterolateral and inferior perfusion defects with variable extents and severities. The area under the ROC curve was estimated from the IO, IO-MM, and CHO test statistics and served as the figure-of-merit. The optimal collimator for the IO had a resolution of 9-11 mm FWHM at 10 cm, which is poorer resolution than typical collimators used for MPS. When the IO-MM and CHO used a geometric or no model of the CDR, the optimal collimator shifted toward higher resolution than that obtained using the IO and the CHO with full CDR modeling. With the optimal collimator, the IO-MM and CHO using geometric modeling gave similar performance to full CDR modeling. Collimators with poorer resolution were optimal when CDR modeling was used. The agreement of rankings between the IO-MM and CHO confirmed that the IO-MM is useful for optimization tasks when model mismatch is present due to its substantially reduced computational burden compared to the CHO.

Comparison of the scanning linear estimator (SLE) and ROI methods for quantitative SPECT imaging

In quantitative emission tomography, tumor activity is typically estimated from calculations on a region of interest (ROI) identified in the reconstructed slices. In these calculations, unpredictable bias arising from the null functions of the imaging system affects ROI estimates. The magnitude of this bias depends upon the tumor size and location. In prior work it has been shown that the scanning linear estimator (SLE), which operates on the raw projection data, is an unbiased estimator of activity when the size and location of the tumor are known. In this work, we performed analytic simulation of SPECT imaging with a parallel-hole medium-energy collimator. Distance-dependent system spatial resolution and non-uniform attenuation were included in the imaging simulation. We compared the task of activity estimation by the ROI and SLE methods for a range of tumor sizes (diameter: 1-3 cm) and activities (contrast ratio: 1-10) added to uniform and non-uniform liver backgrounds. Using the correct value for the tumor shape and location is an idealized approximation to how task estimation would occur clinically. Thus we determined how perturbing this idealized prior knowledge impacted the performance of both techniques. To implement the SLE for the non-uniform background, we used a novel iterative algorithm for pre-whitening stationary noise within a compact region. Estimation task performance was compared using the ensemble mean-squared error (EMSE) as the criterion. The SLE method performed substantially better than the ROI method (i.e. EMSE(SLE) was 23-174 times lower) when the background is uniform and tumor location and size are known accurately. The variance of the SLE increased when a non-uniform liver texture was introduced but the EMSE(SLE) continued to be 5-20 times lower than the ROI method. In summary, SLE outperformed ROI under almost all conditions that we tested.

Characterization of parallel-hole collimator using Monte Carlo Simulation

Objective:

Accuracy of in vivo activity quantification improves after the correction of penetrated and scattered photons. However, accurate assessment is not possible with physical experiment. We have used Monte Carlo Simulation to accurately assess the contribution of penetrated and scattered photons in the photopeak window.

Materials and Methods:

Simulations were performed with Simulation of Imaging Nuclear Detectors Monte Carlo Code. The simulations were set up in such a way that it provides geometric, penetration, and scatter components after each simulation and writes binary images to a data file. These components were analyzed graphically using Microsoft Excel (Microsoft Corporation, USA). Each binary image was imported in software (ImageJ) and logarithmic transformation was applied for visual assessment of image quality, plotting profile across the center of the images and calculating full width at half maximum (FWHM) in horizontal and vertical directions.

Results:

The geometric, penetration, and scatter at 140 keV for low-energy general-purpose were 93.20%, 4.13%, 2.67% respectively. Similarly, geometric, penetration, and scatter at 140 keV for low-energy high-resolution (LEHR), medium-energy general-purpose (MEGP), and high-energy general-purpose (HEGP) collimator were (94.06%, 3.39%, 2.55%), (96.42%, 1.52%, 2.06%), and (96.70%, 1.45%, 1.85%), respectively. For MEGP collimator at 245 keV photon and for HEGP collimator at 364 keV were 89.10%, 7.08%, 3.82% and 67.78%, 18.63%, 13.59%, respectively.

Conclusion:

Low-energy general-purpose and LEHR collimator is best to image 140 keV photon. HEGP can be used for 245 keV and 364 keV; however, correction for penetration and scatter must be applied if one is interested to quantify the in vivo activity of energy 364 keV. Due to heavy penetration and scattering, 511 keV photons should not be imaged with HEGP collimator.

Oxytocin and vasopressin systems in genetic syndromes and neurodevelopmental disorders

Oxytocin (OT) and arginine vasopressin (AVP) are two small, related neuropeptide hormones found in many mammalian species, including humans. Dysregulation of these neuropeptides have been associated with changes in behavior, especially social interactions. We review how the OT and AVP systems have been investigated in Autism Spectrum Disorder (ASD), Prader-Willi Syndrome (PWS), Williams Syndrome (WS) and Fragile X syndrome (FXS). All of these neurodevelopmental disorders (NDD) are marked by social deficits. While PWS, WS and FXS have identified genetic mutations, ASD stems from multiple genes with complex interactions. Animal models of NDD are invaluable for studying the role and relatedness of OT and AVP in the developing brain. We present data from a FXS mouse model affecting the fragile X mental retardation 1 (Fmr1) gene, resulting in decreased OT and AVP staining cells in some brain regions. Reviewing the research about OT and AVP in these NDD suggests that altered OT pathways may be downstream from different etiological factors and perturbations in development. This has implications for ongoing studies of the therapeutic application of OT in NDD. This article is part of a Special Issue entitled Oxytocin and Social Behav.

Quantitative SPECT/CT

Conventional nuclear medical imaging uses radiopharmaceuticals labeled by single-photon emitters such as Tc-99m, I-123, or I-131 in vivo. Classical clinical examples are the study of bone metabolism by bone scintigraphy with the Tc-99m-labeled polyphosphonates or of iodine transport into the thyroid gland using Tc-99m-pertechnetate. With single-photon emission-computed tomography (SPECT), the distribution of these radiopharmaceuticals within the human body is three-dimensionally visualized. Contrary to positron emission tomography (PET), current SPECT technology does not allow the quantification of regional values of radioactivity tissue concentration as SPECT images are grossly compromised by artifacts caused by photon scatter and attenuation. With the advent of hybrid imaging systems combining a SPECT camera with an X-ray computerized (CT) scanner in one gantry, reliable corrections for these artifacts seem possible, allowing truly quantitative SPECT.

The effect of adaptive iterative dose reduction on image quality in 320-detector row CT coronary angiography

OBJECTIVE

To evaluate the effect of adaptive iterative dose reduction (AIDR) on image noise and image quality as compared with standard filtered back projection (FBP) in 320-detector row CT coronary angiography (CTCA).

METHODS

50 patients (14 females, mean age 68 ± 9 years) who underwent CTCA (100 kV or 120 kV, 400-580 mA) within a single heartbeat were enrolled. Studies were reconstructed with FBP and subsequently AIDR. Image noise, vessel contrast and contrast-to-noise ratio (CNR) in the coronary arteries were evaluated. Overall image quality for coronary arteries was assessed using a five-point scale (1, non-diagnostic; 5, excellent).

RESULTS

All the examinations were performed in a single heartbeat. Image noise in the aorta was significantly lower in data sets reconstructed with AIDR than in those reconstructed with FBP (21.4 ± 3.1 HU vs 36.9 ± 4.5 HU; p<0.001). No significant differences were observed between FBP and AIDR for the mean vessel contrast (HU) in the proximal coronary arteries. Consequently, CNRs in the proximal coronary arteries were higher in the AIDR group than in the FBP group (p<0.001). The mean image quality score was improved by AIDR (3.75 ± 0.38 vs 4.24 ± 0.38; p<0.001).

CONCLUSION

The use of AIDR reduces image noise and improves image quality in 320-detector row CTCA.

Recent advances in iterative reconstruction for clinical SPECT/PET and CT

Statistical iterative reconstruction is now widely used in clinical practice and has contributed to significant improvement in image quality in recent years. Although primarily used for reconstruction in emission tomography (both single photon emission computed tomography (SPECT) and positron emission tomography (PET)) there is increasing interest in also applying similar algorithms to x-ray computed tomography (CT). There is increasing complexity in the factors that are included in the reconstruction, a demonstration of the versatility of the approach. Research continues with exploration of methods for further improving reconstruction quality with effective correction for various sources of artefact.

Absolute quantification in SPECT

Single-photon emission computed tomography (SPECT) allows the three-dimensional visualization of radioactivity within the human body and is widely used for clinical purposes. In SPECT, image quality is compromised by several factors including photon attenuation, photon scatter, the partial volume effect, and motion artefacts. These variables also confound the capacity of SPECT to quantify the concentration of radioactivity within given volumes of interest in absolute units, e.g. as kilobecquerels per cubic centimetre. In the last decade, considerable technical progress has been achieved in SPECT image reconstruction, involving, in particular, the development of iterative image reconstruction techniques. Furthermore, hybrid cameras integrating a SPECT camera with an X-ray CT scanner have become commercially available. These systems allow the acquisition of SPECT and CT datasets registered to each other with a high anatomical accuracy. First studies have shown that iterative SPECT image reconstruction techniques incorporating information from SPECT/CT image datasets greatly increase the accuracy of SPECT in quantifying radioactivity concentrations in phantoms and also in humans. This new potential of SPECT may improve not only diagnostic accuracy, but also dosimetry for internal radiotherapy.

Image quality of multiplanar reconstruction of pulmonary CT scans using adaptive statistical iterative reconstruction

OBJECTIVE

We investigated the image quality of multiplanar reconstruction (MPR) using adaptive statistical iterative reconstruction (ASIR).

METHODS

Inflated and fixed lungs were scanned with a garnet detector CT in high-resolution mode (HR mode) or non-high-resolution (HR) mode, and MPR images were then reconstructed. Observers compared 15 MPR images of ASIR (40%) and ASIR (80%) with those of ASIR (0%), and assessed image quality using a visual five-point scale (1, definitely inferior; 5, definitely superior), with particular emphasis on normal pulmonary structures, artefacts, noise and overall image quality.

RESULTS

The mean overall image quality scores in HR mode were 3.67 with ASIR (40%) and 4.97 with ASIR (80%). Those in non-HR mode were 3.27 with ASIR (40%) and 3.90 with ASIR (80%). The mean artefact scores in HR mode were 3.13 with ASIR (40%) and 3.63 with ASIR (80%), but those in non-HR mode were 2.87 with ASIR (40%) and 2.53 with ASIR (80%). The mean scores of the other parameters were greater than 3, whereas those in HR mode were higher than those in non-HR mode. There were significant differences between ASIR (40%) and ASIR (80%) in overall image quality (p<0.01). Contrast medium in the injection syringe was scanned to analyse image quality; ASIR did not suppress the severe artefacts of contrast medium.

CONCLUSION

In general, MPR image quality with ASIR (80%) was superior to that with ASIR (40%). However, there was an increased incidence of artefacts by ASIR when CT images were obtained in non-HR mode.

Adaptive statistical iterative reconstruction technique for pulmonary CT: image quality of the cadaveric lung on standard- and reduced-dose CT

RATIONALE AND OBJECTIVES

To evaluate thin-section computed tomography (CT) images of the lung reconstructed using adaptive statistical iterative reconstruction (ASIR) on standard- and reduced-dose CT.

MATERIALS AND METHODS

Eleven cadaveric lungs were scanned by multidetector-row CT with two different tube currents (standard dose, 400 mA; reduced dose, 10 mA). The degree of ASIR was classified into six different levels: 0% (non-ASIR), 20%, 40%, 60%, 80%, and 100% (maximum-ASIR). The ASIR (20%, 60%, and 100%) images were compared with the ASIR (0%) images and assessed visually by three independent observers for image quality using a 7-point scale. The evaluation items included abnormal CT findings, normal lung structures, and subjective visual noise. The median scores assigned by the three observers were analyzed statistically. Quantitative noise was calculated by measuring the standard deviation in a circular region of interest on each selected image of ASIR (0%-100%).

RESULTS

On standard-dose CT, the overall image quality significantly improved with increasing degree of ASIR (P ≤ .009, Wilcoxon signed-ranks test with Bonferroni correction). As ASIR increased, however, intralobular reticular opacities and peripheral vessels tended to be obscure. On reduced-dose CT, the overall image quality of ASIR (100%) was significantly better than that of ASIR (20%) (P ≤ .009). As ASIR increased, however, intralobular reticular opacities tended to be obscure. Using ASIR significantly reduced subjective and quantitative image noise on both standard- and reduced-dose CT (P < .001, Bonferroni/Dunn's method).

CONCLUSION

ASIR improves the image quality by decreasing image noise. Maximum-ASIR may be needed for improving image quality on highly reduced-dose CT. However, excessive ASIR may obscure subtle shadows.

Diminishing the impact of the partial volume effect in cardiac SPECT perfusion imaging

The partial volume effect (PVE) significantly restricts the absolute quantification of regional myocardial uptake and thereby limits the accuracy of absolute measurement of blood flow and coronary flow reserve by SPECT. The template-projection-reconstruction method has been previously developed for PVE compensation. This method assumes the availability of coregistered high-spatial resolution anatomical information as is now becoming available with commercial dual-modality imaging systems such as SPECT/CTs. The objective of this investigation was to determine the extent to which the impact of the PVE on cardiac perfusion SPECT imaging can be diminished if coregistered high-spatial resolution anatomical information is available. For this investigation the authors introduced an additional parameter into the template-projection-reconstruction compensation equation called the voxel filling fraction (F). This parameter specifies the extent to which structure edge voxels in the emission reconstruction are filled by the structure in question as determined by the higher spatial-resolution imaging modality and the fractional presence of the structure at different states of physiological motion as in combining phases of cardiac motion. During correction the removal of spillover to the cardiac region from the surrounding structures is performed first by using reconstructed templates of neighboring structures (liver, blood pool, lungs) to calculate spillover fractions. This is followed by determining recovery coefficients for all voxels within the heart wall from the reconstruction of the template projections of the left and right ventricles (LV and RV). The emission data are subsequently divided by these recovery coefficients taking into account the filling fraction F. The mathematical cardiac torso phantom was used for investigation correction of PVE for a normal LV distribution, a defect in the inferior wall, and a defect in the anterior wall. PVE correction resulted in a dramatic visual reduction in the impact of extracardiac activity, improved the uniformity of the normally perfused heart wall, and enhanced defect visibility without undue noise amplification. No significant artifacts were seen with PVE correction in the presence of mild (one voxel) misregistration. A statistically significant improvement in the accuracy of the count levels within the normal heart wall was also noted. However, residual spillover of counts from within the myocardium creates a bias in regions of decreased wall counts (perfusion defects/abnormal wall motion) when the anatomical imaging modality does not allow definition of templates for defects present in the heart during emission imaging.

Absolute quantitation of myocardial blood flow with (201)Tl and dynamic SPECT in canine: optimisation and validation of kinetic modelling

PURPOSE

201Tl has been extensively used for myocardial perfusion and viability assessment. Unlike 99mTc-labelled agents, such as 99mTc-sestamibi and 99mTc-tetrofosmine, the regional concentration of 201Tl varies with time. This study is intended to validate a kinetic modelling approach for in vivo quantitative estimation of regional myocardial blood flow (MBF) and volume of distribution of 201Tl using dynamic SPECT.

METHODS

Dynamic SPECT was carried out on 20 normal canines after the intravenous administration of 201Tl using a commercial SPECT system. Seven animals were studied at rest, nine during adenosine infusion, and four after beta-blocker administration. Quantitative images were reconstructed with a previously validated technique, employing OS-EM with attenuation-correction, and transmission-dependent convolution subtraction scatter correction. Measured regional time-activity curves in myocardial segments were fitted to two- and three-compartment models. Regional MBF was defined as the influx rate constant (K(1)) with corrections for the partial volume effect, haematocrit and limited first-pass extraction fraction, and was compared with that determined from radio-labelled microspheres experiments.

RESULTS

Regional time-activity curves responded well to pharmacological stress. Quantitative MBF values were higher with adenosine and decreased after beta-blocker compared to a resting condition. MBFs obtained with SPECT (MBF(SPECT)) correlated well with the MBF values obtained by the radio-labelled microspheres (MBF(MS)) (MBF(SPECT) = -0.067 + 1.042 x MBF(MS), p < 0.001). The three-compartment model provided better fit than the two-compartment model, but the difference in MBF values between the two methods was small and could be accounted for with a simple linear regression.

CONCLUSION

Absolute quantitation of regional MBF, for a wide physiological flow range, appears to be feasible using 201Tl and dynamic SPECT.

Image reconstruction from limited Fourier data

We consider the problem of reconstructing a function f with bounded support S from finitely many values of its Fourier transform F. Although f cannot be band limited since it has bounded support, it is typically the case that f can be modeled as the restriction to S of a sigma-band-limited function, say g. Our reconstruction method is based on such a model for f. Of particular interest is the effect of the choice of sigma > 0 on the resolution.

Automatic estimation of detector radial position for contoured SPECT acquisition using CT images on a SPECT/CT system

An algorithm was developed to estimate noncircular orbit (NCO) single-photon emission computed tomography (SPECT) detector radius on a SPECT/CT imaging system using the CT images, for incorporation into collimator resolution modeling for iterative SPECT reconstruction. Simulated male abdominal (arms up), male head and neck (arms down) and female chest (arms down) anthropomorphic phantom, and ten patient, medium-energy SPECT/CT scans were acquired on a hybrid imaging system. The algorithm simulated inward SPECT detector radial motion and object contour detection at each projection angle, employing the calculated average CT image and a fixed Hounsfield unit (HU) threshold. Calculated radii were compared to the observed true radii, and optimal CT threshold values, corresponding to patient bed and clothing surfaces, were found to be between -970 and -950 HU. The algorithm was constrained by the 45 cm CT field-of-view (FOV), which limited the detected radii to < or = 22.5 cm and led to occasional radius underestimation in the case of object truncation by CT. Two methods incorporating the algorithm were implemented: physical model (PM) and best fit (BF). The PM method computed an offset that produced maximum overlap of calculated and true radii for the phantom scans, and applied that offset as a calculated-to-true radius transformation. For the BF method, the calculated-to-true radius transformation was based upon a linear regression between calculated and true radii. For the PM method, a fixed offset of +2.75 cm provided maximum calculated-to-true radius overlap for the phantom study, which accounted for the camera system's object contour detect sensor surface-to-detector face distance. For the BF method, a linear regression of true versus calculated radius from a reference patient scan was used as a calculated-to-true radius transform. Both methods were applied to ten patient scans. For -970 and -950 HU thresholds, the combined overall average root-mean-square (rms) error in radial position for eight patient scans without truncation were 3.37 cm (12.9%) for PM and 1.99 cm (8.6%) for BF, indicating BF is superior to PM in the absence of truncation. For two patient scans with truncation, the rms error was 3.24 cm (12.2%) for PM and 4.10 cm (18.2%) for BF. The slightly better performance of PM in the case of truncation is anomalous, due to FOV edge truncation artifacts in the CT reconstruction, and thus is suspect. The calculated NCO contour for a patient SPECT/CT scan was used with an iterative reconstruction algorithm that incorporated compensation for system resolution. The resulting image was qualitatively superior to the image obtained by reconstructing the data using the fixed radius stored by the scanner. The result was also superior to the image reconstructed using the iterative algorithm provided with the system, which does not incorporate resolution modeling. These results suggest that, under conditions of no or only mild lateral truncation of the CT scan, the algorithm is capable of providing radius estimates suitable for iterative SPECT reconstruction collimator geometric resolution modeling.

Iterative three-dimensional expectation maximization restoration of single photon emission computed tomography images: Application in striatal imaging

Single photon emission computed tomography imaging suffers from poor spatial resolution and high statistical noise. Consequently, the contrast of small structures is reduced, the visual detection of defects is limited and precise quantification is difficult. To improve the contrast, it is possible to include the spatially variant point spread function of the detection system into the iterative reconstruction algorithm. This kind of method is well known to be effective, but time consuming. We have developed a faster method to account for the spatial resolution loss in three dimensions, based on a postreconstruction restoration method. The method uses two steps. First, a noncorrected iterative ordered subsets expectation maximization (OSEM) reconstruction is performed and, in the second step, a three-dimensional (3D) iterative maximum likelihood expectation maximization (ML-EM) a posteriori spatial restoration of the reconstructed volume is done. In this paper, we compare to the standard OSEM-3D method, in three studies (two in simulation and one from experimental data). In the two first studies, contrast, noise, and visual detection of defects are studied. In the third study, a quantitative analysis is performed from data obtained with an anthropomorphic striatal phantom filled with 123-I. From the simulations, we demonstrate that contrast as a function of noise and lesion detectability are very similar for both OSEM-3D and OSEM-R methods. In the experimental study, we obtained very similar values of activity-quantification ratios for different regions in the brain. The advantage of OSEM-R compared to OSEM-3D is a substantial gain of processing time. This gain depends on several factors. In a typical situation, for a 128 x 128 acquisition of 120 projections, OSEM-R is 13 or 25 times faster than OSEM-3D, depending on the calculation method used in the iterative restoration. In this paper, the OSEM-R method is tested with the approximation of depth independent resolution. For the striatum this approximation is appropriate, but for other clinical situations we will need to include a spatially varying response. Such a response is already included in OSEM-3D.

A Monte Carlo and physical phantom evaluation of quantitative In-111 SPECT

Accurate estimation of the 3D in vivo activity distribution is important for dose estimation in targeted radionuclide therapy (TRT). Although SPECT can potentially provide such estimates, SPECT without compensation for image degrading factors is not quantitatively accurate. In this work, we evaluated quantitative SPECT (QSPECT) reconstruction methods that include compensation for various physical effects. Experimental projection data were obtained using a GE VH/Hawkeye system and an RSD torso phantom. Known activities of In-111 chloride were placed in the lungs, liver, heart, background and two spherical compartments with inner diameters of 22 mm and 34 mm. The 3D NCAT phantom with organ activities based on clinically derived In-111 ibritumomab tiuxetan data was used for the Monte Carlo (MC) simulation studies. Low-noise projection data were simulated using previously validated MC simulation methods. Fifty sets of noisy projections with realistic count levels were generated. Reconstructions were performed using the OS-EM algorithm with various combinations of attenuation (A), scatter (S), geometric response (G), collimator-detector response (D) and partial volume compensation (PVC). The QSPECT images from the various combinations of compensations were evaluated in terms of the accuracy and precision of the estimates of the total activity in each organ. For experimental data, the errors in organ activities for ADS and PVC compensation were less than 6.5% except the smaller sphere (-11.9%). For the noisy simulated data, the errors in organ activity for ADS compensation were less than 5.5% except the lungs (20.9%) and blood vessels (15.2%). Errors for other combinations of compensations were significantly (A, AS) or somewhat (AGS) larger. With added PVC, the error in the organ activities improved slightly except for the lungs (11.5%) and blood vessels (3.6%) where the improvement was more substantial. The standard deviation/mean ratios were all less than 1.5%. We conclude that QSPECT methods with appropriate compensations provided accurate In-111 organ activity estimates. For the collimator used, AGS was almost as good as ADS and may be preferable due to the reduced reconstruction time. PVC was important for small structures such as tumours or for organs in close proximity to regions with high activity. The improved quantitative accuracy from QSPECT methods has the potential for improving organ dose estimations in TRT.

Partial volume effect compensation for quantitative brain SPECT imaging

Partial volume (PV) effects degrade the quantitative accuracy of SPECT brain images. In this paper, we extended a PV compensation (PVC) method originally developed for brain PET, the geometric transfer matrix (GTM) method, to brain SPECT using iterative reconstruction-based compensations. In the GTM method a linear transform between the true regional activities and the measured results was assumed. Elements of the GTM were calculated by projecting and reconstructing maps with uniform regions representing different structures. However, with iterative reconstruction methods, especially when reconstruction-based compensation for detector response was applied, we found that it was important to treat the region maps as a perturbation to the reconstructed image in the estimation of the GTM. This modified method, termed perturbation-based GTM (pGTM) was evaluated using Monte Carlo (MC) simulated and experimentally acquired data. Results showed great improvement of the quantitative accuracy in brain SPECT imaging. For MC simulated data, PVC using pGTM reduced the underestimation of striatal activities from 30% to less than 1.2%. For experimental data, PVC using pGTM reduced the underestimation of striatal activities from 36% to less than 7.8%. The underestimation of the striatum to background activity ratio was also improved from 31% to 2.7%.

Attenuation correction single-photon emission computed tomography myocardial perfusion imaging

Clinicians now rely heavily on the results of single-photon emission computed tomography (SPECT) myocardial perfusion imaging for diagnosing coronary disease and for planning therapy. However, the technique is imperfect for these purposes, mainly because of technical limitations, the most prominent of which is the effect of soft-tissue attenuation on apparent tracer distribution. Providers have attempted to compensate for this by a number of indirect approaches. Recently, validated hardware and software solutions for directly correcting image data for soft-tissue attenuation have become widely available commercially. Optimal application requires an understanding of the technical details that differ somewhat from system to system, the quality control prerequisites, knowledge of the importance of the transmission map quality, and how dedicated SPECT and SPECT-computed tomography systems present different challenges. In addition, the clinical literature is expanding rapidly, including studies on diagnostic accuracy, image appearances, quantitative analysis, appropriate patients for attenuation correction, clinical utility, incremental value in relation to ECG-gating, and risk stratification.

Performance evaluation of OSEM reconstruction algorithm incorporating three-dimensional distance-dependent resolution compensation for brain SPECT: a simulation study

Iterative reconstruction techniques such as an ordered subsets-expectation maximization (OSEM) algorithm can easily incorporated various physical models of attenuation or scatter. We implemented OSEM reconstruction algorithm incorporating compensation for distance-dependent blurring due to the collimator in SPECT. The algorithm was examined by computer simulation to estimate the accuracy for brain perfusion study.

METHODS

The detector response was assumed to be a two-dimensional Gauss function and the width of the function varied linearly with the source-to-detector distance. The attenuation compensation (AC) was also included. To investigate the properties of the algorithm, we performed computer simulations with the point source and digital brain phantoms. In the point source phantom, the uniformity of FWHM for the radial, tangential and longitudinal directions was evaluated on the reconstruction image. As for the brain phantom, quantitative accuracy was estimated by comparing the reconstructed images with the true image by the mean square error (MSE) and the ratio of gray and white matter counts (G/W). Both noise free and noisy simulations were examined.

RESULTS

In the point source simulation, FWHM in radial, tangential and longitudinal directions were 14.7, 14.7 and 15.0 mm at the image center and were 15.9, 9.83 and 10.6 mm at a distance of 15 cm from the center by using FBP, respectively. On the other hand, they were 8.12, 8.12 and 7.83 mm at the image center, and were 7.45, 7.44 and 7.01 mm at 15 cm from the center by OSEM with distance-dependent resolution compensation (DRC). An isotropic and stationary resolution was obtained at any location by OSEM with DRC. The spatial resolution was also improved about 6.5 mm by OSEM with DRC at the image center. In the brain phantom simulation, the blurring at the edge of the brain structure was eliminated by using OSEM with both DRC and AC. The G/W was 2.95 and 2.68 for noise free and noisy cases, respectively, when no compensation was performed. But the values for G/W without and with noise became 3.45 and 3.21 with AC only and were improved to 3.75 and 3.71 with both AC and DRC. The G/W approached the true value (4.00) by using OSEM with both AC and DRC even when there was statistical noise.

CONCLUSION

In conclusion, OSEM reconstruction including the distance-dependent resolution compensation algorithm was reasonably successful in achieving isotropic and stationary resolution and improving the quantitative accuracy for brain perfusion SPECT.

Resolution recovery for list-mode reconstruction in SPECT

The purpose of the study was to evaluate the resolution recovery in the list-mode iterative reconstruction algorithm (LMIRA) for SPECT. In this study we compare the performance of the proposed method with other iterative resolution recovery methods for different noise levels. We developed an iterative reconstruction method which uses list-mode data instead of binned data. The new algorithm makes use of a more accurate model of the collimator structure. We compared the SPECT list-mode reconstruction with MLEM, OSEM and RBI, all including resolution recovery. For the evaluation we used Gaussian shaped sources with different FWHM at three different locations and three noise levels. For these distributions we calculated the reconstructed images for a different number of iterations. The absolute error for the reconstructed images was used to evaluate the performance. The performance of all four methods is comparable for the sources located in the centre of the field of view. For the sources located out of the centre, the error of the list-mode method is significantly lower than that of the other methods. Splitting the system model into a separate object-dependent and detector-dependent module gives us a flexible reconstruction method. With this we can very easily adapt the resolution recovery to different collimator types.

Response to incremental doses of dobutamine early after reperfusion is predictive of the degree of myocardial salvage in dogs with experimental acute myocardial infarction

OBJECTIVES

We sought to determine whether the inotropic response to dobutamine might be useful for estimating the extent of viable myocardium soon after reperfusion.

BACKGROUND

Early identification of viable myocardium in the presence of severe left ventricular dysfunction after reperfusion is important for clinical decision making.

METHODS

Nine open-chest dogs had left anterior descending coronary artery occlusion for 40 to 180 min, followed by gradual reperfusion. The systolic thickening response to incremental dobutamine doses was measured with ultrasonic crystals and regional flow by microspheres.

RESULTS

Dogs were divided into two groups based on triphenyl tetralozium chloride infarct size (group 1: 9.3 +/- 3.0% risk area; group 2: 51.1 +/- 4.8%). In group 2 dogs with larger infarcts, regional flow during peak dobutamine was lower than it was in group 1 in endocardial (1.15 +/- 0.22 vs. 2.64 +/- 0.33 mL x min(-1) x g(-1)) and midwall (1.47 +/- 0.32 vs. 2.92 +/- 0.36 mL x min(-1) x g(-1)) layers, and endocardial flow in group 2 failed to increase from baseline (0.96 +/- 0.07 vs. 1.15 +/- 0.22 mL x min(-1) x g(-1)). Group 1 dogs demonstrated a dose dependent increase in systolic thickening with dobutamine versus a blunted response in group 2. The inotropic response to only 10 microg x kg(-1) x min(-1) of dobutamine was predictive of the degree of myocardial salvage.

CONCLUSIONS

In the early postischemic stunning phase of reperfusion, the inotropic response to dobutamine is predictive of the degree of myocardial salvage and ultimate infarct size. The ability to distinguish between stunned versus necrotic myocardium early after reperfusion was most likely due to the presence of subendocardial flow reserve during dobutamine in dogs with predominantly salvaged myocardium.

Application of distance-dependent resolution compensation and post-reconstruction filtering for myocardial SPECT

Compensation for distance-dependent resolution can be directly incorporated in maximum likelihood reconstruction. Our objective was to examine the effectiveness of this compensation using either the standard expectation maximization (EM) algorithm or an accelerated algorithm based on use of ordered subsets (OSEM). We also investigated the application of post-reconstruction filtering in combination with resolution compensation. Using the MCAT phantom, projections were simulated for 360 degrees data, including attenuation and distance-dependent resolution. Projection data were reconstructed using conventional EM and OSEM with subset size 2 and 4, with/without 3D compensation for detector response (CDR). Also post-reconstruction filtering (PRF) was performed using a 3D Butterworth filter of order 5 with various cutoff frequencies (0.2-1.2 cycles cm(-1)). Image quality and reconstruction accuracy were improved when CDR was included. Image noise was lower with CDR for a given iteration number. PRF with cutoff frequency greater than 0.6 cycles cm(-1) improved noise with no reduction in recovery coefficient for myocardium but the effect was less when CDR was incorporated in the reconstruction. CDR alone provided better results than use of PRF without CDR. Results suggest that using CDR without PRF, and stopping at a small number of iterations, may provide sufficiently good results for myocardial SPECT. Similar behaviour was demonstrated for OSEM.

Quantitative assessment of regional myocardial blood flow with thallium-201 and SPECT

Thallium-201 has been used extensively as a myocardial perfusion agent and to assess myocardial viability. Unlike other 99mTc-labeled agents such as 99mTc-sestamibi and 99mTc-tetrofosmine, the regional concentration of 201Tl varies with time, and its kinetics make it a potential candidate for estimating absolute physiologic parameters with kinetic model analysis. This article outlines a strategy for quantitative assessment of regional myocardial blood flow in man using 201Tl and dynamic single photon emission computed tomography (SPECT). Quantitatively accurate SPECT images that are proportional to the true radioactivity distribution are prerequisites for model-based kinetic analysis. Our technique for quantitative SPECT includes ordered-subset maximum likelihood-expectation maximization (ML-EM) reconstruction with transmission data-based attenuation correction and transmission-dependent convolution subtraction scatter correction. A three-compartment model was found to reproduce the observed regional time-activity curves well, and dog experiments demonstrated that influx rate constant (K1) values estimated from the dynamic SPECT data correlated well with absolute myocardial blood flow determined by in vitro microspheres for a physiologically wide range of flows. Several possible strategies for simplifying the study procedures, without compromising accuracy, are also presented, which should make absolute quantitation of regional myocardial blood flow feasible using 201Tl and a conventional SPECT camera in a clinical setting.

Comparison of frequency-distance relationship and Gaussian-diffusion-based methods of compensation for distance-dependent spatial resolution in SPECT imaging

The goal of this investigation was to compare resolution recovery versus noise level of two methods for compensation of distance-dependent resolution (DDR) in SPECT imaging. The two methods of compensation were restoration filtering based on the frequency-distance relationship (FDR) prior to iterative reconstruction, and modelling DDR in the projector/backprojector pair employed in iterative reconstruction. FDR restoration filtering was computationally faster than modelling the detector response in iterative reconstruction. Using Gaussian diffusion to model the detector response in iterative reconstruction sped up the process by a factor of 2.5 over frequency domain filtering in the projector/backprojector pair. Gaussian diffusion modelling resulted in a better resolution versus noise tradeoff than either FDR restoration filtering or solely modelling attenuation in the projector/backprojector pair of iterative reconstruction. For the pixel size investigated herein (0.317 cm), accounting for DDR in the projector/backprojector pair by Gaussian diffusion, or by applying a blurring function based on the distance from the face of the collimator at each distance, resulted in very similar resolution recovery and slice noise level.