Designed especially for neurobiologists, FluoRender is an interactive tool for multi-channel fluorescence microscopy data visualization and analysis.
Deep brain stimulation
BrainStimulator is a set of networks that are used in SCIRun to perform simulations of brain stimulation such as transcranial direct current stimulation (tDCS) and magnetic transcranial stimulation (TMS).
Developing software tools for science has always been a central vision of the SCI Institute.

Image Analysis

SCI's imaging work addresses fundamental questions in 2D and 3D image processing, including filtering, segmentation, surface reconstruction, and shape analysis. In low-level image processing, this effort has produce new nonparametric methods for modeling image statistics, which have resulted in better algorithms for denoising and reconstruction. Work with particle systems has led to new methods for visualizing and analyzing 3D surfaces. Our work in image processing also includes applications of advanced computing to 3D images, which has resulted in new parallel algorithms and real-time implementations on graphics processing units (GPUs). Application areas include medical image analysis, biological image processing, defense, environmental monitoring, and oil and gas.


ross

Ross Whitaker

Segmentation
sarang

Sarang Joshi

Shape Statistics
Segmentation
Brain Atlasing
tolga

Tolga Tasdizen

Image Processing
Machine Learning
chris

Chris Johnson

Diffusion Tensor Analysis
shireen

Shireen Elhabian

Image Analysis
Computer Vision


Funded Research Projects:



Publications in Image Analysis:


Informative features of local field potential signals in primary visual cortex during natural image stimulation
S.M. Seyedhosseini, S. Shushruth, T. Davis, J.M. Ichida, P.A. House, B. Greger, A. Angelucci, T. Tasdizen. In Journal of Neurophysiology, Vol. 113, No. 5, American Physiological Society, pp. 1520--1532. March, 2015.
DOI: 10.1152/jn.00278.2014

The local field potential (LFP) is of growing importance in neurophysiology as a metric of network activity and as a readout signal for use in brain-machine interfaces. However, there are uncertainties regarding the kind and visual field extent of information carried by LFP signals, as well as the specific features of the LFP signal conveying such information, especially under naturalistic conditions. To address these questions, we recorded LFP responses to natural images in V1 of awake and anesthetized macaques using Utah multielectrode arrays. First, we have shown that it is possible to identify presented natural images from the LFP responses they evoke using trained Gabor wavelet (GW) models. Because GW models were devised to explain the spiking responses of V1 cells, this finding suggests that local spiking activity and LFPs (thought to reflect primarily local synaptic activity) carry similar visual information. Second, models trained on scalar metrics, such as the evoked LFP response range, provide robust image identification, supporting the informative nature of even simple LFP features. Third, image identification is robust only for the first 300 ms following image presentation, and image information is not restricted to any of the spectral bands. This suggests that the short-latency broadband LFP response carries most information during natural scene viewing. Finally, best image identification was achieved by GW models incorporating information at the scale of ∼0.5° in size and trained using four different orientations. This suggests that during natural image viewing, LFPs carry stimulus-specific information at spatial scales corresponding to few orientation columns in macaque V1.



Nonlinear Regression with Logistic Product Basis Networks
M. Sajjadi, M. Seyedhosseini,, T. Tasdizen. In IEEE Signal Processing Letters, Vol. 22, No. 8, IEEE, pp. 1011--1015. Aug, 2015.
DOI: 10.1109/lsp.2014.2380791

We introduce a novel general regression model that is based on a linear combination of a new set of non-local basis functions that forms an effective feature space. We propose a training algorithm that learns all the model parameters simultaneously and offer an initialization scheme for parameters of the basis functions. We show through several experiments that the proposed method offers better coverage for high-dimensional space compared to local Gaussian basis functions and provides competitive performance in comparison to other state-of-the-art regression methods.



Disjunctive normal random forests
M. Seyedhosseini , T. Tasdizen. In Pattern Recognition, Vol. 48, No. 3, Elsevier BV, pp. 976--983. March, 2015.
DOI: 10.1016/j.patcog.2014.08.023

