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

High-resolution diffusion-weighted brain MRI under motion

Magnetic resonance imaging is one of the fastest developing medical imaging techniques. It provides excellent soft tissue contrast and has been a leading tool for neuroradiology and neuroscience research over the last decades. One of the possible MR imaging contrasts is the ability to visualize diffusion processes. The method, referred to as diffusion-weighted imaging, is one of the most common clinical contrasts but is prone to artifacts and is challenging to acquire at high resolutions. This thesis aimed to improve the resolution of diffusion weighted imaging, both in a clinical and in a research context. While diffusion-weighted imaging traditionally has been considered a 2D technique the manuscripts and methods presented here explore 3D diffusion acquisitions with isotropic resolution. Acquiring multiple small 3D volumes, or slabs, which are combined into one full volume has been the method of choice in this work. The first paper presented explores a parallel imaging driven multi-echo EPI readout to enable high resolution with reduced geometric distortions. The work performed on diffusion phase correction lead to an understanding that was used for the subsequent multi-slab papers. The second and third papers introduce the diffusion-weighted 3D multi-slab echo-planar imaging technique and explore its advantages and performance. As the method requires a slightly increased acquisition time the need for prospective motion correction became apparent. The forth paper suggests a new motion navigator using the subcutaneous fat surrounding the skull for rigid body head motion estimation, dubbed FatNav. The spatially sparse representation of the fat signal allowed for high parallel imaging acceleration factors, short acquisition times, and reduced geometric distortions of the navigator. The fifth manuscript presents a combination of the high-resolution 3D multi-slab technique and a modified FatNav module. Unlike our first FatNav implementation, using a single sagittal slab, this modified navigator acquired orthogonal projections of the head using the fat signal alone. The combined use of both presented methods provides a promising start for a fully motion corrected high-resolution diffusion acquisition in a clinical setting

Effect of Gradient Vectors Scheme and Noise Correction on Fractional Anisotropy in Diffusion Tensor Imaging of the Peripheral Nervous System

Diffusion Tensor Imaging (DTI) is a method widely used in research and clinic, especially for imaging and connectivity analysis of the white brain matter. Despite the many possibilities offered by DTI, this method suffers from an inherently low signal-to-noise ratio (SNR), since both the long echo time and the diffusion gradients weaken the signal. The SNR is particularly low at high spatial resolution, e.g. in the DTI of nerves. A low SNR leads to systematic and statistical errors in parameters calculated from the DTI, e.g. fractional anisotropy (FA). A low SNR can be partially compensated by increasing the number of diffusion directions or using methods for a posteriori noise correction. The most robust method for anatomical structures with unknown orientation is to distribute the diffusion gradients evenly in space. However, if the preferred direction of the anatomical structure is known in advance, it may be advantageous to limit the diffusion gradients to a cone centered on the axis of the structure. The aim of this work was to develop a DTI method with high accuracy and reliability for application in peripheral nerves. Two methods to reduce image noise were investigated: (1) A newly developed scheme of diffusion gradient vectors (DGV), where the vectors are restricted to a cone with an aperture angle Theta around the axis of the nerve and (2) different methods for a posteriori noise correction. For this purpose, Monte Carlo simulations were performed based on realistic values for diffusivity, FA and noise obtained from clinical investigations and studies. Furthermore, the methods were tested in a specially designed phantom simulating diffusion in peripheral nerves (FA = 0.65). These investigations were performed on a 3 Tesla whole-body magnetic resonance (MR) scanner. To determine the accuracy and reliability of the DTI using the appropriate measurement or correction procedures, systematic deviations of FA from baseline and the statistical error of FA were measured. The newly developed DGV scheme with limited space coverage was compared with gradient schemes with uniform space coverage (Jones, Downhill Simplex Method (DSM), gradient scheme of the manufacturer) based on their condition number (CN). The study showed that with the newly developed DGV scheme FA can be measured with high accuracy when the angle Theta is at least 45° or 60°. The minimum Theta depends on the number of gradient directions and on FA. Basically, the higher the FA value and the greater the number of gradients, the better the accuracy of the DGV scheme. For N = 30, the DGV allowed an exact determination of FA for the entire FA range (0.4 - 0.8) investigated in this study, if Theta ≥45° was. It could be shown that when using the new DGV scheme, a slight inclination of the investigated structure (≤30°) does not affect the accuracy of FA. CN of the developed DGV-scheme was higher than CN of the Jones-scheme and the DSM-scheme for N = 6; for N≥10 CN of the new DSM-scheme was lower than that of the Jones-scheme. However, it is also not to be expected that a method that concentrates the gradient vectors on a limited segment of space is as insensitive to interference as schemes with uniform gradient distribution. Nevertheless, the CN of the new DGV method was in the same order of magnitude as that of the other methods. A comparison of the different a posteriori correction methods showed that the power image method is the most effective and robust method and compensates for both the systematic and statistical errors of FA. The efficiency of the power image method is independent of the number of diffusion gradients used. In addition, the method works reliably - regardless of the method used for the coil combination (square sum versus adaptive combination). In contrast, both correction factor methods used in this study were less efficient in terms of noise correction; furthermore, the correction efficiency depended on the coil combination method. In conclusion, a combination of the newly developed DGV scheme with the power image method for a posteriori correction allows DTI of peripheral nerves with high SNR, high accuracy and reliability of the calculated parameters (e.g. FA) without the need for additional acquisition time. So far, however, these newly developed and tested methods have not yet been applied in studies or clinical trials

