Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 18 October 2021
Sec. Biomechanics
This article is part of the Research Topic Advances in Brain Mechanics View all 12 articles

Subject-Specific Head Model Generation by Mesh Morphing: A Personalization Framework and Its Applications

  • Division of Neuronic Engineering, Department of Biomedical Engineering and Health Systems, KTH Royal Institute of Technology, Stockholm, Sweden

Finite element (FE) head models have become powerful tools in many fields within neuroscience, especially for studying the biomechanics of traumatic brain injury (TBI). Subject-specific head models accounting for geometric variations among subjects are needed for more reliable predictions. However, the generation of such models suitable for studying TBIs remains a significant challenge and has been a bottleneck hindering personalized simulations. This study presents a personalization framework for generating subject-specific models across the lifespan and for pathological brains with significant anatomical changes by morphing a baseline model. The framework consists of hierarchical multiple feature and multimodality imaging registrations, mesh morphing, and mesh grouping, which is shown to be efficient with a heterogeneous dataset including a newborn, 1-year-old (1Y), 2Y, adult, 92Y, and a hydrocephalus brain. The generated models of the six subjects show competitive personalization accuracy, demonstrating the capacity of the framework for generating subject-specific models with significant anatomical differences. The family of the generated head models allows studying age-dependent and groupwise brain injury mechanisms. The framework for efficient generation of subject-specific FE head models helps to facilitate personalized simulations in many fields of neuroscience.

Introduction

Finite element (FE) head models have become powerful tools to simulate brain stimulations with direct current (tDCS) (Datta et al., 2009; Datta et al., 2012; Huang et al., 2013; Windhoff et al., 2013; Opitz et al., 2015; Alekseichuk et al., 2019; Li et al., 2020; Wang et al., 2020), magnetic (TMS) (Opitz et al., 2013), and ultrasound (TUS) (Legon et al., 2014). Such models are also being used to study the development of neurodegenerative diseases (Fornari et al., 2019; Noël and Kuhl, 2019; Weickenmeier et al., 2019) and biomechanical consequences of neurosurgery (Weickenmeier et al., 2017; von Holst and Li, 2014; Li et al., 2015; Ji et al., 2009; Hu et al., 2007; Miller et al., 2010). In particular, FE head models have been tremendously used to study traumatic brain injuries (TBIs) in the last decades (see reviews (Giudice et al., 2019; Horstemeyer et al., 2019; Madhukar and Ostoja-Starzewski, 2019)). Meshing is a first step in generating FE models by discretizing a continuous domain into a finite number of elements, e.g., tetrahedral or hexahedral elements. Generation of FE head models is often time-consuming and challenging due to its complex geometry, though tetrahedral elements are relatively easier to generate, e.g., automatic pipelines have been reported to efficiently generate tetrahedral head models for simulating brain stimulations (Huang et al., 2013; Windhoff et al., 2013). Such efficiency is partially attributed to the well-developed automatic tetrahedral meshing algorithms within mathematics and computer science (Baker, 2005); it is also because the involved partial differential equations (PDEs) are less computationally demanding, permitting a huge number of tetrahedral elements (up to >10 million (Datta et al., 2009)) to capture anatomical details of the brain. Thus, personalized simulations with anatomically detailed subject-specific head models are largely facilitated in these brain stimulation fields.

In contrast to tetrahedrons, hexahedral elements are much more challenging to generate (Baker, 2005; Shepherd and Johnson, 2008) but are preferred in FE head models intended for studying TBIs (hereafter called head injury models) due to their higher efficiency for simulating the incompressible brain under impact (Li et al., 2021). Furthermore, the involved PDEs in head impacts consist of geometrical, material nonlinearity, and complex contacting algorithms, which are more computationally demanding. As a result, many state-of-the-art head injury models (Zhang et al., 2001; Horgan and Gilchrist, 2003; Kleiven, 2007; Takhounts et al., 2008; Mao et al., 2013a; Sahoo et al., 2014; Ji et al., 2015a; Atsumi et al., 2016) use hexahedrons despite the meshing challenges and have simplified brains to reduce the number of elements for computational efficiency. For example, these models have smoothed out brain surfaces and do not have sulci, and gyri, resulting in fewer elements (often <1 million). While a simplified representation of the brain is a reasonable trade-off for computational efficiency, it’s also partially due to the challenges for meshing techniques (e.g., the blocking technique (Mao et al., 2013b)) to capture the anatomical details. It is worth mentioning that the voxel approach is efficient in generating hexahedrons by converting image voxels to hexahedral elements, either directly or with smoothing algorithms. However, a known concern is a less accurate peak strain/stress predicted from such models, especially on the surfaces due to jaggedness. Nevertheless, careful choice of sufficiently refined mesh and result analysis allow such models to provide valuable insights due to their anatomical accuracy (see discussion in (Li et al., 2021)). Besides the much less developed automatic algorithms for generating hexahedrons (Baker, 2005; Shepherd and Johnson, 2008), a necessity to include falx and tentorium to account for their important structural influence on brain mechanical responses during impact (Ho et al., 2017) poses an additional challenge for subject-specific head injury model generation while both structures are often neglected in head models for simulating tDCS, TMS, and TUS. A detailed analysis of the current meshing challenge for head injury models is found in a previous study (Li et al., 2021).

Therefore, the generation of FE head injury models with anatomical details remains a challenge and has become a bottleneck hindering personalized simulations. FE head models without anatomical details such as sulci and gyri also hinder studying detailed mechanisms at areas of interest, such as chronic traumatic encephalopathy (CTE) with pathologies observed at sulcal depth (McKee et al., 2015). Studies have also shown that the brain size/shape influences brain mechanical responses significantly under impact (Kleiven and von Holst, 2002; Li et al., 2021), suggesting the importance of using personalized models to study the onset of TBI in real life. Along with the many existing adult healthy FE head models, there are only a few elderly, children, and infant models (e.g., (Li et al., 2011; Giordano et al., 2017; Li et al., 2017; Li and Kleiven, 2018; Hajiaghamemar et al., 2019; Li et al., 2019; Zhou et al., 2019; Zhou et al., 2020)). TBIs are influencing all age groups, especially infants and the elderly are overrepresented (Pedersen et al., 2015). Thus, it is imperative to investigate efficient approaches for generating detailed subject-specific head injury models across the lifespan and for pathological brains to understand the injury mechanisms and develop preventions.

This study addresses the challenge of generating subject-specific head injury models with hexahedrons, especially concerns about mesh morphing, which is an efficient approach for generating subject-specific models. The approach has been used in many biomechanics fields on different organs (Couteau et al., 2000; Castellano-Smith et al., 2001; Fernandez et al., 2004; Sigal et al., 2008; Bucki et al., 2010; Bijar et al., 2016; Park et al., 2017), full-body models (Davis et al., 2016; Beillas and Berthet, 2017; Liu et al., 2020), as well as for detailed (Giudice et al., 2020; Giudice et al., 2021; Li et al., 2021; Montanino et al., 2021) and simplified brain models (Hu et al., 2007; Ji et al., 2011; Ji et al., 2015b; Wu et al., 2019). A typical procedure involves image registration (rigid or affine and followed by nonlinear registrations), from which a displacement field representing the geometrical difference between the subject and baseline model is obtained. The displacement field is then applied to morph the baseline model, resulting in a personalized model with updated nodal coordinates while preserving element connections. The displacement field derived from image registrations should generally comply with continuum mechanics conditions on motion, requiring diffeomorphic, non-folding, and one-to-one correspondence to avoid excessive element distortions (Bucki et al., 2010).

In particular, deformable image registration-based mesh morphing has been applied to personalize detailed brain models of healthy subjects (Giudice et al., 2020; Giudice et al., 2021; Li et al., 2021; Montanino et al., 2021). However, despite intensive efforts, inter-subject registration between brains with significant anatomical differences is still challenging within neuroimaging field with limited registration accuracy (Kim et al., 2015). Moreover, when applying image registration for mesh morphing, there is a higher requirement on the smoothness of the obtained displacement field to ensure acceptable element quality in the morphed mesh. Therefore, one major challenge for using mesh morphing to generate subject-specific FE head models is how to design an image registration pipeline that leads to high registration accuracy, meanwhile, not causes excessive element distortions. In a previous study (Li et al., 2021), we proposed a hierarchical image registration pipeline that allows efficient generation of subject-specific head models for healthy adult subjects. But a pipeline that allows morphing a baseline model to subjects with significant anatomical differences is yet to be developed.

