Research Projects

Connecting Theoretical Tomography with the World of Computers


Involved researchers from our group: Richard Huber

Description of the research project:
Computed tomography is a crucial method in medicine, allowing doctors to investigate a patient's body's interior. Tomographic reconstruction results in an ill-posed inverse problem whose solution is non-trivial. While theoretical tomography is commonly understood as operator equations in function spaces, in practice, only finitely many measurements can be made, and finitely many calculations can be performed to process them.

Hence, in practice, tomographic reconstruction is performed on computers using discretizations, i.e., reduced models of finite dimensions.

This naturally raises the question of whether the discrete problems are related to our continuous understanding of computed tomography. One cannot hope for discrete inverse problems to fully describe the infinite-dimensional setting; however, one would expect that when refining discretizations, e.g., by increasing resolutions, the discretizations are more and more representative of the infinite-dimensional situation.

The investigations move in multiple directions that benefit each other. Quite naturally the discrete forward operator should be as representative of the true operator as possible, a goal that leads to development of better discretization schemes. But that alone is not enough to guarantee that the solutions we obtained for the discrete problems, in fact, converge to solutions of the continuous inverse problem. Investigating for which methods and which senses this is the case is another key aspect of this project.

Associated publications:

[HA2025] Richard Huber. Convergence of Ray- and Pixel-Driven Discretization Frameworks in the Strong Operator Topology. 2025
[BH2021] Bredies, Kristian and Huber, Richard. Convergence Analysis of Pixel-Driven Radon and Fanbeam Transforms. SIAM Journal on Numerical Analysis, 59(3):1399-1432, 2021
[HB2025] Richard Huber. A Novel Interpretation of the Radon Transform's Ray and Pixel-Driven Discretizations Under Balanced Resolutions. Scale Space and Variational Methods in Computer Vision, Springer Nature Switzerland, 2025

Characterization of Plausible Tomography Data


Involved researchers from our group: Richard Huber

Description of the research project:
Tomography is a crucial method in medicine, allowing doctors to investigate a patient's body's interior from a sequence of independent projections (views of the body -- e.g., via X-rays -- from different directions). The tomographic reconstruction process is ill-posed, i.e., it can be very strongly impacted by minor corruptions of the used data, resulting in undesired artifacts. However, medical imaging data obtained in practice are unavoidably corrupted, not only by random noise, but also by more systematic errors caused by patients' motion or imperfect detector setups.

Thus, preprocessing the data into a suitable form is a key step of tomographic reconstruction. To that end, detecting and correcting said corruptions prior to reconstruction is crucial. Doing so requires differentiation between plausible and implausible data.

One way to describe the plausibility of data is by checking whether they are in the range of tomographic projection operators describing the measurement process. The level of deviation from said range can indicate the level of corruption. However, how can we check whether data is in the range without doing reconstruction? Data in the range inherently possesses information overlaps between different projections that systematic corruptions will perturb.

Hence, the goal of this research is the characterization of tomographic projection operators' ranges and their related inherent information overlap. We develop general methods for characterizing the range of tomographic projection operators and describe how they can be employed in tomographic preprocessing.

Associated publications:

[HCD2025] Huber, Richard and Clackdoyle, Rolf and Desbat, Laurent. Determination of Range Conditions for General Projection Pair Operators. 2025
[H2023] Huber, Richard and Clackdoyle, Rolf and Desbat, Laurent. Pairwise data consistency conditions for the exponential fanbeam transform. Conference Proceedings for the 17th International Meeting on Fully Three-Dimensional Image Reconstruction in Radiology and Nuclear Medicine, 2023
Poster on Data consistency conditions < >

Structured model learning


Involved researchers from our group: Martin Holler, Richard Huber, Štěpán Zapadlo, Erion Morina

Description of the research project:

The field of learning PDE-based models from data is growing rapidly. It is driven by the need to bridge physical principles with data-driven insights and enables successful advances in learning unknown parameters in partially specified PDEs, discovering entirely new PDE structures, and approximating solution operators with remarkable accuracy. Alongside practical applications, there is increasing attention to theoretical analysis addressing consistency, convergence, and generalization. These are essential for ensuring the reliability and robustness of learned models in scientific and engineering contexts.

