Hybrid scatter correction applied to quantitative holmium-166 SPECT

Ho-166 is a combined beta-gamma emitter of which the betas can be used therapeutically. From the 81 keV gammas of Ho-166, SPECT images can be obtained, which give opportunities to guide Ho-166 therapy. Accurate reconstruction of Ho-166 images is currently hampered by photopeak-scatter in the patient, down-scatter in the detector, collimator and patient caused by the 1.4 MeV photons and by bremsstrahlung. We developed and validated a method for quantitative SPECT of Ho-166 that involves correction for both types of scatter plus non-uniform attenuation correction using attenuation maps. Photopeak-scatter (S) is compensated for by a rapid 3D Monte Carlo (MC) method that is incorporated in ordered subset (OS) reconstruction of the emission data, together with simultaneous correction for attenuation (A) and detector response (D); this method is referred to as OS-ADS. Additionally, for correction of down-scatter, we use a 14 keV wide energy window centred at 118 keV (OS-ADSS). Due to a limited number of available energy windows, the same 118 keV energy window was used for down-scatter correction of the simultaneously acquired Gd-153 transmission data. Validations were performed using physical phantom experiments carried out on a dual-head SPECT system; Gd-153 transmission line sources were used for acquiring attenuation maps. For quantitative comparison of OS-ADS and OS-ADSS, bottles filled with Ho-166 were placed in both a cylindrical phantom and an anthropomorphic thorax phantom. Both OS-ADS and OS-ADSS were compared with an ordered subset reconstruction without any scatter correction (OS-AD). Underestimations of about 20% in the attenuation map were reduced to a few per cent after down-scatter correction. The average deviation from the true activity contained in the bottles was +72% with OS-AD. Using OS-ADS, this average overestimation was reduced to +28% and with OS-ADSS the deviation was further reduced to 16%. With OS-AD and OS-ADS, these numbers were more sensitive to the choice of volumes of interest than with OS-ADSS. For the reconstructed activity distributions, erroneous background activity found with OS-AD was reduced by a factor of approximately 2 by applying OS-ADS and reduced by a factor of approximately 4 by applying OS-ADSS. The combined attenuation, photopeak-scatter and down-scatter correction framework proposed here greatly enhanced the quantitative accuracy of Ho-166 imaging, which is of the uppermost importance for image-guided therapies. It is expected that the method, with adapted window settings, also can be applied to other isotopes with high energy peaks that contaminate the photopeak data, such as I-131 or In-111.     

Some physical factors influencing the accuracy of convolution scatter correction in SPECT

Some important physical factors influencing the accuracy of convolution scatter correction techniques in SPECT are presented. In these techniques scatter correction in the projection relies on filter functions, QF, evaluated by Fourier transforms, from measured scatter functions, Qp, obtained from point spread functions. The spatial resolution has a marginal effect on Qp. Thus a single QF can be used in the scatter correction of SPECT measurements acquired with the low energy high resolution or the low energy general purpose collimators and over a wide range of patient-collimator distances. However, it is necessary to examine the details of the shape of point spread functions during evaluation of Qp. QF is completely described by scatter amplitude AF, slope BF and filter sum SF. SF is obtained by summation of the values of QF occupying a 31 x 31 pixels matrix. Regardless of differences in amplitude and slope, two filter functions are shown to be equivalent in terms of scatter correction ability, whenever their sums are equal. On the basis of filter sum, the observed small influence of ellipticity on QF implies that an average function can be used in scatter correcting SPECT measurements conducted with elliptic objects. SF is shown to increase with a decrease in photon energy and with an increase in window size. Thus, scatter correction by convolution may be severely hampered by photon statistics when SPECT imaging is done with low-energy photons. It is pointless to use unnecessarily large discriminator windows, in the hope of improving photon statistics, since most of the extra events acquired will eventually be subtracted during scatter correction. Regardless of the observed moderate reduction in SF when a lung-equivalent material replaces a portion of a water phantom, further studies are needed to develop a technique that is capable of handling attenuation and scatter corrections simultaneously. Whenever superficial and inner radioactive distributions coexist the observed reduction of SF close to the phantom surface indicates that scatter correction of such distributions has to rely on two distinct filter functions. Corrections based on a surface function produce accurate results in the superficial region, while the central distributions are substantially overestimated. Surface radioactive distributions introduce appreciable errors in the determination of central distributions when corrections are based on central filter function. This function introduces a reduction of about 40% in the measured surface concentration.     