Thus, this study aims at developing a personalization framework capable of generating subject-specific head models across the lifespan and for pathological brains with significant anatomical changes. The framework consists of hierarchical multiple feature and multimodality image registration pipelines, mesh morphing, and mesh grouping. Six subject-specific head models are generated to demonstrate its capacity, including a newborn, 1Y, 2Y, adult, 92Y, and a hydrocephalus brain. The results show that the framework is robust to generate subject-specific models across the lifespan and for pathological brains with significant anatomical changes by morphing a baseline model. This framework helps to facilitate personalized simulations in many fields within neurosciences, especially for studying TBIs in which personalized simulations are hindered due to the meshing challenge.

Materials and Methods

Subjects

Images of the six subjects (Figures 1A–F) are acquired from previously published open-access datasets, except the hydrocephalus brain is from the author’s previous study. The baseline ICBM image (Figure 1G) corresponds to the baseline head model. For detailed preprocessing steps for these images, the readers are referred to the original studies. A brief description is found below and summarized in Table 1.

• Images of a newborn (denoted as 0Y afterward), 1Y, and 2Y are obtained from the UNC Infant 0-1-2 atlases (Shi et al., 2011) constructed based on 95 subjects with complete 0-1-2Y longitudinal scans of T1W and T2W images acquired with a 3T MRI scanner. Each atlas consists of T1W images, tissue probability maps, and anatomical parcellation maps.

• Image of a single subject from the WU-Minn HCP database in the 26–30 age group, including T1W and T2W images, was acquired with a 3T MRI scanner (Van Essen et al., 2013).

• Image of an elderly (92Y) from the Brain Imaging of Normal Subjects (BRAINS) atlas was created from 48 healthy elderly subjects within age group 9193Y as detailed by Dickie et al. (2016). The atlas contains T1W and tissue probability maps.

• Image of a hydrocephalus subject with a mass lesion at the brain stem front is reused from a previous study (Li and von Holst, 2013).

• The 1-mm isotropic ICBM 2009c Nonlinear Symmetric template (Fonov et al., 2009; Fonov et al., 2011) was constructed based on T1W images from 152 subjects between 18.5–43.5Y acquired on a 1.5 T MRI scanner.

FIGURE 1
www.frontiersin.org

FIGURE 1. Image data used in this study. Axial, coronal, and sagittal views of (A) 40-week-old newborn, (B) 1-year-old, (C) 2-year-old, (D) an adult, (E) an elderly of 92-year-old, (F) hydrocephalus brain, and (G) the ICBM baseline (the same length scale applies).

TABLE 1
www.frontiersin.org

TABLE 1. Subjects involved in this study.

Baseline FE Head Model

A previously developed FE head injury model (the ADAPT model) (Li et al., 2021) serves as a baseline in this study, which is morphed to obtain subject-specific head models. The ADAPT model has been generated based on and has the same geometry as the ICBM template. The model includes the brain, skull, meninges, CSF, and superior sagittal sinus (SSS) (Figure 2). The brain is divided into primary structures of cerebral gray matter (GM) (i.e., cerebral cortex), cerebral white matter (WM), corpus callosum (CC), brain stem (BS), cerebellum GM and WM, thalamus, and hippocampus. The cerebrum is further divided into frontal, frontal, parietal, temporal, and occipital lobes; CSF is divided into outer CSF and ventricular system including lateral ventricles and 3rd and 4th ventricles connected by the cerebral aqueduct. Continuous mesh is used between brain components throughout the model. The total number of elements in the head model is 4.4 million hexahedral and 0.54 million quad elements. The minimum Jacobian in the brain is 0.45. The brain is modeled as hyper-viscoelastic material to account for large deformations and strain rate dependence of the tissue. Pia, dura/falx/tentorium are modeled with nonlinear hyperelastic material using simplified rubber/foam based on the average stressstrain experimental data (van Noort et al., 1981; Aimedieu and Grebe, 2004). The model has been validated against experimental data of close to or injury level brainskull relative motion, brain strain, and intracranial pressure. Details of the model development, validation, and capacity to study brain responses under impact are presented earlier (Li et al., 2021).

FIGURE 2
www.frontiersin.org

FIGURE 2. The baseline ADAPT head model with major components illustrated. The meshes are not shown for a better illustration.

Personalization Framework for Subject-specific Head Model Generation

The personalization framework consists of image registration pipelines, mesh morphing, and mesh grouping (Figure 3). Image registration is an essential part of the framework. A complete registration pipeline involves hierarchical registrations with multiple features and multimodality images shown at the lower row of Figure 3. The sum of dense displacement fields obtained from each registration step is used to morph the baseline head model to obtain subject-specific models. Afterward, the WM of the morphed brain is regrouped according to the segmented WM image mask of the subject, resulting in the final subject-specific model. Details of each component of the framework are presented in the following subsections.

FIGURE 3
www.frontiersin.org

FIGURE 3. Overview of the personalization framework.

Registration Pipeline With Multiple Features and Multimodality Imaging

A complete registration pipeline contains five steps (Figure 3 lower row). First, Demons registration is performed between the segmented cranial masks of the baseline ICBM (corresponding to the baseline ADAPT head model) and the subject after being rigidly aligned, resulting in a transformation, i.e., dense displacement field gdemo. Second, Demons registration of features is performed, obtaining gf1,gf2. Third, Dramms registration is performed on T1W images inherited from Demons steps, obtaining gdram. Next, Dramms registration is performed with multimodality images, obtaining gm1,gm2. Finally, brain lesions are handled by more Demons feature registration steps, obtaining gfn. In all registration steps, the subject’s image serves as moving image, and the baseline ICBM image serves as fixed image. Note that features in this study refer to the segmented binary images of anatomical regions such as lateral ventricles, corpus callosum, or lesion. The input images to Demons registration steps are segmented binary masks. Thus, these steps capture local anatomical changes between the moving and fixed images only in size and shape, while subsequent Dramms registrations capture the internal anatomical differences within the volumes of the binary image masks. The practical usage of the registration pipeline is demonstrated in Application of the Framework for Subject-specific Model Generation with six subjects.

In particular, for all Demons steps described above, the diffeomorphic Demons registration algorithm (Vercauteren et al., 2009) implemented in the open-source software Slicer 3D is used. Dramms registration algorithm (Ou et al., 2011) implemented as open-source code by the authors (Dramms version 1.5.1, 2018) is used on MRI images of different modalities. Note that for all the six subjects, a smoothness weight, i.e., the -g option (see DRAMMS Software Manual), is always set to 1.0 in Dramms registration to ensure a smooth displacement field.

Morphing

The sum of dense displacement fields from all registration steps (Eq. 1) represents the anatomical differences between the subject and the baseline ICBM images.

gsubj=gdemo+gf1+gf2+gdram+gm1+gm2 +gfn.(1)

As the baseline, ADAPT model is in the same space as the ICBM image; thus, applying gsubj to the baseline head model (Figure 2) leads to a subject-specific head model of the subject. For this, the following step is performed to morph the nodes of the baseline head model to new positions:

xi=Xi+ui,(2)

where Xi is the nodal coordinate of node i, ui  is the linearly interpolated displacement vector at node n from gsubj, xi is the updated nodal coordinate, together with the same element definitions as the baseline, forming a subject-specific head model.

Grouping of WM

To capture the subject’s WM, the morphed brain elements are regrouped based on the segmented binary image of the subject’s cerebral WM. This is achieved by assigning brain FE elements as WM based on Cartesian coordinates of the segmented WM voxels with the following procedures:

- For each element, all WM voxels inside or intersect to a single element of the brain are identified based on spatial coordinates.

- The eight vertices and one centroid of each voxel (i,j,k) are judged; vertices gain a weight of one if falling inside the element; the centroid gains a weight of two if falling inside the element. Weights of the eight vertices and the centroid of the voxel add up, resulting in a total weighting factor for each voxel wi,j,k

wi,j,k=m=19wi,j,km(3)

- Finally, weights of each voxel belong to the same label (e.g., the segmented binary image with label A) added up, obtaining a final weight factor for each label. The element is grouped to the label with the largest weight.

