Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Regional brain development analysis through registration using anisotropic similarity, a constrained affine transformation

  • Antoine Legouhy ,

    Roles Conceptualization, Data curation, Methodology, Software, Validation, Writing – original draft, Writing – review & editing

    antoine.legouhy@irisa.fr

    Affiliation CNRS, INRIA, INSERM, IRISA UMR 6074, Empenn ERL U-1228, Univ Rennes, Rennes, France

  • Olivier Commowick,

    Roles Methodology, Software, Supervision, Writing – review & editing

    Affiliation CNRS, INRIA, INSERM, IRISA UMR 6074, Empenn ERL U-1228, Univ Rennes, Rennes, France

  • Maïa Proisy,

    Roles Data curation, Funding acquisition, Writing – review & editing

    Affiliations CNRS, INRIA, INSERM, IRISA UMR 6074, Empenn ERL U-1228, Univ Rennes, Rennes, France, Radiology Department, CHU Rennes, Rennes, France

  • François Rousseau,

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation IMT Atlantique, LaTIM U1101 INSERM, UBL, Brest, France

  • Christian Barillot

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation CNRS, INRIA, INSERM, IRISA UMR 6074, Empenn ERL U-1228, Univ Rennes, Rennes, France

Abstract

We propose a novel method to quantify brain growth in 3 arbitrary orthogonal directions of the brain or its sub-regions through linear registration. This is achieved by introducing a 9 degrees of freedom (dof) transformation called anisotropic similarity which is an affine transformation with constrained scaling directions along arbitrarily chosen orthogonal vectors. This gives the opportunity to extract scaling factors describing brain growth along those directions by registering a database of subjects onto a common reference. This information about directional growth brings insights that are not usually available in longitudinal volumetric analysis. The interest of this method is illustrated by studying the anisotropic regional and global brain development of 308 healthy subjects betwen 0 and 19 years old. A gender comparison of those scaling factors is also performed for four age-intervals. We demonstrate through these applications the stability of the method to the chosen reference and its ability to highlight growth differences accros regions and gender.

Introduction

In pediatric image analysis, the study of brain development provides insights in the normal trend of brain evolution and enables early detection of abnormalities. Many types of morphometric measurements based on structural images have been explored and have shown their reliability as biomarkers in clinical use as established in [16]. Evaluated on a database of subjects covering a period of interest, it allows to better model the brain development and to highlight changes in growth, shape, structure, etc. Those measurements can be conducted on geometrical objects of different dimensions. Unidimensional ones such as the bicaudate ratio (minimum intercaudate distance divided by brain width along the same line) have been explored in [7] and [8] but also biparietal, bifrontal and transverse cerebellar diameters in [9]. Surface-based quantities such as cortical surface area in [10], corpus callosum mid-sagittal area in [11] or cortical folding metrics as in [12] have also been studied. However, the vast majority of studies are based on 3D features through the assessment of region of interest (ROI) volumes. Volumetric measures of different regions of the brain have been considered for specific ages or various temporal ranges. A far from exhaustive list is presented in Fig 1.

thumbnail
Fig 1. Ages covered by different brain volumetric studies, group comparison between: a—Normal / mental retardation, b—Male / female, c—Normal / motor disturbances, d—Normal / premature.

[11, 1330].

https://doi.org/10.1371/journal.pone.0214174.g001

Studied regions are very heterogeneous from large areas such as the whole brain itself, cerebellum, lobes or partitions of those to smaller ones such as basal ganglia, hippocampus, thalamus sometimes even separated according to the composition of their tissues (white matter (WM), gray matter (GM), cerebro-spinal fluid (CSF)). Some group comparisons have also been performed mostly between male and female or between preterm and term newborns.

Morphometric measurements can be determined manually. However, this requires the intervention of a medical expert able to select specific landmarks or perform segmentation. These tasks are highly time consuming with a potentially large inter-expert variability. Advances in computational medical imaging allow nowadays the use of semi-automated (requiring some human intervention) or fully-automated techniques. They involve algorithms able to automatically perform operations such as registration and segmentation.

A major drawback of purely volumetric measurements is that they do not provide any information on the shape of the regions or about the anisotropy of their development. In this paper, a new method is proposed that aims at quantifying global and regional brain growth in three arbitrary orthogonal directions of the brain (or ROI) through linear registration. To do so, a transformation called anisotropic similarity is introduced in section 1. It is an affine transformation with scaling directions constrained by orthogonal vectors arbitrarily chosen. A method to estimate, in a 3 dimensional space, the optimal anisotropic similarity for the least squares problem of distances between two sets of paired points is presented in section 2. Those results will then be used to create a registration algorithm based on this transformation. By registering a database of subjects onto a common basis (i.e. an atlas segmented in different ROIs) using anisotropic similarity, we have the opportunity to extract global or regional scaling ratios for all those subjects along arbitrary chosen orthogonal directions.

A direct application is, using the pipeline exposed in section 3, the exploration of regional scaling ratios growth charts along three fixed orthogonal directions through the ages highlighting anisotropic brain development. Resulting curves for whole brain and ROIs (lobes, basal ganglias, cerebellum…) are presented in section 4.2. A comparison of scaling factors from males and females is performed for four different age-intervals between 0 and 19 years old in section 4.3. Finally, the influence of the common reference image on the resulting scaling factors is studied in section 4.4.

Anisotropic similarity registration algorithm as well as other image processing tools used in this paper are publicly available in Anima (open source software for medical image processing: github.com/Inria-Visages/Anima-Public/ - RRID:SCR_017017).

1 Theoretical background

1.1 Generalities about linear transformations and anisotropic similarity

An affine transformation is a composition of a linear map A (N × N matrix) and a translation t (in ) operating on coordinates: y = Ax + t. Using singular value decomposition (SVD) on A, we obtain: (1) where W and V are unitary matrices and D is a positive diagonal matrix. By introducing R = V.Det(V), U = W.Det(W) and S = Det(V)Det(W)D, we get a modified decomposition: (2) where U is a rotation matrix defining the directions of scaling, S is a diagonal scaling matrix and R is a rotation matrix. We define a new linear transformation, hereafter named anisotropic similarity, which is an affine transformation with constrained directions of scaling. In other words, we define an anisotropic similarity transformation as an affine one where U is fixed. Summing up, we have the following in 3D space:

  • An affine transformation has 12 degrees of freedom:
    • a rotation (3 dof): the matrix U determines scaling directions.
    • an anisotropic scaling (3 dof): matrix S.
    • a rotation (3 dof): matrix R.
    • a translation (3 dof): vector t.
  • For an anisotropic similarity, the directions of scaling defined by U are constrained. This leaves 9 dof: 3 for rotation, 3 for scaling and 3 for translation.
  • For a similarity, the scaling part is constrained to have identical values on the diagonal i.e. S = s Id with leading to a linear part of the form sRUT. This leaves 7 dof: 3 for rotation, 1 for scaling and 3 for translation.
  • For a rigid transformation, the scaling part is constrained to identity leading to a linear part of the form RUT which is a rotation matrix since rotation matrices form a group for matrix multiplication. This leaves 6 dof: 3 for rotation and 3 for translation.

