Computation of the Mid-Sagittal Plane in 3D Brain Images

International audienceWe present a new method to automatically compute, reorient, and recenter the mid-sagittal plane in anatomical and functional three-dimensional (3-D) brain images. This iterative approach is composed of two steps. At first, given an initial guess of the mid-sagittal plane (generally, the central plane of the image grid), the computation of local similarity measures between the two sides of the head allows to identify homologous anatomical structures or functional areas, by way of a block matching procedure. The output is a set of point-to-point correspondences: the centers of homologous blocks. Subsequently, we define the mid-sagittal plane as the one best superposing the points on one side and their counterparts on the other side by reflective symmetry. Practically, the computation of the parameters characterizing the plane is performed by a least trimmed squares estimation. Then, the estimated plane is aligned with the center of the image grid, and the whole process is iterated until convergence. The robust estimation technique we use allows normal or abnormal asymmetrical structures or areas to be treated as outliers, and the plane to be mainly computed from the underlying gross symmetry of the brain. The algorithm is fast and accurate, even for strongly tilted heads, and even in presence of high acquisition noise and bias field, as shown on a large set of synthetic data. The algorithm has also been visually evaluated on a large set of real magnetic resonance (MR) images.We present a few results on isotropic as well as anisotropic anatomical (MR and computed tomography) and functional (single photon emission computed tomography and positron emission tomography) real images, for normal and pathological subjects

Computation of the Mid-Sagittal Plane in 3D Medical Images of the Head

We present a new symmetry-based method allowing to compute, reorient and recenter the mid-sagittal plane in anatomical and functional 3D images of the brain. In the literature, there are mainly two definitions of this plane: it is either the plane best fitting the inter-hemispheric fissure of the brain, or the plane best superposing the two sides of the head by reflective symmetry. We use this latter definition in our method, which is composed of two steps. At first, the computation of local similarity measures between the two sides of the brain allows to match homologous anatomical structures or functional areas, by way of a block matching procedure. The output is a set of point-to-point correspondences: the centers of homologous blocks. Subsequently, we define the mid-sagittal plane as the one superposing at best the points in one side of the head and their counterparts in the other side by reflective symmetry. The estimatio- n of the parameters characterizing the plane is performed by a least trimmed squares optimization scheme. Then, the estimated plane is aligned with the center of the image lattice. This method is fully automated, objective and reproducible. Our method tackles the main issue posed by the sym­me­try-ba- sed approach, that often relies on global similarity measures (such as the cross-correlation) between the intensities of the two flipped versions of the 3D image. The estimation of the mid-sagittal plane can be severely biased when normal or abnormal asymmetries hide the underlying symmetry of the brain or the skull. The computation of local measures of symmetry and the use of a robust estimation technique allow to discriminate between symmetrical and asymmetrical areas. Then, the mid-sagittal plane is mainly computed from the underlying gross symmetry of the brain, because its estimation is robust with respect to normal or abnormal asymmetries which are treated as outliers. We show on a large database of synthetic images that we can obtain a subvoxel accuracy in a CPU time of about 3 minutes, for strongly tilted heads, noisy and biased images. We present results on isotropic or anisotropic anatomical (MR, CT), and functional (SPECT and PET) images

Symmetry Identification Using Partial Surface Matching and Tilt Correction in 3D Brain Images

We propose a novel method to automatically compute the symmetry plane and correct the 3D orientation of patient brain images. Many images of the brain are clinically unreadable because of the misalignment of the patient's head in the scanner. We proposed an algorithm that represents the brain volume as a re-parameterized surface point cloud where each location has been parameterized by its elevation (latitude), azimuth (longitude) and radius. The removal of the interior contents of the brain makes this approach perform robustly in the presence of the brain pathologies, e.g. tumor, stroke and bleed. Thus, we decompose the symmetry plane computation problem into a surface matching routine. The search for the best matching surface is implemented in a multi-resolution paradigm so as to decrease computational time considerably. Spatial affine transform then is performed to rotate the 3D brain images and align them within the coordinate system of the scanner. The corrected brain volume is re-sliced such that each planar image represents the brain at the same axial level

Employing Symmetry Features for Automatic Misalignment Correction in Neuroimages

