Skip to main content
  • Research Article
  • Open access
  • Published:

Homotopic non-local regularized reconstruction from sparse positron emission tomography measurements

Abstract

Background

Positron emission tomography scanners collect measurements of a patient’s in vivo radiotracer distribution. The system detects pairs of gamma rays emitted indirectly by a positron-emitting radionuclide (tracer), which is introduced into the body on a biologically active molecule, and the tomograms must be reconstructed from projections. The reconstruction of tomograms from the acquired PET data is an inverse problem that requires regularization. The use of tightly packed discrete detector rings, although improves signal-to-noise ratio, are often associated with high costs of positron emission tomography systems. Thus a sparse reconstruction, which would be capable of overcoming the noise effect while allowing for a reduced number of detectors, would have a great deal to offer.

Methods

In this study, we introduce and investigate the potential of a homotopic non-local regularization reconstruction framework for effectively reconstructing positron emission tomograms from such sparse measurements.

Results

Results obtained using the proposed approach are compared with traditional filtered back-projection as well as expectation maximization reconstruction with total variation regularization.

Conclusions

A new reconstruction method was developed for the purpose of improving the quality of positron emission tomography reconstruction from sparse measurements. We illustrate that promising reconstruction performance can be achieved for the proposed approach even at low sampling fractions, which allows for the use of significantly fewer detectors and have the potential to reduce scanner costs.

Peer Review reports

Background

The concept of positron emission tomography (PET) imaging revolves around the measurement of a patient’s in vivo radiotracer distribution. The system detects pairs of gamma rays emitted indirectly by a positron-emitting radionuclide (tracer), which is introduced into the body via biologically active molecules. As such, PET imaging can be described via a line-integral model of acquisition [1]. PET data are collected as projections in sinograms or listmode format [2]. The raw data collected by a PET scanner are a list of ‘coincidence events’ representing near-simultaneous detection, each coincidence event represents a line in space connecting the two detectors along which the positron emission occurred. While one can employ the inverse Radon transform to recover tomograms from acquired PET data, such an approach can be unstable, particularly when dealing with noise-contaminated data [3]. Therefore, the filtered back projection algorithm, which is a stabilized and discretized version of the inverse Radon transform, is commonly used in practice [3].

There are two common approaches for tomographic reconstruction: i) filtered back projection (FBP) [4], and ii) iterative expectation-maximization (EM) [5-7]. FBP has been widely used to reconstruct tomograms from the projections in clinical settings due to its overall simplicity and low computational complexity. However, shot noise in the raw data is prominent in the reconstructed tomograms and areas of high tracer uptake tend to form streaks across the tomogram.

Iterative algorithms for reconstructing tomograms from acquired PET data take into account the statistical nature of the acquired projection data and incorporate the physical model into reconstruction. These algorithms compute an estimate of the likely distribution of annihilation events that led to the measured data, based on statistical principles [8]. The advantage is a better noise profile and resistance to the artifacts common with FBP, but are typically more computationally intensive.

Recently, wavelet-based methods [3,8-10] have also been proposed for solving this inverse problem in PET reconstruction. Wavelet-based methods have been applied in research literature in the post-processing stage to suppress artifacts in the reconstructed tomograms [3]. Another attractive approach to reconstructing tomograms from PET data using wavelets is to incorporate wavelets into the reconstruction process itself [9], where an L1-regularization constraint is used to enforce sparsity in the wavelet domain. For example, wavelet regularization using exponential-spline wavelets has been shown to produce good reconstruction results [10].

More recently, the concept of compressive sensing [11] has offered great potential for signal recovery with very high accuracy through sparse measurements acquired at very low sampling fractions. Researchers have also begun to explore the use of compressive sensing for image reconstruction in various medical modalities [12]. To date, most sparse medical image reconstruction algorithms have focused on solving the sparse reconstruction problem via L1 minimization.

The main contribution of this paper is to introduce and investigate the potential of a homotopic non-local regularization (HNLR) reconstruction framework for the reconstruction of PET tomograms from sparse measurements, with the aim to achieve high reconstruction quality. The ability to achieve strong reconstruction results at low sampling fractions can be very helpful for reducing the number of detectors needed in a PET system [7], which in effect helps can lead to reduced scanner costs and potentially improve adoption of PET in developing countries.