1.2 Generalities about linear registration

Registration consists in finding an optimal transformation that matches a moving image onto a reference image. This transformation is usually obtained by maximizing a similarity criterion. Many rigid (or linear in general) registration methods have been developed. They can be divided into two families: the ones that try to match geometrical features such as contours or surfaces, and those called iconic that are based on voxel intensities. Some of them use a global similarity measure between the two images such as mutual information in [31] and [32], while others rely on local similarities. Among this second category of approaches, block matching strategies exposed in [33] and [34] have gained in popularity. In those methods, two steps are iterated:

  1. Matching: for a set of blocks established in the reference image, homologous blocks best satisfying a similarity criterion are searched in the moving image.
  2. Aggregation into a global transformation: an optimization is performed in order to find the global transformation minimizing a distance between the sets of blocks and is then applied to the moving image. Usually, the weighted sum of squared euclidean distance is chosen for the cost function.

In order to perform an anisotropic similarity registration using the block-matching method, the two steps mentioned above have to be iterated. The first one (matching) is performed the same way it would be for any regular linear transformation. It outputs two sets of paired points: x and y that are in our case the centers of the homologous blocks. The second step (aggregation onto a global transformation) however is dependent on the type of linear transformation we want to determine leading to an adapted optimization in each case.

This optimization step consists in finding, in the set of eligible transformations, the one that best maps x onto y. Let x = {x1, …, xM} and y = {y1, …, yM} be two sets of M paired points coming from the matching step. For a global transformation with linear part A and translational part t, the least squares problem associated to the matching of x and y consists in the minimization of the following criterion C: (3)

Remark. For the sake of clarity we presented a version with a non-weighted least squares problem but the reasoning is the same with a weighted one.

The optimal translation can be directly obtained from the optimal linear part (independently of the type of linear transformation) and the barycenters of the two sets of points as developed in [35]. Let and , we have then: (4)

Let and be the barycentric coordinates. Since can be directly expressed from , the problem can then be simplified as the minimization of: (5)

In the case of the linear part being affine, there is no constraint. A closed form solution can therefore be easily found as shown in [35]. For rigid and similarity transformations, constraints lead to more complicated lagrangians but a closed form solution can be found as well using unit quaternions in 3D space as a representation of rotations like in [36] and [35].

2 Optimal anisotropic similarity between two sets of paired points

To our knowledge, the optimization procedure in the case of anisotropic similarities has not been considered in the literature. We thus present a method also based on quaternions to find the optimal anisotropic similarity between two sets of paired points. Writing an anisotropic similarity A as its decomposition: A = RSUT where U is fixed, the optimization over variable A therefore becomes an optimization over variables R and S: (6)

Let , . (7)

R can be expressed using quaternions following [36] and [35] and the problem then becomes (see S1.1): (8)

Where q is a unit quaternion and * is the quaternion multiplication. Let p and q be quaternions, and their conjugated quaternions respectively. There is a matricial representation of quaternions allowing to express quaternion product as a matrix product. Matricial quaternions P and Q are defined such that: Qpq = p * q and . (9)

Using those matricial quaternions on and ξi taken as pure quaternions, we have and . Thus, we obtain the following criterion (see S1.2): (10)

For further computation, we denote and B = ∑i Bi. A lagrangian with unit constraint qT q = 1 has then to be added to ensure a unit quaternion: (11)

Let sj = Sjj. The derivatives of this new formulation can then be written as: (12)

Derivative with respect to q depends upon sj and vice versa. Therefore, a direct solution to the problem of minimizing is difficult to find if not impossible. However, separating the problem between S and q leads to an alternate optimization scheme, each having an analytical solution.

Rotation: (13)

Solving this equation amounts finding the eigen vectors of B. More precisely, the global minimum is the one associated to the smallest eigen value of B as shown in [36], [35].

Scaling: (see S1.3) (14)

Now, interestingly the matrices have a quite trivial form. They are all symmetric, only the placing and indexes change (see S1.4). We finally get the following iterative alternate optimization scheme:

  • For a fixed value of , estimate the new optimal rotation quaternion: as the eigenvector with the smallest eigenvalue of B
  • For a fixed value of , estimate the new optimal scaling matrix following:

3 Material and methods

In the previous section, a method to find the optimal anisotropic similarity between two sets of paired points has been depicted. This gives the opportunity to register a database of subjects onto a common reference image using this type of linear transformation to extract scaling factors along chosen orthogonal directions and to study the variation of these factors on different ROIs between populations or among normal subjects.

3.1 Material

308 T1-weighted images of healthy subjects between 0 and 19 years old have been used, coming from three different studies: ASLpedia (section 6.1.1), C-MIND (section 6.1.2) and the Developing Human Connectome Project (dHCP) (section 6.1.3). Details on age repartition among databases and on image characteristics are given in Fig 2.

thumbnail
Fig 2. Repartition of the subjects selected from three studies over age.

https://doi.org/10.1371/journal.pone.0214174.g002

3.2 Methods

We developed a pipeline composed of 5 steps to extract scaling factors for 3 orthogonal directions on ROIs from a database of subjects.

  1. Choice and construction of the common reference image
  2. Segmentation of the common reference image into different ROIs
  3. Choice of the constrained directions of scaling for the anisotropic similarity registration
  4. Anisotropic similarity registration of a database of subjects to each ROIs of the common reference image to extract relative scaling factors
  5. Renormalization of the relative scaling factors to obtain absolute scaling factors

The above numbers associated to the steps are also associated to the subsections numbering below and to the numbers in Fig 3.

thumbnail
Fig 3. Pipeline for the extraction of scaling factors of a database of subjects using anisotropic similarity registration onto an atlas based on them as common reference image.

https://doi.org/10.1371/journal.pone.0214174.g003

3.2.1 Creation of the common reference image.

