A New Imaging Software To Evaluate The Centration Of Capsulorhexis And Intraocular Lens (IOL) After Cataract Surgery

Purpose: IOL centration and stability after cataract surgery is of high interest for cataract surgeons and IOL-producing companies. We present a new imaging software to evaluate the centration of the rhexis and the centration of the IOL after cataract surgery.Methods: We developed, in collaboration with the Biomedical Imaging Group (BIG), EPFL, Lausanne, a new working tool in order to assess precisely outcomes after IOL-implantation, such as ideal capsulorhexis and IOL-centration. The software is a plug-in of ImageJ, a general-purpose image processing and image-analysis package. The specifications of this software are: evaluation of the rhexis-centration and evaluation the position of the IOL in the posterior chamber. The end points are to analyze the quality of the centration of a rhexis after cataract surgery, the deformation of the rhexis with capsular bag retraction and the centration of the IOL after implantation.Results: This software delivers tools to interactively measure the distances between limbus, IOL and capsulorhexis and its changes over time. The user is invited to adjust nodes of three radial curves for the limbus, rhexis and the optic of the IOL. The radial distances of the curves are computed to evaluate the IOL implantation. The user is also able to define patterns for ideal capsulorhexis and optimal IOL-centration. We are going to present examples of calculations after cataract surgery.Conclusions: Evaluation of the centration of the rhexis and of the IOL after cataract surgery is an important end point for optimal IOL implantation after cataract surgery. Especially multifocal or accommodative lenses need a precise position in the bag with a good stability over time. This software is able to evaluate these parameters just after the surgery but also its changes over time. The results of these evaluations can lead to an optimizing of surgical procedures and materials.

Brain network analysis of patients with Temporal Lobe Epilepsy

Patients with Temporal Lobe Epilepsy (TLE) suffer from widespread subtle white matter abnormalities and abnormal functional connectivity extending beyond the affected lobe, as revealed by Diffusion Tensor MR Imaging, volumetric and functional MRI studies. Diffusion Spectrum Imaging (DSI) is a diffusion imaging technique with high angular resolution for improving the mapping of white matter pathways. In this study, we used DSI, connectivity matrices and topological measures to investigate how the alteration in structural connectivity influences whole brain structural networks. Eleven patients with right-sided TLE and hippocampal sclerosis and 18 controls underwent our DSI protocol at 3T. The cortical and subcortical grey matters were parcellated into 86 regions of interest and the connectivity between every region pair was estimated using global tractography and a connectivity matrix (the adjacency matrix of the structural network). We then compared the networks of patients and controls using topological measures. In patients, we found a higher characteristic path length and a lower clustering coefficient compared to controls. Local measures at node level of the clustering and efficiency showed a significant difference after a multiple comparison correction (Bonferroni). These significant nodes were located within as well outside the temporal lobe, and the localisation of most of them was consistent with regions known to be part of epileptic networks in TLE. Our results show altered connectivity patterns that are concordant with the mapping of functional epileptic networks in patients with TLE. Further studies are needed to establish the relevance of these findings for the propagation of epileptic activity, cognitive deficits in medial TLE and outcome of epilepsy surgery in individual patients.

Weighted Shape-based Averaging with Neighborhood Prior Model for Multiple Atlas Fusion-based Medical Image Segmentation

In medical imaging, merging automated segmentations obtained from multiple atlases has become a standard practice for improving the accuracy. In this letter, we propose two new fusion methods: "Global Weighted Shape-Based Averaging" (GWSBA) and "Local Weighted Shape-Based Averaging" (LWSBA). These methods extend the well known Shape-Based Averaging (SBA) by additionally incorporating the similarity information between the reference (i.e., atlas) images and the target image to be segmented. We also propose a new spatially-varying similarity-weighted neighborhood prior model, and an edge-preserving smoothness term that can be used with many of the existing fusion methods. We first present our new Markov Random Field (MRF) based fusion framework that models the above mentioned information. The proposed methods are evaluated in the context of segmentation of lymph nodes in the head and neck 3D CT images, and they resulted in more accurate segmentations compared to the existing SBA.

A convex optimization framework for global tractography

In this article we present a novel approach for diffusion MRI global tractography. Our formulation models the signal in each voxel as a linear combination of fiber-tract basis func- tions, which consist of a comprehensive set of plausible fiber tracts that are locally compatible with the measured MR signal. This large dictionary of candidate fibers is directly estimated from the data and, subsequently, efficient convex optimization techniques are used for recovering the smallest subset globally best fitting the measured signal. Experimen- tal results conducted on a realistic phantom demonstrate that our approach significantly reduces the computational cost of global tractography while still attaining a reconstruction quality at least as good as the state-of-the-art global methods.

