Single-Course Engelsk 5 ECTS

# Scientific Computing for X-Ray Computed Tomography

## Overall Course Objectives

X-ray Computed Tomography (CT) is used routinely in medicine, materials science and many other applications to reconstruct an object’s interior using mathematical methods and numerical algorithms.
This course focuses on the formulation, implementation, and use of standard reconstruction methods for CT such as Filtered Back Projection, Algebraic Iterative Reconstruction methods, and regularization methods. We give a rigorous mathematical description of the CT reconstruction problem, the associated mathematical formulations, and the underlying computational algorithms – supplemented with hands-on MATLAB computer exercises that illustrate these methods. In addition, we give exercises in large-scale reconstruction of real CT data using the Python package Core Imaging Library (CIL). The goal is that participants will get a basic understanding of the formulation, implementation, and use of basic and advanced CT reconstruction algorithms, and thus be able to use them to perform data analysis for their own CT problems. As part of the course, participants will acquire their own X-ray CT data at the DTU 3D Imaging Facility and reconstruct it using the methods from the course.

See course description in Danish

Learning Objectives

• Describe the underlying physics of a CT scanner.
• Formulate the corresponding mathematical model and the Radon transform.
• Formulate the inverse Radon transform and the Filtered Back Projection algorithm.
• Discretize the Radon transform to obtain a system of linear equations.
• Use the singular value decomposition (SVD) to analyze the reconstruction problem.
• Formulate and use algebraic iterative methods that include simple constraints.
• Formulate the convergence behavior of these methods.
• Use the software package Core Imaging Library (CIL) for large-scale problems.
• Formulate variational problems based on Bayesian noise modeling.
• Formulate and use Tikhonov regularization and Total Variation regularization.
• Use modern numerical methods from convex optimization for CT reconstruction.

Course Content

Introduction to CT and some of its applications. The CT-scanner. The Radon transform and its inverse, Filtered Back Projection. Discretization of the CT problem. The Singular Value Decomposition (SVD) and its use for studying the CT problem. Stability and the need for filtering; truncated SVD.
Algebraic iterative reconstruction algorithms – foundations and convergence properties. Their behavior for noisy data; semi-convergence and stopping rules. The use of GPU computing. The software package Core Imaging Library (CIL) and its algebraic reconstruction algorithms.
Noise models, priors and regularization. Variational formulations and Bayesian modeling. Cases: Total Variation and Tikhonov regularization. Introduction to convex optimization and numerical optimization algorithms. Artifacts in reconstructions and model calibration.

Recommended prerequisites

Experience with Matlab programming (e.g, 02631/33) and numerical computations (e.g., 02601).

Teaching Method

Seminars and small reports on the computer projects.

Faculty

Limited number of seats

Minimum: 5.

Please be aware that this course will only be held if the required minimum number of participants is met. You will be informed 8 days before the start of the course, whether the course will be held.

See course in the course database.

## Registration

Language Engelsk 3 weeks Compute DTU Lyngby Campus 02946 PhD Week 2 Week 4 Mon-fri 8:00-17:00
Price 10.600,00 DKK