A spatially varying compton scatter correction for SPECT utilizing the integral Klein-Nishina cross section

An algorithm correcting for the fraction of scattered events in SPECT and planar imaging is proposed. The algorithm utilizes a pixel-based multi-channel analyser for data acquisition. The method was designed to operate on a local level by three subtraction steps: (a) Subtracting a modified Klein-Nishina single scatter distribution, pixel by pixel, from the events obtained experimentally in the upper half of the photo-peak window. (b) Subtracting a mirrored distribution of the unscattered events hence obtained from that of total events in the lower half of the window, thus giving the scatter distribution in this part of the window. (c) Subtracting the sum of the scatter distributions in both window halves from the corresponding sum of total events in order to obtain the unscattered photons within the photo-peak window. The accuracy of the method was validated experimentally, using a new rCBF phantom allowing for imaging in matter corresponding to soft tissue and approximately in air, respectively. After correction for photon scattering and attenuation, the regional difference in SPECT values in soft tissue equivalent matter and in low-density matter (simulating air) was only 1.5 +/- 7.2% (mean +/- 1 SD), thus indicating a high accuracy of the correction method. Provided that an accurate and stable pixel peak-alignment routine is available, the method can be applied using a minimum of three windows.     

A three-dimensional ray-driven attenuation, scatter and geometric response correction technique for SPECT in inhomogeneous media

The qualitative and quantitative accuracy of SPECT images is degraded by physical factors of attenuation, Compton scatter and spatially varying collimator geometric response. This paper presents a 3D ray-tracing technique for modelling attenuation, scatter and geometric response for SPECT imaging in an inhomogeneous attenuating medium. The model is incorporated into a three-dimensional projector-backprojector and used with the maximum-likelihood expectation-maximization algorithm for reconstruction of parallel-beam data. A transmission map is used to define the inhomogeneous attenuating and scattering object being imaged. The attenuation map defines the probability of photon attenuation between the source and the scattering site, the scattering angle at the scattering site and the probability of attenuation of the scattered photon between the scattering site and the detector. The probability of a photon being scattered through a given angle and being detected in the emission energy window is approximated using a Gaussian function. The parameters of this Gaussian function are determined using physical measurements of parallel-beam scatter line spread functions from a non-uniformly attenuating phantom. The 3D ray-tracing scatter projector-backprojector produces the scatter and primary components. Then, a 3D ray-tracing projector-backprojector is used to model the geometric response of the collimator. From Monte Carlo and physical phantom experiments, it is shown that the best results are obtained by simultaneously correcting attenuation, scatter and geometric response, compared with results obtained with only one or two of the three corrections. It is also shown that a 3D scatter model is more accurate than a 2D model. A transmission map is useful for obtaining measurements of attenuation and scatter in SPECT data, which can be used together with a model of the geometric response of the collimator to obtain corrected images with quantitative and diagnostically accurate information.     

Improved scatter correction using adaptive scatter kernel superposition