wA=(i,j,k)Awi,j,k.(4)

Figure 4 shows the regrouped WM elements of the morphed brain enclosed by the reconstructed surface of the segmented WM.

FIGURE 4
www.frontiersin.org

FIGURE 4. Regrouped WM based on subjects’ WM image mask for all the six subjects. The red color shows the WM elements, and the white transparent shows the surfaces reconstructed from the subject’s segmented WM image mask.

Evaluation of Personalization Accuracy

To evaluate registration accuracy, the baseline ICBM image (imgbaseline), which corresponds to the baseline head model, is warped via the inverse of displacement fields from each registration step (gsubj1) (Eq. 5), resulting in a warped image (imgwarped) (Eq. 6).

gsubj1=gfn1(gm21(gm11(gdram1(gf21(gf11(gdemo1))))),(5)
imgwarped= gfn1(gm21(gm11(gdram1(gf21(gf11(gdemo1(imgbaseline))))))).(6)

DICE and 95th percentile Hausdorff distance (HD95) between the imgwarped and subjects’ images are then calculated to evaluate registration performance. As imgwarped corresponds to the personalized subject-specific model, both metrics also reflect the personalization accuracy of the generated subject-specific models.

To calculate DICE and HD95, automated segmentation is performed using the software FreeSurfer (version 7.1.0) with the default brain segmentation pipeline (recon-all) for both the warped ICBM and subjects’ T1W images. The segmented binary masks for the whole brain and local regions of cerebral GM, WM, CC, BS, hippocampus, thalamus, and cerebellum are used for DICE and HD95 calculation. For the cranial mask, the metrics are calculated based on manually segmented cranial by thresholding followed by noise removal. Similarly, one sagittal slice of CC is manually segmented and used to calculate both matrices. The use of manual segmentation for both regions is due to the insufficient quality (not reflecting the actual anatomy) by recon-all for the current dataset. Note that these segmented binary masks are only used for DICE and HD95 calculation, and the quality of the automatic segmentation has no influence on the subject-specific mesh development process.

DICE

DICE is a single metric to measure the spatial overlap between images defined as twice the number of elements common to both sets divided by the sum of the number of elements in each set (Ou et al., 2014):

DICE(A,B)=2|AB||A|+|B|,(7)

where A and B denote the binary segmentation labels, |A| and |B| are the number of voxels in each set, and |AB| is the number of shared voxels by A and B. The DICE value of 0 implies no overlap between both, whereas a DICE coefficient of one indicates perfect overlap between the warped and the target image.

HD95

Hausdorff distance is defined as

HD(C,D)=max(h(C,D),h(D,C)),(8)

where C,D are the two sets of vertices from two segmented images

h(C,D)=max cC  maxdDcd.(9)

The 95th percentile Hausdorff distance (HD95) is used following earlier studies (Ou et al., 2011; Ou et al., 2014). HD95 ranges from 0 to above; a lower value indicates a better registration accuracy between the warped and the target image.

Application of the Framework for Subject-specific Model Generation

A complete registration pipeline is only needed for the most challenging case; fewer registration steps are sufficient for brains with small anatomical differences compared with the baseline. The following three typical subtypes of the pipeline are used to generate subject-specific models for the six subjects.

Type I This is the basic pipeline containing two steps: Demons registration of the cranial mask and Dramms registration of T1W image. This two-step pipeline has been shown to achieve good registration accuracy for six healthy adult subjects (Li et al., 2021). The capacity of this pipeline is further demonstrated with a 2Y brain.

Type II Multiple feature steps are added to the Type I pipeline, allowing align brains with significant anatomical changes. The capacity of this pipeline is demonstrated with a hydrocephalus brain using three feature steps.

Type III Multi-modality imaging registration steps are added to the Type I pipeline to improve brain alignment as demonstrated with an adult brain.

2YO Model Generation via Pipeline Type I: Two Steps

First, T1W images of ICBM and the 2Y brain are segmented to obtain the cranial masks, which are used as input for Demons registration, from which a dense displacement field gdemo is obtained (Figure 5). Next, the baseline ICBM T1W image is warped by gdemo1. The warped ICBM and subject’s T1W images are then skull stripped and serve as input for Dramms registration, obtaining gdram. The two displacement fields add up gsubj_2YO=gdemo+gdram which is used to morph the baseline mesh, obtaining the subject-specific head model. Finally, the baseline ICBM image warped to the subject is obtained via imgwarped_2Y= (gdram1(gdemo1(imgbaseline))), which is compared with the T1W image of the 2Y brain to evaluate personalization accuracy.

FIGURE 5
www.frontiersin.org

FIGURE 5. Type I pipeline applied for personalizing the baseline ADAPT model to a subject-specific model of a 2Y. The pipeline consists of two steps: (i) Demons registration with cranial masks; (ii) Dramms registration with T1W image. The displacement field obtained from each step is visualized on the grid together with the warped baseline ICBM images to show its effect. The final warped ICBM is overlaid with the subject’s image to visualize registration accuracy.

Hydrocephalus Model via Pipeline Type II: Multiple Features

The workflow is similar to the above, but three additional feature steps are added to capture the enlarged LV, deformed CC, and brain lesion, resulting in five dense displacement fields that add up as gsubj_hydro=gdemo+gf1+gf2+gdram+gf3, which is used to morph the baseline mesh and obtain a subject-specific head model (Figure 6). The baseline ICBM image warped to the subject is obtained via imgwarped_hydro=(gdram1(gf21(gf11(gdemo1(imgbaseline))))), which is compared with the T1W image of the hydrocephalus subject to evaluate personalization accuracy.

FIGURE 6
www.frontiersin.org

FIGURE 6. Type II pipeline applied for personalizing the baseline ADAPT model to a subject-specific model of a hydrocephalus brain. The pipeline consists of five steps: (i) Demons registration with cranial masks; (ii) Demons registration with segmented lateral ventricle (LV) mask for capturing the enlarged LV; (iii) Demons registration with segmented CC mask for capturing the CC shape; (iv) Dramms registration with T1W image for capturing local brain anatomy; (v) Demons registration to drag back the skull mesh which is pushed due to the lesion in the cranial mask in step (i). The displacement field obtained from each step is visualized on the grid together with the warped baseline ICBM images and morphed meshes to show its effect. The final warped ICBM is overlaid with the subject’s image to visualize registration accuracy.

Adult Brain Model via Pipeline Type III: Multimodality

The workflow is similar to Type I, but an additional multimodality T2W registration step is performed to further align the LVs resulting in three dense displacement fields added up as gsubj_adult=gdemo+gdram+gm1, which is used to morph the baseline mesh and obtain the subject-specific head model (Figure 7). The baseline ICBM image warped to the subject is obtained via imgwarped_adult= gm11(gdram1(gdemo1(imgbaseline))), which is compared with the T1W image of the hydrocephalus subject to evaluate personalization accuracy.

FIGURE 7
www.frontiersin.org

FIGURE 7. Type III pipeline applied for personalizing the baseline ADAPT model to a subject-specific model of an adult subject. The pipeline consists of three steps: (i) Demons registration with cranial masks; (ii) Dramms registration with T1W image; (iii) Dramms registration with T2W image for further alignment. The displacement field obtained from each step is visualized on the grid together with the warped baseline ICBM images to show its effect. The final warped ICBM is overlaid with the subject’s image to visualize registration accuracy.

Pipeline for the 0Y, 1Y, and the 92Y

The 92Y uses Type II pipeline, similar to the hydrocephalus subject, except only one feature step for LV is used, i.e., gsubj_92Y=gdemo+gf1+gdram. Interestingly, the Dramms registration captures the thinning of CC without the CC feature step as for the hydrocephalus subject. It could be due to the higher quality image of the 92Y than the hydrocephalus brain.

The 1Y could use the same pipeline as the 2Y. However, in this study, an alternative approach is used, using the 2Y as an intermediate step, i.e., align 1Y T1W (as moving image) to that of the 2Y (as fixed image) by Dramms registration, obtaining a displacement field g1Y_to_2Y. The final displacement field used for personalizing the baseline model to the 1Y writes: gsubj_1Y= gsubj_2YO+g1Y_to_2Y, and the warped image is obtained via imgwarped_1Y= (g1Y_to_2Y1(gsubj_2Y1(imgbaseline))).

The 0Y uses the same pipeline as the 1Y by having 2Y as an intermediate step, i.e., align the 0Y T1W (as moving image) to that of the 2Y (as fixed image) by Dramms registration. Since the cerebellum for the 0Y was stripped in the original database (Figure 1A), a paired T1W image of the 2Y with cerebellum stripped (readily available in the database) is used for registration, from which a displacement field g0Y_to_2Y is obtained. The remaining steps for personalization and image warping are the same as for the 1Y described above. Note that although the registered images do not have cerebellum, the obtained displacement field g0Y_to_2Y defined in the entire image space does cover the cerebellum region despite values close to zero. The displacement field when used to morph the baseline model that has cerebellum, resulting in a final subject-specific model of the 0Y with cerebellum included.

Results

Subject-specific Head Models and Element Quality

The generated head models (Figure 8) and cross-sections (Figure 9) demonstrate the capacity of the framework for generating subject-specific head models with significant anatomical differences; all morphed from a baseline model. Especially, the extensively enlarged LVs and the varying shapes of CC in the hydrocephalus and the elderly 92Y brain are captured (Figures 10, 11). The element quality for the models is listed in Table 2, showing that most brain elements (95.9 ± 1.5% on average for the six subjects) have a Jacobian over 0.5, and the minimum Jacobian in all the six head models is above 0.13 (in the hydrocephalus brain). In this study, the mesh quality is considered satisfactory when at least 95% of the elements have a Jacobian over 0.5.

FIGURE 8
www.frontiersin.org

FIGURE 8. Six subject-specific head models generated including the 0Y, 1Y, 2Y, adult, 92Y, and a hydrocephalus brain (the same length scale applies).

FIGURE 9
www.frontiersin.org

FIGURE 9. Six subject-specific models are aligned together, showing the generated models have widely varying intracranial volumes (upper row) and significant anatomical differences as exemplified with lateral ventricles and corpus callosum (lower row).

FIGURE 10
www.frontiersin.org

FIGURE 10. (A) T1W image of the ICBM baseline; (B) ICBM baseline warped to the six subjects; (C) T1W image of the six subjects. Transverse, sagittal, and coronal cross-sections are captured for each brain.

FIGURE 11
www.frontiersin.org

FIGURE 11. T1W image of the subject is overlaid with the segmented binary masks of the warped baseline, including the cranial mask, the brain, and local brain regions of the cerebellum, hippocampus, thalamus, CC, BS, and lateral ventricles (A). Boxplots of DICE and 95HD. The boxplots show the median, minimum, and maximum values shown (B).

TABLE 2
www.frontiersin.org

TABLE 2. Element quality of the baseline ADAPT model and the six subject-specific head models generated by morphing.

Personalization Accuracy

The baseline ICBM image (Figure 10A) is warped to the six subjects. The warped images (Figure 10B) and subjects’ images (Figure 10C) are compared to evaluate registration accuracy. The segmented binary masks of the final warped baseline and subjects are overlaid to further visualize personalization accuracy (Figure 11A). The evaluated masks include cranial, brain, and six local brain regions. The boxplots of the DICE and HD95 are presented in Figure 11, with values listed in Tables 3, 4. The average DICE scores are all >0.9 for the cranial mask, the brain, cerebellum, CC, being 0.97, 0.90, 0.89, and 0.94, respectively. Since the cerebellum for the 0Y subject image has been stripped, the evaluation of registration accuracy is without this region. The average DICE score for LV is 0.80. DICE score can be improved by incorporating multimodality step, e.g., with T2W image that has higher contrast for CSF/LVs. For example, the pipeline for the adult subject adding the T2W multimodality step improves personalization accuracy than previously achieved (see Supplementary Appendix S1). The DICE values are comparable to that achieved in neuroimaging field (Ou et al., 2014) despite the higher requirement on the smoothness of displacement field for satisfactory element quality in the personalized head FE models.

TABLE 3
www.frontiersin.org

TABLE 3. DICE coefficients for the six subjects.

TABLE 4
www.frontiersin.org

TABLE 4. HD95 for the six subjects.

Hydrocephalus and the Elderly Brain: Importance of the Feature Step and the Higher Requirement on Displacement Smoothness for Mesh Morphing

The mesh after each morphing step shown in Figure 12 illustrates the effect of the feature steps, which allow capturing subject’s cranial shape (Figure 12A), enlarged LVs (Figure 12B), CC (Figure 12C), as well as pushing back of the skull mesh (Figure 12E), while local brain structures are captured by Dramms registration (Figure 12D). The meshes from these intermediate steps are morphed from the baseline ADAPT head model with displacement fields obtained via the image registration pipeline shown in Figure 6.

FIGURE 12
www.frontiersin.org

FIGURE 12. Morphed FE meshes after each of the five steps for the hydrocephalus subject.

To further illustrate the importance of the feature registrations, a parametric pipeline without the LV feature step is performed for the 92Y brain, i.e., a complete pipeline writes gsubj_92Y=gdemo+gdram and imgwarped_92Y= gdram1(gdemo1(imgbaseline)) (Figure 13). The results show that the Dramms registration, even using the largest allowable smoothness weight (g = 1), leads to FE mesh with negative Jacobian in some elements. For example, one FE element at the frontal horn of the LVs with a Jacobian (JFE= 0.5) in the baseline mesh (Figure 13A, right upper), when morphed by the parametric pipeline resulting in a negative Jacobian (JFE= -0.05) (Figure 13C, right upper). In contrast, when morphed by the original three-step pipeline, the same element has a positive value (JFE= 0.31) (Figure 13B, right upper). This parametric pipeline also results in lower registration accuracy than the three-step pipeline (Figures 13A,B). This example also demonstrates the higher requirement on the smoothness of displacement than in the neuroimaging field when only Jacobian of the displacement field (Jimg) is of concern (more detailed analysis presented in Supplementary Appendix S2).

FIGURE 13
www.frontiersin.org

FIGURE 13. Parametric pipeline for the 92Y without the LV feature step compared with the default pipeline. The warped baseline image by the parametric pipeline (B) does not capture the enlarged LV compared with that achieved by the default pipeline (A). The parametric pipeline leads to negative Jacobian in some elements in the personalized mesh (right figure), although the Jacobian map (Jimg) of the final obtained displacement field is all positive. One representative axial slice is shown with the minimum and the maximum value of Jimg in the entire brain indicated (C,D).

Discussions

This study presents a personalization framework for the efficient generation of subject-specific head models. The framework consists of hierarchical multiple feature and multimodality imaging registration pipelines, mesh morphing, and mesh grouping. The registration pipeline achieves competitive registration accuracy despite a higher requirement on the smoothness of the displacement field concerning the element quality of the morphed mesh. The Demons feature registration steps capture significant anatomical differences, allowing a good initialization before applying Dramms registration to further capture the inter-subject anatomical details. The Dramms registration step with multimodality imaging further improves brain alignment. As a final step of the framework, mesh grouping of WM according to the subject’s image mask allows incorporating subject-specific WM directly. The framework is successfully applied to subjects across the lifespan and a hydrocephalus brain with significant anatomical changes, achieving competitive personalization accuracy. The results demonstrate that the framework can personalize the baseline head model to brains with significant anatomical differences, resulting in subject-specific models ready for personalized simulations without manual repairing. To the knowledge of the author, this is the first study aligning such a broad scope of brain images suitable for mesh morphing.

The efficiency of the hierarchical two-step pipeline combining Demons and Dramms (Type I) has been previously assessed with six healthy adult subjects that have high-quality T1W images (Li et al., 2021). In this study, an extended pipeline is proposed for obtaining high-quality alignment across heterogeneous data of lifespan and for pathological brains with significant anatomical changes by introducing multiple feature steps as demonstrated with the hydrocephalus (Figure 12) and the 92Y brain (Figure 13). The registration accuracy for these more challenging cases is comparable with the six healthy adults, with average DICE scores for the cerebellum, CC, and brain all above 0.89. Notably, the average DICE score for LV for the six subjects in this study is 0.80, higher than that of the six adult subjects (0.71) (Li et al., 2021). Note that the same adult subject in an early study (Li et al., 2021) (subject ID 771354) is used here by adding T2W multimodality registration step (Type III). T2W images with higher contrast for CSF/LVs improve personalization accuracy compared with previous results (see Supplementary Appendix S1). The mesh grouping step incorporates subject-specific WM directly, which is important for infant models to accurately capture the rapid transition between GM and WM in early infancy. Thus, the promising performance demonstrates the potential of the framework to personalize the baseline model to almost any brains with significant anatomical changes. Besides hydrocephalus, personalized models for brains with other structural changes such as decompressive craniotomy with brain expanded outside the skull (Holst et al., 2012) can also be achieved.

Inter-subject registration between brains with significant anatomical differences is still challenging in neuroimaging field and has limited registration accuracy (Kim et al., 2015); even more challenging is to apply image registration for mesh morphing due to the higher requirement on the smoothness of the obtained displacement fields concerning element quality of the morphed mesh. There is often a trade-off between registration accuracy and element quality, and higher registration accuracy tends to worsen element quality according to the experience with the six subjects in this study. Especially, FE elements become invalid if their Jacobian become negative, which are not accepted by most FE analysis software. While in neuroimaging field, for physically plausible morphing, only positive Jacobian determinant of displacement field is to be ensured, which is often a looser requirement than FE Jacobian (see detailed analysis in Supplementary Appendix S1). Despite the higher requirement, this study achieves competitive registration accuracy compared with that reported in the neuroimaging field (Ou et al., 2014). For example, a previous study reported Jaccard index below 0.6 for all brain regions using popular deformable registration algorithms for inter-subject registration (Ou et al., 2014), while the average Jaccard index (converted from DICE according to Jaccard index = DICE/(2-DICE) (Ou et al., 2014)) for the six subjects is all above 0.66 for all regions in this study.

The applications of the framework show that different pipelines can be used depending on the anatomical differences between the subject and baseline ICBM, as well as the subject’s image quality and available imaging modalities. For brains that are similar to the baseline, Type I pipeline with fewer steps is sufficient, while for brains with significant anatomical differences compared with the baseline, e.g., the hydrocephalus and elderly brain, Type III pipeline is needed to achieve a proper alignment. Furthermore, when T2W images are available, multimodality image allows better alignment of the brain and CSF/LVs. In principle, more multimodality registrations can be performed if available from the subject as the baseline image contains imaging modalities of T1W, T2W, proton density (PD), and tissue probability maps. Besides, more feature steps can be introduced to handle even more challenging cases. Choosing the proper pipeline for a specific case needs trial and error. An overall guideline is to start from Type I then add more registration steps if needed. Note that the multiple feature steps can be combined into one image with multiple binary masks and perform Demons registration at once. However, one feature in each step, as done in this study, tends to be more robust. Furthermore, the framework, though demonstrated with the ADAPT baseline head model, is equally applicable for personalizing other head models as a baseline, e.g., models with tetrahedral elements as commonly used for tDCS, TMS, TUS, as well as smoothed-voxel brain models.

Compared with existing studies registering adult brains, fewer studies align infant brains, which are more challenging partially due to the rapid development of brain anatomy within the first year, especially T1W images are inversed with densities. Not only more challenging for registration algorithms but also the evaluation of performance is also more difficult as most segmentation algorithms are developed based on adult images, such as FreeSurfer. It’s worth noting that the lowest registration accuracy in all brains is for the thalamus in the 0Y; a visual check shows FreeSurfer automatically segmented thalamus not accurate enough. Future studies can employ infant Freesurfer (Zöllei et al., 2020) for more accurate segmentation for infant brain images, thus allows more objective evaluation of personalization accuracy. For adult brain mesh morphing, image registration-based morphing pipelines proposed earlier show promising performance in generating detailed subject-specific head models of healthy adult brains (Giudice et al., 2020; Li et al., 2021) while the framework proposed in this study allows generating models across the lifespan and for brains with significant anatomical changes, which can be used for studying age-specific and groupwise TBIs. Especially, brains with neurological diseases such as hydrocephalus with extensively enlarged LVs mimicking the elderly brain may provide a possible clue for new insights into TBIs. The approach also opens the opportunity for studying how a potentially vulnerable brain, e.g., a hydrocephalus patient, may sustain a TBI injury risk under fall impact, especially hydrocephalus patients who are more prone to fall. Until today, the biomechanics of TBIs in these groups are much understudied, partially due to the meshing challenge.

Compared with the many existing studies of TBIs for healthy adults, the injury mechanisms of infants and children are understudied. There are few child/infant head models (Li et al., 2011; Giordano et al., 2017; Li et al., 2017). In addition to the meshing challenge for adult models, the development of additional unique features of suture and fontanel plays an essential role in head impact response (Li et al., 2017). Previously, mesh morphing has also been used for morphing a baseline infant head model to different ages using radial basis function (RBF) to interpolate the displacement field obtained from land markers the anatomical features of suture and skull surface (Li et al., 2011). Unlike the image registration-based morphing, the RBF approach needs manual indentation of land markers, which is often tedious (Wu et al., 2019). The RBF approach also does not account for brain anatomies. Comparatively, the morphed detailed infant brain models in this study, when combined with the detailed skull and scalp models (Li et al., 2017; Li et al., 2019), will allow studying brain injury biomechanics under impact for infant head model for abusive head trauma with important legal applications for forensic diagnosis. The newborn infant head models may be used for studying delivery-related neurotrauma and studying new intervention approaches for clinical problems.

Some limitations and future works need to be mentioned here. First, the proposed framework allows efficient generation of subject-specific head models with competitive personalization accuracy and satisfactory element quality without mesh repairing. However, the morphing technique involves manual intervention when selecting which morphing pipeline to use. Thus, there could be user-to-user variability based on which pipelines are chosen and concerns regarding repeatability. For example, selecting improper pipelines could result in reduced morphing accuracy, and certain regions may not be morphed accurately if they are not selected as features by the user. Secondly, this framework requires segmentation for the Demons steps, which may need manual effort to ensure accurate segmentation and would require significant time and effort for large-scale studies. Nevertheless, considering the challenge of generating subject-specific head models, this effort is considered acceptable. Note that this morphing technique generates subject-specific models from a geometric perspective only and does not account for subject-specific material properties. Thirdly, the current framework allows generating head models reflecting the subject’s internal brain structures, but the major sulci and gyri lines are not evaluated like most studies in the neuroimaging field. It should also be noted that the framework does not ensure the same characteristic lengths among generated models of different sizes; the infant brains, in general, have smaller elements than that of an adult model. In this regard, the block-based method has an advantage that allows adjusting mesh densities to maintain similar element characteristic lengths (Mao et al., 2013b). Furthermore, Dramms registration algorithm is chosen for registering brain MRI images in this study since it has a clear advantage to align largely different anatomies such as the ventricles in comparison with other popular registration algorithms (Ou et al., 2014). However, other algorithms, such as implemented in ANTs (Avants et al., 2011) and DARTEL (Vercauteren et al., 2009), when used within the current framework, may achieve similar performance, but it is yet to be investigated. Finally, the framework can be extended to include more registration steps and other advanced nonlinear registration algorithms to handle even more challenging cases.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author Contributions

The author confirms being the sole contributor of this work and has approved it for publication.

Funding

The present study was supported by research funds from KTH Royal Institute of Technology (Stockholm), and the Swedish Research Council Grants (2016-04203 and 2020-04724).

Conflict of Interest

The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

Open access funding provided by KTH Royal Institute of Technology. The author thanks the reviewers whose valuable comments and suggestions significantly improved the quality of the manuscript. This work benefits greatly from the open-source imaging dataset and registration algorithms within the neuroimaging field. The WUM HCP data of the adult subject used in this study were provided [in part] by the Human Connectome Project, WU-Minn Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) funded by the 16 NIH Institutes and Centers that support the NIH Blueprint for Neuroscience Research; and by the McDonnell Center for Systems Neuroscience at Washington University. The data handling was enabled by resources in projects [SNIC 2021/5-406] and [SNIC 2021/5-459] provided by the Swedish National Infrastructure for Computing (SNIC) at PDC, partially funded by the Swedish Research Council through grant agreement no. 2018-05973

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2021.706566/full#supplementary-material

References

Aimedieu, P., and Grebe, R. (2004). Tensile Strength of Cranial Pia Mater: Preliminary Results. J. Neurosurg. 100 (1), 111–114. doi:10.3171/jns.2004.100.1.0111

CrossRef Full Text | Google Scholar

Alekseichuk, I., Falchier, A. Y., Linn, G., Xu, T., Milham, M. P., Schroeder, C. E., et al. (2019). Electric Field Dynamics in the Brain during Multi-Electrode Transcranial Electric Stimulation. Nat. Commun. 10 (1), 2573. doi:10.1038/s41467-019-10581-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Atsumi, N., Nakahira, Y., and Iwamoto, M. (2016). Development and Validation of a Head/brain FE Model and Investigation of Influential Factor on the Brain Response during Head Impact. Ijvs 9 (1), 1–23. doi:10.1504/ijvs.2016.077145

CrossRef Full Text | Google Scholar

Avants, B. B., Tustison, N. J., Song, G., Cook, P. A., Klein, A., and Gee, J. C. (2011). A Reproducible Evaluation of ANTs Similarity Metric Performance in Brain Image Registration. Neuroimage 54 (3), 2033–2044. doi:10.1016/j.neuroimage.2010.09.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Baker, T. J. (2005). Mesh Generation: Art or Science? Prog. Aerospace Sci. 41 (1), 29–63. doi:10.1016/j.paerosci.2005.02.002

