s'authentifier
version française rss feed
HAL : ujm-00670025, version 1

Fiche détaillée  Récupérer au format
2011 IEEE Nuclear Science Symposium and Medical Imaging Conference, Valencia : Espagne (2011)
A new representation and projection model for tomography, based on separable B-splines
Fabien Momey 1, 2, Loïc Denis 1, Catherine Mennessier 1, Éric Thiébaut 2, Jean-Marie Becker 1, Laurent Desbat 3
(02/2012)

Data modelization in tomography is a key point for iterative reconstruction. The design of the projector, i.e. the numerical model of projection, is mostly influenced by the representation of the object of interest, decomposed on a discrete basis of functions. Standard projector models are voxel or ray driven; more advanced models such as distance driven, use simple staircase voxels, giving rise to modelization errors due to their anisotropic behaviour. Moreover approximations made at the projection step amplify these errors. Though a more accurate projection could reduce approximation errors, characteristic functions of staircase voxels constitute a too coarse basis for representing a continuous function. As a result, pure modelization errors still hold. Spherically symmetric volume elements (blobs) have already been studied to eradicate such errors, but at the cost of increased complexity, because they require some tuning parameters for adapting them to this use. We propose to use 3D B-splines, which are piecewise polynomials, as basis functions. When the degree of these polynomials is sufficiently high, they are very close from being with a spherical symmetry, i.e. blobs, avoiding projection inconsistencies, while keeping local influence and separability property. B-splines are considered, in sampling theory, as the almost optimal functions for the discretization of a continuous signal, not necessarily band-limited, potentially allowing to reduce the angular sampling of the data without any loss of quality. We show that the projection of B-splines can be approximated rather accurately by a separable function, independent from the angle of projection, easier to integrate on detector pixels. The higher the degree of the used B-splines, the better the quality of the approximation, but also the larger the number of required operations. Thanks to these approximations, a convenient tradeoff between the need of accuracy and a fast calculation can be obtained. This has resulted in the implementation of a more accurate numerical projector, which can deal with a reduced angular sampling without loss of performance. The additional computation cost is also efficiently reduced. We have studied the quality of enhancement involved by this projector on 2D iterative reconstructions of a Shepp-Logan phantom, from a small number of fan beam projections. Reconstructions have been performed by optimization methods, minimizing the squared data residuals with a regularization term, using an efficient Quasi-Newton optimization algorithm.
1 :  LAboratoire Hubert Curien [Saint Etienne] (LHC)
CNRS : UMR5516 – Université Jean Monnet - Saint-Etienne
2 :  Centre de Recherche Astrophysique de Lyon (CRAL)
CNRS : UMR5574 – INSU – Université Claude Bernard - Lyon I – École Normale Supérieure (ENS) - Lyon
3 :  Techniques de l'Ingénierie Médicale et de la Complexité - Informatique, Mathématiques et Applications, Grenoble (TIMC-IMAG)
CNRS : UMR5525 – Université Joseph Fourier - Grenoble I
GMCAO
Informatique/Imagerie médicale

Informatique/Traitement du signal et de l'image

Sciences de l'ingénieur/Traitement du signal et de l'image
tomography – iterative reconstruction – inverse problems – B-splines – projector – model
Liste des fichiers attachés à ce document : 
PDF
momey_nssmic2011.pdf(403.7 KB)
ANNEX
momey_poster_nssmic2011.pdf(1.3 MB)