Generic acquisition protocol for quantitative MRI of the spinal cord

Quantitative spinal cord (SC) magnetic resonance imaging (MRI) presents many challenges, including a lack of standardized imaging protocols. Here we present a prospectively harmonized quantitative MRI protocol, which we refer to as the spine generic protocol, for users of 3T MRI systems from the three main manufacturers: GE, Philips and Siemens. The protocol provides guidance for assessing SC macrostructural and microstructural integrity: T1-weighted and T2-weighted imaging for SC cross-sectional area computation, multi-echo gradient echo for gray matter cross-sectional area, and magnetization transfer and diffusion weighted imaging for assessing white matter microstructure. In a companion paper from the same authors, the spine generic protocol was used to acquire data across 42 centers in 260 healthy subjects. The key details of the spine generic protocol are also available in an open-access document that can be found at https://github.com/spine-generic/protocols. The protocol will serve as a starting point for researchers and clinicians implementing new SC imaging initiatives so that, in the future, inclusion of the SC in neuroimaging protocols will be more common. The protocol could be implemented by any trained MR technician or by a researcher/clinician familiar with MRI acquisition

Towards Picogram Detection of Superparamagnetic Iron-Oxide Particles Using a Gradiometric Receive Coil

Superparamagnetic iron-oxide nanoparticles can be used in a variety of medical applications like vascular or targeted imaging. Magnetic particle imaging (MPI) is a promising tomographic imaging technique that allows visualizing the 3D nanoparticle distribution concentration in a non-invasive manner. The two main strengths of MPI are high temporal resolution and high sensitivity. While the first has been proven in the assessment of dynamic processes like cardiac imaging, it is unknown how far the detection limit of MPI can be lowered. Within this work, we will present a highly sensitive gradiometric receive-coil unit combined with a noise-matching network tailored for the measurement of mice. The setup is capable of detecting 5 ng of iron in vitro at 2.14 sec acquisition time. In terms of iron concentration we are able to detect 156 {\mu}g/L marking the lowest value that has been reported for an MPI scanner so far. In vivo MPI mouse images of a 512 ng bolus at 21.5 ms acquisition time allow for capturing the flow of an intravenously injected tracer through the heart of a mouse. Since it has been rather difficult to compare detection limits across MPI publications we propose guidelines improving the comparability of future MPI studies.Comment: 15 Pages, 7 Figures, V2: Changed the initials of Author Kannan M Krishnan, added two citations, corrected typo

Compressed sensing electron tomography of needle-shaped biological specimens--Potential for improved reconstruction fidelity with reduced dose.