CrossRef Full Text | Google Scholar

Beillas, P., and Berthet, F. (2017). An Investigation of Human Body Model Morphing for the Assessment of Abdomen Responses to Impact against a Population of Test Subjects. Traffic Inj. Prev. 18 (Suppl. 1), S142–S147. doi:10.1080/15389588.2017.1307971

PubMed Abstract | CrossRef Full Text | Google Scholar

Bijar, A., Rohan, P.-Y., Perrier, P., and Payan, Y. (2016). Atlas-Based Automatic Generation of Subject-specific Finite Element Tongue Meshes. Ann. Biomed. Eng. 44 (1), 16–34. doi:10.1007/s10439-015-1497-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Bucki, M., Lobos, C., and Payan, Y. (2010). A Fast and Robust Patient Specific Finite Element Mesh Registration Technique: Application to 60 Clinical Cases. Med. Image Anal. 14 (3), 303–317. doi:10.1016/j.media.2010.02.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Castellano-Smith, A. D., Hartkens, T., Schnabel, J., Hose, D. R., Liu, H., Hall, W. A., et al. (2001). “Constructing Patient Specific Models for Correcting Intraoperative Brain Deformation,” in Medical Image Computing and Computer-Assisted Intervention – MICCAI 2001. Editors W. J. Niessen, and M. A. Viergever (Berlin, Heidelberg2001: Springer).