For genericity, the common reference image has been chosen to be an atlas made from all the subjects using a modified version of the atlas creation algorithm from [37]. The original method consists in the registration of a dataset of subjects onto a common reference image. Then, the inverse of the average of the local displacements is applied to the average of the registered images to obtain a new reference image. By iterating the process, the atlases (reference images), become less and less biased by the choice of the first reference image in terms of local displacements. However, the unbiasing steps does not take into account affine (global) transformations leading to atlases unbiased up to an affine transformation only. Our method, developed in [38] and available in Anima-Scripts (open source scripts for medical image processing: github.com/Inria-Visages/Anima-Scripts-Public/ - RRID:SCR_017072), takes advantage of the log-Euclidean framework developed in [39] and the Baker-Campbell-Hausdorff formula, mentioned in [40] and [41], allowing to average compositions of affine transformations and diffeomorphisms. This adjustment leads to the creation of atlases up to a rigid transformation. This methods includes in addition mechanisms conceived to improve the robustness of the registration algorithms in case of brains of different sizes through clever initialization and linear transformation decompositions.

3.2.2 Segmentation of the common reference image.

The atlas has then been segmented into 21 regions of interest (ROIs): whole brain, hemispheres, frontal lobe, parietal lobe, temporal lobe, occipital lobe, basal ganglias, cerebellum, insulas, ventricules, corpus callusum and brainstem. All structures were also separated in their left and right sides. To do this segmentation, ALBERTs manual ones ([42] and [43], see acknowledgments 6.1.4) have been used: 20 brains segmented into 50 regions manually drawn based on MRI brain scans that we fused to obtain the wider desired regions. The T1 weighted images of those brains have been registered onto our atlas through affine then diffeomorphic registration. The outputs have then been used to transfer all the segmentations onto our atlas which have been then merged using majority voting following [44]. The segmented atlas is shown in Fig 4. It should be noted that the 20 brains associated to ALBERTs data were only used for segmentation of the ROIs. They were not used for the creation of the common reference and, not used either as subjects registered with anisotropic similarity.

thumbnail
Fig 4. Regions of interest (ROIs) segmented and represented on the common reference image (top), chosen directions of scaling for anisotropic similarity registration defined and represented on the common reference image (bottom).

https://doi.org/10.1371/journal.pone.0214174.g004

3.2.3 Choice of the constrained directions of scaling.

The fixed scaling directions (characterized by the column vectors of the matrix U) are chosen on the reference image such that the first direction (blue in Fig 4) is orthogonal to the mid-sagittal plane (determined using [45]) for symmetry reasons. The other two directions are set using principal component analysis (PCA) on the non zero voxels coordinates projected onto the mid-sagittal plane. The second direction (red in Fig 4) corresponds to the principal direction from the PCA while the third (green in Fig 4) corresponds to the secondary one. three orthogonal directions are now chosen: one through iconic considerations and the other two based on purely geometrical features. In our application, the matrix U is the same for all ROIs of the reference image and is defined using the whole brain. However, it is possible to define a different U for each ROI independently. Chosen directions of scaling are shown in Fig 4.

These directions have been chosen as a proof of concept. Depending on the purpose of the study, other choice may appear more relevant (see Discussion Section 5).

3.2.4 Anisotropic similarity registration.

For each ROI, all subjects undergo an anisotropic similarity registration onto the reference image masked by this ROI. This registration is performed in two steps using in each case our block matching algorithm implemented in Anima (open source software for medical image processing: github.com/Inria-Visages/Anima-Public/ - RRID:SCR_017017):

  1. A similarity TB from whole brain subjects onto whole brain common reference is first estimated.
  2. An anisotropic similarity TC initialized from the previous step output is then computed to bring the subjects onto the atlas masked by the current ROI.

The first transformation, a similarity, is computed indirectly during a process of affine registration. Let A be the linear part of an affine transformation TA. We consider the following SVD: A = V DWT with D diagonal positive, V and W unitary matrices. We define TB (the nearest similarity associated to TA) as the transformation with linear part with being the average of the singular values namely the mean of the diagonal of D, and translation part . We chose the initialization to be a similarity since the composition of a similarity TB and an anisotropic similarity TC associated to a matrix U is still an anisotropic similarity associated to the same U: TBTC = (sBRB)(RCSCUT) = (RBRC)(sBSC)UT = RSUT. The whole registration process is summarized in Fig 5.

thumbnail
Fig 5. Two steps registration process: First an affine from which a nearest similarity is deduced, then an anisotropic similarity.

Transformations are composed and represented with arrows from destination to start since the interpolations occurring in the resampling process are done using the backward mapping. The inverse transformation is actually used on each voxel of the output images to determine the position in the input image from which a value is sampled.

https://doi.org/10.1371/journal.pone.0214174.g005

3.2.5 Extraction of absolute scaling factors.

From the output transformations of the registration step, the relative scaling parameters along the three fixed directions are extracted. Those scaling factors are relative to the reference image. We want to normalize them such that they become about equal to 1 at birth and represent, for other ages, how much a region expanded along the chosen directions since birth. To this end, the fact that all dHCP subjects are very young (less than 1 month) is exploited. The normalization step consist in dividing the relative scaling factors of all the dataset by the average of the ones associated to the dHCP subjects considered as the “root” of the brain expansion. Those new scaling factors will now be considered as absolute scaling factors. At this stage, for each subject, an absolute scaling factor has been determined for each ROI. Those absolute scaling factors are used to model the expansion of the brain toward the chosen directions.

4 Experiments and results

4.1 Model selection

Several models are traditionally used to represent growth in biostatistics such as the exponential or Weibull models. The second one has been considered by [30] as the best suited to model brain growth in terms of volume. Our case however is different: it can be viewed as a 3-way unidimensional approach. In our quest to find the function best suited to model growth curves for our data, we decided to consider, as a prior, that the brain expansion is stopping at some point. Therefore, we restricted the spectrum to functions that have a horizontal asymptote at infinity. The selected candidates to model brain growth in the chosen directions are the following:

  • Rational with polynomials of degree 1 as numerator and denominator:
  • Weibull:
  • Gompertz:
  • Exponential: y = a + becx

For each candidate, the optimal coefficients are estimated through nonlinear regression using the Levenberg-Marquardt iterative weighted least squares algorithm from [46]. In this process, weights are chosen to compensate for local gender repartition. For each subject i, a window of width l = 2 years centered on the subject age is considered. Let nf, nm and n be the number of female, male and total subjects respectively in that window. A correction coefficient is applied if i is a female and if i is a male. Let {y1, …, yn} be the observations (i.e. here the obtained scaling factors), be the average of those and be the fitted values.

