{'fraction_non_alphanumeric': 0.03734777484707335, 'fraction_numerical': 0.024187664852124136, 'mean_word_length': 4.84341695359895, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mapping the Human Cortical Surface by Combining Quantitative T1 with Retinotopy†

We combined quantitative relaxation rate (R1= 1/T1) mapping—to measure local myelination—with fMRI-based retinotopy. Gray–white and pial surfaces were reconstructed and used to sample R1 at different cortical depths. Like myelination, R1 decreased from deeper to superficial layers. R1 decreased passing from V1 and MT, to immediately surrounding areas, then to the angular gyrus. High R1 was correlated across the cortex with convex local curvature so the data was first “de-curved”. By overlaying R1 and retinotopic maps, we found that many visual area borders were associated with significant R1 increases including V1, V3A, MT, V6, V6A, V8/VO1, FST, and VIP. Surprisingly, retinotopic MT occupied only the posterior portion of an oval-shaped lateral occipital R1 maximum. R1 maps were reproducible within individuals and comparable between subjects without intensity normalization, enabling multi-center studies of development, aging, and disease progression, and structure/function mapping in other modalities.


Introduction
Cortical areas are best defined and recognized on the basis of multiple converging techniques (Allman and Kaas 1971). Five measures developed during invasive experiments on animals include 1) receptotopic organization (e.g. retinotopy), 2) architectonic features, 3) connection patterns, 4) neurophysiological properties, and 5) effects of localized lesions. These measures have each been adapted to human brains. Studies on the effects of brain damage in humans have a long pedigree, recently supplemented by transcranial magnetic stimulation studies in normal subjects. Cortical-surface-based functional magnetic resonance imaging retinotopy (Engel et al. 1994;Sereno et al. 1995;DeYoe et al. 1996) is now established and often combined with functional studies, which constitute the overwhelming majority of in vivo human neuroimaging experiments. Connection patterns are beginning to be addressed in vivo with diffusionimaging-based fiber-tracking methods. But despite the historical precedence of postmortem human architectonic studies, in vivo human architectonics has perhaps been the least well-adapted of the 5 measures.
Postmortem cortical parcellation remains an difficult problem a century after the publication of Brodmann's classic cytoarchitectonic maps of the cortex in humans and other species. Related early work using cytoarchitectonic and myeloarchitectonic methods soon resulted in a plethora of human cortical area maps (reviews: Zilles and Amunts 2010;Geyer et al. 2011). But human cortical parcellation then went into decline, especially after the blistering critique of Lashley and Clark (1946). Bailey and von Bonin (1951), for example, divided the human neocortex into only 6 areas, even blurring the sharp border between V1 and V2 originally recognized by Meynert (1867).
Despite subsequent advances in mapping cortical areas in non-human primates (Merzenich and Kaas 1980), specifically human cortical parcellation languished (though see Braak (1980) on "pigmento"-architectonics). The growth of human neuroimaging in the 1980s had the effect of elevating Brodmann's map-after a rough translation from Brodmann's 2D summary images to a 3D atlas-to a world standard. Despite its acknowledged shortcomings (e.g. the middle temporal area (MT/V5) did not appear), the moderate resolution of early volume-averaged neuroimaging data did not initially demand better. Two decades later, as modern postmortem human cortical architectonics augmented by immunohistochemistry began to appear, the situation has improved considerably (e.g. Malikovic et al. 2007;Caspers et al. 2012).
However, a persistent stumbling block to combining structural and functional data has been to integrate postmortem parcellations based on high-resolution microscopic measures with retinotopic, connectional, and functional data from in vivo brains. Variability between postmortem and in vivo brains-whether at the single subject or group levelintroduces irreducible uncontrolled systematic variation.
Most previous attempts at in vivo architectonics have only been able to consider few localized regions of the cortex (Clark et al. 1992;Walters et al. 2003;Eickhoff et al. 2005;Sigalovsky et al. 2006;Geyer et al. 2011). One reason is signal-to-noise constraints due to in vivo scan time limits. But a more serious problem is the difficulty of reliably detecting small signal differences (1-4%) between different parts of non-quantitative T 1 -weighted (T 1 w) images that have large artifactual variations in brightness and contrast due to uncorrected inhomogeneities in the local radiofrequency transmit and receive fields. Widely applied post hoc histogram-based normalization methods (e.g. Dale et al. 1999) and T 1 w/T 2 weighted (T 2 w) ratio methods (Glasser and Van Essen 2011) remain sensitive to the effects of uncorrected B1 "transmit" inhomogeneities (flip-angle variation) on brightness and contrast, which affect both within-and between-individual comparisons. Ratio methods must also contend with vessel artifacts and local spatial distortion that differ between the 2 scan types, and unfavorable error propagation inherent in ratio estimation. Finally, differences in local cortical curvature (convex vs. concave) systematically affect the laminar and myeloarchitecture of the cortex, even within a cortical area ( Fig. 1; Smart and McSherry 1986;Annese et al. 2004).
Here, we present a detailed analysis that combines functional and structural data. To overcome the inherent limitations of non-quantitative structural imaging listed above, quantitative T 1 maps were acquired for all subjects, providing a measure of local myeloarchitecture (Schmierer et al. 2004;Draganski et al. 2011). The quantification allows for the direct comparison of datasets acquired on different subjects at different times and different MRI sites. The T 1 maps were then combined with retinotopic mapping across multiple sessions in the same group of subjects. High-resolution cortical-surface reconstructions of the gray-white matter and pial surfaces were made from the structural scans and used to analyze, visualize, and fuse the T 1 and retinotopy data. Since the longitudinal relaxation rate, R 1 (=1/T 1 ), positively correlates with the level of myelination and the brightness of T 1 w scans as commonly used for morphometry, we henceforth refer exclusively to it.

