Computerized Analysis of Magnetic Resonance Images to Study Cerebral Anatomy in Developing Neonates

The study of cerebral anatomy in developing neonates is of great importance for the understanding of brain development during the early period of life. This dissertation therefore focuses on three challenges in the modelling of cerebral anatomy in neonates during brain development. The methods that have been developed all use Magnetic Resonance Images (MRI) as source data. To facilitate study of vascular development in the neonatal period, a set of image analysis algorithms are developed to automatically extract and model cerebral vessel trees. The whole process consists of cerebral vessel tracking from automatically placed seed points, vessel tree generation, and vasculature registration and matching. These algorithms have been tested on clinical Time-of- Flight (TOF) MR angiographic datasets. To facilitate study of the neonatal cortex a complete cerebral cortex segmentation and reconstruction pipeline has been developed. Segmentation of the neonatal cortex is not effectively done by existing algorithms designed for the adult brain because the contrast between grey and white matter is reversed. This causes pixels containing tissue mixtures to be incorrectly labelled by conventional methods. The neonatal cortical segmentation method that has been developed is based on a novel expectation-maximization (EM) method with explicit correction for mislabelled partial volume voxels. Based on the resulting cortical segmentation, an implicit surface evolution technique is adopted for the reconstruction of the cortex in neonates. The performance of the method is investigated by performing a detailed landmark study. To facilitate study of cortical development, a cortical surface registration algorithm for aligning the cortical surface is developed. The method first inflates extracted cortical surfaces and then performs a non-rigid surface registration using free-form deformations (FFDs) to remove residual alignment. Validation experiments using data labelled by an expert observer demonstrate that the method can capture local changes and follow the growth of specific sulcus

Optimizing Magnetic Resonance Imaging for Image-Guided Radiotherapy

Magnetic resonance imaging (MRI) is playing an increasingly important role in image-guided radiotherapy. MRI provides excellent soft tissue contrast, and is flexible in characterizing various tissue properties including relaxation, diffusion and perfusion. This thesis aims at developing new image analysis and reconstruction algorithms to optimize MRI in support of treatment planning, target delineation and treatment response assessment for radiotherapy. First, unlike Computed Tomography (CT) images, MRI cannot provide electron density information necessary for radiation dose calculation. To address this, we developed a synthetic CT generation algorithm that generates pseudo CT images from MRI, based on tissue classification results on MRI for female pelvic patients. To improve tissue classification accuracy, we learnt a pelvic bone shape model from a training dataset, and integrated the shape model into an intensity-based fuzzy c-menas classification scheme. The shape-regularized tissue classification algorithm is capable of differentiating tissues that have significant overlap in MRI intensity distributions. Treatment planning dose calculations using synthetic CT image volumes generated from the tissue classification results show acceptably small variations as compared to CT volumes. As MRI artifacts, such as B1 filed inhomogeneity (bias field) may negatively impact the tissue classification accuracy, we also developed an algorithm that integrates the correction of bias field into the tissue classification scheme. We modified the fuzzy c-means classification by modeling the image intensity as the true intensity corrupted by the multiplicative bias field. A regularization term further ensures the smoothness of the bias field. We solved the optimization problem using a linearized alternating direction method of multipliers (ADMM) method, which is more computational efficient over existing methods. The second part of this thesis looks at a special MR imaging technique, diffusion-weighted MRI (DWI). By acquiring a series of DWI images with a wide range of b-values, high order diffusion analysis can be performed using the DWI image series and new biomarkers for tumor grading, delineation and treatment response evaluation may be extracted. However, DWI suffers from low signal-to-noise ratio at high b-values, and the multi-b-value acquisition makes the total scan time impractical for clinical use. In this thesis, we proposed an accelerated DWI scheme, that sparsely samples k-space and reconstructs images using a model-based algorithm. Specifically, we built a 3D block-Hankel tensor from k-space samples, and modeled both local and global correlations of the high dimensional k-space data as a low-rank property of the tensor. We also added a phase constraint to account for large phase variations across different b-values, and to allow reconstruction from partial Fourier acquisition, which further accelerates the image acquisition. We proposed an ADMM algorithm to solve the constrained image reconstruction problem. Image reconstructions using both simulated and patient data show improved signal-to-noise ratio. As compared to clinically used parallel imaging scheme which achieves a 4-fold acceleration, our method achieves an 8-fold acceleration. Reconstructed images show reduced reconstruction errors as proved on simulated data and similar diffusion parameter mapping results on patient data.PHDElectrical Engineering: SystemsUniversity of Michigan, Horace H. Rackham School of Graduate Studieshttps://deepblue.lib.umich.edu/bitstream/2027.42/143919/1/llliu_1.pd