We develop a novel supervised learning/classification method, called disjunctive normal random forest (DNRF). A DNRF is an ensemble of randomly trained disjunctive normal decision trees (DNDT). To construct a DNDT, we formulate each decision tree in the random forest as a disjunction of rules, which are conjunctions of Boolean functions. We then approximate this disjunction of conjunctions with a differentiable function and approach the learning process as a risk minimization problem that incorporates the classification error into a single global objective function. The minimization problem is solved using gradient descent. DNRFs are able to learn complex decision boundaries and achieve low generalization error. We present experimental results demonstrating the improved performance of DNDTs and DNRFs over conventional decision trees and random forests. We also show the superior performance of DNRFs over state-of-the-art classification methods on benchmark datasets.



Splenium development and early spoken language in human infants
M. R. Swanson, J. J. Wolff, J. T. Elison, H. Gu, H. C. Hazlett, K. Botteron, M. Styner, S. Paterson, G. Gerig, J. Constantino, S. Dager, A. Estes, C. Vachet, J. Piven. In Developmental Science, Wiley Online Library, 2015.
ISSN: 1467-7687
DOI: 10.1111/desc.12360

The association between developmental trajectories of language-related white matter fiber pathways from 6 to 24 months of age and individual differences in language production at 24 months of age was investigated. The splenium of the corpus callosum, a fiber pathway projecting through the posterior hub of the default mode network to occipital visual areas, was examined as well as pathways implicated in language function in the mature brain, including the arcuate fasciculi, uncinate fasciculi, and inferior longitudinal fasciculi. The hypothesis that the development of neural circuitry supporting domain-general orienting skills would relate to later language performance was tested in a large sample of typically developing infants. The present study included 77 infants with diffusion weighted MRI scans at 6, 12 and 24 months and language assessment at 24 months. The rate of change in splenium development varied significantly as a function of language production, such that children with greater change in fractional anisotropy (FA) from 6 to 24 months produced more words at 24 months. Contrary to findings from older children and adults, significant associations between language production and FA in the arcuate, uncinate, or left inferior longitudinal fasciculi were not observed. The current study highlights the importance of tracing brain development trajectories from infancy to fully elucidate emerging brain–behavior associations while also emphasizing the role of the splenium as a key node in the structural network that supports the acquisition of spoken language.



Evaluating Alignment of Shapes by Ensemble Visualization
M. Raj, M. Mirzargar, R. Kirby, R. Whitaker, J. Preston. In IEEE Computer Graphics and Applications, IEEE, 2015.

The visualization of variability in 3D shapes or surfaces, which is a type of ensemble uncertainty visualization for volume data, provides a means of understanding the underlying distribution for a collection or ensemble of surfaces. While ensemble visualization for surfaces is already described in the literature, we conduct an expert-based evaluation in a particular medical imaging application: the construction of atlases or templates from a population of images. In this work, we extend contour boxplots to 3D, allowing us to evaluate it against an enumeration-style visualization of the ensemble members and also other conventional visualizations used by atlas builders, namely examining the atlas image and the corresponding images/data provided as part of the construction process. We present feedback from domain experts on the efficacy of contour boxplots compared to other modalities when used as part of the atlas construction and analysis stages of their work.



Entropy-based particle correspondence for shape populations,
I. OguzI, J. Cates, M. Datar, B. Paniagua, T. Fletcher, C. Vachet, M. Styner, R. Whitaker. In International Journal of Computer Assisted Radiology and Surgery, Springer, pp. 1-12. December, 2015.

Purpose
Statistical shape analysis of anatomical structures plays an important role in many medical image analysis applications such as understanding the structural changes in anatomy in various stages of growth or disease. Establishing accurate correspondence across object populations is essential for such statistical shape analysis studies.