Materials and Methods
Six subjects (ages 22-55, 3 female) with normal or correctedto-normal vision participated in all parts of the study. Experimental protocols were approved by local ethics committees, and participants gave informed and signed written consent.

Structural Imaging
Structural images were acquired on a whole-body Tim Trio system (nominal field strength 3T, actual 2.89T, Siemens Healthcare, Erlangen, Germany), with body transmit and 32-channel receive head coil at the Wellcome Trust Centre for Neuroimaging. Proton densityweighted (PDw) and T 1 w images were acquired using an in-house multi-echo 3D FLASH pulse sequence ; voxel size: 0.8 × 0.8 × 0.81 mm 3 , field of view (FOV) = 256 × 216 × 194 mm, matrix = 320 × 270 × 240, repetition time (TR) = 23.7 ms, excitation flip angle: 6°(PDw) or 28°(T 1 w)). Acquisition was speeded up by 2× GRAPPA parallel imaging in the phase encoding and 6/8 Partial Fourier in the partition direction. To improve image quality, 4 gradient echoes were acquired (echo delay times (TE) = 2.2, 4.75, 7.3, 9.85 ms) after each excitation pulse. Each session consisted of four 10 min 31 s acquisitions (2 PDw and 2 T 1 w) and 2 shorter scans to estimate field inhomogeneities (see below). Quantitative R 1 maps were estimated from the PDw and T 1 w images according to the formalism developed by  including a correction for imperfect RF spoiling (Preibisch and Deichmann 2009). Recent applications of this method have demonstrated its robustness (Helms et al. , 2009Draganski et al. 2011). To correct for the effect of radio frequency (RF) transmit inhomogeneities on R 1 maps, maps of the transmit field B1 + were acquired using a 3D echo-planar imaging (EPI) spin-echo (SE)/stimulated echo (STE) method (Lutti et al. 2010(Lutti et al. , 2012FOV = 256 × 192 × 192 mm 3 , matrix = 64 × 48 × 48, TE SE /TE STE = 39.38/ 72.62 ms, TR = 500 ms, acquisition time 3 min 48 s), which was corrected for off-resonance effects using a B0 fieldmap. For further details, see Supplementary Information.

Functional Imaging and Analysis
Standard T 2 *-weighted EPI scans (3.2-mm isotropic resolution, 128 or 256 volumes each) and a T 1 w alignment scan with the same orientation and slice block center were acquired on a 1.5 whole-body Tim Avanto system (Siemens Healthcare) at the Birkbeck/UCL Centre for NeuroImaging, with body transmit and 32-channel receive head coil (see Supplementary Information for details). Five of 6 subjects completed at least 8 retinotopy scans (4 scans in 1 subject) and 4 ipsilateral mapping scans (2560 volumes per subject) across 3-5 additional sessions. See Supplementary Information for details of surface-based retinotopy analysis.
In-house OpenGL/Xlib software drove a video front-projection direct-view system that stimulated the visual field to an eccentricity of at least 57°of visual angle in all directions from central fixation. Balanced (clockwise/counterclockwise, in/out) phase-encoded polar Overall cortical average R 1 (equivalent T 1 on right y-axis) decreases monotonically from gray-white boundary (depth fraction 0.0) to a slight plateau at middle depths (0.3-0.6, beginning where second derivative in lower graph crosses zero), and then resumes its decrease in superficial layers (0.7-1.0; error bars: ±1 standard error of the mean over subjects). (B) Cross-ROI differences in average cortical R 1 shown as line for 8 depths (from 0.1 near white matter to 0.9 near pia); y-axes, error bars as in (A); ROIs: Angular (angular gyrus), angular-fs (FreeSurfer angular gyrus label), MT low and MT high [non-overlapping low-and high-probability MT labels (Malikovic et al. 2007;Fischl et al. 2008), V1-fs (FreeSurfer V1)]. All matched-paired t-tests on hypothesized differences significant (P < 0.05) except where marked "m" (P < 0.1), "=" (no significant difference), or "-" (difference opposite prediction). (C) Vertex-wise correlation (adjusted R 2 ) of R 1 and curvature as function of depth (error bars as before over subjects). Scatter plot inset at right is from a single subject at depth 0.5. For comparison, left inset shows myelin-stained section of human cortex with reduced myelination in concave areas (compare 2 white arrows, from Annese et al. 2004).
angle (wedge) and eccentricity (ring) stimuli contained medium-luminance-contrast colored checkerboards, optical flow fields, and simultaneous monitor-for-upside-face-among-right-side-up and monitor-for-number-among-letters tasks to maintain continuous peripheral attention (eight 64 second cycles). Ipsilateral field mapping stimuli used low contrast moving versus stationary gratings avoiding the center-of-gaze and the vertical meridian (see Supplementary Information for details).

Cortical-Surface Reconstruction and Sampling of R 1 Values Within Cortical Ribbon
Cortical surfaces were reconstructed with FreeSurfer (v5.0.0) from the aligned (AFNI 3dAllineate, hand-inspected) average of the 2 highresolution T 1 w scans. We initially used quantitative R 1 scans but experienced localized segmentation failures because some boundaries between the pial surface, the CSF, and the skull have different contrast than what is currently assumed by FreeSurfer algorithms (see Discussion).
R 1 data sets were sampled along the normal to each gray-white matter surface vertex in steps of 10% of cortical thickness (thickness estimated in FreeSurfer; Fischl and Dale 2000) and then smoothed tangentially at each depth with a 4-mm full width half maximum (FWHM) 2D kernel. The FreeSurfer estimate of local curvature (smoothed with a 1-mm FWHM 2D kernel) was used as a linear predictor of R 1 values at each depth. Vertex-wise residuals from this regression were used as "de-curved" estimates of R 1 values whose units are directly comparable to raw de-meaned R 1 values. See Supplementary Information for details of R 1 region of interest (ROI) choice.

R 1 as a Function of Cortical Depth
Postmortem studies of cortical myelination demonstrated a consistent pattern of myelination versus depth across most cortical areas. Myelination is highest immediately above the white matter, there are 2 bands of high myelination in deeper layers (inner/deeper and outer/shallower stria of Baillarger), and finally, there is a stepwise reduction in myelination in superficial layers. To investigate whether this overall pattern was apparent in our R 1 data, we sampled the R 1 values for each surface vertex at different fractional cortical depths, and then averaged these profiles across all vertices.
Average R 1 values ( Fig. 1A) show an extremely consistent decrease from deep to superficial layers (small to large depth fractions, z). A moderate plateau in deeper cortical layers begins around depth fraction 0.3, where the second spatial derivative with respect to z (∂ 2 R 1 /∂z 2 ), shown immediately below, crosses zero (an inflection point). The pattern of R 1 with depth closely resembles the overall profile of myelination seen in postmortem histology. It is considerably smoother, not distinguishing the 2 stria of Baillarger, which is a reflection of much coarser sampling units of in vivo human MRI (0.8 mm 3 here) compared with histology (∼0.01 mm 3 ), as well as the averaging of R 1 across the entire cortex of both hemispheres.

R 1 Variation over Probabilistically Defined Cortical ROIs
For an initial verification of regional differences in R 1 as a measure of myelination, we defined 5 ROIs with large expected myelination differences. The first was the FreeSurfer probabilistic V1 label (V1-fs). The second and third were MT high and MT low (non-overlapping high and low probability of cross-subject overlap regions). Finally, we included 2 angular gyrus labels-a smaller one, angular, defined on our subjects (non-visually responsive region just superior to MT) and a larger one, angular-fs, from the FreeSurfer parcellation.
Based on the histological literature, we predicted 1) that R 1 in V1 and the MT ROIs should exceed R 1 in the angular gyrus ROIs at all cortical depths, and that V1 should have the highest R 1 , 2) R 1 in high-probability MT should exceed R 1 in surrounding lower probability MT in deeper cortical layers (depth fractions 0.1-0.6). The systematic R 1 curves across areas at different depths illustrated in Figure 1B confirmed these predictions. All matched-paired t-tests were significant at P < 0.05 except where noted, with "m" indicating marginal (0.05 < P < 0.10) differences, "=" indicating no significant difference, and "-" indicating difference opposite to predicted direction. ROI-based analyses showed that at the level of macroscopic landmarks, our quantitative estimate of R 1 parallels known regional differences in myelination.