Spatial fuzzy c-mean sobel algorithm with grey wolf optimizer for MRI brain image segmentation

Segmentation is the process of dividing the original image into multiple sub regions called segments in such a way that there is no intersection between any two regions. In medical images, the segmentation is hard to obtain due to the intensity similarity among various regions and the presence of noise in medical images. One of the most popular segmentation algorithms is Spatial Fuzzy C-means (SFCM). Although this algorithm has a good performance in medical images, it suffers from two issues. The first problem is lack of a proper strategy for point initialization step, which must be performed either randomly or manually by human. The second problem of SFCM is having inaccurate segmented edges. The goal of this research is to propose a robust medical image segmentation algorithm that overcomes these weaknesses of SFCM for segmenting magnetic resonance imaging (MRI) brain images with less human intervention. First, in order to find the optimum initial points, a histogram based algorithm in conjunction with Grey Wolf Optimizer (H-GWO) is proposed. The proposed H-GWO algorithm finds the approximate initial point values by the proposed histogram based method and then by taking advantage of GWO, which is a soft computing method, the optimum initial values are found. Second, in order to enhance SFCM segmentation process and achieve higher accurate segmented edges, an edge detection algorithm called Sobel was utilized. Therefore, the proposed hybrid SFCM-Sobel algorithm first finds the edges of the original image by Sobel edge detector algorithm and finally extends the edges of SFCM segmented images to the edges that are detected by Sobel. In order to have a robust segmentation algorithm with less human intervention, the H-GWO and SFCM-Sobel segmentation algorithms are integrated to have a semi-automatic robust segmentation algorithm. The results of the proposed H-GWO algorithms show that optimum initial points are achieved and the segmented images of the SFCM-Sobel algorithm have more accurate edges as compared to recent algorithms. Overall, quantitative analysis indicates that better segmentation accuracy is obtained. Therefore, this algorithm can be utilized to capture more accurate segmented in images in the era of medical imaging

Automatic segmentation of the human thigh muscles in magnetic resonance imaging

Advances in magnetic resonance imaging (MRI) and analysis techniques have improved diagnosis and patient treatment pathways. Typically, image analysis requires substantial technical and medical expertise and MR images can su↵er from artefacts, echo and intensity inhomogeneity due to gradient pulse eddy currents and inherent e↵ects of pulse radiation on MRI radio frequency (RF) coils that complicates the analysis. Processing and analysing serial sections of MRI scans to measure tissue volume is an additional challenge as the shapes and the borders between neighbouring tissues change significantly by anatomical location. Medical imaging solutions are needed to avoid laborious manual segmentation of specified regions of interest (ROI) and operator errors. The work set out in this thesis has addressed this challenge with a specific focus on skeletal muscle segmentation of the thigh. The aim was to develop an MRI segmentation framework for the quadriceps muscles, femur and bone marrow. Four contributions of this research include: (1) the development of a semi-automatic segmentation framework for a single transverse-plane image; (2) automatic segmentation of a single transverseplane image; (3) the automatic segmentation of multiple contiguous transverse-plane images from a full MRI thigh scan; and (4) the use of deep learning for MRI thigh quadriceps segmentation. Novel image processing, statistical analysis and machine learning algorithms were developed for all solutions and they were compared against current gold-standard manual segmentation. Frameworks (1) and (3) require minimal input from the user to delineate the muscle border. Overall, the frameworks in (1), (2) and (3) o↵er very good output performance, with respective framework’s mean segmentation accuracy by JSI and processing time of: (1) 0.95 and 17 sec; (2) 0.85 and 22 sec; and (3) 0.93 and 3 sec. For the framework in (4), the ImageNet trained model was customized by replacing the fully-connected layers in its architecture to convolutional layers (hence the name of Fully Convolutional Network (FCN)) and the pre-trained model was transferred for the ROI segmentation task. With the implementation of post-processing for image filtering and morphology to the segmented ROI, we have successfully accomplished a new benchmark for thigh MRI analysis. The mean accuracy and processing time with this framework are 0.9502 (by JSI ) and 0.117 sec per image, respectively