CrossRef Full Text | Google Scholar

Couteau, B., Payan, Y., and Lavallée, S. (2000). The Mesh-Matching Algorithm: an Automatic 3D Mesh Generator for Finite Element Structures. J. Biomech. 33 (8), 1005–1009. doi:10.1016/s0021-9290(00)00055-5

CrossRef Full Text | Google Scholar

Datta, A., Bansal, V., Diaz, J., Patel, J., Reato, D., and Bikson, M. (2009). Gyri-precise Head Model of Transcranial Direct Current Stimulation: Improved Spatial Focality Using a Ring Electrode versus Conventional Rectangular Pad. Brain Stimulation 2 (4), 201–207. 207 e1. doi:10.1016/j.brs.2009.03.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Datta, A., Truong, D., Minhas, P., Parra, L. C., and Bikson, M. (2012). Inter-Individual Variation during Transcranial Direct Current Stimulation and Normalization of Dose Using MRI-Derived Computational Models. Front. Psychiatry 3, 91. doi:10.3389/fpsyt.2012.00091

PubMed Abstract | CrossRef Full Text | Google Scholar

Davis, M. L., Koya, B., Schap, J. M., and Gayzik, F. S. (2016). Development and Full Body Validation of a 5th Percentile Female Finite Element Model. Stapp Car Crash J. 60, 509–544. doi:10.4271/2016-22-0015

PubMed Abstract | CrossRef Full Text | Google Scholar

Dickie, D. A., Job, D. E., Rodriguez, D., Robson, A., Danso, S., Pernet, C., et al. (2016). Brain Imaging of Normal Subjects (BRAINS) Age-specific MRI Atlases from Young Adults to the Very Elderly (v1.0) [Dataset]. University of Edinburgh, Edinburgh Imaging, CCBS, BRAINS Imagebank. doi:10.7488/ds/1369

CrossRef Full Text | Google Scholar

Fernandez, J. W., Mithraratne, P., Thrupp, S. F., Tawhai, M. H., and Hunter, P. J. (2004). Anatomically Based Geometric Modelling of the Musculo-Skeletal System and Other Organs. Biomech. Model. Mechanobiology 2 (3), 139–155. doi:10.1007/s10237-003-0036-1

CrossRef Full Text | Google Scholar

Fonov, V., Evans, A. C., Botteron, K., Almli, C. R., McKinstry, R. C., Collins, D. L., et al. (2011). Unbiased Average Age-Appropriate Atlases for Pediatric Studies. Neuroimage 54 (1), 313–327. doi:10.1016/j.neuroimage.2010.07.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Fonov, V., Evans, A., McKinstry, R., Almli, C., and Collins, D. (2009). Unbiased Nonlinear Average Age-Appropriate Brain Templates from Birth to Adulthood. Neuroimage 47, S102. doi:10.1016/s1053-8119(09)70884-5

CrossRef Full Text | Google Scholar

Fornari, S., Schäfer, A., Jucker, M., Goriely, A., and Kuhl, E. (2019). Prion-like Spreading of Alzheimer's Disease within the Brain's Connectome. J. R. Soc. Interf. 16 (159), 20190356. doi:10.1098/rsif.2019.0356

CrossRef Full Text | Google Scholar

Giordano, C., Li, X., and Kleiven, S. (2017). Performances of the PIPER Scalable Child Human Body Model in Accident Reconstruction. PLoS One 12 (11), e0187916. doi:10.1371/journal.pone.0187916

PubMed Abstract | CrossRef Full Text | Google Scholar

Giudice, J. S., Alshareef, A., Wu, T., Gancayco, C. A., Reynier, K. A., Tustison, N. J., et al. (2020). An Image Registration-Based Morphing Technique for Generating Subject-specific Brain Finite Element Models. Ann. Biomed. Eng. 48 (10), 2412–2424. doi:10.1007/s10439-020-02584-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Giudice, J. S., Alshareef, A., Wu, T. T., Knutsen, A. K., Hiscox, L. V., Johnson, C. L., et al. (2021). Calibration of a Heterogeneous Brain Model Using a Subject-specific Inverse Finite Element Approach. Front. Bioeng. Biotechnol. 9. doi:10.3389/fbioe.2021.664268

CrossRef Full Text | Google Scholar

Giudice, J. S., Zeng, W., Wu, T., Alshareef, A., Shedd, D. F., and Panzer, M. B. (2019). An Analytical Review of the Numerical Methods Used for Finite Element Modeling of Traumatic Brain Injury. Ann. Biomed. Eng. 47 (9), 1855–1872. doi:10.1007/s10439-018-02161-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Hajiaghamemar, M., Lan, I. S., Christian, C. W., Coats, B., and Margulies, S. S. (2019). Infant Skull Fracture Risk for Low Height Falls. Int. J. Leg. Med. 133 (3), 847–862. doi:10.1007/s00414-018-1918-1

CrossRef Full Text | Google Scholar

Ho, J., Zhou, Z., Li, X., and Kleiven, S. (2017). The peculiar Properties of the Falx and Tentorium in Brain Injury Biomechanics. J. Biomech. 60, 243–247. doi:10.1016/j.jbiomech.2017.06.023

CrossRef Full Text | Google Scholar

Holst, H., Li, X., and Kleiven, S. (2012). Increased Strain Levels and Water Content in Brain Tissue after Decompressive Craniotomy. Acta Neurochir 154 (9), 1583–1593. doi:10.1007/s00701-012-1393-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Horgan, T. J., and Gilchrist, M. D. (2003). The Creation of Three-Dimensional Finite Element Models for Simulating Head Impact Biomechanics. Int. J. Crashworthiness 8 (4), 353–366. doi:10.1533/ijcr.2003.0243

