Robust Algorithms for Discrete Tomography
Frank Tabak
Daily supervisor: Martin van Gijzen
Supervisor CWI: Joost Batenburg
Site of the project: CWI
Science Park 123
1098XG Amsterdam
Start of the Project: February 2012
In June 2012 the
Interim Thesis
and a
presentation has been given.
To illustrate the algorithm see the following video's:
blob_hole
en
shapes.avi.
The Master project has been finished in November 2012
by the completion of the
Masters Thesis
and a final
presentation
has been given.
For working address etc. we refer to our
alumnipage.
Summary of the master project:
Tomography methods concentrate on reconstructing objects from multiple projections that are obtained by sending, for example, X-rays through the object. Applications of these methods are, among others, radiology (CT-, MRI- and PET scans), geophysics and material science. The tomographic problems can be formulated as a system of linear equations. Unfortunately, these systems are not symmetric or positive (semi)definite and usually rank deficient.
In material science one is often presented with very small objects (like crystals or nano-structures) that consist of one or a small number of different materials, each with its own density. Scanning these small objects can cause damage to the structure and thus one can only take a very limited amount of projections. Fortunately, one can use the prior knowledge about the object to arrive at a reconstruction of the original object. How to arrive at this reconstruction is studied by the field of discrete tomography (DT).
With every kind of tomography, and thus also with DT, one is faced with noisy data. Because of this noise the reconstruction process becomes more difficult since the system of linear equations becomes inconsistent. The DART (Discrete Algebraic Reconstruction Method) algorithm was developed to solve DT problems. This algorithm deals with noise in a very heuristic method. The goal of this project is to investigate how the problem can be regularized such that it deals with the noise in a more efficient and robust manner.
Left: Original image (256 by 256 pixels) which is scanned using 5 projection angles. Right: Reconstructing of the object from the projections. The number of wrong pixels is 21.
Contact information:
Kees
Vuik
Back to the
home page
or the
Master students page of Kees Vuik