Rigid-body transformation of list-mode projection data for respiratory motion correction in cardiac PET

High-resolution cardiac PET imaging with emphasis on quantification would benefit from eliminating the problem of respiratory movement during data acquisition. Respiratory gating on the basis of list-mode data has been employed previously as one approach to reduce motion effects. However, it results in poor count statistics with degradation of image quality. This work reports on the implementation of a technique to correct for respiratory motion in the area of the heart at no extra cost for count statistics and with the potential to maintain ECG gating, based on rigid-body transformations on list-mode data event-by-event. A motion-corrected data set is obtained by assigning, after pre-correction for detector efficiency and photon attenuation, individual lines-of-response to new detector pairs with consideration of respiratory motion. Parameters of respiratory motion are obtained from a series of gated image sets by means of image registration. Respiration is recorded simultaneously with the list-mode data using an inductive respiration monitor with an elasticized belt at chest level. The accuracy of the technique was assessed with point-source data showing a good correlation between measured and true transformations. The technique was applied on phantom data with simulated respiratory motion, showing successful recovery of tracer distribution and contrast on the motion-corrected images, and on patient data with C15O and 18FDG. Quantitative assessment of preliminary C15O patient data showed improvement in the recovery coefficient at the centre of the left ventricle.     

A method for incorporating organ motion due to breathing into 3D dose calculations.

A method is proposed that incorporates the effects of intratreatment organ motion due to breathing on the dose calculations for the treatment of liver disease. Our method is based on the convolution of a static dose distribution with a probability distribution function (PDF) which describes the nature of the motion. The organ motion due to breathing is assumed here to be one-dimensional (in the superior-inferior direction), and is modeled using a periodic but asymmetric function (more time spent at exhale versus inhale). The dose distribution calculated using convolution-based methods is compared to the static dose distribution using dose difference displays and the effective volume (Veff) of the uninvolved liver, as per a liver dose escalation protocol in use at our institution. The convolution-based calculation is also compared to direct simulations that model individual fractions of a treatment. Analysis shows that incorporation of the organ motion could lead to changes in the dose prescribed for a treatment based on the Veff of the uninvolved liver. Comparison of convolution-based calculations and direct simulation of various worst-case scenarios indicates that a single convolution-based calculation is sufficient to predict the dose distribution for the example treatment plan given.     

A knowledge-based method for reducing attenuation artefacts caused by cardiac appliances in myocardial PET/CT

Attenuation artefacts due to implanted cardiac defibrillator leads have previously been shown to adversely impact cardiac PET/CT imaging. In this study, the severity of the problem is characterized, and an image-based method is described which reduces the resulting artefact in PET. Automatic implantable cardioverter defibrillator (AICD) leads cause a moving-metal artefact in the CT sections from which the PET attenuation correction factors (ACFs) are derived. Fluoroscopic cine images were measured to demonstrate that the defibrillator's highly attenuating distal shocking coil moves rhythmically across distances on the order of 1 cm. Rhythmic motion of this magnitude was created in a phantom with a moving defibrillator lead. A CT study of the phantom showed that the artefact contained regions of incorrect, very high CT values and adjacent regions of incorrect, very low CT values. The study also showed that motion made the artefact more severe. A knowledge-based metal artefact reduction method (MAR) is described that reduces the magnitude of the error in the CT images, without use of the corrupted sinograms. The method modifies the corrupted image through a sequence of artefact detection procedures, morphological operations, adjustments of CT values and three-dimensional filtering. The method treats bone the same as metal. The artefact reduction method is shown to run in a few seconds, and is validated by applying it to a series of phantom studies in which reconstructed PET tracer distribution values are wrong by as much as 60% in regions near the CT artefact when MAR is not applied, but the errors are reduced to about 10% of expected values when MAR is applied. MAR changes PET image values by a few per cent in regions not close to the artefact. The changes can be larger in the vicinity of bone. In patient studies, the PET reconstruction without MAR sometimes results in anomalously high values in the infero-septal wall. Clinical performance of MAR is assessed by two physicians' inspection of images generated in 30 patients with and without MAR. Noticeable image differences are judged in 14 of 28 (50%) observations with AICD leads, and significant clinical impact is judged in 2 of 28 (7%) of those observations. A polar map analysis shows significant differences in 10 of 14 (71%) studies with AICD leads, and 0 of 16 (0%) studies without AICD leads. These results show that the MAR method is successful in reducing the magnitude of the metal artefact without incorrectly altering cases without metal artefact. In spite of profound changes to the CT image from the moving metal, the PET ACF in that study was changed by no more than 20%.     

A spline-regularized minimal residual algorithm for iterative attenuation correction in SPECT.

