The Variation in 3D Face Shapes of Dutch Children for Mask Design

: The use of 3D anthropometric data of children’s heads and faces has great potential in the development of protective gear and medical products that need to provide a close ﬁt in order to function well. Given the lack of detailed data of this kind, the aim of this study is to map the size and shape variation of Dutch children’s heads and faces and investigate possible implications for the design of a ventilation mask. In this study, a dataset of heads and faces of 303 Dutch children aged six months to seven years consisting of traditional measurements and 3D scans were analysed. A principal component analysis (PCA) of facial measurements was performed to map the variation of the children’s face shapes. The ﬁrst principal component describes the overall size, whilst the second principal component captures the more width related variation of the face. After establishing a homology between the 3D scanned face shapes, a second principal component analysis was done on the point coordinates, revealing the most prominent variations in 3D shape within the sample.


Introduction
Anthropometric information is commonly used in the design and evaluation of numerous applications such as workplaces, tools, clothing and wearables. Designers look to translate key body dimensions to relevant product shape and size in order to achieve a good fit. Traditional 1D anthropometric measurements are not able to capture the shape of the human body. However, this morphological information is becoming more important in the development of products such as apparel, backpacks, orthoses and headwear that need to closely fit a part of the body.
In headwear, a good fit is often required in order for the product to function properly. This is especially important in protective gear or medical products as the fit of these products could have a direct impact on the wearer's health and safety [1][2][3][4]. Thus, anthropometric information of the head and face is necessary to describe the variation in size and shape in order to develop a product that fits the user. Previous research has demonstrated the benefit of using 3D anthropometric data in order to understand the morphological variation of the head and face as well as in the improvement of product fit. Various researchers have studied the anthropometric variation of heads and faces based on 3D head scan data [5][6][7] in order to develop new sizing systems [8][9][10], in order to generate representative models [8,[11][12][13] or to improve the fit of certain products [14][15][16][17].
Previous studies all focus on adults, but relatively few studies have been carried out on 3D scan data for children using similar forms of analysis. In addition, the availability of 3D anthropometric data of children's heads and faces for designers is limited. Anthropometric surveys targeting children are a necessity as it has already been determined that children cannot be considered small adults, and thus that one cannot simply Appl. Sci. 2021, 11, 6843 3 of 16 to 7 years were recruited through health centres, primary schools and the university. Of the total population, 17.8% (n = 54) children were considered to be of non-native Dutch origin. This was defined as when the country of origin of either one or both of the child's parents was not the Netherlands. The anthropometric survey was approved by the Human Research Ethics Committee of the Delft University of Technology. Informed consent was obtained from all participants involved in the study. Age categorisation was done according to ISO 15535 [41], which describes the general requirements for establishing anthropometric databases. Age groups were divided as follows: individual age for age group 1 is 0.50 to 1.49 years; for age group 2, it is 1.50 to 2.49; and, for age group 3, it is 2.50 to 3.49, etc. In the survey, both traditional anthropometric measurements as well as 3D image derived measurements were collected as described in Goto et al. [28] through a five-step procedure. First, the following traditional anthropometric head and face dimensions were recorded; head circumference, head length, head height, head width, and face width, as well as more general measurements such as stature and weight. Second, four 3D images captured from different directions per participant were collected using the 3dMD Face system (3dMD Ltd., London, UK). All children were scanned with a neutral face expression. These images were then combined into a 360 degree 3D image of the head. using Artec Studio 9 software (Artec Group, Luxembourg) and the remaining holes in the image were repaired in Geomagic Studio 2013 software (3D Systems, Rock Hill, SC, USA). After that, a total of nineteen landmarks were marked on each 3D face with 3dMD Vultus 2.1 software (3dMD Ltd., London, UK) as shown in Figure 1. The landmarks that were included in the survey could be identified on the 3D image without palpation. The definitions of these landmark locations can be found in [28]. Finally, eight facial dimensions were measured and extracted after the 3D images were aligned according to the Frankfort horizontal plane (Martin and Knussmann, 1988) using MATLAB TM software (The MathWorks, Inc., Natick, MA, USA). An overview of the full procedure is illustrated in Figure 2. Table 1 shows the summarized information of the data including mean and standard deviations for each dimension per age group, males and females combined.