Methods
In this paper, we present an entropy-based correspondence framework for computing point-based correspondence among populations of surfaces in a groupwise manner. This robust framework is parameterization-free and computationally efficient. We review the core principles of this method as well as various extensions to deal effectively with surfaces of complex geometry and application-driven correspondence metrics.

Results
We apply our method to synthetic and biological datasets to illustrate the concepts proposed and compare the performance of our framework to existing techniques.

Conclusions
Through the numerous extensions and variations presented here, we create a very flexible framework that can effectively handle objects of various topologies, multi-object complexes, open surfaces, and objects of complex geometry such as high-curvature regions or extremely thin features.



A Mixture Model for Automatic Diffeomorphic Multi-Atlas Building
M. Zhang, H. Shao, P. T. Fletcher. In MICCAI Workshop, Springer, 2015.

Computing image atlases that are representative of a dataset
is an important first step for statistical analysis of images. Most current approaches estimate a single atlas to represent the average of a large population of images, however, a single atlas is not sufficiently expressive to capture distributions of images with multiple modes. In this paper, we present a mixture model for building diffeomorphic multi-atlases that can represent sub-populations without knowing the category of each observed data point. In our probabilistic model, we treat diffeomorphic image transformations as latent variables, and integrate them out using a Monte Carlo Expectation Maximization (MCEM) algorithm via Hamiltonian Monte Carlo (HMC) sampling. A key benefit of our model is that the mixture modeling inference procedure results in an automatic clustering of the dataset. Using 2D synthetic data generated from known parameters, we demonstrate the ability of our model to successfully recover the multi-atlas and automatically cluster the dataset. We also show the effectiveness of the proposed method in a multi-atlas estimation problem for 3D brain images.



The DTI Challenge: Toward Standardized Evaluation of Diffusion Tensor Imaging Tractography for Neurosurgery
S. Pujol, W. Wells, C. Pierpaoli, C. Brun, J. Gee, G. Cheng, B. Vemuri, O. Commowick, S. Prima, A. Stamm, M. Goubran, A. Khan, T. Peters, P. Neher, K. H. Maier-Hein, Y. Shi, A. Tristan-Vega, G. Veni, R. Whitaker, M. Styner, C.F. Westin, S. Gouttard, I. Norton, L. Chauvin, H. Mamata, G. Gerig, A. Nabavi, A. Golby,, R. Kikinis. In Journal of Neuroimaging, Wiley, August, 2015.
DOI: 10.1111/jon.12283

BACKGROUND AND PURPOSE
Diffusion tensor imaging (DTI) tractography reconstruction of white matter pathways can help guide brain tumor resection. However, DTI tracts are complex mathematical objects and the validity of tractography-derived information in clinical settings has yet to be fully established. To address this issue, we initiated the DTI Challenge, an international working group of clinicians and scientists whose goal was to provide standardized evaluation of tractography methods for neurosurgery. The purpose of this empirical study was to evaluate different tractography techniques in the first DTI Challenge workshop.

METHODS
Eight international teams from leading institutions reconstructed the pyramidal tract in four neurosurgical cases presenting with a glioma near the motor cortex. Tractography methods included deterministic, probabilistic, filtered, and global approaches. Standardized evaluation of the tracts consisted in the qualitative review of the pyramidal pathways by a panel of neurosurgeons and DTI experts and the quantitative evaluation of the degree of agreement among methods.

RESULTS
The evaluation of tractography reconstructions showed a great interalgorithm variability. Although most methods found projections of the pyramidal tract from the medial portion of the motor strip, only a few algorithms could trace the lateral projections from the hand, face, and tongue area. In addition, the structure of disagreement among methods was similar across hemispheres despite the anatomical distortions caused by pathological tissues.

CONCLUSIONS
The DTI Challenge provides a benchmark for the standardized evaluation of tractography methods on neurosurgical data. This study suggests that there are still limitations to the clinical use of tractography for neurosurgical decision making.



