Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2930902

Formats

Article sections

Authors

Related links

Neuroimage. Author manuscript; available in PMC 2011 November 1.

Published in final edited form as:

Published online 2010 June 30. doi: 10.1016/j.neuroimage.2010.06.054

PMCID: PMC2930902

NIHMSID: NIHMS219350

Piotr A. Habas,^{a,}^{b,}^{*} Kio Kim,^{a,}^{b} James M. Corbett-Detig,^{a,}^{b} Francois Rousseau,^{c} Orit A. Glenn,^{b} A. James Barkovich,^{b} and Colin Studholme^{a,}^{b}

The publisher's final edited version of this article is available at Neuroimage

See other articles in PMC that cite the published article.

Modeling and analysis of MR images of the developing human brain is a challenge due to rapid changes in brain morphology and morphometry. We present an approach to the construction of a spatiotemporal atlas of the fetal brain with temporal models of MR intensity, tissue probability and shape changes. This spatiotemporal model is created from a set of reconstructed MR images of fetal subjects with different gestational ages. Groupwise registration of manual segmentations and voxelwise nonlinear modeling allow us to capture the appearance, disappearance and spatial variation of brain structures over time. Applying this model to atlas-based segmentation, we generate age-specific MR templates and tissue probability maps and use them to initialize automatic tissue delineation in new MR images. The choice of model parameters and the final performance are evaluated using clinical MR scans of young fetuses with gestational ages ranging from 20.57 to 24.71 weeks. Experimental results indicate that quadratic temporal models can correctly capture growth-related changes in the fetal brain anatomy and provide improvement in accuracy of atlas-based tissue segmentation.

Magnetic resonance (MR) imaging is an essential tool for the study of early human brain development in utero (Girard et al., 1995; Prayer et al., 2006; Rutherford et al., 2008; Glenn, 2010). Recently developed methods of image reconstruction permit the formation of motion-corrected 3D volumes of the fetal brain from 2D MR scans acquired in utero (Rousseau et al., 2006; Jiang et al., 2007; Kim et al., 2010). Accurate and efficient automatic segmentation of individual tissues (Habas et al., in press) is then the next key step in quantification of fetal brain development. This allows us to extract patterns of normal brain tissue growth in utero and correctly identify cases that may be related to abnormal neurodevelopmental outcomes.

Fetal MRI is usually performed after 20 weeks gestational age (GA) when the main steps of organogenesis are completed. At this stage, the fetal brain consists of a mixture of developing gray matter, developing white matter and other transient structures such as the germinal matrix (Prayer et al., 2006). During embryology and early fetal life, the germinal matrix is a site of production of both neurons and glial cells which then migrate out to their final locations. The volume of the germinal matrix reaches its peak at about 23 weeks GA and decreases subsequently (Battin et al., 1998; Kinoshita et al., 2001).

The key challenge in automated analysis of fetal brain MRI data is rapid changes occurring over very short timescales as illustrated in Fig. 1. These include growth-related changes in the size and shape of the brain, maturation-related changes in MR intensities of developing gray and white matter and the appearance and complete disappearance of the germinal matrix from different brain regions. Therefore, to meaningfully label tissues present within a given MR image of the fetal brain, it is necessary to interpret the underlying anatomy in close relation to its developmental stage.

Axial and coronal views from reconstructed T2w MR images of typical fetal brain anatomies at around 21, 22, 23 and 24 weeks gestational age.

Numerous studies with older anatomies suggested that the use of age-specific atlases can significantly improve the results of automated analysis of brain MRI data. In multiple neonatal segmentation studies (Warfield et al., 2000; Prastawa et al., 2005; Xue et al., 2007; Weisenfeld and Warfield, 2009), specific templates for newborn brains were used to correctly label developing tissues, including myelinated and non-myelinated white matter. Murgasova et al. (2007) showed that the use of age-specific atlases significantly improves the accuracy of the expectation-maximization (EM) segmentation algorithm for brain structures in 1-year-old and 2-year-old children. Yoon et al. (2009) constructed an age-matched intensity template and tissue probability maps from a set of pediatric MR images and demonstrated their beneficial effect on morphometric analysis of brain MR data of 8-year-old children. Similarly, the results of a study with adult brain anatomies (Studholme et al., 2001) indicated the importance of using an age-appropriate template for the group being analyzed.

In our previous work (Habas et al., in press), we focused on an early developmental stage (20.5–22.5 weeks GA) and demonstrated the feasibility of atlas-based segmentation of developing tissues in the human brain from motion-corrected in utero MRI. As the study was based on a relatively homogeneous group of young fetuses, a single probabilistic atlas of tissue distribution was constructed and applied to drive automatic segmentations. Extending this framework to a wider range of gestational ages, however, given the dynamics of changes occurring in the young fetal brain, separate atlases should be constructed and used for different weeks of early brain development. Such an approach is not practical due to potentially insufficient availability of MR scans of subjects with normal brain development at certain gestational ages as well as time and expertise required to perform accurate manual segmentations.

Instead, in this paper, we present an approach to the construction of a complete spatiotemporal atlas of the human fetal brain with temporal parametric models of MR intensity, tissue probability and shape changes. Spatiotemporal registration and modeling have been previously applied for cardiac MRI (Meyer et al., 1996; Heller et al., 2001; Perperidis et al., 2005) and are becoming popular in the domain of brain image analysis (Davis et al., 2007; Pohl et al., 2007). However, building an atlas of brain tissue growth poses significantly different challenges than temporal modeling of soft tissue deformation. Here, we demonstrate that a spatiotemporal atlas of the fetal brain with an appropriate degree of nonlinearity is able to correctly capture the spatial and temporal variation in presence of developing brain tissues as well as their appearance in MR images. We also show that age-specific MR intensity templates and tissue probability maps generated from the spatiotemporal atlas can be applied for improved tissue segmentation in MR images of new fetal subjects.