Based on these statistics, the chosen candidate for the modeling will be the one that best satisfies a criterion quantifying the goodness of fit. This indicator should evaluate the accuracy of the model i.e. how close the model is to the observation while discouraging overfitting. It therefore consists in a tradeoff between accuracy and parsimony. It has been shown in [47] that the coefficient of determination is not, at least when considered alone, an appropriate measure for the goodness of fit in the case of nonlinear model selection. A more adapted goodness of fit for nonlinear model selection is the Akaike information criterion (AIC) developed in [48] and [49]. Based on information theory, it proposes to estimate the information loss induced by each candidate model to represent an unknown process that supposedly generated the data as shown in [50]. This is made possible through the estimation of the Kullback-Leibler divergence related to the maximized log-likelihood. AIC is defined by: , where is the maximum likelihood of the candidate model and the first term penalizes a large number of parameters p. Therefore, the preferred model among the candidates is the one with the lower AIC. Note that AIC of a model taken alone is meaningless, it makes sense only when compared to the one of the other models. This is why it is recommended to consider it along with another statistic that quantifies the error between the model and the data like mean of squared errors (MSE): which is the average of the residuals. A corrected version of the AIC has been developed to avoid overfitting in the case of small sample sizes: . To facilitate the interpretation that can be quite obscure using raw AIC, following [51], it is possible to transform those values into conditional probabilities for each model called Akaike weights. Defined for each model i by , those weights represent the probability for each candidate i to be the best suited in the sense of AIC to model the data among all the candidates.

All the goodness of fit depicted above as well as MSE have been evaluated for each of our candidates to model the scaling factors for each ROI. We present the results of this evaluation in Fig 6. The Gompertz and exponential models are largely below the other two. Even though the Weibull model behaves relatively well, the rational one shows better scores whatever the tested goodness of fit.

thumbnail
Fig 6. Goodness of fits for each candidate to model the ouputed scaling factors averaged in the 3 directions.

Boxplots are performed along the ROIs, ROI IDs are displayed on the left. Akaike weights are computed on mean (blue) and median (red) AICc.

https://doi.org/10.1371/journal.pone.0214174.g006

4.2 Directional growth curves

From the previous sections, scaling factors in each direction for each ROI are now modeled using a rational function with polynomials of degree 1 as numerator and denominator chosen after model selection. Results for all regions studied are presented in Figs 713. The method presented by [52] is used to compute simultaneous 99% confidence intervals for fitted values. The black curve represents the average brain growth computed as the mean of the directional models (Figs 713).

thumbnail
Fig 7. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g007

thumbnail
Fig 8. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g008

thumbnail
Fig 9. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g009

thumbnail
Fig 10. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g010

thumbnail
Fig 11. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g011

thumbnail
Fig 12. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g012

thumbnail
Fig 13. Resulting scaling factors as a function of the age in years for differents ROIs, along direction 1 (blue), 2 (red), 3 (green).

Fitted using rational model together with 99% confidence intervals for fitted values, Black curves represent the average models along the three directions.

https://doi.org/10.1371/journal.pone.0214174.g013

4.3 Male vs female comparison

Gender, like age, is a characteristic of the subjects available in all the studies we considered. The aim of this section is to evaluate if differences in terms of scaling factor can be found between genders. We divided our data into four age-intervals based on the age of the subjects. The first one contains dHCP participants (newborns), the second one is composed of all non-newborn subjects between 0 and 6 years old, the third one between 6 and 12 and the fourth older to 12. Repartition of the subjects in terms of gender, age class and study is shown in Table 1.

thumbnail
Table 1. Repartition of the subjects in term of age class, gender and study.

https://doi.org/10.1371/journal.pone.0214174.t001

For each of these age-intervals, and each of the chosen scaling directions, and each ROI, we performed a test to evaluate if the scaling factors for male subjects are greater than scaling factors for female subjects. Since these data are not normally distributed in those subdivisions, we used two-tailed Wilcoxon-Man-Whitney U-tests. For each of those tests, the null hypothesis H0 is the following: the distribution of the scaling factors between males and females are equal, while the alternative hypothesis H1 states: the distributions of males and females are different. We performed 252 tests in total: 4 age-intervals × 21 ROIs × 3 directions whose results are shown in Fig 14.

thumbnail
Fig 14. Male vs female comparison using Wilcoxon-Man-Whitney U-test and H0: The distribution of the scaling factors of males and females are equal, H1: The distributions of males and females are different.

In color: p-values for H0 rejection for FDR at level 5% (Benjamini and Hochberg method). Numerically: the size of the effect d for each test.

https://doi.org/10.1371/journal.pone.0214174.g014

A type 1 error, or false positive, occurs when H0 is incorrectly rejected. Since we are doing multiple comparisons, rejecting H0 based on the risk of type 1 error α = 5%, may lead in our case to an expected number of false positives greater than 12. Instead of using α, we therefore adopted the false discovery rate (FDR) that controls the proportion of false positives among the tests where H0 has been rejected. Therefore, we stated the acceptance or rejection of H0 based on a FDR at level 5%. This has been done using Benjamini and Hochberg procedure from [53] and corresponds to reject H0 when the p-value is less than 0.0077 (Fig 14). FDR has been preferred to family-wise error rate (FWER), that controls the risk of at least 1 false positive among the whole family of tests, because of the over-conservatism of this last type of procedure leading to poor test power (probability of correctly rejecting H0). Additionally, we calculated, for each test, the effect size d following: (Fig 14), where {Sm} (resp. {Sf}) is the set of scaling factor of males (resp. females) used for the test. We preferred the use of median instead of mean due to the fact that we do not know the distribution of the data a priori and we performed ranksum type tests.

For all the tests that lead to a rejection of the null hypothesis, scaling factors were higher for males both in terms of means and medians. Tests show that scaling factors of males seem higher in the second age class (0-6), brainwise and mainly in temporal and cerebellum areas along the direction 1. This is also notable in the same regions between 6 and 12 years, this time along direction 3. For the older class (12-19), this phenomenon essentially appears brainwise along the direction 3 and in the parietal lobes along direction 1.

4.4 Influence of the common reference

To evaluate the influence of the common reference image, the whole process described previously is reproduced using six different reference images. Those are atlases for six arbitrary chosen time-points t1, …, t6 using subjects among the population depicted in Section 3.1. Time-points are chosen such that five of them cover the period in which the majority of the brain expansion occurs, the last is positioned later, in a stabilized area. Following [38], an atlas for a time-point ti is created using only the subjects with ages in a range around ti (in Fig 15, line 4). The contributions of those subjects are then weighted using a weight function (plotted in Fig 15, line 6) such that subjects closer in age to ti are given more importance. The six outputed atlases (presented Fig 15) are different in terms of shape and intensity, they will be used to challenge the robustness of our method to extract scaling factors.