Performance of an Efficient Image-registration Algorithm in Processing MR Renography Data,
C.C. Conlin, J.L. Zhang, F. Rousset, C. Vachet, Y. Zhao, K.A. Morton, K. Carlston, G. Gerig, V.S. Lee. In J Magnetic Resonance Imaging, July, 2015.
DOI: 10.1002/jmri.25000

PURPOSE:
To evaluate the performance of an edge-based registration technique in correcting for respiratory motion artifacts in magnetic resonance renographic (MRR) data and to examine the efficiency of a semiautomatic software package in processing renographic data from a cohort of clinical patients.

MATERIALS AND METHODS:
The developed software incorporates an image-registration algorithm based on the generalized Hough transform of edge maps. It was used to estimate glomerular filtration rate (GFR), renal plasma flow (RPF), and mean transit time (MTT) from 36 patients who underwent free-breathing MRR at 3T using saturation-recovery turbo-FLASH. The processing time required for each patient was recorded. Renal parameter estimates and model-fitting residues from the software were compared to those from a previously reported technique. Interreader variability in the software was quantified by the standard deviation of parameter estimates among three readers. GFR estimates from our software were also compared to a reference standard from nuclear medicine.

RESULTS:
The time taken to process one patient's data with the software averaged 12 ± 4 minutes. The applied image registration effectively reduced motion artifacts in dynamic images by providing renal tracer-retention curves with significantly smaller fitting residues (P < 0.01) than unregistered data or data registered by the previously reported technique. Interreader variability was less than 10% for all parameters. GFR estimates from the proposed method showed greater concordance with reference values (P < 0.05).

CONCLUSION:
These results suggest that the proposed software can process MRR data efficiently and accurately. Its incorporated registration technique based on the generalized Hough transform effectively reduces respiratory motion artifacts in free-breathing renographic acquisitions. J. Magn. Reson. Imaging 2015.



A Kalman Filtering Perspective for Multiatlas Segmentation
Y. Gao, L. Zhu, J. Cates, R. S. MacLeod, S. Bouix,, A. Tannenbaum. In SIAM J. Imaging Sciences, Vol. 8, No. 2, pp. 1007-1029. 2015.
DOI: 10.1137/130933423

In multiatlas segmentation, one typically registers several atlases to the novel image, and their respective segmented label images are transformed and fused to form the final segmentation. In this work, we provide a new dynamical system perspective for multiatlas segmentation, inspired by the following fact: The transformation that aligns the current atlas to the novel image can be not only computed by direct registration but also inferred from the transformation that aligns the previous atlas to the image together with the transformation between the two atlases. This process is similar to the global positioning system on a vehicle, which gets position by inquiring from the satellite and by employing the previous location and velocity—neither answer in isolation being perfect. To solve this problem, a dynamical system scheme is crucial to combine the two pieces of information; for example, a Kalman filtering scheme is used. Accordingly, in this work, a Kalman multiatlas segmentation is proposed to stabilize the global/affine registration step. The contributions of this work are twofold. First, it provides a new dynamical systematic perspective for standard independent multiatlas registrations, and it is solved by Kalman filtering. Second, with very little extra computation, it can be combined with most existing multiatlas segmentation schemes for better registration/segmentation accuracy.



Modeling Brain Growth and Development
N. Sadeghi, J. H. Gilmore , G. Gerig. In Brain, Vol. 1, pp. 429-436. 2015.
DOI: 10.1016/B978-0-12-397025-1.00314-6

Early brain development is characterized by rapid organization and structuring. Magnetic resonance–diffusion tensor imaging (MR-DTI) provides the possibility of capturing these changes noninvasively by following individuals longitudinally to better understand departures from normal brain development in subjects at risk for mental illness. This article illustrates the modeling of neurodevelopmental trajectories using a recently developed framework. Descriptions include the estimation of normative models for healthy singletons and twins and a statistical framework to predict development at 2 years of age only based on neonatal image data – a capability with excellent potential for preclinical diagnosis and eventual early therapeutic intervention.