In order to construct a spatiotemporal atlas of the fetal brain, multiple time point imaging of the same fetal anatomy is not a viable route due to physical and emotional burden for the mother. It is not ethically appropriate to scan pregnant women with short enough time intervals to capture subtle changes in morphology and morphometry of the developing fetal brain. As a result, we have to use a collection of clinical MR images of different fetuses with different gestational ages. This approach, on the other hand, allows our model to capture also between-subject variability which is important for building statistical priors for the analysis of new subject images. To bring a set of fetal anatomies into spatial correspondence, we take advantage of their accurate manual segmentations and perform groupwise registration of tissue probability maps. After registration, we model the effect of the gestational age on probabilities of developing cortical gray matter, developing white matter, the germinal matrix and ventricles at each voxel of the fetal brain. In a similar fashion, we create spatiotemporal models of MR intensity, shape and global scale changes to build a complete four-dimensional atlas of the fetal brain anatomy. From this spatiotemporal atlas, we generate age-specific MR intensity templates and corresponding tissue probability maps and apply them for automatic segmentation of developing tissues in MR images of new fetal subjects.

Fetal imaging is performed in our institution on a 1.5T scanner (GE Healthcare, Milwaukee, WI) without sedation of the mother or the fetus. For each subject, multiple stacks of single-shot fast spin-echo (SSFSE) T2-weighted (T2w) slice images (pixel size 1mm×1mm, slice thickness ≈3mm) are acquired in the approximately axial, sagittal and coronal planes with respect to the fetal brain. The MR sequence parameters (TR=4500 ms, TE=91 ms) were originally designed for clinical scans and cannot be adjusted for image analysis purposes. After acquisition, image stacks are registered using the slice intersection motion correction (SIMC) technique (Kim et al., 2010) to account for spontaneous fetal movement during scanning and reconstructed into 3D volumes with isotropic voxel size of 0.5 mm. Reconstructed volumes are manually segmented into regions of developing cortical gray matter (GM), developing white matter (WM) also including the deep gray nuclei, the germinal matrix (GMAT), ventricles (VENT) and non-brain structures (NB).

In order to accurately relate data from *S* training anatomies with different gestational ages, we need to define meaningful and consistent mapping between subjects with different shapes and tissue morphologies. Due to limited contrast of T2w MR images used in clinical practice, there is a substantial intensity overlap between tissues in the growing fetal brain, especially between developing cortical gray matter and the germinal matrix as shown in Fig. 1. As this may cause incorrect alignment of developing tissues during intensity-based registration, we take advantage of existing manual segmentations of all subjects and perform groupwise registration of tissue maps rather than MR intensity images. Unlike fully data-driven approaches, this allows us to resolve registration ambiguities between different tissue types that appear with very similar intensities on T2w MR images used in the study.

We are still, however, faced with the problem of temporally inconsistent structural features such as the border between the germinal matrix and developing white matter gradually disappearing with brain growth. This can induce artifactually large deformations during nonlinear registration and lead to extreme warping between subjects with and without the germinal matrix in different brain regions. In our approach, we exclude such temporally inconsistent tissue boundaries from the registration model and rely only on three basic classes that exist at all gestational ages, namely

- developing cortical gray matter (GM),
- developing white matter and the germinal matrix combined (WM+GMAT), and
- “others” including ventricles (VENT), extraventricular cerebrospinal fluid (CSF) and non-brain structures (NB) (e.g., skull, amniotic fluid, etc.).

To provide continuous band-limited images which can be interpolated and transformed without introducing aliasing artifacts as well as to reduce the influence of any voxel errors potentially present in the manual tracings, we generate a set of *S* smooth maps

$${\mathbf{m}}_{s}(\mathbf{x})=\left[\begin{array}{cc}{m}_{s}^{(\text{GM})}(\mathbf{x})\hfill & {m}_{s}^{(\text{WM}+\text{GMAT})}(\mathbf{x})\hfill \end{array}\right]$$

(1)

by blurring of binary regions extracted from the full manual segmentation of each subject *s*. The scalar map components, ${m}_{s}^{(\text{GM})}(\mathbf{x})\phantom{\rule{thinmathspace}{0ex}}\text{and}\phantom{\rule{thinmathspace}{0ex}}{m}_{s}^{(\text{WM}+\text{GMAT})}(\mathbf{x})$, range from 0 to 1 and represent the presence of each class at spatial location **x** as shown in Fig. 2.

Axial views from an MR intensity image *i*_{s}(**x**), its full manual segmentation **p**_{s}(**x**) as well as components ${m}_{s}^{(\text{GM})}(\mathbf{x})\phantom{\rule{thinmathspace}{0ex}}\text{and}\phantom{\rule{thinmathspace}{0ex}}{m}_{s}^{(\text{WM}+\text{GMAT})}(\mathbf{x})$ of a smooth tissue map **m**_{s}(**x**) used for spatial normalization.

To initialize groupwise registration, we select a reference anatomy and perform a global linear registration of scalar maps ${m}_{s}^{(\text{WM}+\text{GMAT})}(\mathbf{x})$ to estimate subjects-specific rotations, translations and scaling parameters **a**_{s} mapping from the reference coordinates to each of the *S* training anatomies. As a global transformation, this affine registration can only approximately align the main tissue boundaries of interest even when both components of **m**_{s}(**x**) are used. Starting from this initial alignment, fine scale differences between brain anatomies including changes in cortical folding are then estimated by a more localized deformable template-free groupwise registration (Studholme and Cardenas, 2004), here using full tissue maps **m**_{s}(**x**) rather than MRI intensities. This template-free approach avoids both the use of a single target brain and the construction of an intermediate blurred target image by forming a collective measure of alignment for all subject anatomies.