CrossRef Full Text | Google Scholar

Horstemeyer, M. F., Panzer, M. B., and Prabhu, R. K. (2019). State-of-the-Art Modeling and Simulation of the Brain's Response to Mechanical Loads. Ann. Biomed. Eng. 47 (9), 1829–1831. doi:10.1007/s10439-019-02351-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Jin, X., Lee, J. B., Zhang, L., Chaudhary, V., Guthikonda, M., et al. (2007). Intraoperative Brain Shift Prediction Using a 3D Inhomogeneous Patient-specific Finite Element Model. Jns 106 (1), 164–169. doi:10.3171/jns.2007.106.1.164

CrossRef Full Text | Google Scholar

Huang, Y., Dmochowski, J. P., Su, Y., Datta, A., Rorden, C., and Parra, L. C. (2013). Automated MRI Segmentation for Individualized Modeling of Current Flow in the Human Head. J. Neural Eng. 10 (6), 066004. doi:10.1088/1741-2560/10/6/066004

CrossRef Full Text | Google Scholar

Ji, S., Ford, J. C., Greenwald, R. M., Beckwith, J. G., Paulsen, K. D., Flashman, L. A., et al. (2011). Automated Subject-specific, Hexahedral Mesh Generation via Image Registration. Finite Elem. Anal. Des. 47 (10), 1178–1185. doi:10.1016/j.finel.2011.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, S., Roberts, D. W., Hartov, A., and Paulsen, K. D. (2009). Brain-skull Contact Boundary Conditions in an Inverse Computational Deformation Model. Med. Image Anal. 13 (4), 659–672. doi:10.1016/j.media.2009.05.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ji, S., Zhao, W., Ford, J. C., Beckwith, J. G., Bolander, R. P., Greenwald, R. M., et al. (2015). Group-wise Evaluation and Comparison of white Matter Fiber Strain and Maximum Principal Strain in Sports-Related Concussion. J. Neurotrauma 32 (7), 441–454. doi:10.1089/neu.2013.3268

CrossRef Full Text | Google Scholar

Ji, S., Zhao, W., Ford, J. C., Beckwith, J. G., Bolander, R. P., Greenwald, R. M., et al. (2015). Group-Wise Evaluation and Comparison of White Matter Fiber Strain and Maximum Principal Strain in Sports-Related Concussion. J. Neurotrauma 32 (7), 441–454. doi:10.1089/neu.2013.3268

CrossRef Full Text | Google Scholar

Kim, M., Wu, G., Wang, Q., Lee, S.-W., and Shen, D. (2015). Improved Image Registration by Sparse Patch-Based Deformation Estimation. Neuroimage 105, 257–268. doi:10.1016/j.neuroimage.2014.10.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleiven, S. (2007). Predictors for Traumatic Brain Injuries Evaluated through Accident Reconstructions. Stapp Car Crash J. 51, 81–114. doi:10.4271/2007-22-0003

PubMed Abstract | CrossRef Full Text | Google Scholar

Kleiven, S., and von Holst, H. (2002). Consequences of Head Size Following Trauma to the Human Head. J. Biomech. 35 (2), 153–160. doi:10.1016/s0021-9290(01)00202-0

CrossRef Full Text | Google Scholar

Legon, W., Sato, T. F., Opitz, A., Mueller, J., Barbour, A., Williams, A., et al. (2014). Transcranial Focused Ultrasound Modulates the Activity of Primary Somatosensory Cortex in Humans. Nat. Neurosci. 17 (2), 322–329. doi:10.1038/nn.3620

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, N., Baldermann, J. C., Kibleur, A., Treu, S., Akram, H., Elias, G. J. B., et al. (2020). A Unified Connectomic Target for Deep Brain Stimulation in Obsessive-Compulsive Disorder. Nat. Commun. 11 (1), 3364. doi:10.1038/s41467-020-16734-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X. G., and von Holst, H. (2013). Quantification of Stretching in the Ventricular Wall and Corpus Callosum and Corticospinal Tracts in Hydrocephalus before and after Ventriculoperitoneal Shunt Operation. J. Appl. Maths. 2013, 350359. Article ID. doi:10.1155/2013/350359

CrossRef Full Text | Google Scholar

Li, X., Holst, H., and Holst, H. (2015). Finite Element Modeling of Decompressive Craniectomy (DC) and its Clinical Validation. Abse 2015, 1–9. doi:10.15764/abse.2015.01001

CrossRef Full Text | Google Scholar

Li, X., and Kleiven, S. (2018). Improved Safety Standards Are Needed to Better Protect Younger Children at Playgrounds. Sci. Rep. 8 (1), 15061. doi:10.1038/s41598-018-33393-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Sandler, H., and Kleiven, S. (2019). Infant Skull Fractures: Accident or Abuse? Forensic Sci. Int. 294, 173–182. doi:10.1016/j.forsciint.2018.11.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Sandler, H., and Kleiven, S. (2017). The Importance of Nonlinear Tissue Modelling in Finite Element Simulations of Infant Head Impacts. Biomech. Model. Mechanobiol 16 (3), 823–840. doi:10.1007/s10237-016-0855-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Zhou, Z., and Kleiven, S. (2021). An Anatomically Detailed and Personalizable Head Injury Model: Significance of Brain and white Matter Tract Morphological Variability on Strain. Biomech. Model. Mechanobiol. 20 (2), 403–431. doi:10.1007/s10237-020-01391-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Hu, J., Reed, M. P., Rupp, J. D., Hoff, C. N., Zhang, J., et al. (2011). Development, Validation, and Application of a Parametric Pediatric Head Finite Element Model for Impact Simulations. Ann. Biomed. Eng. 39 (12), 2984–2997. doi:10.1007/s10439-011-0409-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S., Beillas, P., Ding, L., and Wang, X. (2020). Morphing an Existing Open Source Human Body Model into a Personalized Model for Seating Discomfort Investigation. SAE Int. 2020, 01–0874. doi:10.4271/2020-01-0874

CrossRef Full Text | Google Scholar

Madhukar, A., and Ostoja-Starzewski, M. (2019). Finite Element Methods in Human Head Impact Simulations: A Review. Ann. Biomed. Eng. 47 (9), 1832–1854. doi:10.1007/s10439-019-02205-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Mao, H., Zhang, L., Jiang, B., Genthikatti, V. V., Jin, X., Zhu, F., et al. (2013). Development of a Finite Element Human Head Model Partially Validated with Thirty Five Experimental Cases. J. Biomech. Eng. 135 (11), 111002. doi:10.1115/1.4025101

CrossRef Full Text | Google Scholar

Mao, H., Gao, H., Cao, L., Genthikatti, V. V., and Yang, K. H. (2013). Development of High-Quality Hexahedral Human Brain Meshes Using Feature-Based Multi-Block Approach. Comp. Methods Biomech. Biomed. Eng. 16 (3), 271–279. doi:10.1080/10255842.2011.617005

CrossRef Full Text | Google Scholar

McKee, A. C., Stein, T. D., Kiernan, P. T., and Alvarez, V. E. (2015). The Neuropathology of Chronic Traumatic Encephalopathy. Brain Pathol. 25 (3), 350–364. doi:10.1111/bpa.12248

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, K., Wittek, A., Joldes, G., Horton, A., Dutta-Roy, T., Berger, J., et al. (2010). Modelling Brain Deformations for Computer-Integrated Neurosurgery. Int. J. Numer. Meth. Biomed. Engng. 26 (1), 117–138. doi:10.1002/cnm.1260

CrossRef Full Text | Google Scholar

Montanino, A., Li, X., Zhou, Z., Zeineh, M., Camarillo, D., and Kleiven, S. (2021). Subject-specific Multiscale Analysis of Concussion: from Macroscopic Loads to Molecular-Level Damage. Brain Multiphysics 2, 100027. doi:10.1016/j.brain.2021.100027

CrossRef Full Text | Google Scholar

Noël, L., and Kuhl, E. (2019). Modeling Neurodegeneration in Chronic Traumatic Encephalopathy Using Gradient Damage Models. Comput. Mech. 64 (5), 1375–1387. doi:10.1007/s00466-019-01717-z

CrossRef Full Text | Google Scholar