Our research group focuses on the framework of "structured model learning." Structured model learning builds on approximate physical knowledge and addresses the limitations of overly coarse abstractions by incorporating fine-scale physics learned from data. This approach improves interpretability, accuracy, and generality, enabling the model to handle complex scenarios and external factors beyond the scope of simplified physical principles. Augmenting known physics with data-driven components enables the model to effectively describe phenomena, even in non-ideal or challenging contexts. However, it is crucial to learn only what is necessary to ensure that the augmentation remains efficient and grounded in physical understanding.

Our research group is actively engaged in structured model learning, focusing on both the theoretical foundations and the practical applications. We have investigated the unique identifiability of learned fine-scale physics and established a convergence result with practical relevance (see [HM2024]). Currently, a major focus of our work is learning multi-pool dynamics in magnetic resonance imaging (MRI).

Associated publications:

[HM2024] Holler, Martin and Morina, Erion. On uniqueness in structured model learning. 2024
[AHN2023] Christian Aarset and Martin Holler and Tram T. N. Nguyen. Learning-informed parameter identification in nonlinear time-dependent PDEs. Applied Mathematics & Optimization, 88(3):1-53, 2023
Unique Learning < >

Variational Methods for Inverse Problems in Imaging


Involved researchers from our group: Martin Holler, Richard Huber

Description of the research project:

Variational methods are a state-of-the-art approach for solving inverse problems in image processing and beyond. By employing suitable regularization functionals, variational methods enable provably stable and consistent reconstructions for a wide range of inverse imaging problems and effectively compensate for missing data through appropriate mathematical modeling.

Our research group has made several fundamental contributions to the field of variational methods for inverse problems. These include:

  1. the development and analysis of higher-order, non-smooth regularization functionals [BCHKS2025, BCH2022, BH2020, BHSW2018, HHK2018, BH2014],
  2. the creation of non-smooth models for image and video decompression [BH2015a, BH2015b, BH2015c, BH2012], and
  3. the development of novel models for texture data [CHP2020].

Associated publications:

[BH2012] Kristian Bredies and Martin Holler. A total variation-based JPEG decompression model. SIAM Journal on Imaging Sciences, 5(1):366--393, 2012
[BH2014] Kristian Bredies and Martin Holler. Regularization of linear inverse problems with total generalized variation. Journal of Inverse and Ill-Posed Problems, 22(6):871--913, 2014
[BH2015a] Kristian Bredies and Martin Holler. A TGV-based framework for variational image decompression, zooming and reconstruction. Part I: Analytics. SIAM Journal on Imaging Sciences, 8(4):2814-2850, 2015
[BH2015b] Kristian Bredies and Martin Holler. A TGV-based framework for variational image decompression, zooming and reconstruction. Part II: Numerics. SIAM Journal on Imaging Sciences, 8(4):2851--2886, 2015
[HHK2018] Martin Holler and Richard Huber and Florian Knoll. Coupled regularization with multiple data discrepancies. Inverse Problems, 34(8):084003, 2018
[BHSW2018] Kristian Bredies and Martin Holler and Martin Storath and Andreas Weinmann. Total Generalized Variation for manifold-valued data. SIAM Journal on Imaging Sciences, 11(3):1785-1848, 2018
[CHP2020] Antonin Chambolle and Martin Holler and Thomas Pock. A Convex Variational Model for Learning Convolutional Image Atoms from Incomplete Data. Journal of Mathematical Imaging and Vision, 62(3):417-444, Springer, 2020
[BH2020] Kristian Bredies and Martin Holler. Higher-order total variation approaches and generalisations. Inverse Problems. Topical Review, 36(12):123001, 2020
[BCH2022] Kristian Bredies and Marcello Carioni and Martin Holler. Regularization Graphs -- A unified framework for variational regularization of inverse problems. Inverse Problems, 38(10):105006, 2022
[BCHKS2025] Kristian Bredies, Marcello Carioni, Martin Holler, Yury Korolev and Carola-Bibiane Schönlieb. A sparse optimization approach to infinite infimal convolution regularization. Numerische Mathematik, 157:1-96, 2025
[BH2013] Kristan Bredies and Martin Holler. A TGV Regularized Wavelet Based Zooming Model. Scale Space and Variational Methods in Computer Vision, Springer Berlin Heidelberg, 2013
[BH2015c] Kristian Bredies and Martin Holler. Artifact-Free Variational MPEG Decompression. Scale Space and Variational Methods in Computer Vision, Springer, 2015
[BH2012] Kristian Bredies and Martin Holler. A pointwise characterization of the subdifferential of the total variation functional. 2016