Return to search

On the Application of Compressed Sensing to Magnetic Resonance Imaging / Über die Anwendung von Compressed Sensing in der Magnetresonanztomographie

This thesis investigated the potential of Compressed Sensing (CS) applied to Magnetic Resonance Imaging (MRI). CS is a novel image reconstruction method that emerged from the field of information theory. The framework of CS was first published in technical reports in 2004 by Candès and Donoho. Two years later, the theory of CS was published in a conference abstract and two papers. Candès and Donoho proved that it is possible, with overwhelming probability, to reconstruct a noise-free sparse signal from incomplete frequency samples (e.g., Fourier coefficients). Hereby, it is assumed a priori that the desired signal for reconstruction is sparse. A signal is considered “sparse“ when the number of non-zero elements is significantly smaller than the number of all elements. Sparsity is the most important foundation of CS. When an ideal noise-free signal with few non-zero elements is given, it should be understandably possible to obtain the relevant information from fewer Fourier coefficients than dictated by the Nyquist-Shannon criterion. The theory of CS is based on noise-free sparse signals. As soon as noise is introduced, no exact sparsity can be specified since all elements have signal intensities that are non-zero. However, with the addition of little or moderate noise, an approximate sparsity that can be exploited using the CS framework will still be given. The ability to reconstruct noisy undersampled sparse MRI data using CS has been extensively demonstrated. Although most MR datasets are not sparse in image space, they can be efficiently sparsified by a sparsifying transform. In this thesis, the data are either sparse in the image domain, after Discrete Gradient transformation, or after subtraction of a temporally averaged dataset from the data to be reconstructed (dynamic imaging). The aim of this thesis was to identify possible applications of CS to MRI. Two different algorithms were considered for reconstructing the undersampled sparse data with the CS concept. The Nonlinear Conjugate Gradient based technique with a relaxed data consistency constraint as suggested by Lustig et al. is termed Relaxed DC method. An alternative represents the Gradient or Steepest Descent algorithm with strict data consistency and is, therefore, termed the Strict DC method. Chapter 3 presents simulations illustrating which of these two reconstruction algorithms is best suited to recover undersampled sparse MR datasets. The results lead to the decision for the Strict DC method as reconstruction technique in this thesis. After these simulations, different applications and extensions of CS are demonstrated. Chapter 4 shows how CS benefits spectroscopic 19F imaging at 7 T, allowing a significant reduction of measurement times during in vivo experiments. Furthermore, it allows highly resolved spectroscopic 3D imaging in acceptable measurement times for in vivo applications. Chapter 5 introduces an extension of the Strict DC method called CS-CC (CS on Combined Coils), which allows efficient processing of sparse undersampled multi-coil data. It takes advantage of a concept named “Joint Sparsity“, which exploits the fact that all channels of a coil array detect the same sparse object weighted with the coil sensitivity profiles. The practical use of this new algorithm is demonstrated in dynamic radial cardiac imaging. Accurate reconstructions of cardiac motion in free breathing without ECG triggering were obtained for high undersampling factors. An Iterative GRAPPA algorithm is introduced in Chapter 6 that can recover undersampled data from arbitrary (Non-Cartesian) trajectories and works solely in the Cartesian plane. This characteristic makes the proposed Iterative GRAPPA computationally more efficient than SPIRiT. Iterative GRAPPA was developed in a preceding step to combine parallel imaging with CS. Optimal parameters for Iterative GRAPPA (e.g. number of iterations, GRAPPA kernel size) were determined in phantom experiments and verified by retrospectively undersampling and reconstructing a radial cardiac cine dataset. The synergistic combination of the coil-by-coil Strict DC CS method and Iterative GRAPPA called CS-GRAPPA is presented in Chapter 7. CS-GRAPPA allows accurate reconstruction of undersampled data from even higher acceleration factors than each individual method. It is a formulation equivalent to L1-SPIRiT but computationally more efficient. Additionally, a comparison with CS-CC is given. Interestingly, exploiting joint sparsity in CS-CC is slightly more efficient than the proposed CS-GRAPPA, a hybrid of parallel imaging and CS. The last chapter of this thesis concludes the findings presented in this dissertation. Future applications expected to benefit from CS are discussed and possible synergistic combinations with other existing MR methodologies for accelerated imaging are also contemplated. / In der vorliegenden Arbeit wurde untersucht, welches Potential die Anwendung von Compressed Sensing (CS) in der Magnetresonanztomographie (MRT) hat. CS ist eine neue Bildrekonstruktionsmethode aus der Informationstheorie. Das Grundgerüst für CS wurde zuerst in zwei technischen Berichten von Candès und Donoho aus dem Jahr 2004 vorgestellt. Zwei Jahre später wurde die CS-Theorie in einem Konferenzbeitrag und zwei wissenschaftlichen Artikeln veröffentlicht. Candés und Donoho zeigten, dass es mit überwältigender Wahrscheinlichkeit möglich ist, ein rauschfreies sparses Signal aus unvollständig vorliegender Frequenzinformation zu rekonstruieren. Hierfür ist eine wichtige A-priori-Annahme, dass das gewünschte Signal, welches rekonstruiert werden soll, sparse sein soll. Man spricht von sparsen Signalen, falls die Anzahl der Elemente mit Intensität größer Null signifikant kleiner als die Anzahl aller Elemente ist. Die CS-Theorie basiert auf rauschfreien, sparsen Signalen. Sobald Rauschen auftritt, kann keine exakte Sparsity mehr bestimmt werden, da alle Elemente Signalintensitäten größer Null haben. Falls jedoch nur wenig oder moderates Rauschen hinzugefügt wird ist immer noch näherungsweise eine Sparsity gegeben, die mit Hilfe von CS ausgenutzt werden kann. Die meisten MR-Datensätze sind nicht-sparse im Bildraum, können allerdings durch eine sog. Sparsifizierungstransformation effektiv sparsifiziert werden. In der vorliegenden Arbeit sind die Daten entweder im Bildraum sparse, nach einer Diskreten-Gradienten-Transformation oder nachdem bei dynamischen Daten ein zeitlich gemittelter Datensatz von den zu rekonstruierenden Daten abgezogen worden ist. Das Ziel dieser Arbeit war es, mögliche Anwendungen für CS in der MRT zu identifizieren. Zwei unterschiedliche Algorithmen wurden untersucht, um unterabgetastete sparse Daten mit dem CS-Konzept zu rekonstruieren. Eine Technik, die auf einer Nichtlinearen Methode der Konjugierten Gradienten basiert und eine gelockerte Datenkonsistenzbedingung beinhaltet, wird als Relaxed DC-Methode bezeichnet. Eine Alternative stellt der Gradienten- oder Steilster-Abstieg-Algorithmus dar, der strikte Datenkonsistenz fordert und daher als Strict-DC-Methode bezeichnet wird. Kapitel 3 zeigt Simulationen, die darlegen, dass die Strict-DC-Methode am besten zur Datenrekonstruktion in dieser Arbeit geeignet ist. Kapitel 4 zeigt, in wie fern die spektroskopische 19F-Bildgebung bei 7 T von CS profitieren kann, indem CS eine signifikante Reduktion der Messzeiten bei in vivo Experimenten erlaubt. Desweiteren ermöglicht CS hochaufgelöste spektroskopische 3D-Bildgebung in akzeptablen Messzeiten für in vivo Anwendungen. Kapitel 5 führt eine Erweiterung der Strict-DC-Methode ein, die CS-CC genannt wird, welche eine effiziente Bearbeitung von sparsen unterabgetasteten Multi-Empfänger-Datensätzen erlaubt. Hierbei profitiert CS-CC von einem Konzept namens "Joint Sparsity", welches ausnutzt, dass alle Empfangskanäle eines Spulenarrays dasselbe sparse Objekt detektieren, jeweils gewichtet mit den entsprechenden Spulensensitivitätsprofilen. Der praktische Nutzen dieses neuen Algorithmus wird an einem dynamischen radialen Herzdatensatz verdeutlicht. Akkurate Rekonstruktionen der Herzbewegung in freier Atmung und ohne EKG-Trigger konnten bei hohen Unterabtastfaktoren erreicht werden. Ein Iterativer-GRAPPA-Algorithmus, der unterabgetastete Daten beliebiger (nicht-kartesischer) Trajektorien rekonstruieren kann und ausschließlich auf einem kartesischen Gitter arbeitet, wird in Kapitel 6 vorgestellt. Das vorgeschlagene Iterative GRAPPA ist vom Rechenaufwand her effizienter als SPIRiT und wurde als ein vorhergehender Schritt zur Kombination von Paralleler Bildgebung und Compressed Sensing entwickelt. Optimale Parameter für Iteratives GRAPPA (z.B. Anzahl an Iterationen, GRAPPA-Kern-Größe) wurden in Phantom-Experimenten bestimmt und mittels Rekonstruktionen an einem retrospektiv unterabgetasteten radialen Herzdatensatz verifiziert. Die synergetische Kombination der spulenweise angewendeten Strict-DC-Methode und Iterativem GRAPPA genannt CS-GRAPPA wird in Kapitel 7 präsentiert. CS-GRAPPA erlaubt akkurate Rekonstruktionen unterabgetasteter Daten von höheren Beschleunigungsfaktoren, als mit den jeweiligen Einzelmethoden möglich gewesen wäre. Die Formulierung ist äquivalent zu L1-SPIRiT, allerdings vom Rechenaufwand effizienter. Es wurde zusätzlich ein Vergleich zu CS-CC durchgeführt. Interessanterweise hat sich gezeigt, dass das Ausnutzen der Joint Sparsity in CS-CC etwas effizienter ist als das vorgeschlagene CS-GRAPPA, das ein Hybrid aus Compressed Sensing und Paralleler Bildgebung ist. Im abschließenden Kapitel dieser Dissertation werden die Ergebnisse zusammengefasst und Schlussfolgerungen daraus gezogen. Zukünftige Anwendungen werden diskutiert, die von CS profitieren und mögliche synergetische Kombinationen mit anderen existierenden MR-Methoden für beschleunigte Bildgebung werden angesprochen.

Identiferoai:union.ndltd.org:uni-wuerzburg.de/oai:opus.bibliothek.uni-wuerzburg.de:6121
Date January 2011
CreatorsFischer, André
Source SetsUniversity of Würzburg
LanguageEnglish
Detected LanguageEnglish
Typedoctoralthesis, doc-type:doctoralThesis
Formatapplication/pdf
Rightshttps://creativecommons.org/licenses/by-nc-nd/3.0/de/deed.de, info:eu-repo/semantics/openAccess

Page generated in 0.0028 seconds