To define such a collective measure of alignment for multiple tissue classes, we begin by considering the match between a pair of subject anatomies (s,q) from the training set. The squared difference

$${d}_{(\mathrm{s},\mathrm{q})}(\mathbf{x})={\Vert {\mathbf{m}}_{\mathrm{s}}(\mathbf{x}+{\mathbf{u}}_{\mathrm{s}}(\mathbf{x}))-{\mathbf{m}}_{\mathrm{q}}(\mathbf{x}+{\mathbf{u}}_{\mathrm{q}}(\mathbf{x}))\Vert}_{1}^{2}$$

(2)

between the subjects' tissue maps **m**_{s}(**x**) and **m**_{q}(**x**) transformed to the intermediate average space by subject-specific displacement fields **u**_{s}(**x**) and **u**_{q}(**x**), respectively, is used as the alignment measure of the pair. The derivative of *d*_{(s,q)}(**x**) with respect to **u**_{s}(**x**) is used as the registration force. The collective force *d _{s}*(

$${\mathbf{d}}_{\mathrm{s}}(\mathbf{x})=-{\displaystyle \sum _{q\ne s}{{\mathbf{u}}_{s}(\mathbf{x})}_{{d}_{(s,q)}}}$$

(3)

To regularize the displacement fields **u**_{s}(**x**), we use an elastic deformation model (Bajcsy et al., 1983; Shen and Davatzikos, 2002) in which the tissue difference force **d**_{s}(**x**) is balanced by linear elastic energy

$${\mathbf{d}}_{s}(\mathbf{x})=\mu {2}^{{\mathbf{u}}_{s}}$$

(4)

where μ and λ are constants determining the strength of regularization. The value of β determines the influence of a symmetric penalty term derived from the current average displacement

$$\overline{\mathbf{u}}(\mathbf{x})=\frac{1}{S}{\displaystyle \sum _{s=1}^{S}{\mathbf{u}}_{s}(\mathbf{x})}$$

(5)

which drives the average deformation of all subjects toward zero. An iterative gradient descent is used to solve Eq. 4 over all subjects. This seeks a set of *S* displacement fields **u**_{s}(**x**) that balance the elastic deformation energy with the image dissimilarity while seeking an average anatomy equidistant from each of the *S* subjects. With the parameters λ and μ experimentally set to suitable values, the resulting set of *S* subject-specific displacement fields **u**_{s}(**x**) provides an invertible mapping between the group's average shape and each of the *S* training anatomies.

After registration, for each subject *s* with gestational age *t _{s}*, we use its global linear scaling

The growth-related changes in the size and shape of the fetal brain are captured firstly by the three directional components of the initial affine scaling ${a}_{s}=\left[\begin{array}{ccc}{a}_{s}^{(1)}\hfill & {a}_{s}^{(2)}\hfill & {a}_{s}^{(3)}\hfill \end{array}\right]$ and then by the displacement fields **u**_{s}(**x**) estimated during the elastic groupwise registration.

For each linear scaling component, a parametric temporal model *a*^{(j)}(*t*) is obtained by polynomial fitting to values ${a}_{s}^{(j)}$ corresponding to the gestational ages *t _{s}* of the

$${a}^{(j)}(t)={\displaystyle \sum _{d=0}^{{D}_{a}}{\mathrm{(\mathit{\text{jd}})}}^{{t}^{d}}}$$

(6)

with coefficients ^{(jd)} estimated through least squares fitting to the *S* spatial scaling factors ${a}_{s}^{(j)}$

$${\mathbf{D}}_{a}=[\begin{array}{ccc}\hfill {t}_{1}^{0}\hfill & \hfill {t}_{1}^{1}\hfill & \hfill \hfill {t}_{1}^{{D}_{a}}\hfill \hfill & \hfill {t}_{2}^{0}\hfill & \hfill {t}_{2}^{1}\hfill & \hfill \hfill {t}_{2}^{{D}_{a}}\hfill \hfill & \hfill \hfill \hfill \hfill \hfill & \hfill {t}_{S}^{0}\hfill & \hfill {t}_{S}^{1}\hfill & \hfill \hfill {t}_{S}^{{D}_{a}}\hfill \hfill \hfill & ]\hfill \hfill \end{array}$$

(7)

$$[\begin{array}{c}\hfill {\mathrm{(j0)}}^{}\hfill & \hfill {\mathrm{(j1)}}^{}\hfill & \hfill \hfill & \hfill {\mathrm{({\mathit{\text{jD}}}_{a})}}^{}\hfill \\ ]\\ ={({\mathbf{D}}_{a}^{T}{\mathbf{D}}_{a})}^{-1}{\mathbf{D}}_{a}^{T}\phantom{\rule{thinmathspace}{0ex}}[\begin{array}{c}\hfill {a}_{1}^{(j)}\hfill \\ \hfill {a}_{2}^{(j)}\hfill \\ \hfill \hfill & \hfill {a}_{S}^{(j)}\hfill \end{array}]\end{array}$$

(8)

Temporal modeling of shape deformations is performed separately for each component of the three-dimensional displacement vector **u**(**x**). For each deformation direction *j*, changes in *u*^{(j)}(**x**,*t*) are modeled using a temporal polynomial of degree *D _{u}*