thumbnail
Fig 15. Characteristics of the 6 atlases used as reference image (time is displayed in years).

https://doi.org/10.1371/journal.pone.0214174.g015

The method developed in section 3.2 is used for each of these reference images, on which directions of scaling for the anisotropic similarity registration have been established the same way. A scaling factor s(i, j, d, k) is thus computed for each ROI j of each subject i based on each reference image k along each chosen direction d. To quantify the influence of the reference image on absolute scaling factors, the results, using the six reference images previously depicted, are compared through two approaches:

  1. A pairwise study to evaluate whether or not reference atlases closer to each other in age are more likely to generate closer results.
  2. A study of the standard deviation among results for all reference atlases to evaluate how far they are from the average results.

4.4.1 Study of pairwise distances between scaling factors by reference images in each direction for each ROI.

Our aim is to determine whether or not reference images closer to each other (atlases at shorter temporal distances) are more likely to generate less important absolute differences between their results. We therefore compute the absolute difference of the resulting scaling factors between each pairwise combinations of reference images. Then, those distances are normalized by the average of corresponding scaling factors between the two atlases such that it can be seen as a percentage of it (relative distance). The relative distance between scaling factors from reference atlases k and l is then computed as: (15)

After examination of all the pairwise combinations, the temporal distance between the reference images does not seem to have an impact on the distance of the scaling ratios associated to each other (Fig 16). The highest median of relative distance happens to be between atlases 2 and 5 for right basal ganglia, but does go above 8% of difference.

thumbnail
Fig 16. Relative distances between reference atlas 1 and 2 (top), 1 and 6 (bottom).

Boxplots among subjects for each ROI j, each direction d: boxplot(Dk,l(., j, d)).

https://doi.org/10.1371/journal.pone.0214174.g016

4.4.2 Study of the standard deviation among reference images in each direction for each ROI.

This method gives an average measure of the distances between the results for each atlas and the average results. Those distances are normalized by the average of corresponding scaling factors of all the atlases. The relative standard deviation between scaling from all reference atlases is then computed as: (16)

The graphs (Fig 17) suggest that the method, when applied to large regions such as whole brain and hemispheres, is really robust to reference image change. Occipital lobes and cerebellum however seem to be more vulnerable areas. Those two regions share a common border and we believe that the segmentation process is a crucial step in that case. The cerebellum position indeed varies quickly in early stages of life and our decision to use segmentations based on neonates can be a bit inadequate for this area in particular. We also think that the way we chose to define the constrained directions of scaling (especially those using purely geometrical considerations through PCA on voxel coordinates) may not be the best suited for robustness in those areas. More anatomical features could lead to even smaller influence of the reference image.

thumbnail
Fig 17. Relative standard deviation between reference atlases.

Boxplots among subjects for each ROI j, each direction d: boxplot(D(., j, d)).

https://doi.org/10.1371/journal.pone.0214174.g017

5 Discussion

5.1 Comparison of the results to the literature

Our method aims at characterizing directional growth at the scale of a ROI through the extraction of scaling factors, which are easily interpretable measures. In most longitudinal analyses of brain development, the focus is on volumetry. This type of measurement reflects the global growth of a region but lacks characterization of directional expansion. Other approaches are based on cortical surface curvature as in [12], [54] or [55] to reflect gyri and sulci widening or deepening, giving information about directional expansion at a more local scale.

Values from volumetric measurements cannot be directly compared to our scaling factors: they are three-times unidimensional whereas volume is not. Our anisotropic scaling factors should instead be seen as a complementary information to the volume data. However, there are still some aspects that we can compare:

  • The pace of growth: is the growth monotonous or does it increase or decrease in speed at given ages? Are there regions that seem to evolve at a quicker pace than the others do?
  • The asymptote: at what age does the growth reach a peak and begins to stagnate? Are there regions that seem to reach a higher growth asymptote than others do?

Since the anisotropy is a missing information in volumetric studies, it is appropriate to use the average models of the three directions (black curves in Fig 7) for the comparisons. In [13], [30], [23] and [18], the pace of the whole brain volume curves is quite similar to our findings, showing a strong initial increase followed by a decline around 2 years. In [25], it is also noticeable that, analogously to our results, the cerebellum grows by a much more important factor than the brain as a whole. In [18], basal ganglia grow by a factor comparable to the whole brain which is also consistent with our results.

5.2 Implementation considerations

5.2.1 Directions of scaling.

The choice of the scaling directions, characterized by the matrix U, is crucial. These three directions can be different for each ROI, they are defined on the reference image and must be orthogonal. The brain being a rather symmetrical organ about its mid-sagittal plane, choosing a direction orthogonal to it seems natural. Given that, there is no obvious choice for the other two. Using PCA on non-zeros voxel coordinates of the reference image (atlas of the population in our case) projected on mid-sagittal plane (since the last two directions have to be orthogonal to the first one) seemed to be quite intuitive. This corresponds to the principal sagittal directions of the average model of the population. However, it is based on purely geometric features ignoring iconic or anatomical considerations. This choice must above all be seen as a proof of concept and should not be considered as the one to be necessarily adopted. A future work could be to define, for a given purpose, a criterion of goodness for the chosen directions and, to optimize to find the directions that best satisfy this criterion. A more anatomically-oriented approach could be to ask a medical expert to point, on the reference image, the anterior commissure—posterior commissure (AC-PC) line. This has been well adopted as a standard by the clinical neuroimaging community even though it is mostly a convention for visualization and at the cost of introducing a human interaction or a preprocessing step. This method is therefore very flexible in the choice of the scaling directions and the ROIs, yet it has shown oneself robust concerning the choice of the common reference image.

5.2.2 ROI segmentation.

Although it does not call into question the method itself, there is room for improvements in the way we segmented the ROIs. The main difficulty is to find a method that is reproducible while being adaptable to brains from subjects scanned across a wide range of ages, which induces a large variability in contrast and shape. It has been shown in [56] that the automatic segmentation of 1 year old subjects was more successful when using prior segmentations of 2 years old than when using prior segmentations of adults. In addition, [57] showed that, in the context of multi-atlas segmentation, similarity selection of the atlases increases significantly the accuracy compared to the fusion of random sets of atlases. In the same paper, it has also been established that using an age-based selection gives similar results to similarity selection. Those results support the assertion that segmentation accuracy is improved by using age-adapted atlases. However, the concern of producing segmentations on a longitudinal database that are adapted to the maturation of each brain while being consistent among the subjects is still, to our knowledge, an open question. Our choice to use ALBERTs atlases to segment all the subjects regardless of their age was driven by a desire of consistency/reproducibility to the expense of age-adaptation since we did not find methods offering a satisfying trade-off between the two.