Electron tomography is an invaluable method for 3D cellular imaging. The technique is, however, limited by the specimen geometry, with a loss of resolution due to a restricted tilt range, an increase in specimen thickness with tilt, and a resultant need for subjective and time-consuming manual segmentation. Here we show that 3D reconstructions of needle-shaped biological samples exhibit isotropic resolution, facilitating improved automated segmentation and feature detection. By using scanning transmission electron tomography, with small probe convergence angles, high spatial resolution is maintained over large depths of field and across the tilt range. Moreover, the application of compressed sensing methods to the needle data demonstrates how high fidelity reconstructions may be achieved with far fewer images (and thus greatly reduced dose) than needed by conventional methods. These findings open the door to high fidelity electron tomography over critically relevant length-scales, filling an important gap between existing 3D cellular imaging techniques.The research leading to these results has received funding from the European Union Seventh Framework Programme under Grant Agreement 312483 - ESTEEM2 (Integrated Infrastructure Initiative–I3), as well as from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC grant agreement 291522 - 3DIMAGE. B.W. and E.S. acknowledge financial support from the Deutsche Forschungsgemeinschaft (DFG) within the framework of the SPP 1570 as well as through the Cluster of Excellence “Engineering of Advanced Materials” at the Friedrich-Alexander-Universität ErlangenNürnberg. G.D. and C.D. acknowledge funding from the ERC under grant number 259619 PHOTO EM. B.W. acknowledges the Research Training Group “Disperse Systems for Electronic Applications” (DFG GEPRIS GRK 1161). R.L. acknowledges a Junior Research Fellowship from Clare College.This is the final version of the article. It first appeared from Elsevier via http://dx.doi.org/10.1016/j.ultramic.2015.10.02

Segmentation of images with low-contrast edges

A vast amount of the current research in medical image analysis has aimed to develop improved techniques of image segmentation. Of the existing approaches, active contour methods have proven effective by incorporating edge or region information from the image into a level set formulation. However, complications arise in images containing regions of low-contrast due to noise, occlusions, or partial volume effects, which are often unavoidable in practical applications. Incorporating prior shape information into the segmentation framework provides a more accurate and robust solution by constraining the evolving contour to resemble a target shape. Two methods are presented to incorporate a shape prior into existing active contour segmentation methods, including the edge-based geodesic active contours model and a fast update implementation of the region-based Chan-Vese model. Applying these methods to synthetic and real images demonstrates that an improved result can be obtained for images containing low-contrast edge regions

Tomographic measurement of all orthogonal components of three-dimensional displacement fields within scattering materials using wavelength scanning interferometry

Experimental mechanics is currently contemplating tremendous opportunities of further advancements thanks to a combination of powerful computational techniques and also fullfield non-contact methods to measure displacement and strain fields in a wide variety of materials. Identification techniques, aimed to evaluate material mechanical properties given known loads and measured displacement or strain fields, are bound to benefit from increased data availability (both in density and dimensionality) and efficient inversion methods such as finite element updating (FEU) and the virtual fields method (VFM). They work at their best when provided with dense and multicomponent experimental displacement (or strain) data, i.e. when all orthogonal components of displacements (or all components of the strain tensor) are known at points closely spaced within the volume of the material under study. Although a very challenging requirement, an increasing number of techniques are emerging to provide such data. In this Thesis, a novel wavelength scanning interferometry (WSI) system that provides three dimensional (3-D) displacement fields inside the volume of semi-transparent scattering materials is proposed. Sequences of two-dimensional interferograms are recorded whilst tuning the frequency of a laser at a constant rate. A new approach based on frequency multiplexing is used to encode the interference signal corresponding to multiple illumination directions at different spectral bands. Different optical paths along each illumination direction ensure that the signals corresponding to each sensitivity vector do not overlap in the frequency domain. All the information required to reconstruct the location and the 3-D displacement vector of scattering points within the material is thus recorded simultaneously in a single wavelength scan. By comparing phase data volumes obtained for two successive scans, all orthogonal components of the three dimensional displacement field introduced between scans (e.g. by means of loading or moving the sample under study) are readily obtained with high displacement sensitivity. The fundamental principle that describes the technique is presented in detail, including the correspondence between interference signal frequency and its associated depth within the sample, depth range, depth resolution, transverse resolution and displacement sensitivity. Data processing of the interference signal includes Fourier transformation, noise reduction, re-registration of data volumes, measurement of the illumination and sensitivity vectors from experimental data using a datum surface, phase difference evaluation, 3-D phase unwrapping and 3-D displacement field evaluation. Experiments consisting of controlled rigid body rotations and translations of a phantom were performed to validate the results. Both in-plane and the out-of-plane displacement components were measured for each voxel in the resulting data volume, showing an excellent agreement with the expected 3-D displacement

