A modified fuzzy C means algorithm for shading correction in craniofacial CBCT images

CBCT images suffer from acute shading artifacts primarily due to scatter. Numerous image-domain correction algorithms have been proposed in the literature that use patient-specific planning CT images to estimate shading contributions in CBCT images. However, in the context of radiosurgery applications such as gamma knife, planning images are often acquired through MRI which impedes the use of polynomial fitting approaches for shading correction. We present a new shading correction approach that is independent of planning CT images. Our algorithm is based on the assumption that true CBCT images follow a uniform volumetric intensity distribution per material, and scatter perturbs this uniform texture by contributing cupping and shading artifacts in the image domain. The framework is a combination of fuzzy C-means coupled with a neighborhood regularization term and Otsu's method. Experimental results on artificially simulated craniofacial CBCT images are provided to demonstrate the effectiveness of our algorithm. Spatial non-uniformity is reduced from 16% to 7% in soft tissue and from 44% to 8% in bone regions. With shading-correction, thresholding based segmentation accuracy for bone pixels is improved from 85% to 91% when compared to thresholding without shading-correction. The proposed algorithm is thus practical and qualifies as a plug and play extension into any CBCT reconstruction software for shading correction.Comment: 15 pages, published in CMBEBIH 201

A 3D environment for surgical planning and simulation

The use of Computed Tomography (CT) images and their three-dimensional (3D) reconstruction has spread in the last decade for implantology and surgery. A common use of acquired CT datasets is to be handled by dedicated software that provide a work context to accomplish preoperative planning upon. These software are able to exploit image processing techniques and computer graphics to provide fundamental information needed to work in safety, in order to minimize the surgeon possible error during the surgical operation. However, most of them carry on lacks and flaws, that compromise the precision and additional safety that their use should provide. The research accomplished during my PhD career has concerned the development of an optimized software for surgical preoperative planning. With this purpose, the state of the art has been analyzed, and main deficiencies have been identified. Then, in order to produce practical solutions, those lacks and defects have been contextualized in a medical field in particular: it has been opted for oral implantology, due to the available support of a pool of implantologists. It has emerged that most software systems for oral implantology, that are based on a multi-view approach, often accompanied with a 3D rendered model, are affected by the following problems: unreliability of measurements computed upon misleading views (panoramic one), as well as a not optimized use of the 3D environment, significant planning errors implied by the software work context (incorrect cross-sectional planes), and absence of automatic recognition of fundamental anatomies (as the mandibular canal). Thus, it has been defined a fully 3D approach, and a planning software system in particular, where image processing and computer graphic techniques have been used to create a smooth and user-friendly completely-3D environment to work upon for oral implant planning and simulation. Interpolation of the axial slices is used to produce a continuous radiographic volume and to get an isotropic voxel, in order to achieve a correct work context. Freedom of choosing, arbitrarily, during the planning phase, the best cross-sectional plane for achieving correct measurements is obtained through interpolation and texture generation. Correct orientation of the planned implants is also easily computed, by exploiting a radiological mask with radio-opaque markers, worn by the patient during the CT scan, and reconstructing the cross-sectional images along the preferred directions. The mandibular canal is automatically recognised through an adaptive surface-extracting statistical-segmentation based algorithm developed on purpose. Then, aiming at completing the overall approach, interfacing between the software and an anthropomorphic robot, in order to being able to transfer the planning on a surgical guide, has been achieved through proper coordinates change and exploiting a physical reference frame in the radiological stent. Finally, every software feature has been evaluated and validated, statistically or clinically, and it has resulted that the precision achieved outperforms the one in literature

Upper airways segmentation using principal curvatures

Esta tesis propone una nueva técnica para segmentar las vías aéreas superiores. Esta propuesta permite la extracción de estructuras curvilíneas usando curvaturas principales. La propuesta permite la extracción de éstas estructuras en imágenes 2D y 3D. Entre las principales novedades se encuentra la propuesta de un nuevo criterio de parada en la propagación del algoritmo de realce de contraste (operador multi-escala de tipo sombrero alto). De la misma forma, el criterio de parada propuesto es usado para detener los algoritmos de difusión anisotrópica. Además, un nuevo criterio es propuesto para seleccionar las curvaturas principales que conforman las estructuras curvilíneas, que se basa en los criterios propuestos por Steger, Deng et. al. y Armande et. al. Además, se propone un nuevo algoritmo para realizar la supresión de nomáximos que permite reducir la presencia de discontinuidades en el borde de las estructuras curvilíneas. Para extraer los bordes de las estructuras curvilíneas, se utiliza un algoritmo de enlace que incluye un nuevo criterio de distancia para reducir la aparición de agujeros en la estructura final. Finalmente, con base en los resultados obtenidos, se utiliza un algoritmo morfológico para cerrar los agujeros y se aplica un algoritmo de crecimiento de regiones para obtener la segmentación final de las vías respiratorias superiores.This dissertation proposes a new approach to segment the upper airways. This proposal allows the extraction of curvilinear structures based on the principal curvatures. The proposal allows extracting these structures from 2D and 3D images. Among the main novelties is the proposal of a new stopping criterion to stop the propagation of the contrast enhancement algorithm (multiscale top-hat morphological operator). In the same way, the proposed stopping criterion is used to stop the anisotropic diffusion algorithms. In addition, a new criterion is proposed to select the principal curvatures that make up the curvilinear structures, which is based on the criteria proposed by Steger, Deng et. al. and Armande et. al. Furthermore, a new algorithm to perform the non-maximum suppression that allows reducing the presence of discontinuities in the border of curvilinear structures is proposed. To extract the edges of the curvilinear structures, a linking algorithm is used that includes a new distance criterion to reduce the appearance of gaps in the final structure. Finally, based on the obtained results, a morphological algorithm is used to close the gaps and a region growing algorithm to obtain the final upper airways segmentation is applied.Doctor en IngenieríaDoctorad