Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 19 May 2023
Sec. Biomechanics
This article is part of the Research Topic Understanding Age and Sex-Related Differences in the Biomechanics of Road Traffic Associated Injuries Through Population Diversity Analyses, Volume II View all 5 articles

Personalization of human body models and beyond via image registration

Xiaogai Li
Xiaogai Li1*Qiantailang Yuan&#x;Qiantailang Yuan1Natalia Lindgren&#x;Natalia Lindgren1Qi HuangQi Huang1Madelen FahlstedtMadelen Fahlstedt2Jonas 
sth,Jonas Östh3,4Bengt Pipkorn,Bengt Pipkorn4,5Lotta Jakobsson,Lotta Jakobsson3,4Svein KleivenSvein Kleiven1
  • 1Division of Neuronic Engineering, Department of Biomedical Engineering and Health Systems, KTH Royal Institute of Technology, Huddinge, Sweden
  • 2Mips AB, Täby, Sweden
  • 3Volvo Cars Safety Centre, Gothenburg, Sweden
  • 4Division of Vehicle Safety, Department of Mechanics and Maritime Sciences, Chalmers University of Technology, Gothenburg, Sweden
  • 5Autoliv Research, Vargarda, Sweden

Finite element human body models (HBMs) are becoming increasingly important numerical tools for traffic safety. Developing a validated and reliable HBM from the start requires integrated efforts and continues to be a challenging task. Mesh morphing is an efficient technique to generate personalized HBMs accounting for individual anatomy once a baseline model has been developed. This study presents a new image registration–based mesh morphing method to generate personalized HBMs. The method is demonstrated by morphing four baseline HBMs (SAFER, THUMS, and VIVA+ in both seated and standing postures) into ten subjects with varying heights, body mass indices (BMIs), and sex. The resulting personalized HBMs show comparable element quality to the baseline models. This method enables the comparison of HBMs by morphing them into the same subject, eliminating geometric differences. The method also shows superior geometry correction capabilities, which facilitates converting a seated HBM to a standing one, combined with additional positioning tools. Furthermore, this method can be extended to personalize other models, and the feasibility of morphing vehicle models has been illustrated. In conclusion, this new image registration–based mesh morphing method allows rapid and robust personalization of HBMs, facilitating personalized simulations.

1 Introduction

Finite element (FE) human body models (HBMs) are becoming increasingly important numerical tools in vehicle safety for understanding injury mechanisms and developing prevention strategies (Östh et al., 2015; Boyle et al., 2019; Hu et al., 2019; Jakobsson et al., 2019; Pipkorn et al., 2019; Boyle et al., 2020; Hwang et al., 2020; Tang et al., 2020; Grebonval et al., 2021; Leledakis et al., 2021; Larsson et al., 2022; Bohman et al., 2022; Booth et al., 2022; Corrales et al., 2022; Erlinger et al., 2022; Mishra et al., 2022; Piqueras et al., 2022; Putra et al., 2022; Östh et al., 2022). HBMs have advantages over crash test dummies, such as representing varying loading directions, different anthropometries and sex, and muscle tonus (Jakobsson et al., 2019). Examples of HBMs include the THUMS (Shigeta et al., 2009), GHBMC (Gayzik et al., 2012), VIVA+ (John et al., 2022), SAFER HBM (Pipkorn et al., 2021), and PIPER (Beillas et al., 2016).

Personalized HBMs are required in situations where the individual anatomy is important, such as in the reconstruction of post-mortem human subject (PMHS) impact for HBM validation or accident reconstructions to elucidate injury mechanisms. However, developing a validated and reliable HBM from the start is time-consuming and remains challenging. Mesh morphing is a more efficient alternative for generating personalized HBMs than generating one from the start, which involves moving the nodes of the baseline HBM into a subject by using a displacement field that reflects the anatomical differences between the two. A literature review by Hu et al. (2012) highlighted the capacity of thin-plate spline radial basis function (RBF) for personalizing HBMs to a diverse population. Vavalle et al. (2014) morphed a GHBMC 50th-percentile male model into a 95th-percentile male, while Schoell et al. (2015) morphed the same into a 65-year-old (65YO) male. A study by Jolivet et al. (2015) also demonstrated that mesh morphing can generate HBMs with reasonable mesh quality and geometry accuracy by properly defining geometry targets; Beillas and Berthet (2017) morphed the GHBMC of a 95th male and a 5th female to 52 subjects based on anthropometric data. A series of studies from the University of Michigan have developed methodologies for generating parametric HBMs for a diverse population based on statistical shape models (SSMs) (Shi et al., 2015; Hwang et al., 2016a; Hwang, Hallman et al., 2016b; Hu et al., 2016; Zhang et al., 2017; Hu et al., 2019; Hwang et al., 2020; Tang et al., 2020), and this method has been adopted by John et al. (2022) to morph a 50th-percentile VIVA+ female model from seated to standing postures and a 50th male.

In addition to the abovementioned work, the open-source PIPER software has been developed in a European PIPER project, enabling personalization and positioning of HBMs through kriging (Jolivet et al., 2015; Janak et al., 2018; Janak et al., 2021), an interpolation method to deform geometrical models based on a set of sources and the associated target control points. The software is accompanied by a 6YO child model that is scalable from a 1.5YO to 12YO through metadata files containing landmarks and control points. The PIPER software has successfully been used to personalize and position various HBMs such as GHBMC (Grebonval et al., 2021; Corrales et al., 2022; Erlinger et al., 2022), THUMS (Germanetti et al., 2020), and VIVA (Kleinbach et al., 2018), demonstrating the versatility of the PIPER software beyond positioning the paired PIPER child model (Giordano et al., 2017; Li and Kleiven, 2018).

The abovementioned studies applied RBF or kriging interpolation methods to personalize the HBMs, which according to Janak et al. (2021), leads to an identical interpolation function for the parameters used in these studies. Both RBF and kriging require the identification of landmarks on the skin and skeleton/joints of both the baseline HBM (as source) and the subject (as target), and the landmarks have to correspond between the two. Based on these landmarks, a displacement field is calculated to move the baseline HBM nodes into the target subject. However, identifying the landmarks can be time-consuming and requires manual effort, which has led to the development of methods to simplify the process (Wu et al., 2019). The PIPER software uses a kriging method and requires metadata files with landmarks to personalize and position HBMs (PIPER Software Framework and Application: User guide, 2017), and new metadata have to be developed for new HBMs.

Both the RBF and kriging morphing methods are computationally expensive due to the task of inverting a large matrix. The RBF can handle a few thousand landmarks that define body posture and external surface shape, but for more accurate results with hundreds of thousands of landmarks, it becomes computationally challenging. To address this issue, a regional mesh morphing strategy has been proposed to split the body into parts and perform the morphing sequentially (Zhang et al., 2017; John et al., 2022). Additionally, a workflow has been suggested to efficiently handle the computation by using iterative subsampling and spatial subdivision methods (Janak et al., 2021).

Our study introduces a new method to personalize HBMs using image registration. We aim to overcome previous challenges associated with landmark-based morphing and offer a rapid and effective way to personalize HBMs. Recent studies have shown that image registration–based mesh morphing can be effective for personalizing human brain models (Giudice et al., 2020; Li, 2021; Li et al., 2021). However, it remains to be decided whether this method can be extended to personalize whole-body HBMs, in particular with the inclusion of skeletons. Further investigation is required to explore the feasibility and potential benefits of using image registration–based mesh morphing for whole-body HBMs, which motivates this study. We also explored the application of the method for geometrical correction and evaluated its potential in personalizing vehicle models.

2 Methods

Section 2.1 describes an overview of how we applied the image registration–based morphing method to four baseline HBMs and morphed them into ten subjects to demonstrate its effectiveness. The details of the method are explained in Section 2.2, followed by a demonstration of the usage of the method in Section 2.3. Finally, evaluations of morphed models are presented, which include element quality and runnability tests (Section 2.4).

2.1 HBM personalization: baseline HBMs and subjects

The four baseline HBMs include two models of 50th-percentile males in the standing posture (SAFER HBM and THUMS V4.02) and two models of 50th-percentile females in the seated and standing postures (VIVA+) (illustrated in Figure 1 and listed in Table 1). The anthropometric data of the eight target subjects are in Table 2. In total, ten personalized models were generated, and some were morphed into the same subject (subj1) as shown in the morphing map (Figure 1 and Table 3). The subsequent sections provide further information on the baseline HBMs and body shapes of the target subjects.

FIGURE 1
www.frontiersin.org

FIGURE 1. HBM personalization morphing map. Four baseline HBMs were morphed into eight subjects with varying heights (ranging from 160 to 187 cm) and a BMI of up to 34. SAFER HBM was morphed into six subjects (subj1 to subj6), THUMS was morphed into one subject (subj1), and the VIVA+ female was morphed into its corresponding male for both seated (subj7) and standing (subj8) postures. Additionally, the VIVA+ female standing model was morphed into subj1.

TABLE 1
www.frontiersin.org

TABLE 1. Anthropometric data of the four baseline HBMs.

TABLE 2
www.frontiersin.org

TABLE 2. Anthropometric data for the eight subjects for personalizing baseline HBMs.

TABLE 3
www.frontiersin.org

TABLE 3. List of the ten personalized HBMs from four baseline.

2.1.1 Baseline HBMs

2.1.1.1 SAFER-pedestrian HBM

A pedestrian version of the SAFER HBM was used, which was developed by Lindgren et al. (2023), through positioning and morphing of the SAFER HBM v10 in a seated position (Pipkorn et al., 2021). The following steps were taken: first, the SAFER seated occupant was positioned to a standing posture using the software Oasys PRIMER (Oasys, Solihull, United Kingdom). This led to a standing postured model but with distorted buttocks and a shorter standing body height. This PRIMER-positioned model was then geometrically corrected by morphing it to a 50th-percentile male body shape (see Supplementary Appendix S1), leading to the baseline pedestrian SAFER HBM. This model was then personalized to target subjects. Further details of the development and validation of the pedestrian SAFER HBM are presented in the study by Lindgren et al. (2023).

2.1.1.2 THUMS 50th pedestrian male

THUMS V4.02 pedestrian model was downloaded from the free access data set (https://www.toyota.co.jp/thums/download). The original model’s arm was rotated slightly within the software Oasys PRIMER to create the baseline model, ensuring that the arms were in a similar position to that of the target subjects for further personalization in this study.

2.1.1.3 VIVA+ HBMs in seated and standing postures

The 50th-percentile females in seated and standing postures (labeled as 50F-seated and 50F-standing) were downloaded from an open-source data set (https://openvt.eu/fem/viva/vivaplus/-/tree/main/model) and used as the baseline models to generate personalized HBMs for their corresponding male subjects. Note that the male versions of the VIVA+ model (labeled as 50M-seated and 50M-stand) were also downloaded. These models were previously morphed from their female counterparts using an RBF method (John et al., 2022). The inclusion of the male versions of the HBMs served two purposes: 1) to reverse engineer the body shapes to which the female models were morphed and 2) to compare the morphed models with the current method in this study with the RBF approach (Section 3.5).

2.1.2 Body and skeleton shapes of target subjects

For subjects 1–6, the body shapes of the skin and skeleton were generated based on previously developed statistical skin and skeleton models of SMPL (Loper et al., 2015), SMPLX (Bogo et al., 2016), and OSSO (Keller et al., 2022). Briefly, the subject’s anthropometric information, such as height, weight, and sex, was used as input to generate a SMPLX skin surface of the subject. The SMPLX was then positioned to match the stance of the baseline HBM using the SMPLX built-in tools. The surface model was then converted to SMPL format to allow subsequent generation of the corresponding skeleton of the positioned skin using the OSSO algorithm (OSSO: Obtaining Skeletal Shape from Outside) presented by Keller et al. (2022).

For subjects 7 and 8, the skin and skeletons were reverse-engineered from the VIVA+ male models of the seated and standing positions, respectively, by extracting the skin and skeleton surfaces from the downloaded 50M-seated and 50M-stand HBMs.

2.2 Method for HBM personalization via image registration

The image registration–based morphing method consists of three modules: i) pre-processing, ii) image registration pipeline, and iii) post-processing. Pre-processing involves converting the shape of the human body, that includes the skin and skeleton, into binary images through voxelization. The pipeline uses Demons registration to obtain a displacement field that represents the anatomical differences between the baseline HBM and target subject, which is then used to morph the baseline into a personalized HBM. Post-processing assesses the accuracy of the personalization. Figure 2 illustrates the morphing method applied to personalize THUMS to subject 1. Further details on each module are supplied below.

FIGURE 2
www.frontiersin.org

FIGURE 2. Overview of the image registration–based mesh morphing method for HBM personalization, which included three modules: (i) pre-processing extracts body shapes from the baseline HBM (1a_stl_hbm) and subject (1b_stl_subj) and converts them into binary images (2a_img_hbm and 2b_img_subj); (ii) Demons registration is performed with 2a_img_hbm as the fixed image and 2b_img_subj as the moving image to obtain a displacement field (gdemons). This displacement field, defined on the fixed image space, is used to morph the nodes of the baseline HBM (0_baseline_hbm) into a personalized HBM (4_subject_hbm) that corresponds to the target subject. The image corresponding to the baseline HBM (2a_img_hbm) is then warped by the inverse of the displacement field (gdemons1), resulting in the warped image 2c_img_basewarped2subj, which corresponds to the personalized HBM (4_subject_hbm). Comparing image 2c_img_basewarped2subj with the golden truth image 2a_img_subj allows for quantifying personalization accuracy, as measured by 5_DICE.

2.2.1 Pre-processing to voxelize body shape and skeleton to a binary image

The surface model of the baseline HBM and subject, which includes the skin surface and skeleton, are converted into binary images by 1) determining the minimum and maximum of the xyz coordinates of the polygon mesh, 2) defining the resolution of the output image, 3) voxelizing the inside as 1 and 2 and the outside as 0, and 4) outputting the binary image with flesh image value 1, skeleton 2, and background 0. The voxelization step is done for both the flesh and skeleton using the “Convert model to segmentation node” module in Slicer 3D, which allows the easy conversion of surface models into binary images.

2.2.2 Demons image registration

In this study, a non-linear registration method, the Diffeomorphic Demons algorithm (Vercauteren et al., 2009) implemented in the open-source software 3D Slicer was used. This algorithm was chosen for its ability to handle large anatomical differences, as demonstrated in previous studies (von Holst et al., 2012; Li et al., 2013; von Holst and Li, 2013; Li, 2021; Li, Zhou, and Kleiven, 2021). The binary images of the baseline HBM and subject are first rigidly aligned, then Demons registration is performed, where the subject’s image serves as the moving image and the baseline image as the fixed image. The registration process calculates a displacement field that aligns the two images as accurately as possible. Note that a smooth factor of 2 is chosen for Demons registration which is shown to be efficient for all the tested cases and allows a smooth displacement field capturing the anatomical differences between the HBMs and subjects, while allowing to handle different internal organ shapes between the two.

2.2.3 Mesh morphing

The displacement field gdemons obtained from the abovementioned registration step is defined for every voxel in the fixed image. As the baseline HBM is in the same space as the fixed image, applying gdemons to the baseline HBM leads to a personalized HBM. This is achieved by morphing the nodes of the baseline mesh to new positions using the formula:

xi=Xi+ui,(1)

where Xi is the nodal coordinate of node i, ui is the linearly interpolated displacement vector at node n from gdemons, and xi is the updated nodal coordinate. The personalized HBM is formed by the morphed nodes and same element definitions as the baseline.

2.2.4 Post-process for personalization accuracy: DICE and HD95 distance

The registration accuracy is evaluated by calculating the DICE and 95th-percentile Hausdorff distance (HD95) between a warped image (imgwarped) and the subject’s image. The baseline image (imgbaseline) is warped via the inverse of displacement fields from each registration step (gdemons1) resulting in a warped image (imgwarped)

imgwarped=gdemons1imgbaseline.(2)

As imgwarped corresponds to the personalized model, therefore both metrics also reflect the personalization accuracy.

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).

DICEA,B=2ABA+B,(3)

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; a DICE value of 0 implies no overlap, whereas a DICE coefficient of 1 indicates perfect overlap between the warped and target images.

The Hausdorff distance is defined as

HDC,D=maxhC,D,hD,C,(4)

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

hC,D=maxcCmaxdDcd.(5)

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 0, where a lower value indicates better registration accuracy between the warped and target images. Note that when there is a substantial difference in the internal organ shape between the HBMs and subjects, a decreased registration accuracy is expected. This is because, to maintain valid element quality, Demons registration with a smooth factor of 2 is used to morph the baseline HBM to the subjects, accounting for the overall shape, which can result in lower accuracy.

2.3 Personalization pipelines and demonstration applications

Three typical personalization pipelines (types I, II, and III) are described, along with their applications in generating HBMs a6, b1, and d8.

Type I is a basic pipeline (illustrated in Figure 2) that works well for most cases and has been used to generate nine personalized models except a4. An example usage of Type I is demonstrated by personalizing the SAFER HBM to a subject with a high BMI of 34 (a6) in Section 2.3.1. A skin-only parametric pipeline was also performed to emphasize the importance of the skeleton during morphing.

Type II pipeline adds extra steps to Type I for aligning subjects with significant anatomical differences with the baseline HBM. These extra steps can be a global transformation, multiple Demons registrations, or multiple morph steps focusing on local regions. For example, Type II was used to generate a4 (subject 4 with a height of 160 cm, morphed from a baseline HBM pedestrian SAFER with a height of 178 cm). In this case, a global transformation was used to elongate the subject image by a factor calculated as 178/160, which was then used as the input for the subsequent Demons registration.

Type III “shielding” pipeline is used to intentionally prevent morphing parts of the HBM. The application of Type III is demonstrated by 1) the THUMS model, where the head is “shielded” from morphing, as detailed in Section 2.3.3, and 2) the SAFER model, as detailed in Supplementary Appendix S1.

Both demonstrations highlight the ability of this method for geometrical correction. To further illustrate the geometrical correction capacity, another demonstration is shown for correcting the PIPER 18YO head model, as detailed in Supplementary Appendix S3. We also demonstrate the application of the method for morphing vehicles. As the focus of this study is on HBM morphing, all parts related to vehicle morphing are presented in Supplementary Appendix S2.

2.3.1 Type I pipeline application: morphing HBM to high-BMI subject (a6) and skin-only

Figure 3 displays the process of morphing the baseline SAFER HBM to subject 6 with a BMI of 34, along with the comparison of a parametric pipeline with skin only, while excluding the skeleton. The personalized mesh shows an expansion in the belly region when compared to the baseline HBM (Figure 3E vs. Figure 3D). Overlaying the baseline displacement field (in green frame lines, Figure 3F) with that obtained from the parametric pipeline (Figure 3F) highlights the difference between them and the need for including the skeleton in the registration step, especially for target subjects with high BMIs.

FIGURE 3
www.frontiersin.org

FIGURE 3. Baseline SAFER HBM (A) morphed to subj6 with a BMI of 34 (B), with an overlay that demonstrates that the skeleton differs less than the flesh of a person with high BMI than the baseline (C). The mesh of the baseline (D) and personalized model (E) shows the expansion of the mesh at the belly which is deformed by the displacement field (F) obtained from the Demons registration with the skeleton included. A parametric pipeline is performed using body shape only without a skeleton, resulting in a morphed mesh (E2) by the displacement field (F2). The comparison of the two displacement fields (F and F2) highlights the difference between the two (green frame shows baseline pipeline and red shows parametric pipeline).

2.3.2 Type I pipeline application: morphing female to male (d8)

The VIVA+ female in a standing posture was morphed into a 50th male using Type I pipeline (as shown in Figure 4). The differences between the two before morphing are illustrated in Figure 4C and D, which become minimal after morphing (Figures 4E, F), resulting in a DICE score of 0.96 and HD95 of 5.7 mm for the flesh. The absolute distance errors for both the skin and skeleton are visualized in Figure 4G, showing the skin to be almost perfectly aligned except for the fingers with a larger difference. The average distance error for the skeleton is 2.39 mm (Figure 4G), with the largest error observed at the ribs and hip edge. Registration accuracy in terms of DICE and HD95 for all other subjects is presented in the Results section.

FIGURE 4
www.frontiersin.org

FIGURE 4. VIVA+ female in standing posture (A) was morphed to a subject (50th male) (B), overlaying together shows the difference before morphing in shape (C) and the voxelized image (D) with the former in color and the latter in gray. The image after morphing (F) is used for calculating DICE for personalization accuracy, resulting in DICE 0.96 and HD95 5.7 mm. To depict personalization accuracy, the personalized male subject from the female baseline (E1) is superimposed onto the baseline female (E2). And the absolute error of the nearest distance between the two surfaces was calculated and visualized for both the skin and skeleton surfaces (G).