The rest of the paper is organized as follows. First, the underlying methodology behind the proposed use of a homotopic non-local regularization (HNLR) reconstruction framework for the reconstruction of PET tomograms from sparse measurements is described in Section “Methods”. The experimental results are presented in Section “Results and discussion”. Finally, conclusions are drawn and future work is discussed in Section “Conclusion”.

Methods

The underlying methodology behind the proposed approach to reconstructing PET tomograms from sparse measurements can be explained as follows. The underlying imaging system can be represented with the following linear relationship

$$ P(\alpha,s)=\mathbb{R}{f(x)} $$
((1))

where P is the set of sampled measurements, s is the distance of straight line L from the origin, α is the angle that the normal vector to L makes with the x axis, is the Radon transform, and f is the unknown tomogram. The goal of reconstruction is to use the sampled measurements P (projections in the Radon transform domain through the object) to find the tomogram f in image space.

L p minimization for sparse reconstruction

In PET, the number of projections is related to the number of detectors. Therefore, taking a high number of projections requires tightly packed discrete detector rings [7], which could lead to high scanner costs. Instead, our goal is to measure only a fraction of the projections, hence reducing the number of detectors needed, and still be able to get high quality reconstructions. Consequently, Eq. 1 should be rewritten as:

$$ f(x)=\mathbb{R}^{-1} \left \{ \Phi P(\alpha,s) \right \} $$
((2))

where \(\mathbb {R}^{-1}\) is the inverse Radon transform operator and Φ is a measurement operator defining which of the sites are measured (non-measured projection space indices will be assigned 0 by Φ). Our goal is to reconstruct f(x) from a sparse sampling of P(α,s). With only few measurements made, not all sites are measured, making this problem an ill-posed inverse problem [13], and as such there may exist multiple solutions. A common method to solving this inverse problem is through the use of L2 minimization:

$$ \hat{f}(x) \in \underset{f(x)}{arg\, min} \left \|\, f(x) \right \|_{2}\; \; s.t. \; \mathbb{R}{f(x)}=\Phi P(\alpha,s) $$
((3))

where \(\hat {f}(x)\) is the estimated signal in image space, and \(\hat {P}(\alpha,s)\) is the estimated signal in projection space. However, solving the problem in this manner results in significant artifacts such as aliasing and blurring in practical situations [12].

Recently, scientists realized that many signals, such as PET imagery, possess an inherent sparsity in some sparsifying transform domain [14]. According to the emerging theory of compressive sensing [11,15], a better estimate of f(x) can be obtained by maximizing the sparsity of the signal in the transform domain and enforcing data fidelity in projection space domain. One can formulate the reconstruction problem as a constrained L0 minimization problem,

$$ \hat{f}(x) \in \underset{f(x)}{arg\, min} \left \| \Psi f(x) \right \|_{0}\; \; s.t. \; \mathbb{R}{f(x)}=\Phi P(\alpha,s) $$
((4))

where Ψ is the sparse transform operator. Two commonly used sparse transform operators are the finite differential transformation and the wavelet transformation [16].

Unfortunately, solving the L0 problem is essentially NP hard, and as such is intractable in practice [17]. To address this problem, pioneering work done by Candes [11] and Donoho [15] showed that under certain conditions, L0 and L1 minimization can lead to the same solutions. As such, one can instead solve the L1 minimization as:

$$ \hat{f}(x) \in \underset{f(x)}{arg\, min} \left \| \Psi f(x) \right \|_{1}\; \; s.t. \; \mathbb{R}{f(x)}=\Phi P(\alpha,s) $$
((5))

Theoretically, under certain conditions, solving the L1 problem can get exactly the same solution as solving the L0 problem, although a substantial increase in the number of measurements is required [12,18].

Typically, a total variation (TV) penalty [19] is employed in the L1 minimization framework, which is known to have an edge-preservation effect, and to account for the unavoidable noise in the measurements, the data fidelity constraint is typically replaced by a L2 norm constraint:

$$ \begin{aligned} \hat{f}(x) \in \underset{f(x)}{arg\, min} \left \| \Psi f(x) \right \|_{1}+TV\left(f(x)\right)~~~~ \\ s.t. \; \left \| \mathbb{R}{f(x)}-\Phi P(\alpha,s) \right \|_{2}<\epsilon \end{aligned}\vspace*{10pt} $$
((6))

Homotopic, non-local regularization based reconstruction

There are two fundamental limitations associated with the reconstruction framework described in Eq. 6 for the purpose of reconstructing tomograms [20]. First, the use of a TV penalty can lead to the loss of fine structure at low sampling fractions [12]. Since PET images are characterized by complex functional process variations, the TV approach may not be well-suited for reconstructing these data. Second, as already mentioned, the number of measurements required for solving the L1 problem can be noticeably higher than that for the L0 problem, which is undesirable as it leads to an increased number of detectors needed. Hence, an alternative reconstruction strategy that addresses these two problems is desired.

To achieve the theoretical capability of the constrained L0 minimization approach without a drastic increase in the number of measurements, Trzasko et al. [12] introduced a homotopic L0 minimization framework, which uses an increasing approximation framework to get close to the L0 minimization problem:

$$ \begin{aligned} \hat{f}(x) \in {\lim}_{\sigma \to 0} \underset{f(x)}{arg\, min} \; \rho (|\Psi f(x) |,\sigma) \\ s.t. \; \left \| \mathbb{R}{f(x)}-\Phi P(\alpha,s) \right \|_{2}<\epsilon \end{aligned} $$
((7))

where ρ is the homotopic approximation of the L0 norm, which approaches the L0 norm as σ approaches 0. Experimental results showed that this strategy is able to approach the capabilities of the L0 minimization approach [21], thus addressing the problem pertaining to the number of measurements required.

To address the problem of detail loss from using total variation, we instead integrate the concept of non-local regularization, often used in image processing for improved detail preservation [19,20,22-24], into the homotopic L0 minimization framework [25] for the purpose of tomogram reconstruction from sparse PET measurements, as it is well-suited for handling fine structural details. The proposed homotopic, non-local regularization framework for sparse PET reconstruction is formulated as follows:

$$ \begin{aligned} \hat{f}(x) \in {\lim}_{\sigma \to 0} \underset{f(x)}{arg\, min} \; \eta (|\Phi f(x) |,\sigma) \\ s.t. \; \left \| \mathbb{R}{f(x)}-\Phi P(\alpha,s) \right \|_{2}<\epsilon \end{aligned} $$
((8))

where η denotes the homotopic non-local regularization functional:

$$ \eta(|\Phi f(x)|,\sigma)=\left [\sum\limits_{\textit{x}\in\Omega} \sum\limits_{y\in N(x)}w(x,y,\sigma) \cdot (N(x)-N(y))^{2} \right] $$
((9))

σ controls the approximation degree, and w(x,y,σ) is defined as

$$ w\left [ x,y,\sigma \right ]=exp\left \{ -\frac{(N(x)-N(y)}{2\sigma^{2}} \right \} $$
((10))

where N(x) is a neighborhood around x. Essentially, the idea is to minimize the homotopic non-local regularization function, denoted by η(|Φf(x)|,σ), at decreasing values of approximation degree, denoted by σ, such that η approximates the L0 norm closer and closer. An approximate solution of Eq. 8 can be obtained using an iterative optimization approach, with decreasing values of σ as the number of iterations increases. To reduce the computation complexity, the neighborhood search space is limited to a window search around the pixel to be estimated. The regularization term η(|Ψf(x)|,σ) is designed to efficiently suppress artifacts associated with incomplete measurements while preserving detail and structure, and is enforced via steepest descent [26,27], while the data fidelity term is designed to ensure that the reconstructed tomogram complies with sampled measurements while accounting for some level of noise artifacts. The data fidelity term is enforced at measured projection space indices at each iteration i as follows:

$$ {\fontsize{7.8pt}{9.6pt}\selectfont{\begin{aligned} &P_{i}(\alpha,s)\\ &\,=\,\left\{\!\!\begin{array}{ccccc} \mathbb{R}\,{\hat{f_{i}}(x)} &\! \text{if} &\! \left \| \mathbb{R}{\hat{f_{i}}(x)}-\Phi P(\alpha,s) \right \|_{2}\!\!<\epsilon\\ P(\alpha,s)+\left(\!2H\left(\!\mathbb{R}\,{\hat{f_{i}}(x)}-\Phi P(\alpha,s)\!\right)-1\!\right)\frac{\epsilon}{b} & \!\text{if} & otherwise\\ \end{array} \right. \end{aligned}}} $$
((11))

where H(.) is the Heaviside step function and b is the number of measured projection space indices. The pseudocode for the proposed homotopic non-local regularized reconstruction method is shown in Algorithm 1. Note that a L2 data term is used here for simplicity in the iterative optimization realization shown in Algorithm 1, and more advanced data terms such as the Kullback-Leibler divergence [28,29], which takes better advantage of the Poisson noise statistics of PET, may be used.

figure a

Experiments

Data description

To evaluate the effectiveness of the proposed method, sparse reconstructions were performed using a simulated phantom generated using ASIM [30] at 30% sampling fraction (defined here as the percentage of projection angles from which measurements are obtained). The simulated phantom consists of 4 spots of different sizes and is useful for observing the effects of reduced sampling on spatial resolution. Also, since the simulated phantom generated using ASIM is contaminated by simulated Poisson-distributed noise, it is also useful for observing the performance of the proposed method when faced with noisy data. Since the total number of projection angles used in this study is 180 projection angles, that means at a 30% sampling fraction the number of projection angles used for sparse reconstruction is 54 projection angles. Also, the object covers approximately 70% of the field of view. Furthermore, sparse reconstructions at different sampling fractions were performed using three real PET data sets hosted by Harvard Medical School [31], in accordance with Harvard Medical School ethics procedures. The PET data sets are acquired with Fluorine-18 Deoxyglucose (FDG) as the biologically active molecule. The first PET data set (PET1) is acquired from a 73 year-old man who sought medical attention due to a grand mal seizure and progressive speech difficulties, with a brain biopsy revealing grade II astrocytoma. The second PET data set (PET2) was taken from a 53 year-old man who sought medical attention due to a grand mal seizure, with a brain biopsy revealing grade IV astrocytoma. The third PET data set (PET3) is acquired from a 70 year-old man with mild Alzheimer’s disease about 9 months prior to imaging. The PET data sets are 256×256×23 voxels. The data sets were forward-projected from a previous reconstruction obtained via filtered back-projection (FBP) [2] based on measurements from 100% of the projection angles in order to obtain projection data for evaluation.

Implementation

For comparison purposes, the proposed homotopic non-local regularization (HNLR) algorithm was compared with the standard filtered back-projection (FBP) algorithm [2], as well as the more advanced iterative expectation maximization algorithm with total variation regularization (TVR) algorithm [6]. The standard FBP algorithm was chosen for baseline comparisons as it is widely used in clinical settings for its speed and ease of implementation. The TVR algorithm was chosen as it is one of the more advanced reconstruction approaches available for PET reconstruction [6,7]. To allow for quantitative assessment, the signal-to-noise ratio (SNR) was computed for reconstructed data at different sampling fractions. Furthermore, a qualitative visual assessment is performed on the reconstructed data to investigate the reconstruction performance and the preservation of details achieved using the tested methods. Different sampling fractions were achieved by sampling the projection angles at evenly spaced intervals. For example, to achieve a 50% sampling fraction, every other projection angle from the set of total projection angles was sampled. This methodology was used for all the data sets and the comparison of all three methods. For this study, the parameters of the HNLR algorithm were set as follows as it yielded strong reconstruction performance. The neighborhood N(x) was chosen as a 9×9 neighborhood around x, σ1 and σ t are set as 20% and 1% of the dynamic range in image space, respectively, λ is set as 0.8, ε is set as 2% of the dynamic range in projection space multiplied by the number of measured projection space indices (b), and n is set as 100 as these parameters were found empirically to provide strong reconstruction performance. Note that the parameters of the HNLR algorithm can be tuned based on the PET imaging system the algorithm is incorporated into to obtain optimal performance for clinical data obtained from the PET system.

Results and discussion

In order to perform a comprehensive and systematic assessment of the reconstruction performance of the different methods, the peak signal-to-noise ratio (PSNR) was computed for a wide range of sampling fractions. The PSNR metric was computed as follows:

$$ PSNR=10 \cdot {log}_{10}\left (\frac{max(x)^{2}}{MSE} \right) $$
((12))

and (MSE) was defined as mean squared error between original image and reconstructed image:

$$ MSE=\frac{1}{N}\sum\limits_{\textit{x}\in\Omega}\left(f(x)-\hat{f}(x)\right)^{2} $$
((13))

where f(x) is original image, \(\hat {f}(x)\) is reconstructed image, and N is the number of pixels in each image.

Figure 1 shows PET tomograms for the simulated phantom reconstructed at 30% sampling using the three reconstruction methods. We can see very clearly that the FBP method leads to considerable artifacts in the reconstructed tomogram, with significant streaking artifacts over both the background and the phantom details. Both TVR and HNLR methods were able to significantly suppress these artifacts caused due to sampling at sub-Nyquist sampling fraction. It can also be observed that the details (e.g., dark circles) in the tomogram reconstructed using HNLR appear sharper than that reconstructed using TVR. Furthermore, it can be observed that while the reconstructed PET tomograms produced by the reconstruction methods are able to resolve all 4 spots in the phantom, the resolution of the smallest spot noticeably suffers as a result of the reduction in sampling and illustrates the fundamental tradeoff between spatial resolution and sampling fraction. For quantitative evaluation, the PSNR for the PET tomograms reconstructed using FBP, TVR, and HNLR are calculated as 8.03 db, 38.48 dB, and 38.84 dB, respectively.

Figure 1
figure 1

Reconstructed tomograms from simulated phantom using three different reconstruction methods based on 30% sampling. We can see very clearly that the FBP method leads to considerable artifacts in the reconstructed tomogram, with significant streaking artifacts over both the background and the phantom details. Both TVR and HNLR methods were able to significantly suppress these artifacts caused due to sampling at sub-Nyquist sampling fractions. It can also be observed that the details (e.g., dark circles) in the tomogram reconstructed using HNLR appear sharper than that reconstructed using TVR.

The average PSNR of the reconstructed PET data sets is shown in Figure 2, as a function of the sampling fraction (i.e., the percentage of projections sampled used to perform sparse reconstruction). It can be observed that the proposed HNLR method produces reconstructed PET data with higher PSNR values for all levels of sampling fraction when compared to the FBP method and the TVR method. It is important to note that there is no fundamental tradeoff in terms of spatial resolution as the data fidelity in projection space is preserved in the proposed HNLR method.

Figure 2
figure 2

Average SNR of three different reconstruction methods across PET sets, as a function of sampling fraction. The proposed HNLR method produces reconstructed PET data with higher SNR values for all levels of sampling fraction when compared to the FBP method and the TVR method.

Figures 3, 4 and 5 shows PET tomograms from the PET data sets reconstructed at 30% sampling using the three reconstruction methods. We can see very clearly that the FBP method leads to considerable artifacts, particularly streaking artifacts, in both the background and imaged brain areas, which can affect the visibility of certain functional details in the reconstructed data. TVR provides improved reconstruction performance than FBP, in terms of reduced artifacts, but still has noticeable artifacts in the background area and detail loss in the imaged brain area, which can obscure important functional characteristics in the tomograms. Finally, the reconstructed tomogram produced using the proposed HNLR approach exhibits significantly fewer artifacts, particularly streaking artifacts, in the background area compared to FBP and TVR, while providing better functional detail preservation in the imaged brain areas in the reconstructed tomograms than TVR. The artifact reduction by both the TVR and HNLR methods within the brain region is particularly noticeable in Figure 4, where it is very clear that streaking artifacts (caused by sampling at a sub-Nyquist sampling fraction) is prominently overlaid on the entire brain region in the reconstruction produced using the FBP method (see red arrows), while they are significantly reduced in the reconstructions produced by TVR and HNLR as well as an improvement in contrast. It is also important to highlight that, within the same brain region in the HNLR reconstruction, there is an overall improvement in sharpness when compared to the TVR reconstruction (see green arrows). Figure 4 also shows two line profiles from the PET data sets reconstructed using the methods, and further illustrates the reduction in artifacts and improved contrast achieved by TVR and HNLR when compared to FBP.

Figure 3
figure 3

Reconstructed tomograms from the PET1 data set using three different reconstruction methods based on 30% sampling. The PET1 data set was taken from a 73 year-old man who sought medical attention due to a grand mal seizure and progressive speech difficulties, with a brain biopsy revealing grade II astrocytoma. It can be observed that the reconstructed tomogram using the FBP method exhibits the most significant artifacts (both in the background area and the imaged brain area), when compared to TVR and the proposed HNLR methods. The reconstructed tomogram using the TVR method exhibits significantly reduced but still noticeable artifacts in the background area, and no artifacts in the imaged brain area. However, there is noticeable loss in fine detail in the imaged brain area. Finally, the reconstructed tomogram using the proposed HNLR method exhibits the least amount of artifacts in the background area when compared to FBP and TVR, no artifacts in the imaged brain area, and noticeably better preservation of fine detail in the imaged brain area.

Figure 4
figure 4

Reconstructed tomograms from the PET2 data set using three different reconstruction methods based on 30% sampling. The PET2 data set was taken from a 53 year-old man who sought medical attention due to a grand mal seizure, with a brain biopsy revealing grade IV astrocytoma. As with the reconstructed tomograms from the PET1 data set, it can be observed that the FBP method results in the most significant artifacts in both the background and imaged brain areas when compared to TVR and HNLR, the TVR method results in reduced artifacts but noticeable loss in fine detail in the imaged brain area, and the HNLR method results in the least amount of artifacts and better preservation of fine detail in the imaged brain area. It is very clear that streaking artifacts (caused by sampling at a sub-Nyquist sampling fraction) is prominently overlaid on the entire brain region in the reconstruction produced using the FBP method (see red arrows), while they are significantly reduced in the reconstructions produced by TVR and HNLR. It is also important to highlight that, within the same brain region in the HNLR reconstruction, there is an overall improvement in sharpness when compared to the TVR reconstruction (see green arrows). Two line profiles (indicated by orange lines) from the PET data sets reconstructed using the methods are shown at the bottom, and further illustrates the reduction in artifacts and improved contrast achieved by TVR and HNLR when compared to FBP.

Figure 5
figure 5

Reconstructed tomograms from the PET3 data set using three different reconstruction methods based on 30% sampling. The PET3 data set was acquired from a 70 year-old man with mild Alzheimer’s disease about 9 months prior to imaging. As with the reconstructed tomograms from the PET1 and PET2 data sets, it can be observed that the FBP method results in the most significant artifacts in both the background and imaged brain areas when compared to TVR and HNLR. The TVR method results in reduced artifacts but noticeable loss in functional variation differences in the imaged brain area, and the HNLR method results in the least amount of artifacts and better preservation of functional variation differences in the imaged brain area.

Figure 6 shows the reconstruction results of all three tested methods at different sampling fractions (from 20% to 50%). Each row shows the reconstruction results of same method, while each column shows the reconstruction results acquired at the same sampling fraction, but with different methods. It can be observed that FBP produces results that contain a large amount of artifacts when sampled at low sampling fractions. The amount of artifacts decrease dramatically when using TVR compared to FBP. However, it can be seen that the reconstruction results using TVR have noticeable loss in detail. The third row of Figure 6 shows reconstruction results of the proposed method at three different sampling fractions, which has noticeably fewer artifact residuals compared to FBP and TVR even at a 20% sampling fraction. It is important to note that while the proposed HNLR reconstruction method can help reduce artifacts caused by sampling at sub-Nyquist sampling fractions, it is limited in its performance when dealing with very low sampling fractions, as evident in Figure 6(g) where noticeable artifacts still persists.

Figure 6
figure 6

Another set of reconstructed tomograms from the PET3 data set using three different reconstruction methods (FBP: a-c, TVR: d-f, HNLR: g-i) at 20%, 35% and 50% sampling fractions. It can be observed that the proposed HNLR method can achieve strong reconstruction results even at relatively low 20% sparse sampling fraction, while the reconstruction results of FBP has very strong artifacts and the reconstruction results of TVR show noticeable loss in detail.

Conclusion

The potential use of homotopic non-local regularization (HNLR) reconstruction framework to maintaining high reconstruction quality while reducing sampling fraction in sparse PET reconstruction was studied. The proposed reconstruction framework is designed to address issues associated with artifacts and detail loss with sparse PET reconstruction. A comparative analysis using real PET data sets was performed to compare traditional FBP as well as TVR reconstruction with the proposed HNLR reconstruction method. Based on experimental results, we illustrate that promising reconstruction performance can be achieved using HNLR even at low sampling fractions. Future work includes a more comprehensive evaluative study using larger PET data sets to further validate the performance of the proposed HNLR framework, the investigation of different sparsifying transforms to study their potential for improving reconstruction quality, and the investigation of incorporating alternative representations such as multiscale texture representations and local phase [32,33] for enforcing regularization. Furthermore, we wish to investigate more advanced data terms such as the Kullback-Leibler divergence [28,29], which takes better advantage of the Poisson noise statistics of PET, as well as investigate the application of the HNLR framework for other medical imaging modalities such as diffusion weighted magnetic resonance imaging [34,35] and correlated diffusion imaging [36]. We also aim to investigate methods for automatically optimizing the parameters of the HNLR framework to obtain optimal reconstruction accuracy. Finally, we wish to investigate the effects of employing the HNLR framework for sparse reconstruction on medical image analysis tasks such as medical image annotation and contouring [37], as well as medical image registration [32,33,38-41].

References

  1. Valk P, Bailey D, Townsend D, Maisey M. Positron emission tomography. London: Springer; 2002.

    Google Scholar 

  2. Brinks R, Busch M. Local compensation for respiratory motion in list-mode PET. Springer Proceedings in Physics: Adv Med Imaging. 2007; 114:31–6.

    Article  Google Scholar 

  3. Kalifa J, Laine A, Esser P. Tomographic reconstruction with non-linear diagonal estimators. Proc SPIE. 2000; 4119:1–11.

    Article  Google Scholar 

  4. Wang W, Gindi G. Noise analysis of MAP-EM algorithms for emission tomography. Phys Med Biol. 1997; 42:2215–32.

    Article  CAS  PubMed  Google Scholar 

  5. Wang C, Snyder W, Bilbro G, Santago P. Performance evaluation of filtered backprojection reconstruction and iterative reconstruction methods for PET images. Comput Biol Med. 1998; 28:13–24.

    Article  PubMed  Google Scholar 

  6. Jonsson E, Huang S, Chan T. Total variation regularization in positron emission tomography. Technical report, Los, Angeles, California: Dept Mathematics, Univ California; 1998, pp. 1–25.

    Google Scholar 

  7. Valiollahzadeh S, Chang T, John W Clark J, Mawlawi OR. Image recovery in pet scanners with partial detector rings using compressive sensing. IEEE Nucl Sci Symp Med Imaging Confererence. 2012:3036–9.

  8. Zhou J, Coatrieux J, Bousse A, Shu H, Luo L. A bayesian map-em algorithm for PET image reconstruction using wavelet transform. IEEE Trans on Nucl Sci. 2007; 54:1660–9.

    Article  Google Scholar 

  9. Pustelnik N, Chaux C, Pesquet JC, Sureau FC, Dusch E. Adapted convex optimization algorithm for wavelet-based dynamic PET reconstruction. In: the proceeding of the 10th, International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine: 2009. p. 5–10.

  10. Verhaeghe J, Ville DVD, Khalidov I, Asseler YD, Lemahieu I, Unser M. Dynamic PET reconstruction using wavelet regularization with adapted basis functions. IEEE Trans Med Image. 2008; 27:943–59.

    Article  Google Scholar 

  11. Candës E, Romberg J, Tao T. Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. IEEE Trans Inf Theory. 2006; 52:489–509.

    Article  Google Scholar 

  12. Trzasko J, Manduca A. Highly undersampled magnetic resonance image reconstruction via homotopic l0-minimization. IEEE Trans Med Imaging. 2009; 28:106–21.

    Article  PubMed  Google Scholar 

  13. Fieguth P. Statistical image processing and multidimensional modeling. New York: Springer; 2010.

    Google Scholar 

  14. Hanif A, Mansoor A, Ejaz T. Iterative tomographic image reconstruction by compressive sampling. In: Image Process (ICIP): 2010. p. 4313–6.

  15. Donoho D. Compressive sensing. IEEE Trans Inf Theory. 2006; 52:1289–306.

    Article  Google Scholar 

  16. Lustig M, Donoho D, Pauly J. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magn Reson Med. 2007; 58:1182–95.

    Article  PubMed  Google Scholar 

  17. Natarajan BK. Sparse approximate solutions to linear systems. SIAM J Comput. 1995; 24:227–34.

    Article  Google Scholar 

  18. Wong A, Mishra A, Clausi D, Fieguth P. Sparse reconstruction of breast MRI using homotopic l0 minimization in a regional sparsified domain. Biomed Eng IEEE Trans. 2010; 60:743–52.

    Article  Google Scholar 

  19. Gilboa G, Osher S. Nonlocal operators with applications to image processing. Multiscale Model Simul. 2008; 7:1005–28.

    Article  Google Scholar 

  20. Liang D, Wang H, Ying L. Sense reconstruction with nonlocal TV regularization. Proc IEEE Eng Med Biol Soc. 2009; 2009:1032–5.

    Google Scholar 

  21. Trzasko J, Manduca A, Borisch E. Sparse MRI reconstruction via multiscale L0-Continuation. In: Proc IEEE Workshop Stat Signal Process: 2007. p. 176–80.

  22. Wong A, Mishra A, Fieguth P, Clausi D. A perceptually adaptive approach to image denoising using anisotropic non-local means. In: Proc Int Conference Image Process: 2008. p. 1–4.

  23. Wong A, Orchard J. An adaptive non-local means approach to exemplar-based inpainting. In: Proc Int, Conference Image Process: 2008. p. 2600–3.

  24. Orchard J, Ebrahimi M, Wong A. Efficient nonlocal-means denoising using the svd. In: Proc Int Conference Image Process: 2008. p. 1732–5.

  25. Liu C, Wong A, Fieguth P, Bizheva K, Bie H. Homotopic, non-local sparse reconstruction of optical coherence tomography (OCT) imagery. Opt Express. 1020; 20:0–11.

  26. Rosenbloom P. The method of steepest descent. Proc Symp Appl Math. 1956; 6:127–76.

    Article  Google Scholar 

  27. Wang X. Method of steepest descent and its applications. IEEE Microwave Wireless Components Lett. 2008; 12:24–26.

    Google Scholar 

  28. Pustelnik N, Chaux C, Pesquet JC, Comtat C. Parallel algorithm and hybrid regularization for dynamic pet reconstruction. In: IEEE Nucl Sci Symp Med Imaging Confererence: 2010. p. 2423–7.

  29. Anthoine S, Aujol J, Mélot C, Boursier Y. Some proximal methods for cbct and pet tomography. inverse problems in imaging. Inverse Probl Imaging. 2012; 6:565–98.

    Article  Google Scholar 

  30. Anthoine S, Aujol J, Mélot C, Boursier Y. Asim pet simulation software. 2012. http://depts.washington.edu/asimuw/.

  31. Johnson K, Becker J. The whole brain atlas. 2008. www.med.harvard.edu/AANLIB/home.html.

  32. Wong A, Orchard J. Robust multimodal registration using local phase-coherence representations. J Signal Process Syst. 2009; 54:89–100.

    Article  Google Scholar 

  33. Wong A, Clausi D, Fieguth P. Cpol: Complex phase order likelihood as a similarity measure for mr-ct registration. Med Image Anal. 2010; 14:50–7.

    Article  PubMed  Google Scholar 

  34. Koh D, Padhani A. Diffusion-weighted mri: a new functional clinical technique for tumour imaging. Br J Radiology. 2006; 79:633–5.

    Article  Google Scholar 

  35. Shafiee M, Lui D, Haider S, Cameron A, Wong A, et al. Apparent ultra-high b-value diffusion-weighted image reconstruction via hidden conditional random fields. IEEE Trans Med Imaging. 2015; 1:1–15.

    Article  Google Scholar 

  36. Wong A, Glaister J, Cameron A, Haider M. Correlated diffusion imaging. BMC Med Imaging. 2013; 13:1–7.

    Article  Google Scholar 

  37. Mishra A, Wong A, Zhang W, Fieguth P, Clausi D. Improved interactive medical image segmentation using enhanced intelligent scissors (eis). In: Annu Int Conference IEEE Eng Med Biol Soc: 2008. p. 1–4.

  38. Wong A, Bishop W, Orchard J. Efficient multi-modal least-squares alignment of medical images using quasi-orientation maps. In: Proc Int Conference Image Process, Comput Vision, Pattern Recognit: 2006. p. 74–80.

  39. Wong A, Bishop W. Efficient least squares fusion of mri and ct images using a phase congruency model. Pattern Recognit Lett. 2008; 29:173–80.

    Article  Google Scholar 

  40. Roche A, Malandain G, Pennec X, Ayache N. The correlation ratio as a new similarity measure for multimodal image registration. Lecture Notes Comput Sci: Med Image Comput Comput-Assisted Intervention. 1998; 1496:1115–24.

    Article  Google Scholar 

  41. Wong A, Fieguth P. Fast phase-based registration of multimodal image data. Signal Process. 2009; 89:724–37.

    Article  Google Scholar 

Download references

Acknowledgements

This work was supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada, the Canada Research Chairs program, the Ontario Ministry of Economic Development and Innovation, and the Chinese Council Scholarship. The authors would also like to thank Dr. Keith A. Johnson for the test data.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alexander Wong.

Additional information

Competing interests

The authors have declared that no competing interests exist.

Authors’ contributions

Conceived and designed the HNLR methodology: AW, CL. Performed the experiments: AW, CL, XW. Analyzed the data: AW, CL, XW, PF, HB. Wrote the paper: AW, CL, XW, PF, HB. All authors read and approved the final manuscript.

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wong, A., Liu, C., Wang, X.Y. et al. Homotopic non-local regularized reconstruction from sparse positron emission tomography measurements. BMC Med Imaging 15, 10 (2015). https://doi.org/10.1186/s12880-015-0052-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12880-015-0052-5

Keywords