Study of Computational Image Matching Techniques: Improving Our View of Biomedical Image Data

Image matching techniques are proven to be necessary in various fields of science and engineering, with many new methods and applications introduced over the years. In this PhD thesis, several computational image matching methods are introduced and investigated for improving the analysis of various biomedical image data. These improvements include the use of matching techniques for enhancing visualization of cross-sectional imaging modalities such as Computed Tomography (CT) and Magnetic Resonance Imaging (MRI), denoising of retinal Optical Coherence Tomography (OCT), and high quality 3D reconstruction of surfaces from Scanning Electron Microscope (SEM) images. This work greatly improves the process of data interpretation of image data with far reaching consequences for basic sciences research. The thesis starts with a general notion of the problem of image matching followed by an overview of the topics covered in the thesis. This is followed by introduction and investigation of several applications of image matching/registration in biomdecial image processing: a) registration-based slice interpolation, b) fast mesh-based deformable image registration and c) use of simultaneous rigid registration and Robust Principal Component Analysis (RPCA) for speckle noise reduction of retinal OCT images. Moving towards a different notion of image matching/correspondence, the problem of view synthesis and 3D reconstruction, with a focus on 3D reconstruction of microscopic samples from 2D images captured by SEM, is considered next. Starting from sparse feature-based matching techniques, an extensive analysis is provided for using several well-known feature detector/descriptor techniques, namely ORB, BRIEF, SURF and SIFT, for the problem of multi-view 3D reconstruction. This chapter contains qualitative and quantitative comparisons in order to reveal the shortcomings of the sparse feature-based techniques. This is followed by introduction of a novel framework using sparse-dense matching/correspondence for high quality 3D reconstruction of SEM images. As will be shown, the proposed framework results in better reconstructions when compared with state-of-the-art sparse-feature based techniques. Even though the proposed framework produces satisfactory results, there is room for improvements. These improvements become more necessary when dealing with higher complexity microscopic samples imaged by SEM as well as in cases with large displacements between corresponding points in micrographs. Therefore, based on the proposed framework, a new approach is proposed for high quality 3D reconstruction of microscopic samples. While in case of having simpler microscopic samples the performance of the two proposed techniques are comparable, the new technique results in more truthful reconstruction of highly complex samples. The thesis is concluded with an overview of the thesis and also pointers regarding future directions of the research using both multi-view and photometric techniques for 3D reconstruction of SEM images

Change blindness: eradication of gestalt strategies

Arrays of eight, texture-defined rectangles were used as stimuli in a one-shot change blindness (CB) task where there was a 50% chance that one rectangle would change orientation between two successive presentations separated by an interval. CB was eliminated by cueing the target rectangle in the first stimulus, reduced by cueing in the interval and unaffected by cueing in the second presentation. This supports the idea that a representation was formed that persisted through the interval before being 'overwritten' by the second presentation (Landman et al, 2003 Vision Research 43149–164]. Another possibility is that participants used some kind of grouping or Gestalt strategy. To test this we changed the spatial position of the rectangles in the second presentation by shifting them along imaginary spokes (by ±1 degree) emanating from the central fixation point. There was no significant difference seen in performance between this and the standard task [F(1,4)=2.565, p=0.185]. This may suggest two things: (i) Gestalt grouping is not used as a strategy in these tasks, and (ii) it gives further weight to the argument that objects may be stored and retrieved from a pre-attentional store during this task