Relationship of R 1 to Local Cortical Curvature
The human cortex has deep sulci, but also a complex pattern of concavity and convexity, including many convex regions buried within major sulci. Postmortem studies in humans show that myeloarchitecture varies significantly with local cortical convexity (Annese et al. 2004); more convex regions are thicker and more myelinated, especially in middle and upper layers. We therefore examined the relationship between R 1 and local curvature as a function of cortical depth. Figure 1C shows that there is a moderately strong correlation between local curvature and R 1 at middle cortical depth fractions (up to a maximum average adjusted R 2 = 0.14) that falls off as one approaches the white matter and the pial surface. The bottom right inset scatterplot includes a least-squares fit line for R 1 against vertex-wise curvature; vertices in convex regions of the cortex (often but not exclusively on gyri) have higher R 1 than vertices in more concave regions. Since R 1 varies more at the gray-white matter border and at the pial border than it does within the cortex itself, small errors in tissue segmentation during surface reconstruction can have a large effect on laminar R 1 estimates near those boundaries. The fact that R 1 /curvature correlation is strongest at middle sampling depths suggests that segmentation errors are not responsible for our result.
Because local curvature is also related to local cortical thickness (average vertex-wise adjusted correlation between curvature and thickness: Adjusted R 2 = 0.081), local cortical thickness might explain some of the variation in R 1 . However, the strong correlation between convexity and R 1 survives the addition of vertex-wise thickness as a covariate. To minimize the effects of local curvature in our subsequent analyses of R 1 variation across the cortical sheet, we therefore initially regressed R 1 against curvature, and then used the "de-curved" residual values as our estimates of regional differences in cortical myelination.
Test-Retest Consistency of R 1 Across the Cortex The top 3 rows of Figure 2 show medial and lateral views of the cortical-surface-based cross-subject average ) of R 1 from 6 subjects (2 R 1 scans each), displayed on 1 of the 6 subjects' surfaces. The top row (gray-white matter surface), second row ( pial surface), and third through sixth rows (inflated surface) show same magnification views of one hemisphere. The heat scale overlay in all rows but row 3 shows "de-curved" R 1 values sampled from voxels lying at a cortical depth fraction of 0.5, while row 3 is not "de-curved" (but still de-meaned). The effects of "de-curving" can be seen by comparing rows 3 and 4 (e.g. the thin line posterior to MT disappears). Though regional variations in R 1 are small (<0.04 s −1 ), they are remarkably consistent over repeated measurements (bottom 2 rows). R 1 data were masked on the midline and in the anterior insula-in the first case because single hemisphere surfaces cut into the corpus callosum and thalamus, and in second case because of the difficulty of accurately reconstructing the gray-white matter surface superficial to the claustrum. Small regions of high cortical R 1 in posterior cingulate and anterior insula may thus have been omitted here.
These data show a strong resemblance to Flechsig's (Flechsig 1920) survey of perinatal infant myelogenesis, re-rendered in the top left inset of Figure 2. Flechsig found that early myelinating cortical regions (dark shading) tended to be the most densely stained for myelin in the adult cortex. These regions are clearly visible in our R 1 data and include primary motor (M-I) and somatosensory (S-I) areas, primary and secondary auditory areas, as well as a number of visual areas described in more detail in the next section.