In SPECT, regularization is necessary to avoid divergence of the iterative algorithms used for non uniform attenuation compensation. In this paper, we propose a spline-based regularization method for the minimal residual algorithm. First, the acquisition noise is filtered using a statistical model involving spline smoothing so that the filtered projections belong to a Sobolev space with specific continuity and derivability properties. Then, during the iterative reconstruction procedure, the continuity of the inverse Radon transform between Sobolev spaces is used to design a spline-regularized filtered backprojection method, by which the known regularity properties of the projections determine those of the corresponding reconstructed slices. This ensures that the activity distributions estimated at each iteration present regularity properties, which avoids computational noise amplification, thus stabilizing the iterative process. Analytical and Monte Carlo simulations are used to show that the proposed spline-regularized minimal residual algorithm converges to a satisfactory stable solution in terms of restored activity and homogeneity, using at most 25 iterations, whereas the non regularized version of the algorithm diverges. Choosing the number of iterations is therefore no longer a critical issue for this reconstruction procedure.     

A comparison of MR-based attenuation correction in PET versus SPECT

Attenuation correction (AC) is a critical step in the reconstruction of quantitatively accurate positron emission tomography (PET) and single photon emission computed tomography (SPECT) images. Several groups have proposed magnetic resonance (MR)-based AC algorithms for application in hybrid PET/MR systems. However, none of these approaches have been tested on SPECT data. Since SPECT/MR systems are under active development, it is important to ascertain whether MR-based AC algorithms validated for PET can be applied to SPECT. To investigate this issue, two imaging experiments were performed: one with an anthropomorphic chest phantom and one with two groups of canines. Both groups of canines were imaged from neck to abdomen, one with PET/CT and MR (n = 4) and the other with SPECT/CT and MR (n = 4), while the phantom was imaged with all modalities. The quality of the nuclear medicine reconstructions using MR-based attenuation maps was compared between PET and SPECT on global and local scales. In addition, the sensitivity of these reconstructions to variations in the attenuation map was ascertained. On both scales, it was found that the SPECT reconstructions were of higher fidelity than the PET reconstructions. Further, they were less sensitive to changes to the MR-based attenuation map. Thus, MR-based AC algorithms that have been designed for PET/MR can be expected to demonstrate improved performance when used for SPECT/MR.     

Calculation and validation of the use of effective attenuation coefficient for attenuation correction in In-111 SPECT

Nuclear medicine tracers using 111In as a radiolabel are increasing in their use, especially in the domain of oncologic imaging. In these applications, it often is critical to have the capability of quantifying radionuclide uptake and being able to relate it to the biological properties of the tumor. However, images from single photon emission computed tomography (SPECT) can be degraded by photon attenuation, photon scattering, and collimator blurring; without compensation for these effects, image quality can be degraded, and accurate and precise quantification is impossible. Although attenuation correction for SPECT is becoming more common, most implementations can only model single energy radionuclides such as 99mTc and 123I. Thus, attenuation correction for 111In is challenging because it emits two photons (171 and 245 keV) at nearly equal rates (90.2% and 94% emission probabilities). In this paper, we present a method of calculating a single "effective" attenuation coefficient for the dual-energy emissions of 111In, and that can be used to correct for photon attenuation in radionuclide images acquired with this radionuclide. Using this methodology, we can derive an effective linear attenuation coefficient Micro(eff) and an effective photon energy E(eff) based on the emission probabilities and linear attenuation coefficients of the 111In photons. This approach allows us to treat the emissions from 111In as a single photon with an effective energy of 210 keV. We obtained emission projection data from a tank filled with a uniform solution of 111In. The projection data were reconstructed using an iterative maximum-likelihood algorithm with no attenuation correction, and with attenuation correction assuming photon energies of 171, 245, and 210 keV (the derived E(eff)). The reconstructed tomographic images demonstrate that the use of no attenuation correction, or correction assuming photon energies of 171 or 245 keV introduces inaccuracies into the reconstructed radioactivity distribution when compared against the effective energy method. In summary, this work provides both a theoretical framework and experimental methodology of attenuation correction for the dual-energy emissions from 111In. Although these results are specific to 111In, the foundation could easily be extended to other multiple-energy isotopes.     

Impact of attenuation correction strategies on the quantification of High Resolution Research Tomograph PET studies