An Efficient Segmentation Method for Ultrasound Images based on a Semi-supervised Approach and Patch-based Features

Segmenting ultrasound images is a challenging problemwhere standard unsupervised segmentation methods such asthe well-known Chan-Vese method fail. We propose in thispaper an efficient segmentation method for this class ofimages. Our proposed algorithm is based on asemi-supervised approach (user labels) and the use ofimage patches as data features. We also consider thePearson distance between patches, which has been shown tobe robust w.r.t speckle noise present in ultrasoundimages. Our results on phantom and clinical data show avery high similarity agreement with the ground truthprovided by a medical expert.

Ultra-High Field (7 Tesla) MRI for Gamma Knife Surgery Targeting of the Ventro-Intermediate Nucleus: A Pilot in Vivo Study

Aim: Gamma Knife surgery (GKS) is a non-invasive neurosurgical stereotactic procedure, increasingly used as an alternative to open functional procedures. This includes the targeting of the ventro-intermediate (Vim) nucleus of the thalamus for tremor. We currently perform an indirect targeting, using the "quadrilatere of Guyot," as the Vim nucleus is not visible on current 3 Tesla (T) MRI acquisitions. The primary objective of the current study was to enhance anatomic imaging for Vim GKS using high-field (7 T) MRI, with the aim of refining the visualization and precision of anatomical targeting. Method: Five young healthy subjects (mean age 23 years) were scanned both on 3 and 7 T MRI in Lausanne University Hospital (CHUV) and Center for Biomedical Imaging (CIBM). Classical T1-weighted MPRAGE, T2 CISS sequences (replacing former ventriculography) and diffusion tensor imaging were acquired at 3T. We obtained high-resolution susceptibility weighted images (SWI) at 7T for the visualization of thalamic subparts. SWI was further integrated for the first time into Leksell Gamma Plan® (LGP) software and co-registered with the 3T images. A simulation of targeting of the Vim was done using the "quadrilatere of Guyot" methodology on the 3T images. Furthermore, a correlation with the position of the found target on SWI was performed. The atlas of Morel et al. was used to confirm the findings on a detailed computer analysis outside LGP. Also, 3T and 7T MRI of one patient undergoing GKS Vim thalamotomy, were obtained before and 2 years after the procedure, and studied similarly. Results: The use of SWI provided a superior resolution and improved image contrast within the central gray matter. This allowed visualization and direct delineation of groups of thalamic nuclei in vivo, including the Vim. The position of the target, as assessed with the "quadrilatere of Guyot" method on 3 T, perfectly matched with the supposed one of the Vim on the SWI. Furthermore, a 3-dimensional model of the Vim target area was created on the basis of 3T and 7T images. Conclusion: This is the first report of the integration of SWI high-field MRI into the LGP in healthy subjects and in one patient treated GKS Vim thalamotomy. This approach aims at the improvement of targeting validation and further direct targeting of the Vim in tremor. The anatomical correlation between the direct visualization on 7T and the current targeting methods on 3T seems to show a very good anatomical matching.

How to measure cortical folding from MR images: a step-by-step tutorial to compute local gyrification index.