5.2.3 Normalization step.

We have used a normalization step instead of the raw scaling factors relative to the reference image. Even though it is in a way more abstract, the relative scaling factors can actually be used as such. A subject with abnormalities that would show itself atypical in terms of scaling factors compared to the normal model would be so, with or without normalization. Now it is indeed more tangible for interpretation to have a concrete base like birth. In that case, scaling factors have an intuitive meaning since they quantify how much a region grows along a given direction since birth. In the event that there is no neonatal data to normalize the relative results, the following could be done:

  • If a model can be fitted from the relative scaling factors, it might be possible to extrapolate values for birth and therefore use them for normalization.
  • If not, relative scaling factors may be used or another base should be considered. Taking values in the growth stabilization zone may also be considered as it gives easily interpretable results as well.

That being said, neonatal MRI data are nowadays quite abundant in open access, especially structural images. The dHCP team (see Acknowledgment 6.1.3) announced in September their second release with over 500 neonates. Our normalization can therefore be performed when using this data.

5.2.4 Image quality.

Images used in the experiments are of good quality compared to what can be commonly found in clinical routine. Registration is sensitive to image quality but it is mostly the case in the context of non-linear registration where there are many degrees of freedom and very local displacements are involved. Even if the image quality can have an impact, the sparsity in terms of degrees of freedom and the global aspect of the anisotropic registration make it much more resilient about image size and quality. In addition, the bigger is the size of the ROI is, the less sensitive to image quality the result will be. The registration algorithm being based on block-matching, the ROI must be chosen such that their sizes are at least one order larger than the size of the blocks. A block being a neighborhood of voxels, the poorer the image quality, the larger the ROI must be and too narrow structures should be avoided. Finally, conducting the study at a group scale attenuates the influence of potential registrations errors induced by a poor image quality.

5.3 Clinical applications

We focused on the expansion of structures of a database of healthy subjects but we can also imagine using this method for patients. Intra-individual surveys are also possible, for subjects that had multiple scans through time, to monitor the evolution of a brain sub-region or any part of the body and infer the way it is going to expand.

Brain morphometric measurements have been performed on infants with deformational plagiocephaly and controls in [58]. They found that volume analysis does not allow the distinction between cases and controls. However, their unidimensional measurements based especially on distances between landmarks were more successful. Those nevertheless necessitate point picking on all subjects and may be sensitive to human error. Considering this, it could be interesting to try, in future work, our automatic method (only the computation of directions for each ROI on the common reference is needed) on this kind of pathology. Some chosen directions could for example be colinear to the line segments that have proved to be relevant in this paper. Those measurements might be linked to developmental delays that affect infants with deformational plagiocephaly.

Conclusion

We have presented a method that allows the extraction of regional and global scaling factors along arbitrary chosen orthogonal directions. This is done through linear registration using a 9 dof transformation, anisotropic similarity, which is an affine transformation with constrained directions of scaling.

The main methodological contribution of this paper concerns the resolution of the problem of finding the optimal anisotropic similarity that best matches two sets of paired points. This result has made possible the development of a block-matching registration algorithm based on this transformation.

Given this new type of registration, our second contribution was to map a database of subjects between 0 and 19 years old using anisotropic similarity onto a common reference image on which the constrained directions of scaling of our choosing have been fixed. Registrations have been performed brainwise and ROI wise (lobes, cerebellum, basal ganglias…). For genericity purpose, we chose this reference image to be an atlas made from the subjects. Based on symmetry and geometrical considerations, we defined the same constrained directions of scaling for all ROIs even though it is possible to choose different ones for each. As an output, we obtained for each subject, for each ROI, for each chosen direction a scaling, a scaling factor that we normalized such that it represents an expansion factor from birth.

Those scaling factors have been used to model the anisotropic development of the brain. After model selection, it has been determined that rational function with polynomials of degree 1 as numerator and denominator is the best suited among the tested candidates for that modeling. Curves representing scaling factors as a function of the age for each ROIs, each chosen direction, along with associated confidence intervals have then beeen computed on a combination of four databases.

Tests to determine the influence of gender in those scaling factors have been performed for different age-intervals. Finally, two experiments have been conducted to evaluate the influence of the aforementioned common reference image. The results have shown small relative differences depending on the choice of the reference image leading to the conclusion that the method is robust in that aspect.

6 Acknowledgments

6.1 Databases

6.1.1 ASLpedia.

A retrospective ASL study on routine pediatric MRI performed at Rennes University Hospital Neuropediatric radiology Department.

6.1.2 C-MIND.

Data Repository created by the C-MIND study of Normal Brain Development. This is a multisite, longitudinal study of typically developing children from ages newborn through young adulthood conducted by Cincinnati Children’s Hospital Medical Center and UCLA and supported by the National Institute of Child Health and Human Development (Contract HHSN275200900018C). A listing of the participating sites and a complete listing of the study investigators can be found at: https://research.cchmc.org/c-mind.

6.1.3 The Developing Human Connectome Project (dHCP).

Led by King’s College London, Imperial College London and Oxford University, aims to make major scientific progress by creating the first 4-dimensional connectome of early life.

https://data.developingconnectome.org/.

6.1.4 ALBERTs.

See [42] and [43] for details about segmentations. Copyright Imperial College of Science, Technology and Medicine and Ioannis S. Gousias 2013. All rights reserved. https://brain-development.org/brain-atlases/neonatal-brain-atlases/neonatal-brain-atlas-gousias/.

Supporting information

S1 Appendix. Detailed calculations for optimizing anisotropic similarity between two sets of paired points.

https://doi.org/10.1371/journal.pone.0214174.s001

(PDF)

