Next Article in Journal
A Structure Fidelity Approach for Big Data Collection in Wireless Sensor Networks
Previous Article in Journal
Evaluation of the Leap Motion Controller as a New Contact-Free Pointing Device
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified OMP Algorithm for Exponentially Decaying Signals

by
Krzysztof Kazimierczuk
1,* and
Paweł Kasprzak
2
1
Centre of New Technologies, University of Warsaw, Banacha 2C 02-097, Poland
2
Department of Mathematical Methods in Physics, Faculty of Physics, University of Warsaw, Pasteura 5, Warsaw 02-093, Poland
*
Author to whom correspondence should be addressed.
Sensors 2015, 15(1), 234-247; https://doi.org/10.3390/s150100234
Submission received: 12 September 2014 / Accepted: 15 December 2014 / Published: 24 December 2014
(This article belongs to the Section Sensor Networks)

Abstract

: A group of signal reconstruction methods, referred to as compressed sensing (CS), has recently found a variety of applications in numerous branches of science and technology. However, the condition of the applicability of standard CS algorithms (e.g., orthogonal matching pursuit, OMP), i.e., the existence of the strictly sparse representation of a signal, is rarely met. Thus, dedicated algorithms for solving particular problems have to be developed. In this paper, we introduce a modification of OMP motivated by nuclear magnetic resonance (NMR) application of CS. The algorithm is based on the fact that the NMR spectrum consists of Lorentzian peaks and matches a single Lorentzian peak in each of its iterations. Thus, we propose the name Lorentzian peak matching pursuit (LPMP). We also consider certain modification of the algorithm by introducing the allowed positions of the Lorentzian peaks' centers. Our results show that the LPMP algorithm outperforms other CS algorithms when applied to exponentially decaying signals.

1. Introduction

Compressed sensing (CS) theory provides the methods to solve an underdetermined linear equation y = Ax within a set of sparse vectors x ∈ ℂn. The unconstrained linear equation for an m × n matrix A of a maximum rank may have infinitely many solutions for m < n. However, the number of sparse solutions happens to be finite. Vectors in ℂn having at most l ∈ ℕ non-zero coordinates are referred to as the l-sparse vectors.

All l-sparse vectors fulfilling the equation y = Ax may be found by considering ( n l ) linear sub-problems of solving y = Ax with the constraint x ∈ ℂI, where ℂI ⊂ ℂn denotes the subspace of vectors supported by I ⊂ [1,n]. This provides all l-sparse solutions, but the overall problem was proven to be NP-hard. Remarkably, the sparsest vector solving y = Ax may be found as the solution of the l1-minimization problem:

x ^ = argmin y = A x x 1
where ‖ x1 denotes the l-norm of x. The l-norm convexity and the linear constraint y = Ax allow the application of linear programming algorithms, as long as A has a restricted isometry property; see [1].

A number of CS algorithms has been proposed, including: iterative soft thresholding (IST) [2], orthogonal matching pursuit (OMP) [3], compressive sampling matching (COSAMP) [4], stage-wise orthogonal matching pursuit (StOMP) [5] and iterative re-weighted least squares (IRLS) [6]. Inspired by the greedy CS algorithm, OMP, and motivated by a case of nuclear magnetic resonance (NMR) spectroscopy, we propose an approach that is designed specifically to recover a signal being a combination of exponentially decaying components. The spectrum of such a signal consists of Lorentzian peaks. The algorithm matches a single peak in each of its iterations, and thus, the name Lorentzian peak matching pursuit (LPMP) is proposed. The peaks' centers and heights are found as in the OMP algorithm, whereas the widths of the peaks are determined in a novel way.

2. NMR Spectroscopy and CS

The sources of the free induction decay (FID) signal measured in NMR spectroscopy are magnetic moments of atomic nuclei with non-zero spin, e.g., 1H,13C,15N,19F, that undergo coherent precession when excited by radio-frequency irradiation in the presence of a high static magnetic field [7]. The time dependence of an ideal FID signal f(t) is given by the following linear combination:

f ( t ) = k = 1 N a k f k ( t )
of the exponentially decaying oscillations:
f k ( t ) = { e 2 π γ k t e 2 π ι ω k t t 0 0 t < 0
where N is a number of nuclei groups differing in resonance frequency; ak are amplitudes (the (x0221D) number of nuclei in the k-th group); γk is a relaxation rate for the k-th group; and ωk is an oscillation frequency of the k-th group.

The precession frequencies (and thus, the frequencies of the emitted FID signal) are in the order of MHz, but their band is very narrow, typically a few kHz. In such a narrow range, tens or even hundreds of peaks are present. Often, the spectral resolution (peak width compared to the difference in resonance frequencies) is too low to resolve peaks and, thus, makes an analysis difficult. Fortunately, the peak dispersion can be improved by the acquisition of the multidimensional NMR signal [8] where additional (“indirect”) time dimensions are introduced. Peaks in the multidimensional NMR spectrum correspond to groups of nuclei coupled by some physical interaction. Using the above notation, the formula of the two-dimensional FID is:

f ( t 1 , t 2 ) = k 1 , k 2 = 1 N a k 1 , k 2 f k 1 ( t 1 ) f k 2 ( t 2 )
where t1 is indirectly sampled time and t2 directly sampled time.

The reason for the application of CS in multidimensional NMR spectroscopy is a time-consuming sampling of the indirect dimensions of the FID signal (t1 in Equation (2)). There are three reasons for the lengthy acquisition. Firstly, each measurement point in the indirect dimension takes a few seconds of “real” time. Secondly, the sampling density has to fulfil the Nyquist theorem [9]. Finally, the sampling has to reach far into the indirect time domain to provide sufficient spectral resolution [10]. The latter, together with the Nyquist condition for the sampling density, means that, often, many thousands of sampling points have to be acquired. This makes an experiment unacceptably long (days in the case of three or four dimensions) or completely impossible in the case of higher dimensionality (5D and more). Thus, there is a need for alternative sampling approaches that allow one to save expensive spectrometer time preserving the spectral information. Most of them use the non-uniform sampling, where the major part of the points is not measured, but reconstructed after the experiment based on certain assumptions [11,12]. The CS, where the assumption is spectral sparsity, has been successfully applied in NMR exploiting the algorithms based on convex l1-norm or non-convex lp-norm minimization [13,14]. The history of the application of the “greedy” algorithm, OMP, in NMR is long; in fact, the algorithm was transferred from radio astronomy, where it existed under the name, CLEAN [15]. Direct application of CLEAN in NMR is possible [16], but far from being perfect, because of the aforementioned approximate sparsity of the signal. Thus, various modifications of CLEAN have been proposed [17-19]. Importantly, none of these works refer to OMP formalism and do not explain the modifications by the sparsity criteria. We propose an approach that directly exploits the fact that the NMR spectrum consists of a low number of Lorentzian peaks. Thus, LPMP is based on the sparsity concept, which is more relevant in the case of NMR than classical CS methods.

3. OMP Algorithm

Let us fix the notation. For m,n ∈ ℤ and m < n, we write [m, n] = {m, m + 1,…, n}. Let m, n ∈ ℕ. For matrix A ∈ Mm×n(ℂ), its range and kernel are denoted ran A and ker A, respectively. The Hermitian conjugate of A is denoted by A*. The number of non-zero coordinates of xCn is denoted by ‖x0, and supp x (the support of x) denotes the set of x non-zero coordinates. Note that ‖x0 = # supp x (the cardinality of a set I is denoted by # I). We define χn,l ⊂ ℂn:

χ n , l = { x C n : x 0 l }

OMP algorithm arose in the context of the sparse approximation problem; see [3]. In order to formulate the problem, let us fix a matrix A ⊂ Mm×n(ℂ) satisfying ran A = ℂm. The matrix A is in this context referred to as the dictionary matrix.

Definition 1

Adopting the above notation, we say that xχn,l is an optimal l-sparse approximate representation of a vector y ∈ ℂm if:

y A x 2 = min { y A x 2 : x χ n , l }

We say that y is l-representable, if there exists a vector xχn,l, such that y = Ax.

For a subset I ⊂ [1, n] of cardinality l, we define a vector subspace ℂI ⊂ ℂn:

C I = { x C n : supp ( x ) I }

The restriction of A to ℂI is denoted by AI : ℂI → ℂm.

Definition 2

We say that A ⊂ Mm×n(ℂ) is l-injective if ker AI = {0} for all I ⊂ [1, n], # I = l.

For an l-injective matrix A and I ⊂ [1, n], # I = l, we define:

x 1 = argmin x C I y A I x 2

Since xI is a vector satisfying AIxI = yI, where yI is an orthogonal projection of y onto ran AI, we may see the uniqueness of xI for l-injective A. Let us note that xχn,l is an optimal l-sparse approximate representation of y ∈ ℂm if:

y A x 2 min { I : # I = l } y A x I 2

Although l-injectivity of the matrix A does not exclude a vector x′ ∈ χn,l, x′ ≠ x, such that ‖yAx2 = ‖yAx2, we shall usually choose an optimal l-sparse approximate and denote it by xl.

The OMP algorithm provides a certain approximation of xl ∈ ℂn, which we denote by x l OMP Remarkably, x l OMP is found on the basis of x l 1 OMP Denoting the x l OMP support by I l OMP, the index il ∈ [1, n] is given by:

i l = argmax i [ 1 , n ] | A * ( y A x l 1 ) | i

The OMP vector x l OMP is defined as x l OMP = x I l OMP; see Equation (3) for the definition of x I l OMP. For the coherence, restricted isometry property (RIP) conditions that imply the equality x l OMP = x l, we refer to [3] and [20], respectively. For the recovery limits of the OMP, see [21]. The OMP algorithm (Algorithm 1) is schematically described below.


Algorithm 1 Orthogonal matching pursuit.

  • Input:

    • measurement matrix A ∈ Mm×n(ℂ)

    • measurement vector y ∈ ℂm

    • accuracy parameter ε > 0

  • Output:

    • x ∈ ℂn

  • Initialization:

    • I = ∅, x = 0

  • The main loop:

    whileyAx2ε do

    • I = I { argmax j [ 1 , n ] | A * ( y A x ) j | }

    • x = argmin supp ( z ) I y A z 2

    end while


4. Free Induction Decay Signals in NMR

The free induction decay signal in NMR consists of the combination of exponentially decaying oscillations (see Equations (1) and (2)). The signal can be multidimensional, but for the sake of simplicity, we will limit ourselves to the one-dimensional indirect part of the two-dimensional signal, corresponding to f(t) from Equation (1). Given the damping factor γ ≥ 0 and the oscillation frequency ω ∈ (x0211D), a single decaying oscillation is given by:

f ( t ) = { e 2 π γ t e 2 π l ω t t 0 0 t < 0

Its Fourier transform:

f ^ ( v ) = R f ( t ) e 2 π l v t d t
has the form of complex Lorentzian peak:
f ^ ( v ) = 1 2 π 1 γ ι ( ω v )

The FID signal is sampled at the rate given by the Nyquist theorem for the range cut off ν ∈ [–Ω, Ω]. In order to perform the NMR signal processing, the frequency resolution ΔΩ is introduced, where Δ Ω = Ω 2 N + 1 and N ∈ ℕ. For any l ∈ [–N,N], we get the basic complex oscillation e2πιωlt, where ω l = l Ω 2 N + 1. The NMR signals can be approximated by the finite combination of the basic complex oscillations:

f ( t ) l = N N a l e 2 π l ω l t

The approximation being 2 N + 1 Ω -periodic works only for 0 t 2 N + 1 Ω. Assuming that the series l = N N a l e 2 π l ω l tevaluated at t k = k Ω coincides with sampling points of f, i.e., f(tk) for k ∈ [0,2N], we get:

a l = 1 2 N + 1 k = 0 2 N f ( t k ) e 2 π l ω l t k

Introducing the time unit 1 Ω and the frequency unit Ω 2 N + 1, we may enumerate the time and frequency by the integers [0, 2N] and [–N, N], respectively. Then, the relation between the peak intensities a and the time samples f is given by a = Ff, where F is the Fourier transform matrix:

F k l = e 2 π l k l / 2 N + 1

l ∈ [–N, N] and k ∈ [0, 2N]. In order to model the Lorentzian peaks within our framework, let us introduce the width peak unit Γ and a damping matrix C ∈ M(2N+1)×(2N+1)(ℂ):

C k k = δ k k e 2 π k Γ / Ω
for k, k′ ∈ [0, 2N]. Then, the Lorentzian peak of width Γ centered at l is represented by Ll :
L l = F 1 C F l
where Fl is the l-th column of F. Similarly, the Lorentzian peak of the width jΓ centered at l is represented by:
L l j = F 1 C j F l

Using the above notation, the FID signal is given by:

f = j l b j l F L l j
where bjl ∈ ℂ.

5. Lorentzian Peak Matching Pursuit

The spectrum corresponding to FID signal (Equation (5)) is given by:

f ^ = j l b j l L l j

We shall assume that j ∈ [0, J], which gives the maximal width JΓ and the width resolution Γ. Let K = {k1, k2,…, km} ⊂ [0, 2N] be the sampling scheme, y ∈ ℂm the measured signal yi = f(tki) and A = Fk the matrix consisting of rows of the Fourier matrix F enumerated by the elements of K.

The Lorentzian peak matching pursuit algorithm matches a Lorentzian peak L l j in each of its iterations. In order to match a peak L l j LPMP establishes first the peak's center l. The center l1 ∈ [–N, N] of the first peak is determined by:

l 1 argmax l [ N , N ] | ( A * y ) l |

In order to determine the width of the first peak j1 ∈ [0, J], we find bl1,j ∈ ℂ, such that:

y b l 1 , j A L l 1 j 2 = min { y z A L l 1 j 2 : z C }

Then, j1 ∈ [0, J] is a number satisfying:

y b l 1 , j 1 A L l 1 j 1 2 = min { y b l 1 , j A L l 1 j 2 : j [ 0 , J ] }
and we put f ^ 1 = b l 1 , j 1 L l 1 j 1.

In the second step of the algorithm, we determine the second peak's center l2 by:

l 2 = argmax l [ N , N ] { | A * ( y A f ^ 1 ) l | }

In order to determine the second peak's width j2 ∈ [0, J], we find bl2,j, bl1,j1 ∈ ℂ minimizing:

y z 2 A L l 2 j z 1 A L l 1 j 1 2
where z1,z2 ∈ ℂ Defining j2 E [0, J] as a number satisfying:
y b l 2 , j 2 A L l 2 j 2 b l 1 , j 1 A L l 1 j 1 2 = min { y b l 2 A L l 2 j b l 1 , j 1 A L l 1 j 1 2 : j [ 0 , J ] }
we put f ^ 2 = b l 2 , j 2 L l 2 j 2 + b l 1 , j 1 L l 1 j 1.

We begin the third step of the algorithm by determining the third peak's center l3:

l 3 = argmax l [ N , N ] | A * ( y A f ^ 2 ) l |

Then, j3 and l3 are found as in the previous steps.

Note that in the k-th step of the algorithm, the peak's width jkΓ is found and not changed afterwards. The blk,jk coordinate may, in turn, vary in the n-th step of the algorithm for n > k. The algorithm is performed until the accuracy threshold ‖yAx‖ < ε is achieved; for a different stopping criterion, see Section 5.1.

5.1. Stopping Criterion

If a number of Lorentzian peaks in the spectrum cannot be assumed a priori, then the stopping criterion for the LPMP becomes a crucial issue.

Notably, the dimension of a measurement vector y provides a mathematical bound for the number of iterations in LPMP. The implementation of the LPMP algorithm uses an inversion of a matrix, which ceases to be invertible when the number of peaks exceeds the number of measurements. In what follows, we shall introduce the stopping criteria that further restrict the number of iterations. One commonly used criterion is the threshold ε of accuracy used in Algorithm 2.

Let us denote the result of the k-th LPMP iteration by k. The k + 1 iteration is performed if:

y A f ^ k 2 α y A f ^ k 1 2
where α ∈ [0,1] is an user-defined number. The criterion reflects the idea that the LPMP is continued as long as its actual iteration increases the explanation of the measured data by the factor α with respect to the previous iteration. In our simulations described in the next section, we put α = 0.98.

Alternatively, the noise amplitude σ can be used for a stopping criterion. In this case, LPMP stops on the k-th iteration if |f̂k(lk)| ≤ σ.


Algorithm 2 Lorentzian peak matching pursuit.

  • Input:

    • measurement matrix A ∈ MK×2N+1(ℂ)

    • measurement vector y ∈ ℂK

    • accuracy ε > 0

    • Lorentzian peaks L l i C 2 N + 1 where j ∈ [0, J], l ∈ [−N, N]

  • Output:

    • f̂LPMP∈ ℂ2N+1

  • Initialization:

    • f̂LPMP= 0, j = 0, b = 0, i = 0

  • The main loop:

    whileyAf̂LPMP2ε do

    • i = i + 1

    • l i arg max l [ N , N ] | A * ( y A f ^ LPMP ) l |

    • ( b l i , j , b l i 1 , j j 1 , , b l i , j 1 ) = argmin z C i y A ( z i L l i j + z i 1 L l i 1 j i 1 + z 1 L l 1 j 1 ) 2

    • j i = argmin j [ 0 , J ] y A ( b l i , j L l i j + b l i 1 , j i 1 L l i 1 j i 1 + b l 1 , j 1 L l 1 j 1 ) 2

    • f ^ LPMP = b l i , j i L l i j i + b l i 1 , j i 1 L l i 1 j i 1 + b l 1 , j 1 L l 1 j 1

    end while


5.2. Mask

NMR experiments are often performed in series, with some or all peaks preserving their positions in at least some of the dimensions. This fact has been exploited to restrict the allowed peak frequencies and to improve the results of NUS (non-uniform sampling) reconstruction, e.g., in hyperdimensional spectroscopy [22], the SIFT method [23] or high-dimensional spectra [24]. It is easy to use this information in the LPMP algorithm by introducing the subset Mask ⊂ [−N, N] of admissible peak centers in the frequency domain. In what follows, we give a detailed description of the masked LPMP algorithm (Algorithm 3), taking into an account the stopping criterion described in Section 5.1.


Algorithm 3 Lorentzian peak matching pursuit with a mask.

  • Input:

    • measurement matrix A ∈ MK×2N+1(ℂ)

    • measurement vector y ∈ ℂK

    • stopping parameter α > 0

    • Lorentzian peaks L l i C 2 N + 1 where j ∈ [0, J], l ∈ [−N, N]

    • Mask ⊂ [−N, N]

  • Output:

    • f̂ LPMP ∈ ℂ2N+1

  • Initialization:

    • f̂ LPMP = 0, j = 0,b = 0,i = 0,δ = ‖y2

  • The main loop:

    whileyAf̂LPMP2 ≤ αδ or i = 0 do

    • δ = ‖yAf̂2

    • i = i + 1

    • l i arg max l Mask | A * ( y A f ^ LPMP ) l |

    • ( b l i , j , b l i 1 , j i 1 , , b l 1 , j 1 ) = argmin z C i y A ( z i L l i j + z i 1 L l i 1 j i 1 + + z 1 L l 1 j 1 ) 2

    • j i = argmin j [ 0 , J ] y A ( b l i , j L l i j + b l i 1 , j i 1 L l i 1 j i 1 + b l 1 , j 1 L l 1 j 1 ) 2

    • f ^ LPMP = b l i , j i L l i j i + b l i 1 , j i 1 L l i 1 j i 1 + + b l 1 , j 1 L l 1 j 1

    end while


6. Results and Discussion

In this section, we present the results of the LPMP performance tests. Extensive simulations on different levels of sampling sparsity provided a quantitative characterization of the method. As an example of an application, we have chosen the challenging NOESY (nuclear Overhauser effect spectroscopy) experiment [25] with a high dynamic range of peak intensities.

The parameters of the synthetic spectrum used in the simulations are given in Table 1. The spectral parameters were chosen to simulate various experimental scenarios, i.e., (I) peak widths differ from peak to peak, which corresponds to the differences in relaxation rates in the case of the NMR experiment; (II) the range of intensities is high; (III) two of the peaks overlap; and (IV) the peaks' centers and widths are off-grid.

An example of the LPMP recovery of an NMR spectrum based on 120 sampling points drawn at random from a 256-point full grid is presented in Figure 1. The signal was contaminated with 5% of noise. We performed systemic simulations with a varying sampling level. For each level, 50 different sampling schedules were generated. Denoting the frequency spectrum vector obtained from the full sampling by and the LPMP recovery by LPMP, we compute the recovery error ‖LPMP‖, and then, we average it over 50 sampling distributions. The LPMP performance is depicted in Figure 2 (left), where it is also compared with other popular CS algorithms: IST, OMP, COSAMP and StOMP. For all sampling levels, LPMP is revealed to be superior to other CS algorithms.

In order to test the masked version of LPMP algorithm, we introduced the mask of ±5 frequency points around the peaks' centers. The simulation results are presented in Figure 2 (right). As expected, masking has a more pronounced effect on the low sampling levels.

The qualitative aspects of LPMP reconstruction can be evaluated on the data from the 2D NOESY experiment. Each reconstructed 1D slice of the 2D NOESY spectrum contains one dominating peak accompanied by a group of smaller peaks. Such a case of large intensity range is very challenging for CS methods. In Figure 3, we present the results of the LPMP spectrum recovery obtained on the basis of 200, 225 and 250 random measurements of the corresponding FID signal of length 512. The 2D NOESY spectrum was recorded on a 600-MHz Varian UNITY Inova spectrometer at 25 °C using a sample of human ubiquitin protein (1 mM concentration).

In order to explain the effectiveness of LPMP algorithm in the NMR context, let us emphasize the key difference between the OMP and LPMP methods. OMP tries to obtain the frequency-domain representation having the minimal number of non-zeros, which cannot be exact, even for a single Lorentzian peak. LPMP improves OMP by replacing every determined spectral non-zero with the Lorentzian peak of the appropriate width. Roughly speaking, this corresponds to the concept of sparsity in the sense of a number of Lorentzian peaks.

7. Conclusions

Our results show that LPMP is superior to the generic CS methods when applied to exponentially decaying signals. Comparison with OMP, StOMP, IST and COSAMP revealed that the assumption of the Lorentzian spectral line shape improves the quality of the reconstruction with LPMP. Moreover, its performance can be further improved by limiting (“masking”) the support of a signal. We believe that this method will find applications in NMR spectroscopy. Other fields of possible application may include different kinds of spectroscopy, wireless communications, sonar and radar. In fact, reconstruction of any non-uniformly sampled exponentially-damped signal can benefit from the application of LPMP.

Acknowledgments

The authors thank the Polish Ministry of Science and Higher Education for support with Iuventus Plus Grant No. IP2011 023171 and the National Centre of Science for support with Sonata Bis 2 Grant No. DEC-2012/07/E/ST4/01386. Krzysztof Kazimierczuk thanks the Foundation for Polish Science for the support with the Team program.

Author Contributions

Paweł Kasprzak developed the core of the algorithm and described it. Krzysztof Kazimierczuk adapted it to NMR requirements and wrote the NMR part of the publication.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Candès, E.J.; Romberg, J.K.; Tao, T. Stable signal recovery from incomplete and inaccurate measurements. Commun. Pure Appl. Math. 2006, 59, 1207–1223. [Google Scholar]
  2. Papoulis, A. A new algorithm in spectral analysis and band-limited extrapolation. IEEE T Circuits Syst. 1975, 22, 735–742. [Google Scholar]
  3. Tropp, J.A. Greed is good: Algorithmic results for sparse approximation. IEEE Trans. Inf. Theory 2004, 50, 2231–2242. [Google Scholar]
  4. Needell, D.; Tropp, J. CoSaMP: Iterative signal recovery from incomplete and inaccurate samples. Appl. Comput. Harmonic Anal. 2009, 26, 301–321. [Google Scholar]
  5. Donoho, D.L.; Tsaig, Y.; Drori, I.; Starck, J.L. Sparse solution of underdetermined linear equations by stage-wise orthogonal matching pursuit. IEEE Trans. Inf. Theory 2012, 58, 1094–1121. [Google Scholar]
  6. Candès, E.J.; Wakin, M.B.; Boyd, S.P. Enhancing Sparsity by Reweighted l1 Minimization. J. Fourier Anal. Appl. 2008, 14, 877–905. [Google Scholar]
  7. Keeler, J. Understanding NMR Spectroscopy, 2nd ed.; Wiley: Chichester, UK, 2010; p. 526. [Google Scholar]
  8. Jeener, J. AMPERE International Summer School; Basko Polje, Yugoslavia, 1971. [Google Scholar]
  9. Nyquist, H. Certain topics in telegraph transmission theory. Trans. Am. Inst. Electr. Eng. 1928, 47, 617–644. [Google Scholar]
  10. Szantay, C. NMR and the uncertainty principle: How to and how not to interpret homogeneous line broadening and pulse nonselectivity. IV. (Un?)certainty. Concept. Magn. Reson. A. 2008, 32A, 373–404. [Google Scholar]
  11. Kazimierczuk, K.; Stanek, J.; Zawadzka-Kazimierczuk, A.; Koźmiński, W. Random sampling in multidimensional NMR spectroscopy. Prog. Nucl. Magn. Reson. Spectrosc. 2010, 57, 420–434. [Google Scholar]
  12. Coggins, B.E.; Venters, R.A.; Zhou, P. Radial sampling for fast NMR: Concepts and practices over three decades. Prog. Nucl. Magn. Reson. Spectrosc. 2010, 57, 381–419. [Google Scholar]
  13. Kazimierczuk, K.; Orekhov, V.Y. Accelerated NMR spectroscopy by using compressed sensing. Angew. Chem. Int. Ed. Engl. 2011, 50, 5556–5559. [Google Scholar]
  14. Holland, D.J.; Bostock, M.J.; Gladden, L.F.; Nietlispach, D. Fast multidimensional NMR spectroscopy using compressed sensing. Angew. Chem. Int. Ed. Engl. 2011, 50, 6548–6551. [Google Scholar]
  15. Högbom, J.A. Aperture synthesis with a non-regular distribution of interferometer baselines. Astron. Astrophys. Suppl. 174. 15, 417–426.
  16. Coggins, B.E.; Zhou, P. High resolution 4-D spectroscopy with sparse concentric shell sampling and FFT-CLEAN. J. Biomol. NMR 2008, 42, 225–239. [Google Scholar]
  17. Stanek, J.; Kozmiński, W. Iterative algorithm of discrete Fourier transform for processing randomly sampled NMR data sets. J. Biomol. NMR 2010, 47, 65–77. [Google Scholar]
  18. Coggins, B.E.; Werner-Allen, J.W.; Yan, A.; Zhou, P. Rapid Protein Global Fold Determination Using Ultrasparse Sampling, High-Dynamic Range Artifact Suppression, and Time-Shared NOESY. J. Am. Chem. Soc. 2012, 134, 18619–18630. [Google Scholar]
  19. Kazimierczuk, K.; Zawadzka, A.; Kozmiński, W.; Zhukov, I. Lineshapes and Artifacts in Multidimensional Fourier Transform of Arbitrary Sampled NMR Data Sets. J. Magn. Reson. 2007, 188, 344–356. [Google Scholar]
  20. Zhang, T. Sparse recovery with orthogonal matching pursuit under RIP. IEEE Trans. Inf. Theory 2011, 57, 6215–6221. [Google Scholar]
  21. Wang, J.; Shim, B. On the Recovery Limit of Sparse Signals Using Orthogonal Matching Pursuit. IEEE Trans. Signal Process 2012, 60, 4973–4976. [Google Scholar]
  22. Jaravine, V.A.; Zhuravleva, A.V.; Permi, P.; Ibraghimov, I.; Orekhov, V.Y. Hyperdimensional NMR spectroscopy with nonlinear sampling. J. Am. Chem. Soc. 2008, 130, 3927–3936. [Google Scholar]
  23. Matsuki, Y.; Konuma, T.; Fujiwara, T.; Sugase, K. Boosting Protein Dynamics Studies Using Quantitative Nonuniform Sampling NMR Spectroscopy. J. Phys. Chem. B. 2011, 115, 13740–13745. [Google Scholar]
  24. Kazimierczuk, K.; Zawadzka-Kazimierczuk, A.; Koźmiński, W. Non-uniform frequency domain for optimal exploitation of non-uniform sampling. J. Magn. Reson. 2010, 205, 286–292. [Google Scholar]
  25. Kaiser, R. Use of the Nuclear Overhauser Effect in the Analysis of High-Resolution Nuclear Magnetic Resonance Spectra. J. Chem. Phys. 1963, 39, 2435–2442. [Google Scholar]
Figure 1. Example of the Lorentzian peak matching pursuit (LPMP) recovery. The size of the spectrum is 256. The number of randomly chosen samples is 120. The result of LPMP recovery (blue line) is compared with the original spectrum (red line).
Figure 1. Example of the Lorentzian peak matching pursuit (LPMP) recovery. The size of the spectrum is 256. The number of randomly chosen samples is 120. The result of LPMP recovery (blue line) is compared with the original spectrum (red line).
Sensors 15 00234f1 1024
Figure 2. Quality of the reconstructions of the spectrum from Figure 1. (Left) The comparison of the performance of various algorithms: LPMP (solid line), iterative soft thresholding (IST) (dashed line), stage-wise orthogonal matching pursuit (StOMP) (dotted line), compressive sampling matching (COSAMP) (stared line) and orthogonal matching pursuit (OMP) (squared line). (Right) The performance of the masked LPMP algorithm (solid line) compared with the performance of non-masked LPMP algorithm (dashed line). The absolute recovery errors (in intensity units) are plotted vs. the number of sampling points.
Figure 2. Quality of the reconstructions of the spectrum from Figure 1. (Left) The comparison of the performance of various algorithms: LPMP (solid line), iterative soft thresholding (IST) (dashed line), stage-wise orthogonal matching pursuit (StOMP) (dotted line), compressive sampling matching (COSAMP) (stared line) and orthogonal matching pursuit (OMP) (squared line). (Right) The performance of the masked LPMP algorithm (solid line) compared with the performance of non-masked LPMP algorithm (dashed line). The absolute recovery errors (in intensity units) are plotted vs. the number of sampling points.
Sensors 15 00234f2 1024
Figure 3. 1D indirect dimension slice from 2D NOESY spectrum: LPMP-reconstruction from 200, 225 and 250 samples out of a 512-point full grid compared with a fully sampled spectrum.
Figure 3. 1D indirect dimension slice from 2D NOESY spectrum: LPMP-reconstruction from 200, 225 and 250 samples out of a 512-point full grid compared with a fully sampled spectrum.
Sensors 15 00234f3 1024
Table 1. Peak coordinates for the simulation presented in Figure 1.
Table 1. Peak coordinates for the simulation presented in Figure 1.
Center5.622.767.870.4194.5239.25
width5.25.42.62.841.5

Share and Cite

MDPI and ACS Style

Kazimierczuk, K.; Kasprzak, P. Modified OMP Algorithm for Exponentially Decaying Signals. Sensors 2015, 15, 234-247. https://doi.org/10.3390/s150100234

AMA Style

Kazimierczuk K, Kasprzak P. Modified OMP Algorithm for Exponentially Decaying Signals. Sensors. 2015; 15(1):234-247. https://doi.org/10.3390/s150100234

Chicago/Turabian Style

Kazimierczuk, Krzysztof, and Paweł Kasprzak. 2015. "Modified OMP Algorithm for Exponentially Decaying Signals" Sensors 15, no. 1: 234-247. https://doi.org/10.3390/s150100234

APA Style

Kazimierczuk, K., & Kasprzak, P. (2015). Modified OMP Algorithm for Exponentially Decaying Signals. Sensors, 15(1), 234-247. https://doi.org/10.3390/s150100234

Article Metrics

Back to TopTop