Foetal echocardiographic segmentation

Congenital heart disease affects just under one percentage of all live births [1]. Those defects that manifest themselves as changes to the cardiac chamber volumes are the motivation for the research presented in this thesis. Blood volume measurements in vivo require delineation of the cardiac chambers and manual tracing of foetal cardiac chambers is very time consuming and operator dependent. This thesis presents a multi region based level set snake deformable model applied in both 2D and 3D which can automatically adapt to some extent towards ultrasound noise such as attenuation, speckle and partial occlusion artefacts. The algorithm presented is named Mumford Shah Sarti Collision Detection (MSSCD). The level set methods presented in this thesis have an optional shape prior term for constraining the segmentation by a template registered to the image in the presence of shadowing and heavy noise. When applied to real data in the absence of the template the MSSCD algorithm is initialised from seed primitives placed at the centre of each cardiac chamber. The voxel statistics inside the chamber is determined before evolution. The MSSCD stops at open boundaries between two chambers as the two approaching level set fronts meet. This has significance when determining volumes for all cardiac compartments since cardiac indices assume that each chamber is treated in isolation. Comparison of the segmentation results from the implemented snakes including a previous level set method in the foetal cardiac literature show that in both 2D and 3D on both real and synthetic data, the MSSCD formulation is better suited to these types of data. All the algorithms tested in this thesis are within 2mm error to manually traced segmentation of the foetal cardiac datasets. This corresponds to less than 10% of the length of a foetal heart. In addition to comparison with manual tracings all the amorphous deformable model segmentations in this thesis are validated using a physical phantom. The volume estimation of the phantom by the MSSCD segmentation is to within 13% of the physically determined volume

Segmentation of pelvic structures from preoperative images for surgical planning and guidance

Prostate cancer is one of the most frequently diagnosed malignancies globally and the second leading cause of cancer-related mortality in males in the developed world. In recent decades, many techniques have been proposed for prostate cancer diagnosis and treatment. With the development of imaging technologies such as CT and MRI, image-guided procedures have become increasingly important as a means to improve clinical outcomes. Analysis of the preoperative images and construction of 3D models prior to treatment would help doctors to better localize and visualize the structures of interest, plan the procedure, diagnose disease and guide the surgery or therapy. This requires efficient and robust medical image analysis and segmentation technologies to be developed. The thesis mainly focuses on the development of segmentation techniques in pelvic MRI for image-guided robotic-assisted laparoscopic radical prostatectomy and external-beam radiation therapy. A fully automated multi-atlas framework is proposed for bony pelvis segmentation in MRI, using the guidance of MRI AE-SDM. With the guidance of the AE-SDM, a multi-atlas segmentation algorithm is used to delineate the bony pelvis in a new \ac{MRI} where there is no CT available. The proposed technique outperforms state-of-the-art algorithms for MRI bony pelvis segmentation. With the SDM of pelvis and its segmented surface, an accurate 3D pelvimetry system is designed and implemented to measure a comprehensive set of pelvic geometric parameters for the examination of the relationship between these parameters and the difficulty of robotic-assisted laparoscopic radical prostatectomy. This system can be used in both manual and automated manner with a user-friendly interface. A fully automated and robust multi-atlas based segmentation has also been developed to delineate the prostate in diagnostic MR scans, which have large variation in both intensity and shape of prostate. Two image analysis techniques are proposed, including patch-based label fusion with local appearance-specific atlases and multi-atlas propagation via a manifold graph on a database of both labeled and unlabeled images when limited labeled atlases are available. The proposed techniques can achieve more robust and accurate segmentation results than other multi-atlas based methods. The seminal vesicles are also an interesting structure for therapy planning, particularly for external-beam radiation therapy. As existing methods fail for the very onerous task of segmenting the seminal vesicles, a multi-atlas learning framework via random decision forests with graph cuts refinement has further been proposed to solve this difficult problem. Motivated by the performance of this technique, I further extend the multi-atlas learning to segment the prostate fully automatically using multispectral (T1 and T2-weighted) MR images via hybrid \ac{RF} classifiers and a multi-image graph cuts technique. The proposed method compares favorably to the previously proposed multi-atlas based prostate segmentation. The work in this thesis covers different techniques for pelvic image segmentation in MRI. These techniques have been continually developed and refined, and their application to different specific problems shows ever more promising results.Open Acces

White matter volume assessment in premature infants on MRI at term - computer aided volume analysis

The objective of this study is the development of an automatic segmentation framework for measuring volume changes in the white matter tissue from premature infant MRI data. The early stage of the brain development presents several major computational challenges such as structure and shape variations between patients. Furthermore, a high water content is present in the brain tissue, that leads to inconsistencies and overlapping intensity values across different brain structures. Another problem lies in low-frequency multiplicative intensity variations, which arises from an inhomogeneous magnetic field during the MRI acquisition. Finally, the segmentation is influenced by the partial volume effects which describe voxels that are generated by more than one tissue type. To overcome these challenges, this study is divided into three parts with the intention to locally segment the white matter tissue without the guidance of an atlas. Firstly, a novel brain extraction method is proposed with the aim to remove all non-brain tissue. The data quality can be improved by noise reduction using an anisotropic diffusion filter and intensity variations adjustments throughout the volume. In order to minimise the influence of missing contours and overlapping intensity values between brain and nonbrain tissue, a brain mask is created and applied during the extraction of the brain tissue. Secondly, the low-frequency intensity inhomogeneities are addressed by calculating the bias field which can be separated and corrected using low pass filtering. Finally, the segmentation process is performed by combining probabilistic clustering with classification algorithms. In order to achieve the final segmentation, the algorithm starts with a pre-segmentation procedure which was applied to reduce the intensity inhomogeneities within the white matter tissue. The key element in the segmentation process is the classification of diffused and missing contours as well as the partial volume voxels by performing a voxel reclassification scheme. The white matter segmentation framework was tested using the Dice Similarity Metric, and the numerical evaluation demonstrated precise segmentation results

A comparative evaluation for liver segmentation from spir images and a novel level set method using signed pressure force function

Thesis (Doctoral)--Izmir Institute of Technology, Electronics and Communication Engineering, Izmir, 2013Includes bibliographical references (leaves: 118-135)Text in English; Abstract: Turkish and Englishxv, 145 leavesDeveloping a robust method for liver segmentation from magnetic resonance images is a challenging task due to similar intensity values between adjacent organs, geometrically complex liver structure and injection of contrast media, which causes all tissues to have different gray level values. Several artifacts of pulsation and motion, and partial volume effects also increase difficulties for automatic liver segmentation from magnetic resonance images. In this thesis, we present an overview about liver segmentation methods in magnetic resonance images and show comparative results of seven different liver segmentation approaches chosen from deterministic (K-means based), probabilistic (Gaussian model based), supervised neural network (multilayer perceptron based) and deformable model based (level set) segmentation methods. The results of qualitative and quantitative analysis using sensitivity, specificity and accuracy metrics show that the multilayer perceptron based approach and a level set based approach which uses a distance regularization term and signed pressure force function are reasonable methods for liver segmentation from spectral pre-saturation inversion recovery images. However, the multilayer perceptron based segmentation method requires a higher computational cost. The distance regularization term based automatic level set method is very sensitive to chosen variance of Gaussian function. Our proposed level set based method that uses a novel signed pressure force function, which can control the direction and velocity of the evolving active contour, is faster and solves several problems of other applied methods such as sensitivity to initial contour or variance parameter of the Gaussian kernel in edge stopping functions without using any regularization term