In this study, the quantitative accuracy of different attenuation correction strategies presently available for the High Resolution Research Tomograph (HRRT) was investigated. These attenuation correction methods differ in reconstruction and processing (segmentation) algorithms used for generating a micro-image from measured 2D transmission scans, an intermediate step in the generation of 3D attenuation correction factors. Available methods are maximum-a-posteriori reconstruction (MAP-TR), unweighted OSEM (UW-OSEM) and NEC-TR, which transforms sinogram values back to their noise equivalent counts (NEC) to restore Poisson distribution. All methods can be applied with or without micro-image segmentation. However, for MAP-TR a micro-histogram is a prior during reconstruction. All possible strategies were evaluated using phantoms of various sizes, simulating preclinical and clinical situations. Furthermore, effects of emission contamination of the transmission scan on the accuracy of various attenuation correction strategies were studied. Finally, the accuracy of various attenuation corrections strategies and its relative impact on the reconstructed activity concentration (AC) were evaluated using small animal and human brain studies. For small structures, MAP-TR with human brain priors showed smaller differences in micro-values for transmission scans with and without emission contamination (<8%) than the other methods (<26%). In addition, it showed best agreement with true AC (deviation <4.5%). A specific prior designed to take into account the presence of small animal fixation devices only very slightly improved AC precision to 4.3%. All methods scaled micro-values of a large homogeneous phantom to within 4% of the water peak, but MAP-TR provided most accurate AC after reconstruction. However, for clinical data MAP-TR using the default prior settings overestimated the thickness of the skull, resulting in overestimations of micro-values in regions near the skull and thus in incorrect AC for cortical regions. Using NEC-TR with segmentation or MAP-TR with an adjusted human brain prior showed less overestimation in both skull thickness and AC for these structures and are therefore the recommended methods for human brain studies.     

Optimization of attenuation correction for positron emission tomography studies of thorax and pelvis using count-based transmission scans

The quality of thorax and pelvis transmission scans and therefore of attenuation correction in PET depends on patient thickness and transmission rod source strength. The purpose of the present study was to assess the feasibility of using count-based transmission scans, thereby guaranteeing more consistent image quality and more precise quantification than with fixed transmission scan duration. First, the relation between noise equivalent counts (NEC) of 10 min calibration transmission scans and rod source activity was determined over a period of 1.5 years. Second, the relation between transmission scan counts and uniform phantom diameter was studied numerically, determining the relative contribution of counts from lines of response passing through the phantom as compared with the total number of counts. Finally, the relation between patient weight and transmission scan duration was determined for 35 patients, who were scanned at the level of thorax or pelvis. After installation of new rod sources, the NEC of transmission scans first increased slightly (5%) with decreasing rod source activity and after 3 months decreased with a rate of 2-3% per month. The numerical simulation showed that the number of transmission scan counts from lines of response passing through the phantom increased with phantom diameter up to 7 cm. For phantoms larger than 7 cm, the number of these counts decreased at approximately the same rate as the total number of transmission scan counts. Patient data confirmed that the total number of transmission scan counts decreased with increasing patient weight with about 0.5% kg(-1). It can be concluded that count-based transmission scans compensate for radioactive decay of the rod sources. With count-based transmission scans, rod sources can be used for up to 1.5 years at the cost of a 50% increased transmission scan duration. For phantoms with diameters of more than 7 cm and for patients scanned at the level of thorax or pelvis, use of count-based transmission scans is feasible and results in statistically more consistent transmission scans as compared with fixed transmission scan duration.     

A comparison of attenuation correction methods for quantitative single photon emission computed tomography

Transverse section tomograms of experimental phantoms and patients have been obtained using a GE 400T camera and a filtered back-projection reconstruction technique. These tomograms have been compared with the corresponding sections reconstructed from the same tomographic projection data, but using iterative algorithms with correction for photon attenuation. The comparison assesses the importance of including a correction for attenuation as well as demonstrating how closely a simple geometric attenuation correction, applied to the filtered back-projection reconstruction method, approximates to a more accurate correction incorporated in the computation of line integrals during iterative reconstruction. A comparison is also made between the behaviour of reconstruction algorithms with simulated projection data and real data in terms of convergence properties, and some shortcomings arising from simulation are noted.     

A novel interpolation strategy for estimating subsample speckle motion

Multidimensional, high-resolution ultrasonic imaging of rapidly moving tissue is primarily limited by sparse sampling in the lateral dimension. In order to achieve acceptable spatial resolution and velocity quantization, interpolation of laterally sampled data is necessary. We present a novel method for estimating lateral subsample speckle motion and compare it with traditional interpolation methods. This method, called grid slopes, requires no a priori knowledge and can be applied to data with as few as two samples in the lateral dimension. Computer simulations were performed to compare grid slopes with two conventional interpolation schemes, parabolic fit and cubic spline. Results of computer simulations show that parabolic fit and cubic spline performed poorly at translations greater than 0.5 samples, and translations less than 0.5 samples were subject to an estimation bias. Grid slopes accurately estimated translations between 0 and 1 samples without estimation bias at high signal-to-noise ratios. Given that the grid slopes interpolation technique performs well at high signal-to-noise ratios, one pertinent clinical application might be tissue motion tracking.