$${u}^{(j)}(\mathbf{x},t)={\displaystyle \sum _{d=0}^{{D}_{u}}{\mathrm{(\mathit{\text{jd}})}}^{(}}$$

(9)

with coefficients ^{(jd)}(**x**) found through least squares fitting to the *S* values ${u}_{s}^{(j)}(\mathbf{x})$ extracted from the age-specific displacement fields **u**_{s}(**x**).

$${\mathbf{D}}_{u}=[\begin{array}{ccc}\hfill {t}_{1}^{0}\hfill & \hfill {t}_{1}^{1}\hfill & \hfill \hfill {t}_{1}^{{D}_{u}}\hfill \hfill & \hfill {t}_{2}^{0}\hfill & \hfill {t}_{2}^{1}\hfill & \hfill \hfill {t}_{2}^{{D}_{u}}\hfill \hfill & \hfill \hfill \hfill \hfill \hfill & \hfill {t}_{S}^{0}\hfill & \hfill {t}_{S}^{1}\hfill & \hfill \hfill {t}_{S}^{{D}_{u}}\hfill \hfill \hfill & ]\hfill \hfill \end{array}$$

(10)

$$[\begin{array}{c}\hfill {\mathrm{(j0)}}^{(}\hfill & \hfill {\mathrm{(j1)}}^{(}\hfill & \hfill \hfill & \hfill {\mathrm{({\mathit{\text{jD}}}_{u})}}^{(}\hfill \\ ]\\ ={({\mathbf{D}}_{u}^{T}{\mathbf{D}}_{u})}^{-1}{\mathbf{D}}_{u}^{T}\phantom{\rule{thinmathspace}{0ex}}[\begin{array}{c}\hfill {u}_{1}^{(j)}(\mathbf{x})\hfill \\ \hfill {u}_{2}^{(j)}(\mathbf{x})\hfill \\ \hfill \hfill & \hfill {u}_{S}^{(j)}(\mathbf{x})\hfill \end{array}]\end{array}$$

(11)

The MRI data of each individual are transformed to the group coordinate system allowing the evaluation of age-related intensity changes across the atlas population. The relative contrasts of different tissues are known to vary spatially across the brain even within a given tissue class and over time as the fetal brain develops. As a result, we have to capture not only the average intensity of each tissue class but also the temporal changes at each voxel location. MR intensity changes at voxel **x** are modeled using a temporal polynomial of degree *D _{i}*

$$i(\mathbf{x},t)={\displaystyle \sum _{d=0}^{{D}_{i}}{\mathrm{(d)}}^{}(\mathbf{x}){t}^{d}}$$

(12)

with coefficients ^{(d)}(**x**) found through least squares fitting to the *S* observed intensity values *i _{s}*(

$${\mathbf{D}}_{i}=[\begin{array}{ccc}\hfill {t}_{1}^{0}\hfill & \hfill {t}_{1}^{1}\hfill & \hfill \hfill {t}_{1}^{{D}_{i}}\hfill \hfill & \hfill {t}_{2}^{0}\hfill & \hfill {t}_{2}^{1}\hfill & \hfill \hfill {t}_{2}^{{D}_{i}}\hfill \hfill & \hfill \hfill \hfill \hfill \hfill & \hfill {t}_{S}^{0}\hfill & \hfill {t}_{S}^{1}\hfill & \hfill \hfill {t}_{S}^{{D}_{i}}\hfill \hfill \hfill & ]\hfill \hfill \end{array}$$

(13)

$$[\begin{array}{c}\hfill {\mathrm{(0)}}^{(}\hfill & \hfill {\mathrm{(1)}}^{(}\hfill & \hfill \hfill & \hfill {\mathrm{({D}_{i})}}^{(}\hfill \\ ]\\ ={({\mathbf{D}}_{i}^{T}{\mathbf{D}}_{i})}^{-1}{\mathbf{D}}_{i}^{T}\phantom{\rule{thinmathspace}{0ex}}[\begin{array}{c}\hfill {i}_{1}(\mathbf{x})\hfill \\ \hfill {i}_{2}(\mathbf{x})\hfill \\ \hfill \hfill & \hfill {i}_{S}(\mathbf{x})\hfill \end{array}]\end{array}$$

(14)

The final component of the spatiotemporal atlas is provided by spatially normalized manual tissue segmentations of each subject *s* with gestational age *t _{s}*. From these discrete segmentations ${p}_{s}(x)=[{p}_{s}^{(\text{GM})}(\mathbf{x})\phantom{\rule{thinmathspace}{0ex}}{p}_{s}^{(\text{WM})}(\mathbf{x})\phantom{\rule{thinmathspace}{0ex}}{p}_{s}^{(\text{GMAT})}(\mathbf{x})\phantom{\rule{thinmathspace}{0ex}}{p}_{s}^{(\text{VENT})}\phantom{\rule{thinmathspace}{0ex}}(\mathbf{x})\phantom{\rule{thinmathspace}{0ex}}{p}_{s}^{(\text{NB})}(\mathbf{x})]$, we generate a set of continuous models

Fitting of temporal polynomials *p*(^{1)}(**x**,*t*) (green line) and *p*^{(2)}(**x**,*t*) (red line) to subject-specific probabilities ${p}_{s}^{(k)}(\mathbf{x})$ of two abstract mutually exclusive classes (green and red dots, respectively). (A) Unconstrained modeling directly in the probability **...**

In the multinomial LogOdds representation suitable for medical imaging problems with multiple labels (Pohl et al., 2007), the LogOdds (·) of probability ${p}_{s}^{(k)}(\mathbf{x})$ is the logarithm of the ratio between ${p}_{s}^{(k)}(\mathbf{x})$ and the last entry of **p**_{s}(**x**)

$$\mathrm{({p}_{s}^{(k)}(\mathbf{x}))=\mathit{\text{log}}\left(\frac{{p}_{s}^{(k)}(\mathbf{x})}{{p}_{s}^{(\text{NB})}(\mathbf{x})}\right)={l}_{s}^{(k)}(\mathbf{x})}$$

(15)

where *k*{*GM, WM, GMAT, VENT*}. In the context of fetal brain tissue segmentation, ${p}_{s}^{(k)}(\mathbf{x})$ is the probability that voxel **x** in manual segmentation of subject anatomy *s* is assigned to tissue class *k* whereas ${p}_{s}^{(\text{NB})}(\mathbf{x})$ is the probability of assignment of voxel **x** to non-brain structures. The inverse LogOdds function ^{−1}(·) is the generalized logistic function

$${\mathrm{-1}}^{(}$$

(16)

where

$$Z=1+{\displaystyle \sum _{k\ne \mathit{\text{NB}}}\mathit{\text{exp}}({l}_{s}^{(k)}(\mathbf{x}))}$$

(17)

is a normalizing factor.

After transformation of the age-dependent tissue probabilities **p**_{s}(**x**) to the space of LogOdds, the LogOdds of class *k*{GM, WM, GMAT, VENT} are modeled using a temporal polynomial of degree *D _{p}*

$${l}^{(k)}(\mathbf{x},t)={\displaystyle \sum _{d=0}^{{D}_{p}}{\mathrm{(\mathit{\text{kd}})}}^{(}}$$

(18)

with coefficients ^{(kd)}(**x**) found through least squares fitting to the *S* observed LogOdds values ${l}_{s}^{(k)}(\mathbf{x})$.

$${\mathbf{D}}_{p}=[\begin{array}{ccc}\hfill {t}_{1}^{0}\hfill & \hfill {t}_{1}^{1}\hfill & \hfill \hfill {t}_{1}^{{D}_{p}}\hfill \hfill & \hfill {t}_{2}^{0}\hfill & \hfill {t}_{2}^{1}\hfill & \hfill \hfill {t}_{2}^{{D}_{p}}\hfill \hfill & \hfill \hfill \hfill \hfill \hfill & \hfill {t}_{S}^{0}\hfill & \hfill {t}_{S}^{1}\hfill & \hfill \hfill {t}_{S}^{{D}_{p}}\hfill \hfill \hfill & ]\hfill \hfill \end{array}$$

(19)

$$[\begin{array}{c}\hfill {\mathrm{(k0)}}^{(}\hfill & \hfill {\mathrm{(k1)}}^{(}\hfill & \hfill \hfill & \hfill {\mathrm{({\mathit{\text{kD}}}_{p})}}^{(}\hfill \\ ]\\ ={({\mathbf{D}}_{p}^{T}{\mathbf{D}}_{p})}^{-1}{\mathbf{D}}_{p}^{T}\phantom{\rule{thinmathspace}{0ex}}[\begin{array}{c}\hfill {l}_{1}^{(k)}(\mathbf{x})\hfill \\ \hfill {l}_{2}^{(k)}(\mathbf{x})\hfill \\ \hfill \hfill & \hfill {l}_{S}^{(k)}(\mathbf{x})\hfill \end{array}]\end{array}$$

(20)

A complete age-specific fetal brain anatomycan be synthesized from the estimated model parameters , , and by discretizing the continuous temporal models of global linear scaling, local deformation, MR intensity and tissue probability, respectively, for a given gestational age *t* (Fig. 5). A synthetic age-specific MR image *i*(**x**,*t*) is composed in the space of the atlas from intensity estimates for each voxel **x** at gestational age *t* (Fig. 5A). To obtain a corresponding vector of expected tissue probabilities **p**(**x**,*t*), the continuous temporal models *l*^{(k)}(**x**,*t*) are discretized in the space of LogOdds and the resulting values are transferred to the probability space using the inverse LogOdds function *O*^{−1}(·).

Synthesis of a complete fetal brain anatomy with application to segmentation of new fetal MRI. (A) An age-specific MR intensity template and tissue probability maps are synthesized from model parameters in the space of the spatiotemporal atlas. (B) Age-specific **...**

The inverse of the age-specific shape deformation **u**(**x**,*t*) and global linear scaling **a**(*t*) calculated from respective temporal models need to be applied to both the synthesized MR image *i*(**x**,*t*) and the tissue probability map **p**(**x**,*t*) to transform them to the expected size and shape of the fetal brain at the gestational age *t* (Fig. 5B).

Age-specific MR intensity templates and tissue probability maps generated from the spatiotemporal atlas can be used for tissue segmentation of new subject anatomies as also shown in Fig. 5. First, the MR image of a new subject to be segmented needs to be registered to the age-matched synthetic MR intensity template (Fig. 5C). This is achieved by a sequence of global linear transformations followed by multiple elastic deformations. Then, based on the inverse of the estimated spatial transformation, the age-specific tissue probability map is warped to the space of the new subject (Fig. 5D) and used as a source of spatially varying priors for atlas-based tissue segmentation of the new fetal anatomy.

The study was performed using clinical MR scans of 20 fetal subjects with normal brain development selected for reconstruction and manual tracing from a larger pool collected at University of California San Francisco. The mothers were referred for fetal MRI due to questionable abnormalities on prenatal ultrasound, a prior abnormal pregnancy or recruited as volunteers. The gestational ages *t _{s}* of the fetuses were almost uniformly distributed between 20.57 and 24.71 weeks as shown in Fig. 6.

For each subject, a high-resolution MR volume was formed using the slice intersection motion correction (SIMC) technique (Kim et al., 2010). The reconstructed MR volumes were manually segmented into regions of developing cortical gray matter (GM), developing white matter (WM), the germinal matrix (GMAT), ventricles (VENT) and non-brain structures (NB). The resulting tissue label maps served as reference segmentations for further modeling and validation. The reconstructed 3D images were also corrected for intensity in homogeneity commonly present in fetal MRI data. Low-degree polynomial models of the bias field were estimated using the EM algorithm (Van Leemput et al., 1999) with fixed binary tissue probabilities given by the reference manual segmentations. Finally, the intensities of all images were linearly adjusted based on the average intensity of ventricular CSF.

To bring all 20 fetal anatomies into spatial correspondence, tissue maps **m**_{s}(**x**) were extracted from full manual segmentation of each subject *s* and smoothed with a 2-mm full width at half maximum (FWHM) Gaussian kernel. The ${m}_{s}^{(\text{WM}+\text{GMAT})}(\mathbf{x})$ component of each tissue map was first globally aligned to an average shape anatomy from a previous segmentation study (Habas et al., in press) by maximization of normalized mutual information (Studholme et al., 1999). The obtained 20 age-specific estimates of global linear scaling **a**_{s} in three orthogonal directions (left–right, anterior–posterior and superior–inferior) were used for temporal modeling of changes in the overall size of the fetal brain.

This initial linear registration also provided a starting point for further collective alignment of the full tissue maps **m**_{s}(**x**) into a common average space. To resolve large shape differences between fetal brain anatomies with different gestational ages and keep the mapping diffeomorphic, we used a composition of two elastic deformations. For each of them, we set λ=0 as suggested by Modersitzki (2004) and experimentally chose the value of μ=15 to ensure that folding of the deformations did not occur during registration. The value of β=200 was also selected empirically and the final transformations were not found to be very sensitive to this parameter. Typical registration results for a representative subset of 8 subjects with different gestational ages are presented in Fig. 7.

Tissue maps ${m}_{s}^{(\text{WM}+\text{GMAT})}(\mathbf{x})$ of 8 fetal subjects after rigid registration (top row), global linear registration (middle row) and groupwise multiple elastic registration (bottom row).

After groupwise registration, the 20 age-specific displacement fields **u**_{s}(**x**) were directly used for spatiotemporal modeling of shape changes. Based on **a**_{s} and **u**_{s}(**x**), the full manual segmentations of each subject *s* were transformed to the average space using nearest-neighbor interpolation and provided 20 data points **p**_{s}(**x**) for spatiotemporal modeling for tissue probabilities. The MR intensity images of all subjects were warped to the average space using linear interpolation and provided 20 data points *i _{s}*(

Temporal modeling of global size changes was performed independently for each of the three orthogonal directions as shown in Fig. 8. The ability of polynomial models with different degrees *D _{a}* to capture temporal patterns present in the training data was evaluated by the coefficient of determination,

Temporal modeling of global linear scaling in three orthogonal directions using linear models (*D*_{a}=1, dashed lines) and quadratic models (*D*_{a}=2, solid lines).

For evaluation of temporal models of shape changes, spatial maps of the *R*^{2} coefficient were generated from polynomial shape models *u*(**x**,*t*) with degrees *D _{u}*=1,

Fig. 10 shows spatial maps of the *R*^{2} coefficient obtained from temporal models of intensity changes *i*(**x**,*t*) with degrees *D _{i}*=1,

Smoothed maps of *R*^{2} coefficient for temporal models of intensity changes with *D*_{i}=1, *D*_{i}=2, and *D*_{i}=3. Cooler colors represent lower *R*^{2} values while warmer colors represent higher values of the *R*^{2} coefficient. Maps were automatically masked to the brain **...**

For evaluation of tissue probability models, we performed qualitative assessment of age-specific tissue probability maps generated from temporal models with different degrees *D _{p}* and quantitative analysis of their impact on atlas-based tissue segmentation of new fetal anatomies. Initially, we created a spatiotemporal atlas with the degree of tissue probability model

Average tissue probability maps *p*^{(k)}(**x**) for developing gray matter (GM), developing white matter (WM), the germinal matrix (GMAT) and ventricles (VENT) generated in the average space from a spatiotemporal atlas with *D*_{p}=0 (no temporal dependency).

Then, to model temporal changes in tissue presence at different regions of the fetal brain, we also created spatiotemporal atlases with *D _{p}*=1 and

Age-specific tissue probability maps *p*^{(k)} (**x**,*t*) for developing white matter (WM) and the germinal matrix (GMAT) generated in the average space from a spatiotemporal atlas with *D*_{p}=2.

Based on previously reported results for neonates and young children (Murgasova et al., 2007; Yoon et al., 2009), we expected that the use of age-specific atlases may improve the results of atlas-based tissue segmentation of the fetal brain. To verify this hypothesis, we performed the following leave-one-out validation experiment.

After groupwise registration of the entire study population (20 subjects), each subject was held out once for testing purposes. The remaining 19 subjects were used to build a set of spatiotemporal atlases with *D _{a}*=2,

Table 1 presents summary results of the segmentation experiment. For all degrees of the tissue probability model, the accuracy of automatic atlas-based segmentation for developing cortical gray matter, developing white matter and ventricles (DSC>0.8) is similar to DSC values previously reported in segmentation studies with neonates (Prastawa et al., 2005; Xue et al., 2007) and young children (Murgasova et al., 2007). The inferior average accuracy for segmentation of the germinal matrix with a single atlas (*D _{p}*=0, DSC=0.675) arises mainly from oversegmenting of this structure in older subjects. Fig. 15 shows results of atlas-based tissue segmentation for a fetal subject at 24 weeks GA. For

The impact of the degree of the tissue probability model (*D*_{p}) on automatic atlas-based tissue segmentation of a fetal subject at 24 weeks GA. Arrows indicate areas where the use of age-specific atlases generated from models with *D*_{p}=1 and *D*_{p}=2 improves **...**

Average DSC accuracy of atlas-based EM segmentation with probabilistic atlases generated from spatiotemporal models with *D*_{p}=0, *D*_{p}=1 and *D*_{p}=2 (20 subjects).

The use of age-specific atlases generated from tissue probability models with *D _{p}*=1 and

We presented an approach to modeling of early human brain growth from clinical MRI. Using a set of reconstructed 3D MR images and their manually defined tissue segmentations, we created a spatiotemporal atlas of MR intensity, tissue probability and shape changes in the human fetal brain during early in utero development. From this continuous mathematical model, we can synthesize maps of MR intensity and tissue occurrence for any gestational age covered by the atlas and apply them for improved tissue segmentation in new fetal anatomies.

As multiple time point imaging of the same fetal anatomy is not feasible, the spatiotemporal atlas was built using MR images of different fetuses with different gestational ages. To bring all anatomies into collective alignment in the average space of the group, we applied groupwise registration of tissue maps extracted from fixed manual segmentations. Although this approach did not define full correspondence between images with different gestational ages, it allowed us to consistently align the main tissue boundaries present in the developing fetal brain and resolve matching ambiguities between multiple tissue types that appear with similar intensities on T2w MR images. Extending this framework to a wider range of gestational ages, we plan to replace the single average space with temporarily parameterized average geometry. Another possible direction of future work is to combine the groupwise registration with simultaneous tissue segmentation and atlas building (Ashburner and Friston, 2005; Pohl et al., 2006a).

After registration, we performed parametric modeling of changes in MR intensity, tissue probability and shape of the developing fetal brain. Although additional features such as gender or disease group could be also considered, we used gestational age as the only independent variable for this initial work. For all analyzed properties of the fetal brain anatomy, quadratic temporal models provided a better description of patterns observed for the study population than linear models. For modeling of tissue probabilities, quadratic polynomials allowed us to represent both appearance and disappearance of the germinal matrix from different regions of the brain. As tissue occurrence or MR intensity are not expected to oscillate, the use of higher order temporal models was neither justified nor necessary.

Our morphology model included three main tissue types present in the young fetal brain (developing cortical gray matter, developing white matter, the germinal matrix) and ventricles. Other structures, such as the basal ganglia or the subplate, are difficult to see at this early stage (20–24 weeks GA) but become more visible as the fetuses mature. As a temporary solution, we combined these regions of the fetal brain with developing white matter. In the future, however, they will be fully incorporated into the modeling process as separate and important structures.

Spatiotemporal modeling of all properties was performed in a voxelwise manner, without explicit constraints on spatial smoothness of synthesized maps. Although not an issue for MR intensity and tissue probabilities, this may potentially cause problems in modeling of shape changes. In practice, however, given the spatial smoothness of all subject maps, we found that the shape model provided invertible diffeomorphic mappings for all time points. In the future, we plan to compare our voxelwise approach to alternative techniques such as kernel regression for large deformation models (Davis et al., 2007) and explore extensions to global modeling using ridge regression (Sjostrand et al., 2008).

Finally, we presented a practical application of the spatiotemporal atlas for tissue segmentation in new MR images. As the synthesis of age-matched MR templates reduces the geometric differences that have to be resolved during registration, the anatomy to be segmented can be precisely aligned with the source of tissue priors. The use of age-specific probabilistic atlases improves the overall accuracy of automatic tissue segmentation, as measured by the average DSC values, and allows for correct delineation of transient brain structures that evolve rapidly with gestational age. In clinical applications, accurate quantification of the fetal brain anatomy in relation to its developmental stage may help detect subtle abnormalities that otherwise could be missed.

This work was primarily funded by NIH/NINDS Grant R01 NS 055064. Imaging studies were partially funded by NIH/NINDS Grant K23 NS 052506 and NIH/NCRR UCSF-CTSI Grant UL1 RR 024131. The work of F. Rousseau was supported by the European Research Council under the European Community Seventh Framework Programme (FP7/2007-2013 Grant Agreement 207667).

The authors thank Dr. Julia A. Scott for help with editing of the manuscript.

- Ashburner J, Friston KJ. Unified segmentation. Neuroimage. 2005;26(3):839–851. [PubMed]
- Bajcsy R, Lieberson R, Reivich M. A computerized system for the elastic matching of deformed radiographic images to idealized atlas images. J. Comput. Assist. Tomogr. 1983;7(4):618–625. [PubMed]
- Battin MR, Maalouf EF, Counsell SJ, Herlihy AH, Rutherford MA, Azzopardi D, Edwards AD. Magnetic resonance imaging of the brain in very preterm infants: visualization of the germinal matrix, early myelination, and cortical folding. Pediatrics. 1998;101(6):957–962. [PubMed]
- Davis BC, Fletcher PT, Bullitt E, Joshi SC. Population shape regression from random design data; Proc. International Conference on Computer Vision; 2007. pp. 1–7.
- Dice LR. Measures of the amount of ecologic association between species. Ecology. 1945;26(3):297–302.
- Girard N, Raybaud C, Poncet M. In vivo MR study of brain maturation in normal fetuses. Am. J. Neuroradiol. 1995;16(2):407–413. [PubMed]
- Glenn OA. MR imaging of the fetal brain. Pediatr. Radiol. 2010;40(1):68–81. [PMC free article] [PubMed]
- Habas PA, Kim K, Rousseau F, Glenn OA, Barkovich AJ, Studholme C. Atlas-based segmentation of developing tissues in the human brain with quantitative validation in young fetuses. Hum. Brain Mapp. in press. doi:10.1002/hbm.20935. [PMC free article] [PubMed]
- Heller EN, Staib LH, Dione DP, Constable RT, Shi CQ, Duncan JS, Sinusas AJ. A new method for quantification of spatial and temporal parameters of endocardial motion: evaluation of experimental infarction using magnetic resonance imaging. Can. J. Cardiol. 2001;17(3):309–318. [PubMed]
- Jiang S, Xue H, Glover A, Rutherford M, Rueckert D, Hajnal JV. MRI of moving subjects using multislice snapshot images with volume reconstruction (SVR): application to fetal, neonatal, and adult brain studies. IEEE Trans. Med. Imaging. 2007;26(7):967–980. [PubMed]
- Kim K, Habas PA, Rousseau F, Glenn OA, Barkovich AJ, Studholme C. Intersection based motion correction of multislice MRI for 3-D in utero fetal brain image formation. IEEE Trans. Med. Imaging. 2010;29(1):146–158. [PMC free article] [PubMed]
- Kinoshita Y, Okudera T, Tsuru E, Yokota A. Volumetric analysis of the germinal matrix and lateral ventricles performed using MR images of postmortem fetuses. Am. J. Neuroradiol. 2001;22(2):382–388. [PubMed]
- Meyer FG, Constable RT, Sinusas AJ, Duncan JS. Tracking myocardial deformation using phase contrast MR velocity fields: a stochastic approach. IEEE Trans. Med. Imaging. 1996;15(4):453–465. [PubMed]
- Modersitzki J. Numerical methods for image registration. Oxford Science Publications; 2004.
- Murgasova M, Dyet L, Edwards AD, Rutherford M, Hajnal JV, Rueckert D. Segmentation of brain MRI in young children. Acad. Radiol. 2007;14(11):1350–1366. [PubMed]
- Perperidis D, Mohiaddin RH, Rueckert D. Spatio-temporal free-form registration of cardiac MR image sequences. Med. Image Anal. 2005;9(5):441–456. [PubMed]
- Pohl KM, Fisher J, Grimson WEL, Kikinis R, Wells WM. A Bayesian model for joint segmentation and registration. Neuroimage. 2006a;31(1):228–239. [PubMed]
- Pohl KM, Fisher J, Shenton M, McCarley RW, Grimson WEL, Kikinis R, Wells WM. Logarithm odds maps for shape representation. Medical Image Computing and Computer-Assisted Intervention: LNCS. 2006b;4191:955–963. [PMC free article] [PubMed]
- Pohl KM, Fisher J, Bouix S, Shenton M, McCarley RW, Grimson WEL, Kikinis R, Wells WM. Using the logarithm of odds to define a vector space on probabilistic atlases. Med. Image Anal. 2007;11(5):465–477. [PMC free article] [PubMed]
- Prastawa M, Gilmore JH, Lin W, Gerig G. Automatic segmentation of MR images of the developing newborn brain. Med. Image Anal. 2005;9(5):457–466. [PubMed]
- Prayer D, Kasprian G, Krampl E, Ulm B, Witzani L, Prayer L, Brugger PC. MRI of normal fetal brain development. Eur. J. Radiol. 2006;57(2):199–216. [PubMed]
- Rousseau F, Glenn OA, Iordanova B, Rodriguez-Carranza CE, Vigneron DB, Barkovich AJ, Studholme C. Registration-based approach for reconstruction of high-resolution in utero fetal MR brain images. Acad. Radiol. 2006;13(9):1072–1081. [PubMed]
- Rutherford M, Jiang S, Allsop J, Perkins L, Srinivasan L, Hayat T, Kumar S, Hajnal J. MR imaging methods for assessing fetal brain development. Dev. Neurobiol. 2008;68(6):700–711. [PubMed]
- Shen D, Davatzikos C. HAMMER: hierarchical attribute matching mechanism for elastic registration. IEEE Trans. Med. Imaging. 2002;21(11):1421–1439. [PubMed]
- Sjostrand K, Cardenas VA, Larsen R, Studholme C. A generalization of voxel-wise procedures for high-dimensional statistical inference using ridge regression. Medical Imaging 2008; Image Processing: Proc. SPIE; 2008. p. 69140A.
- Studholme C, Cardenas VA. A template free approach to volumetric spatial normalization of brain anatomy. Pattern Recognit. Lett. 2004;25(10):1191–1202.
- Studholme C, Hill DLG, Hawkes DJ. An overlap invariant entropy measure of 3D medical image alignment. Pattern Recognit. 1999;32(1):71–86.
- Studholme C, Cardenas VA, Weiner MW. Multiscale image and multiscale deformation of brain anatomy for building average brain atlases. Medical Imaging 2001; Image Processing: Proc. SPIE; 2001. pp. 557–568.
- Van Leemput K, Maes F, Vandermeulen D, Suetens P. Automated model-based bias field correction of MR images of the brain. IEEE Trans. Med. Imaging. 1999;18(10):885–896. [PubMed]
- Viola P, Wells WM. Alignment by maximization of mutual information. Int. J. Comput. Vision. 1997;24(2):137–154.
- Warfield S, Kaus M, Jolesz FA, Kikinis R. Adaptive, template moderated, spatially varying statistical classification. Med. Image Anal. 2000;4(1):43–55. [PubMed]
- Weisenfeld NI, Warfield SK. Automatic segmentation of newborn brain MRI. Neuroimage. 2009;47(2):564–572. [PMC free article] [PubMed]
- Xue H, Srinivasan L, Jiang S, Rutherford M, Edwards AD, Rueckert D, Hajnal JV. Automatic segmentation and reconstruction of the cortex from neonatal MRI. Neuroimage. 2007;38(3):461–477. [PubMed]
- Yoon U, Fonov VS, Perusse D, Evans AC. The effect of template choice on morphometric analysis of pediatric brain data. Neuroimage. 2009;45(3):769–777. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |