Implementación Paralela de un Método Iterativo de Reconstrucción de Imagen de TAC y Evaluación Paramétrica en el Grid

Abstract

The SART algorithm s goal is to do an iterative CT reconstruction. This iterative reconstruction obtains an aproximate solution and can be done with a set of only a few projections instead of the full projections needed when using direct methods, this way the dose that the patient takes can be reduced.

In this work we want to do a parallelization of the algorithm in C, for that purpose we will use PETSc and CUDA libraries. After this parallelization we will do an study on a grid, in this study we will launch multiple simultaneous tasks, changing these parameters that can be tuned by the user in order to find those parameters that are optimum for different image qualities, number of projections used, etc.

El algoritmo SART tiene como objetivo la reconstrucción iterativa de imágenes de TAC. Esta reconstrucción iterativa obtiene una solución aproximada y puede realizarse a partir de un conjunto menor de proyecciones que los métodos directos, reduciendo así la dosis necesaria que recibe el paciente.

En este trabajo se realizara una implementación paralela del método propuesto en el lenguaje de programación C, para ello utilizaremos librerías de altas prestaciones para matrices dispersas y CUDA. A continuación, se procederá a realizar un estudio en el grid, en el cual se realizarán multitud de pruebas simultáneas para realizar un barrido de los parámetros que puede ajustar el usuario en este algoritmo; con el objetivo de determinar aquellos valores que son óptimos para diferentes calidades de imágenes a reconstruir, número de proyecciones usadas, etc.