Altered corpus callosum morphology associated with autism over the first 2 years of life
J. J. Wolff, G. Gerig, J. D. Lewis, T. Soda, M. A. Styner, C. Vachet, K. N. Botteron, J. T. Elison, S. R. Dager, A. M. Estes, H. C. Hazlett, R. T. Schultz, L. Zwaigenbaum, J. Piven. In Brain, 2015.
DOI: 10.1093/brain/awv118

Numerous brain imaging studies indicate that the corpus callosum is smaller in older children and adults with autism spectrum disorder. However, there are no published studies examining the morphological development of this connective pathway in infants at-risk for the disorder. Magnetic resonance imaging data were collected from 270 infants at high familial risk for autism spectrum disorder and 108 low-risk controls at 6, 12 and 24 months of age, with 83% of infants contributing two or more data points. Fifty-seven children met criteria for ASD based on clinical-best estimate diagnosis at age 2 years. Corpora callosa were measured for area, length and thickness by automated segmentation. We found significantly increased corpus callosum area and thickness in children with autism spectrum disorder starting at 6 months of age. These differences were particularly robust in the anterior corpus callosum at the 6 and 12 month time points. Regression analysis indicated that radial diffusivity in this region, measured by diffusion tensor imaging, inversely predicted thickness. Measures of area and thickness in the first year of life were correlated with repetitive behaviours at age 2 years. In contrast to work from older children and adults, our findings suggest that the corpus callosum may be larger in infants who go on to develop autism spectrum disorder. This result was apparent with or without adjustment for total brain volume. Although we did not see a significant interaction between group and age, cross-sectional data indicated that area and thickness differences diminish by age 2 years. Regression data incorporating diffusion tensor imaging suggest that microstructural properties of callosal white matter, which includes myelination and axon composition, may explain group differences in morphology.



Investigating maternal brain structure and its relationship to substance use and motivational systems
H. J.V. Rutherford, G. Gerig, S. Gouttard, M. N. Potenza, L. C. Mayes. In Yale Journal of Biology and Medicine, in print, 2015.

Substance use during pregnancy and the postpartum period may have significant implications for both mother and the developing child. However, the neurobiological basis of the impact of substance use on parenting is less well understood. Here we examined the impact of maternal substance use on cortical gray matter (GM) and white matter volumes, and whether this was associated with individual differences in motivational systems of behavioral activation and inhibition. Mothers were included in the substance-using group if any addictive substance was used during pregnancy and/or in the immediate postpartum period (within 3 months of delivery). GM volume was reduced in substance-using mothers compared to non-substance-using mothers, particularly in frontal brain regions. In substance-using mothers, we also found that frontal GM was negatively correlated with levels of behavioral activation (i.e., the motivation to approach rewarding stimuli). This effect was absent in non-substance-using mothers. Taken together, these findings indicate a reduction in GM volume is associated with substance use, and that frontal GM volumetric differences may be related to approach motivation in substance-using mothers.



Bayesian Principal Geodesic Analysis for Estimating Intrinsic Diffeomorphic Image Variability
M. Zhang, P. T. Fletcher. In Medical Image Analysis (accepted), 2015.

In this paper, we present a generative Bayesian approach for estimating the low-dimensional latent space of diffeomorphic shape variability in a population of images. We develop a latent variable model for principal geodesic analysis (PGA) that provides a probabilistic framework for factor analysis in the space of diffeomorphisms. A sparsity prior in the model results in automatic selection of the number of relevant dimensions by driving unnecessary principal geodesics to zero. To infer model parameters, including the image atlas, principal geodesic deformations, and the effective dimensionality, we introduce an expectation maximization (EM) algorithm. We evaluate our proposed model on 2D synthetic data and the 3D OASIS brain database of magnetic resonance images, and show that the automatically selected latent dimensions from our model are able to reconstruct unobserved testing images with lower error than both linear principal component analysis (LPCA) in the image space and tangent space principal component analysis (TPCA) in the diffeomorphism space.