Opitz, A., Legon, W., Rowlands, A., Bickel, W. K., Paulus, W., and Tyler, W. J. (2013). Physiological Observations Validate Finite Element Models for Estimating Subject-specific Electric Field Distributions Induced by Transcranial Magnetic Stimulation of the Human Motor Cortex. Neuroimage 81, 253–264. doi:10.1016/j.neuroimage.2013.04.067

PubMed Abstract | CrossRef Full Text | Google Scholar

Opitz, A., Paulus, W., Will, S., Antunes, A., and Thielscher, A. (2015). Determinants of the Electric Field during Transcranial Direct Current Stimulation. Neuroimage 109, 140–150. doi:10.1016/j.neuroimage.2015.01.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Ou, Y., Akbari, H., Bilello, M., Da, X., and Davatzikos, C. (2014). Comparative Evaluation of Registration Algorithms in Different Brain Databases with Varying Difficulty: Results and Insights. IEEE Trans. Med. Imaging 33 (10), 2039–2065. doi:10.1109/tmi.2014.2330355

PubMed Abstract | CrossRef Full Text | Google Scholar

Ou, Y., Sotiras, A., Paragios, N., and Davatzikos, C. (2011). DRAMMS: Deformable Registration via Attribute Matching and Mutual-Saliency Weighting. Med. Image Anal. 15 (4), 622–639. doi:10.1016/j.media.2010.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, G., Kim, T., Forman, J., Panzer, M. B., and Crandall, J. R. (2017). Prediction of the Structural Response of the Femoral Shaft under Dynamic Loading Using Subject-specific Finite Element Models. Comp. Methods Biomech. Biomed. Eng. 20 (11), 1151–1166. doi:10.1080/10255842.2017.1340459

CrossRef Full Text | Google Scholar

Pedersen, K., Fahlstedt, M., Jacobsson, A., Kleiven, S., and von Holst, H. (2015). A National Survey of Traumatic Brain Injuries Admitted to Hospitals in Sweden from 1987 to 2010. Neuroepidemiology 45 (1), 20–27. doi:10.1159/000381780

PubMed Abstract | CrossRef Full Text | Google Scholar

Sahoo, D., Deck, C., and Willinger, R. (2014). Development and Validation of an Advanced Anisotropic Visco-Hyperelastic Human Brain FE Model. J. Mech. Behav. Biomed. Mater. 33, 24–42. doi:10.1016/j.jmbbm.2013.08.022

CrossRef Full Text | Google Scholar

Shepherd, J. F., and Johnson, C. R. (2008). Hexahedral Mesh Generation Constraints. Eng. Comput. 24 (3), 195–213. doi:10.1007/s00366-008-0091-4

CrossRef Full Text | Google Scholar

Shi, F., Yap, P.-T., Wu, G., Jia, H., Gilmore, J. H., Lin, W., et al. (2011). Infant Brain Atlases from Neonates to 1- and 2-Year-Olds. PLoS One 6 (4), e18746. doi:10.1371/journal.pone.0018746

PubMed Abstract | CrossRef Full Text | Google Scholar

Sigal, I. A., Hardisty, M. R., and Whyne, C. M. (2008). Mesh-morphing Algorithms for Specimen-specific Finite Element Modeling. J. Biomech. 41 (7), 1381–1389. doi:10.1016/j.jbiomech.2008.02.019

CrossRef Full Text | Google Scholar

Takhounts, E. G., Ridella, S. A., Hasija, V., Tannous, R. E., Campbell, J. Q., Malone, D., et al. (2008). Investigation of Traumatic Brain Injuries Using the Next Generation of Simulated Injury Monitor (SIMon) Finite Element Head Model. Stapp Car Crash J. 52, 1–31. doi:10.4271/2008-22-0001

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Essen, D. C., Smith, S. M., Barch, D. M., Behrens, T. E. J., Yacoub, E., Ugurbil, K., et al. (2013). The Wu-Minn Human Connectome Project: an Overview. Neuroimage 80, 62–79. doi:10.1016/j.neuroimage.2013.05.041

PubMed Abstract | CrossRef Full Text | Google Scholar

van Noort, R., Black, M. M., Martin, T. R. P., and Meanley, S. (1981). A Study of the Uniaxial Mechanical Properties of Human Dura Mater Preserved in Glycerol. Biomaterials 2 (1), 41–45. doi:10.1016/0142-9612(81)90086-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Vercauteren, T., Pennec, X., Perchant, A., and Ayache, N. (2009). Diffeomorphic Demons: Efficient Non-parametric Image Registration. Neuroimage 45 (1 Suppl. l), S61–S72. doi:10.1016/j.neuroimage.2008.10.040

PubMed Abstract | CrossRef Full Text | Google Scholar

von Holst, H., and Li, X. (2014). Decompressive Craniectomy (DC) at the Non-injured Side of the Brain Has the Potential to Improve Patient Outcome as Measured with Computational Simulation. Acta Neurochir 156 (10), 1961–1967. doi:10.1007/s00701-014-2195-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, T., Kleiven, S., and Li, X. (2020). Electroosmosis Based Novel Treatment Approach for Cerebral Edema. IEEE Trans Biomed Eng. PP.

Google Scholar

Weickenmeier, J., Jucker, M., Goriely, A., and Kuhl, E. (2019). A Physics-Based Model Explains the Prion-like Features of Neurodegeneration in Alzheimer's Disease, Parkinson's Disease, and Amyotrophic Lateral Sclerosis. J. Mech. Phys. Sol. 124, 264–281. doi:10.1016/j.jmps.2018.10.013

CrossRef Full Text | Google Scholar

Weickenmeier, J., Saez, P., Butler, C. A. M., Young, P. G., Goriely, A., and Kuhl, E. (2017). Bulging Brains. J. Elast 129 (1-2), 197–212. doi:10.1007/978-94-024-1220-8_10

CrossRef Full Text | Google Scholar

Windhoff, M., Opitz, A., and Thielscher, A. (2013). Electric Field Calculations in Brain Stimulation Based on Finite Elements: An Optimized Processing Pipeline for the Generation and Usage of Accurate Individual Head Models. Hum. Brain Mapp. 34 (4), 923–935. doi:10.1002/hbm.21479

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, J., Cai, M., Li, J., Cao, L., Xu, L., Li, N., et al. (2019). Development and Validation of a Semi-automatic Landmark Extraction Method for Mesh Morphing. Med. Eng. Phys. 70, 62–71. doi:10.1016/j.medengphy.2019.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, L., Yang, K. H., Dwarampudi, R., Omori, K., Li, T., Chang, K., et al. (2001). Recent Advances in Brain Injury Research: a New Human Head Model Development and Validation. Stapp Car Crash J. 45, 369–394. doi:10.4271/2001-22-0017

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z., Li, X., and Kleiven, S. (2019). Biomechanics of Acute Subdural Hematoma in the Elderly: A Fluid-Structure Interaction Study. J. Neurotrauma 36 (13), 2099–2108. doi:10.1089/neu.2018.6143

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Z., Li, X., and Kleiven, S. (2020). Evaluation of Brain-Skull Interface Modelling Approaches on the Prediction of Acute Subdural Hematoma in the Elderly. J. Biomech. 105, 109787. doi:10.1016/j.jbiomech.2020.109787

CrossRef Full Text | Google Scholar

Zöllei, L., Iglesias, J. E., Ou, Y., Grant, P. E., and Fischl, B. (2020). Infant FreeSurfer: An Automated Segmentation and Surface Extraction Pipeline for T1-Weighted Neuroimaging Data of Infants 0-2 Years. Neuroimage 218, 116946. doi:10.1016/j.neuroimage.2020.116946

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: finite element modeling, personalized simulation, traumatic brain injury, brain stimulation, neuroimage registration, biomechanics

Citation: Li X (2021) Subject-Specific Head Model Generation by Mesh Morphing: A Personalization Framework and Its Applications. Front. Bioeng. Biotechnol. 9:706566. doi: 10.3389/fbioe.2021.706566

Received: 07 May 2021; Accepted: 03 September 2021;
Published: 18 October 2021.

Edited by:

Silvia Budday, University of Erlangen Nuremberg, Germany

Reviewed by:

Haojie Mao, Western University, Canada
J. Sebastian Giudice, University of Virginia, United States

Copyright © 2021 Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xiaogai Li, xiaogai@kth.se

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.