Catálogo de publicaciones - libros

In this paper, we treat the problem of reducing the false positives (FP) in the automatic detection of colorectal polyps at Computer Aided Detection in Computed Tomography Colonography (CAD-CTC) as a shape-filtering task. From the extracted candidate surface, we obtain a reliable shape distribution function and analyse it in the Fourier domain and use the resulting spectral data to classify the candidate surface as belonging to a polyp or a non-polyp class. The developed shape filtering scheme is computationally efficient (takes approximately 2 seconds per dataset to detect the polyps from the colonic surface) and offers robust polyp detection with an overall false positive rate of 5.44 per dataset at a sensitivity of 100% for polyps greater than 10mm when it was applied to standard and low dose CT data.

We examined the echotexture in ultrasonographic images of the wall of dominant ovulatory follicles in women during natural menstrual cycles and dominant anovulatory follicles which developed in women using oral contraceptives (OC). Ovarian follicles in women are fluid-filled structures in the ovary that contain oocytes (eggs). Dominant follicles are physiologically selected for preferential development and ovulation. Statistically significant differences between the two classes of follicles were observed for two co-occurrence matrix derived texture features and two edge-frequency based texture features which allowed accurate distinction of healthy and atretic follicles of similar diameters. Trend analysis revealed consistent turning points in time series of texture features between 3 and 4 days prior to ovulation coinciding with the time at which follicles are being biologically “prepared” for ovulation.

A new approach to fitting pharmacokinetic models to DCE-MRI data is described. The method relies on fitting individual concentration curves to a small set of basis functions and then making use of a look up table to relate the fitting coefficients to pre-calculated pharmacokinetic parameters. This is significantly faster than traditional non-linear fitting methods. Using simulated data and assuming a Tofts model, the accuracy of this direct approach is compared to the Levenberg-Marquardt algorithm. The effect of signal to noise ratio and the number of basis functions used on the accuracy is investigated. The basis fitting approach is slightly less accurate than the traditional non-linear least squares approach but the ten-fold improvement in speed makes the new technique useful as it can be used to generate pharmacokinetic maps in a clinically acceptable timeframe.

In this paper, we present a systematic investigation of the variability of the human cortical folding using a generic gyral model (GGM). The GGM consists of a fixed number of vertices that can be registered non-linearly to an individual anatomy so that for each individual we have a clearly defined set of landmarks that is spread across the cortex. This allows us to obtain a regionalized estimation of inter-subject variability. Since the GGM is stratified into different levels of depth, it also allows us to estimate variability as a function of depth. As another application of a polygonal line representation underlying the generic gyral model, we present a cortical parcellation scheme that can be used to regionalize cortical measurements.

The X-ray imaging equipment could be used to measure hemodynamic function in addition to visualizing the morphology. The parameters of specific interest are arterial blood flow and velocity. Current monoplane X-ray systems can perform 3D reconstruction of the arterial tree as well as to capture the propagation of the injected contrast agent on a sequence of 2D angiograms. We combine the 2D digital subtraction angiography sequence with the mechanically registered 3D volume of the vessel tree. From 3D vessel tree we extract each vessel and obtain its centerline and cross-section area. We get our velocity estimation from 2D sequence by comparing time-density signals measured at different ends of the projected vessel. From the average velocity and cross-section area we get the average blood flow estimate for each vessel. The algorithm described here is applied to datasets from real neuroradiological studies.

Quantitative analysis of computed tomographic (CT) images of the lungs is becoming increasingly useful in the medical and surgical management of subjects with Chronic Obstructive Pulmonary Disease (COPD). Current methods for the assessment of airway wall work well in idealized models of the airway. We propose a new method for airway wall detection based on phase congruency. This method does not rely on either a specific model of the airway or the point spread function of the scanner. Our results show that our method gives a better localization of the airway wall than ”full width at a half max” and is less sensitive to different reconstruction kernels and radiation doses.

We present a novel method for curved planar reformation (CPR) of spine images obtained by magnetic resonance (MR) imaging. CPR images, created via a transformation from image-based to spine-based coordinate system, follow the structural shape of the spine and allow the whole course of the curved structure to be viewed in a single image. The spine-based coordinate system is defined on the 3D spine curve and on the axial vertebral rotation, both described by polynomial models. The 3D spine curve passes through the centers of vertebral bodies, and the axial vertebral rotation determines the rotation of vertebral spinous processes around the spine. The optimal polynomial parameters are found in an optimization framework, based on image analysis. The method was evaluated on 19 MR images of the spine from 10 patients.

In large multi-center studies it is important to quantify data variations due to differences between sites and over time. Here we investigate inter-site variability in signal to noise ratio (SNR), percent integral uniformity (PIU), width and height using the American College of Radiology (ACR) phantom scans from the NIHPD project. Longitudinal variations are also analyzed. All measurements are fully automated. Our results show that the mean SNR, PIU and the 2 metric values were statistically different across sites. The maximum mean difference in diameter across sites was 2 mm (1.1%), and the maximum mean difference in height was 2.5 mm (1.7%). Over time, an average drift of 0.4 mm per year was observed for the diameter while a drift of 0.5 mm per year was observed for the height. Trends observed over time often depended not only on site, but also on modality and scanner manufacturer.

Evaluating registration algorithms is difficult due to the lack of gold standard in most clinical procedures. The is a real-data based statistical method providing an alternative registration reference through a computationally intensive image database registration procedure. We propose in this paper an efficient implementation of this method through a grid-interfaced workflow enactor enabling the concurrent processing of hundreds of image registrations in a couple of hours only. The performances of two different grid infrastructures were compared. We computed the accuracy of 4 different rigid registration algorithms on longitudinal MRI images of brain tumors. Results showed an average subvoxel accuracy of 0.4 mm and 0.15 degrees in rotation.

Wireless endoscopy is a very recent and at the same time unique technique allowing to visualize and study the occurrence of contractions and to analyze the intestine motility. Feature extraction is essential for getting efficient patterns to detect contractions in wireless video endoscopy of small intestine. We propose a novel method based on anisotropic image filtering and efficient statistical classification of contraction features. In particular, we apply the image gradient tensor for mining informative skeletons from the original image and a sequence of descriptors for capturing the characteristic pattern of contractions. Features extracted from the endoluminal images were evaluated in terms of their discriminatory ability in correct classifying images as either belonging to contractions or not. Classification was performed by means of a support vector machine classifier with a radial basis function kernel. Our classification rates gave sensitivity of the order of 90.84% and specificity of the order of 94.43% respectively. These preliminary results highlight the high efficiency of the selected descriptors and support the feasibility of the proposed method in assisting the automatic detection and analysis of contractions.