2.3.3 Type III pipeline application: “shielding” and geometrical correction (b1)

The THUMS HBM has an unusually prominent chest, which is shown to be corrected through the current morphing method. The basic Type I pipeline used for morphing THUMS to subj1 resulted in an unusual head shape (Figure 2 4_subject_hbm) due to the topological difference between the baseline and subject. To resolve this, the Type III pipeline was used to “shield” the difference in the head, meaning the head was not to be modified during morphing. This involved an additional step of replacing the head surface (scalp and skull) of the subject (Figure 5B) with that of THUMS (Figure 5A). Demons registration was then performed between the baseline HBM and subject (Figure 5C), which resulted in minimal displacement in the head region as indicated by the regular and undistorted frame lines (Figure 5E). The Type III pipeline led to a final personalized model from THUMS to subject 1(b') with a corrected body and the prominent chest pushed down (Figure 5F). The arm was also slightly repositioned during the morphing as indicated by the red arrows (Figures 5C, E).

FIGURE 5
www.frontiersin.org

FIGURE 5. Baseline THUMS body and skeleton surfaces extracted (A) and installed to the head (scalp and skull) surfaces of the subject (B) to “shield” the head during morphing to geometrically morph other body part only with the head unmodified. Before morphing, a prominent chest is seen (C), and after personalization, the difference between the subject (in gray) and the personalized image (in color) is almost invisible (D). The displacement field shown in frame lines (E) illustrates the deformation of the chest, and the arm is slightly repositioned. (F) shows the final personalized HBM from THUMS using the Type III pipeline, which has a similar body to when using Type I pipeline (Figure 2, 4), but with an unmodified head.

2.4 Analysis of personalized models: element quality and runnability test

Element quality was analyzed in the personalized HBMs in the form of Jacobian and aspect ratio, two commonly used indices for assessing element quality as reviewed earlier by Burkhart et al. (2013). Additionally, the minimum length of the elements was evaluated, as it determines the critical time step in the explicit FE dynamics analysis.

As HBMs contain complex contacts, a realizable morphing method should also ensure the maintaining of contacts, which are checked. The runnability of the personalized models was tested by subjecting them to a pedestrian impact to the side by a generic vehicle model developed by Pipkorn et al. (2014) at a velocity of 40 km/h. For details of the HBMs regarding materials, and interaction between organs, the readers are referred to the original study presenting these HBMs as mentioned in Section 2.1.1. The details of the generic vehicle model are found in the study by Pipkorn et al. (2014). LS-Dyna single precision version 13.0 was used to simulate the contact between the HBM and vehicle with a coefficient of 0.2.

3 Results

Personalization accuracy is presented in Section 3.1, followed by an analysis of element quality in Section 3.2 and runnability test results in Section 3.3. The analysis of all three HBMs morphed to the same subject can be found in Section 3.4, and a comparison between the current image registration–based approach and RBF is presented in Section 3.5.

3.1 Personalized HBMs and accuracy

The ten personalized models are visualized in Figure 6. The voxelized images of the baseline HBM and subjects that were used in the registration process are shown in Figure 7. The images show a significant difference between the baseline HBM and subjects before morphing (Figure 7, left upper row), and after morphing, the difference is nearly imperceptible (Figure 7, left lower row), which indicates a high degree of personalization accuracy. This accuracy is further quantified with DICE and HD95 values (Figure 7, right plot). The mean DICE and HD95 for the ten personalized models are 0.94 and 10.63 mm, respectively.

FIGURE 6
www.frontiersin.org

FIGURE 6. Ten personalized models morphed from four baseline HBMs. Each box shows the baseline HBM and its personalized models. Box1 shows SAFER HBM morphed to six subjects (subj1 to subj6). Box2 shows THUMS morphed to subject 1. Box3 and Box4 show the VIVA+ female 50th morphed into a male version in standing and seated postures, respectively.

FIGURE 7
www.frontiersin.org

FIGURE 7. Overlay of voxelized images between baseline HBM and subjects (left upper row) with the gray image representing the subject (golden truth) and the colored image representing the baseline HBM. The same is shown after morphing where the colored images correspond to the personalized HBMs (left lower row). The right plot shows the DICE and HD95 values for the flesh part, with each dot representing one personalized model. It should be noted that the flesh part in this study refers to the total area minus the skeleton part, as indicated by the green image. The DICE values for the skeleton are not included since the values do not represent personalization accuracy due to topology differences between the baseline and HBM in the skeleton.

3.2 Element quality

The ten personalized HBMs have comparable element quality with their respective baseline model in terms of Jacobian, aspect ratio, and element length (Figure 8).

FIGURE 8
www.frontiersin.org

FIGURE 8. Histogram plot of the element quality of the ten personalized HBMs compared with its baseline, which includes Jacobian (first row), aspect ratio (second row), and element length (third row).

3.3 Runnability test of personalized models

Three personalized HBMs (a6, b1', and d8) were selected for the runnability test, representing one baseline each from SAFER, THUMS, and VIVA+ in the standing posture. The models were subjected to a side impact using a generic vehicle buck model, showing that the personalized models are directly runnable (Figure 9) and can withstand such deformation with impact force comparable to a car accident.

FIGURE 9
www.frontiersin.org

FIGURE 9. Runnability test of three personalized models simulating a side impact captured at 0, 60, 120, and 160 ms.

3.4 Different HBMs personalized to the same 50th male

The proposed method facilitates easy personalization of different HBMs into the same subject, thereby eliminating geometrical differences between the models. This is demonstrated by overlaying three personalized models (i.e., a1, b1, and d1) morphed to subj1 from three baseline HBMs: SAFER HBM, THUMS, and VIVA+ female stand (Figure 10). The figures show observable geometrical differences among the HBMs before morphing (Figure 10, left box), while all models conformed to the same subject body shape and skeleton after morphing (Figure 10, right box).

FIGURE 10
www.frontiersin.org

FIGURE 10. Overlay of three baseline HBMs before and after morphing to the same subject. The left box shows the three baseline HBMs before morphing (blue for SAFER, red for THUMS, and pink for VIVA+ standing). The right box shows the three same HBMs after morphing to subj1, which eliminates the geometrical differences among the HBMs.

3.5 Comparison to RBF morphing

A comparison was made between two personalized models generated in this study (c7 and d8) and models obtained with an RBF approach by John et al. (2022). The personalized models obtained using the current morphing method (represented by gray wire in Figure 11) are overlaid with the RBF models (represented by green in Figure 11), and the models were found to be very similar. The difference was quantified in Figures 8G, H (since subject 7 and 8 body shapes were reverse-engineered from the RBF models), with DICE values of 0.96 and 0.97, respectively, for the standing and seated models, and HD95 values of 5.72 mm for the standing model and 3.15 mm for the seated one.

FIGURE 11
www.frontiersin.org

FIGURE 11. VIVA+ male 50th standing morphed from female via image registration in this study is visualized in gray wire compared with the morphed one using the RBF method by John et al. (2022), which is shown in green. (A) Similar comparison is made for the seated model (B).

4 Discussion

This study introduces an image registration–based method for personalization of HBMs and evaluates its performance in terms of personalization accuracy, element quality, and runnability. The ten personalized models achieved a high level of personalization accuracy with a mean DICE of 0.94 and a mean HD95 of 10.63 mm. The element quality of the personalized HBMs was comparable to their respective baseline, and the personalized models were directly runnable, without or only requiring minimal manual repair to avoid potential intersecting contacts in some morphed models. The method also enables easy personalization of different HBMs into the same subject, eliminating the geometric differences. A comparison with the RBF approach showed that the personalized models obtained using the current method are similar to those obtained using the RBF approach.

4.1 Fast and landmark-free HBM personalization method using image registration

The image registration–based method proposed in this study is landmark-free and distinguishes itself from traditional methods like RBF and kriging. Rather than relying on landmarks, this method converts surface models into binary images and employs Demons registration to obtain a displacement field that captures the anatomical difference between the subject and baseline HBM. This displacement field is then used to morph the baseline HBM to a target subject. One of the significant advantages of this method is that the resulting dense displacement field from registration captures detailed anatomy, and the Demons registration is “forgiving,” leading to valid element quality in personalized HBMs, even if there are holes or gaps in the voxelized binary image. By contrast, the RBF method often only uses thousands of landmarks since it involves the computationally heavy task of inverting a large matrix.

In image registration–based mesh morphing methods, there is a trade-off between personalization accuracy and element quality, determined by the chosen registration algorithm through the resulting displacement field. The Demons registration chosen in this study achieves a balance between personalization accuracy and element quality, yielding personalized HBMs with comparable element quality to the baseline. Of the steps, only the preprocessing module requires manual work and is quick to complete, while all other steps are automated and take minutes to complete. The extracted surface models for the baseline HBM can be reused when personalizing to other target subjects, making it easy to personalize baseline HBMs to new subjects.

4.2 Quality of personalized HBMs: element quality and contact

HBMs often contain many (up to thousands in some HBMs) surface contacts defined between body parts, e.g., the skeleton, inner organs, and flesh, to maintain separation. Intersecting contact surfaces will lead to locking contacts, which in turn will lead to unrealistic responses. Morphing HBMs can introduce new intersections, which must be carefully evaluated to ensure accurate results. The current morphing method generally introduces none or very few new intersections in the contact surface, depending on the quality of the baseline HBMs. For example, in the case of d8 (morphing VIVA+ female to male), no new intersections were introduced, whereas for b1 (THUMS morphing to subj1), two new contacting element intersections were introduced, and for the SAFER baseline to subj1, eight new contacting element intersections were introduced (illustrated in Supplementary Figure S2). These new intersections could nevertheless be manually repaired easily by translating the nodal coordinates of the intersecting elements to separate them.

4.3 Pipeline subtypes (Types I, II, and III)

The basic pipeline (Type I) was effective for nine out of ten evaluated cases, which included subjects with significant differences from the baseline, such as morphing a female baseline model to one of male (c7 and d8) and high BMIs (a6). However, it failed when morphing SAFER HBM to subj4, resulting in a distorted foot. To solve this issue, Type II pipeline was used to first globally elongate the voxelized image of subj4 before Demons registration. Type II pipeline can also involve multiple steps of Demons registrations for better accuracy in local regions of interest. Adding more regions, however, increases the risk of the decreased element quality. Type III pipeline is similar to Type I but has a “shielding” function to prevent certain regions, such as the head and foot, from being morphed.

4.4 Geometrical correction capacity facilitates extreme positioning and morphing of head models

Our method also provides easy geometric correction, as demonstrated by the successful correction of the THUMS chest (Figure 5) and SAFER buttocks (Supplementary Figure S1). It also facilitates extreme positioning of HBMs without the need for re-meshing. For example, the current morphing method when combined with the positioning tool Oasys PRIMER converted the SAFER HBM from a seated to standing version without re-meshing (Supplementary Figure S1). Converting a seated HBM to a standing one is challenging and often requires re-meshing or manual repair of the mesh (Peres, 2018; John et al., 2022). In a previous study, positioning of PIPER from seated to standing resulted in poor element quality at the joints, requiring re-meshing (Peres, 2018). A similar challenge was encountered when converting the VIVA+ female seated model to standing using an RBF-based approach, which resulted in problematic elements around the joints requiring manual repair (John et al., 2022).

Our approach for converting the SAFER occupant to a pedestrian (Lindgren et al., 2023) eliminates the need for manual repair. The process involves positioning the SAFER occupant to a pedestrian using the software Oasys PRIMER, similar to that done earlier for the PIPER HBM (Peres, 2018), but instead of re-meshing the flesh and skin, we used image registration to correct the shape and length. This produces a pedestrian model with a similar element quality as the baseline. Our method also allows for easy correction of models at the component level, as demonstrated by the successful correction of an 18YO PIPER head model. Our markerless method can be easily applied to personalize other existing HBMs, and it has a particular strength in generating subject-specific models when imaging data is available, such as full or partial body MRI images of the spine (Booth et al., 2022). The same pipeline used for the PIPER head (Supplementary Figure S5) can be used by segmenting CT/MRI images into binary images first.

4.5 Feasibility for morphing vehicle models

The proposed method is not limited to morphing HBMs and can be used to morph other types of models. To demonstrate this, we morphed an SUV into a sedan and vice versa, using voxelized images of the body and wheels while treating other parts, such as the windshield, as a single entity. The resulting morphed vehicle model is acceptable (Supplementary Figure S4B) and shows the potential of the method. However, the accuracy of the rear windshield of the morphed sedan is not ideal, as the SUV roof is morphed into the rear window. To improve accuracy, the windshield and body can be voxelized separately. It is also important to note that cars have different interiors and body parts, and our method is only suitable for rough estimation and morphing a baseline vehicle into a new car with similar structures.

4.6 Surface models of body shape and skeleton

We used a new open-source data set to obtain outer body shape models based on anthropometric data of height and weight. The skeletons were then automatically embedded using the open-source OSSO algorithm developed by researchers in the computer vision field (Keller et al., 2022). This approach offers an alternative to existing methods, such as the one used by John et al. (2022), where skeletons are embedded into outer body shapes based on bony landmarks using algorithms developed by Reed and Ebert. (2013).

We observed penetrations in the assembled skeleton and skin in some subjects, particularly in the pelvis area for subjects with lower BMI. However, the current method is forgiving, as it did not cause element quality problems in all the evaluated cases. In fact, including the skeleton would more likely lead to poor element quality than skin-only. As an alternative, skin-only body shapes could be used during morphing if relevant. Our results show that morphing using skin-only body shapes leads to acceptable personalized models. In this case, HumanShape.org is a valuable database for obtaining skin-only body shapes. However, including a skeleton would allow for more accurate personalization, especially for subjects with higher BMIs.

4.7 Limitations and future work

The current method is efficient for personalization and can position certain body parts, such as the arms and trunk. However, its positioning capacity is limited, particularly in regions where the source and target images do not overlap, such as in the hands. As a result, subjects’ arms must be positioned similar to baseline HBMs before morphing. The accuracy of personalization depends on the image registration algorithm used. The currently used Demons registration achieves good overall accuracy in the personalized models. However, some areas, such as the skeleton, may not be as accurate, where the rib cage deviates from the subject. This is potentially problematic if conducting a study on rib fractures, as previous research has shown that rib cage shape has a significant impact on fracture risk (Larsson et al., 2022). To address this issue, local morphing of specific regions can be performed to improve accuracy. Finally, the currently chosen Demons registration leads to a balanced personalization accuracy and element quality, but alternative algorithms, such as DRAMMS (Ou et al., 2011), could be explored in future studies.

5 Conclusion

In conclusion, this study introduces a new and efficient landmark-free method for personalizing HBMs using image registration. This automated approach allows for rapid personalization of HBMs to new subjects, incorporating the body shape and skeletons. The resulting personalized HBMs have element quality comparable to the baseline. The method also has strength in geometrical correction and facilitates extreme positioning of HBMs when combined with other positioning software. The proposed method has been demonstrated to be an efficient and robust approach for personalizing a range of HBMs, with potential applications beyond this. It also has particular strengths in geometric correction.

Data availability statement

The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.

Author contributions

XL developed the morphing pipeline and performed morphing, analyzed the results, and wrote the manuscript; QY was responsible for all SMPL-related aspects; NL was responsible for fixing SAFER HBM pedestrian, positioning all baseline HBMs, and performing contact check and runnability simulations; MF contributed to PIPER morphing; QH contributed to vehicle morphing; BP, LJ, and JÖ provided occupant version of the SAFER HBM; JÖ and LJ contributed to contact check; XL, MF, BP, LJ, and SK together obtained Vinnova (no. 2019-03386) funding. All authors contributed to the article and approved the submitted version.

Funding

The authors acknowledge the fundings from Vinnova (no. 2019-03386) and the Swedish Research Council (no. 2020-04724 and 2020-04496) for providing partial financial support to this project. Open access funding was provided by KTH Royal Institute of Technology.

Acknowledgments

This work benefits greatly from the open-source HBMs (VIVA and THUMS), open statistical body and skeleton models (SMPL, SMPLX, and OSSO algorithm), and imaging registration algorithm implemented in open-source software 3D Slicer.

Conflict of interest

MF was employed by Mips AB. JO and LJ were employed by Volvo Cars Safety Centre. BP was employed by Autoliv Research.

The remaining authors declare 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, editors, and 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.

Supplementary material

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

References

Beillas, P., Giordano, C., Alvarez, V., Li, X., Ying, X., Chevalier, M. C., et al. (2016). “Development and performance of the PIPER scalable child human body models,” in 14th international conference on the protection of children in cars (Munich. Germany: Department of Automotive Engineering), 19.

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, S142–S147. doi:10.1080/15389588.2017.1307971

PubMed Abstract | CrossRef Full Text | Google Scholar

Bogo, F., Kanazawa, A., Lassner, C., Gehler, P., Romero, J., and Black, M. J. (2016). 'Keep it SMPL: Automatic estimation of 3D human pose and shape from a single image. Comput. Vis. Eccv. 9909, 561–578.

CrossRef Full Text | Google Scholar

Bohman, K., El-Mobader, S., and Jakobsson, L. (2022). 'Effects of restraint parameters using PIPER 6y in reclined seating during frontal impact. Traffic Inj. Prev. 23, S123–S129. doi:10.1080/15389588.2022.2125304

PubMed Abstract | CrossRef Full Text | Google Scholar

Booth, G. R., Roberts, C. W., Romani, S., Mergler, O., Zhang, H., Pipkorn, B., et al. (2022). Methodology to measure seat belt fit in relation to skeletal geometry using an upright open MRI. Traffic Inj. Prev. 23 (sup1), S199–S201.

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyle, K., Fanta, A., Reed, M. P., Fischer, K., Smith, A., Adler, A., et al. (2020). 'Restraint systems considering occupant diversity and pre-crash posture. Traffic Inj. Prev. 21, S31–S36. doi:10.1080/15389588.2021.1895989

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyle, K. J., Reed, M. P., Zaseck, L. W., and Hu, J. (2019). “A human modelling study on occupant kinematics in highly reclined seats during frontal crashes,” in Proceedings of the international research conference on the Biomechanics of impact (Florence, Italy: IRCOBI).

Google Scholar

Burkhart, T. A., Andrews, D. M., and Dunning, C. E. (2013). 'Finite element modeling mesh quality, energy balance and validation methods: A review with recommendations associated with the modeling of bone tissue. J. Biomechanics 46, 1477–1488. doi:10.1016/j.jbiomech.2013.03.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Corrales, M. A., Bolte, J., Malcolm, S., Pipkorn, B., and Cronin, D. S. (2022). 'Methodology to geometrically age human body models to average and subject-specific anthropometrics, demonstrated using a small stature female model assessed in a side impact. Comput. Methods Biomechanics Biomed. Eng. 2022, 1–12. doi:10.1080/10255842.2022.2112187

CrossRef Full Text | Google Scholar

Erlinger, N., Kofler, D., Heider, E., and Klug, C. (2022). “Effects of boundary conditions and posture on simulations with human body models of braking events,” in International research Council on the Biomechanics of injury 2022, Porto, Portugal, September 14–16 (IRCOBI, International Research Council on biomechanics on Injury), 628–652.

Google Scholar

Evans, A. C.Brain Development Cooperative Group (2006). The NIH MRI study of normal brain development. Neuroimage 30, 184–202. doi:10.1016/j.neuroimage.2005.09.068

PubMed Abstract | CrossRef Full Text | Google Scholar

Forman, J. L., Joodaki, H., Ali, F., Riley, P., Bollapragada, V., Lessley, D., et al. (2015). “Biofidelity corridors for whole-body pedestrian impact with a generic buck,” in IRCOBI conf, Lyon, France, September 9–11, 2015, 356–372.

Google Scholar

Gayzik, F. S., Moreno, D. P., Rhyne, A. C., and Stitzel, J. D. (2012). “Development of a full human body finite element model for blunt injury prediction utilizing a multi-modality medical imaging protocol,” in 12th international LS-DYNA user conference (Dearborn, MI, USA: LS-DYNA), 3–5.

Google Scholar

Germanetti, F. D. F., Cappellino, F., Puppini, R., Belingardi, G., and Scattina, A. (2020). “Effects of non-conventional driving posture on the passive safety of autonomous vehicles - a numerical study using the THUMS human model,” in IRCOBI conference (Online), 20–54.

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, 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, 2412–2424. doi:10.1007/s10439-020-02584-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Grebonval, C., Trosseille, X., Petit, P., Wang, X., and Beillas, P. (2021). 'Effects of seat pan and pelvis angles on the occupant response in a reclined position during a frontal crash. PloS One 16, e0257292. doi:10.1371/journal.pone.0257292

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, J., Fanta, A., Neal, M. .O., Reed, M. .P., and Wang, J. .T. (2016). “Vehicle crash simulations with morphed GHBMC human models of different stature, BMI, and age,” in 4th international digital human modeling conference, Montréal, Québec, Canada, June 15–17, 2016.

Google Scholar

Hu, J., Rupp, J., and Reed, M. (2012). Focusing on vulnerable populations in crashes: Recent advances in finite element human models for injury biomechanics research. J. Automot. Saf. Energy 3, 295–307.

Google Scholar

Hu, J., Zhang, K., Reed, M. P., Wang, J. T., Neal, M., and Lin, C. H. (2019). 'Frontal crash simulations using parametric human models representing a diverse population. Traffic Inj. Prev. 20, S97–S105. doi:10.1080/15389588.2019.1581926

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, S., and Yang, J. (2010). 'Optimization of a reversible hood for protecting a pedestrian's head during car collisions. Accid. Analysis Prev. 42, 1136–1143. doi:10.1016/j.aap.2009.12.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, E., Hallman, J., Klein, K., Rupp, J., Reed, M., and Hu, J. (2016b). Rapid development of diverse human body models for crash simulations through mesh morphing. SAE Technical Paper.

CrossRef Full Text | Google Scholar

Hwang, E., Hu, J., Chen, C., Klein, K. F., Miller, C. S., Reed, M. P., et al. (2016a). 'Development, evaluation, and sensitivity analysis of parametric finite element whole-body human models in side impacts. Stapp Car Crash J. 60, 473–508. doi:10.4271/2016-22-0014

PubMed Abstract | CrossRef Full Text | Google Scholar

Hwang, E., Hu, J., and Reed, M. P. (2020). 'Validating diverse human body models against side impact tests with post-mortem human subjects. J. Biomechanics 98, 109444. doi:10.1016/j.jbiomech.2019.109444

PubMed Abstract | CrossRef Full Text | Google Scholar

Jakobsson, L., Bohman, K., and Svanberg, B. (2019). “Occupant protection for AD cars–the paradigm shift in crash safety?,” in 26th international technical conference on the enhanced safety of vehicles (ESV): Technology: Enabling a safer TomorrowNational highway traffic safety administration (Eindhoven , Netherlands: National Highway Traffic Safety Administration).

Google Scholar

Janak, T., Lafon, Y., Petit, P., and Beillas, P. (2018). “Transformation smoothing to use after positioning of finite element human body models,” in Proceedings of the international research conference on the Biomechanics of impact (Athens, Greece: IRCOBI), 18–33.

Google Scholar

Janak, T., Lafon, Y., Petit, P., and Beillas, P. (2021). 'A method to use kriging with large sets of control points to morph finite element models of the human body. J. Biomechanical Eng. 143, 021013. doi:10.1115/1.4048575

CrossRef Full Text | Google Scholar

John, J., Klug, C., Kranjec, M., Svenning, E., and Iraeus, J. (2022). 'Hello, world! VIVA+: A human body model lineup to evaluate sex-differences in crash protection. Front. Bioeng. Biotechnol. 10, 918904. doi:10.3389/fbioe.2022.918904

PubMed Abstract | CrossRef Full Text | Google Scholar

Jolivet, E., Lafon, Y., Petit, P., and Beillas, P. (2015). 'Comparison of kriging and moving least square methods to change the geometry of human body models. Stapp Car Crash J. 59, 337–357. doi:10.4271/2015-22-0013

PubMed Abstract | CrossRef Full Text | Google Scholar

Keller, M., Zuffi, S., Black, M. J., and Pujades, S. (2022). “OSSO: Obtaining skeletal shape from outside,” in Proceedings of the IEEE/CVF conference on computer vision and pattern recognition, New Orleans, LA, June 19–24, 2022 (IEEE), 20492–20501.

CrossRef Full Text | Google Scholar

Kleinbach, C., Fehr, J., and Brolin, K. (2018). “Simulation of average female rear-end volunteer tests using the active ViVA OpenHBM,” in IRCOBI conference 2018, Athens, Greece, September 12–14, 2018 (IRCOBI), 18–34.

Google Scholar

Klug, C., Feist, F., Raffler, M., Sinz, W., Petit, P., James, E., et al. (2017). “Development of a procedure to compare kinematics of human body models for pedestrian simulations,” in International research Council on the Biomechanics of injury 2017 (IRCOBI, International Research Council on biomechanics on Injury), 508–530.

Google Scholar

Kramlich, T., Langwieder, K., Lang, D., and Hell, W. (2002). “Accident characteristics in car-to-pedestrian impacts,” in Proc. IRCOBI conf, Munich, Germany, September 18–20, 2002 (IRCOBI), 119–130.

Google Scholar

Larsson, K. .J., Iraeus, J., Pipkorn, B., and Holcombe, S. (2022b). Influence of individual ribcage shape variability on occupant rib fracture risk. Front. Bioeng. Biotechnol. 11, 1–18. doi:10.3389/fbioe.2023.1154272

CrossRef Full Text | Google Scholar

Larsson, K. J., Pipkorn, B., Iraeus, J., Forman, J., and Hu, J. (2022a). 'Evaluation of a diverse population of morphed human body models for prediction of vehicle occupant crash kinematics. Comput. Methods Biomechanics Biomed. Eng. 25, 1125–1155. doi:10.1080/10255842.2021.2003790

CrossRef Full Text | Google Scholar

Leledakis, A., Osth, J., Davidsson, J., and Jakobsson, L. (2021). 'The influence of car passengers' sitting postures in intersection crashes. Accid. Anal. Prev. 157, 106170. doi:10.1016/j.aap.2021.106170

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | 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, 15061. doi:10.1038/s41598-018-33393-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., von Holst, H., and Kleiven, S. (2013). Decompressive craniectomy causes a significant strain increase in axonal fiber tracts. J. Clin. Neurosci. 20, 509–513. doi:10.1016/j.jocn.2012.04.019

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. Biomechanics Model. Mechanobiol. 20, 403–431. doi:10.1007/s10237-020-01391-8

CrossRef Full Text | Google Scholar

Lindgren, N., Yuan, Q., Pipkorn, B., Kleiven, S., and Li, X. (2023). Development and Whole-Body Validation of Personalizable Female and Male Pedestrian SAFER Human Body Models (2023). arXiv

Google Scholar

Liu, Y., Wan, X., Xu, W., Shi, L., Bai, Z., and Wang, F. (2022). 'A novel approach to investigate effects of front-end structures on injury response of e-bike riders: Combining Monte Carlo sampling, automatic operation, and data mining. Accid. Analysis Prev. 168, 106599. doi:10.1016/j.aap.2022.106599

CrossRef Full Text | Google Scholar

Loper, M., Mahmood, N., Romero, J., Pons-Moll, G., and Black, M. J. (2015). 'SMPL: A skinned multi-person linear model. Acm Trans. Graph. 34, 1–16. doi:10.1145/2816795.2818013

CrossRef Full Text | Google Scholar

Mishra, E., Mroz, K., Pipkorn, B., and Lubbe, N. (2022). Effects of automated emergency braking and seatbelt pre-pretensioning on occupant injury risks in high-severity frontal crashes. Front. Future Transp. 3, 883951. doi:10.3389/ffutr.2022.883951

CrossRef Full Text | Google Scholar

Östh, J., Brolin, K., Davidsson, J., Pipkorn, B., Jakobsson, L., Törnvall, F., et al. (2015). “Muscle activation strategies in human body models for the development of integrated safety,” in Proceedings of the 24th international techical conference on the enhanced safety of vehicles (Gothenburg, Sweden: ESV), 8–11.

Google Scholar

Östh, J., Larsson, E., and Jakobsson, L. (2022). “Human body model muscle activation influence on crash response,” in Proceedings of the IRCOBI conference, Portugal, September 14–16, 2022 (IRCOBI).

Google Scholar

Ou, Y. M., 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, 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, 622–639. doi:10.1016/j.media.2010.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Peres, J. (2018). “A pedestrian version of the piper 6 year old child model,” in 2018 IRCOBI conference proceedings, Athens, Greece, September 12–14, 2018 (IRCOBI), 247–248.

Google Scholar

PIPER Software Framework and Application: User guide (2017). PIPER manual. PIPER Project Partners. Available at: http://doc.piper-project.org/framework/1.0.1/user-guide.pdf.

Google Scholar

Pipkorn, B., Forsberg, C., Takahashi, Y., Ikeda, M., Fredriksson, R., Svensson, C., et al. (2014). “'Development and component validation of a generic vehicle front buck for pedestrian impact evaluation,” in Proceedings of the international research Council on the Biomechanics of impact, Berlin, Germany, September 10–12, 2014 (IRCOBI), (Online), 718–729.

Google Scholar

Pipkorn, B., Iraeus, J., Björklund, M., Bunketorp, O., and Jakobsson, L. (2019). “Multi-scale validation of a rib fracture prediction method for human body models,” in Proceedings of the IRCOBI conference, Florence, Italy, September 11–13, 2019 (IRCOBI), 175–192.

Google Scholar

Pipkorn, B., Jonas, Ö., Brynskog, E., Larsson, E., Rydqvist, L., Iraeus, J., et al. (2021). “Validation of the SAFER human body model kinematics in far-side impacts,” in Proceedings of the IRCOBI conference, September 8–10, 2021 (IRCOBI), 444–476.

Google Scholar

Piqueras, A., Pipkorn, B., Iraeus, J., Maza-Frechin, M., and Lopez-Valdes, F. J. (2022). 'Assessment of in situ chest deflection of post mortem human subjects (PMHS) and personalized human body models (HBM) in nearside oblique impacts. Traffic Inj. Prev. 23, 181–186. doi:10.1080/15389588.2022.2036731

PubMed Abstract | CrossRef Full Text | Google Scholar

Putra, I. P. A., Iraeus, J., Sato, F., Svensson, M. Y., and Thomson, R. (2022). 'Finite element human body models with active reflexive muscles suitable for sex based whiplash injury prediction. Front. Bioeng. Biotechnol. 10, 968939. doi:10.3389/fbioe.2022.968939

PubMed Abstract | CrossRef Full Text | Google Scholar

Reed, M. P., and Ebert, S. M. (2013). “Elderly occupants: Posture, body shape, and belt fit,” in University of Michigan (Ann Arbor: Transportation Research Institute).

Google Scholar

Schoell, S. L., Weaver, A. A., Urban, J. E., Jones, D. A., Stitzel, J. D., Hwang, E., et al. (2015). 'Development and validation of an older occupant finite element model of a mid-sized male for investigation of age-related injury risk. Stapp Car Crash J. 59, 359–383. doi:10.4271/2015-22-0014

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, X., Cao, L., Reed, M. P., Rupp, J. D., and Hu, J. (2015). 'Effects of obesity on occupant responses in frontal crashes: A simulation analysis using human body models. Comput. methods biomechanics Biomed. Eng. 18, 1280–1292. doi:10.1080/10255842.2014.900544

PubMed Abstract | CrossRef Full Text | Google Scholar

Shigeta, K., Kitagawa, Y., and Yasuki, T. (2009). “Development of next generation human FE model capable of organ injury prediction,” in Proceedings of the 21st annual enhanced safety of vehicles (GERMANY: National Highway Traffic Safety Administration), 15–18.

Google Scholar

Tang, J., Zhou, Q., Nie, B., and Hu, J. (2020). Obesity effects on pedestrian lower extremity injuries in vehicle-to-pedestrian impacts: A numerical investigation using human body models. Traffic Inj. Prev. 21, 569–574. doi:10.1080/15389588.2020.1821195

PubMed Abstract | CrossRef Full Text | Google Scholar

Vavalle, N. A., Schoell, S. L., Weaver, A. A., Stitzel, J. D., and Gayzik, F. S. (2014). 'Application of radial basis function methods in the development of a 95th percentile male seated FEA model. Stapp Car Crash J. 58, 361–384. doi:10.4271/2014-22-0013

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, S61–S72. doi:10.1016/j.neuroimage.2008.10.040

PubMed Abstract | CrossRef Full Text | Google Scholar

von Holst, H., and Li, X. G. (2013). 'Quantification of stretching in the ventricular wall and corpus callosum and corticospinal tracts in hydrocephalus before and after ventriculoperitoneal shunt operation. J. Appl. Math. 2013, 1–10. doi:10.1155/2013/350359

CrossRef Full Text | Google Scholar

von Holst, H., Li, X., and Kleiven, S. (2012). 'Increased strain levels and water content in brain tissue after decompressive craniotomy. Acta Neurochir. 154, 1583–1593. doi:10.1007/s00701-012-1393-2

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, K., Cao, L. B., Fanta, A., Reed, M. P., Neal, M., Wang, J. T., et al. (2017). 'An automated method to morph finite element whole-body human models with a wide range of stature and body shape for both men and women. J. Biomechanics 60, 253–260. doi:10.1016/j.jbiomech.2017.06.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: finite element human body model, image registration, mesh morphing, personalized simulations, traffic safety

Citation: Li X, Yuan Q, Lindgren N, Huang Q, Fahlstedt M, Östh J, Pipkorn B, Jakobsson L and Kleiven S (2023) Personalization of human body models and beyond via image registration. Front. Bioeng. Biotechnol. 11:1169365. doi: 10.3389/fbioe.2023.1169365

Received: 19 February 2023; Accepted: 28 April 2023;
Published: 19 May 2023.

Edited by:

Daniel Nicolella, Southwest Research Institute (SwRI), United States

Reviewed by:

Fang Wang, Changsha University of Science and Technology, China
Amin Komeili, University of Guelph, Canada

Copyright © 2023 Li, Yuan, Lindgren, Huang, Fahlstedt, Östh, Pipkorn, Jakobsson and Kleiven. 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, eGlhb2dhaUBrdGguc2U=

These authors have contributed equally to this work

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.