Accurate scatter correction is required to produce high-quality reconstructions of x-ray cone-beam computed tomography (CBCT) scans. This paper describes new scatter kernel superposition (SKS) algorithms for deconvolving scatter from projection data. The algorithms are designed to improve upon the conventional approach whose accuracy is limited by the use of symmetric kernels that characterize the scatter properties of uniform slabs. To model scatter transport in more realistic objects, nonstationary kernels, whose shapes adapt to local thickness variations in the projection data, are proposed. Two methods are introduced: (1) adaptive scatter kernel superposition (ASKS) requiring spatial domain convolutions and (2) fast adaptive scatter kernel superposition (fASKS) where, through a linearity approximation, convolution is efficiently performed in Fourier space. The conventional SKS algorithm, ASKS, and fASKS, were tested with Monte Carlo simulations and with phantom data acquired on a table-top CBCT system matching the Varian On-Board Imager (OBI). All three models accounted for scatter point-spread broadening due to object thickening, object edge effects, detector scatter properties and an anti-scatter grid. Hounsfield unit (HU) errors in reconstructions of a large pelvis phantom with a measured maximum scatter-to-primary ratio over 200% were reduced from -90 ± 58 HU (mean ± standard deviation) with no scatter correction to 53 ± 82 HU with SKS, to 19 ± 25 HU with fASKS and to 13 ± 21 HU with ASKS. HU accuracies and measured contrast were similarly improved in reconstructions of a body-sized elliptical Catphan phantom. The results show that the adaptive SKS methods offer significant advantages over the conventional scatter deconvolution technique.     

Dual-window scatter correction and energy window setting in cerebral blood flow SPECT: a Monte Carlo study

The image quality in SPECT studies of the regional cerebral blood flow (rCBF) performed with 99mTc-HMPAO is degraded by scattered photons. The finite energy resolution of the gamma camera makes the detection of scattered photons unavoidable, and this is observed in the image as an impaired contrast between grey and white matter structures. In this work, a Monte Carlo simulated SPECT study of a realistic voxel-based brain phantom was used to evaluate the resulting contrast-to-noise ratio for a number of energy window settings, with and without the dual-window scatter correction. Values of the scaling factor k, used to obtain the fraction of scattered photons in the photopeak window, were estimated for each energy window. The use of a narrower, asymmetric, energy discrimination window improved the contrast, with a subsequent increase in statistical noise due to the lower number of counts. The photopeak-window setting giving the best contrast-to-noise ratio was found to be the same whether or not scatter correction was applied. Its value was 17% centred at 142 keV. At the optimum photopeak-window setting, the contrast was improved by using scatter correction, but the contrast-to-noise ratio was made worse.     

Implementation of an iterative scatter correction, the influence of attenuation map quality and their effect on absolute quantitation in SPECT

We investigated the accuracy of qSPECT, a quantitative SPECT reconstruction algorithm we have developed which employs corrections for collimator blurring, photon attenuation and scatter, and provides images in units of absolute radiotracer concentrations (kBq cm(-3)). Using simulated and experimental phantom data with characteristics similar to clinical cardiac perfusion data, we studied the implementation of a scatter correction (SC) as part of an iterative reconstruction protocol. Additionally, with experimental phantom studies we examined the influence of CT-based attenuation maps, relative to those obtained from conventional SPECT transmission scans, on SCs and quantitation. Our results indicate that the qSPECT estimated scatter corrections did not change appreciably after the third iteration of the reconstruction. For the simulated data, qSPECT concentrations agreed with images reconstructed using ideal, scatter-free, simulated data to within 6%. For the experimental data, we observed small systematic differences in the scatter fractions for data using different combinations of SCs and attenuation maps. The SCs were found to be significantly influenced by errors in image coregistration. The reconstructed concentrations using CT-based corrections were more quantitatively accurate than those using attenuation maps from conventional SPECT transmission scans. However, segmenting the attenuation maps from SPECT transmission scans could provide sufficient accuracy for most applications.     

Acceleration of Monte Carlo-based scatter compensation for cardiac SPECT

Single proton emission computed tomography (SPECT) images are degraded by photon scatter making scatter compensation essential for accurate reconstruction. Reconstruction-based scatter compensation with Monte Carlo (MC) modelling of scatter shows promise for accurate scatter correction, but it is normally hampered by long computation times. The aim of this work was to accelerate the MC-based scatter compensation using coarse grid and intermittent scatter modelling. The acceleration methods were compared to un-accelerated implementation using MC-simulated projection data of the mathematical cardiac torso (MCAT) phantom modelling (99m)Tc uptake and clinical myocardial perfusion studies. The results showed that when combined the acceleration methods reduced the reconstruction time for 10 ordered subset expectation maximization (OS-EM) iterations from 56 to 11 min without a significant reduction in image quality indicating that the coarse grid and intermittent scatter modelling are suitable for MC-based scatter compensation in cardiac SPECT.     

