Exact Geosedics and Shortest Paths on Polyhedral Surface

We present two algorithms for computing distances along a non-convex polyhedral surface. The first algorithm computes exact minimal-geodesic distances and the second algorithm combines these distances to compute exact shortest-path distances along the surface. Both algorithms have been extended to compute the exact minimalgeodesic paths and shortest paths. These algorithms have been implemented and validated on surfaces for which the correct solutions are known, in order to verify the accuracy and to measure the run-time performance, which is cubic or less for each algorithm. The exact-distance computations carried out by these algorithms are feasible for large-scale surfaces containing tens of thousands of vertices, and are a necessary component of near-isometric surface flattening methods that accurately transform curved manifolds into flat representations.National Institute for Biomedical Imaging and Bioengineering (R01 EB001550

Multi-Area Visuotopic Map Complexes in Macaque Striate and Extra-striate Cortex

We propose that a simple, closed-form mathematical expression--the Wedge-Dipole mapping--provides a concise approximation to the full-field, two-dimensional topographic structure of macaque V1, V2, and V3. A single map function, which we term a map complex, acts as a simultaneous descriptor of all three areas. Quantitative estimation of the Wedge-Dipole parameters is provided via 2DG data of central-field V1 topography and a publicly available data set of full-field macaque V1 and V2 topography. Good quantitative agreement is obtained between the data and the model presented here. The increasing importance of fMRI-based brain imaging motivates the development of more sophisticated two-dimensional models of cortical visuotopy, in contrast to the one-dimensional approximations that have been in common use. One reason is that topography has traditionally supplied an important aspect of "ground truth", or validation, for brain imaging, suggesting that further development of high-resolution fMRI will be facilitated by this data analysis. In addition, several important insights into the nature of cortical topography follows from this work. The presence of anisotropy in cortical magnification factor is shown to follow mathematically from the shared boundary conditions at the V1-V2 and V2-V3 borders, and therefore may not causally follow from the existence of columnar systems in these areas, as is widely assumed. An application of the Wedge-Dipole model to localizing aspects of visual processing to specific cortical areas--extending previous work in correlating V1 cortical magnification factor to retinal anatomy or visual psychophysics data--is briefly discussed.National Institute of Health/National Institute of Biomedical Imaging and Bioengineering (R01 EB001550

Reliable 3D mapping of ocular dominance columns in humans using GE-EPI fMRI at 7 T

Since the discovery of the BOLD effect, detection of ocular dominance columns (ODCs) in primary visual cortex (V1) served as a benchmark for high-precision functional magnetic resonance imaging (fMRI) (Menon et al., 1997; Dechent and Frahm 2000; Cheng et al., 2001; Yacoub et al., 2007). Although gradient-echo (GE) echo-planar imaging (EPI) is often used at lower field strengths, the applicability for high-resolution fMRI at higher field strengths is still under debate because of its inherent sensitivity to large draining veins (Polimeni et al., 2010). To counteract the loss of specificity, it was recently suggested to only sample far away from the pial surface when using GE-EPI (Nasr et al., 2016; Polimeni et al., 2017). Here, we assessed whether differential ocular dominance responses can be resolved using GE-EPI with different isotropic resolutions (0.8 mm and 1.0 mm) and how the corresponding BOLD signal is distributed across the cortex

Intersubject Regularity in the Intrinsic Shape of Human V1

Previous studies have reported considerable intersubject variability in the three-dimensional geometry of the human primary visual cortex (V1). Here we demonstrate that much of this variability is due to extrinsic geometric features of the cortical folds, and that the intrinsic shape of V1 is similar across individuals. V1 was imaged in ten ex vivo human hemispheres using high-resolution (200 μm) structural magnetic resonance imaging at high field strength (7 T). Manual tracings of the stria of Gennari were used to construct a surface representation, which was computationally flattened into the plane with minimal metric distortion. The instrinsic shape of V1 was determined from the boundary of the planar representation of the stria. An ellipse provided a simple parametric shape model that was a good approximation to the boundary of flattened V1. The aspect ration of the best-fitting ellipse was found to be consistent across subject, with a mean of 1.85 and standard deviation of 0.12. Optimal rigid alignment of size-normalized V1 produced greater overlap than that achieved by previous studies using different registration methods. A shape analysis of published macaque data indicated that the intrinsic shape of macaque V1 is also stereotyped, and similar to the human V1 shape. Previoud measurements of the functional boundary of V1 in human and macaque are in close agreement with these results

Sparsity-Promoting Calibration for GRAPPA Accelerated Parallel MRI Reconstruction

The amount of calibration data needed to produce images of adequate quality can prevent auto-calibrating parallel imaging reconstruction methods like generalized autocalibrating partially parallel acquisitions (GRAPPA) from achieving a high total acceleration factor. To improve the quality of calibration when the number of auto-calibration signal (ACS) lines is restricted, we propose a sparsity-promoting regularized calibration method that finds a GRAPPA kernel consistent with the ACS fit equations that yields jointly sparse reconstructed coil channel images. Several experiments evaluate the performance of the proposed method relative to unregularized and existing regularized calibration methods for both low-quality and underdetermined fits from the ACS lines. These experiments demonstrate that the proposed method, like other regularization methods, is capable of mitigating noise amplification, and in addition, the proposed method is particularly effective at minimizing coherent aliasing artifacts caused by poor kernel calibration in real data. Using the proposed method, we can increase the total achievable acceleration while reducing degradation of the reconstructed image better than existing regularized calibration methods.National Science Foundation (U.S.) (CAREER Grant 0643836)National Institutes of Health (U.S.) (Grant NIH R01 EB007942)National Institutes of Health (U.S.) (Grant NIH R01 EB006847)National Institutes of Health (U.S.) (Grant NIH P41 RR014075)National Institutes of Health (U.S.) (Grant NIH K01 EB011498)National Institutes of Health (U.S.) (Grant NIH F32 EB015914)National Science Foundation (U.S.). Graduate Research Fellowship Progra

Accelerated parallel magnetic resonance imaging reconstruction using joint estimation with a sparse signal model

Accelerating magnetic resonance imaging (MRI) by reducing the number of acquired k-space scan lines benefits conventional MRI significantly by decreasing the time subjects remain in the magnet. In this paper, we formulate a novel method for Joint estimation from Undersampled LinEs in Parallel MRI (JULEP) that simultaneously calibrates the GeneRalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) reconstruction kernel and reconstructs the full multi-channel k-space. We employ a joint sparsity signal model for the channel images in conjunction with observation models for both the acquired data and GRAPPA reconstructed k-space. We demonstrate using real MRI data that JULEP outperforms conventional GRAPPA reconstruction at high levels of undersampling, increasing the peak-signal-to-noise ratio by up to 10 dB.National Science Foundation (U.S.) (CAREER Grant 0643836)National Center for Research Resources (U.S.) (P41 RR014075)National Institutes of Health (U.S.) (NIH R01 EB007942)National Institutes of Health (U.S.) (NIH R01 EB006847)Siemens CorporationNational Science Foundation (U.S.). Graduate Research Fellowship Progra

Risk factors for obstructive sleep apnea syndrome in children: state of the art

The obstructive sleep apnea syndrome (OSAS) represents only part of a large group of pathologies of variable entity called respiratory sleep disorders (RSD) which include simple snoring and increased upper airway resistance syndrome (UARS). Although the etiopathogenesis of adult OSAS is well known, many aspects of this syndrome in children are still debated. Its prevalence is about 2% in children from 2 to 8 years of age, mostly related to the size of the upper airways adenoid tissue. Several risk factors linked to the development of OSAS are typical of the pediatric age. The object of this paper is to analyze the state of the art on this specific topic, discussing its implications in terms of diagnosis and management

High resolution quantitative and functional MRI indicate lower myelination of thin and thick stripes in human secondary visual cortex

The characterization of cortical myelination is essential for the study of structure-function relationships in the human brain. However, knowledge about cortical myelination is largely based on post mortem histology, which generally renders direct comparison to function impossible. The repeating pattern of pale-thin-pale-thick stripes of cytochrome oxidase (CO) activity in the primate secondary visual cortex (V2) is a prominent columnar system which is known to be differentiable by myelin content as well. However, depending on the applied histological method, higher myelination in both thin/thick and pale stripes were found, respectively. We used quantitative magnetic resonance imaging (qMRI) in conjunction with functional magnetic resonance imaging (fMRI) at ultra-high field strength (7T) to localize and study myelination of stripes in several humans at sub-millimeter resolution in vivo. Thin and thick stripes were functionally localized by exploiting their sensitivity to color and binocular disparity, respectively. Resulting functional activation maps showed robust stripe patterns in V2 which enabled further comparison of quantitative relaxation parameters between stripe types. Thereby, we found lower longitudinal relaxation rates (R1) of thin and thick stripes compared to surrounding gray matter in the order of 1-2%, indicating higher myelination of pale stripes. No differences for effective transverse relaxation rates (R2*) were found. The study demonstrates the feasibility to investigate structure-function relationships in living humans within one cortical area at the level of columnar systems using qMRI