Prenatal Drug Exposure Affects Neonatal Brain Functional Connectivity
A. P. Salzwedel, K. M. Grewen, C. Vachet, G. Gerig, W. Lin,, W. Gao. In The Journal of Neuroscience, Vol. 35, No. 14, pp. 5860-5869. April, 2015.
DOI: 10.1523/JNEUROSCI.4333-14.2015

Prenatal drug exposure, particularly prenatal cocaine exposure (PCE), incurs great public and scientific interest because of its associated neurodevelopmental consequences. However, the neural underpinnings of PCE remain essentially uncharted, and existing studies in school-aged children and adolescents are confounded greatly by postnatal environmental factors. In this study, leveraging a large neonate sample (N = 152) and non-invasive resting-state functional magnetic resonance imaging, we compared human infants with PCE comorbid with other drugs (such as nicotine, alcohol, marijuana, and antidepressant) with infants with similar non-cocaine poly drug exposure and drug-free controls. We aimed to characterize the neural correlates of PCE based on functional connectivity measurements of the amygdala and insula at the earliest stage of development. Our results revealed common drug exposure-related connectivity disruptions within the amygdala–frontal, insula–frontal, and insula–sensorimotor circuits. Moreover, a cocaine-specific effect was detected within a subregion of the amygdala–frontal network. This pathway is thought to play an important role in arousal regulation, which has been shown to be irregular in PCE infants and adolescents. These novel results provide the earliest human-based functional delineations of the neural-developmental consequences of prenatal drug exposure and thus open a new window for the advancement of effective strategies aimed at early risk identification and intervention.



Finite-Dimensional Lie Algebras for Fast Diffeomorphic Image Registration
M. Zhang, P. T. Fletcher. In Information Processing in Medical Imaging (IPMI), 2015.

This paper presents a fast geodesic shooting algorithm for diffeomorphic image registration. We first introduce a novel finite-dimensional Lie algebra structure on the space of bandlimited velocity fields. We then show that this space can effectively represent initial velocities for diffeomorphic image registration at much lower dimensions than typically used, with little to no loss in registration accuracy. We then leverage the fact that the geodesic evolution equations, as well as the adjoint Jacobi field equations needed for gradient descent methods, can be computed entirely in this finite-dimensional Lie algebra. The result is a geodesic shooting method for large deformation metric mapping (LDDMM) that is dramatically faster and less memory intensive than state-of-the-art methods. We demonstrate the effectiveness of our model to register 3D brain images and compare its registration accuracy, runtime, and memory consumption with leading LDDMM methods. We also show how our algorithm breaks through the prohibitive time and memory requirements of diffeomorphic atlas building.



Spatio-temporal Image Analysis for Longitudinal and Time-Series Image Data,
S. Durrleman, T.P. Fletcher, G. Gerig, M. Niethammer, X. Pennec (Eds.). In Proceedings of the Third International Workshop, STIA 2014, Image Processing, Computer Vision, Pattern Recognition, and Graphics, Vol. 8682, Springer LNCS, 2015.
ISBN: 978-3-319-14905-9

This book constitutes the thoroughly refereed post-conference proceedings of the Third
International Workshop on Spatio-temporal Image Analysis for Longitudinal and Time-
Series Image Data, STIA 2014, held in conjunction with MICCAI 2014 in Boston, MA, USA, in
September 2014.

The 7 papers presented in this volume were carefully reviewed and selected from 15
submissions. They are organized in topical sections named: longitudinal registration and
shape modeling, longitudinal modeling, reconstruction from longitudinal data, and 4D
image processing.



Cell tracking using particle filters with implicit convex shape model in 4D confocal microscopy images
N. Ramesh, T. Tasdizen. In 2014 IEEE International Conference on Image Processing (ICIP), IEEE, Oct, 2014.
DOI: 10.1109/icip.2014.7025089