Relation of R 1 to Retinotopic Maps
Visual areas were mapped in the same set of 6 subjects using a wide-field direct-view front-projection system. In 2-3 additional scan sessions for each subject, we mapped polar angle, eccentricity, and ipsilateral low-contrast visual field representations (Huk et al. 2002). We then performed cortical-surface-based cross-subject averages of complexvalued (amplitude and phase) data Hagler et al. 2007) and calculated visual field sign from the polar angle and eccentricity averages to identify borders in early areas (Sereno et al. 1995). A detailed pattern of correspondences between R 1 and retinotopic borders emerged when the 2 data sets were fused on the same surface.
The retinotopic and R 1 averages from both hemispheres in Figures 3 and 4 are illustrated in identical lateral, posterior, medial, and inferior poses (inflated surface) to aid comparison and are described together. Iso-R 1 borders (ΔR 1 = 0.020 s −1 ) around R 1 maxima were traced and superimposed on the polar angle data using thin white dashes to aid comparisons (see Supplementary Information for a movie that smoothly varies view and color scale contrast, more clearly illustrating the systematic fine level differences in the data shown in Fig. 4). Retinotopic borders (vertical meridian → circles, horizontal meridian → thick black dashes) are illustrated on both figures as is the posterior boundary of ipsilateral responses to low contrast moving gratings (thick yellow dashes in top row).
V1 (see medial views) was characterized by high R 1 (ΔR 1 ∼0.031 s −1 ), and its retinotopically defined borders in both hemispheres corresponded almost exactly with a sharp drop in R 1 . Superior to V1, V6 was bilaterally identified by a characteristic upper-to-lower field transition (moving superiorly and posteriorly) on the posterior bank of the parieto-occipital sulcus (Pitzalis et al. 2006) with R 1 values similar to V1 (ΔR 1 ∼0.035 s −1 ) with a shallow R 1 valley between it and V1 (see contrast ramp movie in Supplementary Information). Interestingly, Flechsig (1920) identified an early myelinating medial area in an almost identical location (Fig. 2, top inset, medial view). There was an additional shallower maximum of myelination (ΔR 1 ∼0.012 s −1 ) just anterior and superior to V6, containing mostly lower visual fields that we have tentatively labeled V6A (see Gamberini et al. 2011, for finer subdivisions).
In lateral occipital cortex, there was a prominent oval of high R 1 (ΔR 1 ∼0.029 s −1 ). We initially expected this might correspond to MT/V5. However, both polar angle retinotopy (reversal at an upper visual field vertical meridian) and the posterior bound of ipsilateral responses were situated near Figure 2. Cross-subject surface average relaxation rate (R 1 ) test/retest. Top left inset redrawn from Flechsig (1920) shows early myelinating areas in gray. Rows 1-6 show spherical morph cross-subject average ΔR 1 (difference above mean) projected back onto the gray-white, pial, and inflated surfaces of 1 subject. R 1 values in each subject were sampled at a point halfway between the individual's gray-white matter and pial surfaces. Row 3 shows the data before removing variance due to local cortical curvature; that variance is removed in all other rows. Rows 5 and 6 show excellent scan/rescan reliability of average. Lower middle inset shows same magnification lateral occipital cutout with individual subjects' R 1 maxima (yellow contours), superimposed on MT postmortem probability map in purple (Malikovic et al. 2007;Fischl et al. 2008), showing that the individual and average R 1 maxima differ subtly from the postmortem map. See text for midline and anterior insula mask.
the middle of the oval. This region was located within a sulcus ( Supplementary Fig. 2, showing individual subject data on each folded surface; the average was also resampled to each subject and shown as a green outline). We labeled the posterior part MT. To avoid name proliferation, we labeled the anterior and inferior part of this area, fundus of the superior temporal sulcus area (FST), on the basis of a similar retinotopic layout in non-human primates [e.g. macaque monkey, owl monkeys; see also Kolster et al. (2010) and Amano et al. (2009) who have subdivided this region into a larger number of areas with somewhat different names].
In dorsolateral occipital cortex near the midline, we identified a patch of high R 1 (ΔR 1 ∼0.028 s −1 ) that corresponded well with V3A (defined by the first appearance of upper field here moving anteriorly) in the left hemisphere and roughly so in the right. Inferiorly, beginning with the center-of-gaze representation of V1 (see inferior view), there is a large anteriorly extending region of high R 1 values (ΔR 1 ∼0.031 s −1 ) that also covers the adjoining center-of-gaze representations V2, V3/VP, and hV4 in both hemispheres. This maximum extends across parts of at least 4 different visual areas and may be artifactual (see Discussion). Still on the inferior surface, moving medially along the horizontal meridian of hV4 (Wandell et al. 2007), an anterior right-angle bend is visible in a region originally identified as V8 by Hadjikhani et al. (1998). This precisely corresponds to a small maximum of R 1 (ΔR 1 ∼0.021 s −1 ) in both hemispheres.
Back on the lateral surface, there is a prominent maximum of R 1 (ΔR 1 ∼0.020 s −1 ) in both hemispheres at the expected location of the ventral intraparietal area (VIP) . This is connected to the large high R 1 strip in somatosensory and motor cortex via a small isthmus. In our Figure 4. Cross-subject surface average quantitative R 1 (=1/T 1 ) maps. Spherical morph average maps of quantitative R 1 values sampled at 50% of cortical thickness are projected back to the same subject as in Figure 3, and posed and annotated identically. Quantitative R 1 is illustrated as variation from the mean (ΔR 1 ) after removing variance due to local cortical curvature. A histogram of ΔR 1 values from 3 right hemisphere visual areas (V6, MT, and VIP, defined by retinotopy) and 2 right hemisphere non-visual areas (inferior pre-cuneus "default mode" area and non-visual lateral occipital cortex) is shown in the lower middle inset. The non-visual regions of interest are shown as dashed lines in matching blue and green (upper left and middle left). The maxima shown are 3-4% higher than the average R 1 . polar angle mapping, this region contained lower visual field and horizontal meridian responses. The lower field bias may partly reflect head-centered remapping of receptive fields as a result of the lowered gaze position common with direct-view mapping.
Just posterior to VIP, best visible in the posterior view, there is an elongated region of moderately high R 1 values (ΔR 1 ∼0.013 s −1 ) in both hemispheres that extends posteriorly through the region originally identified as the human intraparietal area (LIP) (Sereno et al. 2001) eventually joining up with the prominent V3A maximum. Subsequently, this region has been subdivided into a number of LIP's (or IPS's) and contains several somewhat variable visual maps (e.g. Swisher et al. 2007;Silver and Kastner 2009). Finally, in frontal cortex, there are shallow R 1 maxima in the frontal eye fields, in an area identified as the polysensory zone (Graziano and Gandhi 2000;Sereno and Huang 2006;Huang and Sereno 2007), which appears as an extension off the motor strip, and there is a shallow R 1 maximum in the retinotopic part of the dorsolateral prefrontal cortex (Hagler and Sereno 2006).

Discussion
Cortical-surface-based analysis of quantitative longitudinal relaxation rate, R 1 (=1/T 1 ), showed reliable small (1-4%) regional differences in cortical gray matter likely due to differences in myelination, as verified by test-retest and by statistical analysis in selected regions defined in previous studies ( postmortem-defined V1 and MT, angular gyrus). We compared this in vivo architectonic measure with retinotopic mapping data obtained in separate sessions from the same subjects. This revealed that a number of boundaries visible in the R 1 maps correspond to previously recognized retinotopic borders, demonstrating the feasibility of using R 1 maps to help with cortical parcellation. Though not all adjoining areas have distinguishable R 1 values, a number can now be identified structurally without elaborate functional scans.
The cortex-wide average of R 1 as a function of laminar depth resembled postmortem myeloarchitectonic measures, systematically decreasing from the gray-white matter border to the pial surface with a moderate plateau in middle layers. These results are consistent with findings from localized high-resolution imaging and postmortem imaging (Walters et al. 2003;Eickhoff et al. 2005). The R 1 ranking of different ROIs-from most to least myelinated-is also remarkably consistent over sampling depths.
R 1 was significantly related to local curvature at intermediate depths-with convex parts of the cortex more strongly myelinated-quantitatively confirming observations in myelin-stained sections (Smart and McSherry 1986;Annese et al. 2004). Local curvature predicted R 1 better than a measure of whether a cortical region was on a sulcus or gyrus (the primary measure used for cross-subject surface alignment; Fischl et al. 1999). The strength and generality of this finding across the cortex suggests that there is a fundamental asymmetry between concave and convex bends in the cortex having to do with differences in the way deeper and more superficial layers can be deformed (Xu et al. 2010).
At a macroscopic level, the pattern of R 1 across the cortex showed a remarkable resemblance to postmortem studies of the order of myelogenesis (Flechsig 1920). By comparing these R 1 maps with retinotopic maps, we found that the borders of a number of visual areas defined in the same subjects-including V1, V3A, V6, V6A, V8, VIP, retinotopic dorsolateral prefrontal cortex, and the frontal eye fields-were marked by abrupt increases in R 1 (0.02-0.04 s −1 ). Perhaps our most surprising finding concerned the heavily myelinated oval in the lateral occipital cortex. Previously identified as MT/V5 in a number of studies, retinotopic mapping ( polar angle, posterior boundary of ipsilateral responses) showed that MT proper only accounted for the posterior 1/3 to 1/2 of that oval. This suggests that previous in vivo and postmortem studies in humans may have overestimated the extent of MT. Given that there are several areas with above average myelination anterior and inferior to MT in monkeys (MST and FST-e.g. Bock et al. 2009), this observation is not unprecedented; however, those other areas may be relatively larger or contain additional subdivisions in humans. Somatomotor and auditory cortex were also prominently recognizable (treated elsewhere).

