![MathJax Logo](/templates/jsp/_style2/_tandf/pb2/images/math-jax.gif)
Abstract
Purpose: The purpose of this work was to develop an online dynamic cardiac MRI model to reconstruct image frames from partial acquisition of the Cartesian k-space data, which utilizes structural knowledge of consecutive image frames. Materials and methods: Using an elastic-net model, the proposed algorithm reconstructs dynamic images using both L1 and L2 norm operations. The L1 norm enforces the sparsity of the frame difference, while the L2 norm with motion-adaptive weights catches the internal structure of frame differences. Unlike other online methods such as the Kalman filter (KF) technique, the new model requires no assumption of Gaussian noise, and can faithfully reconstruct the dynamic images within a compressive sensing framework. Results: The proposed method was evaluated using simulated dynamic phantoms with 40 frames of images (128 × 128) and a cardiac MRI cine of 25 frames (256 × 256). Both results showed that the new model offered a better performance than the online KF method in depicting simulated phantom and cardiac dynamics. Conclusion: It is concluded that the proposed imaging model can be used to capture a large variety of objects in motion from highly under-sampled k-space data, and being particularly useful for improving temporal resolution of cardiac MRI.
Introduction
Dynamic MRI is an important imaging technique in clinical applications, such as interventional and functional imaging, that need to capture objects such as heart and intake agent in motion with high spatial and temporal resolutions. Ideally, this can be achieved by acquiring full k-space data of each frame and then reconstructing the image via inverse Fourier transform. In this straightforward method, the temporal resolution hinges on the speed of acquisition, which is limited by physical or physiological constraints. To improve temporal resolution or to speed-up scan procedure, one effective method is to acquire partial k-space data of each frame. In recent years, several strategies [Citation1] have been proposed to reconstruct images from partially acquired k-space data via effective exploitation of spatial and/or temporal correlations of dynamic images. For example, dynamic models have been previously built to exploit temporal correlation between frames. Using the block match algorithm from video compression, k–t FOCUSS (FOCal Underdetermined System Solver) with ME/MC (Motion Estimation/Motion Compensation) [Citation2] estimates the motion by relocating blocks of pixels from a reference frame into current frame and then encodes the sparse residual with k–t FOCUSS [Citation3]; MASTeR (Motion-Adaptive Spatio-Temporal Regularization) [Citation4] constructs a motion-adaptive linear dynamic system and recovers images based on the sparsity of forward and backward differences. However, they both work offline, that is, the acquisition and reconstruction of each frame cannot proceed alternately, because the reconstruction of current frame depends on the information of subsequent frames. As a result, those offline methods limit their usages in interventional applications, such as MRI-guided surgery.
To monitor the motion of objects during scan, several online reconstruction techniques have been proposed. Bilen [Citation5] uses optical flow to estimate the object’s motion, which relocates each pixel from previous frame into current frame, and then recovers images by compressed sensing (CS). However, because both quality reference frame and current frame are not available in the process of imaging as opposed to video compression, the motion estimation by simply relocating pixels from reference frame may incur loss of quality for reconstructed images. To avoid this problem, the Kalman filter (KF) technique has been used to estimate unknown variables recursively from measurements observed over time. In dynamic MRI, the CS scheme is incorporated into the framework of KF to reconstruct the sparse residual [Citation6]. Sümbület [Citation7,8] modelled dynamic MRI as a discrete linear system, and recovered images out of the model frame by frame using KF. However, to make KF computationally tractable, the covariance matrix in KF is simplified by ignoring off-diagonal elements, which can only be possible for non-Cartesian sampling trajectory. For Cartesian trajectory, Feng [Citation9] broke down an image into columns with much lower dimensions along the phase-encoding direction for each read-out position, and then reconstructed each column independently using KF. Hong [Citation10] improved the models proposed in [Citation9] by formulating KF as an optimization programme and then imposing additional sparsity constraints on it. However, in both models of [Citation9] and [Citation10], the system noise of KF model is actually the difference between two consecutive frames, which is usually not Gaussian and hence violates the condition of the KF.
In this work, with the consideration of the fact that the frame difference is caused by object motion and is usually sparse, an elastic-net [Citation11] model is proposed that describes the object dynamics in an online manner. The method reconstructs images using a hybrid optimization procedure that combines the L1 and L2 norms, in which L1 norm is used to promote sparsity and L2 norm with motion-adaptive weights is employed to catch the internal structure of frame differences. The proposed imaging model does not require Gaussian assumption of the system and measurement noise, as does the KF. Simulated phantom and cardiac MRI cine were used to demonstrate improved performance of the proposed method.
Theory
Overview of the discrete linear system model for dynamic MRI
This section briefly reviews the discrete linear system for dynamic MRI. The temporal correlation of consecutive frames of the imaging objects can be described by a discrete linear system. When the k-space is sampled in a Cartesian trajectory, dynamic MRI can be described by(1a)
(1a)
(1b)
(1b)
where is an image column vector assuming the row vector is along the phase-encoding direction,
is the corresponding column vector after applying an inverse Fourier transform to row vectors of acquired k-space data,
is the under-sampled Fourier transform,
is the frame difference,
is the measurement noise of MRI scanners. The subscript
denotes the time index of a frame. For an N by N MR image,
and
are vectors of size N and
is a matrix of K by N, where K is the number of sampled lines. The relationship of
and
is shown in Figure .
Figure 1. The reconstruction of a 3D cine MRI can be broken into the reconstruction of multiple 2D x–t images along the read-out direction. The symbol ⊙ denotes element-wise multiplication.
![Figure 1. The reconstruction of a 3D cine MRI can be broken into the reconstruction of multiple 2D x–t images along the read-out direction. The symbol ⊙ denotes element-wise multiplication.](/cms/asset/d03b1639-4b10-4174-b361-b1981a60e5b7/gipe_a_1141207_f0001_b.gif)
Given and
,
can be reconstructed by a model formulized in terms of
and
with a prior structure. The estimate for
is then obtained by adding reconstructed
to the previous frame
. With different assumptions of the structure of frame difference
, various modalities had been proposed to recover images from Equation (1). For example, reference [Citation12] modelled
as a rank-deficient matrix and images are recovered from the partially sampled k-space data via nuclear norm minimization; references [Citation13,14] observed that
is sparse and thus the dynamic frames can be recovered with the help of CS technique; reference [Citation9] assumed both
and
are ‘Gaussian’ signal, thus KF was applied to recover
for Cartesian trajectory. The KF is a recursive algorithm that estimates the state
of a process in a way that minimizes the variance of the estimation error when
and
represent independent white-noise sequences.[Citation16,19] It incorporates a new measurement
into the a priori estimate
to obtain an improved a posteriori estimate
.[Citation16] Compared with other filters such as Wiener filter [Citation19], the KF is both memory and computation efficient and is well suited to online/real-time execution.
Specifically, given , the KF technique produces a recursive solution to the above discrete linear system, that can be described using the following equations [Citation16]
(2a)
(2a)
(2b)
(2b)
(2c)
(2c)
where is defined to be the estimation error covariance matrix and
be the Kalman gain [Citation16] at the moment
. Substituting Equation (2a) into Equation (2b), we have
(3)
(3)
If the measurement noise covariance matrix is diagonal and the time-invariant, that is
, Equation (3) is actually equivalent to the following equation
(4)
(4)
Because is symmetric and semi-definite, and it can be decomposed as
, so Equation (4) reveals that KF minimizes
with the fidelity constraint
.
Elastic-net imaging model
To further reveal the structure of frame difference, we depict one of frame differences with the most non-zero elements for a cardiac cine, as shown in Figure , indicating that the frame difference reflects the motion of the heart (Figure (a)), and that it is compressible and sparse (Figure (b)).
In dynamic MRI, large similarity between adjacent frames usually exists; in other words, the differences of them are generally small, as demonstrated by Figure (a). So, the following imaging model can be initially set up(5)
(5)
Equation (5) is a mixed-norm model and is known as the elastic-net [Citation11] in the field of statistics. It can also be rewritten to(6)
(6)
In the original elastic-net model, is a constant scalar. That means all elements of
are penalized with the same weights towards zero. However, Figure (a) shows that the elements of
should be penalized adaptively to catch the regional variation of the motion of the imaging object. Let
denote the adaptive weights of
, then Equation (5) can be re-expressed as
(7)
(7)
where is the regularization constant for L2 norm term, and Equation (7) is equivalent to
(8)
(8)
Equation (8) is the proposed imaging model and can be solved using various sparse optimization solvers.
Relationship to existing work
In [Citation10], to enforce the sparsity of the KF solution , Equation (4) is solved with an L1 norm constraint as
(9)
(9)
Equation (9) was then solved using FOCUSS [Citation15] and the results show that this sparsity-enforced KF technique statistically improves the original KF.[Citation9]
Comparing Equation (9) with Equation (7), we show that the difference between the KF technique and the proposed method is that the KF weights using the inverse of
, while the proposed method uses an estimated vector of weights. Specifically, Equation (9) is different from Equation (7) in two aspects:
(1) | Equation (7) does not require | ||||
(2) | The |
Reconstruction
To solve Equation (8), should first be estimated. As stated earlier,
is a vector of weights to penalize the element of
towards zeros. Obviously, stationary voxels should be penalized more heavily than moving voxels. We observe that the variance of moving voxels is usually larger than stationary voxels. Therefore, the variance of one voxel from first frame to previous frame indicates whether the voxel is stationary or moving. So, it is reasonable to use variance of the voxels to derive an estimation for
. To suppress the impact of very early frames, only
frames earlier than the frame being recovered are used to calculate
. In addition, the value of
is normalized to (0, 1] to prevent amplifying the noise. For the first frame and because no variance is available, the CS-MRI [Citation17] is used to reconstruct the image from under-sampled k-space data. Starting from the second frame,
is calculated using variance of the voxels, which will weight adaptively to catch the signal structure from a priori knowledge.
Given , Equation (8) is within the framework of the compressive sensing (CS). Among other requirements, the CS requires an L1 solver. In this work, the Gradient Projection for Sparse Reconstruction with Barzilai and Borwein approach (GPSR-BB) [Citation18] version 5.0 is used to recover the difference out of Equation (8). It is initialized with 0 and stops when the norm of the difference between two consecutive estimates divided by the norm of one of them falls below 10−5.
Methods
The proposed model is implemented in the CS framework and is verified by using a simulated phantom and cardiac MRI cine. In this work, the k-space is under-sampled using a Gaussian random pattern along the phase-encoding direction. Because our model reconstructs current frame by reconstructing the frame difference with previous frame, the reconstruction quality of the first frame will provide a good foundation for subsequent frames. During acquisition, a denser sampling pattern is implemented for the earlier frames and less data was collected in the later ones. Specifically, the first frame was under-sampled by a factor of two, the second frame by four and the rest of the frames by eight. And is the normalized variance of
frames earlier than the frame being recovered because
yields better results overall. The state-of-the-art online method [Citation9] and the offline method k–t FOCUSS [Citation3] were compared with the proposed method, using the same under-sampling patterns. Theoretically, offline methods demonstrate a better performance than online ones because they use the information from all under-sampled data when recovering one frame. All simulations and reconstructions are implemented in MATLAB R2010b (MathWorks, Inc., Natick, MA).
Simulated phantom
A numerical phantom data-set was designed to reflect overall image structure in an axial MRI image of the heart. This data-set was generated using the same method as [Citation9] by constructing a dynamic image series containing three pairs of concentric circles with a fixed radius with only changing grey level, a slowly oscillating radius and a rapidly oscillating radius, respectively. The circle with fixed radius simulates imaging background, and the slowly oscillating circle simulates respiration, while the rapidly oscillating circle simulates heartbeat. These image series are then transformed by the Fourier transform and are under-sampled randomly in the phase-encoding direction. This data-set contains 40 frames of image with the size of 128 × 128. For this simulated data-set, we set ,
for Equation (8).
Cardiac MRI cine
For cardiac cine data, 25 frames of full k-space data were acquired using a 1.5T Philips scanner with a 345 × 270 mm FOV and the matrix size is 256 × 256. The data was obtained using a steady-state free precession sequence with a flip angle of 50 degrees and TR = 3.45 ms. This data-set covered one cardiac phase of the diastole and systole. In this work, the full k-space reconstruction is treated as ‘ground truth’. For this cardiac data-set, we set ,
for Equation (8).
Results
Simulated phantom
Figure (a) shows the first frame of the simulated phantom. The circle at the bottom-right is stationary with only the grey level (concentration) changed; the one near the centre of the FOV is slowly oscillating, while the top-left is rapidly oscillating. We will inspect the area enclosed by dashed rectangle to verify how well the methods catch the rapid oscillation. For this simulated data-set, because the data is generated periodically and fits the assumption of k–t FOCUSS which enforces periodic motion, the k–t FOCUSS outperforms the Kalman and proposed method for most frames.
Figure 3. One frame of simulated phantom (a) with a rectangle ROI and the Gaussian random under-sampling pattern along phase-encoding direction (b).
![Figure 3. One frame of simulated phantom (a) with a rectangle ROI and the Gaussian random under-sampling pattern along phase-encoding direction (b).](/cms/asset/5e37bea9-32b0-4307-bf9f-cef029176d95/gipe_a_1141207_f0003_oc.gif)
Figure shows results of one read-out position crossing the centre of a rapidly oscillating circle. All three methods capture overall dynamics of the motion, while all smooth out the minor grey level change at frames 6 and 31 marked by two arrows in Figure (b). The k–t FOCUSS demonstrates a better performance than the KF and proposed method as expected. However, k–t FOCUSS exhibits periodic pattern near frame 20 in recovered images, because it minimizes the L1 norm of the x-f space using FOCUSS to enforce periodic pattern in x–t space. The absolute error maps show that proposed method slightly out-performs KF, especially in the area of bright inner circle (marked by arrows in Figure (d and e)).
Figure 4. The dashed line in (a) shows one read-out position across the rapidly oscillating circle, (b) is the truth of the x–t space at the read-out position, and (c,d,e) is the recovered x–t space(left half) and their absolute error maps(right half) of the k-t FOCUSS, KF and proposed method.
![Figure 4. The dashed line in (a) shows one read-out position across the rapidly oscillating circle, (b) is the truth of the x–t space at the read-out position, and (c,d,e) is the recovered x–t space(left half) and their absolute error maps(right half) of the k-t FOCUSS, KF and proposed method.](/cms/asset/6caa828b-92a7-44fb-ba45-62a61458b6f6/gipe_a_1141207_f0004_oc.gif)
To evaluate how the proposed method catches slowly oscillating circle and stationary circle with only the grey level changed, we inspected another read-out position covering both slowly oscillating and stationary circles, as shown in Figure . From this figure, we can see that all three methods catch slowly oscillating circle well. However, for stationary circle with changing grey level, these methods fail to catch the variation. The k–t FOCUSS smoothens out almost all grey level changes, especially for frames 6 and 31 marked by two arrows in Figure (c). It can be seen from the absolute error map that proposed method achieved a smaller error than the KF, especially in areas marked by arrows in Figure (d and e).
Figure 5. The dashed line in (a) shows one read-out position across both the slowly oscillating and the stationary circles, (b) is the truth of the x–t space at the read-out position, and (c,d,e) are the recovered x–t space(left half) and their absolute error maps(right half) of k–t FOCUSS, KF and the proposed method.
![Figure 5. The dashed line in (a) shows one read-out position across both the slowly oscillating and the stationary circles, (b) is the truth of the x–t space at the read-out position, and (c,d,e) are the recovered x–t space(left half) and their absolute error maps(right half) of k–t FOCUSS, KF and the proposed method.](/cms/asset/5670d4d2-9c49-4006-b3b6-06e850e44a62/gipe_a_1141207_f0005_oc.gif)
To show the overall error of recovered images, Figure (a) shows the normalized mean squared error () of rapidly oscillating circle as marked by dashed rectangle in Figure (a). It can be seen that proposed method outperforms KF for all frames but 23 and 24. The k–t FOCUSS has a large NMSE for frame 31 because it cannot catch grey level changes of the inner bright circle as shown in the Figure (d).
Cardiac MRI cine
Figure (a) shows one frame of the cardiac cine, in addition to the region of interest (ROI) covering both left and right ventricle cavities.
Figure 7. One frame of cardiac cine (a) with a rectangle ROI and the Gaussian random under-sampling pattern (b).
![Figure 7. One frame of cardiac cine (a) with a rectangle ROI and the Gaussian random under-sampling pattern (b).](/cms/asset/8eceefc6-2282-4af0-88d1-737f5ed6a960/gipe_a_1141207_f0007_oc.gif)
Figure (c) shows that k–t FOCUSS smoothens out some details (marked by black arrows) of the motion, because it enforces periodic pattern in its imaging model. This causes the k–t FOCUSS does not capture the motion of cardiac systole corresponding frames of 9–14. Visually, KF and proposed method had very similar results, except the areas marked by white arrows in Figure (d and e). In addition, the KF had a constant stripe in the position marked by the black arrow in Figure (d), while the proposed method had a similar pattern at the beginning of the imaging (corresponding to a cardiac diastole with a slow motion activity).This gradually disappears as time passes, as shown in Figure (e).This is because the KF updates weights using only the under-sampled Fourier transform and constant covariance matrices of the system and measurement noise, while proposed method corrects weights adaptively, using previous frames to incorporate the motion information.
Figure 8. (a) shows one read-out position across area of rapid cardiac deformation and (b) is the ground truth of the x–t space at the read-out position,(c,d,e) show the recovered x–t space and their error maps of the three methods.
![Figure 8. (a) shows one read-out position across area of rapid cardiac deformation and (b) is the ground truth of the x–t space at the read-out position,(c,d,e) show the recovered x–t space and their error maps of the three methods.](/cms/asset/16bbed72-7078-4d46-964b-c21b6b835559/gipe_a_1141207_f0008_oc.gif)
Figure shows the overall NMSE and SSIM (Structural SIMilarity) [Citation20] of the ROI. It can be seen that the k–t FOCUSS had a smaller NMSE at the beginning and end of imaging (corresponding to a cardiac diastole), while KF and proposed method had smaller errors in the middle of the imaging (corresponding to a cardiac systole). This may also be observed in SSIM of Figure (e). Compared with KF, proposed method had a smaller NMSE for all frames but 13. However, the SSIM reveals that proposed method outperforms the KF gradually as time goes by. This verifies that the proposed model absorbs the motion information gradually to improve the reconstruction.
Figure 9. The NMSE (a) and SSIM (e) of the ROI and the ground truth, recovered images and their error maps of three methods for frames 8(b), 13(c), 24(d), respectively.
![Figure 9. The NMSE (a) and SSIM (e) of the ROI and the ground truth, recovered images and their error maps of three methods for frames 8(b), 13(c), 24(d), respectively.](/cms/asset/03584829-0e6f-47fd-a71d-addadf7e6b92/gipe_a_1141207_f0009_oc.gif)
Figure compares how three methods catch the motion of the heart. The k–t FOCUSS catches the motion better in the early and late stages of cardiac cycle, corresponding to slow motion of cardiac wall; however, this method offers a limited capability in depicting rapid contraction of the ventricle wall (see the 7–13 frames in Figure (b)). The Kalman and proposed method have a very similar error and follow the ground truth well (see Figure (c)). In comparison, without any restrictions in the periodic motion model set-up, proposed method can more accurately capture the motion activity over entire cardiac cycle (see Figure (d)).
Discussion
In this work, a new online model was developed for rapid dynamic MRI using sub-sampling of the Cartesian k-space data. Previously, dynamic MRI has been investigated using the KF method; however, KF requires the system noise to be Gaussian, which is not suitable for modelling dynamic MRI. In dynamic MRI, the system noise is associated with the frame difference, that is, it contains the motion pattern that is usually not in the form of Gaussian distribution. With a basic assumption of measurement noise is time-invariant and white, the KF technique implements a weighted least square optimization to update the weighting matrix. Noticing that the difference of adjacent frames is sparse and the weighting scheme of KF cannot follow the motion of imaging target, we set up an elastic-net model using L1 norm to promote sparsity and weighted L2 norm to catch the motion adaptively.
In the KF method, the weighting matrix is determined by under-sampled Fourier transform, covariance of system and measurement noise, which are chosen to be constant all the time and are updated independently with state variables of the model. While in the proposed method, the weighting matrix is calculated according to the variance of each voxels. As a result, the voxels with larger variance will get reconstructed under sparsity constraints, while voxels with smaller variance will be dominated by previous frame. From comparative results shown above, we can see that proposed method offers a better performance in terms of depicting cardiac dynamics.
The proposed model is solved within the framework of compressive sensing, which benefits from several efficient algorithms. Nonetheless, the Cartesian trajectory enables concurrent reconstruction of all x–t images along read-out direction, which can further reduce the reconstruction time of a frame to only one column by allocating one unit of parallel computing hardware to each read-out position, thus making real-time imaging possible.
The limitation of the proposed model comes from three aspects. First, the model estimates current frame from previous frame using , yielding the frame difference that is not Gaussian. If current frame can be better estimated from previous frames using
, where
is a matrix for motion estimation, to make the frame difference be nearly Gaussian, KF can be applied naturally. However, the motion estimation will incur other possibilities and complexity. Second, it is observed that, the current estimation of adaptive weights by calculating the variance is intuitive and straightforward; however, it may be not robust against the background noise. Third, the proposed model does not accurately represent the real cardiac motion and deformation. And the parameters in the model may contain parametric errors,[Citation21,22] which may be suppressed by the model selection statistically.
Conclusion
This paper presented an online dynamic cardiac MRI model, which can be used to reconstruct image frames from partial acquisition of the k-space data by combining motion-adaptive, sparse constraints. Simulation results showed, compared with online KF method, the proposed method offered a better solution in depicting simulated phantom and cardiac dynamics. Future work will seek better methods for motion estimation and adding a prior structural information for higher spatial and temporal resolution. In addition, the efficiency of the proposed imaging method will be enhanced using hardware-based parallel-computing scheme, which is critical for dynamic MRI studies.
Acknowledgment
This work was supported by National Key Lab of Fundamental Science of Micro/Nano-Device[2014MS01] and Fundamental Research Funds for the Central Universities[CDJZR12090003]
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- Tsao J, Kozerke S. MRI temporal acceleration techniques. J. Magn. Reson. Imaging. 2012;36:543–560.10.1002/jmri.v36.3
- Jung H, Sung K, Nayak KS, et al. k-t FOCUSS: A general compressed sensing framework for high resolution dynamic MRI. Magn. Reson. Med.. 2009;61:103–116.10.1002/mrm.v61:1
- Jung H, Ye JC, Kim EY. Improved k-t BLAST and k-t SENSE using FOCUSS. Phys. Med. Biol. 2007;52:3201–3226.10.1088/0031-9155/52/11/018
- Asif MS, Hamilton L, Brummer M, et al. Motion-adaptive spatio-temporal regularization for accelerated dynamic MRI. Magn. Reson. Med.. 2013;70:800–812.10.1002/mrm.24524
- Bilen C, Wang Y, Selesnick I. Compressed sensing for moving imagery in medical imaging. arXiv:1203.5772.
- Vaswani N. Compressive sensing on Kalman filtered residual. arXiv:0912.1628v3.
- Sümbül U, Santos JM, Pauly JM. A practical acceleration algorithm for real-time imaging. IEEE Trans. Med. Imaging. 2009;28:2042–2051.10.1109/TMI.2009.2030474
- Sümbül U, Santos JM, Pauly JM. Improved time series reconstruction for dynamic magnetic resonance imaging. IEEE Trans. Med. Imaging. 2009;28:1093–1104.10.1109/TMI.2008.2012030
- Feng X, Salerno M, Kramer CM, et al. Kalman filter techniques for accelerated Cartesian dynamic cardiac imaging. Magn. Reson. Med.. 2013;69:1346–1356.10.1002/mrm.v69.5
- Mingjian H, Feng L, Xiaohong Z, et al. Sparsity-enforced Kalman filter technique for dynamic cardiac imaging. Proceedings of the 21stAnnual Meeting ISMRM. Salt Lake City (UT), 2013. Abstract 4698.
- Hui Z, Trevor H. Regularization and variable selection via the elastic net. J. Royal Stat. Soc. Ser. B. 2005;67:301–320.
- Majumdar A, Ward RK. Causal dynamic MRI reconstruction via nuclear norm minimization. Magn. Reson. Imaging. 2012;30:1483–1494.10.1016/j.mri.2012.04.012
- Majumdar A, Ward RK, Aboulnasr T. Compressed sensing based on real-time dynamic MRI reconstruction.IEEE Trans. Med. Imaging. 2012;31:2253–2266.10.1109/TMI.2012.2215921
- Majumdar A. Motion-predicted online dynamic MRI reconstruction from partially sampled k-space data. Magn. Reson. Imaging. 2013;31:1578–1586.10.1016/j.mri.2013.06.003
- Gorodnitsky IF, Rao BD. Sparse signal reconstruction from limited data using FOCUSS: re-weighted minimum norm algorithm. IEEE Trans. Signal Process.;45:600–616.
- Welchand G, Bishop G. An introduction to the Kalman Filter, TR 95-041. Chapel Hill (NC), Department of Computer Science, University of North Carolina.
- Lustig M, Donoho D, Pauly JM. Sparse MRI: the application of compressed sensing for rapid MR imaging. Magn. Reson. Med.. 2007;58:1182–1195.10.1002/(ISSN)1522-2594
- Mario AT, Figueiredo RD, Nowak RD, et al. Gradient projection for sparse reconstruction: Application to compressed sensing and other inverse problems. J. Sel. Top. Signal Process. 2007;1:586–597.
- Sorenson HW. Least-squares estimation: from Gauss to Kalman. IEEE Spectr. 1970;7:63–68.10.1109/MSPEC.1970.5213471
- Wang Z, Bovik AC, Sheikh HR, et al. Image quality assessment: from error visibility to structural similarity. IEEE Trans. Image Process.. 2004;13:600–612.10.1109/TIP.2003.819861
- Peruzzi NJ, Chavarette FR, Balthazar JM, et al. The dynamic behavior of a parametrically excited time-periodic MEMS taking into account parametric errors. J. Vib. Control. 2015;1077546315573913:1–10.
- Balthazar JM, Tusset AM, Bueno AM. TM-AFM nonlinear motion control with robustness analysis to parametric errors in the control signal determination. J. Theor. Appl. Mech. 2014;52:93–106.