DISS. ETH NO. 23491 Effective Joint Regularization in Medical Imaging for Nonsmooth Nonlinear Priors A thesis submitted to attain the degree of DOCTOR OF SCIENCES of ETH ZURICH (Dr. sc. ETH Zurich) presented by Valery Vishnevskiy MSc in Applied Mathematics and System Programming Lomonosov Moscow State University born March 21, 1989 citizen of Russia accepted on the recommendation of Prof. Dr. Gábor Székely, examiner Prof. Dr. Simon Arridge, co-examiner Prof. Dr. Sebastian Kozerke, co-examiner Dr. Christine Tanner, co-examiner 2016 Abstract This thesis describes methods for medical image registration, denoising and anatomy reconstruction. We develop and investigate a broad range of regularization approaches to improve accuracy of the provided estimates. The formulated problems are then solved within an optimization framework based on the method of multipliers. Smoothness priors are conveniently used in medical imaging, as they provide tractable formulation of the optimization problem. However they generally reduce the resolution of the estimate due to the oversmoothing of the initial problem. More flexible regularizations (such as total variation) do allow discontinuities, but at the cost of complicated numerical optimization schemes. In this work we exploit recent advances in mathematical optimization to allow efficient implementation of such flexible regularizations by adopting the method of multipliers for medical imaging tasks. To further improve plausibility of the provided estimates, we develop and present regularizations that are aimed to impose prior knowledge in a joint fashion. We employ and adapt different formulations of the total variation measure to tackle problems of image registration and denoising. Transformation consistency over loops is then used to determine the correctness of registrations in the cases of linear and nonlinear displacement fields, yielding a novel consistency models. A sparsity inducing regularization is employed to account for possibly uncompensated cardiac motion for the denoising of in vivo diffusion weighted images. Specialized regularization was developed to jointly account for the helical myofiber architecture and the uncertainty of the diffusion estimates. All methods have been evaluated on clinical images from different modalities (CT, MRI and DW-MRI) and different anatomies (thorax, abdomen, cardiac, corpus callosum). The accuracy of the resulting estimates was shown to be improved compared to the generic and specialized state of the art methods. Zusammenfassung Die vorliegende Arbeit beschreibt Methoden zur Registrierung und Rauschreduktion sowie Verfahren zur anatomischen Rekonstruktion von medizinischen Bildern. Eine Reihe von Regularisierungsansätzen wird vorgestellt, um die Genauigkeit der verfügbaren Schätzungen zu verbessern. Die Problemstellungen sind im Rahmen der Theorie der mathematischen Optimierung formuliert und werden mit Hilfe des Verfahrens der Multiplikatoren gelöst. Es ist in der medizinischen Bildgebung weit verbreitet, a-priori eine Glattheit der Lösung anzunehmen. Diese Annahme erlaubt es, die resultierenden Probleme effizient zu lösen, führt aber häufig durch Überglättung zu einer qualitativen Verminderung der Ergebnisse. Universellere Regularisierungen (wie zum Beispiel die Total-Variation-Regularisierung) ermöglichen es Diskontinuitäten auf Kosten von Algorithmen mit erhöhter Berechnungskomplexität zu modellieren. In der vorliegenden Arbeit verwenden wir die neuesten Erkenntnisse auf dem Gebiet der mathematischen Optimierungsmethoden mit dem Ziel, solche universellen Regularisierungen in der medizinischen Bildgebung mittels des sogenannten Verfahrens der Multiplikatoren wirksam umzusetzen. Zur weiteren Verbesserung der Vorhersagewahrscheinlichkeit entwerfen wir Regularisierungen, welche Vorkenntnisse pauschal einbeziehen. Wir verwenden verschiedene Formulierungen der Total-Variation-Regularisierung und adaptieren diese für die Rauschreduktion und Registrierung von medizinischen Bildern. Der residuale Fehler von Registrierungs-Zyklen über mehrere Bilder wird verwendet, um Fehler einzelner, paarweise linearer und nicht linearer Registrierungen abzuschätzen, was zu einem neuen Konsistenzmodel führt. Des Weiteren wird eine sparsity-fördernde Regularisierung verwendet, um möglicherweise nicht kompensierte Herzbewegung während der Rauschreduktion von in vivo Diffusions-Kardiogrammen zu berücksichtigen. Eine besondere Regularisierung wurde entwickelt, welche sowohl die Struktur der Herzmuskelfasern als auch die Unsicherheit der Diffusionsschätzungen berücksichtigt. Die vorgeschlagenen Methoden wurden für verschiedene medizinische Bildmodalitäten (CT, MRI and DW-MRI) und anatomische Regionen (Brust, Bauchraum, Herz, Corpus callosum) getestet. Wir zeigen, dass die hier vorgeschlagenen Methoden die Genauigkeit der Ergebnisse im Vergleich sowohl zu universellen Verfahren, als auch zu spezialisierten Methoden deutlich erhöhen.
© Copyright 2024 ExpyDoc