Bayesian frameworks are commonly used in tracking algorithms. An important example is the particle filter, where a stochastic motion model describes the evolution of the state, and the observation model relates the noisy measurements to the state. Particle filters have been used to track the lineage of cells. Propagating the shape model of the cell through the particle filter is beneficial for tracking. We approximate arbitrary shapes of cells with a novel implicit convex function. The importance sampling step of the particle filter is defined using the cost associated with fitting our implicit convex shape model to the observations. Our technique is capable of tracking the lineage of cells for nonmitotic stages. We validate our algorithm by tracking the lineage of retinal and lens cells in zebrafish embryos.



Multiatlas Segmentation as Nonparametric Regression
S.P. Awate, R.T. Whitaker. In IEEE Trans Med Imaging, April, 2014.
PubMed ID: 24802528

This paper proposes a novel theoretical framework to model and analyze the statistical characteristics of a wide range of segmentation methods that incorporate a database of label maps or atlases; such methods are termed as label fusion or multiatlas segmentation. We model these multiatlas segmentation problems as nonparametric regression problems in the high-dimensional space of image patches. We analyze the nonparametric estimator's convergence behavior that characterizes expected segmentation error as a function of the size of the multiatlas database. We show that this error has an analytic form involving several parameters that are fundamental to the specific segmentation problem (determined by the chosen anatomical structure, imaging modality, registration algorithm, and labelfusion algorithm). We describe how to estimate these parameters and show that several human anatomical structures exhibit the trends modeled analytically. We use these parameter estimates to optimize the regression estimator. We show that the expected error for large database sizes is well predicted by models learned on small databases. Thus, a few expert segmentations can help predict the database sizes required to keep the expected error below a specified tolerance level. Such cost-benefit analysis is crucial for deploying clinical multiatlas segmentation systems.



Joint Longitudinal Modeling of Brain Appearance in Multimodal MRI for the Characterization of Early Brain Developmental Processes
A. Vardhan, N. Sadeghi, C. Vachet, J. Piven, G. Gerig. In Spatiotemporal Image Analysis for Longitudinal and Time-Series Image Data (STIA'14) , LNCS. MICCAI'14, Springer Verlag, June, 2014.

Early brain maturational processes such as myelination manifest as changes in the relative appearance of white-gray matter tissue classes in MR images. Imaging modalities such as T1W (T1-Weighted) and T2W (T2-Weighted) MRI each display specific patterns of appearance change associated with distinct neurobiological components of these maturational processes. In this paper we present a framework to jointly model multimodal appearance changes across time for a longitudinal imaging dataset, resulting in quantitative assessment of the patterns of early brain maturation not yet available to clinicians. We measure appearance by quantifying contrast between white and gray matter in terms of the distance between their intensity distributions, a method demonstrated to be relatively stable to interscan variability. A multivariate nonlinear mixed effects (NLME) model is used for joint statistical modeling of this contrast measure across multiple imaging modalities. The multivariate NLME procedure considers correlations between modalities in addition to intra-modal variability. The parameters of the logistic growth function used in NLME modeling provide useful quantitative information about the timing and progression of contrast change in multimodal datasets. Inverted patterns of relative white-gray matter intensity gradient that are observable in T1W scans with respect to T2W scans are characterized by the SIR (Signal Intensity Ratio). The CONTDIR (Contrast Direction) which measures the direction of the gradient at each time point relative to that in the adult-like scan adds a directional attribute to contrast. The major contribution of this paper is a framework for joint multimodal temporal modeling of white-gray matter MRI contrast change and estimation of subject-specific and population growth trajectories. Results confirm qualitative descriptions of growth patterns in pediatric radiology studies and our new quantitative modeling scheme has the potential to advance understanding of variability of brain tissue maturation and to eventually differentiate normal from abnormal growth for early diagnosis of pathology.