Participants and Data Collection
A survey was conducted in order to collect anthropometric data of children's heads and faces [28,40]. A total of 302 Dutch children (128 females, 174 males) aged 6 months to 7 years were recruited through health centres, primary schools and the university. Of the total population, 17.8% (n = 54) children were considered to be of non-native Dutch origin. This was defined as when the country of origin of either one or both of the child's parents was not the Netherlands. The anthropometric survey was approved by the Human Research Ethics Committee (HREC) of the Delft University of Technology in December 2012. Because at that time the HREC did not provide an ethical approval code, the application was re-submitted on the 20th of July 2021 and approved, this time with the following approval number 1736. Informed consent was obtained from all participants involved in the study. Age categorisation was done according to ISO 15535 [41], which describes the general requirements for establishing anthropometric databases. Age groups were divided as follows: individual age for age group 1 is 0.50 to 1.49 years; for age group 2, it is 1.50 to 2.49; and, for age group 3, it is 2.50 to 3.49, etc.
In the survey, both traditional anthropometric measurements as well as 3D image derived measurements were collected as described in Goto et al. [28] through a five-step procedure. First, the following traditional anthropometric head and face dimensions were recorded; head circumference, head length, head height, head width, and face width, as well as more general measurements such as stature and weight. Second, four 3D images captured from different directions per participant were collected using the 3dMD Face system (3dMD Ltd., London, UK). All children were scanned with a neutral face expression. These images were then combined into a 360 degree 3D image of the head. using Artec Studio 9 software (Artec Group, Luxembourg) and the remaining holes in the image were repaired in Geomagic Studio 2013 software (3D Systems, Rock Hill, SC, USA). After that, a total of nineteen landmarks were marked on each 3D face with 3dMD Vultus 2.1 software (3dMD Ltd., London, UK) as shown in Figure 1. The landmarks that were included in the survey could be identified on the 3D image without palpation. The definitions of these landmark locations can be found in [28]. Finally, eight facial dimensions were measured and extracted after the 3D images were aligned according to the Frankfort horizontal plane (Martin and Knussmann, 1988) using MATLAB TM software (The Math-Works, Inc., Natick, MA, USA). An overview of the full procedure is illustrated in Figure  2. Table 1 shows the summarized information of the data including mean and standard deviations for each dimension per age group, males and females combined.      First, an independent t-test was performed to investigate similarities between gender for each dimension per age group. As opposed to anthropometric data of adults, male and female data of children can sometimes be combined depending on the application [20]. When considering a product that needs to fit a certain age range of children, the variability of relevant dimensions has to be taken into account in the design. Because the variability of dimensions of children's heads and faces due to age is greater than the variability due to gender, it is often valid to combine the data of the different genders [20]. A total of 21 dimensions served as input for the t-test. A multiple comparisons correction was applied by employing the Benjamini-Hochberg procedure [42] with a false discovery rate (FDR) set at 0.05. MATLAB TM was used for statistical analysis.

Measurement Based Analysis of Face Variation
The PCA was performed to find important factors that explain the variation of the children's faces (MATLAB TM ). The input variables for the PCA in this study were chosen based on a review of previous studies and through discussion in a panel consisting of four anthropometry experts and ergonomists. The input variables for the PCA were considered relevant to mask design and the selection was based on the method proposed in two previous studies by Zhuang et al. [31] and Amirav et al. [25]. Zhuang et al. [31] identified 10 dimensions that were considered related to respirator fit in a study focussing on respirators for the adult civilian workforce namely, minimum frontal breadth, face width, bigonial breadth, face length, inter-pupillary distance, head breadth, nose protrusion, nose breadth, nasal root breadth and subnasale-sellion length. Amirav et al. [25] considered two facial dimensions relevant in the development of aerosol masks for children, namely, the width of the mouth and sellion-pogonion length. To our knowledge, similar information about dimensions related to fit does not exist for ventilation masks for children. As a result, PCA was applied to the following 9 measurements that were chosen for this study; face length, sellion-pogonion length, nose bridge length, mouth width, nose tip protrusion, nasal root breadth, nose breadth, inter-pupillary distance and face width. The definition of each measurement can be found in Table 2. The principal component (PC) loadings were calculated (MATLAB TM ) to give insight into the influence of different dimensions on the variation. Table 2. Definitions of the nine selected measurements for the measurements based PCA.