Cortical folding (gyrification) is determined during the first months of life, so that adverse events occurring during this period leave traces that will be identifiable at any age. As recently reviewed by Mangin and colleagues(2), several methods exist to quantify different characteristics of gyrification. For instance, sulcal morphometry can be used to measure shape descriptors such as the depth, length or indices of inter-hemispheric asymmetry(3). These geometrical properties have the advantage of being easy to interpret. However, sulcal morphometry tightly relies on the accurate identification of a given set of sulci and hence provides a fragmented description of gyrification. A more fine-grained quantification of gyrification can be achieved with curvature-based measurements, where smoothed absolute mean curvature is typically computed at thousands of points over the cortical surface(4). The curvature is however not straightforward to comprehend, as it remains unclear if there is any direct relationship between the curvedness and a biologically meaningful correlate such as cortical volume or surface. To address the diverse issues raised by the measurement of cortical folding, we previously developed an algorithm to quantify local gyrification with an exquisite spatial resolution and of simple interpretation. Our method is inspired of the Gyrification Index(5), a method originally used in comparative neuroanatomy to evaluate the cortical folding differences across species. In our implementation, which we name local Gyrification Index (lGI(1)), we measure the amount of cortex buried within the sulcal folds as compared with the amount of visible cortex in circular regions of interest. Given that the cortex grows primarily through radial expansion(6), our method was specifically designed to identify early defects of cortical development. In this article, we detail the computation of local Gyrification Index, which is now freely distributed as a part of the FreeSurfer Software (http://surfer.nmr.mgh.harvard.edu/, Martinos Center for Biomedical Imaging, Massachusetts General Hospital). FreeSurfer provides a set of automated reconstruction tools of the brain's cortical surface from structural MRI data. The cortical surface extracted in the native space of the images with sub-millimeter accuracy is then further used for the creation of an outer surface, which will serve as a basis for the lGI calculation. A circular region of interest is then delineated on the outer surface, and its corresponding region of interest on the cortical surface is identified using a matching algorithm as described in our validation study(1). This process is repeatedly iterated with largely overlapping regions of interest, resulting in cortical maps of gyrification for subsequent statistical comparisons (Fig. 1). Of note, another measurement of local gyrification with a similar inspiration was proposed by Toro and colleagues(7), where the folding index at each point is computed as the ratio of the cortical area contained in a sphere divided by the area of a disc with the same radius. The two implementations differ in that the one by Toro et al. is based on Euclidian distances and thus considers discontinuous patches of cortical area, whereas ours uses a strict geodesic algorithm and include only the continuous patch of cortical area opening at the brain surface in a circular region of interest.

MRI-based prostate brachytherapy seed localization.

A magnetic resonance imaging (MRI) pulse sequence and a corresponding image processing algorithm to localize prostate brachytherapy seeds during or after therapy are presented. Inversion-Recovery with ON-resonant water suppression (IRON) is an MRI methodology that generates positive contrast in regions of magnetic field susceptibility, as created by prostate brachytherapy seeds. Phantoms comprising of several materials found in brachytherapy seeds were created to assess the usability of the IRON pulse sequence for imaging seeds. Resulting images show that seed materials are clearly visible with high contrast using IRON, agreeing with theoretical predictions. A seed localization algorithm to process IRON images demonstrates the potential of this imaging technique for seed localization and dosimetry.

A realistic Computed Tomography simulator for small motion analysis of cerebral aneurysms.

This paper describes a realistic simulator for the Computed Tomography (CT) scan process for motion analysis. In fact, we are currently developing a new framework to find small motion from the CT scan. In order to prove the fidelity of this framework, or potentially any other algorithm, we present in this paper a simulator to simulate the whole CT acquisition process with a priori known parameters. In other words, it is a digital phantom for the motion analysis that can be used to compare the results of any related algorithm with the ground-truth realistic analytical model. Such a simulator can be used by the community to test different algorithms in the biomedical imaging domain. The most important features of this simulator are its different considerations to simulate the best the real acquisition process and its generality.

The role of embolic protection devices during carotid stenting prior to cardiac surgery in asymptomatic patients: Empty filters?

Objectives: The purpose of this study was to analyze the debris captured in the distal protection filters used during carotid artery stenting (CAS). Background: CAS is an option available to high-risk patients requiring revascularization. Filters are suggested for optimal stroke prevention during CAS. Methods: From May 2005 to June 2007, filters from 59 asymptomatic patients who underwent CAS were collected and sent to a specialized laboratory for light-microscope and histological analysis. Peri- and postprocedural outcomes were assessed during 1-year follow-up. Results: On the basis of biomedical imaging of the filter debris, the captured material could not be identified as embolized particles from the carotid plaque. On histological analysis the debris consisted mainly of red blood cell aggregates and/ or platelets, occasionally accompanied by granulocytes. We found no consistent histological evidence of embolized particles originating from atherosclerotic plaques. Post-procedure, three neurological events were reported: two (3.4%) transient ischemic attacks (TIA) and one (1.7%) ipsilateral minor stroke. Conclusion: The filters used during CAS in asymptomatic patients planned for cardiac surgery often remained empty. These findings may be explained by assuming that asymptomatic patients feature a different atherosclerotic plaque composition or stabilization through antiplatelet medication. Larger, randomized trials are clearly warranted, especially in the asymptomatic population. © 2012 Wiley Periodicals, Inc.

Vessel centerline tracking and boundary sgementation in coronary MRA with minimal manual interaction.

Magnetic resonance angiography (MRA) provides a noninvasive means to detect the presence, location and severity of atherosclerosis throughout the vascular system. In such studies, and especially those in the coronary arteries, the vessel luminal area is typically measured at multiple cross-sectional locations along the course of the artery. The advent of fast volumetric imaging techniques covering proximal to mid segments of coronary arteries necessitates automatic analysis tools requiring minimal manual interactions to robustly measure cross-sectional area along the three-dimensional track of the arteries in under-sampled and non-isotropic datasets. In this work, we present a modular approach based on level set methods to track the vessel centerline, segment the vessel boundaries, and measure transversal area using two user-selected endpoints in each coronary of interest. Arterial area and vessel length are measured using our method and compared to the standard Soap-Bubble reformatting and analysis tool in in-vivo non-contrast enhanced coronary MRA images.