Generalized dual-energy-window scatter compensation in spatially varying media for SPECT

The detection of scattered photons in the photopeak energy window hinders accurate activity estimation in single-photon-emission computed tomography (SPECT). To compensate for photons scattered in spatially varying media, a framework for generalized dual-energy-window scatter subtraction has been developed. Generalized scatter subtraction factors are introduced, and these factors are decomposed into terms dependent on the uniform (average) and spatially varying components of the source activity distribution. The variation of these factors with projection pixel location and gamma camera position is analysed for a simulated myocardial perfusion study with a 99Tc(m) source radionuclide and a non-uniform thorax model. Monte Carlo methods are used to model photon transport and detection. The application of pixel-dependent scatter subtraction factors for scatter compensation is evaluated in an image reconstruction experiment for this simulated myocardial perfusion study. Generalized matrix inverses with noise-dependent regularization are used for image reconstruction. For this simulation, use of a pixel-dependent scatter subtraction factor and a constant scatter subtraction factor are effective for scatter compensation. Activity estimates within the left ventricular myocardium for these two methods are practically the same as those obtained from image reconstructions where the detection of Compton-scattered photons is included in the system matrix.     

Scatter-to-primary based scatter fractions for transmission-dependent convolution subtraction of SPECT images

In single photon emission computed tomography (SPECT), transmission-dependent convolution subtraction has been shown to be useful when correcting for scattered events. The method is based on convolution subtraction, but includes a matrix of scatter fractions instead of a global scatter fraction. The method can be extended to iteratively improve the scatter estimate, but in this note we show that this requires a modification of the theory to use scatter-to-total scatter fractions for the first iteration only and scatter-to-primary fractions thereafter. To demonstrate this, scatter correction is performed on a Monte Carlo simulated image of a point source of activity in water. The modification of the theory is compared to corrections where the scatter fractions are based on the scatter-to-total ratio, using one and ten iterations. The resulting ratios of subtracted to original counts are compared to the true scatter-to-total ratio of the simulation and the most accurate result is found for our modification of the theory.     

X-ray scatter correction algorithm for cone beam CT imaging

Developing and optimizing an x-ray scatter control and reduction technique is one of the major challenges for cone beam computed tomography (CBCT) because CBCT will be much less immune to scatter than fan-beam CT. X-ray scatter reduces image contrast, increases image noise and introduces reconstruction error into CBCT. To reduce scatter interference, a practical algorithm that is based upon the beam stop array technique and image sequence processing has been developed on a flat panel detector-based CBCT prototype scanner. This paper presents a beam stop array-based scatter correction algorithm and the evaluation results through phantom studies. The results indicate that the beam stop array-based scatter correction algorithm is practical and effective to reduce and correct x-ray scatter for a CBCT imaging task.     

A convolution-subtraction scatter correction method for 3D PET

3D acquisition and reconstruction in positron emission tomography (PET) produce data with improved signal-to-noise ratios compared with conventional 2D slice-oriented methods. However, the sensitivity increase is accompanied by an increase in the number of scattered photons and random coincidences detected. This paper presents a scatter correction technique for 3D PET data where an estimate of the scattered photon distribution is subtracted from the data before reconstruction. The scatter distribution is estimated by iteratively convolving the photopeak projections with a mono-exponential kernel. The method accounts for the 3D acquisition geometry and nature of scatter by performing the scatter estimation on 2D projections. The assumptions of the method have been investigated by measuring the variation in the scatter fraction and the scatter function at different positions in a cylinder. Both parameters were found to vary by up to 50% from the centre to the edge of a large water-filled cylinder. Despite this, in a uniform cylinder containing water with different concentrations of radioactivity, scatter was reduced from 25% in a non-radioactive region to less than 5% using the convolution-subtraction method. In addition, the relative concentration of a cylinder containing an increased concentration, which was underestimated by almost 50% without scatter correction, was within 5% of the true concentration after correction.     

Review: Bioresorbable polymeric stents: current status and future promise