Measurement Definition
Face length Straight-line distance between the sellion (s) and menton (me) landmarks Sellion-pogonion length Straight-line distance between the sellion (s) and pogonion (pg) landmarks Nose bridge length Straight-line distance between the sellion (s) to pronasale (prn) landmarks Mouth width Straight-line distance between the left and right chelion landmarks (ch-ch) Nose tip protrusion Straight-line distance between the the subnasale (sn) and the pronasale (prn) landmarks Nasal root breadth The horizontal breadth of the nosals root spanning from the left nasal root point to the right (nrp-nrp) Nose breadth Straight line distance between the left and right alare landmarks (al-al) Inter-pupillary distance The straight-line distance between the centre of the left and the centre of the right pupil Face width Horizontal breadth between the left and right zygion landmarks (zy-zy)

3D Shape Based Analysis of Face Variation
In order to analyse the shape variation of the 3D face images, a morphological correspondence between each individual face image needs to be realized. Meshes of individual scans contain a varying amount of 3D data points that are distributed in different ways, which makes it impossible to compare them directly. By creating a correspondence, meshes are converted in such a way that the data points are approximately uniformly distributed over the shape, creating meshes that have the same number of data points and connectivity (triangles), with comparable landmark locations. This is referred to as homologous meshes.
In this study, the so-called non-rigid template registration method [43] is used using Wrap 3.4 software (Russian3dscanner, Moscow, Russia). First, a high-quality individual scan was selected to serve as the template mesh. This scan was then processed to be topologically equivalent to a disc (no surface handles and a single boundary at the neck) and with a uniform distribution of vertices. The obtained template was then used to create these homologous meshes by deforming the template mesh towards each individual face scan (target scan) as is illustrated in Figure 3. The annotated landmarks were used to steer the registration by forcing an exact match. Finally, the homologous meshes were spatially aligned (limited to translation and rotation) via partial Generalized Procrustes Analysis [44].

Gender Comparison
The differences in dimensions between gender for each age group are shown in Table  3. Overall, the mean head and face measurements for the male participants were larger than the female participants (120 out of 133 comparisons (90%) throughout all the age categories). However, only some incidental significant differences could be observed (7/133). They were scattered throughout different age groups and related to different dimensions. Therefore, for the remainder of this study, we work with a combined dataset in terms of gender.  A PCA was then conducted on the 3D coordinates of the homological meshes in order to analyse and visualize the 3D shape variation of the children's faces. The principal component analysis was conducted on the x-, y-, and z-values of each vertex of the 3D mesh data using a custom script in Python. The average head mesh was processed in ParaView to extract the relevant facial area through clipping with two interactively positioned planes. The resulting facial region of interest of the average shape was then warped to each individual face through the homology, effectively resulting in corresponding regions of interest for all subjects in the database, forming the input for the facial PCA analysis. Face shape modes that represent the shape variation are visualized along the first 7 Principal Components (PC's).

Gender Comparison
The differences in dimensions between gender for each age group are shown in Table 3. Overall, the mean head and face measurements for the male participants were Appl. Sci. 2021, 11, 6843 7 of 16 larger than the female participants (120 out of 133 comparisons (90%) throughout all the age categories). However, only some incidental significant differences could be observed (7/133). They were scattered throughout different age groups and related to different dimensions. Therefore, for the remainder of this study, we work with a combined dataset in terms of gender.

Face Variation
A PCA was conducted with the following dimensions; face length, sellion-pogonion length, nose bridge length, mouth width, nose tip protrusion, nasal root breadth, nose breadth, inter-pupillary distance and face width. The plot in Figure 4 shows the cumulative variance explained of each component. The threshold of the cumulative variance was set at 90% with each individual PC explaining at least 5% of the variance, resulting in the five PC's that are presented in Table 4. The first five PC's account for 90.39% of the total variation of the sample. The PC loadings of the first principal component (PC1) are all positive and relatively high. This means that each dimension contributes considerably to PC1 and thus captures the overall size and shape of the face. PC1 explains 56.94% of the variation. The PC loadings of the length related dimensions are all negative for the second principal component. This means that 12.22% of the variation is more width related, varying from broader shorter faces with less protruded noses (high PC 2 loading) to longer narrow faces with more protruded noses (low PC 2 loading). The general loadings of PC 2 to 5 are relatively low, which indicates that the differences in face shape are smaller.