A novel method to automatically compute the symmetry plane and to correct the 3D orientation of neuro-images is presented. In acquisition of neuroimaging scans, the lack of perfect alignment of a patient's head makes it challenging to evaluate brain images. By deploying a shape-based criterion, the symmetry plane is defined as a plane that best matches external surface points on one side of the head, with their counterparts on the other side. In our method, the head volume is represented as a re-parameterized surface point cloud, where each location is parameterized by its elevation (latitude), azimuth (longitude), and radius. The search for the best matching surfaces is implemented in a multi-resolution paradigm, and the computation time is significantly decreased. The algorithm was quantitatively evaluated using in both simulated data and in real T1, T2, Flair magnetic resonance patient images. This algorithm is found to be fast (< 10s per MR volume), robust and accurate (< .6 degree of Mean Angular Error), invariant to the acquisition noise, slice thickness, bias field, and pathological asymmetries

Mid-sagittal plane and mid-sagittal surface optimization in brain MRI using a local symmetry measure

This paper describes methods for automatic localization of the mid-sagittal plane (MSP) and mid-sagittal sur-face (MSS). The data used is a subset of the Leukoaraiosis And DISability (LADIS) study consisting of three-dimensional magnetic resonance brain data from 62 elderly subjects (age 66 to 84 years). Traditionally, the mid-sagittal plane is localized by global measures. However, this approach fails when the partitioning plane between the brain hemispheres does not coincide with the symmetry plane of the head. We instead propose to use a sparse set of profiles in the plane normal direction and maximize the local symmetry around these using a general-purpose optimizer. The plane is parameterized by azimuth and elevation angles along with the distance to the origin in the normal direction. This approach leads to solutions confirmed as the optimal MSP in 98 percent of the subjects. Despite the name, the mid-sagittal plane is not always planar, but a curved surface resulting in poor partitioning of the brain hemispheres. To account for this, this paper also investigates an opti-mization strategy which fits a thin-plate spline surface to the brain data using a robust least median of squares estimator. Albeit computationally more expensive, mid-sagittal surface fitting demonstrated convincingly better partitioning of curved brains into cerebral hemispheres. 1

Automatic linear measurements of the fetal brain on MRI with deep neural networks

Timely, accurate and reliable assessment of fetal brain development is essential to reduce short and long-term risks to fetus and mother. Fetal MRI is increasingly used for fetal brain assessment. Three key biometric linear measurements important for fetal brain evaluation are Cerebral Biparietal Diameter (CBD), Bone Biparietal Diameter (BBD), and Trans-Cerebellum Diameter (TCD), obtained manually by expert radiologists on reference slices, which is time consuming and prone to human error. The aim of this study was to develop a fully automatic method computing the CBD, BBD and TCD measurements from fetal brain MRI. The input is fetal brain MRI volumes which may include the fetal body and the mother's abdomen. The outputs are the measurement values and reference slices on which the measurements were computed. The method, which follows the manual measurements principle, consists of five stages: 1) computation of a Region Of Interest that includes the fetal brain with an anisotropic 3D U-Net classifier; 2) reference slice selection with a Convolutional Neural Network; 3) slice-wise fetal brain structures segmentation with a multiclass U-Net classifier; 4) computation of the fetal brain midsagittal line and fetal brain orientation, and; 5) computation of the measurements. Experimental results on 214 volumes for CBD, BBD and TCD measurements yielded a mean $L_1$ difference of 1.55mm, 1.45mm and 1.23mm respectively, and a Bland-Altman 95% confidence interval ($CI_{95}$) of 3.92mm, 3.98mm and 2.25mm respectively. These results are similar to the manual inter-observer variability. The proposed automatic method for computing biometric linear measurements of the fetal brain from MR imaging achieves human level performance. It has the potential of being a useful method for the assessment of fetal brain biometry in normal and pathological cases, and of improving routine clinical practice.Comment: 15 pages, 8 figures, presented in CARS 2020, submitted to IJCAR

Glioma Diagnosis Aid through CNNs and Fuzzy-C Means for MRI

Glioma is a type of brain tumor that causes mortality in many cases. Early diagnosis is an important factor. Typically, it is detected through MRI and then either a treatment is applied, or it is removed through surgery. Deep-learning techniques are becoming popular in medical applications and image-based diagnosis. Convolutional Neural Networks are the preferred architecture for object detection and classification in images. In this paper, we present a study to evaluate the efficiency of using CNNs for diagnosis aids in glioma detection and the improvement of the method when using a clustering method (Fuzzy C-means) for preprocessing the input MRI dataset. Results offered an accuracy improvement from 0.77 to 0.81 when using Fuzzy C-Means.Ministerio de Economía y Competitividad TEC2016-77785-