References

  1. 1. Dubois J, Lefèvre J, Angleys H, Leroy F, Fischer C, Lebenberg J, et al. The dynamics of cortical folding waves and prematurity-related deviations revealed by spatial and spectral analysis of gyrification. NeuroImage. 2019;185(February 2018):934–946. pmid:29522888
  2. 2. Giorgio A, De Stefano N. Clinical use of brain volumetry. Journal of Magnetic Resonance Imaging. 2013;37(1):1–14. pmid:23255412
  3. 3. Durazzo TC, Tosun D, Buckley S, Gazdzinski S, Mon A, Fryer SL, et al. Cortical thickness, surface area, and volume of the brain reward system in alcohol dependence: Relationships to relapse and extended abstinence. Alcoholism: Clinical and Experimental Research. 2011;35(6):1187–1200.
  4. 4. Frisoni GB, Fox NC, Jack CR, Scheltens P, Thompson PM. The clinical use of structural MRI in Alzheimer disease; 2010.
  5. 5. Novak K, Czech T, Prayer D, Dietrich W, Serles W, Lehr S, et al. Individual variations in the sulcal anatomy of the basal temporal lobe and its relevance for epilepsy surgery: an anatomical study performed using magnetic resonance imaging. Journal of Neurosurgery. 2009;96(3):464–473.
  6. 6. Tosun D, Siddarth P, Levitt J, Caplan R. Cortical thickness and sulcal depth: insights on development and psychopathology in paediatric epilepsy. BJPsych Open. 2015;. pmid:27703737
  7. 7. Aylward EH, Schwartz J, Machlin S, Pearlson G. Bicaudate ratio as a measure of caudate volume on MR images. American Journal of Neuroradiology. 1991;12(6):1217–1222. pmid:1763757
  8. 8. Bermel RA, Bakshi R, Tjoa C, Puli SR, Jacobs L. Bicaudate ratio as a magnetic resonance imaging marker of brain atrophy in multiple sclerosis. Archives of Neurology. 2002;59(2):275–280. pmid:11843699
  9. 9. Tich S Nguyen, Anderson PJ, Hunt RW, Lee KJ, Doyle LW, Inder TE. Neurodevelopmental and perinatal correlates of simple brain metrics in very preterm infants. Archives of Pediatrics and Adolescent Medicine. 2011;165(3):216–222.
  10. 10. Winkler AM, Sabuncu MR, Yeo BTT, Fischl B, Greve DN, Kochunov P, et al. Measuring and comparing brain cortical surface area and other areal quantities. NeuroImage. 2012;61(4):1428–1443. pmid:22446492
  11. 11. Giedd JN, Rapoport JL. Structural MRI of Pediatric Brain Development: What Have We Learned and Where Are We Going? Neuron. 2010;67(5):728–734. pmid:20826305
  12. 12. Lefèvre J, Germanaud D, Dubois J, Rousseau F, De MacEdo Santos I, Angleys H, et al. Are developmental trajectories of cortical folding comparable between cross-sectional datasets of fetuses and preterm newborns? Cerebral Cortex. 2016;26(7):3023–3035. pmid:26045567
  13. 13. Hamano K, Iwasaki N, Kawashima K, Takita H. Volumetric quantification of brain volume in children using sequential CT scans. Neuroradiology. 1990;32(4):300–303. pmid:2234389
  14. 14. Pfefferbaum A, Mathalon DH, Sullivan EV, Rawles JM, Zipursky RB, Lim KO. A Quantitative Magnetic Resonance Imaging Study of Changes in Brain Morphology From Infancy to Late Adulthood; 1994.
  15. 15. Reiss AL, Abrams MT, Singer HS, Ross JL, Denckla MB. Brain development, gender and IQ in children. A volumetric imaging study. Brain: a journal of neurology. 1996;119 (Pt 5):1763–74.
  16. 16. Matsumae M, Kikinis R, Mórocz IA, Lorenzo AV, Sándor T, Albert MS, et al. Age-related changes in intracranial compartment volumes in normal adults assessed by magnetic resonance imaging. Journal of Neurosurgery. 1996;84(6):982–991. pmid:8847593
  17. 17. Murphy DG, DeCarli C, McIntosh AR, Daly E, Mentis MJ, Pietrini P, et al. Sex differences in human brain morphometry and metabolism: an in vivo quantitative magnetic resonance imaging and positron emission tomography study on the effect of aging. Archives of general psychiatry. 1996;53(7):585–94. pmid:8660125
  18. 18. Iwasaki N, Hamano K, Okada Y, Horigome Y, Nakayama J, Takeya T, et al. Volumetric quantification of brain development using MRI. Neuroradiology. 1997;39(12):841–846. pmid:9457706
  19. 19. Hüppi PS, Warfield S, Kikinis R, Barnes PD, Zientara GP, Jolesz Fa, et al. Quantitative magnetic resonance imaging of brain development in premature and mature newborns. Annals of neurology. 1998;43(2):224–235. pmid:9485064
  20. 20. Coffey CE. Sex Differences in Brain Aging. Arch Neurol. 1998;. pmid:9482358
  21. 21. Giedd JN, Blumenthal J, Jeffries NO, Castellanos FX, Liu H, Zijdenbos A, et al. Brain development during childhood and adolescence: a longitudinal MRI study. Nature Neuroscience. 1999;2(10):861–863. pmid:10491603
  22. 22. Utsunomiya H, Takano K, Okazaki M, Mitsudome A. Development of the temporal lobe in infants and children: Analysis by MR-based volumetry. American Journal of Neuroradiology. 1999;20(4):717–723. pmid:10319988
  23. 23. Courchesne E, Chisum HJ, Townsend J, Cowles A, Covington J, Egaas B, et al. Normal Brain Development and Aging: Quantitative Analysis at in Vivo MR Imaging in Healthy Volunteers. Radiology. 2000;216(3):672–682. pmid:10966694
  24. 24. Gousias IS, Rueckert D, Heckemann RA, Dyet LE, Boardman JP, Edwards AD, et al. Automatic segmentation of brain MRIs of 2-year-olds into 83 regions of interest. NeuroImage. 2008;40(2):672–684. pmid:18234511
  25. 25. Knickmeyer RC, Gouttard S, Kang C, Evans D, Smith JK, Hamer RM, et al. A Structural MRI Study of Human Brain Development from Birth to 2 Years. Journal of Neuroscience. 2010;28(47):12176–12182.
  26. 26. Kuklisova-Murgasova M, Aljabar P, Srinivasan L, Counsell SJ, Doria V, Serag A, et al. A dynamic 4D probabilistic atlas of the developing brain. NeuroImage. 2011;. pmid:20969966
  27. 27. Lange N. Total and regional brain volumes in a population-based normative sample from 4 to 18 years: The NIH MRI study of normal brain development. Cerebral Cortex. 2012;22(1):1–12.
  28. 28. Choe Ms, Ortiz-Mantilla S, Makris N, Gregas M, Bacic J, Haehn D, et al. Regional Infant Brain Development: An MRI-Based Morphometric Analysis in 3 to 13 Month Olds. Cerebral Cortex. 2013;23(9):2100–2117. pmid:22772652
  29. 29. Makropoulos A, Aljabar P, Wright R, Hüning B, Merchant N, Arichi T, et al. Regional growth and atlasing of the developing human brain. NeuroImage. 2016;125:456–478. pmid:26499811
  30. 30. Peterson M, Warf BC, Schiff SJ. Normative human brain volume growth. Journal of Neurosurgery: Pediatrics. 2018;21(May):1–8.
  31. 31. Wells WM, Viola P, Atsumi H, Nakajima S, Kikinis R. Multi-Modal Volume Registration by Maximisation of Mutual Information. Medical Image Analysis. 1996;.
  32. 32. Maes F, Collignon A, Vandermeulen D, Marchal G, Suetens P. Multimodality image registration by maximization of mutual information. IEEE Transactions on Medical Imaging. 1997;16(2):187–198. pmid:9101328
  33. 33. Ourselin S, Roche A, Prima S, Ayache N. Block Matching: A General Framework to Improve Robustness of Rigid Registration of Medical Images. Springer Berlin Heidelberg; 2000. p. 557–566. Available from: http://link.springer.com/10.1007/978-3-540-40899-4_57.
  34. 34. Commowick O, Wiest-Daessle N, Prima S. Block-matching strategies for rigid registration of multimodal medical images. In: Proceedings—International Symposium on Biomedical Imaging; 2012. p. 700–703.
  35. 35. Pennec X. L’incertitude dans les problèmes de reconnaissance et de recalage—Applications en imagerie médicale et biologie moléculaire. Ecole Polytechnique X; 1996.
  36. 36. Horn BKP. Closed-form solution of absolute orientation using unit quaternions. Journal of the Optical Society of America A. 1987;4(4):629.
  37. 37. Guimond A, Meunier J, Thirion JP. Average brain models: A convergence study. Computer vision and image understanding. 2000;77(2):192–210.
  38. 38. Legouhy A, Commowick O, Rousseau F, Barillot C. Unbiased longitudinal brain atlas creation using robust linear registration and log-Euclidean framework for diffeomorphisms. In: Proceedings—International Symposium on Biomedical Imaging. IEEE; 2019. p. 1038–1041.
  39. 39. Arsigny V, Commowick O, Pennec X, Ayache N. A Log-Euclidean Framework for Statistics on Diffeomorphisms. Springer Berlin Heidelberg; 2006. p. 924–931. Available from: http://link.springer.com/10.1007/11866565_113.
  40. 40. Vercauteren T, Pennec X, Perchant A, Ayache N. Symmetric Log-Domain Diffeomorphic Registration: A Demons-Based Approach. In: Medical Image Computing and Computer Assisted Intervention. Springer Berlin Heidelberg; 2008. p. 754–761. Available from: http://link.springer.com/10.1007/978-3-540-85988-8_90.
  41. 41. Bossa M, Hernandez M, Olmos S. Contributions to 3D diffeomorphic atlas estimation: Application to brain images. Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics). 2007;4791 LNCS(PART 1):667–674.
  42. 42. Gousias IS, Edwards AD, Rutherford MA, Counsell SJ, Hajnal JV, Rueckert D, et al. Magnetic resonance imaging of the newborn brain: Manual segmentation of labelled atlases in term-born and preterm infants. NeuroImage. 2012;62(3):1499–1509. pmid:22713673
  43. 43. Gousias IS, Hammers A, Counsell SJ, Srinivasan L, Rutherford MA, Heckemann RA, et al. Magnetic Resonance Imaging of the Newborn Brain: Automatic Segmentation of Brain Images into 50 Anatomical Regions. PLoS ONE. 2013;8(4). pmid:23565180
  44. 44. Heckemann RA, Hajnal JV, Aljabar P, Rueckert D, Hammers A. Automatic anatomical brain MRI segmentation combining label propagation and decision fusion. NeuroImage. 2006;33(1):115–126. pmid:16860573
  45. 45. Prima S, Ourselin S, Ayache N. Computation of the mid-sagittal plane in 3-D brain images. IEEE Transactions on Medical Imaging. 2002;21(2):122–138. pmid:11929100
  46. 46. Seber GAF, Wild CJ. Nonlinear Regression; 1989.
  47. 47. Spiess AN, Neumeyer N. An evaluation of R2 as an inadequate measure for nonlinear models in pharmacological and biochemical research: A Monte Carlo approach. BMC Pharmacology. 2010;10:1–11.
  48. 48. Akaike H. A new look at the statistical model identification. IEEE Transactions on Automatic Control. 1974;19(6):716–723.
  49. 49. Burnham KP, Anderson DR. Model Selection and Multimodel Inference; 2002.
  50. 50. Burnham KP, Anderson DR. Multimodel inference: Understanding AIC and BIC in model selection. Sociological Methods and Research. 2004;33(2):261–304.
  51. 51. Wagenmakers EJ, Ferrel S. AIC model selection using Akaike weights. Psychonomic Bulletin & Review. 2004;11(1):192–196.
  52. 52. Lane TP, DuMouchel WH. Simultaneous Confidence Intervals in Multiple Regression. The American Statistician. 1994;48(4):315–321.
  53. 53. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society. 1995;57(1):289–300.
  54. 54. Auzias G, De Guio F, Pepe A, Rousseau F, Mangin JF, Girard N, et al. Model-driven parameterization of fetal cortical surfaces. Proceedings—International Symposium on Biomedical Imaging. 2015;2015-July:1260–1263.
  55. 55. Kim SH, Lyu I, Fonov V, Vachet C, Hazlett HC, Smith RG, et al. Development of cortical shape in the human brain from 6 to 24months of age via a novel measure of shape complexity. NeuroImage. 2016;135:163–176. pmid:27150231
  56. 56. Gousias IS, Hammers A, Heckemann RA, Counsell SJ, Dyet LE, Boardman JP, et al. Atlas selection strategy for automatic segmentation of pediatric brain MRIs into 83 ROIs. 2010 IEEE International Conference on Imaging Systems and Techniques, IST 2010—Proceedings. 2010; p. 290–293.
  57. 57. Aljabar P, Heckemann RA, Hammers A, Hajnal JV, Rueckert D. Multi-atlas based segmentation of brain images: Atlas selection and its effect on accuracy. NeuroImage. 2009;46(3):726–738. pmid:19245840
  58. 58. Collett BR, Aylward EH, Berg J, Davidoff C, Norden J, Cunningham ML, et al. Brain volume and shape in infants with deformational plagiocephaly. Child’s Nervous System. 2012;28(7):1083–1090. pmid:22447491