Face Shape Variation
A second PCA was conducted on the face area in order to investigate the shape variation of the children's faces. The PCA was conducted on the vertices of the mesh of the face area. The first seven principal component scores were extracted, which explained at least 90% of the variation of the face shapes of the sample (see Figure 5 for the cumulative proportion of the explained variance). The PC based face shape modes are visualised in Figure 6, together with the colour map, projected on the mean, which visualises the magnitude of displacement for the shape mode per PC.
The first PC describes the variation of the overall size of the face, changing from small to large and accounts for 72% of the total variation. The second PC is related to the width of the face, varying from narrow to wide faces. The third PC is more related to the width of the forehead and the length and angle of the chin whilst the fourth PC shows the variation of the shape of the forehead related to the shape of the jaw. The fifth component shows the variation of the shape of the jaw and the depth of the face and, for the sixth PC, the variation represents the ratio between the lower face height and the head height and variance in nasal tip protrusion. Finally, the seventh PC shows the variation in face and jaw shape and the nose and lip protrusion. The second to seventh PC respectively account for 6%, 4%, 3%, 3%, 2%, and 1% of the variability. Figure 7 shows the scatter plot of the sample along the first and second PC's including the different shape modes.

Face Shape Variation
A second PCA was conducted on the face area in order to investigate the shape variation of the children's faces. The PCA was conducted on the vertices of the mesh of the face area. The first seven principal component scores were extracted, which explained at least 90% of the variation of the face shapes of the sample (see Figure 5 for the cumulative proportion of the explained variance). The PC based face shape modes are visualised in Figure 6, together with the colour map, projected on the mean, which visualises the magnitude of displacement for the shape mode per PC.  The first PC describes the variation of the overall size of the face, changing from small to large and accounts for 72% of the total variation. The second PC is related to the width of the face, varying from narrow to wide faces. The third PC is more related to the width of the forehead and the length and angle of the chin whilst the fourth PC shows the variation of the shape of the forehead related to the shape of the jaw. The fifth component shows the variation of the shape of the jaw and the depth of the face and, for the sixth PC, the variation represents the ratio between the lower face height and the head height and variance in nasal tip protrusion. Finally, the seventh PC shows the variation in face and jaw shape and the nose and lip protrusion. The second to seventh PC respectively account for 6%, 4%, 3%, 3%, 2%, and 1% of the variability. Figure 7 shows the scatter plot of the sample along the first and second PC's including the different shape modes.

Discussion
3D anthropometric data of children's faces are necessary in order to develop head and face gear with a proper fit. However, detailed (3D) anthropometric data of children's heads and faces are still lacking, even though there is a clear demand from the industry [45]. The aim of this study was to map the variation of the size and shape of children's faces, with a focus on the area relevant to mask design, by a multivariate and shape-based approach. An independent t-test was conducted in order to compare differences between gender, and the results showed that only 20 out of the 252 comparisons were considered significantly different, and these differences appeared to be distributed among the dimensions. When considering anthropometric data of children of a certain age range for applications in product design, genders can often be combined because of the greater variability of different dimensions compared to the variability due to gender [20]. Therefore, the PCA was conducted without differentiating between gender.
In addition, the children's dataset was not classified in age groups for the PCA. When designing a product for children within a certain age range, it is often more appropriate to investigate the variability of relevant dimensions irrespective of age [24,46]. There is a large variation in size amongst children within the same age group, which results in an overlap between successive age groups for most of the body dimensions [28]. In addition, the mean values and standard deviations of these body dimensions increase with age, which indicates an increasing differentiation among age groups. The differences in body shape between children of the same age will become larger and larger [46]. Indeed, the scatterplot of the PC scores ( Figure 8) shows that colours representing each age group are scattered throughout the graph, which illustrates the overlap between different age groups for multiple face dimensions.
A PCA was conducted on a selection of relevant dimensions for mask design, in an attempt to describe the morphological trend in the dataset to be used in the development of a ventilation mask. The result of the PCA describes the morphological distribution of the children's faces over an age span of 0.5 to 7 years. The first principal component describes the variation in the overall size of the face. While the PC1 score increases, the overall size of the face also increases. The second component describes the width of the face. Children with a relatively high PC 2 score have short and broad faces, whereas children