Potential Applications
In contrast to T 1 w imaging, quantitative R 1 imaging directly estimates a basic physical property, the longitudinal relaxation rate of protons in a given tissue. Thus, the resulting maps are directly comparable without intensity normalization across individuals, scanners, and time-which especially lends itself to studying cortical areal differences, development, health versus disease, and disease progression. In this respect, quantitative neuroimaging is superior to standard postmortem myelin stains, which, despite their much higher resolution, are more subject to the uncontrollable vagaries of silver impregnation.
As an example, R 1 values in middle cortical layers for the high-probability MT label differ significantly from immediately adjacent regions, with an average R 1 value of ∼0.62 s −1 at a fractional cortical depth 0.5 (Fig. 1B), and an average decurved ΔR 1 of 0.023 s −1 . To investigate whether this highly myelinated region could be identified on a single subject basis, we thresholded each subject's surface-smoothed (4-mm FWHM) de-curved R 1 map at the average ΔR 1 value of 0.023 s −1 and searched for a disconnected supra-threshold R 1 patch within a lateral occipital MT+ search space based on previous human studies (Annese et al. 2005). Such patches were identified in all 6 subjects, and in 11 of 12 hemispheres ( Supplementary Fig. 1).

Methodological Issues and Prospects
Automated in vivo mapping of architectonics in single subjects requires detecting local changes in R 1 on the order of 1% at a high isotropic resolution of ∼800 um, which is made possible by advanced quantitative methods. In particular, we employed high-quality mapping of the B1+ transmit field (Lutti et al. 2010(Lutti et al. , 2012 and highly sensitive multi-echo 3D FLASH R 1 mapping to achieve maximal resolution, accuracy, and precision in the shortest time possible . We optimized the TR, flip angle and RF spoiling phase increment for highest accuracy and precision ) and applied corrections for imperfect spoiling of magnetization coherence pathways (Preibisch and Deichmann 2009). No post hoc brightness intensity normalization was used. The excellent scan-rescan reproducibility (Fig. 2) demonstrates that our R 1 estimates are robust, showing that quantitative in vivo maps can be obtained in only 25 min of scan time (see Supplementary Information for further discussion and description of R 1 mapping).
In recent work done in parallel with ours, cortical myelination has been estimated using a ratio between T 1 w and T 2 w scans (Glasser and Van Essen 2011). The results were broadly similar to ours but focused on large group results. Also, an additional series of post hoc normalization steps were required to visualize this non-quantitative data, which contain uncorrected artifacts in tissue contrast and brightness due to non-uniform B1 + transmit fields and to the effects of dividing pixel values from the 2 different image types that have contrasts sensitive to different aspects of microanatomy as well as different artifacts. The normalization steps are more difficult to reproduce at other centers than would be a hard threshold at cortical depth = 0.5 of R 1 = 0.62 s −1 .
Although local increases in cortical R 1 corresponded well with known patterns of myelination and with retinotopically defined areas, there were some unexpected local increases in R 1 -for example, at the nearly adjoined center-of-gaze regions of V1, V2, V3, and hV4. Because R 1 differences between cortical areas are comparable with R 1 differences between different cortical depths ( Fig. 1), detection of interareal differences requires uniformly accurate estimation of the gray-white and pial surfaces just as much as it requires R 1 uniformity. To improve those estimates, the more uniform R 1 images rather than T 1 w images could be used to reconstruct the pial surface. However, this will require extending the current FreeSurfer pial-surface-finding algorithm, which is finely adapted to T 1 w image intensity relations among tissue types. For example, in quantitative R 1 images, there is no change in image intensity between the superior margin of the cerebellum and inferior margin of the inferior occipital lobewhich leads to a ballooning of the pial surface therewhereas in T 1 w images, the cerebellum and inferior occipital lobe are separated by a thin but distinct hypointensity corresponding to the dura-likely a mixed contrast artifact, but one that FreeSurfer relies on. Another region where the Free-Surfer pial-surface reconstruction using R 1 images fails is near the thin temporal bone. It will be necessary to strategically incorporate quantitative proton density images into the pipeline to resolve this surface-reconstruction problem.
Another problem arises in areas such as S-I and V1 where parts of the cortical ribbon are so thin (∼1.5 mm) that they are spanned by less than two 0.8 mm 3 voxels. The problem of misestimating the cortical depth fraction and hence misestimating R 1 is particularly acute in this case, which might explain the unexpected reduction in R 1 in parafoveal V1 (see medial surface views in Fig. 3).
Quantitative R 1 mapping is fitted to answering a number of questions beyond those introduced here. For example, we are currently using it to characterize auditory and somatomotor fields (Dick et al. 2011) and are combining it with high-angular-resolution diffusion measurements in the cortex. This method could also illuminate structure/function mapping within cortical areas (e.g. in V2 stripes). Finally, the significant individual differences in local quantitative R 1 we have uncovered may have functional correlates.