Titelaufnahme

Titel
Numerical treatment of a Mumford-Shah model for x-ray tomography / vorgelegt von Elena Hoetzl
Verfasser/ VerfasserinHoetzl, Elena
Begutachter / BegutachterinRing, Wolfgang ; Kappel, Franz
Erschienen2009
UmfangVIII, 109 Bl. : Zsfassung ; Ill., graph. Darst.
HochschulschriftGraz, Univ., Diss., 2009
SpracheEnglisch
Bibl. ReferenzOeBB
DokumenttypDissertation
Schlagwörter (DE)Röntgentomographie / Mumford-Shah Funktional / Inverse Probleme / Form Optimierung / Level Set Methode / Geometrische Sensitivitätsanalyse
Schlagwörter (EN)X-ray tomography / Mumford-Shah functional / Inverse problems / Shape optimization / Level set method / Shape sensitivity analysis
Schlagwörter (GND)Mumford-Shah-Funktional / Tomografie / Sensitivitätsanalyse / Mumford-Shah-Funktional / Tomografie / Sensitivitätsanalyse / Online-Publikation
URNurn:nbn:at:at-ubg:1-217 Persistent Identifier (URN)
Zugriffsbeschränkung
 Das Werk ist frei verfügbar
Dateien
Numerical treatment of a Mumford-Shah model for x-ray tomography [1.68 mb]
Links
Nachweis
Klassifikation
Zusammenfassung (Deutsch)

The goal of this work is to identify a density function of a physical body from a given data as the result of X-rays traveling through the body under different angles and offsets. The mathematical relation between parameter and data is described by the Radon transform. More specifically, the aim of the reconstruction is to find the singularity set of the density --- the collection of lines across which the density is discontinuous --- and the smooth density distribution on the complement of the singularity set.Mumford-Shah models are designed to extract simultaneously functional and geometric information for inaccessible parameters from indirect measurements. We propose a piecewise smooth Mumford-Shah model for the inversion of the tomography data. In our approach the functional variable is eliminated by solving a classical variational problem for each fixed geometry. The solution is then inserted in the Mumford-Shah cost functional leading to a geometrical optimization problem for the singularity set. The resulting shape optimization problem is solved using a shape sensitivity calculus and a propagation of shape variables in the level-set form.A specific difficulty poses the solution of the optimality system for the fixed geometry, which has the form of a coupled system of integro-differential equations on variable and irregular domains.A finite difference method based approach for the determination of a piecewise smooth density function as thesolution of the optimality system is presented. A second order accurate finite difference discretization is proposed.Here a standard five-point stencil is used on regular points of the underlying uniform grid and modificationsof the standard stencil are made at points close to the contour. The optimality system is solved iteratively.The Fast Multipole Method is used to solve domain integrals and numerical experiments are presented.

Zusammenfassung (Englisch)

The goal of this work is to identify a density function of a physical body from a given data as the result of X-rays traveling through the body under different angles and offsets. The mathematical relation between parameter and data is described by the Radon transform. More specifically, the aim of the reconstruction is to find the singularity set of the density --- the collection of lines across which the density is discontinuous --- and the smooth density distribution on the complement of the singularity set.Mumford-Shah models are designed to extract simultaneously functional and geometric information for inaccessible parameters from indirect measurements. We propose a piecewise smooth Mumford-Shah model for the inversion of the tomography data. In our approach the functional variable is eliminated by solving a classical variational problem for each fixed geometry. The solution is then inserted in the Mumford-Shah cost functional leading to a geometrical optimization problem for the singularity set. The resulting shape optimization problem is solved using a shape sensitivity calculus and a propagation of shape variables in the level-set form.A specific difficulty poses the solution of the optimality system for the fixed geometry, which has the form of a coupled system of integro-differential equations on variable and irregular domains.A finite difference method based approach for the determination of a piecewise smooth density function as thesolution of the optimality system is presented. A second order accurate finite difference discretization is proposed.Here a standard five-point stencil is used on regular points of the underlying uniform grid and modificationsof the standard stencil are made at points close to the contour. The optimality system is solved iteratively.The Fast Multipole Method is used to solve domain integrals and numerical experiments are presented.