Metal stents and, more recently, polymer-coated metal stents are used to stabilize dissections, eliminate vessel recoil, and guide remodeling after balloon angioplasty and other treatments for arterial disease. Bioresorbable polymeric stents are being developed to improve the biocompatibility and the drug reservoir capacity of metal stents, and to offer a transient alternative to the permanent metallic stent implant. Following a brief review of metal stent technology, the emerging class of expandable, bioresorbable polymeric stents is described, with emphasis on developments in the authors' laboratory.     

Text mining of cancer-related information: Review of current status and future directions

PurposeThis paper reviews the research literature on text mining (TM) with the aim to find out (1) which cancer domains have been the subject of TM efforts, (2) which knowledge resources can support TM of cancer-related information and (3) to what extent systems that rely on knowledge and computational methods can convert text data into useful clinical information. These questions were used to determine the current state of the art in this particular strand of TM and suggest future directions in TM development to support cancer research.MethodsA review of the research on TM of cancer-related information was carried out. A literature search was conducted on the Medline database as well as IEEE Xplore and ACM digital libraries to address the interdisciplinary nature of such research. The search results were supplemented with the literature identified through Google Scholar.ResultsA range of studies have proven the feasibility of TM for extracting structured information from clinical narratives such as those found in pathology or radiology reports. In this article, we provide a critical overview of the current state of the art for TM related to cancer. The review highlighted a strong bias towards symbolic methods, e.g. named entity recognition (NER) based on dictionary lookup and information extraction (IE) relying on pattern matching. The F-measure of NER ranges between 80% and 90%, while that of IE for simple tasks is in the high 90s. To further improve the performance, TM approaches need to deal effectively with idiosyncrasies of the clinical sublanguage such as non-standard abbreviations as well as a high degree of spelling and grammatical errors. This requires a shift from rule-based methods to machine learning following the success of similar trends in biological applications of TM. Machine learning approaches require large training datasets, but clinical narratives are not readily available for TM research due to privacy and confidentiality concerns. This issue remains the main bottleneck for progress in this area. In addition, there is a need for a comprehensive cancer ontology that would enable semantic representation of textual information found in narrative reports.     

A general framework and review of scatter correction methods in cone beam CT. Part 2: scatter estimation approaches

The main components of scatter correction procedures are scatter estimation and a scatter compensation algorithm. This paper completes a previous paper where a general framework for scatter compensation was presented under the prerequisite that a scatter estimation method is already available. In the current paper, the authors give a systematic review of the variety of scatter estimation approaches. Scatter estimation methods are based on measurements, mathematical-physical models, or combinations of both. For completeness they present an overview of measurement-based methods, but the main topic is the theoretically more demanding models, as analytical, Monte-Carlo, and hybrid models. Further classifications are 3D image-based and 2D projection-based approaches. The authors present a system-theoretic framework, which allows to proceed top-down from a general 3D formulation, by successive approximations, to efficient 2D approaches. A widely useful method is the beam-scatter-kernel superposition approach. Together with the review of standard methods, the authors discuss their limitations and how to take into account the issues of object dependency, spatial variance, deformation of scatter kernels, external and internal absorbers. Open questions for further investigations are indicated. Finally, the authors refer on some special issues and applications, such as bow-tie filter, offset detector, truncated data, and dual-source CT.     

Scatter correction in SPECT using non-uniform attenuation data

Quantitative assessment of activity levels with SPECT is difficult because of attenuation and scattering of gamma rays within the object. To study the effect of attenuation and scatter on SPECT quantitation, phantom studies were performed with non-uniform attenuation. Simulated transmission CT data provided information about the distribution of attenuation coefficients within the source. Attenuation correction was performed by an iterative reprojection technique. Scatter correction was done by convolution of the attenuation-corrected image and an appropriate filter. The filter characteristics depended on the attenuation and activity measurement at each pixel. The scatter correction could compensate completely for the 28% scatter component from a line source, and the 61% component from a thick, extended source. Accuracy of regional activity ratios and the linearity of the relationship between true radioactivity and the SPECT measurement were both significantly improved by these corrections. The present method is expected to be valuable for the quantitative assessment of regional activity.