Discussion
3D anthropometric data of children's faces are necessary in order to develop head and face gear with a proper fit. However, detailed (3D) anthropometric data of children's heads and faces are still lacking, even though there is a clear demand from the industry [45]. The aim of this study was to map the variation of the size and shape of children's faces, with a focus on the area relevant to mask design, by a multivariate and shape-based approach. An independent t-test was conducted in order to compare differences between gender, and the results showed that only 20 out of the 252 comparisons were considered significantly different, and these differences appeared to be distributed among the dimensions. When considering anthropometric data of children of a certain age range for applications in product design, genders can often be combined because of the greater variability of different dimensions compared to the variability due to gender [20]. Therefore, the PCA was conducted without differentiating between gender.
In addition, the children's dataset was not classified in age groups for the PCA. When designing a product for children within a certain age range, it is often more appropriate to investigate the variability of relevant dimensions irrespective of age [24,46]. There is a large variation in size amongst children within the same age group, which results in an overlap between successive age groups for most of the body dimensions [28]. In addition, the mean values and standard deviations of these body dimensions increase with age, which indicates an increasing differentiation among age groups. The differences in body shape between children of the same age will become larger and larger [46]. Indeed, the scatterplot of the PC scores ( Figure 8) shows that colours representing each age group are scattered throughout the graph, which illustrates the overlap between different age groups for multiple face dimensions.
A PCA was conducted on a selection of relevant dimensions for mask design, in an attempt to describe the morphological trend in the dataset to be used in the development of a ventilation mask. The result of the PCA describes the morphological distribution of the children's faces over an age span of 0.5 to 7 years. The first principal component describes the variation in the overall size of the face. While the PC1 score increases, the overall size of the face also increases. The second component describes the width of the face. Children with a relatively high PC 2 score have short and broad faces, whereas children with a low PC 2 score have longer and narrower faces. The findings of this study are in agreement with the study of Seo et al. [26] who studied facial dimensions of Korean children for respirator design, where face length and face width were also found to influence the variation the strongest. Interestingly, Zhuang et al. [31] use a similar description for the size categories of adults based on a PCA as part of his study for the development of a new fit test panel for respirators-namely, small, medium, large, short-wide and long-narrow. This suggests that the overall characteristics that describe the face variation are similar for adults; however, there is one clear difference. The relatively high scores for PC 1 in this dataset (0.623-0.907) compared to the PC 1 scores for adults (0.194-0.426) indicate that the overall size of the face varies more in children. This can be due to age differences, resulting in a relatively larger difference between the smallest face and the largest face in the children's dataset.
A shape based PCA offers an even more detailed way of investigating face shapes. By conducting a PCA on the 3D location of the vertices of each 3D face area, the face shape variation of the dataset could be revealed. In this study, the first seven PC's were selected in order to investigate the shape variation. When determining the number of relevant principle components, a selection of components needs to be made that explain a cumulative percentage of the population. In this study, these seven PC's explain 91% of the total variation and, for most applications, a threshold between 70% and 90% is considered sufficient [32,47]. The shape variation along each PC is visualized through the shape modes that were generated for the average, +3 and −3 standard deviations. Similar to the measurement based PCA, the first PC describes the variation in the overall size of the face and the second PC describes the variation in the width of the face that resembles the results of the shape based PCA conducted by Zhuang et al. [37] and Luximon et al. [5] for their study of face shape variations of U.S. civilian workers and Chinese adults, respectively. The shape variations in the remaining PC's are more subtle and contribute less to the overall variation. Apart from the shape modes that are previously presented in Figures 6 and 7, and in order to illustrate the face variation of real participants, Figure 8 shows the measurementbased PC scatterplot with a preliminary selection of 3D face scans of eight participants. A selection was made of five faces that were distributed among the PC 1-axis, including a small (A), close to average (F) and a large (I) face. Secondly, for B and C, two additional A shape based PCA offers an even more detailed way of investigating face shapes. By conducting a PCA on the 3D location of the vertices of each 3D face area, the face shape variation of the dataset could be revealed. In this study, the first seven PC's were selected in order to investigate the shape variation. When determining the number of relevant principle components, a selection of components needs to be made that explain a cumulative percentage of the population. In this study, these seven PC's explain 91% of the total variation and, for most applications, a threshold between 70% and 90% is considered sufficient [32,47]. The shape variation along each PC is visualized through the shape modes that were generated for the average, +3 and −3 standard deviations. Similar to the measurement based PCA, the first PC describes the variation in the overall size of the face and the second PC describes the variation in the width of the face that resembles the results of the shape based PCA conducted by Zhuang et al. [37] and Luximon et al. [5] for their study of face shape variations of U.S. civilian workers and Chinese adults, respectively. The shape variations in the remaining PC's are more subtle and contribute less to the overall variation.
Apart from the shape modes that are previously presented in Figures 6 and 7, and in order to illustrate the face variation of real participants, Figure 8 shows the measurementbased PC scatterplot with a preliminary selection of 3D face scans of eight participants. A selection was made of five faces that were distributed among the PC 1-axis, including a small (A), close to average (F) and a large (I) face. Secondly, for B and C, two additional faces were selected to illustrate the differences along the PC2 axis. Next, to illustrate the variation in shape and the possible implications for the design of a ventilation mask, contours were projected on each of these faces (Rhinoceros ® ). Each contour passes the sellion, pogonion and approximately 10 mm distance from the left and right cheillion landmarks. They were then aligned at the sellion landmark. These points were selected because they could represent the preliminary contour of the rim of the ventilation mask. This illustration shows that the contours indeed not only vary in size but also in shape ( Figure 9). When analysing the contours from the front, we observe that, as the contour size increases, the shape varies but this does not necessarily scale proportionally. From the side, we observe that, as the contour increases in height, the depth of the contour also increases. This suggests that, when developing ventilation masks for children, one should take into account the facial characteristics as well as the size in order to achieve a good fit-for instance, by using a parametric design that is adjusted to different face shapes rather than simply scaling a product to different sizes.
faces were selected to illustrate the differences along the PC2 axis. Next, to illustrate the variation in shape and the possible implications for the design of a ventilation mask, contours were projected on each of these faces (Rhinoceros ® ). Each contour passes the sellion, pogonion and approximately 10 mm distance from the left and right cheillion landmarks. They were then aligned at the sellion landmark. These points were selected because they could represent the preliminary contour of the rim of the ventilation mask. This illustration shows that the contours indeed not only vary in size but also in shape (Figure 9). When analysing the contours from the front, we observe that, as the contour size increases, the shape varies but this does not necessarily scale proportionally. From the side, we observe that, as the contour increases in height, the depth of the contour also increases. This suggests that, when developing ventilation masks for children, one should take into account the facial characteristics as well as the size in order to achieve a good fit-for instance, by using a parametric design that is adjusted to different face shapes rather than simply scaling a product to different sizes. Both the dimension and shape based PCA give insight into how the face shape of Dutch children of different age groups varies. Each method reveals this information concerning the variation in face shapes in a different way. The dimension based PCA gives insight into which product relevant dimension influences the variation of face shape and to which extent. This helps the designer to understand the relationship between the product relevant dimensions and the face shape variation of the children in the dataset. The shape based PCA helps the designer to understand shape variation in a more visual way and shows how the face-scapes vary. The shape modes that are generated based on the analysis can facilitate the design of face related products in a Computer Aided Design Both the dimension and shape based PCA give insight into how the face shape of Dutch children of different age groups varies. Each method reveals this information concerning the variation in face shapes in a different way. The dimension based PCA gives insight into which product relevant dimension influences the variation of face shape and to which extent. This helps the designer to understand the relationship between the product relevant dimensions and the face shape variation of the children in the dataset. The shape based PCA helps the designer to understand shape variation in a more visual way and shows how the face-scapes vary. The shape modes that are generated based on the analysis can facilitate the design of face related products in a Computer Aided Design (CAD) environment. For a designer of head and face related products, both analysis methods can be useful and are complementary because they facilitate in organising and presenting the complex 3D data.
However, in order for the designer to be able to conduct a PCA and get the most out of the data, they need advanced knowledge of statistics, the ability to apply custom algorithms for data processing and to translate the results of these analysis to generate shape modes and specific 3D modelling skills to be able to visualize this data. These are not necessarily all in a designer's repertoire. Nevertheless, the richness of 3D anthropometric data is an advantage when it comes to its versatility in applications, and it can therefore be applied in different phases of the product development process and for different purposes, from product design and sizing to the virtual evaluation of a product's fit using actual face scans of individuals. This shows that, despite the usefulness of 3D anthropometric data for product design, designers are faced with a dilemma. Either they have to go through the complex process of processing, analyzing, translating and visualizing data alone or with experts in order to utilize the richness that 3D Data can provide, or they have to rely on available tools whose functionality may not always align with the designers' objectives. Thus, there is a clear need to provide designers with an intuitive way to access and explore body shape variation related to their product design. One step has been made by presenting current data via the so-called Mannequin tool of the DINED platform, which is accessible through https://dined.io.tudelft.nl/en (accessed on 15 May 2021) [48]. This tool implements a regression model relating the 1D measurements to the 3D PCA scores and thus provides a means to generate and explore the variation of children's 3D head shapes.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to privacy.

Conflicts of Interest:
The authors declare no conflict of interest.